新疆黑蜂基因組選擇信號(hào)分析

何金明，孫智禹，劉玉玲，徐 凱，牛慶生

(吉林省養(yǎng)蜂科學(xué)研究所，吉林 132108)

新疆黑蜂是一個(gè)被列入中國(guó)《中國(guó)畜禽遺傳資源志·蜜蜂志》的本地西方蜜蜂地方品種，中心產(chǎn)區(qū)為阿爾泰山、天山山脈及伊犁河谷地區(qū)，是新疆地區(qū)特有的、與當(dāng)?shù)貧夂颦h(huán)境及蜜源等生態(tài)條件相適應(yīng)的優(yōu)良蜂種，特別適合新疆北疆淺山區(qū)丘陵地帶飼養(yǎng)，具有抗寒、抗逆、抗病、產(chǎn)卵力強(qiáng)、采集力強(qiáng)、節(jié)省飼料及蜂蜜純天然且含天然藥用成分等優(yōu)點(diǎn)[1]。近年來(lái)，隨著飼養(yǎng)方式的改變和引入蜂種的大量應(yīng)用，新疆黑蜂的飼養(yǎng)范圍逐年縮小，飼養(yǎng)量也不斷減少。然而，飼養(yǎng)過(guò)程中的盲目引種導(dǎo)致新疆黑蜂雜化程度非常嚴(yán)重，純種新疆黑蜂數(shù)量急劇下降。2006年，新疆黑蜂被列為《國(guó)家級(jí)畜禽遺傳資源保護(hù)名錄》的瀕危物種之一，對(duì)新疆黑蜂的保護(hù)工作已迫在眉睫[2-3]。

中國(guó)西方蜜蜂(Apismellifera)遺傳資源較為豐富，既有新疆黑蜂、東北黑蜂及浙江漿蜂等地方特色品種，又有意大利蜜蜂、卡尼鄂拉蜂等引入蜂種，不同蜂種在外部形態(tài)、生產(chǎn)性能、抗病力及適應(yīng)性等方面具有較大差異，與各蜂種在進(jìn)化過(guò)程中受到的長(zhǎng)期自然選擇和人工定向選擇密切相關(guān)。選擇作用會(huì)使生物群體在世代傳遞過(guò)程中的某些基因型個(gè)體比例發(fā)生變化，繼而在基因組上留下相應(yīng)的選擇印跡，即選擇信號(hào)(selection signal)[4]。利用基因組上的選擇信號(hào)進(jìn)行群體遺傳學(xué)分析有助于了解種群的進(jìn)化歷史，并對(duì)研究物種的種質(zhì)特性、遺傳育種等問(wèn)題提供幫助。隨著基因組學(xué)和高通量測(cè)序技術(shù)的快速發(fā)展，大范圍的基因組選擇信號(hào)得以快速檢測(cè)，為畜禽種質(zhì)資源的鑒定和保護(hù)利用提供了新的方法[5]。通過(guò)全基因組重測(cè)序技術(shù)檢測(cè)選擇信號(hào)以挖掘動(dòng)物優(yōu)良性狀相關(guān)候選基因已被成功應(yīng)用于牛[6-7]、豬[8]、雞[9]、羊[10-11]等畜禽中，并鑒定出了多個(gè)與繁殖、毛色、生長(zhǎng)等經(jīng)濟(jì)性狀相關(guān)的重要基因。動(dòng)物基因組選擇信號(hào)分析的成功應(yīng)用為揭示蜜蜂遺傳資源進(jìn)化歷史及發(fā)掘其優(yōu)良性狀相關(guān)潛在遺傳學(xué)基礎(chǔ)提供了新的研究方向。當(dāng)前關(guān)于新疆黑蜂的多數(shù)研究?jī)H限于常規(guī)生產(chǎn)性能、繁殖性能、蜂蜜王漿品質(zhì)的測(cè)定與比較，對(duì)其種質(zhì)資源的鑒別僅能根據(jù)形態(tài)特征等傳統(tǒng)方法進(jìn)行判斷，缺乏對(duì)其品種優(yōu)良特性及候選基因的研究，其遺傳進(jìn)化關(guān)系尚不明確[12-13]。鑒于上述情形，本研究利用全基因組重測(cè)序技術(shù)對(duì)新疆黑蜂和其他引入品種進(jìn)行不同蜂種間的遺傳差異檢測(cè)，分析新疆黑蜂的遺傳進(jìn)化關(guān)系，并通過(guò)選擇信號(hào)分析挖掘新疆黑蜂特異性狀的相關(guān)基因，以期為制定更有效的保種措施和進(jìn)一步的蜂種資源利用與開發(fā)提供理論依據(jù)。

1 材料與方法

1.1 材料

供試蜂種為來(lái)源于吉林省養(yǎng)蜂科學(xué)研究所(國(guó)家級(jí)蜜蜂基因庫(kù))的新疆黑蜂(Xinjiang Black bee,XJB)、意大利蜜蜂(Apismelliferaligustica,ITA)、歐洲黑蜂(Apismelliferamellifera,EUR)、卡尼鄂拉蜂(Apismelliferacarnica,CAR)、高加索蜂(Apismelliferacaucasica,CAU)的核心保種群。于2020年5月15日從各蜂種中分別選取10群培育蜂王，2020年6月1日從每群中分別選取1只剛羽化出房的處女蜂王進(jìn)行試驗(yàn)。長(zhǎng)白山中華蜜蜂(Apismelliferacaucasica,A.c.c)采自位于吉林省敦化市的長(zhǎng)白山中華蜜蜂保種場(chǎng)，選取1只未交尾的處女蜂王作為群體遺傳分析的參考群體。收集的蜂王樣本立即投入液氮中，返回實(shí)驗(yàn)室后(約4 h)立即放入-80 ℃保存，用于基因組DNA提取。

1.2 方法

1.2.1 基因組重測(cè)序數(shù)據(jù)獲取 使用天根生化科技(北京)有限公司組織基因組DNA試劑盒提取蜜蜂組織基因組DNA，構(gòu)建500 bp雙末端測(cè)序文庫(kù)，利用Illumina HiSeq PE150平臺(tái)進(jìn)行測(cè)序。

1.2.2 測(cè)序數(shù)據(jù)質(zhì)控、過(guò)濾和處理 為了得到有效且高質(zhì)量測(cè)序數(shù)據(jù)(clean reads)，需要對(duì)Illumina HiSeq PE150平臺(tái)獲得的原始測(cè)序數(shù)據(jù)(raw reads)進(jìn)行嚴(yán)格過(guò)濾。過(guò)濾原則如下：①去除含有接頭序列的reads；②去除單端測(cè)序reads中含N比例>10%的paired reads；③僅保留單端測(cè)序reads中含低質(zhì)量堿基數(shù)<50%的序列。將質(zhì)控后的高質(zhì)量數(shù)據(jù)通過(guò)BWA(參數(shù)：mem-t-k32-M-R)軟件與西方蜜蜂參考基因組(https:∥ftp.ncbi.nlm.nih.gov/genomes/all/GCF/003/254/395/GCF_003254395.2_Amel_HAv3.1/GCF_003254395.2_Amel_HAv3.1_genomic.fna.gz)進(jìn)行比對(duì)，獲得比對(duì)參考基因組后的BAM文件，用SAMTOOLS對(duì)BAM文件進(jìn)行排序、去除重復(fù)，并進(jìn)行SNP和Indel檢測(cè)，獲得高質(zhì)量的遺傳變異位點(diǎn)。使用ANNOVAR軟件對(duì)識(shí)別出的高質(zhì)量遺傳變異進(jìn)行SNP信息注釋。

1.2.3 生物信息學(xué)分析 基于全基因組重測(cè)序技術(shù)檢測(cè)出的高質(zhì)量SNP，采用Mega 6.0在線軟件(http:∥www.megasoftware.net/)的鄰接法(Neighbor-Joining)構(gòu)建系統(tǒng)進(jìn)化樹；使用Plink 1.09在線軟件(https:∥www.cog-genomics.org/plink/1.9/)進(jìn)行主成分分析(principal component analysis,PCA)；利用Admixture 3.0在線軟件(http:∥software.genetics.ucla.edu/admixture/)進(jìn)行群體遺傳結(jié)構(gòu)分析；應(yīng)用Haploview 4.1在線軟件(http:∥www.broad.mit.edu/mpg/haploview/)進(jìn)行連鎖不平衡分析。 通過(guò)Plink 1.09軟件，基于群體遺傳分化指數(shù)(Fst)和核苷酸多樣性比值(θπ)兩種方法篩選排名5%的候選選擇區(qū)域，進(jìn)行選擇信號(hào)的篩選。利用David在線軟件(https:∥david.ncifcrf.gov/)對(duì)提取到的受選擇區(qū)域候選基因進(jìn)行GO功能和KEGG通路富集分析。

2 結(jié) 果

2.1 新疆黑蜂基因組SNP鑒定

新疆黑蜂全基因組重測(cè)序的數(shù)據(jù)經(jīng)質(zhì)控后，最終鑒定出了1 728 216個(gè)SNPs，新疆黑蜂中有73 812個(gè)SNPs位點(diǎn)位于編碼區(qū)，占4.3%；有1 159 436個(gè)SNPs位點(diǎn)位于內(nèi)含子區(qū)，占67%；394 942個(gè)SNPs位點(diǎn)位于基因間區(qū)，占22.85%。新疆黑蜂基因組中存在18 526個(gè)非同義突變位點(diǎn)，可引起氨基酸的改變，可作為該蜂種遺傳特性候選SNP位點(diǎn)。

2.2 系統(tǒng)進(jìn)化樹分析

通過(guò)對(duì)新疆黑蜂、意大利蜜蜂、高加索蜂、歐洲黑蜂、卡尼鄂拉蜂的SNP進(jìn)行系統(tǒng)發(fā)育關(guān)系分析，結(jié)果見圖1。由圖1可知，意大利蜜蜂先與歐洲黑蜂聚在一支，其次是卡尼鄂拉蜂、新疆黑蜂，最后是高加索蜂，與其他蜂種相比，新疆黑蜂與高加索蜂親緣關(guān)系較遠(yuǎn)；新疆黑蜂所有個(gè)體聚為一小支，其群體血統(tǒng)單一，與其他蜂種無(wú)個(gè)體交叉現(xiàn)象。

XJB，新疆黑蜂；ITA，意大利蜜蜂；CAU，高加索蜂；EUR，歐洲黑蜂；CAR，卡尼鄂拉蜂；A.c.c，中華蜜蜂。下同XJB,Xinjiang Black bee;ITA,Apis mellifera ligustica;CAU,Apis mellifera caucasica;EUR,Apis mellifera mellifera;CAR,Apis mellifera carnica;A.c.c,Apis cerana cerana.The same as below圖1 新疆黑蜂及4個(gè)引入蜂種的系統(tǒng)進(jìn)化樹Fig.1 Phylogenetic tree of Xinjiang Black bee and 4 introduced bees

2.3 SNP主成分分析

對(duì)所有試驗(yàn)個(gè)體使用最大的2個(gè)主成分，第一主成分(PC1)和第二主成分(PC2)作圖(圖2)。由圖2可知，第一主成分將高加索蜂(CAU)分開， 第二主成分將其余蜂種分成了新疆黑蜂(XJB)、意大利蜜蜂(ITA)以及歐洲黑蜂(EUR)和卡尼鄂拉蜂(CAR)聚為一簇的群體，新疆黑蜂與其他蜂種之間存在一定的差異。

圖中每一個(gè)點(diǎn)代表一個(gè)樣品，不同的顏色代表不同的蜂種Each point represents a sample,different colors represent different bees圖2 主成分分析結(jié)果Fig.2 Results of principal component analysis

2.4 群體遺傳結(jié)構(gòu)分析

通過(guò)對(duì)5個(gè)蜂種的SNP進(jìn)行群體遺傳結(jié)構(gòu)分析，繪制出了Structure推疊圖(圖3)。由圖3可知，假設(shè)祖先群體亞群數(shù)K=2時(shí)，50個(gè)個(gè)體被劃分為了2個(gè)亞群體(深藍(lán)和淺藍(lán))，高加索蜂、歐洲黑蜂、意大利蜂血統(tǒng)比較純，新疆黑蜂和卡尼鄂拉蜂的個(gè)體有2種組成成分，說(shuō)明其很有可能是由2個(gè)祖先亞群雜交而來(lái)；K=3時(shí)，其他蜂種沒有明顯變化，意大利蜂分離出來(lái)；當(dāng)K>3時(shí)，新疆黑蜂被完全分化出來(lái)，其所屬祖先亞群與意大利蜂、卡尼鄂拉蜂、高加索蜂、歐洲黑蜂完全不同。

2.5 連鎖不平衡分析

由圖4可知，各蜂種的衰減速度是歐洲黑蜂>卡尼鄂拉蜂>高加索蜂>新疆黑蜂>意大利蜜蜂，即意大利蜜蜂的連鎖不平衡衰減速度最小，其次是新疆黑蜂，表明在進(jìn)化過(guò)程中這2個(gè)群體與其他群體相比受馴化程度較高，受選擇強(qiáng)度較大。

圖中分別展示了預(yù)設(shè)為2～6個(gè)群體的結(jié)構(gòu)；每份樣品用一個(gè)垂直線條表示，每個(gè)垂直線條的顏色表示該樣品祖先種群的比例The figure shows the structure of the predetermined 2 to 6 groups,respectively;Each sample is respresented by a vertical line,and line color indicates ancestral population proportion圖3 群體遺傳結(jié)構(gòu)分析結(jié)果Fig.3 Results of population genetic structure analysis

圖4 連鎖不平衡衰減圖Fig.4 Linkage disequilibrium patterns map

2.6 選擇信號(hào)分析

通過(guò)計(jì)算基因組的θπ和Fst值，對(duì)篩選出的前5%區(qū)域的受選擇信號(hào)進(jìn)行分析，將受選擇的基因進(jìn)行GO功能富集分析，結(jié)果顯示，新疆黑蜂受選擇基因主要富集在代謝、繁殖、信號(hào)轉(zhuǎn)導(dǎo)等生物過(guò)程，細(xì)胞組分主要富集在細(xì)胞器、核小體、染色質(zhì)和膜等，分子功能主要富集在有機(jī)環(huán)狀化合物結(jié)合、離子結(jié)合等功能，其中在這些主要富集的通路中共篩選出了細(xì)胞色素P450 314A亞家族多肽1(cytochrome P450 family 314 subfamily A polypeptide 1，Cyp314A1)、毛透明蛋白(Trichohyalin-like)、卷曲螺旋結(jié)構(gòu)域蛋白112(coiled-coil domain-containing protein 112，CCDC112)、山梨醇脫氫酶(Sorbitol dehydrogenase)、絲氨酸蛋白酶抑制劑3(serine protease inhibitor 3，SPI-3)、鐵蛋白(soma ferritin，F(xiàn)er3HCH)、張力蛋白1(tensin-1，TNS1)、磷酸酶與張力蛋白同源物(phosphatase and tensin-like，Pten)、腺苷琥珀酸裂解酶(adenylosuccinate lyase，ADSL)、章魚胺受體β2R(octopamine receptor beta-2R，Octβ2R)、G補(bǔ)綴FHA域血管生成因子1(angiogenic factor with G patch and FHA domains 1，AFFG1)、鳥嘌呤核苷酸交換因子17(rho guanine nucleotide exchange factor 17，ARHGEF17) 12個(gè)候選基因(表1)，這些基因可能與新疆黑蜂的抗寒、產(chǎn)卵力強(qiáng)、抗病等生物學(xué)特性相關(guān)，參與了機(jī)體的繁殖、生長(zhǎng)發(fā)育及免疫等生理活動(dòng)。KEGG通路富集分析結(jié)果表明，卡尼鄂拉蜂與新疆黑蜂基因組相比，新疆黑蜂受選擇基因主要富集于嘌呤代謝、胞吞作用等信號(hào)通路；高加索蜂與新疆黑蜂基因組相比，新疆黑蜂受選擇基因主要富集于胞吞作用、磷酸肌醇代謝等信號(hào)通路；歐洲黑蜂與新疆黑蜂基因組相比，新疆黑蜂受選擇基因主要富集于物質(zhì)代謝等信號(hào)通路，意大利蜂與新疆黑蜂基因組相比，新疆黑蜂受選擇基因主要富集于嘌呤代謝等信號(hào)通路(圖5)。

表1 新疆黑蜂受選擇基因GO功能富集分析

A～D，分別為與卡尼鄂拉蜂、高加索蜂、歐洲黑蜂、意大利蜂相比新疆黑蜂受選擇的KEGG通路富集分析圖A-D,KEGG pathway enrichment analysis of selected genes in Xinjiang Black bees compared with Apis mellifera carnica,Apis mellifera caucasica,Apis mellifera mellifera,Apis mellifera ligustica圖5 新疆黑蜂受選擇基因KEGG通路富集分析Fig.5 KEGG pathway enrichment analysis of selected genes in Xinjiang Black bee

3 討 論

據(jù)報(bào)道，俄國(guó)人于1900年將黑色蜂種帶入了中國(guó)新疆伊犁和阿勒泰地區(qū)，1919年又經(jīng)哈薩克斯坦將黑蜂帶入了新疆布爾津縣，1925-1926年又將黑蜂帶到了新疆伊寧飼養(yǎng)，慢慢發(fā)展到了整個(gè)伊犁地區(qū)，經(jīng)過(guò)百余年的長(zhǎng)期混養(yǎng)、自然雜交和人為選育逐漸形成了適應(yīng)新疆氣候環(huán)境的蜂蜜高產(chǎn)型地方品種[1]。但是關(guān)于該蜂種是何起源未見詳細(xì)記錄。一直以來(lái)，新疆黑蜂被認(rèn)為是歐洲黑蜂或高加索蜂的后代[14-15]，但都沒有相關(guān)的遺傳學(xué)證據(jù)。本研究對(duì)國(guó)家級(jí)蜜蜂基因庫(kù)(吉林)異地保存的新疆黑蜂保種群體進(jìn)行遺傳學(xué)分析，發(fā)現(xiàn)新疆黑蜂與卡尼鄂拉蜂、意大利蜂、歐洲黑蜂親緣關(guān)系較近，而與高加索蜂的親緣關(guān)系較遠(yuǎn)，所以暫不支持新疆黑蜂來(lái)源于高加索蜂。西方蜜蜂的起源進(jìn)化和分類問(wèn)題一直存在爭(zhēng)議，在目前已被確認(rèn)的29個(gè)西方蜜蜂亞種中，根據(jù)形態(tài)學(xué)及線粒體分子數(shù)據(jù)被劃分為4個(gè)主要譜系：A系(非洲)、M系(西歐)、C系(東歐)、O系(中東)[16]。本研究系統(tǒng)進(jìn)化樹結(jié)果與上述分類略有不同，可能是由于分析所用的數(shù)據(jù)集和具體的分析方法不同。此外，結(jié)合系統(tǒng)發(fā)育分析、主成分分析及群體遺傳結(jié)構(gòu)和連鎖不平衡分析結(jié)果表明，新疆黑蜂在經(jīng)歷了長(zhǎng)時(shí)間的人工選擇及對(duì)中國(guó)新疆地區(qū)的氣候環(huán)境產(chǎn)生高度適應(yīng)性后，與其他引入蜂種已經(jīng)產(chǎn)生了明顯的遺傳差異，形成了獨(dú)立的遺傳背景，且在進(jìn)化過(guò)程中受到的選擇強(qiáng)度較大，推測(cè)可能是由于長(zhǎng)期的自然選擇和人工選育造成的。

新疆黑蜂對(duì)中國(guó)新疆地區(qū)的氣候、蜜源等生態(tài)條件逐漸形成了很強(qiáng)的適應(yīng)性。本研究篩選出了毛透明蛋白(Trichohyalin-like)、山梨醇脫氫酶(Sorbitol dehydrogenase)等新疆黑蜂受選擇的候選基因,這些基因已被證明與昆蟲的抗寒及能量代謝有關(guān)。毛透明蛋白是一種富含谷氨酸、丙氨酸、精氨酸、賴氨酸、絲氨酸的結(jié)構(gòu)蛋白，其特征氨基酸也是抗寒物質(zhì)中的關(guān)鍵氨基酸。研究顯示，不同溫度梯度下的煙粉虱MED隱種和Asia Ⅱ 3隱種Trichohyalin基因的表達(dá)存在差異，且該基因在抵御低溫過(guò)程中發(fā)揮了一定的調(diào)控作用[17]，有助于煙粉虱在低溫條件下生存；山梨醇脫氫酶對(duì)昆蟲滯育期能量代謝的調(diào)控起到重要作用，它可以將越冬期的抗寒物質(zhì)山梨醇轉(zhuǎn)化成糖原，進(jìn)而分解成可供昆蟲生長(zhǎng)發(fā)育的能源物質(zhì)[18]。新疆黑蜂可以在惡劣的氣候、地理?xiàng)l件下生存，是中國(guó)唯一能夠在野外生存的西方蜜蜂類型，而這些候選基因的功能將有助于新疆黑蜂抗寒、越冬性能強(qiáng)等生物學(xué)特性的形成。

本研究篩選出的CYP314A1、CCDC112、Octβ2R等新疆黑蜂受選擇的候選基因在昆蟲以及人類繁殖、產(chǎn)卵等生命活動(dòng)中發(fā)揮著重要作用。研究發(fā)現(xiàn)，CYP314A1基因?qū)儆贑YP314家族成員，與昆蟲中誘導(dǎo)卵子發(fā)生和抗藥性形成等生理活動(dòng)密切相關(guān)[19-20]；CCDC112是與高齡婦女產(chǎn)卵有關(guān)的差異蛋白之一[21]；Octβ2R基因可以激活多種信號(hào)通路以誘導(dǎo)產(chǎn)卵，而Octβ2R基因純合突變則導(dǎo)致果蠅雌性個(gè)體不育，它們雖表現(xiàn)出正常的求愛、交配、精子儲(chǔ)存和交配后排斥行為，但最后無(wú)法產(chǎn)卵[22]。新疆黑蜂體內(nèi)CYP314A1、CCDC112、Octβ2R等繁殖候選基因在進(jìn)化過(guò)程中受到選擇可能有利于該蜂種高繁殖能力的形成，因而具有產(chǎn)卵力較強(qiáng)、子脾面積大且春季繁殖快等特點(diǎn)，這些基因可以作為日后研究新疆黑蜂優(yōu)異繁殖性能的候選基因。此外，本研究中還篩選出了一些與生長(zhǎng)發(fā)育相關(guān)的基因，它們也在新疆黑蜂進(jìn)化的過(guò)程中受到了選擇，如與細(xì)胞凋亡、增殖及血管生成有關(guān)的Kruppel樣因子6(Kruppel-like factor 6，KLF6)[23]；通過(guò)誘導(dǎo)細(xì)胞凋亡來(lái)參與抵抗家蠶核型多角體病毒侵染的細(xì)胞色素C(Cytc)[24]；在細(xì)胞黏附、運(yùn)動(dòng)、生長(zhǎng)等細(xì)胞生物學(xué)活動(dòng)中發(fā)揮著重要作用的TNS1基因[25]；協(xié)作參與維持基因組穩(wěn)定性、修復(fù)DNA、參與控制基因重組的Pten基因[26]；與肌肉、生長(zhǎng)、維持細(xì)胞代謝有關(guān)的ADSL基因[27]；高表達(dá)于血管內(nèi)皮細(xì)胞，編碼蛋白具有強(qiáng)力血管生成作用的AFFG1以及對(duì)維持內(nèi)皮間連接起著非常重要作用的ARHGEF17等基因[28-29]。KLF6、TNS1等上述基因的表達(dá)調(diào)控可能在新疆黑蜂的生長(zhǎng)發(fā)育過(guò)程中起到了重要作用，使其更能應(yīng)對(duì)伊犁北疆寒冷多風(fēng)的惡劣生存環(huán)境。

新疆黑蜂在抗病力方面也優(yōu)于其他西方蜜蜂品種，本研究篩選出了SPI-3、Fer3HCH等多個(gè)與昆蟲免疫相關(guān)的基因。已有研究發(fā)現(xiàn)，絲氨酸蛋白酶抑制劑(Serpins)是昆蟲體內(nèi)一類非常重要的蛋白酶抑制劑，昆蟲免疫應(yīng)答的重要組成部分，其家族成員SPI-1和SPI-2基因在正常大鼠肝臟中表達(dá)，但SPI-3基因在正常大鼠中幾乎處于沉默，卻在大鼠肝臟急性炎癥期間以及大鼠腦短暫缺血后變得短暫活躍，表明SPI-3基因可能參與了針對(duì)炎癥應(yīng)激的細(xì)胞防御反應(yīng)[30-31]。昆蟲中的Fer3HCH基因通常存在于液泡系統(tǒng)和血淋巴中，研究認(rèn)為該基因是赤擬谷盜的一種潛在的抗菌效應(yīng)因子，在抵御病原感染方面具有一定作用；另外，果蠅體內(nèi)基因過(guò)表達(dá)試驗(yàn)結(jié)果表明，F(xiàn)er3HCH基因在保護(hù)機(jī)體免受氧化應(yīng)激脅迫等方面發(fā)揮著重要作用[32]。以上與免疫相關(guān)的候選基因在新疆黑蜂體內(nèi)受選擇，有利于該蜂種形成免疫力強(qiáng)的體質(zhì)，從而具有較強(qiáng)的抗病力。

4 結(jié) 論

本研究通過(guò)全基因組重測(cè)序技術(shù)在新疆黑蜂中共鑒定出1 728 216個(gè)SNPs，其中18 526個(gè)非同義突變位點(diǎn)是后續(xù)研究新疆黑蜂遺傳特性候選的SNPs位點(diǎn)；篩選到了CYP314A1、Trichohyalin-like、CCDC112、Sorbitol dehydrogenase、SPI-3、Fer3HCH、TNS1、Pten、ADSL、Octβ2R、AFFG1、ARHGEF17共12個(gè)候選基因，解析了新疆黑蜂的抗寒、產(chǎn)卵力強(qiáng)、抗病等優(yōu)異性狀形成的分子機(jī)理，為新疆黑蜂種質(zhì)資源評(píng)價(jià)、保護(hù)及開展分子標(biāo)記輔助選擇提供了理論依據(jù)。