Host and rumen microbiome contributions to feed efficiency traits in Holstein cows

It is now widely accepted that dairy cow performance is influenced by both the host genome and rumen microbiome composition. The contributions of the genome and the microbiome to the phenotypes of interest are quantified by heritability (h2) and microbiability (m2), respectively. However, if the genome and microbiome are included in the model, then the h2 reflects only the contribution of the direct genetic effects quantified as direct heritability (hd2), and the holobiont effect reflects the joint action of the genome and the microbiome, quantified as the holobiability (ho2). The objectives of this study were to estimate h2, hd2,m2, and ho2 for dry matter intake, milk energy, and residual feed intake; and to evaluate the predictive ability of different models, including genome, microbiome, and their interaction. Data consisted of feed efficiency records, SNP genotype data, and 16S rRNA rumen microbial abundances from 448 mid-lactation Holstein cows from 2 research farms. Three kernel models were fit to each trait: one with only the genomic effect (model G), one with the genomic and microbiome effects (model GM), and one with the genomic, microbiome, and interaction effects (model GMO). The model GMO, or holobiont model, showed the best goodness-of-fit. The hd2 estimates were always 10% to 15% lower than h2 estimates for all traits, suggesting a mediated genetic effect through the rumen microbiome, and m2 estimates were moderate for all traits, and up to 26% for milk energy. The ho2 was greater than the sum of hd2 and m2, suggesting that the genome-by-microbiome interaction had a sizable effect on feed efficiency. Kernel models fitting the rumen microbiome (i.e., models GM and GMO) showed larger predictive correlations and smaller prediction bias than the model G. These findings reveal a moderate contribution of the rumen microbiome to feed efficiency traits in lactating Holstein cows and strongly suggest that the rumen microbiome mediates part of the host genetic effect.

Leveraging microbiome information for animal genetic improvement

There is growing evidence that the microbiome influences host phenotypic variation. Incorporating information about the holobiont - the host and its microbiome - into genomic prediction models may accelerate genetic improvements in farmed animal populations. Importantly, these models must account for the indirect effects of the host genome on microbiome-mediated phenotypes.

Microbiability of milk composition and genetic control of microbiota effects in sheep

Recently, high-dimensional omics data are becoming available in larger quantities, and models have been developed that integrate them with genomics to understand in finer detail the relationship between genotype and phenotype, and thus improve the performance of genetic evaluations. Our objectives are to quantify the effect of the inclusion of microbiome data in the genetic evaluation for dairy traits in sheep, through the estimation of the heritability, microbiability, and how the microbiome effect on dairy traits decomposes into genetic and nongenetic parts. In this study we analyzed milk and rumen samples of 795 Lacaune dairy ewes. We included, as phenotype, dairy traits and milk fatty acids and proteins composition; as omics measurements, 16S rRNA rumen bacterial abundances; and as genotyping, 54K SNP chip for all ewes. Two nested genomic models were used: a first model to predict the individual contributions of the genetic and microbial abundances to phenotypes, and a second model to predict the additive genetic effect of the microbial community. In addition, microbiome-wide association studies for all dairy traits were applied using the 2,059 rumen bacterial abundances, and the genetic correlations between microbiome principal components and dairy traits were estimated. Results showed that in general the inclusion of both genetic and microbiome effect did not improve the fit of the model compared with the model with the genetic effect only. In addition, for all dairy traits the total heritability was equal to the direct heritability after fitting microbiota effects, due to a microbiability being almost zero for most dairy traits and heritability of the microbial community was very close to zero. Microbiome-wide association studies did not show operational taxonomic units with major effect for any of the dairy traits evaluated, and the genetic correlations between the first 5 principal components and dairy traits were low to moderate. So far, we can conclude that, using a substantial data set of 795 Lacaune dairy ewes, rumen bacterial abundances do not provide improved genetic evaluation for dairy traits in sheep.

Genetic and microbiome analysis of feed efficiency in laying hens

Improving feed efficiency is an important target for poultry breeding. Feed efficiency is affected by host genetics and the gut microbiota, but many of the mechanisms remain elusive in laying hens, especially in the late laying period. In this study, we measured feed intake, body weight, and egg mass of 714 hens from a pedigreed line from 69 to 72 wk of age and calculated the residual feed intake (RFI) and feed conversion ratio (FCR). In addition, fecal samples were also collected for 16S ribosomal RNA gene sequencing (V4 region). Genetic analysis was then conducted in DMU packages by using AI-REML with animal model. Moderate heritability estimates for FCR (h² = 0.31) and RFI (h² = 0.52) were observed, suggesting that proper selection programs can directly improve feed efficiency. Genetically, RFI was less correlated with body weight and egg mass than that of FCR. The phenotypic variance explained by gut microbial variance is defined as the microbiability (m²). The microbiability estimates for FCR (m² = 0.03) and RFI (m² = 0.16) suggested the gut microbiota was also involved in the regulation of feed efficiency. In addition, our results showed that the effect of host genetics on fecal microbiota was minor in three aspects: 1) microbial diversity indexes had low heritability estimates, and genera with heritability estimates more than 0.1 accounted for only 1.07% of the tested fecal microbiota; 2) the genetic relationship correlations between host genetics and different microbial distance were very weak, ranging from -0.0057 to -0.0003; 3) the microbial distance between different kinships showed no significant difference. Since the RFI has the highest microbiability, we further screened out three genera, including Anaerosporobacter, Candidatus Stoquefichus, and Fournierella, which were negatively correlated with RFI and played positive roles in improving the feed efficiency. These findings contribute to a great understanding of the genetic background and microbial influences on feed efficiency.

Heritable and Nonheritable Rumen Bacteria Are Associated with Different Characters of Lactation Performance of Dairy Cows

Recent studies have reported that some rumen microbes are heritable. However, it is necessary to clarify the functions and specific contributions of the heritable rumen microbes to cattle phenotypes (microbiability) in comparison with those that are nonheritable. This study aimed to identify the distribution and predicted functions of heritable and nonheritable bacterial taxa at species level in the rumen of dairy cows and their respective contributions to energy-corrected milk yield, protein content and yield, and fat content and yield in milk. Thirty-two heritable and 674 nonheritable bacterial taxa were identified at species level, and the functional analysis revealed that predicted microbial functions for both groups were mainly enriched for energy, amino acid, and ribonucleotide metabolism. The mean microbiability (to reflect a single taxon's contribution) of heritable bacteria was found to range from 0.16% to 0.33% for the different milk traits, whereas the range for nonheritable bacteria was 0.03% to 0.06%. These findings suggest a strong contribution by host genetics in shaping the rumen microbiota, which contribute significantly to milk production traits. Therefore, there is an opportunity to further improve milk production traits through attention to host genetics and the interaction with the rumen microbiota. IMPORTANCE Rumen bacteria produce volatile fatty acids which exert a far-reaching influence on hepatic metabolism, mammary gland metabolism, and animal production. In the current study, 32 heritable and 674 nonheritable bacterial taxa at species level were identified, and shown to have different microbiability (overall community contribution) and mean microbiability (the average of a single taxon's contribution) for lactation performance. The predicted functions of heritable and nonheritable bacterial taxa also differed, suggesting that targeted nutritional and genetic breeding approaches could be used to manipulate them to improve dairy cow performance.

Exploring methods to summarize gut microbiota composition for microbiability estimation and phenotypic prediction in swine

The microbial composition resemblance among individuals in a group can be summarized in a square covariance matrix and fitted in linear models. We investigated eight approaches to create the matrix that quantified the resemblance between animals based on the gut microbiota composition. We aimed to compare the performance of different methods in estimating trait microbiability and predicting growth and body composition traits in three pig breeds. This study included 651 purebred boars from either breed: Duroc (n = 205), Landrace (n = 226), and Large White (n = 220). Growth and body composition traits, including body weight (BW), ultrasound backfat thickness (BF), ultrasound loin depth (LD), and ultrasound intramuscular fat (IMF) content, were measured on live animals at the market weight (156 ± 2.5 days of age). Rectal swabs were taken from each animal at 158 ± 4 days of age and subjected to 16S rRNA gene sequencing. Eight methods were used to create the microbial similarity matrices, including four kernel functions (Linear Kernel, LK; Polynomial Kernel, PK; Gaussian Kernel, GK; Arc-cosine Kernel with one hidden layer, AK1), two dissimilarity methods (Bray-Curtis, BC; Jaccard, JA), and two ordination methods (Metric Multidimensional Scaling, MDS; Detrended Correspondence analysis, DCA). Based on the matrix used, microbiability estimates ranged from 0.07 to 0.21 and 0.12 to 0.53 for Duroc, 0.03 to 0.21 and 0.05 to 0.44 for Landrace, and 0.02 to 0.24 and 0.05 to 0.52 for Large White pigs averaged over traits in the model with sire, pen, and microbiome, and model with the only microbiome, respectively. The GK, JA, BC, and AK1 obtained greater microbiability estimates than the remaining methods across traits and breeds. Predictions were made within each breed group using four-fold cross-validation based on the relatedness of sires in each breed group. The prediction accuracy ranged from 0.03 to 0.18 for BW, 0.08 to 0.31 for BF, 0.21 to 0.48 for LD, and 0.04 to 0.16 for IMF when averaged across breeds. The BC, MDS, LK, and JA achieved better accuracy than other methods in most predictions. Overall, the PK and DCA exhibited the worst performance compared to other microbiability estimation and prediction methods. The current study shows how alternative approaches summarized the resemblance of gut microbiota composition among animals and contributed this information to variance component estimation and phenotypic prediction in swine.

Gut Microbial Composition and Predicted Functions Are Not Associated with Feather Pecking and Antagonistic Behavior in Laying Hens

BACKGROUND

Feather pecking is a well-known problem in layer flocks that causes animal welfare restrictions and contributes to economic losses. Birds' gut microbiota has been linked to feather pecking. This study aims to characterize the microbial communities of two laying hen lines divergently selected for high (HFP) and low (LFP) feather pecking and investigates if the microbiota is associated with feather pecking or agonistic behavior.

METHODS

Besides phenotyping for the behavioral traits, microbial communities from the digesta and mucosa of the ileum and caeca were investigated using target amplicon sequencing and functional predictions. Microbiability was estimated with a microbial mixed linear model.

RESULTS

Ileum digesta showed an increase in the abundance of the genus Lactobacillus in LFP, while Escherichia was abundant in HFP hens. In the caeca digesta and mucosa of the LFP line were more abundant Faecalibacterium and Blautia. Tryptophan metabolism and lysine degradation were higher in both digesta and mucosa of the HFP hens. Linear models revealed that the two lines differ significantly in all behavior traits. Microbiabilities were close to zero and not significant in both lines and for all traits.

CONCLUSIONS

Trait variation was not affected by the gut microbial composition in both selection lines.

Microbiability of meat quality and carcass composition traits in swine

The impact of gut microbiome composition was investigated at different stages of production (weaning, Mid‐test and Off‐test) on meat quality and carcass composition traits of 1,123 three‐way crossbred pigs. Data were analysed using linear mixed models which included the fixed effects of dam line, contemporary group and gender as well as the random effects of pen, animal and microbiome information at different stages. The contribution of the microbiome to all traits was prominent although it varied over time, increasing from weaning to Off‐test for most traits. Microbiability estimates of carcass composition traits were greater than that of meat quality traits. Among all of the traits analysed, belly weight (BEL) had a higher microbiability estimate (0.29 ± 0.04). Adding microbiome information did not affect the estimates of genomic heritability of meat quality traits but affected the estimates of carcass composition traits. Fat depth had a greater decrease (10%) in genomic heritability at Off‐test. High microbial correlations were found among different traits, particularly with traits related to fat deposition with a decrease in the genomic correlation up to 20% for loin weight and BEL. This suggested that genomic correlation was partially contributed by genetic similarity of microbiome composition. The results indicated that better understanding of microbial composition could aid the improvement of complex traits, particularly the carcass composition traits in swine by inclusion of microbiome information in the genetic evaluation process.

The Gut Microbial Architecture of Efficiency Traits in the Domestic Poultry Model Species Japanese Quail (Coturnix japonica) Assessed by Mixed Linear Models

It is well known that mammals and avian gut microbiota compositions are shaped by the host genomes and affect quantitative traits. The microbial architecture describes the impact of the microbiota composition on quantitative trait variation and the number and effect distribution of microbiota features. In the present study the gut microbial architecture of feed-related traits phosphorus and calcium utilization, daily gain, feed intake and feed per gain ratio in the domestic poultry model species Japanese quail were assessed by mixed linear models. The ileum microbiota composition was characterized by 16S rRNA amplicon sequencing techniques of growing individuals. The microbiability of the traits was on a similar level as the narrow sense heritability and was highly significant except for calcium utilization. The animal microbial correlation of the traits was substantial. Microbiome-wide association analyses revealed several traits associated and highly significant microbiota features, both on the bacteria genera as well as on the operational taxonomic unit level. Most features were significant for more than one trait, which explained the high microbial correlations. It can be concluded that the traits are polymicrobial determined with some microbiota features with larger effects and many with small effects. The results are important for the development of hologenomic selection schemes for feed-related traits in avian breeding programs that are targeting the host genome and the metagenome simultaneously.