近年来,长链非编码RNA(lncRNA)作为一类长度超过200个核苷酸的非蛋白编码转录本,已被鉴定为多种生物学过程的关键调控分子
[1]。Hox转录反义RNA(HOTAIR)作为首个被阐明功能机制的lncRNA,在乳腺癌、胰腺癌等多种恶性肿瘤中呈现异常表达,其过表达水平与肿瘤分期、远处转移及总体生存率呈显著相关
[2]。多项研究表明
[3-5],HOTAIR通过调控染色质开放性动态改变转移相关基因表达谱,其诱导的表观遗传状态改变具有可逆性特征,这为靶向干预提供了重要理论基础。HOTAIR可正向调控HER2等致癌基因的表达,同时抑制PTEN等抑癌基因的转录活性,并通过激活PI3K/AKT/mTOR和Wnt/β-catenin等关键致癌信号通路导致肿瘤靶向耐药并与TGF-β1,SIRT1和E-cadherin等蛋白相互作用促进肿瘤恶性进展
[6, 7]。敲除HOTAIR可抑制MCF-7细胞增殖、迁移及紫杉醇耐药性,并通过增强DNA损伤修复缺陷使肿瘤细胞对化疗敏感性提升
[1, 8, 9]。
目前针对HOTAIR的研究仍存在一些局限:首先,对肿瘤研究覆盖面不广,主要集中在消化系统肿瘤
[10]、妇科肿瘤
[11]、肺部肿瘤
[12]以及非肿瘤慢性疾病
[13]等领域,对乳腺癌相关研究不够深入。并且多关注其促转移机制,缺乏遗传变异与治疗敏感性的关联研究
[14];其次,人群中HOTAIR单核苷酸多态性研究较少,全基因组关联研究(GWAS)发现的rs920778多态性虽已证实与乳腺癌发生发展相关
[15]。但其对靶向治疗的影响尚未阐明;再者,HOTAIR调控网络存在显著种族差异,东亚人群尤其是中国人群特异性数据亟待补充。
基于此,本研究拟对中国皖北地区女性乳腺癌患者采用TaqMan探针法进行HOTAIR基因分型,探讨rs920778位点多态性与乳腺癌易感性及抗HER2靶向治疗敏感性的关系。通过整合表观遗传调控与遗传变异双维度,揭示HOTAIR异常表达的分子基础。同时本研究也是首次尝试将HOTAIR遗传变异与抗HER2治疗响应性相关联,以期弥补lncRNA多态性研究在精准治疗领域的不足;为制定适合中国区域人群的个体化诊疗方案提供分子流行病学证据;该研究将为建立基于lncRNA遗传变异的精准诊疗策略提供新思路和依据,特别是针对携带rs920778高危基因型患者的早期干预和靶向治疗方案优化具有较强的临床意义。现报道如下:
1 资料和方法
1.1 患者招募和数据收集
实验组纳入于2022年7月~2024年12月经病理确诊的原发性乳腺癌患者共287例。纳入标准:经空心针穿刺活检或手术标本病理证实为浸润性癌;初诊且未接受系统性治疗(包括化疗、靶向治疗、内分泌治疗及免疫治疗);年龄18~75岁;临床资料完整(包括影像学、病理学及随访数据);ECOG评分≤2分,预期生存期≥6个月;签署知情同意书并自愿参与研究。排除标准:合并其他恶性肿瘤病史(非黑色素瘤皮肤癌或宫颈原位癌除外);存在严重心、肝、肾功能不全(ALT/AST>3倍正常值上限,肌酐清除率<50 mL/min);妊娠或哺乳期女性;既往接受过同侧乳腺放疗或胸壁放疗;患有活动性自身免疫性疾病或长期使用免疫抑制剂;存在精神疾病或依从性差无法配合随访。所有病例均严格遵循美国国家综合癌症网络(NCCN)乳腺癌临床实践指南(2022年第3版)进行TNM分期及分子分型。
同期招募260例健康女性志愿者,通过本院健康体检中心纳入,其年龄分布与实验组匹配(采用频数匹配法,按5岁年龄分层,确保各层比例一致)。纳入标准:乳腺超声及钼靶检查BI-RADS分级≤3级;无恶性肿瘤个人史及家族史(一级亲属);未接受雌孕激素替代治疗;近5年内未参与其他干预性临床研究。排除标准:乳腺超声或钼靶检查BI-RADS分级>3级,或存在可疑恶性病变;恶性肿瘤个人史及一级亲属家族史(乳腺癌、卵巢癌等);当前或既往接受雌孕激素替代治疗(HRT)或避孕药(持续使用≥6个月);合并其他慢性疾病(如系统性红斑狼疮、类风湿关节炎)或长期服用免疫调节药物;近5年内参与过其他干预性临床研究;妊娠或哺乳期女性;基因检测已知携带BRCA1/2、PALB2等乳腺癌易感基因致病突变。
通过医院电子病历系统(EMR)及结构化病例报告表(CRF),系统收集患者的临床病理参数,包括人口统计学特征(年龄、吸烟史、产次、绝经状态)、肿瘤生物学特性(TNM分期、分子亚型、组织学分级、淋巴血管侵犯(LVI)、腋窝淋巴结转移状态(SLNB/ALND证实)及远处转移)。所有受试者在入组前均已签署书面知情同意书,本研究经蚌埠医科大学第一附属医院伦理审查委员会批准(伦理批号:[2021038])。
1.2 DNA提取和基因分型
从每位参与者获得约1.5 mL外周静脉血,并立即储存在-80 ℃下。使用Relax Gene血液DNA系统从这些样品中提取基因组DNA(天根生物科技有限公司)。根据制造商的说明,SNP rs920778使用TaqMan等位基因辨别试验进行分型(TaqMan SNP测定:C9162435_20,Applied Biosystems,Foster City,CA,USA)对SNP rs920778进行基因分型。
1.3 RNA提取和实时PCR
从287名乳腺癌患者中获得新鲜肿瘤组织和血样标本。将肿瘤组织研磨成小片,并在1mL Trizol试剂(Invitrogen,卡尔斯巴德,CA,USA)中裂解,并按照制造商的说明收集全部RNA。然后使用PrimeScriptTM RT试剂盒(Takara)将RNA转录成cDNA,并使用SYBR EXTaq试剂盒(Takara)进行实时PCR。HOTAIR的引物序列由上海生工生物工程股份有限公司(中国上海)设计与合成,如下:正向,5'-GGTAGAAA AAGCAACCACGAAGC-3';反向,5'-ACATAAACCT CTGTCTGTGAGTGCC-3'。GAPDH用作内部对照(正向,5'-GTCAACGGATTTGGTCTGTATT-3';反向,5'-AGTCTTCTGGGTGGCAGTGAT-3'。HOTAIR的相对表达使用2-ΔΔCt法。
1.4 系统治疗与疗效评估
在晚期HER2阳性乳腺癌队列中,所有患者均接受一线THP联合治疗方案:曲妥珠单抗(Trastuzumab)首剂8 mg/kg静脉输注,后续维持剂量6 mg/kg(每3周1次);帕妥珠单抗(Pertuzumab)首剂840 mg静脉输注,后续维持剂量420 mg (每3周1次);联合紫杉类药物(多西他赛75 mg/m²,紫杉醇80 mg/m²或白蛋白结合型紫杉醇260 mg/m²)静脉给药,21 d为1周期。治疗周期第1天进行基线评估,后续每2个治疗周期(42±3 d)采用增强CT/MRI进行RECIST v1.1标准评估,持续至疾病进展(PD)或完成至少3次疗效评估(中位评估次数4,范围3~7)。疗效分类标准为:完全缓解(CR)定义为所有靶病灶消失;部分缓解(PR)指靶病灶直径总和减少≥30%;疾病稳定( SD)为变化介于-30%至+20%;疾病进展(PD)指靶病灶直径总和增加≥20%或出现新病灶。客观缓解率(ORR)计算为(CR+PR)/总病例数×100%。
在早期HER2阳性新辅助治疗队列中,所有患者接受剂量密集型TCbHP方案:曲妥珠单抗(首剂8 mg/kg,维持6 mg/kg,1次/3周)、帕妥珠单抗(首剂840 mg,维持420 mg,1次/3周)、卡铂(AUC 6,1次/3周)联合紫杉类药物(多西他赛75 mg/m²,紫杉醇80 mg/m²或白蛋白结合型紫杉醇260 mg/m²)。每2个治疗周期(42±3 d)通过乳腺MRI进行肿瘤负荷动态监测,直至完成6周期治疗或达到手术指征。术后病理评估采用Miller-Payne分级系统:MP1(癌细胞数量无减少)、MP2(癌细胞数量减少≤30%)、MP3(癌细胞数量减少31%~90%)、MP4(癌细胞数量减少91%~99%,仅存极少量散在癌细胞)、MP5(无浸润癌残留,允许导管原位癌存在)。观察指标包括乳房病理完全缓解(bpCR,定义为ypT0/Tis ypN0/1)和总病理缓解率(tpCR,定义为ypT0/Tis ypN0)。
1.5 统计学分析
所有统计分析均采用SPSS Statistics 22.0 (IBM Corp., Armonk, NY,USA)和Haploview 4.2 (Broad Institute, Cambridge,MA,USA)软件平台完成。Pearson χ²检验比较乳腺癌组与健康对照组间基因型分布及等位基因频率的差异性以及rs920778不同基因型与乳腺癌临床病理特征(包括TNM分期、组织学分级、激素受体状态等)的相关性。多因素logistic回归模型计算比值比(OR)及其95%置信区间(CI)。所有统计检验均为双侧检验,P<0.05为差异具有统计学意义。
2 结果
2.1 乳腺癌组织中rs920778位点多态性与HOTAIR基因的表达的关系
在287名乳腺癌患者中,HOTAIR基因的rs920778位点的TT基因型(
n=60)相比于CC基因型(
n=128),可显著上调HOTAIR表达(
P<0.001)。CT基因型(
n=99)相比于CC基因型,并不能提高HOTAIR的表达(
P=0.330,
图1)。
2.2 乳腺癌的易感性与HOTAIR基因rs920778位点多态性的关系
共547名年龄匹配的女性(287名乳腺癌患者和260名健康对照),两组年龄分别为46.5±8.8岁和44.8±7.2岁,差异无统计学意义(P>0.05),并且在月经状态、婚姻状况、生育状况、吸烟史等方面差异均无统计学意义(P>0.05)。
两组中,HOTAIR基因rs920778多态性的基因型频率均未偏离哈代-温伯格平衡(HWE)(
P>0.05)。T等位基因在乳腺癌患者组中的频率为39.5%,健康对照组为29.7%,其与乳腺癌风险增加显著相关(OR=1.54,95%
CI 1.09~2.19,
P=0.017)。与对照组相比,在乳腺癌患者组中TT基因型的频率明显增高(20.7%
vs 10.0%)。当将TT和CT受试者结合在一起进行分析时,含有T等位基因的患者,罹患乳腺癌风险增加(TT+CT
vs CC:OR=1.83,95%
CI 1.26~2.54,
P<0.001,
表1)。
2.3 HOTAIR基因rs920778位点多态性与乳腺癌临床病理特征的关系
rs920778位点的TT基因型与乳腺癌临床分期偏晚(
P<0.001)、不良组织学分级(
P<0.001)、存在淋巴结转移(
P<0.001)和分子亚型(三阴性,
P<0.001)相关。但与患者诊断时的年龄(
P>0.05)、婚育状况(
P>0.05)、绝经状态(
P>0.05)或脉管浸润(
P>0.05)无关(
表2)。
2.4 HOTAIR基因rs920778位点多态性对HER2阳性晚期乳腺癌治疗疗效的影响
共纳入43例HER2阳性晚期乳腺癌,治疗后评估完全缓解率(CR):23.2%(10/43),部分缓解率(PR):46.5%(20/43),客观有效率(ORR)为69.8%(30/43)。拥有rs920778位点的TT基因型患者治疗疗效显著较差(
P<0.001,
表3)。
2.5 HOTAIR基因rs920778位点多态性对HER2阳性早期乳腺癌新辅助治疗疗效的影响
共纳入48例早期HER2阳性乳腺癌,符合新辅助治疗指证,完成新辅助治疗后全部进行手术并评估,bpCR率为58.3%(28/48),tpCR率为47.9%(23/48),拥有rs 920778位点的TT基因型患者经新辅助治疗后pCR率显著降低(
P=0.018),但bpCR率与tpCR率差异无统计学意义(
P=0.728,
表4)。
3 讨论
HOTAIR是一种位于人类12号染色体HOXC基因簇的长链非编码RNA(lncRNA),通过结合染色质修饰复合体PRC2和LSD1诱导组蛋白H3K27三甲基化(H3K27me3)及降低H3K4me2水平,从而沉默靶基因表达。在乳腺癌中,HOTAIR的高表达与肿瘤转移和不良预后密切相关
[16, 17]。本研究探讨了HOTAIR基因rs920778多态性在中国皖北地区乳腺癌人群中的遗传易感性、临床病理特征及治疗反应调控作用。通过整合基因分型、分子表达谱与临床疗效数据,我们发现TT基因型不仅是乳腺癌发生的独立危险因素,还与肿瘤恶性表型及HER2靶向治疗耐药显著相关。这一发现与我们团队既往研究相互印证
[18],为乳腺癌的分子分型及精准治疗提供了新的遗传学视角。
从分子机制层面,本研究在临床样本中证实rs920778的TT基因型与HOTAIR表达上调存在显著相关性(
P=0.012)。这一现象可能与SNP位点对转录因子结合能力的调控密切相关。rs920778位于HOTAIR基因第5内含子区(hg38 chr12:54,343,774),该区域包含多个转录因子结合基序。生物信息学分析显示,T等位基因可能破坏原有YY1结合位点,同时形成新的SP1结合域
[19]。体外荧光素酶报告基因实验证实,携带T等位基因的启动子片段具有更强的转录活性(约3.2倍)
[20]。此外,表观遗传学研究表明,该SNP可能通过改变染色质三维构象,增强HOTAIR启动子与远端增强子的空间互作
[21]。这些分子机制共同导致TT基因型个体HOTAIR的异常高表达,进而通过其致癌通路促进肿瘤发生
[22]。值得注意的是,HOTAIR过表达可抑制BRCA1介导的同源重组修复
[23][2],导致基因组不稳定性增加,这与三阴性乳腺癌(TNBC)的分子特征高度吻合,可能部分解释了TT基因型与TNBC亚型的强相关性。与Qian等
[24]研究结果一致。
在遗传易感性方面,本研究显示T等位基因携带者(TT+CT)的乳腺癌风险较CC基因型增加1.83倍(
P=0.008),其中TT纯合子的OR值高达2.37(95%
CI 1.53-3.67)。这一结果与东亚人群GWAS数据(OR=1.45)存在一定差异
[15],可能源于地域遗传背景的异质性。皖北地区人群的T等位基因频率(39.5%)高于北京(32.1%)和上海(28.6%)的报道值
[3, 4, 25, 26],提示该SNP可能作为区域特异性遗传标志物。多因素分析进一步显示,TT基因型与晚期肿瘤分期(70% III/IV期病例)、高组织学分级(80% III级)及淋巴结转移(80%)等侵袭性特征显著相关。这些发现与HOTAIR的已知功能高度一致:其可通过募集PRC2复合体介导H3K27me3修饰,沉默PTEN、DKK1等抑癌基因
[27];作为miR-34a分子海绵解除对MET、AXL等促转移基因的抑制
[28];并通过激活Wnt/β-catenin通路促进上皮-间质转化(EMT)
[5, 29]。值得注意的是,TT基因型与脉管侵犯无显著相关性(
P=0.844),提示HOTAIR可能更倾向于调控淋巴结转移的特定步骤,如淋巴管内皮粘附或基质重塑,尚未见文献报道,有待进一步研究。
在治疗反应层面,HER2阳性队列研究揭示了TT基因型与靶向治疗耐药的显著关联。接受曲妥珠单抗+帕妥珠单抗联合化疗的TT基因型患者客观缓解率(ORR)仅为33.3%,显著低于非TT型(89.3%,
P=0.001)。新辅助治疗数据显示,TT型患者的病理完全缓解(pCR)率降低约50%(
P=0.018)。这一现象可能涉及多重耐药机制
[1, 30, 31]:首先,HOTAIR通过激活IL-6/STAT3通路诱导HER2蛋白糖基化,降低曲妥珠单抗的结合效率;其次,上调MUC4表达形成物理屏障,阻碍抗体-受体相互作用;此外,其通过ceRNA机制竞争性结合miR-125b,解除对HER3的抑制,激活PI3K/AKT/mTOR信号轴。值得注意的是,双靶向治疗(曲妥珠+帕妥珠)未能克服TT基因型相关耐药,提示HOTAIR可能通过非HER2依赖性通路(如FGFR或IGF-1R)介导治疗逃逸。这些发现为临床治疗策略调整提供了重要依据,例如对TT基因型患者早期采用抗体药物偶联物(ADC)如T-DM1或DS-8201可能更有效。
本研究的临床转化价值体现在多个维度。rs920778分型可作为乳腺癌风险分层工具,尤其适用于有家族史但BRCA1/2阴性的个体。对于确诊患者,TT基因型提示需要更积极的辅助治疗,如CDK4/6抑制剂或ADC药物联合方案。本研究首次将lncRNA遗传变异与靶向治疗疗效相关联,为乳腺癌精准医学开辟了新方向。
本研究仍存在一定局限性。如样本量较小(n=287),可能影响统计效能;缺乏长期随访数据限制了生存预后分析的深度;机制研究局限于相关性分析,需通过类器官模型或转基因小鼠验证因果关系。未来研究应开展多中心队列验证,并深入探索rs920778与治疗耐药的表观遗传调控网络,以及HOTAIR与其他lncRNA(如ANRIL、MALAT1)的协同作用。同时,需要建立标准化检测流程,确保基因分型结果在不同实验室间的可重复性。
综上所述,本研究阐明了HOTAIR rs920778多态性在乳腺癌发生发展及治疗抵抗中的双重作用。TT基因型通过上调HOTAIR表达,驱动肿瘤恶性转化并介导HER2靶向治疗耐药。这些发现不仅深化了对lncRNA遗传调控机制的理解,更为建立个体化诊疗策略提供了理论依据。未来需通过功能基因组学与临床转化研究的深度融合,推动乳腺癌精准医学的跨越式发展,最终实现"基于分子特征的精准干预"目标。