基于16SrDNA深度测序检测不同大豆根际土壤原核微生物的方法
【技术领域】
[0001 ]本发明属于土壤微生物学技术领域，具体涉及大田种植方式下，基于16SrDNA深度测序检测不同大豆根际土壤原核微生物的方法。
【背景技术】
[0002]根际(Rhizosphere)被定义为紧密围绕着有生命力根系的生物活性区域，其中包含大量微生物，例如细菌、放线菌和真菌等，也包括一些藻类和病毒等，范围一般为距离根表面Icm以内。不同植物甚至同一植物的不同生育期，其根系分泌物会对根际微生物的生长发育产生影响;而根际微生物又参与土壤有机质的分解、腐殖质的形成、养分的转化和循环等多种生理生化过程，尤其是其中的固氮菌群，通过生物固氮作用在常温下把氮气转化为能被植物利用于合成各种氨基酸的氨。土壤根际微生物群落不仅受到植物根系影响，也受土壤类型、农业耕作管理、季节变化、重金属污染、杀虫剂处理等影响，同时也受到检测技术的影响。因此，准确检测不同基因型植物包括转基因作物的根际微生物群落组成、结构、多样性及相对丰度差异，是判断和评价不同基因型植物对土壤微环境影响的重要基础，也是评估转基因作物释放对土壤微生态风险的重要基础。
[0003]传统的用于研究微生物的分离培养法主要有平板培养法和显微镜观察等方法，虽然简单方便，至今仍然是一种不可替代的研究手段，但也存在一些弊端，即上述方法无法研究土壤中大量的且不可培养的微生物。B1LOG微平板分析法适用于检测可培养的且生长速度较快的微生物及其代谢多样性，但是实验结果并不能完全代表原位条件下微生物代谢多样性的结果，样品的处理方法、培养条件以及B1LOG微平板的类型等都会对结果造成一定程度的误差。醌指纹法的优点是简单快速，因此也被广泛地应用于微生物细胞结构多样性的分析中，但是醌指纹法也存在一些不足，例如它不能具体反映到哪个属、种的多样性变化等。磷脂脂肪酸(PLFA)分析法是基于PLFA是活体细胞膜的一个重要组成部分，含量大约占细胞干重的5% ；用PLFA的变化来反映土壤生态系统中微生物的种类组成和数量变化，该方法不需要进行微生物培养，可直接鉴定微生物群落多样性变化，但该方法与PLFA是否完全提取以及操作过程中PLFA是否稳定有着很大的关系。
[0004]土壤微生物的遗传多样性是指土壤生态系统中的微生物在基因水平上所携带的遗传物质和遗传信息的总和;与高等生物相比，微生物的遗传多样性更加突出，不同微生物种群间的遗传物质和遗传信息具有非常显著的差异。土壤微生物遗传多样性的分析方法主要包括两类:第一类是基于杂交的分析方法，第二类是基于PCR的分析方法。杂交分析方法是基于核酸分子的碱基互补配对原理，利用特异性的探针与待测样品的核酸分子进行杂交的过程;杂交技术是通过测定DNA重组率或DNA杂交动力相似度来研究土壤微生物的遗传多样性，用于研究微生物遗传多样性的杂交探针主要包括rRNA基因探针、抗性基因探针和编码代谢酶基因探针等，其中应用比较普遍的杂交技术包括焚光原位杂交(Fluorescence inSitu Hybridizat1n,FISH)技术，肽核酸(peptide nucleic acids，PNA)杂交技术以及限制性片段长度多态性(restrict1n fragment length polymorphism,RFLP)技术等，各有优点，但均是低通量的分析方法。
[0005]运用PCR技术扩增细菌基因组中的16S rDNA或是真菌基因组中的18S rDNA、ITS区域可以对微生物进行分类鉴定，该技术已经普遍应用到探讨微生物的多样性、鉴定微生物的种类以及系统发育关系等领域。基于PCR的分析方法大致分为三阶段:第一，样品的核酸提取和纯化;第二，目的基因的PCR扩增;第三，扩增片段的变性处理与分析检测，包括聚丙烯酰胺变性梯度凝胶电泳(DGGE)、单链构象多态性(SSCP)、末端限制性片段长度多态性分析(T-RFLP)、核糖体基因间隔区分析(RISA)和自动核糖体间隔区基因分析(ARISA)，或者第三，扩增片段的测序与分析，例如16S rDNA高通量测序技术等。其中，实时荧光定量PCR具有很多优点，如高精度、高灵敏度、高特异性、具有实时性等，但同时也存在一些不足之处，例如成本比较高，操作时还要考虑同源和异源DNA背景、寡核苷酸杂交特异性、荧光染料或特异性探针的浓度以及PCR产物的大小等因素，这些因素都会引起结果的偏差。DGGE是一种利用聚丙烯酰胺变性凝胶电泳的方法将具有相同长度但不同序列组成的核酸片段分离开的技术，虽然已经普遍应用到土壤微生物群落多样性分析的研究领域中，但是该技术也存在一些不足之处，例如可重复性差，耗时耗力，而且只能检测土壤中巨量复杂微生物中的很小一部分等。16S rDNA高通量测序技术，也称为16S rDNA深度测序技术，其原理是因为16SrDNA是原核微生物分类学研究中最常用的“分子钟”，其序列包含9个高可变区(hypervariable reg1n)和8个保守区(constant reg1n)，可变区序列因原核微生物物种不同而异，且变异程度与原核微生物的系统进化密切相关即提取样品中微生物的总DNA，然后用16S rRNA基因(S卩16S rDNA)的通用引物进行PCR扩增，获得绝大部分微生物16S rDNA高可变区的扩增产物，构建扩增产物的文库，进行大规模、高通量测序，然后比较分析测序数据，对土壤微生物群落的多样性进行研究。
【发明内容】
[0006]本发明需要解决的问题是克服现有技术检测通量低、规模小、对不可培养微生物的覆盖度低等不足，提供一种基于16SrDNA深度测序检测不同大豆根际土壤原核微生物的方法，为不同基因型大豆包括转基因大豆及其受体品种对根际土壤原核微生物群落影响的准确评估提供了一种新方法。
[0007]本发明的技术方案:
[0008]本发明所述基于16SrDNA深度测序检测不同大豆根际土壤原核微生物的方法，由以下步骤构成:
[0009](I)大田试验地设计，每一种基因型大豆种植于三个或三个以上小区，每个小区面积不小于12m2(4mX3m)，每个小区有两个取样点，每个取样点取两棵或两棵以上植株；
[0010](2)先除去地面杂草和枯枝落叶，铲除大田土壤表面l-2cm的表土，然后将在盛花期或苗期或鼓粒期或成熟期的大豆植株连根从土中取出，用抖落法取根系抖落土作为系统对照，混合均匀并去除根毛，于零下20°C_零下70°C保存，再用捋取法取下紧粘于大豆根的根际土壤，混合均匀并去除根毛，并于零下70°C保存；
[0011](3)从上述各土壤样品中，用可去除土壤中残留的腐殖质及几乎所有其他的PCR抑制因子的PowerSoil DNA Isolat1n Kit(MoB1 Laboratories Inc.,Carlsbad,CA,USA)提取高质量微生物宏基因组DNA;
[0012](4)用双标签融合引物对，其一含P5 Illumina接头序列、8核苷酸的标签以及特异性引物515F(5 ‘ -GTGCCAGCMGCCGCGGTAA-3，)，另一引物含P7 Il Iumina接头序列、8核苷酸的标签以及806R(5‘-GGACTACHVGGGTWTCTAAT-3’)，做PCR扩增所述宏基因组DNA中 16SrDNA的V4高变区片段以构建文库；
[0013](5)对合格的文库用Illumina Miseq平台进行双末端250个核苷酸(英语缩写为PE250)边合成边测序的深度测序，并对边合成边测序的读取序列(英语为READS)进行质量控制；
[0014](6)把纯净的成对READS拼接成TAG，然后聚类成“可操作分类学单元”(英语缩写为OTU)，再做物种组成、结构、多样性及相对丰度差异的显著性分析；
[0015](7)以根系抖落土壤微生物群落的组成、结构及多样性作为系统对照，克服土壤异质性的影响，准确测定根际土壤原核微生物群落组成、结构、多样性、相对丰度，比较不同大豆之间的异同。
[0016]上述系统对照为大田种植方式下的大豆根系抖落土壤，每种基因型大豆的系统对照土壤样品至少3个生物学重复，每种基因型大豆的根际土壤样品至少3个生物学重复。所述的高质量宏基因组DNA用PowerSoil DNA Isolat1n Kit提取。其中的土壤样品的均质化在康宁LSE涡旋混合器上以最大转速2850rpm运转10分钟完成。
[0017]深度测序是在Illumina Miseq平台上，以PE250模式做16S rDNA的V4高变区片段高通量测序，每个系统对照土壤样品DNA至少产出3.8万条有效TAG以用于聚类成0TU;每个根际土壤样品DNA至少产出13万条有效TAG以用于聚类成0TU。
[0018]本发明的有益效果是:16SrDNA深度测序方法具有很多优点:首先是大规模、高通量，可以检测到土壤中的绝大部分不可培养的原核微生物种类，其次相比于宏基组高通量测序方法，费用较低(700?800元/样品DNA)，因而可以更加全面(不同作物品种植前、苗期、花期、成熟期等多个时期；同一时期各品种多点样品)、直接、原位地研究不同作物对于土壤原核微生物群落组成、结构和遗传多样性的影响。准确检测不同基因型植物包括转基因作物的根际微生物群落组成、结构、多样性及相对丰度差异，是判断和评价不同基因型植物对土壤微环境影响的重要基础，也