《Frontiers in Genome Editing》：Evaluation of genetic variation in tumor suppressor miRNA encoding and their target genes in breast cancer; focus on miRNA interaction and expression analysis
编辑推荐：
本研究揭示了位于KRAS基因3′非翻译区（3′UTR）的两种单核苷酸多态性（SNP）——rs712和rs9266——在乳腺癌发生中的关键作用。它们通过破坏肿瘤抑制性微小RNA（miRNA， 包括hsa-let-7c和hsa-miR-181c）的结合，解除对癌基因KRAS的转录后抑制，导致其表达上调。这一调控失衡与患者的不良预后相关，凸显了KRAS及其调控miRNA作为乳腺癌潜在诊断、预后标志物和治疗靶点的重要意义。
引言：探索乳腺癌中的基因调控网络
尽管乳腺癌（Breast Cancer， BC）的治疗手段已取得显著进步，但其复杂的分子机制，特别是疾病进展和治疗抵抗背后的机理，仍需深入探索。表观遗传因素，如组蛋白修饰、DNA甲基化和非编码RNA，已成为影响疾病发展的重要角色。其中，微小RNA（microRNA， miRNA）作为一类重要的非编码RNA，在包括乳腺癌在内的多种癌症中已被广泛研究。miRNA通过其种子序列与信使RNA（mRNA）的3′非翻译区（3′ Untranslated Region， 3′UTR）进行碱基互补配对，从而调控靶基因的表达。miRNA根据其功能可分为肿瘤抑制性miRNA（tumor suppressor miRNAs， tsmiRs）和致癌miRNA（oncogenic miRNAs， oncomiRs）。本研究的核心目的在于评估肿瘤抑制性miRNA编码基因（hsa-let-7c、hsa-miR-34a和hsa-miR-145a）及其靶基因（KRAS、IGFBP6和IGF1R）3′UTR的遗传变异，并阐明重要变异对RNA二级结构、疾病结局和化疗反应的潜在功能影响。
材料与方法：从基因到功能的系统性研究
研究招募了208名经组织病理学确诊的女性乳腺癌患者（年龄20-85岁，平均55.01 ± 11.9岁），她们来自印度旁遮普邦马尔瓦地区。所有参与者均签署了知情同意书。研究通过文献综述，确定了三个关键的肿瘤抑制性miRNA（hsa-let-7c、hsa-miR-34a和hsa-miR-145a）及其对应靶基因（KRAS、IGF1R和IGFBP6）的完整序列作为基因型分析对象。使用酚-氯仿法从外周血样本中提取基因组DNA，通过聚合酶链式反应（PCR）扩增所选miRNA基因的全长序列以及靶基因的3′UTR区域，并采用桑格测序技术进行测序。基因型分布和等位基因频率被计算分析。
测序结果显示，仅在KRAS基因的3′UTR中发现了两个变异：rs712和rs9266，而在分析的miRNA基因及其他靶基因中未检测到变异。为阐明rs9266对miRNA结合的功能性后果，研究将野生型（Wild-Type， WT）和突变型KRAS 3′UTR序列克隆到pMIR报告载体中。在MCF-7细胞系中共转染WT或突变型KRAS构建体以及hsa-let-7c和hsa-miR-181c的模拟物，通过双荧光素酶报告基因实验在转染48小时后测定荧光素酶活性，评估变异对miRNA介导的转录后调控的影响。所有实验均设置生物学三次重复。
研究从13份乳腺癌活检样本中提取RNA，根据rs712和rs9266基因型（2例野生型、9例杂合型、2例纯合突变型）分层，定量分析了KRAS以及已知可结合KRAS 3′UTR的miRNA（hsa-let-7c和hsa-miR-181c）的表达水平。β-肌动蛋白（β-actin）和U6 snRNA分别作为mRNA和miRNA的内参基因，采用2-ΔΔCt方法计算相对表达量。同时，通过Western印迹（Western blotting）和免疫组织化学（Immunohistochemistry， IHC）对来自三种KRAS基因型（WT、HT、HM）的乳腺癌活检组织蛋白提取物进行KRAS蛋白表达的验证。IHC采用半定量评分法评估染色强度（0 = 无， 1 = 弱， 2 = 中等， 3 = 强）和阳性细胞百分比，计算组织学评分（H-score）。
生物信息学分析方面，使用RNAfold网络服务器对WT和变异序列进行KRAS 3′UTR的RNA二级结构预测。生存分析利用Kaplan-Meier绘图仪平台，从GEO数据库提取临床数据集，分析KRAS、hsa-let-7c和hsa-miR-181c表达与乳腺癌患者生存结局（无远处转移生存期、总生存期、无进展生存期、无复发生存期）的关联。此外，利用ROCplot.org评估了KRAS表达在区分化疗应答者和非应答者方面的预测价值。
结果：揭示变异的功能影响与临床意义
对208名患者的人口统计学和临床特征分析为后续研究提供了坚实的队列基础。基因型分布分析显示，rs712和rs9266两个变异在患者群体中广泛存在。携带rs712和rs9266野生型TT/TT基因型的患者在年轻患者组（诊断时<45岁）和老年患者组（诊断时≥45岁）中占比略低，而携带变异基因型（TG+GG / TC+CC）的频率在两组中均较高。不同年龄初产、绝经状态及肿瘤分级患者中的基因型分布也显示出类似趋势，变异基因型普遍占比较高。
双荧光素酶报告基因实验明确证实，rs9266变异显著破坏了肿瘤抑制性miRNA与KRAS 3′UTR的结合，尤其是对hsa-miR-181c的影响最为显著。hsa-miR-181c的存在显著抑制了野生型构建体的荧光素酶活性，证实了其对KRAS的直接靶向作用。相比之下，rs9266突变消除了这种抑制作用，导致突变构建体的荧光素酶活性急剧升高，表明KRAS的表达被强力解除抑制而上调。rs9266对hsa-let-7c结合的影响虽然相对温和，但仍具有统计学意义。这些发现确立了rs9266作为一个关键的功能性变异，破坏了肿瘤抑制性miRNA的调控，从而驱动致癌性KRAS的过表达。
表达分析揭示了基因型驱动的表达模式。KRAS的表达在rs712和rs9266变异携带者中显示出高度显著的、阶梯式的升高：野生型（WT）患者维持最低的KRAS表达，杂合型（HT）患者表现出急剧的、最高的表达升高，而纯合突变型（HM）患者的表达也有显著但略低于HT的升高。相反，肿瘤抑制性miRNA hsa-miR-181c和hsa-let-7c则表现出引人注目的逆向表达模式：野生型个体拥有强健的miRNA表达，这是维持对KRAS严格调控所必需的；而杂合型携带者表现出显著下降，纯合突变型患者则遭受了最严重的miRNA耗竭。这种深刻的负相关关系强调了关键调控轴线的破坏。
免疫组织化学和Western印迹分析有力地证实了KRAS蛋白表达的基因型依赖性调控。野生型基因型携带者的乳腺癌组织显示出微弱的KRAS染色和相应较低水平的蛋白表达，反映了最低的致癌激活。与此形成鲜明对比的是，杂合型基因型表现出最强烈的免疫染色和显著升高的KRAS蛋白丰度，反映了翻译水平的急剧上调。纯合突变型则显示出中等的染色强度和蛋白表达，进一步强化了这些变异对KRAS表达的梯度式影响。
RNA二级结构预测显示，rs712和rs9266变异导致KRAS 3′UTR的结构发生了改变。与野生型相比，变异序列的最小自由能（Minimum Free Energy， MFE）略低，表明稳定性可能略有增加，同时热力学集合的多样性降低，暗示RNA构象的可变性减少。这些结构变化可能影响miRNA的结合以及RNA-蛋白质相互作用。
基于RNA芯片数据的Kaplan-Meier生存曲线分析显示，高KRAS表达与乳腺癌患者的某些不良生存结局相关。尽管与无远处转移生存期、总生存期和无进展生存期的关联不显著，但与无复发生存期（Relapse-Free Survival， RFS）存在显著关联，高KRAS表达预示着更高的复发风险。
相反，基于RNA-seq数据的Kaplan-Meier生存分析揭示，hsa-let-7c和hsa-miR-181c的高表达与乳腺癌患者更好的总生存期显著相关。高水平的hsa-let-7c与改善的总生存期强相关，风险比小于1，表明其具有保护作用。同样，高hsa-miR-181c表达的患者也表现出显著更好的总生存期。这些结果巩固了hsa-let-7c和hsa-miR-181c作为关键保护性miRNA在乳腺癌中的作用。
受试者工作特征曲线（Receiver Operating Characteristic， ROC）分析评估了KRAS表达对乳腺癌患者化疗反应的预测潜力。对于2级肿瘤患者，ROC曲线显示出中等的曲线下面积（Area Under the Curve， AUC），且具有高度显著的统计学意义，表明KRAS表达在预测该亚组化疗反应方面具有一定的准确性。然而，对于3级肿瘤患者，KRAS表达未能显示出具有统计学意义的预测能力，提示其单独作为预测高级别肿瘤化疗反应的生物标志物可能不足。
讨论与结论：意义、局限与未来方向
本研究强调了KRAS基因3′UTR中的rs712和rs9266变异在乳腺癌发展中的重要作用。rs9266变异破坏了miRNA（尤其是hsa-miR-181c）介导的KRAS调控，导致KRAS表达增加。杂合型基因型显示出最高的KRAS表达，其次是纯合突变型，而野生型表达最低。KRAS表达与肿瘤抑制性miRNA hsa-miR-181c和hsa-let-7c呈负相关，突显了它们在癌变中的调节作用。生物信息学分析揭示了这些遗传变异引起的KRAS 3′UTR结构改变，这可能影响miRNA结合和调控动态。生存分析和ROC分析进一步强化了KRAS、hsa-let-7c和hsa-miR-181c的预后意义：KRAS过表达与不良结局相关，而hsa-let-7c和hsa-miR-181c水平升高则与更好的总生存期相关。这些发现确立了KRAS及其调控miRNA作为潜在的诊断、预后和治疗生物标志物。
研究的优势在于采用了整合遗传筛查、功能验证和表达分析的多模式设计。荧光素酶报告基因实验直接证实了miRNA-3′UTR相互作用，增强了结果的可信度。qRT-PCR、IHC和Western印迹结果的一致性进一步增强了结果的真实性。应用KM绘图仪和ROCplot等多种生存分析工具，增强了研究结果的临床转化相关性。
研究也存在一些局限性。首先，由于缺乏健康的对照人群，无法直接建立观察到的KRAS变异与乳腺癌风险之间的相关性。其次，表达结果的普适性受到相对较小的活检样本量的限制。第三，需要进一步研究明确这些变异对下游致癌信号通路（如MAPK/ERK和PI3K/AKT通路）的影响，以建立直接的机制联系。最后，需要在其他种族群体中进行研究，以验证这些发现在更广泛人群中的适用性。
未来研究可结合先进的转录组计算分析和深入的机制探究。例如，应用批量RNA-seq去卷积分析和加权基因共表达网络分析，以识别受变异影响的共表达模块和候选竞争性内源RNA。实验上，可以利用CRISPR/Cas9技术在相关乳腺癌细胞系中构建rs712/rs9266等位基因，结合miRNA扰动和下游磷酸化蛋白质组学分析，明确这些3′UTR变异是否导致MAPK/ERK和PI3K/AKT信号通路的持续激活并改变药物敏感性。此外，探索通过miRNA疗法、反义寡核苷酸或靶向KRAS稳定性的小分子调节剂来恢复miRNA介导的抑制，是值得关注的精准医学新方向。本研究为KRAS 3′UTR变异通过改变肿瘤抑制性miRNA的转录后调控参与乳腺癌病理生理过程提供了基础证据，为未来开发针对KRAS驱动型乳腺癌的治疗和预后策略铺平了道路。
