Dissecting loci that underpin the genetic correlations between production, fertility, and urea traits in Australian Holstein cattle.

Unfavorable genetic correlations between milk production, fertility, and urea traits have been reported. However, knowledge of the genomic regions associated with these unfavorable correlations is limited. Here, we used the correlation scan method to identify and investigate the regions driving or antagonizing the genetic correlations between production vs. fertility, urea vs. fertility, and urea vs. production traits. Driving regions produce an estimate of correlation that is in the same direction as the global correlation. Antagonizing regions produce an estimate in the opposite direction of the global estimates. Our dataset comprised 6567, 4700, and 12,658 Holstein cattle with records of production traits (milk yield, fat yield, and protein yield), fertility (calving interval) and urea traits (milk urea nitrogen and blood urea nitrogen predicted using milk-mid-infrared spectroscopy), respectively. Several regions across the genome drive the correlations between production, fertility, and urea traits. Antagonizing regions were confined to certain parts of the genome and the genes within these regions were mostly involved in preventing metabolic dysregulation, liver reprogramming, metabolism remodeling, and lipid homeostasis. The driving regions were enriched for QTL related to puberty, milk, and health-related traits. Antagonizing regions were mostly related to muscle development, metabolic body weight, and milk traits. In conclusion, we have identified genomic regions of potential importance for dairy cattle breeding. Future studies could investigate the antagonizing regions as potential genomic regions to break the unfavorable correlations and improve milk production as well as fertility and urea traits.

Synchronization of breeding and its impact on genetic parameters and evaluation of female fertility traits.

Achieving an acceptable level of fertility in herds is difficult for many dairy producers because identifying cows in estrus has become challenging owing to poor estrus expression, increased herd size, and lack of time and skilled labor for estrus detection. As a result, synchronization of estrus is often used to manage reproduction. The aims of this study were (1) to identify artificial inseminations (AI) that were performed following synchronization and (2) to assess the effect of synchronization on genetic parameters and evaluation of fertility traits. This study used breeding data collected between 1995 and 2021 from over 4,600 Australian dairy herds that had at least 30 matings per year. Because breeding methods were not reported, the recording pattern of breeding dates showing a large proportion of the total AI being recorded on a single date of the year served as an indicator of synchronization. First, the proportion of AI recorded on each day of the year was calculated for each herd-year. Subsequently, synchronization was defined when a herd with, for instance, only 30 matings in a year, had at least 0.20 or more AI on the same day. As the number of breedings in a herd-year increased, the threshold for classifying AI was continuously reduced from 0.20 to as low as 0.03 under the assumption that mating of many cows on a single date becomes increasingly difficult without synchronization. From the current data, we deduced that 0.11 of all AI were possibly performed following synchronization (i.e., timed AI, TAI). The proportion of AI classified as TAI increased over time and with herd size. Although the deviation from equal numbers of mating on 7 d of the week was not used for classifying AI, 0.44 of AI being categorized as TAI were performed on just 2 d of the week. When data classified as TAI were used for estimating genetic parameters and breeding values, the interval between calving and first service (CFS) was found to be the most affected trait. The phenotypic and additive genetic variance and heritability, as well as variability and reliability of estimated breeding values of bulls and cows for CFS were lower for TAI than for AI performed following detected estrus (i.e., estrus-detected AI, EAI). For calving interval, first service nonreturn rate (FNRR), and successful calving rate to first service, genetic correlations between the same trait measured in TAI and EAI were close to 1, in contrast to 0.55 for CFS. The lower genetic variances and heritabilities for FNRR and calving interval in TAI than in EAI suggests that synchronization reduces the genetic variability of fertility. In conclusion, TAI makes CFS an ineffective measure of fertility. One approach to minimize this effect on genetic evaluations is to identify TAI (using the method described for example) and then set the CFS of these cows as missing records when running multitrait genetic evaluations of fertility traits that include CFS. In the long term, the most practical and accurate way to reduce the effect of synchronization on genetic evaluations is to record TAI along with mating data.

BayesR3 enables fast MCMC blocked processing for largescale multi-trait genomic prediction and QTN mapping analysis.

Bayesian methods, such as BayesR, for predicting the genetic value or risk of individuals from their genotypes, such as Single Nucleotide Polymorphisms (SNP), are often implemented using a Markov Chain Monte Carlo (MCMC) process. However, the generation of Markov chains is computationally slow. We introduce a form of blocked Gibbs sampling for estimating SNP effects from Markov chains that greatly reduces computational time by sampling each SNP effect iteratively n-times from conditional block posteriors. Subsequent iteration over all blocks m-times produces chains of length m × n. We use this strategy to solve large-scale genomic prediction and fine-mapping problems using the Bayesian MCMC mixed-effects genetic model, BayesR3. We validate the method using simulated data, followed by analysis of empirical dairy cattle data using high dimension milk mid infra-red spectra data as an example of "omics" data and show its use to increase the precision of mapping variants affecting milk, fat, and protein yields relative to a univariate analysis of milk, fat, and protein.

Using mid-infrared spectroscopy to increase GWAS power to detect QTL associated with blood urea nitrogen.

Blood urea nitrogen (BUN) is an indicator trait for urinary nitrogen excretion. Measuring BUN level requires a blood sample, which limits the number of records that can be obtained. Alternatively, BUN can be predicted using mid-infrared (MIR) spectroscopy of a milk sample and thus records become available on many more cows through routine milk recording processes. The genetic correlation between MIR predicted BUN (MBUN) and BUN is 0.90. Hence, genetically, BUN and MBUN can be considered as the same trait. The objective of our study was to perform genome-wide association studies (GWAS) for BUN and MBUN, compare these two GWAS and detect quantitative trait loci (QTL) for both traits, and compare the detected QTL with previously reported QTL for milk urea nitrogen (MUN). The dataset used for our analyses included 2098 and 18,120 phenotypes for BUN and MBUN, respectively, and imputed whole-genome sequence data. The GWAS for MBUN was carried out using either the full dataset, the 2098 cows with records for BUN, or 2000 randomly selected cows, so that the dataset size is comparable to that for BUN. The GWAS results for BUN and MBUN were very different, in spite of the strong genetic correlation between the two traits. We detected 12 QTL for MBUN, on bovine chromosomes 2, 3, 9, 11, 12, 14 and X, and one QTL for BUN on chromosome 13. The QTL detected on chromosomes 11, 14 and X overlapped with QTL detected for MUN. The GWAS results were highly sensitive to the subset of records used. Hence, caution is warranted when interpreting GWAS based on small datasets, such as for BUN. MBUN may provide an attractive alternative to perform a more powerful GWAS to detect QTL for BUN.

GWAS and genomic prediction of milk urea nitrogen in Australian and New Zealand dairy cattle.

BACKGROUND: Urinary nitrogen leakage is an environmental concern in dairy cattle. Selection for reduced urinary nitrogen leakage may be done using indicator traits such as milk urea nitrogen (MUN). The result of a previous study indicated that the genetic correlation between MUN in Australia (AUS) and MUN in New Zealand (NZL) was only low to moderate (between 0.14 and 0.58). In this context, an alternative is to select sequence variants based on genome-wide association studies (GWAS) with a view to improve genomic prediction accuracies. A GWAS can also be used to detect quantitative trait loci (QTL) associated with MUN. Therefore, our objectives were to perform within-country GWAS and a meta-GWAS for MUN using records from up to 33,873 dairy cows and imputed whole-genome sequence data, to compare QTL detected in the GWAS for MUN in AUS and NZL, and to use sequence variants selected from the meta-GWAS to improve the prediction accuracy for MUN based on a joint AUS-NZL reference set. RESULTS: Using the meta-GWAS, we detected 14 QTL for MUN, located on chromosomes 1, 6, 11, 14, 19, 22, 26 and the X chromosome. The three most significant QTL encompassed the casein genes on chromosome 6, PAEP on chromosome 11 and DGAT1 on chromosome 14. We selected 50,000 sequence variants that had the same direction of effect for MUN in AUS and MUN in NZL and that were most significant in the meta-analysis for the GWAS. The selected sequence variants yielded a genetic correlation between MUN in AUS and MUN in NZL of 0.95 and substantially increased prediction accuracy in both countries. CONCLUSIONS: Our results demonstrate how the sharing of data between two countries can increase the power of a GWAS and increase the accuracy of genomic prediction using a multi-country reference population and sequence variants selected based on a meta-GWAS.

Multiple Country Approach to Improve the Test-Day Prediction of Dairy Cows' Dry Matter Intake.

We predicted dry matter intake of dairy cows using parity, week of lactation, milk yield, milk mid-infrared (MIR) spectrum, and MIR-based predictions of bodyweight, fat, protein, lactose, and fatty acids content in milk. The dataset comprised 10,711 samples of 534 dairy cows with a geographical diversity (Australia, Canada, Denmark, and Ireland). We set up partial least square (PLS) regressions with different constructs and a one-hidden-layer artificial neural network (ANN) using the highest contribution variables. In the ANN, we replaced the spectra with their projections to the 25 first PLS factors explaining 99% of the spectral variability to reduce the model complexity. Cow-independent 10 × 10-fold cross-validation (CV) achieved the best performance with root mean square errors (RMSECV) of 3.27 ± 0.08 kg for the PLS regression and 3.25 ± 0.13 kg for ANN. Although the available data were significantly different, we also performed a country-independent validation (CIV) to measure the models' performance fairly. We found RMSECIV varying from 3.73 to 6.03 kg for PLS and 3.69 to 5.08 kg for ANN. Ultimately, based on the country-independent validation, we discussed the developed models' performance with those achieved by the National Research Council's equation.

Validation of Dairy Cow Bodyweight Prediction Using Traits Easily Recorded by Dairy Herd Improvement Organizations and Its Potential Improvement Using Feature Selection Algorithms.

Knowing the body weight (BW) of a cow at a specific moment or measuring its changes through time is of interest for management purposes. The current work aimed to validate the feasibility of predicting BW using the day in milk, parity, milk yield, and milk mid-infrared (MIR) spectrum from a multiple-country dataset and reduce the number of predictors to limit the risk of over-fitting and potentially improve its accuracy. The BW modeling procedure involved feature selections and herd-independent validation in identifying the most interesting subsets of predictors and then external validation of the models. From 1849 records collected in 9 herds from 360 Holstein cows, the best performing models achieved a root mean square error (RMSE) for the herd-independent validation between 52 ± 2.34 kg to 56 ± 3.16 kg, including from 5 to 62 predictors. Among these models, three performed remarkably well in external validation using an independent dataset (N = 4067), resulting in RMSE ranging from 52 to 56 kg. The results suggest that multiple optimal BW predictive models coexist due to the high correlations between adjacent spectral points.

Mitochondrial protein gene expression and the oxidative phosphorylation pathway associated with feed efficiency and energy balance in dairy cattle.

Feed efficiency and energy balance are important traits underpinning profitability and environmental sustainability in animal production. They are complex traits, and our understanding of their underlying biology is currently limited. One measure of feed efficiency is residual feed intake (RFI), which is the difference between actual and predicted intake. Variation in RFI among individuals is attributable to the metabolic efficiency of energy utilization. High RFI (H_RFI) animals require more energy per unit of weight gain or milk produced compared with low RFI (L_RFI) animals. Energy balance (EB) is a closely related trait calculated very similarly to RFI. Cellular energy metabolism in mitochondria involves mitochondrial protein (MiP) encoded by both nuclear (NuMiP) and mitochondrial (MtMiP) genomes. We hypothesized that MiP genes are differentially expressed (DE) between H_RFI and L_RFI animal groups and similarly between negative and positive EB groups. Our study aimed to characterize MiP gene expression in white blood cells of H_RFI and L_RFI cows using RNA sequencing to identify genes and biological pathways associated with feed efficiency in dairy cattle. We used the top and bottom 14 cows ranked for RFI and EB out of 109 animals as H_RFI and L_RFI, and positive and negative EB groups, respectively. The gene expression counts across all nuclear and mitochondrial genes for animals in each group were used for differential gene expression analyses, weighted gene correlation network analysis, functional enrichment, and identification of hub genes. Out of 244 DE genes between RFI groups, 38 were MiP genes. The DE genes were enriched for the oxidative phosphorylation (OXPHOS) and ribosome pathways. The DE MiP genes were underexpressed in L_RFI (and negative EB) compared with the H_RFI (and positive EB) groups, suggestive of reduced mitochondrial activity in the L_RFI group. None of the MtMiP genes were among the DE MiP genes between the groups, which suggests a non-rate limiting role of MtMiP genes in feed efficiency and warrants further investigation. The role of MiP, particularly the NuMiP and OXPHOS pathways in RFI, was also supported by our gene correlation network analysis and the hub gene identification. We validated the findings in an independent data set. Overall, our study suggested that differences in feed efficiency in dairy cows may be linked to differences in cellular energy demand. This study broadens our knowledge of the biology of feed efficiency in dairy cattle.