Comparison of parametric, orthogonal, and spline functions to model individual lactation curves for milk yield in Canadian Holsteins

Test day records for milk yield of 57,390 first lactation Canadian Holsteins were analyzed with a linear model that included the fixed effects of herd-test date and days in milk (DIM) interval nested within age and calving season. Residuals from this model were analyzed as a new variable and fitted with a five parameter model, fourth-order Legendre polynomials, with linear, quadratic and cubic spline models with three knots. The fit of the models was rather poor, with about 30%-40% of the curves showing an adjusted R-square lower than 0.20 across all models. Results underline a great difficulty in modelling individual deviations around the mean curve for milk yield. However, the Ali and Schaeffer (5 parameter) model and the fourth-order Legendre polynomials were able to detect two basic shapes of individual deviations among the mean curve. Quadratic and, especially, cubic spline functions had better fitting performances but a poor predictive ability due to their great flexibility that results in an abrupt change of the estimated curve when data are missing. Parametric and orthogonal polynomials seem to be robust and affordable under this standpoint

Non-additive genetic effects for fertility traits in Canadian Holstein cattle (Open Access publication )

The effects of additive, dominance, additive by dominance, additive by additive and dominance by dominance genetic effects on age at first service, non-return rates and interval from calving to first service were estimated. Practical considerations of computing additive and dominance relationships using the genomic relationship matrix are discussed. The final strategy utilized several groups of 1000 animals (heifers or cows) in which all animals had a non-zero dominance relationship with at least one other animal in the group. Direct inversion of relationship matrices was possible within the 1000 animal subsets. Estimates of variances were obtained using Bayesian methodology via Gibbs sampling. Estimated non-additive genetic variances were generally as large as or larger than the additive genetic variance in most cases, except for non-return rates and interval from calving to first service for cows. Non-additive genetic effects appear to be of sizeable magnitude for fertility traits and should be included in models intended for estimating additive genetic merit. However, computing additive and dominance relationships for all possible pairs of individuals is very time consuming in populations of more than 200 000 animals

Fit of different functions to the individual deviations in random regression test day models for milk yield in dairy cattle

The shape of individual deviations of milk yield for dairy cattle from the fixed part of a random regression test day model (RRTDM) was investigated. Data were 53,217 TD records for milk yield of 6,229 first lactation Canadian Holsteins in Ontario. Data were fitted with a model that included the fixed effects of herd-testdate, DIM interval nested within age and season of calving. Residuals of the model were then fitted with the following functions: Ali and Schaeffer 5 parameter model, fourth-order Legendre Polynomials, and cubic spline with three, four or five knots. Result confirm the great variability of shape that can be found when individual lactation are modeled. Cubic splines gave better fitting pe4rformances although together with a marked tendency to yield aberrant estimates at the edge of the lactation trajectory

Genome-wide association study to identify genomic regions and positional candidate genes associated with male fertility in beef cattle

[EN] Fertility plays a key role in the success of calf production, but there is evidence that reproductive efficiency in beef cattle has decreased during the past half-century worldwide. Therefore, identifying animals with superior fertility could significantly impact cow-calf production efficiency. The objective of this research was to identify candidate regions affecting bull fertility in beef cattle and positional candidate genes annotated within these regions. A GWAS using a weighted single-step genomic BLUP approach was performed on 265 crossbred beef bulls to identify markers associated with scrotal circumference (SC) and sperm motility (SM). Eight windows containing 32 positional candidate genes and five windows containing 28 positional candidate genes explained more than 1% of the genetic variance for SC and SM, respectively. These windows were selected to perform gene annotation, QTL enrichment, and functional analyses. Functional candidate gene prioritization analysis revealed 14 prioritized candidate genes for SC of which MAP3K1 and VIP were previously found to play roles in male fertility. A different set of 14 prioritized genes were identified for SM and five were previously identified as regulators of male fertility (SOD2, TCP1, PACRG, SPEF2, PRLR). Significant enrichment results were identified for fertility and body conformation QTLs within the candidate windows. Gene ontology enrichment analysis including biological processes, molecular functions, and cellular components revealed significant GO terms associated with male fertility. The identification of these regions contributes to a better understanding of fertility associated traits and facilitates the discovery of positional candidate genes for future investigation of causal mutations and their implications.SIThe authors acknowledge financial support from the (FDE.13.17) Sustainable Beef and Forage Science Cluster funded by the Canadian Beef Cattle Check-Off, Beef Cattle Research Council (BCRC), Alberta Beef Producers, Alberta Cattle Feeders’ Association, Beef Farmers of Ontario, La Fédération des Productuers de bovins du Québec, and Agriculture and Agri-Food Canada’s Canadian Agricultural Partnership. This study was also supported by the Ontario Ministry of Agriculture, Food, and Rural Affairs (OMAFRA), Ontario Ministry of Research and Innovation, Agriculture and Agri-Food Canada (AAFC), and Natural Sciences and Engineering Research Council of Canada (NSERC) Discovery Grant. Hannah Sweett was supported by the OMAFRA Highly Qualified Personnel Scholarship Program

Identification of unique ROH regions with unfavorable effects on production and fertility traits in Canadian Holsteins

Background: The advent of genomic information and the reduction in the cost of genotyping have led to the use of genomic information to estimate genomic inbreeding as an alternative to pedigree inbreeding. Using genomic measures, effects of genomic inbreeding on production and fertility traits have been observed. However, there have been limited studies on the specific genomic regions causing the observed negative association with the trait of interest. Our aim was to identify unique run of homozygosity (ROH) genotypes present within a given genomic window that display negative associations with production and fertility traits and to quantify the effects of these identified ROH genotypes. Methods: In total, 50,575 genotypes based on a 50K single nucleotide polymorphism (SNP) array and 259,871 pedigree records were available. Of these 50,575 genotypes, 46,430 cows with phenotypic records for production and fertility traits and having a first calving date between 2008 and 2018 were available. Unique ROH genotypes identified using a sliding-window approach were fitted into an animal mixed model as fixed effects to determine their effect on production and fertility traits. Results: In total, 133 and 34 unique ROH genotypes with unfavorable effects were identified for production and fertility traits, respectively, at a 1% genome-wise false discovery rate. Most of these ROH regions were located on bovine chromosomes 8, 13, 14 and 19 for both production and fertility traits. For production traits, the average of all the unfavorably identified unique ROH genotypes effects were estimated to decrease milk yield by 247.30 kg, fat yield by 11.46 kg and protein yield by 8.11 kg. Similarly, for fertility traits, an average 4.81-day extension in first service to conception, a 0.16 increase in number of services, and a - 0.07 incidence in 56-day non-return rate were observed. Furthermore, a ROH region located on bovine chromosome 19 was identified that, when homozygous, had a negative effect on production traits. Signatures of selection proximate to this region have implicated GH1 as a potential candidate gene, which encodes the growth hormone that binds the growth hormone receptor. This observed negative effect could be a consequence of unfavorable alleles in linkage disequilibrium with favorable alleles. Conclusions: ROH genotypes with unfavorable effects on production and fertility traits were identified within and across multiple traits on most chromosomes. These identified ROH genotypes could be included in mate selection programs to minimize their frequency in future generations

Estimates of genetic parameters for rumination time, feed efficiency, and methane production traits in first lactation Holstein cows

The large-scale recording of traits such as feed efficiency (FE) and methane emissions (ME) for use in genetic improvement programs is complex, costly, and time-consuming. Therefore, heritable traits that can be continuously recorded in dairy herds and are correlated with FE and ME traits could provide useful information for genetic evaluation. Rumination time has been suggested to be associated with FE, methane production (MeP; ME in g/d), and production traits at the phenotypic level. Therefore, the objective of this study was to investigate the genetic relationships among rumination time (RT), FE, methane and production traits using 7,358 records from 656 first-lactation Holstein cows. The estimated heritabilities were moderate for RT (0.45 ± 0.14), MeP (0.36 ± 0.12), milk yield (0.40 ± 0.08), fat yield (0.29 ± 0.06), protein yield (0.32 ± 0.07), and energy-corrected milk (0.28 ± 0.07), but were low and nonsignificant for FE (0.15 ± 0.07), which was defined as the residual of the multiple linear regression of DMI on energy-corrected milk and metabolic body weight. A favorable negative genetic correlation was estimated between RT and MeP (−0.53 ± 0.24), whereas a positive favorable correlation was estimated between RT and energy-corrected milk (0.49 ± 0.11). The estimated genetic correlation of RT with FE (−0.01 ± 0.17) was not significantly different from zero but showed a trend of a low correlation with dry matter intake (0.21 ± 0.13). These results indicate that RT is genetically associated with MeP and milk production traits, but high standard errors indicate that further analyses should be conducted to verify these findings when more data for RT, MeP, and FE become available.</p

Parent-offspring genotyped trios unravelling genomic regions with gametic and genotypic epistatic transmission bias on the cattle genome

Several biological mechanisms affecting the sperm and ova fertility and viability at developmental stages of the reproductive cycle resulted in observable transmission ratio distortion (i.e., deviation from Mendelian expectations). Gene-by-gene interactions (or epistasis) could also potentially cause specific transmission ratio distortion patterns at different loci as unfavorable allelic combinations are under-represented, exhibiting deviation from Mendelian proportions. Here, we aimed to detect pairs of loci with epistatic transmission ratio distortion using 283,817 parent-offspring genotyped trios (sire-dam-offspring) of Holstein cattle. Allelic and genotypic parameterization for epistatic transmission ratio distortion were developed and implemented to scan the whole genome. Different epistatic transmission ratio distortion patterns were observed. Using genotypic models, 7, 19 and 6 pairs of genomic regions were found with decisive evidence with additive-by-additive, additive-by-dominance/dominance-by-additive and dominance-by-dominance effects, respectively. Using the allelic transmission ratio distortion model, more insight was gained in understanding the penetrance of single-locus distortions, revealing 17 pairs of SNPs. Scanning for the depletion of individuals carrying pairs of homozygous genotypes for unlinked loci, revealed 56 pairs of SNPs with recessive epistatic transmission ratio distortion patterns. The maximum number of expected homozygous offspring, with none of them observed, was 23. Finally, in this study, we identified candidate genomic regions harboring epistatic interactions with potential biological implications in economically important traits, such as reproduction

The genetic architecture of milk ELISA scores as an indicator of Johne's disease (paratuberculosis) in dairy cattle

peer-reviewedJohne's disease (or paratuberculosis), caused by Mycobacterium avium ssp. paratuberculosis (MAP) infection, is a globally prevalent disease with severe economic and welfare implications. With no effective treatment available, understanding the role of genetics influencing host infection status is essential to develop selection strategies to breed for increased resistance to MAP infection. The main objectives of this study were to estimate genetic parameters for the MAP-specific antibody response using milk ELISA scores in Canadian Holstein cattle as an indicator of resistance to Johne's disease, and to unravel genomic regions and candidate genes significantly associated with MAP infection. After data editing, 168,987 milk ELISA records from 2,306 herds, obtained from CanWest Dairy Herd Improvement, were used for further analyses. Variance and heritability estimates for MAP infection status were determined using univariate linear animal models under 3 scenarios: (a) SCEN1: the complete data set (all herds); (b) SCEN2: herds with at least one suspect or test-positive animal (ELISA optical density ≥0.07); and (c) SCEN3: herds with at least one test-positive animal (ELISA optical density ≥0.11). Heritability estimates were calculated as 0.066, 0.064, and 0.063 for SCEN1, SCEN2, and SCEN3, respectively. The correlations between estimated breeding values for resistance to MAP infection and other economically important traits, when significant, were favorable and of low magnitude. Genome-wide association analyses identified important genomic regions on Bos taurus autosome (BTA)1, BTA7, BTA9, BTA14, BTA15, BTA17, BTA19, and BTA25 showing significant association with MAP infection status. These regions included 2 single nucleotide polymorphisms located 2 kb upstream of positional candidate genes CD86 and WNT9B, which play key roles in host immune response and tissue homeostasis. This study revealed the genetic architecture of MAP infection in Canadian Holstein cattle as measured by milk ELISA scores by estimating genetic parameters along with the identification of genomic regions potentially influencing MAP infection status. These findings will be of significant value toward implementing genetic and genomic evaluations for resistance to MAP infection in Holstein cattle

Association of genetic polymorphisms related to Johne’s disease with estimated breeding values of Holstein sires for milk ELISA test scores

peer-reviewedBackground Johne’s disease (JD) is a chronic intestinal inflammatory disease caused by Mycobacterium avium subsp. paratuberculosis (MAP) infection in ruminants. Since there are currently no effective vaccine or treatment options available to control JD, genetic selection may be an alternative strategy to enhance JD resistance. Numerous Single Nucleotide Polymorphisms (SNPs) have been reported to be associated with MAP infection status based on published genome-wide association and candidate gene studies. The main objective of this study was to validate these SNPs that were previously identified to be associated with JD by testing their effect on Holstein bulls’ estimated breeding values (EBVs) for milk ELISA test scores, an indirect indicator of MAP infection status in cattle. Results Three SNPs, rs41810662, rs41617133 and rs110225854, located on Bos taurus autosomes (BTA) 16, 23 and 26, respectively, were confirmed as significantly associated with Holstein bulls’ EBVs for milk ELISA test score (FDR < 0.01) based on General Quasi Likelihood Scoring analysis (GQLS) analysis. Single-SNP regression analysis identified four SNPs that were associated with sire EBVs (FDR < 0.05). This includes two SNPs that were common with GQLS (rs41810662 and rs41617133), with the other two SNPs being rs110494981 and rs136182707, located on BTA9 and BTA16, respectively. Conclusions The findings of this study validate the association of SNPs with JD MAP infection status and highlight the need to further investigate the genomic regions harboring these SNPs