单细胞元分析揭示帕金森病性别差异中转座元件失调

《Neurobiology of Disease》：A single-cell meta-analysis evidences transposable element dysregulation in sex-based differences in Parkinson's disease

　　

转座元件（TEs）（占人类基因组约45%的可移动遗传元件）最近被列为帕金森病（PD）的潜在贡献因素；然而，它们的作用及性别特异性影响仍知之甚少。在此，研究人员展示了首个跨4个黑质单核RNA测序（snRNA-seq）数据集的整合性TE表达元分析，共包含66例患者，构建了PD中TE失调的细胞类型分辨图谱。研究人员在主要脑细胞类型（即神经元、星形胶质细胞、少突胶质细胞和小胶质细胞）中鉴定到广泛的TE激活，其中L1在神经元中显著上调，HERV在少突胶质细胞中显著上调。性别分层分析揭示了明显的男性和女性偏向性TE特征，表明每种性别中受到独特影响的调控程序，包括小胶质细胞中的MIR元件和神经元中的Alu亚家族。相关性和基因组邻近性分析还发现了与重要PD通路（如神经炎症或髓鞘形成）相关的TE-基因关联。总体而言，该研究将TE定位为PD病理学中潜在的性别调节贡献因素，并提供了一个公共网络资源（PATOSS）用于探索PD相关的TE转录失调。

**论文解读：单细胞元分析揭示帕金森病性别差异中转座元件失调的细胞类型特异性特征**

**一、研究背景与意义**

帕金森病（Parkinson's disease, PD）是仅次于阿尔茨海默病的第二大常见神经退行性疾病，也是最常见的运动障碍疾病。全球约有1200万人罹患PD，其发病率与患病率随人口老龄化显著上升。PD的神经病理学特征主要表现为黑质致密部（substantia nigra pars compacta, SNpc）多巴胺能神经元的进行性丢失，导致一系列运动和非运动症状的出现。尽管PD的病因尚未完全阐明，但已有研究确认年龄、农药暴露、特定遗传变异及性别等因素可增加疾病易感性。值得注意的是，大量研究强调了生物学性别对PD病理生理、症状表现及治疗反应的显著影响：男性PD发病率和患病率更高，发病年龄更早，且更易出现僵硬、冻结步态、认知障碍和排尿功能障碍；而女性患者则以震颤、疼痛、抑郁和运动障碍为主要表现。在病理生理层面，女性患者表现出更高的纹状体多巴胺能活性和较低的氧化应激及小胶质细胞激活水平，这可能与雌激素的神经保护和抗炎作用有关。

转座元件（Transposable Elements, TEs）是可移动的重复DNA序列，约占人类基因组的45%，是基因组完整性和细胞功能的重要调节因子。TEs的转录激活对发育和神经元群体遗传嵌合等关键生物学过程至关重要，但其失调可导致严重的有害后果，包括双链断裂累积、基因组完整性受损，以及诱导"病毒模拟"（viral mimicry）状态，激活先天免疫传感器并触发无菌性炎症反应。在衰老相关疾病如PD中，与衰老相关的异染色质破坏和沉默蛋白丢失可扰乱抑制TE激活的表观遗传机制（主要是DNA甲基化和抑制性组蛋白修饰），从而触发TE异常激活。这种失调表达与基因组不稳定性和神经炎症密切相关，可能促进PD的发生和进展。已有研究将LINE-1（L1）和内源性逆转录病毒（ERV）家族的TE多态性和表达增加与PD病理和临床进展联系起来。

然而，目前对TE在神经退行性过程中贡献的确切机制仍不完全清楚。以往研究多采用疾病模型或死后脑组织的批量测序技术，可能掩盖TE表达的细胞特异性转录动态。少数近期研究通过应用snRNA-seq描述了可能与PD标志性过程（如神经炎症）相关的细胞类型特异性TE模式，但不同数据集间的TE表达模式存在动态和不一致性。此外，TE表达失调中的性别差异在很大程度上仍未得到探索，尚无研究在PD中进行TE的性别分层分析。因此，迫切需要采用元分析等整合策略来提取跨数据集的一致信号，并在单细胞水平上全面理解这些性别差异及其相关生物学机制，这可能为针对TE活性的性别特异性治疗干预开辟新途径。

**二、主要关键技术方法**

研究人员依据PRISMA指南进行了系统综述，从GEO、ENCODE、UCSC Cell Browser、Human Cell Atlas和Single Cell Expression Atlas等公共数据库筛选符合标准的人源PD单细胞转录组数据集，最终纳入4个来自SNpc的snRNA-seq研究（Kamath et al., 2021; A. J. Lee et al., 2023; Martirosyan et al., 2024; Smaji? et al., 2021），共计66例患者（含健康对照）。原始数据经fastqc/multiqc质量控制后，使用10× Genomics Cell Ranger（v7.1.0）比对至GRCh38参考基因组，随后采用soloTE软件推断TE转录激活水平。单个数据集处理使用SingleCellExperiment和DropletUtils构建对象，基于UMI数、检测基因数、线粒体基因百分比和双细胞检测进行质量控制与标准化。TE差异表达分析使用MAST R包构建hurdle模型，以疾病状态和性别的组合作为解释变量，并纳入样本来源作为协变量。元分析采用metafor R包的DerSimonian-Laird随机效应模型，计算合并效应量、I2统计量和Cochran's Q检验，并通过Benjamini & Hochberg校正控制多重检验。TE-基因关联网络基于Pearson相关系数构建，并通过供体水平伪bulk相关元分析和双变量随机效应元分析进行验证。功能富集分析使用topGO R包的weight01算法进行GO生物学过程注释。此外，在独立bulk RNA-seq数据集（GSE168496）中使用STAR比对和Telescope定量进行正交验证。所有结果通过R shiny开发的PATOSS网络平台公开发布。

**三、研究结果**

**3.1 跨个体数据集的差异TE表达具有一致性和细胞类型特异性**

研究人员对4个SNpc snRNA-seq数据集进行重新处理和TE定量分析。soloTE分析鉴定到与TE相关的转录本，占大多数样本总UMI的10-20%。TE表达水平在星形胶质细胞中普遍略高于神经元，大多数细胞群体的中位TE表达维持在15-30%。研究在数据集中所有主要细胞类型中均鉴定到显著TE标志物，例如L1PA12被鉴定为显著的神经元标志物，在约90%的神经元中表达；LTR46（ERV/LTR家族成员）被鉴定为少突胶质细胞标志物。跨数据集比较分析揭示了每种细胞类型的一小组共同标志物（少突胶质细胞8个、星形胶质细胞9个、神经元6个、小胶质细胞4个），可作为主要脑细胞类型的特征性标记。差异表达分析在PDvsC（PD vs Control）、IDF（女性疾病影响）、IDM（男性疾病影响）和SDID（性别差异疾病影响）四种对比中鉴定了差异表达TE（DETE）。值得注意的是，神经元在PDvsC对比中产生最多DETE而小胶质细胞最少，但在SDID对比中这一模式发生反转，小胶质细胞比神经元获得更多DETE。

**3.2 共识TE模式呈现L1和HERV细胞类型特异性失调**

通过细胞类型特异性元分析，共鉴定到537个DETE（IDF 182个、IDM 81个、SDID 30个、PDvsC 244个），较之前同一数据集的基因表达元分析鉴定到的14457个DEG显著减少，反映了TE亚家族数量较少（约1000）及其检测稳定性较低的特点。跨对比和细胞类型的DETE比较分析显示重叠有限，提示TE失调具有高度对比和细胞类型特异性。在L1家族中，PDvsC对比显示所有细胞类型均存在失调亚家族，其中神经元表现出最强信号，39个L1亚家族在PD中显著上调（如L1MC5）。性别特异性信号方面，少突胶质细胞中L1M3c和L1MDa在男性中表现出更强的PD相关效应。HERV家族在少突胶质细胞中显示最高数量的失调亚家族（31个DETE），且在SDID对比中是最常失调的TE，大多数元件在男性中表现出更强的PD相关效应（如LTR12C和Harlequin-int）。独立bulk RNA-seq队列（GSE168496）的正交分析显示候选TE和集合水平分析与元分析结果具有一致的方向趋势。

**3.3 相关的TE-基因关联显示更短的基因组距离并涉及PD相关生物学通路**

研究人员计算DETE与DEG在4种对比中的logFC相关系数，以推断潜在功能关联。距离分析显示，除星形胶质细胞外，所有细胞类型中相关DETE与DEG之间的距离显著短于非相关DETE与DEG之间的距离，其中小胶质细胞和神经元的显著性和效应量最高。基于这些关联构建的TE-基因网络揭示了不同细胞类型的特异性模式：星形胶质细胞网络含146个关联，连接125个DEG与8个DETE亚家族，其中DNA转座子HSMAR2作为主要枢纽（56个DEG），富集的生物学过程包括端粒组织与维持、先天免疫反应和轴突延伸；神经元网络最小（80个关联），含77个DEG与8个DETE亚家族，大多数DETE表现出女性偏向性PD相关效应，富集过程涉及淀粉样β肽加工、泛素化和凋亡过程以及NF-κB信号通路；小胶质细胞网络含126个关联，MIR家族（MIRc和MIR3）和L1MDa LINE为主要贡献者，且均表现出男性偏向性PD相关效应，富集分析显示与前炎症细胞因子产生通路强相关；少突胶质细胞网络最大（658个关联），涉及287个DEG和49个DETE亚家族，富集过程包括DNA代谢加工、基因调控与转录、病毒基因组复制、髓鞘形成和氧化应激反应。供体水平伪bulk相关元分析和双变量随机效应模型验证了优先枢纽关联的稳健性。

**3.4 失调TE更接近PARK基因，主要在神经元中**

研究人员比较了DETE与非显著TE与60个PD最高排名基因（来自OpenTargets数据库）的最短基因组距离。总体而言，除少突胶质细胞外，所有细胞类型中DETE出现在显著更接近大部分PD基因的位置。神经元在全部基因面板上均显示显著结果和更大效应量。核心PD基因（PARK基因）中，神经元中的SNCA（PARK1/4）和DJ-1（PARK7）、星形胶质细胞中的PINK1（PARK6）以及小胶质细胞中的GBA1显示出DETE更接近其基因组位置的明显趋势。进一步分析选择各细胞类型中在PDvsC对比中显著差异表达的PD标志基因，鉴定到小胶质细胞中的LRRK2（PARK8）、VPS35（PARK17）和ATP6V1B2，星形胶质细胞中的COMT和MAPT，以及神经元中的BCHE、TMEM175、DNAJC13和VPS35。

**3.5 PATOSS：探索PD中共识TE转录失调的网络资源**

研究人员开发了PATOSS（PArkinson Transcriptomic Overview of Sex-differences at a Single cell resolution）网络资源平台（https://bioinfo.cipf.es/cbl-atlas-pd/），提供三个主要功能：TE表达探索（动态表格和自定义条形图/热图）、TE-基因网络探索（交互式图形和边缘支持数据表格）、TE富集探索（点图和基因/TE信息表格）。

**四、讨论与结论**

本研究提供的首个PD中TE失调的单细胞元分析图谱，揭示了SNpc中稳健的细胞类型特异性和性别依赖性TE特征。研究发现神经元经历最广泛的TE激活，特别是L1家族中39个亚家族在PD中显著上调，这可能与多巴胺能环境中高水平活性氧（ROS）触发的DNA去甲基化有关。DETE与PD标志基因位点（如神经元中的SNCA和PARK7、星形胶质细胞中的PINK1、小胶质细胞中的GBA1）的邻近性，提示TE可能作为假定转录顺式调控元件影响基因网络，或反映共享表观遗传脆弱性的基因组区域。性别特异性分析揭示了各细胞类型中独特的性别偏向性TE特征：神经元中的Alu亚家族（如AluYc）与女性偏向性PD效应相关，小胶质细胞中的L1MDa和MIR3/MIRc与男性偏向性效应相关，少突胶质细胞中的HERV TE也主要表现出男性偏向性效应。这些发现表明TE不仅是神经退行性变的次要后果，更可能是性别差异性生物学机制的活跃中介者。研究也存在局限性，包括TE异质性、TE-基因共表达网络仅基于4种对比计算的统计功效限制，以及纳入研究中性别代表不均（如GSE157783仅含每条件1例女性样本）可能引入的偏倚。未来需要结合多组学数据和实验模型来充分验证TE激活与PD发病之间的因果关系。结论部分指出，该研究提供了首个PD中TE失调的单细胞元分析，揭示了SNpc中稳健的细胞类型特异性和性别依赖性TE特征，表明TE可能参与PD相关的关键过程，影响疾病相关基因调控和性别特异性病理，将TE确立为PD中有前景的分子参与者，为未来机制和治疗研究提供了宝贵途径。