蘇曉玲,陳道印,曹聯飛,華啟云,趙東緒,羅谷輝
(1.金華市農業科學研究院,浙江 金華 321017;2.浙江省農業科學院 畜牧獸醫研究所,浙江 杭州 310021)
東方蜜蜂(Apiscerana)是我國重要的本土蜂種資源,在維持生態系統平衡和促進農作物增產方面起著重要作用。東方蜜蜂在我國分布廣泛,并在長期自然選擇中對地理和生態環境產生了適應性,形成了不同的生態型[1]。研究證明,東方蜜蜂的形態特征在不同地域表現出明顯的變異[2-5],而且受到海拔和氣候等因素的影響[6-9]。20世紀初東方蜜蜂種群銳減,蜜蜂種群遺傳多樣性下降[10],21世紀以來東方蜜蜂種群數量開始快速增加,但其種質資源多樣性不容樂觀,蜂群跨地區運輸和蜜蜂隨意引種等原因導致其遺傳多樣性降低。為了科學保護和利用東方蜜蜂蜂種資源,種質資源評價和遺傳多樣性研究是當務之急[11],而形態特征測定是蜜蜂種質資源評價最常用和最基礎的手段[12-13]。
金華為浙中丘陵盆地地貌,盆周向盆地中心呈現出中山、低山、丘陵崗地、河谷平原階梯式層狀分布,地理環境非常適宜東方蜜蜂生存,境內東方蜜蜂資源豐富。該地區屬亞熱帶季風氣候,光熱水條件優越,但時空分布不均勻,小氣候多樣,有一定垂直差異。近年來,在政策扶持和項目帶動下,浙江省東方蜜蜂種群數量成倍增長[14],金華市東方蜜蜂種群數量由2012年的1萬余群增長到7萬余群[15]。然而目前還沒有關于金華市東方蜜蜂內部遺傳分化的研究。本研究通過蜜蜂形態特征分析金華東方蜜蜂資源遺傳分化及地理環境對形態的影響,為該地東方蜜蜂資源保護和開發利用提供科學參考依據,也為我國東方蜜蜂資源的研究提供重要基礎資料。
于2019年9—11月,在浙江省金華市8個樣點采集東方蜜蜂(圖1和表1),每個樣點采集5群,共40群。每群采集50只成年工蜂,以75%的酒精浸泡帶回實驗室,低溫保存備用。同時,收集了廣東東方蜜蜂(GD)和吉林東方蜜蜂(JL)各15群作為外群分析樣本。

圖1 浙江金華市8個東方蜜蜂采樣點

表1 金華市東方蜜蜂樣本采集點信息
用帶有圖形采集器的體式顯微鏡(型號:GL-99TI)與電腦形態測量系統進行形態測定。每個樣點解剖40只成年工蜂,每群隨機解剖8只,采集標準形態標記圖像,測定16項與形狀大小相關的形態指標,測定的指標包括吻長、右前翅(翅長、翅寬)、右后足(脛節長、股節長、基跗節長、基跗節寬)、背板(第三背板長、第四背板長、第五背板長、第六背板長)、腹板(第四腹板長、蠟鏡長、第七腹板長、第七腹板寬)和第五背板絨毛帶長。測定方法參照Ruttner1988年提出的測定標準進行[16]。
對東方蜜蜂形態數據利用SPSS23.0軟件進行分析。將各樣點形態指標作單因素方差(ANOVA)分析,根據樣本均值和標準差分析種群形態特征,比較兩兩之間差異顯著性,分別作顯著性水平0.05、0.01檢驗。采用主成分分析和聚類分析研究金華東方蜜蜂與外群東方蜜蜂的形態分化。采用逐步判別和聚類分析研究金華地區東方蜜蜂種群形態分化。采用相關分析法分析不同樣點的形態學特征與其生態環境中的不同生態因子的關系。
金華東方蜜蜂8個樣點及外群形態數據的平均值與標準誤見表2。金華東方蜜蜂16個形態指標均極顯著小于吉林樣本,除翅長、翅寬、基跗節寬和第三背板長與廣東樣本差異不顯著外,其他形態指標與廣東樣本有一定顯著差異。

表2 金華市東方蜜蜂8個樣點及外群形態數據的平均值與標準誤 單位:mm
金華東方蜜蜂內部變異分析表明,16個形態特征之間均存在極顯著差異,武義新宅、婺城北山、蘭溪馬澗樣本的吻長顯著大于其他樣點;婺城塔石樣本具有最小的翅長(8.53±0.12)mm和翅寬(2.90±0.07)mm,且與武義宣平樣本差異極顯著;武義新宅和磐安方前2個樣本具有較大的第三背板長(2.07±0.10和2.07±0.08)mm、第四背板長(1.86±0.05和1.86±0.04)mm、第五背板長(1.81±0.05和1.81±0.05)mm和第六背板長(1.84±0.08和1.79±0.08)mm,且與永康方巖、婺城塔石、武義宣平和蘭溪馬澗樣本差異顯著;永康方巖樣本具有最小的第四背板長(1.77±0.08)mm、第五背板長(1.73±0.07)mm、第六背板長(1.73±0.09)mm和第四腹板長(2.49±0.08)mm,除第四腹板長外,與婺城北山、武義新宅和磐安方前樣本差異顯著;武義新宅樣本蠟鏡長平均值最大(1.26±0.05)mm,與其他樣本差異顯著。
將金華8個樣點和2個外群東方蜜蜂的16個形態遺傳標記數據進行主成分分析,從中提出2個特征值高的因素(>1)。因素1特征根為7.093,解釋總體變異的44.3%,第四腹板長、第七腹板寬、第七腹板長、蠟鏡長、第三背板長、第四背板長、第五背板長、第六背板長、第七背板長、第五背板絨毛帶長、翅長、翅寬和吻長對因素1貢獻較大。因素2特征根1.282,解釋總體變異的8.0%,形態標記脛節長、股節長、基跗節長和基跗節寬對因素2貢獻較大。以上2個因素包括52.4%的數據變異。
對因素1和因素2進行分析。由圖2可知,所有樣本可以分為2個種群:金華東方蜜蜂和廣東樣本為1個地理種群,吉林樣本單獨為1個種群。金華東方蜜蜂的因素1的值比吉林樣本要小,由于因素1所涉及的性狀主要與體型大小有關,因此說明,金華東方蜜蜂個體比吉林樣本個體要小。
對金華8個樣點和2個外群東方蜜蜂的16個形態遺傳標記數據進行聚類分析。由圖3可知,金華8個樣點與廣東樣本聚為一支,吉林樣本單獨聚為一支。結果說明,金華東方蜜蜂與吉林樣本遺傳距離遠,而浙江金華地區東方蜜蜂內部遺傳距離小,金華東方蜜蜂與廣東樣本具有相似的遺傳特征。

圖3 聚類分析結果
根據Wilks’ Lambda統計量檢驗,9個判別效果極顯著的形態標記進入判別函數,按判別函數特征值從大到小依次為吻長、第四背板長、蠟鏡長、股節長、第四腹板長、第七腹板長、第五背板長、第五背板絨毛帶長和第七腹板長。逐步判別分析建立的判別函數共7個,前3個判別函數分別解釋總體變異的39.5%、26.9%和16.7%,累積貢獻率為83.2%。
利用貢獻率最高的2個函數,分析金華地區東方蜜蜂內部的遺傳分化,逐步判別分析散點圖4顯示,函數1橫軸方向樣點明顯分離成兩組,武義新宅、磐安方前、婺城北山、蘭溪馬澗東方蜜蜂樣點的重心較近,婺城長山、武義宣平、婺城塔石、永康方巖東方蜜蜂4個樣點的重心較近。利用判別分析組質心處的函數進行聚類分析,武義新宅、磐安方前、婺城北山、蘭溪馬澗聚為一支,婺城長山、武義宣平、婺城塔石、永康方巖聚為一支。結果說明,金華地區東方蜜蜂內部分化成2個種群。

圖4 金華市東方蜜蜂逐步判別分析
對金華8個樣點東方蜜蜂的形態特征與環境因子進行相關性分析。由表3可知,16項形態特征中,有12項指標與環境因子呈一定相關性,9項指標與海拔的相關系數顯著,其中吻長、翅寬、第四背板長、第六背板長、蠟鏡長、第七腹板長和第七腹板寬呈顯著正相關,脛節長和基跗節長呈顯著負相關。6項指標與年平均氣溫的相關系數極顯著,其中翅寬、蠟鏡長第七腹板寬呈顯著負相關,脛節長、基跗節長和第七腹板長呈顯著正相關;9項指標與年最低氣溫呈極顯著負相關,包括吻長、翅寬、股節長、基跗節寬、第四背板長、第五背板長、第六背板長、蠟鏡長和第七腹板寬。

圖5 金華市東方蜜蜂判別函數重心聚類

表3 金華市東方蜜蜂形態特征與環境因子相關性
金華東方蜜蜂16個形態指標均極顯著小于吉林樣本,表明金華東方蜜蜂體型小于吉林樣本;散點圖和聚類分析結果顯示,金華東方蜜蜂與廣東樣本聚成一支,吉林樣本單獨聚為一支,群體距離與地理距離相符。譚墾等[17-18]研究表明,東北東方蜜蜂受低溫環境影響進化出了較大的體型,本研究與其結果一致。另外,金華東方蜜蜂與吉林、廣東樣本形態數據對比分析也說明了本研究所測數據的準確性和可靠性。
根據國家畜禽遺傳資源志東方蜜蜂類型的劃分標準,我國東方蜜蜂分為9個不同類型,廣東東方蜜蜂屬于典型的華南中蜂類型,而浙江東方蜜蜂分成2個生態型,浙江西部為華中中蜂類型,浙江大部分地區為華南中蜂類型。金華地處浙江省中部,處于2個生態型分布的過渡地帶。本研究表明,金華東方蜜蜂除翅長、翅寬、基跗節寬和第三背板長與廣東樣本差異不顯著外,其他形態指標與廣東樣本顯著差異。由于缺乏華中中蜂形態樣本,本研究未明確金華地區東方蜜蜂與華中中蜂的關系,金華東方蜜蜂與2個生態型之間的關系需要進一步的研究。
金華地區東方蜜蜂變異分析表明,16個形態特征之間均存在極顯著的變異性,說明金華地區東方蜜蜂種群內遺傳變異豐富。判別—聚類分析顯示金華8個樣點分成2個種群:婺城北山、武義新宅、磐安方前與蘭溪馬澗為1個種群;永康方巖、婺城塔石、婺城長山與武義宣平為1個種群。婺城北山、武義新宅與磐安方前3個樣點彼此之間的直線距離均在60 km以上,卻聚為一支;而這3個樣點與婺城長山、永康方巖、婺城塔石樣點之間直線距離僅為15~62 km,這表明地理距離不是影響金華地區東方蜜蜂遺傳變異的主要原因。婺城北山、武義新宅為海拔最高的2個樣點,年均氣溫與年最低氣溫在所有樣點中處于非常低的水平,磐安方前的年均氣溫與最低氣溫在所有樣點中也處于較低的水平。需要注意的是蘭溪馬澗樣點的海拔與氣溫與其他3個樣點并不一致,但仍在同一個分組中,可能是由于其與婺城北山樣點之間的直線距離為8 km,而蜂王與雄蜂的交配距離約7.5 km[19],兩地基因流交流頻繁導致的。永康方巖、婺城塔石、婺城長山與武義宣平相較則海拔較低、氣溫較高。本研究結果表明,海拔與溫度可能是影響金華地區東方蜜蜂遺傳變異的主要原因,這與王桂芝等[6,8,20]的研究結果一致。
在金華東方蜜蜂樣本中,武義新宅、磐安方前和婺城北山3個樣本具有最大的第四背板長、第五背板長和第六背板長,而且武義新宅樣本的蠟鏡長最長,磐安方前樣本第五背板絨毛帶長最長,婺城北山樣本翅寬最寬。工蜂背板長是衡量蜜蜂腹部大小的相對標準,因此,這3個地區的蜜蜂具有較大的腹部。而武義新宅、磐安方前和婺城北山樣本的年均氣溫與年最低氣溫較低,這可能是金華東方蜜蜂對低溫環境進行的適應性進化,這與譚墾等[8]的研究結果一致。
環境因子對金華東方蜜蜂形態特征的影響表現明顯,表明金華東方蜜蜂具有比較明顯的地理變異性。本研究16項形態特征中,有12項指標與環境因子呈顯著相關性,其中翅寬、蠟鏡長和第七腹板寬同時與海拔呈正相關,與年平均氣溫和年最低氣溫呈極顯著負相關,說明金華地區東方蜜蜂這3個指標受環境影響大,海拔越高,溫度越低,這3個指標越大。