人类疟原虫基因组被组织成数千个共表达的线性调控单元

The human malaria parasite genome is configured into thousands of coexpressed linear regulatory units

作者信息Chengqi Wang, Justin Gibbons, Swamy R Adapa, Jenna Oberstaller, Xiangyun Liao, Min Zhang, John H Adams, Rays H Y Jiang
PMID33272860
发布时间2020-09
DOI10.1016/j.jgg.2020.08.005

实验完整度

本文主要为基于公共数据的生物信息学分析,未提供具体实验方法细节,属基因组数据分析研究。

主要模型

恶性疟原虫3D7株 红细胞内期发育同步化寄生虫

重点核对

使用的RNA-Seq和ATAC-Seq数据来自八个时间点(5, 10, 15, 20, 25, 30, 45, 40小时)的恶性疟原虫3D7寄生虫 AP2-I结合位点来自ChIP-Seq数据(Santos等,2017) HMM模型的平均后验概率阈值设为0.99 调控单元长度中位数为4.7kb

摘要

人类疟原虫恶性疟原虫在蚊子和人类截然不同的宿主环境中生存,仅拥有有限数量的转录因子。恶性疟原虫基因组中调控元件的性质及其靶基因仍然难以捉摸。在这里,我们发现这种真核寄生虫利用一种高效的方式来最大化利用遗传和表观遗传调控,在血液感染期间形成调控单元(RUs)。位于同一调控单元中的基因往往随时间具有相同的表达模式,并与调控元件处的开放染色质相关。为了精确定义和量化这些调控单元,我们开发了一种新的隐马尔可夫模型,通过整合表达和表观遗传证据,以全基因组方式捕获调控结构。我们成功鉴定了数千个调控单元,并与先前发现进行了交叉验证。我们发现更多参与红细胞(RBC)入侵的基因与PfAP2-I(AP2-I)转录因子位于同一调控单元中,表明AP2-I负责调控红细胞入侵。我们的研究为紧凑的真核基因组提供了一种调控机制,并为恶性疟原虫红细胞内期的体内转录调控提供了新见解。

实验结论

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

研究问题
恶性疟原虫基因组中调控元件的性质及其靶基因是什么?
核心机制
恶性疟原虫基因组被组织成调控单元,其中基因与开放染色质区域协同表达,AP2-I转录因子调控同一调控单元中的红细胞入侵基因。
主要证据
通过整合RNA-Seq和ATAC-Seq数据,使用隐马尔可夫模型识别出2072个调控单元,并验证了AP2-I结合位点与调控单元内基因表达的关联。
研究意义
该研究为理解紧凑真核基因组的转录调控机制提供了新见解,有助于探索针对疟疾的药物靶点。

研究路径

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

1

数据获取与预处理

获取全基因组染色质可及性和转录组数据,为识别调控单元提供基础。

从公共数据库下载恶性疟原虫3D7株八个时间点的ATAC-Seq和RNA-Seq数据,以及AP2-I的ChIP-Seq数据。

2

计算方向性向量

量化每个基因组元件(基因或ATAC峰)的上下游调控偏向。

对每个基因组元件,计算其与上下游元件的表达和染色质可及性信号的Pearson相关性,形成方向性向量。

3

基于HMM的调控单元识别

利用隐马尔可夫模型推断调控单元的位置。

以方向性向量为观测,构建四状态HMM(上游偏向、无偏向、下游偏向、双向),使用Baum-Welch算法估计参数,后验概率≥0.99的连续状态区域被定义为调控单元。

4

调控单元与已有数据的比较验证

评估识别出的调控单元的可靠性。

将HMM识别的调控单元与简单相关性截断法识别的单元比较,并检查与已知共表达基因簇和eQTL的关联。

5

AP2-I靶基因的重新定义

基于调控单元重新确定AP2-I转录因子的靶基因。

将AP2-I结合位点映射到调控单元,将同一调控单元内的基因预测为AP2-I靶基因,并与以往方法比较。

研究方法

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

关键环节

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

环节核对要点
数据获取
数据集的来源和版本
阅读提示:Methods部分“Public data sets analyzed”
调控单元识别
HMM的状态数、后验概率阈值
阅读提示:Methods部分“DV, RU, and boundary calling”
AP2-I靶基因预测
AP2-I结合位点数据的来源和调用标准
阅读提示:Methods部分“Public data sets analyzed”和结果2.4节