基于16S rDNA測(cè)序技術(shù)分析不同齡期蠶沙中菌群多樣性

鄭天瑤 宿樹蘭 蔡紅蝶

[摘要]該文對(duì)不同齡期烘干前后蠶沙中菌群多樣性進(jìn)行分析與評(píng)價(jià)，以期明確其菌群組成及其豐富度差異，闡明烘干處理對(duì)中藥材蠶沙中菌群的影響，為其功效科學(xué)內(nèi)涵的揭示提供依據(jù)和參考。應(yīng)用高通量測(cè)序技術(shù)測(cè)定蠶沙細(xì)菌的16S rDNAV4 變異區(qū)序列，應(yīng)用Qiime，Mothur，PICRUSt等軟件整理和統(tǒng)計(jì)樣品序列數(shù)目和操作分類單元（OTUs）數(shù)量，分析蠶沙樣品中菌群的組成、豐度、分布、Alpha多樣性、Beta多樣性，菌群差異性及對(duì)群落的代謝功能等進(jìn)行預(yù)測(cè)。該研究獲得用于分析的有效序列數(shù)為259 250；稀疏曲線表明測(cè)序深度充分，OTU 的數(shù)量接近于飽和。蠶沙菌群主要由變形菌門（Proteobacteria，893%）、放線菌門（Actinobacteria，50%）、厚壁菌門（Firmicutes，44%）、擬桿菌門（Bacteroidetes，11%）、藍(lán)細(xì)菌（Cyanobacteria，02%）組成，其優(yōu)勢(shì)菌門為變形菌門。烘干后蠶沙的菌群多樣性及豐富度降低，五齡蠶沙尤為明顯。PICRUSt分析蠶沙菌群對(duì)應(yīng)的基因功能發(fā)現(xiàn)膜轉(zhuǎn)運(yùn)、碳水化合物代謝、氨基酸代謝、脂質(zhì)代謝、核苷酸代謝，細(xì)胞過程與信號(hào)傳導(dǎo)等基因功能豐度較高。烘干處理可顯著降低蠶沙中致病菌的種類和數(shù)目，有利于藥材藥用品質(zhì)，相比于低齡期蠶沙，五齡蠶沙更適宜入藥。Illumina MiSeq 高通量測(cè)序技術(shù)為蠶沙菌群的研究提供了更加準(zhǔn)確、科學(xué)的數(shù)據(jù)資源。

[關(guān)鍵詞]不同齡期； 蠶沙； 菌群多樣性； 16S rDNA高通量測(cè)序

Silkworm excrement bacterial communities diversity in different instars

based on 16S rDNA sequence analysis

ZHENG Tianyao1，2， SU Shulan2*， CAI Hongdie2， DAI Xinxin2， OUYANG Zhen1， DUAN Jinao2*

（1 Jiangsu University， Zhenjiang 212013， China；

2 Jiangsu Collaborative Innovation Center of Chinese Medicinal Resources Industrialization， National and

Local Collaborative Engineering Center of Chinese Medicinal Resources Industrialization and Formulae

Innovative Medicine， Nanjing University of Chinese Medicine， Nanjing 210023， China）

[Abstract]This paper investigated the diversity of the silkworm excrement bacterial communities in different ages before and after drying， aiming to clarify the differences of bacterial communities in composition and bacterial abundance and the influences of drying treatment， and provide scientific basis for the efficacy of scientific connotation and utilization of silkworm excrement Highthroughput sequencing technique was used to measure the sequence of 16S rDNAV4 variable region of bacteria in silkworm excrement QIIME， Mothur and PICRUSt software programs were employed to sort and calculate the number of sequences and operational taxonomic units （OTUs） for each sample Thereafter， the abundance， distribution， alpha diversity index of species， beta diversity and bacterial communities diversity among different sample groups and predicted the bacterial gene functions were analyzed In this study， the numbers of effective sequences for six samples were 259 250； the rarefaction curves showed a sufficient sequencing depth， and the number of OTUs was close to saturation The bacteria in silkworm excrement belonged to the following five phylums： Proteobacteria （893%）， Actinobacteria （50%）， Firmicutes （44%）， Bacteroidetes （11%） and Cyanobacteria （02%） The dominant specie was Cyanobacteria of the total bacteria identified， respectively The abundances and diversities of the silkworm excrement bacterial communities have been reduced after drying treatment， especially the silkworm excrement of the fifth instar PICRUSt analysis was performed to show that abundance of the functional genes such as membrane transport， carbohydrate metabolism， amino acid metabolism， lipid metabolism， nucleotide metabolism， cellular processes and signaling were relatively high The result showed that the drying treatment could decreased the species and numbers of pathogenic bacteria in silkworm excrement obviously and improve the quality of medicinal materials Compared with the lower ages， silkworm excrement of fifth instar seems like to be more suitable for use in medicine Illumina MiSeq highthroughput sequencing system provides a more accurate and scientific data resource for the study of bacteria in silkworm excrementendprint

[Key words]different instars； silkworm excrement； bacterial communities diversity； 16SrDNA highthroughput sequencing

蠶沙，即蠶糞，為蠶蛾科昆蟲家蠶Ombyx mori L幼蟲的干燥糞便，是我國的傳統(tǒng)中藥，在《名醫(yī)別錄》和《本草綱目》均有記載，具有祛風(fēng)除濕、清熱明目、活血定痛的功效。用于主治皮膚不仁、關(guān)節(jié)不遂、肢體麻木、風(fēng)疹瘙癢等癥。隨著近年來研究的不斷深入，關(guān)于蠶沙的化學(xué)成分、藥理作用、臨床應(yīng)用的研究日益增多。現(xiàn)代研究表明，蠶沙具有改善和治療貧血[1]、抗炎、鎮(zhèn)痛、抑菌[2]、保肝作用[3]等多種藥理活性，且對(duì)糖尿病患者血糖的改善也有一定療效[4]。以蠶沙為原料提取的葉綠素、果膠、葉黃素、類胡蘿卜素、植物醇、葉蛋白等，廣泛用于食品、飲料、化妝品、醫(yī)藥等領(lǐng)域[5]，而關(guān)于蠶沙的微生物學(xué)研究較少。鐘楊生等[6]研究了蠶沙堆肥過程中微生物群落數(shù)量的動(dòng)態(tài)變化，史才娟等[7]研究了不同齡期蠶沙中細(xì)菌、真菌、放線菌的數(shù)量變化，但關(guān)于蠶沙中菌群的具體組成目前尚不清楚。

隨著分子生物技術(shù)的發(fā)展，一些非培養(yǎng)型技術(shù)被用于腸道菌群、土壤菌群、植物內(nèi)生菌等微生物群落結(jié)構(gòu)的研究中，例如以16S rDNA為基礎(chǔ)的DGGE[8]，Real timePCR[9]，F(xiàn)LSH[10]。近年來高通量測(cè)序技術(shù)已較為廣泛地應(yīng)用在微生物的研究中，是目前最為先進(jìn)的檢測(cè)技術(shù)，具備通量高、測(cè)序快、準(zhǔn)確度高等特點(diǎn)。

本研究采用Illumina MiSeq第二代測(cè)序方法測(cè)定6種不同齡期烘干前后的蠶沙中菌群組成及豐富度，比較分析烘干加工前后、不同齡期蠶沙的菌群結(jié)構(gòu)的差異性，同時(shí)對(duì)其代謝功能進(jìn)行初步預(yù)測(cè)，以期為蠶沙藥材品質(zhì)及功能主治提供科學(xué)依據(jù)。

1材料與方法

11樣品采集

所有蠶沙樣品均采自鎮(zhèn)江市江蘇科技大學(xué)中國農(nóng)業(yè)科學(xué)院蠶業(yè)研究所。采樣時(shí)，在家蠶排便后，用無菌鑷子收集糞便放入無菌袋中，取部分樣品置于烘箱中50 ℃烘干備用，其余樣品取樣后立即用準(zhǔn)備好的冰袋冷藏，并迅速轉(zhuǎn)移到-80 ℃冰箱保存?zhèn)溆谩?/p>

12試劑與儀器

RS232G紫外分光光度計(jì)；Eppendorf；DYY6C電泳儀；PCR儀2720；凝膠成像系統(tǒng)（BG）；臺(tái)式紫外分析儀（BG）；Pico17離心機(jī)（Thermo）；Agilent，2100；BioTek，F(xiàn)Lx800；Promega，QuantiFluorTBS380；Illumina Miseq高通量測(cè)序儀（Illumina）。

Millipore超純水；瓊脂糖（英杰公司，75510019）；EB溴化乙錠清除液（生工，EX328）；PowerSoil DNA提取試劑盒（美國MoBio，12888）；5×反應(yīng)緩沖液；5×High GC Buffer；脫氧核糖核苷三磷酸（dNTP，10 mmol·L-1）；Axy Prep DNA Gel Extration kit試劑盒（Axygen，APGX500）；TruSeq Nano DNA LT Sample Prep kit試劑盒（FC1214001或FC1214002）；Agencourt AMPure XP Beads kit試劑盒（A63881）；QuantiT PicoGreen dsDNA Assay Kit試劑盒（Invitrogen，P7589）；Agilent High Sensitivity DNA Kit試劑盒（Agilent，50674626）。

13蠶沙樣品DNA的提取

將每種原始糞便樣品解凍、混勻后，各取025 g樣品進(jìn)行DNA的提取。樣品使用MoBio PowerSoil DNA Isolation Kit（12888）試劑盒提取DNA[11]，用08%的瓊脂糖凝膠電泳進(jìn)行分子大小判斷，利用紫外分光光度計(jì)對(duì)DNA進(jìn)行定量，取適量的樣品于離心管中，使用無菌水稀釋樣品至1 mg·L-1 [12]。

1416S rDNAV4區(qū)的PCR擴(kuò)增

以稀釋后的基因組DNA為模板，使用帶標(biāo)簽的16S rDNAV4 區(qū)特異引物520F （5′標(biāo)簽+GCACCTAAYTGGGYDTAAAGNG3′）和802R（5′TACNVGGGTATCTAATCC3′），使用高效和高保真酶（Q5高保真DNA聚合酶；5倍反應(yīng)緩沖液；5倍高GC緩沖液）進(jìn)行PCR擴(kuò)增，確保擴(kuò)增的效率和準(zhǔn)確性。

15PCR 產(chǎn)物的混樣和純化

PCR產(chǎn)物使用2%濃度的瓊脂糖凝膠進(jìn)行電泳檢測(cè)；根據(jù)PCR產(chǎn)物濃度進(jìn)行等濃度混樣，充分混勻后使用2%的瓊脂糖凝膠電泳檢測(cè)PCR產(chǎn)物，切取目的片段并用Axygen凝膠回收試劑盒回收目的片段。

1616S rDNA 文庫構(gòu)建與測(cè)序

161構(gòu)建文庫利用Illumina公司的TruSeq Nano DNA LT Library Prep Kit試劑盒進(jìn)行建庫。首先對(duì)基因片段進(jìn)行末端修復(fù)，將DNA 5′端突出的堿基切除，3′端缺失的堿基補(bǔ)齊，同時(shí)在5′端加上1個(gè)磷酸基團(tuán)。在修復(fù)平整的DNA片段3′端引入單堿基“A”，從而防止DNA片段的自連，同時(shí)保證DNA與3′端有1個(gè)突出T堿基的測(cè)序接頭相連。在連接酶的作用下，孵育含有標(biāo)簽的接頭與DNA片段，使其相連。通過PCR擴(kuò)增已經(jīng)加上接頭的DNA片段，利用BECKMAN AMPure XP beads試劑盒純化PCR體系。利用2%瓊脂糖凝膠電泳來對(duì)文庫做最終的片段選擇與純化。構(gòu)建好的文庫取1 μL，在Agilent Bioanalyzer機(jī)器上用Agilent High Sensitivity DNA Kit對(duì)文庫做2100質(zhì)檢，利用QuantiT PicoGreen dsDNA Assay Kit在Promega QuantiFluor上對(duì)文庫進(jìn)行定量，合格后，使用MiSeq進(jìn)行上機(jī)測(cè)序[1316]。endprint

162上機(jī)測(cè)序?qū)悠稤NA文庫均一化至10 nmol·L-1后等體積混合?；旌玫奈膸熘鸩较♂尪恐?～5 pmol·L-1后用Illumina MiSeq 測(cè)序儀測(cè)序。

17生物信息學(xué)分析

測(cè)序得到的原始數(shù)據(jù)（raw data），存在一定比例的干擾數(shù)據(jù)（dirty data），為了使信息分析的結(jié)果更加準(zhǔn)確、可靠，首先對(duì)原始數(shù)據(jù)進(jìn)行篩查、處理和質(zhì)量控制，得到高質(zhì)量的有效數(shù)據(jù)（effective tags）[17]運(yùn)用QIIME軟件識(shí)別疑問序列。要求序列長度≥150 bp，且不允許存在模糊堿基N，剔除5′端引物錯(cuò)配堿基數(shù)>1的序列和含有連續(xù)相同堿基數(shù)>8的序列。利用QIIME軟件調(diào)用USEARCH檢查并剔除嵌合體序列。最后去除barcode及引物序列，使用UCLUST軟件，按照97%相似性進(jìn)行操作分類單位（operational taxonomic units，OTU）聚類，挑選出每個(gè)OTU的代表序列[18]。利用Greengene 數(shù)據(jù)庫對(duì)代表序列進(jìn)行物種注釋。通過對(duì)OTUs進(jìn)行豐度，Alphadeversity，Betadeversity，LefSe分析以及物種在各個(gè)分類水平上的群落結(jié)果統(tǒng)計(jì)分析，得到微生物群落結(jié)構(gòu)組成[1920]。利用KEGG數(shù)據(jù)庫通過PICRUSt軟件對(duì)樣品的群落代謝功能進(jìn)行預(yù)測(cè)。

2結(jié)果

21測(cè)序結(jié)果的質(zhì)量分析

蠶沙樣品所測(cè)得的有效序列共274 750條，其中高質(zhì)量序列259 250 條，占序列總數(shù)的9436%。長度分布在300～500 bp，長度為450 bp 的序列最多，有20多萬條。從序列長度的分布來看，與16SrDNAV4 區(qū)序列長度吻合（表1）。

22OTU劃分和分類地位鑒定

利用Mothur軟件計(jì)算在97%的相似水平上每個(gè)樣品的OTU數(shù)量，OTU的數(shù)量可以代表樣品物種的豐度[21]。結(jié)果表明，在所有樣品中，OTU數(shù)量最多可達(dá)到369，最少為243。這表明蠶沙樣品的菌群豐富度很高，干品與鮮品之間菌群豐富度差異較大，不同齡期鮮品菌群豐富度明顯高于干品，且樣品均不存在未分類OTU（表2）。在門、綱、目水平上OTU幾乎可以全部被檢測(cè)，在屬、種水平上被檢測(cè)出的OTU的數(shù)量明顯減少（圖1）。

23Alpha多樣性分析

231稀疏曲線利用已測(cè)得16S rDNA序列中已知的各種OTU的相對(duì)比例，模擬輸入序列數(shù)目（小于總的樣品序列條數(shù)）與OTU個(gè)數(shù)產(chǎn)出間的相互關(guān)系。對(duì)基于97%相似度歸類的OTU用 Observed species，Chao1，Shannon 3種方法進(jìn)行稀疏曲線（rarefaction curve）繪制（圖2）。

從稀疏曲線來看，隨著測(cè)序數(shù)量的增加，稀釋曲線斜率逐漸降低，趨向平坦但未進(jìn)入平臺(tái)期，說明再增加測(cè)序數(shù)量也只會(huì)產(chǎn)生少量新的OTUs；從Chao1指數(shù)稀疏曲線來看（圖2C），Chao1指數(shù)在0～2×104的測(cè)序數(shù)量范圍內(nèi)增幅最大，在2×104以后曲線斜率下降，且4LX，5LX 的chao1曲線的增幅十分相似，表明2種蠶沙樣品菌群豐富度相似；從Shannon指數(shù)稀疏曲線來看（圖2B），6種蠶沙樣品的曲線數(shù)值先是直線上升直至趨向平坦，說明測(cè)序量足夠大，能覆蓋樣品中絕大多數(shù)微生物信息，再增加測(cè)序深度蠶沙樣品的菌群多樣性也不會(huì)發(fā)生改變。4LG蠶沙的3種曲線均明顯短于其他樣品組，說明測(cè)序數(shù)量足夠，已達(dá)到全部覆蓋。烘干后蠶沙樣品的菌群的豐富度及多樣性均有明顯下降，五齡蠶沙受烘干的影響較大，烘干處理可明顯降低五齡蠶沙的菌群多樣性與豐富度。

232Alpha多樣性指數(shù)對(duì)OTU豐度矩陣中的全體樣品根據(jù)最低測(cè)序深度統(tǒng)一進(jìn)行隨機(jī)重抽樣（即“序列量拉平處理”），使用QIIME軟件分別對(duì)每個(gè)樣品計(jì)算chao1，ACE，Shannon指數(shù)（表3）。

由表3可知，烘干后蠶沙樣品中chao1，ACE指數(shù)隨著齡期的增加而逐漸降低，說明干蠶沙隨著齡期的增加菌群的豐富度逐漸降低；新鮮蠶沙chao1，

ACE指數(shù)增長趨勢(shì)與其相反，說明隨著齡期的增加新鮮蠶沙的菌群豐富度逐漸升高。Shannon指數(shù)在烘干前、后的蠶沙樣品中均隨著齡期的增大而呈下降趨勢(shì)，說明齡期增大，蠶沙菌群的多樣性有所下降。

24各分類水平的分類學(xué)組成分析

241門水平的群落分類學(xué)組成分析蠶沙主要由變形菌門Proteobacteria、厚壁菌門Firmicutes、放線菌門Actinobacteria、擬桿菌門Bacteroidetes、藍(lán)細(xì)菌Cyanobacteria組成。其中相對(duì)豐度最大的優(yōu)勢(shì)菌門為變形菌門Proteobacteria，高達(dá)90%。烘干處理前后的放線菌門Actinobacteria、厚壁菌門Firmicutes豐度差異較大（圖3）。

242屬水平上的各樣品群落構(gòu)成分析除變形菌門腸桿菌科未分類菌屬（Enterobacteriaceae，759%），蠶沙細(xì)菌主要分布于以下10個(gè)屬，蒼白桿菌屬（Ochrobactrum，41%）、節(jié)細(xì)菌屬（Arthrobacter，37%）、不動(dòng)桿菌屬（Acinetobacter，35%）、腸球菌（Enterococcus，20%）、葡萄球菌屬（Staphylococcus，16%）、農(nóng)桿菌屬（Agrobacterium，07%）、甲基桿菌屬（Methylobacterium，06%）、克雷白氏桿菌屬（Klebsiella，05%）、腸桿菌屬（Enterobacter，04%）、沙雷氏菌屬（Serratia，03%）。3LG蠶沙的節(jié)細(xì)菌屬Enterococcus含量明顯高于其他；4LG，4LX，5LX的蒼白桿菌屬Ochrobactrum含量較高；5LG的腸球菌屬Enterococcus含量較高；4LG，3LG的不動(dòng)桿菌屬Acinetobacter含量較高（圖4）。使用R軟件，對(duì)豐度前50位的屬進(jìn)行聚類分析并繪制熱圖，熱圖中每小格代表所在樣品中某種菌的相對(duì)豐度，顏色越紅代表相對(duì)豐度越高。在屬水平上，4LX與5LX、3LG與3LX、4LG與5LG的菌群內(nèi)部結(jié)構(gòu)較相似，可分別聚為一類（圖5）。endprint

25Beta多樣性分析

基于UniFrac[22]距離上對(duì)屬水平各樣品之間進(jìn)行加權(quán)（weighted）和非加權(quán)（unweighted）的主坐標(biāo)分析（principal coordinate analysis，PCoA），同時(shí)采用基于UniFrac的非加權(quán)組平均法（unweighted pairgroup method with arithmetic means，UPGMA）進(jìn)行樣品聚類。6種樣品微生物群落構(gòu)成的PcoA圖中，坐標(biāo)軸括號(hào)中的百分比代表了對(duì)應(yīng)的主成分所能解釋的原始數(shù)據(jù)中差異的比例（圖6），基于UniFrac的加權(quán)和非加權(quán)的主坐標(biāo)分析（weighted and unweighted UniFrac principal coordinate analysis）其第一主成分和第二主成分的貢獻(xiàn)率分別為4423%，1811%和5535%，3204%。說明在加權(quán)和非加權(quán)主坐標(biāo)分析中所選主成分均可以充分解釋原始數(shù)據(jù)中的差異。通過圖中樣品點(diǎn)距離的遠(yuǎn)近，可觀察個(gè)體或群體間的差異，樣品點(diǎn)越相近說明樣品間微生物群落構(gòu)成越相似。由基于UniFrac的非加權(quán)主坐標(biāo)分析圖（圖6A）和基于UniFrac的加權(quán)主坐標(biāo)分析圖（圖6B）可知，3LG，3LX始終處于同一象限，說明三齡蠶沙樣品在烘干前后菌群多樣性無明顯改變，烘干處理對(duì)其影響較小。這點(diǎn)在UPGMA聚類圖中（圖7）也得到印證。4LX與5LX無法同其他樣品聚為一簇（圖7），這表明二者微生物群落結(jié)構(gòu)相似而與其他樣品相差較大，烘干對(duì)四齡、五齡蠶沙樣品有顯著影響。

26LefSe分析

LefSe分析，根據(jù)分類學(xué)組成對(duì)樣品按照不同的分組條件進(jìn)行線性判別分析，找出對(duì)樣品劃分產(chǎn)生顯著性差異影響的群落[2325]。將蠶沙樣品分為2組，烘干處理前與烘干處理后，以P<05，LOD>2為標(biāo)準(zhǔn)，在各分類水平找出具有顯著性差異的群落。柱形圖表示在2組中的差異菌屬信息，紅色表示在烘干處理蠶沙組中高出，綠色表示在新鮮蠶沙組中高出；所有菌在門、綱、目、科、屬，水平的差異信息用餅形圖表示。經(jīng)LEfSe分析（圖8），干蠶沙中的檸檬酸桿菌屬Citrobacter、愛文氏菌屬Ewingella顯著高于新鮮蠶沙，新鮮蠶沙中的甲基桿菌屬M(fèi)ethylobacterium、甲基桿菌科Methylobacteriaceae、鏈球菌科Streptococcaceae、藍(lán)細(xì)菌（cyanobacteria）、芽孢桿菌屬Bacillus和0319_6G20科菌（0319_6G20）顯著高于干蠶沙。由分析結(jié)果可知，干蠶沙中的條件致病菌要少于新鮮蠶沙，相比于新鮮蠶沙，干蠶沙更適于入藥。另干蠶沙中的檸檬酸桿菌屬可發(fā)酵葡萄糖產(chǎn)酸產(chǎn)氣，艾文氏菌屬可發(fā)酵D阿拉伯糖、D纖維二糖產(chǎn)酸，提示蠶沙的降糖作用可能與蠶沙中的相關(guān)菌群的作用有關(guān)。

27群落代謝功能分析

通過PICRUSt[26]軟件，根據(jù)KEGG數(shù)據(jù)庫中微生物代謝功能的類別對(duì)群落樣本的代謝功能進(jìn)行預(yù)測(cè)。統(tǒng)計(jì)蠶沙菌群相對(duì)豐度大于005的相關(guān)代謝功能（表4），蠶沙菌群涉及的代謝功能主要有膜轉(zhuǎn)運(yùn)、碳水化合物代謝、氨基酸代謝、脂質(zhì)代謝、核苷酸

3討論

本文對(duì)6種不同齡期、烘干處理前后的蠶沙樣品中菌群多樣性進(jìn)行分析與評(píng)價(jià)，對(duì)其16S rDNAV4區(qū)進(jìn)行Illumina MiSeq平臺(tái)測(cè)序，共得到274 750個(gè)有效序列。6種蠶沙樣品的3種稀疏曲線數(shù)值先是上升直至趨向平坦，說明測(cè)序量足夠大，能覆蓋樣品中絕大多數(shù)微生物信息。

蠶沙作為一種動(dòng)物藥，其本身為家蠶的排泄物，相對(duì)與傳統(tǒng)的中藥材，蠶沙中所含有的微生物較多。藥典中關(guān)于中藥制劑、生物制品等有相應(yīng)的微生物限度要求[27]，這也提示對(duì)于蠶沙這類微生物含量較多的藥材其菌群鑒定的必要性。本研究以蠶沙中的菌群結(jié)構(gòu)為切入點(diǎn)，首先明確了蠶沙中菌群的組成和相對(duì)豐度，并從菌群結(jié)構(gòu)的角度對(duì)中藥材蠶沙的烘干處理進(jìn)行了評(píng)價(jià)。從蠶沙樣品在不同分類水平上的OTU數(shù)量，樣品的chao1指數(shù)、ACE指數(shù)、shannon指數(shù)可明顯的看出烘干處理后的蠶沙樣品中菌群的豐度和多樣性都低于新鮮蠶沙，芽孢桿菌、鏈球菌等一些條件致病菌在干蠶沙中也明顯減少，充分說明烘干處理可有利于蠶沙藥用。

在對(duì)蠶沙樣品進(jìn)行分類學(xué)鑒定時(shí)發(fā)現(xiàn)蠶沙主要由變形菌門Proteobacteria、厚壁菌門Firmicutes、放線菌門Actinobacteria、擬桿菌門Bacteroidetes、藍(lán)細(xì)菌Cyanobacteria組成。其優(yōu)勢(shì)菌門為變形菌門。各分類水平的菌群鑒定均顯示3齡蠶沙的菌群多樣性與結(jié)構(gòu)同4齡、5齡蠶沙相差較大，原因是家蠶在1～3齡期為稚蠶，4～5齡期為壯蠶，分屬家蠶生長發(fā)育的不同階段，其腸道菌群本身的差別較大。此外，有報(bào)道[28]通過剖取家蠶腸道內(nèi)容物進(jìn)行16SrDNA測(cè)序，確定了家蠶腸道細(xì)菌的優(yōu)勢(shì)類群為葡萄球菌屬、陰溝腸桿菌（腸桿菌屬）、蠟樣芽胞桿菌（芽孢桿菌屬）。這與本研究所得的蠶沙菌群的優(yōu)勢(shì)類群有所不同，推測(cè)家蠶在食用桑葉的過程中體內(nèi)的微生物系統(tǒng)發(fā)生了改變，可能有新的菌種產(chǎn)生。加之蠶沙中存在消化不完全的桑葉，導(dǎo)致蠶沙的菌群一部分來源于桑葉的內(nèi)生菌，值得進(jìn)一步深入研究。

結(jié)合對(duì)人腸道菌群的相關(guān)研究[29]可發(fā)現(xiàn)蠶沙的菌群結(jié)構(gòu)與人體的腸道菌群結(jié)構(gòu)差異較大，若服用不當(dāng)其含有的條件致病菌極有可能破壞人體的腸道菌群平衡，從而不但達(dá)不到治病的效果反而會(huì)致病。本研究結(jié)果也提示，在蠶沙入藥時(shí)要經(jīng)過適當(dāng)加工處理，鮮品不適于入藥，三齡蠶沙不適于入藥；同時(shí)充分考慮其可能對(duì)機(jī)體微生物生態(tài)環(huán)境產(chǎn)生改變，從而對(duì)機(jī)體產(chǎn)生影響。蠶沙菌群優(yōu)勢(shì)種屬代謝功能的預(yù)測(cè)說明蠶沙對(duì)機(jī)體糖、脂肪、蛋白代謝有一定的影響，表明蠶沙的藥效可能與其菌群的作用有關(guān)，尚需進(jìn)一步深入研究。

本實(shí)驗(yàn)利用Illumina MiSeq 高通量測(cè)序技術(shù)分析了蠶沙菌群多樣性，全面而準(zhǔn)確地分析了蠶沙菌群的種類組成，發(fā)現(xiàn)蠶沙菌群分布相對(duì)集中，主要分布于腸桿菌科（Enterobacteriaceae，80%），以下7個(gè)種屬為蠶沙菌群的優(yōu)勢(shì)種群：變形菌門下的蒼白桿菌屬Ochrobactrum、不動(dòng)桿菌屬Acinetobacter、節(jié)細(xì)菌屬Arthrobacter、甲基桿菌Methylobacterium、腸桿菌屬Enterobacter、克雷白氏桿菌屬Klebsiella及厚壁菌門下的葡萄球菌Staphylococcus，為蠶沙的后續(xù)研究及開發(fā)利用提供了科學(xué)依據(jù)。endprint

[參考文獻(xiàn)]

[1]林庚庭， 魏克民， 梁衛(wèi)青， 等 蠶沙提取物對(duì)再生障礙性貧血小鼠細(xì)胞因子影響的實(shí)驗(yàn)研究[J]. 中國中醫(yī)藥科技， 2008， 15（2）：117.

[2]施文君， 楊云帆， 朱思然， 等 蠶沙抗炎鎮(zhèn)痛作用實(shí)驗(yàn)研究[J]. 亞太傳統(tǒng)醫(yī)藥， 2013， 9（9）：44.

[3]Raghavendra R， Neelagund S， Kuluvar G， et al. Protective effect of partially purified 35 kDa protein from silk worm （Bombyx mori） fecal matter against carbon tetrachloride induced hepatotoxicity and in vitro antiviral properties[J]. Pharmaceutical Biol， 2010， 48（12）：1426.

[4]楊倉良， 齊英杰 動(dòng)物本草[M]北京：中醫(yī)古籍出版社，2001

[5]楊海霞， 朱祥瑞， 房澤民 蠶沙的開發(fā)利用研究進(jìn)展[J]. 蠶桑通報(bào)， 2002， 33（3）： 9.

[6]鐘楊生， 陳芳艷， 劉佳， 等 蠶沙堆肥過程中微生物群落的動(dòng)態(tài)變化[J]. 廣東蠶業(yè)， 2012（4）：21.

[7]史才娟， 陳芳艷， 鐘楊生， 等 不同齡期家蠶蠶沙部分理化和微生物學(xué)性質(zhì)研究[J]. 廣東農(nóng)業(yè)科學(xué)， 2013， 40（3）：109

[8]Reetta M， Satokari1， Elaine E Diversity of Bidobacterium and Lactobacillus spp in breastfed and formulafed infants as assessed by 16S rDNA sequence differences[J]. Microb Ecol Health Dis， 2002， 14： 97.

[9]Penders J， Vink C， Driessen C， et al. Quantification of Bifidobacterium spp Escherichia coli and Clostridium difficile in faecal samples of breastfed and formulafed infants by realtime PCR[J]. FEMS Microbiol Lett， 2005， 243（1）： 141.

[10]Klaassens E S， Vos W M D， Vaughan E E Metaproteomics approach to study the functionality of the microbiota in the human infant gastrointestinal tract[J]. Appl Environ Microbiol， 2007， 73（4）： 1388.

[11]倪桂萍， 王延平， 王華田， 等 楊樹人工林土壤細(xì)菌DNA 的提取與擴(kuò)增[J]. 山東大學(xué)學(xué)報(bào)：理學(xué)版，2013， 48（5）： 23.

[12]陳澤斌， 李冰， 王定康， 等 應(yīng)用Illumina MiSeq高通量測(cè)序技術(shù)分析玉米內(nèi)生細(xì)菌多樣性[J]. 現(xiàn)代食品科技， 2016（2）： 113.

[13]Caporaso J G， Lauber C L， Walters W A， et al. Global patterns of 16s rRNA diversity at a depth of millions of sequences per sample[J]. Proc Natl Acad Sci USA， 2011， 108 （25）： 4516.

[14]Youssef N， Sheik C S， Krumholz L R， et al. Comparison of species richness estimates obtained using nearly complete fragments and simulated pyrosequencinggenerated fragments in 16S rRNA genebased environmental surveys[J]. Appl Environ Microbiol， 2009， 75（16）： 5227.

[15]Hess M， Sczyrba A， Egan R， et al. Metagenomic discovery of biomassdegrading genes and genomes from cow rumen[J]. Science， 2011， 331（6016）： 463.

[16]Luo C， Tsementzi D， Kyrpides N， et al. Direct comparisons of Illumina vs Roche 454 sequencing technologies on the same microbial community DNA sample[J]. PLoS ONE， 2012， 7（2）： e30087

[17]Edgar R C， Haas B J， Clemente J C， et al. UCHIIME improves sensitivity and speed of chimera detection[J]. Bioinformatics， 2011， 27（16）： 2194.endprint

[18]Edgar R C Search and clustering orders of magnitude faster than BLAST[J]. Bioinformatics， 2010， 26（19）： 2460.

[19]Wang Q， Garrity G M， Tiedje J M， et al. Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy[J]. Appl Environ Microbiol， 2007， 73（16）： 5261.

[20]Desantis T Z， Hugenholtz P， Larsen N， et al. Greengenes， a chimerachecked 16S rRNA gene database and workbench compatible with ARB[J]. Appl Environ Microbiol， 2006， 72（7）： 5069

[21]Zhang M L， Zhang M H， Zhang C H， et al. Pattern extraction of structural responses of gut microhiota to rotavirus infection via multivariate statistical analysis of clone library data[J]. FEMS Microhiol Eco1， 2009， 70： 177.

[22]Legiewicz M， Lozupone C， Knight R， et al. Size， constant sequences， and optimal selection[J]. RNA， 2005， 11（11）： 1701

[23]Segata N Metagenomic biomarker discovery and explanation[J]. Genome Biol， 2011， 12（6）：1.

[24]Scher J U， Sczesnak A， Longman R S， et al. Expansion of intestinal Prevotella copri correlates with enhanced susceptibility to arthritis[J]. Elife Sci， 2013， 2（1629）： e01202.

[25]Wang Y， Sheng H F， He Y， et al. Comparison of the levels of bacterial diversity in freshwater， intertidal wetland， and marine sediments by using millions of Illumina tags[J]. Appl Environ Microbiol， 2012， 78（23）： 8264.

[26]Langille M G， Zaneveld J， Caporaso J G， et al. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences[J]. Nat Biotechnol， 2013， 31（9）： 814.

[27]中國藥典一部[S] 2015：6

[28]袁志輝， 藍(lán)希鉗， 楊廷， 等 家蠶腸道細(xì)菌群體調(diào)查與分析[J]. 微生物學(xué)報(bào)， 2006， 46（2）： 285.

[29]劉瑞雪， 李勇超， 張波 腸道菌群微生態(tài)平衡與人體健康的研究進(jìn)展[J]. 食品工業(yè)科技， 2016， 37（6）：383.

[責(zé)任編輯呂冬梅]endprint