Genome-wide association studies for epistatic genetic effects on fertility and reproduction traits in Holstein cattle.

Non-additive genetic effects are well known to play an important role in the phenotypic expression of complex traits, such as fertility and reproduction. In this study, a genome scan was performed using 41,640 single nucleotide polymorphism (SNP) markers to identify genomic regions associated with epistatic (additive-by-additive) effects in fertility and reproduction traits in Holstein cattle. Nine fertility and reproduction traits were analysed on 5825 and 6090 Holstein heifers and cows with phenotypes and genotypes, respectively. The Marginal Epistasis Test (MAPIT) was used to identify SNPs with significant marginal epistatic effects at a chromosome-wise 5% and 10% false discovery rate (FDR) level. The -log10 (p) values were adjusted by the genomic inflation factor (λ) to correct for the potential bias on the p-values and minimize the possible effects of population stratification. After adjustments, MAPIT enabled the identification of genomic regions with significant marginal epistatic effects for heifers on BTA5 for age at first insemination, BTA3 and BTA24 for non-return rate (NRR); BTA16 and BTA28 for gestation length (GL); BTA1, BTA4 and BTA17 for stillbirth (SB). For the cow traits, MAPIT enabled the identification of regions on BTA11 for GL, BTA11 and BTA16 for SB and BTA19 for calf size (CZ). An additional approach for mapping epistasis in a genome-wide association study was also proposed, in which the genome scan was performed using estimates of epistatic values as the input pseudo-phenotypes, computed using single-trait animal models. Significant SNPs were identified at the chromosome-wise 5% and 10% FDR levels for all traits. For the heifer traits, significant regions were found on BTA7 for AFS; BTA12 for NRR; BTA14 and BTA19 for GL; BTA19 for calving ease (CE); BTA5, BTA24, BTA25 and in the X chromosome for SB; BTA23 and in the X chromosome for CZ and in the X chromosome for the number of services (NS). For the cow traits, significant regions were found on BTA29 and in the X chromosome for NRR, BTA11, BTA16 and in the X chromosome for SB, BTA2 for GL, BTA28 for CZ, BTA19 for calving to first insemination, and in the X chromosome for NS and first insemination to conception. The results suggest that the epistatic genetic effects are likely due to many loci with a small effect rather than few loci with a large effect and/or a single SNP marker alone do not capture the epistatic effects well. The genomic architecture of fertility and reproduction traits is complex, and these results should be validated in independent dairy cattle populations and using alternative statistical models.

Genomic prediction of fertility and calving traits in Holstein cattle based on models including epistatic genetic effects.

The goal of this study was to investigate whether the inclusion of genomic information and epistatic (additive by additive) genetic effects would increase the accuracy of predicting phenotypes adjusted for known environmental effects, reduce prediction bias and minimize the confounding between additive and additive by additive epistatic effects on fertility and calving traits in Holstein cattle. Phenotypic and genotypic records were available for 6090 cows. Eight cow traits were assessed including 56-day nonreturn rate (NRR), number of services (NS), calving to first insemination (CTFS), first insemination to conception (FSTC), gestation length (GL), calving ease (CE), stillbirth (SB) and calf size (CZ). Four scenarios were assessed for their ability to predict adjusted phenotypes, which included: (1) traditional pedigree-based Best Linear Unbiased Prediction (P-BLUP) for additive genetic effects (PA); (2) P-BLUP for additive and epistatic (additive by additive) genetic effects (PAE); (3) genomic BLUP (G-BLUP) for additive genetic effects (GA); and (4) G-BLUP for additive and epistatic genetic effects (GAEn, where n = 1-3 depending on the alternative ways to construct the epistatic genomic matrix used). Constructing epistatic relationship matrix as the Hadamard product of the additive genomic relationship matrix (GAE1), which is the usual method and implicitly assumes a model that fits all pairwise interactions between markers twice and includes the interactions of the markers with themselves (dominance). Two additional constructions of the epistatic genomic relationship matrix were compared to test whether removing the double counting of interactions and the interaction of the markers with themselves (GAE2), and removing double counting of interactions between markers, but including the interaction of the markers with themselves (GAE3) would had an impact on the prediction and estimation error correlation (i.e. confounding) between additive and epistatic genetic effects. Fitting epistatic genetic effects explained up to 5.7% of the variance for NRR (GAE3), 7.7% for NS (GAE1), 11.9% for CTFS (GAE3), 11.1% for FSTC (GAE2), 25.7% for GL (GAE1), 2.3% for CE (GAE1), 14.3% for SB (GAE3) and 15.2% for CZ (GAE1). Despite a substantial proportion of variance being explained by epistatic effects for some traits, the prediction accuracies were similar or lower for GAE models compared with pedigree models and genomic models without epistatic effects. Although the prediction accuracy of direct genomic values did not change significantly between the three variations of the epistatic genetic relationship matrix used, removing the interaction of the markers with themselves reduced the confounding between additive and additive by additive epistatic effects. These results suggest that epistatic genetic effects contribute to the variance of some fertility and calving traits in Holstein cattle. However, the inclusion of epistatic genetic effects in the genomic prediction of these traits is complex and warrant further investigation.

Estimation of additive and non-additive genetic effects for fertility and reproduction traits in North American Holstein cattle using genomic information.

Non-additive genetic effects are usually ignored in animal breeding programs due to data structure (e.g., incomplete pedigree), computational limitations and over-parameterization of the models. However, non-additive genetic effects may play an important role in the expression of complex traits in livestock species, such as fertility and reproduction traits. In this study, components of genetic variance for additive and non-additive genetic effects were estimated for a variety of fertility and reproduction traits in Holstein cattle using pedigree and genomic relationship matrices. Four linear models were used: (a) an additive genetic model; (b) a model including both additive and epistatic (additive by additive) genetic effects; (c) a model including both additive and dominance effects; and (d) a full model including additive, epistatic and dominance genetic effects. Nine fertility and reproduction traits were analysed, and models were run separately for heifers (N = 5,825) and cows (N = 6,090). For some traits, a larger proportion of phenotypic variance was explained by non-additive genetic effects compared with additive effects, indicating that epistasis, dominance or a combination thereof is of great importance. Epistatic genetic effects contributed more to the total phenotypic variance than dominance genetic effects. Although these models varied considerably in the partitioning of the components of genetic variance, the models including a non-additive genetic effect did not show a clear advantage over the additive model based on the Akaike information criterion. The partitioning of variance components resulted in a re-ranking of cows based solely on the cows' additive genetic effects between models, indicating that adjusting for non-additive genetic effects could affect selection decisions made in dairy cattle breeding programs. These results suggest that non-additive genetic effects play an important role in some fertility and reproduction traits in Holstein cattle.

Estimation of direct and maternal genetic parameters for individual birth weight, weaning weight, and probe weight in Yorkshire and Landrace pigs.

As a result of selecting for increased litter size, newborn piglets are being born lighter and have a lower chance of survival. Raising fewer pigs to market weight would have a negative impact on the industry and farmer profitability; thus, understanding the genetics of individual growth performance traits will determine whether these traits will play an important role in pig breeding schemes. This study aimed to estimate genetic parameters for individual birth weight (BW), weaning weight (WW), and probe weight (PW) in Canadian-purebred Yorkshire and Landrace pigs. PW is a live weight taken at the time of the ultrasound measurements, when pigs weigh about 100 kg. Data were collected from 2 large and related breeding herds from 2003 to 2015. Four linear animal models were used, which included the following: Model 1-direct additive genetic effect; Model 2-direct additive genetic and maternal genetic effect; Model 3-direct additive genetic and common litter effect; and Model 4-direct additive genetic, maternal genetic, and common litter effect. The model which included all 3 random effects (Model 4) was determined to be the best fit to the data. Low to moderate direct heritability estimates were observed as follows: 0.15 ± 0.03 for BW, 0.04 ± 0.01 for WW, and 0.33 ± 0.03 for PW for the Yorkshire breed; and 0.05 ± 0.01 for BW, 0.01 ± 0.01 for WW, and 0.27 ± 0.03 for PW in the Landrace breed. As expected, the direct heritability estimates increased with age as a result of decreased maternal influence on the trait. Bivariate animal models were also used to estimate genetic and environmental correlations between traits. Strong direct genetic correlations were observed between BW and WW in both breeds. Based on the estimates of genetic parameters, individual BW could be evaluated and considered in breeding programs aiming to increase BW and improve subsequent performance. Different selection emphasis could also be applied on direct and maternal additive genetic effects on BW to optimize the breeding programs and improve selection efficiency.