构建并训练PAL-AI模型预测蛙卵母细胞尾长度变化
建立从mRNA序列预测poly(A)尾长度变化的机器学习模型
基于k-mer特征的线性回归模型和整合卷积与循环神经网络的PAL-AI模型，使用蛙内源mRNA尾长度变化数据进行10折交叉验证训练和测试
关于丁香通
公司信息
个人用户
企业机构
无忧采购 轻松科研
- 移动端
PAL-AI揭示控制卵母细胞成熟过程中poly(A)尾长度的遗传决定因素，并与人类生育力相关
PAL-AI reveals genetic determinants that control poly(A)-tail length during oocyte maturation, with relevance to human fertility.
实验完整度
高包含体外卵母细胞注射模型、体内外功能验证、单核苷酸突变实验及人类基因组数据分析等多层证据
主要模型
非洲爪蟾卵母细胞 小鼠卵母细胞 人类卵母细胞
重点核对
非洲爪蟾卵母细胞注射mRNA文库（N60(LC)-PASmos）及单核苷酸突变文库，孕酮处理0 h和7 h后测尾长 PAL-AI模型在蛙内源mRNA上10折交叉验证，Pearson Rp=0.82 单核苷酸突变效应预测与实验测量比较，10个mRNA中6个Rp=0.78–0.90 小鼠和人类卵母细胞数据集验证，PAL-AI-ms模型Rp分别为0.64和0.66 All of Us Research Program和gnomAD中预测破坏尾延长的变异等位基因频率耗竭分析
摘要
在哺乳动物和其他动物的卵母细胞中，基因调控主要通过poly(A)尾长度的变化介导。在此，我们介绍PAL-AI，一种集成的神经网络机器学习模型，能够准确预测青蛙和哺乳动物成熟卵母细胞中尾长度的变化。我们表明，PAL-AI学习到了已知和先前未知的序列元件及其控制poly(A)尾长度的上下文特征，使其能够预测由3′非翻译区单核苷酸替换引起的尾长度变化。它还预测了由尾长度介导的翻译变化，使我们能够提名对卵母细胞成熟重要的基因。当我们将预测的人类卵母细胞尾长度变化与All of Us Research Program和gnomAD的基因组数据集进行比较时，我们发现预测会破坏尾延长的遗传变异在人群中一直受到负选择，从而将mRNA尾延长与人类女性生育力联系起来。
实验结论
提炼研究问题、关键发现与证据，快速把握文章的核心贡献。
- 研究问题
- 能否开发机器学习模型准确预测蛙和哺乳动物卵母细胞成熟过程中poly(A)尾长度变化，并评估遗传变异对尾长度控制及人类生育力的影响？
- 核心机制
- PAL-AI学习到CPE、PAS及其上下文特征（位置、间距、侧翼核苷酸）和UGU/GUU辅助基序控制胞质多聚腺苷酸化，预测的尾长度变化与翻译效率变化相关。
- 主要证据
- PAL-AI在蛙内源mRNA上解释>67%方差（Rp=0.82）；在注射mRNA文库和单核苷酸突变文库中预测与实测相关（Rp最高0.90）；多物种模型预测鼠和人卵母细胞尾变化；人类遗传变异分析显示破坏尾延长的变异在常见等位基因中耗竭。
- 研究意义
- 将mRNA尾延长与人类女性生育力联系起来，提示破坏尾长度控制的3′UTR变异可能是女性不孕的遗传原因之一，PAL-AI可作为提名因果变异的工具。
研究路径
按研究推进顺序梳理实验设计、验证步骤与关键观察。
体外验证PAL-AI对注射mRNA尾长度变化的预测
检验模型能否预测外源注射mRNA报告基因的尾长度变化
将N60(LC)-PASmos mRNA文库注射到蛙卵母细胞，孕酮处理0 h和7 h后收集RNA，高通量测序测定尾长分布，与PAL-AI预测比较
单核苷酸突变文库验证PAL-AI预测突变效应
评估PAL-AI预测单核苷酸替换对尾长度影响的能力
构建包含10个蛙mRNA 3′UTR最后100 nt所有单核苷酸替换的文库（3000变体+野生型），注射蛙卵母细胞，测序测定尾长变化，与预测比较
多物种PAL-AI模型预测鼠和人卵母细胞尾长度变化
检验PAL-AI能否跨物种预测哺乳动物卵母细胞尾长度变化
使用蛙、鼠、人卵母细胞数据训练PAL-AI-ms模型，在鼠GV-MII和人GV-MII卵母细胞数据集上测试预测性能
分析人类遗传变异预测效应对进化选择和等位基因频率的影响
评估预测破坏尾延长的遗传变异是否在人类群体中受到负选择
使用PAL-AI-ms预测人类3′UTR所有可能单核苷酸突变的尾长度变化，与phyloP保守性评分比较；分析All of Us Research Program和gnomAD中变异等位基因频率与预测效应的关系
研究方法
按研究目的归类文中使用的方法，便于定位所需技术。
验证机制链路
RNA结构预测
产品清单
| 实验环节 | 名称 | 品牌 | 货号 |
|---|---|---|---|
| mRNA文库制备 | KAPA HiFi HotStart试剂盒 | Roche | KK2502 |
| 凝胶纯化 | GeneJet凝胶回收试剂盒 | Thermo Fisher | K0692 |
| 体外转录 | SUPERase·In RNA酶抑制剂 | Thermo Fisher | AM2694 |
| T7 RNA聚合酶 | -- | -- | |
| DNA模板去除 | DNase I | New England Biolabs | M0303S |
| RNA加帽 | 痘苗病毒加帽系统 | New England Biolabs | M2080S |
| RNA纯化脱盐 | Micro Bio-Spin P-30柱 | Bio-Rad | 7326250 |
| 去除2′,3′-环磷酸 | T4多核苷酸激酶 | New England Biolabs | M0201S |
| 卵母细胞培养 | 庆大霉素 | Thermo Fisher | 15750060 |
| 卵母细胞注射 | PLI-100 Plus显微注射仪 | -- | -- |
| 卵母细胞成熟诱导 | 孕酮 | Millipore Sigma | P0130 |
| 卵母细胞裂解 | cOmplete蛋白酶抑制剂混合物 | MilliporeSigma | 11836170001 |
| RNA提取 | TRIzol LS | Thermo Fisher | 10296010 |
| Phasemaker管 | Thermo Fisher | A33248 | |
| mRNA富集 | MyOne链霉亲和素C1磁珠 | Thermo Fisher | 65002 |
| 3′接头连接 | T4 RNA连接酶1 | NEB | M0204S |
| 逆转录 | SuperScript IV逆转录酶 | ThermoFisher | 18090050 |
| PCR产物纯化 | AMPure XP磁珠 | Beckman Coulter | A63881 |
| 单核苷酸突变文库合成 | SurePrint HiFi寡核苷酸 | Agilent Technologies | -- |
关键环节
汇总复现实验时建议重点确认的条件及原文阅读提示。
| 环节 | 核对要点 |
|---|---|
| 蛙卵母细胞注射mRNA文库 | 卵母细胞来源（Xenopus1, 12005）、发育阶段V–VI、注射量4 nl（0.1 pmol/μl）、孕酮浓度10 μg/ml、收集时间0 h和7 h、每组100个卵母细胞 阅读提示：Methods中Reporter mRNA library injection and sample collection小节 |
| 单核苷酸突变文库实验 | 10个mRNA 3′UTR最后100 nt、3000变体+野生型、注射条件同上、测序数据分析中尾长起点位置101 阅读提示：Methods中Preparation of the single-nucleotide mutagenesis library小节及Variant and tail-length analysis小节 |
| PAL-AI模型训练与评估 | 10折交叉验证、输入区域最后2000 nt 3′UTR加编码区、超参数（学习率0.0025169等）、每个时间点至少50个poly(A)标签 阅读提示：Methods中PAL-AI for poly(A) tail-length changes小节及Supplementary Data 1 |
| 人类遗传变异分析 | All of Us Research Program外显子callset v8、gnomAD v4.1、仅分析3′UTR最后100 nt、变异效应定义为与参考等位基因预测尾长变化的差值 阅读提示：Methods中Prediction of poly(A) tail-length changes for human genetic variants小节 |