衰老小鼠前额叶皮层LINE-1表达上调与神经元基因网络下调的关联研究

《Frontiers in Aging Neuroscience》：Increased LINE-1 expression is associated with neuronal gene network downregulation in the aged mouse prefrontal cortex

　　

背景： 衰老与表观遗传改变、神经炎症及神经元功能的进行性衰退相关。逆转座子LINE-1（L1）的去抑制已被认为与衰老及年龄相关疾病有关。然而，在正常衰老大脑中，L1的区域性表达及其与DNA甲基化、基因组L1拷贝数以及衰老相关转录变化之间的关系仍未完全阐明。 方法： 研究人员检测了年轻（17周龄，n=10）和衰老（101周龄，n=7）雄性C57BL/6J小鼠前额叶皮层（PFC）、海马（Hp）和伏隔核（NAc）中的L1表达。L1表达、DNA甲基化和基因组拷贝数分别通过RT-qPCR、焦磷酸测序和基因组qPCR进行定量。通过RNA-seq分析PFC的转录组变化，并随后进行加权基因共表达网络分析（WGCNA）。 结果： 在衰老PFC中，多个活跃L1亚家族的L1表达显著增加，而在其他脑区未观察到显著变化。尽管L1表达升高，但活跃L1启动子的DNA甲基化或基因组L1拷贝数在年龄组之间均未检测到明显变化。WGCNA鉴定了三个主要的衰老相关模块。其中，仅神经元模块与L1表达呈显著负相关，且该模块富集了参与突触信号传导、神经元发育和神经元投射的基因。 结论： L1表达在衰老小鼠PFC中增加，并与神经元基因网络的协调性下调相关，而非与其他主要的衰老相关转录程序相关。这些发现进一步表征了正常脑衰老过程中L1表达增加的分子背景。

论文解读

一、研究背景与意义

衰老是一个伴随多器官功能衰退和疾病易感性增加的渐进性生物学过程。在分子层面，衰老以基因组不稳定、表观遗传失调和慢性炎症为特征，同时伴随染色质组织的改变，包括异染色质松弛以及富含转座元件（TEs）的重复基因组区域内转录抑制的丧失。长散在核元件-1（LINE-1，L1）是哺乳动物基因组中唯一具有自主活性的逆转座子家族，约占小鼠和人类基因组的20%。完整的L1元件包含5′非翻译区（UTR）内部启动子及编码ORF1和ORF2蛋白的两个开放阅读框，这些蛋白通过RNA中间体介导逆转座，导致新的L1拷贝插入基因组。尽管L1活性在体细胞组织中通常被表观遗传和细胞防御机制所抑制，但越来越多的证据表明，在衰老和细胞衰老过程中L1的抑制会被削弱。

已有研究提出了L1可能通过多种机制促进衰老：衰老相关的L1抑制减弱与调控因子改变及DNA甲基化降低有关；在衰老细胞中，L1去抑制导致L1来源的核酸在细胞质中积累，激活先天免疫通路并促进以I型干扰素信号为特征的衰老相关分泌表型（SASP）；此外，L1 RNA还涉及早衰综合征中的异染色质侵蚀和细胞功能障碍。然而，在正常脑衰老过程中，L1激活是否伴随类似机制仍不清楚。近年来的证据表明，脑衰老过程中存在转座元件调控的年龄相关变化，包括ORF1表达增加和细胞类型依赖性的表观遗传轨迹改变，这提示L1失调可能与衰老大脑中的特定分子程序相关。然而，L1表达增加是否与特定的衰老相关转录程序相关联，抑或仅反映更广泛的年龄相关分子变化，仍不明确。基于L1去抑制与衰老及年龄相关神经元改变之间的关联，研究人员假设L1表达在正常脑衰老过程中会以区域特异性方式增加，并与神经元功能相关的协调转录变化相关联，同时进一步检验L1表达变化是否伴随DNA甲基化或基因组拷贝数的改变。

二、研究设计与技术方法

本研究使用雄性C57BL/6J小鼠，分为年轻组（17周龄，n=10）和衰老组（101周龄，n=7），来源于日本Charles River Laboratories。研究人员主要采用了以下关键技术方法：

1. 1.
   RT-qPCR：使用针对L1保守序列（5′UTR、ORF1、ORF2）及L1亚家族特异性引物，定量检测PFC、Hp和NAc三个脑区的L1表达水平及衰老和L1调控相关基因的表达。
2. 2.
   焦磷酸测序（Pyrosequencing）：检测活跃L1亚家族（TfI、GfII和A）5′UTR区域选定CpG位点的DNA甲基化水平。
3. 3.
   基因组qPCR：以5S rRNA为内参，定量检测L1基因组拷贝数。
4. 4.
   RNA-seq：对PFC样本进行链特异性mRNA文库构建和Illumina NovaSeq 6000双端150bp测序，每样本约2000万读段对，比对至GRCm38参考基因组。
5. 5.
   WGCNA：基于TPM矩阵构建基因共表达模块，识别衰老相关模块并分析模块特征基因与L1表达的相关性。
6. 6.
   TEtranscripts和MORE-RNAseq：分别用于定量总L1表达和具有逆转座能力的全长完整L1（rc-L1s）表达。

三、研究结果

3.1 衰老小鼠PFC中L1表达增加

通过RT-qPCR检测L1 mRNA表达，发现衰老组小鼠PFC中L1表达在5′UTR、ORF1和ORF2三个区域均显著高于年轻组，而Hp和NAc中年龄相关增加极为有限。进一步分析L1亚家族表达，PFC中所有检测的亚家族（A、Gf和Tf）在衰老小鼠中均显著增加，而Hp和NAc中未检测到显著变化。

3.2 衰老小鼠脑中衰老和L1调控相关基因的表达

Bdnf在PFC和NAc中显著降低，与其在脑衰老中的已知作用一致。Trex1和Bmi1在三个脑区中均显著升高，Rb1和Yy1在NAc中显著升高。Sirtuin家族方面，PFC和Hp中Sirt3显著增加，NAc中Sirt1、Sirt4和Sirt5显著增加。总体而言，衰老伴随多个L1调控基因和sirtuin家族成员的表达增加，且表达模式因脑区而异。

3.3 PFC中L1的DNA甲基化水平未改变

通过亚硫酸氢盐PCR结合焦磷酸测序检测活跃L1亚家族5′UTR区域选定CpG位点的DNA甲基化水平，衰老小鼠在所检测的CpG位点未观察到显著的甲基化变化。

3.4 前额叶皮层中L1 DNA拷贝数未改变

通过qPCR检测L1基因组拷贝数，年轻和衰老小鼠之间无论是保守L1序列还是活跃L1亚家族均无显著差异。

3.5 RNA-seq揭示与衰老过程中L1表达增加相关的转录组改变

RNA-seq分析确认了RT-qPCR的结果，包括衰老小鼠中Bdnf表达降低以及Bmi1和Sirt3表达增加。TEtranscripts分析确认了衰老小鼠L1表达增加。MORE-RNAseq定量结果显示rc-L1s表达与总L1表达呈正相关（r=0.56，P=0.02）。共鉴定出570个差异表达基因（DEGs），其中422个上调、148个下调。

3.6 衰老相关基因共表达模块的鉴定

WGCNA分析鉴定出20个共表达模块。模块-性状相关分析显示三个模块与衰老显著相关：棕色模块与衰老呈强负相关（r=?0.903），蓝绿色模块（r=0.724）和浅绿色模块（r=0.637）与衰老呈正相关。DEG富集分析证实仅这三个模块显著富集衰老相关DEGs，共包含89.1%的衰老相关DEGs。功能富集分析显示：棕色模块富集神经元和突触功能相关基因（突触组织、树突发育、神经递质信号）；蓝绿色模块富集核糖体和应激反应相关基因（核糖体生物发生、蛋白质量控制、氧化应激、抗病毒反应）；浅绿色模块富集髓鞘形成和胶质细胞激活相关基因（少突胶质细胞分化、轴突包绕、白细胞迁移）。

3.7 L1表达优先与神经元共表达模块相关

以ORF2表达作为L1活性的指标（与TEtranscripts估计的L1表达强相关，r=0.81），分析模块特征基因与L1表达的相关性。在三个衰老相关模块中，仅棕色模块（神经元模块）与ORF2表达呈显著负相关（r=?0.662，FDR校正后P=0.0114），蓝绿色和浅绿色模块均无显著关联。

四、讨论与结论

本研究表明L1表达在衰老小鼠PFC中选择性增加，而Hp和NAc中仅观察到微小变化。通过RNA-seq与WGCNA整合分析，研究人员识别出脑衰老的三个主要转录程序：神经元基因表达的协调下调、细胞应激相关通路的激活以及胶质相关基因表达的增加。值得注意的是，在衰老相关模块中，仅神经元共表达模块与L1表达显著相关。L1表达与神经元模块之间的负相关可能反映了共享的年龄相关调控环境，其中神经元转录程序减少而L1抑制被放松。然而，现有数据无法确定神经元转录变化是导致L1去抑制的原因还是L1活性增加的结果。

L1表达增加并未伴随活跃L1启动子的低甲基化或基因组L1拷贝数的增加，这可能与DNA甲基化检测使用bulk PFC组织有关，细胞类型特异性的甲基化变化可能被稀释或掩盖。此外，检测的CpG位点仅覆盖三个活跃L1亚家族5′UTR的选定区域，不能排除更局部或细胞类型特异性的L1甲基化改变。Trex1、Rb1、Yy1和Bmi1等L1调控基因在衰老PFC中并未出现表达降低，而是维持或上调，提示正常脑衰老与病理性衰老或衰老细胞中L1调控的转录调控存在差异。

研究局限性包括：样本量相对较小（衰老组n=7）限制了相关分析的统计功效；研究主要基于转录组关联，无法建立L1去抑制与神经元功能障碍之间的因果关系；bulk组织分析无法确定L1表达增加的具体细胞类型；无法排除罕见体细胞逆转座事件；仅纳入雄性小鼠，性别差异的影响尚不清楚。

研究结论： L1表达在正常脑衰老过程中的前额叶皮层中增加，并且在脑衰老过程中发生的主要转录变化中，L1表达优先与神经元转录程序相关联。