程曉亮，陳劍波*，王樹(shù)華，韓文儒，穆秀梅，淡江華

(1.山西農(nóng)業(yè)大學(xué) 動(dòng)物醫(yī)學(xué)學(xué)院，山西 太谷 030801；2.山西農(nóng)業(yè)大學(xué) 動(dòng)物科學(xué)學(xué)院，山西 太谷 030801)

養(yǎng)殖環(huán)境質(zhì)量是影響?zhàn)B豬生產(chǎn)的重要因素，豬舍內(nèi)有害氣體及懸浮粒子等環(huán)境污染物影響豬只生長(zhǎng)發(fā)育，誘發(fā)多種疾病，造成嚴(yán)重經(jīng)濟(jì)損失。豬舍內(nèi)空氣中的微生物(細(xì)菌、真菌、病毒等)會(huì)附著在懸浮粒子上，所形成微生物氣溶膠成為豬氣源性疾病的主要誘因。研究證實(shí)，如非洲豬瘟、豬瘟、豬鏈球菌病等可通過(guò)氣溶膠實(shí)現(xiàn)短距傳播，豬口蹄疫、豬偽狂犬病、豬支原體肺炎等也可經(jīng)氣溶膠長(zhǎng)距傳播。此外，如果攜帶人畜共患病病原體的氣溶膠向周邊居民區(qū)擴(kuò)散，會(huì)造成嚴(yán)重的生物安全風(fēng)險(xiǎn)。

微生物氣溶膠作為豬舍內(nèi)氣載污染物的重要來(lái)源，已引起國(guó)內(nèi)外研究人員廣泛關(guān)注。研究表明豬舍內(nèi)微生物氣溶膠和豬糞便中的微生物組成有較高相似性。隨著舍內(nèi)環(huán)境變化，氣溶膠中的微生物會(huì)快速適應(yīng)，產(chǎn)生時(shí)間及空間上的動(dòng)態(tài)變化模式。如豬舍類(lèi)型、飼養(yǎng)模式及季節(jié)變化等會(huì)導(dǎo)致舍內(nèi)空氣中細(xì)菌濃度存在差異。這些研究建立在傳統(tǒng)的分離培養(yǎng)方法之上，但其無(wú)法全面揭示豬舍內(nèi)微生物氣溶膠的豐富度和多樣性。隨著高通量測(cè)序技術(shù)的快速發(fā)展及應(yīng)用，為豬舍空氣樣本中微生物的多樣性檢測(cè)提供了技術(shù)支撐。本研究以山西省某典型規(guī)?；i場(chǎng)為研究對(duì)象，分別于不同季節(jié)對(duì)各類(lèi)型豬舍內(nèi)空氣樣本進(jìn)行采集，利用高通量測(cè)序技術(shù)對(duì)細(xì)菌的群落結(jié)構(gòu)及多樣性比較分析，探索氣溶膠中的細(xì)菌在空間和時(shí)間上的分布、變化特征，以期為豬舍內(nèi)空氣質(zhì)量控制和豬場(chǎng)生物安全防控等方面提供理論依據(jù)。

1 材料與方法

1.1 試驗(yàn)豬場(chǎng)及飼養(yǎng)管理

選擇山西省晉中市介休市某規(guī)?；i場(chǎng)開(kāi)展研究。該豬場(chǎng)存欄基礎(chǔ)母豬480余頭，其飼養(yǎng)規(guī)模、豬舍結(jié)構(gòu)、通風(fēng)模式及管理水平在山西省乃至中國(guó)北方地區(qū)均具有代表性。豬舍采取東西走向，長(zhǎng)約25 m，寬約8 m，高約3 m。保育豬舍共28欄，每欄6～8頭；育肥豬舍共10欄，每欄14～15頭；哺乳母豬舍存欄20頭母豬。半封閉式豬舍，舍頂部均勻分布3個(gè)通風(fēng)口，自然采光及通風(fēng)，冬季地暖保溫。豬只自由飲水，早晚人工清糞。

1.2 樣品采集與處理

分別于夏季(采樣時(shí)舍外溫度為27.3～34.0 ℃)、秋季(采樣時(shí)舍外溫度為18.3～18.7 ℃)、冬季(采樣時(shí)舍外溫度為-5.4～-1.5 ℃)進(jìn)行豬舍內(nèi)氣溶膠的采集和處理。BX、BQ、BD分別代表保育豬舍的夏季、秋季、冬季空氣樣本，F(xiàn)X、FQ、FD分別代表育肥豬舍的夏季、秋季、冬季空氣樣本，MX、MQ、MD分別代表哺乳母豬舍的夏季、秋季、冬季空氣樣本。選用六級(jí)空氣微生物采樣器采集豬舍空氣樣本(JWL-S6型，北京先能技術(shù)開(kāi)發(fā)有限責(zé)任公司，空氣流速率為28.3 L·min，采集時(shí)間為3 min)，選用玻纖膜(天津市領(lǐng)航實(shí)驗(yàn)設(shè)備有限公司)為采樣介質(zhì)，采樣點(diǎn)為豬舍中間過(guò)道的中點(diǎn)位置，距地面高度為1.2 m(豬背高處)，每個(gè)采樣點(diǎn)設(shè)3個(gè)重復(fù)，采樣前后均需對(duì)采樣器進(jìn)行消毒。玻纖膜空氣樣本采集后放入無(wú)菌凍存管中，并迅速投入液氮，帶回實(shí)驗(yàn)室，于-80 ℃冰箱保存。

1.3 DNA提取與測(cè)序

將剪碎的玻纖膜空氣樣本放入裂解液中充分裂解，待樣本裂解完畢后，采用德國(guó)MN公司的NucleoSpin 96 Soil (4×96)試劑盒進(jìn)行細(xì)菌基因組DNA的提取，并用DNA 濃度檢測(cè)儀器測(cè)定基因組DNA的濃度，將符合測(cè)樣要求的細(xì)菌基因組DNA樣本送百邁克生物科技有限公司開(kāi)展高通量測(cè)序。

以送檢DNA樣本作為PCR反應(yīng)模板，以338F：5'-ACTCCTACGGGAGGCAGCAG-3'和806R：5'-GGACTACHVGGGTWTCTAAT-3'作為引物，擴(kuò)增細(xì)菌基因組16S rRNA的v3～v4區(qū)域。將PCR產(chǎn)物純化、定量與均一化處理，形成測(cè)序文庫(kù)。選用Illumina HiSeq 2500平臺(tái)對(duì)質(zhì)檢合格文庫(kù)測(cè)序。

1.4 數(shù)據(jù)處理

優(yōu)化測(cè)序結(jié)果，將優(yōu)化序列(Tags)聚類(lèi)，劃分OTUs(Operational taxonomic units)，揭示豬舍空氣樣本的細(xì)菌群落結(jié)構(gòu)；分別采用Mothur(version v.1.30)軟件和QIIME軟件對(duì)送檢樣本進(jìn)行α多樣性(Alpha diversity)、β多樣性(Beta diversity)分析，比較不同季節(jié)及各類(lèi)型豬舍內(nèi)氣溶膠中細(xì)菌多樣性方面存在的差異。此外，使用SPSS 24.0進(jìn)行檢驗(yàn)。當(dāng)<0.05時(shí)，表示在95%的置信區(qū)間內(nèi)具有統(tǒng)計(jì)學(xué)上的顯著差異。

2 結(jié)果與分析

2.1 測(cè)序數(shù)據(jù)統(tǒng)計(jì)

選用Illumina Miseq高通量測(cè)序技術(shù)，分別測(cè)得不同季節(jié)、各類(lèi)型豬舍內(nèi)空氣樣本的序列數(shù)據(jù)。其中，共產(chǎn)生優(yōu)化序列1 524 972條，每個(gè)樣本至少產(chǎn)生43 145條優(yōu)化序列。以序列相似度97%為標(biāo)準(zhǔn)，將優(yōu)化序列劃分為OTUs，共產(chǎn)生OUTs為1 547個(gè)，不同樣本OTUs范圍為315～1 159個(gè)。

2.2 豬舍空氣中細(xì)菌α多樣性分析

對(duì)豬舍內(nèi)空氣樣本開(kāi)展α多樣性指數(shù)分析，結(jié)果見(jiàn)表1，不同季節(jié)、各類(lèi)型豬舍內(nèi)的空氣樣本均獲得99%以上的高覆蓋度，反映測(cè)序結(jié)果的全面性。在圖1中，隨著樣品序列數(shù)量的增加，Shannon指數(shù)稀釋曲線逐漸趨于平坦，OTUs不再隨著測(cè)序數(shù)量增加而增長(zhǎng)。說(shuō)明測(cè)序數(shù)量能夠代表空氣樣本中的大多數(shù)細(xì)菌信息，反映測(cè)序結(jié)果的真實(shí)性。

圖1 樣本Shannon指數(shù)稀釋曲線

由表1可以看出：各類(lèi)型豬舍空氣中細(xì)菌的豐富度(Chao1指數(shù)、ACE指數(shù))和多樣性(Shannon指數(shù)、Simpson指數(shù))均表現(xiàn)出明顯的季節(jié)性變化趨勢(shì)，冬季>秋季>夏季。其中，育肥豬舍秋、冬兩季更替時(shí)菌群豐富度和多樣性指數(shù)變化顯著，冬季Chao1、ACE指數(shù)極顯著大于夏秋季(<0.01)，Shannon指數(shù)顯著大于夏秋季，Simpson指數(shù)顯著小于夏秋季(<0.05)；哺乳母豬舍夏、秋兩季更替時(shí)菌群豐富度和多樣性指數(shù)變化顯著，秋季Chao1、ACE、Shannon指數(shù)極顯著大于夏季，Simpson指數(shù)極顯著小于夏季(<0.01)。此外，由表2可見(jiàn)，秋季育肥豬舍菌群豐富度和多樣性指數(shù)均極顯著小于哺乳母豬舍(<0.01)；夏季保育豬舍菌群豐富度和多樣性指數(shù)略大于育肥豬舍、哺乳母豬舍，但差異不顯著(>0.05)。

表1 不同季節(jié)豬舍空氣樣本Alpha多樣性指數(shù)統(tǒng)計(jì)

表2 不同類(lèi)型豬舍空氣樣本Alpha多樣性指數(shù)統(tǒng)計(jì)

2.3 豬舍空氣中細(xì)菌群落結(jié)構(gòu)差異分析

為進(jìn)一步研究豬舍類(lèi)型、季節(jié)變化對(duì)豬舍空氣中細(xì)菌氣溶膠群落結(jié)構(gòu)的影響，利用非度量多維尺度分析法(NMDS)對(duì)不同樣本的菌群結(jié)構(gòu)差異進(jìn)行分析(圖2)。在所有豬舍空氣樣本中，冬季各樣本間距離較近，且單獨(dú)聚為一簇，而夏、秋季的樣本相互重疊。這說(shuō)明冬季各類(lèi)型豬舍間空氣樣本菌群結(jié)構(gòu)差異較小，且冬季樣本菌群特征明顯區(qū)分于其他季節(jié)。NMDS 分析的 stress 值為0.028，表明排序結(jié)果可信。

圖2 樣本NMDS分析圖

2.4 豬舍空氣中細(xì)菌群落結(jié)構(gòu)組成分析

在不同季節(jié)，各類(lèi)型豬舍空氣樣本中共獲得細(xì)菌分為23個(gè)門(mén)，547個(gè)屬。

從表3可以看出：在門(mén)水平下，豬舍空氣中細(xì)菌主要由變形菌門(mén)(Proteobacteria)、厚壁菌門(mén)(Firmicutes)、擬桿菌門(mén)(Bacteroidetes)、放線菌門(mén)(Actinobacteria)、酸桿菌門(mén)(Acidobacteria)等組成。優(yōu)勢(shì)菌門(mén)的相對(duì)豐度值均表現(xiàn)出季節(jié)性差異，秋、冬季節(jié)更替時(shí)變化明顯，夏、秋季節(jié)更替時(shí)則差異均不顯著(>0.05)。在保育豬舍，冬季變形菌門(mén)(Proteobacteria)、酸桿菌門(mén)(Acidobacteria)相對(duì)豐度值均顯著大于夏、秋季(<0.05)；芽單胞菌門(mén)(Gemmatimonadetes)相對(duì)豐度值極顯著大于夏、秋季(<0.01)；厚壁菌門(mén)(Firmicutes) 相對(duì)豐度值顯著小于秋季(<0.05)，極顯著小于夏季(<0.01)。在育肥豬舍，冬季變形菌門(mén)(roteobacteria)相對(duì)豐度值顯著大于夏季(<0.05)；酸桿菌門(mén)(Acidobacteria)相對(duì)豐度值顯著大于夏、秋季(<0.05)；芽單胞菌門(mén)(Gemmatimonadetes)相對(duì)豐度值顯著大于秋季(<0.05)，極顯著大于夏季(<0.01)；厚壁菌門(mén)(Firmicutes) 相對(duì)豐度值顯著小于夏、秋季(<0.05)；在哺乳母豬舍，冬季變形菌門(mén)(roteobacteria)顯著大于夏季(<0.05)；酸桿菌門(mén)(Acidobacteria)、綠彎菌門(mén)(Chloroflexi)相對(duì)豐度值均極顯著大于夏、秋季(<0.01)；芽單胞菌門(mén)(Gemmatimonadetes)相對(duì)豐度值顯著大于秋季(<0.05)，極顯著大于夏季(<0.01)；厚壁菌門(mén)(Firmicutes) 相對(duì)豐度值顯著小于夏、秋季(<0.05)。此外，從豬舍類(lèi)型分析(表4)，各類(lèi)型豬舍間細(xì)菌群落結(jié)構(gòu)相似度較高。在秋季，保育豬舍藍(lán)藻細(xì)菌門(mén)(Cyanobacteria)相對(duì)豐度值顯著小于哺乳母豬舍(<0.05)，其他優(yōu)勢(shì)菌門(mén)的相對(duì)豐度值均差異不顯著(>0.05)。

表3 不同季節(jié)豬舍空氣樣本門(mén)水平上的優(yōu)勢(shì)菌群占比

表4 不同類(lèi)型豬舍空氣樣本門(mén)水平上的優(yōu)勢(shì)菌群占比

從屬分類(lèi)水平挑選豐富度最高的前20個(gè)菌群作物種豐富熱圖。從圖3可以看出：豬舍空氣中細(xì)菌群落結(jié)構(gòu)表現(xiàn)出明顯的季節(jié)性差異。在夏、秋季，___、乳酸桿菌屬()、葡萄球菌屬()為優(yōu)勢(shì)菌屬。在冬季，鞘脂單胞菌屬()、____6、___為優(yōu)勢(shì)菌屬。此外，從豬舍類(lèi)型分析，各類(lèi)型豬舍間細(xì)菌群落結(jié)構(gòu)相似度較高。

圖3 屬水平上的物種豐度熱圖

將豬舍空氣中細(xì)菌相對(duì)豐度較高(在屬分類(lèi)水平豐富度最高的前50個(gè)物種，且相對(duì)豐度值大于0.1%)的潛在動(dòng)物病原菌屬進(jìn)行整理，結(jié)果見(jiàn)表5。從表5可以看出：病原菌屬相對(duì)豐度值均表現(xiàn)出明顯的季節(jié)性差異，夏、秋季高于冬季。其中，在保育豬舍，夏季埃希氏-志賀菌屬(-)、不動(dòng)桿菌屬()相對(duì)豐度值顯著大于冬季(<0.05)，秋季棒狀桿菌屬(_1)、鏈球菌屬()相對(duì)豐度值顯著大于冬季(<0.05)；在育肥豬舍，夏季鏈球菌屬()、不動(dòng)桿菌屬()相對(duì)豐度值顯著大于冬季(<0.05)；在哺乳母豬舍，夏季埃希氏-志賀菌屬(-)、鏈球菌屬()、葡萄球菌屬()、棒狀桿菌屬(_1)相對(duì)豐度值顯著大于冬季(<0.05)，秋季埃希氏-志賀菌屬(-)、不動(dòng)桿菌屬()、假單胞菌屬()、梭狀芽孢桿菌屬(___1)相對(duì)豐度值顯著大于冬季(<0.05)。此外，從豬舍類(lèi)型分析(見(jiàn)表6)，在夏季，各類(lèi)型豬舍間葡萄球菌屬()相對(duì)豐度值均維持較高的水平，彼此差異不顯著(>0.05)；在夏季，育肥豬舍埃希氏-志賀菌屬(-)相對(duì)豐度值高于保育豬舍和哺乳母豬舍，但差異不顯著(>0.05)。

表5 不同季節(jié)豬舍空氣樣本優(yōu)勢(shì)菌屬中潛在病原菌分布情況

表6 不同類(lèi)型豬舍空氣樣本優(yōu)勢(shì)菌屬中潛在病原菌分布情況

3 討 論

3.1 豬舍空氣中細(xì)菌多樣性分析

國(guó)內(nèi)前期關(guān)于豬舍內(nèi)微生物氣溶膠豐度變化規(guī)律的研究多基于傳統(tǒng)培養(yǎng)法。有研究表明豬舍類(lèi)型、季節(jié)變化對(duì)氣溶膠中可培養(yǎng)的總需氧菌及致病菌數(shù)量有明顯影響。本人前期研究表明，保育豬舍總需氧菌濃度季節(jié)性變化趨勢(shì)為冬季顯著高于其他季節(jié)，夏季在部分時(shí)間點(diǎn)高于春、秋季。與傳統(tǒng)培養(yǎng)法相比，高通量測(cè)序技術(shù)在分析微生物氣溶膠的多樣性及豐度方面有著天然的優(yōu)勢(shì)。本研究利用高通量測(cè)序技術(shù)獲得豬舍空氣樣本總OTU數(shù)為1 547，不同樣本OTU數(shù)范圍為315～1 159。從季節(jié)變化上分析，豬舍空氣中細(xì)菌豐富度和多樣性指數(shù)表現(xiàn)出明顯的季節(jié)性差異，冬季>秋季>夏季。有研究表明分娩豬舍內(nèi)小氣候變化特別是舍內(nèi)風(fēng)速等環(huán)境因素對(duì)細(xì)菌氣溶膠的豐富度和多樣性有一定的影響。本研究試驗(yàn)豬場(chǎng)位于山西省中部盆地，為典型大陸季風(fēng)氣候，四季分明，秋冬季早晚溫差大。各類(lèi)型豬舍均采取自然通風(fēng)模式。在夏季，豬舍通過(guò)開(kāi)窗通風(fēng)實(shí)現(xiàn)降溫；進(jìn)入秋季，天氣轉(zhuǎn)涼，哺乳母豬舍通風(fēng)次數(shù)明顯減少，舍內(nèi)空氣中細(xì)菌豐富度及多樣性指數(shù)顯著升高；進(jìn)入冬季，育肥豬舍也采取密閉保溫模式，其豐富度和多樣性指數(shù)也顯著升高。此外，從豬舍類(lèi)型上分析，在夏、秋季，各類(lèi)型豬舍間菌群豐富度和多樣性均存在明顯差異。有研究結(jié)果表明，除了通風(fēng)因素外，舍內(nèi)其他環(huán)境因素的變化也會(huì)影響細(xì)菌的多樣性，即在相對(duì)封閉條件下，豬群飼養(yǎng)密度、舍內(nèi)環(huán)境衛(wèi)生等可能是影響豬舍空氣中細(xì)菌豐富度及多樣性的重要因素。

3.2 豬舍空氣中優(yōu)勢(shì)菌群分析

通過(guò)細(xì)菌群落結(jié)構(gòu)分析，在不同季節(jié)，各類(lèi)型豬舍空氣樣本中共獲得23個(gè)門(mén)，547個(gè)屬。在門(mén)水平下，變形菌門(mén)(Proteobacteria)、厚壁菌門(mén)(Firmicutes)、擬桿菌門(mén)(Bacteroidetes)、放線菌門(mén)(Actinobacteria)、酸桿菌門(mén)(Acidobacteria)等為豬舍微生物氣溶膠樣本中的優(yōu)勢(shì)菌門(mén)。這與國(guó)內(nèi)相關(guān)報(bào)道較為一致。在屬水平下，豬舍空氣中的優(yōu)勢(shì)菌屬種類(lèi)及相對(duì)豐度存在季節(jié)性差異，夏、秋季優(yōu)勢(shì)菌屬為_(kāi)__、乳酸桿菌屬()、葡萄球菌屬()等、冬季優(yōu)勢(shì)菌屬為鞘脂單胞菌屬()、____6、___等。張喜慶在豬舍內(nèi)氣溶膠樣本中發(fā)現(xiàn)梭狀芽孢桿菌屬、乳酸桿菌屬、尤微菌屬為主要優(yōu)勢(shì)菌屬，同本研究結(jié)果存在差異，這可能與不同地區(qū)地理、氣候環(huán)境或試驗(yàn)豬場(chǎng)在飼養(yǎng)方式、營(yíng)養(yǎng)配方上存在的差異相關(guān)。此外，不同類(lèi)型豬舍內(nèi)環(huán)境存在差異，導(dǎo)致豬舍內(nèi)菌群結(jié)構(gòu)不同，本研究也發(fā)現(xiàn)，各類(lèi)型豬舍間空氣中優(yōu)勢(shì)菌屬的相對(duì)豐度值存在差異，但大部分差異不顯著。

在豐度水平前50的細(xì)菌群落中發(fā)現(xiàn)文獻(xiàn)中已報(bào)道的病原菌，包括不動(dòng)桿菌屬()、梭狀芽孢桿菌屬(___1)、棒狀桿菌屬(_1)、埃希氏-志賀菌屬(-)、假單胞菌屬()、葡萄球菌屬()、鏈球菌屬()等。病原菌不僅對(duì)畜禽以及養(yǎng)殖場(chǎng)工作人員的健康產(chǎn)生影響，也給養(yǎng)殖場(chǎng)的生物安全防控帶來(lái)巨大挑戰(zhàn)。其中，鏈球菌通過(guò)粘膜感染引起的敗血癥、腦膜炎、關(guān)節(jié)炎甚至急性死亡，同時(shí)繼發(fā)多種呼吸道炎癥，是一種重要的人畜共患病病原。梭狀芽孢桿菌屬中的產(chǎn)氣莢膜梭菌和埃希氏-志賀菌屬中的大腸桿菌是引起仔豬細(xì)菌性腹瀉的常見(jiàn)病原體。葡萄球菌作為一種條件致病菌廣泛存在于畜禽舍空氣內(nèi)、動(dòng)物體表，可誘發(fā)仔豬滲出性皮炎、呼吸道感染、腦膜炎、敗血癥等多種疾病。本研究發(fā)現(xiàn)豬舍空氣樣本中優(yōu)勢(shì)病原菌屬相對(duì)豐度值的季節(jié)性變化趨勢(shì)較一致，表現(xiàn)為夏、秋季明顯高于冬季，這與樣本的菌群豐富度及多樣性變化趨勢(shì)相反。雖然秋、冬季交替后病原菌屬的相對(duì)含量明顯下降，但考慮到冬季冷熱應(yīng)激以及有害氣體、粉塵濃度升高等因素造成的不利影響，豬群相對(duì)易感，此時(shí)仍需加強(qiáng)微生物氣溶膠特別是病原菌的有效控制，切斷疾病的傳播渠道，降低豬場(chǎng)生物安全風(fēng)險(xiǎn)。

4 結(jié) 論

豬舍空氣中細(xì)菌豐富度和多樣性表現(xiàn)出季節(jié)性變化趨勢(shì)，冬季>秋季>夏季；在夏、秋季，豬舍類(lèi)型也是影響細(xì)菌多樣性的重要因素。豬舍空氣中細(xì)菌群落結(jié)構(gòu)表現(xiàn)為季節(jié)性差異，秋、冬季節(jié)更替時(shí)變化明顯，各類(lèi)型豬舍間細(xì)菌群落結(jié)構(gòu)相似度較高。在豬舍空氣樣本的優(yōu)勢(shì)菌屬中發(fā)現(xiàn)大量潛在動(dòng)物病原菌，并具有較明顯的季節(jié)性變化特征，豬場(chǎng)應(yīng)加強(qiáng)相應(yīng)的生物安全措施。