一种基于高通量基因测序建立粪便菌群基因文库的方法
技术领域
本发明属于生物技术领域中细菌的分子生物学检测方法技术领域，特别是涉及一种基于高通量基因测序建立粪便菌群基因文库的方法。
背景技术
近年来肠道微生物与饮食和健康的关系被越来越多的阐释，在很大程度上得益于肠道微生物研究技术的快速发展。肠道微生物的研究大体经历了培养依赖的方法，非培养依赖的传统分子生物学方法，基于测序的高通量组学方法3个阶段。
通过比较3个阶段研究方法的特点与应用，人们可以发现培养依赖的方法在鉴定菌种的同时即可获得相应菌株，方便后续研究，但培养法本身费时费力，肠道微生物以厌氧菌和兼性厌氧菌为主，培养起来更加困难，并且在培养的过程中菌种比例会发生改变，使得其应用存在瓶颈。非培养依赖的传统分子生物学方法可以不经培养直接从样品中提取肠道微生物基因组，利用分子生物学手段进行分离、鉴定和定量，使其结果可以比较准确的反应肠道微生物中高丰度菌种的组成和真实比例。
期刊Journal of Clinical Microbiology的第1122–1124页中，Kevin Jeng公开了沙门氏菌的快速检测，题目为Application of a 16S rRNA PCR–High-Resolution Melt Analysis Assay for Rapid Detection of Salmonella Bacteremia。该发明通过针对沙门氏菌反应性血清型肠炎菌血症患者进行分析，以确定其感染情况。
但是传统分子生物学方法存在着通量低的缺陷，靶向方法如实时定量PCR每次只能研究一种或一类肠道微生物，而非靶向方法如最常用的梯度变形凝胶电泳(Denaturing gradient gel electrophoresis，DGGE)受限于灵敏度，往往只能研究肠道中高丰度的微生物。肠道微生物群落组成复杂，每种细菌都与其他细菌形成复杂的相互关系网络，低丰度的菌种同样扮演者重要的角色，所以传统分子生物学研究结果往往具有片面性。
16SrDNA因其序列在物种间的高度多样性，成为细菌分类学研究的“分子钟”，具9个可变区和9个保守区间隔排列的特征，可变区因细菌而异，且变异程度与细菌的系统发育密切相关。因此通过分析可变区的序列即可得到各细菌的分类学特征。
传统分子生物学方法中有时也应用Sanger法对单一菌种的16SrDNA进行测序和鉴定，但受限于通量低，效率低，不能很好地满足科研、临床的需要。
发明内容
本发明的目的是针对以上要解决的技术问题，提供一种获得覆盖整个微生物群落状况的方法，以为实现快速、便捷及准确检测人肠道菌种组成与研究菌落生态结构建立基础性的菌落数据库。
具体的，为了解决上述问题，本发明提供了一种肠道微生物全部菌种的结构和功能的靶标基因序列库的建立方法，尤其具体的，提供了一种基于高通量基因测序建立粪便菌群基因文库的方法，该方法采用“巢式PCR”的方法对16S rDNA进行富集扩增，最大限度减少了宿主和食物残渣基因组污染，同时将V3和V6组合，进行特异扩增并进行大规模并行测序，获得菌群的靶标基因序列库。利用该序列库，能够将菌群鉴定从属(Genus)级别，提升至种(Species)级别。
具体的，根据上述方法，其中“巢式PCR”的方法包括从总基因组DNA中进行16S rDNA全长富集扩增，扩增所用的引物组由3对PCR引物组成，该3对引物序列如下：
第1对PCR引物序列为：
SEQ ID NO.1：5’-AGAGTTTGATCCTGGCTCA-3’(16S-27F)
SEQ ID NO.2：5’-GGTTACCTTGTTACGACT-3’(16S-1492R)
第2对PCR引物序列为：
SEQ ID NO.3：5’-CTCCTACGGGAGGCAGCA-3’(V3-357F)
SEQ ID NO.4：5’-ATTACCGCGGCTGCTGG-3’(v3-534R)
第3对PCR引物序列为：
SEQ ID NO.5：5’-CAACGCGAAGAACCTTAC-3’(v6-967F)
SEQ ID NO.6：5’-GGGTTGCGCTCGTTG-3’(v6-1100R)
上述引物组中，第1对PCR引物序列扩增细菌16S rDNA的位点27/1492区域，扩增产物长度为1465bps，约15kb。
上述引物组中，第2对PCR引物序列扩增细菌16S rDNA V3的位点357/534区域，扩增产物长度约为177bps；
上述引物组中，第3对PCR引物序列扩增细菌16S rDNA V6的位点967/1100区域，扩增产物长度约为133bps；
为此，本发明还提供了一种基于高通量基因测序建立粪便菌群基因文库的方法，该方法包括如下步骤：
步骤一、提取粪便菌群总基因组DNA；
步骤二、用上述第1对引物对总基因组DNA进行16S rDNA全长富集扩增，并用16S rDNA的V3区和V6区第2对引物和第3对引物分别进行特异扩增；
步骤三、进行大规模并行测序，建立文库。
上述方法中，步骤二的特异扩增的产物分别为约177bps和133bps。
上述方法中，步骤三的大规模并行测序使用Ion Torrent个人化操作基因组测序仪进行大规模并行测序。
上述方法中，步骤一提取粪便菌群总基因组DNA按照按照(TIANGEN公司DP328)试剂盒操作说明书所述的条件提取粪便菌群总基因组DNA进行，优选的，根据肠道菌进行条件优化，优化后的条件如下：
将样品预处理的裂解温度调整为95℃孵育10min，能使较难破壁的革兰氏阳性细菌裂解破壁，使得抽提产物能覆盖到全部微生物，同时使DNA结构完整性保持最佳状态；样品裂解后，在原有一块的基础上，再添加一块抑制剂吸附片InhibitEX，能针对于粪便样品杂质去除PCR抑制作用。
优选的，上述步骤二中，对总基因组DNA进行16S rDNA全长富集扩增的方法为“巢式PCR三步循环法”。
具体的，巢式PCR三步循环法的操作步骤为：
首先将细菌总基因组DNA稀释到浓度为1ng/ul，取1μl为模板作为第一次扩增的模板，以第1对引物进行PCR扩增，得到16S rDNA的PCR产物；第二次扩增，以1ul其扩增富集的16S rDNA的PCR产物稀释50倍为模板，引物为第2对引物和第3对引物，用相同配液体系和反应条件进行二次扩增。
其中巢式PCR三步循环法的反应体系为：PCR为50μl总体系，具体是2×Taq PCR StarMix with Loading Dye 25μl，Primer-F/Primer-R 2μl，浓度稀释到1ng/μl以下的细菌gDNA模板2ul，用灭菌水补足50μl体系。
其中巢式PCR三步循环法的反应条件为：95℃预变性3min，然后95℃变性30sec，58℃退火30sec，72℃延伸30sec，共循环28次，最后72℃再延伸5min。
上述方法中，16S rDNA和V3、V6区PCR产物鉴定方法为：取5μl第一次扩增完后的16S rDNA产物和第二次扩增的V3、V6区PCR产物直接跑1％琼脂糖凝胶电泳。
V3、V6区PCR产物回收方法为：利用Backman品牌AgencourtAMPure磁珠按样品与磁珠体积1.5倍体积纯化回收目的片段。
上述方法中，优选的，步骤三的操作步骤如下：
(1)末端补平DNA、纯化回收目的片段：利用Backman品牌AgencourtAMPure磁珠按样品与磁珠体积1.5倍体积(75μL)纯化回收目的片段；
(2)连接接头，缺口平移、纯化回收目的片段：利用Backman品牌AgencourtAMPure磁珠按样品与磁珠体积1.2倍体积(120μL)纯化回收目的片段；
(3)文库扩增、纯化回收目的片段：利用Backman品牌AgencourtAMPure磁珠按样品与磁珠体积1.2倍体积(120μL)纯化回收目的片段；
(4)文库质控与决定模板稀释因子：
检测浓度：取2μl扩增后的文库，将文库稀释为终浓度100pM的稀释因子，通过Qubit 2.0HS DNA荧光检测试剂进行浓度检测；
检测片段分布：根据Qubit检测浓度将待测样品稀释至1ng/μl，通过High Sensitivity DNA Assay试剂盒芯片检测；
分别将V3、V6等摩尔浓度混合。取100ng作为建库起始量。
(5)模板制备：采用Ion OneTouch 200Template Kit OT2试剂盒和OneTouch仪器进行乳液PCR；采用MyOne beads磁珠和Ion OneTouch ES仪器进行ISPs富集，最终得到3’端带有磁珠的阳性模板文库；
(6)测序及测序结果分析：将制备好的文库在Life technologies公司的Ion PGM sequencer高通量测序仪上进行测序，样本生成测序读长在200bp，深度要求在100X以上，数据量在0.05M reads以上；经FastQC-3.4.1.1分析测序质控，使用QIIME软件，通过与数据库比对，即可获得菌群生态的物种组成、比例。通过获得菌群生态的物种组成、比例判断优势劣势菌情况。
上述方法中，V3、V6区纯化回收产物质控与决定建库起始量为：
(1)检测浓度：取2μl V3、V6区纯化回收产物，通过Qubit 2.0HS DNA荧光检测试剂进行浓度检测；
(2)检测片段分布：根据Qubit检测浓度将待测样品稀释至1ng/μl，通过High Sensitivity DNA Assay试剂盒芯片检测。
本发明还提供了一种从总基因组DNA中进行16S rDNA全长富集扩增所用的引物组，该引物组由3对PCR引物组成，该3对引物序列如下：
第1对PCR引物序列为：
SEQ ID NO.1：5’-AGAGTTTGATCCTGGCTCA-3’
SEQ ID NO.2：5’-GGTTACCTTGTTACGACT-3’
第2对PCR引物序列为：
SEQ ID NO.3：5’-CTCCTACGGGAGGCAGCA-3’
SEQ ID NO.4：5’-ATTACCGCGGCTGCTGG-3’
第3对PCR引物序列为：
SEQ ID NO.5：5’-CAACGCGAAGAACCTTAC-3’
SEQ ID NO.6：5’-GGGTTGCGCTCGTTG-3’。
本发明与现有技术相比具有以下优点：
(1)建库速度快，通量大：检测方法搭配了目前最快的高通量测序平台life品牌的Ion Torrent PGM测序仪，上百个样品并行检测，上机测序时间仅为3个小时，文库满足科研、临床需要；
(2)文库高特异性，高灵敏度：V3和V6组合，菌群从鉴定到“属(Genus)”级别，提升至鉴定到“种(Species)”的级别。同时，选择“巢式PCR”的方法对以16S rDNA进行富集扩增，最大限度减少了宿主和食物残渣基因组污染；
(3)文库建立采用原位活检：以新鲜或冻存的样品原位活检粪便的菌群组成，不需要经过筛选培养，从而保持原来菌群生态组成。
(4)本发明以高通量组学为手段，能够快速、方便、同步、高效、高特异、高灵敏地检测到原来菌群的生态组成，获得覆盖整个微生物群落的文库信息，真正实现了对肠道微生物全部菌种的结构和功能的整体研究。
附图说明
图1为实施例粪便菌群总基因组抽提的检测结果。
图2为实施例粪便菌群16SrDNA扩增的检测结果。
图3为实施例粪便菌群V3、V6扩增的检测结果。
图4为实施例粪便菌群V3扩增纯化回收的2100检测结果。
图5为实施例粪便菌群V6扩增纯化回收的2100检测结果。
图6为实施例粪便菌群V3、V6文库构建后的2100检测结果。
图7为实施例粪便菌群V3、V6文库测序上机后的数据量结果。
图8为实施例粪便菌群V3、V6文库测序上机后片段分布结果。
图9为实施例粪便菌群的微生物百分比热图和柱状图。
图10为实施例粪便菌群的微生物丰富度柱状图。
附图1-8中外文的含义如下：
M指maker 1kb，DNA分子量标准；
bp，碱基对个数
16s rDNA，16s核糖体DNA；
kb，一千个碱基对个数
V3，核糖体16s DNA V3区；
V6，核糖体16s DNA V6区；
FU，电信号值；
Total Bases，总数据量；
M，Mb，兆字节；
Key Signal，关键信号值；
ISP Loading，Ion Sphere Particles Loading，离子球粒子覆盖度；
ISP Density，Ion Sphere Particles Density，离子球粒子密度；
Wells，纳米孔数；
GZ-170-2015.09.08-16s v3v6-200，测序流水号和测序样品命名；
Loading Density，Ion Sphere Particles Density average，离子球离子密度平均值；
(Avg-76％)，本次离子球离子密度测定的平均值是76％，正常范围是在65-85％之间，属于较佳值。
None bp Mean，无平均值；
median，中值；
Mode，众数；
Read Length，测序序列长度；
Read length Histogram测序序列长度直方图；
count，计数值；
附图9中外文的含义如下：
Color Key and Histogram，颜色信号与柱状图；
Value，值；
count，计数点；
AB1.1，携带艰难梭菌AB毒素阳性的粪便样品1
AB1.2，携带艰难梭菌AB毒素阳性的粪便样品2
NC，normal control，空白对照
附图10中外文的含义如下：
AB1-1，携带艰难梭菌AB毒素阳性的粪便样品1
AB1-2，携带艰难梭菌AB毒素阳性的粪便样品2
NC，normal control，空白对照
具体实施方式
下面通过结合附图和具体实施例对本发明的技术方案做进一步详细说明，但本发明并不限于以下实施例。实施例中所用的材料为广州军区总医院提供的一组确诊为携带艰难梭菌AB毒素阳性的患者粪便。
本发明的建立文库的实验流程具体步骤如下。
1.提取粪便菌群总基因组DNA：
严格按照TIANamp Stool DNA Kit(TIANGEN公司DP328)试剂盒操作说明书提取粪便菌群总基因组DNA，并针对肠道菌进行了条件优化：将样品预处理的裂解温度调整为95℃孵育10min，能使较难破壁的革兰氏阳性细菌裂解破壁，保证抽提产物能覆盖到全部微生物，同时使 DNA结构完整性保持最佳状态；样品裂解后，再添加一块抑制剂吸附片InhibitEX，能针对于粪便样品杂质去除PCR抑制作用。具体操作如下：
1.1提取粪便菌群总基因组DNA
(1)称取粪便样本50mg至1.5ml离心管中，并将管置于冰上。
(若是液态样本则转移200μl至离心管中)
(2)向样本中加入1.5ml缓冲液GSL，间歇振荡3min至样本混匀。
(3)95℃孵育10min。
(4)涡旋15sec。12,000rpm离心1min。转移上清液1ml至新的1.5ml离心管。
(5)添加一块抑制剂吸附片InhibitEX，振荡至吸附片彻底打开重悬。室温孵育5min，使吸附片能充分作用。
(6)12,000rpm离心3min。
(7)将上一步所得上清液转移至新的1.5ml离心管，重复步骤6。
(8)转移所得上清液400μl至新的1.5ml离心管，加入20μl ProteinaseK。
(9)加入400μl缓冲液GB，涡旋15sec。
(10)70℃孵育10min。
(11)加入400μl无水乙醇，涡旋混匀。
(12)将上一步所得溶液加入到一个吸附柱中12,000rpm离心30sec，倒掉废液，将吸附柱放入收集管中。
(13)向吸附柱中加入500μl缓冲液GD，12,000rpm离心30sec，倒掉废液，将吸附柱放入收集管中。
(14)向吸附柱CR2中加入600μl漂洗液PW，12,000rpm离心30sec，倒掉废液，吸附柱CR2放入收集管中。
(15)重复操作步骤14。
(16)将吸附柱放回收集管中，12,000rpm离心2min，倒掉废液。将吸附柱置于室温放置数分钟，以彻底晾干吸附材料中残余的漂洗液。
(18)将吸附柱转入一个干净的离心管中，向吸附膜的中间部位悬空滴加50μlpH值7.5的无核酸水，放置在65℃烘箱5min，12,000rpm离心2min，将溶液收集到离心管中。DNA产物应保存在-20℃。
1.2菌群总基因组浓度检测与质控
对抽提后的DNA片段用微量紫外分光光度计和1％琼脂糖凝胶电泳检测浓度与纯度。DNA应在OD260处有显著吸收峰，OD260/OD280比值应为1.8-2.0，浓度为500ng/μl。电泳胶图中在约15kb处条带单一清晰，说明DNA质量合格。细菌总基因组电泳图如图1所示：M：Marker；1-18:菌群总基因组DNA。
2.设计PCR引物富集16SrDNA和捕获V3、V6区
2.1富集16SrDNA和捕获V3、V6区
设计3对PCR引物对抽提出来的总基因组DNA进行16S rDNA全长富集扩增，保证扩增的产物能完全覆盖到高突变区域；然后取1μlPCR产物为模板，用16S rDNA的V3区和V6区特异引物分别进行特异捕获。
第一对PCR引物序列为：
16S-27F:5’-AGAGTTTGATCCTGGCTCA-3’
16S-1492R:5’-GGTTACCTTGTTACGACT-3’
扩增细菌16S rDNA的位点27/1492区域，扩增产物长度为1465bps，约15kb；
第二对PCR引物序列为：
V3-357F:5’-CTCCTACGGGAGGCAGCA-3’
v3-534R:5’-ATTACCGCGGCTGCTGG-3’
扩增细菌16S rDNA V3的位点357/534区域，扩增产物长度约为177bps；
第三对PCR引物序列为：
v6-967F:5’-CAACGCGAAGAACCTTAC-3’
v6-1100R:5’-GGGTTGCGCTCGTTG-3’
扩增细菌16S rDNA V6的位点967/1100区域，扩增产物长度约为133bps；
用上述引物配合配合如下反应体系：PCR为50μl总体系:
首先将细菌总基因组DNA稀释至浓度为1ng/μl，取1μl为模板反应条件为“巢式PCR三步循环法”
第二次扩增，以1ul其扩增富集的16SrDNA的PCR产物稀释50倍为模板，分别用V3和V6区两对引物，用相同配液体系和反应条件进行二次扩增。目的模板的浓度稀释至1ng/μl以下有利于提高PCR的特异性。
2.2 16SrDNA和V3、V6区PCR产物鉴定
取5μl第一次扩增完后的16SrDNA产物和第二次扩增的V3、V6区PCR产物直接跑1％琼脂糖凝胶电泳，如图2、图3所示。图2，扩增产物16S rDNA约为15kb；图3，V3、V6区扩增产物分别为约177bps和154bps。
2.3V3、V6区PCR产物磁珠回收
利用Backman品牌AgencourtAMPure磁珠按样品与磁珠体积1.5倍体积纯化回收目的片段。具体方法如下:
(1)将样品补水至50μl,加入1.5倍体积的磁珠(75μl)，用移液器吹打混匀，室温孵育5分钟；
(2)将样本管子放置在磁力架上，等待2-3分钟直到管内溶液完全澄清。
(3)小心吸取上清液,移除和弃掉上清液。
(4)加500μL新鲜配置的70％的酒精溶液。通过旋转管子清洗磁珠，旋转180°。总共旋转5次以彻底清洗。注：旋转可以利用磁性打散珠，去除底物。
(5)放置样本管子在磁力架上澄清后去除上清液,再旋转离心后去除残留的酒精溶液。
(6)在室温条件下完全风干磁珠5分钟，直到磁珠表面干燥，无液体光泽。
(7)加25μL LowTE到样本管子内并震荡10秒。用枪头上下吹吸溶液多次以保证洗脱完全混匀。
(8)将样本管子放置在磁力架上，待溶液澄清后,将洗脱的溶液(上清液)转移到一个新的.5ml低吸附管子内保存。
2.4V3、V6区纯化回收产物质控与决定建库起始量
(1)Qubit检测浓度
取2μl V3、V6区纯化回收产物，通过Qubit 2.0HS DNA荧光检测试剂进行浓度检测。检测结果：V3产物为53ng/μl、V6产物为25ng/μl。从浓度上看，说明此磁珠回收的产量相对于切胶回收法较高效，同时弥补了切胶法无法去除大片段污染的不足。
(2)安捷伦2100检测片段分布
根据Qubit检测浓度将待测样品稀释至1ng/μl，通过High Sensitivity DNA Assay试剂盒芯片检测。检测结果:V3区主峰为179bp和205bp大小；V6区主峰为154bp大小。从片段分布上看，V3、V6区条带都符合预期大小，对于V3区存在两个主峰，205bp的主峰可能是肠道特有的优势菌的扩增产物。片段大小在160-200之间，有利于构建200base读长的测序DNA文库。安捷伦2100仪器对文库的片段大小检测如图4、5。
3.大规模并行测序
文库构建采用life品牌的DNA文库构建试剂盒，按照200bp扩增子的标准操作流程进行建库，并针对肠道菌16SrDNA的片段特征进行优化。文库构建即为对靶目标序列前后两端加上与引物匹配的位点并为识别不同样品插入约十个碱基的特异性标签双链接头。此标签接头，为后面的扩增反应及测序提供匹配位点和识别信号。
3.1末端补平DNA
(1)在低吸附的1.5ml管内混合并混匀以下试剂成分，然后放置在室温条件下孵育30分钟。
(2)利用Backman品牌Agencourt AMPure磁珠按样品与磁珠体积1.5倍体积(75μL)纯化回收目的片段，经验证此浓度能确保片段选择中最大限度保留目标主峰200base，去除接头二聚体的污染。
3.2连接接头，缺口平移
(1)在连接反应体系中混合并混匀以下试剂成分。
(2)在PCR仪上按照以下程序孵育：
(3)利用Backman品牌Agencourt AMPure磁珠按样品与磁珠体积1.2倍体积(120μL)纯化回收目的片段，经验证此浓度能确保片段选择中最大限度保留加接头后的目标主峰。
3.3文库扩增
(1)在PCR管中混合并混匀以下试剂成分。
(2)设置以下程序进行PCR反应。
(3)利用Backman品牌Agencourt AMPure磁珠按样品与磁珠体积1.2倍体积(120μL)纯化回收目的片段，经验证此浓度能确保片段选择中最大限度保留目标主峰300，同时去除接头二聚体和引物污染。
3.4文库质控与决定模板稀释因子
(1)Qubit检测浓度
取2μl扩增后的文库，通过Qubit 2.0HS DNA荧光检测试剂进行浓度检测。检测结果：V3和V6片段混合文库为1.68ng/μl。从浓度上看，在1-10ng/μl范围，又不过度扩增而失真。(2)安捷伦2100检测片段分布
根据Qubit检测浓度将待测样品稀释至1ng/μl，通过High Sensitivity DNA Assay试剂盒芯片检测。检测结果:V3和V6片段混合文库为237bp、263bp和288bp大小；此为目标序列加上两端接头85bp，实际插入片段为152bp、178bp和203bp，此与建库前加进去的V3和V6区片段模板完全一致，且大小符合测序长度要求，由此说明本批样品文库构建质量合格。根据安捷伦2100提供的50-1000bp范围的摩尔浓度将文库稀释为终浓度100pM的稀释因子。16SV3V6文库安捷伦2100检测结果见图6。
3.5模板制备
文库构建完成后，通过在油包水乳液PCR中序列簇成百万倍扩增，并以此为模板进行上机测序。经验证取6μl的100pM文库稀释因子能降低主峰分散带来的多克隆现象，可获得更多的测序仪数据量。采用Ion OneTouch 200Template Kit OT2试剂盒和OneTouch仪器进行乳液PCR；采用MyOne beads磁珠和Ion OneTouch ES仪器进行ISPs富集，最终得到3’端带有磁珠的阳性模板文库。
3.6Ion Torrent PGM测序
对PGM测序仪初始化前，首先按Ion PGM Sequencing 200kit v2试剂盒清洗标准操作流程对测序仪管道先后用氯水和18MΩ去离子水清洗15min。然后，按照安装W1、W2、W3反应溶液进行初始化调节pH值，校正至pH为7.5-7.6。初始化后，用100％异丙醇和退火缓冲液Annealing buffer清洗314测序芯片，然后放到PGM测序仪上反应孔信号检测。同时，向富集好的文库中加入测序引物、文库内参Control Ion Spheres。然后置PCR仪上95℃，2min；37℃，2min退火变性，加入测序聚合酶，将文库样品加载到测序芯片上。最后，将芯片放到PGM测序仪上，设置500个flow测序循环，反应两个小时，测序完成。
3.7测序结果分析
样品通过上机测序，得到的最原始的测序数据以FASTQ格式的形式储存。
本次测序数据结果报告：314基因芯片实际测得有效区域达76％，最终产生的原始测序数据总量为45.1Mb；测序读长为150bp、179、205bp，均反映插入的目的片段实际长度。(见图7为实施例粪便菌群V3、V6文库测序上机后的数据量结果；图8为实施例粪便菌群V3、V6文库测序上机后片段分布结果)
本次16S rRNA测序得到的FASTQ结果，通过使用QIIME软件进行操作单元分类(没有进行过滤)，得到分别有自己代表序列的操作单元，使用QIIME里面的分类功能，根据Greengene database分类OTU的代表序列，得到其界、门、纲、目、科、属分别是那种微生物。分类分析结果得到29个能分类到属并大于0.001(0.1％)的微生物优势属，并且根据结果绘制了图9微生物百分比热图和柱状图和图10微生物丰富度柱状图。
根据分类分析学的分析结果，数据用更直观的方法去展示。
根据图9，微生物百分比热图，能更好观测到不同样品中的微生物差异情况：
根据样品聚类情况，知道携带艰难梭菌患者样品AB1.1和AB1.2的相似度比样品正常人NC(空白对照)的情况更接近，且可以明显观测微生物百分比含量在不同样品间的差别。
根据图10，微生物丰富度柱状图，能清晰看到在三组样品中33个微生物优势属占样品百分比的情况，更好地表达样品内优势菌属的占百分比的状况。如图10所示，Ruminococcus(约9％-25％)，Corynebacterium(约3％-9％)、Enhydrobacter(约2％-5％)，Clostridium(约2％-7％)，Enterococcus(约1％-2％)分别在3个样品中均为优势的5个菌属，其中本次实验最为关注的Clostridium(梭菌属)在三组样品中的百分比含量分别是NC样为2.65％，AB1-1样为5.61％，AB1-2样为7.04％。
由此可见，艰难梭菌属患者比正常人患者所含比例明显增高，其菌群结构的动态也有稍微差别。
