曾 晶 王子龍 吳小波 顏偉玉 曾志將
(江西農業大學蜜蜂研究所,南昌 330045)
蜂王和工蜂都是由受精卵發育而成,它們的遺傳物質完全相同,但由于幼蟲發育過程中得到的食物和發育空間的不同,發育形成形態、壽命以及其它許多生物學特性迥然不同的蜂王和工蜂,這就是蜜蜂級型分化現象。
蜜蜂級型分化機理,一直是廣大生物學家所關注熱門研究領域。早期的研究者認為:食物數量和成分差異是引起蜜蜂級型分化的主要原因[1]。近幾來的研究表明:DNA甲基化、組蛋白修飾、“royalactin”蛋白質都與蜜蜂級型分化相關[2-4]。
Ku charski等受果蠅實驗啟示,通過檢測dynactin p62甲基化水平的差異來研究西方蜜蜂(Apis mellifera)級型分化,結果發現營養因素可通過影響DNA 甲基化水平調控西方蜜蜂級型分化[2]。Shi等首次發現空間因素也可以通過影響dynactin p62的甲基化來調控西方蜜蜂雌的級型分化[5-6],并且發現:發現飼喂異種蜂王漿可以更顯著降低雌性蜜蜂幼蟲dynactin p62基因甲基化水平[7]。但目前關于蜜蜂DNA甲基化的研究主要集中在西方蜜蜂上,對中華蜜蜂(Apis cerana cerana,簡稱中蜂)的甲基化模式研究很少。正是鑒于此,我們在對中蜂dynactin p62基因克隆及差異表達研究此基礎上[8],開展了發育營養和空間因素對中蜂幼蟲dynactin p62基因表達影響的研究工作,現總結報道如下。
1.1.1 供試蜜蜂
取自江西農業大學蜜蜂研究所飼養的中華蜜蜂。
1.1.2 試劑及儀器
試劑:總RNA提取試劑盒Trizol,RNA酶抑制劑購自全式金公司;DNA marker DL2000,M-MLV反轉錄酶及SYBR GreenⅡ熒光定量試劑購自TaKaRa公司。
儀器:普通離心機(飛鴿 KA-1000 型,上海安亭科學儀器廠公司產品),臺式冷凍離心機(Eppendorf 5810R),移液器(Eppendorf公司產品),普通 PCR儀(Eppendorf Mastercycler personal),MyiQ2 Real-Time PCR System(Bio-Rad公司產品)。
1.2.1 隔王產卵
為了得到基本一致1日齡幼蟲,用隔王框進行隔王產卵。
1.2.2 實驗分組
實驗分為A、B、C三個組,A組是在王臺中,移入1日齡幼蟲,一直飼喂新鮮蜂王漿直至6日齡幼蟲;B組是在王臺中,移入1日齡幼蟲,前3天飼喂新鮮蜂王漿,后3天飼喂大幼蟲食物直至6日齡幼蟲;C組是在工蜂巢房中,移入1日齡幼蟲后,一直飼喂新鮮蜂王漿直至6日齡幼蟲。這樣A組和B組是發育空間一樣,營養因素不同;A組和C組是營養因素一樣,發育空間不同。
1.2.3 人工飼喂幼蟲
吸取0.2ml新鮮蜂王漿(或大幼蟲食物)于王臺或工蜂巢房中,放入培養箱中(溫度35℃,濕度90%)30min,用移蟲針移取1日齡至王臺或工蜂巢房中,幼蟲浮在食物上,每隔12h換一次食物。
1.2.4 樣品采集
每個實驗組取8只發育正常的6日齡幼蟲,每只幼蟲作為1個樣品,用于dynactin p62基因表達量的檢測。取樣后蜜蜂幼蟲樣品迅速放入液氮速凍,-80℃保存,用于提取總RNA。
1.2.5 RNA的提取和cDNA第一鏈的合成
每個樣品用液氮研磨后,用 Trizol 試劑盒進行RNA 的提取。所有操作均按照試劑盒說明書進行,RNA最后溶于30 μL RNA-free的DEPC水中,經瓊脂糖凝膠電泳和紫外分光光度計檢測后,放入(80℃保存。用反轉錄試劑盒對總RNA進行反轉錄,反應體系為50 μL:8 μL總RNA,10 μL Buffer,8 μL dNTP,1.5 μL M-MLV反轉錄酶,3 μL Oligo dT,1 μL RNA酶抑制劑,18.5 μL DEPC水。反轉錄反應條件如下:體系混勻后,42℃反應60 min,75℃ 5 min。反轉錄產物保存于(80℃。
1.2.6 熒光定量PCR
定量引物設計:在中華蜜蜂dynactin p62基因閱讀框區設計多對特異引物用于定量 PCR 反應,根據Livak和Schmittgen的方法,本研究對設計的特異性熒光定量引物進行了篩選,獲得符合要求的特異性引物,用于后續的實時熒光定量 PCR 反應。同時本研究以β-Actin 基因作為內參基因(引物為βact-Q-F和βact-Q-R)。熒光定量所用引物序列(表1)。

表1 中華蜜蜂dynactin p62基因引物
熒光定量 PCR 反應:qRT-PCR反應體系為10 μL:2μL反轉錄產物,目的基因上游和下游引物各0.3 μL,2.4 μL H2O,5 μL SYBR GreenⅡ;PCR擴增程序:94℃預變性3 min;94℃變性30 s,60.9℃和55.8℃(分別對應基因引物和內參引物)退火45 s,72℃延伸40 s,40個循環;72℃延伸10 min;最后以每5 s上升0.5℃的速度從61℃到95℃記錄熔解曲線,每個反轉錄cDNA樣品。反應結束后收集目的基因與內參基因的CT值。
收集目的基因(dynactin p62)與內參基因(β-Actin)的CT值,用qpcR package,R(Hornik 2011)軟件計算擴增效率,計算基因相對表達量(r)的數據分析方法參考Huang等(2012)[9]。
各試驗組基因表達量的差異,采用StatView軟件“ANOVA and t-test”中的“ANOVA or ANCOVA”進行統計分析。
在發育空間因素不變的情況下,雖然A組和B組3日齡后的幼蟲食物營養因素不同,但兩組6日齡幼蟲dynactin p62基因表達量差異不顯著(P>0.05)(圖1)。

圖1 營養因素對中蜂6日齡幼蟲dynactin p62基因表達影響
在發育營養因素不變的情況下,A組和C組幼蟲發育空間因素不同,兩組6日齡幼蟲dynactin p62基因表達量差異顯著(P<0.05)(圖2)。

圖2 空間因素對中蜂6日齡幼蟲dynactin p62基因表達影響
dynactin p62是dynactin蛋白的3個組成部分之一,位于末端,生成的dynactin p62 siRNA能有效干擾外源基因的表達,同時dynactin p62能識別并連接ATP水解酶,并調節其活性。目前蜜蜂的dynactin p62基因的生物功能還不清楚。通過對中蜂dynactin p62基因的克隆和測序可以知道,其基因組DNA序列全長為2 403 bp,cDNA序列全長為1 491 bp,推測蛋白質序列為496個氨基酸,與西方蜜蜂dynactin p62的相似性88%[8]。
dynactin p62基因在不同發育時期中蜂均有表達,剛羽化雌性蜜蜂的表達量顯著高于幼蟲期,且工蜂的表達量顯著高于蜂王;但dynactin p62基因在雄蜂中表達量沒有明顯的規律性。這提示dynactin p62基因可能與蜜蜂級型分化有關[8]。本研究首次研究了發育營養和空間因素對中蜂6日齡幼蟲dynactin p62基因表達影響,發現發育空間因素顯著影響了中蜂6日齡幼蟲dynactin p62基因表達量,且是工蜂巢房中幼蟲表達量高于王臺中幼蟲,這與以前研究結果一致[8]。同時發現營養因素沒有顯著影響中蜂6日齡幼蟲dynactin p62基因表達,這與Weaver 發現的前3.5天食物質量是幼蟲級型分化關鍵時期相吻合[10],但詳細的分子機理還待進一步探討。
[1] 曾志將.蜜蜂生物學[M].北京:中國農業出版社,2007.
[2] Kucharski R,Maleszka J,Foret S,Maleszka R. Nutritional control of reproductive status in honeybees via DNA methylation[J].Science,2008,319:1827-1830.
[3] Spannhoff A,Kim YK,Raynal NJM,Gharibyan V,Su MB,Zhou YY,Li J,Castellano S,Sbardella,G,Issa JPJ. Histone deacetylase inhibitor activity in royal jelly might facilitate caste switching in bees[J]. EMBO reports,2011,12(3):238-243.
[4] Kamakura M. Royalactin induces queen differentiation in honeybees[J]. Nature,2011,473:478- 483.
[5] Shi YY,Huang ZY,Zeng ZJ,Wang ZL,Wu XB,Yan WY.Diet and cell size both affect queen-worker differentiation through DNA Methylation in honey bees(Apis mellifera,Apidae)[J]. PloS ONE,2011,6(4):e18808.
[6] 石元元,曾志將,吳小波,顏偉玉,王子龍. 人工注射 Dnmt3 siRNA對意大利蜜蜂雌蜂發育的影響[J]. 昆蟲學報,2011,54(3):272-278.
[7] 石元元,田柳青,張飛,劉俊峰,顏偉玉,曾志將. 蜂王漿對雌性蜜蜂幼蟲dynactin p62甲基化影響[J]. 中國農業科學,2012,45(23):4909-4915.
[8] 曾晶,潘其忠,王子龍,吳小波,顏偉玉,曾志將. 中華蜜蜂dynactin p62基因的克隆及不同發育階段表達分析[J]. 昆蟲學報,2012,55(10):1132-1141.
[9] Huang Q,Kryger P,Le Conte Y,Moritz RF,. Survival and immune response of drones of a Nosemosis tolerant honey bee strain towards N. ceranae infections[J]. J. Invertebr. Pathol.,2012,109(3):297-302.
[10] Weaver N.Physiology of caste determination[J]. Annu Rev Entomol,1966,11:79-102.