鮑鳴　奚業(yè)文　海思嬈　王晨龍　陳闖將　王希春　李西雷

摘 要：通過PCR擴(kuò)增和一代測(cè)序技術(shù)獲得麥魚 （Rhinogobius sp.） 線粒體基因組全序列，并對(duì)其序列結(jié)構(gòu)進(jìn)行了分析。結(jié)果表明：1） 麥魚線粒體基因組全長(zhǎng)16 511 bp，其堿基G+C的含量為47.7%，A+T的含量為52.3%，呈現(xiàn)AT偏好性；2） 共含有蛋白質(zhì)編碼基因13個(gè)，rRNA基因2個(gè)，tRNA基因22個(gè)和D-loop區(qū)，除ND6、tRNA-Gln、tRNA-Tyr、tRNA-Ala、tRNA-Asn、tRNA-Cys、tRNA-Ser、tRNA-Glu、tRNA-Pro基因在L鏈上編碼之外，其他基因均在H鏈上進(jìn)行編碼；3） 基于線粒體基因組全序列和Cytb基因，分別構(gòu)建了鰕虎魚類中部分屬的進(jìn)化樹，其表明麥魚與波氏吻鰕虎魚、褐吻鰕虎魚、子陵吻鰕虎魚聚在一個(gè)分支上，均屬于吻鰕虎魚屬，并且其與波氏吻鰕虎魚的親緣關(guān)系最近。

關(guān)鍵詞：麥魚（Rhinogobius sp.）；線粒體基因組；序列分析；系統(tǒng)發(fā)育進(jìn)化；物種鑒定

麥魚（Rhinogobius sp.） 屬鱸形目（Perciformes），鰕虎魚亞目（Gobioidei），鰕虎魚科（Gobiidae），吻鰕虎魚屬（Rhinogobius），為安徽省池州市東至縣所特有的小型魚類之一，因盛產(chǎn)于刈麥季節(jié)而定名，與著名的廬山石魚屬于同一種類[1]。迄今對(duì)于麥魚的研究甚少，麥魚準(zhǔn)確的物種鑒定、分類地位及遺傳多樣性尚未確定，現(xiàn)僅有根據(jù)形態(tài)學(xué)特征進(jìn)行的分類研究，其養(yǎng)殖潛力尚待評(píng)估與開發(fā)[1]。

線粒體基因組數(shù)據(jù)特點(diǎn)顯著，如結(jié)構(gòu)緊密、能夠自主復(fù)制、編碼效率高、無組織特異性、無重組現(xiàn)象、進(jìn)化速率快以及遵循嚴(yán)格的母系遺傳等，在系統(tǒng)與進(jìn)化生物學(xué)和群體遺傳學(xué)等許多生物學(xué)研究中已廣泛應(yīng)用[2-9]。目前對(duì)鰕虎魚科多屬種的分類較為混亂，當(dāng)前最被接受的分類是NELSON（2006）提出的9科5亞科分類方法，中國(guó)鰕虎魚分類與此一致[10-11]。根據(jù)分布范圍已經(jīng)證實(shí)過去記載于我國(guó)的櫛鰕虎魚屬種類實(shí)際上應(yīng)歸為吻鰕虎魚屬[12]。在已有的研究中，研究者根據(jù)形態(tài)特點(diǎn)結(jié)合生態(tài)習(xí)性將麥魚鑒定為吻鰕虎魚和洞庭吻鰕虎魚[1]，伍漢霖等[13]、劉蟬馨[14]、羅智等[15]根據(jù)生物及生態(tài)學(xué)特性研究鰕虎魚類的分類地位。曲麗艷[10]通過核基因及線粒體COI基因?qū)?8種鰕虎魚進(jìn)行了系統(tǒng)發(fā)育分析，將傳統(tǒng)分類中歸為鰕虎魚亞科的睛尾蝌蚪鰕虎魚聚到小鰕虎亞科，其他大部分結(jié)果與形態(tài)研究一致。ZHONG等[16]通過線粒體基因?qū)Σㄊ衔泅y虎魚（Rhinogobius Cliffordpopei）進(jìn)行的系統(tǒng)發(fā)育分析表明，波氏吻鰕虎魚與其他吻鰕虎魚屬魚類聚在一個(gè)分支，但吻鰕虎魚屬與縞鰕虎魚屬的進(jìn)化關(guān)系尚不清楚。

本試驗(yàn)使用PCR擴(kuò)增和一代測(cè)序技術(shù)，測(cè)定了麥魚的線粒體基因組全序列，分析了基因的排序及其組成，通過與26種鰕虎魚類線粒體基因組進(jìn)行比對(duì)，構(gòu)建鰕虎魚類的系統(tǒng)進(jìn)化樹，進(jìn)一步鑒定了麥魚的分類地位，補(bǔ)充了鰕虎魚科的線粒體基因組數(shù)據(jù)，也為吻鰕虎魚屬的分子鑒定和系統(tǒng)發(fā)育研究提供了參考基礎(chǔ)。

1 材料與方法

1.1 試驗(yàn)材料

本試驗(yàn)中所用的麥魚樣品均由2022年7月采自安徽省池州市東至縣。采集樣本的肌肉組織以及肝臟組織，去離子水清洗，-80 ℃冰箱中保存?zhèn)溆谩?/p>

1.2 線粒體DNA提取

使用滅菌的剪刀切割30 g左右的肌肉組織于液氮中，-80 ℃保存待用，線粒體DNA嚴(yán)格按照上海諾倫生物醫(yī)藥技術(shù)有限公司細(xì)胞線粒體DNA抽提試劑盒說明書進(jìn)行，1%瓊脂糖凝膠電泳檢測(cè)。

1.3 線粒體基因組測(cè)序

引物設(shè)計(jì)根據(jù)麥魚種在NCBI上的近緣物種的線粒體基因序列，PCR擴(kuò)增程序?yàn)椋?4 ℃預(yù)變性30 s，94 ℃變性30 s，55 ℃退火30 s，72 ℃延伸1 min/kb，30個(gè)循環(huán)后，72 ℃延伸10 min。取10 μL PCR產(chǎn)物于1%瓊脂糖凝膠中電泳檢測(cè)，DNA純化樣品直接測(cè)序。

1.4 序列的拼接與分析

使用BioEdit進(jìn)行序列分析，刪去可信度低的測(cè)序片段，并統(tǒng)計(jì)其長(zhǎng)度，再利用NCBI在線功能BLAST與數(shù)據(jù)庫(kù)里面已經(jīng)報(bào)道過的鰕虎魚亞科進(jìn)行比對(duì)分析，并使用Contig Express進(jìn)行序列的拼接，最終補(bǔ)全序列，將麥魚全線粒體基因組序列提交至GenBank數(shù)據(jù)庫(kù)，獲得序列號(hào)為QDA22533.1—QDA22545.1。

使用NCBI在線分析Open Reading Frame Finder預(yù)測(cè)確定編碼位置。用NCBI中已有的鰕虎魚類線粒體基因組進(jìn)行比對(duì)，注釋蛋白質(zhì)編碼基因（Protein coding genes，PCGs）及D-Loop區(qū)。軟件MEGA7.0分析線粒體基因序列的堿基組成、密碼子使用情況以及堿基偏移程度。

利用在線tRNAscan-SE軟件進(jìn)行麥魚tRNA基因定位，并利用The Mfold Web Server 在線軟件有效預(yù)測(cè)麥魚線粒體tRNA的二級(jí)結(jié)構(gòu)。

1.5 系統(tǒng)進(jìn)化樹的構(gòu)建

使用 MEGA 7.0基于K-2-P模型采用自展法（bootstrap）1 000次檢驗(yàn)分支可信度，構(gòu)建鄰接法（NJ法）分子系統(tǒng)進(jìn)化樹，將在GeneBank中下載的同一屬26種鰕虎魚類的線粒體基因組全序列以及Cytb基因序列與本試驗(yàn)得到的麥魚線粒體基因組全序列（16 511 bp）進(jìn)行多重序列比對(duì)和分析，以了解鰕虎魚類之間的發(fā)育關(guān)系。使用的基因基本信息見表1。

2 結(jié)果與分析

2.1 麥魚線粒體基因組結(jié)構(gòu)與特征

麥魚線粒體基因組全長(zhǎng)16 511 bp，堿基含量分別為A （26.7%） 、T （25.6%） 、G （16.9%） 、C （30.8%） ；G+C含量為47.7%，A+T含量為52.3%，說明本試驗(yàn)中麥魚線粒體基因組呈現(xiàn)AT偏好性。麥魚線粒體全基因組序列中共含37個(gè)基因，分別為PCGs 13個(gè)、tRNA基因22個(gè)、rRNA基因2個(gè)和d-loop區(qū) 1個(gè)（見封三圖1）。ND6基因以及tRNA基因中的tRNA-Gln、tRNA-Tyr、tRNA-Ala、tRNA-Asn、tRNA-Cys、tRNA-Ser、tRNA-Glu和tRNA-Pro均在L-鏈上進(jìn)行編碼，余下的28個(gè)蛋白質(zhì)基因及tRNA基因則均在H-鏈上進(jìn)行編碼（表2） 。

麥魚線粒體基因組結(jié)構(gòu)與眾多魚類基因組結(jié)構(gòu)相同，堿基排列較緊密，但仍存在部分重疊與間隔[3-5]。其基因間隔共有11處，間隔長(zhǎng)度共計(jì)58 bp，其中間隔最大處出現(xiàn)于tRNA-Cys基因與tRNA-Asn基因之間，間隔為32 bp。而基因重疊共計(jì)出現(xiàn)6次，共計(jì)21 bp，長(zhǎng)度不一，其中重疊長(zhǎng)度最多的一處出現(xiàn)在ATP8與ATP6之間以及ND4L基因與ND4基因之間，重疊7 bp。既沒有出現(xiàn)重疊也沒有出現(xiàn)間隔處總計(jì)22處 （表2） 。

2.2 線粒體基因組蛋白質(zhì)編碼基因

本次試驗(yàn)獲得的麥魚線粒體基因組的13個(gè)蛋白質(zhì)編碼基因序列總長(zhǎng)度為11 428 bp，含3 810個(gè)密碼子 （終止密碼子除外） ，在H鏈上進(jìn)行編碼的有12個(gè)基因（ND1，ND2、COX1、COX2、ATP8、ATP6、COX3、ND3、ND4L、ND4、ND5、Cytb） ，僅ND6基因在L鏈上進(jìn)行編碼。蛋白質(zhì)編碼基因擁有兩種起始密碼子（ATG和GTG），除了COX1蛋白質(zhì)編碼基因的起始密碼子為GTG以外，其余12個(gè)蛋白質(zhì)編碼基因的起始密碼子均為ATG。終止密碼子有四種（TAG、TAA、TA-、T--），使用情況相較起始密碼子更為復(fù)雜。其中ND1、COX1、ATP8、ND4L、ND5基因的終止密碼子為TAA；ND2、ND6基因的終止密碼子為TAG；而COX2、COX3、ND3、ND4、Cytb基因則是以不完全的T作為終止密碼子；ATP6基因的終止密碼子為TA-。終止密碼子使用情況與太平洋鱈[3]、黃條鰤[4]相同。

軟件分析蛋白質(zhì)編碼基因密碼子的使用頻率（表3）以及同義密碼子使用頻率的結(jié)果表明：13個(gè)蛋白質(zhì)編碼基因中有35個(gè)偏好密碼子（RSCU密碼子），第三位為C或T的密碼子具有很高的使用頻率，以G為結(jié)尾的密碼子使用頻率最小。第三位堿基以A或者T結(jié)尾的密碼子中，UUA（L）、UCA（S）、UAU（Y）、CAU（H）、AAU（N）、GAU（D）、UGU（C）、CGU（R）、AGU（S）、GGU（G）、GGA（G）密碼子的RSCU小于1，其他密碼子RSCU均大于1。

2.3 線粒體基因組tRNA基因及rRNA基因

麥魚線粒體基因組共計(jì)22個(gè)tRNA基因，其長(zhǎng)度分布于68 bp到76 bp之間，總計(jì)1 556 bp，不均勻地分布于13個(gè)蛋白質(zhì)編碼基因之間。分析結(jié)果顯示，tRNA基因之間存在著相同程度的重疊及不同程度的間隔，tRNA-Ile與tRNA-Gln、tRNA-Thr與tRNA-Pro、tRNA-Gln與tRNA-Met分別都存在相同程度的核苷酸重疊，三組重疊的核苷酸都是1bp。其間隔的長(zhǎng)度從1 bp至32 bp不等，如tRNA-Asn基因與tRNA-Cys基因之間的間隔為32 bp。tRNA-Gln、tRNA-Asn、tRNA-Ala、tRNA-Cys、tRNA-Tyr、tRNA-Glu、tRNA-Ser、tRNA-Pro基因均在L鏈進(jìn)行編碼，其余14個(gè)tRNA基因在H鏈編碼。tRNA的二級(jí)結(jié)構(gòu)預(yù)測(cè)結(jié)果顯示（詳見圖2），tRNA-Val、tRNA-Phe等21個(gè)tRNA基因均可形成明顯的三葉草結(jié)構(gòu)，tRNA-Ser （GCT） 基因缺失二氫尿嘧啶臂 （DHU臂） 。

另外，麥魚線粒體基因組包括rRNA基因2個(gè)，12S rRNA以及16S rRNA基因，均處于H鏈上，位于tRNA-Phe與tRNA-Leu之間，由tRNA-Val基因?qū)烧叻珠_。12S rRNA基因和16S rRNA基因的A+T含量分別為52.8%和52.7%，呈現(xiàn)AT堿基偏好性。

2.4 系統(tǒng)進(jìn)化樹分析

2.4.1 線粒體全基因組NJ進(jìn)化樹分析 系統(tǒng)進(jìn)化樹（見圖3）分為兩個(gè)大支，第一個(gè)分支由大彈涂魚屬、青彈涂魚屬、副平牙鰕虎魚屬、狼牙鰕虎魚屬 、鰻鰕虎魚屬、齒彈涂魚屬、彈涂魚屬、縞鰕虎魚屬、竿鰕虎魚屬、鯔鰕虎魚屬等構(gòu)成；第二個(gè)分支由尖塘鱧屬、烏塘鱧屬、吻鰕虎魚屬構(gòu)成。

由圖3可知本次試驗(yàn)獲得的麥魚線粒體基因組和子陵吻鰕虎魚、褐吻鰕虎魚、波氏吻鰕虎魚聚在較近的分支上，與波氏吻鰕虎魚親緣關(guān)系最近。并進(jìn)一步證實(shí)過去記載于我國(guó)的櫛鰕虎魚屬種類實(shí)際上應(yīng)歸為吻鰕虎魚屬，結(jié)果與前人的分類結(jié)果相契合[17]。

2.4.2 Cytb基因NJ進(jìn)化樹分析 本分析采用當(dāng)前在物種鑒別、遺傳多樣性和群體分化中常用的Cytb基因作為分子標(biāo)記，研究麥魚的系統(tǒng)發(fā)育地位[3，18-20]。該系統(tǒng)進(jìn)化樹（見封三圖4） 主要分為三個(gè)分支，青彈涂魚屬、副平牙鰕虎魚屬、狼鰕虎魚屬、齒彈涂魚屬和彈涂魚屬聚為一個(gè)大支，縞鰕虎魚屬為一個(gè)分支，吻鰕虎魚屬、烏塘鱧屬、竿鰕虎魚屬、尖塘鱧屬和鯔鰕虎魚屬聚為一個(gè)分支，烏塘鱧屬與吻鰕虎魚屬和本次研究的麥魚距離最近，由圖4可知，根據(jù)Cytb基因構(gòu)建的進(jìn)化樹與線粒體全基因組構(gòu)建的進(jìn)化樹結(jié)果幾乎完全一致。

3 討論

通過研究獲得麥魚線粒體基因組全長(zhǎng)16 511 bp，包括22個(gè)tRNA 基因、13個(gè)蛋白編碼基因，2個(gè)rRNA 和1段 D-loop 控制區(qū)，在玻璃缺鰭鯰（Kryptopterus vitreolus）[5]中也有此類結(jié)構(gòu)。其線粒體基因組中，A+T堿基 （52.3%） 大于G+C堿基 （47.7%） 含量，呈現(xiàn)AT偏好性，這與太平洋鱈 （Gadus macrocephalus）[3]、黃條鰤（Seriola aureovittata）[4]、玻璃缺鰭鯰[5]、吻鰕虎魚屬（Rhinogobius）[21]等魚類相似。

在麥魚線粒體基因組中，蛋白質(zhì)編碼基因既存在于H鏈也存在于L鏈上。Brown等[6]報(bào)道，H鏈上的蛋白質(zhì)編碼基因不容易形成有保護(hù)的雙鏈，因該H鏈大多是水解的單鏈。與黃條鰤[4]相似，麥魚線粒體基因組中絕大多數(shù)蛋白質(zhì)編碼基因都位于H鏈上，易于發(fā)生相應(yīng)的水解和氧化，而位于L鏈上的ND6基因會(huì)穩(wěn)定很多，這也體現(xiàn)了ND6基因的差異性和重要性。

在麥魚mtDNA的tRNA基因和rRNA基因中，含有22個(gè)tRNA，它們的長(zhǎng)度分布于68 bp到76 bp之間，總計(jì)1 556 bp，其中絕大部分tRNA的二級(jí)結(jié)構(gòu)都滿足三葉草的形狀，穿插于不同的蛋白質(zhì)編碼基因與rRNA中間。僅tRNA-Ser （GCT） 基因缺失二氫尿嘧啶臂 （DHU臂），這與毛明光等[3]研究的太平洋鱈十分相似。麥魚線粒體基因組中含有2個(gè)rRNA，16S rRNA和12S rRNA，其與相鄰基因并無間隔和重疊。線粒體基因組中含有的重復(fù)序列、部分插入序列和序列等都屬于物種進(jìn)化的特征[7]。

本研究的麥魚與波氏吻鰕虎魚、褐吻鰕虎魚、子陵吻鰕虎魚聚在一個(gè)分支上，均屬于吻鰕虎魚屬，并且其與波氏吻鰕虎魚的親緣關(guān)系最近。此外，吻鰕虎魚屬與中華烏塘鱧、尖塘鱧、斑駁尖塘鱧所聚成的塘鱧屬比較接近，這與曲麗艷[10]、孫希福[18]、金逍逍等[22]、ZHONG等[16]等的分析結(jié)果基本相似。本次試驗(yàn)的麥魚不管是線粒體基因組全序列進(jìn)化樹，還是使用Cytb基因構(gòu)建的進(jìn)化樹分析，都出現(xiàn)了烏塘鱧屬、吻鰕虎魚屬與麥魚的距離最近的結(jié)果，但是對(duì)于其他分析結(jié)果有所不同，縞鰕虎魚屬在線粒體基因組分析中和其他屬歸為一支，在使用Cytb基因進(jìn)行進(jìn)化樹分析時(shí)，單獨(dú)地劃為一個(gè)分支。產(chǎn)生這樣的結(jié)果是因?yàn)榧棺祫?dòng)物線粒體基因的不同區(qū)域的進(jìn)化速率不同，所以線粒體各個(gè)基因的系統(tǒng)發(fā)育信息的表現(xiàn)也會(huì)有所不同。

目前，線粒體基因組中的相關(guān)基因被大量應(yīng)用于系統(tǒng)發(fā)育、物種分類等相關(guān)研究領(lǐng)域[2]。KIM等[23]、黃小林等[8]基于mtDNA中的12S rRNA、ND4L基因構(gòu)建系統(tǒng)進(jìn)化樹研究矛尾刺鰕虎魚（Acanthogobius hasta）、籃子魚（Siganidae Siganus）的進(jìn)化地位。隨著分子生物學(xué)手段的普遍使用，國(guó)內(nèi)外學(xué)者對(duì)鰕虎魚科魚類系統(tǒng)發(fā)育研究也越來越多，AGORRETA等[24]分析了222種歐洲的鰕虎魚目魚類在系統(tǒng)進(jìn)化上的關(guān)系，THACKER等證實(shí)了絲鰭鱚屬是鰕虎魚科的姐妹屬[25]，研究了鰕虎魚科在棘鰭魚目中的系統(tǒng)分布[26]。普遍來說線粒體基因的進(jìn)化速度大于物種的進(jìn)化速度，不同的地理環(huán)境、生活習(xí)慣也都有機(jī)會(huì)導(dǎo)致物種DNA的變異，要解釋清楚還需要進(jìn)一步加深對(duì)于線粒體基因的研究。

參考文獻(xiàn)：

[1] 安徽省池州地區(qū)升金湖水產(chǎn)資源調(diào)查小組.池州地區(qū)的麥魚[J].動(dòng)物學(xué)雜志，1976，4（19）：36-37.

[2] 李雪娟，楊婧，王俊紅，等.線粒體基因組數(shù)據(jù)等分析方法和軟件[J].應(yīng)用昆蟲學(xué)報(bào)，2013，50（1）：298-304.

[3] 毛明光，顧杰，劉瑞婷，等.太平洋鱈線粒體全基因組測(cè)序及結(jié)構(gòu)特征分析[J].水生生物學(xué)報(bào)，2019，43（1）：17-26.

[4] 史寶，柳學(xué)周，劉永山，等.黃條鰤線粒體全基因組測(cè)序及結(jié)構(gòu)特征分析[J].中國(guó)水產(chǎn)科學(xué)，2019，26（3）：405-415.

[5] 孫志鵬，呂偉華，匡友誼，等.玻璃缺鰭鯰線粒體全基因組序列測(cè)定與分析[J].水產(chǎn)學(xué)雜志，2018，31（2）：12-18.

[6] BROWN W M，PRAGER E M，WANG A，et al.Mitochondrial DNA sequences of primates：tempo and mode of evolution[J].Journal of Molecular Evolution，1982，18（4）：225-239.

[7] 高祥剛，鹿志創(chuàng)，田甲申，等.斑海豹線粒體基因組序列分析及比較研究[J].水產(chǎn)科學(xué)，2017，36（4）：5.

[8] 黃小林，楊育凱，林黑著，等.籃子魚屬線粒體基因組序列系統(tǒng)發(fā)育信息分析[J].生物學(xué)雜志，2018，35（5）：33-36.

[9] 楊楊，宋小晶，唐文喬，等.克氏光唇魚線粒體基因組測(cè)定及光唇魚屬的系統(tǒng)發(fā)育分析[J].動(dòng)物學(xué)雜志，2018，53（2）：207-219.

[10] 曲麗艷.基于形態(tài)、線粒體及核基因?qū)|海28種鰕虎魚類的系統(tǒng)進(jìn)化分析[D].上海：上海海洋大學(xué)，2018：42-47.

[11] NELSON J S.Fish of the world[M].4th ed.New York：John Wiley and Sons， 2006：601.

[12] 郭昶暢.中國(guó)沿海石首魚科魚類的鑒定、分類和分子系統(tǒng)發(fā)育研究[D].廈門：廈門大學(xué)，2017.

[13] 伍漢霖，鐘俊生.中國(guó)動(dòng)物志-硬骨魚綱、鱸形目（五）、鰕虎魚亞目[M].北京：科學(xué)出版社，2008.

[14] 劉蟬馨.紋縞鰕虎魚的生物學(xué)和生態(tài)學(xué)[J].水產(chǎn)科學(xué)，1985（4）：13-17.

[15] 羅智，李曉東，白海娟，等.野生和養(yǎng)殖矛尾復(fù)鰕虎魚營(yíng)養(yǎng)組成和形態(tài)學(xué)的比較研究[J].上海水產(chǎn)大學(xué)學(xué)報(bào)，2008，17（2）：182-186.

[16] ZHONG L Q，WANG M H，LI D，et al.Complete mitochondrial genome of freshwater goby Rhinogobius cliffordpopei（Perciformes，Gobiidae）：genome characterization and phylogenetic analysis[J].Genes & Genomics，2018，40：1137-1148.

[17] 李帆，鐘俊生.中國(guó)浙江省吻鰕虎魚屬一新種 （鱸形目：鰕虎魚科） [J]．動(dòng)物學(xué)研究，2007，28（5）：539-544.

[18] 孫希福.基于形態(tài)學(xué)和分子生物學(xué)資料探討中國(guó)沿海10種鰕虎魚類的系統(tǒng)發(fā)育關(guān)系[D].青島：中國(guó)海洋大學(xué)，2009.

[19] 陳雷，程遠(yuǎn)志，王日昕，等.2種鰕虎魚細(xì)胞色素b基因全序列克隆與分析[J].南方水產(chǎn)，2010，6（5）：1673-2227.

[20] AKIHITO，IWATA A，KOBAYASHI T，et al.Evolutionary aspects of gobioid fishes based upon a phylogenetic analysis of mitochondrial cytochrome b genes[J].GENE.2000，259（1-2）：5-15.

[21] WANG D，DAI C X，LI Q，et al.Complete mitochondrial genome and phylogenic analysis of Rhinogobius cliffordpopei（Perciformes，Gobiidae）[J].Taylor & Francis，2019，4（2）：2473-2474.

[22] 金逍逍，孫悅娜，王日昕，等.鰕虎魚類線粒體全基因組序列結(jié)構(gòu)特征分析及系統(tǒng)發(fā)育關(guān)系探討[J].遺傳.2013，（12）：1391-1402.

[23] KIM I C，KWEON H S，Kim Y J，et al.The complete mitochondrial genome of the javeline goby Acanthogobius hasta （Perciformes，Gobiidae） and phylogenetic considerations[J].GENE.2004，336（2）：147-153.

[24] AGORRETA A，MAURO D S，SCHLIEWEN U，et al.Molecular phylogenetics of Gobioidei and phylogenetic placement of European gobies[J].Molecular Phylogenetics and Evolution，2013，69：619-633.

[25] THACKER C E，SATOH T P，KATAYAMA E，et al.Molecular phylogeny of Percomorpha resolves Trichonotus as the sister lineage to Gobioidei （Teleostei：Gobiiformes） and confirms the polyphyly of Trachinoidei[J].Molecular Phylogenet Evolution，2015：172-179.

[26] THACKER C E.Phylogeny of Gobioidei and placement within Acanthomorpha，with a new classification and investigation of diversification and character evolution[J].Copeia，2009（1）：93-104.

The analysis of complete mitochondrial genome sequences and gene organization of Rhinogobius sp.（Perciformes，Gobiidae）

BAO Ming1，XI Yewen1，HAI Sirao2，WANG Chenlong2，CHEN Chuangjiang2，WANG Xichun2，LI Xilei2

（1.Aquatic Technology Promotion Station of Anhui，Hefei 230601，China;2.College of Animal Science and Technology，Anhui Agricultural University，Hefei 230036，China）

Abstract：The PCR and the first-generation sequencing technology were used to obtain the complete mitochondrial genome （mtDNA） sequences of Rhinogobius sp..The mtDNA was annotated and its structure was analyzed.In present study，the circular genome was 16 511 bp in length.The G+C content was 47.7% and the A+T content was 52.3%，which showed a tendency of high AT.The mitochondrial genome encodes 13 proteins and 2 rRNAs，22tRNAs and a D-Loop region.Except for ND6，tRNA-Gln，tRNA-Tyr，tRNA-Ala，tRNA-Asn，tRNA-Cys， tRNA-Ser，tRNA-Glu，tRNA-Pro，all the other genes were encoded on the heavy strand （H strand） .Based on the whole mitochondrial genome sequence and the Cytb gene，the pthlogenetic trees were constructed separately with several species in Gobiidae.It showed that Rhinogobius sp.clustered with R. cliffordpopei，R.brunneus and R.giurinus on a branch，belonging to Rhinogobius and having the closest relationship with R. cliffordpopei.

Key words：Rhinogobius sp.; mitochondrial genome; sequence analysis; phylogenetic evolution; species identification

（收稿日期：2023-05-17）

作者簡(jiǎn)介：鮑鳴（1979.7-），女，研究生，工程師，研究方向：水產(chǎn)品質(zhì)量安全、漁業(yè)環(huán)境監(jiān)測(cè)等。E-mail： baominglj@126.com。