CN107619873B - 基于关联分析和KASP开发waxy1基因内分子标记 - Google Patents
基于关联分析和KASP开发waxy1基因内分子标记 Download PDFInfo
- Publication number
- CN107619873B CN107619873B CN201710769206.4A CN201710769206A CN107619873B CN 107619873 B CN107619873 B CN 107619873B CN 201710769206 A CN201710769206 A CN 201710769206A CN 107619873 B CN107619873 B CN 107619873B
- Authority
- CN
- China
- Prior art keywords
- primer
- waxy
- snp
- indel
- corn
- 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.)
- Expired - Fee Related
Links
- 108090000623 proteins and genes Proteins 0.000 title claims abstract description 73
- 238000010219 correlation analysis Methods 0.000 title abstract description 27
- 238000011161 development Methods 0.000 title description 8
- 240000008042 Zea mays Species 0.000 claims abstract description 137
- 235000002017 Zea mays subsp mays Nutrition 0.000 claims abstract description 137
- 235000005824 Zea mays ssp. parviglumis Nutrition 0.000 claims abstract description 90
- 235000005822 corn Nutrition 0.000 claims abstract description 90
- 239000000463 material Substances 0.000 claims abstract description 57
- 238000009395 breeding Methods 0.000 claims abstract description 13
- 230000001488 breeding effect Effects 0.000 claims abstract description 12
- 235000016383 Zea mays subsp huehuetenangensis Nutrition 0.000 claims description 47
- 235000009973 maize Nutrition 0.000 claims description 47
- 239000003147 molecular marker Substances 0.000 claims description 29
- 238000012217 deletion Methods 0.000 claims description 19
- 230000037430 deletion Effects 0.000 claims description 19
- 101000690100 Homo sapiens U1 small nuclear ribonucleoprotein 70 kDa Proteins 0.000 claims description 14
- 101100029173 Phaeosphaeria nodorum (strain SN15 / ATCC MYA-4574 / FGSC 10173) SNP2 gene Proteins 0.000 claims description 14
- 101100094821 Saccharomyces cerevisiae (strain ATCC 204508 / S288c) SMX2 gene Proteins 0.000 claims description 14
- 102100024121 U1 small nuclear ribonucleoprotein 70 kDa Human genes 0.000 claims description 14
- 238000012216 screening Methods 0.000 claims description 12
- 238000000246 agarose gel electrophoresis Methods 0.000 claims description 8
- 108020005345 3' Untranslated Regions Proteins 0.000 claims description 4
- 230000002441 reversible effect Effects 0.000 claims description 4
- 239000012634 fragment Substances 0.000 claims description 3
- 238000012408 PCR amplification Methods 0.000 claims description 2
- 230000000977 initiatory effect Effects 0.000 claims 2
- 239000000126 substance Substances 0.000 claims 2
- 238000000034 method Methods 0.000 abstract description 30
- 239000003550 marker Substances 0.000 abstract description 21
- 238000012098 association analyses Methods 0.000 abstract description 17
- 238000005516 engineering process Methods 0.000 abstract description 15
- 229920002472 Starch Polymers 0.000 abstract description 12
- 235000019698 starch Nutrition 0.000 abstract description 12
- 239000008107 starch Substances 0.000 abstract description 12
- 230000000694 effects Effects 0.000 abstract description 7
- 239000000499 gel Substances 0.000 abstract description 5
- 230000035772 mutation Effects 0.000 abstract description 5
- 238000001514 detection method Methods 0.000 abstract description 4
- 239000011543 agarose gel Substances 0.000 abstract description 3
- 229920002401 polyacrylamide Polymers 0.000 abstract description 2
- 229920000856 Amylose Polymers 0.000 description 42
- 108020004414 DNA Proteins 0.000 description 22
- 238000004458 analytical method Methods 0.000 description 17
- LFQSCWFLJHTTHZ-UHFFFAOYSA-N Ethanol Chemical compound CCO LFQSCWFLJHTTHZ-UHFFFAOYSA-N 0.000 description 16
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 13
- 238000010586 diagram Methods 0.000 description 10
- 239000006228 supernatant Substances 0.000 description 9
- HEDRZPFGACZZDS-UHFFFAOYSA-N Chloroform Chemical compound ClC(Cl)Cl HEDRZPFGACZZDS-UHFFFAOYSA-N 0.000 description 8
- 238000013461 design Methods 0.000 description 8
- 230000002068 genetic effect Effects 0.000 description 8
- 239000000203 mixture Substances 0.000 description 8
- 238000012163 sequencing technique Methods 0.000 description 8
- 239000003153 chemical reaction reagent Substances 0.000 description 7
- 238000002156 mixing Methods 0.000 description 7
- HEMHJVSKTPXQMS-UHFFFAOYSA-M Sodium hydroxide Chemical compound [OH-].[Na+] HEMHJVSKTPXQMS-UHFFFAOYSA-M 0.000 description 6
- 230000015572 biosynthetic process Effects 0.000 description 6
- 238000003205 genotyping method Methods 0.000 description 6
- 238000003786 synthesis reaction Methods 0.000 description 6
- 238000009826 distribution Methods 0.000 description 5
- 235000019441 ethanol Nutrition 0.000 description 5
- 238000002887 multiple sequence alignment Methods 0.000 description 5
- 239000008188 pellet Substances 0.000 description 5
- 239000000243 solution Substances 0.000 description 5
- 101150084750 1 gene Proteins 0.000 description 4
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 4
- 108091092584 GDNA Proteins 0.000 description 4
- KFZMGEQAYNKOFK-UHFFFAOYSA-N Isopropanol Chemical compound CC(C)O KFZMGEQAYNKOFK-UHFFFAOYSA-N 0.000 description 4
- 230000003321 amplification Effects 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 4
- 238000003199 nucleic acid amplification method Methods 0.000 description 4
- 230000003287 optical effect Effects 0.000 description 4
- 102000004169 proteins and genes Human genes 0.000 description 4
- 238000005215 recombination Methods 0.000 description 4
- 230000006798 recombination Effects 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 108091032973 (ribonucleotides)n+m Proteins 0.000 description 3
- ZCYVEMRRCGMTRW-UHFFFAOYSA-N 7553-56-2 Chemical compound [I] ZCYVEMRRCGMTRW-UHFFFAOYSA-N 0.000 description 3
- 229920000945 Amylopectin Polymers 0.000 description 3
- 208000003643 Callosities Diseases 0.000 description 3
- LZZYPRNAOMGNLH-UHFFFAOYSA-M Cetrimonium bromide Chemical compound [Br-].CCCCCCCCCCCCCCCC[N+](C)(C)C LZZYPRNAOMGNLH-UHFFFAOYSA-M 0.000 description 3
- 206010020649 Hyperkeratosis Diseases 0.000 description 3
- 108010039811 Starch synthase Proteins 0.000 description 3
- 230000009418 agronomic effect Effects 0.000 description 3
- 238000007844 allele-specific PCR Methods 0.000 description 3
- 238000000137 annealing Methods 0.000 description 3
- 210000000349 chromosome Anatomy 0.000 description 3
- 230000002860 competitive effect Effects 0.000 description 3
- 238000004925 denaturation Methods 0.000 description 3
- 230000036425 denaturation Effects 0.000 description 3
- 238000013537 high throughput screening Methods 0.000 description 3
- 239000011630 iodine Substances 0.000 description 3
- 229910052740 iodine Inorganic materials 0.000 description 3
- 239000007788 liquid Substances 0.000 description 3
- TWRXJAOTZQYOKJ-UHFFFAOYSA-L magnesium chloride Substances [Mg+2].[Cl-].[Cl-] TWRXJAOTZQYOKJ-UHFFFAOYSA-L 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 239000000047 product Substances 0.000 description 3
- 238000003753 real-time PCR Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000007447 staining method Methods 0.000 description 3
- 238000012795 verification Methods 0.000 description 3
- HRPVXLWXLXDGHG-UHFFFAOYSA-N Acrylamide Chemical compound NC(=O)C=C HRPVXLWXLXDGHG-UHFFFAOYSA-N 0.000 description 2
- 108700028369 Alleles Proteins 0.000 description 2
- 241000196324 Embryophyta Species 0.000 description 2
- 244000061456 Solanum tuberosum Species 0.000 description 2
- 235000002595 Solanum tuberosum Nutrition 0.000 description 2
- 229910000831 Steel Inorganic materials 0.000 description 2
- 150000001413 amino acids Chemical class 0.000 description 2
- 238000003556 assay Methods 0.000 description 2
- 239000011324 bead Substances 0.000 description 2
- 238000009835 boiling Methods 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 238000000227 grinding Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- PHTQWCKDNZKARW-UHFFFAOYSA-N isoamylol Chemical compound CC(C)CCO PHTQWCKDNZKARW-UHFFFAOYSA-N 0.000 description 2
- 238000009405 line breeding Methods 0.000 description 2
- 229910001629 magnesium chloride Inorganic materials 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 238000005065 mining Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000013642 negative control Substances 0.000 description 2
- 229910052757 nitrogen Inorganic materials 0.000 description 2
- 230000010355 oscillation Effects 0.000 description 2
- 239000000843 powder Substances 0.000 description 2
- 239000002244 precipitate Substances 0.000 description 2
- 238000002360 preparation method Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 108060006613 prolamin Proteins 0.000 description 2
- 238000000746 purification Methods 0.000 description 2
- 238000003757 reverse transcription PCR Methods 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 238000002864 sequence alignment Methods 0.000 description 2
- 239000010959 steel Substances 0.000 description 2
- 108091026890 Coding region Proteins 0.000 description 1
- 108010014303 DNA-directed DNA polymerase Proteins 0.000 description 1
- 102000016928 DNA-directed DNA polymerase Human genes 0.000 description 1
- 102000004190 Enzymes Human genes 0.000 description 1
- 108090000790 Enzymes Proteins 0.000 description 1
- 108700024394 Exon Proteins 0.000 description 1
- 102100034343 Integrase Human genes 0.000 description 1
- KDXKERNSBIXSRK-UHFFFAOYSA-N Lysine Natural products NCCCCC(N)C(O)=O KDXKERNSBIXSRK-UHFFFAOYSA-N 0.000 description 1
- 239000004472 Lysine Substances 0.000 description 1
- 108010092799 RNA-directed DNA polymerase Proteins 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000004913 activation Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 239000007853 buffer solution Substances 0.000 description 1
- 238000010804 cDNA synthesis Methods 0.000 description 1
- 238000005119 centrifugation Methods 0.000 description 1
- 235000013339 cereals Nutrition 0.000 description 1
- 238000000546 chi-square test Methods 0.000 description 1
- 238000010367 cloning Methods 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 239000002299 complementary DNA Substances 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000004090 dissolution Methods 0.000 description 1
- 238000004043 dyeing Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000011536 extraction buffer Substances 0.000 description 1
- 238000002866 fluorescence resonance energy transfer Methods 0.000 description 1
- 238000007716 flux method Methods 0.000 description 1
- 235000013305 food Nutrition 0.000 description 1
- 238000007710 freezing Methods 0.000 description 1
- 230000008014 freezing Effects 0.000 description 1
- 238000001502 gel electrophoresis Methods 0.000 description 1
- 238000012248 genetic selection Methods 0.000 description 1
- 238000012165 high-throughput sequencing Methods 0.000 description 1
- 238000009396 hybridization Methods 0.000 description 1
- 238000009399 inbreeding Methods 0.000 description 1
- 238000002372 labelling Methods 0.000 description 1
- 230000000670 limiting effect Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 108020004999 messenger RNA Proteins 0.000 description 1
- 239000002366 mineral element Substances 0.000 description 1
- 210000003205 muscle Anatomy 0.000 description 1
- 108020004707 nucleic acids Proteins 0.000 description 1
- 102000039446 nucleic acids Human genes 0.000 description 1
- 150000007523 nucleic acids Chemical class 0.000 description 1
- 235000015097 nutrients Nutrition 0.000 description 1
- 235000019629 palatability Nutrition 0.000 description 1
- 238000002264 polyacrylamide gel electrophoresis Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 239000002994 raw material Substances 0.000 description 1
- 239000011535 reaction buffer Substances 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000002829 reductive effect Effects 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 238000007480 sanger sequencing Methods 0.000 description 1
- 238000007789 sealing Methods 0.000 description 1
- 238000009394 selective breeding Methods 0.000 description 1
- 238000007873 sieving Methods 0.000 description 1
- 150000003384 small molecules Chemical class 0.000 description 1
- 238000003756 stirring Methods 0.000 description 1
- -1 sugar Chemical class 0.000 description 1
- 238000005382 thermal cycling Methods 0.000 description 1
- 210000001519 tissue Anatomy 0.000 description 1
- 230000001988 toxicity Effects 0.000 description 1
- 231100000419 toxicity Toxicity 0.000 description 1
- 238000011426 transformation method Methods 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
- 239000013598 vector Substances 0.000 description 1
- 235000013343 vitamin Nutrition 0.000 description 1
- 239000011782 vitamin Substances 0.000 description 1
- 229940088594 vitamin Drugs 0.000 description 1
- 229930003231 vitamin Natural products 0.000 description 1
- 238000005406 washing Methods 0.000 description 1
- 238000009941 weaving Methods 0.000 description 1
Images
Landscapes
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
本发明公开了基于关联分析及KASP开发waxy1基因内分子标记;主要是通过全基因组关联分析和候选基因关联分析技术,找到位于waxy1基因内与直连淀粉含量极显著相关的一个InDel和两个SNP标记。所检测到的InDel标记为中国糯玉米的主要突变类型,本发明基于该InDel开发的标记可通过琼脂糖胶对糯玉米材料进行鉴定,相比以往的聚丙烯酰胺胶鉴定更加方便快捷;同时通过对59份玉米自交系材料的waxy1基因序列区域进行配对连锁不平衡检验,验证了该InDel标记与两个SNP标记紧密连锁,并开发了基于KASP技术的两个SNP的引物,可实现对糯玉米材料快速高效准确的鉴定,对糯玉米分子育种具有促进作用。
Description
技术领域
本发明属于生物技术领域,具体涉及基于全基因组关联分析和KASP技术开发waxy1基因内分子标记。
背景技术
玉米是很多国家广泛种植的主要农作物之一,且是重要的粮食、经济和饲料作物。现代玉米是野生玉米经历近万年的历史进化而来,在其过程中形成了遍布全球的独特地理品系及突变品系。这些品系与野生玉米之间以及品系内部之间特定性状的差异实质上是基因组上的某些位点受到长期的人工定向选择下变异积累而导致的。这些位点通常控制着重要的农艺性状,是育种中需要筛选的重要基因位点。因此,如何高效地鉴定、利用和聚合这些位点是当前育种行业面临的问题和挑战。起源于中国西南地区的糯玉米就是这些优异的突变品系之一(Weatherwax P,Randolph L,Sprague G.History and origin of corn incorn and corn improvement.Academic Press,New York,USA,1955.黄玉碧,荣廷昭.我国糯玉米种质资源的遗传多样性和起源进化.作物杂志,1998,1(998):27-80.),其籽粒中的营养成分十分丰富,糖分、蛋白质、赖氨酸等多种氨基酸含量均高于普通玉米,并含有多种维生素和矿质元素,其淀粉分子量比普通玉米小10多倍,适口性好且更易消化吸收,深受中国及东南亚国家人民的喜爱,同时也是加工支链淀粉的重要原料。
种质资源是培育一切优良品种的基础,因此如何高效的筛选和培育糯玉米自交系是目前研究的热点之一。常规的糯玉米自交系选育基本采用自交系选育方法、二环系选育法和回交转育法(刘纪麟.玉米育种学.北京:中国农业出版社,2002.)。这些方法为了保证隐性waxy1基因不丢失,通常需要育种家观察籽粒表型或通过碘染检法查配子体基因型,但这种方式育种周期长且选择效率低,已经无法完全满足当前玉米生产的需要,同时玉米的大部分农艺性状诸如产量、品质等均为数量性状,而数量性状的表型与基因型之间受环境、亲缘关系及群体结构等多种因素的影响(Yu J,Pressoir G,Briggs W H,et al.A unifiedmixed-model method for association mapping that accounts for multiple levelsof relatedness[J].Nature genetics,2006,38(2):203-208.),仅依赖于表型选择的常规育种手段并不准确。
针对以上问题,近年来分子标记在糯玉米育种中的应用得到了广泛的发展,已经报道的有RAPD标记即随机扩增多态性DNA,缺点是它的稳定性和重复性不高;还有SSR标记,它的数量丰富,并覆盖整个基因组,但密度不够,在选育过程中往往易受重组、选择、环境等变化的影响导致基因组中LD衰退(李春艳,刘建华.分子标记技术在我国糯玉米育种中的应用.广东农业科学,2011,38(15):113-115.)。糯玉米主要是因为第9号染色体waxy1基因发生隐性突变导致胚乳中直链淀粉降低甚至缺失而产生的(Demerec M.A case of pollendimorphism in maize.American Journal of Botany,1924:461-464.),在对waxy1基因编码区的多态性位点进行研究中发现,糯玉米第7个外显子和内含子上有30bp的缺失,或是第外10显个子有15bp缺失,同时发现waxy1基因位点的遗传多样性很低(田孟良,黄玉碧,谭功燮等.西南糯玉米地方品种waxy基因序列多态性分析.作物学报,2008,34(5):729-736.FanL,Quan L,Leng X,et al.Molecular evidence for post-domestication selection inthe Waxy gene of Chinese waxy maize.Molecular breeding,2008,22(3):329.)。在以往对糯玉米地方品种waxy1基因序列多态性分析报道中,主要关注了waxy1基因的第9到14外显子,未从全局上进行连锁分析;同时他们开发的标记均通过聚丙烯酰胺凝胶电泳进行鉴定,虽然可以区分小分子,但丙烯酰胺胶具有毒性,且制胶及染色相对复杂,高通量筛选效率低下。
发明内容
本发明的目的在针对上述现有技术存在的操作复杂、选择效率低、不利于高通量筛选等缺点,提供一种基于关联分析和KASP技术开发waxy1基因内分子标记。本发明利用高通量测序技术和含有糯玉米材料的自然群体,通过全基因组关联分析挖掘潜在的与直链淀粉表型紧密连锁的SNP,以应用于今后的育种,同时通过候选基因关联分析技术我们拟从waxy1基因序列挖掘出更多的多态性位点,为设计更灵活高效经济的SNP检测方法(KASP法)打下基础,从而为糯玉米基础材料的创制提供更多高效准确的鉴定手段。
本发明的目的是通过以下技术方案来实现的:
第一方面,本发明涉及一种基于关联分析及KASP技术开发waxy1基因内分子标记的方法,通过碘染法测定了多份玉米自交系的直链淀粉含量并结合其对应的重测序基因型数据,以waxy1作为候选基因进行关联分析,获得waxy1基因内与直链淀粉含量显著相关的一个InDel和两个SNP分子标记,并通过配对连锁不平衡检验验证了InDel标记与两个SNP标记紧密连锁,并利用KASP技术设计了两个SNP的引物可实现对糯质玉米的鉴定。
优选的,所述关联分析具体为:结合表型数据通过全基因组关联分析筛选一些可能存在相互关系的基因作为候选基因,按表型数据分布再选择一个子群体对候选基因进行测序,分析候选基因序列的多态性进行关联分析;或者利用GWAS或QTL mapping分析所获得的候选区间,对该区间的所有基因再与表型进行关联分析。
优选的,所述KASP技术开发SNP引物具体为:对子群体的候选基因的基因型数据进行配对连锁不平衡检验,获得与InDel标记紧密连锁的两个SNP位点,并设计竞争性等位基因特异性PCR(KASP)引物,对SNP进行基因分型,鉴定糯质玉米。
优选的,所述方法包括如下步骤:
S1、对若干已知直链淀粉含量的马铃薯直链淀粉标样进行淀粉含量测定,获得光密度读数与直链淀粉含量的标准曲线;
S2、对待测玉米自交系群体进行同步骤S1的淀粉含量测定,获得光密度度数,并根据标准曲线计算直链淀粉含量;
S3、将群体自交系的表型数据与相应的SNP数据,通过使用GAPIT软件(版本:2.22)中的压缩混合线性模型(PCA群体结构+Kinship亲缘关系)来进行全基因组关联分析,获得显著的SNP区间内的候选基因waxy1进行后续分析;
S4、克隆多份(59份)子群体waxy1基因gDNA全长,并进行多重序列比对,获得Phylip alignment格式的基因型数据,导入Tassel 5使用GLM进行基于候选基因的关联分析,获得在waxy1基因内与直连淀粉显著相关的多态性位点;
S5、使用Tassel 5对子群体的基因型数据进行配对连锁不平衡检验,分析显著性位点间的连锁不平衡程度,获得与中国糯质玉米主要突变类型(InDel标记)显著相关的两个SNP,针对SNP的侧翼序列信息进行KASP引物的设计和合成(将SNP的侧翼序列信息提交LGC公司(Laboratory of the Government Chemist,Hoddeston,UK)进行KASP引物的设计和合成),使用糯质玉米和非糯质玉米对两个SNP引物适用性进行鉴定。
第二方面,本发明涉及一种玉米糯质的InDel分子标记,是位于waxy1基因第7个外显子末端的30~33bp缺失的InDel。
优选的,所述InDel与糯质表型呈现紧密连锁
第三方面,本发明还涉及一种本发明所述的玉米糯质的InDel分子标记的引物,可以通过琼脂糖凝胶电泳进行检测,wxl2引物效果最佳,其正向引物序列如SEQ ID NO.1所示,反向引物序列如SEQ ID NO.2所示:
wxl2_fwd:TCGTGCTACCTCAAGAGCAACTA SEQ ID NO.1,
wxl2_rev:GAGCATGGAGAACGAAGACGG SEQ ID NO.2。
第四方面,本发明涉及一种本发明所述的InDel分子标记在检测玉米基因型、筛选糯质玉米或分子标记辅助选择育种中的用途。
优选的,在利用InDel分子标记筛选糯质玉米中,针对waxy1基因内的InDel缺失,在两端保守区域设计引物获取扩增片段,直接通过高浓度的琼脂糖凝胶电泳对糯玉米材料进行鉴定。
优选的,设计的引物为如SEQ ID NO.1和2所示序列的wxl2引物。
优选的,利用wxl2引物在低直链玉米材料中可扩增出长度为70bp的含标记InDel的特征条带,在正常玉米材料中可扩增出长度为100~103bp的特征条带。
第五方面,本发明还涉及一种玉米糯质的SNP分子标记,包括位于waxy1基因第14个外显子的SNP1和位于3’UTR区的SNP2。
优选的,所述SNP1和SNP2与糯质表型呈现紧密连锁。
优选的,所述SNP1和SNP2与InDel缺失紧密连锁。(r2>0.9)。
第六方面,本发明涉及一种本发明所述的玉米糯质的SNP分子标记的KASP引物,分子标记SNP1对应的WXT引物包括:如SEQ ID NO.3所示的引物1、如SEQ ID NO.4所示的引物2、如SEQ ID NO.5所示的引物3;分子标记SNP2对应的WXT2引物包括:如SEQ ID NO.6所示的引物4、如SEQ ID NO.7所示的引物5、如SEQ ID NO.8所示的引物6:
SEQ ID NO.3:CAGCCTCGGGGTCGCC,
SEQ ID NO.4:CTCAGCCTCGGGGTCGCT,
SEQ ID NO.5:CGGCGCGATCTCCTCGCCTT,
SEQ ID NO.6:CCACTGGCCACCACTACACTAT,
SEQ ID NO.7:CACTGGCCACCACTACACTAC,
SEQ ID NO.8:TTATGTGATATGGACAAGTRTGTGTA。
第七方面,本发明涉及一种如本发明所述的玉米糯质的SNP分子标记在检测玉米基因型、筛选糯质玉米或分子标记辅助选择育种中的用途。
优选的,在利用SNP分子标记筛选糯质玉米中,利用KASP引物进行PCR扩增,WET引物的SNP T在低直链玉米材料中激发HEX荧光,SNP C在非低直链玉米材料中激发FAM荧光;WET2引物的SNP A在低直链玉米材料中激发FAM荧光,SNP G在非低直链玉米材料中激发HEX荧光。
与现有技术相比,本发明具有如下有益效果:
1、建立了基于全基因组关联分析技术并迅速挖掘与表型紧密连锁的分子标记方法;
2、设计了基于第7个外显子末端的30~33bp缺失的分子标记,可以通过3%的琼脂糖凝胶电泳直接识别,适用于低通量的筛选;
3、分别设计了位于waxy1第14个外显子及3’UTR区的两个SNP(WXT及WXT2)的KASP分子标记,通过荧光定量PCR仪对玉米快速分型,适用于高通量的筛选。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述,本发明的其它特征、目的和优点将会变得更明显:
图1为直链淀粉含量全基因组关联分析示意图;其中,A为是基于混合线性模型(MLM)全基因组关联分析Quantile–quantile plot图;B为直链淀粉全基因组关联分析Manhattan Plot图,x轴代表玉米10条染色体的物理位置,单位为bp;y轴代表相应SNP的Pvalue的log10的负对数,虚线代表P value=5×10-8,虚线以上代表该位点SNP与表型极显著相关;C为Manhattan Plot在最显著SNP位点左右各0.2MB区间内所含基因的示意图,并标注各基因的ID(MaizeGDB release 5b.60);
图2为waxy1紧密连锁SNP标记T在低直链淀粉自交系中的分布示意图;其中,A中左半部分的圆为500自交系中含有SNP T的玉米材料,共26份,右半部分的圆为500自交系中直链淀粉含量低于22%的玉米材料,共27份,其中重叠的玉米材料为15份;B为500份自交系材料在最显著SNP位点的不同等位基因的直连淀粉含量的箱线图;
图3为直链淀粉含量在关联分析群体中的分布示意图;其中,A为503份玉米自交系直链淀粉含量的分布;B为503份玉米自交系直链淀粉含量的描述统计;
图4为基于候选基因waxy1的直链淀粉含量的关联分析及配对连锁不平衡分析示意图;其中,A为关联分析结果,虚线以上(P value<5×10-8)表示与直链淀粉含量紧密关联的多态性位点,其中三角代表InDel,圆形代表SNP,标注相应的ID以及距离ATG的物理位置;B为配对连锁不平衡分析示意图,黑色实线突出显示了InDel与两个SNP的连锁不平衡程度;
图5为waxy1基因gDNA全长的多重序列比对示意图;其中,最左侧显示的为多重序列比对所使用的24份玉米自交系,其中2-13号为糯玉米,14-25号为非糯玉米(以B73为参考序列)。
图6为InDel标记的设计及位置信息示意图;其中wxl1~wxl4为条带差异引物,wxl5为条带有无引物,其反向引物20bp中有10bp在缺失上,因此低直链淀粉材料将无法扩出目的序列;
图7为InDel标记对糯玉米材料的鉴定结构示意图;其中,A上半部分为24份样品直链淀粉含量,下半部分为相应的4种条带差异引物的琼脂糖凝胶电泳结果;B为阴性对照wxl5扩增结果;C为wxl2在市售粘糯玉米杂交种中的验证;
图8为基于KASP技术的waxy1基因SNP引物的基因分型散点图;A为WXT引物的分型情况,其中横坐标为HEX荧光值,纵坐标为FAM荧光值;B为WXT2引物的分型情况,其中横坐标为FAM荧光值,纵坐标为HEX荧光值;
图9为PCR反应条件示意图。
具体实施方式
下面结合实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明,但不以任何形式限制本发明。应当指出的是,对本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干调整和改进。这些都属于本发明的保护范围。
实施例
本实施例涉及基于全基因组关联分析和KASP技术开发waxy1基因内分子标记;具体内容如下:
1.材料和方法:
1.1植物材料
本发明收集来源于中国农业大学的465份温带自交系,38份来自于云南省农业科学院的优质蛋白玉米自交系和5份购于亚马逊(由深圳市飞度时代网络有限公司销售)粘糯玉米杂交种(分别为紫粘玉米,彩粘玉米,白粘玉米,白糯玉米,金糯玉米),其它所有玉米遗传材料来自本实验室。除5份糯玉米杂交种外其余全部遗传材料于2013年11月海南三亚播种自交繁殖。
1.2玉米籽粒表观直链淀粉含量(AAC)测定
标准曲线的制作:称取37度烘箱过夜后的马铃薯直链淀粉标样(Sigma.公司,货号:A0512)100mg。分别加入编号的容量瓶中,然后各加1ml无水乙醇,使淀粉浸润,再各加9ml 1N NaOH,沸水浴10min,冷却后,定容至100ml,充分混匀,按下列表1中的比例吸取至50ml离心管。
表1
离心管编号 | 0 | 1 | 2 | 3 | 4 | 5 | 6 |
直链淀粉(ml) | 0 | 0.1 | 0.4 | 0.5 | 0.75 | 1.0 | 1.25 |
水(ml) | 2.5 | 2.4 | 2.1 | 2.0 | 1.75 | 1.5 | 1.25 |
直链淀粉含量(%) | 0 | 4 | 16 | 20 | 30 | 40 | 50 |
支链淀粉含量(%) | 100 | 96 | 84 | 80 | 70 | 60 | 50 |
随后向各离心管加入25ml ddH20,0.5ml 1N HAc及0.5ml I-KI试剂(0.2%I-KI试剂配制方法:2g KI容于少量ddH20中,再加0.2g I2,待溶解后,用ddH20定容至100ml,混匀。),并定容至50ml,充分混匀,放置10min,测620nm光密度读数。同时,配置空白对照(无淀粉溶液其余一致)。
样品淀粉含量测定:玉米干种子,打样机打碎。过80目筛,37℃烘箱过夜。取25mg样品干粉(称取时要迅速),由于玉米籽粒干粉中含有较多的醇溶蛋白会影响最终的显色,故加入80%的乙醇4ml,振荡混匀,80℃水浴20min(每隔5min振荡混匀一次),4000rpm离心10min,倒掉上清以去除大部分醇溶蛋白。随后加入300μl无水乙醇,振荡悬浮起沉淀并立即加入4.5ml 1N NaOH,振荡摇匀,沸水浴10min,取500μl至50ml离心管,加水25ml,加1ml 1NHAc,加I-KI试剂1ml,定容至50ml,混匀,放置10min。配置空白对照(无淀粉溶液其余一致),备用。分光光度计测620nm光密度读数,并根据标准曲线拟合方程计算直链淀粉含量。
1.3DNA提取和纯化
样品DNA的提取采用经典的CTAB法,取玉米幼苗叶片于2ml离心管中,用液氮冷冻并加钢珠至于研磨仪中粉碎,随后加入预热的CTAB提取缓冲液,摇匀,放入65℃水浴锅中水浴60min-90min,其间每10-15min摇动或搅拌一次;离心管取出水浴锅后,置于室温条件5-10min(勿低于15℃,否则CTAB易出现沉淀);加入等体积的氯仿:异戊醇(24:1)至离心管,密封后放入摇床轻摇15-20min;室温条件12,000rpm离心5min后,小心吸取上层清液至新1.5ml离心管中;加2/3体积的异丙醇或2倍体积的无水乙醇至离心管中,缓慢摇动离心管;室温10,000rpm离心5min,倒掉上清。用1ml的75%酒精洗DNA沉淀两次,每次在室温条件下10,000rpm离心5min后,将乙醇液体小心倾倒出(不要倒出核酸沉淀);在室温中轻微晾干DNA沉淀;加入0.2ml ddH2O溶解DNA沉淀。
1.4RNA提取和纯化
取50-100mg(约一到两个玉米胚乳)16天的玉米籽粒胚乳组织于2ml离心管用液氮冷冻,加入350μl去蛋白溶液和200μl氯仿,并加钢珠用研磨仪打碎。取上清,再加200μl氯仿重复一次。取200μl上清,加入1ml Trizol试剂,振荡混匀室温放置5min。加入200μl上清氯仿,混匀冰上5min。12,000rpm,4℃离心15min。取400μl上清液小心转移到新的1.5ml离心管中,加入等量体积的异丙醇,上下颠倒几次混匀,4℃下放置15min。再12,000rpm,4℃离心15min。小心移去上清液,防止RNA沉淀丢失。再用DEPC水配制70%乙醇洗涤1次,12,000rpm,4℃离心5min。尽可能彻底地吸走上清,将RNA放在室温,待乙醇完全挥发掉,最后沉淀用20μl DEPC水溶解。
1.5cDNA合成
cDNA合成使用SuperScriptIII First-Stand synthesis试剂盒(Life Science公司)。在RNase-free的无菌离心管中将2μg RNA与1μl Oligo(dT)和1μl dNTP混匀,用水补齐到10μL,置于65℃水浴5min。迅速取出,冰浴1min。加入以下组分并混匀:10×反应缓冲液2μL;25mM MgCl2 4μL;0.1M DTT 2μL;RNaseOUT1μL和1μL SuperScriptIII RT。将离心管轻轻混匀,快离后放置于50℃水浴锅中50min,最后置于85℃5min使反转录酶失活后备用。
1.6PCR引物设计和反应体系及条件
引物设计使用NCBI网站BLAST中的Primer-BLAST(可同时实现引物设计和引物的特异性检验)。PCR产物仅用于跑胶的使用LifeFeng公司的2xTaq PCR Master mix(货号:PT202-02)扩增,用于测序的使用TOYOBO公司的KOD FX Neo进行扩增;反应体系参照各试剂说明书,PCR反应条件如图9所示。
引物的标注使用SnapgGene(版本3.1.4)软件中的Features功能。多重序列比对使用Geneious(版本9.0.4)的Align/Assemble功能中Multiple Align的ClustalW比对算法以及SnapgGene(版本3.1.4)软件tools中的Align Multiple Sequence功能。
1.7关联分析
使用GAPIT软件(版本:2.22)中的压缩混合线性模型(PCA群体结构+Kinship亲缘关系)来进行GWAS的分析。Kinship矩阵通过软件emmax-kin(版本:beta-07Mar2010)来计算,参数为-v-h-s-d 10。PCA的计算使用GCTA软件(版本:gcta64),分为两步:第一,估计遗传关系矩阵,参数为--autosome-num 10--autosome--make-grm。第二,将遗传关系矩阵放到PCA计算模块,输出特征值和特征向量,参数为:--grm
-pca 20。特征值最大的3个特征向量用于GWAS的分析。分析结果通过shell脚本进行整理并用perl脚本绘制Manhattan Plot。
基于候选基因的关联分析首先设计6对引物克隆waxy1基因gDNA全长,然后测序结果使用Geneious(版本9.0.4)的Align/Assemble功能中的De Novo Assemble进行拼接,拼接后的序列比对使用Align/Assemble功能中的Multiple Align的MUSCLE算法,将比对结果两端裁齐(根据基因注释信息)存为Phylip alignment格式作为基因型数据,将基因型数据导入Tassel 5(版本5.2.28),并通过Filter中的sites功能将基因型数据中的非差异性位点去除,随后导入表型数据,使用Data中的Intersect Jion将基因型数据和表型数据结合,结合后的文件使用Analysis中的GLM(一般线性模型)进行关联分析,Manhattan Plot使用result中的功能实现。
1.8基于KASP技术SNP引物的设计
使用Tassel 5对waxy1比对的gDNA序列进行Pairwise LD分析,找到两个与waxy1功能性缺失标记紧密关联的SNP(r2>0.9),并将其侧翼序列信息提交LGC公司(Laboratoryof the Government Chemist,Hoddeston,UK)进行KASP引物的设计和合成(表2);PCR体系为10μL,其中5μL DNA(10ngμL–1),KASP Master mix 5μL(LGC Genomics,Hoddeston,UK)(组分包括通用的FRET cassette荧光引物,分别含有FAM荧光基团和HEX荧光基团,ROX内参染料,KlearTaq DNA聚合酶,dNTP和MgCl2。),0.14μL KASP Assay mix(含有3种特异性的引物,包括通用反向引物以及根据SNP设计的分别含有FAM荧光基团和HEX荧光基团的正向引物两个)。PCR程序为94℃热激活15min;94℃变性20s,61~55℃退火和延伸60s,10个循环(touch-down,每循环降低0.6℃);94℃变性20s,55℃退火和延伸60s,26个循环。荧光信号检测在Bio-Rad公司的CFX Connect荧光定量PCR仪上完成。可再次在RT-PCR上进行热循环以提高分型效果,程序为95℃变性20s,57℃退火和延伸1min,8个循环。结果数据导出csv格式,编写R脚本绘制基因分型散点图。
表2 KASP引物的基本信息
2.实验结果:
2.1基于全基因组关联分析技术挖掘直链淀粉表型紧密连锁分子标记
近年来,植物全基因组关联分析(GWAS)及基于候选基因的关联分析技术在解析复杂数量性状方面已得到广泛的利用。本发明通过碘染法测定了503份玉米自交系的直链淀粉含量和其对应的重测序基因型数据进行全基因组关联分析,目的在于挖掘潜在的与直链淀粉表型紧密连锁的SNP及开发新型的分子标记对糯玉米材料进行选育。GWAS分析结果表明(见图1),采用MLM对直链淀粉进行关联分析的模型合理有效(图1B),多个SNP位点与直链淀粉含量呈现显著关联(注:位于虚线以上,P value<5×10-8)(图1A),在各SNP的LD(连锁不平衡)衰减距离内,发现了多个注释基因可能与直链淀粉的调控相关,其中在9号染色体的第23,283,117bp位置(物理位置基于B73 RefGen_v2)的SNP与直链淀粉含量最显著相关,在其LD衰减距离内存在编码颗粒结合型淀粉合成酶(GBSS1)的基因waxy1(图1C),GBSS1已被证明是直链淀粉合成的关键酶,waxy1基因的突变将导致籽粒胚乳中支链淀粉含量几乎为100%,只含少量直链淀粉,籽粒呈糯性。
该位置SNP信息为C或T,其中SNP T与低直链表型呈极显著相关(p<0.01,卡方检验),在直链淀粉含量低于11%的14份材料中,有13份材料含有SNP T,仅一份为SNP C,在所有含有SNP T的26材料中直链淀粉含量低于22%的有15个(图2A),同时在500份玉米自交系中,含有SNP T的材料直链淀粉含量普遍偏低,中位数在15%左右,而含有SNP C的中位数在26%左右。因此该SNP位点可以作为分子标记对低直链玉米进行选育。
2.2基于候选基因的关联分析技术开发基因内选择标记
玉米是雌雄异株作物,其在选育过程中易受重组、选择、环境等变化的影响导致基因组中LD衰退。因此在上述结果中所挖掘的SNP虽然与表型紧密关联,但在基因组中的位置不够理想。基于候选基因的关联分析技术是根据表型筛选一些可能存在相互关系的基因进行关联分析,或者利用GWAS或QTL mapping分析所获得的候选区间,对该区间的所有基因再与表型进行关联分析。该方法直接对候选基因进行测序,分析候选基因序列的多态性进行关联分析,因此可以挖掘基因内的功能性标记,使所开发的分子标记更具适用性。
本发明根据已测定的503份玉米自交系直链淀粉含量的分布情况(见图3),挑选59份自交系样品(见表3),以waxy1作为候选基因进行关联分析,挖掘基因内的分子标记,结果如图4所示,经分析发现,有三处多态性差异与直链淀粉含量呈极显著相关,分别是位于ATG起始(位置参考B73基因组)第1,289~1,321位点一个33bp的缺失,命名为InDel,以及第3,044位点和3242位点的两个SNP,命名为SNP1和SNP2,其中InDel和SNP1的-log10(P)值为13.25,对表型的贡献率高达64.52%;SNP2的-log10(P)值为7.79,对表型的贡献率为44.28%。
表3用于候选基因waxy1关联分析和多重序列比对的59份玉米自交系
根据多重序列比对数据(图5,其中显示了waxy1关联分析所得到的3个显著性位点的具体序列信息及和糯玉米和非糯玉米的对应情况)和直链淀粉含量(见表3)进一步分析发现,其中编号2-13的样品均为糯玉米自交系,14-25为正常玉米自交系,InDel普遍存在于糯玉米材料中,缺失位置信息为7号外显子末端的4bp及其相连的7号内含子的29bp,但与有些材料相比为26bp,由于缺失发生在外显子末端,因此影响了mRNA正常剪接,经cDNA测序发现,根据糯玉米材料的不同,7号内含子被不同程度的保留;SNP1位于第14号外显子中部,糯玉米材料中普遍为T,其它材料为C,但编码的氨基酸未发生改变(由GCC变为GCT,都编码的是Ala);SNP2位于3’UTR区,糯玉米材料中主要为A,其它材料主要G。由配对连锁不平衡分析可知,SNP1和SNP2与InDel紧密连锁(r2>0.9)。
综上所述,该30~33bp的InDel缺失以及SNP1和SNP2与糯质表型呈现紧密连锁,同时它们位于waxy1基因内部,在杂交组合中不会被交换掉,可以作为筛选糯质材料的功能性分子标记,相比其他类标记更具有适用性。
2.3基于waxy1基因的InDel标记的开发及验证
在中国糯玉米wx基因种质资源遗传多样性的研究中,也报道了第7个外显子和其中30bp的缺失并开发了分子标记引物,他们将632bp的扩增条带进行测序与比对,鉴定基因型。通常建立于测序之上的方法比较费时费钱,本发明开发了既能准确又能快速鉴定基因型的分子标记。针对waxy1基因内的30~33bp InDel缺失,在两端保守区域设计引物获取合适大小的扩增片段,可以直接通过高浓度的琼脂糖凝胶电泳对糯玉米材料进行鉴定,相比丙烯酰胺凝胶电泳和Sanger测序更加的快捷高效。本发明共设计了5对引物对InDel标记进行验证,5对引物位置如图6所示,序列信息见表4,其中wxl5为阴性对照(图7B)。根据序列比对数据可知,5对引物基本处在序列保守区域,因此具有材料广泛适用性。
本发明使用了3种不同浓度的琼脂糖凝胶(2%、3%,4%)对InDel标记进行鉴定,比较分析发现,3%的琼脂糖凝胶在较易制备的同时能在较短时间(约20min)内区分条带差异,同时由图7A可知,4对引物的扩增条带在12份低直链玉米材料均比正常玉米材料的条带要小;引物之间相比,wxl1和wxl2的分离效果较好,wxl2的分离效果最好;从测序结果可知wxl2的扩增产物在低直链玉米材料中的长度为70bp,正常玉米材料中的长度为100~103bp。
本发明同时用wxl2引物在503份的关联分析群体中复查,结合图2可知,该群体低于17%的玉米材料有15份,wxl2共检测出13份玉米材料含有InDel缺失且完全包含在15份低直链玉米材料中。
表4 waxy1-InDel引物列表
本发明也对wxl2在市售玉米中的兼容性进行了鉴定,图7C中的2、4、6,8和10泳道为正常玉米自交系,3、5、7,9和11泳道为网购粘糯玉米杂交种,粘糯玉米通常含有较低的直链淀粉含量,由图可知,wxl2具有极好的筛选效果。
2.4基于KASP技术的waxy1基因SNP标记的开发及验证
连锁不平衡(linkage disequlhbrium,LD)是指不同基因座上等位基因间的非随机关联,其可以是两个标记间或两个基因/QTL间或一个基因/QTL与一个标记座位间的非随机关联。其用D表示,取值范围为0-1,但其度量一般用归一化的D’值和r2值来检验,当D’或r2等于0时,说明该基因座之间不存在连锁不平衡;当D’或r2等于1时,说明基因座之间完全连锁不平衡,观察一个基因座就可以获得另一个基因座的全部信息。在r2和D’中,D’较多地反映了重组率的影响,而r2可同时考虑重组率和突变率的影响,因此,更能客观地反映不同基因座基因间的连锁不平衡关系。本发明通过关联分析和Pairwise LD分析检测到了两个SNP(SNP1和SNP2)不仅与糯质表型极显著相关,同时还与waxy1功能性缺失标记(InDel)紧密连锁(r2>0.8)(图4),因此通过鉴定两个SNP也可以对糯玉米材料进行筛选。
KASP技术即竞争性等位基因特异性PCR技术,是基于引物末端碱基的特异匹配激发不同的荧光,仅通过RT-PCR检测荧光即可对SNP进行分型,因此通过设计KASP引物可以更高效地对SNP标记进行鉴定。使用鉴定InDel标记相同的24份玉米自交系,基于KASP技术基因分型情况如下:其中WXT引物(基于SNP1)的SNP T对应12份糯玉米材料,激发HEX荧光,C对应12份非糯玉米材料,激发FAM荧光;WXT2引物(基于SNP2)的SNP A对应12份糯玉米材料,激发FAM荧光,G对应12份非糯玉米材料,激发HEX荧光(表2),从两个SNP引物的荧光信号二维图(图8)中可知,两个SNP引物都能将糯玉米及非糯玉米很好的区分开,但相比WXT(图8A),WXT2(图8B)的分型更加聚合,但WXT与waxy1功能性缺失标记的连锁程度更好(图4)。
综上所述,本发明通过全基因组关联分析及候选基因的关联分析技术开发了wx基因内的三个选择标记:30-33bp InDel和两个SNPs。本发明开发验证的wxl2引物扩增的PCR产物(~100bp)比已经报道的(~632bp)要短,可直接通过3%浓度的琼脂糖凝胶电泳就可以比较~30bp的差异,无需测序或聚丙烯酰胺凝胶即可达到快速分离糯玉米和非糯玉米的目的。该方法对实验设备及化学试剂要求低,常规的分子生物学实验室即可完成,但对于大数量样本的基因型鉴定还是费时费力。因此本发明进一步开发了更灵活经济的SNP检测方法,该方法对仪器试剂的要求高,需要荧光定量PCR仪和KASP引物,但是这种基于竞争性等位基因特异性PCR对SNP基因分型的方法既便宜又快速。本发明验证开发的WXT2引物分型效果最好,可以将糯玉米及非糯玉米很好的区分开,这种高通量的方法将在促进筛选玉米重要农艺性状糯性或低直链淀粉中发挥重要作用。
以上对本发明的具体实施例进行了描述。需要理解的是,本发明并不局限于上述特定实施方式,本领域技术人员可以在权利要求的范围内做出各种变形或修改,这并不影响本发明的实质内容。
SEQUENCE LISTING
<110> 上海交通大学
<120> 基于关联分析和KASP开发waxy1基因内分子标记
<130> DAG28641
<160> 16
<170> PatentIn version 3.5
<210> 1
<211> 23
<212> DNA
<213> Artificial Sequence
<220>
<223> wxl2_fwd
<400> 1
tcgtgctacc tcaagagcaa cta 23
<210> 2
<211> 21
<212> DNA
<213> Artificial Sequence
<220>
<223> wxl2_rev
<400> 2
gagcatggag aacgaagacg g 21
<210> 3
<211> 16
<212> DNA
<213> Artificial Sequence
<220>
<223> WXT Primer_AlleleFAM
<400> 3
cagcctcggg gtcgcc 16
<210> 4
<211> 18
<212> DNA
<213> Artificial Sequence
<220>
<223> WXT Primer_AlleleHEX
<400> 4
ctcagcctcg gggtcgct 18
<210> 5
<211> 20
<212> DNA
<213> Artificial Sequence
<220>
<223> WXT Primer_Common
<400> 5
cggcgcgatc tcctcgcctt 20
<210> 6
<211> 22
<212> DNA
<213> Artificial Sequence
<220>
<223> WXT2 Primer_AlleleFAM
<400> 6
ccactggcca ccactacact at 22
<210> 7
<211> 21
<212> DNA
<213> Artificial Sequence
<220>
<223> WXT2 Primer_AlleleHEX
<400> 7
cactggccac cactacacta c 21
<210> 8
<211> 26
<212> DNA
<213> Artificial Sequence
<220>
<223> WXT2 Primer_Common
<400> 8
ttatgtgata tggacaagtr tgtgta 26
<210> 9
<211> 22
<212> DNA
<213> Artificial Sequence
<220>
<223> wxl1_fwd
<400> 9
tccatccata caaggggagg ac 22
<210> 10
<211> 20
<212> DNA
<213> Artificial Sequence
<220>
<223> wxl1_rev
<400> 10
agcatggaga acgaagacgg 20
<210> 11
<211> 23
<212> DNA
<213> Artificial Sequence
<220>
<223> wxl3_fwd
<400> 11
tacctcaaga gcaactacca gtc 23
<210> 12
<211> 23
<212> DNA
<213> Artificial Sequence
<220>
<223> wxl3_rev
<400> 12
ggagatgttg tggatgcaga aag 23
<210> 13
<211> 22
<212> DNA
<213> Artificial Sequence
<220>
<223> wxl4_fwd
<400> 13
tccatccata caaggggagg ac 22
<210> 14
<211> 22
<212> DNA
<213> Artificial Sequence
<220>
<223> wxl4_rev
<400> 14
cctggtagga gatgttgtgg at 22
<210> 15
<211> 23
<212> DNA
<213> Artificial Sequence
<220>
<223> wxl5_fwd
<400> 15
gcctcaacaa caacccatac ttc 23
<210> 16
<211> 20
<212> DNA
<213> Artificial Sequence
<220>
<223> wxl5_rev
<400> 16
gaaggcaacc tttgcgtccc 20
Claims (7)
1.一种玉米糯质的InDel分子标记,其特征在于,是位于waxy1基因第7个外显子末端的30~33bp缺失的InDel;位置参考B73基因组,是位于ATG起始第1,289~1,321位点一个33bp的InDel。
2.一种如权利要求1所述的玉米糯质的InDel分子标记的引物,其特征在于,可以通过琼脂糖凝胶电泳进行检测,所述引物包括wx12引物,其正向引物序列如SEQ ID NO.1所示,反向引物序列如SEQ ID NO.2所示。
3.一种玉米糯质的SNP分子标记,其特征在于,包括位于waxy1基因第14个外显子的SNP1和位于3’UTR区的SNP2;位置参考B73基因组,是位于ATG起始第3,044位点的SNP1、第3242位点的SNP2。
4.一种如权利要求3所述的玉米糯质的SNP分子标记的KASP引物,其特征在于,分子标记SNP1对应的WXT引物包括:如SEQ ID NO.3所示的引物1、如SEQ ID NO.4所示的引物2、如SEQ ID NO.5所示的引物3;分子标记SNP2对应的WXT2引物包括:如SEQ ID NO.6所示的引物4、如SEQ ID NO.7所示的引物5、如SEQ ID NO.8所示的引物6。
5.一种如权利要求1所述的玉米糯质的InDel分子标记或如权利要求3所述的玉米糯质的SNP分子标记在检测玉米基因型、筛选糯质玉米或分子标记辅助选择育种中的用途。
6.根据权利要求5所述的用途,其特征在于,在利用InDel分子标记筛选糯质玉米中,针对waxy1基因内的InDel缺失,在两端保守区域设计引物获取扩增片段,直接通过高浓度的琼脂糖凝胶电泳对糯玉米材料进行鉴定。
7.根据权利要求5所述的用途,其特征在于,在利用SNP分子标记筛选糯质玉米中,利用KASP引物进行PCR扩增,WXT引物的SNP T在低直链玉米材料中激发HEX荧光,SNP C在非低直链玉米材料中激发FAM荧光;WXT2引物的SNP A在低直链玉米材料中激发FAM荧光,SNP G在非低直链玉米材料中激发HEX荧光。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710769206.4A CN107619873B (zh) | 2017-08-30 | 2017-08-30 | 基于关联分析和KASP开发waxy1基因内分子标记 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710769206.4A CN107619873B (zh) | 2017-08-30 | 2017-08-30 | 基于关联分析和KASP开发waxy1基因内分子标记 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107619873A CN107619873A (zh) | 2018-01-23 |
CN107619873B true CN107619873B (zh) | 2021-10-12 |
Family
ID=61089143
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710769206.4A Expired - Fee Related CN107619873B (zh) | 2017-08-30 | 2017-08-30 | 基于关联分析和KASP开发waxy1基因内分子标记 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107619873B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108456740B (zh) * | 2018-03-02 | 2021-07-27 | 江苏里下河地区农业科学研究所 | 一个水稻稻瘟病抗性位点‘Pi-jx’及其Indel标记引物和育种应用 |
CN109448789A (zh) * | 2018-09-17 | 2019-03-08 | 上海派森诺生物科技股份有限公司 | 一种基于perl语言的种群特异SNP位点的自动化分析方法 |
CN109652582B (zh) * | 2019-01-08 | 2022-07-22 | 上海市农业科学院 | 用于鉴定沪玉糯3号糯玉米的snp分子标记及其应用 |
CN109735653A (zh) * | 2019-03-15 | 2019-05-10 | 广西壮族自治区农业科学院 | 一种水稻蜡质基因Wx的荧光分子标记及其应用 |
CN110157834B (zh) * | 2019-07-01 | 2022-04-19 | 北京市农林科学院 | 与玉米高赖氨酸性状相关的snp位点oe2-5及其应用 |
CN110610744B (zh) * | 2019-09-11 | 2020-10-23 | 华中农业大学 | 一种高效可并行运算且高准确性的全基因组选择方法 |
CN110894539B (zh) * | 2019-12-04 | 2022-01-25 | 云南师范大学 | 鉴定二倍体马铃薯是否自交亲和的方法 |
CN111518940B (zh) * | 2020-05-09 | 2023-03-24 | 华南师范大学 | 一种检测Wx1分子标记及应用于甜糯玉米选育的方法 |
CN113528696A (zh) * | 2021-06-24 | 2021-10-22 | 扬州大学 | 玉米糯性基因的kasp分子标记的开发方法及应用 |
CN114591928B (zh) * | 2022-04-02 | 2023-10-24 | 山东中农天泰种业有限公司 | 一个与玉米籽粒支链淀粉含量相关的dCAPS分子标记 |
CN116467596B (zh) * | 2023-04-11 | 2024-03-26 | 广州国家现代农业产业科技创新中心 | 水稻粒长预测模型的训练方法、形态预测方法及装置 |
CN117144040B (zh) * | 2023-09-06 | 2024-09-17 | 上海市农业科学院 | 一种鲜食玉米基因型分型芯片及其应用 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105603083B (zh) * | 2016-01-29 | 2019-12-17 | 北京市农林科学院 | Ssr与snp两种分子标记配合使用进行玉米辅助高效育种的方法 |
CN105671042A (zh) * | 2016-03-28 | 2016-06-15 | 浙江师范大学 | 稻米直链淀粉含量微控基因ssiii-1分子标记及应用 |
-
2017
- 2017-08-30 CN CN201710769206.4A patent/CN107619873B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN107619873A (zh) | 2018-01-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107619873B (zh) | 基于关联分析和KASP开发waxy1基因内分子标记 | |
CN113957155B (zh) | 与草鱼性状相关的snp分子标记及其应用 | |
CN109628628B (zh) | 水稻抗稻瘟病基因Pi2的SNP标记的开发和应用 | |
CN107227373B (zh) | 一种粳稻抗倒伏基因的snp功能分子标记及应用 | |
CN112195265B (zh) | 一种鉴定辣椒杂交种纯度的snp位点、引物组及应用 | |
CN111471790B (zh) | 与小麦籽粒灌浆速率QTL QGfr.sicau-7D.1紧密连锁的分子标记及应用 | |
CN114717355A (zh) | 一种西瓜全基因组SNP-Panel | |
CN107858447B (zh) | 用于鉴定桃花单瓣/重瓣性状的单核苷酸多态性标记位点、引物对、试剂盒及应用 | |
CN110938635B (zh) | 一种粳稻粒长基因及与其紧密连锁的分子标记和应用 | |
CN109735648B (zh) | 一种筛选不同千粒重小麦的方法及其专用试剂盒 | |
CN106755465B (zh) | 与小麦旗叶长QTL QFll.sicau-2D紧密连锁的分子标记及应用 | |
CN111647677A (zh) | 与小麦籽粒灌浆速率QTL QGfr.sicau-6D紧密连锁的分子标记及应用 | |
CN114717352B (zh) | 水稻耐高温调控基因Hsp70的分子标记及其应用 | |
WO2017133936A1 (en) | Plant breeding using next generation sequencing | |
CN114480718B (zh) | 一种基于kasp技术用于水稻耐高温基因分型的引物组、检测试剂盒及其应用 | |
CN115852033A (zh) | 提高稻米品质的gs3基因和gw5基因的分子标记 | |
CN112143830B (zh) | 水稻剑叶宽调控基因nal1的分子标记及其应用 | |
CN111705157B (zh) | 用于甘蓝自交不亲和系ⅰ类s单元型筛选的pcr标记及引物 | |
CN102443583A (zh) | 外显子保守序列扩增多态性分子标记及其分析方法 | |
CN109652582B (zh) | 用于鉴定沪玉糯3号糯玉米的snp分子标记及其应用 | |
CN109913579A (zh) | 一种大麦磷素高效利用qtl位点的分子标记及应用 | |
CN113736906B (zh) | 用于检测番茄黄萎病抗性的snp位点组合及其应用 | |
CN115044700B (zh) | 水稻抽穗期QTL qHd1功能标记及其应用 | |
CN116622883B (zh) | 一种鉴定水稻h23转化体的特异性探针、引物、试剂盒和方法 | |
CN116397040B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20211012 |