通过体外TFAM基因敲除减少线粒体DNA拷贝数，可以重塑核表观基因组和转录组

《Epigenomics》：Mitochondrial DNA copy number reduction via in vitro TFAM knockout remodels the nuclear epigenome and transcriptome

mtDNA-CN降低通过核基因组甲基化与基因表达调控影响疾病机制。研究利用TFAM基因敲除细胞模型，发现mtDNA-CN减少导致2924个差异甲基化位点、67个差异甲基化区域和102个差异表达基因，涉及GABA受体、神经活性配体信号通路、ABCD1/2基因及细胞信号传导。整合分析表明mtDNA-CN通过代谢物调控（如SAM）和ROS信号影响核表观遗传与转录网络，为年龄相关疾病（如神经退行性疾病、肾脏疾病）的分子机制提供新见解。

　　
线粒体DNA拷贝数（mtDNA-CN）与多种年龄相关疾病及全因死亡率存在显著关联，但其分子机制尚未完全阐明。该研究通过构建TFAM基因敲除的HEK293T细胞模型，系统探究了mtDNA-CN减少对核DNA甲基化及基因表达的影响，揭示了核基因组表观修饰与线粒体功能间的潜在联系。

研究团队采用CRISPR-Cas9技术成功敲除TFAM基因，通过qPCR验证mtDNA-CN在三个独立敲除细胞系中较对照组下降18倍，且该变化在细胞传代过程中保持稳定。为排除脱靶效应，研究者对8个预测的潜在脱靶位点进行测序验证，发现所有目标位点均特异性编辑，未检测到脱靶突变。

在表观基因组分析中，差异甲基化位点达2924个，形成67个差异甲基化区域（DMRs）。值得注意的是，DMRs主要分布在非编码区域，包括增强子、异染色质及开放海区，其中75%区域呈现低甲基化，27%区域甲基化程度升高。这种非均匀的甲基化改变提示mtDNA-CN可能通过代谢中间产物调控表观遗传修饰。

转录组测序发现67个差异表达基因，其中GABA受体亚基基因（如GABRB1）和线粒体相关转运蛋白基因（如ABCD1）显著下调。通过ELMERv2整合分析，确认24个基因-甲基化位点对（Gene-CpG pairs），其中CA2和CA8基因的甲基化改变与表达下调呈现强相关性（Pearson相关系数均＞0.95）。

功能富集分析揭示三大核心通路：首先，GABA能信号通路显著富集，包括GABA受体活性（FDR＜2.3e-3）、神经递质调控（FDR＜3.8e-4）及离子通道复合物（FDR＜2.3e-3）。其次，线粒体-核互作相关通路如ABC转运蛋白活性（FDR＜1.7e-2）和氧化磷酸化调控（FDR＜2.3e-3）被反复提及。第三，非编码调控区域的异常甲基化（如增强子低甲基化、异染色质高甲基化）可能通过改变染色质可及性间接调控基因表达。

机制探讨部分指出，TFAM缺失导致线粒体转录活性下降，mtDNA拷贝数减少可能通过以下途径影响核基因组：1）琥珀酸半醛积累抑制DNA甲基转移酶活性；2）α-酮戊二酸减少影响S-腺苷甲硫氨酸合成；3）线粒体ROS爆发激活Nrf2信号通路，进而改变组蛋白修饰。过表达TFAM实验证实，在48小时干预后，6/8个差异表达基因（如CA2、CA8、RPL22L1）的表达水平可逆性恢复，提示TFAM可能通过调控线粒体代谢产物间接影响核基因表达。

研究创新性体现在首次建立完整的mtDNA-CN→表观修饰→转录调控的三级作用模型。特别发现GABA受体通路作为关键枢纽，其甲基化状态与线粒体功能状态存在双向调控关系：mtDNA-CN减少导致GABRB1基因启动子区低甲基化，激活相关转录因子，促使抑制性突触传递；同时，GABA能神经元功能异常可能通过迷走神经反射进一步影响线粒体功能。这种核-线粒体互作网络可能解释为何mtDNA-CN可作为多系统疾病的生物标志物。

局限性与改进方向：1）细胞模型选择HEK293T系，可能无法完全反映体细胞特异性；2）甲基化与转录调控的因果方向仍需验证；3）样本量较小（n=6），建议后续扩大样本量。未来研究可结合单细胞测序解析时空特异性变化，或利用器官芯片模拟复杂生理环境。

该研究为理解线粒体基因组表观调控提供了重要模型，证实mtDNA-CN通过非编码区甲基化改变影响核基因表达，特别是GABA能信号通路和线粒体-核穿梭机制。这些发现不仅完善了线粒体核基因组互作的理论框架，更为开发基于mtDNA-CN的疾病预警系统及靶向治疗策略（如调节TFAM表达）提供了新思路。