之前我们平台分享过《Nature:宏基因组关联分析》综述,让大家系统的了解这一领域。同时还分享过一篇 《GigaScience:谷子产量与微生物组关联分析》是植物领域中的优秀工作。今天再分享一篇在动物中的研究,供大家学习。值得一提的是此文也是微生物所王军老师早年参与的文章之一。

基于易感自免疫皮肤水疱小鼠全基因组关联基因与微生物互作

Genome-wide mapping of gene–microbiota interactions in susceptibility to autoimmune skin blistering

此文2013年发表于Nature Communication上,截止2018年1月6日Google Scholar统计己被引49次。

摘要

易感慢性感染疾病,通常于自免疫遗传因素,和环境风险因子共同决定。在疾病与健康人群中微生物群落通常存在差异,然后这些差异是原因,还是疾病的结果仍然末知。
本文,基于小鼠获得性大疱性表皮松解模型研究,发现了宿主基因与微生物组互作对疾病发生存在贡献。
我们鉴定到的基因座位对皮肤微生物组的改变和易感皮肤水泡状态存在贡献,并且两者间存在共同的基因。
进一步,把细菌物种作为疾病的共变量,我们鉴定了一个新的病症位点。主要的共变量物种类别丰度下降伴随着疾病风险增加,提供了其在防病中的重要作用。这些找到的益生菌及群落提供了预防和治疗此类疾病的可能手段。

主要结果

图1. Chao1物种丰富度。Chao1指数是基于种水平OTUs估计,其中健康119个,EBA 64个样品,误差线代表95%置信区间,显著性为Wilcoxon rank-sum test.

图2. 监督式PCoA分析Jaccard距离。以疾病状态为限制因子,显著性统计基于permutation test,图中了各组的重心和68%的置信区间椭圆?

图3. QTL定位皮肤微生物。19个小鼠常染色体和X染色体添加1199个SNPs(黑线)、同时鉴定了种和属水平的QTLs并用门水平上色。同时还标注了两篇之前研究中报导的QTLs位点。

图4. EBA感染中基因与微生物互作
(a) 曼哈顿图展示1199个SNP与表型与疾病表型-log10 P-values沿每条染色体上的分布;(b) 添加Staphylococcus spp. (OTUID 173469)丰度作为共因子分析关联;(C) rs6211533位点添加共变量后出现显著位点;(D) 动物在Staphylococcus spp.丰度高低分为两组,发现率也存在显著差异;三种颜色代表三种基因型。

图5. 协变量OTUs间的共相关矩阵。OTUs间的皮尔森相关系数与EBA易感座位显著共变化。对于OTUs属水平无名,则显示分类学最高级别的分类并添加星,只有显著的才显示颜色。

方法

DNA提取和16S测序

使用PowerSoil Kit提取小鼠耳朵提取细菌DNA。扩增使用高变区V1-V2(27F AGAGTTTGATCCTGGCTCAG - 338R TGCTGCCTCCCGTAGGAGT)

物种注释

RDP分类器使用0.8的最小置信值,fixrank参数要求分界到属水平。种水平注释采用blast比对 Greengene并阈值为E-value < 0.001。

QTL分析数据的准备

在每个样品中reads先转换为比便,再进行log10对数变换。在OTUs、属、科、目、纲和门水平进行合并汇总(summary_tax.py)。

微生物组核心鉴定

筛选bin中reads数量大于20,且至少在20个动物中出现的OTUs,找到131个OTUs占总数据量的80%。

基因分型和QTL分析

使用芯片鉴定了261个动物个体的1449个SNP,其中1199个有效SNP。Happy包用于QTL和协变量QTL分析,期待存在座位的加性贡献。每个log10转换的分类学bins作为一个独立的表型,使用happy包的hfit命令去校正笼子、家族效应。在输入文件中添加协变量的矩阵。全基因组显著采用1000次排列表型,采用QQ-plot展示log10 P-value,

Reference

  1. Srinivas, G., et al. (2013). “Genome-wide mapping of gene-microbiota interactions in susceptibility to autoimmune skin blistering.” Nat Commun 4: 2462.
  2. Belheouane, M., et al. (2017). “Improved detection of gene-microbe interactions in the mouse skin microbiota using high-resolution QTL mapping of 16S rRNA transcripts.” Microbiome 5(1): 59.

QTL分析

QTL分析软件 happy,目前是官网 http://www.well.ox.ac.uk/happy 无法打开;Rbioconductor也没有此包。

http://www.rqtl.org/ 此网站有分析QTL分析软件 大全。推荐使用qtl,一直在更新

# 稳定版
install.packages("qtl")
library(qtl)# 开发版
library(devtools) install_github("kbroman/qtl@devel")

猜你喜欢

  • Top热文:1提DNA发Nature 2论文图表规范 3肠道菌群 4生信宝典 5实验vs分析 6我们的未来在哪里
  • 一文读懂:1微生物组 2进化树 3网络图 4图片配色
  • 必备技能:1提问 2搜索 3Endnote 4在线绘图
  • 文献阅读 1geenmedical 2热心肠 3SemanticScholar
  • 扩增子分析:1图表解读 2分析流程 3统计绘图 4美化进化树ggtree
  • 16S预测功能 1KEEG通路PICRUSt 2元素循环FAPROTAX 3表型bugbase 4KO通路Tax4Fun
  • 宏基因组分析:1拼接MEGAHIT 2基因注释Prokka 3Binning 4MetaPhlAn2
  • 科研经验:1云笔记 2云协作 3公众号 4微生物组入门 5成为团队VIP
  • 系列教程:1生信入门Biostar 2微生物组 3宏基因组
  • 生物科普 1生命大跃进 2细胞的暗战 3人体奥秘 4人体内旅行 5神奇细胞

写在后面

宏基因组公众号经过半年发展,原创文章过百,关注人数12000+。为鼓励读者交流、快速解决科研困难,我们建立了“宏基因组”专业讨论群,目前己有国内外100+ PI,900+ 一线科研人员加入。参与讨论,获得专业解答,欢迎分享此文至朋友圈,并扫码加主编好友带你入群,务必备注“姓名-单位-研究方向-职称/年级”。技术问题寻求帮助,首先阅读《如何优雅的提问》学习解决问题思路,仍末解决群内讨论,问题不私聊,帮助同行。

学习扩增子、宏基因组科研思路和分析实战,关注“宏基因组”

点击阅读原文,跳转最新文章目录阅读
https://mp.weixin.qq.com/s/5jQspEvH5_4Xmart22gjMA

NC:自体免疫水泡皮肤病中鉴定基因与微生物组互作(微生物组关联分析MWAS)相关推荐

  1. mLife | 刘星吟等肠道菌群、基因变异和饮食互作介导孤独症的发生

    南京医科大学基础医学院刘星吟教授的观点文章"The interaction of gut microbiota, genetic variation and diet in Autism S ...

  2. 双代号网络图基础算法_从网络图探寻基因互作的蛛丝马迹(5)

    在前面的4期中,我们分别给大家讲解了网络图的构造. STRING 数据库.Cytoscape 软件的安装以及使用,链接如下: 从网络图探寻基因互作的蛛丝马迹(1) [科研猫·绘图]从网络图探寻基因互作 ...

  3. SBIO | 浙大陈云组综述农业中细菌与真菌的互作机制

    点击蓝字|关注我们 导读  ‍ 细菌和真菌是环境微生物组中的主要成员.细菌-真菌互作(Bacterial-fungal interactions,BFIs)及其调控是影响宿主及生态系统健康和功能的重要 ...

  4. 文献(5): 单细胞转录组揭示人类纤维化皮肤病中成纤维细胞的异质性

    一句话评价:小样本量单细胞转录组在冷门疾病领域应用前景尚可,常规套路也能发NC 1. Title 标题可以看出本研究的主要内容:①构建成纤维细胞的图谱/揭示成纤维细胞的异质性:②找到一群(相对于正常样 ...

  5. 热点综述 | 空间多组学技术在肿瘤免疫微环境研究中的应用

    目前空间多组学技术已用于研究几种类型癌症中肿瘤免疫微环境的转录组.蛋白质组和代谢组,从这些方法获得的数据已与免疫组化和多参数分析相结合,以产生癌症进展的标记物.2022年<Journal of ...

  6. MPB:微生物所东秀珠组-​​基于16S rRNA基因和基因组序列对细菌物种的初步鉴定...

    为进一步提高<微生物组实验手册>稿件质量,本项目新增大众评审环节.文章在通过同行评审后,采用公众号推送方式分享全文,任何人均可在线提交修改意见.公众号格式显示略有问题,建议电脑端点击文末阅 ...

  7. metaProdigal:宏基因组序列中的基因和翻译起始位点预测

    文章目录 metaProdigal:宏基因组序列中的基因和翻译起始位点预测 热心肠日报 摘要 动机 Motivation 结果 Results 可用性 Availability 主要结果 表1. 大肠 ...

  8. 易基因|3文聚焦:宏病毒组测序在肠病中的应用研究

    大家好,这里是专注表观组学十余年,领跑多组学科研服务的易基因.本期我们通过3篇高分文章来解读宏病毒组测序在肠病中的应用研究. 01 小鼠结肠炎研究揭示噬菌体在疾病发生中的作用 标题:Murine co ...

  9. RNA 13. SCI 文章中加权基因共表达网络分析之 WGCNA

    WGCNA 分析流程 2008 年发表在 BMC 之后的影响力还是很高的,先后在各大期刊都能看到,但是就其分析的过程来看,还是需要有一定 R 语言的基础才能完整的复现出来文章中的结果,这期就搞出来供大 ...

最新文章

  1. python—多线程之线程之间共享数据(Queue)
  2. 尝试做了个Chrome浏览器的扩展程序,显示extensions的页面
  3. Android之6.0 权限申请封装
  4. nginx php访问日志配置,nginx php-fpm 输出php错误日志的配置方法
  5. python如何用色度表示数值大小_python入门004数字(例程界面很清晰大小也合适)
  6. MySQL安装目录分类_分类 安装 下的文章
  7. mysql存储过程中的异常处理
  8. MQ通道搭建以及连通性检查
  9. 【优化算法】混沌游戏优化 (CGO)【含Matlab源码 1010期】
  10. C1WebChart 图形化处理。
  11. 在linux下运行锐捷客户端,锐捷Linux客户端使用方法(基于ubuntu16.04)
  12. 一个门外汉的产品设计漫谈:从初级到神级的设计
  13. 腾跃英语计算机学院微信公众号,英语四级报名_微信还能这么玩:Geek大学生搭建英语课堂互动系统_沪江英语...
  14. 利用计算机发布调度命令时必须严格遵守,调度命令规范格式(22页)-原创力文档...
  15. 一年级计算机上册计划进度表,一年级上册语文教学计划及进度表
  16. raphael.js菜鸟笔记(二)
  17. 【MSFconsole进阶】evasion(免杀模块):免杀模块使用方法,相关设置
  18. Redis中的Stream的实现Radix Tree源码解读
  19. Android DataBinding 详解
  20. Js isNaN()和Number.isNaN()、isFinite() 和 Number.isFinite()

热门文章

  1. 2019了,互联网研发是不是都快要被淘汰了?
  2. 调用训练好的模型(tensorflow)
  3. winsock select模型实现
  4. SQL与NoSQL区别-规范化
  5. android 动态 修改图标,android-动态更改ActionBar中的图标
  6. 双链表(删除节点操作)
  7. 二叉搜索树-创建最小高度树(递归)
  8. 31 元素滚动scroll系列
  9. Python3多线程threading介绍(转载)
  10. 18岁智商低的表现_吃手是宝宝聪明的信号?婴儿智商高的5个讯号,吃手只是其中一个...