Including microbiome information in a multi-trait genomic evaluation: a case study on longitudinal growth performance in beef cattle.

BACKGROUND: Growth rate is an important component of feed conversion efficiency in cattle and varies across the different stages of the finishing period. The metabolic effect of the rumen microbiome is essential for cattle growth, and investigating the genomic and microbial factors that underlie this temporal variation can help maximize feed conversion efficiency at each growth stage. RESULTS: By analysing longitudinal body weights during the finishing period and genomic and metagenomic data from 359 beef cattle, our study demonstrates that the influence of the host genome on the functional rumen microbiome contributes to the temporal variation in average daily gain (ADG) in different months (ADG1, ADG2, ADG3, ADG4). Five hundred and thirty-three additive log-ratio transformed microbial genes (alr-MG) had non-zero genomic correlations (rg) with at least one ADG-trait (ranging from |0.21| to |0.42|). Only a few alr-MG correlated with more than one ADG-trait, which suggests that a differential host-microbiome determinism underlies ADG at different stages. These alr-MG were involved in ribosomal biosynthesis, energy processes, sulphur and aminoacid metabolism and transport, or lipopolysaccharide signalling, among others. We selected two alternative subsets of 32 alr-MG that had a non-uniform or a uniform rg sign with all the ADG-traits, regardless of the rg magnitude, and used them to develop a microbiome-driven breeding strategy based on alr-MG only, or combined with ADG-traits, which was aimed at shaping the rumen microbiome towards increased ADG at all finishing stages. Combining alr-MG information with ADG records increased prediction accuracy of genomic estimated breeding values (GEBV) by 11 to 22% relative to the direct breeding strategy (using ADG-traits only), whereas using microbiome information, only, achieved lower accuracies (from 7 to 41%). Predicted selection responses varied consistently with accuracies. Restricting alr-MG based on their rg sign (uniform subset) did not yield a gain in the predicted response compared to the non-uniform subset, which is explained by the absence of alr-MG showing non-zero rg at least with more than one of the ADG-traits. CONCLUSIONS: Our work sheds light on the role of the microbial metabolism in the growth trajectory of beef cattle at the genomic level and provides insights into the potential benefits of using microbiome information in future genomic breeding programs to accurately estimate GEBV and increase ADG at each finishing stage in beef cattle.

Identification of intestinal and fecal microbial biomarkers using a porcine social stress model.

Understanding the relationships between social stress and the gastrointestinal microbiota, and how they influence host health and performance is expected to have many scientific and commercial implementations in different species, including identification and improvement of challenges to animal welfare and health. In particular, the study of the stress impact on the gastrointestinal microbiota of pigs may be of interest as a model for human health. A porcine stress model based on repeated regrouping and reduced space allowance during the last 4 weeks of the finishing period was developed to identify stress-induced changes in the gut microbiome composition. The application of the porcine stress model resulted in a significant increase in salivary cortisol concentration over the course of the trial and decreased growth performance and appetite. The applied social stress resulted in 32 bacteria being either enriched (13) or depleted (19) in the intestine and feces. Fecal samples showed a greater number of microbial genera influenced by stress than caecum or colon samples. Our trial revealed that the opportunistic pathogens Treponema and Clostridium were enriched in colonic and fecal samples from stressed pigs. Additionally, genera such as Streptococcus, Parabacteroides, Desulfovibrio, Terrisporobacter, Marvinbryantia, and Romboutsia were found to be enriched in response to social stress. In contrast, the genera Prevotella, Faecalibacterium, Butyricicoccus, Dialister, Alloprevotella, Megasphaera, and Mitsuokella were depleted. These depleted bacteria are of great interest because they synthesize metabolites [e.g., short-chain fatty acids (SCFA), in particular, butyrate] showing beneficial health benefits due to inhibitory effects on pathogenic bacteria in different animal species. Of particular interest are Dialister and Faecalibacterium, as their depletion was identified in a human study to be associated with inferior quality of life and depression. We also revealed that some pigs were more susceptible to pathogens as indicated by large enrichments of opportunistic pathogens of Clostridium, Treponema, Streptococcus and Campylobacter. Generally, our results provide further evidence for the microbiota-gut-brain axis as indicated by an increase in cortisol concentration due to social stress regulated by the hypothalamic-pituitary-adrenal axis, and a change in microbiota composition, particularly of bacteria known to be associated with pathogenicity and mental health diseases.

KOunt: a reproducible KEGG orthologue abundance workflow.

SUMMARY: Accurate gene prediction is essential for successful metagenome analysis. We present KOunt, a Snakemake pipeline, that precisely quantifies KEGG orthologue abundance. AVAILABILITY AND IMPLEMENTATION: KOunt is available on GitHub: https://github.com/WatsonLab/KOunt. The KOunt reference database is available on figshare: https://doi.org/10.6084/m9.figshare.21269715. Test data are available at https://doi.org/10.6084/m9.figshare.22250152 and version 1.2.0 of KOunt at https://doi.org/10.6084/m9.figshare.23607834.

Different microbial genera drive methane emissions in beef cattle fed with two extreme diets.

The ratio of forage to concentrate in cattle feeding has a major influence on the composition of the microbiota in the rumen and on the mass of methane produced. Using methane measurements and microbiota data from 26 cattle we aimed to investigate the relationships between microbial relative abundances and methane emissions, and identify potential biomarkers, in animals fed two extreme diets - a poor quality fresh cut grass diet (GRASS) or a high concentrate total mixed ration (TMR). Direct comparisons of the effects of such extreme diets on the composition of rumen microbiota have rarely been studied. Data were analyzed considering their multivariate and compositional nature. Diet had a relevant effect on methane yield of +10.6 g of methane/kg of dry matter intake for GRASS with respect to TMR, and on the centered log-ratio transformed abundance of 22 microbial genera. When predicting methane yield based on the abundance of 28 and 25 selected microbial genera in GRASS and TMR, respectively, we achieved cross-validation prediction accuracies of 66.5 ± 9% and 85 ± 8%. Only the abundance of Fibrobacter had a consistent negative association with methane yield in both diets, whereas most microbial genera were associated with methane yield in only one of the two diets. This study highlights the stark contrast in the microbiota controlling methane yield between animals fed a high concentrate diet, such as that found on intensive finishing units, and a low-quality grass forage that is often found in extensive grazing systems. This contrast must be taken into consideration when developing strategies to reduce methane emissions by manipulation of the rumen microbial composition.

Microbiome-driven breeding strategy potentially improves beef fatty acid profile benefiting human health and reduces methane emissions.

BACKGROUND: Healthier ruminant products can be achieved by adequate manipulation of the rumen microbiota to increase the flux of beneficial fatty acids reaching host tissues. Genomic selection to modify the microbiome function provides a permanent and accumulative solution, which may have also favourable consequences in other traits of interest (e.g. methane emissions). Possibly due to a lack of data, this strategy has never been explored. RESULTS: This study provides a comprehensive identification of ruminal microbial mechanisms under host genomic influence that directly or indirectly affect the content of unsaturated fatty acids in beef associated with human dietary health benefits C18:3n-3, C20:5n-3, C22:5n-3, C22:6n-3 or cis-9, trans-11 C18:2 and trans-11 C18:1 in relation to hypercholesterolemic saturated fatty acids C12:0, C14:0 and C16:0, referred to as N3 and CLA indices. We first identified that ~27.6% (1002/3633) of the functional core additive log-ratio transformed microbial gene abundances (alr-MG) in the rumen were at least moderately host-genomically influenced (HGFC). Of these, 372 alr-MG were host-genomically correlated with the N3 index (n=290), CLA index (n=66) or with both (n=16), indicating that the HGFC influence on beef fatty acid composition is much more complex than the direct regulation of microbial lipolysis and biohydrogenation of dietary lipids and that N3 index variation is more strongly subjected to variations in the HGFC than CLA. Of these 372 alr-MG, 110 were correlated with the N3 and/or CLA index in the same direction, suggesting the opportunity for enhancement of both indices simultaneously through a microbiome-driven breeding strategy. These microbial genes were involved in microbial protein synthesis (aroF and serA), carbohydrate metabolism and transport (galT, msmX), lipopolysaccharide biosynthesis (kdsA, lpxD, lpxB), or flagellar synthesis (flgB, fliN) in certain genera within the Proteobacteria phyla (e.g. Serratia, Aeromonas). A microbiome-driven breeding strategy based on these microbial mechanisms as sole information criteria resulted in a positive selection response for both indices (1.36±0.24 and 0.79±0.21 sd of N3 and CLA indices, at 2.06 selection intensity). When evaluating the impact of our microbiome-driven breeding strategy to increase N3 and CLA indices on the environmental trait methane emissions (g/kg of dry matter intake), we obtained a correlated mitigation response of -0.41±0.12 sd. CONCLUSION: This research provides insight on the possibility of using the ruminal functional microbiome as information for host genomic selection, which could simultaneously improve several microbiome-driven traits of interest, in this study exemplified with meat quality traits and methane emissions. Video Abstract.

Bovine host genome acts on rumen microbiome function linked to methane emissions.

Our study provides substantial evidence that the host genome affects the comprehensive function of the microbiome in the rumen of bovines. Of 1,107/225/1,141 rumen microbial genera/metagenome assembled uncultured genomes (RUGs)/genes identified from whole metagenomics sequencing, 194/14/337 had significant host genomic effects (heritabilities ranging from 0.13 to 0.61), revealing that substantial variation of the microbiome is under host genomic control. We found 29/22/115 microbial genera/RUGs/genes host-genomically correlated (|0.59| to |0.93|) with emissions of the potent greenhouse gas methane (CH4), highlighting the strength of a common host genomic control of specific microbial processes and CH4. Only one of these microbial genes was directly involved in methanogenesis (cofG), whereas others were involved in providing substrates for archaea (e.g. bcd and pccB), important microbial interspecies communication mechanisms (ABC.PE.P), host-microbiome interaction (TSTA3) and genetic information processes (RP-L35). In our population, selection based on abundances of the 30 most informative microbial genes provided a mitigation potential of 17% of mean CH4 emissions per generation, which is higher than for selection based on measured CH4 using respiration chambers (13%), indicating the high potential of microbiome-driven breeding to cumulatively reduce CH4 emissions and mitigate climate change.

Concepts and Consequences of a Core Gut Microbiota for Animal Growth and Development.

Animal microbiomes are occasionally considered as an extension of host anatomy, physiology, and even their genomic architecture. Their compositions encompass variable and constant portions when examined across multiple hosts. The latter, termed the core microbiome, is viewed as more accommodated to its host environment and suggested to benefit host fitness. Nevertheless, discrepancies in its definitions, characteristics, and importance to its hosts exist across studies. We survey studies that characterize the core microbiome, detail its current definitions and available methods to identify it, and emphasize the crucial need to upgrade and standardize the methodologies among studies. We highlight ruminants as a case study and discussthe link between the core microbiome and host physiology and genetics, as well as potential factors that shape it. We conclude with main directives of action to better understand the host-core microbiome axis and acquire the necessary insights into its controlled modulation.

Compositional Data Analysis of Microbiome and Any-Omics Datasets: A Validation of the Additive Logratio Transformation.

Microbiome and omics datasets are, by their intrinsic biological nature, of high dimensionality, characterized by counts of large numbers of components (microbial genes, operational taxonomic units, RNA transcripts, etc.). These data are generally regarded as compositional since the total number of counts identified within a sample is irrelevant. The central concept in compositional data analysis is the logratio transformation, the simplest being the additive logratios with respect to a fixed reference component. A full set of additive logratios is not isometric, that is they do not reproduce the geometry of all pairwise logratios exactly, but their lack of isometry can be measured by the Procrustes correlation. The reference component can be chosen to maximize the Procrustes correlation between the additive logratio geometry and the exact logratio geometry, and for high-dimensional data there are many potential references. As a secondary criterion, minimizing the variance of the reference component's log-transformed relative abundance values makes the subsequent interpretation of the logratios even easier. On each of three high-dimensional omics datasets the additive logratio transformation was performed, using references that were identified according to the abovementioned criteria. For each dataset the compositional data structure was successfully reproduced, that is the additive logratios were very close to being isometric. The Procrustes correlations achieved for these datasets were 0.9991, 0.9974, and 0.9902, respectively. We thus demonstrate, for high-dimensional compositional data, that additive logratios can provide a valid choice as transformed variables, which (a) are subcompositionally coherent, (b) explain 100% of the total logratio variance and (c) come measurably very close to being isometric. The interpretation of additive logratios is much simpler than the complex isometric alternatives and, when the variance of the log-transformed reference is very low, it is even simpler since each additive logratio can be identified with a corresponding compositional component.

Comprehensive functional core microbiome comparison in genetically obese and lean hosts under the same environment.

Our study provides an exhaustive comparison of the microbiome core functionalities (captured by 3,936 microbial gene abundances) between hosts with divergent genotypes for intramuscular lipid deposition. After 10 generations of divergent selection for intramuscular fat in rabbits and 4.14 phenotypic standard deviations (SD) of selection response, we applied a combination of compositional and multivariate statistical techniques to identify 122 cecum microbial genes with differential abundances between the lines (ranging from -0.75 to +0.73 SD). This work elucidates that microbial biosynthesis lipopolysaccharides, peptidoglycans, lipoproteins, mucin components, and NADH reductases, amongst others, are influenced by the host genetic determination for lipid accretion in muscle. We also differentiated between host-genetically influenced microbial mechanisms regulating lipid deposition in body or intramuscular reservoirs, with only 28 out of 122 MGs commonly contributing to both. Importantly, the results of this study are of relevant interest for the efficient development of strategies fighting obesity.

Identification of Complex Rumen Microbiome Interaction Within Diverse Functional Niches as Mechanisms Affecting the Variation of Methane Emissions in Bovine.

A network analysis including relative abundances of all ruminal microbial genera (archaea, bacteria, fungi, and protists) and their genes was performed to improve our understanding of how the interactions within the ruminal microbiome affects methane emissions (CH4). Metagenomics and CH4 data were available from 63 bovines of a two-breed rotational cross, offered two basal diets. Co-abundance network analysis revealed 10 clusters of functional niches. The most abundant hydrogenotrophic Methanobacteriales with key microbial genes involved in methanogenesis occupied a different functional niche (i.e., "methanogenesis" cluster) than methylotrophic Methanomassiliicoccales (Candidatus Methanomethylophylus) and acetogens (Blautia). Fungi and protists clustered together and other plant fiber degraders like Fibrobacter occupied a seperate cluster. A Partial Least Squares analysis approach to predict CH4 variation in each cluster showed the methanogenesis cluster had the best prediction ability (57.3%). However, the most important explanatory variables in this cluster were genes involved in complex carbohydrate degradation, metabolism of sugars and amino acids and Candidatus Azobacteroides carrying nitrogen fixation genes, but not methanogenic archaea and their genes. The cluster containing Fibrobacter, isolated from other microorganisms, was positively associated with CH4 and explained 49.8% of its variability, showing fermentative advantages compared to other bacteria and fungi in providing substrates (e.g., formate) for methanogenesis. In other clusters, genes with enhancing effect on CH4 were related to lactate and butyrate (Butyrivibrio and Pseudobutyrivibrio) production and simple amino acids metabolism. In comparison, ruminal genes negatively related to CH4 were involved in carbohydrate degradation via lactate and succinate and synthesis of more complex amino acids by γ-Proteobacteria. When analyzing low- and high-methane emitters data in separate networks, competition between methanogens in the methanogenesis cluster was uncovered by a broader diversity of methanogens involved in the three methanogenesis pathways and larger interactions within and between communities in low compared to high emitters. Generally, our results suggest that differences in CH4 are mainly explained by other microbial communities and their activities rather than being only methanogens-driven. Our study provides insight into the interactions of the rumen microbial communities and their genes by uncovering functional niches affecting CH4, which will benefit the development of efficient CH4 mitigation strategies.

Correlated responses to selection for intramuscular fat in several muscles in rabbits.

The aim of this study was to evaluate the effect of six generations of selection for intramuscular fat (IMF) in muscle Longissimus thoracis et lumborum (LTL) at 9 wk. in IMF and fatty acid (FA) composition of muscles with diverse metabolic profile in rabbits. Direct response to selection was 0.33 g of IMF/ 100 g of LTL, around 0.4 SD per generation. A positive correlated response was observed in IMF of Biceps femoris, Supraspinatus and Semimembranosus proprius muscles at 9 wk., representing around 0.2 SD of the trait per generation. Selection affected similarly the FA composition of all muscles at 9 wk., high-IMF line showing greater monounsaturated but lower polyunsaturated FA percentages than low-IMF line, whereas no differences were observed for saturated FA. Traits were also measured at 13 wk. and correlated responses were in the same direction. Our results suggest a common genetic background for IMF and FA composition in muscles with different metabolic profile in rabbits.

Selection for environmental variance of litter size in rabbits.

BACKGROUND: In recent years, there has been an increasing interest in the genetic determination of environmental variance. In the case of litter size, environmental variance can be related to the capacity of animals to adapt to new environmental conditions, which can improve animal welfare. RESULTS: We developed a ten-generation divergent selection experiment on environmental variance. We selected one line of rabbits for litter size homogeneity and one line for litter size heterogeneity by measuring intra-doe phenotypic variance. We proved that environmental variance of litter size is genetically determined and can be modified by selection. Response to selection was 4.5% of the original environmental variance per generation. Litter size was consistently higher in the Low line than in the High line during the entire experiment. CONCLUSIONS: We conclude that environmental variance of litter size is genetically determined based on the results of our divergent selection experiment. This has implications for animal welfare, since animals that cope better with their environment have better welfare than more sensitive animals. We also conclude that selection for reduced environmental variance of litter size does not depress litter size.