ORIGINAL RESEARCH article
Sec. Systems Microbiology
Volume 9 - 2018 | https://doi.org/10.3389/fmicb.2018.03217
Composition of Rumen Bacterial Community in Dairy Cows With Different Levels of Somatic Cell Counts
- MoE Key Laboratory of Molecular Animal Nutrition, College of Animal Sciences, Institute of Diary Science, Zhejiang University, Hangzhou, China
Mastitis is an inflammatory disease, affects the dairy industry and has a severe economic impact. During subclinical mastitis, milk production and milk quality deteriorates. Recently, rumen microbial composition has been linked to rumen health, but few studies have investigated the effect of rumen microbiota on mammary health in cows. This study was undertaken to identify the rumen microbial composition and associated microbial fermentation in cows with different somatic cell counts (SCC), with the speculation that cows with different health statuses of the mammary gland have different rumen bacterial composition and diversity. A total of 319 Holstein dairy cows fed the same diet and under the same management were selected and divided into four groups as SCC1 (N = 175), SCC2 (N = 49), SCC3 (N = 49), and SCC4 (N = 46) with < 200,000, 200,001–500,000, 500,001–1,000,000, and >1,000,000 somatic cells/mL, respectively. Further, 20 cows with the lowest SCC and 20 cows with the highest SCC were identified. The rumen microbial composition was profiled using 16S rRNA sequencing, along with measurement of rumen fermentation parameters and milking performance. Compared to low SCC, cows with high SCC showed poorer milk yield, milk composition, and rumen volatile fatty acids concentration, but higher rumen bacterial diversity. Although the predominant rumen bacterial taxa did not vary among the SCC groups, the relative abundance of phyla SR1 and Actinobacteria, unclassified family Clostridiales and genus Butyrivibrio were significantly different. In addition, Proteobacteria and family Succinivibrionaceae were enriched in cows with low SCC. Our results suggest that specific rumen microbes are altered in cows with high SCC.
Mastitis of dairy cows is an inflammatory disease caused by intramammary infection (IMI), with high incidence and prevalence. It involves a host-pathogen interaction driven by host and bacterial determinants. Subclinical mastitis is characterized by no visible changes in the appearance of milk or the udder, but milk production decreases, milk composition is altered, somatic cell count increases (SCC), and pathogens may be present in the secretion (Batavani et al., 2007). The IMI in dairy cows can cause severe economic losses due to reduced milk production, discarded milk, lower conception rates, premature culling, and treatment costs (Seegers et al., 2003; Halasa et al., 2007). In the dairy industry, SCC have been extensively used to distinguish healthy quarters from quarters with an inflammatory response due to IMI (Tsenkova et al., 2001).
Recently, the fundamental role of naturally occurring diverse gastrointestinal microbes in the well-being of ruminants and humans has been demonstrated (Clemente et al., 2012; Mao et al., 2015). The constituents of microbiota have been shown to interact with each other and with the host immune system in ways that influence the development of disease. In human studies, it has been reported that milk bacteria as well as mastitis can be affected by gut microbiota (Rodriguez, 2014). The oral administration of strains of probiotics isolated from milk has been shown to be an alternative method for the treatment of human mastitis (Arroyo et al., 2010). These studies lead to the speculation that there could be an interaction between gut microbiota and mammary gland, although the mechanism remains under-investigated.
The rumen is considered as one of the major organs in dairy cows, which directly affect milk production and health of the host. Dairy cows depend on rumen microbes to convert human-indigestible feedstock to human-edible products, by providing volatile fatty acids (VFAs) and microbial proteins. In recent years, major efforts in characterizing the rumen microbial composition and function have been made, along with studying the factors that affect rumen microbiota, such as diet (Zhu et al., 2017), lactation stage, rumen development and age (in pre-weaned calves) (Robert et al., 2012). In addition, host-microbial interactions in the rumen have been studied, aiming to maximize the performance of dairy cows, including lactation performance (Xue et al., 2018), feed efficiency, and regulation of methane emission (McCann et al., 2014). Rumen microbial composition has also been linked to the health of rumen (acidosis and subacute acidosis) (Mao et al., 2015). To our knowledge, it is unknown whether the rumen microbiota could be associated with mammary health of the host.
In the current study, it is hypothesized that rumen microbial composition, diversity, and rumen fermentation parameters differed in cows with different health statuses of mammary gland (different levels of SCC). The16S rRNA gene sequencing was performed to investigate the rumen bacterial composition, and microbial fermentation products were characterized in 319 dairy cows with different levels of SCC. The current study provides a better understanding of rumen microbial composition and fermentation in cows with different health statuses of mammary gland.
Materials and Methods
The experimental procedures used in this study were approved by the Animal Care Committee of Zhejiang University (Hangzhou, China) and were in accordance with the University's guidelines for animal research.
A total of 319 multiparous (parity ranged from 2 to 7) Holstein dairy cows in mid-lactation (day in milk [DIM] ranged from 92 to 218) with no clinical signs of disease were selected from 652 lactating cows in a commercial dairy farm (Hangzhou, China). The animals were fed the same diet three times daily at 0630, 1400 and 2000 h, with forage to concentrate ratio of 45:55 (Supplementary Table S1), and had free access to water. Animals were divided into four groups according to the SCC in milk, including 175 cows with SCC ranging from 0 to 200,000 cells/mL (SCC1), 49 cows with SCC ranging from 200,001 to 500,000 cells/mL (SCC2), 49 cows with SCC ranging from 500,001 to 1,000,000 cells/mL (SCC3), and 46 cows with SCC greater than 1,000,000 cells/mL (SCC4). The SCC cut-offs were established based on a previous study on SCC values in healthy quarters (Hand et al., 2012). In order to increase the statistical power when comparing the rumen microbiota and fermentation (Supplementary Table S2), cows with extremely low and high SCC (including 20 cows with the lowest SCC [L_SCC] and 20 cows with the highest SCC [H_SCC]) were selected from the 319 cows based on the extreme phenotype characterization approach (Shabat et al., 2016). A previous study has reported that rumen samples from 16 cows are sufficient to cover the microbial diversity in rumen fed the same diet (Jami and Mizrahi, 2012).
Sample Collection and Analysis
Milk yield of all the 319 cows was recorded for three consecutive days. Milk samples of all the 319 cows were collected on the third day (sampling day) for the measurement of milk protein, fat, lactose, urea nitrogen (MUN) and SCC by infrared analysis (Laporte and Paquin, 1999) using a spectrophotometer (Foss-400; Foss Electric A/S, Hillerod, Denmark).
Rumen content of all the 319 cows was collected using oral stomach tubes (Shen et al., 2012) before the morning feeding on the sampling day. The pH of the rumen fluid was measured immediately using a pH meter (FE20-FiveEasy Plus™; Mettler Toledo Instruments Co. Ltd., Shanghai, China). The rumen content of all the 319 cows were used for analysis of VFAs and ammonia nitrogen (NH3-N). One set of triplicate samples of 1 mL were acidified with 20 μL orthophosphoric acid for analysis of VFAs, one set of triplicate samples of 1 mL were acidified with 6 N HCl and frozen at −20°C for analysis of NH3-N, as described by Hu et al. (2005). The third set of triplicate samples of 1 mL were frozen at −80°C until DNA extraction.
DNA Extraction and Sequencing
Total DNA was extracted from rumen content of all the 319 cows using the bead-beating method (Li et al., 2009). DNA qualities and quantities were measured using the NanoDrop 2000 Spectrophotometer (NanoDrop Technologies, Wilmington, DE, USA). The DNA was amplified using the 341F/806R primer set (341 F: 5′-CCTATYGGGRBGCASCAG-3′, 806R: 5′-GGACTACNNGGGTATCTAAT-3′), which targets the V3-V4 region of the bacterial 16S rRNA gene. The reaction contained 0.5 U of Taq polymerase (TransGen Biotech, Beijing, China) in 25 μL of 10 × PCR buffer, 200 μM of each dNTP, 0.2 μM of each primer and 2 μL of DNA (50 ng/μL). Double distilled water was added to make the volume reach 25 μL. The PCR reactions were performed using Phusion High-Fidelity PCR MasterMix (New England Biolabs LTD., Beijing, China) with the following program, 94°C for 3 min, 35 cycles of 94°C for 45 s, 50°C for 60s and 72°C for 90s, followed by 72°C for 10 min. The PCR products were visualized on 2% agarose gels and purified using the QIAquick Gel Extraction Kit (QIAGEN, Dusseldorf, Germany). Amplicon sequencing was conducted on an Illumina HiSeq platform using the paired-end 2 × 250 bp sequencing (Caporaso et al., 2012).
Analysis of Sequencing Data
Paired-end reads were merged using FLASH (V1.2.7, http://ccb.jhu.edu/software/FLASH/) (Magoč and Salzberg, 2011). Sequences were de-multiplexed and quality-filtered using QIIME (version 2, http://qiime.org/index.html) and bases with quality scores higher than 20 were retained for further analysis (Caporaso et al., 2010). Chimeric sequences were identified and removed using the UCHIME algorithm (UCHIME Algorithm, http://www.drive5.com/usearch/manual/uchime_algo.html) (Edgar et al., 2011). Operational taxonomic units (OTUs) were clustered based on 97% similarity threshold using UPARSE (Uparse v7.0.1001, http://drive5.com/uparse/) (Edgar, 2013), and taxonomy was assigned using the latest Greengenes database (http://greengenes.secondgenome.com Greengenes May 2013 release). The OTU-level alpha diversity of the bacterial communities was determined using various diversity indices (Chao1, Shannon, Simpson and Ace indexes) and calculated using procedures within QIIME (version 2). Jack-knifed beta diversity was calculated based on OTU-level weighted and unweighted Unifrac distances using QIIME (version 2), and visualized by principal coordinates analysis (PCoA) using PAST software (version 3.18, http://folk.uio.no/ohammer/past).
Lactation performance (DIM, parity, milk yield, milk fat, milk protein, lactose, MUN and SCC), total and individual rumen VFA (acetate, propionate, butyrate, isobutyrate, valerate, and isovalerate) concentrations and alpha diversity indices of the four SCC groups were analyzed using one-way ANOVA; and the relative abundances of bacterial taxa among the four SCC groups were compared using Kruskal-Wallis H test (four SCC groups). Scheffe's method was used for multiple comparisons of means among the four SCC groups. Lactation performance, total and individual rumen VFAs concentrations and alpha diversity indices in H_SCC and L_SCC groups were analyzed using t-test; and the relative abundances of bacterial taxa between the two groups were compared using Wilcoxon rank-sum test. Statistical analysis was performed using SPSS (version 22), with statistical significances declared at P < 0.05. The P-values from the Kruskal-Wallis H test and Wilcoxon rank-sum test were adjusted by the false discovery rate (FDR < 0.05) (Benjamini and Hochberg, 1995), with statistical significances declared at FDR-adjusted P < 0.05. Linear discriminant analysis effect size (LEfSe) was used to further compare relative abundances of microbial taxa in H_SCC and L_SCC groups, and significant differences were considered by a linear discriminant analysis (LDA) score > 2 and P < 0.05.
Measurement of SCC, Lactation Performance and Rumen VFAs
The parities of the 319 dairy cows were 3.1 ± 0.07 (mean ± standard error of the mean [SEM]), and the DIMs of the 319 dairy cows were 160 ± 1.87 (Table 1). Milk yield in SCC1 was the highest among the four groups, followed by SCC2 and SCC3, and the lowest in SCC4 (P < 0.01). Milk fat, milk protein and MUN content were significantly higher (P < 0.01), and the lactose content was significantly lower (P < 0.01) in SCC4 compared to the other 3 groups (Table 1). Rumen pH, and concentration of VFAs and NH3-N were not different among the four SCC groups. With the increase in SCC, the acetate to propionate (A: P) ratio increased, with higher A: P ratio in SCC4 (P < 0.05) than in SCC1 (Table 1).
Table 1. The somatic cell counts (SCC), milk yield, milk composition, and rumen fermentation parameters of the lactating dairy cows.
Milk yield (P < 0.01), lactose content (P < 0.01) and MUN (P = 0.03) were significantly higher in L_SCC, whereas the content of milk fat (P = 0.01) and milk protein (P < 0.01) were significantly higher in H_SCC. Molar proportion of propionate (P = 0.03) and valerate (P = 0.04) were significantly higher in L_SCC group than in H_SCC, whereas the A: P ratio was significantly higher in H_SCC (P = 0.02, Table 1).
Rumen Bacterial Communities in Cows With Different Levels of SCC
The amplicon sequencing of rumen samples generated a total of 19,253,662 high-quality sequences across all samples, and an average of 1,826 OTUs per sample, with a Good's coverage of 99.9% across all samples. A total of 26 bacterial phyla were identified across all the samples, with Firmicutes (49.8 ± 0.45%) and Bacteroidetes (33.6 ± 0.48%) being the most abundant and Proteobacteria (9.50 ± 0.38%), Tenericutes (2.16 ± 0.04%) and Spirochaetes (1.42 ± 0.04%) being less abundant (Table 2). At genus level, a total of 415 genera were identified. The bacterial genera with relative abundance > 1% were considered as predominant, accounting for 85.9 ± 0.13% of the total sequences. Among the predominant genera, Prevotella (19.3 ± 0.48%), unclassified_o_Clostridiales (9.16 ± 0.11%), Ruminococcus (8.37 ± 0.10%), unclassified_f_Succinivibrionaceae (8.03 ± 0.38%), and unclassified_f_Ruminococcaceae (7.95 ± 0.12%) were the most abundant.
Table 2. Phylum and genus composition (Relative abundances >1%) of the rumen bacteria in cows with different somatic cell counts (SCC).
The comparison of alpha diversity indices showed a significant difference of Shannon index among the four SCC groups at OTU level (Table 3). The Shannon index was significantly higher in SCC4 (P = 0.04) compared to SCC1. The PCoA based on weighted (Figure 1A) and unweighted Unifrac distances (Supplementary Figure S1A) showed no clear separation among the four SCC groups or between SCC1 and SCC4 (Figure 1B and Supplementary Figure S1B). Specific bacterial genera were unique to the rumen of cows with different levels of SCC (Supplementary Table S3), and included 29, 6, 10 and 8 unique genera in SCC1, SCC2, SCC3 and SCC4 groups, respectively. However, these unique bacterial genera in the different SCC groups were in low abundance, with relative abundance < 0.0002%.
Figure 1. Principal coordinates analysis (PCoA) based on weighted Unifrac distances of OTUs. Samples are indicated by points and colored for different SCC groups. (A) PCoA of the four SCC groups. (B) PCoA of the four SCC groups with SCC1 and SCC4 highlighted. (C) PCoA of L_SCC and H_ SCC groups.
The Shannon index was significantly higher (P = 0.03), and the Simpson index significantly lower in H_SCC (P = 0.04) compared to the L_SCC group (Table 3). The PCoA based weighted (Figure 1C) and unweighted Unifrac distances (Supplementary Figure S1C) showed no clear separation between L_SCC and H_SCC groups. In addition, 64 and 31 unique bacterial genera existed in L_SCC and H_SCC groups, respectively (Supplementary Table S4).
Differential Rumen Bacterial Taxa in Cows With Different Levels of SCC
Further, Kruskal-Wallis H analysis of the relative abundances of phyla revealed that phyla SR1 (FDR-adjusted P = 0.05) and Actinobacteria (FDR-adjusted P < 0.01) were significantly different among the four SCC groups (Figure 2A). At genus level, Butyrivibrio was significantly higher in SCC4 (FDR-adjusted P = 0.02, Figure 2B). Three unclassified taxa belonging to the families Clostridiales (FDR-adjusted P = 0.04), S24-7 (FDR-adjusted P = 0.04) and RF 16 (FDR-adjusted P = 0.03) were significantly different among the four SCC groups.
Figure 2. Comparison of relative abundances of bacterial phyla and genera in the different SCC groups. (A) The 10 most abundant phyla in the four groups. (B) The 20 most abundant genera are presented. *0.01 < P ≤ 0.05, **0.001 < P ≤ 0.01.
The most differentially abundant bacterial taxa in the H_SCC group tested by the LEfSe analysis belonged to the genera unclassified_f_RF16, Paenibacillus, unclassified_c_Deltaproteobacteria, unclassified_f _BS11, unclassified_f _Succinivibrionaceae, Lysinibacillus, unclassified_f _Pirellulaceae, unclassified_o _GMD14H09 and Sutterella (Figure 3). The genera unclassified_f _Succinivibrionaceae, Rhodobacter, Comamonas, Enterococcus, and unclassified_c _Gammaproteobacteria were more abundant in the L_SCC group (Figure 3). The genera unclassified_f_RF16 and unclassified_f_Succinivibrionaceae weighted most to the differences between the communities, with an absolute LDA score factor of approximately 3.
Figure 3. The LDA effect size (LEfSe) analysis of bacterial taxa between L_SCC and H_SCC groups. (A) Cladogram shows significantly enriched bacterial taxa (from phylum to genus level). (B) Bar chart shows LDA score of L_SCC and H_SCC groups. Significant differences are defined as P < 0.05 and LDA score > 2.0.
Animal IMI is the major factor that influences SCC (González-Rodríguez et al., 1995), with high SCC (> 200,000 cells/mL) in dairy cows indicating infection and low milk quality (Pyorala, 2003). The SCC at the cow level is commonly used as a tool to diagnose subclinical mastitis in dairy farms. Although the variation of SCC may exist among mammary quarters within a cow, the SCC at cow level can also reflect the general statuses of mammary quarters (Schukken et al., 2003) and is easier to detect relatively compared with SCC at quarter level. In this study, none of the cows in the cohort showed clinical signs of disease, but the cohort had varied SCC (ranging from 8,000 to 11,963,000 cells/mL), indicating that higher SCC cows (> 200,000 cells/mL, including cows in SCC2, SCC3, and SCC4) are at risk of IMI and suffer from subclinical mastitis.
It has been reported that milk yield and composition were negatively affected by IMI, and the degree of changes depended on the inflammatory response (Pyorala, 2003). This change can be the result of damage to mammary epithelial cells by either mastitis-causing bacteria or by the cow's immune response to infection (Gray et al., 2005). It has also been reported that the ruminal VFAs derived from microbial fermentation are also key factors that directly affect the milk biosynthesis (Flint and Bayer, 2008). The VFA concentrations were significantly different between cows with the lowest SCC (L_SCC group) and highest SCC (H_SCC group), with higher concentrations of propionate and valerate, and a lower A: P ratio in the cows of L_SCC group. Such differences in rumen microbial metabolites suggests an altered trend in rumen microbial fermentation in animals at risk of IMI. As the major microbial fermentation product, propionate is the main glucogenic nutrient in ruminants (Danfær, 1994; Reynolds et al., 2003), and glucose is the main precursor for lactose synthesis. The increase in post-hepatic glucose availability could be a potential regulator of milk yield. In this study, the higher rumen propionate, lactose, and better milk yield in cows with lower SCC compared to cows with higher SCC, suggests that IMI cows have a different rumen microbial fermentation pattern from healthy cows, and the fermentation products affect milk synthesis in the mammary gland. These observations could be associated with differential microbial composition and function in the rumen.
When the alpha diversity indices of rumen microbial community were compared, the Shannon index was significantly higher in H_SCC group, and the Simpson index were significantly lower in the L_SCC group. The differences in the alpha diversity indices between the two groups suggested that the rumen microbial population had higher diversity and less strength of dominance in cows with IMI (Rocchini et al., 2015). Additionally, the comparison of rumen bacterial taxa indicated that the predominant phyla (Firmicutes and Bacteroidetes) along with the well-studied genera were not different in cows with different SCC levels. For example, the main cellulolytic and hemicellulolytic bacteria (Fibrobacter, Ruminococcus, Clostridium and Eubacterium), amylolytic bacteria (Streptococcus and Ruminobacter), proteolytic bacteria (Treponema and Lachnospira) and saccharolytic bacteria (Succinivibrio, Selenomonas, Lactobacillus, and Bifidobacterium) did not vary among cows with different SCC levels. These results suggest that the bacterial communities, which play a fundamental role in rumen functions, are relatively stable in cows with different health statuses of mammary gland.
Specific microbes, notably at phylum level were observed to be variable among cows with different levels of SCC (four groups). The SCC1 group had lower abundance of SR1 and higher abundance of Actinobacteria than the other groups. The phylum SR1 exists in the rumen of many herbivores (Davis et al., 2009), and Actinobacteria are regular, though infrequent, members of the rumen microflora, representing up to 3% of total rumen bacteria (Singh et al., 2012). However, information is limited on the ecology and biology of these two taxa (Sulak et al., 2012; Campbell et al., 2013). Further studies to identify the functions of taxa within SR1 and Actinobacteria in the rumen are required to determine whether these two taxa are linked to IMI in dairy cows. Additionally, when the H_SCC and L_SCC groups were compared, Proteobacteria was significantly enriched in the L_SCC group, and unclassified Succinivibrionaceae, the most abundant genera within Proteobacteria, was the major contributor to this difference. Microbes belonging to Succinivibrionaceae commonly produce succinate, the precursor of propionate (Pope et al., 2011). Although succinate was not measured in this study, propionate was significantly higher in the L_SCC group, suggesting that this taxon could contribute to the higher propionate in the L_SCC group.
Traditionally, the mammary gland is considered a sterile environment, and bacterial cells in milk are considered to be the result of environmental contamination. However, the detection of live bacterial cells and/or DNA from anaerobic species that are usually related to the gut environment in milk has fueled scientific debate on the origin of milk-associated bacteria (Jost et al., 2014). Young et al. (2015) investigated the microbial composition of milk, blood, and feces of healthy lactating cows, and found that OTUs belonging to Ruminococcus, Bifidobacterium, and Peptostreptococcaceae were present in all three samples from the same animal, suggesting the existence of endogenous entero-mammary pathways during lactation in dairy cows. In a recent study comparing structural and functional features of the gut microbiome between mastitic and healthy cows, Ma et al. (2016) reported that the changes of fecal microbial community of mastitis cows was similar to that of the milk, characterized by a general increase in the mastitis pathogens (Enterococcus, Streptococcus, and Staphylococcus) and deprivation of Lactobacillus and its members. The above studies indicated that the gut microbiota could be associated with different health statuses of mammary gland. Since the gastrointestinal cross-talk in ruminants (including different segments, such as rumen, small and large intestine), along with the relationship between gastrointestinal microbiota and health statuses of mammary gland in ruminants has not been studied in detail, it is therefore recommended that both gut and rumen microbiota are taken into consideration and to investigate the relationship between mammary gland health and gastrointestinal microbiota in further studies.
Although the cows with the lowest and highest SCC in the cohort were selected and compared to get adequate power to detect the rumen microbial components, the power analysis based on the microbiome outcomes still indicated a low power (Supplementary Table S2). Due to the animals in this study showed no clinical signs of mastitis, it is speculated that the rumen microbial community may differ more evidently between healthy cows and cows with clinical mastitis. Such speculations need to be confirmed by further studies comparing rumen microbiota between healthy cows and cows with clinical mastitis.
Compared to cows with lower SCC (71,460 ± 3,890 cells/mL), cows with higher SCC (3,107,610 ± 368,100 cells/mL) had lower milk yield, poorer milk composition, and reduced rumen VFA concentrations, but higher rumen bacterial diversity. Although the beta diversity and the predominant rumen bacterial taxa did not vary among cows with different SCC levels, the relative abundance of phyla SR1 and Actinobacteria, unclassified family Clostridiales and genus Butyrivibrio were significantly different among the SCC groups. Proteobacteria and family Succinivibrionaceae were enriched in cows with extremely low SCC (L_SCC group). The rumen microbiota and rumen fermentations are different in the context of different statuses of udder health, however, this association behind the results needs to be further investigated. The results in this study improve our understanding of the rumen microbiota under different health statuses of the mammary gland.
YZ, MX, and JL conceived and designed the study. YZ and MX performed the collection of samples, analyzed the sequencing data, interpreted the data, prepared the figures and tables, and drafted the manuscript. JL helped interpret the data and revise the paper. All authors have read and approved the final manuscript.
This study was financially supported by the China Agriculture Research System (No. CARS-36).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The authors acknowledge the members of the Institute of Dairy Science at Zhejiang University (Hangzhou, China) for their assistance with sampling and analysis.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2018.03217/full#supplementary-material
Arroyo, R., Martín, V., Maldonado, A., Jiménez, E., Fernández, L., and Rodríguez, J. M. (2010). Treatment of infectious mastitis during lactation: antibiotics versus oral administration of lactobacilli isolated from breast milk. Clin. Infect. Dis. 50, 1551–1558. doi: 10.1086/652763
Batavani, R., Asri, S., and Naebzadeh, H. (2007). The effect of subclinical mastitis on milk composition in dairy cows. Iran. J. Vet. Res. 8, 205–211. doi: 10.22099/ijvr.2007.925
Benjamini, Y., and Hochberg, Y. (1995). Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Stat. Soc. Series B. Stat. Methodol. 57, 289–300.
Campbell, J. H., O'Donoghue, P., Campbell, A. G., Schwientek, P., Sczyrba, A., Woyke, T., et al. (2013). UGA is an additional glycine codon in uncultured SR1 bacteria from the human microbiota. Proc. Natl. Acad. Sci. U.S.A. 110, 5540–5545. doi: 10.1073/pnas.1303090110
Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336. doi: 10.1038/nmeth.f.303
Caporaso, J. G., Lauber, C. L., Walters, W. A., Berg-Lyons, D., Huntley, J., Fierer, N., et al. (2012). Ultra-high-throughput microbial community analysis on the illumina HiSeq and MiSeq platforms. ISME J. 6, 1621–1624. doi: 10.1038/ismej.2012.8
Clemente, J. C., Ursell, L. K., Parfrey, L. W., and Knight, R. (2012). The impact of the gut microbiota on human health: an integrative view. Cell 148, 1258–1270. doi: 10.1016/j.cell.2012.01.035
Danfær, A. (1994). Nutrient metabolism and utilization in the liver. Livest. Prod. Sci. 39, 115–127. doi: 10.1016/0301-6226(94)90163-5
Davis, J. P., Youssef, N. H., and Elshahed, M. S. (2009). Assessment of the diversity, abundance, and ecological distribution of members of candidate division SR1 reveals a high level of phylogenetic diversity but limited morphotypic diversity. Appl. Environ. Microbiol. 75, 4139–4148. doi: 10.1128/AEM.00137-09
Edgar, R. C. (2013). UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 10, 996–1000. doi: 10.1038/nmeth.2604
Edgar, R. C., Haas, B. J., Clemente, J. C., Quince, C., and Knight, R. (2011). UCHIME improves sensitivity and speed of chimera detection. Bioinformatics 27, 2194–2200. doi: 10.1093/bioinformatics/btr381
Flint, H. J., and Bayer, E. A. (2008). Plant cell wall breakdown by anaerobic microorganisms from the mammalian digestive tract. Ann. N. Y. Acad. Sci. 1125, 280–288. doi: 10.1196/annals.1419.022
González-Rodríguez, M. C., Gonzalo, C., San Primitivo, F., and Cármenes, P. (1995). Relationship between somatic cell count and intramammary infection of the half udder in dairy ewes. J. Dairy Sci. 78, 2753–2759. doi: 10.3168/jds.S0022-0302(95)76906-5
Gray, C., Strandberg, Y., Donaldson, L., and Tellam, R. K. (2005). Bovine mammary epithelial cells, initiators of innate immune responses to mastitis. Aust. J. Exp. Agr. 45, 757–761. doi: 10.1071/EA05046
Halasa, T., Huijps, K., Osteras, O., and Hogeveen, H. (2007). Economic effects of bovine mastitis and mastitis management: a review. Vet. Q. 29, 18–31. doi: 10.1080/01652176.2007.9695224
Hand, K. J., Godkin, A., and Kelton, D. F. (2012). Milk production and somatic cell counts: a cow-level analysis. J. Dairy Sci. 95, 1358–1362. doi: 10.3168/jds.2011-4927
Hu, W. L., Liu, J. X., Ye, J. A., Wu, Y. M., and Guo, Y. Q. (2005). Effect of tea saponin on rumen fermentation in vitro. Anim. Feed Sci. Technol. 120, 333–339. doi: 10.1016/j.anifeedsci.2005.02.029
Jami, E., and Mizrahi, I. (2012). Composition and similarity of bovine rumen microbiota across individual animals. PLoS ONE 7:e33306. doi: 10.1371/journal.pone.0033306
Jost, T., Lacroix, C., Braegger, C. P., Rochat, F., and Chassard, C. (2014). Vertical mother-neonate transfer of maternal gut bacteria via breastfeeding. Environ. Microbiol. 16, 2891–2904. doi: 10.1111/1462-2920.12238
Laporte, M.-F., and Paquin, P. (1999). Near-infrared analysis of fat, protein and casein in cow's milk. J. Agr. Food Chem. 47, 2600–2605. doi: 10.1021/jf980929r
Li, M., Penner, G. B., Hernandez-Sanabria, E., Oba, M., and Guan, L. L. (2009). Effects of sampling location and time, and host animal on assessment of bacterial diversity and fermentation parameters in the bovine rumen. J. Appl. Microbiol. 107, 1924–1934. doi: 10.1111/j.1365-2672.2009.04376.x
Ma, C., Zhao, J., Xi, X., Ding, J., Wang, H., Zhang, H., et al. (2016). Bovine mastitis may be associated with the deprivation of gut Lactobacillus. Benef. Microbes 7, 95–102. doi: 10.3920/BM2015.0048
Magoč, T., and Salzberg, S. L. (2011). FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics 27, 2957–2963. doi: 10.1093/bioinformatics/btr507
Mao, S. Y., Huo, W. J., and Zhu, W. Y. (2015). Microbiome-metabolome analysis reveals unhealthy alterations in the composition and metabolism of ruminal microbiota with increasing dietary grain in a goat model. Environ. Microbiol. 18, 525–41. doi: 10.1111/1462-2920.12724
McCann, J. C., Wickersham, T. A., and Loor, J. J. (2014). High-throughput methods redefine the rumen microbiome and its relationship with nutrition and metabolism. Bioinform. Biol. Insights 8, 109–125. doi: 10.4137/BBI.S15389
Pope, P. B., Smith, W., Denman, S. E., Tringe, S. G., Barry, K., and Hugenholtz, P. (2011). Isolation of Succinivibrionaceae implicated in low methane emissions from tammar wallabies. Science 333, 646–648. doi: 10.1126/science.1205760
Pyorala, S. (2003). Indicators of inflammation in the diagnosis of mastitis. Vet. Res. 34, 565–578. doi: 10.1051/vetres:2003026
Reynolds, C. K., Aikman, P. C., Lupoli, B., Humphries, D. J., and Beever, D. E. (2003). Splanchnic metabolism of dairy cows during the transition from late gestation through early lactation. J. Dairy Sci. 86, 1201–1217. doi: 10.3168/jds.S0022-0302(03)73704-7
Robert, W. L., Erin, E. C., Congjun, L., Ransom, L. B., and Michael, E. S. (2012). Characterization of the rumen microbiota of pre-ruminant calves using metagenomic tools. Environ. Microbiol. 14,129–139. doi: 10.1111/j.1462-2920.2011.02543.x
Rocchini, D., Hernandez-Stefanoni, J. L., and He, K. S. (2015). Advancing species diversity estimate by remotely sensed proxies: a conceptual review. Ecol. Inform. 25, 22–28. doi: 10.1016/j.ecoinf.2014.10.006
Rodriguez, J. M. (2014). The origin of human milk bacteria: is there a bacterial entero-mammary pathway during late pregnancy and lactation? Adv. Nutr. 5, 779–784. doi: 10.3945/an.114.007229
Schukken, Y. H., Wilson, D. J., Welcome, F., Tikofsky, L. G., and Gonzalez, R. N. (2003). Monitoring udder health and milk quality using somatic cell counts. Vet. Res. 34, 579–596. doi: 10.1051/vetres:2003028
Seegers, H., Fourichon, C., and Beaudeau, F. (2003). Production effects related to mastitis and mastitis economics in dairy cattle herds. Vet. Res. 34, 475–491. doi: 10.1051/vetres:2003027
Shabat, S. K. B., Sasson, G., Doron-Faigenboim, A., Durman, T., Yaacoby, S., Miller, M. E. B., et al. (2016). Specific microbiome-dependent mechanisms underlie the energy harvest efficiency of ruminants. ISME J. 10, 2958–2972. doi: 10.1038/ismej.2016.62
Shen, J. S., Chai, Z., Song, L. J., Liu, J. X., and Wu, Y. M. (2012). Insertion depth of oral stomach tubes may affect the fermentation parameters of ruminal fluid collected in dairy cows. J. Dairy Sci. 95, 5978–5984. doi: 10.3168/jds.2012-5499
Singh, K. M., Tripathi, A. K., Pandya, P. R., Parnerkar, S., Rank, D. N., Kothari, R. K., et al. (2012). Methanogen diversity in the rumen of Indian Surti buffalo (Bubalus bubalis), assessed by 16S rDNA analysis. Res. Vet. Sci. 92, 451–455. doi: 10.1016/j.rvsc.2011.03.022
Sulak, M., Sikorova, L., Jankuvova, J., Javorsky, P., and Pristas, P. (2012). Variability of actinobacteria, a minor component of rumen microflora. Folia. Microbiol. (Praha) 57, 351–353. doi: 10.1007/s12223-012-0140-7
Tsenkova, R., Atanassova, S., Kawano, S., and Toyoda, K. (2001). Somatic cell count determination in cow's milk by near-infrared spectroscopy: a new diagnostic tool. J. Anim. Sci. 79, 2550–2557. doi: 10.2527/2001.79102550x
Xue, M. Y., Sun, H. Z., Wu, X. H., Guan, L. L., and Liu, J. X. (2018). Assessment of rumen microbiota from a large dairy cattle cohort reveals the pan and core bacteriomes contributing to varied phenotypes. Appl. Environ. Microbiol. 84, e00970–e00918. doi: 10.1128/AEM.00970-18
Young, W., Hine, B. C., Wallace, O. A. M., Callaghan, M., and Bibiloni, R. (2015). Transfer of intestinal bacterial components to mammary secretions in the cow. PeerJ. 3:e888. doi: 10.7717/peerj.888
Zhu, Z. G., Kristensen, L., Difford, G. F., Poulsen, M., Noel, S. J., Ai-Soud, et al. (2017). Changes in rumen bacterial and archaeal communities over the transition period in primiparous Holstein dairy cows. J. Dairy Sci. 101, 9847–9862. doi: 10.3168/jds.2017-14366
Keywords: somatic cell counts, subclinical mastitis, milk production, rumen fermentation, rumen microbiota
Citation: Zhong Y, Xue M and Liu J (2018) Composition of Rumen Bacterial Community in Dairy Cows With Different Levels of Somatic Cell Counts. Front. Microbiol. 9:3217. doi: 10.3389/fmicb.2018.03217
Received: 18 October 2018; Accepted: 11 December 2018;
Published: 24 December 2018.
Edited by:Garret Suen, University of Wisconsin-Madison, United States
Reviewed by:Stephanie A. Metzger, University of Wisconsin-Madison, United States
Erika Korzune Ganda, Cornell University, United States
Copyright © 2018 Zhong, Xue and Liu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jianxin Liu, firstname.lastname@example.org
†These authors have contributed equally to this work