周文才 李應(yīng) 姚丹 萬煒 詹洪平 黎華君 賀興江 冉曜琦 于瀛龍 韋小平

DOI：10.3969/j.issn.2095-1191.2023.06.023

摘要：【目的】明確貴州省東方蜜蜂群體的遺傳結(jié)構(gòu)和遺傳多樣性，掌握其遺傳現(xiàn)狀，為貴州省東方蜜蜂的遺傳資源評(píng)價(jià)及保護(hù)利用提供參考依據(jù)?！痉椒ā炕诰€粒體tRNAleu～COII和COI片段對(duì)貴州省內(nèi)19個(gè)采樣點(diǎn)的1282群東方蜜蜂進(jìn)行遺傳分析，計(jì)算單倍型多樣性（Hd）、平均核苷酸差異數(shù)（K）和核苷酸多樣性（π）等遺傳多樣性指數(shù)，采用中性檢驗(yàn)（Tajima?s D和Fu?s Fs）評(píng)估群體擴(kuò)張情況，通過計(jì)算遺傳分化系數(shù)（Fst）、構(gòu)建系統(tǒng)發(fā)育進(jìn)化樹和分子方差分析評(píng)估其遺傳分化，并開展環(huán)境因子相關(guān)分析以明確遺傳分化的作用因子?！窘Y(jié)果】貴州東方蜜蜂在線粒體tRNAleu～COII片段上的Hd介于0.5421～0.9167，π介于0.00174～0.00527，K介于0.611～1.854，其中臺(tái)江、雷山、麻江、羅甸、紫云、晴隆和威寧等7個(gè)采樣點(diǎn)發(fā)生過群體擴(kuò)張；在線粒體COI片段上的Hd介于0.5710～0.9109，π介于0.00109～0.00261，K介于0.874～2.084，其中松桃、江口、印江、臺(tái)江、石阡、雷山、羅甸、紫云、黔西、晴隆和威寧等11個(gè)采樣點(diǎn)發(fā)生過群體擴(kuò)張。貴州東、西部東方蜜蜂群體間存在中度遺傳分化（在tRNAleu～COII片段上的Fst=0.04±0.02，在COI片段上的Fst=0.07±0.04），貴州西部東方蜜蜂群體內(nèi)也存在中度遺傳分化（在tRNAleu～COII片段上的Fst=0.04±0.02，在COI片段上的Fst=0.07±0.03）。相關(guān)分析結(jié)果顯示，貴州西部東方蜜蜂組群的遺傳分化程度與海拔呈弱相關(guān)，但與溫度、濕度等環(huán)境因子呈顯著正相關(guān)?！窘Y(jié)論】貴州省東方蜜蜂在線粒體tRNAleu～COII和COI片段上的遺傳多樣性處于較高水平，主要?dú)w功于蜂群數(shù)量大和飼養(yǎng)方式較原始；貴州東方蜜蜂整體上呈現(xiàn)中度遺傳分化，且西部東方蜜蜂組群內(nèi)部的遺傳分化程度較東部東方蜜蜂組群更高，主要與溫度和濕度等環(huán)境因子密切相關(guān)，是其為適應(yīng)所處環(huán)境而進(jìn)化產(chǎn)生，而非海拔因素造成。

關(guān)鍵詞：東方蜜蜂；tRNAleu～COII片段；COI片段；遺傳多樣性；遺傳分化；環(huán)境因子

中圖分類號(hào)：S891.8? ? ? ? ? ? ? ? ? ? ? ? ? 文獻(xiàn)標(biāo)志碼：A 文章編號(hào)：2095-1191（2023）06-1819-10

Population genetic analysis of Apis cerana in Guizhou based on mitochondrial tRNAleu-COII and COI

ZHOU Wen-cai， LI Ying， YAO Dan， WAN Wei， ZHAN Hong-ping， LI Hua-jun，

HE Xing-jiang， RAN Yao-qi， YU Ying-long*， WEI Xiao-ping*

（Institute of Integrated Agriculture Development，Guizhou Academy of Agricultural Sciences，

Guiyang，Guizhou? 550006，China）

Abstract：【Objective】The purpose of the study was to clarify the population genetic structure and genetic diversity of Apis cerana in Guizhou， to understand its genetic status， and to provide a reference basis for the evaluation， conservation and utilization of Apis cerana genetic resources in Guizhou. 【Method】Genetic analysis was conducted on 1282 populations of A. cerana from 19 sampling sites in Guizhou based on mitochondrial tRNAleu-COII and COI fragments. Genetic diversity indexes such as haplotype diversity （Hd）， the average nucleotide difference （K） and nucleotide diversity （π） were calculated. Tajima’s D and Fu’s Fs neutral tests were carried out to evaluate population expansion. Genetic differentiation was assessed by calculating the genetic differentiation coefficient （Fst）， constructing phylogenetic evolutionary trees and molecular variance analysis， and environmental factor correlation analysis was carried out to clarify the role of factors in genetic differentiation. 【Result】The Hd of A. cerana in Guizhou on the mitochondrial tRNAleu-COII fragment ranged from 0.5421 to 0.9167， π ranged from 0.00174 to 0.00527， and K ranged from 0.611 to 1.854， with colony expansions occurring in seven sampling sites， including Taijiang， Leishan， Majiang， Luodian， Ziyun， Qinglong， and Weining. On the mitochondrial COI fragment， Hd ranged from 0.5710 to 0.9109， π ranged from 0.00109 to 0.00261， and K ranged from 0.874 to 2.084， with population expansion occurring at 11 sampling sites， including Songtao， Jiangkou， Yinjiang， Taijiang， Shiqian， Leishan， Luodian， Ziyun， Qianxi， Qinglong， and Weining. There existed moderate genetic differentiation （on tRNAleu-COII fragment， Fst=0.04±0.02， on COI fragment Fst=0.07±0.04） between A. cerana populations in the eastern and western of Guizhou， and moderate genetic differentiation （on tRNAleu-COII fragment， Fst=0.04±0.02， on COI fragment， Fst=0.07±0.03） within A. cerana populations in the western part of Guizhou. The results of correlation analysis showed that genetic differentiation of A. cerana colonies in western Guizhou was weakly correlated with altitude， but significantly positively correlated with environmental factors such as temperature and humidity. 【Conclusion】The genetic diversity of A. cerana in Guizhou is at a high level in mitochondrial tRNAleu-COII and COI fragments， which is mainly attributed to the huge population size and primitive rearing methods. Guizhou A. cerana as a whole show moderate genetic differentiation， and the degree of genetic differentiation within the western A. cerana populations is higher than that of the eastern A. cerana， which is mainly closely related to the environmental factors， such as temperature and humidity， and arises from their evolution to adapt to the environment they live in rather than due to altitudinal factors.

Key words： Apis cerana； tRNAleu-COII fragment； COI fragment； genetic diversity； genetic differentiation； environmental factor

Foundation items： Guizhou Science and Technology Plan Project （QKHJC〔2019〕1453）；Construction Project of China Agriculture Research System （CARS-44-SYZ-22）；Guizhou Academy of Agricultural Sciences Project （QNKYQNKJJJ〔2020〕08，QNKYZZZY〔2022〕16）

0 引言

【研究意義】東方蜜蜂（Apis cerana）是我國本土蜂種，廣泛分布在除內(nèi)蒙古和新疆外的我國其他地區(qū)，具有善于利用零星蜜粉源、耐低溫、抗螨等優(yōu)良特性，在維持生物多樣性、維護(hù)生態(tài)系統(tǒng)平衡及糧食安全等方面發(fā)揮著重要作用（楊冠煌，2009）。貴州地處云貴高原東段，屬亞熱帶濕潤季風(fēng)氣候，自東向西海拔逐漸升高，省內(nèi)有東方蜜蜂51萬群，占全國總量的10%，是我國東方蜜蜂的主產(chǎn)區(qū)之一（韋小平等，2020）。開展貴州省本土東方蜜蜂群體遺傳分析，既是掌握蜂種資源的基礎(chǔ)性工作，也是蜂種資源保護(hù)與利用的前提?！厩叭搜芯窟M(jìn)展】《中國畜禽遺傳資源志 蜜蜂志》將貴州省東方蜜蜂劃分為云貴高原中蜂和華中中蜂，其中，云貴高原中蜂分布于貴州省西部，華中中蜂分布在貴州省東部（國家畜禽遺傳資源委員會(huì)，2011）。至今，國內(nèi)已有諸多學(xué)者基于形態(tài)、微衛(wèi)星遺傳分子標(biāo)記等開展了貴州東方蜜蜂群體遺傳分析，但得出的相關(guān)研究結(jié)果存在明顯差異。徐祖蔭（1986，2009）利用7個(gè)形態(tài)標(biāo)記均值將貴州東方蜜蜂劃分為東部的山地生態(tài)型和西部的云貴高原生態(tài)型，以織金—冊亨、赫章—盤縣為界限；楊冠煌（2001）基于12個(gè)形態(tài)標(biāo)記均值將東方蜜蜂劃分為華中生態(tài)型和云貴高原生態(tài)型；駱群等（2015）利用32個(gè)形態(tài)標(biāo)記進(jìn)行多元化統(tǒng)計(jì)分析，但未發(fā)現(xiàn)貴州省東方蜜蜂發(fā)生種群分化；于瀛龍等（2017）利用31個(gè)微衛(wèi)星位點(diǎn)對(duì)貴州省東方蜜蜂進(jìn)行分析，也未發(fā)現(xiàn)貴州省東方蜜蜂發(fā)生遺傳分化。此外，有學(xué)者通過線粒體tRNAleu～COII、COI及COII等片段分別開展了我國西藏、廣西、重慶、浙江、云南等地的東方蜜蜂遺傳多樣性研究（姬聰慧等，2017；劉意秋等，2018；曹聯(lián)飛等，2021；周姝婧等，2021），證實(shí)山脈、高原等自然隔離條件，以及人為干擾造成的基因交流等是塑造東方蜜蜂群體遺傳結(jié)構(gòu)的主要因素，東方蜜蜂不同地理群體間遺傳分化明顯，且蘊(yùn)藏著豐富的遺傳多樣性，其研究結(jié)論對(duì)東方蜜蜂遺傳資源保護(hù)與合理開發(fā)利用具有基礎(chǔ)性指導(dǎo)意義，但針對(duì)于貴州省本土東方蜜蜂群體的研究僅趙金玉等（2022）對(duì)79只貴州省東方蜜蜂進(jìn)行線粒體COI～COII片段分析?！颈狙芯壳腥朦c(diǎn)】線粒體tRNAleu～COII和COI片段是目前廣泛應(yīng)用于蜜蜂遺傳多樣性研究的有效手段，但鮮見應(yīng)用于貴州省東方蜜蜂的系統(tǒng)研究?！緮M解決的關(guān)鍵問題】基于2個(gè)線粒體片段tRNAleu～COII和COI對(duì)貴州全省范圍內(nèi)的東方蜜蜂開展大樣本量群體遺傳分析，以期明確貴州省東方蜜蜂群體的遺傳結(jié)構(gòu)和遺傳多樣性，掌握其遺傳現(xiàn)狀，為貴州省東方蜜蜂的遺傳資源評(píng)價(jià)及保護(hù)利用提供參考依據(jù)。

1 材料與方法

1. 1 樣本采集

東方蜜蜂樣本采自貴州省19個(gè)主要縣（市），共1282群（表1），樣本量介于10～176群/采樣點(diǎn)。東方蜜蜂樣本均采自健康蜂群，為減少引種、育王等對(duì)代表性的影響，蜂群以半野生的飼養(yǎng)蜂群為主。采集后的樣本以無水乙醇處理后-20 ℃保存?zhèn)溆谩?/p>

1. 2 試驗(yàn)方法

每個(gè)蜂群取1只工蜂，利用UNIQ-10柱式試劑盒、SanTaq PCR Mix擴(kuò)增試劑盒［生工生物工程（上海）股份有限公司］從工蜂胸部提取基因組。為提高結(jié)果可靠性，本研究擴(kuò)增2個(gè)經(jīng)典的線粒體片段tRNAleu～COII（1～352 bp）和COI（1～799 bp），擴(kuò)增條件參照Garnery等（1992）、Zhao等（2014）的研究方法。PCR擴(kuò)增產(chǎn)物送至生工生物工程（上海）股份有限公司分別進(jìn)行單向、雙向測序，測序結(jié)果經(jīng)人工檢查，若存在嚴(yán)重雙峰現(xiàn)象則需重新測序。

1. 3 統(tǒng)計(jì)分析

利用ClustalX進(jìn)行測序序列比對(duì)，并將序列數(shù)據(jù)上傳至NCBI以確定其單倍型（Larkin et al.，2007）；通過DNAsp 5.0計(jì)算單倍型多樣性（Hd）、平均核苷酸差異數(shù)（K）和核苷酸多樣性（π）（Librado and Rozas，2009），以及Tajima?s D（Tajima，1989）和Fu?s Fs（Fu，1997）；使用Arlequin 3.11計(jì)算不同采樣點(diǎn)間的遺傳分化系數(shù)（Fst），并進(jìn)行分子方差分析（Excoffier and Lischer，2010；張桂寧等，2022）；利用MEGA 7.0的鄰接法（Neighbor-joining，NJ）基于Kimura雙參數(shù)（Kimura 2-parameter，K2P）構(gòu)建系統(tǒng)發(fā)育進(jìn)化樹（Kumar et al.，2016），外群選用Apis mellifera（GenBank登錄號(hào)KX908209）的對(duì)應(yīng)片段。采用GenAlex 6.5進(jìn)行遺傳分化與海拔等環(huán)境因子的相關(guān)分析（Peakall and Smouse，2012），其中，海拔數(shù)據(jù)來源于實(shí)地測量，其余主要環(huán)境因子數(shù)據(jù)來源于中國氣象信息網(wǎng)站，包括：月平均氣壓、月平均氣溫、月平均最高氣溫、月平均最低氣溫、每月最高氣溫≥30 ℃日數(shù)、每月最高氣溫≥35 ℃日數(shù)、每月最低氣溫≤2 ℃日數(shù)、每月最低氣溫≤0 ℃日數(shù)、月平均相對(duì)濕度及月平均降水量。

2 結(jié)果與分析

2. 1 貴州東方蜜蜂遺傳多樣性與歷史動(dòng)態(tài)

在貴州東方蜜蜂線粒體tRNAleu～COII片段上共發(fā)現(xiàn)98種單倍型（圖1），其中28種為新發(fā)現(xiàn)單倍型，對(duì)應(yīng)的GenBank登錄號(hào)分別為MT607430～MT607438、MT607519～MT607520、MZ970560～MZ970561、OL333585～ OL333588、OL333590～OL333598及OL989889～OL989890。各采樣點(diǎn)的單倍型種類介于3～27種，Hd介于0.5421～0.9167，π介于0.00174～0.00527，K介于0.611～1.854，其中臺(tái)江、雷山、麻江、羅甸、紫云、晴隆和威寧等7 個(gè)采樣點(diǎn)發(fā)生過群體擴(kuò)張，中性檢驗(yàn)（Tajima?s D，F(xiàn)u?s Fs）呈顯著性負(fù)值（表2）。

在貴州東方蜜蜂線粒體COI片段上共發(fā)現(xiàn)143 種單倍型（圖2），其中122種為新發(fā)現(xiàn)單倍型，對(duì)應(yīng)的GenBank登錄號(hào)分別為MT607439～MT607518、OK001729～OK001739及OL308023～OL308053。各采樣點(diǎn)的單倍型種類介于5～36種，Hd介于0.5710～0.9109，π介于0.00109～0.00261，K介于0.874～2.084，其中松桃、江口、印江、臺(tái)江、石阡、雷山、羅甸、紫云、黔西、晴隆和威寧等11個(gè)采樣點(diǎn)發(fā)生過群體擴(kuò)張，中性檢驗(yàn)（Tajima?s D，F(xiàn)u?s Fs）呈顯著性負(fù)值（表3）。

2. 2 貴州東方蜜蜂遺傳分化情況

貴州東方蜜蜂整體上呈現(xiàn)出中度遺傳分化（表4）。在線粒體tRNAleu～COII片段上的Fst介于 -0.01～0.15，平均為0.05±0.03；在線粒體COI片段上的Fst介于-0.01～0.18，平均為0.06±0.04。貴州東部（松桃、江口、沿河、印江、臺(tái)江、石阡、雷山等7個(gè)采樣點(diǎn)）東方蜜蜂群體與貴州西部（紫云、黔西、大方、晴隆、赫章、盤州、威寧等7個(gè)采樣點(diǎn)）東方蜜蜂群體間存在中度遺傳分化，在線粒體tRNAleu～COII片段上的Fst為0～0.09，平均為0.04±0.02；在線粒體COI片段上的Fst為0～0.18，平均為0.07±0.04。此外，貴州東部的松桃、江口、沿河、印江、臺(tái)江、石阡、雷山等7個(gè)采樣點(diǎn)東方蜜蜂在線粒體tRNAleu～COII片段上的Fst為0～0.09，平均為0.03±0.02；在線粒體COI片段上的Fst為0～0.13，平均為0.03±0.04。貴州西部的紫云、黔西、大方、晴隆、赫章、盤州、威寧等7 個(gè)采樣點(diǎn)東方蜜蜂在線粒體tRNAleu～COII片段上的Fst為0.01～0.07，平均為0.04±0.02；在線粒體COI片段上的Fst為0.03～0.15，平均為0.07±0.03。其他采樣點(diǎn)東方蜜蜂在線粒體tRNAleu～COII片段上的Fst介于0.05～0.15，平均為0.10±0.04；在線粒體COI片段上的Fst介于-0.01～0.09，平均為0.05±0.03。

分子方差分析時(shí)將貴州東方蜜蜂視為同一組群，基于線粒體tRNAleu～COII片段的分析結(jié)果顯示，群體間遺傳變異占總變異的4.76%，群體內(nèi)遺傳變異占總變異的95.24%（表5）；基于線粒體COI片段的分析結(jié)果顯示，群體間遺傳變異占總變異的5.18%，群體內(nèi)遺傳變異占總變異的94.82%（表6）?？梢姡F州東方蜜蜂不同群體間存在遺傳分化。

按照地理區(qū)域，將貴州東方蜜蜂設(shè)為貴州東部（松桃、江口、沿河、印江、臺(tái)江、石阡、雷山等7個(gè)采樣點(diǎn)）和貴州西部（紫云、黔西、大方、晴隆、赫章、盤州、威寧等7個(gè)采樣點(diǎn)）2個(gè)組群進(jìn)行分析，基于線粒體tRNAleu～COII片段的分析結(jié)果顯示，組群間遺傳變異占總變異的0.25%，組群內(nèi)群體間遺傳變異占總變異的3.84%，群體內(nèi)遺傳變異占總變異的95.91%（表7）；基于線粒體COI片段的分析結(jié)果顯示，組群間遺傳變異占總變異的1.45%，組群內(nèi)群體間遺傳變異占總變異的4.51%，群體內(nèi)遺傳變異占總變異的94.04%（表8）?？梢?，貴州東方蜜蜂遺傳分化還存在于貴州東西部之間。

2. 3 貴州東方蜜蜂遺傳分化與環(huán)境因子的相關(guān)性

對(duì)存在中度遺傳分化的貴州西部7個(gè)采樣點(diǎn)東方蜜蜂進(jìn)行Mantel檢驗(yàn)分析，結(jié)果顯示：線粒體COI片段的遺傳分化程度與月平均氣壓（R2=0.1390，P=0.04）、月平均最高氣溫（R2=0.2515，P=0.01）、每月最低氣溫≤0 ℃日數(shù)（R2=0.1413，P=0.03）及月平均相對(duì)濕度（R2=0.1146，P=0.04）呈顯著正相關(guān)（圖3）。線粒體tRNAleu～COII片段的遺傳分化程度與月平均降水量（R2=0.1622，P=0.03）呈顯著正相關(guān)（圖4）。線粒體tRNAleu～COII和COI片段的遺傳分化程度與海拔均呈弱相關(guān)性（tRNAleu～COII-R2=0.0191，P=0.10；COI-R2=0.0819，P=0.03）。

2. 4 貴州東方蜜蜂線粒體tRNAleu～COII和COI單倍型的系統(tǒng)發(fā)育進(jìn)化分析

基于線粒體tRNAleu～COII和COI片段各單倍型構(gòu)建的系統(tǒng)發(fā)育進(jìn)化樹（圖5和圖6）均顯示，貴州東方蜜蜂各單倍型聚類緊密，單倍型種類間的親緣關(guān)系較近，從單倍型種類角度上并未體現(xiàn)出明顯的地理發(fā)生關(guān)系。

3 討論

我國東方蜜蜂遺傳多樣性分析結(jié)果顯示，東方蜜蜂線粒體tRNAleu～COII的Hd介于0.171～0.911（周姝婧等，2012，2021），線粒體COI片段的Hd介于0.067～0.871（Zhao et al.，2014；劉意秋等，2018）。對(duì)比發(fā)現(xiàn)，貴州東方蜜蜂在線粒體tRNAleu～COII和COI片段上表現(xiàn)出較高的遺傳多樣性水平，究其原因可能是：（1）貴州省森林資源豐富，東方蜜蜂蜂群數(shù)量龐大。貴州省現(xiàn)存的東方蜜蜂蜂群數(shù)量約占全國的10%（韋小平等，2020），龐大的蜂群數(shù)量為其遺傳多樣性的形成奠定了基礎(chǔ)。2020年貴州省森林覆蓋率為61.51%，野生蜜源植物有243 種（韋小平等，2020），為蜂群的生存繁衍提供了物質(zhì)基礎(chǔ)。本研究的中性檢驗(yàn)（Tajima?s D，F(xiàn)u?s Fs）結(jié)果也顯示貴州東方蜜蜂蜂群體發(fā)生過擴(kuò)張，說明自然環(huán)境適宜蜜蜂群體數(shù)量增長。（2）貴州東方蜜蜂養(yǎng)殖仍保留有土法飼養(yǎng)，蜂群遺傳結(jié)構(gòu)受人為干擾的影響較小，尤其在自然保護(hù)區(qū)（江口和石阡）、西部山區(qū)（威寧）、南部石漠化地區(qū)（羅甸）等區(qū)域，土法飼養(yǎng)蜂群比例更高。土法飼養(yǎng)方式下蜂群操作不便，蜂農(nóng)全年開箱操作保持在1～2次，雖然養(yǎng)蜂生產(chǎn)效率低，但減少了對(duì)蜂群的干擾，有效避免了蜜蜂種群遺傳結(jié)構(gòu)受人為干擾，在自然保護(hù)區(qū)更是嚴(yán)格限制外來蜂群的進(jìn)入，使得東方蜜蜂遺傳多樣性極少受到人為因素的影響。

貴州東方蜜蜂整體上呈現(xiàn)中度遺傳分化，且西部東方蜜蜂組群內(nèi)部的遺傳分化程度較東部東方蜜蜂組群更高。相關(guān)分析發(fā)現(xiàn)，貴州西部東方蜜蜂組群的遺傳分化程度與海拔呈弱相關(guān)，但與溫度、濕度等環(huán)境因子呈顯著正相關(guān)，推測貴州西部東方蜜蜂的遺傳分化是其為適應(yīng)所處環(huán)境而進(jìn)化產(chǎn)生，而非海拔因素造成。溫度、濕度等環(huán)境因子直接影響蜜蜂出巢采集及婚飛等巢外活動(dòng)，同時(shí)影響蜜源植物的分布。前人研究發(fā)現(xiàn)，濕度明顯影響蜂王的正常交配，隨著降雨量的增加，與蜂王交配的雄蜂數(shù)量呈下降趨勢（Neves et al.，2011；Polatto et al.，2014）；溫度與西方蜜蜂不同亞種的分布范圍密切相關(guān)，西方蜜蜂中卡尼鄂拉亞種和馬其頓亞種的分布范圍以年均溫9 ℃為分界線（Coroian et al.，2014）。本研究中，致使貴州西部東方蜜蜂呈中度遺傳分化的主要環(huán)境因子可能是溫度和濕度。

本研究結(jié)果顯示，貴州東西部東方蜜蜂組群的遺傳分化以石阡—雷山、黔西—紫云為界限，在前人研究中貴州東方蜜蜂的遺傳分化則以織金—冊亨、赫章—盤縣為界限（徐祖蔭，1986；國家畜禽遺傳資源委員會(huì)，2011）。究其原因可能是：（1）貴州呈自東向西逐漸升高的地勢，無明顯的地理隔離因素，蜜蜂遺傳分化現(xiàn)象與所處棲息地的小環(huán)境相關(guān)，很難清晰界定遺傳分化的界限。（2）東方蜜蜂種群分布處于動(dòng)態(tài)變化。蜜蜂種群的大小及分布等與溫度和濕度等自然環(huán)境、蜜粉源等食物來源、人為干擾等因素密切相關(guān)，因此種群的分布界限也處于動(dòng)態(tài)變化過程中。

4 結(jié)論

貴州省東方蜜蜂在線粒體tRNAleu～COII和COI片段上的遺傳多樣性處于較高水平，主要?dú)w功于蜂群數(shù)量大和飼養(yǎng)方式較原始；貴州東方蜜蜂整體上呈現(xiàn)中度遺傳分化，且西部東方蜜蜂組群內(nèi)部的遺傳分化程度較東部東方蜜蜂組群更高，主要與溫度和濕度等環(huán)境因子密切相關(guān)，是其為適應(yīng)所處環(huán)境而進(jìn)化產(chǎn)生，而非海拔因素造成。

參考文獻(xiàn)：

曹聯(lián)飛，蘇曉玲，陳道印，趙東緒，華啟云. 2021. 浙江省金華市中華蜜蜂線粒體DNA遺傳多樣性研究［J］. 環(huán)境昆蟲學(xué)報(bào)，43（4）：986-991. ［Cao L F，Su X L，Chen D Y，Zhao D X，Hua Q Y. 2021. Genetic diversity of Apis cerana cerana based on mitochondrial DNA in Jinhua，Zhejiang，China［J］. Journal of Environmental Entomology，43（4）：986-991.］ doi：10.3969/j.issn.1674-0858.2021. 04.18.

國家畜禽遺傳資源委員會(huì). 2011. 中國畜禽遺傳資源志 蜜蜂志［M］. 北京：中國農(nóng)業(yè)出版社. ［China National Commission of Animal Genetic Resources. 2011. Animal genetic resources in China：Bees［M］. Beijing：China Agriculture Press.］

姬聰慧，戴榮國，羅文華，任勤，劉佳霖，王瑞生，曹聯(lián)飛. 2017. 重慶大巴山地區(qū)中蜂遺傳多樣性分析［J］. 黑龍江畜牧獸醫(yī)，（9）：232-233. ［Ji C H，Dai R G，Luo W H，Ren Q，Liu J L，Wang R S，Cao L F. 2017. Genetic diversity of Apis cerana cerana fabricius in the Daba mountain of Chongqing［J］. Heilongjiang Animal Science and Veterinary Medicine，（9）：232-233.］ doi：10.13881/j.cnki.hljxmsy.2017.0820.

劉意秋，周丹銀，龔雪陽，馮毅楠，和紹禹，趙文正. 2018. 線粒體DNA分析青藏高原東方蜜蜂（Apis cerana）種群遺傳結(jié)構(gòu)［J］. 山西農(nóng)業(yè)大學(xué)學(xué)報(bào)（自然科學(xué)版），38（11）：60-66. ［Liu Y Q，Zhou D Y，Gong X Y，F(xiàn)eng Y N，He S Y，Zhao W Z. 2018. MtDNA variation of Apis cerana population in the Qinghai-Tibet Plateau［J］. Journal of Shanxi Agricultural University （Natural Science Edition），38（11）：60-66.］ doi：10.13842/j.cnki.issn1671-8151.201807015.

駱群，周姝婧，徐新建，朱翔杰，經(jīng)翩翩，周冰峰. 2015. 貴州東方蜜蜂形態(tài)遺傳分析［J］. 福建農(nóng)林大學(xué)學(xué)報(bào)（自然科學(xué)版），44（3）：298-302. ［Luo Q，Zhou S J，Xu X J，Zhu X J，Jing P P，Zhou B F. 2015. Morphometric genetic diffe-rentiation of Apis cerana from Guizhou Province［J］. Journal of Fujian Agriculture and Forestry University （Natural Science Edition），44（3）：298-302.］ doi：10.13323/j.cnki.j.fafu（nat.sci.）.2015.03.014.

韋小平，賀興江，向華，任建峰，鄭鋼，龍立炎，吳建勇. 2020. 貴州蜂產(chǎn)業(yè)發(fā)展現(xiàn)狀與趨勢［J］. 貴州農(nóng)業(yè)科學(xué)，48（10）：87-92. ［Wei X P，He X J，Xiang H，Ren J F，Zheng G，Long L Y，Wu J Y. 2020. Development status and trend of apiculture in Guizhou Province［J］. Guizhou Agricultural Sciences，48（10）：87-92.］ doi：10.3969/j.issn. 1001-3601.2020.10.019.

徐祖蔭. 1986. 貴州省中蜂資源調(diào)查及其開發(fā)利用研究——貴州中蜂不同生態(tài)類型的劃分、評(píng)價(jià)及開發(fā)利用意見［J］. 貴州畜牧獸醫(yī)科技，（1）：13-19. ［Xu Z Y. 1986. Investigation and development and utilization of Apis cerana resources in Guizhou Province—Division，evaluation and development and utilization of different ecological types of Apis cerana in Guizhou Province［J］. Guizhou Journal of Animal Husbandry and Veterinary Medicine，（1）：13-19.］

徐祖蔭. 2009. 蜂海求索：徐祖蔭養(yǎng)蜂論文集［M］. 貴陽：貴州科技出版社. ［Xu Z Y. 2009. Apiculture research—Xu Zuyin essays collection［M］. Guiyang：Guizhou Science and Technology Publishing House.］

楊冠煌. 2001. 中華蜜蜂［M］. 北京：中國農(nóng)業(yè)科技出版社. ［Yang G H. 2001. Chinese honeybee［M］. Beijing：China Agricultural Science and Technology Press.］

楊冠煌. 2009. 中華蜜蜂在我國森林生態(tài)系統(tǒng)中的作用［J］. 中國蜂業(yè)，60（4）：5-7. ［Yang G H. 2009. The effect of Apis cerana cerana on forest ecosystems［J］. Apiculture of China，60（4）：5-7.］ doi：10.3969/j.issn.0412-4367.2009.04.001.

于瀛龍，周姝婧，徐新建，朱翔杰，楊凱杰，陳道印，周冰峰. 2017. 貴州省東方蜜蜂微衛(wèi)星DNA遺傳分化與遺傳多樣性分析［J］. 福建農(nóng)林大學(xué)學(xué)報(bào)（自然科學(xué)版），46（3）：323-328. ［Yu Y L，Zhou S J，Xu X J，Zhu X J，Yang K J，Chen D Y，Zhou B F. 2017. Genetic diversity and genetic differentiation of Apis cerana in Guizhou Province of southwest China［J］. Journal of Fujian Agriculture and Forestry University （Natural Science Edition），46（3）：323-328.］ doi：10.13323/j.cnki.j.fafu（nat.sci.）.2017.03.015.

張桂寧，方弟安，薛向平，毛成誠，彭云鑫. 2022. 基于線粒體COII基因序列的長江下游翹嘴群體遺傳多樣性分析［J］. 南方農(nóng)業(yè)學(xué)報(bào)，53（7）：2025-2032. ［Zhang G N，F(xiàn)ang D A，Xue X P，Mao C C，Peng Y X. 2022. Genetic diversity analysis of Culter alburnus populations in the lower reach of Yangtze River based on COII gene sequences［J］. Journal of Southern Agriculture，53（7）：2025-2032.］ doi：10.3969/j.issn.2095-1191.2022.07.026.

趙金玉，龍見坤，David Tembong Cham，陳祥盛，常志敏. 2022. 貴州省東方蜜蜂Apis cerana線粒體DNA遺傳多樣性研究［J］. 環(huán)境昆蟲學(xué)報(bào)，44（3）：595-605. ［Zhao J Y，Long J K，David T C，Chen X S，Chang Z M. 2022. Mitochondrial DNA diversity of Apis cerana in Guizhou Province［J］. Journal of Environmental Entomology，44（3）：595-605.］ doi：10.3969/j.issn.1674-0858.2022.03.10.

周姝婧，徐新建，朱翔杰，高景林，周冰峰. 2012. 海南中華蜜蜂線粒體DNA的遺傳多樣性［J］. 福建農(nóng)林大學(xué)學(xué)報(bào)（自然科學(xué)版），41（2）：170-175. ［Zhou S J，Xu X J，Zhu X J，Gao J L，Zhou B F. 2012. Genetic diversity of Apis cerana cerana in Hainan based on mitochondrial DNA［J］. Journal of Fujian Agriculture and Forestry University （Natural Science Edition），41（2）：170-175.］ doi：10.3969/j.issn.1671-5470.2012.02.013.

周姝婧，朱翔杰，徐新建，胡軍軍，于瀛龍，朱詩謠，熊長安，周冰峰. 2021. 廣西東方蜜蜂遺傳多樣性分析［J］. 應(yīng)用昆蟲學(xué)報(bào)，58（3）：672-684. ［Zhou S J，Zhu X J，Xu X J，Hu J J，Yu Y L，Zhu S Y，Xiong C A，Zhou B F. 2021. Genetic diversity of Apis cerana in Guangxi，China［J］. Chinese Journal of Applied Entomology，58（3）：672-684.］ doi：10.7679/j.issn.2095-1353.2021.068.

Coroian C O，Mu?oz I，Schlüns E A，Paniti-Teleky O R，Erler S，F(xiàn)urdui E M，M?rghita? L A，Dezmirean D S，Schlüns H，de la Rúa P D L，Moritz R F A. 2014. Climate rather than geography separates two European honeybee subspecies［J］. Molecular Ecology，23（9）：2353-2361. doi：10.1111/mec.12731.

Excoffier L，Lischer H E. 2010. Arlequin suite ver 3.5：A new series of programs to perform population genetics analyses under Linux and Windows［J］. Molecular Ecology Resources，10（3）：564-567. doi：10.1111/j.1755-0998.2010. 02847.x.

Fu Y X. 1997. Statistical tests of neutrality of mutations against population growth，hitchhiking and background selection［J］. Genetics，147（2）：915-925. doi：10.1093/genetics/147.2.915.

Garnery L，Cornuet J M，Solignac M. 1992. Evolutionary history of the honey bee Apis mellifera inferred from mitochondrial DNA analysis［J］. Molecular Ecology，1（3）：145-154. doi：10.1111/j.1365-294X.1992.tb00170.x.

Kumar S，Stecher G，Tamura K. 2016. MEGA7：Molecular evolutionary genetics analysis version 7.0 for bigger datasets［J］. Molecular Biology and Evolution，33（7）：1870-1874. doi：10.1093/molbev/msw054.

Larkin M A，Blackshields G，Brown N P，Chenna R，McGettigan P A，McWilliam H，Valentin F，Wallace I M，Wilm A，Lopez R，Thompson J D，Gibson T J，Higgins D G. 2007. ClustalW and ClustalX version 2.0［J］. Bioinformatics，23（21）：2947-2948. doi：10.1093/bioinformatics/btm404.

Librado P，Rozas J. 2009. DnaSP v5：A software for comprehensive analysis of DNA polymorphism data［J］. Bioinformatics，25（11）：1451-1452. doi：10.1093/bioinforma-tics/btp187.

Neves E F，F(xiàn)aita M R，Gaia L D O，Alves V V，Antonialli W F. 2011. Influence of climate factors on flight activity of drones of Apis mellifera （Hymenoptera：Apidae）［J］. Sociobiology，57（1）：107-114.

Peakall R，Smouse P E. 2012. GenAlEx 6.5：Genetic analysis in excel. Population genetic software for teaching and research—An update［J］. Bioinformatics，28（19）：2537-2539. doi：10.1093/bioinformatics/bts460.

Polatto L P，Chaud-Netto J，Alves-Junior V V. 2014. Influence of abiotic factors and floral resource availability on daily foraging activity of bees［J］. Journal of Insect Behavior，27（5）：593-612. doi：10.1007/s10905-014-9452-6.

Tajima F. 1989. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism［J］. Genetics，123（3）：585-595. doi：10.1093/genetics/123.3.585.

Zhao W Z，Tan K，Zhou D Y，Wang M，Cheng C F，Yu Z Y，Miao Y W，He S Y. 2014. Phylogeographic analysis of Apis cerana populations on Hainan Island and southern mainland China，based on mitochondrial DNA sequences［J］. Apidologie，45（1）：21-33. doi：10.1007/s13592-013-0223-8.

（責(zé)任編輯 蘭宗寶）

收稿日期：2022-07-12

基金項(xiàng)目：貴州省科技計(jì)劃項(xiàng)目（黔科合基礎(chǔ)〔2019〕1453號(hào)）；國家現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系建設(shè)專項(xiàng)（CARS-44-SYZ-22）；貴州省農(nóng)業(yè)科學(xué)院科技項(xiàng)目（黔農(nóng)科院青年科技基金〔2020〕08號(hào)，黔農(nóng)科院種質(zhì)資源〔2022〕16號(hào)）

通訊作者：于瀛龍（1989-），http：//orcid.org/0000-0002-2187-1687，博士，主要從事蜜蜂生物學(xué)研究工作，E-mail：yuyinglong2012@163.com；韋小平（1977-），http：//orcid.org/0000-0002-4999-0952，博士，研究員，主要從事蜜蜂生物學(xué)與生態(tài)學(xué)研究工作，E-mail：xiaopinggzaas@126.com

第一作者：周文才（1992-），http：//orcid.org/0000-0002-5424-3142，主要從事蜜蜂生物學(xué)研究工作，E-mail：2629779134@qq.com