Hanne Honerlagen1
Henry Reyer1Dierck Segelke2Carolin Beatrix Maria Müller3Marie Christin Prahl3
Siriluck Ponsuksili1
Nares Trakooljul1
Norbert Reinsch4Björn Kuhla3
Klaus Wimmers1,5*- 1Research Institute for Farm Animal Biology (FBN), Institute of Genome Biology, Dummerstorf, Germany
- 2IT-Solutions for Animal Production, Vereinigte Informationssysteme Tierhaltung w.V. (vit), Verden, Germany
- 3Research Institute for Farm Animal Biology (FBN), Institute of Nutritional Physiology “Oskar Kellner”, Dummerstorf, Germany
- 4Research Institute for Farm Animal Biology (FBN), Institute of Genetics and Biometry, Dummerstorf, Germany
- 5Faculty of Agricultural and Environmental Sciences, University of Rostock, Rostock, Germany
Efforts to reduce nitrogen (N) emissions are currently based on the optimization of dietary- N supply at average herd N requirements. The implementation of the considerable individual differences and predispositions in N- use efficiency and N- excretion in breeding programs is hampered by the difficulty of data collection. Cow individual milk urea (MU) concentration has been proposed as an easy-to-measure surrogate trait, but recent studies questioned its predictive power. Therefore, a deeper understanding of the biological mechanisms underlying predisposed higher (HMUg) or lower (LMUg) MU concentration in dairy cows is needed. Considering the complex N- metabolism in ruminants, the distinction between HMUg and LMUg could be based on differences in (i) the rumen microbial community, (ii) the host-specific transcription processes in the rumen villi, and (iii) the host–microbe interaction in the rumen. Therefore, rumen fluid and rumen epithelial samples from 10 HMUg and 10 LMUg cows were analyzed by 16S sequencing and HiSeq sequencing. In addition, the effect of dietary-N reduction on ruminal shifts was investigated in a second step. In total, 10 differentially abundant genera (DAG) were identified between HMUg and LMUg cows, elucidating greater abundances of ureolytic Succinivibrionaceae_UCG-002 and Ruminococcaceae_unclassified in LMUg animals and enhanced occurrences of Butyvibrio in HMUg cows. Differential expression analysis revealed genes of the bovine Major Histocompatibility Complex (BOLA genes) as well as MX1, ISG15, and PRSS2 displaying candidates of MU predisposition that further attributed to enhanced immune system activities in LMUg cows. A number of significant correlations between microbial genera and host transcript abundances were uncovered, including strikingly positive correlations of BOLA-DRA transcripts with Roseburia and Lachnospiraceae family abundances that might constitute particularly prominent microbial–host interplays of MU predisposition. The reduction of feed-N was followed by 18 DAG in HMUg and 19 DAG in LMUg, depicting pronounced interest on Shuttleworthia, which displayed controversial adaption in HMUg and LMUg cows. Lowering feed-N further elicited massive downregulation of immune response and energy metabolism pathways in LMUg. Considering breeding selection strategies, this study attributed information content to MU about predisposed ruminal N-utilization in Holstein–Friesians.
Introduction
Nitrogen (N) excretions of dairy cows contribute substantially to N-deposition in the environment (Castillo et al., 2000). Considering the relevance for climate change (nitrous oxide (N2O) atmospheric emissions) and groundwater quality (nitrate () contamination), the dairy sector is demanded to reduce its N-emissions in the near future (Uwizeye et al., 2020). The primary source of N-emissions from dairy cows is represented by unutilized N, which goes into the cow via feed crude protein (CP) and is excreted via urine, milk and feces (Bryant, 1970; Abdoun et al., 2006). In this context, the ruminal conversion rate of feed N into valuable microbial protein, which is the result of rumen microbial growth, mainly determines the number of N-emissions (Tan et al., 2021).
In detail, the cow acquires N via the dietary CP. A small part of CP passes through the rumen as undegradable protein (UDP) and is absorbed in the small intestine or excreted via feces. The majority of dietary N from CP is provided to the rumen microbes for microbial growth. However, the rate of microbial N-utilization depends on various factors such as ruminal pH, dietary energy content, the absorption capacity of the rumen epithelium and the microbial community. Thus, a part of the ruminal N-influx is incorporated into microbial protein, whereas the remaining part is microbially converted to non-protein-N (NPN) in forms of ammonia (NH3), ammonium-ions (), peptides and amino acids (Bryant, 1970; Hartinger et al., 2018; Cholewińska et al., 2020). The cow absorbs amino acids and peptides deriving from the microbial protein in the small intestine whereas the NPN-compounds are absorbed by the rumen epithelium. The amino acids and peptides can be utilized by the cow as valuable N-sources, while absorbed NH3 and is transported with the blood stream to the liver. There, the ammonia is detoxified, yielding urea in an energy-intensive process. Urea is secreted back into the blood cycle from which it is eliminated via urine (urinary urea, UU) and milk (milk urea, MU), or secreted into the rumen either across the rumen epithelium or via saliva (Burgos et al., 2007; Spek et al., 2012). The recycled ruminal urea is converted by epithelium-bound, ureolytic rumen bacteria to CO2 and NH3. The latter, in turn, can be utilized by NH3-assimilating microbes for the metabolism of microbial protein (Jin et al., 2018).
In symbiosis with their microbes, dairy cows may utilize on average 25% of the dietary N (Calsamiglia et al., 2010). However, substantial variation in the N-utilization efficiency, which is defined as the ratio of gram N in product to gram N intake between individual cows has been observed, ranging from 15 to 40% (Calsamiglia et al., 2010; Spek et al., 2012; Gardiner et al., 2016; Chen et al., 2021). In recent years various studies have focused on the reduction of dietary CP levels to minimize N-emissions at herd level, but neither the cow-specific potential has been exploited intensely nor has the predisposition for individual N-use efficiency been considered for breeding progress yet (Cantalapiedra-Hijar et al., 2020; Bergen, 2021; Jahnel et al., 2021; Te Pas et al., 2021). One reason is that data collection for the cow-individual N-use efficiency is not feasible in huge cow cohorts. Generally, neither the cow individual N-uptake nor the content of N in urine and feces is measurable in practice operating dairy farms. Thus, in the last years the cow-specific MU concentration has been considered as a substitute trait (Burgos et al., 2007, 2010; König et al., 2008; Honerlagen et al., 2021; Jahnel et al., 2021). MU values are available to breeding companies at a large scale due to the results of monthly milk records. Furthermore, MU is moderately heritable and only weak or even no genetic correlations to further milk performance parameters (i.e., milk protein yield, milk fat yield) were evidenced (Wood et al., 2003; Miglior et al., 2007), which would, in principle, support MU as a trait for breeding selection. Due to the phenotypic correlation between MU, UU and blood urea concentration (Spek et al., 2012; Guliński et al., 2016) it has been hypothesized that breeding selection on lower MU concentration could lead to lower N-urine concentration and thus to lower N-emissions from dairy farming (Beatson et al., 2019). In addition, lower MU concentrations were also suggested as biomarkers for increased N-utilization efficiency (Nousiainen et al., 2004; Lavery and Ferris, 2021), being predominantly explained by enhanced ruminal N-utilization (Guliński et al., 2016). However, up to now the validity of MU as a substitute trait for cow-individual N-use efficiency and N-excretion is relying on the phenotypic correlation between MU and UU, but its utility has been questioned by the results of two recent studies (Correa-Luna et al., 2021; Müller et al., 2021). Thus, it became evident, that a deeper knowledge of the biological mechanisms underlying predisposed MU concentrations is indispensable before MU can be seriously proposed for breeding selection strategies.
Considering the large impact of ruminal processes on the N-metabolism in dairy cows, this study aims to elucidate ruminal mechanisms that distinguish between cows with predisposed higher (HMUg) and lower (LMUg) MU concentration to contribute to the validation of MU as substitutive trait for N-utilization efficiency and N-excretion. In detail, this study focuses on differences in (i) the rumen microbial community, (ii) the host-specific transcriptional processes in the rumen villi and (iii) the host–microbe interaction in the rumen between HMUg and LMUg cows.
Moreover, assumed that LMUg cows may possess lower N-input requirements as a consequence of more efficient N-utilization, which would result in an N-oversupply if fed a ration according to the recommendations, this study further investigated the effect of reduced dietary N-supply on the microbial community and on host specific rumen tissue processes in HMUg and LMUg cows.
Materials and methods
Animal husbandry and sampling were carried out according to the guidelines of the German Animal Protection Law. All protocols were approved by the Institute's Animal Welfare Commission. The experimental protocol is in strict compliance with the German Animal Welfare Legislation, and has been approved by the Ethics Committee of the Federal State of Mecklenburg, Western Pomerania, Germany (Landesamt für Landwirtschaft, Lebensmittelsicherheit und Fischerei; LALLF M-V7221.3-2-019/19) and is in accordance with the ARRIVE guidelines.
Cow population and sample collection
The predisposition of cows for a higher or lower MU concentration (HMUg and LMUg) was determined by deregressed proofs of an estimated breeding value for MU (EBVMU). The EBVMUs were estimated by vit Verden (Vereinigte Informationssysteme Tierhaltung, Verden, Germany) applying the established model for somatic cell score in milk [Vereinigte Informationssysteme Tierhaltung (VIT)., (2020, access 20.01.2020)]. The estimation was based on a milk record dataset of around 8 million lactation curves of milk urea values from Holstein cows in their first to third parities. The deregression of EBVMUs is a pre step for genomic breeding value estimation. The HMUg was determined as predisposition by positive values, LMUg by negative values. Based on this categorization, 10 HMUg and 10 LMUg cows were housed in pairs, comprising one HMUg and one LMUg cow, respectiveley. The phenotypes of MU concentration were obtained from the recent five milk records depicting HMUg = 278.1 ± 15.8 (mean ± SD) mg/L and LMUg = 181.54 ± 15.08 (mean ± SD) mg/L with comparable milk yields between the groups.
These cow groups were fed an isocaloric diet [10.1 ± 0.2 MJ metabolisable energy/kg dry matter (DM)], with either recommended (normal) CP content (NP with 157 ± 2.68 (mean ± SD) g CP/kg DM) or low CP content [LP with 139 ± 8.37 (mean ± SD) g CP/kg DM], resulting in a 4 ×5 balanced design (NPxHMUg, NPxLMUg, LPxHMUg, LPxLMUg; each n = 5). The experiment was conducted in 10 trials, with each trial comprising one cow pair (HMUg and LMUg), which was fed either the NP or the LP diet. The LP and NP rations were offered twice-daily for ad libitum intake as a total mixed ration. A detailed description of the diet, the feeding scheme as well as phenotypic analyses of N-metabolites (i.e., ammonia concentration in the rumen fluids, plasma urea, and glutamine concentration) is documented in Müller et al. (2021). After 2 weeks, the cows were slaughtered 4 h after morning milking and feeding by exsanguination after captive bolt stunning. Samples of rumen villi from the saccus dorsalis containing all rumen villus cell layers were collected, i.e., the basal membrane cells, which connect villi to blood stream and muscle tissue. The samples were rinsed, snap frozen in liquid nitrogen and stored at −80°C until RNA isolation. In addition, the rumen fluids (50 ml) were collected as a homogeneous mixture of the entire rumen content and stored at −20° until isolation of microbial DNA.
Rumen fluids: Microbial DNA extraction, 16S rRNA amplicon sequencing, and data preparation
The rumen fluid samples (50 ml) were mixed and microbial DNA was extracted with the PowerLyzer PowerSoil DNA isolation kit (QIAGEN, Hilden, Germany) following manufacture's recommendations, but utilizing 800 μl rumen fluid sample (instead of 250 μl) and two additional heat incubation steps (step 1: 70°C, 10 min; step 2: 95°C, 10 min) before the bead beating procedure. The quantity of the extracted DNA was determined on NanoDrop ND-2000 (Thermo Fisher Scientific, Dreieich, Germany) and afterwards diluted to 10 ng/μl per sample. The diluted DNA extracts were utilized to amplify the V4 hyper-variable region of the 16S rRNA gene. The specific primer enclosed the V4-specific sequence, the sequencing flow cell adapter and a specific barcode for each primer as previously described (Kozich et al., 2013). The amplification was conducted in duplicate by polymerase chain reaction (PCR) with initial denaturation at 95°C for 2 min, followed by 35 cycles at 95°C for 30 s, 50°C for 60 s, and 72°C for 90 s, and a final extension for 10 min at 72°C using the 5 Prime HotMasterMix (5 Prime, Hamburg, Germany). The PCR products (16S amplicons) were checked on agarose gel and subsequently purified, normalized and pooled for sequencing using the SequalPrep normalization plate kit (Thermo Fisher Scientific, Dreieich, Germany) The sequencing of the amplicons was conducted on HiSeq2500 (Illumina, San Diego, CA) with 250 bp paired-end reads.
The sequence data were then trimmed for the adapter sequences and filtered with mothur software for sequences that exceeded 275 bp in length, which contained ambiguous base calls or excessively long homopolymers (>8 bases), as well as putative chimeras were removed [version 1.44.1 (Schloss et al., 2009)].The sequences were further globally aligned to the Silva reference database (release 138; https://www.arb-silva.de/, access 10.6.20) and subsequently aggregated into operational taxonomic units (OTU) considering a sequence identity of ≥97%. Taxonomic annotations for OTUs were obtained from the Silva database (release 138). Considering stratification by dispersion, the count data was rarefied to the sample with the lowest read depth (745,235 reads) as suggested by Mcmurdie and Holmes (2014).
Rumen fluids: Statistical analysis of 16S rRNA sequencing data
Initially, very low abundant OTUs with <10 counts in all 20 cows were pre-discarded from further analyses. The microbial community in rumen fluids was analyzed with inverse Simpson diversity index (alpha diversity) at OTU level utilizing the R packages “vegan” (Oksanen et al., 2013) and “agricolae” (De Mendiburu Delgado, 2009). Means of Simpson indices of the predisposition groups (HMUg and LMUg) and the diet groups (NP and LP in HMUg and LMUg cows, respectively) were tested for significant differences by applying Student's t-test. The subsequent analyses were conducted at genus level, excluding genera with <10 counts in 5 or more cows to avoid potential noise by barely observed reads. The remaining dataset was tested for significantly differentially abundant genera (DAG) between the groups (HMUg–LMUg; NPxHMUg–LPxHMUg; NPxLMUg–LPxLMUg) utilizing Wald test in the DESeq2 R package (Love et al., 2014). For the contrast of predisposition (HMUg–LMUg) the trial was implemented as further fixed effect in the statistical model, accounting for different CP levels in the feed ration. For the differences between the diets, the contrasts within each predisposition group were analyzed and expressed as fold changes obtained from DESeq2 delogarithmized log2fold changes. Genera with p-values below 0.05 were considered as DAG.
Rumen tissue: RNA extraction, library preparation, sequencing, and data preparation
The total RNA from homogenized rumen villi was isolated using TRI reagent (Sigma–Aldrich, Taufkirchen, Germany), followed by DNaseI treatment (Roche, Mannheim, Germany) and a subsequent column-based purification step utilizing the NucleoSpin RNA Kit (Macherey-Nagel, Düren, Germany). The RNA quantity was measured on NanoDrop ND-2000 (Thermo Fisher Scientific) and further assessed for integrity by application on agarose gel and via 2100 Bioanalyzer system (Agilent, Santa Clara, CA). The RNA library preparation was performed using the TruSeq Stranded mRNA kit according to the manufacture's recommendation (Illumina). The resulting libraries were quality checked on the 2100 Bioanalyzer system and quantified using the Invitrogen Qubit dsDNA HS kit (Thermo Fisher Scientific). Libraries were paired-end sequenced for 2 ×71 cycles on HiSeq 2500 (Illumina). The raw sequencing reads (fastq) were quality-checked utilizing FastQC, version 0.11.8 (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/). The data preprocessing was conducted with Trim Galore v.0.6.5 (https://www.bioinformatics.babraham.ac.uk/projects/trim_galore/). The low-quality reads (mean Q-score <20) and short reads (<20 bp) were filtered out. Adapter-like sequences at the 3'-end of sequence reads were trimmed. The remaining high quality paired-end reads comprised a read depth of 20,324,785 ± 3,234,163 (mean ± SD) per sample and were further aligned to the current reference genome, Bos_taurus.ARS-UCD1.2 (access 01. 08.2020) with an average alignment rate of 97.86% ± 0.01% (SD) using Hisat2, version 2.2.0 (Kim et al., 2015, 2019; Pertea et al., 2016). The number of reads uniquely mapped to each gene were extracted from the HISAT2 mapping results using HTSeq, version 0.12.4 (Anders et al., 2015).
Rumen tissue: Statistical analysis of RNAseq data
The RNA dataset was analyzed for significantly differentially expressed genes (DEG) between groups—corresponding to the 16S microbiota sequence data (HMUg–LMUg; NPxHMUg–LPxHMUg; and NPxLMUg–LPxLMUg)—by utilizing Wald test in DESeq2 package (Love et al., 2014). Genes with <10 counts in 10 or more cows were excluded from analysis, leading to a dataset of 12,627 genes. The trial was considered as further fixed effect for HMUg–LMUg contrast. Genes with a Benjamini–Hochberg adjusted p-values <0.05 were considered as DEG. The DEG identifiers were extracted to investigate the biological function of the respective gene using the databases GeneCards (http://www.genecards.org), Ensembl (http://www.ensembl.org), and the published literature. Furthermore, genes revealing p-values below 0.01 (corresponding adjusted p-values <0.6) were imported into Ingenuity Pathway Analysis (IPA; Ingenuity® Systems, http://www.ingenuity.com) for enrichment analysis between the groups (HMUg–LMUg; NPxHMUg–LPxHMUg; NPxLMUg–LPxLMUg). Pathways with an adjusted p-value <0.05 were defined as significantly enriched between the groups. Significant pathway inhibition or activation was defined by Z-scores (significant inhibition: Z-score < -2; significant activation: Z-score >2) in IPA. The Z-scores between −2 and 0 or 0 and 2 were considered indicators of inhibition or activation, respectively. The pathways were further assigned to general functions after evaluation in Targeted Explorer (IPA). The pathways of cancer and disease were excluded.
Integration of microbial and transcriptome data: Statistical analysis
Datasets of rumen microbial genera (n = 229 microbes) and gene expression (n = 12,627 genes) were initially transformed by variance-stabilizing transformation (DESeq2 R package). Subsequently, sparse Partial Least Squares Discriminant Analysis (sPLS-DA) was applied on both datasets. The sPLS-DA is incorporated in R package “mixOmics” (version 6.6.2, available at: http://mixomics.org) (Rohart et al., 2017) and achieves feature selection and dimension reduction simultaneously in dependence on a specific trait (Chun and Keleş, 2010). In this case, sPLS-DA was utilized to select a key subset of microbial genera and a key subset of genes that discriminate the microbial community and the rumen tissue transcriptome between HMUg and LMUg cows on two components each. The discrimination by the generated subsets was visualized in two plots.
The data subsets were further analyzed for microbial abundance–gene transcript correlations by utilizing Pearson statistics. Microbe–gene pairs exhibiting p-values below 0.01 were considered as significantly correlated. The results were visualized in a correlation map using R package “pheatmap” (version 1.0.12, available at: https://cran.r-project.org/web/packages/pheatmap) (Kolde and Kolde, 2015). The correlation map identified two dominant clusters (I and II). The genes attributing to the respective clusters were subsequently analyzed for significantly enriched pathways between HMUg and LMUg predisposition in IPA (adjusted p <0.05).
Results
The microbial community in rumen fluids comprised 7,697 OTUs after filtering, which attributed to 24 phyla, 45 orders, 101 classes, 168 families, and 350 genera across all 20 cows. The filtered transcriptome dataset constituted 12,767 genes.
Microbial analyses
The microbial composition ranked at inverse Simpson values of 87.51 ± 39.13 (mean ± SD) in HMUg cows and 80.95 ± 34.24 (mean ± SD) in LMUg cows, indicating no significant differences between the predispositions (Figure 1A). With lower dietary CP level, alpha diversity numerically decreased, but not significantly, neither in HMUg nor in LMUg cows (Figures 1B,C).
Figure 1. Alpha diversity of the microbial community in rumen fluids of (A) HMUg and LMUg cows and after lowering dietary N in HMUg (B) and LMUg (C) cows, respectively. Groups are displayed on the x-axis; inverse Simpson values scale the y-axis.
The DESeq2 analysis at genus level revealed 10 DAG between HMUg and LMUg cows, which covered relative abundances between 0.0016 and 2.3668% across all cows (Table 1). Notably, the six most abundant genera (Prevotella, Rikenellaceae_RC9_gut_group, Methanobrevibacter, Christensenellaceae_R-7_group, F082_ge, and NK4A214_group) captured together 51.6% of total rumen abundances and did not differ significantly between HMUg and LMUg cows. Within the DAG, Unclassified Ruminococcae exhibited the broadest relative abundance (2.3668%) and manifested a 1.39-fold higher presence in LMUg compared to HMUg animals, whereas Butyrivibrio was significantly more abundant in HMUg compared to LMUg cows. Moreover, Lachnospiraceae_UCG-010 showed significantly higher occurrence in HMUg compared to LMUg, while Succinivibrionaceae_UCG-002 displayed a 2.24-fold higher occurrence in LMUg than in HMUg animals. However, highest fold changes were observed for Acetobacter and Monoglobus, with Acetobacter depicting 3.25-fold more abundance in HMUg and Monoglobus being 2.33-times more prevalent in LMUg animals.
Table 1. Significantly differentially abundant genera (DAG) in the rumen of HMUg and LMUg cows.
For the comparison of NP and LP diets, 18 DAG were identified in HMUg and 19 DAG in LMUg cows, including four DAG that were previously identified in HMUg–LMUg comparison (Table 2). The reduction of dietary CP content was generally accompanied by high fold changes of genera abundances. For instance, MVP-15_ge abundances decreased 133-fold in HMUg, whereas Sharpea abundances increased 24-fold in HMUg and 166-fold in LMUg with LP feeding. Incertae sedis was less abundant in both, HMUg and LMUg animals, when cows were exposed to the LP diet. Furthermore, the abundances of Shuttleworthia increased in HMUg, but decreased in LMUg cows with the LP diet. Following dietary CP reduction, HMUg animals exhibited a further increase of Acidaminococcus abundances and of three genera belonging to phylum Desulfobacterota (Desulfovibrio, Desulfobulbus, Desulfobulbaceae_unclassified) with mean relative abundances between 0.0015 and 0.0993%. In LMUg animals, the significant increases of Streptococcus and Lachnospiraceae_NK3A20_group abundances were observed on the LP compared to NP ration. On the contrary, Alphaproteobacteria_unclassified and Succinivibrionaceae_unclassified abundances decreased remarkably under conditions of LP feeding in LMUg cows at 4.05- and 13.3-times, respectively. Interestingly, Acetobacter, whose abundance was significantly lower in LMUg than in HMUg cows (Table 1), showed a further 21-fold decrease in LMUg cows upon the reduction of the dietary CP content (Table 2).
Table 2. Significantly differentially abundant genera (DAG) in the rumen comparing normal crude protein (NP) and low crude protein (LP) diet offered to HMUg and LMUg cows.
Transcriptome analyses
In the filtered transcriptome dataset, 28 genes were identified as differentially expressed in the rumen between HMUg and LMUg cows (Supplementary material 1, p. 1). The highest fold changes were observed for Matrix Metallopeptidase 3 (MMP3) and Transition Protein 2 (TNP2), with higher transcript abundances in HMUg and LMUg, respectively. Furthermore, 166 genes passed the relaxed significance threshold (p <0.01) for enrichment analysis (Supplementary material 1, p. 2) and revealed 13 significantly enriched pathways that contributed almost exclusively (9 out of 13) to the immune system (Table 3). Antigen Presentation Pathway revealed the most prominent adjusted p-value. In addition, Interferon Signaling and Neuroinflammation Signaling pointed as two further candidates, displaying significant inhibition in HMUg. Whereas, Ubiquitin Like Modifier 15 (ISG15) and MX Dynamin Like GTPase 1 (MX1) were included in the former, MMP3 and Prostaglandin-Endoperoxide Synthase 2 (PTGS2) contributed to the latter. Interestingly, these four genes also ranked under the top 25 genes in terms of absolute fold changes among all annotated genes in the filtered dataset (Supplementary material 1). However, Serine Protease 2 (PRSS2) and UL16 Binding Protein 17 (ULBP17) displayed the most prominent fold changes.
Table 3. Significantly enriched pathways in rumen tissue comparing HMUg and LMUg cows.
Considering the reduction of dietary CP in HMUg cows, no DEGs were observed, but 68 genes passed the relaxed significance threshold for pathway analysis. This list culminated in the identification of Intrinsic Prothrombin Activation Pathway and MSP-RON Signaling Pathway as significantly enriched (Table 4; Supplementary material 2, p. 1). In contrast, LMUg animals showed 1,220 DEGs and 139 significantly enriched pathways when dietary CP content was reduced (Table 4; Supplementary material 2, p. 2, 3). Whereas, Oxidative Phosphorylation and Semaphorin Neuronal Repulsive Signaling Pathway were significantly activated when LMUg cows were fed the LP ration, additional 79 pathways were identified to be significantly inhibited, mainly representing pathways of immune response and energy metabolism as well as pathways that contribute to cell cycle and cell maintenance (Table 4).
Table 4. Significantly enriched pathways between normal protein (NP) and low protein (LP) diet in rumen tissue of HMUg and LMUg cows.
Integration of microbial and transcriptome data
The multivariate discriminant analysis (sPLS-DA) revealed a microbial subset comprising 50 genera in the first component and 5 genera in the second component separating HMUg and LMUg cows' rumen microbial communities (Figure 2A; Supplementary material 3, p. 1). Two genera contributed to both components, giving a subset of 53 genera that distinguished HMUg and LMUg cows. Eight genera of the subset were also identified as DAG in HMUg–LMUg contrast. The transcriptome subset included 100 genes, comprising 50 genes in each component (Figure 2B; Supplementary material 3, p. 2). Twelve genes thereof were uncovered as DEG between HMUg and LMUg.
Figure 2. Sparse partial least square-discriminant analysis of HMUg and LMUg cows, driven by (A) 53 key microbial genera in the rumen fluid and by (B) 100 key genes in the rumen tissue.
The correlation analysis between the microbial and the gene subset separating HMUg and LMUg cows detected 157 significantly correlated microbial genera–gene pairs. The two pronounced clusters I and II were generated by two groups of microbes (A and B) and two groups of genes (I and II) that correlated contrarily (Figure 3; Supplementary material 3). In detail, Cluster I depicted gene Group I (27 genes) being positively correlated to microbial Group A (24 genera) and negatively correlated to microbial Group B (29 genera). Cluster II on the other hand showed negative correlations between gene Group II (73 genes) and microbial Group A, but positive correlations between gene Group II and microbial Group B.
Figure 3. Correlation map of ruminal host-microbiota interaction depicting the relationship of microbial subset abundances and gene subset transcripts that separate HMUg and LMUg rumen profiles. Subsets derived from sPLS-DA analysis; estimated correlations are based on the microbial abundances and gene transcripts in HMUg and LMUg cows; *indicates significant correlation (p <0.01).
The subsequent pathway analysis of gene Groups I and II revealed significant enrichment of Urea Cycle, Arginine Degradation and Citrulline Metabolism for Group I, which was mainly determined by Arginase 1 (ARG1) expression (Table 5). The ARG1 expression was negatively correlated with UCG-005 abundances (r = −0.64, p <0.01; Figure 3; Supplementary material 4). Furthermore, the differential expression of Glutamine-fructose-6-phosphate transaminase 1 (GFPT1) in Group I contributed to the significant enrichment of UDP-N-acetyl-D-glucosamine biosynthesis II (Table 5). Moreover, GFPT1 was negatively correlated with UCG-010_ge (r = −0.68, p <0.01), Lachnospiraceae_UCG-006 (r = −0.56, p = 0.011), and Roseburia abundances (r = −0.56, p = 0.01) (Figure 3; Supplementary material 4). A higher expression level of GFTP1 was accompanied by increased abundances of Lachnospiraceae_UCG-010 (r = 0.52, p = 0.02), Acetobacter (r = 0.56, p = 0.01), and Shuttleworthia (r = 0.55, p = 0.01), which had been pre-identified as DAG in the comparison of HMUg and LMUg cows (Figure 3; Table 1; Supplementary material 4).
Table 5. Significantly enriched pathways between HMUg and LMUg cows based on gene Groups I and II derived from ruminal host-microbiota interaction analysis (Figure 3).
Gene Group II revealed significant enrichment of four pathways that correspond to the immune system, including B Cell Development, IL-4 Signaling and Antigen Presentation Pathway, in overlap to the analysis of ruminal DEGs (Tables 3, 5). Particularly, genes of the Major Histocompatibility Complex (BOLA-B, BOLA-DOB, and BOLA-DRA) as well as Interferon Regulatory Factor 4 (IRF4) contributed to the significant enrichment of immune pathways in Group II (Table 5). BOLA-DRA transcripts were significantly less abundant in HMUg compared to LMUg animals (DEG) and revealed positive correlations with abundances of Roseburia (r = 0.53, p = 0.02) and Lachnospiraceae_FE2018_group (r = 0.52, p = 0.02) (Figure 3; Supplementary materials 1, 4). In addition, the presence of Gamma-Glutamyltransferase 7 (GGT7) and Leukotriene C4 Synthase (LTC4S) in Group II caused the significant enrichment of Leukotriene Biosynthesis Pathway (Table 5). The transcript abundance of both genes was highly correlated with the abundance of Prevotellaceae_UCG-004 (GGT7: r = −0.67, p <0.01, LTC4S: r = −0.52, p = 0.02) (Figure 3; Supplementary material 4). Moreover, GGT7 transcripts showed further a prominent negative correlation to Lachnospiraceae_UCG-010 abundances (r = −0.05, p = 0.02), whereas LTC4S transcripts negatively correlated with the occurrence of Coxiella (r = −0.58, p <0.01).
Furthermore, UCG-010_ge genus attracted interest due to its particularly high correlations with transcript abundances of KH Domain Containing 4 (KHDC4; r = −0.72, p = 0.0003), Ribosomal RNA Processing 9 (RRP9; r = 0.71, p <0.0005), Lipase Maturation Factor 1 (LMF1; r = 0.73, p <0.0003) and DNA Polymerase Delta 4 (POLD4; r = 0.73, p = 0.0002). Furthermore, POLD4 expression showed a strikingly high negative correlation with the abundance of Shuttleworthia (r = −0.74, p = 0.00017), whereas KHDC4 transcripts showed a prominent correlation with Acetobacter occurrences (r = 0.69, p = 0.0008).
Discussion
In regard to the underlying biological principles of N-metabolism in ruminants, it was hypothesized that the predisposed variance in MU concentration is—in complement to inherent post-ruminal processes—mainly attributed to the ruminal mechanisms that differentiate between HMUg and LMUg cows and thus lead to N-utilization and N-excretion variance in genetically divergent cow cohorts. Specifically, the rumen microbial N-utilization, the absorption and diffusion processes of the rumen epithelium and the ruminal interaction between microbes and host epithelium were thought to characterize the ruminal make up of HMUg and LMUg cows.
Microbial analyses
Although no differences in the overall microbial diversity have been uncovered between HMUg and LMUg cows, some DAGs were identified. These DAG might (i) influence the amount of ruminal NPN in forms of NH3 or , (ii) determine the utilization of blood urea by ureolytic activity, or (iii) influence the permeability for the transport of NH3, and urea across the rumen epithelium.
Specifically, LMUg animals displayed a significantly higher abundance of Succinivibrionaceae_UCG-002, which is an ureolytic genus in the rumen. Ureolytic bacteria metabolize urea to CO2 and NH3. The latter can be further utilized by other ruminal genera for microbial growth, which subsequently reduces the ruminal NPN-pool (Rosendahl, 2015). Moreover, the conversion of urea into NH3 increases the urea concentration gradient between blood and the rumen lumen. This massively alters the urea transport rate, reduces blood urea concentration and promotes N-fixation into microbial protein (Rosendahl, 2015; Jin et al., 2016). Interestingly, significantly lower blood urea concentrations and lower blood urea pool sizes were observed in cows displaying the LMU compared to the HMU phenotype in a recent study (Müller et al., 2021). Furthermore, unclassified genera of the Succinivibrionaceae family have already been suspected as being more abundant in ruminants with increased N-use efficiency (Jin et al., 2016).
LMUg animals additionally exhibited significantly higher abundances of unclassified Ruminococcaceae. Various Ruminococcaceae genera correlated with the N-metabolism in goats, including Ruminococcus_2 displaying significantly higher abundances in animals with higher N utilization efficiency (Wang et al., 2019). Moreover, the Ruminococcaceae family has been suggested to play a fundamental role in amino acid and protein metabolism in the bovine rumen, and the Ruminococcaceae NK4A214 group has been attributed a role to N-recycling in cattle (Pacífico et al., 2021). Hu et al. (2019) found a positive correlation between the Ruminococcus abundances and the gap junctions in the rumen epithelium of yaks and postulated an absorptive response of the rumen epithelium to Ruminococcus metabolites. Some Ruminococcaceae genera were also reported to have urease activity (Patra and Aschenbach, 2018), and thus it might be conceivable that unclassified Ruminococcaceae increases the urea concentration gradient across the rumen wall and support the diffusion of urea in LMUg animals.
Furthermore, HMUg animals showed higher abundances of Butyvibrio, which is one of the most relevant butyrate producers in the rumen (Meehan and Beiko, 2014; Henderson et al., 2015; Kong et al., 2020). Interestingly, two Butyvibrio species, namely, B. fibrisolvens and B. proteoclasticus, were postulated as key species involved in protein degradation in the bovine rumen (Wallace and Brammall, 1985; Attwood and Reilly, 1995; Henderson et al., 2015). Due to their high proteolytic activities, they belong to the hyper-ammonia producing bacteria species (HAB) (Hartinger et al., 2018). Also, HAB species are known to negatively influence the N-efficiency of ruminants, since they cause steep increases of NH3, which cannot be directly utilized by the rumen microbes and thus leads to enhanced NH3 efflux into the bloodstream (Bento et al., 2015; Hartinger et al., 2018). Based on phenotypic categorization into HMU and LMU cow groups, HMU cows were indeed shown to exhibit numerically higher ruminal NH3 concentrations than LMU cows (Müller et al., 2021). Furthermore, higher Butyrivibrio abundances were found to be correlated with lower N-recycling efficiency in beef cattle (Alves et al., 2020).
The influence of different CP levels on the rumen microbial community has been broadly evidenced (Lapierre et al., 2005; Aguiar et al., 2014; Patra and Aschenbach, 2018; Müller et al., 2021). In general, it is known that lower CP levels in the diets of dairy cows reduce rumen fluid NH3 and , enhance the relative blood urea reflux and faster N-assimilation and N-fixation into microbial protein.
In both HMUg and LMUg animals, the reduction of dietary CP content induced significant changes in microbial abundances, which might be caused by microbial N-use competition and the general dietary nutrient availability. In this context, Sharpea abundances increased simultaneously with the decline of Incertae_sedis abundances in both, HMUg and LMUg cows, when feed-CP content was lowered. Sharpea is known to digest short-chain carbohydrates (Dias et al., 2017; Trabi et al., 2020), while Incertae_sedis (aggregated genus of the Ethanoligenenaceae and Ruminococcaceae families) primarily converts cellulose (Suen et al., 2011; Bailoni et al., 2021). The metabolic processes and the resulting microbial growth of cellulose digesters are generally slower than the growth of microorganisms, which metabolize the short-chain sugars. However, both microorganisms need to assimilate N components for their microbial growth (Russell et al., 2009). It might be conceivable that Sharpea depressed the growth of Incertae_sedis due to comparatively faster N-assimilation, when CP was lowered. Shuttleworthia attracted further interest since its abundance increased in HMUg, whereas significant abundance decreases were identified in the LMUg group. Several studies reported feed induced abundance alterations of Shuttleworthia, which were related to the dietary starch–fiber ratio and the energy level of the diet (Plaizier et al., 2017; Kotz et al., 2020; Tun et al., 2020). Furthermore, Zhang et al. (2019) postulated significantly different abundances between high and low yielding dairy breeds and identified Shuttleworthia as a genus displaying high relative abundances, which is in common with the present study. These findings promote Shuttleworthia as highly vulnerable and strongly adaptable genus, which might be of importance in characterizing the different ruminal adaption mechanisms in LMUg and HMUg on low CP levels.
Moreover, the abundances of Acidaminococcus and three genera of Desulfobacterota phylum (Desulfovibrio, Desulfobulbus, and Desulfobulbaceae_unclassified) increased in HMUg cows with the reduction of dietary CP content. Interestingly, several species of these genera have been previously categorized as HAB (Attwood and Reilly, 1995; Eschenlauer et al., 2002; Loubinoux et al., 2002). The growth of HAB is known to be suppressed by high NH3 and concentrations in a self-regulating manner (Sales et al., 2000). Interestingly, phenotypic HMU animals displayed significant decreases in previously high ruminal NH3 and concentrations following CP reduction, whereas in LMU phenotypes only moderate ruminal NH3 and declines were observed (Müller et al., 2021). Since genera abundances of the Desulfobacterota phylum appeared to be significantly increased exclusively in HMUg, it is conceivable that their growth was inhibited by high ruminal NH3 and levels in HMUg under NP diet feeding, and that the reduction of dietary N might have had a positive effect on their growth.
In LMUg cows, dietary CP reduction led to increased abundances of Streptococcus, which is one of the most important proteolytic members in the rumen fluid (Attwood and Reilly, 1995; Ling and Armstead, 1995). Streptococcus utilizes rapidly assimilable energy sources (starch), is known to possess ureolytic activity and is highly adaptable to various N-sources, which could provide a competitive advantage under N scarcity (Jin et al., 2016; Hartinger et al., 2018; Patra and Aschenbach, 2018).
Under conditions of LP diet feeding, LMUg cows further exhibited a significant abundance decline of Acetobacter, which had already been identified as significantly less abundant in LMUg compared to HMUg cows in this study (contrast HMUg–LMUg). The further massive abundance decline in LMUg cows from NP to LP feeding suggests Acetobacter an adaptable capacity to different N levels in the rumen, which may contribute to differential N-utilization in LMUg cows.
Transcriptome
While the metabolic activity of the microbial community substantially affects the ruminal protein-N to NPN ratio, the gene expression in the rumen villi impact the ruminal milieu (e.g., active and passive transport processes) and thus the composition of the microbial community (e.g., immune modulation) (Penner et al., 2011; Steele et al., 2016). The pathway analysis of DEGs in rumen villi revealed the enrichment of immune response pathways, which were activated in LMUg compared to HMUg cows. Interestingly, it was recently reported in humans that the intestinal mucosa was able to respond to the microbial products with the targeted expression of receptors, immune genes and signaling molecules to maintain intestinal homeostasis (Tannock and Liu, 2020). In addition, the concept of “immune selection” was described as a way of maintaining host–microbe relationships in the long term as a result of host immune activity (Peterson et al., 2007). In this context, Bessman and Sonnenberg (2016) elucidated the essential role of Antigen Presentation and the Major Histocompatibility Complex, which dictated the host T-cell reactivity to biochemical products of commensal bacteria by the expression of innate lymphoid cells. Notably, the highest Z-score with significant activation in LMUg compared to HMUg cows was found for the Antigen Presentation Pathway in this study. Furthermore, the pathways of IL-4 Signaling, B Cell Development and Glucocorticoid Receptor identified in LMUg–HMUg contrast were enhanced, predominantly by Bovine Leukocyte Antigens (BOLA), which are prominent drivers of the bovine Major Histocompatibility Complex (Ellis and Ballingall, 1999; Lewin et al., 1999). Moreover, Zhang et al. (2021) reported a simultaneous slight activation of the Antigen Presentation Pathway with a significantly higher expression of PRRS2 in the rumen tissue of calves with comparatively higher N-efficiency. In fact, PRSS2 showed a severely higher expression in LMUg compared to HMUg cows. Thus, both antigen presentation and PRRS2 might constitute key candidates distinguishing ruminal N-utilization of LMUg and HMUg. Also, PRSS2 has further been associated with feed efficiency and performance traits in cattle (Abo-Ismail et al., 2013), which might be of interest since MX Dynamin Like GTPase 1 (MX1) and ISG15 Ubiquitin Like Modifier (ISG15) were similarly found to be higher expressed in more feed-efficient heifers, which simultaneously possessed more active immune systems (Paradis et al., 2015). All of these three genes displayed remarkably higher expression levels in LMUg compared to HMUg cows in this study driving the enrichment of Interferon and Neuroinflammation Signaling.
Among the significantly DEG, Matrix Metallopeptidase 3 (MMP3) and TNP2 exhibited the broadest expression differences between cow groups. Metalloproteases were inferred to influence the destruction of tissues by remodeling the extracellular matrix (Coussens et al., 2003; Abendaño et al., 2013) and MMP3 in specific has been recognized to influence the remodeling of adipose tissue in cattle (Jeong et al., 2017). It might be conceivable that MMP3 also attributes to structure processes in the rumen tissue and thus contribute to different absorption and diffusion capacities between HMUg and LMUg epithelia, regarding ruminal NH3, , and urea fluxes. Furthermore, TNP2 has been reported in the context of male fertility (Pasquariello, 2015), but up to now, its function in the rumen of cattle as well as its possible role in N-metabolism remains uncertain. However, given the relationship of fertility traits and immune system activities in dairy cows, the differential expression of TNP2 demands further investigation whether LMUg and HMUg cows differ with respect to fertility parameters (Chebel et al., 2004; Hansen et al., 2004; Moore et al., 2005; König et al., 2008; Weiss et al., 2009; Hurley, 2014).
Although the reduction in dietary CP content did not trigger significantly differential gene expression in rumen tissue of HMUg animals, LMUg cows exhibited severe adaptations in their transcriptome, as characterized by more than 1,000 DEG and 139 significantly enriched pathways. The pathways identified contributed almost exclusively to the immune system, the energy metabolism and the cell cycle. Considering the significantly higher plasma urea pool size during NP compared to LP feeding in cows phenotyped for LMU (Müller et al., 2021), and the general role of immune cells in the rumen wall in protecting the host against infiltration of toxic compounds, the enrichment of these pathways may indicate the necessity for the rumen epithelia to balance the N-gradient and maintain immune homeostasis (Bessman and Sonnenberg, 2016; Tannock and Liu, 2020; Müller et al., 2021). In regard to the severe energy requirements of the immune system (Gleeson et al., 2004; Kvidera et al., 2017), it might further be speculated that the observed effects on energy metabolism pathways in LMUg cows fed the LP ration result from an altered energy demand to stabilize the ruminal N balance.
Relationship between microbial genera and transcriptome subsets
The correlation analysis between microbial and host gene subsets reflecting the rumen microbiome—host interaction, which distinguished the rumen profiles of HMUg and LMUg cows, identified two main clusters. The pathway analysis of Gene Group I revealed the enrichment of Urea Cycle and Arginine Metabolism, primarily due to increased expression of Arginase 1 (ARG1) in LMUg compared to HMUg cows. Arginase is an enzyme that catalyzes the last reaction step in the urea cycle by facilitating the synthesis of arginine into ornithine and urea. In general, the urea cycle is performed predominantly by hepatocytes in the liver (Emmanuel, 1980; Razmi et al., 2005). However, Beck et al. (2009) stated the expression of ARG1 also in the rumen epithelium of cattle, which implies catabolism of N-metabolites via the urea cycle in the rumen wall. Similarly, Emmanuel (1980) identified low ARG1 expression in the rumen epithelium of sheep. Considering the conversion of NH3 into urea, the prominent negative correlation between ARG1 transcripts and UCG-005 abundances calls for further research. Although UCG-005 has not been mentioned directly in the context of N-efficiency in ruminants yet, Amat et al. (2021) predicted UCG-005 a negative impact on Methanobrevibacter abundances in yearling heifers, which is known to contain various N-fixating strains (Poehlein et al., 2018).
Further insights into a possible differentiation of HMUg and LMUg were displayed by Gene Group II, which revealed the enrichment of four immune system pathways, particularly attributed to the expression of BOLA genes. Interestingly, BOLA polymorphisms were not only reported to determine inherent and acquired immune defenses (Rupp and Boichard, 2003), they were also capable to modulate the composition of colostrum microbiota in dairy cows (Derakhshani et al., 2018). Furthermore, BOLA gene expression is known to influence the symbiotic microbial community in the gastrointestinal tract of humans (De Palma et al., 2010), mice (Toivanen et al., 2001), and fish (Bolnick et al., 2014). In the present study, BOLA-DRA transcripts were significantly more abundant in LMUg compared to HMUg cows and showed further striking positive correlations to Roseburia abundances. The contribution of Roseburia to immune responses and to the N-metabolism in cattle, goats, and rabbits has been postulated by various authors (Wang et al., 2011; Alves et al., 2020; Sun et al., 2020). The observed higher abundances of Roseburia and BOLA-DRA transcripts in LMUg compared to HMUg cows might picture a microbial–host interplay that contributes to the distinction of ruminal patterns between HMUg and LMUg cows by affecting immune response and N-metabolism-associated processes.
In addition to the correlation with Roseburia, BOLA-DRA transcripts revealed a severe positive correlation to Lachnospiraceae_FE2018_group abundances, which belongs—like as Roseburia—to the Lachnospiraceae family. Lachnospiraceae_FE2018_group had been mentioned in terms of intestinal barrier function, antioxidant balance, and intestinal inflammation in the gut immune system of broilers (Kong et al., 2020; Liu et al., 2021). Moreover, Wang et al. (2021) postulated a significant positive correlation between Lachnospiraceae family and intestinal inflammation in dairy cows. The Lachnospiraceae family was also found to be significantly more abundant in beef steers with low N-retention efficiency and higher urinary N-excretion (Alves et al., 2020). Furthermore, high abundances of Lachnospiraceae were uncovered in the colon and caecum of goats with low N-utilizing phenotype (Wang et al., 2019). Additionally, Lachnospiraceae occurrences were found to be correlated with a low feed conversion rate and low feed efficiency in beef steers (Hernandez-Sanabria et al., 2010; Carberry et al., 2012).
Conclusion
This study focused on the ruminal background of Holstein cows with predisposed higher or lower MU concentration (HMUg–LMUg) as well as on their adaptation to low CP diets. Compared to HMUg predisposition, LMUg cows displayed higher occurrences of ureolytic genera, such as Succinivibrionaceae_UCG-002 and Ruminococcaceae_unclassified, which might cause lower blood urea concentrations and lower blood urea pool sizes in LMU phenotypes. Similarly, the higher ruminal NH3 concentrations in HMU phenotypes might be attributed to the higher occurrences of high ammonia-producing species hosted by HMUg cows. If enhanced immune responses that were uncovered in LMUg compared to HMUg cows' rumen epithelia influence the epithelial barrier for NPN molecules and thus drive MU phenotype distinction can only be speculated at this point. However, the downregulation of immune responses and energy metabolism pathways in LMUg cows fed a low CP diet might have indicated energy wasting efforts of the rumen epithelia to maintain ruminal N-balance when LMUg cows were exposed to diets with normal CP content. Considering the important role of Arginase in the urea cycle, the observed interplay between ARG1 expression and UCG-005 occurrences is specifically proposed for future research on microbe–host interactions. The ruminal patterns identified for HMUg and LMUg cows contribute to a deeper insight into MU predisposition of Holsteins and attribute to the optimization of N-utilization in tandem with the reduction of N-emissions on dairy farms by future breeding selection strategies.
Data availability statement
The data presented in the study are deposited in the ArrayExpress Microarray Database at EBI (RNA datasets) and in the BioSample Database (16S Amplicon data) repository, accession numbers E-MTAB-9901 (RNA datasets) and PRJNA856508 (16S Amplicon data).
Ethics statement
The animal study was reviewed and approved by Ethics Committee of the State of Mecklenburg-Western Pomerania (State Office for Agriculture, Food Safety and Fisheries; LALLF M-V7221.3-2-019/19). Written informed consent was obtained from the owners for the participation of their animals in this study.
Author contributions
CM, MP, and BK organized and carried out animal husbandry. HH and HR jointly collected the samples, performed statistical analysis, and interpreted the data. HH conducted laboratory work and wrote the manuscript. NT performed RNA and 16S sequencing. SP supported data analysis. DS assisted in data interpretation. NR, BK, and KW conceptualized and supervised the study. All authors reviewed the final manuscript and provided critical feedback.
Funding
As a part of the BlueCow Project (FKZ: 281B101516) this study has received funding from the Federal Ministry of Food and Agriculture.
Conflict of interest
Author DS was employed by IT-Solutions for Animal Production, Vereinigte Informationssysteme Tierhaltung w.V. (vit).
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/fmicb.2022.939711/full#supplementary-material
Supplementary material 1. p. 1: Significantly differentially expressed genes (DEG) between HMUg and LMUg cows. Log2fold changes refer to HMUg. p. 2: Gene set for enrichment analysis (p-values are below 0.01, and the corresponding q-values are below 0.6) between HMUg and LMUg cows. Log2fold changes refer to HMUg.
Supplementary material 2. p. 1: Gene set for enrichment analysis (p-values are below 0.01, and the corresponding q-values are below 0.6) between NP and LP diet in HMUg cows. Log2fold changes refer to HMUg. p. 2: Gene set for enrichment analysis (p-values are below 0.01, and the corresponding q-values are below 0.6) between NP and LP diet in LMUg cows. Log2fold changes refer to NP. p. 3: Significantly enriched pathways (adjusted p-values are below 0.05) between NP and LP diet in LMUg cows.
Supplementary material 3. p. 1: Microbial subset (53 genera) generated by sparse Partial Least Squares Discriminant Analysis (sPLS-DA) separating HMUg and LMUg cows. p. 2: Gene subset (100 genes) generated by Partial Least Squares Discriminant Analysis (sPLS-DA) separating HMUg and LMUg cows.
Supplementary material 4. p. 1: The p-values of the correlation analysis between microbial genera abundances and gene transcripts separating HMUg and LMUg cows (Figure 3). p. 2: Correlation coefficients (r-values) of the correlation analysis between microbial genera abundances and gene transcripts separating HMUg and LMUg cows (Figure 3).
References
Abdoun, K., Stumpff, F., and Martens, H. (2006). Ammonia and urea transport across the rumen epithelium: a review. Anim. Health Res. Rev. 7, 43–59. doi: 10.1017/S1466252307001156
Abendaño, N., Juste, R. A., and Alonso-Hearn, M. (2013). Anti-inflammatory and antiapoptotic responses to infection: a common denominator of human and bovine macrophages infected with Mycobacterium avium subsp. paratuberculosis. BioMed. Res. Int. 2013:7. doi: 10.1155/2013/908348
Abo-Ismail, M., Kelly, M., Squires, E., Swanson, K., Bauck, S., and Miller, S. (2013). Identification of single nucleotide polymorphisms in genes involved in digestive and metabolic processes associated with feed efficiency and performance traits in beef cattle. J. Animal Sci. 91, 2512–2529. doi: 10.2527/jas.2012-5756
Aguiar, S. C. D., Paula, E. M. D., Yoshimura, E. H., Santos, W. B. R. D., Machado, E., Valero, M. V., et al. (2014). Effects of phenolic compounds in propolis on digestive and ruminal parameters in dairy cows. Rev. Brasil. Zootecn. 43, 197–206. doi: 10.1590/S1516-35982014000400006
Alves, K. L. G. C., Granja-Salcedo, Y. T., Messana, J. D., De Souza, V. C., Ganga, M. J. G., Colovate, P. H. D., et al. (2020). Rumen bacterial diversity in relation to nitrogen retention in beef cattle. Anaerobe 67, 102316. doi: 10.1016/j.anaerobe.2020.102316
Amat, S., Holman, D. B., Schmidt, K., Menezes, A. C. B., Baumgaertner, F., Winders, T., et al. (2021). The nasopharyngeal, ruminal, and vaginal microbiota and the core taxa shared across these microbiomes in virgin yearling heifers exposed to divergent in utero nutrition during their first trimester of gestation and in pregnant beef heifers in response to mineral supplementation. bioRxiv. 9, 2011. doi: 10.3390/microorganisms9102011
Anders, S., Pyl, P. T., and Huber, W. (2015). HTSeq—a Python framework to work with high-throughput sequencing data. Bioinformatics 31, 166–169. doi: 10.1093/bioinformatics/btu638
Attwood, G., and Reilly, K. (1995). Identification of proteolytic rumen bacteria isolated from New Zealand cattle. J. Appl. Bacteriol. 79, 22–29. doi: 10.1111/j.1365-2672.1995.tb03119.x
Bailoni, L., Carraro, L., Cardin, M., and Cardazzo, B. (2021). Active rumen bacterial and protozoal communities revealed by RNA-based amplicon sequencing on dairy cows fed different diets at three physiological stages. Microorganisms 9, 754. doi: 10.3390/microorganisms9040754
Beatson, P., Meier, S., Cullen, N., and Eding, H. (2019). Genetic variation in milk urea nitrogen concentration of dairy cattle and its implications for reducing urinary nitrogen excretion. Animal 13, 2164–2171. doi: 10.1017/S1751731119000235
Beck, D., Jiang, H., and Zhang, L. (2009). Elucidating the evolutionary relationships among bos taurus digestive organs using unigene expression data. Int. J. Evolut. Biol. 2010, 803142. doi: 10.4061/2009/803142
Bento, C. B. P., De Azevedo, A. C., Detmann, E., and Mantovani, H. C. (2015). Biochemical and genetic diversity of carbohydrate-fermenting and obligate amino acid-fermenting hyper-ammonia-producing bacteria from Nellore steers fed tropical forages and supplemented with casein. BMC Microbiol. 15, 28. doi: 10.1186/s12866-015-0369-9
Bergen, W. G. (2021). “Amino acids in beef cattle nutrition and production,” in Amino Acids in Nutrition and Health (New York, NY: Springer), 29–42. doi: 10.1007/978-3-030-54462-1_3
Bessman, N. J., and Sonnenberg, G. F. (2016). Emerging roles for antigen presentation in establishing host–microbiome symbiosis. Immunol. Rev. 272, 139–150. doi: 10.1111/imr.12425
Bolnick, D. I., Snowberg, L. K., Caporaso, J. G., Lauber, C., Knight, R., and Stutz, W. E. (2014). Major H istocompatibility C omplex class II b polymorphism influences gut microbiota composition and diversity. Molecular ecology 23, 4831–4845. doi: 10.1111/mec.12846
Bryant, M. (1970). Normal flora—rumen bacteria. Am. J. Clin. Nutr. 23, 1440–1450. doi: 10.1093/ajcn/23.11.1440
Burgos, S. A., Embertson, N. M., Zhao, Y., Mitloehner, F. M., Depeters, E. J., and Fadel, J. G. (2010). Prediction of ammonia emission from dairy cattle manure based on milk urea nitrogen: relation of milk urea nitrogen to ammonia emissions. J. Dairy Sci. 93, 2377–2386. doi: 10.3168/jds.2009-2415
Burgos, S. A., Fadel, J. G., and Depeters, E. J. (2007). Prediction of ammonia emission from dairy cattle manure based on milk urea nitrogen: relation of milk urea nitrogen to urine urea nitrogen excretion. J. Dairy Sci. 90, 5499–5508. doi: 10.3168/jds.2007-0299
Calsamiglia, S., Ferret, A., Reynolds, C. K., Kristensen, N. B., and Van Vuuren, A. M. (2010). Strategies for optimizing nitrogen use by ruminants. Animal 4, 1184–1196. doi: 10.1017/S1751731110000911
Cantalapiedra-Hijar, G., Ortigues-Marty, I., Sepchat, B., Titgemeyer, E., and Bahloul, L. (2020). Methionine-balanced diets improve cattle performance in fattening young bulls fed high-forage diets through changes in nitrogen metabolism. Br. J. Nutr. 124, 273–285. doi: 10.1017/S0007114520001154
Carberry, C. A., Kenny, D. A., Han, S., Mccabe, M. S., and Waters, S. M. (2012). Effect of phenotypic residual feed intake and dietary forage content on the rumen microbial community of beef cattle. Appl. Environ. Microbiol. 78, 4949–4958. doi: 10.1128/AEM.07759-11
Castillo, A., Kebreab, E., Beever, D., and France, J. (2000). A review of efficiency of nitrogen utilisation in lactating dairy cows and its relationship with environmental pollution. J. Animal Feed Sci. 9, 1–32. doi: 10.22358/jafs/68025/2000
Chebel, R. C., Santos, J. E., Reynolds, J. P., Cerri, R. L., Juchem, S. O., and Overton, M. (2004). Factors affecting conception rate after artificial insemination and pregnancy loss in lactating dairy cows. Animal Reprod. Sci. 84, 239–255. doi: 10.1016/j.anireprosci.2003.12.012
Chen, Y., Vanderick, S., Mota, R., Grelet, C., Gengler, N., and Consortium, G. (2021). Estimation of genetic parameters for predicted nitrogen use efficiency and losses in early lactation of Holstein cows. J. Dairy Sci. 104, 4413–4423. doi: 10.3168/jds.2020-18849
Cholewińska, P., Czyz, K., Nowakowski, P., and Wyrostek, A. (2020). The microbiome of the digestive system of ruminants: a review. Animal Health Res. Rev. 21, 3–14. doi: 10.1017/S1466252319000069
Chun, H., and Keleş, S. (2010). Sparse partial least squares regression for simultaneous dimension reduction and variable selection. J. R. Stat. Soc. Ser. B 72, 3–25. doi: 10.1111/j.1467-9868.2009.00723.x
Correa-Luna, M., Donaghy, D., Kemp, P., Schutz, M., and López-Villalobos, N. (2021). Nitrogen use efficiency and excretion in grazing cows with high and low milk urea nitrogen breeding values. Sustainability 13, 9827. doi: 10.3390/su13179827
Coussens, P. M., Colvin, C. J., Rosa, G. J., Perez Laspiur, J., and Elftman, M. D. (2003). Evidence for a novel gene expression program in peripheral blood mononuclear cells from Mycobacterium avium subsp. paratuberculosis-infected cattle. Infect. Immun. 71, 6487–6498. doi: 10.1128/IAI.71.11.6487-6498.2003
De Mendiburu Delgado, F. (2009). Una herramienta de análisis estadístico para la investigaciónagrícola. Lima: Universidad Nacional de Ingeniería.
De Palma, G., Capilla, A., Nadal, I., Nova, E., Pozo, T., Varea, V., et al. (2010). Interplay between human leukocyte antigen genes and the microbial colonization process of the newborn intestine. Curr. Issues Mol. Biol. 12, 1–10. doi: 10.21775/cimb.012.001
Derakhshani, H., Plaizier, J. C., De Buck, J., Barkema, H. W., and Khafipour, E. (2018). Association of bovine major histocompatibility complex (BoLA) gene polymorphism with colostrum and milk microbiota of dairy cows during the first week of lactation. Microbiome 6, 1–18. doi: 10.1186/s40168-018-0586-1
Dias, J., Marcondes, M. I., Noronha, M. F., Resende, R. T., Machado, F. S., Mantovani, H. C., et al. (2017). Effect of pre-weaning diet on the ruminal archaeal, bacterial, and fungal communities of dairy calves. Front. Microbiol. 8, 1553. doi: 10.3389/fmicb.2017.01553
Ellis, S. A., and Ballingall, K. T. (1999). Cattle MHC: evolution in action? Immunol. Rev. 167, 159–168. doi: 10.1111/j.1600-065X.1999.tb01389.x
Emmanuel, B. (1980). Urea cycle enzymes in tissues (liver, rumen epithelium, heart, kidney, lung and spleen) of sheep (Ovis aries). Comp. Biochem. Physiol. B Comp. Biochem. 65, 693–697. doi: 10.1016/0305-0491(80)90181-9
Eschenlauer, S., Mckain, N., Walker, N., Mcewan, N., Newbold, C., and Wallace, R. (2002). Ammonia production by ruminal microorganisms and enumeration, isolation, and characterization of bacteria capable of growth on peptides and amino acids from the sheep rumen. Appl. Environ. Microbiol. 68, 4925–4931. doi: 10.1128/AEM.68.10.4925-4931.2002
Gardiner, C., Clough, T. J., Cameron, K. C., Di, H. J., Edwards, G., and De Klein, C. (2016). Potential for forage diet manipulation in New Zealand pasture ecosystems to mitigate ruminant urine derived N2O emissions: a review. N. Zeal. J. Agric. Res. 59, 301–317. doi: 10.1080/00288233.2016.1190386
Gleeson, M., Nieman, D. C., and Pedersen, B. K. (2004). Exercise, nutrition and immune function. J. Sports Sci. 22, 115–125. doi: 10.1080/0264041031000140590
Guliński, P., Salamończyk, E., and Młynek, K. (2016). Improving nitrogen use efficiency of dairy cows in relation to urea in milk: a review. Anim. Sci. Pap. Rep. 34, 5–24.
Hansen, P. J., Soto, P., and Natzke, R. P. (2004). Mastitis and fertility in cattle–possible involvement of inflammation or immune activation in embryonic mortality. Am. J. Reprod. Immunol. 51, 294–301. doi: 10.1111/j.1600-0897.2004.00160.x
Hartinger, T., Gresner, N., and Südekum, K.-H. (2018). Does intra-ruminal nitrogen recycling waste valuable resources? A review of major players and their manipulation. J. Anim. Sci. Biotechnol. 9, 1–21. doi: 10.1186/s40104-018-0249-x
Henderson, G., Cox, F., Ganesh, S., Jonker, A., Young, W., and Janssen, P. H. (2015). Rumen microbial community composition varies with diet and host, but a core microbiome is found across a wide geographical range. Sci. Rep. 5, 1–15. doi: 10.1038/srep14567
Hernandez-Sanabria, E., Guan, L. L., Goonewardene, L. A., Li, M., Mujibi, D. F., Stothard, P., et al. (2010). Correlation of particular bacterial PCR-denaturing gradient gel electrophoresis patterns with bovine ruminal fermentation parameters and feed efficiency traits. Appl. Environ. Microbiol. 76, 6338–6350. doi: 10.1128/AEM.01052-10
Honerlagen, H., Reyer, H., Oster, M., Ponsuksili, S., Trakooljul, N., Kuhla, B., et al. (2021). Identification of genomic regions influencing N-metabolism and N-excretion in lactating holstein-friesians. Front. Genet. 12, 699550. doi: 10.3389/fgene.2021.699550
Hu, R., Zou, H., Wang, Z., Cao, B., Peng, Q., Jing, X., et al. (2019). Nutritional interventions improved rumen functions and promoted compensatory growth of growth-retarded yaks as revealed by integrated transcripts and microbiome analyses. Front. Microbiol. 10, 318. doi: 10.3389/fmicb.2019.00318
Hurley, D.J. (2014). “The relationship of immunity and reproduction in dairy cows,” in Proceedings of the 25th Annual Florida Ruminant Nutrition Symposium (Gainesville, FL: University of Florida), 1.
Jahnel, R. E., Wittenburg, D., Mayer, M., Taeubert, H., and Reinsch, N. (2021). Estimation of genetic parameters of milk urea content. Zuchtungskunde 93, 157–168.
Jeong, J. Y., Ibrahim, M., Kim, M. J., So, K., Jeong, Y. D., Park, S., et al. (2017). Comparisons of extracellular matrix-related gene expression levels in different adipose tissues from Korean cattle. Livestock Sci. 198, 138–146. doi: 10.1016/j.livsci.2017.02.020
Jin, D., Zhao, S., Wang, P., Zheng, N., Bu, D., Beckers, Y., et al. (2016). Insights into abundant rumen ureolytic bacterial community using rumen simulation system. Front. Microbiol. 7, 1006. doi: 10.3389/fmicb.2016.01006
Jin, D., Zhao, S., Zheng, N., Beckers, Y., and Wang, J. (2018). Urea metabolism and regulation by rumen bacterial urease in ruminants: a review. Ann. Anim. Sci. 18, 303. doi: 10.1515/aoas-2017-0028
Kim, D., Langmead, B., and Salzberg, S. L. (2015). HISAT: a fast spliced aligner with low memory requirements. Nat. Methods 12, 357–360. doi: 10.1038/nmeth.3317
Kim, D., Paggi, J. M., Park, C., Bennett, C., and Salzberg, S. L. (2019). Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat. Biotechnol. 37, 907–915. doi: 10.1038/s41587-019-0201-4
Kong, F., Gao, Y., Tang, M., Fu, T., Diao, Q., Bi, Y., et al. (2020). Effects of dietary rumen–protected Lys levels on rumen fermentation and bacterial community composition in Holstein heifers. Appl. Microbiol. Biotechnol. 104, 6623–6634. doi: 10.1007/s00253-020-10684-y
König, S., Chang, Y., Borstel, U., Gianola, D., and Simianer, H. (2008). Genetic and phenotypic relationships among milk urea nitrogen, fertility, and milk yield in Holstein cows. J. Dairy Sci. 91, 4372–4382. doi: 10.3168/jds.2008-1236
Kotz, A., Azevedo, P. A., Khafipour, E., and Plaizier, J. (2020). Effects of the dietary grain content on rumen and fecal microbiota of dairy cows. Can. J. Anim. Sci. 101, 274–286. doi: 10.1139/cjas-2020-0122
Kozich, J. J., Westcott, S. L., Baxter, N. T., Highlander, S. K., and Schloss, P. D. (2013). Development of a dual-index sequencing strategy and curation pipeline for analyzing amplicon sequence data on the MiSeq Illumina sequencing platform. Appl. Environ. Microbiol. 79, 5112–5120. doi: 10.1128/AEM.01043-13
Kvidera, S., Horst, E., Abuajamieh, M., Mayorga, E., Fernandez, M. S., and Baumgard, L. (2017). Glucose requirements of an activated immune system in lactating Holstein cows. J. Dairy Sci. 100, 2360–2374. doi: 10.3168/jds.2016-12001
Lapierre, H., Berthiaume, R., Raggio, G., Thivierge, M., Doepel, L., Pacheco, D., et al. (2005). The route of absorbed nitrogen into milk protein. Anim. Sci. 80, 11–22. doi: 10.1079/ASC41330011
Lavery, A., and Ferris, C. P. (2021). Proxy measures and novel strategies for estimating nitrogen utilisation efficiency in dairy cattle. Animals 11, 343. doi: 10.3390/ani11020343
Lewin, H. A., Russell, G. C., and Glass, E. J. (1999). Comparative organization and function of the major histocompatibility complex of domesticated cattle. Immunol. Rev. 167, 145–158. doi: 10.1111/j.1600-065X.1999.tb01388.x
Ling, J., and Armstead, I. (1995). The in vitro uptake and metabolism of peptides and amino acids by five species of rumen bacteria. J. Appl. Bacteriol. 78, 116–124. doi: 10.1111/j.1365-2672.1995.tb02831.x
Liu, W.-C., Guo, Y., An, L.-L., and Zhao, Z.-H. (2021). Protective effects of dietary betaine on intestinal barrier function and cecal microbial community in indigenous broiler chickens exposed to high temperature environment. Environ. Sci. Pollut. Res. 28, 10860–10871. doi: 10.1007/s11356-020-11326-6
Loubinoux, J., Valente, F. M., Pereira, I. A., Costa, A., Grimont, P. A., and Le Faou, A. E. (2002). Reclassification of the only species of the genus Desulfomonas, Desulfomonas pigra, as Desulfovibrio piger comb. nov. Int. J. Syst. Evolut. Microbiol. 52, 1305–1308.: doi: 10.1099/00207713-52-4-1305
Love, M. I., Huber, W., and Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15, 1–21. doi: 10.1186/s13059-014-0550-8
Mcmurdie, P. J., and Holmes, S. (2014). Waste not, want not: why rarefying microbiome data is inadmissible. PLoS Comput. Biol. 10, e1003531. doi: 10.1371/journal.pcbi.1003531
Meehan, C. J., and Beiko, R. G. (2014). A phylogenomic view of ecological specialization in the Lachnospiraceae, a family of digestive tract-associated bacteria. Genome Biol. Evolut. 6, 703–713. doi: 10.1093/gbe/evu050
Miglior, F., Sewalem, A., Jamrozik, J., Bohmanova, J., Lefebvre, D. M., and Moore, R. K. (2007). Genetic analysis of milk urea nitrogen and lactose and their relationships with other production traits in Canadian Holstein cattle. J. Dairy Sci. 90, 2468–2479. doi: 10.3168/jds.2006-487
Moore, D. A., Overton, M. W., Chebel, R. C., Truscott, M. L., and Bondurant, R. H. (2005). Evaluation of factors that affect embryonic loss in dairy cattle. J. Am. Vet. Med. Assoc. 226, 1112–1118. doi: 10.2460/javma.2005.226.1112
Müller, C. B. M., Görs, S., Derno, M., Tuchscherer, A., Wimmers, K., Zeyner, A., et al. (2021). Differences between Holstein dairy cows in renal clearance rate of urea affect milk urea concentration and the relationship between milk urea and urinary nitrogen excretion. Sci. Total Environ. 755, 143–198. doi: 10.1016/j.scitotenv.2020.143198
Nousiainen, J., Shingfield, K., and Huhtanen, P. (2004). Evaluation of milk urea nitrogen as a diagnostic of protein feeding. J. Dairy Sci. 87, 386–398. doi: 10.3168/jds.S0022-0302(04)73178-1
Oksanen, J., Blanchet, F. G., Kindt, R., Legendre, P., Minchin, P. R., O'hara, R., et al. (2013). Package ‘vegan'. Commun. Ecol. 2, 1–295. Available online at: http://cran.r-project.org/
Pacífico, C., Petri, R. M., Ricci, S., Mickdam, E., Wetzels, S. U., Neubauer, V., et al. (2021). Unveiling the bovine epimural microbiota composition and putative function. Microorganisms 9, 342. doi: 10.3390/microorganisms9020342
Paradis, F., Yue, S., Grant, J., Stothard, P., Basarab, J., and Fitzsimmons, C. (2015). Transcriptomic analysis by RNA sequencing reveals that hepatic interferon-induced genes may be associated with feed efficiency in beef heifers. J. Animal Sci. 93, 3331–3341. doi: 10.2527/jas.2015-8975
Pasquariello, R. (2015). Molecular Characterization of Factors Which Can Influence the Reproductive Success of Cattle. Milano : Universitá degli Studi di Milano, Dipartimento di Scienze Agrarie e Ambientali.
Patra, A. K., and Aschenbach, J. R. (2018). Ureases in the gastrointestinal tracts of ruminant and monogastric animals and their implication in urea-N/ammonia metabolism: a review. J. Adv. Res. 13, 39–50. doi: 10.1016/j.jare.2018.02.005
Penner, G., Steele, M., Aschenbach, J., and Mcbride, B. (2011). Ruminant nutrition symposium: molecular adaptation of ruminal epithelia to highly fermentable diets. J. Anim. Sci. 89, 1108–1119. doi: 10.2527/jas.2010-3378
Pertea, M., Kim, D., Pertea, G. M., Leek, J. T., and Salzberg, S. L. (2016). Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown. Nat. Protocols 11, 1650–1667. doi: 10.1038/nprot.2016.095
Peterson, D. A., Mcnulty, N. P., Guruge, J. L., and Gordon, J. I. (2007). IgA response to symbiotic bacteria as a mediator of gut homeostasis. Cell Host Microbe 2, 328–339. doi: 10.1016/j.chom.2007.09.013
Plaizier, J. C., Li, S., Tun, H. M., and Khafipour, E. (2017). Nutritional models of experimentally-induced subacute ruminal acidosis (SARA) differ in their impact on rumen and hindgut bacterial communities in dairy cows. Front. Microbiol. 7, 2128. doi: 10.3389/fmicb.2016.02128
Poehlein, A., Schneider, D., Soh, M., Daniel, R., and Seedorf, H. (2018). Comparative genomic analysis of members of the genera Methanosphaera and Methanobrevibacter reveals distinct clades with specific potential metabolic functions. Archaea 2018. doi: 10.1155/2018/7609847
Razmi, N., Ali Jelodar, G., Nazifi, S., and Dehghani, A. (2005). Arginase status in cattle reproductive system. Veterinarski arhiv 75, 31–38.
Rohart, F., Gautier, B., Singh, A., and Cao, K.-A. (2017). mixOmics: An R package for ‘omics feature selection and multiple data integration. PLoS Comput. Biol. 13, e1005752. doi: 10.1371/journal.pcbi.1005752
Rosendahl, J. (2015). Untersuchung über die beteiligung von nicht-selektiven kationenkanälen am ammoniaktransport über das pansenepithel des rindes (Dissertation). Berlin: Freie Universität Berlin, Institut für Veterinär-Physiologie des Fachbereichs Veterinärmedizin.
Rupp, R., and Boichard, D. (2003). Genetics of resistance to mastitis in dairy cattle. Vet. Res. 34, 671–688. doi: 10.1051/vetres:2003020
Russell, J. B., Muck, R. E., and Weimer, P. J. (2009). Quantitative analysis of cellulose degradation and growth of cellulolytic bacteria in the rumen. FEMS Microbiol. Ecol. 67, 183–197. doi: 10.1111/j.1574-6941.2008.00633.x
Sales, M., Lucas, F., and Blanchart, G. (2000). Effects of ammonia and amino acids on the growth and proteolytic activity of three species of rumen bacteria: prevotella albensis, Butyrivibrio fibrisolvens, and Streptococcus bovis. Curr. Microbiol. 40, 380–386. doi: 10.1007/s002840010074
Schloss, P. D., Westcott, S. L., Ryabin, T., Hall, J. R., Hartmann, M., Hollister, E. B., et al. (2009). Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl. Environ. Microbiol. 75, 7537–7541. doi: 10.1128/AEM.01541-09
Spek, J. W., Dijkstra, J., Van Duinkerken, G., and Bannink, A. (2012). A review of factors influencing milk urea concentration and its relationship with urinary urea excretion in lactating dairy cattle. J. Agric. Sci. 151, 407–423. doi: 10.1017/S0021859612000561
Steele, M. A., Penner, G. B., and Chaucheyras-Durand, F. (2016). Development and physiology of the rumen and the lower gut: targets for improving gut health. J. Dairy Sci. 99, 4955–4966. doi: 10.3168/jds.2015-10351
Suen, G., Stevenson, D. M., Bruce, D. C., Chertkov, O., Copeland, A., Cheng, J.-F., et al. (2011). Complete genome of the cellulolytic ruminal bacterium Ruminococcus albus . Am. Soc. Microbiol. 193, 5574–5575. doi: 10.1128/JB.05621-11
Sun, H. Z., Zhou, M., Wang, O., Chen, Y., Liu, J. X., and Guan, L. L. (2020). Multi-omics reveals functional genomic and metabolic mechanisms of milk production and quality in dairy cows. Bioinformatics 36, 2530–2537. doi: 10.1093/bioinformatics/btz951
Tan, P., Liu, H., Zhao, J., Gu, X., Wei, X., Zhang, X., et al. (2021). Amino acids metabolism by rumen microorganisms: nutrition and ecology strategies to reduce nitrogen emissions from the inside to the outside. Sci. Total Environ. 800, 149–596. doi: 10.1016/j.scitotenv.2021.149596
Tannock, G. W., and Liu, Y. (2020). Guided dietary fibre intake as a means of directing short-chain fatty acid production by the gut microbiota. J. R. Soc. N. Zeal. 50, 434–455. doi: 10.1080/03036758.2019.1657471
Te Pas, M. F., Veldkamp, T., De Haas, Y., Bannink, A., and Ellen, E. D. (2021). Adaptation of livestock to new diets using feed components without competition with human edible protein sources—a review of the possibilities and recommendations. Animals 11, 2293. doi: 10.3390/ani11082293
Toivanen, P., Vaahtovuo, J., and Eerola, E. (2001). Influence of major histocompatibility complex on bacterial composition of fecal flora. Infect. Immun. 69, 2372–2377. doi: 10.1128/IAI.69.4.2372-2377.2001
Trabi, E. B., Seddik, H.-E., Xie, F., Wang, X., Liu, J., and Mao, S. (2020). Effect of pelleted high-grain total mixed ration on rumen morphology, epithelium-associated microbiota and gene expression of proinflammatory cytokines and tight junction proteins in Hu sheep. Anim. Feed Sci. Technol. 263, 114453. doi: 10.1016/j.anifeedsci.2020.114453
Tun, H. M., Li, S., Yoon, I., Meale, S. J., Azevedo, P. A., Khafipour, E., et al. (2020). Saccharomyces cerevisiae fermentation products (SCFP) stabilize the ruminal microbiota of lactating dairy cows during periods of a depressed rumen pH. BMC Vet. Res. 16, 237. doi: 10.1186/s12917-020-02437-w
Uwizeye, A., De Boer, I. J. M., Opio, C. I., Schulte, R. P. O., Falcucci, A., Tempio, G., et al. (2020). Nitrogen emissions along global livestock supply chains. Nat. Food 1, 437–446. doi: 10.1038/s43016-020-0113-y
Vereinigte Informationssysteme Tierhaltung (VIT). (2020). Beschreibung der Zuchtwertschätzung für alle Schätzmerkmale bei den Milchrinderrassen für die vit mit der Zuchtwertschätzung beauftragt ist. Available online at: https://www.vit.de/fileadmin/DE/Zuchtwertschaetzung/Zws_Bes_deu.pdf. (accessed January 01, 2020).
Wallace, R., and Brammall, M. L. (1985). The role of different species of bacteria in the hydrolysis of protein in the rumen. Microbiology 131, 821–832. doi: 10.1099/00221287-131-4-821
Wang, L., Liu, K., Wang, Z., Bai, X., Peng, Q., and Jin, L. (2019). Bacterial community diversity associated with different utilization efficiencies of nitrogen in the gastrointestinal tract of goats. Front. Microbiol. 10, 239. doi: 10.3389/fmicb.2019.00239
Wang, M.-Z., Shen, J.-Z., Liu, Y., Dong, S.-H., Yu, L.-H., and Wang, H.-R. (2011). Effects of dietary protein sources on bacterial and protozoal community structure in goat rumen. Chin. J. Anim. Sci. 3:1–12. doi: 10.1371/journal.pone.0221815
Wang, Y., Nan, X., Zhao, Y., Jiang, L., Wang, M., Wang, H., et al. (2021). Rumen microbiome structure and metabolites activity in dairy cows with clinical and subclinical mastitis. J. Anim. Sci. Biotechnol. 12, 1–21. doi: 10.1186/s40104-020-00543-1
Weiss, G., Goldsmith, L. T., Taylor, R. N., Bellet, D., and Taylor, H. S. (2009). Inflammation in reproductive disorders. Reprod. Sci. 16, 216–229. doi: 10.1177/1933719108330087
Wood, G., Boettcher, P., Jamrozik, J., Jansen, G., and Kelton, D. (2003). Estimation of genetic parameters for concentrations of milk urea nitrogen. J. Dairy Sci. 86, 2462–2469. doi: 10.3168/jds.S0022-0302(03)73840-5
Zhang, N., Teng, Z., Li, P., Fu, T., Lian, H., Wang, L., et al. (2021). Oscillating dietary crude protein concentrations increase N retention of calves by affecting urea-N recycling and nitrogen metabolism of rumen bacteria and epithelium. PLoS ONE 16, e0257417. doi: 10.1371/journal.pone.0257417
Keywords: rumen microbiome, host gene expression, microbe–host relationship, dairy cow, milk urea
Citation: Honerlagen H, Reyer H, Segelke D, Müller CBM, Prahl MC, Ponsuksili S, Trakooljul N, Reinsch N, Kuhla B and Wimmers K (2022) Ruminal background of predisposed milk urea (MU) concentration in Holsteins. Front. Microbiol. 13:939711. doi: 10.3389/fmicb.2022.939711
Received: 09 May 2022; Accepted: 26 July 2022;
Published: 13 September 2022.
Edited by:
Alejandro Belanche, University of Zaragoza, SpainReviewed by:
Francesco Rubino, Queen's University Belfast, United KingdomJuan Manuel Palma-Hidalgo, Scotland's Rural College, United Kingdom
Copyright © 2022 Honerlagen, Reyer, Segelke, Müller, Prahl, Ponsuksili, Trakooljul, Reinsch, Kuhla and Wimmers. 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: Klaus Wimmers, wimmers@fbn-dummerstorf.de