周姝婧, 朱翔杰, 徐新建, 蘇顯冰, 蔡宗兵, 周冰峰

[1.福建農(nóng)林大學(xué)動(dòng)物科學(xué)學(xué)院(蜂學(xué)學(xué)院)，福建 福州 350002；2.福建農(nóng)林大學(xué)蜜蜂研究所，福建 福州 350002]

遺傳多樣性對(duì)生物的生存和發(fā)展起到重要作用，是種群適應(yīng)環(huán)境變化的基礎(chǔ)，并將長(zhǎng)期及動(dòng)態(tài)地影響進(jìn)化潛力及適應(yīng)不同環(huán)境的能力.合理、有效地評(píng)估種群遺傳多樣性水平是種群遺傳學(xué)、遺傳多樣性、保護(hù)生物學(xué)等學(xué)科的重要工作[1-2].微衛(wèi)星標(biāo)記是相對(duì)中性的共顯性標(biāo)記，具有多態(tài)性高、擴(kuò)增穩(wěn)定和高效、易于擴(kuò)增，且價(jià)格較低等優(yōu)點(diǎn).在分子生物學(xué)實(shí)驗(yàn)技術(shù)快速發(fā)展的背景下，微衛(wèi)星標(biāo)記仍廣泛用于動(dòng)植物遺傳多樣性研究，種群遺傳分化及相關(guān)的分化機(jī)制研究，種群的遺傳結(jié)構(gòu)與特征、品種鑒定研究，以及評(píng)估外來(lái)種群對(duì)本土種群基因滲透程度等研究[3-5].

蜜蜂遺傳多樣性研究是進(jìn)行蜜蜂遺傳資源保護(hù)和利用的基礎(chǔ)，微衛(wèi)星標(biāo)記在蜜蜂遺傳多樣性領(lǐng)域上的應(yīng)用非常普遍.樣本量對(duì)蜜蜂微衛(wèi)星DNA遺傳多樣性參數(shù)的影響尚未見(jiàn)報(bào)道，該領(lǐng)域研究所需的樣本量始終處于含糊狀態(tài).對(duì)1998—2021年公開(kāi)發(fā)表的論文進(jìn)行樣本量統(tǒng)計(jì)的結(jié)果顯示，蜜蜂微衛(wèi)星DNA遺傳多樣性研究使用的樣本量從個(gè)位數(shù)到百群以上，多為11～60群.其中，樣本量為10群以下約占6.4%，11～20群約占14.4%，21～30群約占18.3%，31～40群約占10.2%，41～50群約占16.0%，51～60群約占17.3%，61～70群約占7.3%，71～80群約占3.5%，81～90群約占2.0%，91～100群約占0.9%，100群以上約占3.8%[6-44].隨著對(duì)微衛(wèi)星DNA遺傳多樣性認(rèn)識(shí)和研究的深入，研究者意識(shí)到樣本量不足可能無(wú)法充分體現(xiàn)樣本代表區(qū)域的遺傳多樣性水平，因此，樣本量小于30群的情況從2010年前的41%下降至近5年的12%.在實(shí)際樣本采集中，常常有許多困難使得樣本難以采集.因此，選擇適當(dāng)大小的樣本量，對(duì)于提高遺傳多樣性研究的效率有著非常重要的意義.

本研究對(duì)5個(gè)180群樣本量的東方蜜蜂種群，設(shè)置15個(gè)樣本量梯度進(jìn)行重復(fù)抽樣，探討在遺傳多樣性研究中，各遺傳多樣性常用參數(shù)在不同樣本量時(shí)與總體的關(guān)系，以及樣本量對(duì)遺傳多樣性參數(shù)檢測(cè)穩(wěn)定性的影響.探索蜜蜂遺傳多樣性研究的最適宜樣本量，可為其他生物的種群遺傳學(xué)研究提供采集策略借鑒.

1 材料與方法

1.1 樣本的采集與保存

在東方蜜蜂分布面積大、蜂群保有量高的陜西省、重慶市、江西省、廣西壯族自治區(qū)和保有量不高的臺(tái)灣省開(kāi)展樣本采集.每個(gè)省份為一個(gè)種群，每個(gè)種群的樣本量為180群.供試樣本以野生蜂群或半野生原始飼養(yǎng)蜂群為主，采樣蜂場(chǎng)無(wú)引種、選育行為，可代表當(dāng)?shù)孛鄯溥z傳結(jié)構(gòu)和多樣性.

1.2 基因組提取與微衛(wèi)星DNA擴(kuò)增

每個(gè)蜂群取1只工蜂提取全基因組進(jìn)行微衛(wèi)星DNA擴(kuò)增.微衛(wèi)星標(biāo)記利用ExTaq Hot Start Version擴(kuò)增試劑盒(大連寶生生物工程有限公司)進(jìn)行36個(gè)微衛(wèi)星標(biāo)記擴(kuò)增試驗(yàn)，引物設(shè)計(jì)、PCR擴(kuò)增體系和擴(kuò)增程序參考文獻(xiàn)[11,14]的方法進(jìn)行.成功擴(kuò)增出目的條帶的樣品，送至北京金唯智公司進(jìn)行基因分型.

1.3 數(shù)據(jù)處理

1.3.1 抽樣與微衛(wèi)星DNA遺傳多樣性參數(shù)的計(jì)算 試驗(yàn)所用的5個(gè)種群，均以10群為間隔，設(shè)置15個(gè)梯度樣本量(10、20、30、40、50、60、70、80、90、100、110、120、130、140、150群)，采用簡(jiǎn)單隨機(jī)抽樣法直接抽樣，每個(gè)梯度樣本量抽樣10次.

計(jì)算各種群微衛(wèi)星DNA遺傳多樣性參數(shù).使用Mstool軟件計(jì)算等位基因數(shù)(Na)、期望雜合度(He)、觀察雜合度(Ho)和多態(tài)信息含量(PIC)；使用FSTAT 2.9.3軟件計(jì)算稀疏化后的種群等位基因豐富度(Rs)；基于MsTools軟件得到的數(shù)據(jù)再利用Popgene1.31軟件計(jì)算香農(nóng)指數(shù)(I)和有效等位基因數(shù)(Ne).

1.3.2 統(tǒng)計(jì)分析 使用SPSS 24軟件對(duì)各種群微衛(wèi)星DNA遺傳多樣性參數(shù)與樣本量進(jìn)行相關(guān)性分析，并進(jìn)行回歸分析和曲線擬合，探討樣本量與各參數(shù)的關(guān)系；使用單因素方差分析5個(gè)種群15個(gè)梯度樣本量及其總樣本量的參數(shù)，參數(shù)兩兩之間的顯著性采用多重比較Tukey方法檢驗(yàn).

2 結(jié)果與分析

2.1 樣本量與遺傳變異豐富度參數(shù)的關(guān)系

等位基因數(shù)和等位基因豐富度是反映種群遺傳變異豐富度的參數(shù).5個(gè)種群的等位基因數(shù)和等位基因豐富度均易受樣本量的影響，與樣本量呈極顯著正相關(guān)關(guān)系，相關(guān)系數(shù)(R2)為0.902～0.958.回歸分析和曲線擬合顯示，等位基因數(shù)、等位基因豐富度與樣本量呈顯著對(duì)數(shù)增長(zhǎng)關(guān)系(圖1).等位基因數(shù)和等位基因豐富度在樣本量為10群時(shí)，分別達(dá)到各種群總樣本量的43%～59%和26%～58%；60群時(shí)，分別達(dá)到各種群總樣本量的74%～89%和72%～89%；150群時(shí)，分別達(dá)到各種群總樣本量的91%～100%和91%～99%.

圖1 樣本量與等位基因數(shù)、等位基因豐富度的相互關(guān)系

方差分析顯示，多樣性越高，能反映總體多樣性水平所需的樣本量就越大.多樣性最高的江西種群和廣西種群，樣本量為150群時(shí)的等位基因數(shù)、等位基因豐富度與總樣本量的差異顯著；陜西種群的等位基因豐富度在樣本量為140～150群時(shí)與總樣本量的差異不顯著；重慶種群的等位基因數(shù)在樣本量為130群以上時(shí)與總樣本量的差異不顯著，等位基因數(shù)在樣本量為120群以上時(shí)與總樣本量的差異不顯著；臺(tái)灣種群的等位基因數(shù)和等位基因豐富度在樣本量為150群以上時(shí)與總樣本量的差異不顯著.

2.2 樣本量與雜合度的關(guān)系

期望雜合度和觀察雜合度隨著樣本量的增加，檢測(cè)值會(huì)趨近總樣本量檢測(cè)值，更穩(wěn)定、準(zhǔn)確地反應(yīng)整個(gè)種群的遺傳雜合度(圖2A、2C).

樣本量為10群時(shí)，觀察雜合度占總樣本量的87%～111%，5個(gè)種群標(biāo)準(zhǔn)差均值為0.020 1；樣本量為150群時(shí)，觀察雜合度占總樣本量的99%～102%，標(biāo)準(zhǔn)差均值降至0.002 2(圖2B).標(biāo)準(zhǔn)差方差分析顯示，樣本量為70群以上的觀察雜合度標(biāo)準(zhǔn)差與150群的差異不顯著.

樣本量為10群時(shí)，期望雜合度占總樣本量的93%～109%，5個(gè)種群標(biāo)準(zhǔn)差均值為0.013 7；樣本量為150群時(shí)，期望雜合度占總樣本量的99%～102%，5個(gè)種群標(biāo)準(zhǔn)差均值為0.001 3(圖2D).標(biāo)準(zhǔn)差方差分析顯示，樣本量為40群以上的期望雜合度標(biāo)準(zhǔn)差與150群的差異不顯著.

圖2 樣本量與雜合度、雜合度標(biāo)準(zhǔn)差的相互關(guān)系

2.3 樣本量與有效等位基因數(shù)的關(guān)系

有效等位基因數(shù)與樣本量呈極顯著正相關(guān)關(guān)系，陜西、重慶、江西、廣西4個(gè)種群的有效等位基因數(shù)與樣本量的R2均在0.600 0以上，臺(tái)灣種群有效等位基因數(shù)與樣本量的R2為0.200 0.回歸分析和曲線擬合顯示，有效等位基因數(shù)與樣本量呈對(duì)數(shù)增長(zhǎng)關(guān)系(圖3A).方差分析顯示，陜西種群樣本量為30群以上時(shí)的有效等位基因數(shù)與總樣本量的差異不顯著，重慶、江西種群樣本量為40群以上時(shí)的有效等位基因數(shù)與總樣本量的差異不顯著，廣西種群樣本量為50群以上時(shí)的有效等位基因數(shù)與總樣本量的差異不顯著，臺(tái)灣種群為10群以上的有效等位基因數(shù)與總樣本量的差異不顯著.有效等位基因數(shù)標(biāo)準(zhǔn)差在10～40群時(shí)下降明顯，40～140群時(shí)下降較為平緩(圖3B).標(biāo)準(zhǔn)差方差分析顯示，樣本量為40群以上的有效等位基因數(shù)與150群的差異不顯著.

圖3 樣本量與有效等位基因數(shù)、有效等位基因數(shù)標(biāo)準(zhǔn)差的相互關(guān)系

2.4 樣本量與多態(tài)信息含量的關(guān)系

5個(gè)種群的多態(tài)信息含量與樣本量呈極顯著正相關(guān)關(guān)系，R2為0.482～0.660.回歸分析和曲線擬合顯示，多態(tài)信息含量與樣本量呈對(duì)數(shù)增長(zhǎng)關(guān)系(圖4A).方差分析顯示，多態(tài)信息含量最高的江西種群，樣本量為60群以上時(shí)的多態(tài)信息含量與總樣本量的差異不顯著；其次是陜西、重慶、廣西種群，樣本量為30群以上時(shí)的多態(tài)信息含量與總樣本量的差異不顯著；多態(tài)信息含量最低的臺(tái)灣種群，樣本量為20群以上時(shí)的多態(tài)信息含量與總樣本量的差異不顯著.標(biāo)準(zhǔn)差方差分析顯示，樣本量為80群以上時(shí)的多態(tài)信息含量標(biāo)準(zhǔn)差與150群的差異不顯著(圖4B).

圖4 樣本量與多態(tài)信息含量、多態(tài)信息含量標(biāo)準(zhǔn)差的相互關(guān)系

2.5 樣本量與香農(nóng)指數(shù)的關(guān)系

5個(gè)種群的香農(nóng)指數(shù)與樣本量均呈極顯著正相關(guān)關(guān)系，R2為0.699～0.762.回歸分析和曲線擬合顯示，香農(nóng)指數(shù)與樣本量呈對(duì)數(shù)增長(zhǎng)關(guān)系(圖5A).方差分析顯示，樣本量為60群以上時(shí)的香農(nóng)指數(shù)與總樣本量的差異不顯著，均達(dá)到總樣本量的95%以上.標(biāo)準(zhǔn)差方差分析顯示，樣本量為70～140群的香農(nóng)指數(shù)標(biāo)準(zhǔn)差與150群的差異不顯著(圖5B).

圖5 樣本量與香農(nóng)指數(shù)、香農(nóng)指數(shù)標(biāo)準(zhǔn)差的相互關(guān)系

3 討論

遺傳多樣性研究在使用微衛(wèi)星標(biāo)記等共顯性分子標(biāo)記時(shí)，常用7種遺傳多樣性參數(shù).由于每個(gè)參數(shù)反映的遺傳多樣性信息不同，能反映總體遺傳多樣性所需要的最小樣本量也不同.本次在東方蜜蜂遺傳多樣性的研究中，反映種群遺傳變異豐富度的等位基因數(shù)和等位基因豐富度等參數(shù)檢測(cè)的穩(wěn)定性，均隨著樣本量的增加而增強(qiáng)，樣本量為60群時(shí)能檢測(cè)到70%的等位基因豐富度，樣本量為150群時(shí)能檢測(cè)到91%以上的等位基因豐富度；反映種群遺傳變異均勻度的期望雜合度最小樣本量為40群，觀察雜合度的最小樣本量為70群，有效等位基因數(shù)的最小樣本量為40群，多態(tài)信息含量最小樣本量為60～80群；同時(shí)反映遺傳豐富度和均勻度的香農(nóng)指數(shù)最小樣本量為60～70群.因此，建議在充分保證樣本代表性的前提下，研究蜜蜂遺傳多樣性的樣本量最少60群，最好在80群以上.重點(diǎn)關(guān)注等位基因數(shù)和等位基因豐富度等參數(shù)的研究，樣本量應(yīng)在150群以上；重點(diǎn)關(guān)注期望雜合度和有效等位基因數(shù)的研究，樣本量可在40群以上.

保護(hù)生物學(xué)和保育遺傳學(xué)領(lǐng)域，在沒(méi)有人為遷入的情況下，等位基因的豐富性比等位基因的均勻性更重要，等位基因數(shù)會(huì)決定生物在自然選擇下的生存極限[45].等位基因數(shù)以及稀疏化后的等位基因豐富度的參數(shù)模型就是從豐富度的角度解釋種群遺傳變異程度[46]，樣本量增加使得種群中的更多基因包括稀有基因被檢測(cè)到，因此在變異數(shù)量的檢測(cè)上受樣本量的影響極大.但研究中發(fā)現(xiàn)，這類標(biāo)記除了直接受樣本量的影響外，還可能受到均勻度的間接影響.臺(tái)灣種群的特點(diǎn)是，種群數(shù)量較少，分布在臺(tái)灣全島，遺傳多樣性和遺傳結(jié)構(gòu)受人工育王、蜂場(chǎng)轉(zhuǎn)地移動(dòng)等人為因素的影響較小，在豐富度較低的同時(shí)均勻度也低，并存在大量的稀有等位基因.因此在研究豐富度相關(guān)參數(shù)時(shí)，臺(tái)灣種群的樣本量需要達(dá)到150群時(shí)才能代表總體.雖然這類標(biāo)記在了解種群的等位基因數(shù)，設(shè)計(jì)、實(shí)施方案保護(hù)本土資源時(shí)非常重要，但仍需要與均勻度參數(shù)結(jié)合在一起來(lái)判斷遺傳多樣性水平.

種群遺傳變異均勻度的參數(shù)在樣本量增大時(shí)，偏離總體均勻度的程度會(huì)減小.期望雜合度、觀察雜合度、有效等位基因數(shù)、多態(tài)信息含量都與種群遺傳變異均勻度有關(guān)，是從等位基因比例的角度說(shuō)明種群內(nèi)遺傳變異的分布.雜合度是最常用的反映遺傳多樣性的參數(shù)，其中，期望雜合度可用在單倍體、二倍體和多倍體的研究中，因此被認(rèn)為是遺傳多樣性的最適參數(shù)之一，已有研究認(rèn)為樣本量對(duì)雜合度幾乎沒(méi)有影響[47-48].但本次對(duì)東方蜜蜂5個(gè)種群15個(gè)梯度樣本量抽樣10次的數(shù)據(jù)進(jìn)行分析，發(fā)現(xiàn)樣本量的增加使雜合度檢測(cè)的穩(wěn)定性增強(qiáng)，更能反映總體雜合度水平.樣本量為70群時(shí)的檢測(cè)值波動(dòng)小，能穩(wěn)定反應(yīng)種群整體觀察雜合度水平；樣本量為40群時(shí)就能穩(wěn)定反應(yīng)種群整體期望雜合度水平.有效等位基因數(shù)也是從均勻度角度觀察遺傳多樣性水平，是種群純合體頻率的倒數(shù)，能很好地區(qū)分同樣的等位基因但頻率不同的種群.在種群雜合體較多時(shí)，江西、廣西、陜西、重慶種群的有效等位基因數(shù)均大于2.5，樣本量為40群時(shí)就能代表總體的有效等位基因數(shù)；而臺(tái)灣種群的遺傳多樣性較低，樣本量為10群時(shí)即可代表總體的有效等位數(shù).因奠基者效應(yīng)、小種群、近交等因素導(dǎo)致純合子較多時(shí)，樣本量適當(dāng)減少也能體現(xiàn)總體的有效等位基因數(shù).

香農(nóng)指數(shù)綜合考慮了等位基因數(shù)和等位基因頻率，是較為全面的反映遺傳多樣性指標(biāo).本研究中，香農(nóng)指數(shù)在樣本量為70群時(shí)即可代表總體多樣性水平.香農(nóng)指數(shù)除了受樣本量的影響外，參與研究的標(biāo)記個(gè)數(shù)會(huì)影響結(jié)果.各種群的香農(nóng)指數(shù)需要在同樣的標(biāo)記中才能有效比較.

東方蜜蜂遺傳多樣性研究，除了考慮樣本量外，還需要根據(jù)本物種特點(diǎn)制定合適的取樣方式.其他動(dòng)物或昆蟲(chóng)的研究，采樣策略需要考慮季節(jié)和性比對(duì)遺傳多樣性的影響[17].由于蜜蜂的生殖特點(diǎn)，同一蜂群的遺傳背景是由蜂王決定的，自然狀態(tài)下的蜂群一年以上才會(huì)換王，因此，季節(jié)對(duì)蜜蜂種群遺傳結(jié)構(gòu)的影響不大；性比對(duì)蜜蜂遺傳多樣性的影響更小.蜜蜂是完全社會(huì)性昆蟲(chóng)，雄性蜜蜂只在特定季節(jié)出現(xiàn)，且數(shù)量較少，因此，蜜蜂遺傳多樣性的研究樣本為雌性工蜂.蜜蜂是經(jīng)濟(jì)昆蟲(chóng)，飼養(yǎng)數(shù)量很大，人為因素的干擾會(huì)影響對(duì)當(dāng)?shù)孛鄯溥z傳多樣性評(píng)估的準(zhǔn)確性.合理的采樣策略，需要避免人為因素的干擾，包括引種、人工育王、育種，得到的蜜蜂樣本才能準(zhǔn)確地代表當(dāng)?shù)氐淖匀环N群，可以準(zhǔn)確地評(píng)估種群遺傳多樣性和遺傳結(jié)構(gòu)，以便設(shè)計(jì)后續(xù)的保護(hù)和利用方案.