鄒邦興,張宇靖,卿圣環(huán),蔣亞明,張友洪,肖文福,周安蓮

(四川省農(nóng)業(yè)科學(xué)院蠶業(yè)研究所,四川 南充 637000)

【研究意義】人工飼料養(yǎng)蠶始于20世紀(jì)30年代日本[1],而我國則從20世紀(jì)60年代開始研究[2],截止目前,人工飼料養(yǎng)蠶技術(shù)已在浙江、江蘇、廣西、四川、山東等18個(gè)蠶桑省區(qū)示范推廣。和傳統(tǒng)桑葉養(yǎng)蠶相比,人工飼料養(yǎng)蠶具有省力高效、安全健康和工廠化操作等優(yōu)勢,是蠶桑生產(chǎn)上的一次技術(shù)革新。隨著人工飼料養(yǎng)殖技術(shù)的推廣,關(guān)于飼料養(yǎng)蠶蠶沙的處理和利用成為行業(yè)面臨的難題,區(qū)別于傳統(tǒng)桑葉育,人工飼料養(yǎng)蠶需要從各方面進(jìn)行成本控制以達(dá)到節(jié)本增效的目的,而當(dāng)前人工飼料養(yǎng)蠶過程中產(chǎn)生的蠶沙主要用于水產(chǎn)品養(yǎng)殖或者直接丟棄。蠶沙作為蠶桑生產(chǎn)過程中的大宗副產(chǎn)物,攜帶大量病原菌,易導(dǎo)致蠶病爆發(fā)和造成環(huán)境污染[3]。陳春等[4]研究發(fā)現(xiàn),每生產(chǎn)100 kg蠶繭可得蠶沙約600 kg,根據(jù)2021年全國蠶繭產(chǎn)量7.172×105t計(jì)算[5],年產(chǎn)蠶沙約4.30×106t,約占2022年全國化肥總產(chǎn)量的7.9%。蠶沙中富含葉綠素[6]、β-胡蘿卜素[6-8]和氮磷鉀[3-4]等元素,因此極具醫(yī)藥和肥料開發(fā)價(jià)值,如何將其合理利用,已成為行業(yè)人員共同面臨的問題。 【前人研究進(jìn)展】作為生物界最龐大、最復(fù)雜的微生態(tài)系統(tǒng),腸道微生物本身及其代謝產(chǎn)物與人類健康和環(huán)境保護(hù)息息相關(guān),因此,研究微生物群落結(jié)構(gòu)與宿主和環(huán)境之間的關(guān)系尤為重要。隨著家蠶人工飼料養(yǎng)蠶技術(shù)的推廣,對飼料和桑葉不同模式飼養(yǎng)下家蠶腸道微生物群落結(jié)構(gòu)的研究已成為熱點(diǎn)。桑基魚塘是我國最早的農(nóng)業(yè)循環(huán)利用模式,也是蠶沙利用最成功和最多的典型實(shí)例[9],但蠶沙中攜帶大量病原微生物,如直接利用可能導(dǎo)致病害發(fā)生和環(huán)境污染,因此,針對蠶沙中大量病原微生物存在,通過對其腸道和蠶沙微生物的研究有助于蠶沙的安全高效利用。荊文輝等[10]研究發(fā)現(xiàn),以日本配方飼料(NOSAN)和自制配方飼料飼養(yǎng)家蠶,家蠶中腸的細(xì)菌菌群存在顯著差異,NOSAN飼料育家蠶中腸的菌群豐度較高,而桑葉飼育組與自制飼料飼育組的菌群豐度相當(dāng)。李剛等[11]以蓖麻蠶品種B7為對象,分別喂食蓖麻葉和臭椿葉,腸道微生物檢測結(jié)果顯示,蓖麻葉育和臭椿葉育的蓖麻蠶腸道優(yōu)勢菌表現(xiàn)相似,其菌群結(jié)構(gòu)和組成與其食性存在一定關(guān)系。史才娟等[12]通過對各齡期蠶沙中的微生物群落進(jìn)行檢測,發(fā)現(xiàn)各齡期蠶沙微生物群落在數(shù)量上有差異,但不顯著,細(xì)菌以3齡蠶的微生物數(shù)量最少,2齡與4齡接近,真菌幾個(gè)齡期間數(shù)量差異不大,放線菌以2齡最多,3齡最少?！颈狙芯壳腥朦c(diǎn)】前期研究主要集中在家蠶腸道微生物群落結(jié)構(gòu)組成,關(guān)于蠶沙微生物群落結(jié)構(gòu)的研究較少。蠶沙作為養(yǎng)蠶過程中的副產(chǎn)物,其產(chǎn)量高,對環(huán)境影響大,因此針對蠶沙的處理一直以來是行業(yè)研究熱點(diǎn)和難點(diǎn)?！緮M解決的關(guān)鍵問題】采用Illumina NovaSeq二代測序方法對全齡飼料育和全齡桑葉育不同模式下5齡蠶沙進(jìn)行細(xì)菌群落組成和豐度檢測,解析不同飼養(yǎng)模式下蠶沙微生物群落結(jié)構(gòu)組成,為后期蠶沙的無害化處理及綜合利用提供科學(xué)參考。

1 材料與方法

1.1 實(shí)驗(yàn)設(shè)計(jì)

共設(shè)計(jì)全齡飼料育和全齡桑葉育2組,所用蠶品種為“9211·9215×川58·川62”(又名“川蠶30號”),飼養(yǎng)階段在四川省農(nóng)業(yè)科學(xué)院蠶業(yè)研究所完成。其中全齡飼料育處理溫濕度為1～2齡溫度29～30 ℃,濕度90%;3齡溫度28～29 ℃,濕度85%,4齡溫度26～27 ℃,濕度80%,5齡溫度25～26 ℃,濕度80%。桑葉育處理采用自然溫濕度,其中1～2齡期溫度不低于28 ℃。

1.2 材料采集

2022年秋季收集2種飼養(yǎng)模式下5齡盛食期蠶沙,剔除雜質(zhì),然后每個(gè)處理5個(gè)重復(fù)送樣,飼料育組樣本編號為T1～T5(記作T,作為處理組),桑葉育組樣本編號為CK1～CK5(記作CK,作為對照組),然后以干冰保存送至上海派森諾生物科技股份有限公司。

1.3 蠶沙微生物檢測

DNA提取、16S rDNA區(qū)域擴(kuò)增、序列測定和文庫構(gòu)建均由上海派森諾生物科技股份有限公司完成。稱取0.2 g蠶沙,提取微生物DNA,利用瓊脂糖凝膠電泳進(jìn)行DNA樣品檢測。然后對其16S V3～V4區(qū)域進(jìn)行擴(kuò)增,擴(kuò)增引物為通用引物F:5′-ACTCCTACGGGAGGCAGCA-3′R:5′-GGACTACHVGGGTWTCTAAT-3′。用2%瓊脂糖凝膠電泳檢測PCR 產(chǎn)物,目的片段用 Axygen 凝膠回收試劑盒回收并純化。純化后產(chǎn)物采用Illumina NovaSeq平臺(tái)對群落DNA片段進(jìn)行雙端(Paired-end)測序,同時(shí)采用DADA2和Vsearch方法進(jìn)行去引物、質(zhì)量過濾、去噪、拼接和去嵌合體等步驟以獲得高質(zhì)量的ASV和OTU序列。然后采用QIIME2的classify-sklearn算法對每個(gè)ASVs的特征序列或每個(gè)OTU的代表序列在QIIME2軟件中使用默認(rèn)參數(shù)進(jìn)行物種注釋。

1.4 數(shù)據(jù)分析

采用WPS 2020進(jìn)行數(shù)據(jù)整理,GraphPad Prism 9進(jìn)行雙尾T檢驗(yàn),Origin 2022進(jìn)行作圖。

2 結(jié)果與分析

2.1 蠶沙微生物DNA測序結(jié)果

從表1可知,桑葉育對照組蠶沙細(xì)菌微生物測序在原始數(shù)據(jù)(Input)、有效序列(Denoised)、拼接后序列(Merged)、去除嵌合體后的高質(zhì)量序列(Non-chimeric)和去除singleton后的序列(Non-singleton)均比全齡飼料育處理組高,且存在顯著差異(P<0.05),但在去除低質(zhì)量序列(Filtered)上兩者差異不顯著(P>0.05)。

2.2 蠶沙微生物的OTU聚類分析

通過高通量測序得到桑葉育對照組和全齡飼料育處理組蠶沙微生物的有效序列分別為41 691條和38 604條,使用QIIME2軟件在97%相似水平上聚類,通過注釋獲得2組樣本的OTU數(shù)目,分別為全齡飼料育處理組1983個(gè)、桑葉育對照組529個(gè),2組共有的OTU為133個(gè)(圖1)。說明全齡飼料育較桑葉育可以增加蠶沙微生物豐度。

圖1 不同飼養(yǎng)模式蠶沙微生物OTU聚類Fig.1 OTU cluster of silkworm sand microorganisms in different feeding modes

2.3 蠶沙微生物Alpha多樣性分析

2.3.1 Alpha多樣性 從圖2可知,全齡飼料育組和桑葉育對照組文庫覆蓋率在99%以上,表明所建基因文庫能較完整、真實(shí)地反映樣本蠶沙細(xì)菌群落結(jié)構(gòu)。Chao1指數(shù)和Observed_species指數(shù)反映樣本中細(xì)菌群落的物種豐富度,其值越大,代表樣細(xì)菌群落中的物種總數(shù)越多;Shannon指數(shù)和Simpson指數(shù)代表樣本細(xì)菌群落的物種多樣性,其值越大,代表樣細(xì)菌群落的物種多樣性越大;Pielou_e用于反映物種個(gè)體數(shù)目在群落中分配的均勻程指數(shù)。蠶沙細(xì)菌群落包含的物種數(shù)目在Chao1指數(shù)、Observed_species指數(shù)和均勻度(Pielou_e)指數(shù)的變化趨勢一致,均是全齡飼料育處理組數(shù)據(jù)高于桑葉育對照組,且呈顯著差異(P<0.05);在細(xì)菌群落多樣性方面,Shannon指數(shù)和Simpson指數(shù)變化趨勢一致,均是全齡飼料育處理組高于桑葉育對照組,且呈極顯著差異(P<0.01)。說明,全齡飼料育組蠶沙細(xì)菌群落的物種數(shù)目和物種多樣性顯著高于桑葉育對照組。

2.3.2 蠶沙微生物Chao1和Shannon稀疏曲線 為驗(yàn)證上述所得Alpha多樣性測序數(shù)據(jù)量是否足以反應(yīng)蠶沙微生物多樣性,利用已測得樣本各種OTU數(shù)繪制稀疏曲線,通過該曲線可以在相同測序深度下,比較樣本中OTU數(shù)的多少,從而在一定程度上衡量每個(gè)樣本的物種數(shù)目和多樣性水平。從圖3可知,隨著抽樣水平增加,蠶沙細(xì)菌微生物的Chao1指數(shù)和Shannon指數(shù)的稀疏曲線呈快速增長到趨于平衡,說明測序量足夠大,已基本覆蓋樣品中所有物種,同時(shí)在同一抽樣水平上全齡飼料育處理組樣本中OTU均高于桑葉育對照組,說明全齡飼料育處理組蠶沙微生物在物種數(shù)目和多樣性方面均較桑葉育對照組高。

2.4 蠶沙微生物群落結(jié)構(gòu)

2.4.1 不同飼養(yǎng)模式蠶沙微生物物種組成分類單元數(shù) 通過對抽平后的OTU進(jìn)行統(tǒng)計(jì),獲得每個(gè)樣本中的微生物群落在域、門、綱、目、科、屬、種7個(gè)分類水平各自含有的分類單元的數(shù)量,統(tǒng)計(jì)的數(shù)目中已去除包含unclassified、uncultured、uncultivated、unknown、metagenome等字段的分類單元(表2)。在門、綱、目和科分類水平,桑葉育對照組的分類單元數(shù)量均較全齡飼料育處理組高;在屬(Genus)水平上,兩者一致;在種分類水平上,全齡飼料育處理組的分類單元數(shù)量較桑葉育對照組高。說明,隨著分類單元逐漸細(xì)化,桑葉育對照組和全齡飼料育處理組蠶沙在細(xì)菌群落各分類水平的單元數(shù)量呈先上升后下降趨勢,其中桑葉育對照組在科分類水平達(dá)到最高值,全齡飼料育處理組則在屬分類水平達(dá)到最高值。

表1 樣本測序量統(tǒng)計(jì)

圖2 不同飼養(yǎng)模式下蠶沙細(xì)菌Alpha多樣性Fig.2 Bacterial diversity of silkworm excrement under different feeding modes

圖3 蠶沙細(xì)菌微生物Chao1和Shannon稀疏曲線Fig.3 Chao1 and Shannon sparse curve of bacteria in silkworm excrement

2.4.2 不同飼養(yǎng)模式蠶沙微生物物種組成分析 通過對去除singleton后的特征序列進(jìn)行統(tǒng)計(jì),實(shí)現(xiàn)桑葉育對照組和全齡飼料育處理組蠶沙在門、綱、目、科、屬、種6個(gè)分類水平上的組成分布。在門水平(圖4-A)上,桑葉育對照組蠶沙細(xì)菌群落中相對豐度排前3的依次為藍(lán)細(xì)菌(Cyanobacteria,55.38%)、厚壁菌門(Firmicutes,32.50%)和變形菌門(Proteobacteria,11.55%);全齡飼料育處理組蠶沙細(xì)菌群落中相對豐度排前3依次為變形菌門(Proteobacteria,68.18%)、厚壁菌門(Firmicutes,24.10%)和放線菌門(Actinobacteria,7.28%)。其中厚壁菌門和變形菌門為桑葉育對照組和全齡飼料育處理組蠶沙細(xì)菌優(yōu)勢菌門。

在目水平(圖4-B)上,桑葉育對照組蠶沙細(xì)菌群落中相對豐度大于1%的有腸桿菌目(Enterobacteriales)、木霉菌目(streptophyta)、乳酸桿菌目(Lactobacillales)、芽孢桿菌目(Bacillales)、假單胞菌目(Pseudomonadales)、伯克氏菌目(Burkholderiales)和立克次氏體目(Rickettsiales),其中相對豐度排前3的依次為木霉菌目(55.24%)、乳酸桿菌目(28.22%)和芽孢桿菌目(4.15%)。全齡飼料育處理組蠶沙細(xì)菌群落中相對豐度大于1%的有腸桿菌目(Enterobacteriales)、乳酸桿菌目(Lactobacillales)、芽孢桿菌目(Bacillales)、假單胞菌目(Pseudomonadales)、放線菌目(Actinomycetales)、伯克氏菌目(Burkholderiales)和其他(Others),其中相對豐度排前3的依次為腸桿菌目(58.90%)、乳酸桿菌目(15.14%)和芽孢桿菌目(8.63%)。乳酸桿菌目和芽孢桿菌目均為桑葉育對照組和全齡飼料育處理組蠶沙細(xì)菌目水平優(yōu)勢菌目。

表2 不同飼養(yǎng)模式下蠶沙細(xì)菌群落在各分類水平組成

A、B、C、D分別代表在蠶沙細(xì)菌群落在門、綱、科、屬分類水平;橫坐標(biāo)表示樣品分類,縱坐標(biāo)表示在各分類水平下相對豐度(豐度值前10)。A, B, C and D represent the classification level of the bacterial community in the silkworm excrement at the level of phylum, class, family and genus. The horizontal axis represents sample classification, and the vertical axis represents relative abundance at each classification level (top 10 abundance values).

在科水平(圖4-C)上,桑葉育對照組蠶沙細(xì)菌群落中相對豐度大于1%的有腸桿菌科(Enterobacteriaceae)、腸球菌科(Enterococcaceae)、莫拉氏菌科(Moraxellaceae)、小桿菌科(Exiguobacteraceae)、Mitochondria、草酸桿菌科(Oxalobacteraceae)和其他(Others),其中相對豐度排前3的依次為其他(57.94%)、腸球菌科(27.58%)和小桿菌科(3.88%)。全齡飼料育處理組蠶沙細(xì)菌群落中相對豐度大于1%的有腸桿菌科(Enterobacteriaceae)、腸球菌科(Enterococcaceae)、莫拉氏菌科(Moraxellaceae)、細(xì)球菌科 (Micrococcaceae)、小桿菌科(Exiguobacteraceae)、叢毛單胞菌科(Comamonadaceae)其他(Others),其中相對豐度排前3的依次為腸桿菌科(58.90%)、腸球菌科(14.04%)和莫拉氏菌科(7.07%)。腸球菌科均為桑葉育對照組和全齡飼料育處理組蠶沙細(xì)菌科分類水平優(yōu)勢菌科。

在屬水平(圖4-D)上,桑葉育對照組蠶沙細(xì)菌群落中相對豐度大于1%的有腸球菌屬(Enterococcus)、不動(dòng)桿菌屬(Acinetobacter)、微小桿菌屬(Exiguobacterium)、蠅子草屬(Silene)、青枯菌屬(Ralstonia)和其他(Others),其中相對豐度排前3的依次為其他(60.21%)、腸球菌屬(27.47%)和微小桿菌屬(3.58%)。全齡飼料育處理組蠶沙細(xì)菌群落中相對豐度大于1%的有腸球菌屬(Enterococcus)、不動(dòng)桿菌屬(Acinetobacter)、節(jié)桿菌屬(Arthrobacter)微小桿菌屬(Exiguobacterium)、葡萄球菌屬(Staphylococcus)和其他,其中相對豐度排前3的依次為其他(67.74%)、腸球菌屬(13.28%)和不動(dòng)桿菌屬(6.79%)。腸球菌屬均為桑葉育對照組和全齡飼料育處理組蠶沙細(xì)菌屬分類水平優(yōu)勢菌屬。

2.4.3 蠶沙微生物物種Beta多樣性和物種組成 PCoA分析可以直觀觀察樣本內(nèi)或樣本間的差異,在坐標(biāo)軸上的投影距離越近,表明樣本在相應(yīng)維度中的群落組成越相似。如圖5-A所示,基于Bray-Curtis距離、置信度95%計(jì)算聚類,桑葉育對照組和全齡飼料育處理組之間投影距離較大,說明2種飼育模式下家蠶蠶沙細(xì)菌物種組成存在差異。為進(jìn)一步比較樣本間的物種組成差異,默認(rèn)使用平均豐度前20位的科的豐度數(shù)據(jù)繪制熱圖,實(shí)現(xiàn)2組蠶沙微生物細(xì)菌的物種豐度分布趨勢的展示,如圖5-B所示,在科分類水平上,桑葉育對照組和全齡飼料育處理組樣品間菌群結(jié)構(gòu)組成存在一定差異。其中乳桿菌科(Lactobacillaceae)、甲基桿菌科(Methylobacteriaceae)、小桿菌科(Exiguobacteraceae)、橙單胞菌科(Aurantimonadaceae)、Mitochondria、鏈球菌科(Streptococcaceae)、草酸桿菌科(Oxalobacteraceae)和鞘脂單胞菌科(Sphingomonadaceae)在桑葉育對照組中相對豐度較全齡飼料育處理組中相對豐度高;肉桿菌科(Carnobacteriaceae)、微桿菌科(Microbacteriaceae)、芽孢桿菌科(Bacillaceae)、腸桿菌科(Enterobacteriaceae)、細(xì)球菌科(Micrococcaceae)、葡萄球菌科(Staphylococcaceae)、慢生根瘤菌科(Bradyrhizobiaceae)、動(dòng)球菌科(Planococcaceae)、莫拉氏菌科(Moraxellaceae)、叢毛單胞菌科(Comamonadaceae)和假單胞菌科(Pseudomonadaceae)在全齡飼料育處理組中相對豐度較桑葉育對照組中相對豐度高;腸球菌科(Enterococcaceae)在2組中相對豐度相對一致。

圖5 不同飼育模式蠶沙細(xì)菌Beta多樣性(A)和在科水平豐度聚類(B)Fig.5 Beta diversity (A) and cluster of abundance at family level(B) of silkworm sand bacteria in different feeding modes

3 討 論

蠶沙的安全利用依然是產(chǎn)業(yè)發(fā)展過程中面臨的難題之一。蠶沙富含多種氨基酸和礦質(zhì)元素,能夠?yàn)橥寥捞峁┓柿?因此,蠶沙的肥料化利用成為當(dāng)前蠶沙綜合利用中前景較廣和用量最多的方向。但蠶沙中攜帶大量病原微生物,如直接還田利用會(huì)造成土壤水體污染和給蠶桑生產(chǎn)帶來安全隱患[13],因此在蠶沙肥料化利用過程中需要經(jīng)過發(fā)酵和熟化滅菌處理才能用于土壤肥培管理,增加土壤有機(jī)質(zhì)[13-16],除了為主要農(nóng)作物[17]和經(jīng)濟(jì)作物[18-22]提供肥力外,蠶沙還可用于連作障礙土壤修復(fù)[23]、土壤碳庫平衡[24]和沼氣循環(huán)利用[25-26]等。

近年來,隨著工廠化全齡人工飼料養(yǎng)蠶技術(shù)的推廣,關(guān)于飼料養(yǎng)蠶蠶沙的綜合利用也逐漸成為熱點(diǎn)問題,而當(dāng)前還未見關(guān)于全齡飼料育蠶沙微生物群落結(jié)構(gòu)的報(bào)道。本研究通過對同一家蠶品種采用不同模式飼養(yǎng)(即全齡飼料育和全齡桑葉育),對其五齡蠶沙進(jìn)行細(xì)菌多樣性分析,并比較2種飼養(yǎng)模式下家蠶蠶沙細(xì)菌群落結(jié)構(gòu)組成,飼料育較桑葉育能夠顯著提高蠶沙細(xì)菌群落結(jié)構(gòu)的OTU個(gè)數(shù)、物種數(shù)目和物種多樣性,其原因可能與飼料組成成分有關(guān),對照組采用桑葉喂食家蠶,其蠶沙細(xì)菌群落結(jié)構(gòu)主要受桑葉細(xì)菌群落結(jié)構(gòu)影響,全齡飼料育處理組采用人工配制飼料,該飼料主要由桑葉、玉米、大豆、礦物質(zhì)和少量微量元素組成,從食物來源看,全齡飼料育處理組家蠶取食的食物中細(xì)菌菌落結(jié)構(gòu)較對照豐富,這可能是導(dǎo)致全齡飼料育處理組蠶沙細(xì)菌群落在物種豐富度和物種多樣性方面較桑葉育對照組高的原因之一,具體情況有待進(jìn)一步驗(yàn)證。此外還可能與養(yǎng)蠶環(huán)境、桑葉品種、葉片老熟和桑葉儲(chǔ)存環(huán)境和儲(chǔ)存時(shí)長存在一定關(guān)系。

全齡飼料育處理組和桑葉育對照組之間的細(xì)菌物種組成在相對豐度前20的科分類水平存在一定差異,其中乳桿菌科(Lactobacillaceae)等8種科水平細(xì)菌在桑葉育對照組中相對豐度較全齡飼料育處理組高,肉桿菌科(Carnobacteriaceae)等11科水平細(xì)菌在全齡飼料育處理組中相對豐度較桑葉育對照高,而腸球菌科(Enterococcaceae)在兩組中相對豐度相對一致。腸球菌為需氧或兼性厭氧型,廣泛分布于自然環(huán)境中,革蘭氏陽性,適宜在高鹽(6.5%NaCl)、高堿(pH 9.6)和10～45 ℃環(huán)境下生長,極易形成耐藥性,是造成新生兒敗血癥的第三病因,因此,在蠶沙的無害化處理過程中應(yīng)重視腸球菌的滅活處理,避免造成環(huán)境污染和影響人類健康。同時(shí)樣品內(nèi)部間也存在一定差異,其原因有待進(jìn)一步研究。在細(xì)菌物種分類方面,在門水平上,厚壁菌門和變形菌門均為桑葉育對照組和全齡飼料育處理組蠶沙細(xì)菌門分類水平優(yōu)勢菌,其相對豐度分別為32.50%、11.55%和68.18%、24.10%。在目水平上,乳酸桿菌目均為桑葉育對照組和全齡飼料育處理組蠶沙細(xì)菌目分類水平優(yōu)勢菌,其相對豐度分別為22.82%和15.14%。在科水平上,腸球菌科均為桑葉育對照組和全齡飼料育處理組蠶沙細(xì)菌科分類水平優(yōu)勢菌,其相對豐度分別27.58%和14.04%。 在屬水平上,腸球菌屬均為桑葉育對照組和全齡飼料育處理組蠶沙細(xì)菌屬分類水平優(yōu)勢菌,其相對豐度分別為27.47%和13.28%。物種分類學(xué)結(jié)果顯示,不同飼養(yǎng)模式下家蠶蠶沙細(xì)菌群落結(jié)構(gòu)在門、目、科、屬水平上均存在差異(P<0.05),說明不同飼料來源對蠶沙微生物的物種組成和相對豐度有一定影響。

4 結(jié) 論

全齡桑葉育和全齡飼料育不同飼育模式下蠶沙微生物群落物種組成及相對豐度存在差異,飼料育蠶沙細(xì)菌群落結(jié)構(gòu)的物種豐度和物種多樣性顯著高于桑葉育;在物種組成結(jié)構(gòu)單元上,桑葉育在門、綱、目和科分類單元數(shù)量均較飼料育高,在屬水平上,兩者一致,在種分類單元數(shù)量上飼料育較桑葉育高。