植物根系基因组与数据分析

这篇具有很好参考价值的文章主要介绍了植物根系基因组与数据分析。希望对大家有所帮助。如果存在错误或未考虑完全的地方,请大家不吝赐教,您也可以点击"举报违法"按钮提交疑问。

1.背景

这段内容主要是关于植物对干旱胁迫的反应,并介绍了生活在植物体内外以及根际的真菌和细菌的作用。然而,目前对这些真菌和细菌的稳定性了解甚少。作者通过调查微生物群落组成和微生物相关性的方法,对农业系统中真菌和细菌对干旱的抗性和恢复力进行了研究。研究结果显示,在微生物群落组成层面上,真菌比细菌更抗旱,但在重新灌溉缓解干旱压力后,真菌的恢复能力较低。然而,在所有相关性和共存网络方面,结果更加复杂。总体上,干旱会破坏细菌、真菌以及细菌与真菌之间显著正相关的微生物网络。令人惊讶的是,干旱加强了根际真菌与叶片细菌之间以及根际中丛枝菌根真菌相关性网络的结构。此外,研究还支持了应激梯度假说,因为干旱增加了正相关的相对频率

2.假设

分析H1和H2假设在微生物群落组成层面上的测试结果。

  • 假设1(H1):认为真菌群落比细菌群落更抗旱;
  • 假设2(H2),认为真菌群落在重新灌水后的恢复能力比细菌群落更低。

3.实验

为了测试这两个假设,研究人员从同一片叶子、根系、根际和土壤样本中收集的高粱植株中获得了真菌和细菌的数据集,并进行了混合处理。这些植株分为两组,一组作为对照组进行正常灌溉,另一组在开花前经历了干旱处理,然后在开花后定期进行灌水。通过对这些数据集的分析,研究人员测试了H1和H2假设在微生物群落组成层面上的有效性。

4.专业知识

抗扰性(resistance)和复原力(resilience)

分析中的两个重要指标:抗扰性(resistance)和复原力(resilience)

  • 抗扰性指的是社区在受到干扰后保持不变的程度
  • 复原力则是社区在受到干扰后恢复到原始状态的速度。

真菌与细菌间的相互作用

微生物群落中的关联对生态系统稳定性的重要性。这些关联可以是正向或负向的,强弱不一,并涉及到细菌、真菌之间的相互作用。
有研究表明正向关联的比例与演替、施肥和栖息地等生态因素相关。

5.分析工具与特征

分析工具:Bray-Curtis dissimilarities
特征:根、根际、土壤和叶子等四个部分的细菌(红色和浅红色)和真菌(蓝色和浅蓝色)群落的差异。

6.R语言使用Bray-Curtis

1.安装包vegan,可以使用镜像源加速

代码

 install.packages("vegan", repos="https://mirrors.tuna.tsinghua.edu.cn/CRAN/")

运行案例:

植物根系基因组与数据分析,生物信息学,数据分析,数据挖掘

2.导入包vegan

代码

library(vegan)

运行案例

植物根系基因组与数据分析,生物信息学,数据分析,数据挖掘

3.创建一个示例数据集

代码

data <- data.frame(
 Sample1 = c(5,0,0,0),
 Sample2 = c(0,10,0,0),
 Sample3 = c(0,0,15,0),
 Sample4 = c(0,0,0,20)
 )

运行案例

植物根系基因组与数据分析,生物信息学,数据分析,数据挖掘

3.计算Bray-Curtis距离

代码

bray_curtis_distance <- vegdist(data,method="bray")

运行案例

植物根系基因组与数据分析,生物信息学,数据分析,数据挖掘

4.打印距离矩阵

代码

print(bray_curtis_distance)

运行案例

植物根系基因组与数据分析,生物信息学,数据分析,数据挖掘

7.Bray-Curtis距离介绍

Bray-Curtis距离是以该统计指标的提出者J. Roger Bray和John T. Curtis的名字命名的,主要基于OTUs的计数统计,比较两个群落微生物的组成差异。与unifrac距离,包含的信息完全不一样;相比于jaccard距离,Bray-Curtis则包含了OTUs丰度信息。

bray-curtis距离
其中,S_(A,i)和S_(B,i)表示第i个OTU分别在A群落和B群落中的计数。min表示取两者最小值。
植物根系基因组与数据分析,生物信息学,数据分析,数据挖掘
举例说明:
群落A和群落B的OTU统计如下表,

community OTU1 OTU2 OTU3 OTU4 OTU5
A 10 8 4 1 1
B 7 3 8 4 0

min(S_(A,i),S_(B,i)) = 7+3+4+1+0 = 15

sum(S_(A,i)) = 10+8+4+1+1 = 24
sum(S_(B,i)) = 7+3+4+8+4+0 = 22

D = 1 - 2*15/(24+22) = 0.3478

当然,D值越小越好,表示二者组成差异小。

8.问题与总结

  • 从文章来理解:共存网络分析主要关注微生物个群在扰动下的共振变化,即仅关注显著的正相关关系。
  • 需要加快理解数据的意义以及现有的分析手段

内容引自:Co-occurrence networks reveal more complexity than community composition in resistance and resilience of microbial communities文章来源地址https://www.toymoban.com/news/detail-693759.html

到了这里,关于植物根系基因组与数据分析的文章就介绍完了。如果您还想了解更多内容,请在右上角搜索TOY模板网以前的文章或继续浏览下面的相关文章,希望大家以后多多支持TOY模板网!

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处: 如若内容造成侵权/违法违规/事实不符,请点击违法举报进行投诉反馈,一经查实,立即删除!

领支付宝红包 赞助服务器费用

相关文章

  • liftOver 不同版本基因组文件相互转化

    大家好,我是邓飞。前一段时间有小伙伴在星球提问:想将不同版本的SNP数据合并,不想重新call snp,想把绵羊的V2和V4版本的数据合并,具体来说,是V2转为V4然后与V4合并。 我建议用 liftOver 软件进行处理,并许诺写篇博客介绍一下。 还有小伙伴想把1.2的参考基因组,变为

    2024年02月07日
    浏览(31)
  • 生信步骤|原核生物基因组注释--Prokka

    全基因组注释是鉴定生物基因组特征的过程。Prokka是一个适用于原核生物基因组注释工具,可以注释细菌,古菌和病毒基因组。Prokka在预测基因组CDS区域时采用了多种数据库,内置的三个核心数据库包括ISfinder数据库,NCBI细菌抗性数据库和UniprotKB数据库。 此外,prokka内置基因

    2024年02月04日
    浏览(81)
  • 易基因:人类大脑的单细胞DNA甲基化和3D基因组结构|Science

    大家好,这里是专注表观组学十余年,领跑多组学科研服务的易基因。 高通通量表观基因组分析技术可用于阐明大脑中细胞复杂性的基因调控程序。5\\\'-甲基胞嘧啶 (5mCs)是哺乳动物基因组中最常见的修饰碱基,大多数5mCs发生在胞嘧啶-鸟嘌呤二核苷酸(CpGs)上。CG差异甲基化区

    2024年04月17日
    浏览(37)
  • 基因组组装: 3D-DNA 染色体挂载

    本文将介绍基因组组装过程中,如何利用 HiC 测序数据,进行染色体级别基因组的组装。该过程主要利用 Juicer [1] 和 3D-DNA [2] 进行,有关第一步 Juicer 的过程,已经下方的文章中介绍了,本文主要介绍第二步: 3D-DNA 的安装与使用。 目前基因组组装的主要流程是,利用二代或者

    2024年02月13日
    浏览(33)
  • MetaHipMer2 - MHM2超算系统宏基因组短读长序列组装神器的介绍和使用

    berkeleylab / mhm2 / Downloads — Bitbucket 文章: Terabase-scale metagenome coassembly with MetaHipMer | Scientific Reports MetaHipMer (MHM) 是一种从头开始的宏基因组短读组装器。这是版本 2 (MHM2),完全用 UPC++、CUDA 和 HIP 编写,可以在单服务器和多节点超级计算机上高效运行,可以扩展以共同组装 te

    2024年01月19日
    浏览(40)
  • GEO生信数据挖掘(七)差异基因分析

    上节,我们使用结核病基因数据,做了一个数据预处理的实操案例。例子中结核类型,包括结核,潜隐进展,对照和潜隐,四个类别。本节延续上个数据,进行了差异分析。 加载数据 构建差异比较矩阵 计算差异基因指标 #绘制前40个基因在不同样本之间的热图 差异基因分析

    2024年02月07日
    浏览(46)
  • 使用R语言绘制富集条形图,轻松分析基因表达数据

    富集分析(enrichment analysis)是一种生物信息学方法,它可以帮助我们识别基因或其他的生物实体在某个特定的类别中过度表示的趋势。通俗来说,富集分析通过将基因分类到特定的集合中,然后根据基因在集合中的分布和总体分布的比较,来寻找哪些集合与特定的生物过程、

    2024年02月11日
    浏览(44)
  • 数据结构课程设计C/C++版--植物百科数据的管理与分析

    注意:评测不通过请重置代码仓库,重新评测 第1关:增加植物信息 任务描述 本关任务:已知plant.txt的路径为\\\"data_edit/plant.txt\\\",从plant.txt中读取植物的基本信息,创建一个植物信息的链表,基于该链表,实现植物基本信息的增加功能。 编程要求 根据提示,在右侧编辑器补充

    2024年04月11日
    浏览(184)
  • 《数据结构》课程设计(C/C++版):植物百科数据的管理与分析

    目录 第1关:增加植物信息 第2关:删除植物信息 第3关:修改植物信息 第4关:基于顺序表的顺序查找 第5关:基于链表的顺序查找 第6关:基于顺序表的折半查找 第7关:基于二叉排序树的查找 第8关:基于开放地址法的散列查找 第9关:基于链地址法的散列查找 第10关:基于

    2024年04月09日
    浏览(200)
  • GEO生信数据挖掘(六)实践案例——四分类结核病基因数据预处理分析

    前面五节,我们使用阿尔兹海默症数据做了一个数据预处理案例,包括如下内容: GEO生信数据挖掘(一)数据集下载和初步观察 GEO生信数据挖掘(二)下载基因芯片平台文件及注释 GEO生信数据挖掘(三)芯片探针ID与基因名映射处理 GEO生信数据挖掘(四)数据清洗(离群值

    2024年02月07日
    浏览(54)

觉得文章有用就打赏一下文章作者

支付宝扫一扫打赏

博客赞助

微信扫一扫打赏

请作者喝杯咖啡吧~博客赞助

支付宝扫一扫领取红包,优惠每天领

二维码1

领取红包

二维码2

领红包