• 《中国科学引文数据库(CSCD)》来源期刊
  • 中国科技期刊引证报告(核心版)期刊
  • 《中文核心期刊要目总览》核心期刊
  • RCCSE中国核心学术期刊

基于mtDNA和cpDNA序列的甘薯栽培种及近缘野生种分析

王崇, 王连军, 田小海, 雷剑, 柴沙沙, 焦春海, 杨新笋

王崇, 王连军, 田小海, 等. 基于mtDNA和cpDNA序列的甘薯栽培种及近缘野生种分析[J]. 华南农业大学学报, 2021, 42(4): 25-32. DOI: 10.7671/j.issn.1001-411X.202011026
引用本文: 王崇, 王连军, 田小海, 等. 基于mtDNA和cpDNA序列的甘薯栽培种及近缘野生种分析[J]. 华南农业大学学报, 2021, 42(4): 25-32. DOI: 10.7671/j.issn.1001-411X.202011026
WANG Chong, WANG Lianjun, TIAN Xiaohai, et al. Analyses of Ipomoea batatas cultivated species and wild relatives based on mtDNA and cpDNA sequences[J]. Journal of South China Agricultural University, 2021, 42(4): 25-32. DOI: 10.7671/j.issn.1001-411X.202011026
Citation: WANG Chong, WANG Lianjun, TIAN Xiaohai, et al. Analyses of Ipomoea batatas cultivated species and wild relatives based on mtDNA and cpDNA sequences[J]. Journal of South China Agricultural University, 2021, 42(4): 25-32. DOI: 10.7671/j.issn.1001-411X.202011026

基于mtDNA和cpDNA序列的甘薯栽培种及近缘野生种分析

基金项目: 国家重点研发计划(2019YFD1001300,2019YFD1001304,2019YFD1001305);粮食作物种质创新与遗传改良湖北省重点实验室开放课题(国家现代甘薯产业技术体系建设项目)(CARS-11-C-15);湖北省技术创新专项(重大项目)(2017ABA149);湖北省农业科学院特色学科、湖北省农业科技创新中心资助项目(2007-620-001-03);湖北省农业科学院青年拔尖人才培养计划
详细信息
    作者简介:

    王崇(1994—),男,博士研究生,E-mail: wangchong199409@163.com

    王连军(1980—),男,副研究员,博士,E-mail: wanglianjun10@163.com;†表示同等贡献

    通讯作者:

    焦春海(1962—),男,研究员,博士,E-mail: jiaoch@hotmail.com

    杨新笋(1967—),男,研究员,博士,E-mail: yangxins013@163.com

    †表示同等贡献

  • 中图分类号: S531

Analyses of Ipomoea batatas cultivated species and wild relatives based on mtDNA and cpDNA sequences

  • 摘要:
    目的 

    利用线粒体DNA(Mitochondrial DNA,mtDNA)序列和叶绿体DNA(Chloroplast DNA,cpDNA)matK序列对甘薯Ipomoea batatas栽培种及近缘野生种进行分子鉴定和亲缘关系分析,为甘薯栽培种和近缘野生种的种质鉴定、保护及开发利用提供理论依据。

    方法 

    以3个甘薯栽培种及8个近缘野生种为材料,采用CTAB法提取基因组DNA,通过PCR扩增mtDNA序列和cpDNA matK序列,使用DnaSP 6.0对序列进行核苷酸多态性、单倍型多样性等特征分析,并基于邻接法构建3个甘薯栽培种及8个近缘野生种的系统发育进化树。

    结果 

    5个mtDNA序列和cpDNA matK序列经测序、比对、拼接后,长度为6 713 bp,GC占比在47.79%~48.31%之间。合并序列的单倍型数量、核苷酸多态性、变异位点数量、单一突变位点数量、简约信息位点数量和插入/缺失位点数量分别为9、0.003 25、69、39、30和111。中性检验显示,合并序列差异不显著(P>0.10),遵循中性进化模型。3个甘薯栽培种及8个近缘野生种间的遗传距离在0.000 00~0.005 84之间,平均遗传距离0.003 26,遗传多样性较低;按照亲缘关系被分为2大类,各大类内部亲缘关系较近。

    结论 

    本研究采用的序列可对甘薯栽培种和近缘野生种进行准确的鉴定区分,为甘薯近缘野生种的进化和利用提供参考和理论指导。

    Abstract:
    Objective 

    To conduct molecular identification and genetic relationship analysis of Ipomoea batatas cultivated species and wild relatives based on mitochondrial DNA (mtDNA) and chloroplast DNA (cpDNA) matK sequences, and provide theoretical bases for germplasm identification, protection, development and utilization.

    Method 

    Three cultivated species and eight wild relatives were used as materials, from which total DNA was extracted by the CTAB method. Their mtDNA and cpDNA matK sequences were amplified by PCR. DnaSP 6.0 was used to analyze nucleotide diversity, haplotype diversity and other characteristics. The phylogenetic tree of three cultivated species and eight wild relatives was constructed based on the neighbor-joining method.

    Result 

    The length of five mtDNA regions and one cpDNA region was 6 713 bp after sequencing, alignment and splicing, the GC proportion was 47.79%−48.31%, and the haplotype number, nucleotide diversity, variable site number, singleton variable site number, parsimony informative site number, insertion/deletion site number were 9, 0.003 25, 69, 39, 30, 111, respectively. The neutrality test showed there was no significant difference between Tajima’sD values at the level of P>0.10, which indicated that variation of those regions followed neutral theory of molecular evolution. The genetic distances among three cultivated species and eight wild relatives ranged from 0.000 00 to 0.005 84, with an average genetic distance of 0.003 26, indicating low genetic diversity. The phylogenetic tree divided the 11 species into two categories with close genetic relationship within the category.

    Conclusion 

    The sequences used in this study can accurately identify I. batatas cultivated species and wild relatives, and provide references and theoretical guidance for the evolution and utilization of I. batatas wild relatives.

  • 甘薯Ipomoea batatas是旋花科Convolvulaceae番薯薯Ipomoea双子叶植物,主要起源于南美洲等热带、亚热带地区,是世界第七大粮食作物、中国第五大粮食作物[1]。甘薯富含蛋白质、脂肪、多糖、维生素、胡萝卜素、花青素、氨基酸、钙、磷、钾等多种物质,营养价值丰富[2-3]。甘薯抗旱性强,耐贫瘠能力强,具有极强的适应性,在中国具有悠久的栽培历史[4]。甘薯是无性繁殖作物,染色体数目多,遗传背景复杂,在遗传上高度杂合,中国的甘薯栽培种94%以上具有‘南瑞苕’和‘胜利百号’的遗传信息,造成中国甘薯遗传背景狭窄,严重制约甘薯品种的遗传改良[5]。甘薯近缘野生种种质资源丰富,在世界范围内分布有600~700种,具有丰富的遗传多样性和潜在的基因资源,如高淀粉含量、抗病虫、抗病毒、抗旱等优良特质,可为改良甘薯品种和拓宽甘薯基因库提供重要的基因资源[6-7]

    伴随现代生物技术不断发展、更新,DNA分子标记技术在植物研究领域中广泛应用于植物种质鉴定、遗传多样性分析、指纹图谱构建、亲缘关系鉴定和遗传图谱构建等[8]。Yang等[9]使用SSR标记分析了380个甘薯品种的遗传多样性和遗传结构;苏一钧等[10]筛选出30对多态性丰富的SSR引物,对303份甘薯地方材料进行了遗传多样性和群体结构分析,发现303份甘薯材料间亲缘关系较近,遗传多样性丰富;季志仙等[11]选用6个ISSR标记对17份甘薯品种进行分析,显示同类型品种间的遗传相似性较高,亲缘关系密切。SSR和ISSR等分子标记技术主要是基于植物基因组DNA开发,在植物亲缘鉴定和遗传背景研究中具有极好的应用价值。随着现代生物技术的发展,基于线粒体DNA(Mitochondrial DNA,mtDNA)和叶绿体DNA(Chloroplast DNA,cpDNA)开发的分子标记技术也越来越广泛地应用于植物遗传分析[12-13]

    mtDNA具有细胞质遗传特性,一些mtDNA参与编码具有重要生物学功能的蛋白或亚基,如核糖体和生物氧化链上某些重要酶的部分亚基,还与细胞生命周期、真核细胞抗药性有关[14]。在植物细胞中,mtDNA参与调控子叶大小、叶片形状、植物细胞质雄性不育,还影响拟南芥、玉米和黄瓜等作物的叶绿体结构以及光合作用[15-16]。mtDNA具有进化速度快、缺少重组以及严格母系遗传等特点,且线粒体基因组较小、易于测序,已经成为研究植物遗传进化和系统发育的理想工具[17]。cpDNA具有单倍型、母系遗传、缺乏杂合细胞质、不发生等位基因重组等特性,非编码区不受选择影响,核酸位点表现出丰富的多态性和变异,与mtDNA序列广泛用于研究植物的系统发育、种间或种内的亲缘关系鉴定等研究[18]。Godbout等[19]采用nad3、nad5等mtDNA序列对北美短叶松进行了系统发育研究,结果表明mtDNA可反映北美短叶松的地理关系;Hu等[20]使用mtDNA非编码区序列分析薯蓣属的遗传背景,结果表明线粒体基因重组时存在丰富的变异位点,可用于鉴别区分不同薯蓣品种;马丽等[21]使用matK基因序列对59个石榴品种进行了鉴定和分类,结果表明matK可用于石榴的系统发育分析。

    本研究通过对3个甘薯栽培种及8个近缘野生种进行mtDNA和cpDNA matK序列分析,综合分析3个甘薯栽培种及8个近缘野生种的亲缘关系和遗传多样性,为甘薯近缘野生种的保护及利用提供参考,为甘薯遗传育种提供一定的理论依据。

    试验所用材料为3个甘薯栽培种(‘浙紫薯3号’‘鄂薯6号’‘皖薯7号’)及8个近缘野生种(何鲁牵牛I. holubii、‘黄茑萝’I. hederifoliavar. lutea、空心菜I. aquatica Forsk、‘月光花’I. alba、‘黄色朝颜’I. obscura Keniak、变色牵牛I. indica、旋转牵牛I. digitata huge caudex、树牵牛I. carnea),种植于湖北省农业科学院粮食作物研究所。长势旺盛时进行样品采集。

    采取试验样品新鲜、幼嫩叶片后,参考Yang等[9]的方法,使用改良CTAB法提取甘薯栽培种及近缘野生种的基因组DNA。使用分光光度计测量DNA的浓度,用15 g/L琼脂糖凝胶电泳检测DNA的质量。将DNA质量浓度稀释到70~80 ng/µL,在−20 ℃冰箱中保存,备用。

    本研究选用6对mtDNA和1对cpDNA通用引物(表1),引物由武汉天一辉远生物科技有限公司合成。将7对引物分别在试验材料中进行扩增。PCR反应扩增体系:10× PAGE Buffer(Mg2+)2.5 µL,dNTPs(10 mmol/L)2.0 µL,正向引物和反向引物各1.0 µL,基因组DNA(70~80 ng/uL)1.0 µL,EasyTaq DNA Polymerase 0.5 µL,ddH2O补至总体积为25 µL。PCR扩增程序:94 ℃预变性5 min;94 ℃变性30 s,56~57 ℃退火30 s,72 ℃延伸1 min,35个循环;72 ℃延伸5 min;在4 ℃条件下保存5 min。使用15 g/L的琼脂糖凝胶电泳检测PCR扩增结果,以单一、清晰、明亮条带作为检测标准,将检测合格的扩增产物送到武汉天一辉远生物科技有限公司进行双向测序。

    表  1  试验采用引物信息
    Table  1.  The information of primers used in this experiment
    类别
    Species
    基因序列
    Gene sequence
    引物
    Primer
    片段大小/bp
    Fragment size
    退火温度/℃
    Annealing temperature
    参考文献
    Reference
    mtDNA ccb203 F: ASGTTCTACGGACCGATGCC 500 57 [22]
    R: CACGGGGAGGGAGCRGGCGA
    ccb256 F: GGAAGTTAGCAAAGTTAGAC 520 57 [22]
    R: TTGTTCTTAACAGCGATGGC
    nad2/1-2 F: TTTTCTTCCTCATTCTKATTT 1200 57 [22]
    R: CCACTCTATTGTCCACTTCTA
    nad2/4-5 F: TTCATATAGAATCCATGTCC 1800 57 [22]
    R: CTATTTGTTCTTCGCCGCTT
    nad5/4-5 F: CCAATTTTTGGGCCAATTCC 1400 57 [23]
    R: CATTGCAAAGGCATAATGAT
    nad7/1-2 F: ACCTCAACATCCTGCTGCTC 1200 57 [23]
    R: CGATCAGAATAAGGTAAAGC
    cpDNA matK F: CGTACAGTACTTTTGTGTTTACGAG 1500 57 [24]
    R: ACCCAGTCCATCTGGAAATCTTGGTTC
    下载: 导出CSV 
    | 显示表格

    用软件DNAstar对测序结果进行手动校对,删除两端冗余碱基,将正向序列和反向序列进行拼接,使用软件MEGA X[25]对拼接好的序列进行比对,使用软件DnaSP 6.0[26]对序列进行多态性(π)和单倍型(Hd)等序列特征分析。使用MEGA X软件采用邻接距离矩阵法(Neighbor-joining distance matrix method)构建系统发育进化树,使用自展法(Bootstrap)进行可信度检测,重复1000次。使用软件NetWork的中介邻接网络(Median-joining network,MJ)算法构建单倍型关联图。

    对6个mtDNA序列和cpDNA matK序列的测序结果进行多态性分析,mtDNA序列ccb256在3个甘薯栽培种及8个近缘野生种中无碱基差异,其余6个序列表现出多态性。对表现出多态性的序列进行分析,将ccb203、nad2/1-2、nad2/4-5、nad5/4-5、nad7/1-2和matK在11个材料中的扩增序列进行对位排列后,上述序列的长度分别为494、1 247、1 503、1 665、964和840 bp。6个序列扩增片段GC碱基占比分别为46.52%~47.17%、53.90%~54.69%、48.02%~48.54%、45.93%~46.50%、56.39%~57.23%和32.90%~34.09%。在5个线粒体序列中,序列nad5/4-5的核酸多态性最高(0.003 95),含有20个变异位点、10个单一突变位点、10个简约信息位点、6个插入/缺失位点。ccb203、nad2/1-2、nad2/4-5和nad7/1-2序列的核苷酸多态性、变异位点数、单一突变位点数、简约信息位点数和插入/缺失位点数分别为0.003 20、0.000 17、0.003 43、0.001 86,5、1、12、4,2、1、3、0,3、0、9、4,6、11、63、24。matK序列的核苷酸多态性、变异位点数、单一突变位点数比mtDNA序列的高,简约信息位点数和插入/缺失位点数比mtDNA序列的低。6个序列合并后,序列长度为6 713 bp,共有变异位点69个,其中单一突变位点39个、简约信息位点30个,插入/缺失位点111个(表2)。

    表  2  试验材料的序列多态性信息
    Table  2.  Sequence polymorphism information of experiment materials
    基因序列
    Gene
    sequence
    长度/
    bp
    Length
    GC占
    比/%
    GC
    proportion
    核苷酸多态性
    Nucleotide
    diversity
    变异位点数
    Variable
    site
    number
    单一突变位点数
    Singleton
    variable
    site number
    简约信息位点数
    Parsimony
    informative
    site number
    插入/缺失位点数
    Insertion/
    deletion
    site number
    ccb203 494 46.52~47.17 0.003 20 5 2 3 6
    nad2/1-2 1 247 53.90~54.69 0.000 17 1 1 0 11
    nad2/4-5 1 503 48.02~48.54 0.003 43 12 3 9 63
    nad5/4-5 1 665 45.93~46.50 0.003 95 20 10 10 6
    nad7/1-2 964 56.39~57.23 0.001 86 4 0 4 24
    matK 840 32.90~34.09 0.007 11 27 23 4 1
    合并 Total 6 713 47.79~48.31 0.003 25 69 39 30 111
    下载: 导出CSV 
    | 显示表格

    ccb203、nad2/1-2、nad2/4-5、nad5/4-5、nad7/1-2和matK序列中,单倍型数量最多的是cpDNA序列matK,9个;ccb203、nad2/1-2、nad2/4-5、nad5/4-5和nad7/1-2的单倍型数量分别是5、2、8、8和5。6个序列的单倍型多样性分别为0.818、0.182、0.927、0.927、0.818和0.945,matK区域的单倍型多样性最高。单倍型多样性方差和单倍型多样性标准差最大的区域是nad2/1-2,分别为0.020 61和0.144。6个序列合并后,单倍型数量、单倍型多样性、单倍型多样性方差和单倍型多样性标准差分别是9、0.945、0.004 34和0.066(表3)。

    表  3  试验材料的单倍型多样性
    Table  3.  Haplotype diversity of experiment materials
    基因序列
    Gene
    sequence
    单倍型数量
    Haplotype
    number
    单倍型多样性
    Haplotype
    diversity
    单倍型多样性方差
    Variance of
    haplotype diversity
    单倍型多样性标准差
    Standard deviation of
    haplotype diversity
    ccb203 5 0.818 0.006 82 0.083
    nad2/1-2 2 0.182 0.020 61 0.144
    nad2/4-5 8 0.927 0.004 42 0.066
    nad5/4-5 8 0.927 0.004 42 0.066
    nad7/1-2 5 0.818 0.006 82 0.083
    matK 9 0.945 0.004 34 0.066
    合并 Total 9 0.945 0.004 34 0.066
    下载: 导出CSV 
    | 显示表格

    对6个序列进行Tajima’s D测验和Fu and Li’s D*/F*测验,在5个mtDNA序列中,序列nad7/1-2的Tajima’s D值、Fu and Li’s D*值和Fu and Li’s F*值最大,序列nad2/1-2的Tajima’s D值、Fu and Li’s D*值和Fu and Li’s F*值最小,且5个mtDNA序列的Tajima’s D值、Fu and Li’s D*值和Fu and Li’s F*值均差异不显著(P>0.10),符合中性进化模型。cpDNAmatK的Tajima’s D值、Fu and Li’s D*值和Fu and Li’s F*值分别为−1.627 84、−2.009 96和−2.168 40,在0.05<P<0.10的情况下差异不显著,符合中性进化模型。6个序列合并后,Tajima’sD值、Fu and Li’s D*值和Fu and Li’s F*值分别为−0.639 26、−0.796 77和−0.858 74,差异不显著(P>0.10),符合中性进化模型(表4)。

    对3个甘薯栽培种及8个近缘野生种进行遗传多样性分析,11份试验材料间的遗传距离均在0.000 00~0.005 84之间,‘月光花’与‘黄色朝颜’2个品种间的遗传距离最大,为0.005 84,3个栽培种间的遗传距离为0,3个栽培种及8个近缘野生种间的平均遗传距离为0.003 26 (表5)。

    表  4  Tajima’s D测验和Fu and Li’s D*/F*测验
    Table  4.  Tajima’s D test and Fu and Li’s D*/F* test
    基因序列 Gene sequence Tajima’s D Fu and Li’s D* Fu and Li’s F* P
    ccb203 −0.321 97 −0.083 18 −0.161 10 > 0.10
    nad2/1-2 −1.128 50 −1.289 46 −1.399 19 > 0.10
    nad2/4-5 0.500 00 0.574 40 0.628 71 > 0.10
    nad5/4-5 −0.356 59 −0.514 26 −0.537 17 > 0.10
    nad7/1-2 1.018 28 1.214 66 1.312 69 > 0.10
    matK −1.627 84 −2.009 96 −2.168 40 > 0.05, < 0.10
    合并Total −0.639 26 −0.796 77 −0.858 74 > 0.10
    下载: 导出CSV 
    | 显示表格
    表  5  11份试验材料的遗传距离1)
    Table  5.  Genetic distance among 11 experimental materials
    样品编号
    Sample No.
    1 2 3 4 5 6 7 8 9 10
    1
    2 0.003 47
    3 0.002 99 0.003 63
    4 0.003 78 0.001 42 0.004 42
    5 0.002 99 0.005 21 0.004 73 0.005 84
    6 0.003 15 0.001 10 0.003 78 0.000 94 0.004 89
    7 0.002 52 0.003 15 0.003 63 0.003 78 0.004 26 0.003 15
    8 0.001 89 0.003 15 0.002 68 0.004 10 0.003 63 0.003 47 0.002 52
    9 0.003 47 0.003 63 0.003 15 0.003 78 0.004 89 0.003 15 0.003 15 0.003 15
    10 0.003 47 0.003 63 0.003 15 0.003 78 0.004 89 0.003 15 0.003 15 0.003 15 0.000 00
    11 0.003 47 0.003 63 0.003 15 0.003 78 0.004 89 0.003 15 0.003 15 0.003 15 0.000 00 0.000 00
     1) 1:何鲁牵牛,2:‘黄茑萝’,3:空心菜,4:‘月光花’,5:‘黄色朝颜’,6:变色牵牛,7:旋转牵牛,8:树牵牛,9:‘浙紫薯3号’,10:‘鄂薯6号’,11:‘皖薯7号’
     1) 1: Ipomoea holubii, 2: Ipomoea hederifolia var. lutea, 3: Ipomoea aquatica Forsk, 4: Ipomoea alba, 5: Ipomoea obscura Keniak, 6: Ipomoea indica, 7: Ipomoea digitata huge caudex, 8: Ipomoea carnea, 9: Ipomoea batatas ‘Zhezishu 3’, 10: Ipomoea batatas ‘Eshu 6’, 11: Ipomoea batatas ‘Wanshu 7’
    下载: 导出CSV 
    | 显示表格

    基于邻接法构建甘薯栽培种及近缘野生种的系统发育进化树,可将3个甘薯栽培种及8个近缘野生种分为2大类(Ⅰ和Ⅱ) (图1),第一大类包括7个近缘野生种,第二大类包括3个栽培种和1个近缘野生种。第一大类又可分为2类,何鲁牵牛、‘黄色朝颜’、树牵牛和旋转牵牛被归为一类,‘黄茑萝’‘月光花’和变色牵牛被归为一类;第二大类可分为2类,第一类为空心菜,第二类为‘浙紫薯3号’‘鄂薯6号’和‘皖薯7号’(图1)。

    图  1  邻接法构建系统发育进化树
    图中数字表示自展值百分比/%
    Figure  1.  The phylogenetic tree constructed by neighbor-joining method
    The number in the figure indicates the bootstrap value percentage/%

    使用软件NetWork中介邻接网络算法对3个甘薯栽培种及8个近缘野生种的单倍型绘制网络图,11份材料表现出9种单倍型,图中无明显主体单倍型(图2)。单倍型H_1~H_8分别对应何鲁牵牛、‘黄茑萝’、空心菜、‘月光花’‘黄色朝颜’、变色牵牛、旋转牵牛和树牵牛,单倍型H_9对应‘紫薯3号’‘鄂薯6号’和‘皖薯7号’。单倍型之间呈线状排列,无主体单倍型,单倍型H_1与H_5有共同节点,单倍型H_3与H_8有共同节点,单倍型H_4与H_6有共同节点。中介邻接网络算法分析结果与邻接聚类分析结果一致,甘薯栽培种与野生种之间存在界限。

    图  2  单倍型中介邻接网络图
    H_1~ H_9:单倍型,mv:可能存在的原始单倍型
    Figure  2.  Median-joining network for haplotypes
    H_1−H_9: Haplotype, mv: The possible original haplotype

    甘薯近缘野生种是重要的育种资源,基于mtDNA和cpDNA序列特征对甘薯近缘野生种进行分析,可为甘薯育种提供新的资源。序列特征分析在甘薯及近缘野生种分析中被广泛应用,刘峥等[27]以及俞立璇等[28]使用ITS序列对甘薯栽培种和近缘野生种进行序列特征分析,显示ITS序列可用于旋花科的系统演化关系分析,为甘薯野生种的利用提供了理论依据。matK序列是叶绿体中进化速率较快的基因之一,存在一定程度的分化,由于叶绿体基因组较为保守,其进化速率远低于核基因的进化速率[29]。低进化速率限制了matK及其他叶绿体基因在低级分类层面(如属、亚属)中的应用,matK序列不能应用于甘薯属内间的亲缘关系研究。线粒体基因在系统学研究中应用较为广泛,如nad1基因应用于探讨稻族Oryzeae的系统发育[30]nad5序列应用于秦艽组植物的分子系统发育分析及物种鉴定[31]nad7/1-2序列用于分析突尼斯柑橘砧木种质资源的遗传多样性[32]。与核基因相比,线粒体基因如nad1序列在较高等级类群的系统发育重建中有更高的信息量,变异水平更适于在高等级层面的系统发育重建。相较于核基因和叶绿体基因,线粒体基因变异速率低,结构不稳定,但筛选出合适的线粒体序列片段仍能为植物的系统学研究和亲缘关系鉴定提供新的证据。在植物系统学研究中,采用mtDNA对甘薯近缘野生种进行分析的报道较少,通过mtDNA和cpDNA对甘薯近缘野生种进行分析,可为研究旋花科的系统演化提供新的方法和指导。

    本研究选用6个mtDNA序列和1个cpDNA序列的引物在3个甘薯栽培种及8个近缘野生种中进行扩增,结果显示,在11个材料中,mtDNA序列ccb256无多态性,表现出多态性的序列占比为85.7%,表明基于mtDNA和cpDNA开发的通用引物在甘薯栽培种及近缘野生种的多态性分析中具有良好的多态性。具有多态性的5个mtDNA序列中,nad2/1-2和nad7/1-2序列中的GC占比大于AT占比,ccb203、nad2/4-5和nad5/4-5序列中的GC占比小于AT占比,表明nad2/1-2和nad7/1-2序列更为稳定。matK序列的GC占比小于AT占比,matK序列的核苷酸多样性高于mtDNA序列的,与Yang等[33]基于核基因组和叶绿体基因组分析中国橡树的试验结果基本保持一致。对mtDNA序列的变异位点进行分析,产生的变异位点以单一突变位点和插入/缺失突变为主,序列nad5/4-5的单一突变位点最多(10个),序列nad2/4-5插入/缺失位点在6个序列中最多(63个)。与mtDNA序列相比,cpDNA序列matK的单一突变位点(23个)比mtDNA序列的多,而插入/缺失位点(1个)远少于mtDNA序列的。对6个序列的单倍型多样性进行分析,mtDNA序列中nad2/1-2序列的单倍型数量和单倍型多样性最小(2、0.182)。matK序列的单倍型数量(9)和单倍型多样性(0.945)比5个mtDNA序列的高,表明相较于mtDNA序列,cpDNA序列反映的序列特征信息更为丰富。对甘薯近缘野生种进行分析时,mtDNA和cpDNA序列结合获得的结果更为准确。

    中性检验(Neutrality test)是通过统计学的检验方法分析DNA数据,进而检验自然选择[34]。本研究采用Tajima’s D和Fu and Li’s D*/F*检验,结果显示基于mtDNA序列分析的Tajima’s D值、Fu and Li’s D*值和Fu and Li’s F*值均差异不显著(P>0.10),符合中性进化模型;matK序列的Tajima’s D值、Fu and Li’s D*值和Fu and Li’s F*在0.05<P<0.10时差异不显著,符合中性进化模型;当6个序列合并后,合并序列差异不显著(P>0.10),符合中性进化模型。中性检验的结果表明,甘薯近缘野生种在遗传进化时表现稳定,符合中性进化模型,没有出现种群扩张的情况,单倍型网络无中心再度证明甘薯近缘野生种在遗传进化时无种群扩张的情况。本研究3个甘薯栽培种及8个近缘野生种间的遗传距离在0.000 00~0.005 84之间,平均遗传距离为0.003 26,表明3个甘薯栽培种及8个近缘野生种的遗传距离小,彼此间亲缘关系近,遗传多样性低,如‘浙紫薯3号’‘鄂薯6号’和‘皖薯7号’之间的遗传距离为0,亲缘关系密切。使用5个mtDNA序列和1个cpDNAmatK序列对甘薯栽培种和近缘野生种进行分析时,栽培种之间的遗传距离为0,栽培种与近缘野生种的遗传距离和近缘野生种之间的遗传距离在0.000 94~0.005 84之间,说明mtDNA序列和matK序列在同属种间的系统发育和种内材料的亲缘关系分析中具有一定的参考价值。使用邻接法构建3个栽培种及8个近缘野生种的系统进化树,11个材料被分为2大类:7个近缘野生种归为Ⅰ类,3个栽培种和空心菜归为Ⅱ类。甘薯栽培种和近缘野生种聚类在不同组别中,表明mtDNA序列和cpDNA序列可用于甘薯栽培种和近缘野生种的关系鉴定。空心菜和3个甘薯栽培种归为一类,空心菜与3个栽培种的遗传距离一致,四者之间的遗传距离小于空心菜与‘月光花’‘黄色朝颜’、变色牵牛和旋转牵牛之间的遗传距离,大于空心菜与树牵牛之间的遗传距离,推测空心菜是旋花科植物进化的重要研究对象。系统发育树中,‘黄色朝颜’与3个甘薯栽培种的亲缘关系较远,遗传距离较大,‘月光花’与‘黄色朝颜’的亲缘关系较远,遗传距离较大,‘黄色朝颜’与‘月光花’可用于甘薯育种的亲本选配。

    3个甘薯栽培种及8个近缘野生种的mtDNA和cpDNA序列分析显示,mtDNA序列和cpDNA序列可用于甘薯栽培种及近缘野生种的遗传多样性和亲缘关系研究,为旋花科植物的分类和系统发育研究提供指导,还可为甘薯种质资源的开发利用和筛选提供理论依据。

  • 图  1   邻接法构建系统发育进化树

    图中数字表示自展值百分比/%

    Figure  1.   The phylogenetic tree constructed by neighbor-joining method

    The number in the figure indicates the bootstrap value percentage/%

    图  2   单倍型中介邻接网络图

    H_1~ H_9:单倍型,mv:可能存在的原始单倍型

    Figure  2.   Median-joining network for haplotypes

    H_1−H_9: Haplotype, mv: The possible original haplotype

    表  1   试验采用引物信息

    Table  1   The information of primers used in this experiment

    类别
    Species
    基因序列
    Gene sequence
    引物
    Primer
    片段大小/bp
    Fragment size
    退火温度/℃
    Annealing temperature
    参考文献
    Reference
    mtDNA ccb203 F: ASGTTCTACGGACCGATGCC 500 57 [22]
    R: CACGGGGAGGGAGCRGGCGA
    ccb256 F: GGAAGTTAGCAAAGTTAGAC 520 57 [22]
    R: TTGTTCTTAACAGCGATGGC
    nad2/1-2 F: TTTTCTTCCTCATTCTKATTT 1200 57 [22]
    R: CCACTCTATTGTCCACTTCTA
    nad2/4-5 F: TTCATATAGAATCCATGTCC 1800 57 [22]
    R: CTATTTGTTCTTCGCCGCTT
    nad5/4-5 F: CCAATTTTTGGGCCAATTCC 1400 57 [23]
    R: CATTGCAAAGGCATAATGAT
    nad7/1-2 F: ACCTCAACATCCTGCTGCTC 1200 57 [23]
    R: CGATCAGAATAAGGTAAAGC
    cpDNA matK F: CGTACAGTACTTTTGTGTTTACGAG 1500 57 [24]
    R: ACCCAGTCCATCTGGAAATCTTGGTTC
    下载: 导出CSV

    表  2   试验材料的序列多态性信息

    Table  2   Sequence polymorphism information of experiment materials

    基因序列
    Gene
    sequence
    长度/
    bp
    Length
    GC占
    比/%
    GC
    proportion
    核苷酸多态性
    Nucleotide
    diversity
    变异位点数
    Variable
    site
    number
    单一突变位点数
    Singleton
    variable
    site number
    简约信息位点数
    Parsimony
    informative
    site number
    插入/缺失位点数
    Insertion/
    deletion
    site number
    ccb203 494 46.52~47.17 0.003 20 5 2 3 6
    nad2/1-2 1 247 53.90~54.69 0.000 17 1 1 0 11
    nad2/4-5 1 503 48.02~48.54 0.003 43 12 3 9 63
    nad5/4-5 1 665 45.93~46.50 0.003 95 20 10 10 6
    nad7/1-2 964 56.39~57.23 0.001 86 4 0 4 24
    matK 840 32.90~34.09 0.007 11 27 23 4 1
    合并 Total 6 713 47.79~48.31 0.003 25 69 39 30 111
    下载: 导出CSV

    表  3   试验材料的单倍型多样性

    Table  3   Haplotype diversity of experiment materials

    基因序列
    Gene
    sequence
    单倍型数量
    Haplotype
    number
    单倍型多样性
    Haplotype
    diversity
    单倍型多样性方差
    Variance of
    haplotype diversity
    单倍型多样性标准差
    Standard deviation of
    haplotype diversity
    ccb203 5 0.818 0.006 82 0.083
    nad2/1-2 2 0.182 0.020 61 0.144
    nad2/4-5 8 0.927 0.004 42 0.066
    nad5/4-5 8 0.927 0.004 42 0.066
    nad7/1-2 5 0.818 0.006 82 0.083
    matK 9 0.945 0.004 34 0.066
    合并 Total 9 0.945 0.004 34 0.066
    下载: 导出CSV

    表  4   Tajima’s D测验和Fu and Li’s D*/F*测验

    Table  4   Tajima’s D test and Fu and Li’s D*/F* test

    基因序列 Gene sequence Tajima’s D Fu and Li’s D* Fu and Li’s F* P
    ccb203 −0.321 97 −0.083 18 −0.161 10 > 0.10
    nad2/1-2 −1.128 50 −1.289 46 −1.399 19 > 0.10
    nad2/4-5 0.500 00 0.574 40 0.628 71 > 0.10
    nad5/4-5 −0.356 59 −0.514 26 −0.537 17 > 0.10
    nad7/1-2 1.018 28 1.214 66 1.312 69 > 0.10
    matK −1.627 84 −2.009 96 −2.168 40 > 0.05, < 0.10
    合并Total −0.639 26 −0.796 77 −0.858 74 > 0.10
    下载: 导出CSV

    表  5   11份试验材料的遗传距离1)

    Table  5   Genetic distance among 11 experimental materials

    样品编号
    Sample No.
    1 2 3 4 5 6 7 8 9 10
    1
    2 0.003 47
    3 0.002 99 0.003 63
    4 0.003 78 0.001 42 0.004 42
    5 0.002 99 0.005 21 0.004 73 0.005 84
    6 0.003 15 0.001 10 0.003 78 0.000 94 0.004 89
    7 0.002 52 0.003 15 0.003 63 0.003 78 0.004 26 0.003 15
    8 0.001 89 0.003 15 0.002 68 0.004 10 0.003 63 0.003 47 0.002 52
    9 0.003 47 0.003 63 0.003 15 0.003 78 0.004 89 0.003 15 0.003 15 0.003 15
    10 0.003 47 0.003 63 0.003 15 0.003 78 0.004 89 0.003 15 0.003 15 0.003 15 0.000 00
    11 0.003 47 0.003 63 0.003 15 0.003 78 0.004 89 0.003 15 0.003 15 0.003 15 0.000 00 0.000 00
     1) 1:何鲁牵牛,2:‘黄茑萝’,3:空心菜,4:‘月光花’,5:‘黄色朝颜’,6:变色牵牛,7:旋转牵牛,8:树牵牛,9:‘浙紫薯3号’,10:‘鄂薯6号’,11:‘皖薯7号’
     1) 1: Ipomoea holubii, 2: Ipomoea hederifolia var. lutea, 3: Ipomoea aquatica Forsk, 4: Ipomoea alba, 5: Ipomoea obscura Keniak, 6: Ipomoea indica, 7: Ipomoea digitata huge caudex, 8: Ipomoea carnea, 9: Ipomoea batatas ‘Zhezishu 3’, 10: Ipomoea batatas ‘Eshu 6’, 11: Ipomoea batatas ‘Wanshu 7’
    下载: 导出CSV
  • [1]

    SHEKHAR S, MISHRA D, BURAGOHAIN A K, et al. Comparative analysis of phytochemicals and nutrient availability in two contrasting cultivars of sweet potato (Ipomoea batatas L.)[J]. Food Chemistry, 2015, 173: 957-965. doi: 10.1016/j.foodchem.2014.09.172

    [2]

    MOHANRAJ R, SIVASANKAR S. Sweet potato (Ipomoea batatas [L. ] Lam): A valuable medicinal food: A rievew[J]. Journal of Medicinal Food, 2014, 17(7): 733-741. doi: 10.1089/jmf.2013.2818

    [3]

    YANG Z, ZHU P, KANG H, et al. High-throughput deep sequencing reveals the important role that microRNAs play in the salt response in sweet potato (Ipomoea batatas L.)[J]. BMC Genomics, 2020, 21(1). doi: 10.1186/s12864-020-6567-3.

    [4] 谢一芝, 郭小丁, 贾赵东, 等. 中国食用甘薯育种现状及展望[J]. 江苏农业学报, 2018, 34(6): 1419-1424. doi: 10.3969/j.issn.1000-4440.2018.06.030
    [5] 胡玲, 李强, 王欣, 等. 甘薯地方品种和育成品种的遗传多样性[J]. 江苏农业学报, 2010, 26(5): 925-935. doi: 10.3969/j.issn.1000-4440.2010.05.006
    [6] 李强, 刘庆昌, 马代夫. 甘薯近缘野生种研究利用现状及展望[J]. 分子植物育种, 2006, 4(6S): 105-110.
    [7] 曹清河, 张安, 李鹏, 等. 甘薯近缘野生种的抗病性鉴定与新型种间杂种的获得[J]. 植物遗传资源学报, 2009, 10(2): 224-229.
    [8]

    GARRIDO-CARDENAS J A, MESA-VALLE C, MANZANO-AGUGLIARO F. Trends in plant research using molecular markers[J]. Planta, 2018, 247(3): 543-557. doi: 10.1007/s00425-017-2829-y

    [9]

    YANG X S, SU W J, WANG L J, et al. Molecular diversity and genetic structure of 380 sweetpotato accessions as revealed by SSR markers[J]. Journal of Integrative Agriculture, 2015, 14(4): 633-641. doi: 10.1016/S2095-3119(14)60794-2

    [10] 苏一钧, 王娇, 戴习彬, 等. 303份甘薯地方种SSR遗传多样性与群体结构分析[J]. 植物遗传资源学报, 2018, 19(2): 243-251.
    [11] 季志仙, 王美兴, 范宏环, 等. 基于ISSR指纹的甘薯食用品种的遗传多样性分析[J]. 核农学报, 2014, 28(7): 1197-1202. doi: 10.11869/j.issn.100-8551.2014.07.1197
    [12] 王崇, 王连军, 苏文瑾, 等. 基于cpSSR标记的甘薯品种亲缘关系及遗传多样性分析[J]. 分子植物育种, 2020, 18(5): 1687-1696.
    [13]

    LEE K J, LEE G A, LEE J R, et al. Genetic diversity of sweet potato (Ipomoea batatas L. Lam) germplasms collected worldwide using chloroplast SSR markers[J]. Agronomy, 2019, 9(11): 725-740. doi: 10.3390/agronomy9110725

    [14]

    GUALBERTO J M, NEWTON K J. Plant mitochondrial genomes: Dynamics and mechanisms of mutation[M]//MERCHANT S S. Annual Review of Plant Biology: Vol 68. Palo alto: Annual Reviews. 2017: 225-252.

    [15]

    PELLETIER G, BUDAR F. The molecular biology of cytoplasmically inherited male sterility and prospects for its engineering[J]. Current Opinion in Biotechnology, 2007, 18(2): 121-125. doi: 10.1016/j.copbio.2006.12.002

    [16]

    LILLY J W, BARTOSZEWSKI G, MALEPSZY S, et al. A major deletion in the cucumber mitochondrial genome sorts with the MSC phenotype[J]. Current Genetics, 2001, 40(2): 144-151. doi: 10.1007/s002940100238

    [17]

    QIN Y J, BUAHOM N, KROSCH M N, et al. Genetic diversity and population structure in Bactrocera correcta (Diptera: Tephritidae) inferred from mtDNA cox1 and microsatellite markers[J]. Scientific Reports, 2016, 6: 38476. doi: 10.1038/srep38476.

    [18]

    BOUILLÉ M, SENNEVILLE S, BOUSQUET J. Discordant mtDNA and cpDNA phylogenies indicate geographic speciation and reticulation as driving factors for the diversification of the genus Picea[J]. Tree Genetics & Genomes, 2011, 7(3): 469-484.

    [19]

    GODBOUT J, JARAMILLO-CORREA J P, BEAULIEU J, et al. A mitochondrial DNA minisatellite reveals the postglacial history of jack pine (Pinus banksiana), a broad-range North American conifer[J]. Molecular Ecology, 2005, 14(11): 3497-3512. doi: 10.1111/j.1365-294X.2005.02674.x

    [20]

    HU D, LUO Z. Polymorphisms of amplified mitochondrial DNA non-coding regions in Diospyros spp.[J]. Scientia Horticulturae, 2006, 109(3): 275-281. doi: 10.1016/j.scienta.2006.02.027

    [21] 马丽, 周玉亮, 郝兆祥, 等. 石榴种质资源的matK基因序列分析[J]. 分子植物育种, 2020, 18(16): 5274-5279.
    [22]

    DUMINIL J, PEMONGE M H, PETIT R J. A set of 35 consensus primer pairs amplifying genes and introns of plant mitochondrial DNA[J]. Molecular Ecology Notes, 2002, 2(4): 428-430. doi: 10.1046/j.1471-8286.2002.00263.x

    [23]

    DUMOLIN-LAPEGUE S, PEMONGE M H, PETIT R J. An enlarged set of consensus primers for the study of organelle DNA in plants[J]. Molecular Ecology, 1997, 6(4): 393-397. doi: 10.1046/j.1365-294X.1997.00193.x

    [24]

    SUN X Q, ZHU Y J, GUO J L, et al. DNA barcoding the Dioscorea in China, a vital group in the evolution of monocotyledon: Use of matK gene for species discrimination[J]. PLoS One, 2012, 7(2): e32057. doi: 10.1371/journal.pone.0032057

    [25]

    KUMAR S, STECHER G, LI M, et al. MEGA X: Molecular evolutionary genetics analysis across computing platforms[J]. Molecular Biology and Evolution, 2018, 35(6): 1547-1549. doi: 10.1093/molbev/msy096

    [26]

    ROZAS J, FERRER-MATA A, SÁNCHEZ-DELBARRIO J C, et al. DnaSP 6: DNA sequence polymorphism analysis of large data sets[J]. Molecular Biology and Evolution, 2017, 34(12): 3299-3302. doi: 10.1093/molbev/msx248

    [27] 刘峥, 张汉尧. 旋花科植物ITS序列分析[J]. 西部林业科学, 2012, 41(4): 70-74. doi: 10.3969/j.issn.1672-8246.2012.04.012
    [28] 俞立璇, 刘美艳, 曹清河, 等. 栽培种甘薯及其近缘野生种nrDNA ITS序列分析[J]. 植物科学学报, 2014, 32(1): 40-49.
    [29]

    JIN D P, LEE J H, XU B, et al. Phylogeography of East Asian Lespedeza buergeri (Fabaceae) based on chloroplast and nuclear ribosomal DNA sequence variations[J]. Journal of Plant Research, 2016, 129(5): 793-805. doi: 10.1007/s10265-016-0831-2

    [30] 郭亚龙, 葛颂. 线粒体nad1基因内含子在稻族系统学研究中的价值: 兼论Porteresia的系统位置[J]. 植物分类学报, 2004, 42(4): 333-344.
    [31] 陆佳妮, 赵志礼, 倪梁红, 等. 线粒体nad1/b-c及nad5/d-e在秦艽组植物中物种鉴定意义的评价[J]. 药物学报, 2019, 54(1): 166-172.
    [32]

    SNOUSSI H, DUVAL M F, GARCIA-LOR A, et al. Assessment of the genetic diversity of the Tunisian citrus rootstock germplasm[J]. BMC Genetics, 2012: 13. doi: 10.1186/1471-2156-13-16.

    [33]

    YANG J, VÁZQUEZ L, CHEN X D, et al. Development of chloroplast and nuclear DNA markers for Chinese Oaks (Quercus Subgenus Quercus) and assessment of their utility as DNA barcodes[J]. Frontiers in Plant Science, 2017: 8. doi: 10.3389/fpls.2017.00816.

    [34]

    KORNELIUSSEN T S, MOLTKE I, ALBRECHTSEN A, et al. Calculation of Tajima’s D and other neutrality test statistics from low depth next-generation sequencing data[J]. BMC Bioinformatics, 2013: 14. doi: 10.1186/1471-2105-14-289.

  • 期刊类型引用(0)

    其他类型引用(1)

图(2)  /  表(5)
计量
  • 文章访问数:  493
  • HTML全文浏览量:  10
  • PDF下载量:  671
  • 被引次数: 1
出版历程
  • 收稿日期:  2020-11-22
  • 网络出版日期:  2023-05-17
  • 刊出日期:  2021-07-09

目录

/

返回文章
返回