高原林蛙线粒体全基因组序列测定及系统发育分析.pdf
《高原林蛙线粒体全基因组序列测定及系统发育分析.pdf》由会员分享,可在线阅读,更多相关《高原林蛙线粒体全基因组序列测定及系统发育分析.pdf(8页珍藏版)》请在咨信网上搜索。
1、doi:10.3969/j.issn.2095-1736.2023.04.054高原林蛙线粒体全基因组序列测定及系统发育分析张湑泽1,2,董 豹1,哈金强1,赵 英1,魏生楠1,韦迎婷1(1.青海民族大学 生态环境与资源学院,西宁 810007;2.青海省特色经济植物高值化利用重点实验室,西宁 810007)摘 要 利用二代测序技术 Illumina 平台对高原林蛙(Rana kukunoris)线粒体进行低覆盖全基因组测序,对线粒体基因组进行获取、拼接与注释,获得高原林蛙线粒体全基因组序列,序列全长 21 913 bp,主要包括 tRNA 基因(22 个)、rRNA 基因(2 个)、CDS(
2、13 个)和 D-loop(1 个)4 部分。除轻链中的 NAD6 外,其他蛋白质编码基因都在同一条链上。在编码基因中,共发现 3 个重叠基因:COX1/trnS2、ATP6/ATP8 和 NAD4L/NAD4。高原林蛙线粒体基因组具有 AT 偏好性,蛋白质编码基因中氨基酸使用最频繁的是 Leu 和 Ser,tRNA 基因都有经典的三叶草结构。KaKs分析显示,高原林蛙线粒体基因组编码基因未检测到正选择作用。基于高原林蛙线粒体全基因组构建的系统发育树表明,高原林蛙与中国林蛙亲缘关系最近。研究结果将为高原林蛙的起源、分化和遗传多样性研究提供数据支持。关键词 高原林蛙;线粒体全基因组;序列分析;系
3、统发育;物种分化中图分类号 Q951+.3文献标识码 A文章编号 2095-1736(2023)04-0054-08Complete mitochondrial genome sequence and phylogenetic analysis of Rana kukunorisZHANG Xuze1 2 DONG Bao1 HA Jinqiang1 ZHAO Ying1 WEI Shengnan1 WEI Yingting1 1.College of Eco-Environment and Resources Qinghai Minzu University Xining 810007 Ch
4、ina 2.Qinghai Provincial Key Laboratory of High Value Utilization of Characteristic Economic Plants Xining 810007 China Abstract The high-throughput sequencing of illumina platform was used to sequence the whole mitochondrial genome of Rana kukuno-ris.Then the mitochondrial genome was assembled and
5、annotated.The total length of the sequence was 21 913 bp which mainly in-cluded four parts tRNA gene 22 rRNA gene 2 CDS 13 and D-loop 1.Except NAD6 in the light chain other protein cod-ing genes were on the same chain.Of the coding genes three overlapping genes were found COX1/trnS2 ATP6/ATP8 and NA
6、D4L/NAD4.Mitochondrial genome of Rana kukunoris had AT preference.Leu and Ser were the most frequently used amino acids in protein coding genes and all tRNA had the classic clover structure.KaKs analysis showed that the coding genes of mitochondrial genome were not significantly positively selected
7、in Rana kukunoris.The phylogenetic tree constructed based on the whole mitochondrial genome showed that Rana kukunoris was closely related to Rana chensinensis.It provided the genetic information of Rana kukunoris mitochondri-al genome and aid future investigation of phylogenetics speciation and evo
8、lution of Rana genus.Keywords Rana kukunoris mitochondrial genome sequence analysis phylogeny species differentiation 高原林蛙(Rana kukunoris),隶属于无尾目(Anu-ra)蛙科(Ranidae)林蛙属(Rana),是青藏高原特有种1。高原林蛙曾被认为是中国林蛙的亚种,属于中国林蛙复合体。谢锋等2将中国西北地区的林蛙种群分为两个有效种,即中国林蛙和高原林蛙,以 Nikol skii 等于 1918 年在青海湖采集的高原林蛙(Rana amu-45第 40 卷第 4
9、 期2023 年 8 月生 物 学 杂 志JOURNAL OF BIOLOGY Vol.40 No.4Aug.2023收稿日期:2022-04-23;最后修回日期:2022-06-30基金项目:青海省自然科学基金项目(2020-ZJ-965Q);青海民族大学校级规划项目(2021XJGH15);青海省动物生态基因组学重点实验室开放课题项目(QHEG-2019-04)作者简介:张湑泽,硕士,副教授,研究方向为高原两栖动物多样性与进化,E-mail:zxz1904 rensis kukunoris Nikol skii 1918)标本为模式标本。江建平等3采用 Cytb 基因对高原林蛙不同居群之间
10、的遗传多样性进行分析,基于分子生物学证据确定高原林蛙是林蛙属单独物种,染色体 2n=24。Chen 等4发现高原林蛙存在两性异型现象;Zhou 等5认为青藏高原隆升和气候环境变化促进中国林蛙复合体物种形成,且中国林蛙部分群体与高原林蛙存在渐渗杂交现象;基于 Cytb 基因的高原林蛙群体进化分析,发现青藏高原东北部高原林蛙存在多个冰期避难所,并在祁连山南北形成两个高度分化的世系6;米志平等7研究了高原林蛙皮肤结构,发现其皮肤形态特征与其他物种不同,推测可能因为生活在高海拔地区,生存环境气温低,水分不易蒸发;张湑泽等8对青藏高原东北部高原林蛙皮肤结构观察发现,高原林蛙对青藏高原环境有明显的适应性;
11、Yang 等9通过林蛙属 3 个物种线粒体全基因组序列,分析蛙科的进化关系。Wang等10对来自贺兰山的高原林蛙的线粒体全基因组开展测序,进行系统发育分析。但以上研究还未对高原林蛙线粒体基因组具体编码信息、基因结构、遗传多样性和分子进化开展分析。线粒体 DNA 不仅具有结构简单、无重复、母系遗传和进化速率快等特点11-13,还具有很高的专一性和独特性,是动物系统发育分析和物种鉴定的重要工具14。探究线粒体结构可以为物种起源、分化和系统发育分析提供良好的参考。如今越来越多的学者通过对线粒体基因组的研究分析物种间的进化关系15-17。本研究基于高原林蛙线粒体全基因组,对其组成和结构进行分析,以期对
12、高原林蛙遗传多样性、分类地位和物种保护进行深入研究。1 材料与方法1.1 材料高原林蛙样本于 2020 年 6 月采自青藏高原东北部青海湖周边黑马河镇(363726.63 N,1000720.82 E,3 194 m)。样本采集后保存于 95%乙醇中,采用标准苯酚/氯仿法从高原林蛙肌肉组织中提取基因组DNA,采用 Illumina NovaSeq 平台进行测序,测序读长为 PE150,测序获得 5G 数据量,使用 fastp(version 0.20.0,https: 据 进 行 过 滤,获 得 Clean Data 总 碱 基 数 为6 972 984 300 bp,用 SPAdes 组装成
13、一个完整的高原林蛙线粒体基因组。线粒体基因组全长为 21 913 bp,该序列已经提交 NCBI 基因数据库(GenBank 登录号:MZ043820)。另有下载自 GenBank 的 25 份其他物种的线粒体全基因组序列,并参照 NCBI 上所提交的高原林蛙线粒体全基因组序列(MN733918 和 KU246049)开展分析。实验样本信息见表 1。表 1 用于系统进化分析的物种线粒体 DNA 序列信息Table 1 Mitochondrial DNA sequence of species forphylogenetic analysisGenBank登录号GenBankaccession
14、No.线粒体基因组大小Mt DNA size/bp物种SpeciesNC02352918 808中国林蛙 Rana chensinensisKF84051622 255昆嵛林蛙 Rana kunyuensisMK48311820 120峨眉林蛙 Rana omeimontisKP22228117 343北美洲木蛙 Rana sylvaticaMW52698918 291大别山林蛙 Rana dabieshanensisKT58807119 253桓仁林蛙 Rana huanrensisMF37034820 571黑龙江林蛙 Rana amurensisNC02352818 864东北林蛙 Ra
15、na dybowskiiKM59055015 765韩国林蛙 Rana coreanaKP01311017 805加利福尼亚红腿蛙 Rana draytoniiMN21868719 205镇海林蛙 Rana zhenhaiensisKU24604818 591昭觉林蛙 Rana chaochiaoensisMN73391816 644高原林蛙 Rana kukunorisKU24604918 863高原林蛙 Rana kukunorisKX68610817 681牛蛙 Rana catesbeianaNC04222616 061欧洲林蛙 Rana temporariaKU72030017 21
16、1比利斯牛蛙 Rana pyrenaicaMZ57136517 873越南趾沟蛙 Rana johnsiKP01309617 504佛罗里达沼泽蛙 Rana okaloosaeMT78212119 207桑植蛙 Rana sangzhiensisNC03232817 656中亚林蛙 Rana asiaticaMW00906717 370上野蛙 Rana uenoiMZ68052919 395寒露林蛙 Rana hanluicaMZ68052817 833长肢林蛙 Rana longicrusOL46732117 779武夷林蛙 Rana wuyiensisKP31748217 837高山倭蛙
17、Nanorana parkeriHQ32423217 660倭蛙 Nanorana pleskei1.2 方法1.2.1 全序列特征分析采用 SPAdes v 3.10.1(http:cab.spbu.ru/soft-ware/spades/)软件组装线粒体基因组,组装不依赖参考基因组。使用 Perl 脚本进行计算,包括碱基组成与偏移分析、密码子偏好性。使用 OGDRAW(https:chlorobox.mpimp-golm.mpg.de/OGDraw.html)制作线粒体基因组图谱。使用 Mitos2(http:mitos2.bioinf.uni-leipzig.de)对组装好的序列进行注释
18、,获得 tRNA二级结构。55第 40 卷第 4 期2023 年 8 月生 物 学 杂 志JOURNAL OF BIOLOGY Vol.40 No.4Aug.20231.2.2 高原林蛙线粒体基因组编码基因 KaKs 分析使用 mafft v 7.310(https:mafft.cbrc.jp/align-ment/software/)软件进行基因序列的比对18,使用KaKsCalculator v 2.0(https: 值19。1.2.3 遗传距离分析基于 Kimura 2-parameter 碱基替代模型,使用软件 MEGA11.020计算高原林蛙与其他林蛙属之间的遗传距离。1.2.4 系
19、统进化分析基于测序获得的高原林蛙线粒体全基因组序列和GenBank 上下载的 25 条林蛙属物种线粒体全基因组序列以及 2 条倭蛙属线粒体全基因组序列开展系统发育分析,物种间序列用 MAFFT 软件(v7.427,auto 模式)进行多序列比对,将比对好的数据用 BioEdit 编辑21,使用 RAxMLv8.2.1022(https:cme.h-its.org/exelixis/software.html)和 jModeltest 2.1.723软件进行分析,选用 GTR+GAMMA 模型,构建最大似然进化树,Bootstrap=1 000。2 结果与分析2.1 高原林蛙线粒体全基因组序列分
20、析高原林蛙线粒体全基因组序列总长 21 913 bp,包括 22 个 tRNA 基因、2 个 rRNA 基因、13 个蛋白质编码基因和 1 个 D-loop 区。其中,蛋白质编码基因序列全长11 287 bp,tRNAs 总长 1 530 bp,rRNAs 总长 2 503 bp。高原林蛙线粒体全基因组序列中 A、T、G、C 的碱基含量分别为 28.41%、30.17%、14.28%和 27.13%(图 1),A+T 的碱基含量为 58.58%,G+C 的碱基含量为 41.42%,表现出明显的 A+T 偏好性。不同物种在其基因组上各有差异,表 2 所示为高原林蛙线粒体 DNA 上的 D-loo
21、p 区、RNA 基因、蛋白质编码基因所在位置、蛋白质编码基因的密码子、起始密码子和终止密码子。表 2 高原林蛙线粒体基因组结构Table 2 Mitochondrial genome structure of Rana kukunoris基因Gene位置Location长度Length/bp结构区位置Strand密码子Codon反密码子Anticodon起始密码子Star codon终止密码子Stop codonD-loop16 3236 323+trnL16 3246 39572+TAGtrnT6 3986 46770+TGTtrnP6 4686 53669-TGGtrnF6 5386 60
22、770+GAArrnS6 6087 535928-trnV7 5377 60569+TACrrnL7 6089 1821 575+trnL29 1829 25473+TAANAD19 25510 212958+ATT/TAAATTTAAtrnI10 21310 28371+GATtrnQ10 28410 35471-TTGtrnM10 35410 42269+CATNAD210 42311 4551 033-ATG/TATATGTATtrnW11 45611 52570+TCAtrnA11 52611 59570-TGCtrnN11 59611 66873-GTTtrnC11 69711 76
23、165-GCAtrnY11 76211 82867-GTACOX111 83313 3831 551+ATA/AGGATAAGGtrnS213 37513 44571-TGAtrnD13 44713 51567+GTCCOX213 51614 203688+ATG/CCTATGCCTtrnK14 20414 27269+TTTATP814 27414 435162+ATG/TAAATGTAAATP614 42915 110682+ATG/CTTATGCTTCOX315 11115 894784+ATG/CATATGCATtrnG15 89515 96268+TCCNAD315 96316 30
24、2340+ATG/AGTATGAGT65第 40 卷第 4 期2023 年 8 月生 物 学 杂 志JOURNAL OF BIOLOGY Vol.40 No.4Aug.2023续表 2(Continued Table 2)基因Gene位置Location长度Length/bp结构区位置Strand密码子Codon反密码子Anticodon起始密码子Star codon终止密码子Stop codontrnR16 30316 37169+TCGNAD4L16 37216 656285+GTG/TAAGTGTAANAD416 65018 0091 360+ATG/TCTATGTCTtrnH18 01
25、018 07869+GTGtrnS118 07918 14567+GCTNAD518 14819 9531 806+ATG/AGGATGAGGNAD620 20420 698495-ATG/CATATGCATtrnE20 69920 76769-TTCCOB20 77121 9131 143+ATG/TAAATGTAA图 1 高原林蛙线粒体基因组序列结构图谱Figure 1 Sequence map of mitochondrial genome of Rana kukunoris2.2 高原林蛙蛋白质编码基因分析高原林蛙线粒体基因组共有13个蛋白质编码基因,包括 1 个 Cytb 基因,2
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 高原 线粒体 基因组 序列 测定 系统发育 分析
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。