生殖衰老小鼠卵巢中可变剪接模式的失调

Dysregulation of alternative splicing patterns in the ovaries of reproductively aged mice

作者信息Adnan T Alsamaraee, Vanessa L Correll, Julius O Nyalwidhe, Pavla Brachova, Nehemiah S Alvarez
PMID41630016
发布时间2026-02
DOI10.1093/reprod/xaag019

实验完整度

包含直接RNA测序、Bayesian统计模型、scRNA-seq验证、质谱蛋白组学验证和AlphaFold结构建模等多层级证据,并进行了功能验证。

主要模型

年轻(6周龄)C57BL/6小鼠卵巢 生殖衰老(14月龄)C57BL/6小鼠卵巢 PMSG刺激模型 PMSG+hCG刺激模型

重点核对

小鼠年龄:6周龄(年轻)与14月龄(生殖衰老)C57BL/6 促性腺激素刺激方案:PMSG(5 IU)或PMSG(5 IU)+hCG(5 IU),hCG在PMSG后46小时注射 卵巢采集时间:PMSG注射后46小时或hCG注射后16小时 生物重复数:每组3只小鼠(直接RNA测序),scRNA-seq每组4只 统计阈值:95%最高后验密度区间不含零且后验概率>0.95

摘要

雌性生殖衰老的特征是卵巢功能的进行性衰退,然而驱动这些变化的分子机制仍未完全阐明。在这里,我们使用长读长直接RNA测序来绘制小鼠卵巢在生殖年龄期间的转录本异构体变化。在控制促性腺激素刺激后比较年轻和年老小鼠,我们鉴定了广泛的剪接变化,包括外显子使用、剪接位点选择和转录本边界的改变。老年卵巢表现出增加的异构体多样性,倾向于远端起始和终止位点,并且外显子跳跃和内含子保留事件显著增加。许多这些年龄偏倚的剪接事件改变了开放阅读框,引入了提前终止密码子,或破坏了保守的蛋白结构域。值得注意的是,几个参与呼吸链的线粒体基因受到影响。我们重点介绍了Ndufs4,一个线粒体复合物I亚基,作为衰老促进缺乏典型蛋白家族(Pfam)结构域的短异构体剪接的案例。结构建模表明这种剪接变体可能损害复合物I功能,导致活性氧产生增加。我们的数据提示剪接与衰老卵巢中的线粒体功能障碍之间存在机制联系。这些发现支持卵巢生理学中“剪接-能量-衰老轴”模型,其中线粒体功能下降和适应性或适应不良的剪接变化相互交织。我们的研究揭示可变剪接不仅仅是衰老的副产品,而是一个动态的、转录组范围的调控层,可能影响卵巢寿命。这些见解为研究生殖衰老中的转录后机制开辟了新途径,并强调在卵巢衰退模型中需要考虑异构体水平的调控。

实验结论

提炼研究问题、关键发现与证据,快速把握文章的核心贡献。

研究问题
生殖衰老过程中,卵巢转录组的剪接模式如何发生改变,这些改变是否与线粒体功能障碍相关?
核心机制
衰老促进卵巢中广泛的剪接异常,包括外显子跳跃、内含子保留增加,以及转录本边界向远端偏移,其中线粒体复合物I亚基Ndufs4的剪接变体(缺失Pfam结构域)可能导致复合物I功能受损和ROS增加,构成“剪接-能量-衰老轴”。
主要证据
通过直接RNA测序比较年轻和衰老卵巢,发现大量年龄偏倚的剪接事件;通过Pfam结构域分析和AlphaFold建模,揭示Ndufs4短异构体可能损害复合物I功能;通过scRNA-seq和质谱验证了异构体的细胞类型特异性和翻译。
研究意义
本研究揭示可变剪接在卵巢衰老中的系统性重编程,支持“剪接-能量-衰老轴”模型,为理解卵巢衰老的转录后调控机制提供新视角,并可能为改善年龄相关生育力下降提供新靶点。

研究路径

按研究推进顺序梳理实验设计、验证步骤与关键观察。

1

卵巢组织收集与RNA提取

获得年轻和衰老小鼠卵巢的RNA用于直接RNA测序。

对6周龄和14月龄C57BL/6小鼠进行PMSG或PMSG+hCG刺激后,收集卵巢,使用TRIzol提取总RNA。

2

直接RNA测序文库构建与测序

获得全长转录本序列,用于分析异构体。

使用Oxford Nanopore Direct RNA Sequencing Kit (SQK-RNA002)构建文库,在MinION Mk1B上测序。

3

剪接事件注释与定量

识别和定量不同年龄卵巢间的可变剪接事件。

使用FLAIR比对和整理转录本,然后用自定义Python模块注释剪接事件,并计算percent spliced (PS)值。

4

剪接差异的统计分析

确定在年轻和衰老卵巢之间具有显著差异的剪接事件。

采用贝叶斯beta-binomial模型比较异构体使用比例,计算log2 fold change和95%最高后验密度区间,以及后验概率。

5

剪接事件对蛋白结构域和功能的影响分析

评估年龄偏倚的剪接事件对蛋白结构域的影响。

将在silico翻译的异构体与Pfam数据库比对,分析蛋白结构域的变化,并对Ndufs4进行AlphaFold结构建模。

6

正交验证

验证差异剪接异构体的存在和翻译。

使用scRNA-seq验证异构体的细胞类型特异性表达,并使用质谱检测蛋白肽段。

研究方法

按研究目的归类文中使用的方法,便于定位所需技术。

组织处理和形态学分析
转录组测序
剪接事件分析
数据分析与统计

产品清单

实验环节名称品牌货号
TRIzolInvitrogen--
直接RNA测序试剂盒Oxford Nanopore TechnologiesSQK-RNA002
磁性mRNA分离试剂盒NEBS1550S
T4 DNA连接酶NEBM0202
Superscript IV逆转录酶Thermo Fisher Scientific18090010
RNAClean XP磁珠BeckmanA63987
MinION Mk1B测序仪Oxford Nanopore Technologies--
Chromium Next GEM单细胞3'试剂盒v3.110X Genomics--
地衣芽孢杆菌蛋白酶Creative EnzymesNATE-0633
ReadyCount绿色/红色活力染料InvitrogenA49905
Pierce C18枪头Thermo Scientific87784
孕马血清促性腺激素Lee Biosolutions493-10

关键环节

汇总复现实验时建议重点确认的条件及原文阅读提示。

环节核对要点
动物模型和刺激
小鼠品系(C57BL/6)、年龄(6周龄和14月龄)、PMSG剂量(5 IU)、hCG剂量(5 IU)、给药途径(腹腔注射)、时间(PMSG后46小时,hCG后16小时)
阅读提示:Materials and methods: Ovary collection and RNA isolation
直接RNA测序
RNA投入量(100 ng poly(A) RNA)、试剂盒版本(SQK-RNA002)、测序平台(MinION Mk1B)、flow cell型号(FLO-MIN106D)
阅读提示:Materials and methods: Direct RNA-sequencing of young and aged ovaries
剪接分析
参考基因组版本(GRCm38)、基因注释版本(Ensembl release 102)、FLAIR设置(严格剪接检查)、贝叶斯模型先验(gamma shape=5, scale=1)、显著性阈值(95% HPD区间不含0,后验概率>0.95)
阅读提示:Materials and methods: Alignment and isoform quantification; Bayesian model for differential isoform usage
蛋白结构域分析
Pfam数据库版本(release 35.0)、hmmsearch参数
阅读提示:Materials and methods: Protein domain annotation
scRNA-seq验证
细胞数(目标5000细胞/文库)、文库试剂盒(10X Genomics v3.1)、测序平台(NextSeq 1000/2000)、读长设置(210,10,10,90)
阅读提示:Materials and methods: scRNA-seq; Paired-end scRNA-seq isoform usage validation
质谱验证
卵巢数量(每组10个)、蛋白酶(trypsin Promega V511A)、LC-MS/MS仪器(EASY nLC 1200, Orbitrap Eclipse)、梯度时间(140分钟)、MaxQuant版本和设置
阅读提示:Materials and methods: Sample preparation for mass spectrometry; Mass spectrometry; MS data analysis