CN111647676A - 利用转录组学鉴定紫花苜蓿与根瘤菌结瘤专一性的方法 - Google Patents

利用转录组学鉴定紫花苜蓿与根瘤菌结瘤专一性的方法 Download PDF

Info

Publication number
CN111647676A
CN111647676A CN202010130201.9A CN202010130201A CN111647676A CN 111647676 A CN111647676 A CN 111647676A CN 202010130201 A CN202010130201 A CN 202010130201A CN 111647676 A CN111647676 A CN 111647676A
Authority
CN
China
Prior art keywords
alfalfa
specificity
nodulation
rhizobia
transcriptomics
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.)
Withdrawn
Application number
CN202010130201.9A
Other languages
English (en)
Inventor
康文娟
师尚礼
吴芳
陈永岗
周彤
尹国丽
祁娟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Gansu Agricultural University
Original Assignee
Gansu Agricultural University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Gansu Agricultural University filed Critical Gansu Agricultural University
Priority to CN202010130201.9A priority Critical patent/CN111647676A/zh
Publication of CN111647676A publication Critical patent/CN111647676A/zh
Withdrawn legal-status Critical Current

Links

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/6876—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes
    • C12Q1/6888—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for detection or identification of organisms
    • C12Q1/6895—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for detection or identification of organisms for plants, fungi or algae
    • 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
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B25/00—ICT specially adapted for hybridisation; ICT specially adapted for gene or protein expression
    • G16B25/10—Gene or protein expression profiling; Expression-ratio estimation or normalisation
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B30/00—ICT specially adapted for sequence analysis involving nucleotides or amino acids
    • G16B30/10—Sequence alignment; Homology search
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B50/00—ICT programming tools or database systems specially adapted for bioinformatics
    • G16B50/30—Data warehousing; Computing architectures
    • 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
    • C12Q2600/00—Oligonucleotides characterized by their use
    • C12Q2600/13—Plant traits
    • 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
    • C12Q2600/00—Oligonucleotides characterized by their use
    • C12Q2600/158—Expression markers

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Proteomics, Peptides & Aminoacids (AREA)
  • Organic Chemistry (AREA)
  • Biotechnology (AREA)
  • Analytical Chemistry (AREA)
  • Biophysics (AREA)
  • General Health & Medical Sciences (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Zoology (AREA)
  • Theoretical Computer Science (AREA)
  • Wood Science & Technology (AREA)
  • Genetics & Genomics (AREA)
  • Evolutionary Biology (AREA)
  • Molecular Biology (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Medical Informatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • General Engineering & Computer Science (AREA)
  • Immunology (AREA)
  • Microbiology (AREA)
  • Biochemistry (AREA)
  • Botany (AREA)
  • Mycology (AREA)
  • Bioethics (AREA)
  • Databases & Information Systems (AREA)
  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)

Abstract

本发明公开了利用转录组学鉴定紫花苜蓿与根瘤菌结瘤专一性的方法，是利用转录组学手段，测定苜蓿品种接种不同根瘤菌株后的差异基因，并通过差异基因数量权重，鉴定苜蓿品种与根瘤菌株的结瘤专一性。本发明利用转录组学技术直接鉴定苜蓿品种与根瘤菌株的结瘤专一性，摆脱了对共生表型指标的依赖，避免由表型指标引起的共生效应评价不精准的问题，提高了共生效应鉴定的精准性和促进生长的高效性，结果稳定可靠，易于应用，对苜蓿与根瘤菌的精准共生具有重要的意义。

Description

利用转录组学鉴定紫花苜蓿与根瘤菌结瘤专一性的方法
技术领域
本发明属于农业领域，具体涉及利用转录组学鉴定紫花苜蓿与根瘤菌结瘤专一性的方法。
背景技术
豆科植物与根瘤菌的共生固氮体系是自然界效率最高的固氮体系。农业生产中所栽培的豆科植物品种繁多，而不同豆科植物品种与根瘤菌株间存在共生匹配专一性，即某一根瘤菌株仅能与特定豆科植物基因型高效结瘤固氮。苜蓿品种与根瘤菌株的促生效应具有明显的对苜蓿品种和根瘤菌株专一性共生的依赖性，不同根瘤菌株与苜蓿品种共生根瘤对植株生长产生的增产效应差异明显。因此，快速、精准的鉴定共生专一性，以实现根瘤菌株与苜蓿品种的高效匹配，是提高共生系统促进生长效应的重要途径。
目前，对于苜蓿-根瘤菌共生专一性效应的鉴定缺乏精准和快速简便的方法。通常采用田间苜蓿品种接种根瘤菌株，测定共生表型指标的方法筛选高效的共生组合，且在共生表型指标的选择和测定方法上缺乏统一规范和标准，存在受环境条件影响大、结果可靠性低等问题；同时，虽然占瘤率、結瘤数量、固氮酶活性、固氮量、株高和生物量等指标被大多数人所采用，但是由它们衡量的共生效应往往与实际结果存在不一致性，缺乏有效的指导意义。在分子技术应用方面，有关根瘤菌共生效应的研究，主要集中在结瘤因子与固氮酶合成及调控关键基因等方面，但尚未应用到苜蓿-根瘤菌的专一性效应方面。而转录组学技术已经被用于研究紫花苜蓿与根瘤菌的互作，不仅可以用于寻找差异基因和精确分析苜蓿品种与根瘤菌株专一性高效共生的分子机理，还可以利用转录组学技术鉴定苜蓿品种与根瘤菌株的专一性共生效应，可以从基因簇、注释蛋白和代谢通路等角度更加直观的体现苜蓿品种和根瘤菌株共生互作的有效性，为精准构建苜蓿与根瘤菌高效共生体、提高苜蓿产量和质量提供路径。
发明内容
本发明提供了一种利用转录组学鉴定紫花苜蓿与根瘤菌结瘤专一性的方法，有效的解决了上述背景技术中存在的问题。
本发明采用的技术方案如下：
利用转录组学鉴定紫花苜蓿与根瘤菌结瘤专一性的方法，包括以下步骤：
S1、挑选健康紫花苜蓿品种种子数粒，经处理后种植培养至发芽，形成紫花苜蓿幼苗；
S2、制备根瘤菌液，将根瘤菌株经培养处理后调制成根瘤菌液；
S3、接种根瘤菌株，将制备的根瘤菌液加入紫花苜蓿幼苗根部；
S4、样品采集，将收获的紫花苜蓿植株处理后随机选取若干株，采集紫花苜蓿植株的毛根处理后保存；
S5、总RNA提取，采用总RNA提取试剂盒根据Trizol法提取紫花苜蓿植株的毛根组织总RNA；
S6、对总RNA的质量检测，综合评估样品的浓度、纯度和完整性；
S7、cDNA文库构建及上机检测，以提取和检测的总RNA样品为材料，用带有Oligo(dT)的磁珠富集真核生物mRNA；
S8、转录本De novo组装和基本注释，对检测数据处理得到高质量的clean reads，并进行转录本De novo拼接组装和基本注释；
S9、对clean reads进行基因差异表达分析；
S10、通过基因差异表达分析，实现对专一性共生结瘤效应确定。
所述S4中，收获的紫花苜蓿植株处理方式和紫花苜蓿植株的毛根处理方式为：根瘤菌液接种45d后，收获紫花苜蓿植株并清洗干净，用滤纸吸干水分，随机选取，采集苜蓿植株的毛根并混匀，用锡箔纸包裹，并在泡沫盒内迅速用液氮冷冻，保存备用。
所述S4中，样品的采集部位包括毛根和根瘤。
所述S5还包括：
S501、向毛根组织和根瘤的混合物中加入液氮，迅速研磨至粉末，加入SL 和β-巯基乙醇，涡旋震荡混匀并离心处理；
S502、将上清液转至过滤柱CS中，离心处理，转移上清液至新的RNase-Free 离心管中；
S503、加入无水乙醇，混匀后转入吸附柱CR3中，离心处理，弃废液；
S504、将去蛋白液RW1加入吸附柱CR3中，离心处理，弃废液；
S505、将DNase I储存液加入新的RNase-Free离心管中，再加入RDD溶液，轻柔混匀；
S506、将DNase I工作液加入吸附柱CR3中央，室温静置；
S507、将去蛋白液RW1加入吸附柱CR3中，离心处理，弃废液；
S508、将漂洗液RW加入吸附柱CR3中，离心处理，弃废液，重复处理；
S509、将S508中的CR3吸附柱离心处理后，转移至新的RNase-Free离心管中，将RNase-Free ddH2O悬空滴入吸附膜中部，室温静置，并离心处理，收集RNA溶液；
S5010、取适量RNA进行浓度测定并进行琼脂糖电泳检测，并保存。
所述S6总RNA的质量检测包括采用Nanodrop超微量紫外分光光度计和2％琼脂糖凝胶电泳定量检测RNA条带，再利用Agilent 2100生物分析仪进行RNA 样品质量检测，综合评估样品的浓度、纯度和完整性，其余保存备用。
所述S7中，cDNA文库构建及上机检测包括：
S701、以短片段mRNA为模板，用六碱基随机引物合成第一条cDNA链；
S702、利用dNTPs、缓冲液、DNA polymerase I和RNase H合成第二条cDNA 链；
S703、通过QiaQuick PCR试剂盒纯化、EB缓冲液洗脱、末端修复、加poly (A)和连接测序接头；
S704、采用琼脂糖凝胶电泳选择长度为150～200bp的片段；
S705、采用Phusion高保真度DNA聚合酶、通用PCR引物和index(X)引物进行PCR扩增，并用Agilent 2100生物分析仪检测扩增产物质量；建好的cDNA 测序文库用IlluminaNovaSeq 6000平台进行测序。
所述S8中，转录本De novo组装和基本注释包括：
S801、利用FastQC(v0.11.5)对测序得到的原始数据进行质控，滤去含 adapter、N比例>10％以及低质量的reads，获得高质量的clean reads；
S802、通过碱基组成、质量值分布图和饱和度分析检测数据质量情况；
S803、使用Trinity(v2.2.0)组装软件对过滤后的高质量clean reads进行转录本de novo拼接组装，后进行转录本聚类，取最长转录本为Unigene，最后通过N50数值和序列长度来评估组装结果质量；
S804、采用Blast+(v2.4.0)将转录本注释到Nr、COG/KOG和KEGG数据库；
S805、利用Blast2Go(v2.3.5)软件进行GO注释。
所述S9中，基因差异表达分析包括采用EdgeR软件(v3.14.0)进行差异基因分析，基因的表达量用FPKM值表示。错误发现率≤0.05且|log2FC|>1为筛选差异基因的条件，以表达差异倍数≥2或≤0.5表示表达上调或下调，在转录组比较组合中选取差异基因；通过比对GO数据库对差异基因进行功能注释和显著性富集分析；通过比对KEGG数据库，对差异基因参与的细胞代谢途径及其产物功能进行系统分析。
所述S10中，专一性共生结瘤效应确定方式为：多个根瘤菌株接种于同一个苜蓿品种，对比每两两组合间的差异基因数量，小于或等于苜蓿总基因组的 0.5％时，确定为此菌株与该苜蓿品种专一性弱，表明此根瘤菌株对该苜蓿品种促进生长能力弱；差异基因数量大于或等于苜蓿总基因组的12.2％时，确定为专一性强，表明此根瘤菌株对该苜蓿品种促进生长能力强。
本发明的有益效果为：
1)利用转录组学鉴定专一性共生结瘤效应，摆脱了对共生表型指标的依赖，避免了受环境条件影响而导致结果可靠性低的问题。
2)转录组学技术成熟，操作规范，标准化程度高，结果稳定可靠，易于普及应用。
具体实施方式
下面通过具体的实施例进一步的说明本发明的技术方案：
实施例1
利用转录组学鉴定紫花苜蓿与根瘤菌结瘤专一性的方法，包括以下步骤：
1)苜蓿品种：甘农3号紫花苜蓿(M.sativa cv.Gannong No.3)、甘农9号紫花苜蓿(M.sativa cv.Gannong No.9)、陇中紫花苜蓿(M.sativa cv.Longzhong)和清水紫花苜蓿(M.sativa cv.Qingshui)。
2)根瘤菌株：6株中华根瘤菌(S.meliloti)：LL2、WLP2、QL2、LL1、 G3L3和LP3。
3)培育幼苗：挑选子粒饱满、健康的紫花苜蓿品种种子数粒，在生物安全柜内用碘伏溶液(有效碘含量0.45％～0.55％(W·V-1))震荡灭菌3min，用无菌水冲洗5～6次，于清水琼脂中28℃催芽48h。将河沙洗净，过筛(2-mm)， pH调至中性后105℃烘干，121℃灭菌6h，冷却。在直径13.2cm、高10cm 的塑料盆中装入河沙450g，并均匀植入60粒种子发芽，深度为2cm。之后将塑料小盆放入29cm×20cm×9.5cm的塑料面盆中，每个面盆内放2个小盆，并浇灌500mL灭菌水。培养条件为：光照时间16h·d-1，有光照时温度21～ 25℃，无光照时温度16～20℃，相对湿度45±5％。在第7d浇灌500mL Hoagland 有氮营养液。
4)制备菌液：将步骤2)中的根瘤菌株TY平板活化后转接入50mL TY液体培养基，28℃、180rpm·min-1振荡培养至菌液光密度OD600nm值为0.5～1， 4000rpm·min-1离心5min，抛去上清液后用等体积的无菌水冲洗菌体，并用涡旋振荡器将其打散，最后用无菌水调制成OD600nm为0.5的根瘤菌液。
5)接种根瘤菌株：在幼苗生长至第15d左右时(第1片真叶出现)，用移液器将制备好的菌液加入幼苗根部，每小盆60mL，以不含根瘤菌的TY液体培养基为对照。根瘤菌接种苜蓿品种的组合包括：菌株LL2+甘农9号苜蓿、菌株WLP2+甘农9号苜蓿、菌株LL2+甘农3号苜蓿、菌株QL2+甘农3号苜蓿、菌株WLP2+清水苜蓿、菌株LL1+清水苜蓿、菌株G3L3+陇中苜蓿、菌株LP3+陇中苜蓿。每个接种处理4个小盆作为重复。接种后每隔7d在塑料面盆中浇灌Hoagland无氮营养液500mL。其余时间用灭菌蒸馏水补充水分。
6)样品采集：接种45d后(植株已生长60d)，收获紫花苜蓿植株并清洗干净，用滤纸吸干水分。每盆随机选取20株，采集苜蓿植株的毛根(包括根瘤)并混匀，用锡箔纸包裹，并在泡沫盒内迅速用液氮冷冻，–80℃冰箱保存备用。
7)总RNA提取：采用总RNA提取试剂盒根据Trizol法提取根组织总RNA。 (i)向100mg毛根组织和根瘤的混合物中加入液氮，迅速研磨至粉末，加入 SL 475μL和β-巯基乙醇25μL，涡旋震荡混匀，12000rpm离心2min；(ii) 将上清液转至过滤柱CS中，12000rpm离心2min，转移上清液至新的RNase-Free 离心管中；(i i i)加入无水乙醇(0.4倍上清体积)，混匀后转入吸附柱CR3 中，12000rpm离心15sec，弃废液；(iv)将350μL去蛋白液RW1加入吸附柱CR3中，12000rpm离心15sec，弃废液；(v)将10μLDNase I储存液加入新的RNase-Free离心管中，再加入RDD溶液70μL，轻柔混匀；(vi)将 80μLDNase I工作液加入吸附柱CR3中央，室温静置15min；(vii)将350μ L去蛋白液RW1加入吸附柱CR3中，12000rpm离心15sec，弃废液；(viii) 将500μL漂洗液RW加入吸附柱CR3中，12000rpm离心15sec，弃废液；(ix) 重复步骤(viii)；(x)将步骤(ix)中的CR3吸附柱12000rpm离心2min，转移至新的RNase-Free离心管中，将30～50μLRNase-Free ddH2O悬空滴入吸附膜中部，室温静置2min，12000rpm离心1min，收集RNA溶液；(xi)取适量RNA进行浓度测定并进行琼脂糖电泳检测，其余保存于–80℃冰箱。
8)总RNA的质量检测：采用Nanodrop超微量紫外分光光度计和2％琼脂糖凝胶电泳定量检测RNA条带，再利用Agilent 2100生物分析仪进行RNA样品质量检测，综合评估样品的浓度、纯度和完整性，其余–80℃冰箱保存备用。
9)cDNA文库构建及上机检测：以步骤7)、8)提取和检测的总RNA样品为材料，用带有Oligo(dT)的磁珠富集真核生物mRNA。首先，以短片段mRNA 为模板，用六碱基随机引物(Random hexamers)合成第一条cDNA链；其次，利用dNTPs、缓冲液、DNA polymerase I和RNase H合成第二条cDNA链，通过QiaQuick PCR试剂盒纯化、EB缓冲液洗脱、末端修复、加poly(A)和连接测序接头，采用琼脂糖凝胶电泳选择长度为150～200bp的片段。最后采用Phusion 高保真度DNA聚合酶、通用PCR引物和index(X)引物进行PCR扩增，并用Agilent2100生物分析仪检测扩增产物质量。建好的cDNA测序文库用Illumina NovaSeq 6000平台进行测序。
10)转录本De novo组装和基本注释：利用FastQC(v0.11.5)对测序得到的原始数据进行质控，滤去含adapter、N比例>10％以及低质量的reads，获得高质量的clean reads；并通过碱基组成、质量值分布图和饱和度分析检测数据质量情况。使用Trinity(v2.2.0)组装软件对过滤后的高质量clean reads进行转录本de novo拼接组装，然后进行转录本聚类，取最长转录本为Unigene，最后通过N50数值和序列长度来评估组装结果质量。采用Blast+(v2.4.0)将转录本注释到Nr(Non-redundant protein sequences from NCBI)、COG/KOG (Cluster of Orthologous Groups of proteins)和KEGG(Kyoto encyclopediaof genes and genomes)数据库；利用Blast2Go(v2.3.5)软件进行GO(Gene ontology)注释。
11)基因差异表达分析：利用RSEM软件包(v1.2.31)将全长reads直接比对到参考unigenes上，用FPKM值表示基因表达水平，对基因的表达量进行定量分析。采用EdgeR软件(v3.14.0)进行差异基因分析，错误发现率(False discover rate，FDR)≤0.05且|log2FC|>1为筛选差异基因的条件，以表达差异倍数(Fold change，FC)≥2或≤0.5表示表达上调或下调，在转录组比较组合中选取差异基因。FDR值越小，表达差异越显著。通过比对GO数据库对差异基因进行功能注释和显著性富集分析；通过比对KEGG数据库，对差异基因参与的细胞代谢途径及其产物功能进行系统分析。
12)专一性共生结瘤效应确定：陇中苜蓿上接种菌株组合G3L3 vs.LP3的差异基因数量占苜蓿总基因组的12.2％，这2株菌与陇中苜蓿专一性强；甘农3 号苜蓿上接种菌株组合LL2 vs.QL2的差异基因数量占总基因组的1.7％，其中 QL2仅与甘农3号苜蓿专一性高效共生，LL2与甘农3号和甘农9号均高效共生；清水苜蓿上接种菌株组合WLP2 vs.LL1的差异基因数量占总基因组的2.7％，其中LL1仅与清水苜蓿专一性高效共生，WLP2与清水和甘农9号苜蓿均高效共生；甘农9号苜蓿上接种菌株组合LL2 vs.WLP2的差异基因数量占总基因组的0.5％，这2株菌与甘农9号苜蓿和其它苜蓿品种均共生，专一性弱。
本发明采用转录组学方法，以4个M.sativa品种和6株S.meliloti菌为基础，直接鉴定紫花苜蓿与根瘤菌的结瘤专一性。根瘤菌接种处理对比组合间的差异基因越多，此菌株与该品种的专一性越强。转录组结果与共生效应指标分析结果完全一致，并为结瘤专一性效应提供了合理的解释。
此发明直接采用转录组学鉴定紫花苜蓿与根瘤菌的结瘤专一性，是从植物品种和同一个种内根瘤菌株的角度对专一性共生效应进行的更深层次的分析，极大程度上简化了高效共生体系的发掘和专一性根瘤菌株的鉴定过程。
由于所采用的试验和分析方法简便合理，易于操作。因此本发明具有很好的实用性和可靠性。
对于本领域技术人员而言，显然本发明不限于上述示范性实施例的细节，而且在不背离本发明的精神或基本特征的情况下，能够以其他的具体形式实现本发明。因此，无论从哪一点来看，均应将实施例看作是示范性的，而且是非限制性的，本发明的范围由所附权利要求而不是上述说明限定，因此旨在将落在权利要求的等同要件的含义和范围内的所有变化囊括在本发明内。
此外，应当理解，虽然本说明书按照实施方式加以描述，但并非每个实施方式仅包含一个独立的技术方案，说明书的这种叙述方式仅仅是为清楚起见，本领域技术人员应当将说明书作为一个整体，各实施例中的技术方案也可以经适当组合，形成本领域技术人员可以理解的其他实施方式。

Claims (9)

1.利用转录组学鉴定紫花苜蓿与根瘤菌结瘤专一性的方法，其特征在于，包括以下步骤：
S1、挑选健康紫花苜蓿品种种子数粒，经处理后种植培养至发芽，形成紫花苜蓿幼苗；
S2、制备根瘤菌液，将根瘤菌株经培养处理后调制成根瘤菌液；
S3、接种根瘤菌株，将制备的根瘤菌液加入紫花苜蓿幼苗根部；
S4、样品采集，将收获的紫花苜蓿植株处理后随机选取若干株，采集紫花苜蓿植株的毛根处理后保存；
S5、总RNA提取，采用总RNA提取试剂盒根据Trizol法提取紫花苜蓿植株的毛根组织总RNA；
S6、对总RNA的质量检测，综合评估样品的浓度、纯度和完整性；
S7、cDNA文库构建及上机检测，以提取和检测的总RNA样品为材料，用带有Oligo(dT)的磁珠富集真核生物mRNA；
S8、转录本De novo组装和基本注释，对检测数据处理得到高质量的clean reads，并进行转录本De novo拼接组装和基本注释；
S9、对clean reads进行基因差异表达分析；
S10、通过基因差异表达分析，实现对专一性共生结瘤效应确定。
2.根据权利要求1所述的利用转录组学鉴定紫花苜蓿与根瘤菌结瘤专一性的方法，其特征在于：所述S4中，收获的紫花苜蓿植株处理方式和紫花苜蓿植株的毛根处理方式为：根瘤菌液接种45d后，收获紫花苜蓿植株并清洗干净，用滤纸吸干水分，随机选取，采集苜蓿植株的毛根并混匀，用锡箔纸包裹，并在泡沫盒内迅速用液氮冷冻，保存备用。
3.根据权利要求1所述的利用转录组学鉴定紫花苜蓿与根瘤菌结瘤专一性的方法，其特征在于：所述S4中，样品的采集部位包括毛根和根瘤。
4.根据权利要求1所述的利用转录组学鉴定紫花苜蓿与根瘤菌结瘤专一性的方法，其特征在于：所述S5还包括：
S501、向毛根组织和根瘤的混合物中加入液氮，迅速研磨至粉末，加入SL和β-巯基乙醇，涡旋震荡混匀并离心处理；
S502、将上清液转至过滤柱CS中，离心处理，转移上清液至新的RNase-Free离心管中；
S503、加入无水乙醇，混匀后转入吸附柱CR3中，离心处理，弃废液；
S504、将去蛋白液RW1加入吸附柱CR3中，离心处理，弃废液；
S505、将DNase I储存液加入新的RNase-Free离心管中，再加入RDD溶液，轻柔混匀；
S506、将DNase I工作液加入吸附柱CR3中央，室温静置；
S507、将去蛋白液RW1加入吸附柱CR3中，离心处理，弃废液；
S508、将漂洗液RW加入吸附柱CR3中，离心处理，弃废液，重复处理；
S509、将S508中的CR3吸附柱离心处理后，转移至新的RNase-Free离心管中，将RNase-Free ddH2O悬空滴入吸附膜中部，室温静置，并离心处理，收集RNA溶液；
S5010、取适量RNA进行浓度测定并进行琼脂糖电泳检测，并保存。
5.根据权利要求1所述的利用转录组学鉴定紫花苜蓿与根瘤菌结瘤专一性的方法，其特征在于：所述S6总RNA的质量检测包括采用Nanodrop超微量紫外分光光度计和2％琼脂糖凝胶电泳定量检测RNA条带，再利用Agilent 2100生物分析仪进行RNA样品质量检测，综合评估样品的浓度、纯度和完整性，其余保存备用。
6.根据权利要求1所述的利用转录组学鉴定紫花苜蓿与根瘤菌结瘤专一性的方法，其特征在于：所述S7中，cDNA文库构建及上机检测包括：
S701、以短片段mRNA为模板，用六碱基随机引物合成第一条cDNA链；
S702、利用dNTPs、缓冲液、DNA polymerase I和RNase H合成第二条cDNA链；
S703、通过QiaQuick PCR试剂盒纯化、EB缓冲液洗脱、末端修复、加poly(A)和连接测序接头；
S704、采用琼脂糖凝胶电泳选择长度为150～200bp的片段；
S705、采用Phusion高保真度DNA聚合酶、通用PCR引物和index(X)引物进行PCR扩增，并用Agilent 2100生物分析仪检测扩增产物质量；建好的cDNA测序文库用Illumina NovaSeq6000平台进行测序。
7.根据权利要求1所述的利用转录组学鉴定紫花苜蓿与根瘤菌结瘤专一性的方法，其特征在于：所述S8中，转录本De novo组装和基本注释包括：
S801、利用FastQC(v0.11.5)对测序得到的原始数据进行质控，滤去含adapter、N比例>10％以及低质量的reads，获得高质量的clean reads；
S802、通过碱基组成、质量值分布图和饱和度分析检测数据质量情况；
S803、使用Trinity(v2.2.0)组装软件对过滤后的高质量clean reads进行转录本denovo拼接组装，后进行转录本聚类，取最长转录本为Unigene，最后通过N50数值和序列长度来评估组装结果质量；
S804、采用Blast+(v2.4.0)将转录本注释到Nr、COG/KOG和KEGG数据库；
S805、利用Blast2Go(v2.3.5)软件进行GO注释。
8.根据权利要求1所述的利用转录组学鉴定紫花苜蓿与根瘤菌结瘤专一性的方法，其特征在于：所述S9中，基因差异表达分析包括采用EdgeR软件(v3.14.0)进行差异基因分析，基因的表达量用FPKM值表示。错误发现率≤0.05且|log2FC|>1为筛选差异基因的条件，以表达差异倍数≥2或≤0.5表示表达上调或下调，在转录组比较组合中选取差异基因；通过比对GO数据库对差异基因进行功能注释和显著性富集分析；通过比对KEGG数据库，对差异基因参与的细胞代谢途径及其产物功能进行系统分析。
9.根据权利要求1所述的利用转录组学鉴定紫花苜蓿与根瘤菌结瘤专一性的方法，其特征在于：所述S10中，专一性共生结瘤效应确定方式为：多个根瘤菌株接种于同一个苜蓿品种，对比每两两组合间的差异基因数量，小于或等于苜蓿总基因组的0.5％时，确定为此菌株与该苜蓿品种专一性弱，表明此根瘤菌株对该苜蓿品种促进生长能力弱；差异基因数量大于或等于苜蓿总基因组的12.2％时，确定为专一性强，表明此根瘤菌株对该苜蓿品种促进生长能力强。
CN202010130201.9A 2020-02-28 2020-02-28 利用转录组学鉴定紫花苜蓿与根瘤菌结瘤专一性的方法 Withdrawn CN111647676A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010130201.9A CN111647676A (zh) 2020-02-28 2020-02-28 利用转录组学鉴定紫花苜蓿与根瘤菌结瘤专一性的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010130201.9A CN111647676A (zh) 2020-02-28 2020-02-28 利用转录组学鉴定紫花苜蓿与根瘤菌结瘤专一性的方法

Publications (1)

Publication Number Publication Date
CN111647676A true CN111647676A (zh) 2020-09-11

Family

ID=72344711

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010130201.9A Withdrawn CN111647676A (zh) 2020-02-28 2020-02-28 利用转录组学鉴定紫花苜蓿与根瘤菌结瘤专一性的方法

Country Status (1)

Country Link
CN (1) CN111647676A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113517024A (zh) * 2021-04-25 2021-10-19 北京果壳生物科技有限公司 一种基于ONT全长转录组测序数据的denovo分析方法
CN114317753A (zh) * 2021-12-30 2022-04-12 北京迈基诺基因科技股份有限公司 眼肿瘤融合基因的检测模型及构建方法和检测方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109593655A (zh) * 2019-01-02 2019-04-09 甘肃农业大学 苜蓿品种专一性共生根瘤菌的筛选方法
CN109734785A (zh) * 2019-03-13 2019-05-10 深圳市伊山万瑞科技有限公司 一种来源于苦参的蛋白质及其编码基因与应用
CN109943496A (zh) * 2019-01-02 2019-06-28 甘肃农业大学 一株根瘤菌寻配专一性高效共生苜蓿品种的方法
CN111066593A (zh) * 2020-01-17 2020-04-28 华中农业大学 酸性磷酸酶在苜蓿根瘤和菌根共生中的功能研究方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109593655A (zh) * 2019-01-02 2019-04-09 甘肃农业大学 苜蓿品种专一性共生根瘤菌的筛选方法
CN109943496A (zh) * 2019-01-02 2019-06-28 甘肃农业大学 一株根瘤菌寻配专一性高效共生苜蓿品种的方法
CN109734785A (zh) * 2019-03-13 2019-05-10 深圳市伊山万瑞科技有限公司 一种来源于苦参的蛋白质及其编码基因与应用
CN111066593A (zh) * 2020-01-17 2020-04-28 华中农业大学 酸性磷酸酶在苜蓿根瘤和菌根共生中的功能研究方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
康文娟: "花苜蓿根瘤菌生物型划分及其转录组学分析", 《中国博士学位论文全文数据库（农业科技辑）》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113517024A (zh) * 2021-04-25 2021-10-19 北京果壳生物科技有限公司 一种基于ONT全长转录组测序数据的denovo分析方法
CN114317753A (zh) * 2021-12-30 2022-04-12 北京迈基诺基因科技股份有限公司 眼肿瘤融合基因的检测模型及构建方法和检测方法

Similar Documents

Publication Publication Date Title
BR112014029988B1 (pt) Método para produzir uma planta aprimorada e para seleção genômica de associações de planta-simbionte
CN106755478A (zh) 一种鉴定嘎拉苹果后代植株的ssr分子标记iv及其应用
Moura et al. Validation of reference genes from Eucalyptus spp. under different stress conditions
CN111647676A (zh) 利用转录组学鉴定紫花苜蓿与根瘤菌结瘤专一性的方法
CN112981000B (zh) 与棉花抗旱性相关的InDel分子标记及其应用
CN115700282B (zh) 一种玉米抗旱基因Zm00001d030678及其应用
CN112063738A (zh) 与毛白杨生长相关基因的分子标记、筛选方法、试剂盒及应用
Wong et al. RNA-seq data of Ganoderma boninense at axenic culture condition and under in planta pathogen-oil palm (Elaeis guineensis Jacq.) interaction
CN106755018A (zh) 利用棉花GhPEPC基因创造棉花高油材料的方法
CN103937790A (zh) 一种与油菜种子硫苷含量紧密相关的分子标记及应用
CN108823242A (zh) 一种转化玉米cDNA文库改良水稻品种的方法
CN109599147B (zh) 镉胁迫条件下毛白杨中的一种差异表达的NACs转录因子及其获取方法
CN115896128B (zh) 一种烟草硝酸盐转运蛋白NtNPF6.13、编码基因及其应用
CN114410816B (zh) 一种适用木薯抗病研究的内参基因的筛选方法及应用
CN102559668B (zh) 与玉米抗丝黑穗病基因连锁的SNP位点、基于该位点的分子标记LSdCAP3及其应用
CN107868840A (zh) 一种亚麻中与全生长日数相关联的ssr分子标记和应用
CN114015756A (zh) 凤丹耐涝差异表达基因及基于转录组测序挖掘的方法
CN111549042A (zh) 核酸分子在植物转基因、分子育种、病害防治和分子标记物中的用途
CN110885817A (zh) 一种田菁线粒体及其dna提取的有效方法
CN108546704A (zh) 一种无功能的OsC1基因及其应用
CN115852031B (zh) 与甜瓜种子大小性状共分离的snp标记及其应用
CN116064919B (zh) 甘蔗割手密种不同组织荧光定量内参基因的筛选和应用
CN115820915B (zh) 甘蔗割手密种非生物胁迫条件下的荧光定量内参基因的筛选和应用
CN113373254B (zh) 一组与陆地棉耐盐qtl连锁的分子标记及应用
CN115961081B (zh) 与玉米抗茎基腐病基因位点qFCR9紧密连锁的分子标记及其应用

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
WW01 Invention patent application withdrawn after publication
WW01 Invention patent application withdrawn after publication

Application publication date: 20200911