The fecal microbiota of Holstein cows is heritable and genetically correlated to dairy performances.

The fecal microbiota of ruminants constitutes a diversified community that has been phenotypically associated with a variety of host phenotypes, such as production and health. To gain a better understanding of the complex and interconnected factors that drive the fecal bacterial community, we have aimed to estimate the genetic parameters of the diversity and composition of the fecal microbiota, including heritabilities, genetic correlations among taxa, and genetic correlations between fecal microbiota features and host phenotypes. To achieve this, we analyzed a large population of 1,875 Holstein cows originating from 144 French commercial herds and routinely recorded for production, somatic cell score, and fertility traits. Fecal samples were collected from the animals and subjected to 16S rRNA gene sequencing, with reads classified into Amplicon Sequence Variants (ASVs). The estimated α- and ß-diversity indices (i.e., Observed Richness, Shannon index, Bray-Curtis and Jaccard dissimilarity matrices) and the abundances of ASVs, genera, families and phyla, normalized by centered-log ratio (CLR), were considered as phenotypes. Genetic parameters were calculated using either univariate or bivariate animal models. Heritabilities estimates, ranging from 0.08 to 0.31 for taxa abundances and ß-diversity indices, highlight the influence of the host genetics on the composition of the fecal microbiota. Furthermore, genetic correlations estimated within the microbial community and between microbiota features and host traits reveal the complex networks linking all components of the fecal microbiota together and to their host, thus strengthening the holobiont concept. By estimating the heritabilities of microbiota-associated phenotypes, our study quantifies the impact of the host genetics on the fecal microbiota composition. In addition, genetic correlations between taxonomic groups and between taxa abundances and host performance suggest potential applications for selective breeding to improve host traits or promote a healthier microbiota.

Use of meta-analyses and joint analyses to select variants in whole genome sequences for genomic evaluation: An application in milk production of French dairy cattle breeds.

As a result of the 1000 Bull Genome Project, it has become possible to impute millions of variants, with many of these potentially causative for traits of interest, for thousands of animals that have been genotyped with medium-density chips. This enormous source of data opens up very interesting possibilities for the inclusion of these variants in genomic evaluations. However, for computational reasons, it is not possible to include all variants in genomic evaluation procedures. One potential approach could be to select the most relevant variants based on the results of genome-wide association studies (GWAS); however, the identification of causative mutations is still difficult with this method, partly because of weak imputation accuracy for rare variants. To address this problem, this study assesses the ability of different approaches based on multi-breed GWAS (joint and meta-analyses) to identify single-nucleotide polymorphisms (SNP) for use in genomic evaluation in the 3 main French dairy cattle breeds. A total of 6,262 Holstein bulls, 2,434 Montbéliarde bulls, and 2,175 Normande bulls with daughter yield deviations for 5 milk production traits were imputed for 27 million variants. Within-breed and joint (including all 3 breeds) GWAS were performed and 3 models of meta-analysis were tested: fixed effect, random effect, and Z-score. Comparison of the results of within- and multi-breed GWAS showed that most of the quantitative trait loci identified using within-breed approaches were also found with multi-breed methods. However, the most significant variants identified in each region differed depending on the method used. To determine which approach highlighted the most predictive SNP for each trait, we used a marker-assisted best unbiased linear prediction model to evaluate lists of SNP generated by the different GWAS methods; each list contained between 25 and 2,000 candidate variants per trait, which were identified using a single within- or multi-breed GWAS approach. Among all the multi-breed methods tested in this study, variant selection based on meta-analysis (fixed effect) resulted in the most-accurate genomic evaluation (+1 to +3 points compared with other multi-breed approaches). However, the accuracies of genomic evaluation were always better when variants were selected using the results of within-breed GWAS. As has generally been found in studies of quantitative trait loci, these results suggest that part of the genetic variance of milk production traits is breed specific in Holstein, Montbéliarde, and Normande cattle.

Genomic evaluation of regional dairy cattle breeds in single-breed and multibreed contexts.

An important prerequisite for high prediction accuracy in genomic prediction is the availability of a large training population, which allows accurate marker effect estimation. This requirement is not fulfilled in case of regional breeds with a limited number of breeding animals. We assessed the efficiency of the current French routine genomic evaluation procedure in four regional breeds (Abondance, Tarentaise, French Simmental and Vosgienne) as well as the potential benefits when the training populations consisting of males and females of these breeds are merged to form a multibreed training population. Genomic evaluation was 5-11% more accurate than a pedigree-based BLUP in three of the four breeds, while the numerically smallest breed showed a < 1% increase in accuracy. Multibreed genomic evaluation was beneficial for two breeds (Abondance and French Simmental) with maximum gains of 5 and 8% in correlation coefficients between yield deviations and genomic estimated breeding values, when compared to the single-breed genomic evaluation results. Inflation of genomic evaluation of young candidates was also reduced. Our results indicate that genomic selection can be effective in regional breeds as well. Here, we provide empirical evidence proving that genetic distance between breeds is only one of the factors affecting the efficiency of multibreed genomic evaluation.

Efficiency of multi-breed genomic selection for dairy cattle breeds with different sizes of reference population.

Single-breed genomic selection (GS) based on medium single nucleotide polymorphism (SNP) density (~50,000; 50K) is now routinely implemented in several large cattle breeds. However, building large enough reference populations remains a challenge for many medium or small breeds. The high-density BovineHD BeadChip (HD chip; Illumina Inc., San Diego, CA) containing 777,609 SNP developed in 2010 is characterized by short-distance linkage disequilibrium expected to be maintained across breeds. Therefore, combining reference populations can be envisioned. A population of 1,869 influential ancestors from 3 dairy breeds (Holstein, Montbéliarde, and Normande) was genotyped with the HD chip. Using this sample, 50K genotypes were imputed within breed to high-density genotypes, leading to a large HD reference population. This population was used to develop a multi-breed genomic evaluation. The goal of this paper was to investigate the gain of multi-breed genomic evaluation for a small breed. The advantage of using a large breed (Normande in the present study) to mimic a small breed is the large potential validation population to compare alternative genomic selection approaches more reliably. In the Normande breed, 3 training sets were defined with 1,597, 404, and 198 bulls, and a unique validation set included the 394 youngest bulls. For each training set, estimated breeding values (EBV) were computed using pedigree-based BLUP, single-breed BayesC, or multi-breed BayesC for which the reference population was formed by any of the Normande training data sets and 4,989 Holstein and 1,788 Montbéliarde bulls. Phenotypes were standardized by within-breed genetic standard deviation, the proportion of polygenic variance was set to 30%, and the estimated number of SNP with a nonzero effect was about 7,000. The 2 genomic selection (GS) approaches were performed using either the 50K or HD genotypes. The correlations between EBV and observed daughter yield deviations (DYD) were computed for 6 traits and using the different prediction approaches. Compared with pedigree-based BLUP, the average gain in accuracy with GS in small populations was 0.057 for the single-breed and 0.086 for multi-breed approach. This gain was up to 0.193 and 0.209, respectively, with the large reference population. Improvement of EBV prediction due to the multi-breed evaluation was higher for animals not closely related to the reference population. In the case of a breed with a small reference population size, the increase in correlation due to multi-breed GS was 0.141 for bulls without their sire in reference population compared with 0.016 for bulls with their sire in reference population. These results demonstrate that multi-breed GS can contribute to increase genomic evaluation accuracy in small breeds.

Characterization of bovine vaginal microbiota using 16S rRNA sequencing: associations with host fertility, longevity, health, and production.

Due to their potential impact on the host's phenotype, organ-specific microbiotas are receiving increasing attention in several animal species, including cattle. Specifically, the vaginal microbiota of ruminants is attracting growing interest, due to its predicted critical role on cows' reproductive functions in livestock contexts. Notably, fertility disorders represent a leading cause for culling, and additional research would help to fill relevant knowledge gaps. In the present study, we aimed to characterize the vaginal microbiota of a large cohort of 1171 female dairy cattle from 19 commercial herds in Northern France. Vaginal samples were collected using a swab and the composition of the microbiota was determined through 16S rRNA sequencing targeting the V3-V4 hypervariable regions. Initial analyses allowed us to define the core bacterial vaginal microbiota, comprising all the taxa observed in more than 90% of the animals. Consequently, four phyla, 16 families, 14 genera and a single amplicon sequence variant (ASV) met the criteria, suggesting a high diversity of bacterial vaginal microbiota within the studied population. This variability was partially attributed to various environmental factors such as the herd, sampling season, parity, and lactation stage. Next, we identified numerous significant associations between the diversity and composition of the vaginal microbiota and several traits related to host's production and reproduction performance, as well as reproductive tract health. Specifically, 169 genera were associated with at least one trait, with 69% of them significantly associated with multiple traits. Among these, the abundances of Negativibacillus and Ruminobacter were positively correlated with the cows' performances (i.e., longevity, production performances). Other genera showed mixed relationships with the phenotypes, such as Leptotrichia being overabundant in cows with improved fertility records and reproductive tract health, but also in cows with lower production levels. Overall, the numerous associations underscored the complex interactions between the vaginal microbiota and its host. Given the large number of samples collected from commercial farms and the diversity of the phenotypes considered, this study marks an initial step towards a better understanding of the intimate relationship between the vaginal microbiota and the dairy cow's phenotypes.

Short communication: Bifidobacterium abundance in the faecal microbiota is strongly associated with milk traits in dairy cattle.

The performance of dairy cows is influenced by the microbial communities hosted within their digestive tract. While the rumen microbiota has long been associated with host phenotypes, the impact of the faecal microbiota remains elusive. In this study, we collected 697 faecal samples from commercial Holstein cows and analysed them with 16S rRNA gene analyses. For each animal, routinely recorded data, i.e., milk yield, fat yield, protein yield, fat content, protein content, and an aggregate production trait (pINEL) based on the French economic dairy index, were available to assess the links between the faecal microbiota and host production. Our findings revealed a strong and significant association between the structure of the bacterial and prokaryote community (ß-diversity) and dairy production. In addition, differential abundance analyses identified 48 genera whose abundances were significantly associated with pINEL, milk, fat and protein yield. Among these genera, the increased abundance of Bifidobacterium, and particularly an amplicon sequence variant with a 16S rRNA V3-V4 gene region identical to B. globosum and B. pseudolongum, was found to be the most important for high-yielding animals. Bifidobacterium seemed to be a potential key member of the bovine faecal microbiota that should be further investigated. Conversely, the p-1088-a5 gut group genus was found more abundant in low-productive cows. In conclusion, this study demonstrates significant associations between the faecal microbiota and the performance of dairy cows at the whole lactation scale. A better understanding of the physiology of the gut microbiota could help to improve dairy cow production.

Application of Bayesian least absolute shrinkage and selection operator (LASSO) and BayesCπ methods for genomic selection in French Holstein and Montbéliarde breeds.

Recently, the amount of available single nucleotide polymorphism (SNP) marker data has considerably increased in dairy cattle breeds, both for research purposes and for application in commercial breeding and selection programs. Bayesian methods are currently used in the genomic evaluation of dairy cattle to handle very large sets of explanatory variables with a limited number of observations. In this study, we applied 2 bayesian methods, BayesCπ and bayesian least absolute shrinkage and selection operator (LASSO), to 2 genotyped and phenotyped reference populations consisting of 3,940 Holstein bulls and 1,172 Montbéliarde bulls with approximately 40,000 polymorphic SNP. We compared the accuracy of the bayesian methods for the prediction of 3 traits (milk yield, fat content, and conception rate) with pedigree-based BLUP, genomic BLUP, partial least squares (PLS) regression, and sparse PLS regression, a variable selection PLS variant. The results showed that the correlations between observed and predicted phenotypes were similar in BayesCπ (including or not pedigree information) and bayesian LASSO for most of the traits and whatever the breed. In the Holstein breed, bayesian methods led to higher correlations than other approaches for fat content and were similar to genomic BLUP for milk yield and to genomic BLUP and PLS regression for the conception rate. In the Montbéliarde breed, no method dominated the others, except BayesCπ for fat content. The better performances of the bayesian methods for fat content in Holstein and Montbéliarde breeds are probably due to the effect of the DGAT1 gene. The SNP identified by the BayesCπ, bayesian LASSO, and sparse PLS regression methods, based on their effect on the different traits of interest, were located at almost the same position on the genome. As the bayesian methods resulted in regressions of direct genomic values on daughter trait deviations closer to 1 than for the other methods tested in this study, bayesian methods are suggested for genomic evaluations of French dairy cattle.

A comparison of partial least squares (PLS) and sparse PLS regressions in genomic selection in French dairy cattle.

Genomic selection involves computing a prediction equation from the estimated effects of a large number of DNA markers based on a limited number of genotyped animals with phenotypes. The number of observations is much smaller than the number of independent variables, and the challenge is to find methods that perform well in this context. Partial least squares regression (PLS) and sparse PLS were used with a reference population of 3,940 genotyped and phenotyped French Holstein bulls and 39,738 polymorphic single nucleotide polymorphism markers. Partial least squares regression reduces the number of variables by projecting independent variables onto latent structures. Sparse PLS combines variable selection and modeling in a one-step procedure. Correlations between observed phenotypes and phenotypes predicted by PLS and sparse PLS were similar, but sparse PLS highlighted some genome regions more clearly. Both PLS and sparse PLS were more accurate than pedigree-based BLUP and generally provided lower correlations between observed and predicted phenotypes than did genomic BLUP. Furthermore, PLS and sparse PLS required similar computing time to genomic BLUP for the study of 6 traits.

Issues in association mapping with high-density SNP data and diverse family structures.

Genetic association studies have the potential to identify causative genetic variants with small effects in complex diseases, but it is not at all clear which study designs best balance power with sample size, especially when taking into account the difficulty of obtaining a sample of the necessary structure. The 14 contributions from the Genetic Analysis Workshop 15 group 3 used data sets with rheumatoid arthritis as primary phenotype from problem 2 (real data) and Problem 3 (simulated data) to investigate design and analysis problems that arise in candidate-gene, candidate-region, and genome-wide association studies. We identified three major themes that were addressed by multiple groups: (1) comparing family-based and case-control study designs with each other and with hybrid designs incorporating both related and unrelated individuals; (2) exploring and comparing techniques of combining information from multiple, correlated single-nucleotide polymorphisms; and (3) comparing analyses that select the model(s) of best fit with the ultimate aim of detecting the joint effects of several unlinked single-nucleotide polymorphisms. These contributions achieved some success in improving upon existing methods. For example, tests using related cases and unrelated controls can achieve higher power than the tests using unrelated cases and unrelated controls. Aside from these successes, the group 3 contributions highlight some interesting areas for future research.