食用黄花菜SSR标记开发及指纹图谱库构建

刘新星,欧巧明*,罗俊杰,李忠旺,石有太,崔文娟

(甘肃省农业科学院生物技术研究所,甘肃 兰州,730070)

摘 要 利用转录组测序技术开发SSR标记,旨在为食用黄花菜种质资源鉴定、评价及遗传多样性分析提供技术支撑。测序共得到104 177个Unigene,总长度为79 106 833 bp, 平均长度为759 bp,利用MISA软件对Unigene进行搜索,共检测到15 890个SSR,其中单核苷酸、二核苷酸和三核苷酸重复占比较高,分别为18.57%、27.23%和44.26%。设计合成SSR引物114对,通过聚丙烯酰氨凝胶电泳最终筛选出20对核心引物,20对引物在43份黄花菜种质中共检测到88个等位变异,每个位点的等位变异数为3~8个,平均为4.4个,多态信息含量(polymorphism information content,PIC)值变幅为0.23~0.77,平均为0.499 5。依据扩增产物的片段大小确定了每个SSR位点的等位变异以及对应的参照品种,建立了基于SSR的DNA指纹鉴定技术体系,并构建了43份黄花菜种质的DNA指纹数据。该研究开发的SSR标记具有丰富的多态性。

关键词 黄花菜;SSR标记;指纹图谱

黄花菜俗称“金针菜”,在植物分类学上属百合科萱草属(Hemerocaliss L.)植物[1],是一种药食同源的特色蔬菜,其花蕾可食用,具有丰富的营养价值,与香菇、木耳、冬笋并称为“四大素山珍”[2]。据报道,黄花菜的挥发性成分多达58种[3-4],含有大量的醇、酯等香气物质,赋予黄花菜特有的美味。近年来有研究表明,食用黄花菜还有抗氧化、抑制纤维原细胞增生及阻止癌细胞生殖等作用[5]

随着分子生物学的发展,DNA 分子标记鉴定因快速准确、不受环境影响等优点逐渐发挥重要作用,满足品种快速准确鉴定和指纹图谱构建的要求[6-7]。其中SSR标记具有共显性、高度重复性和丰富的多态性等优点,通过PCR快速检测就能达到鉴定目的。目前基于SSR分子标记的玉米、大豆、小麦、油菜等主要农作物的指纹图谱库已构建完成[7-11]。关于黄花菜品种DNA指纹图谱及分子标识研究,陈士林院士[12]针对药用黄花菜开展过DNA条形码分子鉴定研究。

在农业科技迅速发展和农产品市场竞争日趋剧烈的今天,真正的特色农产品优势发挥离不开现代农业技术。食用黄花菜作为特色加工型农产品,利用现代农业技术快速发展优良种质、产地标识及分子身份证构建研究已迫在眉睫。本研究拟利用转录组测序技术开发黄花菜的SSR标记,构建黄花菜DNA指纹图谱库,预期研究结果为食用黄花菜种质资源鉴定、评价及遗传多样性分析提供技术,进一步为品种更新及产业扶贫提供科技支撑。

1 材料与方法

1.1 实验材料

供试43份食用黄花菜种质资源由甘肃省农科院黄花菜种质资源圃提供,种质名称及来源见表1。

1.2 实验方法

1.2.1 EST-SSR引物设计与筛选

采用二代测序技术对黄花菜进行转录组测序,并使用MISA软件对所有测序所得的转录本进行SSR位点搜索。搜索标准为:SSR重复最低长度为18 bp,单核苷酸、二核苷酸、三核苷酸、四核苷酸、五核苷酸和六核苷酸重复最低次数为12、6、5、5、4和4次,间隔不完全重复的SSR位点不列为搜索对象。结合SSR位点两端保守的特点,利用Primer 3.0软件进行批量引物设计,共设计引物11 363对,根据基因注释分类,挑选出氨基酸代谢和萜类聚酮化合物代谢相关的Unigene中的SSR引物进行合成,共计89对。另从文献中搜索得到25对萱草属通用的SSR引物[13-14],合计114对引物进行筛选鉴定。

表1 黄花菜种质名称及来源
Table 1 The names and original of Hemerocallis citrina

编码名称来源编码名称来源1丁庄大菜江苏23安民花1号湖南祁东2渠县花四川24祁东四月花湖南祁东3福建红花福建25祁东财桥花湖南祁东4祁东中期花湖南祁东26浙江仙居花浙江仙居5重庆红花重庆27黄花型未知6猛子花湖南祁东28黑紫花型未知7祁东白花菜-1湖南祁东29红花型未知8荆州花湖北30山西大同花山西大同9大乌嘴江苏31邵东采集湖南邵阳10大同黄花菜山西大同32沐阳采集江苏宿迁11倒箭花湖南祁东33潼阳采集江苏宿迁12高箭中期花湖南祁东34未知1未知13龙游红花菜浙江35未知2未知14淮阳花河南淮阳36未知3未知15四川武平早四川37未知4未知16沙苑金针陕西38未知5未知17雪不谢苗湖南祁东39金蕾1号甘肃庆阳18祁珍花湖南祁东40航天2号甘肃庆阳19四月白湖南祁东41小花黄花甘肃庆阳20棒槌花湖南祁东42高亭黄花甘肃庆阳21茄子花湖南祁东43航天选甘肃庆阳22祁东白花菜-2湖南祁东

引物设计的原则为EST序列长度大于200 bp,SSR序列的开始和结束位置分别距5’和3’端不少于20 bp,引物长度18~24 bp,退火温度Tm值40~60 ℃,上游和下游引物的Tm值相差不大于5 ℃,GC含量40%~60%,且上游和下游引物的GC含量相差不要太大,产物预期长度100~500 bp。对各条引物的最佳退火温度进行试验,以每条引物的解链温度(Tm值)为参考,上下各浮动5 ℃,通过梯度PCR仪自动生成8个温度梯度,以扩增条带多且明亮、背景清晰者为最适温度选择。

1.2.2 样品DNA提取

样品由硅胶颗粒迅速干燥后带回,每个样品取约50 mg,球磨仪研磨后用天根的植物基因组DNA提取试剂盒进行样品DNA的提取,提取完成后用1%的琼脂糖凝胶电泳进行样品DNA质量的检测,并用超微量分光光度计测定其浓度和纯度,将DNA稀释至20~50 ng/μL,4 ℃保存备用。

1.2.3 PCR扩增

PCR扩增体系总反应体积15 μL,包括2×Taq Master Mix 7.5 μL,DNase-Free Water 4.5 μL,Forward Primer 1 μL,Reverse Primer 1 μL,模板DNA 1 μL。反应程序:95 ℃热启动3 min,95 ℃变性45 s,48~65 ℃ 退火30 s,72 ℃延伸1 min,30个循环后,72 ℃终延伸5 min后结束。

1.2.4 非变性聚丙烯酰氨凝胶电泳与生物分析仪检测

利用1.5%的琼脂糖凝胶和6%的非变性聚丙烯酰胺凝胶进行引物的初筛和多态性筛选,待筛选出核心引物后,使用核心引物对所有供试材料进行PCR扩增,扩增产物稀释10倍后上样至LabChipGx Touch大分子生物分析仪进行自动化检测,软件自动收集电泳数据。

1.2.5 数据分析

利用Excel软件统计分析数据,NTSYS软件的UPGMA进行聚类分析。利用 PIC(polymorphism information content)评价引物多态性,PIC计算如公式(1)所示:

(1)

式中:Pij表示标记i的第j个等位基因在群体中的频率[15]

2 结果与分析

2.1 黄花菜转录组SSR位点分布频率与基序特征

采用Illumina Hiseq技术,从黄花菜样品中测序共得到104 177个Unigene,总长度为79 106 833 bp, 平均长度为759 bp,GC含量为41.99%。利用MISA软件对Unigene进行搜索,共检测到15 890个SSR,分布于12 492个Unigene,占总序列的11.99%,含有2个或2个以上SSR的Unigene有2 545个,复合SSR有1 316个,占总序列的8.28%。从黄花菜EST-SSR中基元分布频率情况(表2)看,单核苷酸到六核苷酸基元类型均有分布,其中三核苷酸重复占比最高,为44.26%,二核苷酸和单核苷酸次之,分别为27.23%和18.57%,四、五、六核苷酸占比较低,不足5%。从重复次数看, 黄花菜转录组SSR重复基序的重复次数分布在4~115次, 其中4~8次的SSR个数占总数的69.93%;其次为9~12次的SSR, 共有2 095个, 占总数的13.18%;13次以上重复的有 2 699个,占16.99%(表2)。其中,单核苷酸重复次数以≥13次为主,占单核苷酸总重复次数的76.37%,二核苷酸重复次数主要以6~7次重复为主,占二核苷酸总重复次数的53.78%,三核苷酸重复中以5~6次重复为主,占三核苷酸总重复次数的72.64%。

从黄花菜基元类型分布图(图1)看,单核苷酸重复基元类型以(A/T)n为主,占单核苷酸重复总数的96.387%;二核苷酸重复基元类型以(AG/CT)n和(AT/AT)n为主,分别占比64.16%和26.25%;三核苷酸重复基元类型以(AAG/CTT)n和(AGG/CTT)n较多,占比分别为26.87%和15.23%,其余还有8种基元类型,占比从2.55%到10.93%均有分布;四、五、六核苷酸重复基元类型以(AAAT/ATTT)n、(AAAAT/ATTTT)n和(AACCCT/AGGGTT)n最多。

表2 黄花菜EST-SSR的分布频率
Table 2 The distribution of the SSR motifs in Hemerocalli scitrina

基序类型EST-SSR不同重复数/次456789101112≥13总计单核苷酸6972 2532 950二核苷酸1 4089196263392412122313514 327三核苷酸3 4841 625815654881309053947 033四核苷酸264981911463405五核苷酸2951071921424六核苷酸611853714211751总计 9063 9403 1871 7691 2944323773059812 69915 890 分布频率/%5.7024.8020.0611.138.142.712.371.926.1716.99

图1 黄花菜EST-SSR基元类型分布
Fig.1 The distribution of EST-SSR types in Hemerocallis citrina

2.2 EST-SSR核心引物筛选及多态性分析

以2份黄花菜种质(祁东白花菜和金蕾1号)为材料,采用PCR扩增琼脂糖凝胶电泳检测,对114对SSR引物进行初步筛选,共选出扩增条带清晰、目标条带大小相一致的引物70对,有效引物比率为60.87%;选择6份不同地方来源的黄花菜作为材料,利用6%非变性聚丙烯酰氨凝胶电泳对70对初选引物进行多态性引物筛选,共筛选出条带清晰、多态性高、扩增稳定的引物33对,多态性引物扩增率为26.09%;进一步选用供试43份黄花菜种质材料,对上述筛选出的33对引物进行PCR扩增、6%非变性聚丙烯酰胺凝胶电泳检测,最终筛选出20对品种类型间多态性高、条带清晰易读取的SSR引物作为黄花菜指纹图谱构建的核心引物(表3)。20对引物中二核苷酸重复4对,三核苷酸重复13对,四、五、六核苷酸重复各占1对,重复次数主要以5~6次为主。

20对核心引物的PIC值变幅为0.23~0.77,平均为0.499 5。其中,引物HH42的PIC值最高,HH78的PIC值最低。20对核心引物中有12对引物为高度多态位点(PIC>0.50),占比60%,7对引物为中度多态位点(0.25≤PIC≤0.50),占比35%,1对引物为低度多态位点(PIC<0.25),占比5%,表明筛选出的20对核心引物在参试品种间具有丰富的多态性。

2.3 指纹图谱库的构建

20对核心引物对43份黄花菜种质进行扩增检测,统计各黄花菜种质在20个SSR位点上的扩增片段大小,用片段大小命名各等位变异,建立包含43份黄花菜种质的DNA指纹图谱库(表4),筛选出安民花1号、祁东四月花、棒槌花、高箭中期花、雪不谢苗、祁东中期花、四川武平旱共7个黄花菜品种的特异标记。在20个位点中共检测到88个等位变异,每个位点的等位变异数为3~8个,平均为4.4个。

2.4 聚类分析

根据NTSYS软件genemic similiarity计算,43份黄花菜种质材料间的遗传距离变化为0~0.813 8,遗传相似度变化范围为0.443 2~1,从聚类图(图2)中看出,其中有5组(I-V)材料,每组材料间相似度为1,组内无法区分,分别是:I组为渠县花、山西大同花、未知2和未知4;II组为淮阳花、沙苑金针和浙江仙居花;III组为祁珍花和未知黄花型;IV组为猛子花、沐阳采集种质和金蕾1号;V组为重庆红花和龙游红花菜。除去以上14种种质资源,剩余的29份种质均可以进行区分,且以上5组材料组间仍然可以进行区分。

表3 20对EST-SSR引物信息
Table 3 The primer information for 20 pairs of EST-SSR

引物编码引物序列基元类型扩增片段长度PICHHC32FRGGAAGAACAAGCGTGTCAAGAAAGGATCTGATCTGCGAAGATGTCT(ACT)51460.61HHC35FRGAATACTTACTCGTCGGGTCTGGGCTCCAGATCTGCTCAAATTGAC(TTTAT)51400.58HHC36FRAGTATAGGAGGGAGGGTAAGTGGAAGCACAGCCGTTTTCTCAATAC(TGC)51490.49HHC40FRTTATTCCGCGCTTATTATGGGGAGGAGAGAATGAAAGGTGGGGATT(AT)61390.66HHC42FRTACGGGCTTGAGTTGATCTGAAAACCGGGCCTGGTATATACAAAAA(TC)71210.77HHC43FRCAGGGCTTCTACGAGGATGTCATTTCCCAGCTTTAAGATTCGCG(GTG)61450.47HHC45FRAGCTTGAACCCGAAAACTCACTACCTCATCCTCTTCAGAACTACGG(ACC)51320.62HHC52FRCAACGACGACACAAACAAACCTAGAGGAAGAGGAGTGGGATCTGAT(CAT)51470.60HHC53FRAAGTTGTGCAAACCACCCAACACAGTGACGGAAAGATCAACAGT(ACC)51540.57HHC65FRTTGATCTCCTGTGTCACAACCATTCCGTAAATTCAAACATGTAGGAACA(GA)101600.57HHC72FRCTTTACATTCCGACAAACCCAGGCACCTCGTGGTCTGATTATGTGA(TA)61160.76HHC78FRAGTTCAAGGAAGTGACGGAACTGGAGAGTCCATCCTACTGCTGTTC(GCA)61020.23HHC80FRTAGCACTGATCAGTAACGATCCGTACTCCAGGATCAACAGTATGCC(ATG)51320.40HHC90FRGAATTTTGTGCCAGCCATATCGACATGCATGCCTACCAATTCATGT(TCT)91500.62HHC93FRCAACAATGCCATCCCCATCATTTAGAACCAGCCTTGTACTGTTCTC(CCGCTG)41570.56HHC95FRAGCGAAAAGAAAATACACCCAATCTTCCTCCACAATCTGAAACAAACG(TTG)81310.50HHC99FRGATCGAGAATTGAGAAGGCAAGGTGTAGTGAACTTATTTAGTGCTAGC(GTAC)41470.56HHC101FRAAACCGAACTGTACTGCAATTGCATCTGAGCCATTCTAATGCTTCG(AAC)51520.27HHC106FRGCCAAAAGAGAGCTCCTACTTCTAACCATCGTCCAGTATTGGTGAT(CTC)61210.26HHC109FRAGTGGCTCAAAATCAACAACACCTAGTTGGGGACGCTTTAGAAGAC(CTA)61000.45

图2 43份黄花菜种质聚类图
Fig.2 The cluster analysis dendrogram of 43 varietiesin in Hemerocallis citrine

3 讨论

种质资源是植物遗传改良的物质基础,然而,日益庞大的种质数量增加了资源收集保存和遗传改良工作量,导致工作效率低下[15]。通过分子标记鉴定,筛选淘汰同名异物、同物异名的冗余种质,提高种质资源的鉴定保存效率。与其他分子标记相比,SSR 分子标记在基因组上分布均匀、数量丰富、稳定性好且技术成熟,可从DNA水平直接鉴定品种间遗传差异[16]。国际植物新品种保护联盟把SSR和SNP作为作物品种DNA 指纹库构建的分子标记[17],并已经在遗传多样性分析、品种鉴定、群体结构及图谱构建等方面得到应用[18-20]。李益等[21]利用21对SSR核心引物,构建了包含500份柑橘品种的 DNA 指纹图谱库,筛选出部分品种的特异标记,可以鉴别221份柑橘品种和87个品种组合。分子标记在黄花菜上的应用十分滞后,主要有以下几方面的进展:郑家祯等[22]利用SCoT分子标记技术对28份国内栽培的黄花菜资源19个形态学特征进行了遗传多样性分析,李森等[23]对85份萱草属植物开花表型性状进行多样性分析,黎海利等[24]借助 AFLP标记对35份萱草野生种和栽培品种进行亲缘关系研究,但是这些研究很少结合形态学和分子标记进行比对分析,且未专门针对国内食用黄花菜的品种资源进行分析。

本研究所选的20对核心引物能将43份地方种质中的29份有效区分,未能区分开的品种有5组共14个。通过分析14个品种的亲本、植物学性状、来源等相关信息发现,这些品种的特征特性较为接近,例如“重庆红花菜”和“龙游红花菜”这一组,2个品种在开花期花朵颜色、大小、叶形等植物学性状完全相同。“渠县花、山西大同花、未知2和未知4”这组可能存在同一个品种在不同地方因品种名称不一样而被重复收集的情况,也就是同种异名的现象。由于黄花菜育苗仍以农民自己传统育种为主,有些品种可能是在地方品种的基础上提纯复壮而来,两者间的遗传背景一致,因此利用有限数量的分子标记仍难以鉴别。

4 结论

本研究筛选出20对扩增稳定、条带清晰、多态性高的核心引物,依据扩增产物的片段大小确定了每个SSR位点的等位变异以及对应的参照品种,建立了基于EST-SSR的DNA指纹鉴定技术体系,并构建了43份黄花菜种质的DNA指纹数据,可以应用于黄花菜种质鉴定和遗传多样性评价。

表4 黄花菜品种基因型数据
Table 4 Genotype data of d Hemerocallis citrina

引物名称等位变异/bp参照品种名称参照品种基因型数据HHC-32158152149146143祁东中期花、猛子花、金蕾1号福建红花、重庆红花、祁东白花菜-1、倒箭花、龙游红花菜、四川武平早、小花黄花、高亭黄花丁庄大菜、雪不谢苗、祁东白花菜-2渠县花、荆州花、大同黄花菜、淮阳花、沙苑金针、祁东四月花、浙江仙居花、航天2号高箭中期花、四月白安民花1号、祁东财桥花158/143152/146149/146146/143149/143143/143HHC-35145140135130福建红花、重庆红花、祁东白花菜-1、龙游红花菜丁庄大菜、渠县花、祁东中期花、猛子花、荆州花、大同黄花菜、倒箭花、淮阳花、四川武平早、沙苑金针、祁珍花、茄子花、祁东四月花、祁东财桥花、浙江仙居花、金蕾1号、小花黄花、高亭黄花大乌嘴、雪不谢苗、棒槌花高箭中期花、航天2号四月白、安民花1号、145/140140/140140/135135/135130/130HHC-36158149146福建红花、重庆红花、祁东白花菜-1、龙游红花菜丁庄大菜、渠县花、猛子花、荆州花、大乌嘴、大同黄花菜、倒箭花、高箭中期花、淮阳花、四川武平早、沙苑金针、雪不谢苗、四月白、茄子花、祁东白花菜-2、祁东四月花、浙江仙居花、金蕾1号、航天2号、小花黄花、高亭黄花祁东中期花、祁珍花、棒槌花、安民花1号、祁东财桥花158/146149/149149/146HHC-40143141139渠县花、淮阳花、大同黄花菜、沙苑金针、浙江仙居花、航天2号荆州花、倒箭花四川武平早、棒槌花、茄子花、祁东财桥花祁东中期花、大乌嘴、高箭中期花、四月白、祁东四月花、小花黄花丁庄大菜、猛子花、祁珍花、祁东白花菜-2、金蕾1号、高亭黄花福建红花、重庆红花、祁东白花菜-1、龙游红花菜、雪不谢苗、安民花1号143/141143/143143/139141/141141/139139/139HHC-42143137129125121福建红花、重庆红花、祁东白花菜-1、龙游红花菜、雪不谢苗安民花1号祁东中期花、猛子花、荆州花、大乌嘴、高箭中期花、祁珍花、四月白、茄子花、金蕾1号四川武平早、棒槌花祁东四月花淮阳花、沙苑金针、祁东白花菜-2、浙江仙居花丁庄大菜、渠县花、大同黄花菜、倒箭花、航天2号、小花黄花、高亭黄花143/143143/125137/129137/125129/129125/125121/121HHC-43148145142139136安民花1号祁东中期花、四川武平旱、航天2号猛子花、雪不谢苗、金蕾1号丁庄大菜、渠县花、荆州花、大乌嘴、大同黄花菜、倒箭花、高箭中期花、淮阳花、沙苑金针、祁东白花菜-2、祁东四月花、祁东财桥花、浙江仙居花、小花黄花、高亭黄花棒槌花祁珍花、四月白、茄子花福建红花、重庆红花、祁东白花菜-1、龙游红花菜148/139145/145145/142145/139142/139139/139145/136HHC-45135132祁东中期花、猛子花、高箭中期花、金蕾1号福建红花、重庆红花、祁东白花菜-1、龙游红花菜、祁东白花菜-2、安民花1号丁庄大菜、渠县花、大乌嘴、大同黄花菜、倒箭花、淮阳花、四川武平早、沙苑金针、雪不谢苗、四月白、棒槌花、茄子花、祁东四月花、祁东财桥花、浙江仙居花、航天2号、小花黄花、高亭黄花135/135135/132132/132HHC-52153150147141安民花1号猛子花、高箭中期花、四月白、祁东财桥花、金蕾1号丁庄大菜、祁东中期花、倒箭花、四川武平早、祁珍花、棒槌花、祁东白花菜-2、祁东四月花、小花黄花、高亭黄花渠县花、荆州花、大乌嘴、大同黄花菜、淮阳花、沙苑金针、茄子花、浙江仙居花、航天2号福建红花、重庆红花、祁东白花菜-1、龙游红花菜、雪不谢苗153/153153/147150/147147/147141/141HHC-53157154151148丁庄大菜、四月白、祁东白花菜-2渠县花、祁东中期花、荆州花、大乌嘴、大同黄花菜、淮阳花、四川武平早、沙苑金针、棒槌花、安民花1号、浙江仙居花、航天2号、小花黄花猛子花、祁东四月花、金蕾1号福建红花、重庆红花、祁东白花菜-1、龙游红花菜、雪不谢苗倒箭花、祁珍花、茄子花、祁东财桥花、高亭黄花157/154154/154154/151151/151154/148

续表4

引物名称等位变异/bp参照品种名称参照品种基因型数据HHC-65158156154150丁庄大菜、四月白、祁东白花菜-2渠县花、祁东中期花、荆州花、大乌嘴、大同黄花菜、倒箭花、高箭中期花、淮阳花、四川武平早、沙苑金针、棒槌花、安民花1号、浙江仙居花、航天2号、小花黄花猛子花、祁东四月花、金蕾1号福建红花、重庆红花、祁东白花菜-1、龙游红花菜、雪不谢苗倒箭花、祁珍花、茄子花、祁东财桥花、高亭黄花158/156156/156156/154154/154156/150HHC-72120118116114高箭中期花雪不谢苗丁庄大菜、福建红花、重庆红花、祁东白花菜-1、荆州花、龙游红花菜祁东中期花、茄子花、安民花1号、祁东四月花四川武平早、祁东白花菜-2、祁东财桥花渠县花、猛子花、大乌嘴、大同黄花菜、倒箭花、淮阳花、沙苑金针、棒槌花、浙江仙居花、金蕾1号、高亭黄航天2号、小花黄花120/120120/118120/116118/118118/116116/116114/114HHC-78110106102祁东四月花福建红花、重庆红花、祁东白花菜-1、龙游红花菜、雪不谢苗、安花1号丁庄大菜、渠县花、祁东中期花、猛子花、荆州花、大乌嘴、大同黄花菜、倒箭花、高箭中期花、淮阳花、四川武平旱、沙苑金针、祁珍花、四月白、棒槌花、茄子花、祁东白花菜-2、祁东财桥花、浙江仙居花、金蕾1号、航天2号、小花黄花、高亭黄花110/102106/102102/102HHC-80144138132四川武平早、小花黄花猛子花、金蕾1号、航天2号雪不谢苗、安民花1号丁庄大菜、渠县花、福建红花、重庆红花、祁东白花菜-1、荆州花、大乌嘴、大同黄花菜、倒箭花、龙游红花菜、淮阳花、沙苑金针、四月白、棒槌花、茄子花、祁东白花菜-2、祁东四月花、浙江仙居花祁东中期花、猛子花、高箭中期花、祁珍花、高亭黄花144/138144/132138/138138/132132/132HHC-90159156153150147144141祁东中期花丁庄大菜、四月白、茄子花、祁东白花菜-2、祁东财桥花渠县花、大乌嘴、大同黄花菜、淮阳花、沙苑金针、浙江仙居花航天2号、小花黄花猛子花、荆州花、倒箭花、高箭中期花、祁珍花、祁东四月花、金蕾1号、高亭黄花棒槌花福建红花、重庆红花、祁东白花菜-1、龙游红花菜、雪不谢苗四川武平早安民花1号159/156156/153156/150156/147156/144153/153153/144150/144150/141HHC-93175157151猛子花、金蕾1号丁庄大菜、渠县花、祁东中期花、荆州花、大乌嘴、大同黄花菜、倒箭花、高箭中期花、淮阳花、四川武平早、沙苑金针、祁珍花、四月白、棒槌花、茄子花、祁东白花菜-2、祁东四月花、祁东财桥花、浙江仙居花、航天2号、小花黄花、高亭黄花雪不谢苗福建红花、重庆红花、祁东白花菜-1、龙游红花菜、安民花1号、175/157157/157157/151151/151HHC-95131128125122丁庄大菜、渠县花、猛子花、大乌嘴、大同黄花菜、高箭中期花、淮阳花、四川武平早、沙苑金针、祁珍花、四月白、棒槌花、祁东白花菜-2、祁东四月花、浙江仙居花、金蕾1号、航天2号、小花黄花祁东中期花、祁东白花菜-1、荆州花、倒箭花、茄子花、祁东财桥花安民花1号、高亭黄花福建红花、重庆红花、龙游红花菜雪不谢苗祁东四月花131/131131/128131/125128/128125/122131/122HHC-99155151147福建红花、重庆红花、祁东白花菜-1、龙游红花菜茄子花、安民花1号丁庄大菜、渠县花、祁东中期花、猛子花、荆州花、大乌嘴、大同黄花菜、倒箭花、高箭中期花、淮阳花、四川武平早、沙苑金针、雪不谢苗、祁珍花、四月白、棒槌花、祁东白花菜-2、祁东四月花、祁东财桥花、浙江仙居花、金蕾1号、航天2号、小花黄花、高亭黄花155/151155/147147/147HHC-101146福建红花、重庆红花、祁东白花菜-1、龙游红花菜、雪不谢苗146/146

续表4

引物名称等位变异/bp参照品种名称参照品种基因型数据HHC-106130127121118雪不谢苗安民花1号、航天2号丁庄大菜、猛子花、荆州花、大乌嘴、倒剪花、四川武平旱、棒槌花、祁东白花菜-2、祁东四月花、祁东财桥花渠县花、祁东中期花、大同黄花菜、高箭中期花、淮阳花、沙苑金针、祁珍花、四月白、茄子花、浙江仙居花、高亭黄花福建红花、重庆红花、祁东白花菜-1、龙游红花菜130/130127/127127/121121/121121/118HHC-10910310097猛子花、大乌嘴、雪不谢苗、金蕾1号丁庄大菜、渠县花、祁东中期花、大同黄花菜、淮阳花、四川武平旱、沙苑金针、祁珍花、四月白、棒槌花、茄子花、祁东白花菜-2、祁东四月花、浙江仙居花、航天2号、小花黄花、高亭黄花福建红花、重庆红花、祁东白花菜-1、高箭中期花、龙游红花菜、安民花1号、祁东财桥花103/103100/100100/97

参考文献

[1] 朱自学.黄花菜良种选育与繁殖技术[J].农业与技术,2007,27(2):132-134.

ZHU Z X.Breeding and propagation techniques of fine varieties of daylily[J].Agriculture andTechnology,2007,27(2):132-134.

[2] 邓放明,尹华,李精华,等.黄花菜应用研究现状与产业化开发对策[J].湖南农业大学学报(自然科学版),2003,29 (6):529-532.

DENG F M, YIN H, LI J H, et al.On latest application and countermeasure for industrizationexploitation of daylily flower[J].Journal of Hunan Agricultural University(Natural Sciences),2003,29 (6):529-532.

[3] DU B J, TANG X S, LIU F, et al.Antidepressant-like effects of the hydroalcoholic extracts of Hemerocallis citrina and its potential active components[J].BMC Complementary and Alternative Medicine,2014,14:326-336.

[4] 翟俊乐,田欢,李孟秋,等.黄花菜抗抑郁作用有效成分的筛选[J].中国食品添加剂,2015(10):93-97.

ZHAI J L,TIAN H, LI M Q, et al.Screen of active anti-depression ingredients from daylily[J].China Food Additives, 2015(10):93-96.

[5] 毛建兰.黄花菜的营养价值及加工技术综述[J].安徽农业科学,2008(3):1 197-1 198.

MAO JL.Summary of nutritional value and processing technology of daylily[J].Journal of Anhui Agriculture Science, 2008(3):1 197-1 198.

[6] 朱岩芳.作物品种分子标记鉴定及指纹图谱构建研究[D].杭州:浙江大学,2013.

ZHU Y F.Research on cultivar identification and DNA fingerprinting of crops based on molecular markers[D].Hangzhou:Zhejiang University, 2013.

[7] 王凤格,杨扬,易红梅,等.中国玉米审定品种标准SSR 指纹库的构建[J].中国农业科学,2017,50(1):1-14.

WANG F G, YANG Y, YI H M, et al.Construction of an SSR-based standard fingerprint database for corn variety authorized in China[J].Scientia Agricultura Sinica, 2017, 50(1):1-14.

[8] 李志江,牛江帅,李忍,等.黑龙江省种植大豆品种遗传多样性分析及与性状关联SSR标记筛选[J].中国粮油学报,2021,36(4):37-44.

LI Z J, NIU J S, LI REN, et al.Genetic diversity analysis of soybean varieties and selection of SSR markers associated with soybeantraits in Heilongjiang Province[J].Journal of the Chinese Cereals and Oils Association,2021,36(4):37-44.

[9] 代畅.利用EST-SSR分子标记构建小麦遗传图谱[D].南充:西华师范大学,2018.

DAI C.Construction of wheat genetic map using EST-SSR molecular markers[D].Nanchong:China West Normal University, 2018.

[10] 王瑞,张福耀,程庆军,等.利用SSR 荧光标记构建20个高粱品种指纹图谱[J].作物学报, 2015, 41(4):658-665.

WANG R, ZHANG F Y, CHENG Q J, et al.Establishment of 20 sorghum hybrids fingerprints using SSR fluorescent marker[J].Acta AgronomicaSinica, 2015, 41(4):658-665.

[11] 何俊平,张书芬,王建平,等.甘蓝型油菜育种亲本材料遗传多样性的SSR分析[J].西南农业学报,2015,28(6):2 374-2 380.

HE J P, ZHANG S F, WANG J P, et al.Genetic diversity of breeding parents of Brassica napus by SSR markers[J].Southwest China Journal of Agricultural Sciences, 2015,28(6):2 374-2 380.

[12] 陈士林.中国药典中药材DNA条形码标准序列[M].北京:科学出版社,2015.

CHEN S L.DNA Barcode Standard Sequence of Chinese Herbal Medicines in Chinese Pharmacopoeia[M].Beijing:Science Press,2015.

[13] 朱华芳,罗玉兰,胡永红,等.萱草属部分种和园艺品种的SSR多态性分析[J].上海交通大学学报(农业科学版),2009,27(2):143-148.

ZHU H F, LUO Y L, HU Y H, et al.Analysis of a section of Hemerocallis species and varieties by SSR- PCR[J].Journal of Shanghai Jiaotong University:Agricultural Science, 2009,27(2):143-148.

[14] 任阳.萱草属植物的形态学性状与分子标记研究[D].长春:吉林农业大学,2017.

REN Y.The morphological characteristics and molecular basis of plant development of Hemerocallis[D].ChangChun:Jilin Agricultural University,2017.

[15] 董美超,岳建强,李进学,等.基于荧光SSR标记的鳄梨种质资源遗传多样性[J].分子植物育种,2020,18(16):5 403-5 410.

DONG M C,YUE J Q,LI J X, et al.Geneticdiversity analysis of germplasm resources of Perseaamericana mill.based on SSR fluorescentmarker[J].Molecular Plant Breeding, 2020,18(16):5 403-5 410.

[16] 吴涛,陈少瑜,肖良俊,等.基于SSR标记的云南省核桃种质资源遗传多样性及核心种质构建[J].植物遗传资源学报,2020,21(3):767-774.

WU T,CHEN S N,XIAO L J,et al.Genetic diversity analysis and core collection construction of walnut germplasm in Yunnan Province[J].Journal of Plant Genatic Resources, 2020,21(3):767-774.

[17] YARED B, HUSSEIN S, MARK L,et al.Genetic diversity and population structure of bread wheat genotypes determined via phenotypic and SSR marker analyses under drought-stress conditions[J].Journal of Crop Improvement,2021,35(3):302-325.

[18] WANG L Y,LI S S,WANG T Y,et al.Genomic SSR and EST‐SSR markers for phylogenetic and pedigree reconstructions —A comparison in sea buckthorn[J].Plant Breeding,2021,140(1):167-183.

[19] CHEN J H,DONG S J,ZHANG X,et al.Genetic diversity of Prunussibirica L.superior accessions based on the SSR markers developed using restriction-site associated DNA sequencing[J].Genetic Resources and Crop Evolution,2020,68(2):615-628.

[20] 马文瑞,严密,刘业伟,等.酿酒葡萄品种遗传多样性的简单序列重复分析和分子身份证的建立[J].食品与发酵工业,2018,44(9):232-238.

MA W R,YAN M,LIU Y W,etal.Analysis of genetic relationship of Vitisvinifera resources by SSR markers[J].Food and Fermentation Industries,2018,44(9):232-238.

[21] 李益,马先锋,唐浩,等.柑橘品种鉴定的SSR标记开发和指纹图谱库构建[J].中国农业科学,2018,51(15):149-159.

LI Y,MA X F,TANG H,et al.SSR Markers screening for identification of citrus cultivar and construction of DNA fingerprintinglibrary[J].ScientiaAgricultura Sinica,2018, 51(15):149-159.

[22] 郑家祯,李和平,赖正锋,等.国内菜用黄花菜种质资源遗传多样性分析[J].福建农业学报,2018,33(10):1 030-1 038.

ZHENG J Z, LI H P, LAI Z F, et al.Genetic diversity of vegetable daylily germplasms in China[J].Fujian Journal of Agricultural Sciences, 2018,33(10):1 030-1 038.

[23] 李森,史青青,侯非凡,等.萱草属种质资源主要观赏性状的多样性分析[J].山西农业大学学报(自然科学版),2016,36(9): 619-627.

LI S, SHI Q Q, HOU F F, et al.Diversity of the main ornamental traits for Henerocallis spp.accessions[J].Shanxi Agricultral University (Natural Science Edition), 2016,36(9):619-627.

[24] 黎海利,董丽.萱草部分野生种和栽培品种亲缘关系的AFLP分析[J].园艺学报,2009,36 (8):1 203-1 208.

LI H L, DONG L.Relationship of Hemerocallis spp.Wild species and cultivars by AFLPmarker[J].Acta Horticulture, 2009,36(8):1 203-1 208.

The development of SSR markers and construction of fingerprint library for edible daylily

LIU Xinxing,OU Qiaoming*,LUO Junjie,LI Zhongwang,SHI Youtai,CUI Wenjuan

(Biotechnology Institute, Gansu Academy of Agricultural Sciences, Lanzhou 730070, China)

ABSTRACT The SSR markers were developed by transcriptome sequencing technology to provide technical support for the identification, evaluation and genetic diversity analysis of edible daylily germplasm resources. The results showed that a total of 104 177 unigenes containing 79 106 833 bpin sequence, with an average length of 759 bp was obtained. The Unigene was searched using MISA software, and a total of 15 890 SSRs were detected which including single nucleotides, dinucleotides and trinucleotides accounting for a relatively high proportion of 18.57%, 27.23% and 44.26%, respectively. A total of 114 pairs of SSR primers were used for the screening, and finally, 20 pairs of primers were screened out by polyacrylamide gel electrophoresis. A total of 88 allelic variants were detected in 43 daylily germplasms, and the number of allelic variants at each locus was 3-8. The polymorphism information content (PIC) was ranged from 0.23 to 0.77. Based on the fragment size of the amplified product, the allelic variation of each SSR locus and the corresponding reference variety were determined, and a DNA fingerprint identification technology system was constructed. The SSR markers developed in this study are rich in polymorphism.

Key words daylily; SSR marker; fingerprint

DOI:10.13995/j.cnki.11-1802/ts.028102

引用格式:刘新星,欧巧明,罗俊杰,等.食用黄花菜SSR标记开发及指纹图谱库构建[J].食品与发酵工业,2022,48(12):188-195.LIU Xinxing,OU Qiaoming,LUO Junjie, et al.The development of SSR markers and construction of fingerprint library for edible daylily[J].Food and Fermentation Industries,2022,48(12):188-195.

第一作者:助理研究员(欧巧明为通信作者,E-mail:ouqioaming@163.com)

基金项目:甘肃省农业科学院中青年基金项目(2019GAAS39;2019GAAS43);甘肃省农业科学院科技成果转化项目(2019GAAS-CGZH03)

收稿日期:2021-05-26,改回日期:2021-08-17