Impacts of preweaning colostrum feeding practices and health measures on dairy cow production, while accounting for genetic potential

Calf management and health are essential for setting up the foundation of a productive cow. The objectives of this study were to estimate the impact of preweaning practices on milk production parameters while accounting for an animal's genetic potential in New Brunswick, Canada. A retrospective cohort study was performed on 220 heifer calves from eight herds born in 2014-2015. Preweaning practices and health data were recorded by producers and reviewed by the herd veterinarian for each calf. The herd veterinarian also visited the farms to collect serum samples from calves and frozen colostrum samples. The production outcomes assessed were milk, protein and fat yields, standardized to 305 d for the first lactation (L1) and a combined group of lactations two and three (L2 + 3). The genomic potential was determined as genomic parent averages (GPA) for the associated production parameters. Analysis was performed with multivariable linear (L1) and linear mixed (L2 + 3) regression models. In L1, for every 1.0 kg increase in weaning weight, milk, protein, and fat yield increased by 25.5, 0.82, and 1.01 kg, respectively (P < 0.006). Colostrum feeding time (CFT) positively impacted L1 milk and protein production, with feeding between 1-2 h of life producing the greatest estimates of 626 kg of milk and 18.2 kg of protein yield (P < 0.007), compared to earlier or later CFT. Fat yield production was decreased by 80.5 kg (P < 0.006) in L1 when evaluating animals that developed a preweaning disease and were not treated with antibiotics compared to healthy untreated animals. Impacts on L2 + 3 were similar across all production outcomes, with a positive interaction effect of CFT and weaning weight. Compared to CFT < 1 h, the later CFT groups of 1-2 h and > 2 h produced greater yield outcomes of 68.2 to 72.6 kg for milk (P < 0.006), 2.06 to 2.15 kg for protein (P < 0.005), and 1.8 to 1.9 kg for fat (P < 0.045) for every 1 kg increase of weaning weight, respectively. The fit of all models was significantly improved with the inclusion of GPA. These results indicate that colostrum management and preweaning health measures impacted production parameters as adults. The inclusion of GPA significantly improved the accuracy of the models, indicating that this can be an important parameter to include in future studies.

Genomic Characterization and Initial Insight into Mastitis-Associated SNP Profiles of Local Latvian Bos taurus Breeds

Latvia has two local Bos taurus breeds-Latvian Brown (LBG) and Latvian Blue (LZG)-characterized by a good adaptation to the local climate, longevity, and high fat and protein contents in milk. Since these are desired traits in the dairy industry, this study investigated the genetic background of the LBG and LZG breeds and identified the genetic factors associated with mastitis. Blood and semen samples were acquired, and whole genome sequencing was then performed to acquire a genomic sequence with at least 35× or 10× coverage. The heterozygosity, nucleotide diversity, and LD analysis indicated that LBG and LZG cows have similar levels of genetic diversity compared to those of other breeds. An analysis of the population structure revealed that each breed clustered together, but the overall differentiation between the breeds was small. The highest genetic variance was observed in the LZG breed compared with the LBG breed. Our results show that SNP rs721295390 is associated with mastitis in the LBG breed, and SNPs rs383806754, chr29:43998719CG>C, and rs462030680 are associated with mastitis in the LZG breed. This study shows that local Latvian LBG and LZG breeds have a pronounced genetic differentiation, with each one suggesting its own mastitis-associated SNP profile.

Multibreed genomic prediction using summary statistics and a breed-origin-of-alleles approach

Because of an increasing interest in crossbreeding between dairy breeds in dairy cattle herds, farmers are requesting breeding values for crossbred animals. However, genomically enhanced breeding values are difficult to predict in crossbred populations because the genetic make-up of crossbred individuals is unlikely to follow the same pattern as for purebreds. Furthermore, sharing genotype and phenotype information between breed populations are not always possible, which means that genetic merit (GM) for crossbred animals may be predicted without the information needed from some pure breeds, resulting in low prediction accuracy. This simulation study investigated the consequences of using summary statistics from single-breed genomic predictions for some or all pure breeds in two- and three-breed rotational crosses, rather than their raw data. A genomic prediction model taking into account the breed-origin of alleles (BOA) was considered. Because of a high genomic correlation between the breeds simulated (0.62-0.87), the prediction accuracies using the BOA approach were similar to a joint model, assuming homogeneous SNP effects for these breeds. Having a reference population with summary statistics available from all pure breeds and full phenotype and genotype information from crossbreds yielded almost as high prediction accuracies (0.720-0.768) as having a reference population with full information from all pure breeds and crossbreds (0.753-0.789). Lacking information from the pure breeds yielded much lower prediction accuracies (0.590-0.676). Furthermore, including crossbred animals in a combined reference population also benefitted prediction accuracies in the purebred animals, especially for the smallest breed population.

Genomics in animal breeding from the perspectives of matrices and molecules

BACKGROUND

This paper describes genomics from two perspectives that are in use in animal breeding and genetics: a statistical perspective concentrating on models for estimating breeding values, and a sequence perspective concentrating on the function of DNA molecules.

MAIN BODY

This paper reviews the development of genomics in animal breeding and speculates on its future from these two perspectives. From the statistical perspective, genomic data are large sets of markers of ancestry; animal breeding makes use of them while remaining agnostic about their function. From the sequence perspective, genomic data are a source of causative variants; what animal breeding needs is to identify and make use of them.

CONCLUSION

The statistical perspective, in the form of genomic selection, is the more applicable in contemporary breeding. Animal genomics researchers using from the sequence perspective are still working towards this the isolation of causative variants, equipped with new technologies but continuing a decades-long line of research.

Blood Transcriptome Analysis of Beef Cow with Different Parity Revealed Candidate Genes and Gene Networks Regulating the Postpartum Diseases

Maternal parity is an important physiological factor influencing beef cow reproductive performance. However, there are few studies on the influence of different calving periods on early growth and postpartum diseases. Here, we conducted blood transcriptomic analysis on cows of different parities for gene discovery. We used Short Time Series Expression Miner (STEM) analysis to determine gene expression levels in cows of various parities and divided multiple parities into three main periods (nulliparous, primiparous, and multiparous) for subsequent analysis. Furthermore, the top 15,000 genes with the lowest median absolute deviation (MAD) were used to build a co-expression network using weighted correlation network analysis (WGCNA), and six independent modules were identified. Combing with Exon Wide Selection Signature (EWSS) and protein-protein interaction (PPI) analysis revealed that TPCN2, KIF22, MICAL3, RUNX2, PDE4A, TESK2, GPM6A, POLR1A, and KLHL6 involved in early growth and postpartum diseases. The GO and KEGG enrichment showed that the Parathyroid hormone synthesis, secretion, and action pathway and stem cell differentiation function-related pathways were enriched. Collectively, our study revealed candidate genes and gene networks regulating the early growth and postpartum diseases and provided new insights into the potential mechanism of reproduction advantages of different parity selection.

Breed of origin of alleles and genomic predictions for crossbred dairy cows

BACKGROUND

In dairy cattle, genomic selection has been implemented successfully for purebred populations, but, to date, genomic estimated breeding values (GEBV) for crossbred cows are rarely available, although they are valuable for rotational crossbreeding schemes that are promoted as efficient strategies. An attractive approach to provide GEBV for crossbreds is to use estimated marker effects from the genetic evaluation of purebreds. The effects of each marker allele in crossbreds can depend on the breed of origin of the allele (BOA), thus applying marker effects based on BOA could result in more accurate GEBV than applying only proportional contribution of the purebreds. Application of BOA models in rotational crossbreeding requires methods for detecting BOA, but the existing methods have not been developed for rotational crossbreeding. Therefore, the aims of this study were to develop and test methods for detecting BOA in a rotational crossbreeding system, and to investigate methods for calculating GEBV for crossbred cows using estimated marker effects from purebreds.

RESULTS

For detecting BOA in crossbred cows from rotational crossbreeding for which pedigree is recorded, we developed the AllOr method based on the comparison of haplotypes in overlapping windows. To calculate the GEBV of crossbred cows, two models were compared: a BOA model where marker effects estimated from purebreds are combined based on the detected BOA; and a breed proportion model where marker effects are combined based on estimated breed proportions. The methods were tested on simulated data that mimic the first four generations of rotational crossbreeding between Holstein, Jersey and Red Dairy Cattle. The AllOr method detected BOA correctly for 99.6% of the marker alleles across the four crossbred generations. The reliability of GEBV was higher with the BOA model than with the breed proportion model for the four generations of crossbreeding, with the largest difference observed in the first generation.

CONCLUSIONS

In rotational crossbreeding for which pedigree is recorded, BOA can be accurately detected using the AllOr method. Combining marker effects estimated from purebreds to predict the breeding value of crossbreds based on BOA is a promising approach to provide GEBV for crossbred dairy cows.

Investigation of β-hydroxybutyrate in early lactation of Simmental cows: Genetic parameters and genomic predictions

Genomic information allows for a more accurate calculation of relationships among animals than the pedigree information, leading to an increase in accuracy of breeding values. Here, we used pedigree-based and single-step genomic approaches to estimate variance components and breeding values for β-hydroxybutyrate milk content (BHB). Additionally, we performed a genome-wide association study (GWAS) to depict its genetic architecture. BHB concentrations within the first 90 days of lactation, estimated from milk medium infrared spectra, were available for 30,461 cows (70,984 records). Genotypes at 42,152 loci were available for 9,123 animals. Low heritabilities were found for BHB using pedigree-based (0.09 ± 0.01) and genomic (0.10 ± 0.01) approaches. Genetic correlation between BHB and milk traits ranged from -0.27 ± 0.06 (BHB and protein percentage) to 0.13 ± 0.07 (BHB and fat-to-protein ratio) using pedigree and from -0.26 ± 0.05 (BHB and protein percentage) to 0.13 ± 0.06 (BHB and fat-to-protein ratio) using genomics. Breeding values were validated for 344 genotyped cows using linear regression method. The genomic EBV (GEBV) had greater accuracy (0.51 vs. 0.45) and regression coefficient (0.98 vs. 0.95) compared to EBV. The correlation between two subsequent evaluations, without and with phenotypes for validation cows, was 0.85 for GEBV and 0.82 for EBV. Predictive ability (correlation between (G)EBV and adjusted phenotypes) was greater when genomic information was used (0.38) than in the pedigree-based approach (0.31). Validation statistics in the pairwise two-trait models (milk yield, fat and protein percentage, urea, fat/protein ratio, lactose and logarithmic transformation of somatic cells count) were very similar to the ones highlighted for the single-trait model. The GWAS allowed discovering four significant markers located on BTA20 (57.5-58.2 Mb), where the ANKH gene is mapped. This gene has been associated with lactose, alpha-lactalbumin and BHB. Results of this study confirmed the usefulness of genomic information to provide more accurate variance components and breeding values, and important insights about the genomic determination of BHB milk content.

Genetic consequences of terminal crossbreeding, genomic test, sexed semen, and beef semen in dairy herds

The development of breeding tools, such as genomic selection and sexed semen, has progressed rapidly in dairy cattle breeding during the past decades. In combination with beef semen, these tools are adopted increasingly at herd level. Dairy crossbreeding is emerging, but the economic and genetic consequences of combining it with the other breeding tools are relatively unknown. We investigated 5 different sexed semen schemes where 0, 50, and 90% of the heifers; 50% of the heifers + 25% of the first-parity cows; and 90% of the heifers + 45% of the first-parity cows were bred to sexed semen. The 5 schemes were combined in scenarios managing pure-breeding or terminal crossbreeding, including genomic testing of all newborn heifers or no testing, and keeping Swedish Red or Swedish Holstein as an initial breed. Thus, 40 scenarios were simulated, combining 2 stochastic simulation models: SimHerd Crossbred (operational returns) and ADAM (genetic returns). The sum of operational and genetic returns equaled the total economic return. Beef semen was used in all scenarios to limit the surplus of replacement heifers. Terminal crossbreeding implied having a nucleus of purebred females, where some were inseminated with semen of the opposite breed. The F₁ crossbred females were inseminated with beef semen. The reproductive performance played a role in improving the benefit of any of the tools. The most considerable total economic returns were achieved when all 4 breeding tools were combined. For Swedish Holstein, the highest total economic return compared with a pure-breeding scenario, without sexed semen and genomic test, was achieved when 90% sexed semen was used in heifers and 45% sexed semen was used for first-parity cows combined with genomic test and crossbreeding (+€58, 33% crossbreds in the herd). The highest total economic return for Swedish Red compared with a pure-breeding scenario, without sexed semen and genomic test, was achieved when 90% sexed semen was used in heifers combined with genomic test and crossbreeding (+€94, 46% crossbreds in the herd). Terminal crossbreeding resulted in lower genetic returns across the herd compared with the corresponding pure-breeding scenarios but was compensated by a higher operational return.

Genomic information allows for more accurate breeding values for milkability in dual-purpose Italian Simmental cattle

Milkability is a trait related to the milking efficiency of an animal, and it is a component of the herd profitability. Due to its economic importance, milkability is currently included in the selection index of the Italian Simmental cattle breed with a weight of 7.5%. This lowly heritable trait is measured on a subjective scale from 1 to 3 (1 = slow, 3 = fast), and genetic evaluations are performed by pedigree-based BLUP. Genomic information is now available for some animals in the Italian Simmental population, and its inclusion in the genetic evaluation system could increase accuracy of breeding values and genetic progress for milkability. The aim of this study was to test the feasibility and advantages of having a genomic evaluation for this trait in the Italian Simmental population. Phenotypes were available for 131,308 cows. A total of 9,526 animals had genotypes for 42,152 loci; among the genotyped animals, 2,455 were cows with phenotypes, and the other were their relatives. The youngest cows with both phenotypes and genotypes (n = 900) were identified as selection candidates. Variance components and heritability were estimated using pedigree information, whereas genetic and genomic evaluations were carried out using BLUP and single-step genomic BLUP (ssGBLUP), respectively. In addition, a weighted ssGBLUP was assessed using genomic regions from a genome-wide association study. Evaluation models were validated using theoretical and realized accuracies. The estimated heritability for milkability was 0.12 ± 0.01. The mean theoretical accuracies for selection candidates were 0.43 ± 0.08 (BLUP) and 0.53 ± 0.06 (ssGBLUP). The mean realized accuracies based on linear regression statistics were 0.29 (BLUP) and 0.40 (ssGBLUP). No genomic regions were significantly associated with milkability, thus no improvements in accuracy were observed when using weighted ssGBLUP. Results indicated that genomic information could improve the accuracy of breeding values and increase genetic progress for milkability in Italian Simmental.