基于16SrRNA部分序列探讨12种鲹科鱼类的分子系统进化关系

  • 格式:pdf
  • 大小:584.46 KB
  • 文档页数:8

第32卷第6期2008年11月水 产 学 报J OURNAL OF FISH ERIES OF CHINAVol.32,No.6Nov.,2008收稿日期:2007209222资助项目:国家自然科学基金(30771652);华南农业大学校长基金(K 06141)作者简介:郑文娟(1982-),女,浙江湖州人,硕士,从事鱼类分子遗传的研究,E 2mail :zhengwj 37965@ 通讯作者:邹记兴,E 2mail :zoujixing @文章编号:1000-0615(2008)06-0847-08基于16S r RNA 部分序列探讨12种鲹科鱼类的分子系统进化关系郑文娟1,2, 朱世华1,2, 邹记兴3, 杨迎春1,2, 沈锡权1(1.宁波大学应用海洋生物技术教育部重点实验室,浙江宁波 315211;2.宁波大学生命科学与生物工程学院生物系,浙江宁波 315211;3.华南农业大学动物科学学院,广东广州 510642)摘要:通过PCR 扩增获得了中国海域的鲹科(Carangidae )8属9种的线粒体16S rRNA 序列片段约598bp 碱基,结合来自GenBank 的3种鲹科鱼类的相应片段序列,并以大斑石鲈为外群,生成供系统发育分析的序列矩阵,利用M EGA version 310软件分析序列的碱基组成、差异百分比和转换/颠换值等,应用最大简约法和邻接法构建系统树。

结果显示:(1)鲹科鱼类的16S rRNA 序列片段生成的序列矩阵中发现有碱基的插入和缺失,共有146bp 变异位点,转换/颠换值为2117,表明基因序列的突变未达到饱和,碱基平均差异为8122%;(2)支持鲹科下设四个亚科(鲹亚科,亚科,鲳鲹亚科,鲹亚科)阶元的分类系统;(3)鲹亚科鲹属下不宜设亚属分类阶元;(4)及达副叶鲹与丽叶鲹亲缘关系近,16S rRNA 序列片段碱基只有1107%的差异,未达到分属水平。

关键词:鲹科;线粒体DNA ;16S rRNA 序列;系统发育中图分类号:S 917 文献标识码:A 鲹科(Carangidae )属鲈形目(Perciformes ),是重要的海水经济鱼类,广泛分布在世界的三大洋(印度洋、太平洋、大西洋),尤其是热带或亚热带海域。

目前,多数学者认为全世界鲹科鱼类分成四个亚科(鲹亚科、鲳鲹亚科、亚科和鲹亚科),32属约140种[1-2],中国分布有21属58种[3],也有学者认为有74种[4]。

基于形态的传统分类与基于分子遗传水平构建的鲹科四个亚科的系统关系存在争议[5-10],同时鲹科鱼类由于同种异名和异种同名的现象比较多,如:及达副叶鲹、丽叶鲹等鲹科鱼类的属种分类关系问题就存在较大争议[2,5,11]。

目前,鲹科鱼类亚科之间及属种之间的系统关系仍不明确,国外对鲹科的分子分类及系统关系问题的研究报道也很少[10,12-14]。

随着现代分子生物学技术高速发展,通过获得物种特定遗传标记的大量数据,构建系统进化树,可以对传统分类提供重要佐证以及补充或修正。

线粒体DNA (mitochondrial DNA ,mtDNA )作为分子标记已成为分子系统学研究中的热点[15-19]。

而mtDNA 的16S rRNA 序列是非编码蛋白质基因,由于未受密码子编码的选择压力影响,其大部分区域发生的变异为中性突变,加之mtDNA 进化速度明显大于细胞核DNA ,因此利用线粒体rRNA 序列同源性比较研究物种间的系统进化关系具有一定的优越性[20],在鱼类中,线粒体16S rRNA 序列也被广泛用于研究不同阶元的系统发育关系[21-22]。

本研究测定了分布于我国的9种鲹科鱼类的16S rRNA 序列片段,并结合从GenBank 中下载的3种鲹科鱼类和作为外类群大斑石鲈(Pom adasys m aculates)的16S rRNA序列相应片段进行了比较分析,再运用最大简约法(M P)和邻接法(NJ)来构建分子系统树,以此探讨鲹科鱼类的分子系统发育。

1 材料与方法1.1 材料研究直接测定了中国鲹科鱼类的8属9种、并结合从GenBank下载3种鲹科的16S rRNA 序列相应片段,总共形成10属12种鲹科鱼类一起进行分析,以同为鲈形目石鲈科的大斑石鲈作为外类群(表1)。

1.2 方法基因组DNA提取和PCR扩增 基因组DNA从保存于95%乙醇的肌肉样品中提取,采用苯酚/氯仿抽提法[18]。

用于PCR扩增的引物为16SAR(5′-GCC T GT T TA TCA AAA ACA T-3′)和16SBR(5′-CCG GTC T GA ACT CA G A TC ACG T-3′)[23]。

PCR反应的模板DNA约为100ng,反应体系总体积为50μL,其中10×Ex T aq Buffer 5μL,dN TPs2μL(各215mmol・L-1),引物各1μL(20μmol・L-1),ExT aq酶013μL(2U)。

PCR反应条件为:94℃预变性4min;94℃变性45s,52℃退火45s,72℃延伸1min,34个循环,最后72℃延伸10min。

每次反应设立不含DNA模板的空白对照,扩增产物经1%的琼脂糖凝胶电泳检测。

表1 鲹科鱼类材料和数据来源T ab.1 The origins of tissue samples of C arangid ae f ishes and DNA d ata亚科subfamilies属genera种species登录号GenBank Accession鲹亚科Caranginae丝鲹属A lectis长吻丝鲹A.i ndica,RüppellEF6132653副叶鲹属A lepes及达副叶鲹A.d jedaba,Forssk…lEF6132693丽叶鲹A.klei nii,BlochEF6132643鲹属Caranx黑尻鲹C.melam p y gus,CuvierAP004445若鲹属Carangoi des甲若鲹C.armat us,RüppellAP004444长吻若鲹C.chrysophrys,CuvierEF6132673圆鲹属Decapterus蓝圆鲹D.maruadsi,Temminck and SchlegelEF6132663细鲹属Selaroi des金带细鲹S.leptolepis,CuvierEF6132703竹鱼属T rachurus竹鱼T.j aponicus,Temminck and SchlegelAP003091亚科Seriolinae条属Zonicht hys黑纹条Z.ni g rof asciata,RüppellEF6132713鲳鲹亚科Trachinotinae鲳鲹属T rachinot us狮鼻鲳鲹T.blochii,LacepèdeEF6132633䲠鲹亚科Chorineminae䲠鲹属Scomberoi des康氏似鲹S.commersonianus Lacepède EF6132683石鲈科Pomadasyidae石鲈属Pomadasys大斑石鲈P.maculates,BlochA F247443 注:3为本实验材料Notes:3The sequences were sequenced in t his study848水 产 学 报32卷  PCR产物电泳割胶回收,用凝胶纯化试剂盒(上海申能)纯化,与p UCm2T载体(上海生工)连接,转化感受态细胞D Hα5,筛选阳性克隆子,提取克隆质粒,用M13引物在自动测序仪(Applied Bio systems3730,上海英俊)正反双向测序,以保证所测序列的准确性。

序列分析 测序后的序列,删除引物序列,与GenBank下载的序列一起用Clustal W排序,生成供系统发育分析的矩阵。

利用M EGA version310软件分析序列的碱基组成和差异百分比、变异位点、简约信息位点数、转换/颠换值。

系统分析用最大简约法(MP)和邻接法(NJ)。

系统分析与分子进化分析使用软件包PAU P Version410b10[24]。

最大简约法(M P)分析使用启发式搜寻(heuristic search),参数设置为:1000 rando m stepwise addition sequence replicates,构树方法采用t ree bisection reco nnection(TBR),所有数据未加权。

邻接法分析使用K imura22 parameter。

系统树分支的置信度采用重复抽样分析(Boot st rap analysis)的方法,重复抽样的次数为1000次[25]。

2 结果2.1 鲹科鱼类16S r RNA序列特征16S rRNA序列经PCR扩增、回收、克隆,本实验测定的9种鲹科鱼类得到长598bp的碱基序列,与GenBank下载的3种鲹科鱼类16S rRNA序列一起读入Clustal W中排序,生成供系统发育分析的序列矩阵中序列长度共有为608 bp碱基,其间存在插入和缺失现象,主要集中在245~280bp和370~380bp两个区域,此外共有146bp变异位点,其中简约信息位点为78个bp。

所得到的16S rRNA序列中T、C、A、G碱基平均含量分别为2314%、2611%、2816%、2119%,其中A+T的含量(5210%)高于G+C含量(4810%),但总体上看来碱基组成的差别并不大,这在其它鱼类也出现类似情况[23,26]。

除去外类群,所分析的鲹科鱼类碱基差异为1107%~14100%,平均差异为8122%。

其中及达副叶鲹与丽叶鲹之间碱基差异最小,只有1107%,从表2中我们可以看出鲹亚科的康氏似鲹与其他鲹科鱼类的碱基差异都很大。

表2 鲹科和外群的16S rRNA部分序列的碱基差异百分比(下三角)和转换/颠换(上三角)T ab.2 Percentage divergence(below diagonal)and number of transitions/transversions(above diagonal)for partial16S rRNA sequences in C arangid ae and the outgroup12345678910111213 120/819/823/817/419/116/626/1114/1029/2226/1936/2942/3325.136/018/427/1227/921/1022/520/1029/2232/1942/2740/2734.941.0716/426/1227/923/1023/520/1030/2227/1939/2741/2745.714.013.6323/1019/919/1221/520/1236/2228/1938/2740/2953.827.267.066.099/520/827/1320/1227/2227/1736/2548/3363.646.686.685.122.5123/722/1024/923/2329/2034/2843/3474.005.706.095.695.135.5229/138/632/2030/1941/2941/3186.874.955.154.767.455.897.8533/1128/2328/1841/2639/3094.365.505.505.885.886.092.518.2736/1833/1741/2544/29 109.619.619.8111.069.208.599.839.6010.2639/1743/2946/31 118.419.638.618.818.229.219.228.629.4410.6938/2248/26 1212.4813.3712.7212.5011.6711.8513.5612.9512.7414.0011.4937/38 1314.6212.9313.1513.3515.9715.0513.9813.3414.2315.0914.4914.57 Codes1-13: A.i ndica, A.d jedaba, A.kleinii, C.melam p y gus, C.armat us, C.chrysop hrys, D.maruadsi,S.leptolepis,T. j aponicus,Z.ni grof asciata,T.blochii,mersonianus,P.maculates 序列中的转换明显比颠换多,转换/颠换(Ts/Tv)比值平均为2117(K imura22parameter),通常认为Ts/Tv大于210时,基因序列的突变未达到饱和[27]。