唐氏综合征胎盘差异miRNA表达谱分析：基于全转录组测序分析技术

唐氏综合征（DS）是最常见的染色体非整倍体，发病率约为1∶700活产。该综合征包括80多种影响所有主要器官的临床表现，以及不同程度的认知功能障碍和颅面异常，这是唐氏综合征患者最常见的临床特征。同时，唐氏综合征常伴发生先天性心脏缺陷、胃肠道异常、儿童白血病或早发性阿尔茨海默病等。唐氏综合征的表型可能是基因失衡所引起，由于21号染色体基因剂量增加了50%且存在与其他基因相互不平衡的作用。目前，唐氏综合征在临床缺乏治疗的有效手段，主要通过产前筛查和诊断的方式减少唐氏患儿的出生。

目前常用于唐氏儿的筛查方法是一种联合检测，包括从母体血液中提取的几种生化标记物和胎儿超声标记物，但检出率仅60%～70%。无创胎儿DNA检查对于唐氏综合征检出率较高，但检查费用高，技术操作复杂，检查周期长，不易在基层医疗机构推广。筛查的阳性结果需要通过侵入性产前诊断检测绒毛膜绒毛取样（CVS）或羊膜穿刺术来确认。然而，这两种手术对操作者的临床经验要求高，且仍有1%导致流产的机率。因此，建立一种安全、价格合理又特异性的唐氏综合征产前筛查手段仍是产前诊断领域探索的方向。

miRNA是真核生物中常见的一类由内源基因编码的长度约为20～25个核苷酸的非编码单链RNA分子，参与转录后基因表达的调控。继妊娠者血浆中胎儿游离DNA被发现后，有文献报道母体外周血中还存在胎儿游离RNA。胎盘是胎儿特有的组织，筛出胎盘特异表达miRNA会产生新的胎儿特异性标志物，对其检测也更具针对性，可以提高检测的特异性和灵敏性，但目前最具特异性的miRNA分子尚有待探讨。

本研究重点比较了正常胎盘和唐氏综合征胎盘中miRNA的表达；通过使用阵列技术对82个miRNA全转录组进行了分析；使用生物信息学工具分析选定的miRNA的靶基因，以绘制潜在的受影响的功能通路。

1 资料和方法

1.1 一般资料

选取同期昆明市妇幼保健院产科产前诊断确诊唐氏综合征孕妇（=3，年龄分别为28，35，47岁）及正常对照孕妇胎盘（=3，年龄分别为29，30，31岁）。本研究通过医院伦理委员会审核（审批号：2021-5-7），孕妇均被提供遗传咨询并签署知情同意书。

1.2 唐氏综合征产前诊断的细胞遗传学和高通量测序分析

使用标准程序对胎儿羊水样本进行染色体分析。超声引导下羊膜腔穿刺术，抽羊水20 mL，离心后弃上清，加入培养基，原位法进行羊水细胞培养，37 ℃培养7～9 d。用GTG显带法对中期染色体进行染色，每个病例检查30个分裂像，分析5个分裂像。采用德国QIAGEN公司DNeasy Blood and Tissue Kit（69506）试剂盒提取羊水细胞DNA。采用染色体CNV 检测试剂盒（KR0040）进行DNA 文库构建。在Illumina Next-Seq500测序平台对DNA文库进行测序。

1.3 标本处理

采用SPSS 23.0进行统计学分析，多组间比较使用Chi-square（N×N）检验分析，两组间比较使用-test分析。对差异表达基因的筛查以|log2(Foldchange)|＞1且＜0.05为阈值。＜0.05为差异有统计学意义。

1.4 转录组测序及数据分析

采用高通量测序技术。选取唐氏综合征胎盘和正常胎盘组织各3例，对子面、母面、绒毛层组织进行混样，充分研磨后TRIzol法提取总RNA。miRNA测序文库制备采用TruSeq Small RNA Sample Prep Kits（Illumina）试剂盒。文库制备工作完成后，对构建好的文库使用Illumina Hiseq2000/2500进行测序。实验流程按照说明书的标准步骤执行。对比分析唐氏综合征胎盘与正常胎盘组织中的差异表达miRNA，以差异倍数≥1.5且＜0.05及检测丰度较高作为显著上调和显著下调的标准。通过miRBase（https://www.mirbase.org/）查询目标miRNA的序列信息及所在染色体。对差异表达miRNA的潜在靶基因进行GO和KEGG-Pathway功能富集分析。

通过使用miRWalk、Targetscan、miRDB软件工具对上调miRNA的靶基因进行预测并通过Venny软件分析，发现上调miRNA的共同靶基因为BTBD3、AUTS2、ALS2等。目前关于BTBD3的功能在很大程度上是未知的，Fagerberg等报道BTBD3基因特异性在胎盘组织和脑组织中存在高表达（图3E）。另有研究发现，在新生儿发育过程中，BTBD3可以将树突导向活动的轴突终末，并调节初级感觉皮层中活动依赖的树突修剪，从而促进神经回路的形成。在人类和啮齿类动物中，BTBD3在皮层-纹状体-丘脑-皮层回路中高表达，包括中背丘脑、前扣带皮层和海马体。在这个回路中，BTBD3的表达在出生后早期迅速增加，在人类和小鼠中，在儿童和青少年时期达到高峰。多篇报道显示BTBD3与神经细胞发育有关。另一靶基因AUST2（自闭症易感性候选基因2）已成为一个与广泛的神经心理疾病相关的关键基因，如自闭症谱系障碍、智力障碍、精神分裂症和癫痫。AUTS2位于染色体7q11.22上，该区域与染色体断点的易感性相关。它共跨越1.2 Mb，包括19个外显子，其中外显子1-6位于5'末端，被大内含子分开，外显子7-13紧凑排列在基因的3'末端。最初在2002年被报道为一对同卵双胞胎因染色体新发平衡易位使AUTS2 基因功能被中断而导致自闭症。AUTS2已被证明参与多种神经发育过程：在细胞核中，它作为神经发育的关键转录调控因子，而在细胞质中，它通过控制细胞骨架重排参与大脑皮质发生，包括神经元迁移和神经发生。在产后，AUTS2调节兴奋性突触的数量，以维持神经回路中兴奋和抑制之间的平衡。在本次唐氏综合征胎盘组织转录组基因测序中显示AUST2、BTBD3 和ALS2 转录组中表达下调。在AUST2与BTBD3均参与神经细胞的发育，鉴于唐氏综合征患儿均存在智力障碍的问题，其是否与BTBD3、AUST2基因表达下调相关是我们下一步将继续研究的方向。

1.5 统计学分析

胎盘娩出后2 h内完成组织采样。以脐带附着点为中心，在对称的四个方位（ABCD），子面、母面、绒毛层分别取1 cm样本，放入含RNA-later 试剂的EP 管中4 ℃过夜，后置于液氮中速冻，之后置于-80 ℃冰箱保存。

2 结果

2.1 唐氏综合征产前诊断结果分析

已有研究显示，通过检测母体外周血中胎儿游离RNA，可将唐氏综合征诊断的敏感度和特异性分别提高至90.0%和96.5%。胎儿来源的游离RNA在母体外周血中以亚细胞小体形式存在，不易被RNA酶降解，且在妊娠期母体血液中含量丰富，有望成为无创途径获得胎儿遗传物质的另一个候选标记物。近年来孕妇血液特异的miRNA检测有望成为一种新型的无创产前检测手段逐渐受到重视。Karaca等对怀有唐氏综合征胎儿的孕妇的羊水miRNA进行研究发现21号染色体相关性miRNA-125b-2，miRNA-155 和miRNA-3156 在怀有唐氏综合征胎儿的孕妇羊水中明显高于正常妊娠或未妊娠女性。我们的检测结果（表1）中显示miR-125b-2-3p在唐氏综合征胎盘中的表达量高于正常胎盘，与报道一致。

她笑后又接着说：“你们男人啊，是不是一看见好看的女人，就像老猫子嗅到腥一样，就是不要了老命，也恨不得一口把她含在嘴里才过得？”

2.2 全转录组测序及质量分析

通过对所有差异表达miRNA 的靶基因GO 和Pathway分析发现，GO富集分析结果表明，生物过程模块多数为参与RNA代谢过程及转录调控，其余与感觉器官发育、眼睛发育等相关，RNA转录过程的调控推测与唐氏综合征患儿各类代谢产物异常表达有关。参与的细胞组分方面，多数富集在细胞膜及细胞质相关组分，其余在整个细胞器、细胞核均有涉及。分子功能近年成为在肿瘤、代谢病及遗传病的研究重点，在本文中，唐氏综合征的分子功能富集主要显示在蛋白质结合和金属离子结合方面，且在各类酶活性中均有富集。研究报道在唐氏综合征中神经生长因子的代谢途径发生了改变：通过水解酶、金属蛋白酶等降低神经生长因子的数量。其余富集到转录调控过程等。

2.3 唐氏综合征胎盘中miRNA差异表达情况

通过对所有差异表达miRNA 的靶基因GO 和Pathway分析发现，GO富集分析结果表明，差异表达miRNA的靶基因主要参与的分子功能包括蛋白结合、金属离子结合、转移酶活性、核苷酸结合、DNA结合、水解酶活性、ATP结合、RNA结合、DNA依赖性RNA聚合酶Ⅱ启动子转录调控、同种蛋白质结合；参与的细胞组分包括质膜部分、胞质部分、细胞核部分、膜的整体成分、胞质部分、核腔、核质部分、突触、细胞连接、细胞器、细胞器腔、膜封闭腔、质膜部分、核仁、突触部分、细胞内非膜细胞器、非膜细胞器、核外围、核基质；参与的生物过程包括转录调控、RNA 代谢过程调控、DNA 依赖性RNA聚合酶Ⅱ启动子转录调控、眼睛发育、感觉器官发育、转录负调控、依赖DNA、磷酸化、磷代谢过程、磷酸盐代谢过程、mRNA代谢过程。KEGG富集分析结果表明差异表达miRNA靶基因主要参与的信号通路包括肿瘤相关信号通路、PI3K-Akt信号通路、Ras信号通路、Rap1信号通路、细胞骨架调控信号通路、嘌呤代谢相关信号通路、P53相关信号通路（图4）。

2.4 差异表达miRNA的靶基因功能富集分析结果

实验筛选出82个差异表达的miRNA。与正常胎盘组织相比，唐氏综合征胎盘组织中有29个miRNA表达上调（变化倍数≥2，＜0.05），15个miRNA表达下调（变化倍数≥2，＜0.05）；根据表达丰度（表达丰度300以上）共筛选出4个表达上调的miRNA，6个表达下调的miRNA。差异表达的miRNA分布在不同的染色体上，其中包括一个21号染色体衍生的miRNA（表1）。差异表达miRNA的火山图、整体层次聚类图显示这些miRNA在两组中表达有显著差异（图3A～C）。通过miRWalk、Targetscan、miRDB软件预测靶基因后进行分析，miRNA-323a-3p、miRNA-199b-5p 和miRNA-125b-2-3p 共同靶基因为ALS2（鸟嘌呤核苷酸交换因子）。miRNA-323a-3p 和miRNA-199b-5p 共同靶基因为BTBD3、AUTS2等（图3D）。对95个人的组织样本进行RNA-seq测序，靶基因BTBD3在27个不同组织中的表达情况（图3E，数据来自Bioproject：PRJEB4337。https://www.ncbi.nlm.nih.gov/gene/22903/）。

3 讨论

唐氏综合征的发生与母亲的怀孕年龄有关，高龄孕妇唐氏综合征患儿的发病率更高，卵子老化也与唐氏综合征的发生有着不可分离的关系。据报道，30岁以下的孕妇孕育唐氏综合征胎儿风险为1/1000，35岁以上孕妇风险率1/270，40岁以上孕妇风险率为1/100，45岁以上孕妇风险率为1/20。95%左右的唐氏综合征患者父母均为正常且家族中也没有唐氏综合征，其发生是受精卵早期细胞分裂错误或生殖细胞（精子或卵子）分裂错误导致的。仅有5%以下的唐氏综合征的发生与父母染色体结构异常（如易位）有关。所以理论上讲，不管是否有家族史，所有怀孕的人都应该进行唐氏筛查，提高孕妇及家属对于唐氏综合征系列产前筛查的重视程度，使更多孕妇主动接受唐氏综合征系列产前筛查，进行唐氏综合征患儿的早期诊断，具有非常重要的现实意义。

唐氏综合征病例产前诊断羊水细胞染色体核型和高通量测序结果一致，均为21三体（图1）。

截止至2018年5月31日，“体育舞蹈圈”公众号粉丝已达到238 805人，并覆盖全国包括台湾在内的34个省、市、自治区、直辖市，目前粉丝数量仍处于健康增长状态。此外，根据线上问卷调查结果显示，如表1，“体育舞蹈圈”微信公众号粉丝群体大致可分为体育舞蹈初学者、等级考试参与者、体育舞蹈在职教师三类群体，且三者需求均以学习体育舞蹈教材及配套等级组合视频为主。

测序结果经base calling软件CASAVA分析后Q值（碱基质量值）显示由于测序一起造成的某个位置碱基错误率小于1%，测序质量可靠。在样本间重复性分析中使用Pearson相关性分析、PCA分析（主成分分析）来判断，显示数据可信度良好（图2）。

平顺县隶属于山西省东南部的长治市，位于太行山西侧，黄土高原东部边缘，处于我国由华北平原向第二阶梯过度的前沿，平均海拔高于1000m的地区占全县总面积的70%以上（图1）。全县地势南高北低、自东南向西北倾斜。境内河流则有浊漳河，东部的虹霓河、南部宑底河系集东南部山地雨水形成，流域面积均较小[10]。县境内地貌以山地为主，平原仅在浊漳河两岸有少量分布[11]。

在本次分析结果中，KEGG富集分析结果表明差异表达miRNA靶基因主要参与的信号通路有肿瘤相关信号通路、PI3K-Akt信号通路等。PI3K-Akt信号通路作为脑组织受体，推测其在唐氏综合征患儿脑发育过程起到重要作用。近年来，关于唐氏综合征生物标记物及发病的分子机制已有一定研究，但尚缺乏更多临床试验验证，本次的检测可能为唐氏综合征发病的分子机制提供了新论据。

工业革命之前，技术进步非常缓慢，一个人一生几乎感受不到技术的进步，沿着同样的技术路径劳动和生活，可谓天经地义。工业革命以后，尤其是20世纪下半叶以来，一个人可以感受到技术的明显进步。1980年代以来，每十年甚至每年都可以感受到技术的明显进步。

本研究所纳入的临床试验普遍质量较低，缺乏高质量、标准化的RCT。由于各研究间针刺所选穴位不同、服用的西药和中成药不同，且各研究间所纳入患者的病情严重程度不一等因素的影响，造成分析过程中各研究间存在较大异质性。在治疗前后NIH-CPSI评分变化量与白细胞降低有效率作为结局指标的分析中，由于纳入的相关分析较少，得出的结果可能可靠性不强。故在今后要多进行相关的临床研究，尤其是与中药汤剂、中药坐浴和微波治疗等方面的对比研究，且要进行高质量、标准化的RCT；尽量多注重治疗后随访的结果，为针刺治疗CP远期疗效提供可靠证据。

综上所述，可能提示miRNA的过度表达属于影响胎盘形成模式及对胎儿发育的所有后果的重要因素。唐氏胎儿在早期发育时胎盘组织中miRNA表达出现上调，进而抑制胎儿神经发育相关靶基因的表达。我们下一步进行血清学验证miRNA表达情况，可用于弥补血清学筛查在唐氏综合征筛查的特异性不足的问题，提高孕早期唐氏综合征筛查的准确性。