大口黑鲈脂代谢相关基因NPY、UCP2、LPL、HL克隆与分子进化分析
MOLECULAR CLONING AND EVOLUTIONAL ANALYSIS OF NPY,UCP2,LPL AND HL GENE OF LARGEMOUTH BASS(MICROPTERUS SALMOIDES)
-
摘要: 采用RT-PCR和RACE方法克隆了大口黑鲈NPY基因cDNA全序列及UCP2、LPL、HL基因cDNA核心片段。序列分析结果表明,大口黑鲈NPY基因cDNA全序列长664 bp,其中5′端非翻译区(5′-UTR)长53 bp,3′端非翻译区(3′-UTR)长311 bp,开放阅读框(ORF)长300 bp,编码99个氨基酸,即前体NPY。大口黑鲈前体NPY包括三个部分,28个氨基酸组成的信号肽、36氨基酸组成的成熟NPY以及32个氨基酸组成的由Gly-Lys-Arg指示的NPY C端肽(CPON)。大口黑鲈UCP2、LPL、HL基因cDNA核心片段长度分别为737 bp、509 bp和666 bp,各自编码245个氨基酸、169个氨基酸和222个氨基酸。将4个基因的氨基酸序列分别与其他物种的氨基酸序列进行同源性比较,并通过MEGA3构建系统树,对这4个脂代谢相关基因的分子进化特征进行了探讨。Abstract: Largemouth bass(Micropterus salmoides) is a piscivorous teleost fish,found only in the fresh waters of North America.It is a good game fish in the United States,and a highly valued cultivated species in China.In order to study its lipid metabolism as well as its evolutional relationship with other percoids,we used RT-PCR method to get the partial cDNA sequences of NPY,UCP2,LPL and HL genes,and deduced their corresponding amino acid sequences.Rapid amplification of the cDNA ends(RACE) PCR was used to amplify the 3′-and 5′-end of NPY cDNA.The complete NPY cDNA sequence of largemouth bass was 664bp in length,composed of an open reading frame(ORF) of 300bp,encoding a protein of 99 amino acid residues.A search of the BLAST neuropeptide database revealed that the deduced neuropeptide sequence of largemouth bass NPY was very similar to the NPY of other species,such as European sea bass(Dicentrarchus labrax,93﹪ identity),orange-spotted grouper(Epinephelus coioides,91﹪ identity),house mouse(Mus musculus,86﹪ identity),and human(Homo sapiens,87﹪ identity).The obtained UCP2 core sequence was 737bp in length,which contains an open reading frame coding 245 amino acids.The length of the obtained LPL core sequence was 509bp,coding 169 amino acids.The length of the obtained HL core sequence was 666bp,coding 222 amino acids.We made the amino acid sequence comparisons among the four obtained genes of largemouth bass and other species,and got four phylogenetic trees by NJ method using software MEGA3.The phylogenetic tree built on the basis of conservative NPY gene,showed a right relationship of fish with mollusk,amphibia,aves and mammals,but it could not present accurate relationship of fish species.The phylogenetic tree built by less conservative UCP2 gene showed that largemouth bass and other percoids were farther away from mammals than the more primitive cyprinids.In order to make a further determination of the evolutional position of largemouth bass,we used the least conservative genes(LPL and HL) of the four to build phylogenetic trees,and found that largemouth bass was the closest to Chinese perch(Siniperca chuatsi),which was distributed only in the fresh waters of East Asia.These results support the close relationship of largemouth bass and Chinese perch.
-
Keywords:
- Largemouth bass /
- NPY /
- UCP2 /
- LPL /
- HL /
- Molecular cloning /
- Phylogenetic analysis
-
太平洋鳕(Gadus macrocephalus)又称大头鳕, 隶属于鳕形目(Gadiformes)、鳕科(Gadidae)、鳕属(Gadus), 属冷水性底层鱼类, 是重要的海洋经济鱼类[1—3]。近年来, 由于过度捕捞及环境污染等问题, 黄渤海域的太平洋鳕产量显著下降, 美国、加拿大等国也对太平洋鳕的捕捞量开始了严格的限制[4, 5]。
目前, 对太平洋鳕的研究主要集中在人工繁殖和免疫学[6—10]。此外, 尹伟力等[11]设计了鳕物种特异性PCR鉴定方法以区分包括太平洋鳕在内的鳕制品与其他鱼类制品。1987年, Grant等[12]的研究显示, 日本海和黄渤海不同海域的太平洋鳕在某些基因位点存在明显差异, 推测太平洋鳕虽然属于洄游性鱼类, 但这种洄游是一种短距离洄游, 造成了较小范围内的地理隔离。由于样品数量较少, Grant等[12]无法明确划分地理隔离的界限。另一方面, 太平洋鳕溯源以及与其他鳕科鱼类的进化关系仍是值得探讨的问题之一。脊椎动物线粒体DNA因具有自主复制、结构紧密、编码效率高、无重组现象、无组织特异性、进化速率快以及遵循严格的母系遗传等特点, 被广泛地应用于系统发育、物种分类和群体遗传进化等领域中[13]。本研究通过第二代测序技术获得太平洋鳕线粒体全基因组序列, 并对该线粒体全基因组序列进行分析, 为进一步研究鳕系统进化关系、鳕鱼种属间分类提供依据, 并且为系统研究因地理隔离而引发太平洋鳕线粒体变异等问题提供重要数据支持。
1. 材料与方法
1.1 实验材料
本实验所用的太平洋鳕捕自大连旅顺海域。解剖并采集太平洋鳕肌肉、肝脏组织后用去离子水冲洗, 液氮速冻, 转移至–80℃冰箱冷冻保存。
1.2 线粒体DNA提取
组织DNA提取按照北京艾德莱生物科技有限公司DNA快速提取试剂盒说明进行。采用NV 3000微量分光光度计检测DNA浓度, 1%琼脂糖凝胶电泳检测DNA的质量和完整性。
1.3 太平洋鳕线粒体基因组测序
首先根据近缘物种序列设计简并引物(表 1), 使用LA Taq polymerase进行PCR扩增。PCR程序为: 94℃ 30s, 50℃退火30s, 72℃延伸1min/kb。直接用PCR产物在ABI 3730全自动测序仪中测序, 使用DNAstar软件对测序结果进行拼接组装, 并人工分析调整, 最终补全序列。将太平洋鳕线粒体全基因序列上传至NCBI数据库, 获得GenBank序列号为KY296294。
表 1 太平洋鳕线粒体基因组全序列扩增使用的引物序列Table 1. Sequences of primers used in amplification of the complete mitochondrial genome in G. macrocephalus编号
Number引物序列
Primer sequences (5ʹ—3ʹ)基因区域
Gene region1 TGGTCCTGACTTTACTATC/
GTACGCTTACCATGTTACG12S rRNA 2 CCTGAAGGAGGATTTAGCAG/
CATTTCACAGGCAACCAGC12S rRNA/
16S rRNA3 CTACTCGTACCTTTTGCATC/
GAACTCAGATCACGTAGGAC16S rRNA 4 GGTATAACCGATCAACGAAC/
CGATTCGTCTCTGCTAGTG16S rRNA/ND1 5 CATATGAAGTAAGTCTTGG/
GAGAGGTAAGGTAGATGCTGND1/ND2 6 CCTAAGTGAATAATTCTTC/
TCTAGGGAAGAATTGGCTAGND2/COX1 7 CAGTTGGGATAGATGTAGAC/
CATGGAAGTGTAGTAACTCCOX1/COX2 8 CAATCACATGACCACTCTG/
TGCTCATGTAAACATGAAGtRNA-Ser/ATP8 9 GTCGAAGCAGTTCCTCTAG/
AAGCTCAGAAGAAGCCTAGCOX2/COX3 10 GCTTTAGGAACAGTCCTTC/
CCTCCTTGAAGTCATTCGTCOX3/ND3 11 GACTACGAAAAGCTATCTC/
ACTATTGAGCCTGCAACAGND3/ND4 12 CAGCACAGGAACTTTATCG/
GCTCAGGGTGATATCTAGAGND4/ND5 13 GTAACCTTTACATCTCTGCC/
TCTTTAGGTGGAGGACTACND5/ND6 14 CAAGCAACATACAACAAGG/
AGAGCCATAATAGAGACCTCND5/Cytb 15 CTCAACTTCTAACAGGACT/
CTCTATACTGAGCTACTAGGCytb 16 CGTTGCAGATATACTAGTTC/
GTCAAGCTTTCACTTATGGCCytb/12S rRNA 1.4 基因组组装和注释
使用GENEIOUS R8软件进行线粒体基因组注释。利用软件MITOS Web Server (http://mitos.bioinf.uni-leipzig.de/index.py)[14]进行转运RNA(tRNA)及核糖体RNA(rRNA)注释, 基于已有的鳕科鱼类线粒体基因组进行比对, 对蛋白质编码基因(Protein coding genes, PCGs)以及D-Loop区进行注释。利用MEGA 7.0软件进行线粒体DNA序列的碱基组成特征、密码子使用情况以及碱基偏移度分析。
1.5 基因二级结构预测
使用在线tRNAscan-SE软件对太平洋鳕线粒体tRNA基因进行定位, 并使用The mfold Web Server (http://unafold.rna.albany.edu)在线软件对线粒体tRNA的二级结构进行预测。
1.6 系统进化树构建
为了分析鳕科鱼类的系统发育关系, 从GenBank中下载相关的线粒体全序列或部分序列, 以斑马鱼(Danio rerio, NC002333)和红鳍东方鲀(Takifugu rubripes, NC004299)等线粒体基因组全序列为外群, 使用MEGA 7.0软件构建进化树, 分别构建11种鳕科鱼类线粒体基因组全序列和Cytb基因序列系统进化树。本研究所使用的基因组具体信息见表 2。
表 2 线粒体系统发育分析所用基因组Table 2. The complete mitochondrial genomes used in phylo-genetic analysis物种
Species序列长度
Length
(bp)GenBank登录号
GenBank
accession No.大西洋鳕 Gadus morhua 16696 NC 002081 格陵兰鳕 Gadus ogac 15534 NC 012323 黑线鳕 Melanogrammus aeglefinus 16582 AM 489717 牙鳕 Merlangius merlangus 16565 NC 007395 杜氏银大眼鳕 Gadiculus argenteus thori 16713 AP 018148 阿拉斯加狭鳕 Gadus chalcogrammus 16570 AB 182308 挪威狭鳕 Theragra finnmarchica1 16571 AM 489718 Theragra finnmarchica2 16570 AM 489719 蓝鳕 Micromesistius poutassou 16573 NC 015102 北鳕 Boreogadus saida 16745 AM 919428 北极鳕 Arctogadus glacialis 16640 NC 010122 江鳕 Lota lota1 16527 AP 004412 Lota lota2 16575 KC 844053 Lota lota3 16547 KM 201364 红鳍东方鲀 Takifugu rubripes 16449 AP 006045 斜带石斑鱼 Epinephelus coioides 16418 EU 043376 斑马鱼 Danio rerio 16596 NC 002333 大菱鲆 Scophthalmus maximus 17583 NC 013183 2. 结果
2.1 太平洋鳕线粒体基因组结构与特征
本研究采用第二代基因测序技术, 获得太平洋鳕线粒体全基因组序列, 其长度为16569 bp。与大多数脊椎动物线粒体组成相似, 太平洋鳕线粒体全基因组序列包含了13个PCGs、22个tRNA基因、2个rRNA基因和1个D-Loop (图 1)。PCGs、rRNA基因、tRNA基因和D-Loop区分别占整个线粒体长度的69.26%、15.79%、9.30%和5.26%。太平洋鳕线粒体包括D-Loop区在内的各基因大小、位点、排列顺序与大多数脊椎动物相同, 但D-Loop区序列长度有较大差异, 同源性相对较低。除tRNA-Gln、tRNA-Trp、tRNA-Ala、tRNA-Asn、tRNA-Cys、tRNA-Ser、tRNA-Glu和tRNA-Pro以及NAD6基因在H-链编码外其余基因及tRNA均在L-链编码(表 3)。
![]()
图 1 太平洋鳕线粒体基因组结构Figure 1. Schematic Map of the mitochondrial genome of G. macrocephalus表 3 太平洋鳕线粒体基因注释结果Table 3. The result of the complete mitochondrial genome annotation for G. macrocephalus基因
Gene定位Position 大小
Size (bp)密码子Codons 反密码子
Anti-codon编码链
Strand间隔/重叠
Space/Overlap起始位点Start 终止为点End 起始位点Start 终止位点Stop tRNA-Phe 1 68 68 GAA L 12S rRNA 69 1020 952 L tRNA-Val 1021 1092 72 TAC L 16S rRNA 1093 2757 1665 L tRNA-Leu 2758 2831 74 TAA L Nad1 2832 3806 975 ATG TAG L tRNA-Ile 3808 3877 70 GAT L 1 tRNA-Gln 3877 3947 71 TTG H –1 tRNA-Met 3947 4016 70 CAT L –1 Nad2 4017 5063 1047 ATG TAG L tRNA-Trp 5062 5130 69 TCA H –2 tRNA-Ala 5132 5200 69 TGC H 1 tRNA-Asn 5203 5275 73 GTT H 2 tRNA-Cys 5308 5374 67 GCA H 32 tRNA-Tyr 5375 5441 67 GTA L Cox1 5443 6993 1551 GTG TAG L 1 tRNA-Ser 6996 7066 71 TGA L 2 tRNA-Asp 7070 7137 68 GTC L 3 Cox2 7151 7849 699 ATG AGA L 13 tRNA-Lys 7842 7915 74 TTT L –8 Atp8 7917 8084 168 ATG TAA L 1 Atp6 8075 8758 684 ATG TAA L –10 Cox3 8758 9543 786 ATG TAA L –1 tRNA-Gly 9543 9612 70 TCC L –1 Nad3 9613 9963 351 ATG TAG L tRNA-Arg 9962 10030 69 TCG L –2 Nad4L 10031 10327 297 ATG TAA L Nad4 10321 11706 1386 ATG AGA L –7 tRNA1-His 11702 11770 69 GTG L –5 tRNA-Ser 11771 11838 68 TGA H tRNA-Leu 11841 11913 73 TAG L 2 Nad5 11914 13752 1839 ATG TAA L Nad6 13748 14269 522 ATG AGG H –5 tRNA-Glu 14270 14338 69 TTC H Cytb 14343 15483 1161 ATG T-- L 4 tRNA-Thr 15484 15554 70 TGT L tRNA-Pro 15629 15698 70 TGG H 74 2.2 核苷酸组成
太平洋鳕线粒体全基因组的A+T含量达到57.3%, 明显高于C+G的含量, 呈现AT偏向性, 并且本文所选物种均有不同程度的AT偏向(表 4)。太平洋鳕线粒体全基因组碱基含量由高到低依次为: T(29.4%)>A(27.9%)>C(25.9%)>G(16.9%), A+T-skew%碱基组成和G+C-skew% 碱基组成分别为−0.026和−0.211。另外, 本文中所选鳕科鱼类除蓝鳕和生活在淡水中的江鳕外, 其余种类线粒体基因组碱基组成均显示出较弱的AT负偏斜现象(A+T-skew%<0), 而外群中除大菱鲆外均显示AT正偏斜现象。文中各鳕科鱼类线粒体基因组中胸腺嘧啶(T)含量最高的要属黑线鳕(30.6%), 而最低的为蓝鳕(27.1%)(表 4)。
表 4 太平洋鳕及其他物种线粒体基因组碱基组成比较Table 4. Comparison of base composition in mitochondrial genome between G. macrocephalus and other species物种Species 长度Length 碱基组成Composition (%) A+T-skew (%) G+C-skew (%) A T G C A+T C+G 太平洋鳕G. macrocephalus 16569 27.9 29.4 16.9 25.9 57.3 42.7 –0.026 –0.211 大西洋鳕G. morhua 16696 28.1 29.6 16.8 25.6 57.6 42.4 –0.027 –0.209 格林兰鳕G. ogac 15534 27.8 29.2 17.0 26.1 57.0 43.0 –0.025 –0.211 黑线鳕M. aeglefinus 16582 28.5 30.6 16.3 24.6 59.1 40.9 –0.034 –0.204 牙鳕M. merlangus 16565 28.8 30.1 16.2 24.8 59.0 41.0 –0.022 –0.210 杜氏银大眼鳕G. argenteus thori 16713 27.2 27.3 16.6 28.9 54.5 45.5 –0.004 –0.269 阿拉斯加狭鳕G. chalcogrammus 16570 28.1 29.5 16.7 25.7 57.6 42.4 –0.024 –0.212 挪威狭鳕T. finnmarchica1 16571 28.1 29.5 16.7 25.7 57.6 42.4 –0.024 –0.212 挪威狭鳕T. finnmarchica2 16570 28.1 29.5 16.7 25.7 57.6 42.4 –0.023 –0.212 蓝鳕M. poutassou 16573 27.4 27.1 17.2 28.2 54.5 45.5 0.006 –0.242 北鳕B. saida 16745 28.1 29.6 16.7 25.5 57.7 42.3 –0.026 –0.208 北极鳕A. glacialis 16640 28.1 29.9 16.6 25.3 58.1 41.9 –0.031 –0.208 江鳕L. lota1 16527 28.4 27.6 16.3 27.7 56.0 44.0 0.014 –0.259 江鳕L. lota2 16575 28.3 27.6 16.4 27.8 55.9 44.1 0.013 –0.258 江鳕L. lota3 16547 28.3 27.5 16.4 27.8 55.8 44.2 0.014 –0.258 红鳍东方鲀T. rubripes 16449 29.8 25.9 15.2 29.0 55.7 44.3 0.071 –0.311 斜带石斑鱼E. coioides 16418 28.7 26.5 15.8 28.9 55.2 44.8 0.039 –0.293 斑马鱼D. rerio 16596 31.9 28.1 16.0 23.9 60.1 39.9 0.063 –0.198 大菱鲆S. maximus 17583 27.0 28.7 17.1 27.2 55.8 44.2 –0.030 –0.228 2.3 蛋白质编码基因
太平洋鳕线粒体基因组13个蛋白质编码基因中的12个基因(NAD1、NAD2、COX1、COX2、ATP8、ATP6、COX3、NAD3、NAD4L、NAD4、NAD5和Cytb)在L链上编码蛋白, 仅NAD6在H链上编码。编码基因拥有2种典型的起始密码子(GTG和ATG), 除COX1基因使用GTG作为起始密码子, 其余所有蛋白质编码基因均采用ATG作为起始密码子。终止密码子的使用情况则有较大的变化。太平洋鳕线粒体有5种终止密码子(TAG、AGA、TAA、AGG、T--)(表 3), 其中TAG为NAD1、NAD2、COX1和NAD3基因的终止密码子; AGA为COX2、NAD4基因的终止密码子; TAA为ATP8、ATP6、COX3、NAD4L和NAD5基因的终止密码子; AGG为NAD6基因的终止密码子; Cytb基因则使用不完全的T作为终止密码子。其中AGG与AGA在哺乳动物线粒体中作为终止密码子使用更为常见。
13个蛋白编码基因序列总长度为11444 bp, 除了终止密码子外共有3815个密码子。使用MEGA 7.0软件分析密码子平均使用频率和相对同义密码子平均使用频率(表 5)。太平洋鳕线粒体基因组13个蛋白质编码基因中存在32个偏好密码子(RSCU密码子), 其中第三位点为A或T(U)的密码子普遍具有较高的使用频率。除了CAU(H)、CGU(R)、AGU(S)、GGU(G)和CCA(P)、AGA(*)密码子之外, 密码子第三位碱基为A或T的密码子RSCU均大于1, 密码子第三位点对A、T的偏好性与蛋白质编码基因密码子的第三位点对A、T偏向性相一致。
表 5 太平洋鳕13个编码蛋白编码基因密码子使用频率Table 5. Total codon average usage in the thirteen protein-coding genes for G. macrocephalusCodon Count RSCU Codon Count RSCU Codon Count RSCU Codon Count RSCU UUU(F) 133 1.14 UCU(S) 68 1.58 UAU(Y) 64 1.17 UGU(C) 17 1.26 UUC(F) 100 0.86 UCC(S) 48 1.12 UAC(Y) 45 0.83 UGC(C) 10 0.74 UUA(L) 158 1.45 UCA(S) 73 1.70 UAA(*) 5 1.67 UGA(W) 103 1.72 UUG(L) 25 0.23 UCG(S) 7 0.16 UAG(*) 4 1.33 UGG(W) 17 0.28 CUU(L) 182 1.67 CCU(P) 86 1.59 CAU(H) 49 0.92 CGU(R) 15 0.80 CUC(L) 100 0.92 CCC(P) 78 1.44 CAC(H) 57 1.08 CGC(R) 14 0.75 CUA(L) 150 1.37 CCA(P) 44 0.81 CAA(Q) 73 1.55 CGA(R) 29 1.55 CUG(L) 40 0.37 CCG(P) 9 0.17 CAG(Q) 21 0.45 CGG(R) 17 0.91 AUU(I) 178 1.41 ACU(T) 84 1.14 AAU(N) 65 1.17 AGU(S) 24 0.56 AUC(I) 75 0.59 ACC(T) 85 1.16 AAC(N) 46 0.83 AGC(S) 38 0.88 AUA(M) 134 1.51 ACA(T) 116 1.58 AAA(K) 56 1.58 AGA(*) 2 0.67 AUG(M) 44 0.49 ACG(T) 9 0.12 AAG(K) 15 0.42 AGG(*) 1 0.33 GUU(V) 73 1.28 GCU(A) 105 1.19 GAU(D) 46 1.16 GGU(G) 51 0.83 GUC(V) 42 0.73 GCC(A) 117 1.33 GAC(D) 33 0.84 GGC(G) 81 1.32 GUA(V) 87 1.52 GCA(A) 118 1.34 GAA(E) 68 1.40 GGA(G) 71 1.16 GUG(V) 27 0.47 GCG(A) 12 0.14 GAG(E) 29 0.60 GGG(G) 42 0.69 注: 黑体为偏好密码子Note: Bold shows preferred codons 2.4 tRNA基因
对太平洋鳕线粒体22个tRNA基因的定位以及二级结构进行预测分析, 结果表明22个tRNA基因分布于13个蛋白质编码基因之间, 大小从67到74 bp不等, 总长度为1540 bp(表 3)。通过比对发现, 太平洋鳕相邻tRNA基因之间存在相互折叠的现象, 如tRNA-Ile和tRNA-Gln、tRNA-Gln和tRNA-Met、tRNA-Thr和tRNA-Pro之间, 相互折叠的核苷酸数目不定, 1—2个较为多见, 最长可达74 bp。22个tRNA中13个由L-链编码, 其他9个则由H-链编码(表 3)。所有的tRNA基因序列的A+T碱基含量低于控制区, 高于蛋白质编码区和rRNA区。tRNA二级结构预测结果显示, tRNA-Phe、tRNA-Val、tRNA-Leu等21个tRNA基因均形成典型的三叶草结构, 而tRNA-Ser(GCT)基因缺失二氢尿嘧啶臂(DHU臂) (图 2)。
![]()
图 2 太平洋鳕tRNA-Ser(GCT)二级结构预测图Figure 2. The second structures of the tRNA-Ser(GCT) genes in G. macrocephalus mitogenome2.5 非编码区
线粒体基因组中的非编码区是不参与编码基因, 而是调节线粒体DNA复制和转录的片段, 主要分布于L-链复制起始区和各个tRNA基因之间, 其中位于tRNA和tRNA之间的称为D-Loop区。根据结构组成和分布位置的不同, 非编码区可分为基因间隔序列区和基因序列重叠区。在太平洋鳕线粒体中共寻找到11处重叠区, 序列总长度达43 bp, 最大重叠碱基数为10 bp, 位于ATP6基因和ATP8基因之间, 最小重叠碱基数为1 bp。非编码区有12个基因间隔序列, 长度在1—74 bp不等, 其中最长的间隔位于tRNA-Thr基因和tRNA-Pro基因之间, 长达74 bp。在tRNA-Asn基因和tRNA-Cys基因之间寻找到能够启动L-链复制的序列(OL), 该序列由14 bp的茎和13 bp的环组成, 并且在5ʹ端有一段短的序列5ʹ-GCCGG-3ʹ(图 3)。
![]()
图 3 L-链复制起始区茎环结构Figure 3. Stem-loop structure of L-strand replication initiation region富含A+T的控制区在调控线粒体DNA的复制和转录中起到重要的作用。本研究中太平洋鳕控制区显示出相对不稳定性, 与哺乳动物相比只能鉴定到一个与CSB-2相对应的保守序列CSB(图 4)。在太平洋鳕D-Loop区发现与终止结合序列区(Terminal associated sequences, TAS)共有序列相似的序列且包含保守的5个核苷酸5ʹ-TACAT-3ʹ, 该序列可以形成发夹结构, 是D环DNA终止的候选位点。另外太平洋鳕D环含有一个17 bp的嘧啶序列。
![]()
图 4 太平洋鳕D-Loop区序列方框: TAS*, 下划线: Py17, 双下划线: CSBFigure 4. Schematic map characterizing of the control region of G. macrocephalusBox: TAS*, Underline: Py17, Double underline: CSB2.6 系统发育树分析
以斑马鱼和红鳍东方鲀等几种常见经济鱼类作为外群, 采用最大似然法(Maximum Likelihood, ML), 构建几种鳕科鱼类的系统进化树(图 5)。蓝鳕、北鳕、杜氏银大眼鳕、江鳕与斜带石斑鱼、斑马鱼、大菱鲆等经济鱼种聚为一大支, 而红鳍东方鲀与其他海水性鳕鱼聚为一大支。传统观点所认为的3种真鳕(太平洋鳕、格陵兰鳕、大西洋鳕)并没有汇聚成一个分支, 北极鳕与太平洋鳕和大西洋鳕两种真鳕汇成一支, 而格陵兰鳕则与阿拉斯加狭鳕以及挪威狭鳕汇成一个大支。另外江鳕虽然作为鳕科鱼类, 但进化分析结果显示, 江鳕、蓝鳕和杜氏银大眼鳕3种鳕科鱼类与其他鳕科鱼类进化关系相对疏远。
![]()
图 5 基于线粒体基因组全序列构建的ML进化树Figure 5. Phylogenetic tree based on the ML analysis of the whole mitochondrial genomeCytb基因因为存在于所有鱼类线粒体中, 且进化速度适中, 在一定的进化尺度内不会受到饱和效应的影响, 常被作为分析物种种间和属间差异的依据[15]。以Cytb基因采用邻接法构建的进化树则更符合传统分类学的观点(图 6)。在淡水环境中生活的江鳕与其他鳕科鱼类进化关系较为疏远, 但仍然与其他鳕科鱼类聚为一大支。
![]()
图 6 采用邻接法构建的鳕科鱼类Cytb基因进化树Figure 6. Phylogenetic tree based on the NJ analysis of the Cytb genes3. 讨论
通过研究太平洋鳕线粒体基因组碱基组成发现其基因组成与其他典型脊椎动物相似, 并且呈明显的AT偏向性, 其他鱼类中也同样报道了此类的AT偏向性, 只是含量因物种的差异而有所区别[16]。连总强等[16]指出大多数鱼类L-链编码了tRNA-Gln、tRNA-Ala、tRNA-Asn、tRNA-Cys、tRNA-Tyr、tRNA-Ser(UCN)、tRNA-Glu和tRNA-Pro 8个tRNA以及ND6基因, 其余tRNA均在H-链编码, 但本研究中太平洋鳕H-链只编码了tRNA-Gln、tRNA-Trp、tRNA-Ala、tRNA-Asn、tRNA-Cys、tRNA-Ser、tRNA-Glu和tRNA-Pro以及NAD6基因, 其他tRNA均在L-链编码, 似乎恰好与大多数鱼类不同。此外, tRNA-Thr和tRNA-Pro之间出现74 bp的基因间隔区, 在其他鱼类中同样存在这种现象, Christoph等[17]在鲅鲷(Petrochromis trewavasae)线粒体中发现最长为4 bp的折叠区, 与本文相比只在折叠的碱基数目方面存在差异。此间隔区未发现与H链或L链复制相关的起始序列。在tRNA-Asn基因和tRNA-Cys基因之间发现的OL序列主要功能是调节L-链的复制[18], 此序列在哺乳动物和其他鱼类线粒体控制区中都存在[19, 20]。
太平洋鳕线粒体基因组蛋白质编码基因中, 起始密码子与大多数鱼类一致, 均使用2种密码子(ATG和GTG), 但是终止密码子则与大多数硬骨鱼不同[16], 除了TAA、TAG、TAA和T--四种终止密码子外, 出现了以AGA和AGG为终止密码子的情况。其中, COX2和NAD4基因以AGA为终止密码子, NAD6基因以AGG为终止密码子。这2种密码子多在哺乳动物线粒体中使用, 而太平洋鳕线粒体中使用这两种密码子预示着太平洋鳕在进化过程中与哺乳动物进化方式更为相似。在碱基偏移方面, 除了生活在南半球的蓝鳕以及淡水类江鳕, 其他鱼类均表现出弱AT负偏移, 并且作为外群物种中的大菱鲆同样表现出弱AT偏移, 但生活在热带和亚热带海域的斜带石斑鱼却没有这种现象, 推测A+T-skew%的偏移很可能与物种生活的地理位置以及环境温度有关, 甚至可能涉及蛋白质编码的方向。Mclean等[21]通过检测9种真细菌中GC和AT偏斜情况, 发现AT偏斜与转录复制的方向有关, 这或许可以解释上文中提及到的太平洋鳕大量tRNA和基因在L-链编码的原因。同时根据Weber等[22]对黑腹果蝇(Drosophila melanogaster)内含子中AT偏斜差异的研究, 发现其可能控制着物种之间的物理距离。本文所涉及的鳕科鱼类的地理隔离是否也是由AT偏斜控制的尚需探究。
线粒体12S rRNA和16S rRNA的功能主要决定于其二级结构。太平洋鳕rRNA的序列十分保守, 与其他已报道的硬骨鱼类十分相似, 如兰州鲇[16]。本研究中22个tRNA基因中除了tRNA-Ser(GCT)外, 其余均形成典型的三叶草结构。tRNA缺少DHU臂, 在DHU臂的位置上形成一个单环结构。此种情况在鱼类线粒体中较为常见。Cheng等[23]研究证明这种缺失DHU臂的tRNA可以调整结构形态, 并不会影响其进入核糖体以及携带并转运氨基酸等功能。
鱼类线粒体控制区包含: 终止序列区(TAS)、中央保守区、H-链复制起始区(OH)等保守序列[24], 而本文中太平洋鳕线粒体控制区中仅含有一个与CSB-2相对应的保守序列CSB和一个与TAS功能类似的序列TAS*(图 4)。TAS*虽然在结构上与其他鱼类不同, 但此序列可以形成发夹结构同时含有保守的五核苷酸序列5ʹ-TACAT-3ʹ, 是D-Loop区DNA终止的候选位点[25]。太平洋鳕D-Loop区中的17 bp嘧啶序列, 是线粒体单链结合蛋白(Single strand DNA-binding protein, mtSSB)假定结合位点。这种蛋白特异性结合单链嘧啶区域, 被认为在DNA复制调节过程中起作用。Johansen等[25]在大西洋鳕中分别发现了17和27 bp的嘧啶序列, 以及位于D-Loop区5ʹ端的4个完美重复序列。太平洋鳕控制区没有发现存在于大西洋鳕控制区的27 bp的嘧啶序列以及5ʹ端的4个完美重复序列。在大西洋鳕线粒体控制区发现的重复序列只是DNA重复复制的结果, 据此推测虽然太平洋鳕线粒体基因组在控制区与其他鱼类存在较大差异, 但是在具有关键功能的序列上通常是保守的或具有相似的代替序列, 似乎太平洋鳕在进化的过程中更加倾向于功能上的完整, 而非结构上的相似。这些D-Loop序列的独特性是否是导致太平洋鳕与其他鳕鱼地理隔绝的另一个原因, 目前尚没有明确的观点。
线粒体基因组相关基因或结构已经被广泛应用于系统发育、物种分类等不同研究领域[26]。本研究以最大似然法构建线粒体全基因组进化树(图 5)并没有显示出与传统分类学观点相似的结果。其中被认为是真鳕的太平洋鳕、大西洋鳕与格陵兰鳕并没有出现在同一支中, 生存环境更加苛刻的北极鳕则取代了格陵兰鳕的位置, 并且同为江鳕的3个不同样本在该进化树中也在不同分支中。以线粒体全基因组构建的进化树可信度较差, 原因可能是这些鱼类线粒体中非编码区存在大量的变异。太平洋鳕存在长达74 bp的基因间隔区和差异较大的D-Loop区。而无论以何种方法构建进化树, 其原理都是检测各个碱基的替换率, 而非编码区的变异则会将线粒体整体进化水平夸大。鉴于太平洋鳕和其他鳕之间, 尤其是与大西洋鳕线粒体基因组有着非常高的同源性, 所以采用Cytb基因作为进化树分析的依据(图 6)。分析结果显示, 与大西洋鳕相比, 太平洋鳕与格陵兰鳕在系统进化上更为亲近, 这似乎也能解释为什么太平洋鳕和大西洋鳕在D-Loop区上有如此巨大的差异。
线粒体基因组进化速度快于物种进化速度, 另外不同地理分布有可能造成动物体DNA的变异, 如太平洋鳕和大西洋鳕2个物种, DNA序列同源性很高, 单独比较分析某一个编码基因时差别更小, 但是它们往往在生活习性等方面有着较大的区别。目前的技术发展水平很难分析这些线粒体基因组, 尤其是控制区上的微小变化所导致的结果, 进一步分析线粒体基因组的变异不仅可以分析不同物种之间的遗传关系, 还可为解决因不同地理分布而导致物种之间的变异等问题提供数据支持。
-
[1] Yang F. Animal nutrition ( second edition). Beijing: China Ag2riculture Press. 2003, 76-88 [杨凤. 动物营养学(第二版).北京: 中国农业出版社. 2003, 76-88]
[2] Feng J, J ia G. Studies on the fatty liver diseases resulted fromdifferent lip id levels in sciaenop s ocellatus diets [J]. ActaHydrobiologica S inica, 2005, 29 (1) : 61-64 [冯健, 贾刚. 饵料中不同脂肪水平诱导红姑鱼脂肪肝病的研究. 水生生物学报, 2005, 29 (1) : 61-64]
[3] Feng J, Qin ZB. Effects of four dietary lip ids on the performanceof growth and body composition in pacific salmon, oncorhynchuy [J]. Acta Hydrobiologica S inica, 2006, 30 (3) : 256-261 [冯健, 覃志彪. 4 种不同脂肪源对太平洋鲑生长和体组成的影响. 水生生物学报, 2006, 30 (3) : 256-261]
[4] Allen JM,Adrian T E, Tatemoto K, et al. Two novel related pep2tides, neuropep tide Y (NPY) and pep tide YY ( PYY) inhibitthe contraction of the electrically stimulated mouse vas deferens [J]. N europeptides, 1982, 3 (2) : 71-77
[5] Inui A. Neuropep tide Y feeding recep tors: are multip le subtypesinvolved [J] ? Trend Pharm acol Sci, 1999, 20 (2) : 43-46
[6] Lopez2PatinoM A, Guijarro A I, Isorna E, et al. Neuropep tide Yhas a stimulatory action on feeding in goldfish (Carassius aura2tus) [J]. Eur J Pham acol, 1999, 377: 147-153
[7] Narnaware Y K, Peyon P P, Lin X, et al. Regulation of food in2take by neuropep tide Y (NPY) in goldfish [J]. Am J PhysiolRegul Integr Com p Physiol, 2000, 279 (3) : 1025-1034
[8] Fleury C, Neverova M, Collins S, et al. Uncoup ling p rotein22: anovel gene linked to obesity and hyperinsulinemia [J]. N at Gen2et, 1997, 15 (3) : 269-272
[9] Zhou Y T, ShimabukuroM, Koyama K, et al. Induction by lep tinof uncoup ling p rotein22 and enzymes of fatty acid oxidation [J].Proc N atl Acad Sci USA, 1997, 94 (12) : 6386-6390
[10] Auwerx J,Leroy P, Schoonjans K. Lipoprotein lipase: recent contri2butions from molecular biology [J]. Crit Rev Clin Lab Sci, 1992,29 (3-4) : 243-268
[11] Zechner R. The tissue2specific exp ression of lipop rotein lipase:imp lications for energy and lipop rotein metabolism [J]. CurrOpin L ipidol, 1997, 8 (2) : 77-88
[12] Kirchgessner T G, Chuat J C, Heinzmann C, et al. Organization ofthe human lipop rotein lipase gene and evolution of the lipase genefamily [J]. Proc N atl Acad Sci USA, 1989, 86 ( 24 ) :9647-9651
[13] Dugi K A,Vaisman B L, SakaiN, et al. Adenovirus2mediated ex2p ression of hepatic lipase in LCAT transgenic mice [J]. L ipidRes, 1997, 38 (9) : 1822-1832
[14] Marques2Vidal P, Azema C, Collet X, et al. Hepatic lipase p ro2motes the up take of HDL esterified cholesterol by the perfused ratliver: a study using reconstituted HDL particles of defined phos2pholip id composition [J]. L ipid Res, 1994, 35 (3) : 373-384
[15] LiaoW Q, Liang X F, Wang L, et al. Molecular cloning andcharacterization of alpha2class glutathione S2transferase gene fromthe liver of silver carp, bighead carp, and other major Chinesefreshwater fishes [J]. B iochem M ol Toxicol, 2006, 20 ( 3 ) :114-126
[16] M Bienengraeber, K S Echtay, M Klingenberg. H + transport byuncoup ling p rotein (UCP21 ) is dependent on a histidine pair,absent in UCP22 and UCP23 [J]. B iochem istry, 1998, 37: 3-8
[17] Nagase I, Yoshida T, SaitoM. Up2regulation of uncoup ling p ro2teins by beta2adrenergic stimulation in L6 myotubes [J]. FEBSLett, 2001, 494 (3) : 175-180
[18] BossO, Hagen T, LowellB B. Uncoup ling p roteins 2 and 3: po2tential regulators ofmitochondrial energymetabolism [J]. D iabe2tes, 2000, 49 (2) : 143-156
[19] Li L X, Skorpen F, Egeberg K. Uncoup ling p rotein22 partici2pates in cellular defense against oxidative stress in clonal beta2cells [J]. B iochem B iophys Res Comm un, 2001, 282 (1) : 273-277
[20] Klingenberg M, Echtay K S. Uncoup ling p roteins: the issuesfrom a biochemist point of view [J]. B iochim B iophys Acta,2001, 1504 (1) : 128-143
[21] Liang X F, Lin X T, Huang F, et al. The liver uncoup ling p ro2tein 2 gene of the red sea bream ( Pagrusm ajor) with reference toits physiological function [J]. Acta Zoologica S inica, 2003, 49(1) : 110-117 [梁旭方, 林小涛, 黄芬, 尾形博. 真鲷肝脏解偶联蛋白2 (UCP2) 基因及其功能的探讨. 动物学报,2003, 49 (1) : 110-117]
下载:
计量
- 文章访问数: 1053
- HTML全文浏览量: 0
- PDF下载量: 541