KRAS 3′UTR基因變異(rs712, rs9266)損害腫瘤抑制性miRNA結(jié)合并促進(jìn)乳腺癌發(fā)生:一項(xiàng)聚焦miRNA相互作用與表達(dá)分析的研究
《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′非翻譯區(qū)(3′UTR)的兩種單核苷酸多態(tài)性(SNP)——rs712和rs9266——在乳腺癌發(fā)生中的關(guān)鍵作用。它們通過破壞腫瘤抑制性微小RNA(miRNA, 包括hsa-let-7c和hsa-miR-181c)的結(jié)合,解除對癌基因KRAS的轉(zhuǎn)錄后抑制,導(dǎo)致其表達(dá)上調(diào)。這一調(diào)控失衡與患者的不良預(yù)后相關(guān),凸顯了KRAS及其調(diào)控miRNA作為乳腺癌潛在診斷、預(yù)后標(biāo)志物和治療靶點(diǎn)的重要意義。
引言:探索乳腺癌中的基因調(diào)控網(wǎng)絡(luò)
盡管乳腺癌(Breast Cancer, BC)的治療手段已取得顯著進(jìn)步,但其復(fù)雜的分子機(jī)制,特別是疾病進(jìn)展和治療抵抗背后的機(jī)理,仍需深入探索。表觀遺傳因素,如組蛋白修飾、DNA甲基化和非編碼RNA,已成為影響疾病發(fā)展的重要角色。其中,微小RNA(microRNA, miRNA)作為一類重要的非編碼RNA,在包括乳腺癌在內(nèi)的多種癌癥中已被廣泛研究。miRNA通過其種子序列與信使RNA(mRNA)的3′非翻譯區(qū)(3′ Untranslated Region, 3′UTR)進(jìn)行堿基互補(bǔ)配對,從而調(diào)控靶基因的表達(dá)。miRNA根據(jù)其功能可分為腫瘤抑制性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二級結(jié)構(gòu)、疾病結(jié)局和化療反應(yīng)的潛在功能影響。
材料與方法:從基因到功能的系統(tǒng)性研究
研究招募了208名經(jīng)組織病理學(xué)確診的女性乳腺癌患者(年齡20-85歲,平均55.01 ± 11.9歲),她們來自印度旁遮普邦馬爾瓦地區(qū)。所有參與者均簽署了知情同意書。研究通過文獻(xiàn)綜述,確定了三個(gè)關(guān)鍵的腫瘤抑制性miRNA(hsa-let-7c、hsa-miR-34a和hsa-miR-145a)及其對應(yīng)靶基因(KRAS、IGF1R和IGFBP6)的完整序列作為基因型分析對象。使用酚-氯仿法從外周血樣本中提取基因組DNA,通過聚合酶鏈?zhǔn)椒磻?yīng)(PCR)擴(kuò)增所選miRNA基因的全長序列以及靶基因的3′UTR區(qū)域,并采用桑格測序技術(shù)進(jìn)行測序;蛐头植己偷任换蝾l率被計(jì)算分析。
測序結(jié)果顯示,僅在KRAS基因的3′UTR中發(fā)現(xiàn)了兩個(gè)變異:rs712和rs9266,而在分析的miRNA基因及其他靶基因中未檢測到變異。為闡明rs9266對miRNA結(jié)合的功能性后果,研究將野生型(Wild-Type, WT)和突變型KRAS 3′UTR序列克隆到pMIR報(bào)告載體中。在MCF-7細(xì)胞系中共轉(zhuǎn)染W(wǎng)T或突變型KRAS構(gòu)建體以及hsa-let-7c和hsa-miR-181c的模擬物,通過雙熒光素酶報(bào)告基因?qū)嶒?yàn)在轉(zhuǎn)染48小時(shí)后測定熒光素酶活性,評估變異對miRNA介導(dǎo)的轉(zhuǎn)錄后調(diào)控的影響。所有實(shí)驗(yàn)均設(shè)置生物學(xué)三次重復(fù)。
研究從13份乳腺癌活檢樣本中提取RNA,根據(jù)rs712和rs9266基因型(2例野生型、9例雜合型、2例純合突變型)分層,定量分析了KRAS以及已知可結(jié)合KRAS 3′UTR的miRNA(hsa-let-7c和hsa-miR-181c)的表達(dá)水平。β-肌動蛋白(β-actin)和U6 snRNA分別作為mRNA和miRNA的內(nèi)參基因,采用2-ΔΔCt方法計(jì)算相對表達(dá)量。同時(shí),通過Western印跡(Western blotting)和免疫組織化學(xué)(Immunohistochemistry, IHC)對來自三種KRAS基因型(WT、HT、HM)的乳腺癌活檢組織蛋白提取物進(jìn)行KRAS蛋白表達(dá)的驗(yàn)證。IHC采用半定量評分法評估染色強(qiáng)度(0 = 無, 1 = 弱, 2 = 中等, 3 = 強(qiáng))和陽性細(xì)胞百分比,計(jì)算組織學(xué)評分(H-score)。
生物信息學(xué)分析方面,使用RNAfold網(wǎng)絡(luò)服務(wù)器對WT和變異序列進(jìn)行KRAS 3′UTR的RNA二級結(jié)構(gòu)預(yù)測。生存分析利用Kaplan-Meier繪圖儀平臺,從GEO數(shù)據(jù)庫提取臨床數(shù)據(jù)集,分析KRAS、hsa-let-7c和hsa-miR-181c表達(dá)與乳腺癌患者生存結(jié)局(無遠(yuǎn)處轉(zhuǎn)移生存期、總生存期、無進(jìn)展生存期、無復(fù)發(fā)生存期)的關(guān)聯(lián)。此外,利用ROCplot.org評估了KRAS表達(dá)在區(qū)分化療應(yīng)答者和非應(yīng)答者方面的預(yù)測價(jià)值。
結(jié)果:揭示變異的功能影響與臨床意義
對208名患者的人口統(tǒng)計(jì)學(xué)和臨床特征分析為后續(xù)研究提供了堅(jiān)實(shí)的隊(duì)列基礎(chǔ);蛐头植挤治鲲@示,rs712和rs9266兩個(gè)變異在患者群體中廣泛存在。攜帶rs712和rs9266野生型TT/TT基因型的患者在年輕患者組(診斷時(shí)<45歲)和老年患者組(診斷時(shí)≥45歲)中占比略低,而攜帶變異基因型(TG+GG / TC+CC)的頻率在兩組中均較高。不同年齡初產(chǎn)、絕經(jīng)狀態(tài)及腫瘤分級患者中的基因型分布也顯示出類似趨勢,變異基因型普遍占比較高。
雙熒光素酶報(bào)告基因?qū)嶒?yàn)明確證實(shí),rs9266變異顯著破壞了腫瘤抑制性miRNA與KRAS 3′UTR的結(jié)合,尤其是對hsa-miR-181c的影響最為顯著。hsa-miR-181c的存在顯著抑制了野生型構(gòu)建體的熒光素酶活性,證實(shí)了其對KRAS的直接靶向作用。相比之下,rs9266突變消除了這種抑制作用,導(dǎo)致突變構(gòu)建體的熒光素酶活性急劇升高,表明KRAS的表達(dá)被強(qiáng)力解除抑制而上調(diào)。rs9266對hsa-let-7c結(jié)合的影響雖然相對溫和,但仍具有統(tǒng)計(jì)學(xué)意義。這些發(fā)現(xiàn)確立了rs9266作為一個(gè)關(guān)鍵的功能性變異,破壞了腫瘤抑制性miRNA的調(diào)控,從而驅(qū)動致癌性KRAS的過表達(dá)。
表達(dá)分析揭示了基因型驅(qū)動的表達(dá)模式。KRAS的表達(dá)在rs712和rs9266變異攜帶者中顯示出高度顯著的、階梯式的升高:野生型(WT)患者維持最低的KRAS表達(dá),雜合型(HT)患者表現(xiàn)出急劇的、最高的表達(dá)升高,而純合突變型(HM)患者的表達(dá)也有顯著但略低于HT的升高。相反,腫瘤抑制性miRNA hsa-miR-181c和hsa-let-7c則表現(xiàn)出引人注目的逆向表達(dá)模式:野生型個(gè)體擁有強(qiáng)健的miRNA表達(dá),這是維持對KRAS嚴(yán)格調(diào)控所必需的;而雜合型攜帶者表現(xiàn)出顯著下降,純合突變型患者則遭受了最嚴(yán)重的miRNA耗竭。這種深刻的負(fù)相關(guān)關(guān)系強(qiáng)調(diào)了關(guān)鍵調(diào)控軸線的破壞。
免疫組織化學(xué)和Western印跡分析有力地證實(shí)了KRAS蛋白表達(dá)的基因型依賴性調(diào)控。野生型基因型攜帶者的乳腺癌組織顯示出微弱的KRAS染色和相應(yīng)較低水平的蛋白表達(dá),反映了最低的致癌激活。與此形成鮮明對比的是,雜合型基因型表現(xiàn)出最強(qiáng)烈的免疫染色和顯著升高的KRAS蛋白豐度,反映了翻譯水平的急劇上調(diào)。純合突變型則顯示出中等的染色強(qiáng)度和蛋白表達(dá),進(jìn)一步強(qiáng)化了這些變異對KRAS表達(dá)的梯度式影響。
RNA二級結(jié)構(gòu)預(yù)測顯示,rs712和rs9266變異導(dǎo)致KRAS 3′UTR的結(jié)構(gòu)發(fā)生了改變。與野生型相比,變異序列的最小自由能(Minimum Free Energy, MFE)略低,表明穩(wěn)定性可能略有增加,同時(shí)熱力學(xué)集合的多樣性降低,暗示RNA構(gòu)象的可變性減少。這些結(jié)構(gòu)變化可能影響miRNA的結(jié)合以及RNA-蛋白質(zhì)相互作用。
基于RNA芯片數(shù)據(jù)的Kaplan-Meier生存曲線分析顯示,高KRAS表達(dá)與乳腺癌患者的某些不良生存結(jié)局相關(guān)。盡管與無遠(yuǎn)處轉(zhuǎn)移生存期、總生存期和無進(jìn)展生存期的關(guān)聯(lián)不顯著,但與無復(fù)發(fā)生存期(Relapse-Free Survival, RFS)存在顯著關(guān)聯(lián),高KRAS表達(dá)預(yù)示著更高的復(fù)發(fā)風(fēng)險(xiǎn)。
相反,基于RNA-seq數(shù)據(jù)的Kaplan-Meier生存分析揭示,hsa-let-7c和hsa-miR-181c的高表達(dá)與乳腺癌患者更好的總生存期顯著相關(guān)。高水平的hsa-let-7c與改善的總生存期強(qiáng)相關(guān),風(fēng)險(xiǎn)比小于1,表明其具有保護(hù)作用。同樣,高h(yuǎn)sa-miR-181c表達(dá)的患者也表現(xiàn)出顯著更好的總生存期。這些結(jié)果鞏固了hsa-let-7c和hsa-miR-181c作為關(guān)鍵保護(hù)性miRNA在乳腺癌中的作用。
受試者工作特征曲線(Receiver Operating Characteristic, ROC)分析評估了KRAS表達(dá)對乳腺癌患者化療反應(yīng)的預(yù)測潛力。對于2級腫瘤患者,ROC曲線顯示出中等的曲線下面積(Area Under the Curve, AUC),且具有高度顯著的統(tǒng)計(jì)學(xué)意義,表明KRAS表達(dá)在預(yù)測該亞組化療反應(yīng)方面具有一定的準(zhǔn)確性。然而,對于3級腫瘤患者,KRAS表達(dá)未能顯示出具有統(tǒng)計(jì)學(xué)意義的預(yù)測能力,提示其單獨(dú)作為預(yù)測高級別腫瘤化療反應(yīng)的生物標(biāo)志物可能不足。
討論與結(jié)論:意義、局限與未來方向
本研究強(qiáng)調(diào)了KRAS基因3′UTR中的rs712和rs9266變異在乳腺癌發(fā)展中的重要作用。rs9266變異破壞了miRNA(尤其是hsa-miR-181c)介導(dǎo)的KRAS調(diào)控,導(dǎo)致KRAS表達(dá)增加。雜合型基因型顯示出最高的KRAS表達(dá),其次是純合突變型,而野生型表達(dá)最低。KRAS表達(dá)與腫瘤抑制性miRNA hsa-miR-181c和hsa-let-7c呈負(fù)相關(guān),突顯了它們在癌變中的調(diào)節(jié)作用。生物信息學(xué)分析揭示了這些遺傳變異引起的KRAS 3′UTR結(jié)構(gòu)改變,這可能影響miRNA結(jié)合和調(diào)控動態(tài)。生存分析和ROC分析進(jìn)一步強(qiáng)化了KRAS、hsa-let-7c和hsa-miR-181c的預(yù)后意義:KRAS過表達(dá)與不良結(jié)局相關(guān),而hsa-let-7c和hsa-miR-181c水平升高則與更好的總生存期相關(guān)。這些發(fā)現(xiàn)確立了KRAS及其調(diào)控miRNA作為潛在的診斷、預(yù)后和治療生物標(biāo)志物。
研究的優(yōu)勢在于采用了整合遺傳篩查、功能驗(yàn)證和表達(dá)分析的多模式設(shè)計(jì)。熒光素酶報(bào)告基因?qū)嶒?yàn)直接證實(shí)了miRNA-3′UTR相互作用,增強(qiáng)了結(jié)果的可信度。qRT-PCR、IHC和Western印跡結(jié)果的一致性進(jìn)一步增強(qiáng)了結(jié)果的真實(shí)性。應(yīng)用KM繪圖儀和ROCplot等多種生存分析工具,增強(qiáng)了研究結(jié)果的臨床轉(zhuǎn)化相關(guān)性。
研究也存在一些局限性。首先,由于缺乏健康的對照人群,無法直接建立觀察到的KRAS變異與乳腺癌風(fēng)險(xiǎn)之間的相關(guān)性。其次,表達(dá)結(jié)果的普適性受到相對較小的活檢樣本量的限制。第三,需要進(jìn)一步研究明確這些變異對下游致癌信號通路(如MAPK/ERK和PI3K/AKT通路)的影響,以建立直接的機(jī)制聯(lián)系。最后,需要在其他種族群體中進(jìn)行研究,以驗(yàn)證這些發(fā)現(xiàn)在更廣泛人群中的適用性。
未來研究可結(jié)合先進(jìn)的轉(zhuǎn)錄組計(jì)算分析和深入的機(jī)制探究。例如,應(yīng)用批量RNA-seq去卷積分析和加權(quán)基因共表達(dá)網(wǎng)絡(luò)分析,以識別受變異影響的共表達(dá)模塊和候選競爭性內(nèi)源RNA。實(shí)驗(yàn)上,可以利用CRISPR/Cas9技術(shù)在相關(guān)乳腺癌細(xì)胞系中構(gòu)建rs712/rs9266等位基因,結(jié)合miRNA擾動和下游磷酸化蛋白質(zhì)組學(xué)分析,明確這些3′UTR變異是否導(dǎo)致MAPK/ERK和PI3K/AKT信號通路的持續(xù)激活并改變藥物敏感性。此外,探索通過miRNA療法、反義寡核苷酸或靶向KRAS穩(wěn)定性的小分子調(diào)節(jié)劑來恢復(fù)miRNA介導(dǎo)的抑制,是值得關(guān)注的精準(zhǔn)醫(yī)學(xué)新方向。本研究為KRAS 3′UTR變異通過改變腫瘤抑制性miRNA的轉(zhuǎn)錄后調(diào)控參與乳腺癌病理生理過程提供了基礎(chǔ)證據(jù),為未來開發(fā)針對KRAS驅(qū)動型乳腺癌的治療和預(yù)后策略鋪平了道路。