劉 波 斯木吉德 王純潔 阿日查 張 劍 陳 浩 何利娜 張 晨 崔銀雪 敖日格樂*

(1.內(nèi)蒙古農(nóng)業(yè)大學(xué) 動(dòng)物科學(xué)學(xué)院，呼和浩特 010018；2.內(nèi)蒙古農(nóng)業(yè)大學(xué) 獸醫(yī)學(xué)院，呼和浩特 010018)

肉牛的體況評(píng)分是評(píng)價(jià)一頭肉牛個(gè)體能量存儲(chǔ)情況以及肉牛能量平衡程度的重要方法，是目前廣泛應(yīng)用的評(píng)估肉牛體能儲(chǔ)備是否平衡的實(shí)用方法，是肉牛飼養(yǎng)管理中不可缺少的工具[1]。體況對(duì)母牛的繁殖機(jī)能會(huì)產(chǎn)生重要影響，在整個(gè)繁殖周期中維持肉牛適宜的體況，既能保持其良好的繁殖性能，又能降低飼養(yǎng)成本，取得較好的經(jīng)濟(jì)效益[2]。周貴等[3]研究表明，為取得理想的繁殖性能，繁殖母牛的體況應(yīng)保持在2.5～3分(5分制)或5～6分(9分制)的范圍內(nèi)。

妊娠初期在整個(gè)妊娠期內(nèi)有重要作用，有研究表明妊娠初期陰道菌群發(fā)生紊亂會(huì)顯著增加早產(chǎn)和流產(chǎn)風(fēng)險(xiǎn)，妊娠初期的陰道菌群狀態(tài)可以預(yù)測(cè)妊娠結(jié)局[4-5]。姜宏星等[6]研究表明奶牛產(chǎn)前產(chǎn)后的體況直接影響本胎次乃至終生的繁殖水平，體況過肥或偏瘦均會(huì)造成奶牛繁殖力下降。周穎等[7]研究表明母豬在妊娠期保持中等體況繁殖性能較好。李元?jiǎng)偟萚8]研究發(fā)現(xiàn)在妊娠初期保持母貂中等體況有助于發(fā)揮其繁殖性能。因此在妊娠初期保持良好的體況是很重要的。牛的陰道中存在大量種類繁多的微生物及多種常駐菌。正常情況下，這些常駐菌群與宿主、環(huán)境之間相互協(xié)調(diào)、相互作用、相互制約，共同維持陰道的平衡狀態(tài)[9]。這種平衡狀態(tài)與牛的生殖健康緊密相關(guān)，一旦遭到破壞，如外界環(huán)境、自身激素水平和免疫力的改變，就會(huì)引發(fā)一些繁殖障礙疾病[10]。劉暢[9]研究發(fā)現(xiàn)健康、體況良好的肉牛陰道細(xì)菌有氣球菌、綠色氣球菌、睡眠嗜血桿菌、多動(dòng)物鏈球菌、嗜冷桿菌、大腸桿菌和鞘氨醇單胞菌，黃體期和卵泡期陰道菌群結(jié)構(gòu)有差異。Mateus等[16]研究指出，與未妊娠的奶牛相比，妊娠期的奶牛往往表現(xiàn)出較低的陰道菌群Alpha多樣性。Gunay等[11]研究發(fā)現(xiàn)狗卵泡期陰道細(xì)菌數(shù)量顯著高于黃體期。

然而到目前為止，關(guān)于牛微生物區(qū)系的研究主要集中在胃腸道，而對(duì)于不同體況放牧肉用繁殖母牛妊娠初期陰道菌群結(jié)構(gòu)的研究尚未見報(bào)道。對(duì)于肉牛養(yǎng)殖來說，保持繁殖母牛擁有較好的繁殖力，對(duì)擴(kuò)大種群以及提高養(yǎng)殖效益具有重要意義。而擁有較好繁殖能力的前提是母牛在妊娠期間保持良好的體況，如果過肥或者偏瘦則會(huì)引起難產(chǎn)、出生犢牛體質(zhì)瘦弱、成活率低等問題，同時(shí)良好的陰道菌群結(jié)構(gòu)可以降低繁殖障礙疾病的發(fā)生率。因此本研究在妊娠初期通過對(duì)不同體況肉牛陰道菌群結(jié)構(gòu)的研究，了解繁殖母牛陰道菌群組成以及是否有致病菌異常增加，對(duì)發(fā)展肉牛健康養(yǎng)殖、降低繁殖障礙疾病發(fā)病率和提高繁殖力具有重要意義。

1 材料與方法

1.1 試驗(yàn)動(dòng)物

在內(nèi)蒙古自治區(qū)錫林郭勒盟烏拉蓋草原選擇相同放牧條件下2～3 胎身體健康的妊娠初期西門塔爾母牛18 頭。采用9分制評(píng)分標(biāo)準(zhǔn)[12]將其分為3組，每組6 頭分別為FC(肥胖組),7～8分、MC(標(biāo)準(zhǔn)組)，5～6分、TC(偏瘦組)，3～4分。

1.2 樣品采集

在樣品采集之前，先用清水將牛外陰部清洗干凈，然后用25%碘溶液對(duì)外陰部進(jìn)行消毒，避免糞便污染。使用本實(shí)驗(yàn)室自制取樣器采集樣品。將棉拭子置于陰道2/3處1 min，使棉拭子與陰道分泌物充分接觸，將棉拭子取出后放在1.5 ml離心管中，迅速放于液氮中并送回實(shí)驗(yàn)室，以便細(xì)菌DNA提取。

1.3 指標(biāo)測(cè)定及數(shù)據(jù)分析

將粘有陰道分泌物的棉簽置于離心管中放置于干冰上送往北京諾和生物信息科技有限公司。對(duì)341F(5′-CCTAYGGGRBGCASCAG-3′)和806R(5′-GGACTACNNGGGTATCTAAT-3′)的16S rRNA V3～V4區(qū)域進(jìn)行擴(kuò)增測(cè)序。詳細(xì)步驟見參考文獻(xiàn)[13]。利用Barcode區(qū)分樣本序列，在對(duì)各組樣本序列進(jìn)行過濾。利用Uparse軟件(Uparse V7.0.1001)對(duì)所有有效數(shù)據(jù)進(jìn)行OTUs聚類和物種分析，所有聚類均在默認(rèn)為97%的相似度下進(jìn)行。使用R軟件(Version 2.15.3)繪制稀釋曲線、Rank abundance曲線，使用Qiime軟件(Version 1.9.1)進(jìn)行Alpha多樣性及菌群結(jié)構(gòu)分析，Alpha 多樣性指數(shù)包括豐富度指數(shù)(Chao1 指數(shù)、ACE指數(shù))、多樣性指數(shù)(Simpson指數(shù)、Shannon 指數(shù))和覆蓋率(Good’s coverage 指數(shù))。利用Beta多樣性分析評(píng)價(jià)樣品在物種復(fù)雜性方面的差異，對(duì)OTU數(shù)據(jù)進(jìn)行主坐標(biāo)分析(Principle coordination analysis，PCoA), 使用層次聚類方法非加權(quán)組平均法(Unweighted pair-group method with arithmetic means，UPGMA)構(gòu)建系統(tǒng)發(fā)育樹，對(duì)距離矩陣進(jìn)行解釋，結(jié)合細(xì)菌群落柱狀圖，研究不同組間或樣品間的菌群相似性和分類學(xué)差異性[14]。

數(shù)據(jù)采用Excel 2016整理，采用SPSS 20.0軟件進(jìn)行方差分析，并以P<0.05進(jìn)行顯著性檢驗(yàn)，結(jié)果均以平均值±標(biāo)準(zhǔn)差表示。

2 結(jié)果與分析

2.1 測(cè)序結(jié)果及α多樣性分析

通過Illumina HiSeq 測(cè)序平臺(tái)測(cè)序?qū)﹃幍兰?xì)菌16S rRNA基因V3-V4高變區(qū)進(jìn)行測(cè)序分析，18 份樣品平均每個(gè)樣品測(cè)得70 567條有效序列，經(jīng)過質(zhì)控平均得到 68 309 條序列，質(zhì)控有效率達(dá)96.84%。以97%的一致性將序列聚類成為OTUs，共得到2 677 個(gè)OTUs，對(duì)OTUs序列與Silva132數(shù)據(jù)庫(kù)進(jìn)行物種注釋。注釋結(jié)果中，共有37.58%的OTU注釋到門水平；32.84%的OTU注釋到科水平；24.36%的OTU注釋到屬水平。

通過繪制稀疏曲線來評(píng)判測(cè)序數(shù)據(jù)量是否足以反映樣品中的物種多樣性(圖1)，隨著測(cè)序深度的增加，在測(cè)序量達(dá)到32 000時(shí)，曲線趨向平坦，說明取樣的數(shù)量合理，測(cè)序深度已經(jīng)基本覆蓋到樣品中所有的物種，且能夠很好地反映出所有樣品的細(xì)菌群落結(jié)構(gòu)和多樣性。

表1顯示了不同體況陰道菌群Alpha多樣性指數(shù)，直接反映微生物群落的豐度及多樣性。Shannon值越大，生物多樣性越高；Simpson值越低，表明生物多樣性越高；Chao1指數(shù)和ACE指數(shù)是估計(jì)物種總數(shù)的常用指數(shù)，值越大，生物多樣性越高。Good’s coverage表示樣本的測(cè)序深度。標(biāo)準(zhǔn)組陰道分泌物微生物群落的Alpha多樣性指數(shù)的Shannon指數(shù)、Chao1指數(shù)和ACE指數(shù)均顯著低于肥胖組和偏瘦組(P<0.05)，Simpson指數(shù)差異不顯著(P>0.05)，表明標(biāo)準(zhǔn)組體況較好，機(jī)體內(nèi)分泌正常、陰道內(nèi)環(huán)境趨于穩(wěn)定，物種多樣性處于較低水平。各組陰道菌群的Good’s coverage指數(shù)均在99.4%以上，顯示樣本測(cè)序質(zhì)量較高，未被檢測(cè)到的可能性極低，足以檢測(cè)陰道分泌物樣本中的微生物群落信息。

不同顏色代表不同樣品Different samples were represented by different colors圖1 18份樣品的稀釋曲線Fig.1 Dilution curve of 18 samples

表1 不同體況妊娠初期母牛陰道菌群α多樣性指數(shù)Table 1 Alpha diversity index of vaginal flora of cows in early pregnancy with different body conditions

2.2 β多樣性分析

β多樣性分析是對(duì)不同樣本的微生物群落構(gòu)成進(jìn)行比較分析。主坐標(biāo)分析反映了各樣本中菌群構(gòu)成的相似性。基于Weighted Unifrac距離進(jìn)行PCoA分析(圖2)，并選取貢獻(xiàn)率最大的主坐標(biāo)組合進(jìn)行作圖展示。如果樣本距離越接近，表示物種組成結(jié)構(gòu)越相似，因此群落結(jié)構(gòu)相似度高的樣本傾向于聚集在一起，群落差異較大的樣本則會(huì)遠(yuǎn)遠(yuǎn)分開。由圖2可知，第一主成分(PC1)和第二主成分(PC2)的貢獻(xiàn)率分別為總變量的58.86%和18.89%，肥胖組(FC)、標(biāo)準(zhǔn)組(MC)、偏瘦組(TC)幾乎聚在一起，差異不顯著(P>0.05)。

橫縱坐標(biāo)分別為第一和第二主成分，百分比表示主成分對(duì)樣品差異的解釋率；圖中的每個(gè)點(diǎn)表示一個(gè)樣品，同一組的樣品用同一種顏色表示The abscissa and ordinate were the first and second principal components respectively, and the percentage represents the interpretation rate of the principal component to the sample difference; each point in the diagram represents a sample, and the samples of the same group were shown in the same color圖2 各組樣品基于Weighted Unifrac距離PCoA分析Fig.2 Weighted Unifrac distance PCoA analysis for each group of samples

2.3 陰道菌群構(gòu)成及豐度分布

2.3.1門水平陰道菌群分布豐度

由表2可知，不同體況妊娠初期繁殖母牛陰道分泌物中10大優(yōu)勢(shì)菌門分別是厚壁菌門(Firmicutes)、擬桿菌門(Bacteroidetes)、變形菌門(Proteobacteria)、無壁菌門(Tenericutes)、放線菌門(Actinobacteria)、梭桿菌門(Fusobacteria)、疣微菌門(Verrucomicrobia)、未知菌(Unidentified_Bacteria)、螺旋體菌門(Spirochaetes)和Melainabacteria。厚壁菌門在肥胖組、標(biāo)準(zhǔn)組、偏瘦組的豐度分別為54.16%、56.38%和57.03%，擬桿菌門在3組的豐度分別為24.56%、30.04%和22.37%，變形菌門在3組的豐度分別為3.86%、3.14%和5.83%。而在非優(yōu)勢(shì)菌門中，無壁菌門在肥胖組、標(biāo)準(zhǔn)組、偏瘦組的豐度分別為1.34%、2.06%和1.72%；放線菌門在3組的豐度分別為1.63%、1.67%和1.79%；梭桿菌門在3組的豐度分別為1.80%、1.25%和1.23%；疣微菌門在3組的豐度分別為0.58%、0.35%和0.70%；未知菌在3組的豐度分別為0.50%、0.36%和0.89%；螺旋體菌門在3組的豐度分別為0.31%、0.25%和0.22%； Melainabacteria在3組的豐度分別為0.26%、0.15%和0.20%。結(jié)果表明，陰道菌群隨著體況不同而發(fā)生變化，偏瘦組變形菌門分布豐度高于標(biāo)準(zhǔn)組和肥胖組，擬桿菌門的分布豐度標(biāo)準(zhǔn)組高于肥胖組和偏瘦組，但差異不顯著(P>0.05)；其他菌門分布豐度無明顯變化，差異不顯著(P>0.05)。

表2 門水平上陰道菌群豐度Table 2 Abundance of vaginal microbiota at phylum level %

2.3.2屬水平陰道菌群分布豐度

在屬水平上，由于屬的種類多且部分屬的相對(duì)豐度較低，因此僅對(duì)排名前10的屬進(jìn)行分析。由表3 可知，不同體況妊娠初期繁殖母牛陰道分泌物中10大優(yōu)勢(shì)菌屬分別是擬桿菌屬(Bacteroides)、阿利斯佩斯菌屬(Alistipes)、氣球菌屬(Aerococcus)、費(fèi)氏菌屬(Faecalibacterium)、鏈桿菌屬(Streptobacillus)、鏈球菌屬(Streptococcus)、棒桿菌屬(Corynebacterium)、嗜組織菌屬(Histophilus)、脲原體菌屬(Ureaplasma)和支原體屬(Mycoplasma)。擬桿菌屬在肥胖組、標(biāo)準(zhǔn)組、偏瘦組的豐度分別為3.79%、4.26%和3.38%；阿利斯佩斯菌屬在3組的豐度分別為2.78%、2.71%和2.15%；氣球菌屬在3組的豐度分別為1.15%、0.85%和1.12%；費(fèi)氏菌屬在3組的豐度分別為0.73%、0.73%和0.33%；鏈桿菌屬在3組的豐度分別為0.63%、0.80%和1.13%；鏈球菌屬在3組的豐度分別為0.58%、0.17%和0.48%；棒桿菌屬在3組的豐度分別為0.51%、0.48%和0.29%；嗜組織菌屬在3組的豐度分別為0.03%、0.99%和0.27%；脲原體屬在3組的豐度分別為0.17%、0.20%和0.14%；支原體屬在3組的豐度分別為0.04%、0.09%和0.01%。結(jié)果表明，隨著體況的不同，標(biāo)準(zhǔn)組嗜組織菌屬分布豐度顯著高于其他2組(P<0.05)、擬桿菌屬分布豐度高于其他2組，氣球菌屬分布豐度低于其他2組，偏瘦組鏈桿菌屬分布豐度高于其他2組，但差異均不顯著(P>0.05)。

表3 屬水平上陰道菌群豐度Table 3 Abundance of vaginal microbiota at genus level %

3 討 論

本研究利用16S rRNA高通量測(cè)序技術(shù)，對(duì)不同體況妊娠初期繁殖母牛的陰道菌群的豐度和多樣性進(jìn)行了比較。結(jié)果表明，體況的不同對(duì)母牛陰道微生物區(qū)系結(jié)構(gòu)有一定影響，菌群豐度各不相同。通過Alpha多樣性指數(shù)分析，標(biāo)準(zhǔn)組母牛陰道菌群多樣性低于肥胖組和偏瘦組。從PCoA分析結(jié)果來看，標(biāo)準(zhǔn)組和肥胖組菌群相似性較高，但和偏瘦組差異不大。Wang等[15]研究表明，患有子宮內(nèi)膜炎奶牛陰道Alpha多樣性指數(shù)顯著高于健康奶牛，這與子宮內(nèi)膜炎的發(fā)生過程一致，即由優(yōu)勢(shì)菌支配的微生物群落結(jié)構(gòu)被破壞，病原菌數(shù)量增加，子宮內(nèi)膜炎發(fā)病的一個(gè)重要因素是陰道內(nèi)的細(xì)菌擴(kuò)增進(jìn)入子宮。而標(biāo)準(zhǔn)組體況母牛機(jī)體狀況良好，機(jī)體為妊娠提供了較好的環(huán)境保障，陰道內(nèi)環(huán)境相對(duì)其他兩組較為穩(wěn)定，因此菌群結(jié)構(gòu)維持在較低水平。通過對(duì)門水平上菌群結(jié)構(gòu)分析發(fā)現(xiàn)，標(biāo)準(zhǔn)組陰道主要菌門為變形菌門、厚壁菌門、擬桿菌門。Laguardia-Nascimento等[16]研究指出不同品種妊娠和非妊娠牛陰道菌群主要由厚壁菌門(40%～50%)、擬桿菌門(15%～25%)、變形菌門(5%～25%)和未知菌門(10%～20%)組成。這與本研究基本吻合，擬桿菌門和無壁菌門有助于子宮內(nèi)膜的恢復(fù)[17]。牛易受到各種病原體引起的與生殖道有關(guān)的疾病的影響，如類桿菌屬、支原體、嗜組桿菌屬、梭桿菌屬、鏈球菌屬等[18,19-20]。球菌屬通常定植在哺乳動(dòng)物的皮膚和粘膜上，也存在于口腔、腸道、泌尿生殖道和上呼吸道[21]。支原體屬和脲原體屬?gòu)V泛定植于動(dòng)物泌尿生殖道黏膜上，并不會(huì)引起任何癥狀，而當(dāng)機(jī)體免疫力低下時(shí)會(huì)引發(fā)炎癥[22]。Funke等[23]研究認(rèn)為棒狀桿菌是人類正常皮膚菌群的一部分。傳統(tǒng)培養(yǎng)方法研究發(fā)現(xiàn)，健康牛陰道菌群主要由鏈球菌屬、葡萄球菌屬、腸桿菌屬、乳桿菌屬及少量真菌組成[24-25]。侯榮倫[26]研究健康奶牛生殖道菌群發(fā)現(xiàn)乳酸菌是主要優(yōu)勢(shì)菌，且不孕癥奶牛生殖道內(nèi)乳酸菌數(shù)量低于健康奶牛。在本試驗(yàn)中，雖然檢測(cè)出乳酸桿菌屬，但其含量在3個(gè)組中均屬于偏低水平，未在圖表中展示。這與劉暢[9]的研究結(jié)果一致，肉牛陰道中存在少量的乳酸桿菌屬。標(biāo)準(zhǔn)組嗜組織菌豐度顯著高于肥胖組和偏瘦組，其他菌屬差異不大。雖然有研究表明該菌可能是引起奶牛子宮內(nèi)膜炎的因素之一且病因不明確[27]，但本試驗(yàn)所選母牛均為未患繁殖障礙疾病牛，只在體況上進(jìn)行區(qū)分，因此可能是其他因素導(dǎo)致該菌豐度增加，這也說明了陰道中存在少量致病菌并不足以致病，但可以使養(yǎng)殖者及時(shí)了解母牛陰道菌群中有害菌是否有升高的趨勢(shì)。因?yàn)榕５耐怅幉课挥诟亻T的正下方，經(jīng)常會(huì)被排出的糞便所覆蓋，而外陰部的開口可能會(huì)使糞便中的少量致病菌進(jìn)入，有研究證明，許多生殖道疾病均是由糞便中微生物引起[28-29]。本試驗(yàn)中檢測(cè)出的少量費(fèi)氏桿菌屬和阿利斯佩斯菌屬也有可能是在樣品采集之前已經(jīng)通過糞便進(jìn)入陰道[30-31]。由于本試驗(yàn)結(jié)果中肥胖組、標(biāo)準(zhǔn)組和偏瘦組陰道菌群的主要優(yōu)勢(shì)菌門和優(yōu)勢(shì)菌屬組間差異并不顯著，可能是因?yàn)楸驹囼?yàn)所選母牛均為體質(zhì)健康并處于妊娠初期階段，而體況對(duì)母牛在妊娠初期階段陰道內(nèi)環(huán)境穩(wěn)態(tài)影響不大進(jìn)而導(dǎo)致主要優(yōu)勢(shì)菌門和優(yōu)勢(shì)菌屬豐度處于正常水平差異不顯著，但在妊娠初期保持母牛陰道內(nèi)環(huán)境穩(wěn)定仍是十分重要的，可在之后的試驗(yàn)中對(duì)妊娠中后期不同體況母牛陰道菌群進(jìn)行測(cè)定，進(jìn)一步探索體況對(duì)陰道菌群的影響。

研究表明，陰道微生物群落通過生物屏障、產(chǎn)生乳酸、細(xì)菌素和過氧化氫等多種機(jī)制在維持生殖道健康方面起著至關(guān)重要的作用[32-33]。機(jī)體會(huì)受到多重因素的影響從而導(dǎo)致陰道菌群結(jié)構(gòu)發(fā)生變化，這些因素包括基因、自身的防御屏障系統(tǒng)、年齡、品種、免疫功能、激素分泌等[7]。體況過肥或過瘦會(huì)影響血液離子、激素、營(yíng)養(yǎng)物質(zhì)代謝水平，降低機(jī)體免疫力和繁殖性能[34]，從而對(duì)妊娠初期母牛陰道微生物穩(wěn)態(tài)產(chǎn)生影響。而陰道微生物菌群失調(diào)可能是引發(fā)陰道炎和生殖道感染的重要因素，從而損害生殖性能，造成重大的經(jīng)濟(jì)損失[35]。另外致病性細(xì)菌感染可引發(fā)子宮內(nèi)膜炎最終導(dǎo)致胚胎早期死亡，嚴(yán)重阻礙繁殖母牛繁殖力。盡管抗生素已經(jīng)被廣泛應(yīng)用于雌性動(dòng)物繁殖疾病，但細(xì)菌耐藥性和抗生素殘留等問題正變得越來越突出[36]。Kavyashree等[37]研究發(fā)現(xiàn)對(duì)陰道菌群鑒定篩查可以有效調(diào)控奶牛生殖健康。因此，可通過對(duì)體況過肥和偏瘦的妊娠母牛分群管理，對(duì)于體況過肥的母?？蛇m當(dāng)縮短其放牧采食時(shí)間，減少其能量攝入，逐步使體況恢復(fù)到標(biāo)準(zhǔn)體況范圍，而對(duì)于體況偏瘦的母牛則延長(zhǎng)其放牧?xí)r間或者收牧后進(jìn)行適量補(bǔ)飼，進(jìn)而增重到標(biāo)準(zhǔn)體況，進(jìn)而使繁殖母牛在妊娠初期避免過肥或過瘦、保持良好的體況，同時(shí)明確微生物區(qū)系在生殖道健康中的作用和了解妊娠初期菌群結(jié)構(gòu)對(duì)防治生殖道感染具有重要意義，可通過提早對(duì)不良體況繁殖母牛進(jìn)行陰道微生物檢查，盡早發(fā)現(xiàn)可能引起繁殖疾病的菌群變化，及早預(yù)防，降低繁殖障礙疾病發(fā)生的概率，減少損失。

4 結(jié) 論

利用基于微生物16S rRNA基因測(cè)序技術(shù)，探索了不同體況妊娠初期繁殖母牛陰道菌群變化情況。本研究發(fā)現(xiàn)，標(biāo)準(zhǔn)體況母牛陰道微生物多樣性和豐度最低，厚壁菌門、擬桿菌門、變形菌門分別為肥胖組、標(biāo)準(zhǔn)組和偏瘦組主要的優(yōu)勢(shì)菌門，擬桿菌屬均為肥胖組、標(biāo)準(zhǔn)組和偏瘦組主要的優(yōu)勢(shì)菌屬，不同體況母牛陰道微生物群落組成與結(jié)構(gòu)均有變化并存在少量致病菌，但是差異不顯著。本研究為使母牛在妊娠初期保持良好體況并有利于繁殖性能的發(fā)揮，同時(shí)為養(yǎng)殖者及時(shí)了解牛陰道菌群隨體況變化提供了重要基礎(chǔ)數(shù)據(jù)。