全 文 :北方园艺2016(16):103~107 ·生物技术·
第一作者简介:赵博(1981-),女,河南南阳人,博士,副研究员,研究
方向为植物细胞学与分子系统发育。E-mail:caljone@163.com.
责任作者:唐文秀(1973-),女,本科,副研究员,现主要从事珍稀濒
危植物保护等研究工作。E-mail:tangwx99@163.com.
基金项目:广西自然科学基金资助项目(2012GXNSFBA053075)。
收稿日期:2016-04-26
DOI:10.11937/bfyy.201616027
基于叶绿体ndhA基因内含子序列的DNA
条形码在秋海棠属物种鉴定中的应用
赵 博1,李 景 剑2,毛 世 忠1,唐 文 秀1
(1.中国科学院 广西植物研究所,广西 桂林541006;2.华南农业大学 林学与风景园林学院,广东 广州510642)
摘 要:以18种秋海棠属植物为试材,基于叶绿体ndhA基因内含子片段,利用Codon Code
Aligner软件拼接序列后,通过MEGA软件进行比对并计算19个秋海棠属植物种间及种内遗传
距离,最后利用邻接法构建分子系统树,研究了该片段作为DNA条形码对秋海棠属植物进行物
种鉴定的可行性。结果表明:秋海棠属19个种类共59个个体的ndhA基因内含子序列长度为
1 182bp,在所考察的候选秋海棠属植物中具有最大的种间变异和较小的种内变异,且二者存在
极显著差异。在系统树中,秋海棠属植物每一物种能够分别形成各自独立的分支,表现出良好单
系性。基于ndhA基因内含子的DNA条形码在识别秋海棠属植物物种方面和传统形态学基本
一致。该研究表明,以ndhA基因内含子作为秋海棠属植物DNA条形码进行物种鉴定具有一定
的可行性。
关键词:秋海棠属;ndhA基因内含子;DNA条形码;物种鉴定
中图分类号:Q 948 文献标识码:A 文章编号:1001-0009(2016)16-0103-05
秋海棠属(Begonia L.)为秋海棠科(Begoniaceae)一
年或多年生草本、灌木或小树枝状植物,全世界约1 500~
1 600种[1]。在中国约有173种,分布于长江以南的地
区,常见于云南、广西和贵州[2-3]。由于秋海棠属植物在
习性、生境、性别、形态等方面均表现出很高的多样性,
品种繁多,其中大部分具有极高的观赏价值,在国内外
花卉市场颇受人们的喜爱[4]。另外,秋海棠属植物也是
我国少数民族中常见的民族药,具有清热解毒、散瘀消
肿的功效[5-6]。由此可见,秋海棠属植物中蕴藏着丰富
的观赏、药用资源有待进一步的开发利用。然而,
秋海棠属植物性状受环境影响较大,即使同一种类的不
同居群,在某些性状上的变异也比较显著,很多近缘种
类难以鉴定[7]。同时该属大多是多汁草本植物,具有肉
质的茎、娇嫩的花和果实,传统方法压制成标本后,在叶
片颜色、花纹、毛被颜色和花被片颜色等方面都有很大
程度的失真,果实和子房的形态解剖学特征等是属下分
组不可缺少的重要鉴定特征,在多汁的果实被压干以后
容易变形,给传统的分类鉴定带来较大的困难[8-9]。因
此,为秋海棠属植物寻找快速、高效的分类鉴定方法一
直都是研究的热点。
DNA条形码(DNA barcoding)是通过使用短的标准
DNA片段,对物种进行快速、准确的识别和鉴定[10-13]。
DNA条形码不仅是传统植物分类与鉴定的强有力补
充,而且能突破对经验的过度依赖,可以在较短时间内
建成易于利用的应用系统,在物种分类和鉴定方面展示
出了强大的生命力[14]。植物DNA条形码片段主要在叶
Abstract:Taking chestnut as test material,GBSS gene was cloned from the first strand of Castanea molssina seed’s
cDNA through RT-PCR methods,and the bioinformatics characters were analysed.The results showed that a ful length
cDNA of Cgbss(GenBank accession number:KU162945)was obtained.The length of Cgbss gene was 2 085bp,which
contained an open reading frame encoding 611amino acids.The bioinformatics analysis indicated that the pI and MW of
animo acids encoded by Cgbss were predicted to be 8.36and 67.580 9kDa,respectively,and the protein had one GT1_
glycogen_synthase domain.And the fragment had 85%homology identical to the already reported sequence of GBSS I.
Keywords:Castanea molssina;granule-bound starch synthase;clone;sequence analysis
301
·生物技术· 北方园艺2016(16):103~107
绿体基因组上进行选择,此外还有核基因ITS[15-17]。因
为用于DNA条形码研究的标准基因,一方面应该足够
保守,能够利用通用引物进行大范围的扩增;另一方面
应该有足够的变异来区分不同物种DNA序列,从而进
行鉴定[18]。而目前为止,植物DNA条形码尚处于对所
提议的各片段比较和评价阶段,还未获得统一的标
准[19]。
目前已有关于利用DNA条形码对秋海棠属植物进
行物种鉴定的相关研究报道,但是这些研究结果表明,
叶绿体基因rbcL、matK和psbA-trnH种内和种间变异
都较小,都不足以区分秋海棠属植物,只有ITS/ITS2种
内和种间变异大,可考虑作为秋海棠属DNA条形码鉴
定的候选片段[8,20-21]。因此,在鉴定秋海棠属植物时还
需选择更多的有效基因。该试验应用植物DNA条形码
片段ndhA基因内含子区对中国秋海棠属植物19个种
类共59个个体进行了研究,评价其PCR扩增和测序成
功率,种内种间差异及物种鉴定能力,探讨该基因片段
是否适合作为秋海棠属的候选DNA条形码。
1 材料与方法
1.1 试验材料
采集分布于中国广西、云南和广东地区的秋海棠属
植物18种共58个个体(表1),标本存放于中国科学院广
西植物研究所,同时从GenBank数据库中下载Begonia
oxyloba序列1条(GenBank ID:JF756335),共计19种59
个个体。
1.2 试验方法
1.2.1 秋海棠DNA的提取 由于秋海棠属植物的
DNA较难提取,在提取总DNA时,用液氮将叶片彻底
磨成粉末状后,加入CTAB-free+2%PVP溶液将粉末混
匀,然后放入-20℃冰冻10min后,离心弃上清液,再
采用植物基因DNA提取试剂盒(Tiangen Biotech Co.,
China)按说明书操作步骤提取样品总DNA。通过紫外
分光光度计和1% 琼脂糖凝胶电泳检测DNA的质量和
浓度,最后将浓度调整为30~50ng·μL-1,于-20℃贮
存备用。
表1 供试18种秋海棠样品情况
Table 1 Begonia oxylobasamples for test
样品名称
Samples name
采集地点
Location
GenBank号
GenBank No.
样品名称
Samples name
采集地点
Location
GenBank号
GenBank No.
Begonia jingxiensis 广西桂林植物园、广西靖西
KT599066
KT599067
KT599068
Begonia pseudodaxinensis 广西桂林植物园、广西大新
KT599047
KT599048
KT599049
Begonia filiformis 广西桂林植物园、广西龙州
KT599077
KT599078
KT599079
Begonia grandis 广西桂林植物园
KT599072
KT599073
Begonia umbraculifolia 广西桂林植物园、广西大新
KT599045
KT599046
Begonia chingi 广西桂林植物园
KT599091
KT599092
KT599093
Begonia leprosa 广西桂林植物园、广西大尧山
KT599063
KT599064
KT599065
Begonia fimbristipula 广西桂林、广东鼎湖山、云南昆明植物研究所
KT599039
KT599040
KT599041
KT599042
KT599043
Begonia daxinensis 广西桂林植物园、广西大新
KT599088
KT599089
KT599090
Begonia palmata 广西金秀
KT599054
KT599055
KT599056
Begonia picturata 广西桂林植物园、广西靖西
KT599050
KT599051
KT599052
KT599053
Begonia longifolia 云南
KT599060
KT599061
KT599062
Begonia luochengensis 广西桂林植物园、广西罗成
KT599057
KT599058
KT599059
Begonia handeli 广西桂林植物园、广西金秀
KT599069
KT599070
KT599071
Begonia fangi 广西桂林植物园、广西龙州
KT599080
KT599081
KT599082
Begonia cathayana 广西桂林植物园、广西上思
KT599094
KT599095
KT599096
Begonia gigaphylla 广西桂林植物园、广西龙州
KT599074
KT599075
KT599076
Begonia edulis 广西桂林植物园、广西靖西
KT599083
KT599084
KT599085
KT599086
KT599087
401
北方园艺2016(16):103~107 ·生物技术·
1.2.2 序列扩增 选用ndhA基因内含子区通用引物
ndhAx1(GCY CAA TCW ATT AGT TATG AAA
TACC)和ndhAx2(GGT TGA CGC CAM ARA TTC
CA)。PCR扩增反应体系(25μL):2×Taq PCR Mix
12.5μL,引物(2.5μmol·L
-1)各1.0μL、DNA模板2μL,
最后用ddH2O补足至25μL。PCR扩增程序:94℃预变
性5min,94℃变性30s,57℃退火45s,68℃延伸2min,
循环40次,72℃延伸7min,4℃保存。PCR产物用1%
琼脂糖凝胶电泳检测后由上海生工生物公司完成测序
工作。
1.2.3 序列分析 测序所获得的ndhA基因内含子区
域峰图经序列拼接软件Codon Code Aligner(Codon Code
Co.,USA)校对拼接,并依据NCBI数据库上已测ndhA
基因内含子序列对所测序列去除低质量序列并结合人
工校对,确定每一位点准确无误。利用MEGA 5.1软件
进行多序列比对、查错并计算变异位点。利用Kimura
两参数模型(Kimura-2-parameter,K2P)计算样品间两两
遗传距离(pairwise distances)、种内距离(distances within
species)和种间距离(distances between species)。采用邻
接距离法(Neighbor-Joining)构建系统发育树,选用K2P
替代模型,进行1 000次自检。遗传距离的计算和系统
发育树的构建均采用MEGA 5.0软件。
2 结果与分析
分子生物学研究需要制备用于PCR反应的核酸,
前期的研究工作发现,秋海棠属植物的核酸非常难以提
取,推测可能是该属植物富含多糖、多酚、酯类和酸类等
次生代谢产物,导致传统的CTAB法和试剂盒法无法提
取到可用于扩增的DNA。课题组应用改良试剂盒法提
取的秋海棠属植物DNA可以满足于SSR引物开发、系
统发育和DNA条形码等分子标记研究。DNA条形码
技术的核心为PCR扩增,而延伸温度是影响PCR扩增
的重要参数,在DNA聚合酶允许的温度范围内,选择合
适的延伸温度可相对提高DNA聚合酶的活性,从而促
进引物和模板间的结合,提高PCR扩增的效率。该试验
对ndhA基因内含子序列的扩增过程中发现,当延伸温
度为72℃时,ndhA基因内含子的引物对于秋海棠属植
物的扩增能力较差,有些物种未能成功扩增出条带;延
伸温度改为68℃时扩增时,58个个体的电泳条带明亮
单一,且测序成功率为100% (图1),所以确定大DNA
条形码引物PCR反应的最适延伸温度为68℃。
58个个体中,ndhA基因内含子的PCR扩增成功率
和测序成功率(测序后获得高质量的序列且正反向序列
顺利拼接即判定为成功)为100%。所得的58条序列比
对后,ndhA基因内含子序列长度为1 113~1 130bp,平
均长度为1 182bp,种间存在101个变异位点。种间和
种内遗传距离的大小是进行物种鉴别的主要标准。结
图1 ndhA基因PCR产物电泳结果
Fig.1 DNA electrophoresis results of ndhA PCR products
果显示,在所考察的候选秋海棠属植物中具有最大的种
间变异和较小的种内变异,种间遗传距离为0.019 40±
0.000 19,种内遗传距离为0.000 20±0.000 06,二者存
在极显著差异。
图2 基于ndhA基因内含子序列的秋海棠属聚类分析
Fig.2 Cluster analysis of Begoniabased on ndhA sequence
通过对全部共59条ndhA基因内含子序列构建邻
接(NJ)树(图2),不同种之间能够明显分开,18个种类
501
·生物技术· 北方园艺2016(16):103~107
的所有个体都形成单系,除了B.gigaphyla的分支支持
率没有高于50%,其它秋海棠的支持率均高于50%,物
种分辨率较高。
3 讨论
该试验采用改良的植物试剂盒方法提取秋海棠属
植物DNA,相较于常规的试剂盒法,仅在试验前期增加
了一个步骤,即在叶片彻底磨成粉末状后,加入CTAB-
free+2%PVP溶液预处理,再按常规的试剂盒提取流程
提取。该方法提取的DNA纯度较高,浓度也较大,说明
该方法适合于提取秋海棠属植物的DNA。
生物条形码联盟(CBOL)认为理想的DNA条形码
应该符合以下3个标准:1)在种内的变异应足够小,同时
在种间有明显的遗传变异和分化可将物种区分开来;2)
所用的片段尽量短,便于PCR扩增,而且一个反应就能
完成测序,尤其是对存在DNA降解的材料(如:炮制过
的药材、存放已久的蜡叶标本);3)两端连接相对保守的
区域,便于设计通用引物。与动物学领域相比,在植物
中的DNA条形码的研究进展相对缓慢,主要有2个方
面原因:1)植物线粒体基因组进化速率较慢,遗传分化
小,因此动物中的标准片段COI不适用于植物;2)系统
学研究中常用的片段(如trnL-F)变异较小,不适合用作
条形码片段[19]。由于有些植物的染色体是多倍体,核基
因组通常具有多拷贝的特性,引物通用性差。因此,植
物中最可能的条形码还是从叶绿体基因组中选择[22-23]。
虽然叶绿体基因相对保守,但仍然包含许多变异区域,
同时叶绿体基因组有其自身的优势:单亲遗传避免了基
因重组;植物个体中均有大量的叶绿体,即使DNA高度
降解也容易扩增。生物条形码联盟建议的植物叶绿体
条形码片段有matK、psbA-trnH、rbcL。
在有关于秋海棠属植物的DNA条形码研究中,焦
丽娟等[8]的研究结果表明,3个叶绿体片段(rbcL、matK
和psbA-trnH)可能是由于进化速率慢,谱系分化不完全
或是秋海棠属自身物种分化时间短等生物学原因导致
种内、种间距离较小,不存在遗传差异鸿沟(barcoding
gap),物种正确鉴定率较低,不适合作为秋海棠属的标
准DNA鉴定条形码。在秋海棠属中,只有核基因ITS
片段种内变异和种间变异比较大,存在遗传差异鸿沟,
单片段的物种正确鉴定率达到100%,可考虑作为秋海
棠属DNA条形码鉴定的标准条码。但是秋海棠属植物
含有不规则的多倍体,ITS片段有多拷贝现象。因此,在
鉴定秋海棠属植物时还需开发有效的叶绿体基因片段。
近年来,ndhA基因内含子片段在一些类群中已被
证实具有很好的鉴定效果。该研究中,基于ndhA基因
内含子区的遗传距离分析表明,种内距离小于种间距
离,存在遗传差异鸿沟,物种分辨率较高,只有B.
gigaphyla物种的分支支持率没有高于50%(图2)。虽
然ndhA基因内含子的各种支持率与ITS相比稍有降
低,但ndhA基因内含子和其它3个片段(rbcL、matK和
psbA-trnH)联合分析的物种正确鉴定率可能会提高,那
么ndhA基因内含子可以作为秋海棠属植物DNA条形
码的一个有效的补充片段。
综上所述,该研究通过对秋海棠属19个物种59个
个体的ndhA基因内含子序列分析,发现这些类群中
ndhA基因内含子序列能够较好的区分各物种,18个种
类的所有个体都形成单系类群,对秋海棠属物种而言,
是一个有效的DNA条形码序列。该研究为探讨以
ndhA基因内含子序列作为DNA条形码对秋海棠属植
物进行物种识别的可行性提供了一定的参考。
参考文献
[1] HUGHES M,HOLLINGSWORTH P M.Population genetic divergence
corresponds with species level biodiversity patterns in the large genus Begonia
[J].Molecular Ecology,2008,17(11):2643-2651.
[2] KU S M,LIU Y,PENG C.Four new species of Begoniasect.Coelocen-
trum(Begoniaceae)from limestone areas in Guangxi,China[J].Bot Stud(Tai-
pei,Taiwan),2006,47:207-222.
[3] de WILDE J J F E,van VALKENBURG J L C H.Begonia sosefiana
(Begoniaceae):A new species in section Loasibegonia from Gabon[J].
Blumea,2005,50:467-471.
[4] THOMAS D C,HUGHES M,PHUTTHAI T,et al.A non-coding
plastid DNA phylogeny of Asian Begonia(Begoniaceae):Evidence for mor-
phological homoplasy and sectional polyphyly[J].Molecular Phylogenetics and
Evolution.2011,60(3):428-444.
[5] 张少平,赖正锋,吴水金,等.药食同源植物紫背天葵研究现状与展
望[J].中国农学通报,2014,30(4):58-61.
[6] 王红珊,曹毅敏,李国豪,等.紫背天葵提取物对糖尿病肾病大鼠的
作用[J].中国生化药物杂志,2012,33(3):272-274.
[7] FORREST L L,HUGHES M,HOLLINGSWORTH P M.Phylogeny
of Begoniausing nuclear ribosomal sequence data and morphological characters
[J].Systematic Botany,2005,30(3):671-682.
[8] 焦丽娟,税玉民.中国秋海棠属(秋海棠科)植物的DNA条形码评价
[J].植物分类与资源学报,2013,35(6):715-724.
[9] 张嵘梅,陈文红,税玉民,等.中国秋海棠属植物的叶表皮特征及其
分类学意义[J].云南植物研究,2008,30(6):665-678.
[10]HEBERT P D N,CYWINSKA A,BALL S L,et al.Biological identifications
through DNA barcodes[J].Proceedings of the Royal Society of London Series
B-Biological Sciences,2003,270(1512):313-321.
[11]CHEN S L,YAO H,HAN J P,et al.Validation of the ITS2region as a
novel DNA barcode for identifying medicinal plant species[J].PloS One,
2010,5(1):e8613.
[12]彭居俐,王绪祯,王丁,等.基于线粒体CO1基因序列的DNA条形码
在鲤科鲌属鱼类物种鉴定中的应用[J].水生生物学报,2009,33(2):271-276.
601
北方园艺2016(16):103~107 ·生物技术·
[13]张辉,姚辉,崔丽娜,等.基于CO1条形码序列的《中国药典》动物药
材鉴定研究[J].世界科学技术-中医药现代化,2013,15(3):371-380.
[14]LI D Z,LIU J Q,CHEN Z D,et al.Plant DNA barcoding in China[J].
Journal of Systematica and Evolution,2011,49(3):165-168.
[15]CHASE M W,SALAMIN N,WILKINSON M,et al.Land plants and
DNA barcodes:Short-term and long-term goals[J].Philos Trans R Soc Lond
B Biol Sci,2005,360:1889-1895.
[16]KRESS W J,WURDACK K J,ZIMMER E A,et al.Use of DNA bar-
codes to identify flowering plants[J].Proc Natl Acad Sci,USA,2005,102:
8369-8374.
[17]KRESS W J,ERICKSON D L.A two-locus global DNA barcode for
land plants:The coding rbcL gene complements the non-coding trnH-psbA
spacer region[J].PLoS One,2007,2(6):e508.
[18]FAZEKAS A J,BURGESS K S,KESANGAKURTI P R,et al.Multiple
multilocus DNA barcodes from the plastid genome discriminate plant species
equaly wel[J].PLoS One,2008,3:e2802.
[19]宁淑萍,颜海飞,郝刚,等.植物DNA条形码研究进展[J].生物多样
性,2008,16(5):417-425.
[20]THOMAS D C,HUGHES M,PHUTTHAI T,et al.A non-coding
plastid DNA phylogeny of Asian Begonia(Begoniaceae):Evidence for mor-
phological homoplasy and sectional polyphyly[J].Molecular Phylogenetics and
Evolution,2011,60:428-444.
[21]CHUNG K F,LEONG W C,RUBITE R R,et al.Phylogenetic analyses
of Begoniasect.Coelocentrumand alied limestone species of China shed light
on the evolution of Sino-Vietnamese karst flora[J].Botanical Studies,2014,
55:1.
[22]COWAN R S,CHASE M W,KRESS W J,et al.300000species to identify:
Problems,progress,and prospects in DNA barcoding of land plants[J].Taxon,
2006,55,611-616.
[23]NEWMASTER S G,FAZEKAS A J,RAGUPATHY S.DNA barcoding
in land plants:Evaluation of rbcL in a multigene tiered approach[J].Canadian
Journal of Botany,2006,84:335-341.
Application of DNA Barcoding Based on the Chloroplast ndhA
Intron Region in Classification of Begonia(Begoniaceae)
ZHAO Bo1,LI Jingjian2,MAO Shizhong1,TANG Wenxiu1
(1.Guangxi Institute of Botany,Chinese Academy of Sciences,Guilin,Guangxi 541006;2.Colege of Forestry and Landscape Architecture,
South China Agricultural University,Guangzhou,Guangdong 510642)
Abstract:Taking 58samples from 18species of Begonia as materials,chloroplast ndhA intron region was used as a
popular marker of DNA barcoding for species identification to evaluate its species resolution rates in the genus Begonia.
The sequences of the ndhA intron region were assembled using the Codon Code Aligner,and then aligned and the genetic
distances were computed using MEGA 5.0in accordance with the kimura 2-parameter(K2P)model.The neighbor-
joining(NJ)phylogenetic trees were constructed.The results showed that the ndhA intron region was about 1 182bp in
length,and exhibited the highest inter-specific divergence,and this was significantly higher than the intra-specific
variation.The ndhA intron region had reasonable high identification eficiency among al tested samples.The NJ trees
revealed that al individuals of each species formed a strong monophyletic group.Species identification with ndhA intron
region in cluster supported the taxonomy based on morphological characters.The conclusion suggested that chloroplast
ndhA intron region was a valid DNA barcoding gene for species identification in genus Begonia.
Keywords:Begonia;ndhA intron region;DNA barcoding;species identification
701