Summary
这篇手稿描述了一个详细的标准化协议高通量 16S rRNA 扩增子测序。该协议引入了一个集成的, 统一的, 可行的, 廉价的协议, 从粪便样本收集, 通过数据分析。该协议允许分析大量的样本, 严格的标准和几个控制。
Abstract
人体肠道菌群在保护细胞免遭伤害、处理能量和养分以及促进免疫力方面起着重要作用。从被认为是健康的微生物组成 (dysbiosis) 的偏差可能会损害导致病理状况的重要功能。最近和正在进行的研究工作已被用于描述微生物组成与人类健康和疾病之间的关联。
高通量测序技术的进展使肠道微生物成分的特征得以表征。这些方法包括 16S rRNA 扩增子测序和猎枪测序。16S rRNA 扩增子测序用于分析分类成分, 而猎枪测序则提供有关基因预测和功能注释的附加信息。使用 16S rRNA 基因可变区的目标测序方法的优点是它的成本与猎枪测序相比大大降低。在 16S rRNA 基因中的序列差异被用作微生物指纹, 用来识别和量化个体样本中不同的分类群。
主要的国际努力已登记了 16S rRNA 扩增子测序标准。然而, 一些研究报告了一个常见的变异来源, 由间歇效应引起。为了最小化这种效果, 必须实现示例收集、处理和排序的统一协议。该协议提出了从粪便样本采集开始到数据分析的广泛应用协议的集成。该协议包括无列的直接 pcr 方法, 能够同时处理和 DNA 提取大量粪便样本, 以及 V4 区域的 PCR 扩增。此外, 该协议描述分析管道, 并提供一个脚本使用最新版本的 QIIME (QIIME 2 版 2017.7. 0 和 DADA2)。这一步步协议的目的是引导那些感兴趣的人开始使用 16S rRNA 扩增子测序的健壮, 繁殖, 易于使用, 详细的方式。
Introduction
为了更好地了解微生物的多样性和丰度, 已作出了集中努力, 作为在健康和病理条件下捕捉个体差异和相似性的另一个方面。年龄2,3, 地理4, 生活方式5,6, 疾病5被证明与肠道微生物组的组成有关, 但许多条件和人口尚未完全特点.最近据报道, 微生物菌群可以被修改为治疗性应用7,8,9。因此, 进一步了解各种生理条件和微生物组成之间的关系, 是迈向优化未来潜在的修改的第一步。
传统的微生物培养方法受低产量的限制10,11, 并被概念化为一个二进制状态, 其中细菌是存在于肠道或没有。高通量 DNA 测序已经彻底改变了微生物生态学, 使微生物群落的所有成员都能被捕获。但是, 序列读取长度和质量仍然是精确分类分配12的重要障碍。此外, 基于高通量的实验可能会受到批处理效果的影响, 其中测量受非生物或非科学变量13的影响。近年来, 已经建立了一些研究人类微生物菌群的方案, 包括美国肠道项目、美国 (美国) 人类微生物学项目和联合王国 (英国) MetaHIT 项目。这些倡议产生了大量的数据, 由于它们的做法缺乏连贯性, 因此难以比较。各种国际项目, 如国际人类微生物菌群联合会、国际人类微生物菌群标准项目以及国家标准和技术研究所 (NIST) 试图解决其中的一些问题14, 并制定了微生物菌群测量标准, 这将有助于实现可靠的生殖结果。这里描述的是一些广泛使用的方法15、16的集成协议, 用于 16S rRNA 高通量测序 (十六年代序列), 从粪便样本收集通过数据分析开始。该议定书描述了一种无柱 PCR 方法, 最初设计用于直接提取植物 DNA16, 以便在相对较短的时间内同时处理大量的粪便样本, 并以高质量的扩增 DNA 为靶向常规测序平台上微生物变量 V4 区域的测序。该议定书的目的是引导感兴趣的科学家开始使用 16S rRNA 扩增子测序, 在健壮, 生殖, 易于使用, 详细的方式, 使用重要的控制。具有引导和详细的分步协议可能会减少批处理效果, 从而允许在实验室之间进行更可比较的排序结果。
Subscription Required. Please recommend JoVE to your librarian.
Protocol
该研究的伦理批准是由示巴地方研究伦理学委员会批准的, 所有方法都是按照相关的指导方针和条例进行的。该议定书收到了当地道德审查委员会的病人同意例外, 因为使用的粪便材料已经作为临床检查的一部分提交给微生物学核心, 而没有其他人的年龄以外的可识别的病人信息,性别和微生物结果。从健康志愿者那里获得书面、知情的同意, 机构审查委员会批准了这项研究。其中一些示例已包含在以前的分析1中。
1. 样品处理
- 收集一个大约5毫米2涂片 (大致大小的铅笔橡皮擦) 从一个新鲜的粪便样本与无菌拭子在试管 (见材料表)。将含有粪便样品的所有拭子贮存在24小时内-80 摄氏度。粪便拭子可以留在那里, 直到进一步处理。
2. DNA 提取
- 在室温下解冻萃取和稀释液 (请参阅材料表)。
- 将粪便拭子转移到空的2毫升收集管 (请参阅材料表)。通过使用干净的剪刀来调整拭子棒的尺寸, 以使管道闭合时具有最小的交叉污染。加入 250 ul 的萃取液, 每个收集管含有粪便拭子和漩涡混合。
- 在沸腾的水浴中加热样品10分钟 (95–100°c).添加 250 ul 稀释溶液的每个样品和涡流混合。
- 储存2毫升的管含有提取的 DNA 和拭子在4摄氏度。
3. PCR 和图书馆准备
对于步骤3.1 和 3.2, 在提供干净、模板和扩增子自由环境的 PCR 工作站中工作。
- 根据其条码标记底漆 (表 1)。将双蒸馏水 (DDW) 中的每一底漆稀释至50微米浓度, 并贮存在-20 摄氏度。
- 用96井板进行 PCR 反应。每个盘子可以包含32不同的样品, 由32个不同的索引引物标记。在室温下解冻前底漆和32反向底漆, 稀释至5微米。
- 准备 pcr 反应混合100反应 (每个井的最终体积将是 20 ul) 混合 100 ul 的5微米前底漆, 1 毫升 2X PCR 主混合, 和 400 ul DDW。把这个 PCR 混合在每个井 (总共96口井) 15 ul。将每5微米反向索引底漆的 1 ul 添加到3个不同的井 (triplicates 中32种不同的引物, 总共提供96个井)。
- 在预 PCR 专用区, 这意味着一个干净的长凳, 是模板和扩增子免费, 添加 4 ul 的每一个提取的 dna 样本到反应混合物 (每个提取的 dna 样本被放大32个样本每 96-井板)。
- 运行 PCR 与以下设置: 初始变性94°c 3 分钟;其次35周期变性在94°c 为1分钟, 退火在55°c 为1分钟和延长在72°c 为1分钟;和最后的延伸在72°c 为10分钟。
- 在一个不同的长凳上, 这将被定义为后 pcr 专用区, 结合每三个 pcr 反应成一个单一的容积管 (60 ul 每样品)。
- 为了评估 pcr amplicons 的质量, 在溴溴染色的1% 琼脂糖凝胶上, 对每个联合 pcr 反应运行 4 ul。在紫外波长为 260 nm 的情况下, 正 amplicons 将出现在预期的 375–425 bp 波段大小。只有这些 amplicons 将包含在后续步骤中。
注: 每样样品都有三份放大, 这意味着每个样品在3种不同的 PCR 反应中被放大。不要放大。
4. 图书馆量化和清洁
- 为了获得所有 PCR 样品的摩尔浓度池, 用双链 DNA (dsDNA 定量试剂) 荧光灯核酸染色 (参见材料表) 量化每个扩增子, 适合于同时量化大型样本量。
注: 由于扩增子的尺寸范围是模棱两可的, 毫微克的实际用量用于加载摩尔浓度。 - 将每个样品的 500 ng 组合成一个单一的无菌管。漩涡混合。
- 在一个溴溴染色的1% 琼脂糖凝胶上运行 200 ul 的汇集图书馆。根据制造商的说明 (请参阅材料表) 并洗脱 80 ul DDW 中的汇集库, 从凝胶中提取 375–425 bp 带。
注: 池库应选择大小, 以减少来自宿主 DNA 的非特定放大产品。 - 根据制造商的说明, 使用高度敏感的 dsDNA 检测套件测量最终库浓度 (请参阅材料表)。根据制造商的说明, 使用对 DNA 库的高度敏感的分离和分析工具包测量准确的库大小 (请参阅材料表)。
5. 测序
- 根据顺序机协议, 将池库稀释到下午7点, 并添加20% 控制库 (参见材料表)。
- 对于排序, 使用自定义的 V4 放大引物免费的读引物 (请参见表 1)。
- 根据制造商的规范, 使用第三个自定义的测序底漆, 在附加周期中读取条形码 (请参见表 1), 并在每个方向上生成175个基的配对端读取。
- 根据制造商的协议, 运行顺序机并获取 FASTQ 文件。
6. 数据处理
- 在 QIIME 2 版 2017.7. 017中实现的数据精选管线中, 缝合并处理重叠的配对端 FASTQ 文件。根据示例特定的条形码解读取。
- 使用 DADA218进行质量控制和序列变体 (SVs) 检测。将读取从 3 ' 端和15个基点从 5 ' 端截断13个基元, 丢弃超过2个预期错误的读取。使用协商一致方法识别和删除合体-合体分别在样本中检测到, 并且在样本的足够部分中发现嵌合体的序列被删除。
- 使用朴素的贝叶斯拟合分类器执行 SVs 分类分类, 在 2013年8月99% 标识 Greengenes 数据库6上进行培训, 用于175长时间读取和正/反向底漆集。
- Rarefy 所有样本的深度为2146个序列。
- 在稀薄的样品上使用加权 UniFrac 测量β多样性 (样品多样性19,20), 以避免样本大小的影响。
- 使用所得到的距离矩阵来执行主坐标分析 (PCoA)。
注意: 数据处理脚本和映射文件作为补充材料提供 (辅助材料 1和2)。
Subscription Required. Please recommend JoVE to your librarian.
Representative Results
该协议的示意图说明显示在图 1中。
我们有前瞻性收集的粪便样本, 从住院病人疑似传染性腹泻。这些样本是在2月至2015年5月期间提交给示巴医疗中心的临床微生物学实验室的, 如前所述的1。粪便标本在临床微生物学实验室进行常规微生物培养, 并在十六年代平行的范围内进行广泛的高通量。此外, 健康成人的粪便样本进行了比较。
在此分析中, 重点是数据的质量控制, 并显示从 QIIME217获得的输出的代表性结果。首先, 回顾了6例阴性对照样品的测序结果, 包括: i. 无模板 pcr, ii.) 萃取稀释液中清洁无菌拭子的 pcr 方法及 iii.) 混合萃取和稀释的 pcr 技术解.表 2显示所有负控制样本都有≤118总读数 (平均29个序列), 主要来自支原体分类, 而所有其他样本分析显示平均总读数为 10622.57 (±1211.34 标准错误), 没有支原体读取通过了主要样品的质量控制。
十六个样本, 每起源于相同的粪便样本, 但通过不同的图书馆准备与不同的反向条形码底漆被用作积极的控制。这16个样本包括8个, 其中有阳性的区域性为弯曲杆菌、沙门氏菌或志贺菌, 临床微生物学实验室报告, 8 报告为有负文化。门级别的区域图显示在图 2中。来自同一粪便样本的样本, 但通过不同的图书馆准备不同的反向条形码引物显示非常一致的相对丰度 (在一级分类学相对丰度重复度之间的壁炉试验数值模拟 p 值 = 1 x 10-04基于9999复制)。
主坐标分析 (PCoA) 降低了输入数据的维数, 用于 UniFrac 矩阵直观地探索样本分离和相似性。在图 3A中, 来自同一原始粪便样本但经过不同库准备的序列样本颜色相同。的确, 那些位于相对地严密在 PCoA 剧情 (壁炉测试在重复 PC1 PC2 价值之间模拟 p 价值 = 1e-04 基于9999复制)。此外, 本研究中样本之间的平均加权 unifrac β多样性距离为 0.706, 而两个重复项之间的平均距离为 0.235 (魏氏秩和测试 p 值 = 4.205 x 10-11)。有趣的是, 对于弯曲杆菌、沙门氏菌或志贺氏菌enteropathogens 的区域性阳性样本, 与样品相比, PC1 值显著不同 (魏氏秩和测试 p 值 = 0.011)。消极文化结果。图 3B显示了相同的 PCoA 图, 但现在样本被 Proteobacteria 的相对丰度所着色。高 Proteobacteria 丰度的样品有显著不同的 PC1 值 (长矛相关 r =-0.533, p 值 = 0.000975)。这些结果与以前对这一数据的分析相一致, 在这里, 住院病人从 Proteobacteria1的分类中有了深刻的增加。在这里我们表明, PC1 与 Proteobacteria 丰度相关, 以文化正样本聚类为主的 PC1 轴的一侧和文化阴性样本聚类主要在另一边, 连同健康的样本。
图 1: 从粪便样本到相对微生物丰度的流程图.1)示例处理: 左侧, 粪便样本在无菌拭子管中收集。正确的, 拭子棒尺寸通过切割来调整, 以使2ml 管的闭合和存储。2) dna 提取: 通过直接 PCR 无柱法提取 dna。3)库准备: 左侧, 提取的 DNA 的 4 ul 用正向/反向索引引物放大。每个反向底漆包含12基本条码。每个样品都有三份放大。一个96井板包含32不同的条形码。对, 将每一个条码三个组合成一个体积管。一个96井板包含96不同的条形码。4)库量化: 上板、阳性样品筛选、溴溴化染色1% 琼脂糖凝胶 amplicons 检测。下面板, 量化的积极 amplicons, 以达到摩尔浓度 500 ng 扩增子从每个样本到一个单一的图书馆池。5)库净化: 上部面板, 汇集库是凝胶提取的大小选择和纯化。下部面板, 测量库的精确大小和浓度。6)排序: 池库的高通量排序。7)对 16S rRNA-扩增子微生物生态学的生物信息学管线进行了测序数据分析。请单击此处查看此图的较大版本.
图 2: 来自同一粪便样本的样本之间的一致性相对丰度, 但经过了不同的库准备.每样的分类相对丰度均按样品显示。相对丰度显示的16粪便样本, 从住院病人的疑似感染性腹泻, 每个通过不同的图书馆准备 (重复1和 2) 与不同的反向条形码引物和3健康控制。请单击此处查看此图的较大版本.
图 3: 系统发育多样性显示了来自同一粪便样本的样本的一致结果, 但经过了不同的库准备.在 UniFrac 矩阵上, 利用主坐标分析 (PCoA) 降低输入数据的维数, 直观地探索样本分离和相似性。两个样本之间的距离代表了它们的微生物组成的差异。(A) 加权 UniFrac PCoA 图。每个点表示一个单序列采样。来自同一原始粪便样本的样本颜色相同。临床微生物学报告的住院病人阳性粪便培养样本1是用正方形标记的, 住院的患者与文化阴性在三角1, 并且健康成人从单独排序运行1标记为填充圆圈。(B) 与 A 中相同的 PCoA, 但根据所示, 现在样品的 Proteobacteria 相对丰度是有色的。请单击此处查看此图的较大版本.
| 底漆名称 | 底漆序列 | 底漆说明 |
| 前进底漆 | AATGATACGGCGACCACCGAGATCTACAC-TATGGTAATT-GT-GTGCCAGCMGCCGCGGTAA | 5 ' 适配器序列向前底漆垫向前底漆链接器向前底漆 |
| 反向索引底漆 | CAAGCAGAAGACGGCATACGAGAT-XXXXXXXXXXXX-AGTCAGTCAG-CC-GACTACHVGGGTWTCTAAT | 3 ' 适配器序列的反向补充-Golay 条码-反向底漆垫-反向底漆链接器-反向底漆 |
| 阅读1测序底漆 | TATGGTAATT-GT-GTGCCAGCMGCCGCGGTAA | 向前底漆垫向前底漆链接器向前底漆 |
| 阅读2测序底漆 | AGTCAGTCAG-CC-GGACTACHVGGGTWTCTAAT | 反向底漆垫-反向底漆链接器-反向底漆 |
| 索引序列底漆 | ATTAGAWACCCBDGTAGTCC-GG-CTGACTGACT | 反向底漆反向补体反向引物链接器反向补片反向底漆垫 |
表 1: 底漆清单。
| 示例 | 读帖子数 初始质量控制在 QIIME |
| 萃取稀释液中清洁拭子的 PCR 研究-1 | 118 |
| 萃取稀释液中清洁拭子的 PCR 研究-2 | 39 |
| 无模板的 PCR-1 | 0 |
| 无模板的 PCR-2 | 0 |
| PCR 在萃取和稀释溶液混合物-1 | 16 |
| PCR 在萃取和稀释溶液混合物-2 | 0 |
| 样品 (平均标准误差) | 10622.57 @ 1211.34 |
表 2: 阴性对照和粪便样本读取次数。
辅助材料 1: 数据处理脚本.请单击此处下载此文件.
辅助材料 2: 映射文件.请单击此处下载此文件.
补充材料 3: 底漆方案.请单击此处下载此文件.
Subscription Required. Please recommend JoVE to your librarian.
Discussion
16S rRNA-扩增子和 metagenomics 猎枪测序在临床微生物学应用中得到了普及21, 2223.这些技术有利于提高捕获可培养和非可培养类群的能力, 提供有关致病接种的相对丰度的数据, 以及它们识别更精确的 polymicrobial 传染性的能力。指纹24。微生物研究领域的进展产生了大量的数据, 由于各种方法缺乏一致性, 因此难以比较。最近几年, 一些财团试图解决其中的一些问题。该议定书遵循类似的图书馆准备工作, 作为地球微生物群项目 (EMP)。然而, 添加的是一个详细的分步协议的 DNA 提取从粪便样品通过分析管道。
此前, 16S rRNA 测序已被用于描述健康非住院患者粪便样本的微生物组成模式, 以及疑似感染性腹泻的住院病人, 以及传统的文化结果。分析结果表明, 住院患者在 Proteobacteria 的分类中有很大的增加1。这里描述的是一个综合协议的广泛使用的方法 16S rRNA 基因扩增子排序使用直接 PCR 无柱方法, 最初设计用于提取植物DNA16。这种方法可以有效和快速地提取高质量的扩增 DNA 库, 针对 16S rRNA 基因的 V4 变量区域, 从大量适合高通量 16S rRNA 测序的样本。
由于这是一种基于 DNA 的测序方法, 因此获得了单个粪便样本中需氧生物和厌氧微生物群落的数据。该协议采用了最初设计为15的引物, 针对 16S rRNA 基因的 V4 区域。重要的是, 反向放大底漆包含十二个基本条码序列, 它支持在每个车道15中共用多达2167种不同的样本。正向放大底漆在适配器区域中包括九个额外的底座, 用于支持序列平台上的配对端排序 ( 表 1和辅助材料 3)。此方法启用了处理一个库, 由250个粪便样本组成, 在一条车道上排序, 其运行深度为 10622.57, 1211.34 (平均标准错误) 读取每个样本的成本为每样本30美元。
为确保质量控制, 建议使用以下负面控制;i) pcr 无 dna 模板, ii) pcr 的 dna 提取从清洁消毒拭子, 和 iii) pcr 的混合萃取和稀释溶液。至于积极的控制, 我们建议包括来自同一粪便样本的样本, 但通过不同的图书馆准备与不同的反向条形码底漆和样品从以往的运行作为内部质量控制。值得注意的是, 其他可选的积极控制是由 NIST 提供的微生物测量标准。相对于实际粪便样本 (表 1), 阴性对照组的平均读数覆盖率明显低于, 主要由支原体分类。我们进一步显示了使用样本所获得的测序结果的一致性, 它们分别来源于相同的粪便材料, 但经过了不同的库准备 (图 3)。这一结果还证实, 在使用我们的集成方法时, 样品之间没有交叉污染。
在本协议中, 我们引入了一个统一的、可行的协议来分析大量的样本。该议定书的目的是指导感兴趣的科学家开始使用 16S rRNA-扩增子测序的健壮, 生殖, 易于使用, 廉价, 详细的方式从粪便收集到数据分析, 使用重要的控制。使用此标准化的详细指导协议, 可以最小化批处理效果, 并允许在不同的实验室之间进行更可比较的排序结果。
Subscription Required. Please recommend JoVE to your librarian.
Disclosures
作者没有什么可透露的。
Acknowledgments
这项工作的一部分是由核心项目 (赠款 41/11)、以色列科学基金会 (908/15 号赠款) 和欧洲克罗恩病组织 (ECCO) 支持。
Materials
| Name | Company | Catalog Number | Comments |
| Primers | Integrated DNA Technologies (IDT) | ||
| Extraction solution | Sigma-Aldrich | E7526 | |
| Dilution solution | Sigma-Aldrich | D5688 | |
| Kapa HiFi HotStart ReadyMix PCR Kit | KAPABIOSYSTEMS | KK2601 | PCR Master mix |
| Quant-iT PicoGreen dsDNA Reagent kit | Invitrogen | P7589 | dsDNA quantify reagent |
| MinElute Gel extraction kit | Qiagen | 28606 | |
| Agarose | Amresco | 0710-250G | |
| Ultra Pure Water Dnase and Rnase Free | Biological Industries | 01-866-1A | |
| Qubit dsDNA HS assay kit | Molecular probes | Q32854 | dsDNA detecting kit |
| High Sensitivity D1000 | Agilent Technologies | Screen Tape 5067-5582 | separation and analysis |
| Screen Tape Assay | Agilent Technologies | Reagents 5067-5583 | for DNA libraries |
| PhiX Control v3 | Illumina | 15017666 | control library |
| MiSeq Reagent Kit v2 (500 cycle) | Illumina | MS-102-2003 | |
| Ethidium Bromide | Amresco | E406-10mL-TAM | |
| 2 mL collection tubes | SARSTEDT | 72.695.400 | Safe Seal collection tubes |
| Plastic stick swab in PP test tube | STERILE INTERIOR | 23117 | |
| Name | Company | Catalog Number | Comments |
| Equipment | |||
| PCR Machine | Applied Biosystems | 2720 Thermal Cycler | |
| Sequncing Machine | Illumina | Miseq | |
| PCR workstation | Biosan | UV-cleaner | |
| scissors | |||
| vortexer | Scientific Industries | Vortex-Genie 2 |
References
- Braun, T., et al. Fecal microbial characterization of hospitalized patients with suspected infectious diarrhea shows significant dysbiosis. Scientific Reports. 7, 1088 (2017).
- De Filippo, C., et al. Impact of diet in shaping gut microbiota revealed by a comparative study in children from Europe and rural Africa. Proceedings of the National Academy of Sciences of the United States of America. 107, 14691-14696 (2010).
- Yatsunenko, T., et al. Human gut microbiome viewed across age and geography. Nature. 486, 222-227 (2012).
- Rodriguez, J. M., et al. The composition of the gut microbiota throughout life, with an emphasis on early life. Microbial Ecology in Health and Disease. 26, 26050 (2015).
- Turnbaugh, P. J., et al. An obesity-associated gut microbiome with increased capacity for energy harvest. Nature. 444, 1027-1031 (2006).
- McDonald, D., et al. An improved Greengenes taxonomy with explicit ranks for ecological and evolutionary analyses of bacteria and archaea. The ISME Journal. 6, 610-618 (2012).
- Bakken, J. S., et al. Treating Clostridium difficile infection with fecal microbiota transplantation. Clinical Gastroenterology and Hepatology. 9, 1044-1049 (2011).
- Khoruts, A., Dicksved, J., Jansson, J. K., Sadowsky, M. J. Changes in the composition of the human fecal microbiome after bacteriotherapy for recurrent Clostridium difficile-associated diarrhea. Journal of Clinical Gastroenterology. 44, 354-360 (2010).
- Nood, E. V., et al. Duodenal infusion of donor feces for recurrent Clostridium difficile. The New England Journal of Medicine. 368, 407-415 (2013).
- Koplan, J. P., Benfari Ferraro, M. J., Fineberg, H. V., Rosenberg, M. L. Value of stool cultures. The Lancet. 316, 413-416 (1980).
- Platts-Mills, J. A., Liu, J., Houpt, E. R. New concepts in diagnostics for infectious diarrhea. Mucosal Immunology. 6, 876-885 (2013).
- Bokulich, N. A., et al. Quality-filtering vastly improves diversity estimates from Illumina amplicon sequencing. Nature Methods. 10, 57-59 (2013).
- Leek, J. T., et al. Tackling the widespread and critical impact of batch effects in high-throughput data. Nature Reviews Genetics. 11, (2010).
- Dubilier, N., McFall-Ngai, M., Zhao, L. Create a global microbiome effort. Nature. 526, 631-634 (2015).
- Caporaso, J. G., et al. Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. PNAS. 108, 4516-4522 (2011).
- Flores, G. E., Henley, J. B., Fierer, N. A direct PCR approach to accelerate analyses of human-associated microbial communities. PloS One. 7, (2012).
- Caporaso, J. G., et al. QIIME allows analysis of high-throughput community sequencing data. Nature Methods. 7, 335-336 (2010).
- Callahan, B. J., et al. DADA2: High-resolution sample inference from Illumina amplicon data. Nature Methods. 13, 581-583 (2016).
- Lozupone, C., Knight, R. UniFrac: A new phylogenetic method for comparing microbial communities. Applied and Environmental Microbiology. 71, 8228-8235 (2005).
- Lozupone, C., Lladser, M. E., Knights, D., Stombaugh, J., Knight, R. UniFrac: An effective distance metric for microbial community comparison. The ISME Journal. 5, 169-172 (2011).
- Srinivasan, R., et al. Use of 16S rRNA gene for identification of a broad range of clinically relevant bacterial pathogens. PloS One. 10, e0117617 (2015).
- Hiergeist, A., Gläsner, J., Reischl, U., Gessner, A. Analyses of intestinal microbiota: culture versus sequencing. ILAR Journal. 56, 228-240 (2015).
- Didelot, X., Bowden, R., Wilson, D. J., Peto, T. E. A., Crook, D. W. Transforming clinical microbiology with bacterial genome sequencing. Nature Reviews Genetics. 13, 601-612 (2012).
- Salipante, S. J., et al. Rapid 16S rRNA next-generation sequencing of polymicrobial clinical samples for diagnosis of complex bacterial infections. PloS One. 8, e65226 (2013).
- Caporaso, J. G., et al. Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms. The ISME Journal. 6, 1621-1624 (2012).
