戴 陽(yáng)， 李志紅＊， 柳麗君， 吳佳教， 鄧裕亮

（1．中國(guó)農(nóng)業(yè)大學(xué)農(nóng)學(xué)與生物技術(shù)學(xué)院，北京 100193； 2．廣東檢驗(yàn)檢疫技術(shù)中心，廣州 510623；3．西雙版納出入境檢驗(yàn)檢疫局，景洪 666100）

番石榴果實(shí)蠅［Bactrocer a correcta （Bezzi）］隸屬于雙翅目（Diptera），實(shí)蠅科（Tephritidae），寡鬃實(shí)蠅亞科（Dacinae），寡鬃實(shí)蠅族（Dacini），果實(shí)蠅屬（Bactr ocer a），是為害水果的一類(lèi)重要經(jīng)濟(jì)實(shí)蠅。番石榴果實(shí)蠅主要分布在亞洲的印度、巴基斯坦、泰國(guó)、尼泊爾、斯里蘭卡、緬甸、越南［1］，目前，該蟲(chóng)在我國(guó)主要分布在臺(tái)灣和云南的局部地區(qū)［2－3］。2007年5月28日，農(nóng)業(yè)部與國(guó)家質(zhì)量監(jiān)督檢驗(yàn)檢疫總局公布《中華人民共和國(guó)進(jìn)境植物檢疫性有害生物名錄》，將果實(shí)蠅屬列為檢疫性有害生物，而番石榴果實(shí)蠅是果實(shí)蠅屬中的重要種類(lèi)。

共生菌是存在于昆蟲(chóng)體內(nèi)特定部位的微生物，可分為初級(jí)共生菌和次級(jí)共生菌。初級(jí)共生菌在宿主體內(nèi)垂直傳遞并與宿主協(xié)同進(jìn)化，次級(jí)共生菌同時(shí)存在垂直傳遞和水平傳遞［4］。共生菌通過(guò)生物合成為宿主提供必需的氨基酸等營(yíng)養(yǎng)物質(zhì)，幫助宿主降解有害物質(zhì)［5－6］，在進(jìn)化中趨向同寄主建立共生關(guān)系，許多共生菌還存在與宿主協(xié)同進(jìn)化的現(xiàn)象［7－8］。隨著動(dòng)物胃腸道微生態(tài)理論的發(fā)展，細(xì)菌和昆蟲(chóng)宿主間廣泛存在的互利互惠關(guān)系不斷被發(fā)現(xiàn)［9］。近年來(lái)，國(guó)外研究者針對(duì)油橄欖實(shí)蠅［Bactrocer a oleae（Gmelin）］、地中海實(shí)蠅［Cer atitis capitata （Wiedemann）］等多種實(shí)蠅的共生菌展開(kāi)研究。Kounatidis用16s r DNA文庫(kù)構(gòu)建技術(shù)，通過(guò)油橄欖實(shí)蠅DNA提取、通用引物的PCR擴(kuò)增及產(chǎn)物純化、克隆轉(zhuǎn)化、測(cè)序及比對(duì)，發(fā)現(xiàn)油橄欖實(shí)蠅的優(yōu)勢(shì)共生菌為醋酸桿菌［10］，Daniela等通過(guò)解剖觀察，發(fā)現(xiàn)地中海實(shí)蠅的主要共生菌為聚團(tuán)腸桿菌［11］，Capuzzo等人利用16s r RNA基因序列研究了油橄欖實(shí)蠅的共生菌，發(fā)現(xiàn)一種與宿主協(xié)同進(jìn)化的共生菌并將其命名為“Cantidatus Er winia Dacicola”［12］。

本研究采用16s r DNA文庫(kù)構(gòu)建技術(shù)及克隆測(cè)序方法，檢測(cè)了番石榴果實(shí)蠅室內(nèi)種群共生菌的多樣性，現(xiàn)將結(jié)果報(bào)道如下。

1 材料和方法

1．1 試驗(yàn)材料及主要生化試劑

1.1.1 試驗(yàn)材料

本研究所用的番石榴果實(shí)蠅樣品由廣東檢驗(yàn)檢疫技術(shù)中心提供，為室內(nèi)飼養(yǎng)。取雌蟲(chóng)、雄蟲(chóng)各1頭用于共生菌16s DNA文庫(kù)的構(gòu)建。樣品浸存于無(wú)水乙醇中，－20℃保存?zhèn)溆谩?/p>

1．1．2 主要試劑

血液／細(xì)胞／組織基因組DNA提取試劑盒（離心柱型）、10×Taq Buffer（含15 mmol／L Mg Cl2）、2．5 mmol／L d NTPs、2．5 U／μL Taq DNA 聚合酶、D2000 DNA分子量標(biāo)準(zhǔn)、瓊脂糖凝膠回收試劑盒均購(gòu)自北京天根生化科技公司，p MD19－T載體購(gòu)自寶生物工程（大連）有限公司，DH5α感受態(tài)細(xì)胞購(gòu)自北京全式金生物技術(shù)有限公司，PCR引物由北京奧科生物技術(shù)公司合成。

1．2 16s r DNA文庫(kù)的建立

1．2．1 單頭實(shí)蠅基因組DNA提取

用無(wú)菌水將實(shí)蠅樣品快速洗凈后，采用 “血液／細(xì)胞／組織基因組DNA提取試劑盒（離心柱型）”提取單頭實(shí)蠅成蟲(chóng)樣品的基因組DNA，瓊脂糖凝膠電泳檢測(cè)提取質(zhì)量。

1．2．2 番石榴果實(shí)蠅共生菌16s r DNA序列擴(kuò)增

參考Willam等的研究［13］，引物采用真細(xì)菌16S r DNA正向通用引物27 F；反向通用引物1492R．，該對(duì)通用引物可以從廣泛的真細(xì)菌中擴(kuò)增出16S r DNA片段。50μL PCR反應(yīng)體系包括：10×PCR緩沖液5．0μL（含15 mmol／L Mg Cl2），正反向引物（10μmol／L）各1．0μL，d NTPs（各2．5 mmol／L）4．0μL，DNA模板2．5μL，Taq酶1．2 U，加無(wú)菌水至50μL；熱循環(huán)程序?yàn)?5℃5 min，接著進(jìn)行35個(gè)循環(huán)（94℃3 min，55℃30 s，72℃1 min），最后在72℃反應(yīng)10 min后結(jié)束。PCR擴(kuò)增產(chǎn)物為各種真細(xì)菌16s r DNA片段的混合物。

表1 PCR引物信息

1．2．3 擴(kuò)增產(chǎn)物純化及克隆

PCR產(chǎn)物用1%瓊脂糖凝膠電泳分離，在紫外凝膠成像儀下將目的條帶切下，采用瓊脂糖凝膠回收試劑盒回收、純化。純化產(chǎn)物與p MD19－T載體于16℃連接過(guò)夜，轉(zhuǎn)化大腸桿菌DH5α感受態(tài)細(xì)胞，藍(lán)白斑篩選陽(yáng)性克隆。

1．2．4 菌落PCR驗(yàn)證

采用菌落PCR驗(yàn)證白色菌落是否為插入了目的片段的陽(yáng)性克隆。25μL PCR反應(yīng)體系包括：10×PCR緩沖液2．5μL（含15 mmol／L MgCl2），正反向引物（10μmol／L）各 0．2μL，d NTP（各2．5 mmol／L）2．0μL，Taq酶0．6U，加無(wú)菌水至25μL，挑取單個(gè)菌落作為DNA模板加入；熱循環(huán)程序?yàn)椋?5℃5 min，接著進(jìn)行38個(gè)循環(huán)（94℃1 min，55℃1 min，72℃2 min），最后在72℃反應(yīng)10 min后結(jié)束。PCR產(chǎn)物用1．5%的瓊脂糖凝膠電泳分離，EB染色10 min，紫外成像系統(tǒng)檢測(cè)。

1．2．5 限制性內(nèi)切酶分析及測(cè)序

限制性內(nèi)切酶為MspⅠ，反應(yīng)體系為10μL，包括10×緩沖液1μL，PCR產(chǎn)物9μL，內(nèi)切酶4 U，加去離子水至10μL，37℃溫浴過(guò)夜。酶切產(chǎn)物采用2%瓊脂糖冰浴電泳分離。

對(duì)酶切圖譜進(jìn)行鑒定、比較，將酶切圖譜一致的歸為同一個(gè)分類(lèi)操作單元OTU（operational taxonomic unit），挑取同一OTU的代表性克隆子送北京奧科公司進(jìn)行雙向測(cè)序。將不同反應(yīng)所獲得的序列進(jìn)行拼接，獲得16s r DNA片段序列（約1．5 kb）。將拼接序列進(jìn)行BLAST比對(duì)（http：∥www．ncbi．nl m．nih．gov／BLAST）。

2 結(jié)果及分析

2．1 番石榴果實(shí)蠅共生菌總DNA的PCR擴(kuò)增

樣品總DNA的PCR擴(kuò)增見(jiàn)圖1，通過(guò)PCR擴(kuò)增獲得的條帶約為1．5 kb，且無(wú)非特異性條帶。得到的擴(kuò)增產(chǎn)物經(jīng)純化后用于16s r DNA文庫(kù)的建立。

圖1 番石榴果實(shí)蠅共生菌的16s r DNA擴(kuò)增圖譜

2．2 番石榴果實(shí)蠅雄蟲(chóng)共生菌及其分析

將酶切圖譜進(jìn)行比較、分組，共獲得11個(gè)OTU，如圖2a所示，每一個(gè)類(lèi)型選取一個(gè)代表性序列進(jìn)行克隆測(cè)序并參考Kounatidis等人的方法統(tǒng)計(jì)豐度［10］，共 獲 得 11 條 16s r DNA 序 列，長(zhǎng) 度 為1．5 kb左右。通過(guò)BLAST N搜索Gen Bank數(shù)據(jù)庫(kù)，選擇數(shù)據(jù)庫(kù)中相似性最高且至少大于95%的序列（表2）。相似度分別從95%到99%。從表2可以看出，番石榴果實(shí)蠅雄蟲(chóng)共生菌包括假單胞菌屬（Pseudomonas）、腸桿菌屬（Enterobacter）、不動(dòng)桿菌屬 （Acinetobacter）、寡 養(yǎng) 單 胞 菌 屬 （Stenotr ophomonas）、脂環(huán)酸芽胞桿菌屬（Alicyclobacill us）、黃桿菌 科 （Flavobacteriaceae）、乳 球 菌屬 （Lactococcus）、噬酸菌屬（Acidovor ax）、Orbus菌屬以及檸檬酸桿菌屬（Citrobacter）。其中假單胞菌有2種，其他共生菌各1種。在這10個(gè)屬（科）中，乳球菌屬和假單胞菌屬占了50%左右，其中乳球菌約為23．6%，假單胞菌約為26．4%。

2．3 番石榴果實(shí)蠅雌蟲(chóng)共生菌及其分析

將番石榴果實(shí)蠅雌蟲(chóng)共生菌文庫(kù)中的13個(gè)OTU選取代表性序列進(jìn)行測(cè)序，去掉空載體及無(wú)效序列，共得13條16s r DNA序列。通過(guò)BLAST N搜索GenBank數(shù)據(jù)庫(kù)，選擇數(shù)據(jù)庫(kù)中相似性最高的序列（表3），相似度從91%到99%。從表3可以看出，番石榴實(shí)蠅雌蟲(chóng)共生菌屬于假單胞菌屬（Pseudomonas）、腸桿菌屬（Enterobacter）、不動(dòng)桿菌屬（Acinetobacter）、叢毛單胞菌屬（Comamonas）、蒼白桿菌屬（Ochrobactr um）、寡養(yǎng)單胞菌屬（Stenotrophomonas）、脫硫弧菌屬（Desul f ovibrio）、醋桿菌屬（Acetobacter）。其中腸桿菌屬有5種，醋桿菌屬有3種，其他屬的細(xì)菌各1種。腸桿菌屬是優(yōu)勢(shì)共生菌，約占33．9%。

圖2 共生菌文庫(kù)中部分克隆子RFLP結(jié)果

表2 番石榴果實(shí)蠅（雄）廣州種群共生菌組成和相對(duì)豐度1）

表3 番石榴果實(shí)蠅（雌）廣州種群共生菌組成和相對(duì)豐度

3 討論

本研究在番石榴果實(shí)蠅雄蟲(chóng)中鑒定了11種細(xì)菌，隸屬于假單胞菌屬、腸桿菌屬、不動(dòng)桿菌屬、寡養(yǎng)單胞菌屬、脂環(huán)酸芽胞桿菌屬、黃桿菌科、乳球菌屬、噬酸菌屬、Orbus菌屬以及檸檬酸桿菌屬，表明了番石榴果實(shí)蠅中存在多種共生菌。番石榴果實(shí)蠅雄蟲(chóng)體內(nèi)豐度較高的假單胞菌在地中海實(shí)蠅和油橄欖實(shí)蠅中也曾被發(fā)現(xiàn)過(guò)［10，14］，但并非優(yōu)勢(shì)菌。而另外一種豐度較高的乳球菌則未有文獻(xiàn)報(bào)道。另外，墨西哥實(shí)蠅中也曾發(fā)現(xiàn)過(guò)檸檬酸桿菌［15－16］。

本研究在番石榴果實(shí)蠅雌蟲(chóng)中檢測(cè)到了13種細(xì)菌，隸屬于假單胞菌屬、腸桿菌屬、不動(dòng)桿菌屬、叢毛單胞菌屬、蒼白桿菌屬、寡養(yǎng)單胞菌屬、脫硫弧菌屬、醋桿菌屬。其中腸桿菌屬和醋酸桿菌為優(yōu)勢(shì)共生菌。

研究表明初級(jí)共生菌與次級(jí)共生菌相比GC含量較低，這有利于基因的快速進(jìn)化［17］，如尖音庫(kù)蚊（Culex pipiens）感染的Wol bachia，其16s r DNA序列的GC含量為46．6%［18］，阮永明檢測(cè)的B型煙粉虱初級(jí)共生菌16s r DNA序列，GC含量為48．2%，次級(jí)共生菌為54．1%［19］；李正西研究桃蚜初生共生菌，發(fā)現(xiàn)其16s r DNA序列GC含量為49．5%，而次生共生菌為55．5%［20］。本研究鑒定的細(xì)菌序列中，GC含量最低為51．0%，多數(shù)接近自由生活的細(xì)菌。據(jù)此可以推測(cè)本研究鑒定的細(xì)菌均為次級(jí)共生菌。

腸桿菌科在實(shí)蠅共生菌中具有重要地位，Daniela在地中海實(shí)蠅中發(fā)現(xiàn)了克雷伯氏菌和腸桿菌［11］，繞實(shí)蠅屬、按實(shí)蠅屬和果實(shí)蠅屬中都有腸桿菌存在［12，16，21］。

次級(jí)共生菌部分來(lái)自于垂直傳遞，部分來(lái)自于不同寄主的水平傳遞。受環(huán)境影響，昆蟲(chóng)的次級(jí)共生菌不如初級(jí)共生菌穩(wěn)定，容易發(fā)生改變。如腸桿菌科細(xì)菌多存在昆蟲(chóng)腸道中，隨取食而攝取，隨糞便而排出［20］。但這并不能排除與宿主協(xié)同進(jìn)化的共生菌的存在。

本研究在雄蟲(chóng)中發(fā)現(xiàn)豐度較高的綠膿假單胞菌和豐度較低的腸桿菌屬細(xì)菌，而在雌蟲(chóng)中這一結(jié)果正好相反。Behar等在地中海實(shí)蠅中同樣也分離出綠膿假單胞菌P．a(chǎn)er uginosa。該菌具有致病性，能降低宿主昆蟲(chóng)的壽命［14］。腸桿菌科細(xì)菌具有拮抗致病菌、降解有毒物質(zhì)、增加宿主適應(yīng)性的功能，另外，腸桿菌還能夠促進(jìn)宿主的碳氮代謝循環(huán)［22］。假單胞菌和腸桿菌互為消長(zhǎng)，這與兩類(lèi)細(xì)菌的生理學(xué)功能相吻合。食物對(duì)昆蟲(chóng)腸道的微生物影響比較復(fù)雜，不同的食物可能會(huì)影響共生菌的組成［23］。有人認(rèn)為在野生狀態(tài)下，昆蟲(chóng)攝取的果實(shí)中含有的次生代謝物質(zhì)會(huì)具有一定的殺菌作用，可能會(huì)影響腸道共生菌，尤其傾向選擇某些有解毒功能的細(xì)菌［24］。蟑螂在取食低蛋白、高纖維食物時(shí)，前腸的鏈球菌（Streptococci）和乳桿菌（Lactobacilli）數(shù)量減少，高蛋白的食物則導(dǎo)致H2和CO2產(chǎn)量減少以及G＋C細(xì)菌增多［25］。因此筆者估計(jì)飼養(yǎng)種群與野外種群的共生菌多樣性可能存在差異。

目前針對(duì)昆蟲(chóng)共生菌的分離鑒定研究方法主要有3種，即傳統(tǒng)的分離培養(yǎng)鑒定法、顯微觀察法和以變性梯度凝膠電泳（PCR－DGGE）和16s r DNA文庫(kù)構(gòu)建為主的分子生物學(xué)方法。相比傳統(tǒng)的微生物學(xué)方法，分子生物學(xué)方法能全面獲得共生菌的信息，而且不受細(xì)菌是否可培養(yǎng)的限制。本研究采用的16s r DNA文庫(kù)構(gòu)建和測(cè)序方法，不僅獲得序列信息比PCR－DGGE更為豐富，還能評(píng)估各細(xì)菌的豐度，且不需要昂貴的儀器設(shè)備［26］。

研究結(jié)果為探索番石榴果實(shí)蠅共生菌的生理功能及其與宿主的協(xié)同進(jìn)化關(guān)系提供了基礎(chǔ)。在明確共生菌多樣性的基礎(chǔ)上，結(jié)合相關(guān)文獻(xiàn)報(bào)道，選取對(duì)番石榴果實(shí)蠅入侵起正作用的共生菌開(kāi)展研究。針對(duì)可培養(yǎng)細(xì)菌，通過(guò)高溫或抗生素去除共生菌，然后通過(guò)飼喂或顯微注射導(dǎo)入，通過(guò)生物學(xué)試驗(yàn)明確共生菌對(duì)寄主的生理功能，可以明確共生菌對(duì)番石榴果實(shí)蠅的生物學(xué)功能，并進(jìn)一步揭示番石榴果實(shí)蠅和共生菌協(xié)同、快速入侵的機(jī)制。

［1］ 梁廣勤，楊國(guó)海，梁帆，等．亞太地區(qū)寡毛實(shí)蠅［M］．廣州：廣東科技出版社，1996．

［2］ 汪興鑒，趙明珠．中國(guó)寡鬃實(shí)蠅屬記述（雙翅目：實(shí)蠅科）［J］．動(dòng)物分類(lèi)學(xué)報(bào)，1989（2）：209－219．

［3］ Foote R H，Blanc F L，Norr bo m A L．Handbook of the fr uit flies（Diptera：Tephritidae）of America Nort h of Mexico［M］．USA：Co mstock，1993．

［4］ 譚周進(jìn)，肖啟明，謝丙炎，等．昆蟲(chóng)內(nèi)共生菌研究概況［J］．微生物學(xué)通報(bào)，2005，32（4）：140－143．

［5］ Ker kut G A，Gilbert L I．Comprehensive insect physiology，biochemistr y and phar macology．Vol．41［M］．El msfor d（NY）：Pergamon Press，1985．

［6］ Lauzon C R，Potter S E，Pr okopy R J．Degradation and detoxification of the dihydr ochalcone phloridzin by Enter obacter agglomer ans，a bacteriu m associated with the apple pest，Rhago－letis pomonell a （Walsh）（Diptera：Tephritidae）［J］．Environmental Entomology，2003，32（5）：953－962．

［7］ Sun X，Cui L，Li Z，Diversity and phylogeny of Wolbachia infecting Bactrocer a dorsalis（Diptera：Tephritidae）populations from China［J］．Environ mental Entomology，2007，36（5）：1283－1289．

［8］ 褚棟，劉國(guó)霞，陶云荔，等．煙粉虱復(fù)合種內(nèi)共生菌多樣性及其生物學(xué)意義［J］．昆蟲(chóng)學(xué)報(bào)，2006，49（4）：687－694．

［9］ Br u mmel T，Ching A，Ser oude L，et al．Dr osophil a lifespan enhancement by exogenous bacteria［J］．PNAS，2004，101（35）：12974－12979．

［10］Kounatidis I，Crotti E，Sapountzis P，et al．Acetobacter tropicalis is a major sy mbiontof the olive fr uit fly（Bactrocer a oleae）［J］．Applied and Envir on mental Micr obiology，2009，75（10）：3281－3288．

［11］Daniela M，Marco R，Ro mano D，et al．Bacteria associated with the esophageal bulb of the medfly Ceratitis capitata （Diptera：Tephritidae）［J］．Current Microbiology，2002，44：120－124．

［12］Capuzzo C，F(xiàn)irrao G，Mazzon L，et al．“Candidatus Er winia dacicola”，a coevolved sy mbiotic bacteriu m of the olive fly Bactrocera oleae （Gmelin）［J］．International Journal of Systematic and Evolutionar y Micr obiology，2005，55（4）：1641－1647．

［13］Weisburg W G，Barns S M，Pelletier D A，et al．16S ribosomal DNA a mplification for phylogenetic study［J］．Jour nal of Bacteriology，1991（173）：697．

［14］Behar A，Yuval B，Jur kevitch E．Gut bacterial communities in the Mediterranean fr uit fly（Cer atitis capitata）and their i mpacton host longevity［J］．Journal of Insect Physiology，2008，54（9）：1377－1383．

［15］Robacker D C，Bartelt R J．Chemicals attractive to Mexican fr uit fly fro m Klebsiell a pneu moniae and Citr obacter f reundii cult ures sampled by solid－phase microextraction［J］．Jour nal of Chemical Ecology，1997，23（12）：2897－2915．

［16］Kuzina L V，Peloquin J J，Vacek D C，et al．Isolation and identification of bacteria associated with adult laboratory Mexican fr uit flies，Anastrepha l udens（Diptera：Tephritidae）［J］．Current Microbiology，2001，42（4）：290－294．

［17］Nancy A，Paul B．Bacterial endosy mbionts in ani mals［J］．Ecology and Industrial Microbiology，2000，3：270－275．

［18］O’Neill S L，Giordano R，Colbert A M，et al．16Sr RNA phylogenetic analysis of the bacterial endosy mbionts associated with cytoplas mic inco mpatibility in insects［J］．PNAS，1992，89（7）：2699－2702．

［19］阮永明，劉樹(shù)生．浙江B型與非B型（China－ZHJ－1）煙粉虱種群共生細(xì)菌的檢測(cè)及系統(tǒng)發(fā)育分析［J］．昆蟲(chóng)學(xué)報(bào)，2005，48（6）：859－865．

［20］李正西，李定旭．桃蚜自然種群初級(jí)和次級(jí)共生菌的分子鑒定［J］．昆蟲(chóng)學(xué)報(bào)，2005，48（5）：810－814．

［21］Howard D J，Bush G L．Infl uence of bacteria on lar val sur vival and develop ment in Rhagoletis（Diptera：Tephritidae）［J］．Annals of the Entomological Society of America，1989，82：633－640．

［22］Yuichi H，Tetsuhiko S，Haji me I．Cloning，sequence analysis and expression in Escherichia coli if the gene encoding a uricase from the yeast－like sy mbionts of bro wn plant hoppers，Nil apar vata l ugens［J］．Insect Biochemistr y and Molecular Biology，2000，30：173－182．

［23］相輝，李木旺，趙勇，等．家蠶幼蟲(chóng)中腸細(xì)菌群落多樣性的PCR－DGGE和16s r DNA文庫(kù)序列分析［J］．昆蟲(chóng)學(xué)報(bào)，2007，50（3）：222－233．

［24］相輝，黃勇平．腸道微生物與昆蟲(chóng)的共生關(guān)系［J］．昆蟲(chóng)知識(shí)，2008，45（5）：687－693．

［25］Santo Domingo J W，Kauf man M G，Klug M J，et al．Infl uence of dieton the str uct ure and f unction of the bacterial hindgut co mmunity of crickets［J］．Molecular Ecology，1998（7）：761－767．

［26］蔣燕，鮑寶龍，謝彩霞，等．16S r DNA文庫(kù)法分析草魚(yú)肝胰臟和膽汁細(xì)菌群落組成［J］．上海海洋大學(xué)學(xué)報(bào)，2010，19（4）：440－446．