- 1Department of Animal Science and Technology, National Taiwan University, Taipei City, Taiwan
- 2Department of Veterinary Medicine, College of Veterinary Medicine, National Chung Hsing University, Taichung City, Taiwan
- 3Department of Wood Based Materials and Design, National Chiayi University, Chiayi City, Taiwan
- 4Graduate Institute of Biomedical and Pharmaceutical Science, Fu-Jen Catholic University, New Taipei City, Taiwan
- 5Biotools Co. Ltd., New Taipei City, Taiwan
Introduction: There are differences in the gut microbiome and metabolome when the host undergoes different physical or pathological conditions. However, the inter-relationship of microbiome and metabolome biomarkers to potentially promote the health of dairy cows needs to be studied. Further, the development of next-generation probiotics for dairy cattle health promotion has not been demonstrated.
Objective: In the present study, we identified the microbiome and metabolome biomarkers associated with healthy cows.
Methods: We analyzed the relationships of the ruminal microorganism profile and metabolites between healthy and mastitis lactating dairy cows. The roles of bacterial biomarker were further verified by in vitro fermentation and cow-to-mouse fecal microbiota transplantation (FMT).
Results: Two species, Ruminococcus flavefaciens and Bifidobacterium longum subsp. longum, and six rumen metabolites were positively correlated with healthy cows by Spearman’s correlation analysis. Through in vitro ruminal fermentation, inoculating R. flavefaciens and B. longum subsp. longum showed the upregulation of the levels of putrescine, xanthurenic acid, and pyridoxal in the mastitis ruminal fluid, which confirmed the inter-relationships between these microbiota and metabolites associated with healthy cows. Further, we verified the role of R. flavefaciens and B. longum subsp. longum in promoting health by FMT. The administration of R. flavefaciens and B. longum subsp. longum reduced the death rate and recovered the bodyweight loss of germ-free mice caused by FMT mastitis feces.
Discussion: We provided evidence that the bacterial biomarkers alter downstream metabolites. This could indirectly indicate that the two bacterial biomarkers have the potential to be used as next-generation probiotics for dairy cattle, although it needs more evidence to support our hypothesis. Two species, R. flavefaciens and B. longum subsp. longum, with three metabolites, putrescine, xanthurenic acid, and pyridoxal, identified in the ruminal fluid, may point to a new health-promoting and disease-preventing approach for dairy cattle.
1. Introduction
Bovine rumen possesses a highly diverse population of microorganisms, including bacteria, protozoa, archaea, and fungi, which degrade and ferment the plant materials into digestible compounds (1). Proteobacteria, Firmicutes, and Bacteroidetes are the dominant phyla of the Kingdom Eubacteria in dairy cattle rumen (2, 3). However, the bacterial abundance in the digestive tract may fluctuate because of age, nutrients, and other factors related to lifestyle (4–6) or immune status (7). Changes in the gut microbiota play an influential role in the performance of the animals (3, 8) as well as contribute to the development of diseases including mastitis (7, 9). Gut dysbiosis also plays an important role in the health of dairy cows due to the difference and diversity of the gastrointestinal microbiota that competes for nutrients, regulate the immune system, and produce metabolites (10). Different microbiota found in feces (9) and ruminal fluid (7) of healthy and mastitis cows, and induction of mastitis in germ-free (GF) mice by fecal microbiota transplantation (FMT) suggested that bovine mastitis is not necessary a local infection of the mammary glands (9). Many evidences have shown that the stress factors such as high concentrate feeding or heat stress, could disturb the rumen microbiota and upregulate the level of lipopolysaccharide (LPS) (11), resulting in changes in the permeability of the rumen epithelial layer (12). The rumen-derived LPS could enter the mammary gland via blood circulation and further impair the blood-milk barrier, leading to inflammation of the mammary gland in cows (12, 13). Additionally, Zhao et al. (13) suggested that ruminal dysbiosis-derived low-grade endotoxemia could cause mastitis and worsen pathogen-induced mastitis by damaging host anti-inflammatory enzymes.
The metabolites, derived from fermentation by rumen microorganisms, are considered as a downstream outcome, illustrating the interaction between microorganisms, hosts, and the microenvironment. Metabolomics has been employed to evaluate the quality of milk (14) and search for new biomarkers for disorders (7). Chuang et al. (7) identified seven rumen fluid metabolites that changed between healthy and mastitis cows, which could be used as potential biomarkers for the diagnosis of mastitis.
The probiotics currently available to farm animals are generally limited to a narrow range of organisms. Characterizing the gut microbiota and metabolites is a novel preventive or therapeutic approach for the development of next-generation probiotics (15, 16). Therefore, the comprehensive description of the ruminal microbiota and metabolome and their roles in health and disease is crucial. Although the ruminal microbial and metabolomic structure in lactating dairy cows with mastitis has been studied (7), the inter-relationship of the microbiome and metabolome biomarkers in promoting health in dairy cows has never been confirmed. The role of bacterial biomarkers in health promotion and disease prevention also remains unknown. Thus, this study first identified the microbiome and metabolome biomarkers associated with healthy cows by evaluating the relationships among the ruminal microbial profile, metabolites, and mastitis outcomes. We then verified the inter-relationship of the microbiome and metabolome biomarkers and the role of bacterial biomarkers in health promotion and disease prevention through in vitro ruminal fermentation and cow-to-mouse fecal microbiota transplantation, respectively.
2. Materials and methods
2.1. Potential ruminal microorganism biomarkers and related metabolites for healthy dairy cows
2.1.1. Animals and sample collection
Thirty lactating Holstein dairy cows with 120–240 milk production days and an average age of 3.53 ± 0.67 years from a commercial farm were involved in the present study. All cows were under the same management, receiving total mixed ration (TMR) feeding and water ad libitum, and milked twice per day. After the outcomes of veterinary diagnosis, raw milk test with California mastitis test (CMT) and somatic cell counts (SCC), and serum proinflammatory cytokines, 15 healthy cows and 15 cows with clinical mastitis were selected for microbiomic and metabolomic analysis. Cows with one quarter milk showed positive reaction by CMT, SCC ≥ 1,000,000 cells/mL, and elevated cytokines in serum were defined as mastitis cows. On the other hand, cows with negative CMT reaction, SCC < 200,000 cells/mL, and no specific cytokines were included in health group (17). Milk, ruminal fluid, and blood samples were collected 2 h after morning feeding (4 h after morning milking), according to previously described methods (7).
2.1.2. Analysis of somatic cell counts and N-acetyl-β-D-glucosaminidase in milk
The California mastitis test kit (ImmuCell Corp., Portland, ME, United States) was used to analyze milk CMT reaction on the farm and was followed the manufacturer’s instructions. The SCC of quarter milk samples was conducted by a Fossomatic FC instrument (Foss Electric, Hillerød, Denmark). Milk N-acetyl-β-D-glucosaminidase (NAGase) activity was measured with a fluoro-optical method described by Kalmus et al. (18).
2.1.3. Analysis of serum cytokines
The commercial enzyme-linked immunosorbent assay kits (Bovine TNF-alpha and IL-6 DuoSet ELISA kit, R&D Systems, Minneapolis, MN, United States) were used to measure the levels of tumour necrosis factor (TNF)-α and interleukin (IL)-6.
2.1.4. Microbiome analysis
Total genomic DNA was extracted from ruminal fluid samples using the bead-beating method (19). The microbiome analysis adopted the method described by Chuang et al. (7) using the Illumina HiSeq 2,500 PE250 platform (20–25). The representative sequence for each operational taxonomic units (OTUs) was analyzed through taxonomic annotation (26, 27) and determined the alpha diversity (Chao1 richness estimator and Shannon’s diversity index). We used partial least squares discriminant analysis (PLS-DA), and the linear discriminant analysis (LDA) effect size (LEfSe) algorithm to analyze the data (28). The false discovery rate (FDR) was used to carry out multiple testing for the correction of the p value using the Benjamini–Hochberg procedure.
2.1.5. Metabolite analysis
The ruminal fluid sample preparation and metabolite analysis adopted the method described by Chuang et al. (7). The orthogonal PLS-DA (oPLS-DA) model with MetaboAnalyst 5.01 (29) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways using the KEGG database were also analyzed.
2.2. Verification of potential biomarkers related to the health of dairy cows
2.2.1. Bacterial preparation
Bifidobacterium longum subsp. longum (BCRC 14664) and Ruminococcus flavefaciens (DSM 25089) were purchased from the Bioresource Collection and Research Center (BCRC, Food Industry Research and Development Institute, Hsinchu, Taiwan) and Deutsche Samulung von Mikroorganisem und Zelkultruen (DSMZ, Braunscheig, Germany), respectively. B. longum subsp. longum was activated three times using the Lactobacilli MRS broth (Lactobacilli de Man, Rogosa and Sharpe broth, Acumedia, Lansing, MI, United States) at 37°C before subsequent analysis. R. flavefaciens was cultured under an anaerobic environment at 37°C using the medium formulated by Wang (30).
2.2.2. In vitro ruminal fermentation
After the outcomes of veterinary diagnosis, CMT, SCC, NAGase, and serum proinflammatory cytokines, 5 healthy cows (H group) and 5 cows with clinical mastitis (M group) under a similar milk production stage (the first 40- to 60 day lactation) and an average age of 2.75 ± 0.71 years old from same commercial farm mention above were selected as rumen fluid donors. After the morning feeding, 1,500 mL of rumen fluid was obtained as described above. The samples were filtered with four layers of cheesecloth and placed in a flask with thermal insulation (39°C) before being transported to the laboratory. The artificial saliva, prepared according to the description of Menke and Steingass (31), was combined with rumen fluid in a non-oxygen atmosphere. The 40 mL mixture was filled into a 100 mL serum bottle with CO2 containing 0.4 mg of feed subtract (fresh TMR prepared by the farm) and bacterial culture. Both the H and M groups of rumen fluid were further divided into 4 sub-groups defined as followed: A, with 1 mL of sterile ddH2O as the control; B, with 1 mL of B. longum subsp. longum (106 CFU/mL) bacterial culture; C, with 1 mL of R. flavefaciens (106 CFU/mL) bacterial culture; D, with 1 mL of each bacterial culture. For fermentation, the bottles were capped and incubated in a shaking incubator (120 rpm) at 39°C. The fermented fluid was collected at 0, 3, and 12 h.
2.2.3. Qualitative metabolites
The ruminal fluid samples were centrifuged at 13,400 × g for 15 min. The supernatants were analyzed using a Shimadzu LC-20A high-performance liquid chromatography (HPLC) system (Shimadzu, Kyoto, Japan) coupled to a linear ion trap-Orbitrap mass spectrometer (LTQ Orbitrap Velos, Thermo Fisher Scientific, Waltham, MA, United States). The standards of metabolites were used for the qualitative analysis.
2.2.4. Fecal microbiota transplantation
Fresh fecal samples from 15 mastitis and 15 healthy cows which were the same as Section 2.1.1. mentioned were, respectively, collected. The preparation procedure followed the method described by Ma et al. (9). GF mice were obtained and housed according to animal care regulations in the germ-free animal facility at the Animal Resource Center, National Taiwan University (Taipei, Taiwan). A total of 15 female adult (8 week-old) C57BL/6 J mice were randomly divided into three groups, which received 0.3 mL fecal supernatant from (i) healthy cows (Control group), (ii) mastitis cows (Mastitis group), or (iii) mastitis cows, plus 108 CFU per day of B. longum subsp. longum and R. flavefaciens administration (M + BR group) for 4 weeks. The animals had measured bodyweight per week to determine the changes in body weight during the experiment period. The three groups of mice were caged in different gnotobiotic isolators after FMT to prevent cross-contamination.
2.3. Statistical analysis
All phenotypic and next-generation sequencing (NGS) data were analyzed with a nonparametric Mann–Whitney U test to identify significant differences between groups. Spearman’s correlation analysis was used to conduct the correlation between the relative abundance of biomarkers and metabolites. Statistical Analysis System v9.4 (SAS Institute Inc., Cary, NC, United States) and R software were used for all statistical analysis.
3. Results
3.1. Healthy status of tested dairy cows
Since the healthy status is crucial for this study, the mastitis cows were selected not only by veterinary diagnosis and CMT, but also by the milk SCC and NAGase as well as serum proinflammatory cytokines. The selected 15 mastitis cows demonstrated significantly higher milk SCC (p < 0.05) and NAGase (p < 0.05) than those of the 15 healthy counterparts (Figure 1). The serum IL-6 in the mastitis cows was also higher than that of healthy cows, which provided a solid foundation for the current study.
3.2. Beta diversity illustrates the dissimilar gut microbiota harbored in ruminants between healthy and mastitis cows
After verifying the healthy status of the cows, we analyzed the ruminal microbiota by NGS. 16S rRNA analysis revealed a total of 778,270 and 794,717 effective tags from 1,200,157 and 1,233,089 raw paired-end reads from healthy and mastitis groups, respectively. The Venn diagram in Figure 2A showed that 2,307 OTUs were identical between the groups, with 75 and 95 unique OTUs for the healthy and mastitis groups, respectively. The alpha diversity (Chao1 richness estimator and Shannon’s diversity index) revealed no significant difference (p > 0.05) between the two groups (Supplementary Figure S1). The top 10 dominant taxa at the genus level, which covered 58% of the total genus level results, were identical between groups but with different proportions (Supplementary Figure S2). Additional beta diversity analysis separated the healthy and mastitis groups using the PLS-DA plot (Figure 2B). PLS1 and PLS2 explained 6.59 and 5.78%, respectively, of the variation in gut microbiota composition, illustrating the dissimilar gut microbiota harbored in ruminants. The predicted phenotypes showed that the healthy group possessed higher relative abundance in stress-tolerant, anaerobic, and Gram-positive bacteria, and lower abundance in Gram-negative and potentially pathogenic bacteria compared to the mastitis counterpart (Figure 2C).
Figure 2. Ruminal bacteria and archaea composition identified by 16S rRNA sequencing of healthy and mastitis cows. (A) Venn diagram illustrating 2,307 operational taxonomic units (OTUs) of core microbiota identified both in healthy and mastitis cows. (B) Partial least squares discriminant analysis (PLS-DA) plot based on the relative abundance of OTUs indicates a significantly different composition of healthy versus mastitis cows. Ellipses represent 95% confidence intervals for each group. (C) Predicted phenotypes. (D) Significant differential biomarkers were identified using the LEfSe algorithm. (E) Spearman’s correlation test between mastitis markers and gut microbial biomarkers at the species level. Each cell was colored corresponding to the Spearman’s correlation results. Significant difference: *p < 0.05 and **p < 0.01. (F) Significant relative abundance of differential biomarkers.
3.3. Identification of the critical ruminal bacterial biomarkers
The beta diversity analysis and predicted phenotypes revealed the dissimilar ruminal bacteria existing between the two groups. Thus, the critical taxa associated with healthy and mastitis groups were then analyzed using the LEfSe algorithm with LDA > 2.5 as the bacterial biomarkers. The results identified 37 influential taxonomic clades, including 17 genera and 3 species (Figure 2D). The most impacted taxa in the healthy group were 12 genera (Ruminococcaceae UCG 014, Eubacterium coprostanoligenes group, Eubacterium ruminantium group, Ruminococcus 1, Syntrophococcus, Dialister, Pseudobutyrivibrio, Desulfovibrio, Lachnoclostridium 12, Ruminococcaceae UCG 007, Peptostreptococcus, Mitsuokella), and 2 species (R. flavefaciens and B. longum subsp. longum). Four genera (Prevotella 1, Prevotellaceae UCG001, Prevotellaceae UCG003, Fibrobacter) and one species (unidentified rumen bacterium RNF82) were the critical taxa in the mastitis group.
3.4. Correlation of mastitis parameters with the bacterial biomarkers
After identifying the bacterial biomarkers in both groups, we illustrated the correlation of mastitis parameters (SCC, NAGase, IL-6) with the bacterial biomarkers at the species level (Figure 2E). The species enriched in the healthy group, R. flavefaciens and B. longum subsp. longum, were negatively correlated with the levels of SCC and NAGase. Conversely, the species enriched in the mastitis group, the unidentified rumen bacterium RNF82, demonstrated positive correlations with the levels of SCC and NAGase. The relative abundance of the bacterial biomarkers related to the mastitis group was paralleled with the above findings. The mastitis group demonstrated a significantly lower relative abundance in the R. flavefaciens and B. longum subsp. longum (p < 0.05) (Figure 2F).
Phylogenetic investigation of communities by the reconstruction of unobserved states (PICRUSt) was applied to investigate the mastitis-associated functional profiles of microbiome communities (Figure 3). In the mastitis group, 13 functions were enriched, including xenobiotics biodegradation and metabolism, energy metabolism, amino acid metabolism by cysteine and methionine metabolism, neurodegenerative and infectious diseases, replication and repair, transcription and translation. The remaining 18 pathways were depleted, including nucleotide metabolism, lipid metabolism, glycan biosynthesis and metabolism, environmental adaptation, and vitamin B-related metabolic pathways.
Figure 3. Comparison of the relative abundance of the PICRUSt functional prediction of the ruminal microbiota between healthy and mastitis groups. The results are presented as mean ± SEM (n = 15). Distinct gene categories were selected according to significant differences in the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway level 2 (Mann–Whitney U test, p < 0.01).
3.5. Metabolomics analysis revealed a dissimilar ruminal metabolite composition
In total, 1,181 compounds were identified through blasting, matching the mzCloud online database among 12,709 practicable peaks. A volcano plot showing the metabolite profile of the statistical significance (VIP > 1, FDR adjusted p < 0.05 threshold) against fold change revealed the ruminal metabolites with significant differences between the two groups (Figure 4A). Additionally, the oPLS-DA plot showed a clear separation between the healthy and mastitis groups, suggesting a dissimilar ruminal metabolite conformation (Figure 4B). Among identified metabolites, 65 metabolites were significantly different between the two groups (VIP > 1, FDR adjusted p < 0.05), of which, 29 were significantly lower in the mastitis group than those in the healthy group (Supplementary Table S1). Further analysis of the metabolic pathway using KEGG (the second level) identified 21 (methionine, putrescine, proline, piperidine, 5-hydroxyindoleacetic acid, N-(2-phenylethyl)-acetamide, 1-pyrroline, 3-acetamidopropanal, nervonyl carnitine, asparaginyl-alanine, triphenylsilanol, 4-aminophenylalanine, 2-phenylbutyramide, pyrophaeophorbide, linoleoyl ethanolamide, 6-pentadecyl salicylic acid, tyr-OEt, 1-(3-aminopropyl)-4-aminobutanal, hexamethylene bisacetamide, 1,5-diphenylcarbohydrazide, and ansamitocin P3) and 16 (carnitine, (alpha)-JWH 073 N-(3-hydroxybutyl) metabolite-d5, (6aR, 11aR)-3-hydroxy-8,9-dimethoxypterocarpan, 2-ethylpyrazine, erucamide, 2-amino-octadecanoic acid, Corey PG-lactone diol, juvenile hormone I, 2,3-dinor-8-iso-PGF2a, 2-acetylpyrazine, hydrocortamate, 9-hydroperoxy-10E, 12-octadecadienoic acid, ethyl 2-furanpropionate, oleoyl ethyl amide, osmundalactone, and 9-decynoic acid) metabolites involved in amino acid and lipid metabolism, respectively. Other pathways such as the metabolism of cofactors and vitamin, nucleotide metabolism, xenobiotic biodegradation metabolism, nicotinate metabolism, carbohydrate metabolism, energy metabolism, nervous system, and replication and repair were also identified.
Figure 4. Different compositions of ruminal metabolites of healthy and mastitis cows. (A) Volcano plot of 1,181 compounds with log-transformed adjusted p-values and fold change. The green and red dots indicate significantly higher metabolites in the healthy and mastitis groups, respectively. (B) Orthogonal PLS-DA (oPLS-DA) plot based on the 1,181 compounds indicates significantly different metabolite compositions of the healthy and inflammatory groups. Ellipses represent 95% confidence intervals for each group. Every dot represents a single individual cow. (C) Reciprocal interrelationships between the ruminal microbiota and metabolome by Spearman’s correlation test shown at the species level. Orange and blue colors indicate positive and negative correlation coefficients, respectively. Symbols indicate the significant correlation between metabolites and biomarkers (*p < 0.05 and **p < 0.01). (D) The significant relative peak area and (E) concentration of metabolites positively correlated with Ruminococcus flavefaciens and Bifidobacterium longum subsp. longum.
3.6. Correlation of the ruminal microbiota and metabolome
The reciprocal inter-relationships between 20 ruminal microbial biomarkers and 65 metabolites, which were significantly different between the two groups, were analyzed with the Spearman’s correlation test. All the metabolites were significantly correlated with some of the critical bacterial biomarkers (p < 0.05) (Supplementary Figure S3). Among them, 3-methylcytosine, 1-(3-aminopropyl)-4-aminobutanal, putrescine, pyridoxal, xanthurenic acid, and Tyr-OEt, are related to amino acid metabolism, replication and repair, and metabolism of cofactors and vitamin (Supplementary Table S1) were significantly positively correlated with the bacterial species biomarkers in the healthy group, R. flavefaciens and B. longum subsp. longum (p < 0.05) (Figure 4C). The unidentified rumen bacterium RNF82, the bacterial biomarkers in the mastitis group, were negatively correlated with cotinine N-oxide, and Corey PG-lactone diol (p < 0.05) and positively correlated with 2-amino-octadecanoic acid, erucamide, and ethosuximide M7 (Figure 4C). By further quantifying the ruminal metabolites, the higher relative peak area of three metabolites, xanthurenic acid, pyridoxal, and putrescine in the healthy group (Figure 4D) demonstrated significantly higher concentrates in the ruminal fluid samples compared with the mastitis counterpart, verifying the ruminal metabolomic finding (Figure 4E).
3.7. Verification of the health-promoting effect of Ruminococcus flavefaciens and Bifidobacterium longum subsp. longum by FMT GF mice
To verify the healthy promoting effect of R. flavefaciens and B. longum subsp. longum, fecal microbiota from the 15 mastitis and 15 healthy cows were, respectively, pooled and inoculated into GF mice. The results showed that the mice transplanting gut microbiota from mastitis cows could decrease the survival rate and bodyweight gain compared with that from healthy cows (Figure 5). Administration of R. flavefaciens and B. longum subsp. longum could increase body weight gain and reduce mortality rate.
Figure 5. Survival rate and body weight gain after the fecal microbiota transplantation (FMT) test and bacterial biomarker supplements.
3.8. Verification of inter-relationships between microbiota and metabolome in vitro
We further investigated the inter-relationships between microbiota and metabolome by in vitro ruminal fermentation. First, the ruminal fluid pH value during a 12 h in vitro fermentation was above 6.0 in both the healthy and mastitis group. The mean pH was 6.72 (HA), 6.84 (HB), 6.83 (HC), 6.87 (HD), 6.95 (MA), 6.91 (MB), 7.00 (MC), and 7.02 (MD). The basic composition of the ruminal fluid in all groups had no significant change except NH3-N (Table 1). R. flavefaciens and B. longum subsp. longum addition demonstrated the trend to upregulate the levels of putrescine, xanthurenic acid, and pyridoxal in the mastitis ruminal fluid (Figure 6), which confirmed the inter-relationships between microbiota and metabolome. Additionally, the levels of total volatile fatty acids (VFA) were upregulated after inoculation with R. flavefaciens (HC and MC groups) and R. flavefaciens + B. longum subsp. longum (HD and MD groups) (Table 2). As supplement of B. longum subsp. longum, the relative amount of lactate would increase (HB and MB groups).
Table 1. Supplement of Ruminococcus flavefacians and Bifidobacterium longum subsp. longum on the digestibility of dietary nutrients in in-vitro fermentation.
Figure 6. Effect of bacterial biomarker supplements on (A) putrescine, (B) xanthurenic acid, and (C) pyridoxal concentration during the 12 h in vitro fermentation.
Table 2. Effects of bacterial biomarker supplements on the concentration of volatile fatty acids (VFA) and lactic acid in healthy and mastitis groups after 12 h of in vitro fermentation.
4. Discussion
In the present study, we verified the inter-relationship of microbiome and metabolome biomarkers to potentially promote the health of dairy cows. First, a stable and resilient core microbiota in the ruminal fluid of the lactating cow w/wo mastitis was observed based on the results of the Venn diagram, Shannon, Chao1, and the relative abundance of the taxa at the different levels. The genera Prevotella 1, Ruminococcaceae NK4A214 group, and Christensenellaceae R7 group, the top three predominant genera in both healthy and mastitis groups, were also the abundant genera in the rumens of lactating cows, dry period cows (32), yak (33), and mastitis cows (7). These three genera play important roles in protein degradation and lipid biohydrogenation (3), SCFA production by the breakdown of fibrous plants (34), and microbial inhibiting activities (35). The Ruminococcaceae NK4A214 group was also positively correlated to milk total solids in lactating cows (36).
Although the bovine rumen demonstrated a highly similar microbial composition, the differences in taxa and individual bacterial abundance still existed between the lactating cows w/wo mastitis, which could effectively distinguish these two groups using the PLS plot and unweighted UniFrac. The ruminal microbiota of mastitis cows was characterized by high Gram-negative and potentially pathogenic bacteria. KEGG pathways with higher abundance in mastitis cows, including xenobiotics acid biodegradation and metabolism, energy metabolism, replication and repair, transcription, translation, and infectious disease, can be related to the inflammation and mucosa repairing in cows (37). The genera Prevotella (Prevotella 1, Prevotellaceae UCG001, Prevotellaceae UCG003) and Fibrobacter were the biomarkers associated with mastitis cows. Prevotella with diverse isoforms is crucial for ruminal fermentation (38); however, the high abundance of this genus was associated with high-grain feed (39) as well as acidosis (40). The studies in humans and animals have connected the increased abundance of Prevotella species at mucosal sites with localized and systemic inflammation disease due to enhancing T helper type 17 (Th17)-mediated mucosal inflammation via IL-8, IL-6, and CCL20 stimulation (41). The increasing serum IL-6 in mastitis cows in this study may partially support the ruminal dysbiosis leading to systemic inflammatory effect.
Conversely, the ruminal microbiota associated with healthy cows was categorized by high-stress tolerant, Gram-positive, and SCFA producing bacteria. This finding was later confirmed by in vitro ruminal fermentation. Increasing VFA and SCFA with healthy ruminal liquid were due to the high SCFA producing bacteria. The reduction of SCFA producing bacteria in ruminal fluid (7) and feces (9) of mastitis cows has been reported. Higher KEGG pathway in lipid and carbohydrate metabolisms in healthy cows also suggested reduced lipid and carbohydrate metabolic activities of the gut microbiota in mastitis cows. Our finding was in line with a study indicating that carbon metabolism was less abundant in mastitis cows (9). Downregulation of carbohydrate metabolism may alter the glucose and carbohydrate balance in the body (42), which affects the energy for maintenance, growth, and production in farm animals (43).
The biomarkers identified in the healthy group, including the genera Ruminococcus, Eubacterium, Lachnoclostridium, and Pseudobutyrivibrio, known as cellulose and fiber degraders (44), associated with high-yield cows (45). The genus Syntrophococcus has been reported to utilize sugars and H2-CO2-using methanogens as electron donors to produce acetate and as an electron acceptor, respectively (46). At the species level, R. flavefaciens could modify the abundance of other cellulolytic bacterial populations (47) and improve the feed efficacy for ruminants (48). Another species, B. longum subsp. longum, a biomarker in the healthy group, was also reported to stabilize gut microbiota and improve the intestinal environment (49). Both species recognized in the healthy cow could be potential probiotics to promote animal health, which warrants further investigation.
From the in vitro ruminal fermentation, additional inoculation of R. flavefaciens and B. longum subsp. longum could significantly impact the levels of NH3-N, total VFA, and VFA profiles in both healthy and mastitis groups. NH3-N is the main nitrogen source used by microbes to synthesize amino acid and peptide bonds for growth (50). Ruminobacter spp. is a hyper-ammonia producing (HAP) bacteria (51). The increase in NH3-N could be explained by the inoculation of R. flavefaciens, which leads to an increased population of HAP bacteria and deaminase activity. Total VFA and VFA profiles are important products of the bacterial fermentation activity in the rumen, which have emerged as key regulators in intestinal and energy homeostasis regulation (52). R. flavefaciens participates in the butyrate metabolic pathway (53, 54). Upregulating the VFA concentration has been reported in repeated ruminal dosing of R. flavefaciens in dairy cows (47). The increase in VFA suggested an upregulating deamination activity (55) by the addition of R. flavefaciens. The increase in lactate is expected with additional B. longum subsp. longum, a lactic acid producer. The administration of lactic acid bacterial probiotics is thought to help rumen microbiota adapt to the presence of lactic acid (56) and prevent lactate accumulation in the rumen (57). Nevertheless, the increase in VFA did not have much physiological impact as the rumen pH due to shifting the microbiota to lactate-consuming bacteria.
FMT verified the inter-relationship among gut dysbiosis, systemic inflammatory effect, and health-promoting ability of two microbial biomarkers, R. flavefaciens and B. longum subsp. longum. Ma et al. (9) found that FMT from diseased cows caused mastitis-like symptoms in mice by shifting the murine intestinal microbiota. Although the acute inflammation led to mice mortality after FMT with mastitis feces, which could not provide solid evidence between mastitis and the dysbiosis of ruminal microbiota, the findings confirmed the impact of gut microbiota as one potential parameter affecting dairy cow health. The increase in the survival rate after FMT of mastitis feces with R. flavefaciens and B. longum subsp. longum also supported the potential efficacy of microbial biomarkers as probiotic treatment, which may point to a new health-promoting and disease-preventing approach.
Besides the microbiota, bacterial products, in turn of metabolome, were also a key factor involved in bovine health and systemic disease outcomes. Three metabolites (i.e., putrescine, xanthurenic acid, and pyridoxal) were verified by further in vitro ruminal fermentation and HPLC qualitative analysis, which were positively correlated with two species biomarkers, R. flavefaciens and B. longum subsp. longum. Putrescine, a biogenic amine produced from the decarboxylation of amino acids by decarboxylase in certain intestinal microorganisms (58), has been found in rumen fluid of healthy animals, which was in line with our findings. Putrescine has antioxidant and anti-inflammatory attributes (59) and is involved in the growth of tissues and organs (60). Xanthurenic acid, a non-indolic catabolite of tryptophan and a metabolite of the kynurenine pathway (61), demonstrated profound effects on the gut microbial composition, host-microbiome interface, and host immune system–intestinal microbiota interactions. Pyridoxal is one of the natural forms available of vitamin B6, supplied either in the diet or by rumen or intestinal symbiosis for bovine species. Vitamin B6 participates in DNA, RNA and protein synthesis. Mastitis is associated with a vitamin B metabolism disorder in intestinal microbiota (9). The upregulated putrescine, xanthurenic acid, and pyridoxal in ruminal fluid after adding R. flavefaciens and B. longum subsp. longum suggested a health effect on modulating intestinal homeostasis and damage repair.
5. Conclusion
Although the bovine rumen possesses a strong core microbial composition, we proved that minor microbiota shifting caused by mastitis could affect the health of dairy cows. This influence is not only because of the rumen microbiota but the downstream microbiome produced by microbiota also plays an important role in health. To the best of our knowledge, this study is the first to verify the inter-relationship of microbiome and metabolome biomarkers for the potential to promote health in dairy cows. Two species, R. flavefaciens and B. longum subsp. longum, with three metabolites, putrescine, xanthurenic acid, and pyridoxal, were identified in the ruminal fluid, which may point to a new direction to promote health and prevent disease in dairy cattle.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.
Ethics statement
The animal studies were approved by Institutional Animal Care and Use Committee of National Taiwan University (IACUC approval no: NTU-107-EL-00221). The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent was obtained from the owners for the participation of their animals in this study.
Author contributions
S-TC and M-JC conceived and designed the experiments. J-CH, S-TC, K-YL, Y-TH, and S-TH performed the experiments. J-CH, K-YL, Y-TH, and S-HC analyzed the data. J-CH, S-TC, and M-JC wrote and revised the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This research was funded by the Ministry of Science and Technology of Taiwan (MOST 105-2313-B-002-041-MY3 and MOST 109-2321-B-002-054).
Conflict of interest
S-HC was employed by Biotools Co. Ltd.
The remaining 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.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2023.1228086/full#supplementary-material
Footnotes
References
1. Flint, HJ. The rumen microbial ecosystem—some recent development. Trends Microbiol. (1997) 5:483–8. doi: 10.1016/S0966-842X(97)01159-1
2. Jami, E, Israel, A, Kotser, A, and Mizrahi, I. Exploring the bovine rumen bacterial community from birth to adulthood. ISME J. (2013) 7:1069–79. doi: 10.1038/ismej.2013.2
3. Jami, E, and Mizrahi, I. Composition and similarity of bovine rumen microbiota across individual animals. PLoS One. (2012) 7:e33306. doi: 10.1371/journal.pone.0033306
4. Magne, F, Gotteland, M, Gauthier, L, Zazueta, A, Pesoa, S, Navarrete, P, et al. The Firmicutes/Bacteroidetes ratio: a relevant marker of gut dysbiosis in obese patients? Nutrients. (2020) 12:1474. doi: 10.3390/nu12051474
5. Mariat, D, Firmesse, O, Levenez, F, Guimarăes, V, Sokol, H, Doré, J, et al. The Firmicutes/Bacteroidetes ratio of the human microbiota changes with age. BMC Microbiol. (2009) 9:123. doi: 10.1186/1471-2180-9-123
6. Wu, D, Vinitchaikul, P, Deng, M, Zhang, G, Sun, L, Gou, X, et al. Host and altitude factors affect rumen bacteria in cattle. Braz J Microbiol. (2020) 51:1573–83. doi: 10.1007/s42770-020-00380-4
7. Chuang, ST, Li, KY, Tu, PW, Ho, ST, Hsu, CC, Hsieh, JC, et al. Investigating the reciprocal inter-relationships among the ruminal microbiota, metabolome, and mastitis in early lactating Holstein dairy cows. Animals. (2021) 11:3108. doi: 10.3390/ani11113108
8. Min, BR, Gurung, N, Shange, R, and Solaiman, S. Potential role of rumen microbiota in altering average daily gain and feed efficiency in meat goats fed simple and mixed pastures using bacterial tag-encoded FLX amplicon pyrosequencing. J Anim Sci. (2019) 97:3523–34. doi: 10.1093/jas/skz193
9. Ma, C, Sun, Z, Zeng, B, Huang, S, Zhoa, J, Zhang, Y, et al. Cow-to-mouse fecal transplantations suggest intestinal microbiome as one cause of mastitis. Microbiome. (2018) 6:200. doi: 10.1186/s40168-018-0578-1
10. Rainard, P, and Foucras, G. A critical appraisal of probiotics for mastitis control. Front Vet Sci. (2018) 5:251. doi: 10.3389/fvets.2018.00251
11. Khafipour, E, Krause, DO, and Plaizier, JC. A grain-based subacute ruminal acidosis challenge causes translocation of lipopolysaccharide and triggers inflammation. J Dairy Sci. (2009) 92:1060–70. doi: 10.3168/jds.2008-1389
12. Hu, X, Li, S, Mu, R, Guo, J, Zhao, C, Cao, Y, et al. The rumen microbiota contributes to the development of mastitis in dairy cows. Microbiol Spectr. (2022) 10:e0251221. doi: 10.1128/spectrum.02512-21
13. Zhao, C, Hu, X, Bao, L, Wu, K, Zhao, Y, Xiang, K, et al. Gut dysbiosis induces the development of mastitis through a reduction in host anti-inflammatory enzyme activity by endotoxemia. Microbiome. (2022) 10:205. doi: 10.1186/s40168-022-01402-z
14. Rocchetti, G, and O’Callaghan, TF. Application of metabolomics to assess milk quality and traceability. Curr Opin Food Sci. (2021) 40:168–78. doi: 10.1016/j.cofs.2021.04.005
15. Hu, X, Li, S, Fu, Y, and Zhang, N. Targeting gut microbiota as a possible therapy for mastitis. Eur J Clin Microbiol Infect Dis. (2019) 38:1409–23. doi: 10.1007/s10096-019-03549-4
16. O’Toole, PW, Marchesi, JR, and Hill, C. Next-generation probiotics: the spectrum from probiotics to live biotherapeutics. Nat Microbiol. (2017) 2:17057. doi: 10.1038/nmicrobiol.2017.57
17. Kimura, M. A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J Mol Evol. (1980) 16:111–20. doi: 10.1007/BF01731581
18. Kalmus, P, Simojoki, H, Pyörälä, S, Taponen, S, Holopainen, J, and Orro, T. Milk haptoglobin, milk amyloid a, and N-acetyl-β-D-glucosaminidase activity in bovines with naturally occurring clinical mastitis diagnosed with a quantitative PCR test. J Dairy Sci. (2013) 96:3662–70. doi: 10.3168/jds.2012-6177
19. Li, SW, Watanabe, K, Hsu, CC, Chao, SH, Yang, ZH, Lin, YJ, et al. Bacterial composition and diversity in breast milk samples from mothers living in Taiwan and mainland China. Front Microbiol. (2017) 8:965. doi: 10.3389/fmicb.2017.00965
20. Bokulich, NA, Subramanian, S, Faith, JJ, Gevers, D, Gordon, JI, Knight, R, et al. Quality-filtering vastly improves diversity estimates from Illumina amplicon sequencing. Nat Methods. (2013) 10:57–9. doi: 10.1038/nmeth.2276
21. Caporaso, JG, Kuczynski, J, Stombaugh, J, Bittinger, K, Bushman, FD, Costello, EK, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. (2010) 7:335–6. doi: 10.1038/nmeth.f.303
22. Edgar, RC. UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat Methods. (2013) 10:996–8. doi: 10.1038/nmeth.2604
23. Edgar, RC, Haas, BJ, Clemente, JC, Quince, C, and Knight, R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics. (2011) 27:2194–200. doi: 10.1093/bioinformatics/btr381
24. Haas, BJ, Gevers, D, Earl, AM, Feldgarden, M, Ward, DV, Giannoukos, G, et al. Chimeric 16S rRNA sequence formation and detection in sanger and 454-pyrosequenced PCR amplicons. Genome Res. (2011) 21:494–504. doi: 10.1101/gr.112730.110
25. Magoč, T, and Salzberg, SL. FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics. (2011) 27:2957–63. doi: 10.1093/bioinformatics/btr507
26. Quast, C, Pruesse, E, Yilmaz, P, Gerken, J, Schweer, T, Yarza, P, et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucl Acids Res. (2013) 41:D590–6. doi: 10.1093/nar/gks1219
27. Wang, Q, Garrity, GM, Tiedje, JM, and Cole, JR. Naïve Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol. (2007) 73:5261–7. doi: 10.1128/AEM.00062-07
28. Segata, N, Izard, J, Waldron, L, Gevers, D, Miropolsky, L, Garrett, WS, et al. Metagenomic biomarker discovery and explanation. Genome Biol. (2011) 12:R60. doi: 10.1186/gb-2011-12-6-r60
29. Pang, Z, Chong, J, Zhou, G, de Lima Morais, DA, Chang, L, Barrette, M, et al. MetaboAnalyst 5.0: narrowing the gap between raw spectra and functional insights. Nucl Acids Res. (2021) 49:W388–96. doi: 10.1093/nar/gkab382
30. Wang, H.T. (2004). Production and utilization of cellulytic enzyme and protease from rumen bacteria. Doctoral dissertation. Taipei, National Taiwan University
31. Menke, KH, and Steingass, H. Estimation of the energetic feed value obtained from chemical analysis and in vitro gas production using rumen fluid. Anim Res Devel. (1988) 28:7–55.
32. Xue, M, Sun, H, Wu, X, Guan, LL, and Liu, J. Assessment of rumen microbiota from a large dairy cattle cohort reveals the pan and core bacteriomes contributing to varied phenotypes. Appl Environ Microbiol. (2018) 84:e00970–18. doi: 10.1128/AEM.00970-18
33. Xin, J, Chai, Z, Zhang, C, Zhang, Q, Zhu, Y, Cao, H, et al. Comparing the microbial community in four stomach of dairy cattle, yellow cattle and three yak herds in Qinghai-Tibetan plateau. Front Microbiol. (2019) 10:1547. doi: 10.3389/fmicb.2019.01547
34. Langille, MG, Zaneveld, J, Caporaso, JG, McDonald, D, Knights, D, Reyes, JA, et al. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat Biotechnol. (2013) 31:814–21. doi: 10.1038/nbt.2676
35. Stevenson, DM, and Weimer, PJ. Dominance of prevotella and low abundance of classical ruminal bacterial species in the bovine rumen revealed by relative quantification real-time PCR. Appl Microbiol Biotechnol. (2007) 75:165–74. doi: 10.1007/s00253-006-0802-y
36. Liu, K, Zhang, Y, Huang, G, Zheng, N, Zhao, S, and Wang, J. Ruminal bacterial community is associated with the variations of total milk solid content in Holstein lactating cows. Anim Nutr. (2022) 9:175–83. doi: 10.1016/j.aninu.2021.12.005
37. da Rocha Lapa, F, da Silva, MD, de Almeida Cabrini, D, and Santos, AR. Anti-inflammatory effects of purine nucleosides, adenosine and inosine, in a mouse model of pleurisy: evidence for the role of adenosine A2 receptors. Purinergic Signal. (2012) 8:693–704. doi: 10.1007/s11302-012-9299-2
38. Hu, Y, He, Y, Gao, S, Liao, Z, Lai, T, Zhou, H, et al. The effect of a diet based on rice straw co-fermented with probiotics and enzymes versus a fresh corn Stover-based diet on the rumen bacterial community and metabolites of beef cattle. Sci Rep. (2020) 10:10721. doi: 10.1038/s41598-020-67716-w
39. Li, F, and Guan, LL. Metatranscriptomic profiling reveals linkages between the active rumen microbiome and feed efficiency in beef cattle. Appl Environ Microbiol. (2017) 83:e00061–17. doi: 10.1128/AEM.00061-17
40. Grilli, DJ, Fliegerová, K, Kopečný, J, Lama, SP, Egea, V, Sohaefer, N, et al. Analysis of the rumen bacterial diversity of goats during shift from forage to concentrate diet. Anaerobe. (2016) 42:17–26. doi: 10.1016/j.anaerobe.2016.07.002
41. Larsen, JM. The immune response to Prevotella bacteria in chronic inflammatory disease. Immunology. (2017) 151:363–74. doi: 10.1111/imm.12760
42. Abbas, Z, Sammad, A, Hu, L, Fang, H, Xu, Q, and Wang, Y. Glucose metabolism and dynamics of facilitative glucose transporters (GLUTs) under the influence of heat stress in dairy cattle. Meta. (2020) 10:312. doi: 10.3390/metabo10080312
43. Nafikov, RA, and Beitz, DC. Carbohydrate and lipid metabolism in farm animals. J Nutr. (2007) 137:702–5. doi: 10.1093/jn/137.3.702
44. Flint, HJ, Duncan, SH, Scott, KP, and Louis, P. Interactions and competition within the microbial community of the human colon: links between diet and health. Environ Microbiol. (2007) 9:1101–11. doi: 10.1111/j.1462-2920.2007.01281.x
45. Tong, J, Zhang, H, Yang, D, Zhang, Y, Xiong, B, and Jiang, L. Illumina sequencing analysis of the ruminal microbiota in high-yield and low-yield lactating dairy cows. PLoS One. (2018) 13:e0198225. doi: 10.1371/journal.pone.0198225
46. Krumholz, LR, and Bryant, MP. Syntrophococcus sucromutans sp. nov. gen. Nov. uses carbohydrates as electron donors and formate, methoxymonobenzenoids or methanobrevibacter as electron acceptor systems. Arch Microbiol. (1986) 143:313–8. doi: 10.1007/BF00412795
47. Chiquette, J, Talbot, G, Markwell, F, Nili, N, and Forster, RJ. Repeated ruminal dosing of Ruminococcus flavefaciens NJ along with a probiotic mixture in forage or concentrate-fed dairy cows: effect on ruminal fermentation, cellulolytic populations and in sacco digestibility. Can J Anim Sci. (2007) 87:237–49. doi: 10.4141/A06-066
48. Hassan, A, Gado, H, Anele, UY, Berasain, MAM, and Salem, AZM. Influence of dietary probiotic inclusion on growth performance, nutrient utilization, ruminal fermentation activities and methane production in growing lambs. Anim Biotechnol. (2020) 31:365–72. doi: 10.1080/10495398.2019.1604380
49. Wong, CB, Odamaki, T, and Xiao, J. Beneficial effects of Bifidobacterium longum subsp. longum BB536 on human health: modulation of gut microbiome as the principal action. J Functional Foods. (2019) 54:506–19. doi: 10.1016/j.jff.2019.02.002
50. Cherdthong, A, and Wanapat, M. Rumen microbes and microbial protein synthesis in Thai native beef cattle fed with feed blocks supplemented with a urea-calcium sulphate mixture. Arch Anim Nutr. (2013) 67:448–60. doi: 10.1080/1745039X.2013.857080
51. Wallace, RJ, McEwan, NR, McIntosh, FM, Teferedegne, B, and Newbold, CJ. Natural product as manipulators of rumen fermentation. Asian-Aust J Anim Sci. (2002) 15:1458–68. doi: 10.5713/ajas.2002.1458
52. Fang, S, Chen, X, Ye, X, Zhou, L, Xue, S, and Gan, Q. Effects of gut microbiome and short-chain fatty acids (SCFAs) on finishing weight of meat rabbits. Front Microbiol. (2020) 11:1835. doi: 10.3389/fmicb.2020.01835
53. Carberry, CA, Kenny, DA, Han, S, McCabe, MS, and Waters, SM. Effect of phenotypic residual feed intake and dietary forage content on the rumen microbial community of beef cattle. Appl Environ Microbiol. (2012) 78:4949–58. doi: 10.1128/AEM.07759-11
54. Izuddin, WI, Loh, TC, Samsudin, AA, Foo, HL, Humam, AM, and Shazali, N. Effects of postbiotic supplementation on growth performance, ruminal fermentation and microbial profile, blood metabolite and GHR, IGF-1 and MCT-1 gene expression in post-weaning lambs. BMC Vet Res. (2019) 15:315. doi: 10.1186/s12917-019-2064-9
55. Hungate, R.E. (1966). Ruminal functions related to rumen microbial activity, The rumen and its microbes (NY: Academic), 148–205
56. Ghorbani, GR, Morgavi, DP, Beauchemin, KA, and Leedle, JA. Effects of bacterial direct-fed microbials on ruminal fermentation, blood variables, and the microbial populations of feedlot cattle. J Anim Sci. (2002) 80:1977–85. doi: 10.2527/2002.8071977x
57. Goto, H, Qadis, AQ, Kim, YH, Ikuta, K, Ichijo, T, and Sato, S. Effects of a bacterial probiotic on ruminal pH and volatile fatty acids during subacute ruminal acidosis (SARA) in cattle. J Vet Med Sci. (2016) 78:1595–600. doi: 10.1292/jvms.16-0211
58. Hill, KJ, and Mangan, JL. The formation and distribution of methylamine in the ruminant digestive tract. Biochem J. (1964) 93:39–45. doi: 10.1042/bj0930039
59. Lagishetty, CV, and Naik, SR. Polyamines: potential anti-inflammatory agents and their possible mechanism of action. Indian J Pharmacol. (2008) 40:121–5. doi: 10.4103/0253-7613.42305
60. Bardócz, S, Grant, G, Brown, DS, Ralph, A, and Pusztai, A. Polyamines in food—implications for growth and health. J Nutr Biochem. (1993) 4:66–71. doi: 10.1016/0955-2863(93)90001-D
Keywords: Ruminococcus flavefaciens, Bifidobacterium longum subsp. longum, metabolites, biomarkers, dairy cows
Citation: Hsieh J-C, Chuang S-T, Hsu Y-T, Ho S-T, Li K-Y, Chou S-H and Chen M-J (2023) In vitro ruminal fermentation and cow-to-mouse fecal transplantations verify the inter-relationship of microbiome and metabolome biomarkers: potential to promote health in dairy cows. Front. Vet. Sci. 10:1228086. doi: 10.3389/fvets.2023.1228086
Edited by:
Wen-Chao Liu, Guangdong Ocean University, ChinaReviewed by:
Cheng Wang, Zhejiang University, ChinaJeehwan Choe, Korea National College of Agriculture and Fisheries – KNCAF, Republic of Korea
Copyright © 2023 Hsieh, Chuang, Hsu, Ho, Li, Chou and Chen. 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: Ming-Ju Chen, cmj@ntu.edu.tw
†These authors have contributed equally to this work