免费文献传递   相关文献

基于RLCKVⅡ基因序列的蜘蛛抱蛋属系统发育分析



全 文 :  文章编号:0427-7104(2013)04-0436-06
收稿日期:2013-01-22
基金项目:国家自然科学基金(30970199)资助项目
作者简介:黄丹妮(1987—),女,硕士研究生;王玉国,男,博士,通讯联系人,E-mail:wangyg@fudan.edu.cn.
基于RLCKVII基因序列的蜘蛛抱蛋属
系统发育分析
黄丹妮,宋志平,王玉国,陈家宽
(复旦大学 生物多样性和生态工程教育部重点实验室,上海200433)
摘 要:蜘蛛抱蛋属(Aspidistra Ker-Gawl.)是单子叶植物中同一属内物种多样性最为丰富的类群之一.由于
缺乏分子序列方面的研究资料,该属种间系统发育关系一直都不清楚.本研究对20种具有代表性的蜘蛛抱蛋属
植物的RLCKVII基因进行了测序分析.结果表明,采用分子系统学证据所构建的该属系统发育关系不支持先前
基于形态证据所建立的分类系统;蜘蛛抱蛋属植物复杂的平行进化关系很可能是该属分化早期频繁种间杂交的
结果.RLCKVII基因适用于禾本科、棕榈科和天门冬科等类群,可作为研究单子叶植物系统发育关系的备选
基因.
关键词:蜘蛛抱蛋属;系统发育;RLCKVII;辐射进化
中图分类号:Q 941        文献标志码:A
蜘蛛抱蛋属(Aspidistra Ker-Gawl.)是广义百合科(Liliaceae)铃兰族(Convalarieae)中种数最多的
属,最新分子系统发育证据将其置于天门冬科(Asparagaceae)酒瓶兰亚科(Nolinoideae)[1].该属现有104
种[2],主要集中分布于中国与越南,我国共有蜘蛛抱蛋属植物69种,其中56种为我国特有.蜘蛛抱蛋属植
物不但种类丰富、特有种多,而且植株形态、果实特征、尤其花的各部分结构复杂多样[3],是探讨植物复杂
性状发生、适应性辐射进化的理想材料.然而,蜘蛛抱蛋属植物的系统发育关系尚不十分清楚,严重制约了
相关研究的深入.
李光照等曾根据植物形态演化趋势,对蜘蛛抱蛋属提出了一个3组20系的分类系统[4].目前在
GenBank中仅可以查到少数蜘蛛抱蛋属植物的叶绿体基因序列,都是在研究其他类群的系统发育时测定
的,而专门针对蜘蛛抱蛋属植物系统学方面的研究还非常少.唐绍清等[5]和黄海等[6]曾对几种蜘蛛抱蛋属
植物的ITS序列进行了分析,但这些研究涉及的种数和克隆序列数都很少,且未设定外类群,难以全面反
映蜘蛛抱蛋属植物属内种间的系统发育关系,也没有对李光照等提出的分类系统进行评价.
在被子植物的系统发育研究中,ITS(核糖体转录间隔区)和cpDNA(叶绿体DNA)非编码区序列在低
分类阶元类群分析中得到了最为广泛的应用[7].与它们相比,理想的核基因序列可在分析时使系统发育信
息最大化,同时使同源基因的数目最小化[8,9].在进行种内或种间水平分析时,往往选择含有长的内含子
或基因间隔区的核基因序列以提高分辨率.编码丝氨酸/苏氨酸蛋白激酶家族 RLCKVII(Serine/
Threonine Protein Kinase RLCKVII)基因具备上述理想核基因的潜力,该基因的序列与水稻第2号染色
体的LOC_Os02g30900序列同源,通过与葱、芭蕉表达序列标签(Expressed Sequence Tag)和水稻BAC
(Bacterial Artificial Chromosome)序列进行比对,可以设计得到保守的内含子搜索引物(Conserved
Intron-Scanning Primers,CISPs),而且该基因片段在单子叶植物水稻(Oryza)、芭蕉(Musa)、葱
(Allium)、高粱(Sorghum)和棕榈(Palm)中能够获得单一扩增,在水稻中扩增得到的片段长度为
1 422bp[10],并包含较长的内含子片段.RLCKVII基因现已被成功应用于棕榈科(Palmae)系统分类和进
化研究[11,12],表明其在单子叶植物研究中可能的应用前景.
第52卷 第4期
2013年8月
复 旦 学 报 (自然科学版)
Journal of Fudan University(Natural Science)
Vol.52No.4
Aug.2013
本研究参照李光照等的分类系统,共选取了20个具有代表性的蜘蛛抱蛋属植物,对它们的RLCKVII
基因序列进行测定分析,并从GenBank数据库下载禾本科和棕榈科的同源序列作为外类群,以此构建蜘
蛛抱蛋属植物的系统发育树,进一步探讨该属植物的进化历史.
1 材料与方法
1.1 材料
用于本研究的20个植物样品均来自广西桂林植物研究所植物园(采自模式产地),涵盖了李光照分类
系统的3组12系的蜘蛛抱蛋属植物(表1).凭证标本藏于广西植物研究所标本馆(IBK)及复旦大学标本
室(FUS).供提取总DNA的叶片材料均采自活体,取生长旺盛的鲜叶5~10g,立即放于100g硅胶中进
行干燥,常温条件下保存.
1.2 方法
1.2.1 总DNA提取
采用CTAB法[13]或DNeasy Plant Mini DNA extraction Kits(天根生化科技有限公司)提取总DNA,
具体操作步骤参照Song等[14]的描述.
1.2.2 PCR扩增、电泳检测及纯化
PCR反应在PTC-200扩增仪(Bio-Rad公司)上进行,反应程序为:94℃预变性4min,38个循环,每
个循环94℃1min,52℃2min,72℃3min,最后72℃延伸10min.采用的50μL反应体系,包含:5μL
10×Buffer,6μL 25mmol/L MgCl2,1.6μL 10mmol/L dNTP,2.5μL 10μmol/L引物,50ng DNA模
板,0.5单位 Taq酶和30μL三蒸水.RLCKVII特异基因引物序列为:正向引物 RLCKVII-F(5-
GGAGATCAACGTCTTCTTGT-3)和反向引物RLCKVII-R(5-GCATACTCAGGAGCACAATA-3),
由上海生工生物工程有限公司合成.dNTP购于博彩生物科技有限公司,其余试剂购于 TaKaRa公司.
PCR产物取5μL于1%琼脂糖凝胶电泳检测,使用DL2000Marker(TaKaRa)进行片段长度标记.采用北
京博大泰克生物有限公司出品的B型小量DNA片段快速割胶回收试剂盒,回收目标条带.
1.2.3 PCR产物的克隆、测序及数据分析
克隆具体操作步骤参照《分子克隆实验指南》[15]的实验流程,克隆载体来自TaKaRa公司的pMD18-
T vector试剂盒,感受态细胞E.coli Top10购自天根生化科技有限公司.克隆菌液样品送往上海美季生
物技术有限公司,利用PE ABIPRISM TM-377型 DNA 自动测序仪对克隆载体上的 M13F(~47)和
M13R(~48)通用引物进行测序,得到的双向序列采用SeqManⅡ软件进行拼接,获得的单一序列和查自
GenBank的外类群同源序列用Clustal X程序进行对位排列[16]后再手工校对.所得数据用PAUP*4.0b
软件构建[17]邻接树和最大简约树,保留支持率50%以上的节点.
2 结果与分析
经克隆测序得到101条蜘蛛抱蛋属植物RLCKVII基因序列,每个样品至少2个单一克隆.合并查自
GenBank的34条外类群序列(表1,见第438页),共计对135条同源序列进行了对位排列.由于序列涉及
单子叶植物的3个科,不同来源的同源基因内含子序列变异较大,因此对位排列引入了较多的空位,最终
序列长度为2 023bp,其中变异位点1 742个,信息位点818个,占序列总长的40.4%;蜘蛛抱蛋属植物序
列中变异位点568个,信息位点375个,占总长的18.5%.每个蜘蛛抱蛋属植物样品有1~3个序列类型,
拥有两个以上类型的样品其不同序列类型间相似度均大于99%.
根据被子植物系统发育小组(Angiosperm Phylogeny Group,APG)所构建的单子叶植物系统树,将
棕榈科同源序列设为系统树中最基部的类群,得到的邻接树和最大简约树拓扑结构一致,本文以最大简约
树作为代表(图1(a),见第439页).从图1(a)可以看出,蜘蛛抱蛋属、禾本科和棕榈科分别为单系类群,支
持率均为100%.在蜘蛛抱蛋属分支内部,有18个种的多个克隆各自聚类,说明RLCKVII基因序列在蜘
蛛抱蛋属植物中分辨率较高.虽然该基因在不同种蜘蛛抱蛋中具有区分度高的特征位点,但该基因整体核
酸替代速率较低,种间进化关系仍较模糊,呈现辐射进化式样,不能支持基于形态特征建立的蜘蛛抱蛋属
734 第4期 黄丹妮等:基于RLCKVII基因序列的蜘蛛抱蛋属系统发育分析
植物分类系统.此外,研究涉及的20个蜘蛛抱蛋属植物中有6个种的不同克隆序列出现在系统树的多个
位置,例如长瓣蜘蛛抱蛋(A.longipetala)共有7个克隆,分为3种类型,聚于系统树的3个不同位置.而
南宁蜘蛛抱蛋(A.sp.)的克隆序列网结于长梗蜘蛛抱蛋(A.longipedunculata)的分支之中.
表1 本研究所用蜘蛛抱蛋属和外类群样品、凭证标本和序列来源
Tab.1 Locality,voucher and sequence source information for samples of Aspidistraand outgroups
种  名 凭 证 标 本 序 列 来 源
小花蜘蛛抱蛋A.minutiflora  APRG 035(IBK) 本研究
啮边蜘蛛抱蛋A.marginella  APRG 001(IBK) 本研究
盈江蜘蛛抱蛋A.yingjiangensis  APRG 071(IBK) 本研究
乐山蜘蛛抱蛋A.leshanensis  APRG 036(IBK) 本研究
隆安蜘蛛抱蛋A.longanensis  APRG 008(IBK) 本研究
罗甸蜘蛛抱蛋A.luodianensis  APRG 011(IBK) 本研究
西林蜘蛛抱蛋A.xilinensis  Wang Y004(FUS) 本研究
黄花蜘蛛抱蛋A.flaviflora  Wang Y013(FUS) 本研究
线叶蜘蛛抱蛋A.linearifolia  APRG 056(IBK) 本研究
广东蜘蛛抱蛋A.lurida  APRG 018(IBK) 本研究
辐花蜘蛛抱蛋A.subrotata  APRG 053(IBK) 本研究
十字蜘蛛抱蛋A.cruciformis  Wang Y111(FUS) 本研究
糙果蜘蛛抱蛋A.muricata  APRG 033(IBK) 本研究
海南蜘蛛抱蛋A.hainanensis  APRG 042(IBK) 本研究
南宁蜘蛛抱蛋A.sp. Wang Y120(FUS) 本研究
巨型蜘蛛抱蛋A.longiloba  APRG 041(IBK) 本研究
峨眉蜘蛛抱蛋A.omeiensis  Wang Y127(FUS) 本研究
长瓣蜘蛛抱蛋A.longipetala  Chen 93190(IBK) 本研究
长梗蜘蛛抱蛋A.longipedunculata  Wang Y023(FUS) 本研究
粽叶草A.oblanceifolia  Wang Y005(FUS) 本研究
大麦 Hordeum vulgare  GenBank  AK371224
玉米Zea mays GenBank  BT041945
水稻Oryza sativa  GenBank  AP003266,NC008394
高粱Sorghum bicolor  GenBank  NC012878
华盛顿棕榈Washingtonia filifera  GenBank  HQ720959,DQ897848
Pritchardia lanaiensis  GenBank  JF904975,JF904977,JF904976,JF904978
Pritchardia flynnii  GenBank  JF904954,JF904956,JF904957
Pritchardia waialealeana  GenBank  JF905022,JF905020
Pritchardia mitiaroana  GenBank  JF904999
Pritchardia martii  GenBank  JF904993
Pritchardia glabrata  GenBank  JF904961
Pholidocarpus macrocarpus  GenBank  HQ720926,HQ720927
Brahea dulcis  GenBank  HQ720851
Saribus rotundifolius  GenBank  HQ720949
Saribus jeanneneyi GenBank  HQ720947
山棕榈Trachycarpus martianus  GenBank  HQ720956
Maxburretia rupicola  GenBank  HQ720924
Livistoninae sp. GenBank  JF292932,JF292933,JF292935,JF292937,JF292940
Johannesteijsmannia perakensis  GenBank  HQ720888
Licuala bacularia  GenBank  HQ720890
糖棕Borassus flabellifer  GenBank  HQ720847
834 复 旦 学 报(自然科学版)    第52卷
图1 基于RLCKVII基因序列构建的系统发育树:(a)最大简约树;(b)棕榈科分支
Fig.1 Phylogenetic tree based on RLCKVII gene:(a)maximum parsimony trees;(b)the branch of Palmae
注:括号中数字表示克隆数;初生组、三线组、四线组为李光照分类系统之组名;南宁蜘蛛抱蛋未定种.
3 讨 论
李光照等曾基于柱头表面纹饰等形态特征对蜘蛛抱蛋属进行过较全面的分类研究,并提出了以初生
组为原始形态,进而逐渐演变出三线组和四线组各种蜘蛛抱蛋属植物的演化趋势[4],但本研究采用
RLCKVII基因序列构建的系统发育树并不支持以上结论.在RLCKVII基因系统树上,来自初生组、三线
组和四线组的蜘蛛抱蛋属植物都没有构成单系类群,除了6个种有不同类型的拷贝之外,大部分物种呈彼
此独立的平行关系,没有明显的发生先后,特别是被认为柱头表面纹饰最为简单的长梗蜘蛛抱蛋,并未出
现在蜘蛛抱蛋属分支更基部的位置.而西林蜘蛛抱蛋(A.xilinensis)和黄花蜘蛛抱蛋(A.flaviflora)所
934 第4期 黄丹妮等:基于RLCKVII基因序列的蜘蛛抱蛋属系统发育分析
构成的分支却以100%的bootstrap支持率处于该属的基部.与蜘蛛抱蛋属不同的是,在RLCKVII基因树
中,棕榈科金棕属各物种之间的关系得以明显分辨,说明蜘蛛抱蛋属的物种间平行进化可能是快速辐射进
化的结果,而非该基因的分辨率低所造成.考虑到序列中的内含子片段较长(对位排列后为1 686bp),而
内含子的进化速率相对较快,进一步说明RLCKVII基因在蜘蛛抱蛋中的低核酸变异速率可能是分化时
间较短所造成.有研究证实,与蜘蛛抱蛋属近缘的舞鹤草属(铃兰族)分化估计在中新世中期之后,距今约
8.3个百万年[18],因此蜘蛛抱蛋属的辐射进化的物种形成模式也可能是一个近期事件.
RLCKVII基因树,在禾本科和棕榈科中能够很好地区分各属,且进化关系明确(图1(b)).事实上,该
基因在棕榈科植物进化的研究中[11,12]被证实为单拷贝或低拷贝核基因.据此推测,在蜘蛛抱蛋属6个物
种中检测到的多个克隆类型并非来自同一家族的不同基因,而是同一基因的不同拷贝.洪德元等[19]最早
进行蜘蛛抱蛋属植物细胞分类学研究,后由黄锦岭等[20]和王任翔等[21]补充,证实蜘蛛抱蛋属植物的染色
体倍性有2x和4x两种类型.本研究没有发现RLCKVII基因的拷贝数与染色体倍性的直接关系,例如四
倍体线叶蜘蛛抱蛋(A.linearifolia)的3个克隆聚在一处,四倍体十字蜘蛛抱蛋(A.cruciformis)的5个
克隆分于两处,另有二倍体罗甸蜘蛛抱蛋(A.luodianensis)的5个克隆也分散在两处(图1(a)).
Yamashita和Tamura[22]通过分析铃兰族系统发育树不同分支的核型和染色体数目,得出x=19为原始
数目,而x=18是从x=19演化而来,说明该属植物很可能是古多倍体.在被子植物中,约70%的种类在
进化史中曾发生过一次或多次多倍化的过程[23],我们推测在蜘蛛抱蛋进化历史中很可能也曾发生过多倍
化事件,并伴随着二倍化过程,该属分化早期频繁的种间杂交使该属植物核基因呈现出复杂的辐射进化
关系.
本研究中,RLCKVII基因无法区分南宁蜘蛛抱蛋与长梗蜘蛛抱蛋,在系统树中,南宁蜘蛛抱蛋的6
个克隆网结于长梗蜘蛛抱蛋之中.由于南宁蜘蛛抱蛋属于未定种,其与长梗蜘蛛抱蛋很可能是同一个
物种.除此之外,每个蜘蛛抱蛋属植物都有各自独特的变异,故RLCKVII基因能够区分本研究涉及的
所有明确的物种.据此,RLCKVII基因在单子叶植物的3个不同的科(天门冬科,禾本科和棕榈科)的
物种鉴定上得到了很好的应用,表明其作为备选核基因,在单子叶植物系统发育研究中可能会有广泛
的应用价值.
参考文献:
[1] Chase M W,Reveal J L,Fay M F.A subfamilial classification for the expanded asparagalean families
Amarylidaceae,Asparagaceae and Xanthorrhoeaceae[J].Botanical Journal of the Linnean Society,
2009,161(2):132-136.
[2] Tilich H J.An updated and improved determaination key for Aspidistra Ker-Gawl.(Ruscaceae,
Monocotyledons)[J].Feddes Repertorium,2008,119(5-6):449-462.
[3] 李光照.蜘蛛抱蛋属植物[M].1版.南宁:广西科学技术出版社,2004:33-41.
[4] 李光照,郎楷永,王任翔,等.中国蜘蛛抱蛋属植物形态演化趋势及新分类系统[J].广西植物,2000,
20(3):201-207.
[5] 唐绍清,苏何玲,李光照.8种蜘蛛抱蛋属植物ITS区序列及其系统学意义[J].广西师范大学学报,2003,
21(4):86-90.
[6] 黄 海,王安仙,刘 杨,等.贵州地区蜘蛛抱蛋属植物ITS1区序列分析及其亲缘关系研究[J].贵州科
学,2010,28(2):50-52.
[7] Bortiri E,Oh S,Jiang J,Baggett S,et al.Phylogeny and systematic of the genus Prunus L.(Rosaceae)
as determined by sequence analysis of ITS and the chloroplast trnL-trnF spacer DNA[J].Systematic
Botany,2001,26(4):797-807.
[8] Oh S,Potter D.Phylogenetic utility of the second intron of LEAFYin Neillia and Stephanandra
(Rosaceae)and implications for the origin of Stephanandra[J].Molecular Phylogenetics and Evolution,
2003,29(2):203-215.
[9] Sang T.Utility of low-copy nuclear gene sequences in plant phylogenetics[J].Critical Reviews in
Biochemistry and Molecular Biology,2002,37(3):121-147.
044 复 旦 学 报(自然科学版)    第52卷
[10] Lohithaswa H C,Feltus F A,Singh H P,et al.Leveraging the rice genome sequence for monocot
comparative and translational genomics[J].Theoretical and Applied Genetics,2007,115(2):237-243.
[11] Bacon C D,McKenna M J,Simmons M P,et al.Evaluating multiple criteria for species delimitation:an
empirical example using Hawaian palms(Arecaceae:Pritchardia)[J].BMC Evolutionary Biology,2012,
12:23-40.
[12] Bacon CD,Bader W J,Simmons M P.Miocene dispersal drives island radiations in the palm tribe
Trachycarpeae(Arecaceae)[J].Systematic Biology,2012,61(3):426-442.
[13] Doyle J J,Dolye J L.A rapid DNA isolation procedure for smal quantities of fresh leaf tissue[J].
Phytochemical Bulletin,1987,19:11-15.
[14] Song Z P,Xu X,Wang B,et al.Genetic diversity in the northernmost Oryza rufipogon populations
estimated by SSR markers[J].Theoretical and Applied Genetics,2003,107(8):1492-1499.
[15] 萨姆布鲁克J,弗里奇 EF,曼尼阿蒂斯 T.分子克隆实验指南[M].2版.黄培堂译.北京:科学出版
社,1996.
[16] Thompson J D,Gibson T J,Plewniak F,et al.The Clustal X windows interface:Flexible strategies for
multiple sequence alignment aided by quality analysis rools[J].Nucleic Acids Research,1997,25(24):
4876-4882.
[17] Swofford D L.PAUP4.0phylogenetic analysis using parsimony(and other methods)[M].Sunderland,
Massachusetts:Sinauer Associates,1999.
[18] Meng Y,Wen J,Nie Z L,et al.Phylogeny and biogeographic diversification of Maianthemum
(Ruscaceae:Polygonatae)[J].Molecular Phylogenetics and Evolution,2008,49(2):424-434.
[19] 洪德元,郎楷永,张志宪.蜘蛛抱蛋属(百合科)的细胞分类学研究[J].植物分类学报,1986,24(5):
353-361.
[20] 黄锦岭,马黎明,洪德元.蜘蛛抱蛋属的细胞分类学研究Ⅱ[J].植物分类学报,1997,35(1):14-23.
[21] 王任翔,李光照,郎楷永,等.中国蜘蛛抱蛋属的细胞分类学研究Ⅱ[J].广西植物,2000,20(2):138-143.
[22] Yamashita J,Tamura M N.Phylogenetic analyses and chromosome evolution in Convalarieae(Ruscaceae
sensu lato),with some taxonomic treatments[J].Journal of Plant Research,2004,117(5):363-370.
[23] Masterson J.Stomatal size in fossil plants:Evidence for polyploidy in majority of angiosperms[J].
Science,1994,264(5157):421-424.
Phylogenetic Analysis of AspidistraInferred from
Sequences of RLCKVII Gene
HUANG Dan-ni,SONG Zhi-ping,WANG Yu-guo,CHEN Jia-kuan
(Ministry of Education Key Laboratory for Biodiversity Science and Ecological Engineering,
Fudan University,Shanghai 200433,China)
Abstract:Aspidistra Ker-Gawl.(Asparagaceae,subfamily Nolinoideae),with high species diversity,is one of the
biggest genera in monocotyledon.However,phylogenetic relationships among species in Aspidistra remain
unknown,due to little molecular sequence analysis.In this study,molecular phylogeny of Aspidistra was
conducted based on the sequence variations of nuclear DNA RLCKVII for 20representative species of Aspidistra.
The results showed that the molecular phylogenetic relationships among Aspidistra species were not consistent
wel with the infrageneric classification based on the morphological characters.The complexity of paralel
phylogenetic relationships indicates that Aspidistra might have experienced frequent episodes of interspecific
hybridization at its initial stages of diversification.The RLCKVII gene can be used for species identification in the
families Poaceae,Palmae and Asparagaceae,suggesting that RLCKVII is a useful candidate gene for phylogenetic
analysis of monocotyledon.
Keywords:Aspidistra;phylogeny;RLCKVII;radiation evolution
144 第4期 黄丹妮等:基于RLCKVII基因序列的蜘蛛抱蛋属系统发育分析