鐘 勇，馬福歡，藍(lán)翊文，何成偉，白華菊

1憑祥海關(guān)，廣西 憑祥 532600；2友誼關(guān)海關(guān)，廣西 憑祥 532600

胸喙亞目昆蟲(包括粉蚧、其他蚧蟲、蚜蟲、木虱和粉虱)以植物汁液為食，植物汁液營(yíng)養(yǎng)上極不平衡，雖然富含主要以蔗糖存在的碳水化合物，但脂質(zhì)和蛋白質(zhì)含量非常少，因此，它們均需內(nèi)共生菌提供必需氨基酸和其他營(yíng)養(yǎng)物質(zhì)來(lái)保證正常的生長(zhǎng)發(fā)育(Baumannetal.,2000; Moran & Telang,1998；Munsonetal.,1992)。研究發(fā)現(xiàn)，粉蚧體內(nèi)含有初生內(nèi)共生菌(Fukatsu & Nikoh,2000; Moran & Telang,1998；Munsonetal.,1992)和次生內(nèi)共生菌(Kanthetietal.,1996；von Dohlenetal.,2001)。在研究昆蟲內(nèi)共生菌的功能及其與宿主間的相互關(guān)系時(shí)，傳統(tǒng)的方法是通過去除內(nèi)共生菌或者利用營(yíng)養(yǎng)缺失實(shí)驗(yàn)來(lái)觀察宿主生理生化的變化。然而，內(nèi)共生菌不能體外培養(yǎng)，這阻礙了進(jìn)一步深入了解共生菌在宿主體內(nèi)所起的作用(楊義婷等，2014)。隨著PCR技術(shù)的發(fā)展和測(cè)序技術(shù)的進(jìn)步，轉(zhuǎn)錄組學(xué)、蛋白質(zhì)組學(xué)、基因組學(xué)等技術(shù)的發(fā)展，為昆蟲內(nèi)共生菌研究提供了便利；最新的基因組學(xué)工具，特別是新一代高通量測(cè)序工具和分析軟件，使得內(nèi)共生菌的大規(guī)模宏基因組測(cè)序和價(jià)格低廉的單個(gè)菌基因組測(cè)序成為可能，為明確共生菌功能提供了直接證據(jù)，同時(shí)也促使從種群和群落水平研究細(xì)菌、真菌等與昆蟲宿主共生關(guān)系成為可能，而這些都將為最終發(fā)現(xiàn)共生菌與宿主共生關(guān)系的潛在機(jī)制奠定基礎(chǔ)。

扶桑綿粉蚧PhenacoccussolenopsisTinsley隸屬半翅目Hemiptera蚧總科Coccoidea粉蚧科Pseudococcidae綿粉蚧亞科Phenacoccinae綿粉蚧屬Phenacoccus，是一種重要的入侵害蟲，為害棉花GossypiumhirsutumL.和其他重要作物，現(xiàn)已幾乎遍布全球(Fuchsetal.，1991；Wang & Yang,2010)。目前，國(guó)內(nèi)外對(duì)其研究多集中在生物學(xué)特性方面，內(nèi)共生菌的研究相對(duì)較少(黃芳等，2015)。這些報(bào)道中，如扶桑綿粉蚧入侵后寄主范圍不斷擴(kuò)大(陳哲等，2012)、出現(xiàn)體色變異型種群(朱藝勇等，2011；Zhaoetal.,2014)等可能均有內(nèi)共生菌參與。同時(shí)，不同的研究對(duì)扶桑綿粉蚧繁殖方式的報(bào)道有較大差異(關(guān)鑫等，2011；唐遠(yuǎn)等，2013; Hodgsonetal.,2008；Huangetal.,2013；Vennilaetal.,2010)，其原因可能是不同種群內(nèi)共生菌感染情況不同，尤其是與生殖有關(guān)的內(nèi)共生菌的感染情況。因此，扶桑綿粉蚧內(nèi)共生菌多樣性的研究，將為共生菌與粉蚧的相互作用研究提供基礎(chǔ)，也為研究扶桑綿粉蚧的入侵機(jī)制提供新的研究思路。

1 材料與方法

1.1 試蟲來(lái)源

扶桑綿粉蚧為國(guó)家熱帶和亞熱帶水果檢疫重點(diǎn)實(shí)驗(yàn)室(憑祥)使用南瓜Cucurbitamoschata(Duch.ex Lam.)Duch.ex Poiret飼養(yǎng)的室內(nèi)種群，原寄主為木槿HibiscusrosasinensisL.(廣西憑祥)。

1.2 高通量測(cè)序

收集試蟲，分為2個(gè)重復(fù)，每重復(fù)雌成蟲30頭；委托深圳華大基因股份有限公司對(duì)整蟲樣本開展微生物宏基因組測(cè)序。

1.3 測(cè)序數(shù)據(jù)處理

原始測(cè)序數(shù)據(jù)6.28×109reads。過濾掉低質(zhì)量序列、接頭序列等雜質(zhì)，從而得到6.12×109的高質(zhì)量Clean Data。組裝過濾后的reads，以得到更長(zhǎng)的contig。然后將reads分別比對(duì)回每個(gè)組裝結(jié)果以評(píng)估組裝效果，最后綜合考慮N50和比對(duì)率情況，選擇一個(gè)最佳k-mer和其對(duì)應(yīng)的組裝結(jié)果。最終組裝獲得1.60×105contigs，組裝長(zhǎng)度2.56×108bp，最大contig長(zhǎng)度89956 bp，最小contig長(zhǎng)度300 bp，組裝平均長(zhǎng)度1600 bp，reads利用率68.32%。

1.4 內(nèi)共生菌豐度分析

使用MEGAN(version 4.6)對(duì)nr比對(duì)結(jié)果進(jìn)行處理。MEGAN根據(jù)LCA算法對(duì)BLAST比對(duì)結(jié)果根據(jù)NCBI物種分類樹進(jìn)行重新梳理，從而得到每條序列的物種注釋信息。將同一物種注釋的基因豐度累加，得到該物種在樣品中的物種豐度。

1.5 內(nèi)共生菌基因分析

對(duì)每個(gè)樣品使用MetaGeneMark進(jìn)行基因預(yù)測(cè)，構(gòu)建高質(zhì)量基因集，將reads比對(duì)到高質(zhì)量基因集計(jì)算基因豐度。對(duì)KEGG(89.1)、COG等數(shù)據(jù)庫(kù)的分析，先將基因注釋到此類數(shù)據(jù)庫(kù)?；赗eads與基因集的比對(duì)結(jié)果,獲取該類功能數(shù)據(jù)庫(kù)的定量(豐度)和定性(基因個(gè)數(shù))信息。

2 結(jié)果與分析

2.1 扶桑綿粉蚧內(nèi)共生菌種類及其豐度

根據(jù)注釋結(jié)果，樣本共鑒定獲得29個(gè)門，47個(gè)綱，105個(gè)目，178個(gè)科，245個(gè)屬，299個(gè)種。在界分類水平上，扶桑綿粉蚧內(nèi)共生物測(cè)序共注釋了古細(xì)菌Archaea、細(xì)菌Bacteria、真核生物Eukaryota和病毒Viruses；細(xì)菌為優(yōu)勢(shì)種群，相對(duì)豐度達(dá)85.16%。在門分類水平上，共注釋了子囊菌門Ascomycota、擔(dān)子菌門Basidiomycota、Candidatus Magasanikbacteria等29個(gè)門；其中變形菌門Proteobacteria為優(yōu)勢(shì)種群(相對(duì)豐度81.36%)。在綱分類水平上，共注釋了α-變形菌綱Alphaproteobacteria、β-變形菌綱Betaproteobacteria、衣原體綱Chlamydiia等47個(gè)綱；γ-變形菌綱和β-變形菌綱為優(yōu)勢(shì)種群，相對(duì)豐度遠(yuǎn)超其他種群(分別為50.04%和27.68%)。在目分類水平上，共注釋了交替單孢菌目Alteromonadales、伯克氏菌目Burkholderiales、彎曲菌目Campylobacterales等105個(gè)目；腸桿菌目Enterobacterales為優(yōu)勢(shì)種群(相對(duì)豐度42.65%)。在科分類水平上，共注釋了酸桿菌科Acetobacteraceae、芽孢桿菌科Bacillaceae、彎曲菌科Campylobacteraceae、Debaryomycetaceae等178個(gè)科；腸桿菌科Enterobacteriaceae為優(yōu)勢(shì)種群(相對(duì)豐度33.48%)。在屬分類水平上，共注釋了醋酸桿菌屬Acetobacter、芽孢桿菌屬Bacillus、Caballeronia等245個(gè)屬；伯克氏菌屬Burkholderia為優(yōu)勢(shì)種群(相對(duì)豐度13.12%)；CandidatusTremblaya相對(duì)豐度3.69%排名第二。在種分類水平上，共注釋了Acetobactermalorum、脆弱類桿菌屬Bacteroidesfragilis、Caballeroniagrimmiae等299個(gè)種；相對(duì)豐度為2.58%的CandidatusTremblayaphenacola為優(yōu)勢(shì)種群(表1)。

表1 扶桑綿粉蚧內(nèi)共生菌豐度(前10)

2.2 扶桑綿粉蚧內(nèi)共生菌基因分析

樣本所測(cè)序列共預(yù)測(cè)了36505個(gè)基因，其中77.27%的基因分布在200～499 nt的范圍內(nèi)(圖1)。基因集KEGG注釋了776個(gè)基因，其中代謝相關(guān)基因最多(384個(gè))，占49.48%；其次為遺傳信息處理基因(207個(gè))，占26.68%；最少的為環(huán)境信息處理基因(8個(gè))，僅占1.03%。代謝相關(guān)基因除了綜合性代謝全譜基因(157個(gè))外，最多的為氨基酸代謝基因(64個(gè))，其次為碳水化合物代謝基因(36個(gè))、輔酶和維生素代謝基因(35個(gè))(圖2)，分別占28.19%、15.86%和15.42%?；蚣疌OG注釋了398個(gè)基因，其中翻譯、核糖體結(jié)構(gòu)和生物發(fā)生的基因最多(110個(gè))，占27.64%；其次為氨基酸轉(zhuǎn)運(yùn)和代謝基因(49個(gè))(圖3)，占12.31%。

圖1 扶桑綿粉蚧內(nèi)共生菌基因長(zhǎng)度分布圖

圖2 扶桑綿粉蚧內(nèi)共生菌KEGG通道注釋圖

圖3 扶桑綿粉蚧內(nèi)共生菌COG通道注釋圖

3 討論

本研究采用宏基因組測(cè)序技術(shù)對(duì)扶桑綿粉蚧內(nèi)共生菌多樣性及群落結(jié)構(gòu)組成進(jìn)行了研究。其中，變形菌門、γ和β變形菌，以及β-變形菌的CandidatusTremblaya是其優(yōu)勢(shì)種群，該結(jié)果與其他昆蟲的報(bào)道相似(毛迎新等，2018；Konoetal.,2008；McCutcheon & von Dohlen，2011；Sergioetal.,2015)。

很多變形菌門細(xì)菌與營(yíng)養(yǎng)代謝相關(guān)，這些種群豐度占優(yōu)勢(shì)的內(nèi)共生菌很可能在扶桑綿粉蚧的營(yíng)養(yǎng)代謝中起很大作用。CandidatusTremblaya是一種與綿粉蚧亞科和粉蚧亞科有共生關(guān)系的內(nèi)共生菌，且多數(shù)綿粉蚧亞科粉蚧具有獨(dú)特的內(nèi)共生菌嵌套體CandidatusTremblayaphenacola(Sergioetal.,2015)。實(shí)驗(yàn)結(jié)果證明,在扶桑綿粉蚧體內(nèi)，CandidatusTremblayaphenacola是種的分類水平的優(yōu)勢(shì)種群。此外，同屬的CandidatusBaumanniacicadellinicola和CandidatusEntotheonellasp.相對(duì)豐度排名分別占第6和第10。

研究發(fā)現(xiàn)，在基因KEGG注釋時(shí)，扶桑綿粉蚧基因主要為代謝相關(guān)基因和遺傳信息處理基因，代謝相關(guān)功能基因中，又以氨基酸代謝基因?yàn)橹?，?shù)量遠(yuǎn)遠(yuǎn)高于其他代謝基因。在橘臀紋粉蚧Planococcuscitri(Risso)體內(nèi)，嵌套菌CandidatusMoranellaendobia和CandidatusTremblayaprinceps共同參與了其體內(nèi)的氨基酸代謝(López-Madrigaletal.,2013；McCutcheon & von Dohlen，2011)。然而，同樣的功能在PhenacoccusavenaeBorchsenius體內(nèi)卻由CandidatusTremblayaphenacola一種共生菌完成(López-Madrigaletal.,2014)。為此，可推測(cè)扶桑綿粉蚧通過刺吸式口器取食富含糖分卻缺乏必需氨基酸的植物韌皮部汁液，內(nèi)共生菌CandidatusTremblaya為其提供了必需氨基酸等營(yíng)養(yǎng)物質(zhì)。同時(shí)，研究檢測(cè)到Drosophilaananassae內(nèi)共生菌Wolbachia也是扶桑綿粉蚧體內(nèi)的優(yōu)勢(shì)種群之一(相對(duì)豐度1.04%)。Wolbachia廣泛存在于節(jié)肢動(dòng)物中，對(duì)宿主生殖有重要影響。由于內(nèi)共生菌在昆蟲體內(nèi)的生存環(huán)境相對(duì)穩(wěn)定，與自由生活的親緣菌相比，基因組更小，其中環(huán)境信息處理可由其寄主提供(López-Madrigaletal.,2015)。因此,扶桑綿粉蚧內(nèi)共生菌基因注釋時(shí)環(huán)境信息處理基因所占比例也很小，僅占1.03%。

本研究共預(yù)測(cè)了36505個(gè)基因，而注釋的基因相對(duì)較少(KEGG，776個(gè)；COG，398個(gè))。這可能與研究目標(biāo)和選取的數(shù)據(jù)庫(kù)有關(guān)。而注釋結(jié)果出現(xiàn)B.ubonensis等種類，可能與基因和物種注釋中存在的偏好性有關(guān)。在基因和物種注釋時(shí)，本研究選取E值<0.00001的object序列作為該基因或者物種的注釋信息，因此存在一定的偏好性，注釋豐度較低的基因和物種可能與實(shí)際有一定差異。

昆蟲內(nèi)共生菌能為宿主昆蟲提供營(yíng)養(yǎng)物質(zhì)、消化食物、解毒及提高免疫活性等功能，還影響宿主的生殖方式、適合度、運(yùn)動(dòng)行為及抗藥性(魏曉瑩等，2019)。本研究借助高通量測(cè)序和生物信息學(xué)技術(shù)，初步確定了扶桑綿粉蚧內(nèi)共生菌菌群的群落結(jié)構(gòu)和優(yōu)勢(shì)菌群，并初步分析了這些內(nèi)共生菌的功能基因情況。但這些內(nèi)共生菌對(duì)扶桑綿粉蚧的生物學(xué)特性有哪些影響、在其體內(nèi)的分布以及與其入侵的關(guān)系還需進(jìn)一步研究。