太湖地區與皖南地區中華蜜蜂微衛星DNA遺傳多樣性分析

吉挺 沈芳 孟祥金 王帥 李文明

摘要：利用篩選后的2對熒光標記微衛星引物，研究太湖地區與皖南山區中華蜜蜂（以下簡稱中蜂）間的遺傳多樣性并分析其遺傳分化。研究檢測判定了60個中蜂個體的基因型，計算2個群體的優勢等位基因頻率（i）、期望雜合度（He）、多態信息含量（IC）、群體內近交系數（Fis），并分析中蜂群體內與群體間的遺傳變異，采取配對試驗均數差異t檢驗比較太湖中蜂與皖南中蜂群體的遺傳差異。結果表明，太湖地區與皖南山區2個中蜂群體間的i、He、IC、Fis差異均不顯著（i=0356>005，He=039>005，IC=0260>005，Fis=0428>005）。可見，太湖地區與皖南地區兩地中蜂群體存在一定的遺傳多樣性，但遺傳分化程度不大。本研究結果對江蘇與安徽兩地中華蜜蜂品種的選育和保護也具有一定的指導意義。

關鍵詞：太湖地區；皖南地區；中華蜜蜂；微衛星DNA標記；遺傳多樣性

中圖分類號： S898文獻標志碼： A[HK]

文章編號：002-302（204）2-002-05

中華蜜蜂（Apis cerana cerana）簡稱中蜂，指分布于我國境內東方蜜蜂地理亞種的總稱，是我國特有的遺傳資源，家養歷史悠久，2006年被列入國家畜禽保護名錄中。江蘇省自古是我國野生中蜂的聚居區，也是最早引入新法飼養的地區之一[2]。近幾十年來，由于西方蜜蜂大量引入及自然條件的破壞，江蘇省中蜂遺傳資源已消失迨盡，目前僅在環太湖地區略有分布，種質資源保護迫在眉睫[3]。同時，隨著近幾年環太湖地區設施農業的快速發展，蜜蜂授粉尤其是中蜂授粉的需求越來越多，現有的中蜂資源已不能滿足農業生產的需求。皖南地區地形復雜，植被豐富，適宜中蜂棲息，中蜂資源十分豐富[4]，而且皖南地區與太湖地區相隔較近，沒有任何地理隔離，所以從皖南地區引入蜂群可能是解決太湖地區中蜂資源缺乏的有效方法。但是引入后是否會改變太湖地區原有中蜂的遺傳結構，繼而對當地種質資源保護造成負面影響是應該考慮的首要問題，所以需要對2個地區現有中蜂資源的遺傳結構與遺傳多樣性進行比較分析。

微衛星是目前普遍使用的分子遺傳標記，具有多態性好、共顯性、易于鑒定、檢測重復性好、在基因組分布廣泛等優點。聯合國糧農組織（FAO）在其持續發展和管理動物遺傳資源的戰略計劃中，推薦將微衛星標記作為優先考慮的分析工具[5-6]。近幾年來，微衛星標記已廣泛應用于我國各地蜜蜂種質資源的分析中，但是目前普遍采用聚丙烯酰胺電泳加放射顯影或銀染的方法，不僅費時費力效率低，而且誤差較大[7-9]。本研究利用篩選后的2對熒光標記微衛星引物對太湖地區與皖南地區的中蜂進行遺傳多樣性比較和遺傳分化分析，旨在初步評價2個地區中蜂群體遺傳結構的差異，繼而探討江蘇省太湖地區異地引入蜂群的可能性。

材料與方法

試驗材料

太湖地區中蜂（DS）來自于江蘇省蘇州市吳中區東山鎮，皖南山區中蜂（X）來自于安徽省宣城市涇縣山區，每個地區各隨機選取30群，蜂群盡量以野生種群為主，每群隨機采集成年工蜂0只，用無水乙醇溶液浸泡保存。中蜂群體的采集信息見表。

參照吉挺所建立的方法[0]提取中蜂基因組DNA，用微量紫外可見光度計（NanoDrop ND-000）測定其含量與純度，-20 ℃ 保存備用。

2微衛星引物的篩選及CR條件的優化

根據筆者所在課題組對中華蜜蜂轉錄組測序結果所獲得的3 000多個微衛星位點進行篩選，每條染色體選擇2～3對微衛星引物，共初步選取54對微衛星引物，送往上海生工生物工程有限公司合成。根據GenBank和相應文獻提供的微衛星引物對所選54對引物進行優化，進一步篩選出 30對擴增效果較好的微衛星引物，具體所選微衛星引物信息見表2。

由于CR反應中影響因素較多，試驗分別設計了以DNA模板濃度、Mg2濃度、Taq酶用量、退火溫度為因素的梯度試驗，最終確定CR的反應體系：0×Buffer 20 μL，25 mmol/L MgCl2 0 μL，0 mmol/L dNTs 05 μL，0 pmol/μL 上游引物 0 μL，0 pmol/μL 下游引物 0 μL，5 U/μL Taq DNA聚

配制3%瓊脂糖凝膠，點樣6 μL，20 V電泳20 min，檢查產物的有無。如有產物，配制8%聚丙烯酰胺凝膠80 mL體系，00 V預電泳5 min，點樣0 μL，20 V電泳3 h。進行硝酸銀染色，直至出現清晰的等位基因條帶。

4熒光引物的篩選和組合

根據聚丙烯酰胺凝膠電泳的結果，盡量選擇不同染色體上的標記，淘汰掉無法擴增以及無多態性的引物，最終選出相對較理想的2對多態性較豐富的微衛星引物：7-CL278Contig（）、5-CL4Contig（γ）、39-CL308Contig6（X）、7-CL462Contig5（S）、8-CL470Contig3（T）、37-CL293Contig（W）、29-CL229Contig33（K）、28-CL650Contig3（B）、33-CL360Contig4（L）、-CL229Contig27（M）、4-CL549Contig3（R）、-CL33Contig4（Q）。每對引物的上游 5′ 端分別用熒光染料FAM（藍色）、HEX（綠色）、TAMRA（黃色）進行標記，獲得的熒光標記引物避光保存。按照微衛星引物的堿基長度進行組合，為了便于區分得精確一些，組合的引物長度至少相差 20 bp，獲得的5個熒光標記微衛星引物組合信息見表3。

CR熒光產物STR分型送至上海生工生物工程有限公司進行，使用ABI-3730XL DNA Analyzer全自動測序儀檢測。上樣的總體積為35 μL，其中的上樣液Hi-Di Formamide 0 μL、GS-500 Size Standard 05 μL，另外的3 μL為混合CR產物。將Hi-Di Formamide 和GS-500 Size Standard 混合均勻后與等量的同一個體的CR產物進行混合，然后進行個體編號，依次加樣到96孔板中，接著變性和檢測。檢測結束后，使用GeneMapper 40軟件自動生成獨立的圖譜文件，讀取片段長度、峰值和峰面積，判斷純合子和雜合子（單峰為純合子，雙峰為雜合子），最后導出Excel數據表格。endprint

6遺傳多樣性分析公式

6等位基因頻率

[Z]i=（2niinijnij2…nijn）/（2N）。

其中：i為第i個等位基因頻率；nii為第i個等位基因純合個數；jn為與i共顯的第n個等位基因；nij為含i與jn共顯的等位基因個體數；N為群體中個體數。

62雜合度

[Z]He=-∑[DD（]ki=-[DD）]2i。

其中：k為等位基因數；i為第i個等位基因頻率。

63多態信息含量

[Z]IC=-∑[DD（]ki=[DD）]2i-∑[DD（]k-i=[DD）]∑[DD（]kj=i[DD）]22i2j=2∑[DD（]k-i=[DD）]∑[DD（]kj=i[DD）]ij（-ij）。

其中：k為等位基因數；i為第i個等位基因頻率；j為第j個等位基因頻率。

64群體內近交系數

[Z]Fis=（Hs-Ho）/Hs。

其中：Ho為觀察雜合度；Hs為平均期望雜合度。

65統計分析

利用Microsatellite-Toolkit軟件根據GeneMapper 40軟件導出的Excel表格來計算等位基因頻率與期望雜合度，根據Botestein等的公式設計、計算群體的多態信息含量，再根據FSTAT程序計算群體內近交系數。

2結果與分析

2DNA提取結果

DNA提取出來以后，將其稀釋到00 ng/μL，部分瓊脂糖凝膠用紫外透視成像系統檢測，其結果如圖所示：純度比較高，無拖尾現象。

[FK（W0][TTtif][FK）]

22熒光產物STR分型結果

使用ABI-3730XL DNA Analyzer全自動測序儀進行檢測掃板，得到的部分微衛星引物的擴增結果與STR分型結果見圖2。圖2所示的是太湖地區中蜂DS0個體在B位點的擴增結果，基因型為32/34（雜合子）；皖南山區中蜂X03個體在T位點的擴增結果，基因型為32/32（純合子）。

23太湖地區與皖南山區中蜂優勢等位基因的比較

由表4可見，利用配對試驗均數差異雙尾t檢驗對2個群體優勢等位基因頻率進行處理，結果i=0356>005，即太湖地區與皖南地區中蜂群體在這2對微衛星引物標記上檢測出來的優勢等位基因差異不顯著。

24群體間相關參數分析

由表5可知，太湖地區中蜂群體的He、IC、Fis平均為055、059、0053，均略高于皖南山區中蜂群體（0473、0423、025）。對2個群體的He、IC、Fis進行t檢驗，結果顯示，2個群體的3個參數差異均不顯著（[HT5”]He=039、IC=0260、Fis=0428）。

3結論與討論

從2006年開始，陸續有研究者利用西方蜜蜂微衛星標記

對我國中蜂遺傳資源及遺傳多樣性進行分析，但普遍方法是在CR擴增后采用聚丙烯酰胺電泳加放射顯影或銀染的方法，通過人為觀察主觀地對多態性條帶進行分析，這種方法最令人質疑的在于很難將相鄰近的2個條帶準確區分開來，從而造成判型誤差[7-8]。目前利用熒光引物擴增的微衛星位點具有較高的分辨率和靈敏度，對擴增產物在單堿基水平和表達量方面也具有較高的準確性。本試驗利用ABI基因測序分析儀對熒光標記DNA片段進行檢測，并利用分子量內標對DNA片段長度進行計算，使2個地區的中蜂STR分型方法較傳統方法更精確，更能真實反映2個群體內的遺傳多樣性及群體間的遺傳分化。

本研究在筆者所在課題組對中華蜜蜂轉錄組測序結果所獲得的 3 000 多個微衛星位點基礎上篩選出2個微衛星位點，盡量保證其在多條染色體或連鎖群中都有分布，因此得出的數據有一定的代表性和可比性。群體中頻率最高的等位基因是該物種中最原始、最保守的等位基因，其余等位基因是進化[CM（25]過程中由該等位基因突變形成的，因此微衛星的多態性能

2K2]夠反映物種的進化歷史[2]。本研究發現2個群體的優勢等位基因頻率間存在差異，但差異不顯著（i=0356>005），表明太湖地區中蜂群體（DS）與皖南地區（X）各自維持著一定的種質特性，但兩者之間遺傳分化較小。這可能與這2個地區相隔距離較近、兩者基因流動值較高有關，表明從皖南山區引入部分中蜂蜂群，擴大太湖地區中蜂種群數量，進行蜂蜜生產與授粉，不會對當地中蜂遺傳資源產生較大影響。

多態信息含量是衡量基因片段多態性較好的指標，當IC>05時，該座位為高度多態性座位；當025

期望雜合度別稱基因多樣度，是一個度量群體遺傳變異的最適參數[4]。本研究所選用的微衛星標記中，、γ、T、R、Q、X、W等7個標記有較高的多態性，群體遺傳多樣性豐富，具有較高的選擇潛力；而B標記在2個群體間表現出完全不同的多態性，其在皖南山區中蜂群體中為純合子，而在太湖地區群體中表現為高度多態（He=0875），今后可將其作為候選遺傳標記區別2個地理群體。

本研究通過對太湖地區與皖南山區中蜂群體進行遺傳多樣性檢測，對2個群體主要遺傳參數進行比較，結果顯示，2個群體間的He、IC、Fis差異均不顯著，其中太湖地區中蜂群體的He、IC還略高于皖南山區，表明太湖地區中蜂遺傳多樣性仍較豐富，這可能與太湖地區蜜源植物資源豐富，適宜中蜂生存，同時該地區與種群數量較多的浙北和皖南山區沒有明顯的地理隔離，造成相互基因流動有一定的關系。這也說明在分析中蜂群體遺傳多樣性時，種群數量與其擁有的遺傳多樣性高低沒有直接聯系。

綜上所述，本研究利用篩選后的2個微衛星DNA熒光標記分析了太湖地區與皖南山區中蜂群體的遺傳多樣性與群體間的差異，結果表明，2個地理群體遺傳多樣性均較豐富，具有較高的保種價值，且兩者間遺傳分化不明顯，完全可能從皖南地區引入中蜂蜂群擴大太湖地區種群數量。

[HS23][HT85H]參考文獻：[HT8SS]

[ZK（#]國家畜禽遺傳資源委員會組 中國畜禽遺傳資源志——蜜蜂志[M] 北京：中國農業出版社，20：2-20

[2][2]吉挺，陳晶，殷靈，等 江蘇省野生及家養中華蜜蜂微衛星DNA遺傳多樣性分析 中國畜牧獸醫，2007，34（2）：39-4

[3]何明 江蘇的中蜂資源亟待保護 蜜蜂雜志，2007，27（3）：6

[4]黃思思 皖南山區中華蜜蜂遺傳多樣性研究[D] 合肥：安徽農業大學，20：-5

[5]Crooijmans R ，Groen A F，van Kampen A ，et al Microsatellite polymorphism in commericial broiler and layers lines estimated using pooled blood samples oultry Science，996，7（7）：904-909

[6]Romanov M N，Weigend S Genetic diversity in chicken populations based on microsatellite markers[C]//roceedings of the Conference “From ay Lush to Genomics：Vidions for Animal Breeding and Genetics”999：7

[7]徐新建，周姝婧，朱翔杰，等 海南島中華蜜蜂遺傳多樣性的微衛星DNA分析 昆蟲學報，203，56（5）：554-560

[8]于瀛龍，周姝婧，徐新建，等 長白山中華蜜蜂（Apis cerana cerana）遺傳多樣性分析 福建農林大學學報：自然科學版，203，42（6）：643-647

[9]陳晶，陳國宏，吉挺，等 沂蒙山中華蜜蜂微衛星DNA遺傳多樣性分析 中國蜂業，2008，59（4）：-3[ZK）]

[0][ZK（#]吉挺 中國東方蜜蜂資源遺傳多樣性分析[D] 揚州：揚州大學，2009：35-48

殷玲 東方蜜蜂抗螨相關基因的篩選及初步驗證[D] 揚州：揚州大學，203：3-32

[2]張揚，黃正洋，陳陽，等 采用微衛星DNA標記分析巢湖鴨異地小群保種效果 中國畜牧雜志，202，48（）：7-

[3][2]Vanhala T，Tuiskula-Haavisto M，Elo K，et al Evaluation of genetic variability and genetic distances between eight chicken lines using microsatellite markers oultry Science，998，77（6）：783-790

[4]Ott Analysis of human genetic linkage[M] Baltimore：ohns Hopkins University ress，200[ZK）][H][FL）]endprint