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

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统计如下表,

communityOTU1OTU2OTU3OTU4OTU5
A108411
B73840

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

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.hqwc.cn/news/96086.html

如若内容造成侵权/违法违规/事实不符,请联系编程知识网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

如何实现的手机实景自动直播,都有哪些功能呢?

手机实景自动直播最近真的太火了&#xff0c;全程只需要一部手机&#xff0c;就能完成24小时直播带货&#xff0c;不需要真人出镜&#xff0c;不需要场地&#xff0c;不需要搭建直播间&#xff0c;只需要一部手机就可以了。真人语音讲解&#xff0c;真人智能回复&#xff0c;实…

力扣:随即指针138. 复制带随机指针的链表

复制带随机指针的链表 OJ链接 分析&#xff1a; 该题的大致题意就是有一个带随机指针的链表&#xff0c;复制这个链表但是不能指向原链表的节点&#xff0c;所以每一个节点都要复制一遍 大神思路&#xff1a; ps:我是学来的 上代码&#xff1a; struct Node* copyRandomList(s…

不同路径【动态规划】

不同路径 一个机器人位于一个 m x n 网格的左上角 &#xff08;起始点在下图中标记为 “Start” &#xff09;。 机器人每次只能向下或者向右移动一步。机器人试图达到网格的右下角&#xff08;在下图中标记为 “Finish” &#xff09;。 问总共有多少条不同的路径&#xff1f;…

四轴飞行器的电池研究(MatlabSimulink仿真)

&#x1f4a5;&#x1f4a5;&#x1f49e;&#x1f49e;欢迎来到本博客❤️❤️&#x1f4a5;&#x1f4a5; &#x1f3c6;博主优势&#xff1a;&#x1f31e;&#x1f31e;&#x1f31e;博客内容尽量做到思维缜密&#xff0c;逻辑清晰&#xff0c;为了方便读者。 ⛳️座右铭&a…

QT(8.30)常用类与组件,实现登录界面

1.作业&#xff1a; 完成一个登录界面(图片未附带): 头文件: #ifndef WIDGET_H #define WIDGET_H#include <QWidget>#include <QLineEdit>//行编辑器#include<QIcon>//图标#include<QLabel>//标签#include<QPushButton>//按钮#include<QIc…

Python|小游戏之猫捉老鼠!!!

最近闲(mang)来(dao)无(fei)事(qi)&#xff0c;喜欢研究一些小游戏&#xff0c;本篇文章我主要介绍使用 turtle 写的一个很简单的猫捉老鼠的小游戏&#xff0c;主要是通过鼠标控制老鼠(Tom)的移动&#xff0c;躲避通过电脑控制的猫(Jerry)的追捕。 游戏主体思考逻辑&#xff1…

WoW GM

当年黑翼被人黑G以后&#xff0c;后来我就自己开团&#xff0c;今天整理电脑还发现截图。。。。

第2章 Linux多进程开发 2.19-2.23 与信号相关的函数

软件中断&#xff1a; 信号名称&#xff08;宏不会变&#xff09;&#xff0c;编号在不同的版本可能会变&#xff1a; core保存进程异常退出的信息&#xff1a; core.c: #include <stdio.h> #include <string.h>int main() {char * buf; // 指针 没有指…

手机怎么剪视频?分享一些剪辑工具和注意事项

视频剪辑是一种将多个视频片段进行剪切、合并和编辑的技术&#xff0c;它可以帮助我们制作出精彩的视频作品。如今&#xff0c;随着智能手机的普及&#xff0c;我们可以随时随地使用手机进行视频剪辑。本文将为大家介绍一些手机剪辑工具和注意事项&#xff0c;帮助大家更好地进…

代码随想录训练营第四十三天|1049. 最后一块石头的重量 II、 494. 目标和、 474.一和零

1049. 最后一块石头的重量 II 力扣题目链接(opens new window) 题目难度&#xff1a;中等 有一堆石头&#xff0c;每块石头的重量都是正整数。 每一回合&#xff0c;从中选出任意两块石头&#xff0c;然后将它们一起粉碎。假设石头的重量分别为 x 和 y&#xff0c;且 x &l…

MATLAB制图代码【第二版】

MATLAB制图代码【第二版】 文档描述 Code describtion: This code is version 2 used for processing the data from the simulation and experiment. Time : 2023.9.3 Author: PEZHANG 这是在第一版基础上&#xff0c;迭代出的第二版MATLAB制图代码&#xff0c;第二版的特点是…

纵行科技与山鹰绿能达成合作,提供物联网资产管理数据服务

近日&#xff0c;纵行科技与山鹰绿能宣布双方达成深度合作关系&#xff0c;纵行科技将为山鹰绿能提供专业的物联网技术服务&#xff0c;使用物联网技术帮助山鹰绿能对循环包装载具等资产进行在线管理和数字化运营。 据悉&#xff0c;山鹰绿能是一家由山鹰国际控股的全资子公司…