(陳二虎 侯秋莉

(南京財(cái)經(jīng)大學(xué)食品科學(xué)與工程學(xué)院；江蘇省現(xiàn)代糧食流通與安全協(xié)同創(chuàng)新中心；江蘇高校糧油質(zhì)量安全控制及深加工重點(diǎn)實(shí)驗(yàn)室1，南京 210023)

(揚(yáng)州大學(xué)園藝與植物保護(hù)學(xué)院2，揚(yáng)州 225009)

表皮蛋白(CPs)是昆蟲表皮的重要結(jié)構(gòu)物質(zhì)，與幾丁質(zhì)相互交聯(lián)形成復(fù)雜穩(wěn)定的結(jié)構(gòu)，以維持表皮彈性和其他物理性質(zhì)[1]。大量研究已經(jīng)證實(shí)CPs在昆蟲生長(zhǎng)發(fā)育、表皮骨化、爬行能力、體型塑造、先天免疫、藥劑抗性等方面均發(fā)揮著重要的生理功能[2-6]。此外，CPs只存在于昆蟲等無(wú)脊椎動(dòng)物中，且不同昆蟲間序列差異大，因此是開發(fā)高靶向性殺蟲劑的理想靶標(biāo)。米象(Sitophilusoryzae)是當(dāng)前世界范圍內(nèi)危害最為嚴(yán)重的儲(chǔ)糧害蟲之一，其日益突出的磷化氫抗性問(wèn)題嚴(yán)重威脅糧食安全[7]。針對(duì)上述現(xiàn)狀，挖掘新的防效好、特異性強(qiáng)、安全性高的殺蟲靶標(biāo)則顯得尤為亟需。

伴隨高通量測(cè)序技術(shù)和生物信息學(xué)的快速發(fā)展，迄今為止已從不同昆蟲的基因組或轉(zhuǎn)錄組數(shù)據(jù)庫(kù)中鑒定獲得大量CPs，包括赤擬谷盜(Triboliumcastaneum)(151個(gè))[8]、馬鈴薯甲蟲(Leptinotarsadecemlineata)(175個(gè))[9]、家蠶(Bombyxmori)(226個(gè))[10]、黑腹果蠅(Drosophilamelanogaster)(174個(gè))[11]、東亞飛蝗(Locustamigratoria)(81個(gè))[12]、等。依據(jù)氨基酸序列的保守結(jié)構(gòu)域，昆蟲CPs被劃分為10余個(gè)家族類群，包括CPR[13,14]、CPAP1(Cuticular Proteins Analogous to Peritrophin 1)[5]、CPAP3(Cuticular Proteins Analogous to Peritrophin 3)、CPLCA、CPLCG、CPLCP、CPLCW[15]、Tweedle、CPF[16,17]、CPFL、CPCFC、CPG[18,19]等。

鑒于昆蟲CPs數(shù)量眾多(約占昆蟲基因數(shù)量的1%)，而且不同昆蟲類群間表皮蛋白基因的數(shù)量和種類差異顯著，因此昆蟲CPs的鑒定、序列分析和家族分類是開展后續(xù)功能研究的重要前提。目前，雖然米象全基因組數(shù)據(jù)已公布，但是目前有關(guān)米象表皮蛋白的研究仍相對(duì)滯后，米象CPs的基因信息尚屬未知。因此，本研究擬基于米象基因組數(shù)據(jù)庫(kù)，開展CPs的鑒定、序列特征及其系統(tǒng)發(fā)育分析，研究結(jié)果可為后續(xù)的功能研究奠定基礎(chǔ)，進(jìn)而為米象的防治提供新的思路和方法，具有重要的理論和應(yīng)用參考價(jià)值。

1 材料與方法

1.1 昆蟲基因信息數(shù)據(jù)來(lái)源

米象基因組數(shù)據(jù)源于National Center for Biotechnology Information(NCBI)公共數(shù)據(jù)庫(kù)(NCBI Assembly: GCF_002938485.1)。其他昆蟲CP基因信息均從NCBI的GnenBank數(shù)據(jù)庫(kù)(http://www.ncbi.nlm.nih.gov/)或者VectorBase數(shù)據(jù)庫(kù)(https://www.vectorbase.org/)下載。

1.2 米象表皮蛋白基因的全基因組鑒定

基于米象基因組數(shù)據(jù)信息，通過(guò)同源序列比對(duì)方法鑒定獲得米象CPs序列(E-value < 1×10-5)。CPR家族以擴(kuò)展保守基序R&R(pfam00379)作為參考序列進(jìn)行BLASTP比對(duì)[20, 21]，運(yùn)用CuticleDB網(wǎng)站http://bioInformatics2.biol.uoa.gr/CuticleDB/index.jsp將CPR家族基因劃分為RR-1和RR-2亞家族[11]；關(guān)于CPAP家族，基于GenBank數(shù)據(jù)庫(kù)下載的赤擬谷盜和馬鈴薯甲蟲CPAP1和CPAP3序列，通過(guò)與米象基因組比對(duì)鑒定得到相應(yīng)CP家族基因[22]；基于昆蟲保守的Tweedle基序(PF03103)識(shí)別米象同源的Tweedle家族[3]；通過(guò)與黑腹果蠅同源序列比對(duì)獲得米象CPLCG和CPLCP家族序列[15]；以CPF家族蛋白特征序列中最保守部分(VSxYSKAVDTPFSSVRKxDxRIVNxA)鑒定獲得米象CPF家族基因[17]；利用赤擬谷盜和馬鈴薯甲蟲CPs序列，經(jīng)同源比對(duì)獲得米象CPFL家族基因序列[9]；依據(jù)保守YPAGVNPAACPNYPYCD基序識(shí)別米象CPCFC家族[19]。

1.3 信號(hào)肽及保守結(jié)構(gòu)域分析

使用SignalP 4.0(http://www.cbs.dtu.dk/services/SignalP-4.0/)[23]預(yù)測(cè)米象CPs信號(hào)肽。分別提取CPs家族保守結(jié)構(gòu)域序列，運(yùn)用ClustalW軟件進(jìn)行多重序列比對(duì)，并結(jié)合在線工具WebLogo(http://weblogo.berkeley.edu/logo.cgi)繪制保守結(jié)構(gòu)域LOGO圖，分析CPs各家族序列特征及保守氨基酸的出現(xiàn)頻率。

1.4 系統(tǒng)發(fā)育分析

分別獲取昆蟲(米象、赤擬谷盜、馬鈴薯甲蟲、黑腹果蠅、家蠶等)CPR家族保守R&R基序，CPAP1和CPAP3家族保守ChtBD2結(jié)構(gòu)域，以及Tweedle、CPLCG、CPLCP、CPCFC、CPF和CPFL家族不包含信號(hào)肽的全長(zhǎng)氨基酸序列。運(yùn)用ClustalW對(duì)上述對(duì)應(yīng)氨基酸序列進(jìn)行多重比對(duì)。利用Mega 6.0軟件中的鄰接法(Neighbor-Joining)分別開展昆蟲不同家族CPs氨基酸序列的系統(tǒng)發(fā)育樹構(gòu)建，其中各樹分支都進(jìn)行1 000次重復(fù)抽樣檢驗(yàn)。

2 結(jié)果與分析

2.1 米象CP基因的全基因組鑒定

基于米象基因組共鑒定獲得135個(gè)CP基因，隸屬9個(gè)家族：CPR(RR1，RR2和CPRNC)、Tweedle、CPAP1、CPAP3、CPLCG、CPLCP、CPF、CPFL和CPCFC，其中有85.93%(116個(gè))的CPs包含信號(hào)肽。通過(guò)與鞘翅目昆蟲赤擬谷盜和馬鈴薯甲蟲的CPs比較分析發(fā)現(xiàn)，不同昆蟲間CPs的數(shù)量存在較大差異，馬鈴薯甲蟲(175個(gè))CPs個(gè)數(shù)最多，赤擬谷盜(151個(gè))次之，米象(135個(gè))最少，其中不同昆蟲間尤以CPR家族基因的數(shù)量差異最為明顯(表1)。此外，本研究發(fā)現(xiàn)馬鈴薯甲蟲的CPs缺少CPLC家族，米象則無(wú)CPLCA家族基因(表1)。

表1 三種鞘翅目昆蟲表皮蛋白基因匯總表

2.2 CPR家族基因

本研究共鑒定獲得含保守R&R基序的CPR家族基因93個(gè)，占米象CP基因總數(shù)的68.89%。依據(jù)CuticleDB運(yùn)算，米象CPR家族分別包含43個(gè)RR1亞家族和45個(gè)RR2亞家族成員，另有5條基因序列被定義為CPRNC(未分類CPR)(表1)。此外，系統(tǒng)發(fā)育分析結(jié)果顯示RR1、RR2和CPRNC各自聚為一支，證實(shí)上述CPR亞家族分類的準(zhǔn)確性(圖1a)。有趣的是，CuticleDB數(shù)據(jù)庫(kù)分類結(jié)果顯示XP_030753085.1蛋白隸屬于RR1亞家族，而進(jìn)化樹分析結(jié)果則與CPRNC蛋白聚為一支，暗示傳統(tǒng)CPR亞家族分類方法與系統(tǒng)發(fā)育關(guān)系存在不一致性(圖1a)。米象RR1和RR2家族蛋白的R&R保守基序分別如圖1b和圖1c所示，其中RR1家族蛋白擁有8個(gè)保守氨基酸識(shí)別位點(diǎn)(YTADENGF)，RR2家族蛋白則含有2個(gè)保守識(shí)別序列(EERDGDVVKG；3個(gè)G-x(3)-VV重復(fù)序列)。為進(jìn)一步探究昆蟲RR1和RR2蛋白的進(jìn)化關(guān)系，本研究以米象、赤擬谷盜、黑腹果蠅和家蠶的R&R基序進(jìn)行系統(tǒng)發(fā)育分析，分別構(gòu)建不同昆蟲RR1和RR2家族蛋白系統(tǒng)發(fā)育樹(圖2和圖3)，結(jié)果表明上述不同昆蟲均形成相應(yīng)物種特異性基因簇，米象RR1蛋白包含4個(gè)基因聚類(圖2)，RR2蛋白則擁有8個(gè)獨(dú)立基因簇(圖3)。

注：a 米象CPR家族蛋白系統(tǒng)發(fā)育分析(鄰接法)；b 米象RR1蛋白R(shí)&R基序分析；c 米象RR2蛋白R(shí)&R基序分析。

注：米象(So)，赤擬谷盜(Tc)，黑腹果蠅(Dm)，家蠶(Bm)；黑括號(hào)表示基因聚為一簇。

2.3 CPAP1和CPAP3家族基因

在米象基因組中共鑒定獲得13個(gè)CPAP1和7個(gè)CPAP3家族蛋白基因(表1)。系統(tǒng)發(fā)育分析結(jié)果表明，不同昆蟲CPAP1和CPAP3家族蛋白各自形成兩個(gè)不同分支的基因簇，且米象與赤擬谷盜親緣關(guān)系較高，兩個(gè)物種的CPAP基因在進(jìn)化樹中均成對(duì)出現(xiàn)(圖4a)。蛋白序列結(jié)構(gòu)分析結(jié)果顯示，米象CPAP1和CPAP3家族蛋白序列分別擁有1個(gè)和3個(gè)保守ChtBD2型幾丁質(zhì)結(jié)合域，且多數(shù)蛋白均含有1個(gè)N末端信號(hào)肽(除SoCPAP1-B、F、J和M)(圖4b和圖4c)。此外，米象CPAP1序列長(zhǎng)度差異較大，氨基酸長(zhǎng)度介于170和169 3 aa之間(圖4b)；CPAP3氨基酸序列長(zhǎng)度的變異則較小，其范圍為230～281 aa(圖4c)。與此同時(shí)，米象CPAP1和CPAP3蛋白家族間的幾丁質(zhì)結(jié)合域數(shù)量和位置存在顯著區(qū)別，即SoCPAP3蛋白的3個(gè)幾丁質(zhì)結(jié)合域主要靠近N端，而SoCPAP1蛋白幾丁質(zhì)結(jié)合域在N端、C端和中間位置均有分布(圖4b和圖4c)。

注：a 昆蟲CPAP1和CPAP3家族蛋白系統(tǒng)發(fā)育分析。米象(So)，赤擬谷盜(Tc)，黑腹果蠅(Dm)，地中海實(shí)蠅Ceratitis capitata(Cc)，家蠶(Bm)，煙草天蛾Manduca sexta(Ms)；b CPAP1和c CPAP3蛋白結(jié)構(gòu)分析。

2.4 Tweedle家族基因

在米象基因組中共鑒定得到6個(gè)Tweedle蛋白基因，根據(jù)米象、赤擬谷盜等11種不同昆蟲的氨基酸序列構(gòu)建Tweedle家族蛋白的系統(tǒng)進(jìn)化樹。結(jié)果顯示，昆蟲Tweedle家族蛋白可形成三個(gè)并列基因簇，分別包括蠅蟲類(地中海實(shí)蠅和黑腹果蠅)、蚊子類(尖音庫(kù)蚊、埃及伊蚊和岡比亞按蚊)的特異性Tweedle基因聚類分支，以及由11種昆蟲共同組成的Tweedle基因混合簇(圖5a)。此外，氨基酸序列分析結(jié)果顯示，米象Tweedle家族蛋白擁有4個(gè)保守氨基酸區(qū)域，且每個(gè)Tweedle家族成員中均存在一個(gè)內(nèi)部重復(fù)結(jié)構(gòu)(圖5b)。

注：a 不同昆蟲Tweedle家族蛋白進(jìn)化關(guān)系分析，包括米象(So)、赤擬谷盜(Tc)、馬鈴薯甲蟲Leptinotarsa decemlineata(Ld)、黑腹果蠅(Dm)、地中海實(shí)蠅Ceratitis capitata(Cc)、岡比亞按蚊Anopheles gambiae(Ag)、埃及伊蚊Aedes aegypti(Aa)、尖音庫(kù)蚊Culex pipiens(Cp)、家蠶(Bm)、麗蠅蛹集金小蜂Nasonia vitripennis(Nv)和意大利蜜蜂Apis mellifera(Am)；b 米象Tweedle蛋白保守序列分析。

2.5 低復(fù)雜度CP家族基因

在米象基因組中共鑒定得到米象低復(fù)雜度CP家族基因7個(gè)，包括1個(gè)CPLCG基因和6個(gè)CPLCP基因，未鑒定到CPLCA基因(表1)。不同昆蟲CPLCG和CPLCP蛋白家族分別聚類形成單獨(dú)分支(圖6a)，并且聚集形成4個(gè)蚊類基因簇(圖6a)，表明該家族蛋白數(shù)量在蚊子類群中得到顯著擴(kuò)展。氨基酸保守序列分析顯示，CPLCP蛋白含有高密度脯氨酸PV和PY重復(fù)序列(圖6b)。

注：a 不同昆蟲CPLCG和CPLCP家族蛋白進(jìn)化關(guān)系分析，包括米象(So)、赤擬谷盜(Tc)、岡比亞按蚊(Ag)、埃及伊蚊(Aa)、尖音庫(kù)蚊(Cp)、黑腹果蠅(Dm)、家蠶(Bm)、麗蠅蛹集金小蜂(Nv)和意大利蜜蜂(Am)；b 米象CPLCP蛋白保守序列分析。

2.6 CPF和CPFL家族基因

在米象基因組共鑒定得到1個(gè)CPF和5個(gè)CPFL家族蛋白基因(表1)。依據(jù)不同昆蟲CPF和CPFL的蛋白序列構(gòu)建進(jìn)化樹，結(jié)果顯示這兩個(gè)家族蛋白分別聚集形成兩個(gè)并列分支(圖7a)。序列分析發(fā)現(xiàn)，鞘翅目昆蟲米象、赤擬谷盜和馬鈴薯甲蟲CPF蛋白分別含有保守的44個(gè)氨基酸和C末端序列(圖7b和圖7c)，而CPFL家族蛋白則僅擁有與CPF蛋白同源性較高的C末端保守序列(圖7d)。

注：a 昆蟲CPF和CPFL家族蛋白進(jìn)化關(guān)系分析，包括米象(So)、赤擬谷盜(Tc)、馬鈴薯甲蟲(Ld)、黑腹果蠅(Dm)和地中海實(shí)蠅(Cc)；(b和c)米象、赤擬谷盜和馬鈴薯甲蟲CPF保守結(jié)構(gòu)域分析；d 米象、赤擬谷盜和馬鈴薯甲蟲CPFL保守結(jié)構(gòu)域分析。

2.7 CPCFC家族基因

本研究在米象基因組中共注釋得到3個(gè)CPCFC家族蛋白基因(表1)。為明確CPCFC蛋白的系統(tǒng)發(fā)育關(guān)系，共選擇9種雙翅目和3種鞘翅目昆蟲的CPCFC蛋白序列開展系統(tǒng)發(fā)育分析，結(jié)果表明這兩大類昆蟲的CPCFC蛋白各自聚集形成一簇(圖8a)。米象、赤擬谷盜和馬鈴薯甲蟲)3種鞘翅目昆蟲的CPCFC氨基酸序列分析結(jié)果顯示，該家族蛋白具有共同的保守結(jié)構(gòu)域，即兩個(gè)重復(fù)的C-X5-C基序(圖8b)。

注：a 昆蟲CPCFC家族蛋白進(jìn)化關(guān)系分析，包括米象(So)、赤擬谷盜(Tc)、馬鈴薯甲蟲(Ld)、桔小實(shí)蠅Bactrocera dorsalis(Bd)、瓜實(shí)蠅Bactrocera cucurbitae(Bc)、地中海實(shí)蠅(Cc)、黑腹果蠅(Dm)、埃及伊蚊(Aa)、岡比亞按蚊(Ag)、不吉按蚊Anopheles funestus(Af)、中華按蚊Anopheles sinensis(As)、四斑按蚊An.quadrimaculatus(Aq)；b 米象、赤擬谷盜和馬鈴薯甲蟲CPCFC保守結(jié)構(gòu)域分析。

3 討論

作為昆蟲表皮重要結(jié)構(gòu)物質(zhì)，CPs不僅種類數(shù)量眾多，而且在昆蟲生長(zhǎng)發(fā)育和環(huán)境適應(yīng)等方面均發(fā)揮關(guān)鍵作用，是開發(fā)新型害蟲控制劑的理想靶標(biāo)[24, 25]。本研究運(yùn)用生物信息學(xué)手段從米象基因組數(shù)據(jù)中共鑒定得到135個(gè)CPs基因，隸屬于9個(gè)家族，包括CPR、Tweedle、CPAP1、CPAP3、CPLCG、CPLCP、CPF、CPFL和CPCFC。

本研究發(fā)現(xiàn)同屬鞘翅目昆蟲的米象、赤擬谷盜和馬鈴薯甲蟲的CPs數(shù)量差異明顯。已有研究證實(shí)CPs基因的存在可能與昆蟲應(yīng)對(duì)外界復(fù)雜環(huán)境密切相關(guān)，因此造成物種間CPs數(shù)量和種類巨大差別的可能原因是不同種類昆蟲所處生存環(huán)境的差異[26,27]。由于米象和赤擬谷盜均為儲(chǔ)糧害蟲，其活動(dòng)范圍主要集中在環(huán)境穩(wěn)定的糧倉(cāng)內(nèi)，相比而言馬鈴薯甲蟲的生存環(huán)境則更為復(fù)雜，因此其擁有CPs數(shù)量也較多。此外，與赤擬谷盜相比，米象幼蟲期和蛹期均于稻谷內(nèi)部進(jìn)行發(fā)育，其相對(duì)穩(wěn)定和受保護(hù)的生長(zhǎng)環(huán)境可能導(dǎo)致該昆蟲進(jìn)化出的CPs數(shù)量少于赤擬谷盜。

米象CPs另外一個(gè)重要特點(diǎn)是諸多CPR家族蛋白在系統(tǒng)發(fā)育上呈現(xiàn)聚類現(xiàn)象，即在進(jìn)化樹上形成許多相應(yīng)物種的特異性CPs基因簇，這意味著聚集的CPs基因具有較高的序列相似度和協(xié)調(diào)一致的進(jìn)化模式。前人亦分別在家蠶、斜紋夜蛾、馬尾松毛蟲、岡比亞按蚊等昆蟲中發(fā)現(xiàn)類似的物種特異性CPR基因聚類，經(jīng)證實(shí)這些基因均來(lái)源于染色體上的同一個(gè)基因簇，暗示其可能是從同一個(gè)祖先基因進(jìn)化而來(lái)，且在各物種中分別進(jìn)行擴(kuò)增，結(jié)果進(jìn)一步證明CPs基因的復(fù)制是在各昆蟲類群中獨(dú)立發(fā)生[10, 15, 27-29]。

除CPRs外，CPAPs是昆蟲第二大類幾丁質(zhì)結(jié)合蛋白，其依據(jù)氨基酸序列所包含幾丁質(zhì)結(jié)合域(ChtBD2)數(shù)量，被進(jìn)一步分為CPAP1(包含1個(gè)ChtBD2)和CPAP3(3個(gè)ChtBD2)兩個(gè)家族[22]。然而，本研究發(fā)現(xiàn)昆蟲CPAPs與CPR蛋白之間的系統(tǒng)發(fā)育關(guān)系存在巨大差異，即不同昆蟲CPAP1和CPAP3同一家族成員間的親緣關(guān)系要強(qiáng)于同一物種各個(gè)蛋白間的關(guān)聯(lián)度，顯示出這些蛋白是依據(jù)某些特定保守結(jié)構(gòu)域進(jìn)化而來(lái)，暗指CPAPs發(fā)展成為多基因家族的過(guò)程要早于不同種類昆蟲的分化過(guò)程。之前的研究發(fā)現(xiàn)昆蟲CPAP1家族基因的數(shù)量在不同物種間變異較大(9-17個(gè))[30]，本論文中米象較赤擬谷盜缺少CPAP1-D和CPAP1-E基因。

與其他CPs相比，不同種類昆蟲間CPAP3家族基因的數(shù)量、種類及其氨基酸序列結(jié)構(gòu)最為保守[30,31]，并且已有研究證實(shí)CPAP3蛋白與幾丁質(zhì)均有較高的結(jié)合能力[32]，暗示該家族蛋白可能通過(guò)與表皮中的幾丁質(zhì)結(jié)合從而在節(jié)肢動(dòng)物的生長(zhǎng)發(fā)育過(guò)程中發(fā)揮重要作用。例如，果蠅CPAP3-A和CPAP3-E蛋白通過(guò)與幾丁質(zhì)相互作用來(lái)控制昆蟲表皮完整性、表皮外骨骼物理性質(zhì)以及幼蟲體型[33, 34]。赤擬谷盜CPAP3家族蛋白則承擔(dān)著多樣的生理功能，利用RNA干擾技術(shù)有效沉默該家族不同基因會(huì)對(duì)昆蟲鞘翅、表皮角質(zhì)層、腸道、卵巢、脂肪體、足、胚胎等造成嚴(yán)重發(fā)育缺陷，從而導(dǎo)致死亡表型[5]。這些研究結(jié)果證明昆蟲特定的關(guān)鍵CPs，如CPAP3家族基因可成為基于dsRNA控制米象的潛在靶位點(diǎn)，后續(xù)可進(jìn)一步針對(duì)該家族基因開展功能研究，驗(yàn)證其在米象生理發(fā)育過(guò)程中的重要功能。

4 結(jié)論

表皮覆蓋于蟲體外表面，承擔(dān)昆蟲外骨骼的功能，在昆蟲生長(zhǎng)發(fā)育過(guò)程中發(fā)揮重要作用，并可抵御外界不良環(huán)境，如病原體、殺蟲劑等，是昆蟲的一道天然保護(hù)屏障。幾丁質(zhì)和不同種類的CPs的比例是決定昆蟲表皮結(jié)構(gòu)、特性和機(jī)械性能的關(guān)鍵因素。因此，本論文針對(duì)米象CPs的全基因組鑒定研究不僅可以豐富昆蟲CP超基因家族數(shù)據(jù)，還可為開展米象CP基因功能及生物學(xué)研究提供基礎(chǔ)的基因序列信息，最終為米象等儲(chǔ)糧害蟲的新型防控策略奠定重要基礎(chǔ)。