CN107815502B - 一种鉴定菊花耐涝性的dCAPS标记开发及应用 - Google Patents

一种鉴定菊花耐涝性的dCAPS标记开发及应用 Download PDF

Info

Publication number
CN107815502B
CN107815502B CN201711176758.0A CN201711176758A CN107815502B CN 107815502 B CN107815502 B CN 107815502B CN 201711176758 A CN201711176758 A CN 201711176758A CN 107815502 B CN107815502 B CN 107815502B
Authority
CN
China
Prior art keywords
chrysanthemum
waterlogging
marker
dcaps
tolerance
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.)
Active
Application number
CN201711176758.0A
Other languages
English (en)
Other versions
CN107815502A (zh
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.)
Nanjing Agricultural University
Original Assignee
Nanjing 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 Nanjing Agricultural University filed Critical Nanjing Agricultural University
Priority to CN201711176758.0A priority Critical patent/CN107815502B/zh
Publication of CN107815502A publication Critical patent/CN107815502A/zh
Application granted granted Critical
Publication of CN107815502B publication Critical patent/CN107815502B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • 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/156—Polymorphic or mutational markers
    • 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/172—Haplotypes

Landscapes

  • Chemical & Material Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Engineering & Computer Science (AREA)
  • Organic Chemistry (AREA)
  • Proteomics, Peptides & Aminoacids (AREA)
  • Health & Medical Sciences (AREA)
  • Biotechnology (AREA)
  • Zoology (AREA)
  • Wood Science & Technology (AREA)
  • Immunology (AREA)
  • Mycology (AREA)
  • Microbiology (AREA)
  • Molecular Biology (AREA)
  • Botany (AREA)
  • Biophysics (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

本发明涉及一种鉴定菊花耐涝性的dCAPS标记开发方法及应用，属于生物技术领域，该方法包括以下几个步骤：a、通过GWAS方法挖掘与菊花耐涝性显著关联的SNP位点；b、SNP突变位点特异性酶切位点分析及dCAPS引物设计；c、结合基因型与表型数据预期酶切扩增多态性；d、dCAPS标记在自然耐涝极端群体的验证；e、在F1后代耐涝极端群体中对WT‑dCAPS1标记进行进一步验证。本发明开发了一个与菊花耐涝性状共分离的dCAPS共显性标记，命名为WT‑dCAPS1，在两个群体的平均鉴定准确率为78.9％，初步认为可以应用于菊花耐涝性分子标记辅助育种，大大缩短育种周期，对于培育耐涝性菊花新品种具有重要的理论和实践意义。

Description

一种鉴定菊花耐涝性的dCAPS标记开发及应用
技术领域
本发明属于生物技术领域，涉及一种鉴定菊花耐涝性的dCAPS标记开发方法及应用，可用于菊花耐涝性分子标记辅助育种。
背景技术
涝渍是危害世界菊花生产的严重自然灾害，培育耐涝性菊花新品种是解决这一问题的直接有效方法。因此，开发与菊花耐涝性相关的特异分子标记具有重要意义。目前，RAPD、AFLP、SRAP、ISSR、SSR等传统的分子标记已经成功应用在菊花种质资源鉴定、遗传多样性分析、连锁图谱构建和关联分析中[1-5]。随着现代分子生物学的发展及基因组测序计划的进行，SNP标记应运而生。SNP(Singie Nucieotide Poiymorphism)，即单核苷酸多态性，是指基因组水平上由单核苷酸的变异所引起的DNA序列多态性，主要包括单个碱基的插入、缺失、转换和颠倒等。SNP在基因组中分布广泛、数量多，遗传稳定性高、遗传分析重现性好且准确性高，易于快速且高通量地进行基因型分型[6]。基于以上优点，SNP已经被广泛地应用于多个物种上的全基因组关联分析(genome wide association study，GWAS)中[7-9]。在菊花上，Chong等人[10]对199份菊花资源利用简化基因组测序SLAF(specific-locusamplified fragment sequencing)技术，检测到＞46万个SNP，并据此分析了菊花的进化关系。
酶切扩增多态性序列(cleaved amplified polymorphic sequence，CAPS)标记是利用特异引物PCR与限制性酶切相结合而产生的一种检测SNP位点的DNA标记[11]。其原理是SNP突变位点位于某种内切酶的识别序列上，经过特异引物PCR扩增后，用相应内切酶对PCR产物进行酶切，从而产生酶切片段的多态性，这样通过琼脂糖或聚丙烯酰胺凝胶电泳、染色即可显示。但是SNP恰好位于某种限制性酶切位点的这种情况还比较少，于是Michaels和Amasino[12]及Neff[13]在CAPS标记基础上，开发了人工引入错配碱基的衍生酶切扩增多态性标记(derived cleaved amplified polymorphic sequence，dCAPS)，从而几乎能够将所有SNP位点转换成可以电泳检测的标记。无论是CAPS还是dCAPS，均具有共显性、位点特异性、操作简单、检测快速、成本低廉且不依赖于精密仪器设备等特点，可用于用于植物基因分型、定位、遗传多样性分析、品种鉴定等方面。
菊花属于高度杂合的异源多倍体，复杂的遗传背景限制了其基因组研究的发展。目前应用在菊花上的分子标记大多为显性标记，且较难直接应用在菊花育种工作中。此前，本发明人已经公开了一种采用盆栽模拟淹水法模拟菊花涝害处理以及利用隶属函数法评价菊花耐涝性的专利(106446596A)[14]；此外本发明人在2016年公开了199份菊花资源的SNP标记[11]，本专利在以上基础上，首先利用GWAS检测到与菊花耐涝性显著关联的SNP，据此成功开发了一个共显性dCAPS标记，并在具有极端耐涝性差异的两个群体进行验证，结果表明该标记可以用于指导菊花耐涝性分子标记辅助育种。
[1]Zhang F,Chen S,Chen F,Fang W,Li F.Apreliminary genetic linkage mapof chrysanthemum(Chrysanthemum morifolium)cultivars using RAPD,ISSR and AFLPmarkers.Scientia Horticulturae,2010,125(3):422-428
[2]Zhang F,Chen S,Chen F,Fang W,Chen Y,Li F.SRAP-based mapping andQTL detection for inflorescence-related traits in chrysanthemum(Dendranthemamorifolium).Molecular Breeding,2011,27(1):11-23
[3]Wang C,Zhang F,Guan Z,Chen S,Jiang J,Fang W,Chen F.Inheritance andmolecular markers for aphid(Macrosiphoniella sanbourni)resistance inchrysanthemum(Chrysanthemum morifolium Ramat.).Scientia Horticulturae,2014,180:220-226
[4]Peng H,Zhang F,Jiang J,Chen S,Fang W,Guan Z,Chen F.Identificationof quantitative trait loci for branching traits of spray cutchrysanthemum.Euphytica,2015,202(3):385-392
[5]Li P,Zhang F,Chen S,Jiang J,Wang H,Su J,Fang W,Guan Z,ChenF.Genetic diversity,population structure and association analysis in cutchrysanthemum(Chrysanthemum morifolium Ramat.).Mol Genet Genomics,2016,291(3):1117–1125
[6]唐立群，肖层林，王伟平.SNP分子标记的研究及其应用进展[J].中国农学通报，2012，28(12)：154-158
[7]陈广凤，陈建省，田纪春.小麦株高相关性状与SNP标记全基因组关联分析[J].作物学报，2015，41(10)：1500-1509
[8]Nimmakayala P,Abburi VL,Saminathan T,Almeida A,Davenport B,Davidson J,C.V.Chandra MR,Gerald H,Andreas E,Doil C,StommelJ,Umesh KR.Genome-wide divergence and linkage disequilibrium analyses for Capsicum baccatumrevealed by genome-anchored single nucleotide polymorphisms.
Frontiers in plant science,2016,7.doi:10.3389/fpls.2016.01646
[9]Zhou Q,Zhou C,Zheng W,Mason AS,FanS,Wu C,Fu D,Huang Y.Genome-wideSNP markers based on SLAF-seq uncover breeding traces in rapeseed(Brassicanapus L.).Frontiers in plant science,2017,8.doi:10.3389/fpls.2017.00648
[10]Chong X,Zhang F,Wu Y,Yang X,Zhao N,Wang H,Guan Z,Fang W,Chen F.ASNP-enabled assessment of genetic diversity,evolutionary relationships andthe identification of candidate genes in chrysanthemum.Genome biology andevolution,2016,8(12):3661-3671
[11]Lee GA,Koh HJ,Chung HK,Dixit A,Chung JW,Ma KH,Lee SY,Lee JR,LeeGS,Gwag JG.Development of SNP-based CAPS and dCAPS markers in eight differentgenes involved in starch biosynthesis in rice.Molecular Breeding,2009,24(1):93–101.
[12]Michaels SD,Amasino RM.A robust method for detecting single-nucleotide changes as polymorphic markers by PCR.Plant Journal,1998,14:381–385
[13]Neff MM,Neff JD,Chory J,Pepper AE.dCAPS,a simpletechnique for thegenetic analysis of single nucleotide polymorphisms:experimental applicationsin Arabidopsis thaliana genetics.Plant Journal,1998,14:387–392
[14]陈发棣，苏江硕，张飞，王海滨，管志勇，房伟民，廖园.一种与菊花耐涝性显著相关的分子标记及其鉴定方法和应用[P].中国专利：106446596A，2017-02-22(审中—实审)
发明内容
为了克服现有技术的不足，本发明通过全基因组关联分析(GWAS)方法筛选出与菊花耐涝性显著关联的SNP位点，据此开发了一个快速、简单易行的共显性dCAPS标记，命名为WT-dCAPS1。我们在两个群体的极端耐涝和不耐涝个体中对该标记进行验证，平均准确率为78.9％，初步认为可用于菊花耐涝性分子标记辅助选择育种，大大缩短育种周期，从而提高育种效率。
本发明的目的通过以下技术手段获得：
一种鉴定菊花耐涝性的dCAPS标记的筛选方法，包括以下步骤：
a、通过GWAS方法挖掘与菊花耐涝性显著关联的SNP位点
选取多份来源不同且无直接亲缘关系的菊花品种，采用盆栽模拟淹水法(106446596A)对该群体进行了三次耐涝性鉴定，获得表型数据；
基于已有的SNP标记数据分别采用ADMIXTURE软件和SPAGeDi软件对该群体进行分析，获得群体结构Q矩阵和亲缘关系系数K矩阵；在TASSEL 5.0软件中，根据三次耐涝鉴定试验表型数据的平均值，利用一般线性模型(general linear model，GLM)和(mixed linearmodel，MLM)两种分析模型进行全基因关联分析(GWAS)，以两种模型均能检测到的显著SNP位点(P<1E-3)为与菊花耐涝性紧密关联的位点，同时得出该标记的贡献率R2；
b、SNP突变位点特异性酶切位点分析及dCAPS引物设计
(1)从步骤a所得的与菊花耐涝性紧密关联的位点中，选择P值最小，且贡献率较大的SNP位点进行后续dCAPS标记开发；
(2)利用在线酶切识别软件dCAPS Finder 2.0(http://helix.wustl.edu/dcaps/dcaps.html)查找步骤(1)中所选择的SNP位点突变所引起的限制性内切酶信息，选用合适的内切酶，且人工引入错配碱基；
(3)利用Primer Primer 5.0软件设计PCR引物，引物设计标准为18-25bp；GC含量40-60％；退火温度为55℃-65℃，且上下游引物的GC含量和退火温度尽量保持接近；无引物二聚体及发卡结构；
c、结合基因型与表型数据预期酶切扩增多态性
对步骤b中选择的SNP位点在该群体进行基因型分析，对选择的SNP位点的基因型分别计算平均表型值，并用Student’s t检验方法分析耐涝差异显著性，结合基因型与表型值预期酶切扩增多态性；
d、dCAPS标记在自然耐涝极端群体的验证
(1)提取多份耐涝菊花品种材料和多份不耐涝菊花品种材料的基因组DNA，用上述步骤b中设计的特异引物对耐涝菊花品种材料和不耐涝菊花品种材料进行PCR扩增，获得目的片段全长；
(2)将扩增产物用步骤b中选用的内切酶进行酶切，并对酶切产物进行电泳检测、银染法染色；酶切扩增多态性分析，如果检测结果与步骤c中的预期相符，则将该SNP位点作为鉴定菊花耐涝性的dCAPS标记，将该标记命名为WT-dCAPS1标记；
e、在F1后代耐涝极端群体中对WT-dCAPS1标记进行进一步验证
选取耐涝性菊花品种和不耐涝性菊花品种进行杂交，在F1杂交群体中选择极端耐涝株系和极端不耐涝株系进一步验证所筛选的WT-dCAPS1标记。
步骤b的(1)中选择的P值最小且贡献率较大的SNP位点为标签Marker6619的75位置处的SNP位点(Marker6619-75)。
步骤b的(2)中选用识别序列为GCTAGC的内切酶NheI，在Marker6619-74处引入错配碱基(T→C)，步骤b的(3)中设计的引物序列分别为：F 5’-ATGCACAAGGCATTGGTCTTGC-3’；R 5’-CCTCAAAGTTTTGGCTATTCTCC-3’。
步骤c中对Marker6619-75位点进行基因型分析，共有C和Y两种基因型，Y代表CT兼并碱基，根据内切酶NheI的识别序列可知，若该SNP位点为C，不能被酶切，若为Y，则可以部分酶切；对Marker6619-75位点的两种基因型分别计算平均表型值，并用Student’s t检验方法分析耐涝差异显著性，显示C基因型的材料耐涝性显著高于Y基因型(P<0.01)，即不耐涝的菊花品种的PCR产物更有可能被酶切。
步骤d的(2)中将扩增产物用内切酶NheI进行酶切，耐涝菊花品种材料的扩增产物不能酶切；不耐涝菊花品种材料的扩增产物酶切产物电泳检测结果包含两条主带170bp和149bp(部分酶切的产物包括170bp，149bp和21bp三个片段，其中21bp的条带太小，电泳检测时不能被检测到)。
上述方法筛选的鉴定菊花耐涝性WT-dCAPS标记。
上述的鉴定菊花耐涝性WT-dCAPS标记在耐涝性菊花育种中的应用。
一种菊花耐涝性dCAPS分子标记辅助选择育种方法，该dCAPS分子标记为Marker6619-75，选用识别序列为GCTAGC的内切酶NheI，设计上下游引物序列：F 5’-ATGCACAAGGCATTGGTCTTGC-3’和R 5’-CCTCAAAGTTTTGGCTATTCTCC-3’，以菊花的基因组DNA为模板，以所设计的上下游引物进行PCR扩增，获得目的片段全长，共170bp；将扩增产物用内切酶NheI进行酶切，如果酶切产物只包含一条170bp的主带，则待测菊花为耐涝菊花品种；若酶切产物的电泳检测结果包含两条主带170bp和149bp，则待测菊花为不耐涝菊花品种。
本发明的有益效果：
本发明建立了一种鉴定菊花耐涝性的dCAPS标记开发体系，首先使用全基因组关联分析方法挖掘和菊花耐涝性显著关联的SNP位点，然后将获得的SNP进一步转化成1个dCAPS标记，并在自然群体和F1后代群体中的极端耐涝和不耐涝个体中进行验证。与目前技术相比，其优点是：
(1)使用在基因组中分布广泛、数量多，遗传稳定性高且遗传分析重现性好、准确性高的SNP标记进行GWAS分析。此外，我们进行了三次表型鉴定试验，并以GLM和MLM两种模型下均能检测到的显著SNP位点(P<1E-3)为与菊花耐涝性紧密关联的标记，在一定程度上较少了外界因素对关联分析结果的影响。
(2)选择了最有可能和耐涝性显著关联的一个SNP进行dCAPS标记开发，结合该位点的基因型和表型数据首先分析得到预期的酶切情况，再和实际获得的酶切扩增多态性进行比较，提高了标记开发的准确性。
(3)对开发的dCAPS标记在两个群体中进行验证，增加了标记开发的准确性。此外，本发明开发的是和菊花耐涝性共分离的功能性dCAPS共显性标记，初步认为可以应用于菊花耐涝性分子标记辅助育种，这样不用进行田间统计即可快速筛选菊花耐涝性，既避免了环境因子对表型鉴定的影响又大大缩短了育种周期，对于培育耐涝性菊花新品种具有重要的理论和实践意义。
附图说明
图1为三次耐涝性表型鉴定数据的箱图。
图2为WT-dCAPS1的PCR扩增片段序列及引物设计图解。
其中，双下划线处的“C”碱基为人工引入的错配碱基，大括号内的碱基[T/C]表示SNP突变位点，箭头表示酶切位点。
图3为Marker6619-75在关联分析群体中两种基因型的耐涝性隶属函数值(MFVW)箱图。
其中，标号为1的箱子表示C基因型；标号为2箱子表示Y基因型；箱子上方的不同大、小写字母表示两者的表型值在P<0.01和P<0.05水平上均具有显著性差异。
图4为开发的WT-dCAPS1标记在自然群体(A)和F1后代群体(B)中的验证。
其中，1-13为13个耐涝品种，2-26为13个不耐涝品种；1为菊花品种‘南农雪峰’，14为菊花品种‘蒙白’，F1群体为‘南农雪峰’和‘蒙白’杂交所得；a-m为13个耐涝株系，n-z为13个不耐涝株系。
具体实施方式
下面结合实施例对本发明做进一步说明，下列实施例中未注明具体条件的实验方法，通常按照本领域的公知手段。
(一)挖掘与菊花耐涝性显著关联的SNP
(1)试验材料和耐涝表型数据的获得：选取不同来源且无直接关系的88份菊花资源作为试验材料，所有材料保存于“中国菊花种质资源保存中心”，本领域技术人员可以从“中国菊花种质资源保存中心”中获得上述种质。采用盆栽模拟淹水法对该自然群体在菊花苗期进行三次耐涝性鉴定试验，淹水处理保持超过水面3cm，淹水3d后统计萎蔫指数，8d后根据形态变化进行打分并统计黄叶率，采用隶属函数法对以上三个指标进行每个品种的菊花耐涝性隶属函数值(MFVW)的计算，以此作为各品种最终的耐涝性鉴定数据，其中MFVW值越大，说明耐涝性越好。
利用SPSS 20.0软件对三次耐涝性鉴定数据进行方差分析和表型箱图绘制，发现三次试验之间无显著性差异(图1)，故将三次试验的MFVW平均值用于后续的GWAS分析(表1)。
表1供试的88个菊花品种及其MFVW
Figure BDA0001478392120000061
Figure BDA0001478392120000071
(2)分子数据的获得：本发明所采用的SNP标记为现有技术已公开的数据。首先对468,521个SNP原始数据在88份资源中进行筛选，筛选标准为：SNP完整度大于50％且最小等位基因频率(Minor Allele Frequency，MAF)大于5％。过滤后获得92,811个有效SNP。基于此，通过ADMIXTURE软件分析样品的群体结构，假设样品的分群数(K值)为1-10，进行聚类，对聚类结果进行交叉验证，根据交叉验证错误率的谷值确定该群体的最优分群数为3，因此在K＝3时生成的Q-矩阵用于后续GWAS分析；使用SPAGeDi软件对该自然群体两两个体间的亲缘关系(relative kinship)进行估计。亲缘关系本身是定义两特定材料之间的遗传相似度与任意材料之间的遗传相似度的相对值，因此当结果出现两材料之间的亲缘关系值小于0时，则直接定义为0，据此建立亲缘关系K-矩阵。
(3)与菊花耐涝性显著关联的SNP的确定：根据不考虑或者考虑(2)中群体结构和亲缘关系对于关联分析结果的影响，采用GLM和MLM两种分析模型，结合(1)中所获得表型数据通过TASSEL 5.0进行GWAS分析。以两种模型均能检测到的显著SNP位点(P<1E-3)为与菊花耐涝性紧密关联的标记，同时得出该标记的贡献率R2。共检测到10个符合该要求的SNP(表2)。进一步分析发现标签Marker6619的75位置处的SNP位点(即Marker6619-75)在两个模型下的P值均最小，且在GLM模型下的贡献率也最大，为20.67％，在MLM模型下的贡献率为16.90％，因此选择该SNP位点进行后续dCAPS标记开发。
表2与菊花耐涝性显著相关的SNP(P<0.001)
Figure BDA0001478392120000081
(二)SNP突变位点特异性酶切位点分析及dCAPS引物设计
(1)特异性酶切位点分析：根据(一)中获得的与耐涝性最为显著关联的Marker6619-75位点，提取SLAF标签Marker6619的序列(双末端测序各100bp)：CACCATTTCTTCTTCAGTGGACATTTCCGAACTATCTATTTCTAACGTCTCAATGCACAAGGCATTGGTCTTGT[Y/C]AGCAGCACTAGCAGAAAGACTTGCANNNNNNNNNNTGAATTCTGATGTATTCACAATCTTTAACTTCTTCTGTGGAGAATAGCCAAAACTTTGAGGTAAATTAAATTTGTCTTTATCATATTCATAACCATTTCG。利用在线酶切识别软件dCAPS Finder 2.0(http://helix.wustl.edu/dcaps/dcaps.html)查找Marker6619-75突变所引起的限制性内切酶信息，选用内切酶NheI，该酶的酶切识别序列为GCTAGC，需要在Marker6619-74位置处人工引入错配碱基(T→C)。
(2)基因型与表型分析：对(1)中的Marker6619-75位点在GWAS分析群体中进行基因型分析，共有C和Y两种基因型，其中Y代表CT兼并碱基，若某材料的基因型为Y表示该材料在该位点为杂合状态；根据上述(1)中NheI的识别序列可知，若该SNP位点为C，不能被酶切，若为Y，则可以部分酶切；对Marker6619-75位点的两种基因型分别计算平均表型值，发现C基因型的平均MFVW为0.62±0.15，Y基因型的平均MFVW为0.29±0.11。用Student’s t检验方法分析耐涝差异显著性，发现C基因型的材料耐涝性显著高于Y基因型(P<0.01)，综合上述(1)和(2)得出不耐涝的菊花品种的PCR产物更有可能被酶切。
(3)引物设计：利用Primer Primer 5.0软件设计PCR引物，引物设计标准为长度18-25bp；GC含量40-60％；退火温度为55℃-65℃，且上下游引物的GC含量和退火温度尽量保持接近；无引物二聚体及发卡结构；上游引物必需包含错配碱基而决不能包含突变碱基，据此设计的上游引物序列为5’-ATGCACAAGGCATTGGTCTTGC-3’，下游引物序列为5’-CCTCAAAGTTTTGGCTATTCTCC-3’。引物设计详情见下表3。
表3WT_dCAPS1标记引物信息
Figure BDA0001478392120000091
(三)PCR扩增及产物测序
(1)采用CTAB法提取13份耐涝菊花品种材料和13份不耐涝菊花品种材料的基因组DNA(表4)，用核酸仪检测DNA的浓度和纯度后，用ddH2O稀释至300ng/μL备用。
表4用于WT-dCAPS1验证的26份菊花品种
Figure BDA0001478392120000092
Figure BDA0001478392120000101
(2)用上述步骤(二)中设计的特异引物对(1)中的26份材料进行PCR扩增。扩增体系为25μL，其中包括：10×Buffer(含Mg2+)2.5μL，dNTP(2.5mM each)2.0μL，前后引物各1μL，TaqDNA聚合酶1U，模板DNA 1μL，ddH2O 17.3μL。PCR扩增程序为：94℃预变性3min；94℃变性30s，57℃退火30s，72℃延伸30s，共35个循环；72℃延伸7min。对PCR产物采用1.5％琼脂糖凝胶电泳检测，获得单一条带，说明引物特异性好。
(3)由于SLAF技术是针对双末端各100bp进行测序，中间的一段序列为未知，因此对(2)中获得的PCR产物进行测序，获得目的片段全长，共170bp：ATGCACAAGGCATTGGTCTTG|C[Y/C]AGCAGCACTAGCAGAAAGACTTGCATTAGTAATGTTTTCCCCCGGTTCATTGTGTTTGGCAAATGCATCTATCAGAGGACTGCTTCTGAATTCTGATGTATTCACAATCTTTAACTTCTTCTGTGGAGAATAGCCAAAACTTTGAGG。NheI的识别序列为GCTAGC，竖线位置为NheI酶切位点，再根据(二)中基因型与表型的分析，推断耐涝菊花品种不能被酶切，而不耐涝菊花品种能够被部分酶切，且酶切产物分别为分子量大小为149bp和21bp的两个条带。值得注意的是，该步骤的PCR产物测序并非必要步骤，但是通过对产物测序既可以验证引物设计的准确性，也能因此推断酶切条带的大小。
(四)dCAPS标记在自然耐涝极端群体的验证
对(三)中获得的PCR产物进行37℃酶切3h。酶切反应体系为50μL，其中包括：10×NEBuffer5μL，PCR产物10μL，ddH2O 34μL，NheI内切酶0.1μL。所用试剂购自New EnglandBiolabs公司，酶切结束后在80℃下热失活20mim。用10％非变性聚丙烯酰胺胶对酶切产物进行电泳检测、银染法染色。酶切扩增多态性分析发现：13份耐涝菊花品种材料中11份不能酶切(84.6％)，而13份不耐涝菊花品种材料中9份能够部分酶切(69.2％)，酶切条带大小符合(三)中的预期(由于另一条21bp的酶切条带太小，电泳时沉到电泳液中不能显示)，所开发的和菊花耐涝性状共分离的WT-dCAPS1标记平均准确率为76.9％(图4A)。
(五)dCAPS标记在F1后代耐涝极端群体的进一步验证
由(四)可知耐涝菊花品种‘南农雪峰’不可以被酶切，而不耐涝菊花品种‘蒙白’可以被部分酶切(图4A)。因此本发明在‘南农雪峰’ב蒙白’的F1杂交群体中(162个株系)选择13个极端耐涝株系和13个极端不耐涝株系进一步验证开发的WT-dCAPS1标记和菊花耐涝性状共分离(表5)。DNA提取、PCR扩增、酶切反应、电泳检测等均与上述26个自然品种的验证方法相同。酶切扩增多态性分析发现：13个耐涝株系中10份不能酶切(76.9％)，而13份不耐涝菊花品种材料中11份能够部分酶切(84.6％)，平均准确率为80.8％(图4B)，从而进一步验证了本发明开发的鉴定菊花耐涝性的dCAPS标记的可行性。
表5用于WT-dCAPS1验证的26个‘南农雪峰’ב蒙白’F1后代株系
Figure BDA0001478392120000111
以上是为了阐述本发明的发明原理而叙述的本发明的示例性实施方式和具体的操作过程，然而本发明的保护范围不仅限于此，本领域技术人员在不脱离本发明实质情况下，可以做出各种改进和变更，这些改进和变更也属于本发明的保护。
序列表
<110> 南京农业大学
<120> 一种鉴定菊花耐涝性的dCAPS标记开发及应
<160> 3
<170> SIPOSequenceListing 1.0
<210> 1
<211> 22
<212> DNA
<213> 人工序列(Artificial Sequence)
<400> 1
atgcacaagg cattggtctt gc 22
<210> 2
<211> 23
<212> DNA
<213> 人工序列(Artificial Sequence)
<400> 2
cctcaaagtt ttggctattc tcc 23
<210> 3
<211> 210
<212> DNA
<213> 人工序列(Artificial Sequence)
<400> 3
caccatttct tcttcagtgg acatttccga actatctatt tctaacgtct caatgcacaa 60
ggcattggtc ttgtcagcag cactagcaga aagacttgca nnnnnnnnnn tgaattctga 120
tgtattcaca atctttaact tcttctgtgg agaatagcca aaactttgag gtaaattaaa 180
tttgtcttta tcatattcat aaccatttcg 210

Claims (6)

1.一种鉴定菊花耐涝性的dCAPS标记的筛选方法，其特征在于包括以下步骤：
a、通过GWAS方法挖掘与菊花耐涝性显著关联的SNP位点
选取多份来源不同且无直接亲缘关系的菊花品种，采用盆栽模拟淹水法对该群体进行了三次耐涝性鉴定，获得表型数据；基于已有的SNP标记数据分别采用ADMIXTURE软件和SPAGeDi软件对该群体进行分析，获得群体结构Q矩阵和亲缘关系系数K矩阵；在TASSEL 5.0软件中，根据三次耐涝鉴定试验表型数据的平均值，利用GLM和MLM两种分析模型进行全基因关联分析，以两种模型均能检测到的P<1E-3的显著SNP位点为与菊花耐涝性紧密关联的位点，同时得出该标记的贡献率R 2；
b、SNP突变位点特异性酶切位点分析及dCAPS引物设计
（1）从步骤a所得的与菊花耐涝性紧密关联的位点中，选择P值最小，且贡献率较大的SNP位点进行后续dCAPS标记开发；
（2）利用在线酶切识别软件dCAPS Finder 2.0查找步骤（1）中所选择的SNP位点突变所引起的限制性内切酶信息，选用合适的内切酶，且人工引入错配碱基；
（3）利用Primer Primer 5.0 软件设计PCR引物，引物设计标准为18-25bp；GC含量40-60%；退火温度为55℃-65℃，且上下游引物的GC含量和退火温度尽量保持接近；无引物二聚体及发卡结构；
步骤b的（1）中选择的P值最小且贡献率较大的SNP位点为标签Marker6619的75 bp位置处的SNP位点（Marker6619-75）；
步骤b的（2）中选用识别序列为GCTAGC的内切酶NheI，在Marker6619-74处引入错配碱基（T→C），步骤b的（3）中设计的引物序列分别为：
F 5’-ATGCACAAGGCATTGGTCTTGC-3’；
R 5’-CCTCAAAGTTTTGGCTATTCTCC-3’；
标签Marker6619的序列为：
CACCATTTCTTCTTCAGTGGACATTTCCGAACTATCTATTTCTAACGTCTCAATGCACAAGGCATTGGTCTTGT[Y/C]AGCAGCACTAGCAGAAAGACTTGCANNNNNNNNNNTGAATTCTGATGTATTCACAATCTTTAACTTCTTCTGTGGAGAATAGCCAAAACTTTGAGGTAAATTAAATTTGTCTTTATCATATTCATAACCATTTCG；
c、结合基因型与表型数据预期酶切扩增多态性
对步骤b中选择的SNP位点在该群体进行基因型分析，对选择的SNP位点的基因型分别计算平均表型值，并用Student’s t检验方法分析耐涝差异显著性，结合基因型与表型值预期酶切扩增多态性；
d、dCAPS标记在自然耐涝极端群体的验证
（1）提取多份耐涝菊花品种材料和多份不耐涝菊花品种材料的基因组DNA，用上述步骤b中设计的特异引物对耐涝菊花品种材料和不耐涝菊花品种材料进行PCR扩增，获得目的片段全长；
（2）将扩增产物用步骤b中选用的内切酶进行酶切，并对酶切产物进行电泳检测、银染法染色；酶切扩增多态性分析，如果检测结果与步骤c中的预期相符，则将该SNP位点作为鉴定菊花耐涝性的dCAPS标记，将该标记命名为WT-dCAPS1标记；
e、在F1后代耐涝极端群体中对WT-dCAPS1标记进行进一步验证
选取耐涝性菊花品种和不耐涝性菊花品种进行杂交，在F1杂交群体中选择极端耐涝株系和极端不耐涝株系进一步验证所筛选的WT-dCAPS1标记。
2.根据权利要求1所述的方法，其特征在于：步骤c中对Marker6619-75位点进行基因型分析，共有C和Y两种基因型，Y代表CT兼并碱基，根据内切酶NheI的识别序列可知，若该SNP位点为C，不能被酶切，若为Y，则可以部分酶切；对Marker6619-75位点的两种基因型分别计算平均表型值，并用Student’s t检验方法分析耐涝差异显著性，显示C基因型的材料耐涝性显著高于Y基因型，即不耐涝的菊花品种的PCR产物更有可能被酶切。
3.根据权利要求1所述的方法，其特征在于：步骤d的（2）中将扩增产物用内切酶NheI进行酶切，耐涝菊花品种材料的扩增产物不能酶切；不耐涝菊花品种材料的扩增产物酶切产物的电泳检测结果包含两条主带170 bp和 149bp。
4.权利要求1~3中任一所述的方法筛选的鉴定菊花耐涝性WT-dCAPS标记。
5.权利要求4所述的鉴定菊花耐涝性WT-dCAPS标记在耐涝性菊花育种中的应用。
6.一种菊花耐涝性dCAPS分子标记辅助选择育种方法，其特征在于：该dCAPS分子标记为Marker6619-75，选用识别序列为GCTAGC的内切酶NheI，设计上下游引物序列：F 5’-ATGCACAAGGCATTGGTCTTGC-3’和R 5’-CCTCAAAGTTTTGGCTATTCTCC-3’，以菊花的基因组DNA为模板，以所设计的上下游引物进行PCR扩增，获得目的片段全长，共170bp；将扩增产物用内切酶NheI进行酶切，如果酶切产物只包含一条170bp的主带，则待测菊花为耐涝菊花品种；若酶切产物电泳检测结果包含两条主带170bp和 149bp，则待测菊花为不耐涝菊花品种。
CN201711176758.0A 2017-11-23 2017-11-23 一种鉴定菊花耐涝性的dCAPS标记开发及应用 Active CN107815502B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711176758.0A CN107815502B (zh) 2017-11-23 2017-11-23 一种鉴定菊花耐涝性的dCAPS标记开发及应用

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711176758.0A CN107815502B (zh) 2017-11-23 2017-11-23 一种鉴定菊花耐涝性的dCAPS标记开发及应用

Publications (2)

Publication Number Publication Date
CN107815502A CN107815502A (zh) 2018-03-20
CN107815502B true CN107815502B (zh) 2021-06-08

Family

ID=61609621

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711176758.0A Active CN107815502B (zh) 2017-11-23 2017-11-23 一种鉴定菊花耐涝性的dCAPS标记开发及应用

Country Status (1)

Country Link
CN (1) CN107815502B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108796053A (zh) * 2018-06-29 2018-11-13 四川农业大学 一种拟南芥基因突变体的鉴定方法
CN116564407B (zh) * 2023-04-10 2024-03-15 南京农业大学 一种基于全基因组选择高效预测菊花花期的方法
CN116469466B (zh) * 2023-04-11 2024-02-09 南京农业大学 一种高效预测菊花耐涝性的方法及其应用

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101387626A (zh) * 2008-09-02 2009-03-18 南京农业大学 一种菊花耐涝性的评价鉴定方法
CN104498475A (zh) * 2014-11-26 2015-04-08 南京农业大学 一种菊花抗旱性关联分子标记的获得方法及应用
CN105112403A (zh) * 2015-09-15 2015-12-02 南京农业大学 菊花耐盐性关联分子标记及其获得方法和应用
CN106446596A (zh) * 2016-07-15 2017-02-22 南京农业大学 一种与菊花耐涝性显著相关的分子标记及其鉴定方法和应用

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101387626A (zh) * 2008-09-02 2009-03-18 南京农业大学 一种菊花耐涝性的评价鉴定方法
CN104498475A (zh) * 2014-11-26 2015-04-08 南京农业大学 一种菊花抗旱性关联分子标记的获得方法及应用
CN105112403A (zh) * 2015-09-15 2015-12-02 南京农业大学 菊花耐盐性关联分子标记及其获得方法和应用
CN106446596A (zh) * 2016-07-15 2017-02-22 南京农业大学 一种与菊花耐涝性显著相关的分子标记及其鉴定方法和应用

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Combining ability, heterosis, genetic distance and their intercorrelations for waterlogging tolerance traits in chrysanthemum;Su, Jiangshuo et al;《EUPHYTICA》;20170228;第213卷(第2期);摘要 *
涝渍胁迫下5种菊花近缘种属植物生理特性;尹冬梅等;《生态学报》;20090430;第29卷(第4期);第2143-2148页 *
菊花及其近缘种属植物耐涝评价体系建立及耐涝性鉴定;尹冬梅等;《植物遗传资源学报》;20091231;第10卷(第3期);第399-404页 *

Also Published As

Publication number Publication date
CN107815502A (zh) 2018-03-20

Similar Documents

Publication Publication Date Title
Zhu et al. Multilocus analysis of nucleotide variation of Oryza sativa and its wild relatives: severe bottleneck during domestication of rice
Wu et al. SNP-based pool genotyping and haplotype analysis accelerate fine-mapping of the wheat genomic region containing stripe rust resistance gene Yr26
Kujur et al. Functionally relevant microsatellite markers from chickpea transcription factor genes for efficient genotyping applications and trait association mapping
Robbins et al. Mapping and linkage disequilibrium analysis with a genome-wide collection of SNPs that detect polymorphism in cultivated tomato
CN108715902B (zh) 梅花垂枝性状snp分子标记及其应用
Wanchana et al. A rapid construction of a physical contig across a 4.5 cM region for rice grain aroma facilitates marker enrichment for positional cloning
Gonzalo et al. Re-evaluation of the role of Indian germplasm as center of melon diversification based on genotyping-by-sequencing analysis
CN107815502B (zh) 一种鉴定菊花耐涝性的dCAPS标记开发及应用
Roychowdhury et al. Molecular marker-assisted technologies for crop improvement
CN111979349B (zh) 控制莲花色性状的主效qtl、snp分子标记及其检测引物和应用
Lee et al. Development of expressed sequence tag derived-simple sequence repeats in the genus Lilium
US20220010325A1 (en) Quantitative trait loci (qtl) associated with shattering-resistant capsules in sesame and uses thereof
Zhou et al. The development of 204 novel EST-SSRs and their use for genetic diversity analyses in cultivated alfalfa
CN109628627A (zh) 水稻广谱抗稻瘟病基因Pigm的SNP标记开发和应用
Kadirvel et al. Genetic markers, trait mapping and marker-assisted selection in plant breeding
US20170022574A1 (en) Molecular markers associated with haploid induction in zea mays
BR112017026015B1 (pt) Métodos para selecionar uma planta de milho com resistência à podridão de colmo causada por antracnose e método para introgredir um alelo qtl associado com a resistência à podridão de colmo causada por antracnose em uma planta de milho
CN102586238A (zh) 稻瘟病抗性基因Pi1功能特异性分子标记Pi1FNP及其方法与应用
TWI635182B (zh) 用於早期鑑別木瓜性別及性別相關性狀之分子標記及其應用
CN109762925B (zh) 水稻抗稻瘟病广谱基因pi9的KASP分子标记及其检测方法和应用
CN108588261B (zh) 一种鉴定位于萝卜R02染色体上的晚抽薹QTL的InDel引物及其应用
CN107365873B (zh) 与谷子叶鞘色性状连锁的分子标记及其应用
CN110540987A (zh) 一种水稻耐低温基因cold1新功能标记的设计及其检测方法
Singh et al. Molecular markers exploited in crop improvement practices
Gaafar et al. Bacterial artificial chromosome-derived molecular markers for early bolting in sugar beet

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
GR01 Patent grant
GR01 Patent grant