[发明专利]一种构建水平基因转移网络的方法

说明书

本发明公开了基因转移技术领域的一种构建水平基因转移网络的方法，包括以下步骤：S1：基于检测到的HGT事件构建初始物种HGT网络；S2：根据边的权重，对边进行过滤；S3：合并初始物种HGT网络的节点；HGT网络可通过一定的方法消除测序量对HGT网络的影响，使结果更加准确：测序量会影响HGT网络，但HGT网络构建之后，可通过在每个样本中只保留若干条最可靠边的手段来消除测序量的影响，其中边的可靠性是通过在两个分类群之间发生HGT中支持reads的总数来量化的，在实验中若对样本使用本方法，可以在每个样本中生成不同分类级别的HGT网络，来进行不同水平的分析，并且HGT网络可以反映微生物之间的基因交流和动态变化，可以提取更多生物学特征和信息。

技术领域

本发明涉及基因转移技术领域，具体为一种构建水平基因转移网络的方法。

背景技术

人体微生物组成了一个复杂的生态系统，其中的成员有着复杂的相互作用。在分析和理解人体微生物组的复杂系统的有力工具中，基于网络的数学模型可用于分析和量化这类系统。一种方法是计算每对微生物丰度的相关系数，形成丰度网络。正相关和负相关分别表明这两种微生物可能存在合作或竞争相互作用。

另一种方法是基于微生物之间的代谢交换来构建人类肠道微生物组的全局代谢相互作用网络。代谢网络可以分基因组(DNA层次)、代谢途径及生化反应网络(蛋白质层次)、代谢流(物流层次)、代谢生理(微生物细胞层次)等层次来讨论。

第三种即为分子生态网络(molecular ecological networks)。它由基因间的皮尔逊相关系数(PCC)构建。

现有方法不能捕捉到物种间的基因交流和动态变化，为此，我们提出一种构建水平基因转移网络的方法。

发明内容

本发明的目的在于提供一种构建水平基因转移网络的方法，以解决上述背景技术中提出的问题。

为实现上述目的，本发明提供如下技术方案：一种构建水平基因转移网络的方法，包括以下步骤：

S1：基于检测到的HGT事件构建初始物种HGT网络；

S2：根据边的权重，对边进行过滤；

S3：合并初始物种HGT网络的节点；

S4：建立HGT事件网络。

优选的，所述S1中检测得到得HGT事件包含五个信息：donor基因组，donor基因组上断点位置，receptor基因组，receptor基因组上断点位置，支持该事件的测序读数数目。

优选的，所述S1中每个断点为一个节点，节点间的边即为支持该断点对的HGT事件，边的权重即为支持该断点对的读数数目。

优选的，所述S2中根据贝叶斯模型，把权重低于一定cutoff的边删去。

优选的，把边删去后，对节点进行过滤。

优选的，计算每个节点的度，并把度数为0的点删去，得到过滤的初始物种HGT网络。

优选的，所述S3中根据不同的层次合并节点，把他们对应的边和权重也相应的合并，得到不同水平的物种HGT网络，对节点进行种、属、科等不同分类层次的标注，并对具有相同分类单元的节点进行合并，通过这些方式，这样就可以在每个样本中生成具有不同分类级别的物种HGT网络。

优选的，所述S3中合并方法方为两种：

第一种方法对基因组划分成等长的window，把属于相同window的节点合并。

第二种方法根据分类水平合并节点。

优选的，对HGT事件定义为不同的分类层次，得到具有不同分类级别的HGT事件网络。