越冬白蠟蟲微生物多樣性分析

孫 濤，王雪慶，趙遵嶺，于淑慧，楊 璞

(中國林業科學研究院資源昆蟲研究所，國家林業局資源昆蟲培育與利用重點實驗室，云南 昆明 650224)

越冬白蠟蟲微生物多樣性分析

孫 濤，王雪慶，趙遵嶺，于淑慧，楊 璞*

(中國林業科學研究院資源昆蟲研究所，國家林業局資源昆蟲培育與利用重點實驗室，云南 昆明 650224)

目的了解白蠟蟲越冬時體內微生物的多樣性，比較昆明與長春越冬蟲體內微生物種類和數量的差異，為了解白蠟蟲低溫適應機制提供有用信息。方法采用Miseq高通量測序方法對昆明越冬雌成蟲(KM)和長春越冬雌成蟲(CC)體內細菌16s rRNA及真菌ITS基因進行測序，利用Usearch 軟件進行細菌和真菌OTU(Operational Taxonomic Unit)劃分，并利用Mothur軟件對OTU進行分類學分析和多樣性指數分析。結果細菌KM和CC樣品中共得到344個OTU，真菌KM和CC樣品中共得到230個OTU。細菌共鑒定到15門、137屬，在KM中乳球菌屬占主要比例(29.78%)，而在CC中立克次氏體占主要比例(55.77%)，真菌中共鑒定到6門、83屬，在KM中僅鑒定到41個屬，而CC中鑒定到75個屬，其中線蟲草科無法歸類的屬在2個樣品中均為優勢菌群，但CC中含量遠低于KM。結論昆明與長春越冬蟲體內細菌及真菌種類和數量均不相同，與昆明越冬蟲不同，需要應對極端低溫的長春越冬蟲體內立克次氏體為優勢菌群。

白蠟蟲；群落組成；越冬蟲；Miseq測序；立克次氏體

白蠟蟲(EriceruspelaChavannes)是我國傳統的資源昆蟲，其2齡雄幼蟲分泌的白蠟具有重要的經濟價值，被廣泛應用于醫藥、食品、航天等領域[1-2]。白蠟蟲分布廣泛，在我國其自然種群主要分布在18°N到42°N地區[3]。李健等在最冷月日均最低氣溫-18.0℃的吉林省四平市發現白蠟蟲依然能夠生長發育繁衍后代[4]，說明白蠟蟲對極端低溫環境具有很強的適應能力。白蠟蟲以雌成蟲蟲態越冬，對白蠟蟲越冬過程中生理狀態等方面的研究將有助于了解其低溫適應機制，并對擴大白蠟產區具有參考意義。

共生微生物在昆蟲生長發育、繁衍進化以及環境適應等方面具有重要作用。例如，在對美洲大蠊(PeriplanetaAmericanaL.)內共生菌的研究中發現其共生菌中含有可抑制病原細菌生長的共生菌，并且存在著能夠降解多糖的真菌[5]。研究發現豌豆蚜(AcyrthosiphonpisumHarris)的共生菌布赫納氏菌屬(Buchnera)能夠幫助其抵御極端溫度[6]。研究表明共生真菌也具有抵御極端溫度侵害的作用，如在對德州切葉蟻(AttatexanaBuckley)進行研究時發現，其所攜帶的Attamyces共生菌株能夠幫助切葉蟻抵御低溫侵害提高適應性[7]。在對白蠟蟲低溫適應機制的研究中，發現長春和昆明越冬蟲在代謝等方面存在很多差異[8-9]，我們推測白蠟蟲內共生微生物也可能在越冬以及抵御極端低溫過程中具有重要作用。

本研究將對白蠟蟲越冬蟲態體內共生微生物進行16s rRNA和ITS基因高通量測序，研究白蠟蟲越冬過程體內微生物多樣性，并將昆明與長春越冬雌成蟲體內微生物進行比較分析，以了解在極端低溫條件下越冬雌成蟲體內微生物種類及數量的變化，分析微生物和白蠟蟲越冬及抗凍可能的相關性。

1 材料和方法

1.1 試驗材料

所用白蠟蟲分別為昆明和長春越冬雌成蟲。昆明越冬雌成蟲(KM)采自云南昆明(102.73°E；25.04°N)；長春越冬雌成蟲(CC)采自于吉林省長春市(125.32°E；43.90°N)。采集時將白蠟蟲從女貞樹上迅速剝離放入液氮中儲存。

1.2 越冬雌成蟲內微生物總DNA提取

按照DNeasy Blood & Tissue Kit試劑盒(QIAGEN, Germany)方法提取越冬雌成蟲內微生物總DNA，將提取出的微生物總DNA進行1%的瓊脂糖凝膠電泳檢測，并用超微量紫外分光光度計(Thermo, 美國)測定DNA濃度。將抽提后的總DNA置于-20℃保存。

1.3 PCR擴增白蠟蟲越冬雌成蟲微生物的16S和ITS基因

以提取的越冬雌成蟲內微生物總DNA為模板，使用通用引物338F和806R擴增細菌16s rRNA基因片段[10]，使用通用引物1737 F和2043 R擴增真菌ITS基因片段[11]。PCR擴增體系如下：總DNA 0.4 μL、上下游引物各0.4 μL、5×FastPfu Buffer 4 μL、2.5 mmol·L-1dNTPs 2 μL、FastPfu Polymerase(TaKaRa, 日本) 0.4 μL以及雙蒸無菌水10.3 μL。擴增程序為：94℃預變性4 min；94℃變性30 s，55℃退火30 s，72℃延伸45 s，28個循環；72℃延伸10 min。

將PCR產物進行1%的瓊脂糖凝膠電泳檢測，使用Axygen DNA膠回收試劑盒(Axygen, 美國)進行回收?；厥债a物使用QuantiFluorTM-ST定量系統(Promega, 美國)定量，按照上機測序要求跟據測序量計算樣本上機量。

1.4 Miseq測序

首先將“Y”字形接頭連接到片段上，通過磁珠篩選將接頭自連片段剔除，然后利用PCR技術擴增已加入“Y”字形接頭的目的片段，經氫氧化鈉變性后得到單鏈DNA片段上機測序文庫。使用Illumina MiSeq(Illumina, 美國)測序儀對獲得的上機測序文庫進行測序分析。

1.5 Miseq測序結果數據分析

根據reads的重疊關系，使用FLASH軟件(Version 3+；http://gnu.org/licenses/gpl.html)將Miseq雙端測序得到的序列拼接得到一條序列，同時對reads質量以及拼接效果進行質控和有效數據的過濾，將接頭和無重復單序列剔除獲得有效序列。利用Usearch 軟件 (Vsesion 7.1；http://drive5.com/uparse/)在97%相似度水平下對有效序列進行OTU(Operational Taxonomic Unit)劃分，并將聚類過程中產生的嵌合體剔除。

根據OTU劃分結果，利用RDP Classifier軟件(Version2.2；http://sourceforge.net/projects/rdp-classifier/)進行分類學分析，在域、界、門、綱、目、科、屬、種等不同分類階元統計其群落組成。利用Mothur軟件(Version v.1.30.1；http://www.mothur.org/wiki/Schloss_SOP#Alpha_diversity)進行Chao(Chao指數)、Ace(Ace指數)、Shannon(香農指數)以及Simpson(辛普森指數)等多樣性指數分析，其中，Ace指數和Chao指數代表了群落的豐富程度，香農指數和辛普森指數表明了群落多樣性，同時辛普森指數也可以呈現優勢種的集中程度。

利用Mothur軟件 (Version v.1.30.1；http://www.mothur.org/wiki) 和R 軟件進行稀釋曲線(rarefaction curve)和Shannon-Wiener曲線(Shannon-Wienercurve)分析，其中，稀釋曲線用于比較測序數量不同的樣本中物種的豐富程度，而Shannon-Wiener曲線則用于反映各樣本在不同測序量時的微生物多樣性。除此之外兩種曲線均被用于說明樣本測序數量的合理性。

2 結果

2.1 OTU數目統計

細菌16S rRNA：KM內共檢測到12 538條reads，CC內共檢測到27 575條reads；真菌ITS：KM內共檢測到14 483條reads， CC內共檢測到21 650條reads。KM和CC兩個樣品細菌一共檢測到344個OTU：包含KM 177個和CC 167個；真菌共有230個OTU：包含KM 168個和CC 62個。

KM和CC的細菌、真菌稀釋曲線均趨于平坦(圖1：A，B)，繼續增加測序量將不會產生更多的OTU，說明數據量合理可以用于其他分析。同時，Shannon-Wiener曲線也均趨于平坦(圖1：C，D)，進一步說明數據量足夠大，可以反映樣本中絕大部分微生物信息。

圖1 KM和CC稀釋曲線及Shannon-Wiener曲線Fig. 1 Rarefaction curve and Shannon-Wiener curve of KM and CCA：細菌稀釋曲線Rarefaction curve of bacteria；B：真菌稀釋曲線Rarefaction curve of bacteria fungi；C：細菌Shannon-Wiener曲線Shannon-Wiener curve of bacteria；D：真菌Shannon-Wiener曲線Shannon-Wiener curve of fungi

2.2 昆明長春越冬雌成蟲內細菌種類及變化

基于Miseq測序結果，在細菌中共鑒定到15門28綱54目84科137屬。

鑒定到的門水平中主要有：變形菌門、厚壁菌門、擬桿菌門、放線菌門，這4個門的總和在KM和CC中分別占98.29%和98.86%，是越冬雌成蟲內細菌的主要成分，其中變形菌門在越冬雌成蟲內細菌中占主要比例。將KM和CC中這4個門所占比例進行對比后發現，除了變形菌門從41.57%(KM)升高至72.25%(CC)，其它3個門在CC中所占比例均下降(圖2:A)。

在KM中共鑒定到115個屬細菌，其中乳球菌屬(Lactococcus)占主要比例(29.78%)，其次為假單胞菌屬(Pseudomonas)占15.05%、鞘氨醇單胞菌屬(Sphingomonas)占5.38%、薄層菌屬(Hymenobacter)占5.31%。而在CC中共鑒定到125個屬細菌，其中立克次氏體(Rickettsia)占主要比例(55.77%)，其次為乳球菌屬(Lactococcus)占17.96%、假單胞菌屬(Pseudomonas)占8.09%。對比KM與CC發現，占主要比例的菌屬發生了改變，KM中的大多數較高豐度的菌屬在CC中所占比例較低(圖2:B)，立克次氏體在KM中僅占0.03%。

2.3 昆明長春越冬雌成蟲內真菌種類及變化

基于Miseq測序結果，在真菌中共鑒定到6門20綱35目46科83屬。

鑒定到的門水平中主要有：子囊菌門、擔子菌門，這2個門的總和在KM和CC中分別占90.42%和99.37%，其中子囊菌門在越冬雌成蟲內真菌中占比最高，為78.45% (圖2:C)。

在KM中僅鑒定到41個屬真菌，其中線蟲草科無法歸類的屬(unclassified Ophiocordycipitaceae)占主要比例(31.06%)，其次為微座孢屬(Microstroma)占13.96%、紅酵母屬(Rhodotorula)占13.14%、莖點霉屬(Phoma)占10.38%、橫斷孢屬(Strelitziana)占8.45%、枝孢霉屬(Cladosporium)占5.26%。而在CC中共鑒定到75個屬真菌，其中線蟲草科無法歸類的屬(unclassified Ophiocordycipitaceae)占主要比例(14.75%)，其次為枝氯霉屬(Ramichloridium)占15.76%、枝孢霉屬(Cladosporium)占9.11%、莖點霉屬(Phoma)占9.03%。對比KM與CC發現，線蟲草科無法歸類的屬、微座孢屬、紅酵母屬等菌屬在CC所占比例低于KM，但CC中枝孢霉屬所占比例高于KM，而且新鑒定到了諸如枝氯霉屬等34個菌屬(圖2：D)。

圖2 KM和CC微生物群落組成 此圖為文章證據材料，根據占比多少不同已在文中有說明Fig. 2 Composition of microorganism communities in the KM and CCA：KM和CC細菌門水平群落組成Composition of bacterial communities at the phylum；B：KM和CC細菌屬水平群落組成Composition of bacterial communities at the genus；C：KM和CC真菌門水平群落組成Composition of fungal communities at the phylum；D：KM和CC真菌屬水平群落組成Composition of fungal communities at the genus

2.4 昆明長春越冬雌成蟲內細菌和真菌多樣性分析

在兩細菌樣本中，CC的Ace指數和Chao指數均稍高于KM，表明在CC細菌群落中其群落豐富度相對高于KM；CC的香農指數低于KM，而辛普森指數高于KM，表明其群落多樣性相對低于KM，辛普森指數高于KM同時也反映了CC優勢種的集中程度要高于KM(表1)。

在兩真菌樣品中，CC的Ace指數、Chao指數和香農指數均高于KM，表明在CC真菌群落中其群落豐富度以及多樣性相對高于KM，而辛普森指數低于KM，說明CC優勢種的集中程度要低于KM(表1)。

表1 KM和CC內細菌和真菌多樣性指數分析

3 討論

昆蟲在生理上主要是通過調節含水量、控制冰核物質、改變代謝途徑以及積累抗凍蛋白等耐寒物質來提高自身適應能力降低低溫傷害[12]。許多研究發現，微生物在昆蟲適應低溫等不良環境脅迫時起到了重要作用[6-7,12-13]。前期在對白蠟蟲耐寒性的研究中，已經鑒定到了抗凍蛋白候選基因并對其中部分基因進行了功能驗證[8]。本研究基于Miseq高通量測序技術對白蠟蟲越冬雌成蟲體內微生物進行研究。

對比昆明長春白蠟蟲越冬雌成蟲體內微生物種類和數量發現，無論是細菌還是真菌都發生了很大的變化，群落豐富度相對提高，尤其是真菌樣品，研究發現有34個屬在CC中被鑒定到而KM中沒有。

在昆明長春白蠟蟲越冬雌成蟲體內均鑒定到了立克次氏體，并且其在長春樣品CC中所占比例高達55.77%。立克次氏體能夠廣泛侵染植物、脊椎動物以及無脊椎動物等。立克次氏體在昆蟲營養、對環境適應性等方面起到了一定的積極作用[14]。在Chen等2000年的研究中發現立克次氏體能夠在不適宜溫度條件下提高豌豆蚜繁殖能力[15]，在對煙粉虱(BemisiatabaciGennadius)中立克次氏體的研究中發現立克次氏體能夠提高寄主對于極端溫度的耐受能力[16]。在對白蠟蟲立克次氏體的研究中發現 ，自然環境下并不是所有白蠟蟲都被立克次氏體寄生[17]，但是本研究的樣品KM和CC中均檢測到立克次氏體的存在，而且長春樣品CC中所占比例很高，因此推測可能是和立克次氏體能夠提高寄主對環境適應性的特點有關。

假單胞菌屬是目前國內外研究發現的具有生物冰核活性的細菌屬，其可以通過誘導結冰降低寄主體溫[18]。研究發現，一些昆蟲能夠在其受到低溫脅迫時將具有生物冰核活性的物質排出體外，降低自身過冷卻點，抵御極端低溫的侵害[19]。從昆明長春白蠟蟲越冬雌成蟲體內微生物高通量測序結果中，我們發現相比昆明越冬雌成蟲樣品KM，假單胞菌屬在CC中所占比例降低，可能是白蠟蟲為了適應更低的溫度環境，很可能將這些具有生物冰核活性的假單胞菌從自身排出，以降低過冷卻點。

除了寄主體內細菌具有提高寄主對于不良環境的耐受能力，有的真菌也具有這種特征[7]。在對松樹蜂(SirexnoctilioFabricius)與其共生真菌關系的研究中發現，其共生真菌能夠協助松樹蜂卵抵御不良環境保障其成功孵化[20]?；诒狙芯扛咄繙y序結果，有34個真菌屬只在CC中出現，可能是受KM和CC所處地理位置以及環境條件不同影響越冬雌成蟲生理和代謝狀態改變所致。對于KM和CC內真菌的差異是否與極端低溫環境有關，還需要進一步研究。

4 結論

基于Miseq高通量測序結果，我們發現長春和昆明越冬雌成蟲體內微生物存在顯著差異，在對長春越冬雌成蟲微生物檢測中發現，有很多能夠協助昆蟲適應不利環境因素的微生物種群相比昆明樣品變化很大，諸如立克次氏體、假單胞菌屬等。

[1] Yang P, Chen X M, Liu W W,etal.Transcriptome analysis of sexually dimorphic Chinese white wax scale insects reveals key differences in developmental programs and transcription factor expression[J].Scientific Reports,2015 (5):8141.

[2] 孫 濤,王雪慶,趙遵嶺,等.白蠟蟲熱激蛋白序列分析[J].林業科學研究, 2016,29(5):778-783.

[3] 張長海.白蠟蟲在我國的地理分布[J].林業科學研究, 1991, 4(2):192-196.

[4] 李 健,張喜儉,程玉林,等.白蠟蟲生物學特性及防治對策[J].吉林農業:學術版,2010(8):170-170.

[5] 阮傳清,劉 波,吳珍泉,等.美洲大蠊內生菌的分離鑒定及酶活性、抑菌作用測定[J].微生物學雜志,2015，35(5):36-42.

[6] Dunbar H E,Wilson A C,Ferguson N R,etal.Aphid thermal tolerance is governed by a point mutation in bacterial symbionts[J] .Plos Biology，2007,5(5)：e96

[7] Mueller U G,Mikheyev A S,Hong E,etal.Evolution of cold-tolerant fungal symbionts permits winter fungiculture by leafcutter ants at the northern frontier of a tropical ant-fungus symbiosis[J].Proceedings of the National Academy of Sciences,2011,108(10):4053-4056.

[8] 楊 璞,陳曉鳴,謝正華,等.白蠟蟲不同蟲態及不同地理種群越冬雌成蟲過冷卻點的研究[J].林業科學研究,2011,24(5):585-590.

[9] Yu S H,Yang Pu,Sun T,etal.Transcriptomic and proteomic analyses on the supercooling ability and mining of antifreeze proteins of the Chinese white wax scale insect[J].Insect Science,2016,23(3):430-437.

[10] Xu N,Tan G,Wang H,etal.Effect of biochar additions to soil on nitrogen leaching,microbial biomass and bacterial community structure[J].European Journal of Soil Biology,2016,74:1-8.

[11] Bokulich N A,Mills D A.Improved internal transcribed spacer (ITS) primer selection enables quantitative, ultra-high-throughput fungal community profiling[J].Applied & Environmental Microbiology,2013,79(8):2519-2526.

[12] 馮宇倩,王錦林,宗世祥.昆蟲越冬蟲態及耐寒策略概述[J].中國農學通報,2014,30(9):22-25.

[13] 馬巨法,程家安.冰核活性微生物對昆蟲耐寒性的影響[J].植物保護,1996，22(2):35-38.

[14] 梁士可,王方海.稻飛虱共生菌的研究進展[C]. 第十七次全國環境微生物學學術研討會.2014.

[15] Chen D Q,Montllor C B,Purcell A H.Fitness effects of two facultative endosymbiotic bacteria on the pea aphid,Acyrthosiphon pisum,and the blue alfalfa aphid,A. kondoi[J].Entomologia Experimentalis Et Applicata,2000,95(3):315-323.

[16] Brumin M,Kontsedalov S,Ghanim M. Rickettsia,influences thermotolerance in the whitefly Bemisia tabaci, B biotype[J].Insect Science,2011,18(1):57-66.

[17] 劉魏魏,楊 璞,阮永明,等.白蠟蟲體內殺雄菌屬(Arsenophonus)共生菌的分子檢測[J].微生物學報,2012,52(8):1002-1010.

[18] 杜鵬瑞,杜 睿,路則棟,等.西藏那曲地區夏季降水中生物冰核的分布特征[J].中國環境科學,2015，35(7):1998-2006.

[19] Zachariassen&Amp K E,Husby J A.Antifreeze effect of thermal hysteresis agents protects highly supercooled insects[J]. Nature,1982,298(5877):865-867.

[20]李大鵬,石 娟,駱有慶.松樹蜂與其共生真菌的互利共生關系[J].昆蟲學報,2015,58(9):1019-1029.

MicroorganismDiversityAnalysisofOverwinteringChineseWhiteWaxScaleInsect

SUNTao,WANGXue-qing,ZHAOZun-ling,YUShu-hui,YANGPu

(Research Institute of Resources Insects, Chinese Academy of Forestry, Key Laboratory of Cultivating and Utilization of Resources Insects of State Forestry Administration, Kunming 650224, Yunnan, China)

ObjectiveThe study aimed to investigate the microorganism diversity and compare the difference of microorganism between overwintering individuals in Kunming and Changchun.MethodThe bacteria 16S rRNA and the fungi internal transcribed space (ITS) of overwintering individuals in Kunming (KM) and Changchun (CC) were sequencing by MiSeq high-throughput sequenced method. The bacterial and fungal Operational Taxonomic Units (OTUs) were obtained by the Usearch software, and the Mothur software was used to calculate and analyze the taxonomy and alpha diversity.ResultAccording to 344 OTUs obtained from KM and CC bacteria samples, 15 phyla and 137 genera were identified. Meanwhile, 6 phyla and 83 genera were identified from 230 OTUs obtained from KM and CC fungi samples. In KM the dominant bacteria,Lactococcusaccounted for 29.78%. And in CCRickettsiawas dominant (55.77%). UnclassifiedOphiocordycipitaceaewas dominant both in KM and CC, but its proportion in CC was much lower than in KM.ConclusionIn this study, lots of differences were found between Kunming and Changchun overwintering individuals on the bacteria and fungi category and quantity. In comparison of KM bacteria sample,Rickettsiabecame the dominant genus in CC.

Ericeruspela; community composition; overwintering individuals; Miseq sequencing;rickettsia

10.13275/j.cnki.lykxyj.2017.06.018

2016-10-07

國家863計劃2014AA021801，云南省應用基礎研究重點項目2013FA052，林業公益性行業科研專項201504302、201304808，國家自然科學基金31572337.

孫 濤(1991—)，男，碩士研究生，主要從事昆蟲分子生物學研究.

* 通訊作者：楊 璞.

S899.1

A

1001-1498(2017)06-1009-06

張 玲)