于瀛龍，周姝婧，徐新建，朱翔杰，周冰峰

(福建農(nóng)林大學(xué)蜂學(xué)學(xué)院，福建福州350002)

遺傳多樣性是物種適應(yīng)環(huán)境變化的必備條件，是進(jìn)化的物質(zhì)基礎(chǔ)，反映了物種的進(jìn)化潛力［1］.遺傳多樣性和遺傳結(jié)構(gòu)是種群遺傳學(xué)、進(jìn)化生物學(xué)和保育遺傳學(xué)等生命科學(xué)的重要內(nèi)容［1－3］.在我國中華蜜蜂遺傳多樣性和遺傳結(jié)構(gòu)的研究中，云南［4］、福建［5－6］、海南［7］、重慶［8］、廣東［9］、華東［10］、東北［11］等地區(qū)中華蜜蜂的研究已開展［12－13］.通過19群樣本的測定，在長白山中華蜜蜂遺傳多樣性的研究中已發(fā)現(xiàn)了5種mtDNA單倍型種類［11，14－15］;嘗試用2群中華蜜蜂樣本的csd基因和用12群中華蜜蜂樣本的RAPD標(biāo)記，對長白山中華蜜蜂的遺傳多樣性和遺傳結(jié)構(gòu)進(jìn)行探索［16－18］.長白山地區(qū)中華蜜蜂mtDNA標(biāo)記、微衛(wèi)星標(biāo)記的遺傳多樣性與遺傳結(jié)構(gòu)的分析還需深入研究.

長白山是東北中華蜜蜂的主要分布區(qū)，是中華蜜蜂分布的東北部邊緣，具有獨特的生物學(xué)特性［19－20］.長白山中華蜜蜂在長期的進(jìn)化過程中形成了群勢強(qiáng)、耐寒力高、分蜂性弱等適應(yīng)當(dāng)?shù)丨h(huán)境的特性［21］，可在最低氣溫－40℃的環(huán)境自然越冬，越冬期長達(dá)半年［11，20］.長白山中華蜜蜂在體色、體長、前翅翅脈、肘脈指數(shù)等主要形態(tài)標(biāo)記表現(xiàn)出特殊性.相比于其他地區(qū)中華蜜蜂，長白山中華蜜蜂個體較大，體色偏黑［21－22］，肘脈指數(shù)較高，為 4.09－5.99［21，23－24］，第三徑中橫脈出現(xiàn)突出的頻率較高，為 30%－80%［22，24］.

本研究在長白山的腹地，中華蜜蜂主產(chǎn)區(qū)采集110群樣本，采用5對微衛(wèi)星引物和mtDNA tRNAleu-COII 357 bp片段，分析該地區(qū)中華蜜蜂的遺傳多樣性與遺傳結(jié)構(gòu).本研究為長白山地區(qū)中華蜜蜂資源保護(hù)和利用提供重要理論依據(jù).

1 材料與方法

1.1 試驗材料

中華蜜蜂樣本于2008年8月采自吉林省長白山區(qū)的安圖縣池北區(qū)、二道白河鎮(zhèn)、松江鎮(zhèn)地區(qū)，總計樣本數(shù)量110群蜜蜂(表1).

表1 長白山中華蜜蜂樣本信息Table 1 Information of samples from Changbai Mountains

1.2 試驗方法

每群取1只工蜂，利用UNIQ10柱式基因組DNA抽提試劑盒(生工生物上海股份有限公司)提取基因組.

PCR擴(kuò)增微衛(wèi)星位點的引物來源于NCBI網(wǎng)站(表2)，反應(yīng)條件為95℃預(yù)變性5 min;95℃變性50 s，退火30 s(退火溫度參照表2)，72℃延伸30 s，30個循環(huán);72℃延伸10 min;4℃保存.微衛(wèi)星PCR產(chǎn)物經(jīng)8%非變性聚丙烯酰胺凝膠電泳檢測.

表2 微衛(wèi)星位點信息Table 2 Information of 5 microsatellite loci

PCR擴(kuò)增mtDNA tRNAleu-COII片段，上游引物為E2:5'-GGCAGAATAAGTGCATTG-3'，下游引物H2:5'-CAATATCATTGATGACC-3'［25］.反應(yīng)條件為95 ℃預(yù)變性5 min;95 ℃變性30 s，50 ℃退火30 s，72 ℃ 延伸30 s，30個循環(huán);72℃延伸8 min;8℃保存.mtDNA PCR產(chǎn)物送至上海英駿生物技術(shù)有限公司進(jìn)行測序.為保證結(jié)果的準(zhǔn)確性，測序結(jié)果中個體數(shù)量在3只以下的單倍型進(jìn)行重復(fù)測序.

1.3 統(tǒng)計方法

根據(jù)8%非變性聚丙烯酰胺電泳結(jié)果進(jìn)行微衛(wèi)星位點的等位基因判讀，利用Ms tools軟件計算期望雜合度(expected heterozygosity，He)、觀察雜合度(observed heterozygosity，Ho)、多態(tài)信息含量(polymorphism information content，PIC)與等位基因頻率.利用Popgene 1.31軟件計算等位基因數(shù)(observed number of alleles，Na)、有效等位基因數(shù)(effective number of alleles，Ne)與香農(nóng)指數(shù)(Shannon's information index，I).

mtDNA測序序列以非編碼區(qū)“AAAATTT”序列為開始，截取357 bp的序列片段進(jìn)行分析，利用ClustalX 1.83軟件對其進(jìn)行比對，統(tǒng)計單倍型種類，利用MEGA 4.0軟件計算序列的堿基組成，利用DnaSP 5.00.07軟件計算單倍型多樣度(haplotype diversity，Hd)，核苷酸多樣度(nucleotide diversity，Pi)與平均核苷酸差異數(shù)(average number of nucleotide differences，K).

2 結(jié)果與分析

2.1 長白山中華蜜蜂遺傳多樣性

2.1.1 微衛(wèi)星位點的遺傳多樣性 長白山中華蜜蜂在A107、AP043、AP085、AP313、AT101位點整體等位基因較少，遺傳多樣性水平低.每個微衛(wèi)星位點發(fā)現(xiàn)1－4個等位基因，平均等位基因數(shù)為2.6000±1.1402，平均有效等位基因數(shù)為1.4181±0.6487.長白山中華蜜蜂5個微衛(wèi)星位點的期望雜合度為0－0.6107，平均期望雜合度為0.2110±0.2333.除A107位點的期望雜合度為0.6107外，其余4個位點的期望雜合度為0－0.1780.5 個位點的觀察雜合度為0－0.1827，平均觀察雜合度為0.0645 ±0.0718.5 個位點的平均香農(nóng)指數(shù)為 0.4031 ±0.4202(表3).

表3 長白山中華蜜蜂微衛(wèi)星DNA的遺傳多樣性參數(shù)Table 3 Genetic diversity parameters of Apis cerana cerana from Changbai Mountains by microsatellite loci

2.1.2 線粒體DNA遺傳多樣性 長白山中華蜜蜂在mtDNA tRNAleu-COII 357 bp片段T、C、A、G 4種核苷酸的平均比例分別為44.8%、9.5%、40.6%和5.0%.長白山中華蜜蜂共檢測出9種單倍型，H1單倍型(Japan1型，AB078733)是文獻(xiàn)報道中分布較為廣泛的類型，相比于H1單倍型，1種單倍型在1個位點上發(fā)生缺失，7種單倍型分別在1個位點發(fā)生突變(表4)，其中單一多態(tài)位點4個，簡約信息位點3個.

表4 長白山中華蜜蜂mtDNA tRNAleu-COII片段單倍型種類及變異信息1)Table 4 Variable sites and haplotypes of Apis cerana cerana from Changbai Mountains in mtDNA tRNAleu-COII region

長白山中華蜜蜂mtDNA tRNAleu-COII 357 bp片段遺傳多樣性偏低，單倍型多樣度為0.249±0.054、核苷酸多樣度為0.00074 ±0.00017、平均核苷酸差異為0.263.

本研究在長白山地區(qū)新發(fā)現(xiàn) H2、H3、H4、H5、H7、H8、H9等7種單倍型.與 NCBI網(wǎng)站進(jìn)行比對，H3(登錄號:KF673779)、H5(登錄號:KF673782)、H7(登錄號:KF673780)、H8(登錄號:KF673781)等4種單倍型為中華蜜蜂中首次發(fā)現(xiàn)(表4)，H2、H4、H9等3種單倍型為在長白山地區(qū)發(fā)現(xiàn)的新記錄.

2.2 長白山中華蜜蜂遺傳結(jié)構(gòu)

2.2.1 微衛(wèi)星遺傳結(jié)構(gòu) 長白山中華蜜蜂5個微衛(wèi)星位點的遺傳結(jié)構(gòu)中，均以一種等位基因為主.AP043、AP085、AP313、AT101等4個位點的主要主等位基因，占所檢測到的等位基因90%以上.AP313位點只有1種等位基因(342 bp).A107有4種等位基因，其中165 bp占所檢測到的等位基因50%以上(表5).2.2.2 線粒體遺傳結(jié)構(gòu) 長白山中華蜜蜂mtDNA tRNAleu-COII 357 bp片段的遺傳結(jié)構(gòu)較單一.9種單倍型中，以H1(Japan1型)單倍型為主，占83.64%.其次為H9單倍型，占6.36%.其余7種單倍型頻率為0.91%－2.73%(表 4).

表5 長白山中華蜜蜂的微衛(wèi)星遺傳結(jié)構(gòu)Table 5 Genetic structure of Apis cerana cerana based on microsatellite in Changbai Mountains

3 討論

(1)在現(xiàn)有中華蜜蜂遺傳多樣性研究中，中華蜜蜂 A107、AP043、AP085、AP313、AT101 這 5 個微衛(wèi)星位點的 PIC 值主要范圍分別是0.5－0.6、0.2－0.5、0.7－0.85、0.2－0.4 和 0.2－0.5［5－6，26］，長白山中華蜜蜂的 PIC 值分別為 0.5541、0.1186、0.1319、0和0.1644，均低于PIC值的主要分布范圍或處于較低水平.期望雜合度也體現(xiàn)出相似的規(guī)律［5－6，26］.目前，較少使用mtDNA tRNAleu-COII片段進(jìn)行中華蜜蜂遺傳多樣性研究，僅有海南中華蜜蜂的相關(guān)研究中運用了該片段.海南中華蜜蜂 mtDNA tRNAleu-COII片段發(fā)現(xiàn)單倍型多樣度、核苷酸多樣度、平均核苷 酸 差異分別為 0.742 ± 0.017、0.00354 ±0.00013和1.248［7］，而長白山中華蜜蜂的這3個指標(biāo)均遠(yuǎn)低于海南中華蜜蜂的研究結(jié)果.所以微衛(wèi)星和線粒體分子遺傳標(biāo)記分析，均說明長白山中華蜜蜂遺傳多樣性水平處于非常低的水平.

(2)長白山植被良好、具有大宗蜜源以及豐富的輔助蜜粉源，是中華蜜蜂優(yōu)良的棲息地.20世紀(jì)前長白山野生中華蜜蜂資源極豐富，此地蜂蜜是皇家貢品［27］.由于西方蜜蜂引進(jìn)后引起種間競爭，導(dǎo)致長白山中華蜜蜂種群數(shù)量下降;外來中華蜜蜂基因滲入，使得當(dāng)?shù)刂腥A蜜蜂面對環(huán)境變化、疾病等惡劣條件的耐受性下降;經(jīng)濟(jì)較發(fā)達(dá)的鄰國對野生蜂蜜的高價需求，促進(jìn)毀巢取蜜的行為，毀滅性減少了當(dāng)?shù)匾吧腥A蜜蜂的種群數(shù)量;社會經(jīng)濟(jì)的發(fā)展，養(yǎng)蜂人的數(shù)量減少，養(yǎng)蜂技術(shù)下降;人工育王技術(shù)應(yīng)用不當(dāng)，以上原因均使得長白山中華蜜蜂種群遺傳多樣性降低.長白山中華蜜蜂形成分散的小種群，遺傳變異低于其它地區(qū)的中華蜜蜂，并表現(xiàn)出較為單一的遺傳結(jié)構(gòu).

(3)長白山中華蜜蜂是分布在我國東北地區(qū)的中華蜜蜂代表，是耐寒型中華蜜蜂的寶貴資源.原長白山中華蜜蜂廣泛分布在黑龍江省的小興安嶺以南、吉林省和遼寧省的長白山脈，具有維持強(qiáng)群、耐寒力強(qiáng)等優(yōu)良蜂種特性.但近年來長白山中華蜜蜂數(shù)量大幅度減少，已成為中華蜜蜂中數(shù)量最少的生態(tài)型，處于瀕危維持狀態(tài)［9，28］.長白山中華蜜蜂存在著滅絕的風(fēng)險.生命資源的滅絕將會永久消失，因此搶救性保護(hù)長白山中華蜜蜂刻不容緩.建議在長白山中華蜜蜂分布較多、生態(tài)條件較好的地方設(shè)立中華蜜蜂保護(hù)區(qū);在保護(hù)區(qū)內(nèi)政府用補(bǔ)貼、技術(shù)支持等形式鼓勵當(dāng)?shù)剞r(nóng)民原始飼養(yǎng)本土的中華蜜蜂，擴(kuò)大長白山中華蜜蜂的種群數(shù)量;禁止西方蜜蜂和轉(zhuǎn)地的中華蜜蜂進(jìn)入，禁止從保護(hù)區(qū)外購買蜂種或蜂王，保持長白山中華蜜蜂的遺傳特征;在保護(hù)區(qū)內(nèi)適當(dāng)?shù)挠杏媱澋姆淙航涣鳎贿M(jìn)行人工選育，保持長白山中華蜜蜂的遺傳多樣性.致謝:樣本采集中得到了楊明福等同志的幫助.本科生夏鳳枝參與部分實驗.謹(jǐn)此致謝!

［1］法蘭克漢，巴盧，布里斯科.保育遺傳學(xué)導(dǎo)論［M］.北京:科學(xué)出版社，2005:30.

［2］沈銀柱，黃占景，王正詢，等.進(jìn)化生物學(xué)［M］.北京:高等教育出版社，2008:112－113.

［3］徐晉麟，趙耕春.基礎(chǔ)遺傳學(xué)［M］.北京:高等教育出版社，2009:282－284.

［4］殷玲，吉挺，陳國宏.利用微衛(wèi)星標(biāo)記分析云南省6個東方蜜蜂(Apis cerena)群體遺傳多樣性及遺傳分化［J］.西南農(nóng)業(yè)學(xué)報，2011，24(2):772－778.

［5］朱翔杰，徐新建，周妹婧，等.武夷山自然保護(hù)區(qū)中華蜜蜂微衛(wèi)星DNA遺傳分析［J］.福建農(nóng)業(yè)學(xué)報，2011，26(6):935－940.

［6］朱翔杰，周冰峰，吳顯達(dá)，等.福建中華蜜蜂微衛(wèi)星標(biāo)記的遺傳多樣性分析［J］.福建農(nóng)林大學(xué)學(xué)報:自然科學(xué)版，2011，40(4):407－411.

［7］周姝婧，徐新建，朱翔杰，等.海南中華蜜蜂線粒體DNA的遺傳多樣性［J］.福建農(nóng)林大學(xué)學(xué)報:自然科學(xué)版，2012，41(2):170－175.

［8］曹聯(lián)飛，姬聰慧，任勤，等.重慶金佛山地區(qū)東方蜜蜂遺傳多樣性研究［J］.西南師范大學(xué)學(xué)報，2010，35(5):162－166.

［9］成純富，周丹銀，劉意秋，等.基于mtDNA Cytb基因序列分析廣東省東方蜜蜂遺傳多樣性［J］.蜜蜂雜志，2012(9):4－6.

［10］吉挺，殷玲，劉敏，等.華東地區(qū)中華蜜蜂六地理種群的遺傳多樣性及遺傳分化［J］.昆蟲學(xué)報，2009，52(4):413－419.

［11］李志勇，薛運波，趙惠燕.長白山區(qū)中華蜜蜂mtDNA非編碼區(qū)至COⅠ基因PCR及序列分析［J］.湖北農(nóng)業(yè)科學(xué)，2009，48(5):1051－1054.

［12］YIN L，JI T，LIU M，et al.Genetic diversity and relationship between genetic distance and geographical distance of 6 Apis cerana cerana populations in China［J］.Research Journal of Animal Sciences，2008，2(6):183－187.

［13］SMITH D R，VILLAFUERTE L，OTIS G，et al.Biogeography of Apis cerana F.and A.nigrocincta:insight from mtDNA studies［J］.Apidologie，2000，31(2):265－279.

［14］TAN K，WARRIT N，SMITH D R.Mitochondrial DNA diversity of Chinese Apis cerana［J］.Apidologie，2007，38(3):238－246.

［15］姜玉鎖，趙慧婷，姜俊兵，等.中國境內(nèi)不同地理型東方蜜蜂線粒體tRNAleu-COII基因多態(tài)性研究［J］.中國農(nóng)業(yè)科學(xué)，2007，40(7):1532－1542.

［16］劉志勇，王子龍，王歡，等.中華蜜蜂csd多態(tài)性分析［J］.中國農(nóng)業(yè)科學(xué)，2011，44(23):4911－4917.

［17］薛運波，李興安，葛鳳晨，等.長白山中華蜜蜂基因組DNA多態(tài)性的研究［J］.中國農(nóng)業(yè)科學(xué)，2007(40):426－432.

［18］薛運波，李興安，葛鳳展，等.長白山中華蜜蜂基因組CBS-700(bp)DNA同源片段克隆及其Southern雜交的研究［J］.中國蜂業(yè)，2006，57(10):8－11，17.

［19］楊冠煌.中華蜜蜂［M］.北京:中國農(nóng)業(yè)科技出版社，2001:13.

［20］龔一飛，張其康.蜜蜂分類與進(jìn)化［M］.福建:科學(xué)技術(shù)出版社，2000:21.

［21］楊冠煌，李桂仙.中華蜜蜂分布和變異(五)—東北地區(qū)中華蜜蜂資源［J］.中國養(yǎng)蜂，1983(5):19.

［22］楊冠煌，許少玉，匡邦郁，等.東方蜜蜂Apis cerana Fab在我國的分布及其亞種分化［J］.云南農(nóng)業(yè)大學(xué)學(xué)報，1986，12(1):89－92.

［23］張連江，樊賢，譚墾，等.長白山東方蜜蜂形態(tài)特征研究［J］.云南農(nóng)業(yè)大學(xué)學(xué)報，2006，21(4):511－516.

［24］葛鳳晨，歷延芳，柏建民，等.長白山中華蜜蜂分布及其生產(chǎn)效率的研究［J］.中國養(yǎng)蜂，2002，53(6):4－7.

［25］GARNERY L，COMUET J M，SOLIGANAC M.Evolutionary history of the honey bee Apis mellifera inferred from mitochondrial DNA analysis［J］.Moleculor Ecology，1992，1(3):145－154.

［26］徐新建，周姝婧，朱翔杰，等.海南島中華蜜蜂遺傳多樣性的微衛(wèi)星DNA分析［J］.昆蟲學(xué)報，2013，56(5):554－560.

［27］葛鳳晨，陳東海，歷延芳，等.長白蜜蜂文化研究［M］.長春:吉林科學(xué)技術(shù)出版社，2006:1－21.

［28］國家畜禽遺傳資源委員.中國畜禽遺傳資源志:蜜蜂志［M］.北京:中國農(nóng)業(yè)出版社，2011:45－49.