利用荧光标记酶切基因进行微生物群落结构分析的方法.pdf

一种生物技术领域的利用荧光标记酶切基因进行微生物群落结构分析的方法，包括如下步骤：提取样本的基因组总DNA，扩增微生物系统发育标记基因；利用限制性内切酶对扩增得到的标记基因进行酶切，得到带有粘性末端的限制性片段；对所有限制性片段进行末端荧光标记；运用DNA遗传分析仪分离并检测已被末端荧光标记的限制性片段，得到样本的微生物组成信息图谱；对信息图谱进行多元变量统计学分析，获得微生物群落结构信息特征与变化规律。本发明的方法利用了毛细管电泳的高通量、高灵敏度和高重复性等优点，能够有效全面地获得微生物群落结构信息。

1、  一种利用荧光标记酶切基因进行微生物群落结构分析的方法，其特征在于，包括如下步骤：
步骤一，提取样本的基因组总DNA，扩增微生物系统发育标记基因；
步骤二，利用限制性内切酶对扩增得到的标记基因进行酶切，得到带有粘性末端的限制性片段；
步骤三，对所有限制性片段进行末端荧光标记；
步骤四，运用DNA遗传分析仪分离并检测已被末端荧光标记的限制性片段，得到样本的微生物组成信息图谱；
步骤五，对信息图谱进行多元变量统计学分析，获得微生物群落结构信息特征与变化规律。

2、  根据权利要求1所述的利用荧光标记酶切基因进行微生物群落结构分析的方法，其特征是，步骤二中，所述酶切具体为，采用4核苷酸识别位点的限制性内切酶对标记基因进行酶切，获得具有粘性末端的酶切片段。

3、  根据权利要求1所述的利用荧光标记酶切基因进行微生物群落结构分析的方法，其特征是，步骤三中，所述末端荧光标记具体为，使用DNA polymeraseKlenow酶进行荧光标记，标记所用的底物是与酶切后得到的粘性末端碱基互补配对的脱氧单核苷酸。

4、  根据权利要求1所述的利用荧光标记酶切基因进行微生物群落结构分析的方法，其特征是，步骤四中，所述分离并检测具体为，利用DNA遗传分析仪进行检测。

5、  根据权利要求1所述的利用荧光标记酶切基因进行微生物群落结构分析的方法，其特征是，步骤五中，所述多元变量统计学分析为主成分分析。

利用荧光标记酶切基因进行微生物群落结构分析的方法
技术领域
本发明涉及一种生物技术领域的分析方法，具体是一种利用荧光标记酶切基因进行微生物群落结构分析的方法。
背景技术
微生物是地球上为数最为众多的一类生物，通常以复杂群落的形式存在，并在维系生态平衡、提供自然资源、调节宿主新陈代谢等方面发挥着重要的作用。目前，已有大量微生物被分离、培养和鉴定。但是由于培养方法有着费时费力的缺点以及技术上的局限性，近年来已有大量分子生物学方法被开发并应用于研究复杂微生物群落的结构，并对研究其功能提供指导。常用的微生物分子生物学技术包括基于分子标签(如SSU rRNA基因、ITS序列等)的分子指纹图谱技术，如变性/温度梯度凝胶电泳(DGGE/TGGE)、单链构像多态性(SSCP)、rDNA限制性片段多态性(ADRDA)、末端限制性酶切片段长度多态性(T-RFLP)等；以及用于全面描述微生物群落结构和物种进化关系的克隆文库分析等。DNA指纹图技术能够把从环境样品中PCR扩增得到的微生物群落中所有种群的特征性基因片段以指纹图谱的形式展现出来，从而快速灵敏地检测微生态系统的动态变化，展现微生物群落结构的多样性。
ADRDA和T-RFLP是两种目前最常用的基于限制性酶切的DNA指纹图谱技术。ARDRA技术是通过PCR得到rDNA的基因序列混合物，用限制性内切酶进行消化，并通过凝胶电泳获得图谱，对得到的所有限制性片段多样性进行全面分析。但基于传统凝胶电泳技术使得其具有分辨率低、重复性差以及不便于数据统计分析等缺点，故目前多用于对环境中分离得到的菌株进行属、种水平的分类和菌株的区分，以及对rDNA克隆文库中的克隆进行分型和评价，仅有少数研究将其运用于分析废水、土壤等微生物群落结构；T-RFLP技术使用5’端用荧光标记的引物来扩增目标基因，进行限制性酶切得到酶切片段，然后用DNA测序仪对带有荧光的末端限制性片段进行检测，从而获得代表微生物群落结构多样性的指纹图谱。依赖于毛细管电泳技术的T-RFLP技术具有高灵敏度，高通量及高重复性的优点，且便于多样性图谱的数字化输出和后续分析，故被广泛应用于包括土壤、水体、特殊地质层、生物降解池以及人和动物的肠道菌群在内的各种环境中微生物群落结构的比较分析和动态监测。但因其只对基因的末端限制性片段进行分析，故在对复杂微生物群落进行分析的时候具有一定的局限性。它仅检测末端限制性片段长度的多态性，因而在图谱中，同样长度的末端限制性片段可能来自几种不同的细菌，有的在进化关系中甚至相去甚远。同时，自然环境中的微生物群落多样性往往过于复杂，故T-RFLP图谱所反映的群落组成很难具体到种属水平，因此，在低分类单元水平上的多样性无法体现出来，但这种多样性恰恰是很多复杂样品的特征。运用多种限制性内切酶对同一样品进行检测，并将数据综合进行分析，可以得到微生物群落组成的更丰富的信息，这一手段已被广泛运用于弥补上述不足。但是，运用多种酶切不但费时费力，而且在数据处理的过程中容易引入误差。因此，结合上述两种方法的优点，发展一种能够准确和全面地对微生物群落结构进行分析的方法，具有重要的意义。
发明内容
本发明的目的在于克服现有技术的局限，提供一种利用荧光标记酶切基因进行微生物群落结构分析的方法。本发明的方法利用了毛细管电泳的高通量、高灵敏度和高重复性等优点，能够有效获得全面的微生物群落结构信息。
本发明是通过以下技术方案实现的，具体步骤如下：
步骤一，提取样本的基因组总DNA，扩增微生物系统发育标记基因；
步骤二，利用限制性内切酶对扩增得到的标记基因进行酶切，得到带有粘性末端的限制性片段；
步骤三，对所有限制性片段进行末端荧光标记；
步骤四，运用DNA遗传分析仪分离并检测已被末端荧光标记的限制性片段，得到样本的微生物组成信息图谱；
步骤五，对信息图谱进行多元变量统计学分析，获得微生物群落结构信息特征与变化规律。
步骤二中，所述酶切具体为，采用4核苷酸识别位点的限制性内切酶对标记基因进行酶切，获得具有粘性末端的酶切片段。
步骤三中，所述末端荧光标记具体为，使用DNA polymerase Klenow酶进行荧光标记，标记所用的底物是与酶切后得到的粘性末端碱基互补配对的脱氧单核苷酸。
步骤四中，所述分离并检测具体为，利用DNA遗传分析仪进行检测。
步骤五中，所述多元变量统计学分析为主成分分析。
本发明涉及的是一种分析微生物群落结构的DNA指纹图谱方法，具体是选取具有微生物分类特征的特定基因片段作为研究对象，使用限制性内切酶对其进行酶切，并对所有的限制性片段进行荧光标记，通过DNA片段荧光分析仪器对所有片段进行检测和分析获得DNA多样性指纹图谱，用于比较研究多样本的微生物群落的结构。
本发明具有如下的有益效果：本发明的方法利用了毛细管电泳的高通量、高灵敏度和高重复性等优点，能够有效全面地获得微生物群落结构信息。
附图说明
图1DNA样品以ADRDA和本发明方法分析结果的比较；
图2第一组模拟微生物群落用本发明方法分析的结果；
图3第二组模拟微生物群落用本发明方法分析的结果；
图4动物肠道微生物群落用本发明方法和T-RFLP法分析结果的比较。
具体实施方式
下面对本发明的实施例做详细说明：本实施例在以本发明技术方案为前提下进行实施，给出了详细的实施方式和具体的操作过程，但本发明的保护范围不限于下述的实施例。
实施例1
下面通过对模拟微生物群落进行分析来进一步描述本发明
1、16S rRNA基因代表序列的选取
本实施例提取样品的基因组总DNA，并选取16S rRNA基因(或其它微生物系统发育标记基因，如gyrB、rnpB等)作为细菌多样性的标记基因，采用合适的通用引物扩增该标记基因片段。其中16S rRNA基因是目前分子生态学中最常用的系统进化标记基因，包含丰富的进化信息，且具有强大的数据库作为后续分析的支撑。本实例以16S rRNA基因为代表。
本实施例所选取的16S rRNA基因代表序列选自李旻等人报道的中国人肠道菌群中的代表序列，共12条序列，记为A～F和H～M，其中5条来自拟杆菌门，4条来自厚壁菌门，3条来自变形菌门。该12条序列的Genbank登录号依次为EF403826，EF403827，EF403833，EF403836，EF403840，EF403957，EF404131，EF404244，EF404279，EF404297，EF404313，EF404522；同时这12条序列已经在《Li M，Wang B，Zhang M，Rantalainen M，Wang S，Zhou H，Zhang Y，Shen J，Pang X，Zhang M，Wei H，Chen Y，Lu H，Zuo J，Su M，Qiu Y，Jia W，Xiao C，Smith LM，Yang S，Holmes E，Tang H，Zhao G，Nicholson JK，Li L，Zhao L.Symbiotic gut microbes modulate human metabolic phenotypes.Proc NatlAcad Sci USA.2008Feb 12；105(6)：2117-22》(共生调节的人体肠道微生物的代谢表型，美国科学院院报)中公开。将这些序列克隆在常规TA克隆载体pGEMT质粒上，记为PA～PF，PH～PM。具体方法为：将携带有质粒的E.coli DH5α细胞于LB培养基中培养后，用试剂盒抽取质粒DNA作为以下研究对象。
2、对单序列进行片段分析
将质粒DNA以16S rRNA基因全长通用引物27f(5’-AGAGTTTGATCCTGGCTCAG)和1492r(5’-GGTTACCTTGTTACGACTT)进行PCR扩增。采用具有4核苷酸识别位点的限制性内切酶MspI(5’-C^CGG)对16S rRNA基因扩增产物进行酶切，得到带有粘性末端的限制性片段，这些片段可作为样品的微生物群落多样性的指标，当一个样品的多样性较高时可能会产生更多种类的酶切片段；产生粘性末端保证了后续标记步骤的进行。在Klenow酶的作用下，将荧光标记的Cy5-dCTP(或其它与酶切片段粘性末端核苷酸碱基互补的荧光标记脱氧三磷酸核苷酸)通过碱基互补连接到粘性末端，使DNA片段进行荧光标记。
运用Beckman CEQ8000 DNA遗传分析仪(或其它可对荧光标记的DNA片段进行分离和片段大小测定的仪器)分离检测经酶切和荧光标记的限制性片段，得到DNA片段大小和丰度组成的信息图谱。
对前述12条序列用本方法进行分析；同时根据这12条16S rRNA基因的序列信息和内切酶MspI的酶切位点(5’-C^CGG)对生成片段进行模拟，得到相应限制性片段的数量和长度的理论值；并对这12条序列进行ARDRA分析。将实验得到的片段信息与模拟结果、ARDRA分析结果作比较。
根据模拟酶切，所有片段都位于0～800bp区间内。考虑到测序仪本身的检测局限，以及去除引物二聚体的影响，本实施例选取70～600bp区间内的片段，比较模拟结果与实验结果得知，在模拟中出现的该区间的69个片段全都能够被检测到，来自各代表序列的片段数目与模拟值一致。ARDRA分析得到片段数目也与本方法检测结果一致(图1a)。图1为DNA样品以ADRDA和本实施例分析结果的比较图，图中，a单条16S rRNA基因片段样本，b动物肠道微生物群落DNA样本。
作为比ARDRA更精细的检测方法，本实施例能够较为准确地测量出片段的长度。经比较发现，片段长度的测量值与模拟值的结果一致。来自不同序列的模拟长度相同的片段，其长度测量值的差别不超过1bp；对于长度在100bp以上的片段，62％的片段其测量值与模拟值之间的误差为±1bp以内，92％为±2bp以内。而在长度在100bp以内的片段，误差范围较大，仅有42％的片段其测量值与模拟值之间的误差为±2bp以内。基于此，为了减小检测误差，准确表现群落结构，本实施例100～600bp范围内的片段信息是准确可靠的，并在后面的分析中，选取本区间的片段用于数据分析。
3、构建模拟微生物群落并对结构进行分析
用前述含有单条序列的质粒构建两组模拟微生物群落。
第一组由PA和PE组成，编号为P1～P8八个组合，其构建方法是向100ng质粒PE中依次加入质粒PA，其质量为10ng，5ng，2ng，1.25ng，1ng，0.5ng，0.2ng，0.1ng，换算成质粒拷贝数即为向拷贝数为4.13*10¹⁰的质粒PE中依次加入拷贝数为4.13*10⁹，2.08*10⁹，8.26*10⁸，5.17*10⁸，4.13*10⁸，2.08*10⁸，8.26*10⁷，4.13*10⁷的质粒PA，其中质粒PA所占百分比依次为：9.09％，4.76％，1.96％，1.23％，0.99％，0.50％，0.20％，0.10％，比例依次减小；
第二组由质粒PA，PB，PE，PF，PL按照不同比例混合构建，记为M1～M5，共五个模拟群落组合，具体组成见表1。
用本实施例前述单个质粒分析实验步骤检测两组模拟微生物群落结构，得到模拟微生物群落DNA酶切片段的长度和丰度值。根据模拟群落中16S rRNA基因序列信息和酶切位点对生成片段进行模拟，得到相应限制性片段的数量和长度的理论值，并与实际检测结果作比较；以来自同一质粒的片段的平均丰度来表示该质粒在群落中含量的检测值，计算其在混合群落中的比例，并将其与各质粒的实际比例进行比较，从而确定本实施例对微生物群落结构组成分析的准确性。
表1各质粒在模拟微生物群落M1～M5中所占的比例