短管赤眼蜂Argonaute蛋白基因家族鑒定及表達分析

楊志強, 王麗艷, 梁洪宇, 劉 旭, 白雪萍, 趙長江, 張海燕

(黑龍江八一農墾大學農學院, 黑龍江省普通高等學校寒地作物種質改良與栽培重點實驗室, 大慶 163319)

短管赤眼蜂Argonaute蛋白基因家族鑒定及表達分析

楊志強#, 王麗艷#, 梁洪宇, 劉 旭, 白雪萍, 趙長江, 張海燕*

(黑龍江八一農墾大學農學院, 黑龍江省普通高等學校寒地作物種質改良與栽培重點實驗室, 大慶 163319)

動物和植物Argonaute蛋白在RNA誘導的基因沉默中發揮重要作用。本研究采用生物信息學方法在全基因組水平鑒定短管赤眼蜂Argonaute蛋白基因家族,鑒定并命名6個TpAGO蛋白。蛋白結構及系統進化分析將TpAGO家族蛋白分為PIWI和Ago兩個亞家族。TpAGO蛋白二級結構的數量比例差異不大，分布卻存在一定差異,其中α-螺旋分布差異尤為顯著。TpAGO蛋白基因結構表現明顯的多樣性,基因外顯子數量差異明顯,從1到20個數量不等;Ka/Ks分析表明TpAGO2和TpAGO5,TpAGO3和TpAGO4兩對基因間進化動力為正選擇作用,其他組合間表現為純化選擇。EST表達分析發現TpAGO基因在昆蟲的卵期、幼蟲期、蛹期和成蟲期均有不同程度的表達,而且靶向EST在同一物種中具有單一性別(雌性或雄性)表達特性,表明短管赤眼蜂TpAGO蛋白基因可參與赤眼蜂性別及發育調控。該研究不僅有助于揭示赤眼蜂產雌孤雌生殖的表觀遺傳機理,也有助于推動赤眼蜂等天敵的生防應用。

短管赤眼蜂; Argonaute; 基因家族; 生物信息學

害蟲生物防治是有害生物治理中最成功、最節約和環境安全的方法。理論研究的核心是明確各類天敵的特性及其對害蟲的控制功能,以及天敵和害蟲之間相互作用規律[1]。赤眼蜂屬膜翅目,是一種卵寄生性天敵昆蟲,目前已成為害蟲生物防治中應用最多,范圍最廣的害蟲天敵之一。其中,感染共生菌Wolbachia的赤眼蜂具有更好的應用前景,因其營產雌孤雌生殖方式致使赤眼蜂后代完全為雌性,同比大大提高了天敵赤眼蜂的寄生能力,即有效增強了天敵對害蟲的控制能力,提高了生物防治效果。盡管已有大量研究表明,赤眼蜂的產雌孤雌生殖方式是由Wolbachia感染所決定的[2-4],但是Wolbachia與赤眼蜂互作調控產雌孤雌生殖的分子機理至今未見明確報道。

MicroRNAs(miRNA)是內源非編碼小RNA,在細胞生物學過程中發揮重要的調控作用,涉及發育、分化、凋亡和免疫等過程,值得一提的是,miRNA也參與Wolbachia與宿主互作調控[5]。那么,小RNA沉默復合體(RNA-induced silencing complex,RISC)的重要組成部分Argonaute蛋白(Argonaute protein)在Wolbachia與宿主互作中的功能如何？Argonaute 蛋白是一個高度保守的家族,含有PIWI和PAZ兩個共同的結構域,存在于絕大部分真核生物中。Argonaute蛋白在各個物種間的數量差距很大,其中人類Homosapiens含有8種[6],擬南芥Arabidopsisthaliana中有10種[7]、果蠅Drosophilamelanogaster中有 5種、家蠶Bombyxmori4種[8],秀麗隱桿線蟲Caenorhabditiselegans27種[9]。但是,迄今為止,關于赤眼蜂Argonaute蛋白及其功能研究尚未見報道。

本研究利用生物信息學方法,在短管赤眼蜂Trichogrammapretiosum基因組水平分離鑒定Argonaute蛋白基因,在對基因結構、蛋白結構和系統進化分析的同時,進行了基于EST的表達譜分析,探討Argonaute蛋白基因的功能。該研究有助于揭示Argonaute蛋白及其協調小RNA途徑在赤眼蜂生長發育和性別調控中的作用,為在表觀遺傳學水平解析赤眼蜂產雌孤雌生殖機理奠定理論基礎。

1 材料與方法

1.1 基因家族成員鑒定

利用已發表的家蠶、果蠅Argonaute蛋白序列為種子序列,通過BLASTp在NCBI數據庫(http:∥www.ncbi.nlm.nih.gov/)中搜索短管赤眼蜂基因組(JARR00000000.1);同時,以蛋白保守結構域PIWI和PAZ為種子在NCBI數據庫搜索短管赤眼蜂基因組,分離鑒定短管赤眼蜂Argonaute蛋白。然后,使用BioEdit軟件構建本地數據庫,對上述兩種方法獲得的短管赤眼蜂Argonaute蛋白進行篩選去冗余并命名TpAGO。采用Expasy(http:∥web.expasy.org/protparam/)分析蛋白分子式、分子量和等電點等數據;采用WoLF PSORT(http:∥www.genscript.com/tools/wolf-psort)對蛋白進行亞細胞定位預測。

1.2 氨基酸比對及系統進化樹的構建

使用DNAMAN對短管赤眼蜂TpAGO蛋白進行比對作圖;使用MEGA7.0對不同昆蟲AGO蛋白采用鄰接法比對并生成系統進化樹,進行1 000次Bootstrap檢驗。

1.3 蛋白高級結構分析

利用SOPMA 在線分析軟件(https:∥npsa-prabi.ibcp.fr/cgi-bin/npsa_automat.pl?page=npsa_sopma.html)分析蛋白二級結構;使用Expasy(https:∥swissmodel.expasy.org/interactive)分析蛋白三級結構。

1.4 基因結構分析

短管赤眼蜂TpAGO全長基因提交GSDS軟件(http:∥gsds.cbi.pku.edu.cn/)繪制基因結構圖。

1.5 堿基替換率分析

將短管赤眼蜂的TpAGO基因的編碼區(去終止子)通過MEGA7.0的ClustalW比對,然后使用DNAsp軟件計算TpAGO基因的非同義替換率Ka、同義替換率Ks以及他們的比值(Ka/Ks)。

1.6 基因表達分析

使用短管赤眼蜂TpAGO基因編碼序列在NCBI的EST數據庫進行BLASTn分析;同時使用TpAGO氨基酸序列在NCBI的EST數據庫進行tBLASTn分析,獲得TpAGO基因的EST表達數據。

2 結果與分析

2.1 短管赤眼蜂TpAGO蛋白基因的鑒定

通過兩種方法共獲得191條序列,一方面采用DNAMAN進行序列分析,去除冗余序列,另一方面利用NCBI數據庫對候選蛋白進行結構域分析,最終獲得6個非冗余的短管赤眼蜂AGO蛋白,分別命名為TpAGO1至TpAGO6(表1)。TpAGO蛋白編碼基因長度介于2 364 bp(TpAGO2.1)至3 144 bp(TpAGO4)之間。其中,4個TpAGO蛋白預測為細胞質或細胞核定位,包括TpAGO1、TpAGO2、TpAGO3和TpAGO4;2個TpAGO蛋白預測為線粒體或細胞核定位,分別為TpAGO5和TpAGO6。

表1短管赤眼蜂TpAGO蛋白基因特征1)

Table1CharactersofTpAGOinTrichogrammapretiosum

蛋白名稱ProteinnameNCBI登錄號Proteinregistrationno．開放閱讀框/bpOpenreadingframe氨基酸Aminoacids分子式Molecularformula分子量/kDMolecularweight亞細胞定位SubcellularlocalizationTpAGO1XP＿014222947．12748915C4534H7194N1298O1357S36102．76Nuc/CytTpAGO2XP＿014229707．12364787C4032H6347N1135O1163S3590．45Cyt/NucTpAGO3XP＿014222980．12727908C4553H7349N1275O1346S38102．71Cyt/NucTpAGO4XP＿014235216．131441047C5198H8245N1523O1550S38118．09Nuc/CytTpAGO5XP＿014225266．12691896C4478H7078N1272O1282S49100．82Mit/NucTpAGO6XP＿014228304．12856951C4720H7451N1343O1377S47106．55Mit/Nuc

1) Nuc:細胞核; Cyt:細胞質; Mit:線粒體。

Nuc:cell nucleus, Cyt:cytoplasm, Mit:mitochondria.

2.2 短管赤眼蜂AGO蛋白結構與系統進化分析

使用DNAMAN對鑒定獲得的6種TpAGO蛋白進行結構域分析,均具有ArgoN、PAZ和PIWI保守結構域,氨基酸組成變異較大(圖1)。通過鄰接法對鑒定的6個短管赤眼蜂、4個家蠶、5個黑腹果蠅、3個麗蠅蛹集金小蜂Nasoniavitripennis的Argonaute蛋白氨基酸序列進行比對構建系統進化樹(圖2)。通過系統進化分析發現,供試的18個Argonaute蛋白可分為兩大亞家族,分別為Ago亞家族和PIWI亞家族。PIWI亞家族蛋白涵蓋4個物種的10個成員,含有典型的PAZ、PIWI和ArgoN結構域,包括短管赤眼蜂TpAGO1、TpAGO2和TpAGO3,果蠅DmPIWI、DmAUB和DmAGO3,家蠶BmPIWI、BmAGO3,麗蠅蛹集金小蜂NvAUB和NvAGO3。其中,短管赤眼蜂TpAGO1和TpAGO2與NvAUB、DmPIWI和DmAUB、BmPIWI親緣關系較近,TpAGO3與NvAGO3、DmAGO3和BmAGO3親緣關系較近。Ago亞家族蛋白涵蓋4個物種的8個成員,含有典型的PAZ、PIWI、ArgoN和Argol1結構域(圖1),包括短管赤眼蜂TpAGO4、TpAGO5和TpAGO6,果蠅DmAGO1和DmAGO2,家蠶BmAGO1和BmAGO2,麗蠅蛹集金小蜂。其中,短管赤眼蜂TpAGO5和TpAGO6與DmAGO1和BmAGO1親緣關系較近;TpAGO4與NvAGO2、DmAGO2和BmAGO2親緣關系較近。

圖1 短管赤眼蜂TpAGO蛋白序列同源性及保守功能區分析Fig.1 Homology and conserved domains of TpAGO proteins

圖2 昆蟲Argonaute蛋白系統進化樹Fig.2 Phylogenetic tree of insect Argonaute proteins

2.3 短管赤眼蜂TpAGO基因結構與進化動力分析

對6個短管赤眼蜂TpAGO基因結構進行分析(圖3)發現,AGO蛋白基因結構表現明顯的多樣性,基因外顯子數量差異明顯,從1到20個數量不等,與蛋白氨基酸結構和亞家族分類無明顯對應關系。其中,TpAGO1含有10個外顯子,TpAGO2有8個外顯子,TpAGO3含有11個外顯子,TpAGO4含有20個外顯子,TpAGO6只有1個外顯子。值得注意的是TpAGO5含有12個外顯子,且內含子跨度較大,最大內含子跨度為32 623 bp,而且內含子中間有許多未知序列,該基因結構或基因組序列需進一步確證。此外6個TpAGO基因終止密碼子也不完全相同,除TpAGO6終止密碼子為tag外,其余5個TpAGO基因的終止子均為taa。

另外,通過對短管赤眼蜂TpAGO基因編碼區序列堿基替代率分析(表2)發現,Ka/Ks比值高于1的組合有TpAGO2和TpAGO5,TpAGO3和TpAGO4,表明這兩對基因間(或者兩對基因所在兩個亞家族間)進化動力為正選擇作用。其他組合的Ka/Ks比值均小于1,表現為純化選擇。但是AGO基因家族成員間Ka/Ks比值整體較高,最低的為0.732 1,表明該家族基因可能經受較大的選擇壓力,可能與基因結構的多樣性有一定關聯。

表2 短管赤眼蜂TpAGO基因的Ka/Ks比值

2.4 短管赤眼蜂TpAGO蛋白高級結構分析

通過軟件對6個TpAGO蛋白的二級和三級結構進行分析(表3,圖4,5),4種二級結構α-螺旋、延伸鏈、β-折疊和無規則卷曲在AGO蛋白間數量比例差異不大,分布卻存在一定差異(圖4),但是AGO蛋白的三級結構基本相似(圖5),相差不大。供試的6個AGO蛋白中α-螺旋氨基酸數量比例在28.37%～40.64%之間;延伸鏈的氨基酸數量比例在17.86%～25.16%之間;β-折疊的氨基酸數量比例在7.59%～10.29%之間;無規則卷曲的氨基酸數量比例在30.07%～44.70%之間。而且,均以α-螺旋和無規則卷曲兩種二級結構為主。其中,α-螺旋分布在AGO蛋白間存在一定差異(圖4),PIWI亞家族蛋白(TpAGO1、TpAGO2和TpAGO3)的N端α-螺旋分布明顯高于Ago亞家族蛋白(TpAGO4、TpAGO5和TpAGO6)。上述結構分析結果與蛋白進化存在明顯的相關性,表明兩個亞家族成員在功能進化上存在一定分歧。

表3短管赤眼蜂TpAGO蛋白二級結構

Table3SecondarystructureoftheTpAGOproteins

二級結構氨基酸數∥百分比/% Number∥PercentageofAAα?螺旋Alphahelix延伸鏈Extendedstrandβ?折疊Betaturn無規則卷曲RandomcoilTpAGO6296∥31．13195∥20．5076∥7．99384∥40．38TpAGO5276∥30．80196∥21．8868∥7．59356∥39．73TpAGO4297∥28．37187∥17．8695∥9．07468∥44．70TpAGO3369∥40．64181∥19．9385∥9．36273∥30．07TpAGO2264∥33．55198∥25．1681∥10．29244∥31TpAGO1323∥35．30192∥20．9876∥8．31324∥35．41

圖4 短管赤眼蜂TpAGO蛋白二級結構Fig.4 Secondary structure of the TpAGO proteins

2.5 短管赤眼蜂TpAGO基因在發育階段的表達

通過TpAGO基因在NCBI的EST數據庫比對分析,EST靶向序列以BLASTn比對中E值小于2e-13或tBLASTn比對中E值小于1e-67為標準,獲得短管赤眼蜂TpAGO基因表達信息(表4),結果表明TpAGO基因在昆蟲的卵期、幼蟲期、蛹期和成蟲期均有不同程度的表達,而且靶向EST在同一物種中具有單一性別(雌性或雄性)表達特性。其中TpAGO1基因靶向EST來源包括麗蠅蛹集金小蜂雄性蛹,埃及伊蚊Aedesaegypti雌成蟲,擬果蠅雌成蟲,褐飛虱Nilaparvatalugens幼蟲和胡蜂Nasoniagiraulti幼蟲。TpAGO2基因靶向EST來源包括擬果蠅雌成蟲,埃及伊蚊雌成蟲,麗蠅蛹集金小蜂蛹和成蟲,胡峰蛹和成蟲、埃及伊蚊幼蟲,褐飛虱幼蟲和玉帶鳳蝶Papiliopolytes蛹。TpAGO3基因靶向EST來源包括麗蠅蛹集金小蜂雌、雄蛹,埃及伊蚊雌成蟲,茶足柄瘤蚜繭蜂Lysiphlebustestaceipes成蟲。TpAGO4基因靶向EST來源包括麗蠅蛹集金小蜂雄蛹,麗蠅蛹集金小蜂雌成蟲,赤擬谷盜Triboliumcastaneum雌性卵,埃及伊蚊雌成蟲。TpAGO5基因靶向EST來源包括雄麗蠅蛹集金小蜂雄性蛹和成蟲,豌豆長管蚜Acyrthosiphonpisum雄性若蟲,桃蚜Myzuspersicae雌成蟲,黑腹果蠅和螺旋蠅Cochliomyiahominivorax卵,二斑葉螨Tetranychusurticae卵,食糞金龜Onthophagustaurus蛹,刺舌蠅Glossinamorsitans成蟲。TpAGO6基因靶向EST來源包括麗蠅蛹集金小蜂雄性蛹和成蟲,家蠶雌性蛹,桃蚜雌成蟲,果蠅卵,豌豆長管蚜若蟲,食糞金龜蛹,玉帶鳳蝶成蟲、三角渦蟲Dugesiajaponica、赤擬谷盜、刺舌蠅和岡比亞按蚊Anophelesgambiae成蟲。

表4短管赤眼蜂TpAGO基因EST表達譜1)

Table4ExpressionprofilesofTpAGOgenesbyESTdata

續表4Table4(Continued)

性別 Sex雄性Male雌性Female混合Mix發育階段 Developmentalstage卵期Egg幼蟲期Larva蛹期Pupa成蟲AdultTpAGO6N．vGE426508．1，N．vGE354623．1M．pEE264736．1，P．pGW819200．1，B．mFY018570．1G．mDV618432．1，D．mAI455659．2，D．mAA390688．1，D．mBP544019．1D．mAI455659．2，D．mAA390688．1，D．mBP544019．1A．pCV842652．1N．vGE354623．1，O．tFG542481．1，B．mFY018570．1N．vGE426508．1，D．jFY948757．1，T．cDT777867．1，G．mDV618432．1，M．pEE264736．1，P．pGW819200．1，A．gBM584762．1

1) 斜體表示使用BLASTn對EST數據庫進行比對,其余為tBLASTn對EST數據庫比對。N.v:麗蠅蛹集金小蜂; N.l:褐飛虱; N.g:胡蜂; A.a:埃及伊蚊; A.p:豌豆長管蚜;A.g:岡比亞按蚊;D.s:擬果蠅; D.j:三角渦蟲; D.m:黑腹果蠅; L.t:茶足柄瘤蚜繭蜂; P.p:玉帶鳳蝶; C.h:螺旋蠅; M.p:桃蚜; T.c:赤擬谷盜; T.u:二斑葉螨;G.m:刺舌蠅; O.t:食糞金龜; B.m:家蠶。

Genes in italic are analyzed by BLASTn in the EST database, and the remaining by tBLASTn. N.v:Nasoniavitripennis; N.l:Nilaparvatalugens; N.g:Nasoniagiraulti; A.a:Aedesaegypti; A.p:Acyrthosiphonpisum; A.g:Anophelesgambiae; D.s:Drosophilasimulans; D.j:Dugesiajaponica; D.m:Drosophilamelanogaster; L.t:Lysiphlebustestaceipes; P.p:Papiliopolytes; C.h:Cochliomyiahominivorax; M.p:Myzuspersicae; T.c:Triboliumcastaneum; T.u:Tetranychusurticae; G.m:Glossinamorsitansmorsitans; O.t:Onthophagustaurus; B.m:Bombyxmori.

3 結論與討論

昆蟲中小分子RNA可分為3種,siRNA(small interfering RNA)、miRNA(microRNA)及piRNA(PIWI-interacting RNA)。在siRNA通路中,RNA誘導沉默復合體(RNA-induced silencing complex,RISC)中的AGO-2與siRNA雙鏈體(siRNA duplex)的引導鏈(guide strand)結合,同時切割siRNA雙鏈體的另一條過客鏈(passenger stand)使其脫離復合體;然后AGO-2蛋白通過siRNA引導鏈連接序列互補的靶mRNA,促成靶基因的沉默。在miRNA通路中,AGO-1蛋白與miRNA(引導鏈)/miRNA*(過客鏈)結合,分離miRNA*[10];然后,與靶基因互補的miRNA能夠促使AGO-1對靶標mRNA切割[11]。在piRNA通路中,PIWI亞家族成員參與調控[12-13],piRNA中多為逆轉錄轉座子和轉座子的反義鏈,PIWI和Aub蛋白與它們結合,引導并切割正義的逆轉錄轉座子轉錄產生本來的正義piRNA[14],AGO-3結合正義piRNA,并切割來源于不同逆轉錄轉座子分化的反義轉錄本[15]。

Argonaute 蛋白家族是一個高度保守的家族,該蛋白以基因家族的形式存在于絕大部分真核生物中。果蠅含有5 個特征性的Argonaute 蛋白:DmPiwi、DmAubergine、DmAgo-l、DmAgo-2 和DmAgo-3,其中DmPIWI、DmAubergine 和DmAgo-3 都屬于PIWI 亞家族[16]。而DmPIWI和DmAub蛋白的突變會影響果蠅胚胎發育。PIWI亞家族蛋白在果蠅中影響轉基因調節的共抑制作用[17],參與生殖干細胞的發育調節[18]。DmAub蛋白可以通過調節oskar和gurken的翻譯及極細胞的形成來影響胚胎發育,DmAub蛋白也可以通過睪丸中的stellate調節抑制來體現它的必要性[19]。DmAGO1和DmAGO2基因表達貫穿果蠅整個胚胎發育時期,其中DmAGO1對于胚胎的發育系統有正調節作用。而且,DmAGO1與DmAGO2也可以在生殖系統外表達,具有更多的功能。另外,家蠶中的Argonaute 家族已經被鑒定出來,包括BmAgo-1、BmAgo-2、BmAgo-3 和BmPIWI,且在數量和序列上都與其他昆蟲是保守的[7]。前兩者屬于Argonaute 亞家族,分別在siRNA/miRNA 通路中發揮作用[20-21],后兩者屬于PIWI亞家族,在家蠶的生殖器官中表達量較高[20]。上述研究表明,Argonaute蛋白極有可能參與赤眼蜂生長發育和性別調控,也有可能參與Wolbachia感染對宿主赤眼蜂產雌孤雌生殖的分子調控。

本研究首次在短管赤眼蜂全基因組水平分離Argonaute蛋白,鑒定的6個TpAGO蛋白基因都含有典型的PAZ和PIWI結構域,與果蠅和家蠶的Argonaute蛋白基因結構基本一致。短管赤眼蜂TpAGO1和TpAGO2與NvAUB,DmPIWI和DmAUB,BmPIWI親緣關系較近,推測赤眼蜂TpAGO1-2可能與果蠅和家蠶AGO蛋白功能相近,參與短管赤眼蜂生殖干細胞的形成,以及在赤眼蜂piRNA通路中起重要作用。TpAGO3與NvAGO3、DmAGO3和BmAGO3親緣關系較近,推測赤眼蜂TpAGO3可能與家蠶AGO3蛋白功能相近,在幼蟲期到成蟲期的表達量明顯不同,可能參與了幼蟲組織的降解或蛹/成蟲組織的生成等發育調節活動。短管赤眼蜂TpAGO5和TpAGO6與DmAGO1和BmAGO1親緣關系較近,推測赤眼蜂TpAGO5-6可能與果蠅和家蠶AGO1蛋白功能相近,參與胚胎發育的調節。TpAGO4與NvAGO2,DmAGO2和BmAGO2親緣關系較近,推測與果蠅DmAGO2功能相似,參與胚胎發育調節。本研究基于EST數據庫揭示的短管赤眼蜂TpAGO基因表達結果發現, 供試6個TpAGO基因可能參與赤眼蜂性別調控,其中,TpAGO4-6不僅在赤眼蜂卵期高表達,與親緣聚類結果相近,而且在蛹期和成蟲期也具有較高表達;TpAGO1-3多在蛹期和成蟲期有較高表達,在卵期沒有表達。上述結果表明赤眼蜂6個TpAGO有可能都參與分裂期性別調控,但是不同亞家族在生長發育階段具有明確的分工。

另外,TpAGO蛋白基因結構表現明顯的多樣性,基因外顯子數量差異明顯,從1到20個數量不等,可能與該家族基因Ka/Ks值偏高,最低的為0.732 1,存在較高的進化選擇壓力有關。其中,TpAGO2和TpAGO5,TpAGO3和TpAGO4兩對基因間Ka/Ks值大于1,進化動力為正選擇作用,其他組合間表現為純化選擇。值得一提的是,TpAGO4與其他的TpAGO蛋白基因相差較大,如編碼蛋白長度最長、分子量、外顯子數最多,在系統進化樹上與果蠅AGO2同源,可能參與siRNA通路對胚胎發育的調控。

[1] 陳學新. 21世紀我國害蟲生物防治研究的進展、問題與展望[J]. 應用昆蟲學報, 2010, 47(4):615-625.

[2] 張海燕, 張瑩, 叢斌,等. 外界生態因子對感染Wolbachia的松毛蟲赤眼蜂生殖穩定性影響[J]. 中國農業科學, 2009, 42(7):2366-2372.

[3] 陳茜, 王麗艷, 楊志強,等. 溫度通過影響Wolbachia滴度調控赤眼蜂生殖方式[J]. 昆蟲學報, 2016, 59(4):464-471.

[4] Werren J, Baldo L, Clark M.Wolbachia: master manipulators of invertebrate biology[J]. Nature Reviews Microbiology, 2008, 6(10):741-751.

[5] Zhang Guangmei, Hussain M, O’Neill S L, et al.Wolbachiauses a host microRNA to regulate transcripts of a methyltransferase, contributing to dengue virus inhibition inAedesaegypti[J]. Proceedings of the National Academy of Sciences, 2013, 110(25):10276-10281.

[6] Sasaki T, Shiohama A, Minoshima S, et al. Identification of eight members of the Argonaute family in the human genome[J]. Genomics, 2003, 82(3):323-330.

[7] Vaucheret H. Plant Argonautes [J]. Trends in Plant Science, 2008, 13(7):350-358.

[8] Wang Genhong, Jiang Liang, Zhu Li, et al. Characterization of Argonaute family members in the silkworm,Bombyxmori[J]. Insect Science, 2013, 20(1):78-91.

[9] Yigit E, Batista P J, Bei Y, et al. Analysis of theC.elegansArgonaute family reveals that distinct Argonautes act sequentially during RNAi [J]. Cell, 2006, 127(4):747-757.

[10] Forstemann K, Horwich M D, Wee L, et al.DrosophilamicroRNAs are sorted into functionally distinct argonaute complexes after production by Dicer-1 [J]. Cell, 2007, 130(2):287-297.

[11] Jinek M, Doudna J A.A three-dimensional view of the molecular machinery of RNA interference [J].Nature,2009, 457(7228): 405-412.

[12] Saito K, Nishida K M, Mori T, et al. Specific association of PIWI with rasiRNAs derived from retrotransposon and heterochromatic regions in theDrosophilagenome [J]. Genes & Development, 2006, 20(16):2214-2222.

[13] Yin Hang, Lin Haifan. An epigenetic activation role of PIWI and a PIWI-associated piRNA inDrosophilamelanogaster[J]. Nature, 2007, 450(7167):304-308.

[14] Gunawardane L S, Saito K, Nishida K M, et al. A slicer-mediated mechanism for repeat-associated siRNA 5′ end formation inDrosophila[J]. Science, 2007, 315(5818):1587-1590.

[15] Zamore P D.RNA silencing: genomic defence with a slice of pi[J]. Nature, 2007, 446(7138):864-865.

[16] Williams R W, Rubin G M. ARGONAUTE1 is required for efficient RNA interference inDrosophilaembryos [J]. Proceedings of the National Academy of Sciences of the United States of America, 2002, 99(10):6889-6894.

[17] Pal-Bhadra M, Bhadra U, Birchler J A.RNAi related mechanisms affect both transcriptional and posttranscriptional transgene silencing inDrosophila[J]. Molecular Cell, 2002, 9(2):315-327.

[18] Cox D N, Chao A, Baker J, et al. A novel class of evolutionarily conserved genes defined byPIWIare essential for stem cell self-renewal [J].Genes & Development,1998,12:3715-3727.

[19] Schmidt A, Palumbo G, Bozzetti M P, et al. Genetic and molecular characterization of stingv, a gene involved in crystal formation and meiotic drive in the male germ line ofDrosophilamelanogaster[J]. Genetics, 1999, 151(2):749-760.

[20] Kawaoka S, Minami K, Katsuma S, et al. Developmentally synchronized expression of twoBombyxmoriPiwi subfamily genes, SIWI and BmAGO3 in germ-line cells [J]. Biochemical and Biophysical Research Communications,2008,367(4):755-760.

[21] Zhu L, Masaki Y, Tatsuke T.A MC motif in silkworm Argonaute 1 is indispensible for translation repression[J]. Insect Molecular Biology, 2013,22(3):320-330.

(責任編輯： 田 喆)

IdentificationandexpressionanalysisofArgonauteproteingenefamilyinTrichogrammapretiosum

Yang Zhiqiang, Wang Liyan, Liang Hongyu, Liu Xu, Bai Xueping, Zhao Changjiang, Zhang Haiyan

(CollegeofAgronomy,KeyLaboratoryofCropGermplasmImprovementandCultivationinColdRegionsofHeilongjiangProvincialHigher-EducationInstitutions,HeilongjiangBayiAgriculturalUniversity,Daqing163319,China)

The Argonaute family play an important role in gene silencing of RNA induction in plants and animals. Six Argonaute proteins were identified and named in genome-wide level inTrichogrammapretiosumby using the bioinformatics methods. Analysis results of protein structure and phylogenetic analysis showed that TpAGO family proteins were divided into two subgroups, PIWI and Ago, and the numbers of secondary structures of TpAGO protein were similar, while the distribution was different among TpAGOs, especially the helix structure was differently distributed. The gene structure of TpAGOs showed obvious diversity, and the exon numbers varied significantly, ranging from 1 to 20. The evolutionary powers of two pairs of genes were positive selection, including TpAGO2 and TpAGO5, TpAGO3 and TpAGO4, and other combinations were driven by purification options. Expressional profiles by ESTs showed that TpAGO was expressed in different developmental stages, including egg, larval, pupal and adult stages, and worked in single gender (female or male) in the same species, which suggested that TpAGOs were involved in sex and development regulation ofT.pretiosum. This research helps to reveal the genetic mechanism of thelytoky pattern and promote the biocontrol application ofTrichogramma.

Trichogrammapretiosum; Argonaute; gene family; bioinformatics analysis

2017-02-05

2017-02-22

國家自然科學基金(31301714);高等學校博士學科點專項科研基金(20132305120001);黑龍江省自然科學基金(C2015039);黑龍江省大學生創新創業訓練項目(201610223039)

* 通信作者 E-mail:zhanghy51@126.com

#為并列第一作者

S 476.3

A

10.3969/j.issn.0529-1542.2017.06.006