Effect of Genetic and Environmental Factors on Twin Pregnancy in Primiparous Dairy Cows.

Twin pregnancies are highly undesirable in dairy cattle; they compromise the health and wellbeing of a cow and dramatically impair the farm economy. Recently, a genomic prediction for twin pregnancies has been developed. The objective of this study was to assess environmental and management risk factors affecting the incidence of twin pregnancies in high-producing dairy cows in their first lactation, with a special emphasis placed on the genomic prediction values of twin pregnancy. Our study population of primiparous cows proved valuable in identifying factors other than genomic predictive values that influence twin pregnancy rates. The odds ratio for twin pregnancies was 0.85 (p < 0.0001) for each unit of a prediction value increase, 3.5 (p = 0.023) for cows becoming pregnant during the negative photoperiod, and 0.33 (p = 0.016) for cows producing ≥42 kg of milk at AI, compared with the remaining cows who produced <42 kg of milk. As a general conclusion, the practical implication of our findings is that genomic prediction values can identify the risk of twin pregnancy at a herd level. Given the cumulative effect of genomic selection, selecting animals with a reduced genetic risk of twin pregnancies can contribute to reducing the incidence of twin pregnancies in dairy herds.

Genomic Prediction for Abortion in Lactating Holstein Dairy Cows.

Abortion in dairy cattle causes great economic losses due to reduced animal health, increase in culling rates, reduction in calf production, and milk yield, among others. Although the etiology of abortions can be of various origins, previous research has shown a genetic component. The objectives of this study were to (1) describe the development of the genomic prediction for cow abortions in lactating Holstein dairy cattle based on producer-recorded data and ssGBLUP methodology and (2) evaluate the efficacy of genomic predictions for cow abortions in commercial herds of US Holstein cows using data from herds that do not contribute phenotypic information to the evaluation. We hypothesized that cows with greater genomic predictions for cow abortions (Z_Abort STA) would have a reduced incidence of abortion. Phenotypic data on abortions, pedigree, and genotypes were collected directly from commercial dairy producers upon obtaining their permission. Abortion was defined as the loss of a confirmed pregnancy after 42 and prior to 260 days of gestation, treated as a binary outcome (0, 1), and analyzed using a threshold model. Data from a different subset of animals were used to test the efficacy of the prediction. The additive genetic variance for the cow abortion trait (Z_Abort) was 0.1235 and heritability was 0.0773. For all animals with genotypes (n = 1,662,251), mean reliability was 42%, and genomic predicted transmitting abilities (gPTAs) ranged from −8.8 to 12.4. Z_Abort had a positive correlation with cow and calf health traits and reproductive traits, and a negative correlation with production traits. Z_Abort effectively identified cows with a greater or lesser risk of abortion (16.6% vs. 11.0% for the worst and best genomics groups, respectively; p < 0.0001). The inclusion of cow abortion genomic predictions in a multi-trait selection index would allow dairy producers and consultants to reduce the incidence of abortion and to select high-producing, healthier, and more profitable cows.

Genomic Prediction for Twin Pregnancies.

Twinning is a multifactorial trait influenced by both genetic and environmental factors that can negatively impact animal welfare and economic sustainability on commercial dairy operations. To date, using genetic selection as a tool for reducing twinning rates on commercial dairies has been proposed, but not yet implemented. In response to this market need, Zoetis (Kalamazoo, MI, USA) has developed a genomic prediction for twin pregnancies, and included it in a comprehensive multitrait selection index. The objectives of this study were to (1) describe a genetic evaluation for twinning in Holstein cattle, (2) demonstrate the efficacy of the predictions, (3) propose strategies to reduce twin pregnancies using this information. Data were retrieved from commercial dairies and provided directly by producers upon obtaining their permission. The twin pregnancies trait (TWIN) was defined as a pregnancy resulting in birth or abortion of twin calves, classified as a binary (0,1) event, and analysed using a threshold animal model. Predictions for a subset of cows were compared to their on-farm twin records. The heritability for twin pregnancies was 0.088, and genomic predicted transmitting abilities ((g)PTAs) ranged from -7.45-20.79. Genetic correlations between TWIN and other traits were low, meaning that improvement for TWIN will not negatively impact improvement for other traits. TWIN was effectively demonstrated to identify cows most and least likely to experience a twin pregnancy in a given lactation, regardless of reproductive protocol used. Effective inclusion of the prediction in a multitrait selection index offers producers a comprehensive tool to inform selection and management decisions. When combined with sound management practices, this presents a compelling opportunity for dairy producers to proactively reduce the incidence of twin pregnancies on commercial dairy operations.

Validation of genomic predictions for a lifetime merit selection index for the US dairy industry.

Selection indices are a critical component of many breeding programs. A common purpose of a selection index is to predict an animal's genetic potential for total economic merit. The objective of this study was to evaluate retrospectively whether a specific selection index comprising genomically-enhanced predicted transmitting abilities had the ability to predict observed lifetime profit in US Holstein animals. The selection index evaluated was dairy wellness profit (DWP$). In total, 2,185 animals were included in this study. Index values were used to rank and assign animals to quartiles (genetic groups: worst 25%, 26-50%, 51-75%, and best 25%). Generalized linear mixed effects models were applied to estimate the associations between index quartile and defined economic outcomes. Similar analyses were conducted to estimate associations between index quartile and observed phenotype to characterize the extent to which profitability outcomes were driven by economically relevant production and health traits. Differences in lifetime profit and annuity value between the best and worst genetic groups for DWP$ were $811 (±297) and $232 (±88), respectively. Significant differences were also observed between top and bottom quartiles for milk production (8,077 kg), fat production (336 kg), protein production (264 kg), live calves (0.5), time spent in the lactating herd (6.6 mo), and cow mortality (8.4%). Additionally, differences in disease incidence were significant between the best and worst DWP$ quartiles for metritis (5.2%), mastitis (14.9%), and lameness (15.9%). The observed results of this study demonstrated the ability of DWP$ predictions to predict lifetime profitability of Holstein animals and its potential utility as a tool to guide selection and breeding programs. Improving DWP$ through genetic selection, when combined with good management practices, provides an opportunity for dairy producers to improve overall herd profitability.

Autozygosity islands and ROH patterns in Nellore lineages: evidence of selection for functionally important traits.

BACKGROUND: The aim of this study was to assess genome-wide autozygosity in a Nellore cattle population and to characterize ROH patterns and autozygosity islands that may have occurred due to selection within its lineages. It attempts also to compare estimates of inbreeding calculated from ROH (FROH), genomic relationship matrix (FGRM), and pedigree-based coefficient (FPED). RESULTS: The average number of ROH per animal was 55.15 ± 13.01 with an average size of 3.24 Mb. The Nellore genome is composed mostly by a high number of shorter segments accounting for 78% of all ROH, although the proportion of the genome covered by them was relatively small. The genome autozygosity proportion indicates moderate to high inbreeding levels for classical standards, with an average value of 7.15% (178.70 Mb). The average of FPED and FROH, and their correlations (- 0.05 to 0.26) were low. Estimates of correlation between FGRM-FPED was zero, while the correlation (- 0.01 to - 0.07) between FGRM-FROH decreased as a function of ROH length, except for FROH > 8Mb (- 0.03). Overall, inbreeding coefficients were not high for the genotyped animals. Autozygosity islands were evident across the genome (n = 62) and their genomic location did not largely differ within lineages. Enriched terms (p < 0.01) associated with defense response to bacteria (GO:0042742), immune complex reaction (GO:0045647), pregnancy-associated glycoproteins genes (GO:0030163), and organism growth (GO:0040014) were described within the autozygotic islands. CONCLUSIONS: Low FPED-FROH correlation estimates indicate that FPED is not the most suitable method for capturing ancient inbreeding when the pedigree does not extend back many generations and FROH should be used instead. Enriched terms (p < 0.01) suggest a strong selection for immune response. Non-overlapping islands within the lineages greatly explain the mechanism underlying selection for functionally important traits in Nellore cattle.

Genome-wide association study provides insights into genes related with horn development in Nelore beef cattle.

The causal mutation for polledness in Nelore (Bos taurus indicus) breed seems to have appeared first in Brazil in 1957. The expression of the polled trait is known to be ruled by a few groups of alleles in taurine breeds; however, the genetic basis of this trait in indicine cattle is still unclear. The aim of this study was to identify genomic regions associated with the hornless trait in a commercial Nelore population. A total of 107,294 animals had phenotypes recorded and 2,238 were genotyped/imputed for 777k SNP. The weighted single-step approach for genome-wide association study (WssGWAS) was used to estimate the SNP effects and variances accounted for by 1 Mb sliding SNP windows. A centromeric region of chromosome 1 with 3.11 Mb size (BTA1: 878,631-3,987,104 bp) was found to be associated with hornless in the studied population. A total of 28 protein-coding genes are mapped in this region, including the taurine Polled locus and the IFNAR1, IFNAR2, IFNGR2, KRTAP11-1, MIS18A, OLIG1, OLIG2, and SOD1 genes, which expression can be related to the horn formation as described in literature. The functional enrichment analysis by DAVID tool revealed cytokine-cytokine receptor interaction, JAK-STAT signaling, natural killer cell mediated cytotoxicity, and osteoclast differentiation pathways as significant (P < 0.05). In addition, a runs of homozygosity (ROH) analysis identified a ROH island in polled animals with 2.47 Mb inside the region identified by WssGWAS. Polledness in Nelore cattle is associated with one region in the genome with 3.1 Mb size in chromosome 1. Several genes are harbored in this region, and they may act together in the determination of the polled/horned phenotype. Fine mapping the locus responsible for polled trait in Nelore breed and the identification of the molecular mechanisms regulating the horn growth deserve further investigation.

Validation of genomic predictions for wellness traits in US Holstein cows.

The objective of this study was to evaluate the efficacy of wellness trait genetic predictions in commercial herds of US Holstein cows from herds that do not contribute phenotypic information to the evaluation. Tissue samples for DNA extraction were collected from more than 3,400 randomly selected pregnant Holstein females in 11 herds and 2 age groups (69% nulliparous, 31% primiparous) approximately 30 to 60 d before their expected calving date. Lactation records from cows that calved between September 1, 2015, and December 31, 2015, were included in the analysis. Genomically enhanced predicted transmitting abilities for the wellness traits of retained placenta, metritis, ketosis, displaced abomasum, mastitis, and lameness were estimated by the Zoetis genetic evaluation and converted into standardized transmitting abilities. Mean reliabilities of the animals in the study ranged between 45 and 47% for each of the 6 traits. Animals were ranked by their standardized transmitting abilities within herd and age group then assigned to 1 of 4 groups of percentile-based genetic groups of equal size. Adverse health events, including retained placenta, metritis, ketosis, displaced abomasum, mastitis, and lameness, were collected from on-farm herd management software, and animal phenotype was coded as either healthy (0), diseased (1), or excluded for each of the 6 outcomes of interest. Statistical analysis was performed using a generalized linear mixed model with genetic group, age group, and lactation as fixed effects, whereas herd and animal nested within herd were set as random effects. Results of the analysis indicated that the wellness trait predictions were associated with differences in phenotypic disease incidence between the worst and best genetic groups. The difference between the worst and best genetic groups in recorded disease incidence was 2.9% for retained placenta, 10.8% for metritis, 1.1% for displaced abomasum, 1.7% for ketosis, 7.4% for mastitis, and 3.9% for lameness. Odds ratio estimates between the highest and lowest genetic groups ranged from 1.6 (lameness) to 17.1 (displaced abomasum) for the 6 traits analyzed. These results indicate that wellness trait information of young calves and heifers can be used to effectively predict meaningful differences in future health performance. Improving wellness traits through direct genetic selection presents a compelling opportunity for dairy producers to help reduce disease incidence and improve profitability when coupled with sound management practices.