西氏熊蜂線粒體基因組結(jié)構(gòu)特征及系統(tǒng)發(fā)育分析

摘要：本研究利用二代測(cè)序技術(shù)對(duì)西氏熊蜂（Bombus sichelii）的線粒體基因組進(jìn)行了測(cè)序、組裝和注釋，并對(duì)其線粒體基因組的組成、結(jié)構(gòu)特征、基因重排及系統(tǒng)發(fā)育等方面進(jìn)行了生物信息學(xué)分析。結(jié)果表明：西氏熊蜂線粒體基因組長度為15 711 bp，包含13個(gè)蛋白編碼基因、22個(gè)tRNA基因和2個(gè)rRNA基因；線粒體基因組全序列、蛋白編碼基因的串聯(lián)序列、tRNA串聯(lián)序列、rRNA串聯(lián)序列的A+T含量均大于80%，表現(xiàn)出明顯的堿基偏倚；13 個(gè)蛋白編碼基因的起始密碼子均為ATN，終止密碼子為TAA，TAG或T --。tRNA二級(jí)結(jié)構(gòu)預(yù)測(cè)發(fā)現(xiàn)，trnS1與其他昆蟲具有相同的特征，即缺少DHU臂和DHU環(huán)，其余21個(gè)tRNA均能形成典型的三葉草結(jié)構(gòu)；系統(tǒng)發(fā)育進(jìn)化分析結(jié)果高度支持西氏熊蜂與火紅熊蜂（B. pyrosoma）隸屬于黑熊蜂亞屬（Melanobombus）。本研究首次獲得西氏熊蜂線粒體基因組序列，初步鑒定了其線粒體基因組的特征，為進(jìn)一步對(duì)熊蜂線粒體基因組學(xué)和系統(tǒng)發(fā)育關(guān)系的研究奠定了基礎(chǔ)。

關(guān)鍵詞：西氏熊蜂；線粒體基因組；生物信息學(xué)；系統(tǒng)發(fā)育分析

中圖分類號(hào)：S891；Q963""" 文獻(xiàn)標(biāo)識(shí)碼：A"""" 文章編號(hào)：1007-0435（2024）08-2408-11

The Mitochondrial Genome Characteristics and Phylogenetic Analysis

of Bombus sichelii

MA Xiao-xuan， LIANG Cheng-bo， LIU Dao-xin， SUN Guo， YAN Jing-yan*

（School of Agriculture and Animal Husbandry， Qinghai University， Xining， Qinghai Province 810016， China）

Abstract：In this study，the mitochondrial genome of Bombus sichelii was sequenced by next-generation sequencing technology，then was assembled and annotated. And the composition，structural characteristics，gene rearrangement and phylogeny of the mitochondrial genome were analyzed by bioinformatics. The results showed that the mitochondrial genome of B. sichelii was 15 711 bp in length，containing 13 protein-coding genes，22 tRNA genes and 2 rRNA genes. The A+T content of the whole mitochondrial genome sequence，the tandem sequence of protein-coding genes，the tandem sequence of tRNA，and the tandem sequence of rRNA were all greater than 80%，showing obvious base bias. The initiation codon of 13 protein-coding genes was ATN，and the stop codons was TAA，TAG，or T--. The secondary structure prediction of tRNAs showed that trnS1 had the same characteristics as other insects. That is，it lacked DHU arm and DHU ring，and the other 21 tRNAs could form a typical clover structure. Phylogenetic findings support the affiliation of B. sichelii and B. pyrosoma to Melanobombus. In this study，the mitochondrial genome sequence of B. sichelii was obtained for the first time，and the characteristics of the mitochondrial genome of Bumblebee were preliminarily identified，which laid a foundation for further research on the relationship between mitochondrial genomics and phylogeny of bumblebees.

Key words：Bombus sichelii;Mitochondrial genome;Bioinformatics;Phylogenetic analysis

熊蜂（Bombus spp.）是膜翅目（Hymenoptera）、細(xì)腰亞目（Apocrita）、蜜蜂科（Apidae）、熊蜂屬（Bombus）昆蟲的總稱［1］。熊蜂是眾多野生植物的重要傳粉者，在高海拔生態(tài)系統(tǒng)中，鳥類和蝴蝶的傳粉作用逐漸減弱，熊蜂則由于其對(duì)高海拔環(huán)境的高度適應(yīng)性、較強(qiáng)的飛行能力和高的花粉傳遞效率成為主導(dǎo)的傳粉者，對(duì)維持高海拔生態(tài)系統(tǒng)的平衡和穩(wěn)定至關(guān)重要［2-4］。但由于棲息地喪失、環(huán)境污染和氣候變化等多種因素的影響，許多生境內(nèi)的熊蜂種類及數(shù)量顯著降低，目前已經(jīng)有36種熊蜂被列入國際自然保護(hù)聯(lián)盟瀕危物種紅色名錄［5-8］。我國在近些年的調(diào)查中發(fā)現(xiàn)，黑足熊蜂（B. atripes）已經(jīng)從云南省內(nèi)消失［9］，所以熊蜂物種多樣性的保護(hù)已迫在眉睫。

線粒體基因組（Mitochondrial genome，mtDNA）是在線粒體內(nèi)發(fā)現(xiàn)的小型環(huán)狀染色體，具有基因進(jìn)化速率快、母系遺傳、高拷貝數(shù)、結(jié)構(gòu)簡單及易獲得等特點(diǎn)，已廣泛應(yīng)用于物種鑒定、種群遺傳結(jié)構(gòu)、系統(tǒng)發(fā)育和系統(tǒng)地理學(xué)等研究領(lǐng)域［10-11］。熊蜂的傳統(tǒng)形態(tài)學(xué)鑒定主要以體色圖案為依據(jù)，然而，一些不同的熊蜂物種常常出現(xiàn)相似的毛色配置，而同種內(nèi)不同級(jí)型或不同區(qū)域的個(gè)體間毛色也會(huì)存在較大變異，這些都給熊蜂的形態(tài)分類和鑒定帶來了較大的困難［12］。隨著分子生物學(xué)技術(shù)的快速發(fā)展，基于線粒體基因序列等分子標(biāo)記的熊蜂分類學(xué)和遺傳學(xué)研究日益受到研究人員的關(guān)注。其中，COI，Cytb，16S rRNA和微衛(wèi)星標(biāo)記等是熊蜂中普遍使用的分子標(biāo)記，Bertsch［13］和劉蘋等［14］利用COI基因明確了明亮熊蜂復(fù)合種的分類地位及其分布特性；Pirounakis等［15］利用Cytb基因分析發(fā)現(xiàn)B. pascuorum在阿爾卑斯山的13個(gè)不同地理分區(qū)中存在遺傳變異。但由于單一基因所攜帶的遺傳信息較少，所以多基因聯(lián)合分析在熊蜂研究中日益增加。例如Williams等［16］利用16S rRNA和COI 2種分子標(biāo)記，對(duì)山熊蜂亞屬（Alpinobombus）中基于形態(tài)學(xué)鑒定存在問題的物種進(jìn)行了綜合鑒定；Cwjas等利用線粒體基因標(biāo)記和微衛(wèi)星標(biāo)記，研究了伊比利亞半島上地熊蜂（B. terrestris）的時(shí)空格局［17］。隨著高通量測(cè)序技術(shù)的發(fā)展，線粒體基因組序列的獲得變得越來越容易，2023年通過對(duì)山東和遼寧2個(gè)地區(qū)的日光蜂（Aphelinus mali）［18］進(jìn)行線粒體基因組測(cè)序分析，發(fā)現(xiàn)2個(gè)支系之間僅存在歷史地理隔離，尚未出現(xiàn)種群間的生殖隔離，證明了2個(gè)支系仍然屬于同一物種。所以在分子技術(shù)不斷發(fā)展的情況下，利用分子手段比利用形態(tài)學(xué)鑒定熊蜂更加高效準(zhǔn)確。近年來，通過眾多學(xué)者不斷豐富線粒體基因的研究，線粒體基因在物種鑒定、遺傳進(jìn)化等方面的應(yīng)用價(jià)值正逐漸增強(qiáng)。

西氏熊蜂隸屬于黑熊蜂亞屬，是目前眾多學(xué)者普遍認(rèn)同的分類結(jié)果，且在我國大部分省份都有分布［2］。2023年胥勛露等［7］在對(duì)青藏高原東北部邊緣的調(diào)查研究中，發(fā)現(xiàn)西氏熊蜂是此地區(qū)的常見種，廣泛分布在 3100～3600 m的高海拔地區(qū)，并且西氏熊蜂的訪花植物種類涉及范圍達(dá)到5科11種，其中對(duì)菊科（Asteraceae）植物、豆科（Fabaceae）植物等具有明顯的訪問偏好。近年來，對(duì)熊蜂的研究工作主要是集中在形態(tài)學(xué)鑒定、蜜源植物的調(diào)查、種質(zhì)資源的保護(hù)及人工繁育等方面，但對(duì)線粒體基因組方面的分子研究尚不完備，其中黑熊蜂亞屬僅有火紅熊蜂（B. pyrosoma）和B. lapidarius的線粒體基因組數(shù)據(jù)被記錄。鑒于此，本研究首先補(bǔ)充西氏熊蜂線粒體基因組數(shù)據(jù)，即采用二代測(cè)序技術(shù)對(duì)西氏熊蜂的線粒體基因組測(cè)序、組裝、注釋，然后對(duì)其線粒體基因組的組成和結(jié)構(gòu)特征等方面進(jìn)行分析，豐富熊蜂線粒體基因組的基礎(chǔ)研究;最后采用最大似然法（Maximum likelihood，ML）和鄰接法（Neighbor-Joining，NJ）構(gòu)建基于13個(gè)亞屬熊蜂線粒體基因組全序列的進(jìn)化樹，探討西氏熊蜂在熊蜂屬內(nèi)的分類地位及各熊蜂亞屬間的進(jìn)化關(guān)系，從而為后續(xù)的群體遺傳學(xué)研究提供重要的基礎(chǔ)信息。

1 材料與方法

1.1 標(biāo)本采集

西氏熊蜂樣品2021年7月在青海省玉樹市（33°11′ N，97°24′ E）采集，浸泡于無水乙醇中并置于-80℃冰箱內(nèi)，保存在青海大學(xué)農(nóng)牧學(xué)院昆蟲學(xué)實(shí)驗(yàn)室。

1.2 線粒體基因組的測(cè)序、組裝及注釋

西氏熊蜂樣品送至北京諾禾致源生物科技有限公司構(gòu)建350 bp的小片段測(cè)序文庫并在Illumina NovaSeq平臺(tái)進(jìn)行高通量測(cè)序，測(cè)序讀長為PE150。利用Getorganelle軟件對(duì)公司測(cè)得的clean reads進(jìn)行組裝，獲得西氏熊蜂的線粒體基因組序列［19］。采用MITOS在線服務(wù)器（http：//mitos.Bioinf.uni-leipzig.de）對(duì)西氏熊蜂線粒體基因組序列進(jìn)行功能注釋。將注釋結(jié)果導(dǎo)入Geneous 9.0.2對(duì)比參考序列，手動(dòng)調(diào)整蛋白編碼基因起始密碼子和終止密碼子位置，遵循盡可能少的基因間隔和基因重疊且起始終止密碼完整的原則。

1.3 西氏熊蜂的生物信息學(xué)分析

采用在線網(wǎng)站CGView Server（http：//cgview.ca/）繪制西氏熊蜂線粒體基因組圈圖。利用Phylosuite v1.2.2提取并分別計(jì)算西氏熊蜂線粒體基因組中各編碼基因的堿基組成、AT偏倚（AT-skew）和GC偏倚（GC-skew），其中，堿基偏倚的計(jì)算公式為：AT偏倚= （A-T） / （A+T）、GC偏倚= （G-C） / （G+C）。利用Phylosuite v1.2.2和Rstudio 4.3.0分析并繪制蛋白編碼基因的相對(duì)同義密碼子使用度柱狀圖（Relative synonymous codon usage，RSCU）［20］。采用在線軟件MITOS和ARWEN（http：//130.235.244.92/ARWEN/）對(duì)tRNA基因進(jìn)行查找及二級(jí)結(jié)構(gòu)的預(yù)測(cè)［21］。對(duì)西氏熊蜂和其余11種熊蜂及果蠅（Drosophila yakuba）的線粒體基因組序列，利用在線網(wǎng)站ITOL（https：//itol.embl.de）繪制13條序列中37種基因的排列，分析不同基因的重排。

1.4 西氏熊蜂系統(tǒng)發(fā)育分析

為研究西氏熊蜂與其他熊蜂的系統(tǒng)發(fā)育關(guān)系，選取已報(bào)道的熊蜂屬13個(gè)亞屬的20條序列作為內(nèi)群，蜜蜂屬的3種昆蟲為外群（表1），使用MEGA 10軟件中的鄰接法和最大似然法構(gòu)建系統(tǒng)進(jìn)化樹［22］。

2 結(jié)果與分析

2.1 西氏熊蜂線粒體基因組的結(jié)構(gòu)特征分析

西氏熊蜂線粒體基因組長度為15 711 bp，包含37個(gè)基因，分別為13個(gè)蛋白編碼基因（protein-coding genes，PCGs）、22個(gè)tRNA基因和2個(gè)rRNA基因（圖1）。其中有9個(gè)PCGs和13個(gè)tRNA基因在主要編碼鏈（H鏈）上；4個(gè)PCGs、9個(gè)tRNA基因和2個(gè)rRNA基因在次要編碼鏈（L鏈）上。

西氏熊蜂的線粒體相鄰基因間存在長度不等的基因間隔和基因重疊，本研究共發(fā)現(xiàn)21個(gè)基因間隔和11個(gè)基因重疊，trnY和trnW間基因間隔最長（198 bp），nad4和nad4L間間隔最短（1 bp）；在相鄰基因的基因重疊中，trnL1和rrnL間基因重疊最長，為45 bp（表2）。

2.2 西氏熊蜂線粒體基因組的堿基組成分析

西氏熊蜂線粒體基因組全序列、蛋白編碼基因、tRNA基因和rRNA基因4個(gè)數(shù)據(jù)集的A+T含量分別為85.9%，83.9%，87.2%和86.7%，其中13個(gè)蛋白編碼基因的A+T含量分別在77%～90.9%之間，其中cox1最低，為77%，nad6最高，達(dá)到90.9%（表3）。對(duì)每個(gè)蛋白編碼基因的AT偏倚分析中發(fā)現(xiàn)，所有蛋白編碼基因都呈現(xiàn)T偏倚，但在GC偏倚分析中4個(gè)蛋白編碼基因（atp6，cox2，nad2，nad5）呈現(xiàn)C偏倚，cox3和cytb無堿基偏倚，其余7個(gè)蛋白編碼基因呈現(xiàn)G偏倚。對(duì)4個(gè)數(shù)據(jù)集的堿基偏倚分析發(fā)現(xiàn)，13個(gè)PCGs串聯(lián)數(shù)據(jù)集與線粒體基因組全序列的堿基偏倚結(jié)果一致，均對(duì)T堿基和C堿基具有明顯的堿基偏倚；22個(gè)tRNAs和2個(gè)rRNAs數(shù)據(jù)集的分析結(jié)果卻表明對(duì)A堿基、G堿基具有明顯的堿基偏倚。所以與其他數(shù)據(jù)集相比，13個(gè)PCGs的串聯(lián)數(shù)據(jù)集與線粒體基因組的組成特征更加一致。

2.3 西氏熊蜂蛋白編碼基因的密碼子使用分析

西氏熊蜂13個(gè)蛋白編碼基因的串聯(lián)序列總長度為11 022 bp，占全基因組的71.15%，共編碼了3 674個(gè)氨基酸。氨基酸的百分含量由高到低排列如下：Leu＞Ile＞Met＞Phe＞Ser＞Asn＞Tyr＞Lys＞Gly＞Thr＞Val＞Pro＞Trp＞Glu＞His＞Asp＞Ala＞Gln＞Arg＞Cys。其中氨基酸占比最高的是亮氨酸（Leu），為15.16%，半胱氨酸（Cys）占比最少，為0.73%。在相對(duì)同義密碼子使用度分析中發(fā)現(xiàn)，所有氨基酸對(duì)第3位堿基是A/U的密碼子具有明顯使用偏好，如精氨酸（Arg）偏好使用CGA、脯氨酸（Pro）偏好使用CCA、天冬氨酰（Asn）偏好使用AAU、纈氨酸（Val）偏好使用GUU等（圖2）。其中13個(gè)蛋白編碼基因的起始密碼子均為ATN，除nad4以不完整的終止密碼子T結(jié)尾，其余12個(gè)蛋白編碼基因均以完整的終止密碼子（TAA/TAG）結(jié)尾。

2.4 西氏熊蜂tRNA二級(jí)結(jié)構(gòu)預(yù)測(cè)分析

西氏熊蜂的線粒體基因組注釋得到22個(gè)tRNA基因，其中13個(gè)位于H鏈，9個(gè)位于L鏈，單個(gè)tRNA基因的序列長度在59～74 bp之間。tRNA基因的二級(jí)結(jié)構(gòu)預(yù)測(cè)結(jié)果如圖3所示，其中trnS1因缺少1個(gè)二氫尿嘧啶臂（DHU臂）和二氫尿嘧啶環(huán)（DHU環(huán)）不能形成典型的三葉草結(jié)構(gòu)，其余tRNA均能構(gòu)成典型的三葉草結(jié)構(gòu)，但trnF，trnH，trnN和trnV基因的二級(jí)結(jié)構(gòu)中缺少TφC環(huán)。22個(gè)tRNA中共發(fā)現(xiàn)12次G-U堿基錯(cuò)配，分別在trnA，trnC，trnD，trnH，trnI，trnK，trnP，trnL1，trnQ和trnV中各出現(xiàn)1次，trnF中出現(xiàn)2次。

2.5 線粒體基因組的基因重排分析

熊蜂與果蠅的線粒體基因排列相比具有明顯不同（圖4），其基因重排主要發(fā)生在以下6個(gè)區(qū)域：trnK移位的trnD-trnK基因簇；trnA移位的trnQ-trnA-trnM-trnI基因簇；trnS1移位的trnE-trnS1-trnF基因簇；trnT移位的trnP-trnT基因簇；trnI移位的trnA-trnM-trnI基因簇；trnW移位的trnC-trnY-trnW基因簇。

不同熊蜂亞屬與果蠅線粒體基因重排的分析結(jié)果如下：

（1）在trnR的移位分析發(fā)現(xiàn)，胸熊蜂亞屬的trnR移位到trnN的下游，形成trnN-trnR基因簇；巨熊蜂亞屬的trnR移位到trnQ的下游，形成trnQ-trnR基因簇，其他5個(gè)亞屬則為trnR-trnN基因簇。（2）trnT移位到nad6上游，形成trnP-trnT基因簇，這種基因重排在熊蜂亞屬、黑熊蜂亞屬、阿熊蜂亞屬、西熊蜂亞屬共4個(gè)亞屬中均有出現(xiàn)，其他亞屬未見此種基因重排。（3）在trnA的移位分析中共發(fā)現(xiàn)2種不同的移位現(xiàn)象，一種是在胸熊蜂亞屬、假熊蜂亞屬、黑熊蜂亞屬、西熊蜂亞屬的西伯熊蜂和阿熊蜂亞屬的克什米爾熊蜂中發(fā)現(xiàn)的trnA移位到trnQ下游，形成trnQ-trnA-trnM-trnI基因簇。另一種在熊蜂亞屬、阿熊蜂亞屬的短頭熊蜂和西熊蜂亞屬的亞西伯熊蜂中發(fā)現(xiàn)的trnA移位到trnM下游，形成trnQ-trnM-trnA-trnI基因簇。（4）在trnI的移位分析中發(fā)現(xiàn)，除假熊蜂亞屬外其他6個(gè)亞屬的trnI均移位到nad2上游，形成trnA-trnM-trnI基因簇，而假熊蜂亞屬的trnI卻移位到trnM上游，形成trnA-trnI-trnM基因簇。

本研究除阿熊蜂亞屬和西熊蜂亞屬內(nèi)具有trnA的基因重排，阿熊蜂亞屬內(nèi)的短頭熊蜂具有nad2的缺失，其他各亞屬內(nèi)均無線粒體基因重排。西氏熊蜂與同樣隸屬于黑熊蜂亞屬的火紅熊蜂的線粒體基因排列順序一致。通過上述的比較分析，發(fā)現(xiàn)熊蜂各亞屬間具有較為頻繁的基因重排，各亞屬內(nèi)的基因重排不頻繁。

2.6 線粒體基因組系統(tǒng)發(fā)育分析

本研究利用最大似然法和鄰接法構(gòu)建基于20種熊蜂線粒體基因組全序列的系統(tǒng)發(fā)育樹（圖5、圖6）。2種方法構(gòu)建的系統(tǒng)發(fā)育樹均高度支持西氏熊蜂與火紅熊蜂共同隸屬于黑熊蜂亞屬，且支持率為100%。2個(gè)系統(tǒng)發(fā)育樹的拓?fù)浣Y(jié)構(gòu)對(duì)各亞屬之間親緣關(guān)系的劃分結(jié)果不完全相同，但兩者均證明假熊蜂亞屬（Mendacibombus）和巴熊蜂亞屬（Confusibombus）是最先分化的一支。在2個(gè)系統(tǒng)發(fā)育樹上可以看出，遠(yuǎn)離根部的分支支持率高，對(duì)亞屬的分類具有100%的支持率；距離根部近的分支節(jié)點(diǎn)支持率低，普遍在80%以下。其中地熊蜂亞屬在ML發(fā)育樹中與巨熊蜂亞屬（Megabombus）的親緣關(guān)系更近，在NJ發(fā)育樹中與凱熊蜂亞屬（Kallobombus）的親緣關(guān)系更近，由此可見利用線粒體基因組序列對(duì)亞屬間的親緣關(guān)系劃分尚且不夠準(zhǔn)確。但2個(gè)系統(tǒng)發(fā)育樹的結(jié)果均可以看出短臉型的熊蜂亞屬（Bombus）、黑熊蜂亞屬（Melanobombus）、阿熊蜂亞屬（Alpigenobombus）、火熊蜂亞屬（Pyrobombus）、西熊蜂亞屬（Sibiricytbombus）和卡熊蜂亞屬（Cullumanobombus）之間的親緣關(guān)系更近，與長臉型的巨熊蜂亞屬（Megabombus）、胸熊蜂亞屬（Thoracytbombus）和擬熊蜂亞屬（Psithyrus）等亞屬親緣關(guān)系相對(duì)較遠(yuǎn)，這與前人根據(jù)頭部的形態(tài)學(xué)特征［23］和單基因分子鑒定的結(jié)果一致。

3 討論

西氏熊蜂首次完成線粒體基因組測(cè)序及功能注釋，其長度為15 711 bp。序列組成分析結(jié)果表明，西氏熊蜂線粒體基因組存在明顯的堿基偏倚，一方面，與前人所發(fā)現(xiàn)的昆蟲線粒體基因組具有高A+T含量的特點(diǎn)相一致，西氏熊蜂的線粒體基因組全序列、蛋白編碼基因、tRNAs基因、rRNAs基因的A+T含量均大于80%；另一方面，在相對(duì)同義密碼子的分析中也發(fā)現(xiàn)第3位堿基是A/U的密碼子具有明顯的高使用現(xiàn)象［10］。西氏熊蜂線粒體蛋白編碼基因的起始密碼子為典型的ATN，終止密碼子主要是TAA和TAG，僅nad4是以不完整的終止密碼子T結(jié)尾，這種使用不完整的終止密碼子在其他膜翅目昆蟲中也時(shí)有發(fā)生［24］。西氏熊蜂線粒體基因組的蛋白編碼序列全長11 022 bp，共編碼氨基酸3647個(gè)，其中頻率最高的依次為亮氨酸（Leu，15.16%）、異亮氨酸（Ile，13.15%）、蛋氨酸（Met，10.94%）和苯丙氨酸（Phe，10.83%），而頻率最低的依次為丙氨酸（Ala，1.42%）、谷氨酰胺（Gln，1.22%）、精氨酸（Arg，1.01%）和半胱氨酸（Cys，0.73%），這種豐度排序與其他膜翅目昆蟲中基本一致［25］。

對(duì)西氏熊蜂線粒體基因組的tRNA二級(jí)結(jié)構(gòu)預(yù)測(cè)發(fā)現(xiàn)，除了trnS1缺少DHU臂和DHU環(huán)，其他都呈現(xiàn)典型的三葉草結(jié)構(gòu)，其中trnS1缺臂的現(xiàn)象在后生動(dòng)物線粒體基因組中普遍存在［26］。本研究在22個(gè)tRNA中共發(fā)現(xiàn)12次堿基錯(cuò)配，這種錯(cuò)配在tRNA的折疊過程中是被允許的，且在其他后生動(dòng)物的線粒體中也普遍存在［27］。本研究所發(fā)現(xiàn)的堿基錯(cuò)配為G-U錯(cuò)配，有學(xué)者認(rèn)為這種錯(cuò)配可在RNA編輯中校正，并不影響tRNA基因的轉(zhuǎn)運(yùn)功能［28］。

果蠅線粒體基因組的基因排列順序一直被認(rèn)為是昆蟲線粒體基因組的原始排列方式［29］。早在上個(gè)世紀(jì)，研究人員就已經(jīng)發(fā)現(xiàn)，膜翅目昆蟲的線粒體基因組比其他目的昆蟲存在更為頻繁的基因重排，且細(xì)腰亞目的基因重排比廣腰亞目更多［30］。通過將西氏熊蜂線粒體基因組的基因排列與果蠅比較，結(jié)果發(fā)現(xiàn)明顯的基因重排，且均為tRNA的移位，蛋白編碼基因無基因重排。這與前人發(fā)現(xiàn)膜翅目線粒體基因組的基因重排主要是tRNA的重排結(jié)果一致［31］。而trnR在巨熊蜂亞屬內(nèi)的遠(yuǎn)距離移位尚不能確定這是否在亞屬內(nèi)普遍發(fā)生，在其他亞屬的研究中也并沒有此類發(fā)現(xiàn)。

本研究中，NJ和ML系統(tǒng)發(fā)育樹的結(jié)果均顯示，西氏熊蜂與火紅熊蜂聚為一支，且與形態(tài)學(xué)鑒定結(jié)果一致均隸屬于黑熊蜂亞屬。由于在2個(gè)系統(tǒng)發(fā)育樹中，各亞屬的親緣關(guān)系不完全一致，但兩棵樹對(duì)遠(yuǎn)離根部的分支聚類均具有高的支持率（100%），即利用線粒體基因組序列對(duì)當(dāng)前物種亞屬的劃分是具有高支持率的，符合Willims等［32］基于形態(tài)學(xué)和線粒體基因片段對(duì)各亞屬的分類及修訂結(jié)果。在2個(gè)系統(tǒng)發(fā)育樹中，距離根部近的分支節(jié)點(diǎn)普遍在80%以下，其結(jié)果對(duì)探討亞屬間親緣關(guān)系的準(zhǔn)確度不高，這或許是由于熊蜂物種的線粒體基因組序列不充足導(dǎo)致，所以對(duì)熊蜂屬內(nèi)更多物種的線粒體基因組進(jìn)行測(cè)序，補(bǔ)充其線粒體基因組數(shù)據(jù)，對(duì)后續(xù)探索各亞屬間的親緣關(guān)系具有重要意義。

4 結(jié)論

本研究獲得了西氏熊蜂的線粒體基因組序列，明確了其線粒體基因組組成和結(jié)構(gòu)特征，即線粒體基因組包含37個(gè)基因的基本組成，對(duì)A/T明顯的堿基偏倚等結(jié)果。西氏熊蜂線粒體基因組與已報(bào)道的膜翅目中大部分昆蟲的線粒體基因組相似，且與蜜蜂科昆蟲線粒體基因組的一般特征一致。西氏熊蜂的研究結(jié)果進(jìn)一步豐富了熊蜂屬的線粒體基因組數(shù)據(jù)，對(duì)西氏熊蜂的分子鑒定和種群遺傳學(xué)研究以及熊蜂屬系統(tǒng)發(fā)育關(guān)系的深入探討等具有重要意義。

參考文獻(xiàn)

［1］ 安建東，黃家興，WILLIAMS P H，等. 河北地區(qū)熊蜂物種多樣性與蜂群繁育特性［J］. 應(yīng)用生態(tài)學(xué)報(bào)，2010，21（6）：1542-1550

［2］ 黃家興，安建東. 中國熊蜂多樣性、人工利用與保護(hù)策略［J］. 生物多樣性，2018，26（5）：486-497

［3］ 蒙艷華，徐環(huán)李，陳軒，等. 塔落巖黃芪主要傳粉蜂的傳粉效率研究［J］. 生物多樣性，2007，15（6）：633-638

［4］ 何亞平，劉建全. 青藏高原高山植物麻花艽的傳粉生態(tài)學(xué)研究［J］. 生態(tài)學(xué)報(bào)，2004，24（2）：215-220

［5］ 梁程博，王久利，孫國，等. 三江源地區(qū)熊蜂物種多樣性研究［J］. 草地學(xué)報(bào)，2022，30（8）：2126-2134

［6］ 張廣碩，丁桂玲，安建東. 熊蜂主要病蟲害的發(fā)生和危害現(xiàn)狀［J］. 應(yīng)用昆蟲學(xué)報(bào)，2021，58（2）：277-288

［7］ 胥勛露，常雪梅，曾桂，等. 甘南草原和若爾蓋濕地的熊蜂物種資源調(diào)查［J］. 西華師范大學(xué)學(xué)報(bào)（自然科學(xué)版），2024，45（3）：239-245

［8］ 劉青青，赫苗花，李飛，等. 不同坡向下高寒礦區(qū)混播草地的恢復(fù)效果評(píng)估［J］. 草地學(xué)報(bào)，2023，31（9）：2834-2842

［9］ 梁鋮，張學(xué)文，黃家興，等. 云南熊蜂地理區(qū)劃及物種多樣性分析［J］. 應(yīng)用昆蟲學(xué)報(bào)，2018，55（6）：1045-1053

［10］朱國淵，張永科，孔祥東，等. 二斑素瓢蟲線粒體基因組全序列測(cè)定和分析［J］. 南方農(nóng)業(yè)學(xué)報(bào)，2023，54（5）：1293-1302

［11］ZHANG J Y，ZHANG L P，YU D N，et al. Complete mitochondrial genomes of Nanorana taihangnica and N.yunnanensis （Anura：Dicroglossidae） with novel gene arrangements and phylogenetic relationship of Dicroglossidae［J］. BMC Evolutionary Biology，2018，18（1）：26-39

［12］黃家興. 華北地區(qū)熊蜂屬（Hymenoptera：Apidae）系統(tǒng)發(fā)育的初步研究［D］. 北京：中國農(nóng)業(yè)科學(xué)院，2006：2

［13］BERTSCH A. Barcoding cryptic bumblebee taxa：B. lucorum，B. crytarum and B. magnus，a case study［J］. Beitrage Zur Entomol，2009，59（2）：287-310

［14］劉蘋，黃家興，吳杰，等. 中國明亮熊蜂復(fù)合種的分子鑒定及分布特性［J］. 昆蟲學(xué)報(bào)，2014，57（2）：235-243

［15］PIROUNAKIS K，KOULIANOS S，SCHMID-HEMPEL P. Genetic variation among European populations of Bombus pascuorum （Hymenoptera：Apidae） from mitochondrial DNA sequence data［J］. European Journal of Entomology，2013，95（1）：27-33

［16］WILLIAMS P H，BEREZIN M V，CANNINGS S G，et al. The arctic and alpine bumblebees of the subgenus Alpinobombus revised from integrative assessment of species’ gene coalescents and morphology （Hymenoptera，Apidae，Bombus）［J］.Zootaxa，2019，4625（1）：1-68

［17］CEJAS D，DE LA RU′A P，ORNOSA C，et al. Spatial and temporal patterns of genetic diversity in Bombus terrestris populations of the Iberian Peninsula and their conservation implications［J］. Scientific Reports，2021，11（1）：22471

［18］SU M，DU L J，ALI M Y，et al. Differences in morphology，mitochondrial genomes，and reproductive compatibility between two clades of parasitic wasps Aphelinus mali （Hymenoptera：Aphelindae） in China［J］. Plos One，2023，18（2）：e0279663

［19］郭冉昊，趙淑文，米福貴，等. 直立型半野生大豆葉綠體基因組分析［J］. 草地學(xué)報(bào)，2023，31（11）：3334-3342

［20］郝新艷，趙淑文，劉嘉偉，等. 雜花苜蓿葉綠體基因組特征及系統(tǒng)發(fā)育分析［J］. 草地學(xué)報(bào)，2023，31（6）：1665-1672

［21］LIN Y J，CAI L N，ZHAO Y Y，et al. Novel mitochondrial gene rearrangement and intergenic regions exist in the mitochondrial genomes from four newly established families of praying mantises （Insecta：Mantodea）［J］. Insects，2022，13（7）：564-564

［22］李瑤瑤，劉云國，劉凌霄，等. 高原鰍屬魚類線粒體全基因組序列結(jié)構(gòu)特征及其系統(tǒng)發(fā)育信息［J］. 煙臺(tái)大學(xué)學(xué)報(bào)（自然科學(xué)與工程版），2021，34（1）：93-99

［23］CAMERON A S，HINES M H，WILLIAMS H P. A comprehensive phylogeny of the bumble bees （Bombus）［J］. Biological Journal of the Linnean Society，2007，91（1）：161-188

［24］趙婉清，李巧，李瑩，等. 獵蝽科昆蟲線粒體基因組的比較研究［J］. 江蘇農(nóng)業(yè)科學(xué)，2020，50（12）：34-41

［25］彭艷，陳斌，李廷景. 黃側(cè)異腹胡蜂線粒體基因組全序列測(cè)定和分析［J］. 昆蟲學(xué)報(bào)，2017，60（4）：464-474

［26］WOLSTENHOLME D R. Animal mitochondrial DNA：structure and evolution［J］. International Review of Cytology，1992（141）：173-216

［27］陳建琴，張振華，邢雨輝，等. 光滑異裝蟹（Heteropanope glabra）線粒體基因組測(cè)定及其在系統(tǒng)發(fā)生研究中的意義［J］. 南京師大學(xué)報(bào)（自然科學(xué)版），2018，41（4）：114-120

［28］韋春夢(mèng)，陳亞珍，劉靜，等. 云南瘤蚱線粒體基因組特征與蚱總科系統(tǒng)發(fā)育分析［J］. 南方農(nóng)業(yè)學(xué)報(bào)，2023，54（5）：1303-1316

［29］CLARY D O，WOLSTENHOOLME D R. The mitochondrial DNA molecular of Drosophila yakuba：nucleotide sequence，gene organization，and genetic code［J］. Journal of Molecular Evolution，1985，22（3）：252-271

［30］JERMIIN L S，CROZIER R H. The cytochrome b region in the mitochondrial DNA of the ant Tetraponera rufoniger：sequence divergence in Hymenoptera may be associated with nucleotide content［J］. Journal of Molecular Evolution，1994，38（3）：282-294

［31］KORKMAZ E M，AYDEMIR H B，TEMEL B，et al. Mitogenome evolution in Cephini （Hymenoptera：Cephidae）：Evidence for parallel adaptive evolution［J］. Biochemical Systematics and Ecology，2017（71）：137-146

［32］WILLIAMS P H，CAMERON S A，HINES H M，et al. A simplified subgeneric classification of the bumblebee （genus Bombus）［J］. Apidologie，2008，39（1）：46-74

（責(zé)任編輯 閔芝智）