单细胞中转录适应的患病率及基因调控约束

Prevalence of and gene regulatory constraints on transcriptional adaptation in single cells

作者信息Ian A Mellis, Madeline E Melzer, Nicholas Bodkin, Yogesh Goyal
PMID39135102
发布时间2024-08-12
DOI10.1186/s13059-024-03351-2

实验完整度

中

整合多组学公共数据集分析、调控子表达稳健性分析、随机数学模型模拟及机器学习分析,但缺乏功能性实验验证(如基因操作及表型检测),主要为计算和建模研究。

主要模型

小鼠胚胎干细胞(mESC) 人源肿瘤细胞系(Perturb-seq数据) 基因调控网络随机数学模型

重点核对

Bulk RNA-seq数据中对照和CRISPR敲除样本的重复数及实验设计(Methods: Selection of CRISPR-Cas9 transcriptomics datasets) 旁系同源物鉴定方法及Ensembl注释版本(Methods: Identification of paralogs) 差异表达阈值:adjusted p-value ≤ 0.05 和 log2 fold-change ≥ 0.5 Perturb-seq数据质控过滤条件:细胞数、UMI数及guide RNA数量(Methods: Perturb-seq reanalysis) 数学模型参数范围及模拟条件(如Gillespie模拟、时间步长)

摘要

细胞和组织通过多种分子机制对遗传扰动具有显著的适应能力。无义突变诱导的转录补偿作为转录适应的一种形式,最近成为一种机制,其中基因的无义突变触发相关基因的上调,可能在细胞和有机体水平上赋予稳健性。然而,除了少数发育背景和精选基因集外,尚未对哺乳动物细胞类型和条件进行全面的全基因组研究。固有随机补偿基因网络的调控效应如何影响单细胞中的表型外显率仍不清楚。我们分析现有的批量及单细胞转录组数据集,以揭示哺乳动物系统中跨不同背景和细胞类型转录适应的普遍性。我们在批量及混合单细胞遗传扰动数据集中对转录因子靶标集进行调控子基因表达分析。我们的结果显示,表现出转录适应的转录因子调控子的表达比没有表现出转录适应的转录因子调控子更具稳健性。最小补偿基因网络的随机数学建模定性地重现了转录适应的几个方面,包括旁系同源物的上调和突变的稳健性。结合网络特征的机器学习分析,我们的框架为哪些调控步骤对转录适应最重要提供了可能的解释。我们的整合方法确定了几种假定的候选基因——表现可能转录适应的基因——用于实验验证,并提供了一个正式的定量框架来测试和精炼转录适应模型。

实验结论

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

研究问题
转录适应(无义诱导的转录补偿)在哺乳动物细胞中的普遍性及其基因调控约束是什么?
核心机制
无义突变产生提前终止密码子,触发无义介导的mRNA降解,进而导致旁系同源基因表达上调(转录适应),维持下游调控子表达的稳健性。
主要证据
对批量RNA-seq数据中74个CRISPR敲除靶标分析发现16个靶基因的旁系同源物上调;Perturb-seq数据中85个旁系同源物-靶基因对进入top-100;转录因子调节子的表达在具有转录适应的靶标中波动更小;随机数学模型模拟重现了旁系同源物上调和表达分布稳健性。
研究意义
该研究提供了一个分析转录适应普遍性的框架,提示在解释CRISPR敲除表型时需考虑转录适应的影响,并有助于设计实验验证和预测哪些基因可能发生转录适应。

研究路径

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

1

公共CRISPR-Cas9转录组数据集筛选与差异表达分析

确定转录适应在批量转录组数据中的普遍性。

从GEO/GREIN筛选公开的CRISPR-Cas9敲除伴随RNA-seq数据,获取对照和敲除样本的转录组,使用DESeq2计算差异表达。

2

旁系同源物鉴定及上调频率分析

评估敲除后旁系同源物上调是否显著高于随机。

使用Ensembl注释鉴定敲除靶标的旁系同源物,计算旁系同源物上调比例,并通过引导抽样构建零分布评估显著性。

3

跨条件一致性分析

验证同一基因的旁系同源物上调在不同条件下是否一致。

使用Lackner et al. 2021的mESC数据,比较在2i培养和分化条件下旁系同源物上调频率及倍数变化的关联性。

4

Perturb-seq单细胞CRISPR筛选分析

在单细胞水平验证旁系同源物上调。

重新分析Frangieh et al. 2021的Perturb-seq数据,比较靶向guide与对照guide处理细胞中旁系同源物表达阳性率或平均表达的差异。

5

转录因子调控子稳健性分析

评估转录适应是否缓冲下游调控子的表达变化。

鉴定转录因子及其旁系同源物,查询DoRothEA数据库中的共同调控子,比较敲除后调控子基因的表达变化幅度。

6

随机数学模型模拟

通过建模重现转录适应中旁系同源物上调及下游效应物表达分布的稳健性。

构建基于电报模型的随机基因调控网络,模拟基因A突变后旁系同源物A'上调,并分析下游效应物B的表达分布形状变化。

研究方法

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

产品清单

实验环节名称品牌货号
PyDESeq2----
clusterProfiler----
BiomaRt----
MATLABMathWorks--
R语言----

关键环节

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

环节核对要点
数据集选择与差异表达分析
CRISPR靶基因、对照设计、重复数、测序深度、表达过滤标准、差异表达阈值
阅读提示:Methods: Selection of CRISPR-Cas9 transcriptomics datasets; Differential gene expression assessment
旁系同源物鉴定
旁系同源物注释来源(Ensembl版本)、序列同源性标准
阅读提示:Methods: Identification of paralogs
Perturb-seq分析
细胞数量、UMI阈值、guide RNA数量、对照处理、表达量化和百分阳性计算
阅读提示:Methods: Perturb-seq-based single-cell gene expression reanalysis
调控子稳健性分析
转录因子定义(AnimalTFDB3版本)、调控子数据库(DoRothEA)、置信度等级、差异表达阈值
阅读提示:Methods: Regulon robustness analysis
随机数学模型模拟
参数范围、模型类型(激活/抑制/多旁系同源物)、模拟算法(Gillespie)、时间步长、基因型(野生型/杂合/纯合)
阅读提示:Methods: Simulations; Parameters; Networks