No hints at glyphosate-induced ruminal dysbiosis in cows

Glyphosate-based herbicides are among the most used non-selective herbicides worldwide and inhibit synthesis of aromatic amino acids in plants, bacteria, and fungi. Given the broad usage, controversies concerning potential effects of glyphosate on health and especially on gut microbiomes arose. For cattle, it has been proposed based on in vitro data that glyphosate has detrimental effects on the ruminal microbiome, which manifest as a specific inhibition of bacteria involved in fiber degradation and as an enrichment of specific pathogens. In the present study, glyphosate effects on the ruminal microbiome were analyzed in vivo using glyphosate contaminated feedstuffs with strong differences in dietary fiber and dietary energy content in order to reproduce the proposed detrimental glyphosate effects on the rumen microbiome. While significant impact of dietary factors on the ruminal microbiome and its products are pointed out, no adverse glyphosate effects on ruminal microbiome composition, diversity, and microbial metabolites are observed.


INTRODUCTION
Glyphosate is one of the most used active substances in common non-selective herbicides worldwide and was considered an advantageous herbicide until the first glyphosate-resistant weeds appeared 1 . Supported by introduction of glyphosate-resistant genetically modified crops in the 1990s, glyphosate applications and concerns about applications steadily increased [2][3][4] . Beneath concerns about human exposure to glyphosate, especially the continuous exposure of livestock to glyphosate is under observation. Usage of glyphosate-based herbicides has been common in Germany for an extended period of time despite increasing strict EU regulations and more than 25% of utilized agricultural areas in Germany had an application of glyphosate-based herbicides in 2009 5 . However, soybean meal from glyphosateresistant soybeans is considered the main source of glyphosate contaminations in German cattle feed, while residues from local crops are negligible 6 . Since glyphosate contaminations in feed production and consequently also ingestion by livestock animals are common, potential exposure and effects of glyphosate on livestock have been evaluated frequently [7][8][9] . For dairy cows, the main proportion of consumed glyphosate, which ranges in average daily exposures from 0.08 to 6.7 mg per cow, is excreted via feces and urine, while only negligible amounts of glyphosate are excreted via milk 6 . Health and performance of dairy cows did not respond to glyphosate contaminations in feed [10][11][12] . Glyphosate inhibits synthesis of aromatic amino acids through inhibition 5-enolpyruvyl-shikimic acid-3-phosphate synthase in plants and susceptible bacteria [13][14][15] . Accordingly, a potential impact of glyphosate on microbiota in digestive organs, which are exposed on regular basis and directly linked to animal health, is in the focus of interest. Recommended glyphosate application has only limited effects on soil microbiomes and different microbial degradation pathways for glyphosate are known [16][17][18] . It was also pointed out that glyphosate has a limited short-term effect on gut microbiota upon ingestion due to sufficient amino acid supply in the gastrointestinal tract 19 . Rumen microbial composition is associated with feed composition, feed efficiency phenotypes, and performance effects, while microbiome disturbances are considered a basis for detrimental effects like pathogen enrichment [20][21][22][23] . In this context variations of concentrate feed proportion (CFP) and dietary fiber content are of special interest, since the chemical composition of diets defines the ruminal fermentative pattern and especially high CFP can lead to adverse effects including subacute ruminal acidosis and ruminal dysbiosis, which is characterized by strong shifts in microbiome composition [24][25][26] . Different studies analyzed a potential introduction of disturbances to the ruminal microbiome by glyphosate in vitro in the context of high CFP challenges and varying dietary fiber contents. The results were conflicting and ranged from no significant glyphosate effects in mixed cultures to proposing an enrichment of Clostridium (C.) botulinum and an increased neurotoxin gene expression [27][28][29][30][31] . Subsequently, such enrichment of C. botulinum, coupled to increased neurotoxin gene expression in the gut of dairy cows, has been controversially discussed in the context of a highly disputed chronic form of botulism [32][33][34][35][36] . However, data concerning the effects of glyphosate in the context of dietary challenges on the ruminal microbiome and putative consequences for animal health from in vivo trials in "real-life" scenarios is still limited.
To address the lack of knowledge about the influences of glyphosate on dairy cows in general and in order to evaluate discrepancies between the different studies, a broad general health screening of dairy cows under realistic housing conditions with different CFP and dietary fiber content in feedstuffs was conducted. Within this screening, influences of glyphosate on animal performance, liver histology, liver gene expression, and functionality of blood cells were tested and occurrence of significant effects was rejected [10][11][12] . To achieve this, cows were assigned to 4 groups receiving diets with a high CFP (60%, HC) or a low CFP (30%, LC) and with (GLY, glyphosate exposure per cow: >73 mg d −1 ; >112 μg kg(bodyweight) −1 d −1 ) or without (CON, glyphosate exposure per cow: <1 mg d −1 ; <1.3 μg kg(bodyweight) −1 d −1 ) glyphosate contaminations 10 . Next to the detailed 1 screening of animal health, especially the question whether glyphosate detrimentally influences the ruminal microbiome in an interactive manner with the dietary composition was of interest. To address this question, rumen fluid samples were collected throughout the trial and subjected to 16S rDNA genotyping, analyses of bacterial endotoxin (lipopolysaccharide = LPS) content, and microbial fermentation products. Finally, a screening for abundance C. botulinum neurotoxin (BoNT) genes was conducted at the beginning and at the end of the trial in order to address the question of whether glyphosate influences a putative enrichment of this pathogen in the context of challenging CFP and varying amounts of dietary fiber.

RESULTS
Analysis of microbial compounds and products in the rumen Throughout the trial, ruminal pH was altered, which was reflected by a significant time effect (p t < 0.01, Fig. 1a). Ruminal ammonia was significantly influenced by CFP (p CFP < 0.05) and changes over time differed between groups with different CFPs (p CFP*t < 0.05, Fig. 1b). The ruminal LPS concentration was significantly influenced by time (p t < 0.01) and CFP (p CFP <0.01). Changes in LPS concentrations over time also clearly differed between groups with different CFPs (p CFP*t < 0.01; Fig. 1c). All three parameters were not significantly influenced by glyphosate effects. Total short-chain fatty acid (SCFA) concentrations in the rumen were significantly varying over time (p t < 0.01) and differed significantly between CON and GLY groups (p GLY < 0.05, Fig. 1d), while no significant differences in variations over time were observed between experimental groups. For the relative proportions of the major SCFAs acetate, propionate, and butyrate, a significant effect of CFP (p CFP < 0.01) was observable and changes over time significantly differed between groups with different CFPs (p CFP*t < 0.01; Fig. 1e-g). Furthermore, a significant alteration over time on the proportions of acetate and propionate was observable (p t < 0.01). In contrast to this, no effect of glyphosate on major SCFA proportions was observable. Similarly, the low relative proportions of valerate and the branched-chain SCFA isobutyrate were influenced by CFP (p CFP < 0.01) and time (p t < 0.01) changes in proportions over time differed between groups with different CFPs (p CFP*t < 0.01, Supplementary Fig. 1a, b). Finally, for the branched-chain SCFA isovalerate, significant influences of glyphosate (p GLY < 0.05) and time (p t < 0.01) were observed, while alterations over time differed between GLY and CON groups (p GLY*t < 0.05; Supplementary Fig. 1c).

Sequencing and data preparation
To evaluate the influences of GLY and CFP on the ruminal microbiome, a 16S rDNA genotyping was conducted. The individual sequencing libraries had an average size of 178,330 reads. After quality filtering an average of 117,601 reads per library remained. Normalization for the size of the smallest library retained 28% of the total initially sequenced reads and 18,518 unique sequences, which were represented as 50,546 remaining reads per individual library (see Supplementary Table 1).

Rumen alpha-diversity
In order to compare animal-individual microbiome compositions, analyses of microbial alpha-diversity, evenness, and dominance were conducted. As measurements for diversity, the number of detected operational taxonomic units (OTUs), the Shannon index, and phylogenetic diversity based on Faith's metric were significantly influenced by CFP (p CFP < 0.01), time (p t < 0.05), and changes over time significantly differed between the groups with different CFPs (p CFP*t < 0.01; Fig. 2a-c). Pielou evenness was similarly influenced by CFP (p CFP < 0.01), time (p t < 0.05), and displayed different alterations for groups with different CFPs over time (p CFP*t < 0.01; Fig. 2d). Berger-Parker dominance was influenced by CFP (p CFP < 0.05) and variations over time differed between groups with different CFP (p CFP*t < 0.05; Fig. 2e). Glyphosate did not influence any of these measures.

Rumen beta-diversity
To assess whether major differences in microbial composition between animals were influenced by glyphosate exposure or different CFPs, a distance matrix based on weighted UniFrac distance was calculated and subjected to an Adonis test (Table 1) 37 . In the Adonis test, significant influences of CFP (p < 0.01) and of an interaction between CFP and time (p = 0.01) on distance between samples were observed, while glyphosate and its interactions with other experimental factors did not significantly influence the distance between samples (Table 1). However, only a limited proportion of the distances between samples could be explained through CFP or its interaction with time (R² = 0.10 for CFP and R² = 0.05 for CFP*t). The observed influence of CFP on distance between samples is supported by the results from permANOVA, permdisp, and anosim. In permANOVA and anosim, differences between groups with different CFP were significantly higher than differences within groups (p < 0.01 for both tests, Supplementary Table 2), while differences between groups were not significantly influenced by glyphosate (p > 0.05, Supplementary Table 2). Significant effects caused through homogeneity of multivariate dispersions were not observed in permdisp (p > 0.05, Supplementary Table 2) for any of the significant group differences detected in permANOVA. In order to visualize distances in microbiome composition between the individual samples, a principal coordinate analysis (PCoA) based on the weighted UniFrac distance matrix was conducted. Assigning color and shape-size to experimental groups made apparent, that animals with different CFPs had bigger distances between each other than animals fed with the same CFP (Fig. 3a). Presence of glyphosate contaminations in the diet in contrast apparently did not influence the distance between samples (Fig. 3a). Coloration by average daily glyphosate exposure in weeks 8 and 16 of the experiment did not indicate clustering of samples based on the amount of exposure (Fig. 3b). Coloration based on total ruminal SCFA concentration resulted in a color gradient in the visualization, which indicated smaller distances of the ruminal microbiome of animals with similar concentrations of total SCFAs than of animals with different concentrations (Fig. 3c). Coloration by ruminal LPS concentration indicated, that the majority of animals had low ruminal LPS concentrations and that most animals with a high ruminal LPS concentration were located in one region in the PCoA (Fig. 3d). Comparison to the experimental groups ( Fig. 3a) indicated, that high ruminal LPS content was coupled to a high CFP and independent of glyphosate exposure. When using colorations based on the major ruminal SCFA proportions, color gradients indicated, that the ruminal microbiomes of animals with similar proportions of SCFAs hat lower distances to each other than the ruminal microbiomes of animals with different proportions of SCFAs ( Fig. 3e-g). A comparison to the experimental groups (Fig. 3a) indicated a glyphosate-independent CFP-dependency of these gradients.

Taxonomic classification
In addition to diversity measures, the taxonomic composition of the rumen microbiota was analyzed. In all experimental groups, the majority of ruminal bacteria were classified as members of the phyla Bacteroidetes and Firmicutes, while Tenericutes, Actinobacteria, and Spirochaetes also showed relatively high abundances. Rare phyla (<1%) and unclassified reads made up for 3.6% (GLY HC in week 8) to 6.2% (GLY LC in week 8) of the total sequenced reads. At the beginning of the trial, the ratio between Bacteroidetes and Firmicutes was at 1.5 for the pooled herd. It stayed in a range between 1.3 and 1.5 for the LC groups throughout the trial, while for the HC groups a decrease to values between 0.9 and 1.1 was observed (Fig. 4). The phylum Bacteroidetes, was completely assigned to the order Bacteroidales and the Prevotellaceae were the dominant family within this order and also within the complete microbiome. Average proportions of total classified reads per group that were assigned to the Prevotellaceae ranged from 38.7% for the GLY HC group in week 16 to 45.7% for the CON LC group in week 16 (Fig. 4  Prevotella 7 was associated with a high CFP, the other genera were visually associated with a low CFP. Further Bacteroidales belonged to uncultured groups from the Rikenellaceae, to the uncultured Bacteroidales F082 and RF16 groups or were considered rare taxa (<1% in all conditions tested). Like Prevotella 1, these taxa were visually associated with a low CFP (Fig. 4, Supplementary Data 1). The Firmicutes were represented by taxa from the order Clostridiales, which accounted for 29.4% to 34.7% of total classified reads and the Erysipelochtrichales, Selenomonadales, and Lactobacillales, which each had high abundance (>1%) in at least one group (Fig. 4, Supplementary Data 1). In groups with a high CFP the proportion of Erysipelochtrichales, represented by uncultured Erysiletrichaceae and Sharpea, was strongly increased in comparison to groups with a low CFP (2.2-3.5% in LC; 9.6-12.1% in HC), while the proportion of Clostridiales did not show a strong response with respect to the high abundance. Within the Clostridiales, however, an enrichment of Lachnospiraceae and a    Association of microbiome data with glyphosate and CFP In order to find microbial taxons that are associated with varying CFPs and glyphosate contaminations in the diet, data of all experimental groups were pooled for week 0, while for weeks 8 and 16 data were pooled according to diet irrespective of time in the trial resulting in groups Week 0, CON HC, CON LC, GLY HC, and GLY LC. The pooled dataset was subjected to PLS-DA and the groups formed three clusters based on the CFP: Week 0 (40% CFP), GLY HC together with CON HC (60% CFP), and GLY LC together with CON LC (30% CFP), while overlaps of GLY groups with their respective CON groups were observed (Fig. 5a). For comparison, separations based on the environmental factor temperaturehumidity index (THI, Fig. 5b) or based on the animal-specific factor year of birth ( Fig. 5c) could be observed in separate PLS-DA analyses. When having a look at the microbial taxons relevant for separation in the PLS-DA analysis based on CFP and glyphosate contaminations (VIP score>1 in at least one of 9 components, Supplementary Data 2), significant changes in the abundance of 9 out of the top 10 microbial taxons (abundance >0.75% in at least one group) were associated with CFP or changes in abundance over time differed between groups with different CFP (Fig. 6a).
Ruminococcaceae UCG 14, the remaining genus in the top 10 microbial taxons, showed significant changes over time independent of glyphosate or CFP (p t < 0.01; Fig. 6a). Differing changes over time between GLY and CON groups or significant influences of glyphosate within the top 10 microbial taxons were observed for Lactobacillus (p GLY < 0.01; p GLY*t < 0.01) and Prevotella 9 (p GLY*t < 0.05; Fig. 6a). When considering further microbial taxons with VIP scores >1 and with an abundance >0.2% in at least one experimental group, 12 out of 14 taxons were significantly influenced by CFP or abundances changed over time in a CFPdependent manner. Among these CFP-responsive taxons, Weissella was influenced by glyphosate (p GLY < 0.05) and changes over time differed between CON and GLY groups (P GLY*t < 0.01). The remaining two taxons were not influenced by CFP. Kandleria was not influenced by any experimental factor, while the Eubacterium ventriosum group was significantly influenced by time and glyphosate (Fig. 6b, Supplementary Table 3). In summary, 21 of 24 microbial taxons relevant for separation in PLS-DA based on experimental factors were significantly influenced by CFP or changes over time were differing based on CFP (Fig. 6a,   Supplementary Table 3). One taxon was not significantly influenced by any experimental factor and abundances of Lactobacillus, Prevotella 9, Weissella, and the Eubacterium ventriosum group were influenced by glyphosate or combinations of glyphosate and further experimental factors (Fig. 6a, b, Supplementary Table 3).
Association of microbiome data with feed and rumen data To get further insights, whether exposure to glyphosate and intake of dietary components correlate with abundance of specific rumen microbial taxons and in how far abundance of specific taxons correlates with rumen microbial products, PLS analyses were conducted. The correlation coefficients between glyphosate exposure and abundance of microbial taxons were weak (−0.3 < r < 0.3) and did not influence clustering, while correlation coefficients for intake of other feed components led to clustering of ruminal bacterial taxons into four major clusters (Fig. 7a). The first cluster, represented by Ruminococcus, Fibrobacter, and Treponema, showed positive correlations with dietary fiber and mineral intake, while displaying negative correlations with intake of dietary sugar, starch, fat, and protein (Fig. 7a, cluster I). A second smaller cluster, represented by Acetobacter, showed correlations with dietary fiber and sugar intake, while negatively correlating with dietary mineral intake (Fig. 7a, cluster II). The third cluster including different Prevotellaceae and Succiniclasticum was negatively correlated to sugar, starch, fat, and protein intake (Fig. 7a, cluster III) and a fourth cluster, represented by Weissella and Catenisphaera was negatively correlated to intake of dietary fiber components, while being positively correlated to intake of dietary starch (Fig. 7a, cluster IV). The PLS analysis of rumen microbial products and bacterial taxons in the rumen led to clustering of ruminal microbes into 4 main clusters (Fig. 7b). In general, genera like Ruminococcus, that were associated with dietary fiber intake, were positively correlated with ruminal acetate and butyrate proportions, while genera like Weissella, positively correlated to dietary starch intake, were also positively correlated to ruminal propionate and valerate proportions (Fig. 7a,  b). Interestingly, Weissella also displayed a positive correlation to ruminal LPS, while belonging to the Gram-positive bacteria. While correlations of microbial genera with feed component intake and microbial products was relatively strong (r up to 0.8 and −0.8), correlations with glyphosate exposure, as mentioned, ruminal pH, pH of feces, and ruminal ammonia were comparably weak (−0.3 ≤ r ≤ 0.3).

Detection of C. botulinum neurotoxin gene presence
Since precise detection of pathogenic Clostridia via 16S rRNA genotyping is not possible due to sequence identities with nonpathogenic Clostridia, detection of C. botulinum neurotoxin (BoNT) genes in feces and rumen fluid at the beginning and at the end of the trial was used to gain more detailed information about these microorganisms. At the beginning of the trial, before glyphosate exposure, bont/A-presence was detected in 6.4% of all animals (n = 4), which is within the borders of average detection among German dairy cattle herds (Fig. 8) 36 . These detections were spread among the groups CON LC (n = 1), which was the only sample with detection in rumen fluid, CON HC (n = 2) and GLY HC (n = 1). At the end of the trial in week 16 bont/E-presence was detected in 11.5% of all animals (n = 7), which is slightly above average detection rate for bont-positive animals in German dairy cattle herds according to literature (Fig. 8) 36 . Detection was distributed over feces samples of all groups and detection rate was highest in CON HC (n = 3) followed by GLY LC (n = 2), CON LC (n = 1), and GLY HC (n = 1).  In general, numbers of bont-positive animals were too low to draw statistically valid conclusions. However, all animals remained healthy throughout the trial and no signs of intestinal neurotoxinrelated effects were observed.

DISCUSSION
Glyphosate has been suspected to have deleterious effects on specific microorganisms in the ruminal microbiome possibly leading to pathogen enrichment and to an inhibited ruminal fiber degradation 27,28 . Moreover, it has been proposed that these effects introduce a dysbiotic state to the ruminal microbiome in interaction with the feed composition 28 . Accordingly, putative disturbances to the rumen microbiome introduced by glyphosate in this study will be discussed in the context of varying CFP and dietary fiber content. On the one hand, possible CFP-related and dietary fiber-related glyphosate effects will be evaluated, while on the other hand the diet-related effects will serve as a comparison and point out the validity of the results of this study. SCFAs, ammonia, pH, and LPS are determinants of microbial metabolism and microbial composition in the rumen, and their concentrations are directly linked to the composition of the diet. High CFP diets are linked to rapid utilization of easy-fermentable carbohydrates by the ruminal microbiome, which can lead to temporary drops in ruminal pH based on acidic fermentation products 38 . These drops can cause subacute ruminal acidosis (SARA), whose susceptibility is varying among individual animals 39 . SARA challenges the ruminal microbiome and can induce lysis of Gram-negative bacteria for example, which in turn is reflected in LPS in the rumen fluid 40,41 . The data in this study suggests increased LPS in rumen fluid in the groups fed with a high CFP, which could be attributed to individual animals with strongly increased LPS concentrations. The absence of changes in ruminal pH in comparison to an expected CFP-induced pH decrease can be explained by the fact that samples were taken in the morning when the cows had a period without fresh feed due to milking and additionally remained without food supply throughout sample collection 10 . Accordingly, rapid drops in pH were not detected. However, pH levels and the number of animals with marked increases in LPS were evenly spread among CON and GLY groups, which indicates that glyphosate had no influence on a putative lysis of ruminal Gram-negative bacteria or ruminal pH in general independent of dietary fiber content or challenge by a high CFP. Remarkably, total ruminal SCFA concentrations seemed to be dropping in the GLY LC group after glyphosate exposure. However, when taking the non-normalized data into account, it becomes apparent that the observed effect is biased by significantly higher total ruminal SCFA concentrations in the group GLY LC (mean 56.2 mmol L −1 ) in comparison to the other groups (pooled mean 45.2 mmol L −1 ) prior to the trial. The stronger decrease in total ruminal SCFA concentrations in this group rather reflects an adjustment to the total SCFA levels of the other groups than a glyphosate effect (mean 35.3 mmol L −1 for GLY LC in week 16, mean 34.3 mmol L −1 for the pooled other groups in week 16). The proportions of SCFAs were displaying a pattern related to dietary fiber and CFP and were not responsive to glyphosate, which was expected based on studies in sheep and in vitro 8,29 . Isovalerate was seemingly influenced by glyphosate and displayed lower proportions in the GLY groups, which is contradicting with existing literature, where isovalerate proportions were increasing upon glyphosate exposure in an in vitro rumen simulation 29 . The isovalerate proportion was highest for group CON LC in week 16 (average 1.79%), while the other groups had lower proportions (average 1.29%). However, isovalerate has a small proportion among total SCFAs which is linked to changes in the major SCFAs and observed variations are within a range that was observed for animals from this herd independent of experimental conditions previously 42,43   significances based on these variations are most likely incidental findings and not reproducible. The increased amount of ruminal ammonia in the groups with a low CFP displays the balance of increased crude protein content in the trial diets, differences in urea recycling, and decreased energy supply inhibiting microbial metabolism 10 . It has been proposed that detrimental glyphosate effects on microbiomes rather occur in fiber-rich diets based on in vitro data 28 . However, neither the data of the alpha-diversity nor beta-diversity analyses indicated an effect of glyphosate on microbial diversity in the rumen. In contrast to this, reduced alphadiversity, evenness, and an increased dominance were related to high CFP. Likewise, in this trial animals fed with an elevated CFP also significantly differed from animals with a low CFP in betadiversity analyses. The composition of the main taxons in the analyzed rumen microbiomes in this study included Bacteroidetes, Firmicutes, Tenericutes, and Spirochaetes on phylum level and Prevotella, Ruminococcus, Butyvibrio, Shuttleworthia, and Clostridium on genus level, which is in agreement with microbiome analyses performed on the herd of the animals in this trial previously and with literature on core rumen microbiomes of dairy cattle 44,45 . A decreased proportion of Bacteroidetes was associated with increased CFP, the proportion of Firmicutes was increased upon high CFP, while no effects of glyphosate were observable. However, a detailed analysis of the data on genus level indicated, that a large proportion of the identified microbial taxons have not been cultured and identification was based on data from previous 16S rDNA genotyping studies of ruminal microbiota. Furthermore, discrepancies between genotype and phenotype have been discussed in the context of ruminal microbiota 44 . Accordingly, it was proposed previously, that characterization of ruminal microbiome should focus on functional properties rather than on genotyping based on 16S rDNA genotyping 44,46 . Despite these challenges related to discrepancies between genotyping and existing phenotypes in the analyzed microbiomes, PLS-DA pointed out a complete overlap of GLY and CON treatment groups based on their CFP. In comparison, seasonal changes in environmental temperature or differences in the year of birth could be detected by PLS-DA. This confirms the general observation, that the ruminal microbiome is shaped by various factors and strongly indicates, that independent of dietary fiber content or challenging conditions for the ruminal microbiome by a high CFP, detrimental influences of glyphosate on the microbiome as a functional unit are most likely negligible in vivo. When considering individual microbial genera, as expected, most relevant microbial taxons responsible for group separation in PLS-DA were statistically associated with the CFP in the diet, while changes related to glyphosate were observed only for 4 taxons. However, when analyzing the abundances outside of the statistical model, it becomes apparent, that glyphosate effects for Weissella and Lactobacillus should be attributed to an increase of Lactobacillaceae in the CON HC group in week 8 rather than specifically to the glyphosate treatment and might reflect short-term drops in ruminal pH 47 . Also, for the taxons from the Eubacterium ventriosum group and Prevotella 9, glyphosate-based effects are possibly statistical artifacts, as abundances clearly follow comparable CFP-related patterns in GLY and in CON groups. Moreover, a physiological relevance is unlikely due to limited abundance of these taxons. In agreement with observed separation based on CFP in PLS-DA, correlations of microbial taxons with intake of specific feed components and concentrations of specific ruminal metabolites could be observed. For example, Weissella abundance was positively correlated to intake of starch and sugar, as well as to the concentrations of propionate, valerate, and LPS in the rumen. Weissella belongs to the hetero-fermentative lactic acid bacteria and produces lactic acid and acetate through rapid fermentation 48 . Accordingly, it seems logical that Weissella occurs in positive correlation with LPS, since rapid fermentation is linked to rumen acidification and consequently lysis of bacteria, as discussed above. Such CFP-dependent effects were observed for multiple bacterial genera in the rumen and were coherent in the context of rumen metabolism. In contrast, for glyphosate, no coherent effects on bacterial genera could be observed and correlations between glyphosate exposure and abundance of bacterial genera remained elusive (−0.3 < r < 0.3). The detection of BoNT gene presence in the fecal and rumen fluid samples was limited to one or two animals per group and sampling and there were no clear indications for influences by CFP or GLY. The proportions of bont-positive animals on herd level were reflecting expectations based on literature 36 . As no signs of disease were observed and no increase of bont-positive animals was observed despite a highly increased glyphosate exposure in the GLY groups in comparison to the average daily exposure of cattle in Germany, a glyphosate-based effect on pathogenic Clostridia can be rejected for this study despite a relatively low number of animals used in the trial. When summarizing the results of the analyses of rumen microbial products, the composition of the ruminal microbiome, the associations of rumen microbial taxa with feed and microbial products and finally the detection of fecal pathogenic Clostridia it becomes obvious, that the principal diet composition unsurprisingly had a strong impact on the ruminal microbiome and its fermentation products, while glyphosate apparently had no relevant effects. This is in agreement with in vitro observations using mixed cultures of microorganisms from ruminal fluid showing that glyphosate had no adverse effects on simulated ruminal microbiomes 29,30 . The absence of long-term glyphosate effects in the gastrointestinal tract is not surprising, since excess of amino acids is available and synthesis of aromatic amino acids can possibly be bypassed through transport 19 . However, it was pointed out recently, that Salmonella isolates from before broad usage of glyphosate are more sensitive toward glyphosate than modern isolates 49 . Given that typical ruminal bacteria have been exposed to low concentrations of glyphosate on regular basis, it might be possible that a common ruminal microbiome is adapted to glyphosate exposure. Consequently, under exposure to "reallife" scenario concentrations of glyphosate, effects might not occur in vivo due to previous adaptations, while studies using approaches with highly elevated glyphosate concentrations in vitro could still display effects due to the nature of study design. This, however, is primarily of academic interest, since in reality most cattle herds worldwide would be carrying adapted microbiomes and are not be exposed to highly elevated ruminal glyphosate concentrations used in in vitro trials.

Ethical statement
The experiment was accomplished in accordance with the German Animal Welfare Act approved by the LAVES (Lower Saxony State Office for Consumer Protection and Food Safety, Germany) at the experimental station of the Institute of Animal Nutrition, Friedrich-Loeffler-Institut (FLI), Braunschweig, Germany 10 .

Experimental feed and set-up
The feed for this trial was produced at the experimental station of the Institute of Animal Nutrition, Friedrich-Loeffler-Institut (FLI), Braunschweig, Germany 10 . In short, subsets of cultivated plants intended for feedstuff production were treated with glyphosate pre-harvest according to the legal regulations. Afterward total mixed rations (TMR) with high (60% CFP, 14% crude fiber (=HC)) and low (30% CFP, 21% crude fiber (=LC)) CFP were prepared as glyphosate contaminated (GLY HC, GLY LC) and glyphosate free (CON HC, CON LC) rations 10 . Sixty-four Holstein cows were assigned to the resulting four groups considering an even distribution of age, performance, and body condition, and 61 of these cows were used for data analyses in the trial 10 . Two cows were removed from the trial due to non-glyphosate-related diseases and suffered from abomasal displacement or general peritonitis. Another cow became dry in trial week 11 and was excluded from the trial 10 . Cows were kept in separated free stall-barns and fed the experimental diets over a period of 16 weeks.

Sample collection
All rumen fluid samples were collected after milking in the morning at the beginning of the trial and after 8 and 16 weeks of feeding the experimental diets by using the oro-ruminal probe 50 . The procedure included discarding~200-300 mL of rumen fluid after initial suction before the sample of interest was collected to minimize saliva contamination. After collection, samples for sequencing and measurements of ruminal SCFAs were stored at −20°C until further processing. Samples for measurements of ammonia and ruminal LPS were centrifuged and heat treated before storage at −80°C. Further animal parameters were documented as described 10 .

Measurement of ruminal LPS content
LPS concentration in rumen fluid samples was determined using the Limulus amebocyte lysate assay (Kinetic-QCL, Lonza, Walkersville, MD, USA) 51 .

Detection of C. botulinum neurotoxin genes
The detection of C. botulinum neurotoxin genes and expression in rumen fluid and fecal samples was conducted by PCR-based methods 36 . In short PCR-based detection of C. botulinum neurotoxin genes A to F was performed by PCR and supported by real-time-PCR-based detection 54,55 .

DNA extraction and 16S rDNA sequencing
Rumen fluid samples were thawed at 4°C overnight and 25 ml sample material was centrifuged for 5 min to remove feed debris and protozoans. Supernatant was collected and processed using a Mikro-Dismembrator S (Sartorius, Goettingen, Germany) for 2 min at a frequency of 2000 min −1 in vials cooled by liquid nitrogen. DNA was isolated and purified using the QIAmp DNA Mini Kit (Quiagen, Hilden, Germany; following the manufacturer's instructions). 16S rDNA was amplified using HotStarTaq DNA polymerase (Qiagen, Hilden, Germany; following the manufacturer's instructions) and derivatives of primers Com1/Com2-ph with barcode and adapter sequences for IonTorrent PGM (Thermo Fisher Scientific, Waltham, MA, USA, primers see Supplementary

Data preparation
After adapter removal and demultiplexing were performed with 454 GENOME SEQUENCER SOFTWARE SUITE 2.6 (Roche, Mannheim, Germany), sequences were imported into QIIME2 37 . Quality filtering, chimera removal, and denoising were applied by usage of the DADA2 plugin with parameters for trimming 15 nt from the 5′ end and a cutoff length of 265 nt resulting in 250 nt read fragments. Libraries were rarefied to the size of the smallest library.

Beta-diversity analysis
Beta-diversity analysis was conducted using the QIIME2 diversity plugin using weighted UniFrac as distance metric 61 . Group differences were visualized using PCoA and the QIIME2 emperor plugin. Tests for significant group differences were conducted using the QIIME2 plugins ADONIS, PERMANOVA, PERMDISP, and ANOSIM. Differences were considered significant with p < 0.05 and highly significant with p < 0.01.

Taxonomic classification
Taxonomic classification was performed against full-length SILVA SSU 132 REF NR 99 (release 132) using the consensus-BLAST plugin in QIIME2 with default parameters 62 . Donut plot visualizations of relative proportions were created using RSTUDIO

Statistical analyses
Unless explicitly stated otherwise, data were collected for each animal in each group (CON HC n = 16, CON LC n = 16, GLY HC n = 15, GLY LC n = 14) at every time point (week 0, week 8, and week 16) and considered for the analyses. For statistical analyses of ruminal pH, ammonia concentrations, total ruminal SCFA concentrations, ruminal SCFA proportions, ruminal LPS, alpha-diversity measures, and abundance of rumen microbial taxons, data were normalized for mean values over all groups in week 0 to account for differences before the beginning of the trial and afterward analyzed with linear mixed effect models using RSTUDIO version 1.1.456 (RSTUDIO Team, 2016) with package NLME. Significance was tested via one-tailed ANOVA. Time (t; weeks of experiment), treatment (GLY or CON), and CFP (HC or LC diet), as well as their interactions (CFP*t, CFP*GLY, GLY*t, CFP*GLY*t) were applied as fixed factors with random intercepts for individual animals. For each variable, the covariance structure was chosen based on the smallest Akaike information criterion 63 . Effects were considered significant with p < 0.05 and highly significant with p < 0.01. Relevant rumen microbial taxons for group differences were determined via sPLS-DA using the RSTUDIO  Taxons with a VIP score >1 in at least one projection component and with an abundance >0.2% of total reads were considered relevant for detailed analyses. Correlations between feed intake and rumen microbial taxons, as well as between rumen microbial taxons and rumen microbial products were determined using sPLS from the R package MIXOMICS with 50-fold validation. For correlations, a robust Pearson's-like −0.3 > r > 0.3 was considered in order to retain glyphosate exposure in the dataset 65 .