CN107904296A - 一种基于深度测序检测玉米根际土壤微生物的方法及系统 - Google Patents
一种基于深度测序检测玉米根际土壤微生物的方法及系统 Download PDFInfo
- Publication number
- CN107904296A CN107904296A CN201711257675.4A CN201711257675A CN107904296A CN 107904296 A CN107904296 A CN 107904296A CN 201711257675 A CN201711257675 A CN 201711257675A CN 107904296 A CN107904296 A CN 107904296A
- Authority
- CN
- China
- Prior art keywords
- rhizosphere
- soil
- crops
- edaphon
- deep sequencing
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 41
- 238000012350 deep sequencing Methods 0.000 title claims abstract description 35
- 238000001514 detection method Methods 0.000 title claims abstract description 27
- 239000002689 soil Substances 0.000 claims abstract description 78
- 244000005700 microbiome Species 0.000 claims abstract description 43
- 238000004458 analytical method Methods 0.000 claims abstract description 29
- 230000009261 transgenic effect Effects 0.000 claims abstract description 20
- 108020004465 16S ribosomal RNA Proteins 0.000 claims abstract description 19
- 239000000203 mixture Substances 0.000 claims abstract description 19
- 241000209140 Triticum Species 0.000 claims abstract description 17
- 235000021307 Triticum Nutrition 0.000 claims abstract description 17
- 235000005824 Zea mays ssp. parviglumis Nutrition 0.000 claims abstract description 16
- 235000002017 Zea mays subsp mays Nutrition 0.000 claims abstract description 16
- 235000005822 corn Nutrition 0.000 claims abstract description 16
- 239000003513 alkali Substances 0.000 claims abstract description 15
- 241000894007 species Species 0.000 claims abstract description 11
- 241000196324 Embryophyta Species 0.000 claims description 23
- 239000002773 nucleotide Substances 0.000 claims description 14
- 125000003729 nucleotide group Chemical group 0.000 claims description 14
- 238000012408 PCR amplification Methods 0.000 claims description 11
- 239000012634 fragment Substances 0.000 claims description 8
- 230000004927 fusion Effects 0.000 claims description 8
- 230000009514 concussion Effects 0.000 claims description 6
- 238000007399 DNA isolation Methods 0.000 claims description 5
- 239000002585 base Substances 0.000 claims description 5
- 238000002474 experimental method Methods 0.000 claims description 4
- 239000007788 liquid Substances 0.000 claims description 4
- 239000000843 powder Substances 0.000 claims description 4
- 238000012360 testing method Methods 0.000 claims description 4
- 238000003908 quality control method Methods 0.000 claims description 3
- 238000012163 sequencing technique Methods 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 2
- 241000209149 Zea Species 0.000 claims 2
- 240000008042 Zea mays Species 0.000 abstract description 14
- 238000005516 engineering process Methods 0.000 abstract description 9
- 208000003643 Callosities Diseases 0.000 abstract description 6
- 206010020649 Hyperkeratosis Diseases 0.000 abstract description 6
- 244000068988 Glycine max Species 0.000 abstract description 2
- 235000010469 Glycine max Nutrition 0.000 abstract description 2
- 230000007613 environmental effect Effects 0.000 abstract description 2
- 108020004414 DNA Proteins 0.000 description 18
- 239000000523 sample Substances 0.000 description 18
- 210000004027 cell Anatomy 0.000 description 10
- 230000008901 benefit Effects 0.000 description 5
- 244000037671 genetically modified crops Species 0.000 description 5
- 239000010977 jade Substances 0.000 description 5
- AZQWKYJCGOJGHM-UHFFFAOYSA-N 1,4-benzoquinone Chemical compound O=C1C=CC(=O)C=C1 AZQWKYJCGOJGHM-UHFFFAOYSA-N 0.000 description 4
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 4
- 230000008859 change Effects 0.000 description 4
- 238000011156 evaluation Methods 0.000 description 4
- 230000000813 microbial effect Effects 0.000 description 4
- 238000000513 principal component analysis Methods 0.000 description 4
- 108090000623 proteins and genes Proteins 0.000 description 4
- 241000894006 Bacteria Species 0.000 description 3
- 108091007491 NSP3 Papain-like protease domains Proteins 0.000 description 3
- 238000012165 high-throughput sequencing Methods 0.000 description 3
- QGZKDVFQNNGYKY-UHFFFAOYSA-N Ammonia Chemical compound N QGZKDVFQNNGYKY-UHFFFAOYSA-N 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 238000010835 comparative analysis Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 230000004907 flux Effects 0.000 description 2
- 238000001502 gel electrophoresis Methods 0.000 description 2
- 239000003864 humus Substances 0.000 description 2
- 238000011065 in-situ storage Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000002906 microbiologic effect Effects 0.000 description 2
- 229910052757 nitrogen Inorganic materials 0.000 description 2
- 238000004321 preservation Methods 0.000 description 2
- 125000006850 spacer group Chemical group 0.000 description 2
- 238000011282 treatment Methods 0.000 description 2
- 108020004463 18S ribosomal RNA Proteins 0.000 description 1
- HRPVXLWXLXDGHG-UHFFFAOYSA-N Acrylamide Chemical compound NC(=O)C=C HRPVXLWXLXDGHG-UHFFFAOYSA-N 0.000 description 1
- 241000186046 Actinomyces Species 0.000 description 1
- 108010043325 Aryl-alcohol dehydrogenase Proteins 0.000 description 1
- 108010049668 Betaine-Aldehyde Dehydrogenase Proteins 0.000 description 1
- 241000195493 Cryptophyta Species 0.000 description 1
- 238000007400 DNA extraction Methods 0.000 description 1
- 108700005088 Fungal Genes Proteins 0.000 description 1
- 241000233866 Fungi Species 0.000 description 1
- 108091028043 Nucleic acid sequence Proteins 0.000 description 1
- 108091034117 Oligonucleotide Proteins 0.000 description 1
- 240000007594 Oryza sativa Species 0.000 description 1
- 235000007164 Oryza sativa Nutrition 0.000 description 1
- 241000425347 Phyla <beetle> Species 0.000 description 1
- 108020001027 Ribosomal DNA Proteins 0.000 description 1
- 238000012300 Sequence Analysis Methods 0.000 description 1
- 241000700605 Viruses Species 0.000 description 1
- 239000002253 acid Substances 0.000 description 1
- 230000000735 allogeneic effect Effects 0.000 description 1
- 150000001413 amino acids Chemical class 0.000 description 1
- 229910021529 ammonia Inorganic materials 0.000 description 1
- 238000012098 association analyses Methods 0.000 description 1
- 230000001580 bacterial effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004178 biological nitrogen fixation Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 239000003398 denaturant Substances 0.000 description 1
- 238000004925 denaturation Methods 0.000 description 1
- 230000036425 denaturation Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 235000014113 dietary fatty acids Nutrition 0.000 description 1
- 239000000686 essence Substances 0.000 description 1
- 238000013401 experimental design Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 210000000416 exudates and transudate Anatomy 0.000 description 1
- 229930195729 fatty acid Natural products 0.000 description 1
- 239000000194 fatty acid Substances 0.000 description 1
- 150000004665 fatty acids Chemical class 0.000 description 1
- 239000007850 fluorescent dye Substances 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 239000011544 gradient gel Substances 0.000 description 1
- 229910001385 heavy metal Inorganic materials 0.000 description 1
- 238000009396 hybridization Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000004060 metabolic process Effects 0.000 description 1
- 230000007269 microbial metabolism Effects 0.000 description 1
- 244000005706 microflora Species 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 150000007523 nucleic acids Chemical group 0.000 description 1
- 235000015097 nutrients Nutrition 0.000 description 1
- 239000000575 pesticide Substances 0.000 description 1
- 238000007747 plating Methods 0.000 description 1
- 229920002401 polyacrylamide Polymers 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000003753 real-time PCR Methods 0.000 description 1
- 238000007894 restriction fragment length polymorphism technique Methods 0.000 description 1
- 210000003705 ribosome Anatomy 0.000 description 1
- 235000009566 rice Nutrition 0.000 description 1
- 230000001932 seasonal effect Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 210000001082 somatic cell Anatomy 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
- 238000003971 tillage Methods 0.000 description 1
Classifications
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
- C12Q1/6869—Methods for sequencing
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Organic Chemistry (AREA)
- Zoology (AREA)
- Wood Science & Technology (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Microbiology (AREA)
- Immunology (AREA)
- Biotechnology (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Analytical Chemistry (AREA)
- Physics & Mathematics (AREA)
- Biochemistry (AREA)
- Bioinformatics & Cheminformatics (AREA)
- General Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Genetics & Genomics (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
本发明属于土壤微生物学技术领域，公开了一种基于深度测序检测玉米根际土壤微生物的方法及系统，采集根际土壤；从土壤中提取高质量微生物基因组DNA；用双标签引物做PCR扩增DNA中16S rDNA V4区,构建文库；对合格文库进行测序，得到纯净READS；5.拼接，每样品产出至少5.0万条有效TAG用于聚类成可操作分类单元；6.物种组成、结构、多样性及相对丰度差异显著性分析；7.准确测定转基因耐盐碱玉米根际土壤原核微生物群落组成、结构、多样性、相对丰度。本发明受国家转基因重大专项课题“转基因玉米小麦大豆环境安全评价技术(2016ZX08011003)”支持。
Description
技术领域
本发明属于土壤微生物学技术领域，尤其涉及一种基于深度测序检测玉米根际土壤微生物的方法及系统。
背景技术
根际(Rhizosphere)被定义为紧密围绕着有生命力根系的生物活性区域，其中包含大量微生物，例如细菌、放线菌和真菌等，也包括一些藻类和病毒等，范围一般为距离根表面1cm以内。不同植物甚至同一植物的不同生育期，其根系分泌物会对根际微生物的生长发育产生影响；而根际微生物又参与土壤有机质的分解、腐殖质的形成、养分的转化和循环等多种生理生化过程，尤其是其中的固氮菌群，通过生物固氮作用在常温下把氮气转化为能被植物利用于合成各种氨基酸的氨。土壤根际微生物群落不仅受到植物根系影响，也受土壤类型、农业耕作管理、季节变化、重金属污染、杀虫剂处理等影响，同时也受到检测技术的影响。因此，准确的获得根际土土样，进而准确检测不同基因型植物包括转基因作物的根际微生物群落组成、结构、多样性及相对丰度差异，是判断和评价不同基因型植物对土壤微环境影响的重要基础，也是评估转基因作物释放对土壤微生态风险的重要基础。
传统的用于分析微生物的分离培养法主要有平板培养法和显微镜观察等方法，但存在很多弊端，即上述方法无法分析土壤中大量的且不可培养的微生物。 BIOLOG微平板分析法适用于检测可培养的且生长速度较快的微生物及其代谢多样性，但是实验结果并不能完全代表原位条件下微生物代谢多样性的结果，样品的处理方法、培养条件以及BIOLOG微平板的类型等都会对结果造成一定程度的误差。醌指纹法的优点是简单快速，因此也被广泛地应用于微生物细胞结构多样性的分析中，但是醌指纹法不能具体反映到哪个属、种的多样性变化等。磷脂脂肪酸(PLFA)分析法是基于PLFA是活体细胞膜的一个重要组成部分，含量大约占细胞干重的5％；用PLFA的变化来反映土壤生态系统中微生物的种类组成和数量变化，该方法不需要进行微生物培养，可直接鉴定微生物群落多样性变化，但该方法与PLFA是否完全提取以及操作过程中PLFA是否稳定有着很大的关系。
而运用PCR技术扩增细菌基因组中的16S rDNA或是真菌基因组中的18S rDNA、ITS区域可以对微生物进行分类鉴定，该技术已经普遍应用到探讨微生物的多样性、鉴定微生物的种类以及系统发育关系等领域。基于PCR的分析方法大致分为三阶段：第一，样品的核酸提取和纯化；第二，目的基因的PCR扩增；第三，扩增片段的变性处理与分析检测，包括聚丙烯酰胺变性梯度凝胶电泳(DGGE)、单链构象多态性(SSCP)、末端限制性片段长度多态性分析 (T-RFLP)、核糖体基因间隔区分析(RISA)和自动核糖体间隔区基因分析 (ARISA)，或者第三，扩增片段的测序与分析，例如16S rDNA高通量测序技术等。其中，实时荧光定量PCR具有很多优点，如高精度、高灵敏度、高特异性、具有实时性等，但同时也存在一些不足之处，例如成本比较高，操作时还要考虑同源和异源DNA背景、寡核苷酸杂交特异性、荧光染料或特异性探针的浓度以及PCR产物的大小等因素，这些因素都会引起结果的偏差。DGGE是一种利用聚丙烯酰胺变性凝胶电泳的方法将具有相同长度但不同序列组成的核酸片段分离开的技术，虽然已经普遍应用到土壤微生物群落多样性分析的分析领域中，但是该技术也存在一些不足之处，例如可重复性差，耗时耗力，而且只能检测土壤中巨量复杂微生物中的很小一部分等。
综上所述，现有技术存在的问题是：
现有技术检测通量低、规模小、对不可培养微生物的覆盖度低；不能准确检测不同基因型植物包括转基因作物的根际微生物群落组成、结构、多样性及相对丰度差异，同时传统的利用抖落法获得的土壤根际土实际上包含了大量的非根际土，导致后续鉴定检测结果不准确。
发明内容
针对现有技术存在的问题，本发明提供了一种基于深度测序检测玉米根际土壤微生物的方法及系统。本发明是基于涡旋震荡分离根际土，进而基于16SrDNA 深度测序检测不同玉米根际土壤原核微生物的方法，为转基因耐盐碱玉米及其受体品种对根际土壤原核微生物群落影响的准确评估提供了一种新方法。
本发明是这样实现的，一种基于深度测序检测玉米根际土壤微生物的方法，所述基于深度测序检测玉米根际土壤微生物的方法包括：
利用涡旋震荡准确分离玉米根际土，提取DNA；
用双标签融合引物对做PCR扩增转基因耐盐碱玉米种植区土壤样品中提取的微生物宏基因组DNA中的16S rDNA V4区，构建文库；
对合格文库进行深度测序，得到纯净READS；
拼接成TAG，土壤样品产出至少5.0万条有效TAG用于聚类成可操作分类单元OTU。
进一步，所述引物对，其中一含P5Illumina接头序列、8核苷酸的标签以及特异性引物515F，515F的核苷酸序列是SEQ ID NO.1；
另一引物含P7Illumina接头序列、8核苷酸的标签以及特异性引物806R， 806R的核苷酸序列是SEQ ID NO.2。
进一步，对合格的文库用IlluminaHiseq PE250进行深度测序，并对读取序列READS，剔除含有Primer/Adaptor的READS；剔除含有过多non-ATCG碱基 N的READS；剔除测序质量较低的碱基数占的比例过高的READS，获得纯净的成对READS。
进一步，把纯净的成对READS拼接成TAG，然后聚类成可操作分类学单元OTU后，再做物种组成、结构、多样性及相对丰度差异的显著性分析。根据得到的群落丰度数据，进行稀有频率数据的多重假设检验和假发现率(FDR) 分析进行差异的显著性分析。
进一步，做物种组成、结构、多样性及相对丰度差异的显著性分析后，还进行：确定根际土壤原核微生物群落组成、结构、多样性、相对丰度，利用方差分析比较不同玉米取样时期之间的异同。
进一步，用双标签融合引物对做PCR扩增转基因耐盐碱玉米种植区土壤样品中提取的微生物宏基因组DNA中，宏基因组DNA用PowerSoil DNAIsolation Kit提取，其中，土壤样品的均质化在涡旋混合器上以转速3000rpm运转15分钟。
进一步，所述的深度测序在Illumina Miseq平台上，以PE250模式做16S rDNA的V4高变区片段高通量测序。
进一步，用双标签融合引物对做PCR扩增转基因耐盐碱玉米种植区土壤样品中提取的微生物宏基因组DNA中的16S rDNA V4区，之前需进行：
将转基因耐盐碱玉米种植于大田试验地三个或三个以上小区，每个小区面积不小于20m2，每个小区有至少3个取样点，每个取样点取至少3棵转基因耐盐碱玉米植株；
先除去地面杂草和枯枝落叶，铲除大田土壤表面1cm-2cm的表土，然后将苗期、拔节期、抽丝散粉期、成熟期的玉米植株连根从土中取出，用抖落法去除根系抖落土，将细根放入装有25mlPBS液的50mL离心管中，利用涡旋震荡将根际土从根上分离，混合均匀，并于-80℃保存。
本发明的另一目的在于提供一种基于深度测序检测玉米根际土壤微生物的系统。
本发明的优点及积极效果为：
本发明提供的根际土获取方式，能极大的提高根际土获取准确性，利用漩涡震荡能完全将根际土从根上分离下来(图2)；
本发明提供的大田种植方式下的玉米根际土壤中，每种土壤样品至少3个生物学重复；每种根际土壤样品至少3个生物学重复。所述的高质量宏基因组DNA用PowerSoil DNAIsolation Kit提取，其中，土壤样品的均质化在涡旋混合器上以最大转速3000rpm运转15分钟完成。深度测序是在Illumina Miseq 平台上，以PE250模式做16S rDNA的V4V5高变区片段高通量测序，每个系统土壤样品DNA至少产出5.0万条有效TAG以用于聚类成OTU。
本发明的有益效果是：16S rDNA深度测序方法具有很多优点：首先是大规模、高通量，可以检测到土壤中的绝大部分不可培养的原核微生物种类，其次相比于宏基组高通量测序方法，费用较低(300～500元/样品DNA)，因而可以更加全面(不同作物品种植前、苗期、花期、成熟期等多个时期；同一时期各品种多点样品)、直接、原位地分析不同作物对于土壤原核微生物群落组成、结构和遗传多样性的影响。准确检测不同基因型植物包括转基因作物的根际微生物群落组成、结构、多样性及相对丰度差异，是判断和评价不同基因型植物对土壤微环境影响的重要基础，也是评估转基因作物释放对土壤微生态风险的重要基础，本发明提供的一种基于16SrDNAV4V5区深度测序基于深度测序检测玉米根际土壤微生物的方法，为准确分离土壤根际土、
转基因耐盐碱玉米及其受体品种对根际土壤原核微生物群落影响的准确评估提供了一种新方法。该方法检测全面，操作简便，成本低,结果可靠。
本发明受国家转基因重大专项课题“转基因玉米小麦大豆环境安全评价技术(2016ZX08011003)”支持。
附图说明
图1是本发明实施例提供的基于深度测序检测玉米根际土壤微生物的方法流程图。
图2是本发明实施例提供的转基因玉根际土壤(SO)高质量微生物宏基因组DNA的检测结果。
图3是本发明实施例提供的转基因玉根际土壤中原核微生物的主成分分析(PCA)结果。
图4是本发明实施例提供的转基因玉根际土壤中原核微生物基于BETA多样性WEIGHTED-UNIFRAC的主关联分析(PCoA)结果。
图5是本发明实施例提供的转基因玉根际土壤中超过99％的主要细菌类别，在门(phyla)分类水平上相对丰度的差异分析。其中*,p<0.05；**,p< 0.01。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
现有技术检测通量低、规模小、对不可培养微生物的覆盖度低；不能准确检测不同基因型植物包括转基因作物的根际微生物群落组成、结构、多样性及相对丰度差异。
16S rDNA高通量测序技术，也称为16S rDNA深度测序技术，其原理是因为16SrDNA是原核微生物分类学分析中最常用的“分子钟”，其序列包含9 个高可变区(hypervariable region)和8个保守区(constant region)，可变区序列因原核微生物物种不同而异，且变异程度与原核微生物的系统进化密切相关即提取样品中微生物的总DNA，然后用16S rRNA基因(即16S rDNA)的通用引物进行PCR扩增，获得绝大部分微生物16S rDNA高可变区的扩增产物，构建扩增产物的文库，进行大规模、高通量测序，然后比较分析测序数据，对土壤微生物群落的多样性进行分析。
下面结合具体实施例对本发明作详细描述。
如图1所示，本发明实施例提供的基于深度测序检测玉米根际土壤微生物的方法，包括：
S101：大田试验地设计，转基因耐盐碱玉米种植于三个或三个以上小区，每个小区面积不小于20m2，每个小区有5个取样点，每个取样点取3棵植株。
S102：先除去地面杂草和枯枝落叶，铲除大田土壤表面1-2cm的表土，然后将苗期、拔节期、抽丝散粉期、成熟期的玉米植株连根从土中取出，用抖落法去除根系抖落土，将细根放入装有25mlPBS液的50mL离心管中，利用涡旋震荡将根际土从根上分离，混合均匀，并于-80℃保存。
S103：从上述各土壤样品中，采用MoBio公司PowerSoil@DNAIsolation Kit 提取高质量微生物宏基因组DNA。
S104：用双标签融合引物对，其一含P5Illumina接头序列、8核苷酸的标签以及特异性引物515F，515F的核苷酸序列是5‘- GTGCCAGCMGCCGCGGTAA-3’，另一引物含P7Illumina接头序列、8核苷酸的标签以及特异性引物806R，806R的核苷酸序列是 5‘-GGACTACHVGGGTWTCTAAT-3’，做PCR扩增所述宏基因组DNA中16S rDNA的V5区片段以构建文库。
S105：对合格的文库用IlluminaHiseq PE250进行深度测序，并对读取序列 READS进行质量控制。
S106：把纯净的成对READS拼接成TAG，然后聚类成“可操作分类学单元”，英语缩写为OTU，再做物种组成、结构、多样性及相对丰度差异的显著性分析。
S107：确测定根际土壤原核微生物群落组成、结构、多样性、相对丰度，比较不同玉米取样时期之间的异同。
所述的系统对照为大田种植方式下的玉米根际土壤，每种土壤样品至少3 个生物学重复；每种根际土壤样品至少3个生物学重复。
所述的高质量宏基因组DNA用PowerSoil DNA Isolation Kit提取，其中，土壤样品的均质化在涡旋混合器上以最大转速3000rpm运转15分钟完成。
所述的深度测序是在Illumina Miseq平台上，以PE250模式做16S rDNA的 V5高变区片段高通量测序，每个系统土壤样品DNA至少产出5.0万条有效TAG 以用于聚类成OTU。
下面结合具体实施例对本发明作进一步描述。
1、玉米材料：郑58；BZ-136含BADH耐盐碱基因玉米。
2、田间试验设计：
大田试验地点在哈尔滨市，每一种玉米种植于三个或三个以上小区，每个小区面积不小于20m2(4m×5m)。每个小区有两个取样点，每个取样点取两棵或两棵以上植株。
3、样品采集：
先除去地面杂草和枯枝落叶，铲除大田土壤表面1-2cm的表土，然后将苗期、拔节期、抽丝散粉期、成熟期的玉米植株连根从土中取出，用抖落法去除根系抖落土，将细根放入装有25mlPBS液的50mL离心管中，利用涡旋震荡将根际土从根上分离，混合均匀，并于-80℃保存；
4、土壤样品中微生物宏基因组DNA的提取：
从上述各土壤样品中，用可去除土壤中残留的腐殖质及几乎所有其他的PCR 抑制因子的PowerSoil DNA Isolation Kit(MoBio Laboratories Inc.,Carlsbad,CA, USA)提取高质量微生物宏基因组DNA(图2)。
5、宏基因组DNA中16SrDNA的V4高变区文库的构建
用双标签融合引物对，其一含P5Illumina接头序列、8核苷酸的标签以及特异性引物515F(SEQ ID NO.1：5‘-GTGCCAGCMGCCGCGGTAA-3’)，另一引物含P7Illumina接头序列、8核苷酸的标签以及806R(SEQ ID NO.2：5 ‘-GGACTACHVGGGTWTCTAAT-3’)，做PCR扩增所述宏基因组DNA中16SrDNA的 V4高变区片段以构建文库。
6、合格文库的高通量测序
对合格的文库用IlluminaHiseq PE250进行深度测序，并对读取序列READS 进行质量控制。
7、生物信息学内容分析
把纯净的成对READS拼接成TAG，然后聚类成“可操作分类学单元”(OTU)；然后做ALPHA多样性、结构、组成及在不同分类水平上相对丰度差异的显著性分析等。
8、数据整合与比较分析
以根际土壤(SO)原核微生物群落的ALPHA多样性、结构(主成分分析，图3； BETA多样性，图4)、组成及在门分类水平上的相对丰度差异的显著性分析(图 5)、在纲分类水平上的相对丰度差异的显著性分析作为系统对照，克服土壤异质性的影响，准确测定根际土壤原核微生物群落的ALPHA多样性、结构(主成分分析，图3；BETA多样性，图4)、组成及在门分类水平上的相对丰度差异的显著性分析(图5)、在纲分类水平上的相对丰度差异的显著性分析以及在转基因玉米根际土壤中主要的物种显著性差异统计，比较转基因耐盐碱玉米根际土壤原核微生物群落之间的异同。
本发明实施例提供的转基因玉米根际土壤中超过95％的主要细菌类别，在纲(classes)分类水平上相对丰度的差异分析中。p<0.05。
本发明实施例提供的转基因玉米根际土壤中主要的物种显著性差异统计中。其中p<0.05。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。
序列表
<110> 东北农业大学
<120> 一种基于深度测序检测玉米根际土壤微生物的方法及系统
<160> 2
<170> SIPOSequenceListing 1.0
<210> 1
<211> 19
<212> DNA
<213> Artificial Sequence
<400> 1
gtgccagcmg ccgcggtaa 19
<210> 2
<211> 20
<212> DNA
<213> Artificial Sequence
<400> 2
ggactachvg ggtwtctaat 20
Claims (9)
1.一种基于深度测序检测玉米根际土壤微生物的方法，其特征在于，所述基于深度测序检测玉米根际土壤微生物的方法包括：
用双标签融合引物对做PCR扩增转基因耐盐碱玉米种植区土壤样品中提取的微生物宏基因组DNA中的16S rDNA V4区，构建文库；
对合格文库进行深度测序，得到纯净READS；
拼接成TAG，土壤样品产出至少5.0万条有效TAG用于聚类成可操作分类单元OTU。
2.如权利要求1所述的基于深度测序检测玉米根际土壤微生物的方法，其特征在于，所述引物对，其中一含P5Illumina接头序列、8核苷酸的标签以及特异性引物515F，515F的核苷酸序列是SEQ ID NO.1；
另一引物含P7Illumina接头序列、8核苷酸的标签以及特异性引物806R，806R的核苷酸序列是SEQ ID NO.2。
3.如权利要求1所述的基于深度测序检测玉米根际土壤微生物的方法，其特征在于，对合格的文库用IlluminaHiseq PE250进行深度测序，并对读取序列READS进行质量控制:
剔除含有Primer/Adaptor的READS；
剔除含有过多non-ATCG碱基N的READS；
剔除测序质量较低的碱基数占的比例过高的READS。
4.如权利要求1所述的基于深度测序检测玉米根际土壤微生物的方法，其特征在于，
把纯净的成对READS拼接成TAG，然后聚类成可操作分类学单元OTU后，再做物种组成、结构、多样性及相对丰度差异的显著性分析；
根据得到的群落丰度数据，进行稀有频率数据的多重假设检验和假发现率FDR分析，评估丰度差异的显著性。
5.如权利要求4所述的基于深度测序检测玉米根际土壤微生物的方法，其特征在于，做物种组成、结构、多样性及相对丰度差异的显著性分析后，还进行：
确定根际土壤原核微生物群落组成、结构、多样性、相对丰度，利用方差分析比较不同玉米取样时期上述指标之间的异同。
6.如权利要求1所述的基于深度测序检测玉米根际土壤微生物的方法，其特征在于，用双标签融合引物对做PCR扩增转基因耐盐碱玉米种植区土壤样品中提取的微生物宏基因组DNA中，宏基因组DNA用PowerSoil DNA Isolation Kit提取，其中，土壤样品的均质化在涡旋混合器上以转速3000rpm运转15分钟。
7.如权利要求1所述的基于深度测序检测玉米根际土壤微生物的方法，其特征在于，所述的深度测序在Illumina Miseq平台上，以PE250模式做16S rDNA的V4高变区片段高通量测序。
8.如权利要求1所述的基于深度测序检测玉米根际土壤微生物的方法，其特征在于，用双标签融合引物对做PCR扩增转基因耐盐碱玉米种植区土壤样品中提取的微生物宏基因组DNA中的16S rDNA V4区，之前需进行：
将转基因耐盐碱玉米种植于大田试验地三个或三个以上小区，每个小区面积不小于20m2，每个小区有至少3个取样点，每个取样点取至少3棵转基因耐盐碱玉米植株；
先除去地面杂草和枯枝落叶，铲除大田土壤表面1cm-2cm的表土，然后将苗期、拔节期、抽丝散粉期、成熟期的玉米植株连根从土中取出，用抖落法去除根系抖落土，将细根放入装有25mlPBS液的50mL离心管中，利用涡旋震荡将根际土从根上分离，混合均匀，并于-80℃保存。
9.一种如权利要求1所述基于深度测序检测玉米根际土壤微生物的方法的基于深度测序检测玉米根际土壤微生物的系统。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711257675.4A CN107904296A (zh) | 2017-12-04 | 2017-12-04 | 一种基于深度测序检测玉米根际土壤微生物的方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711257675.4A CN107904296A (zh) | 2017-12-04 | 2017-12-04 | 一种基于深度测序检测玉米根际土壤微生物的方法及系统 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107904296A true CN107904296A (zh) | 2018-04-13 |
Family
ID=61853882
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711257675.4A Pending CN107904296A (zh) | 2017-12-04 | 2017-12-04 | 一种基于深度测序检测玉米根际土壤微生物的方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107904296A (zh) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109797438A (zh) * | 2019-01-17 | 2019-05-24 | 武汉康测科技有限公司 | 一种用于16S rDNA可变区定量测序文库构建的接头元件及文库构建方法 |
| CN110894527A (zh) * | 2019-11-22 | 2020-03-20 | 东莞市林业科学研究所(广东省林业区域性试验东莞中心、东莞市林业科技推广站) | 一种树木根腐病的评价方法 |
| CN112481409A (zh) * | 2020-12-17 | 2021-03-12 | 东北农业大学 | 一种与玉米耐盐碱性相关的dCAPS分子标记及其应用 |
| CN112575106A (zh) * | 2020-12-17 | 2021-03-30 | 东北农业大学 | 基于转录因子基因ZmNAC89开发的与玉米耐盐碱相关的分子标记及其应用 |
| CN112725481A (zh) * | 2021-01-19 | 2021-04-30 | 贵州省蚕业研究所(贵州省辣椒研究所) | 微生物Pajaroellobacter的新用途 |
| CN113658637A (zh) * | 2021-04-19 | 2021-11-16 | 河南农业大学 | 土壤中玉米秸秆碳同化关键微生物的识别方法 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105525025A (zh) * | 2016-02-17 | 2016-04-27 | 南京大学 | 基于16SrDNA深度测序检测不同大豆根际土壤原核微生物的方法 |
-
2017
- 2017-12-04 CN CN201711257675.4A patent/CN107904296A/zh active Pending
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105525025A (zh) * | 2016-02-17 | 2016-04-27 | 南京大学 | 基于16SrDNA深度测序检测不同大豆根际土壤原核微生物的方法 |
Non-Patent Citations (2)
| Title |
|---|
| JASON A. PEIFFER等: "Diversity and heritability of the maize rhizosphere microbiome under field conditions", 《PNAS》 * |
| 邸宏等: "转BcWRKY1基因玉米耐盐碱性评价及对土壤生态环境影响", 《中国作物学会-2015年学术年会论文摘要集》 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109797438A (zh) * | 2019-01-17 | 2019-05-24 | 武汉康测科技有限公司 | 一种用于16S rDNA可变区定量测序文库构建的接头元件及文库构建方法 |
| CN110894527A (zh) * | 2019-11-22 | 2020-03-20 | 东莞市林业科学研究所(广东省林业区域性试验东莞中心、东莞市林业科技推广站) | 一种树木根腐病的评价方法 |
| CN110894527B (zh) * | 2019-11-22 | 2023-09-19 | 东莞市林业科学研究所(广东省林业区域性试验东莞中心、东莞市林业科技推广站) | 一种树木根腐病的评价方法 |
| CN112481409A (zh) * | 2020-12-17 | 2021-03-12 | 东北农业大学 | 一种与玉米耐盐碱性相关的dCAPS分子标记及其应用 |
| CN112575106A (zh) * | 2020-12-17 | 2021-03-30 | 东北农业大学 | 基于转录因子基因ZmNAC89开发的与玉米耐盐碱相关的分子标记及其应用 |
| CN112725481A (zh) * | 2021-01-19 | 2021-04-30 | 贵州省蚕业研究所(贵州省辣椒研究所) | 微生物Pajaroellobacter的新用途 |
| CN113658637A (zh) * | 2021-04-19 | 2021-11-16 | 河南农业大学 | 土壤中玉米秸秆碳同化关键微生物的识别方法 |
| CN113658637B (zh) * | 2021-04-19 | 2024-02-13 | 河南农业大学 | 土壤中玉米秸秆碳同化关键微生物的识别方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107904296A (zh) | 一种基于深度测序检测玉米根际土壤微生物的方法及系统 | |
| CN105525025A (zh) | 基于16SrDNA深度测序检测不同大豆根际土壤原核微生物的方法 | |
| CN108342465A (zh) | 一种基于高通量测序检测作物根际原核微生物的方法 | |
| Nübel et al. | Quantifying microbial diversity: morphotypes, 16S rRNA genes, and carotenoids of oxygenic phototrophs in microbial mats | |
| CN104711361B (zh) | 快速鉴定西瓜新品种红和平杂交种子纯度的方法以及采用的引物和试剂盒 | |
| CN108486266B (zh) | 玉米叶绿体基因组的分子标记及在品种鉴定中的应用 | |
| CN101736085A (zh) | 牛支原体环介导等温扩增检测方法 | |
| Lee et al. | Polyphasic identification of cyanobacterial isolates from Australia | |
| Joo et al. | Monitoring of phytoplankton community structure using terminal restriction fragment length polymorphism (T-RFLP) | |
| CN112251503A (zh) | 一种基于16S rDNA全长高通量测序检测根际细菌群落的方法 | |
| CN104293778A (zh) | 兰属微卫星标记的建立方法、核心指纹标记库与试剂盒 | |
| CN104673884A (zh) | 利用全基因组和est数据开发多态性est-ssr标记的方法 | |
| Kowalska et al. | DNA barcoding–A new device in phycologist's toolbox | |
| CN106498048A (zh) | 一种与大豆结瘤数相关的qtl、snp分子标记及应用 | |
| CN106521004A (zh) | 与萝卜核雄性不育基因连锁的Indel标记及其应用 | |
| CN108384879A (zh) | 一种用于西瓜杂交品种纯度鉴定的ssr引物及方法 | |
| CN106701967A (zh) | 调控玉米叶夹角主效qtl的分子标记及其应用方法 | |
| Rolland et al. | Molecular identification of broomrape species from a single seed by high resolution melting analysis | |
| Wang et al. | Studying safe storage time of orange peel (Citrus reticulata) using high‐throughput sequencing and conventional pure culture | |
| CN103388030B (zh) | 一种花生不同生长时期根际土壤细菌群落结构的分析方法 | |
| Zhu et al. | Molecular specificity and detection for Pseudanabaena (cyanobacteria) species based on rbcLX sequences | |
| CN109811044A (zh) | 接种菌剂玉米的根际原核微生物多样性检测方法 | |
| CN106676176A (zh) | 一种利用多重pcr对四倍体紫花苜蓿进行ssr分析的方法 | |
| CN106957899A (zh) | 一种利用高通量测序快速分析蓝莓栽培地土壤细菌群落的方法 | |
| CN110669835A (zh) | 基于高通量测序的蕨麻根际微生物多样性及菌群解析方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180413 |