寧夏部分地區(qū)新生犢牛腹瀉的病原調(diào)查及主要病原的遺傳進(jìn)化分析

郭仕輝，趙清梅，陳紹淑，張浩東，田新岳，張津慎，許立華，馬 云，余永濤

(1.寧夏大學(xué) 農(nóng)學(xué)院，銀川 750021；2.北方民族大學(xué) 生物科學(xué)與工程學(xué)院，銀川 750021；3.國家民委發(fā)酵釀造工程生物技術(shù)重點(diǎn)實(shí)驗(yàn)室，銀川 750021)

犢牛腹瀉是犢牛發(fā)生的以消化不良、腹瀉為主要臨床特征的消化道疾病綜合征。犢牛腹瀉多發(fā)生于斷奶前的新生哺乳犢牛，故通常將哺乳期犢牛腹瀉稱為新生犢牛腹瀉(Neonatal calf diarrhea，NCD)[1]。NCD是造成斷奶前犢牛死亡的最主要原因之一。據(jù)報(bào)道，美國NCD的死亡率為57%，韓國為53.4%[2]，給奶牛和肉牛養(yǎng)殖業(yè)帶來了巨大的經(jīng)濟(jì)損失。此外，NCD還會造成后備母牛生長緩慢、生產(chǎn)性能和免疫力低下，飼養(yǎng)和治療成本增加，嚴(yán)重制約牛場的良種選育和牛群擴(kuò)繁[3]。研究表明，中國犢牛腹瀉的發(fā)病率為 9.6%～60.0%，死亡率大約為1%，但大多研究未對NCD和斷奶后犢牛腹瀉進(jìn)行嚴(yán)格區(qū)分，NCD的發(fā)病情況和死亡率未明確報(bào)道[4]。哺乳期犢牛和斷奶后犢牛在生理機(jī)能和消化道解剖結(jié)構(gòu)、飼養(yǎng)管理方式等方面有很大不同，這也決定了NCD和斷奶后犢牛腹瀉的發(fā)病特點(diǎn)、病原流行特點(diǎn)會有所差別，對NCD進(jìn)行深入系統(tǒng)地研究有助于犢牛腹瀉的早期防控[5]。

腸道病原感染是造成犢牛腹瀉的主要原因，流行病學(xué)研究表明牛輪狀病毒(Bovine Rotavirus，BRV)、牛冠狀病毒(Bovine Coronavirus，BCoV)、K99型大腸桿菌和隱孢子蟲(Cryptosporidium)是引起犢牛腹瀉的最常見病原[6]。此外，沙門氏菌(Salmonella)、產(chǎn)氣莢膜梭菌(Clostridiumperfringen)、牛病毒性腹瀉病毒(Bovine viral diarrhea virus，BVDV)、十二指腸賈第鞭毛蟲(Giardiaduodenalis)等感染也可導(dǎo)致犢牛腹瀉[7]。近年研究發(fā)現(xiàn)，在美國、韓國、埃及等數(shù)十個國家腹瀉犢牛的糞便中檢出或分離到牛諾如病毒(Bovine Norovirus，BNoV)和紐布病毒(Nebovirus, NeV)，這兩種病毒可造成腸道黏膜上皮細(xì)胞的損傷，可能也是造成犢牛腹瀉的重要病原[8-10]。自2018年以來，中國多地也檢測到BNoV，但BNoV在寧夏地區(qū)的流行情況尚不清楚[11-13]。目前寧夏地區(qū)關(guān)于犢牛腹瀉相關(guān)病原的流行病學(xué)研究較少，王文佳等[14]、高?；鄣萚15]對寧夏地區(qū)腹瀉犢牛主要病原進(jìn)行了流行病學(xué)調(diào)查，但沒有明確斷奶前犢牛的病原檢出情況，且未對BCoV、BNoV等其他相關(guān)病原進(jìn)行調(diào)查。

為了解寧夏部分地區(qū)NCD相關(guān)病原的流行情況及主要病原的遺傳進(jìn)化規(guī)律，本研究擬對寧夏地區(qū)主要奶牛養(yǎng)殖區(qū)域的部分中大型規(guī)?；鲋?月齡內(nèi)的斷奶前犢牛進(jìn)行腹瀉相關(guān)病原檢測及主要病原的遺傳進(jìn)化關(guān)系研究，以期為寧夏地區(qū)犢牛腹瀉的早期防控提供參考。

1 材料與方法

1.1 材 料

1.1.1 試驗(yàn)材料 2019年7月至2019年9月對寧夏銀川市、吳忠市、石嘴山市主要奶牛養(yǎng)殖園區(qū)的9個中大型規(guī)?；膛?月齡內(nèi)的斷奶前犢牛糞便樣本進(jìn)行隨機(jī)采集，共采集腹瀉犢牛糞便106份和健康犢牛糞便89份。詳見表1。

表1 采樣情況

根據(jù)威斯康星大學(xué)麥迪遜分校的糞便評分標(biāo)準(zhǔn)對犢牛是否發(fā)生腹瀉進(jìn)行評判[16]。糞便成型、質(zhì)地良好為0分，為健康犢牛。半成型、糞便稀軟為1分，為健康犢牛。糞便稀軟、不成型為2分，為腹瀉犢牛。糞便完全成液體狀為3分，為腹瀉犢牛。

1.1.2 主要試劑 RNA提取試劑盒(RNAiso Plus)、反轉(zhuǎn)錄試劑盒(PrimescriptTM1st Strand NCDNA Synthesis Kit)、PremixTaqTM、瓊脂糖凝膠回收試劑盒(TaKaRa MiniBEST Agarose Gel DNA Extraction Kit Ver.4.0)、DL 1000 DNA Marker均購自寶日醫(yī)生物技術(shù)(北京)有限公司；糞便DNA提取試劑盒(Stool DNA Kit)購自O(shè)MEGA公司。

1.1.3 主要儀器 凝膠成像系統(tǒng)(Bio-Rad，ChemiDoc XRS+)、PCR儀(Bio-Rad，S1000 Thermal Cycler)、臺式高速冷凍離心機(jī)(Eppendorf，5430)。

1.2 方 法

1.2.1 樣本采集 首先對待采集犢牛進(jìn)行日齡等背景信息統(tǒng)計(jì)，然后進(jìn)行糞便樣本采集。采用Gomez等[16]報(bào)道的腸道侵入法對犢牛直腸糞便進(jìn)行采集。用酒精棉球?qū)Ω亻T周邊區(qū)域進(jìn)行清潔和消毒，佩戴一次性無菌橡膠手套并將食指伸入直腸內(nèi)輕輕刺激腸壁直至犢牛排出糞便于掌心。將糞便轉(zhuǎn)移至一次性無菌糞便收集杯并置于冰盒中，帶回實(shí)驗(yàn)室于-80 ℃保存，用于病原檢測。

1.2.2 PCR檢測 首先按照糞便DNA提取試劑盒說明書提取DNA，然后分別采用K99型大腸桿菌和隱孢子蟲的特異性引物進(jìn)行PCR擴(kuò)增[17-18]。K99型大腸桿菌 PCR擴(kuò)增體系為Premix 11.25 μL，DNA 1 μL，上下游引物各0.5 μL，dd H2O 11.75 μL。反應(yīng)條件為94 ℃預(yù)變性5 min；94 ℃變性30 s，58 ℃退火30 s，72 ℃延伸50 s，共35個循環(huán)；72 ℃延伸10 min，4 ℃保存。隱孢子蟲檢測采用巢氏PCR，第一輪擴(kuò)增體系和反應(yīng)體系同K99型大腸桿菌，退火溫度為55 ℃；第二輪擴(kuò)增模板為第一輪PCR產(chǎn)物，擴(kuò)增體系和反應(yīng)體系同K99型大腸桿菌，退火溫度為58 ℃。

另將糞便樣本與PBS溶液按1∶5的體積比充分渦旋混勻，-80 ℃冰箱反復(fù)凍融3次，3 000 r/min離心10 min，棄沉淀，上清12 000 r/min離心30 min，離心后取上清用于RNA提取。根據(jù)RNA提取試劑盒(RNAiso Plus)操作程序?qū)μ幚順颖具M(jìn)行RNA提取。根據(jù)PrimescriptTM反轉(zhuǎn)錄試劑盒操作程序?qū)⑻崛〉目俁NA反轉(zhuǎn)錄為cDNA。分別參照BRVVP6基因、BCoVN蛋白[19]基因、BNoVRdRp基因[12]、NeVRdRp基因[11]和BVDV5′UTR基因[20]設(shè)計(jì)并合成引物。以cDNA為模板進(jìn)行PCR擴(kuò)增。最后取PCR產(chǎn)物經(jīng)瓊脂糖凝膠電泳和凝膠成像系統(tǒng)進(jìn)行檢測和分析。引物信息詳見表2。

表2 引物信息

1.2.3 主要病原的遺傳進(jìn)化分析 根據(jù)檢測結(jié)果，對BRV、BCoV、BNoV和隱孢子蟲等主要病原進(jìn)行遺傳進(jìn)化分析。將PCR產(chǎn)物經(jīng)瓊脂糖凝膠回收試劑盒回收，送至上海生工生物有限公司雙向測序。所得序列采用SeqMan(DNAstar 7.0)軟件進(jìn)行序列拼接，在NCBI數(shù)據(jù)庫(https://blast.ncbi.nlm.nih.gov)中進(jìn)行BLAST比對。應(yīng)用DNAstar軟件中的MegAlign程序?qū)LAST比對中序列一致性較高的同源序列進(jìn)行Align分析，根據(jù)Align結(jié)果，應(yīng)用Mega 7.0軟件采用Neighbor-Joining法構(gòu)建系統(tǒng)發(fā)育樹。

1.2.4 統(tǒng)計(jì)分析 應(yīng)用EXCEL 2016軟件對病原檢測結(jié)果進(jìn)行統(tǒng)計(jì)分析并繪圖。

2 結(jié)果與分析

2.1 腹瀉情況統(tǒng)計(jì)

對隨機(jī)采集的106份腹瀉犢牛日齡信息進(jìn)行統(tǒng)計(jì)，結(jié)果如表3所示。11～20日齡的腹瀉犢牛占比最高(33.0%)，其次是1～10日齡(31.1%)和21～30日齡(17.9%)，其他日齡的腹瀉犢牛占比較低，均小于10%。

表3 腹瀉犢牛日齡信息分布情況

2.2 病原檢測

2.2.1 主要病原檢出結(jié)果 對PCR檢測結(jié)果進(jìn)行統(tǒng)計(jì)，結(jié)果如表4所示。腹瀉犢牛BRV、BCoV、BNoV、NeV、K99型大腸桿菌和隱孢子蟲檢出率分別為19.8%、9.43%、12.26%、0.94%、2.83%和6.60%，BVDV檢出率為0。健康犢牛BCoV、BNoV、NeV、K99型大腸桿菌和隱孢子蟲檢出率分別為5.62%、5.62%、1.12%、1.12%和3.37%，BRV和BVDV檢出率為0。

表4 主要病原檢出結(jié)果

2.2.2 主要病原感染情況 各牛場腹瀉犢牛主要病原檢出結(jié)果如圖1所示。檢出2種以上病原的牛場有7個，檢出3種以上病原的牛場有5個，檢出4種以上病原的牛場有2個。BRV、BCoV、BNoV、NeV、BVDV、K99型大腸桿菌和隱孢子蟲的場群流行率分別為77.8%(7/9)、44.5%(4/9)、55.6%(5/9)、11.1%(1/9)、0(0/9)、11.1%(1/9)和33.3%(3/9)。其中BRV牛場檢出率最高，其次是BNoV、BCoV和隱孢子蟲，最后是NeV和K99型大腸桿菌，BVDV未檢出。BCoV檢出主要集中在牛場F4和F7，BNoV檢出主要集中在牛場F5和F8，隱孢子蟲檢出主要集中在牛場F7、F5和F1，NeV僅在F3牛場檢出，K99型大腸桿菌僅在F7牛場檢出。

圖1 9個牛場腹瀉犢牛NCD相關(guān)病原檢出情況

對腹瀉犢牛中BRV、BCoV、BNoV和隱孢子蟲等4種檢出率較高的病原進(jìn)行混合感染情況統(tǒng)計(jì)，結(jié)果如表5所示，單一病原感染的犢牛占 30.2%，2種病原混合感染的犢牛占7.55%、3 種病原混合感染的犢牛占0.94%，未出現(xiàn)4種病原混合感染的情況。腹瀉犢牛主要為單一病原感染，其中BRV(15.1%)、BNoV(10.4%)和BCoV(4.72%)單一感染率較高，隱孢子蟲在混合感染患牛中檢出率較高。

表5 腹瀉犢?；旌细腥厩闆r

2.2.3 不同日齡犢牛的病原檢出結(jié)果 將腹瀉犢牛中檢出率較高的BRV、BCoV、BNoV和隱孢子蟲共4種主要病原檢出結(jié)果按照日齡進(jìn)行分類，結(jié)果如表6所示。感染BRV、BCoV和隱孢子蟲的腹瀉犢牛主要在30日齡以內(nèi)，感染BNoV的腹瀉犢牛主要在20日齡以內(nèi)。

表6 不同日齡腹瀉犢牛中主要病原檢出情況

2.3 NCD主要病原的遺傳進(jìn)化分析

2.3.1 BRV 對21個BRV陽性樣品的VP6基因PCR產(chǎn)物進(jìn)行測序。測序結(jié)果表明，獲得的21個序列中有10個核苷酸序列不一致，分別為D2、D6、D8、F3、F8、K12、O2、O6、P2和Y7。VP6基因BLAST分析結(jié)果表明，所有序列都與GeneBank中已公布的輪狀病毒A群VP6基因同源性最高，核苷酸序列一致性為88.8%～99.3%。遺傳進(jìn)化分析結(jié)果如圖2所示，獲得的10個BRVVP6基因與GeneBank中已公布的BRV A群遺傳進(jìn)化距離最近。其中D2、D6、D8、F3、F8、K12、O2、O6和P2與中國株(Bo/HB-BD/CH)聚在同一個分支上，Y7株與美國株(Hu/2011825363/USA)單獨(dú)聚為一支。

圖2 基于BRV VP6基因序列構(gòu)建的系統(tǒng)發(fā)育樹

2.3.2 BCoV 對15個BCoV陽性樣品的N蛋白基因PCR產(chǎn)物進(jìn)行測序，結(jié)果表明，獲得的15個序列中有9個核苷酸序列不一致，分別為C1、D5、D7、E1、E2、I5、O3、O9和O10。BLAST分析結(jié)果表明，所有序列都與GeneBank中已公布的牛冠狀病毒N蛋白基因同源性最高，核苷酸序列一致性為98.9%～99.9%。遺傳進(jìn)化分析結(jié)果如圖3所示，獲得的9個BCoV N蛋白基因全部與GeneBank中已公布的BCoV N蛋白基因遺傳進(jìn)化距離較近。9個序列與已報(bào)道的中國株(Bo/SC1/CH、Bo/LN1/CH)和美國株(Bo/4-17-08/USA)聚在一個大的分支，但C1、D5、D7、E1、O3、O9和O10在該分支中單獨(dú)聚在一個小的分 支內(nèi)。

圖3 基于BCoV N蛋白區(qū)域構(gòu)建的系統(tǒng)發(fā)育樹

2.3.3 BNoV 對18個BNoV陽性樣品的RdRp基因PCR產(chǎn)物進(jìn)行測序，結(jié)果表明，獲得的18個序列中有12個核苷酸序列不一致，分別為A6、F7、F9、F11、F12、I1、I6、I7、I8、I10，J5和J6。BLAST分析結(jié)果表明，所有序列都與GeneBank中已公布的BNoV GⅢ.2型RdRp基因的同源性最高，核苷酸序列一致性為87.6%～ 98.4%。遺傳進(jìn)化分析結(jié)果如圖4所示，獲得的12個BNoVRdRp基因與GeneBank中已公布的GⅢ.2型BNoV遺傳進(jìn)化距離最近。其中A6、F7、F9、F11、F12、I1、I6、I7、I8和I10等與中國株(Bo/SCZ-18/CH)等聚在同一個分支上，J5與韓國株(Bo/MG-16/KOR)聚為一支，J6單獨(dú)聚為一支，但與中國株(Bo/MI-48/CH)遺傳進(jìn)化距離較近。

圖4 基于BNoV RdRp區(qū)域構(gòu)建的系統(tǒng)發(fā)育樹

2.3.4 隱孢子蟲 對10個陽性樣本的隱孢子蟲SSU區(qū)域PCR產(chǎn)物進(jìn)行測序，測序結(jié)果表明，獲得的10個序列中有6個核苷酸序列不一致，分別為D1、D5、F2、H4、W2和W9。BLAST分析結(jié)果表明，D1、D5、F2和W9序列與GeneBank中已公布的微小隱孢子蟲(Cryptosporidiumparvum)SSU基因的同源性最高，核苷酸序列一致性為95.7%～98.4%；H4、W2與GeneBank中已公布的牛隱孢子蟲(Cryptosporidiumbovis)SSU基因的同源性最高，核苷酸序列一致性為98.5%和99.5%。遺傳進(jìn)化分析結(jié)果如圖5所示，D1、D5、F2和W9樣本序列與GeneBank中已公布的微小隱孢子蟲遺傳進(jìn)化距離最近，H4、W2與GeneBank中已公布的牛隱孢子蟲遺傳進(jìn)化距離最近。其中D1、D5、F2、W9等與微小隱孢子蟲日本株(Ca/NCD01/JAN)、伊朗株(Bo/CPOA5/IRA)、中國株(Bo/245/CH)單獨(dú)聚為一支；W2與牛隱孢子蟲中國株(Bo/C90/CH、Ya/DR-2/CH)聚為一支；H4與牛隱孢子蟲土耳其株(Bo/AF_CR07/TUR)、中國株(Bo/ECOST23057/CH、Bo/151/CH和Bo/GS-15/CH)、孟加拉國株(Bo/GC2/BAN)聚為一支。

圖5 基于隱孢子蟲SSU區(qū)域的系統(tǒng)發(fā)育樹

3 討 論

3.1 新生犢牛腹瀉的日齡特點(diǎn)

新生犢牛由于胃腸道結(jié)構(gòu)發(fā)育不成熟、腸道物理屏障和生物屏障不完整、機(jī)體免疫機(jī)能不完善等原因極易受到病原體侵襲，進(jìn)而發(fā)生腹瀉[21]。新生犢牛日糧以牛奶為主，在犢牛出生后，如果初乳飼喂不及時、常乳消毒不徹底或飲水不清潔等都可引起NCD。此外，犢牛攝入的牛奶可為感染性病原(如BRV等)提供良好的腸道生存環(huán)境，有利于感染性病原的生存[22]。據(jù)美國乳業(yè)國家動物健康監(jiān)測系統(tǒng)報(bào)告稱，NCD多發(fā)生于1月齡以下犢牛[23]。Foster等[5]研究表明，腹瀉多發(fā)于斷奶前犢牛，且以1月齡以內(nèi)犢牛最為常見。本研究中腹瀉犢牛主要集中在1月齡內(nèi)，日齡越小腹瀉犢牛占比越高，與Foster等[5]研究結(jié)果一致。

3.2 寧夏NCD主要病原的流行特點(diǎn)

王文佳等[14]研究表明，寧夏地區(qū)腹瀉犢牛的BRV、K99型大腸桿菌、隱孢子蟲和賈第鞭毛蟲等的檢出率分別為32.28%、11.02%、25.20%和4.72%，腹瀉犢牛以單一病原感染為主。李復(fù)煌[24]研究表明，北京地區(qū)犢牛腹瀉感染類型以單一感染為主，輪狀病毒和隱孢子蟲是最主要的混合感染類型。這與本研究檢測結(jié)果部分一致。國內(nèi)外普遍認(rèn)為引起犢牛腹瀉的主要病原為BRV、BCoV、K99型大腸桿菌和隱孢子蟲，而本研究結(jié)果表明寧夏地區(qū)引起NCD的主要病原為BRV、BNoV、BCoV和隱孢子蟲，檢出率分別為 19.8%、12.26%、9.43%和6.6%。BNoV是近年來新發(fā)現(xiàn)的一種病原，Cho等[21]研究表明，BNoV可侵入空腸、回腸和結(jié)腸遠(yuǎn)端上皮細(xì)胞，引起腸道上皮隱窩和小腸絨毛細(xì)胞萎縮、壞死，減少腸道可吸收面積，進(jìn)而引起腹瀉。王玥琳等[12-13]研究表明，腹瀉糞便樣本中BNoV的檢出率顯著高于健康糞便樣本，中國四川、遼寧、河南、山東、陜西等地腹瀉犢牛BNoV的檢出率分別為 34.78%、75%、11.11%、9.09%、0。本研究中，半數(shù)以上的牧場中均檢出BNoV，且檢出率較高，表明BNoV已經(jīng)在寧夏奶牛場中廣泛流行。

中國犢牛腹瀉發(fā)病率為9.6%～60.0%[4]，但已有的研究多未對斷奶前和斷奶后的犢牛腹瀉進(jìn)行嚴(yán)格區(qū)分，斷奶前犢牛腹瀉的病原感染情況并未明確。實(shí)際生產(chǎn)中，斷奶前犢牛與斷奶后犢牛的飼養(yǎng)方式和養(yǎng)殖環(huán)境明顯不同，斷奶前犢牛均在獨(dú)立牛棚單獨(dú)飼喂，接觸病原的機(jī)會相對較低，且初乳飼喂為新生犢牛的早期階段提供了有效的被動免疫。而斷奶后犢牛一般混群飼喂，接觸病原機(jī)會增多，腸道病原組成可能更為復(fù)雜。因此，斷奶前與斷奶后犢牛腸道感染的病原類型和感染類型可能會有所不同。據(jù)報(bào)道，BRV和BCoV主要感染3周齡以內(nèi)的犢牛[25]，隱孢子蟲主要感染15～21 d犢牛[26]，BNoV主要感染2月齡內(nèi)日齡較小犢牛[27]，與本研究統(tǒng)計(jì)結(jié)果基本 相符。

3.3 NCD主要病原分子特征

根據(jù)輪狀病毒VP6內(nèi)殼多肽的特異性，輪狀病毒被分為A、B、C、D、E、F和G等7個血清型，犢牛感染的輪狀病毒主要為A型[28]。VP6基因的遺傳進(jìn)化分析結(jié)果表明，本研究檢測到的輪狀病毒均屬于A型，大部分毒株與中國河北株(Bo/HB-BD/CH)聚在一支，提示本研究得到的大部分寧夏流行株與河北地區(qū)流行株同源。輪狀病毒可跨種感染，BRV可感染人，并與人輪狀病毒發(fā)生重組[29-30]，本研究中的Y7株與人糞便分離到的美國株(Hu/2011825363/USA)聚為一支，提示該地區(qū)存在可感染人的BRV流行株，因此，加強(qiáng)BRV的研究具有一定的公共衛(wèi)生學(xué)意義。

BCoV中編碼核衣殼蛋白的N蛋白基因高度保守，常作為BCoV檢測的靶基因[31]。N蛋白的遺傳進(jìn)化分析結(jié)果表明，本研究檢測到的核苷酸序列為BCoV N蛋白基因序列，與中國四川報(bào)道株(Bo/SC1/CH)和美國株(Bo/4-17-08/USA)遺傳距離較近[32]，提示寧夏地區(qū)BCoV流行株可能存在國內(nèi)和國際兩種傳播途徑，應(yīng)當(dāng)加強(qiáng)本地區(qū)BCoV流行趨勢的監(jiān)控。

諾如病毒可分為GⅠ、GⅡ、GⅢ、GⅣ、GⅤ和GⅥ基因型，已報(bào)到感染牛的諾如病毒為GⅢ型，該型不感染人[33]。GⅢ型BNoV又可分為GⅢ.1型和GⅢ.2型，臨床上以GⅢ.2型為主，但感染GⅢ.1型BNoV的病牛臨床癥狀更為明顯[25]。通過RdRp基因的遺傳進(jìn)化分析可對諾如病毒進(jìn)行分型鑒定[34]。BNoVRdRp基因的遺傳進(jìn)化分析結(jié)果表明，本研究中檢出的BNoV全部為GⅢ.2型，與GeneBank中收錄的BNoV的RdRp片段核苷酸一致性為87.6%～98.4%。王玥琳[13]在河南、四川等地區(qū)檢出的BNoV大部分為GⅢ.2型，小部分為GⅢ.1型，本研究中尚未檢出GⅢ.1型BNoV。由于BNoV為寧夏地區(qū)首次檢出，且檢出率較高，未檢出GⅢ.1型BNoV，提示應(yīng)該進(jìn)一步擴(kuò)大該地區(qū)BNoV流行株的研究，這將對NCD的防控提供參考依據(jù)。

隱孢子蟲是引起犢牛腹瀉的最主要寄生蟲，中國已報(bào)道的可引起犢牛腹瀉的隱孢子蟲包括瑞氏隱孢子蟲(C.ryanae)、牛隱孢子蟲、微小隱孢子蟲、安氏隱孢子蟲(C.andersoni)和人隱孢子蟲(C.hominis)等，其中牛隱孢子蟲是優(yōu)勢種，其次是微小隱孢子蟲[35]。高?；鄣萚15]報(bào)道寧夏地區(qū)感染奶犢牛的隱孢子蟲主要為牛隱孢子蟲、微小隱孢子蟲和瑞士隱孢子蟲，而本研究中檢測到的隱孢子蟲為牛隱孢子蟲和微小隱孢子蟲，其中微小隱孢子蟲檢出率較高。但高?；踇15]的研究中未對犢牛進(jìn)行嚴(yán)格區(qū)分，而斷奶后犢牛為混群飼喂，感染隱孢子蟲種類可能更多，所以本研究中檢測到的隱孢子蟲種類較少。據(jù)Li等[36]報(bào)道，微小隱孢子蟲的檢出率與犢牛日齡有關(guān)，斷奶前新生犢牛糞便中微小隱孢子蟲的檢出率高于其他種，與本研究檢測結(jié)果一致。但本研究得到的6株隱孢子蟲流行株聚類較為分散，表明該地區(qū)內(nèi)隱孢子蟲具有豐富的遺傳多樣性，提示其演化方向較廣，需要進(jìn)一步擴(kuò)大探究。

4 結(jié) 論

本研究對寧夏部分地區(qū)主要奶牛養(yǎng)殖園區(qū)斷奶前新生犢牛進(jìn)行腹瀉的主要相關(guān)病原檢測，對主要病原的分布和感染特點(diǎn)、遺傳進(jìn)化特征進(jìn)行分析。流行病學(xué)調(diào)查結(jié)果表明，腹瀉主要發(fā)生于30日齡內(nèi)犢牛，BRV、BCoV、BNoV和隱孢子蟲是造成該地區(qū)斷奶前犢牛腹瀉的主要病原，感染形式以單一病原感染為主，BRV和BNoV檢出率相對較高，其中BNoV在寧夏地區(qū)首次檢出。遺傳進(jìn)化分析結(jié)果表明，寧夏部分地區(qū)新生犢牛檢出的BRV、BCoV與國內(nèi)已報(bào)到的毒株基本一致，BNoV主要為GⅢ.2型，隱孢子蟲主要為微小隱孢子蟲，研究結(jié)果為寧夏地區(qū)新生犢牛腹瀉的防控提供重要參考。