Adding gene transcripts into genomic prediction improves accuracy and reveals sampling time dependence.

Recent developments allowed generating multiple high-quality 'omics' data that could increase the predictive performance of genomic prediction for phenotypes and genetic merit in animals and plants. Here, we have assessed the performance of parametric and nonparametric models that leverage transcriptomics in genomic prediction for 13 complex traits recorded in 478 animals from an outbred mouse population. Parametric models were implemented using the best linear unbiased prediction, while nonparametric models were implemented using the gradient boosting machine algorithm. We also propose a new model named GTCBLUP that aims to remove between-omics-layer covariance from predictors, whereas its counterpart GTBLUP does not do that. While gradient boosting machine models captured more phenotypic variation, their predictive performance did not exceed the best linear unbiased prediction models for most traits. Models leveraging gene transcripts captured higher proportions of the phenotypic variance for almost all traits when these were measured closer to the moment of measuring gene transcripts in the liver. In most cases, the combination of layers was not able to outperform the best single-omics models to predict phenotypes. Using only gene transcripts, the gradient boosting machine model was able to outperform best linear unbiased prediction for most traits except body weight, but the same pattern was not observed when using both single nucleotide polymorphism genotypes and gene transcripts. Although the GTCBLUP model was not able to produce the most accurate phenotypic predictions, it showed the highest accuracies for breeding values for 9 out of 13 traits. We recommend using the GTBLUP model for prediction of phenotypes and using the GTCBLUP for prediction of breeding values.

Prediction performance of linear models and gradient boosting machine on complex phenotypes in outbred mice.

We compared the performance of linear (GBLUP, BayesB, and elastic net) methods to a nonparametric tree-based ensemble (gradient boosting machine) method for genomic prediction of complex traits in mice. The dataset used contained genotypes for 50,112 SNP markers and phenotypes for 835 animals from 6 generations. Traits analyzed were bone mineral density, body weight at 10, 15, and 20 weeks, fat percentage, circulating cholesterol, glucose, insulin, triglycerides, and urine creatinine. The youngest generation was used as a validation subset, and predictions were based on all older generations. Model performance was evaluated by comparing predictions for animals in the validation subset against their adjusted phenotypes. Linear models outperformed gradient boosting machine for 7 out of 10 traits. For bone mineral density, cholesterol, and glucose, the gradient boosting machine model showed better prediction accuracy and lower relative root mean squared error than the linear models. Interestingly, for these 3 traits, there is evidence of a relevant portion of phenotypic variance being explained by epistatic effects. Using a subset of top markers selected from a gradient boosting machine model helped for some of the traits to improve the accuracy of prediction when these were fitted into linear and gradient boosting machine models. Our results indicate that gradient boosting machine is more strongly affected by data size and decreased connectedness between reference and validation sets than the linear models. Although the linear models outperformed gradient boosting machine for the polygenic traits, our results suggest that gradient boosting machine is a competitive method to predict complex traits with assumed epistatic effects.

Dual-purpose Guzerá cattle exhibit high dairy performance under heat stress.

The present study evaluated the heat stress response pattern of dual-purpose Guzerá cattle for test-day (TD) milk yield records of first lactation and estimated genetic parameters and trends related to heat stress. A total of 31,435 TD records from 4,486 first lactations of Guzerá cows, collected between 1986 and 2012, were analysed. Two random regression models considered days in milk (DIM) and/or temperature × humidity-dependent (THI) covariate. Impacts of -0.037, -0.019 and -0.006 kg/day/THI for initial and intermediate stages of lactation were observed when considering the mean maximum daily temperature and humidity to calculate THI. Heritability estimates ranged from 0.16 to 0.35 throughout lactation and THI values, suggesting the possibility to expect gains from selection for such trait. The variable trajectory of breeding values for dual-purpose Guzerá sires in response to changes in THI values confirms that the genotype × environment interaction due to heat stress can have some effect on TD milk yield. Despite the high dairy performance of Guzerá cattle under heat stress, estimated genetic trends showed a progressive reduction in heat tolerance. Therefore, new strategies should be adopted to prevent negative impacts of heat stress over milk production in Guzerá animals in future.

Fine mapping of genomic regions associated with female fertility in Nellore beef cattle based on sequence variants from segregating sires.

BACKGROUND: Impaired fertility in cattle limits the efficiency of livestock production systems. Unraveling the genetic architecture of fertility traits would facilitate their improvement by selection. In this study, we characterized SNP chip haplotypes at QTL blocks then used whole-genome sequencing to fine map genomic regions associated with reproduction in a population of Nellore (Bos indicus) heifers. METHODS: The dataset comprised of 1337 heifers genotyped using a GeneSeek® Genomic Profiler panel (74677 SNPs), representing the daughters from 78 sires. After performing marker quality control, 64800 SNPs were retained. Haplotypes carried by each sire at six previously identified QTL on BTAs 5, 14 and 18 for heifer pregnancy and BTAs 8, 11 and 22 for antral follicle count were constructed using findhap software. The significance of the contrasts between the effects of every two paternally-inherited haplotype alleles were used to identify sires that were heterozygous at each QTL. Whole-genome sequencing data localized to the haplotypes from six sires and 20 other ancestors were used to identify sequence variants that were concordant with the haplotype contrasts. Enrichment analyses were applied to these variants using KEGG and MeSH libraries. RESULTS: A total of six (BTA 5), six (BTA 14) and five (BTA 18) sires were heterozygous for heifer pregnancy QTL whereas six (BTA 8), fourteen (BTA 11), and five (BTA 22) sires were heterozygous for number of antral follicles' QTL. Due to inadequate representation of many haplotype alleles in the sequenced animals, fine mapping analysis could only be reliably performed for the QTL on BTA 5 and 14, which had 641 and 3733 concordant candidate sequence variants, respectively. The KEGG "Circadian rhythm" and "Neurotrophin signaling pathway" were significantly associated with the genes in the QTL on BTA 5 whereas 32 MeSH terms were associated with the QTL on BTA 14. Among the concordant sequence variants, 0.2% and 0.3% were classified as missense variants for BTAs 5 and 14, respectively, highlighting the genes MTERF2, RTMB, ENSBTAG00000037306 (miRNA), ENSBTAG00000040351, PRKDC, and RGS20. The potential causal mutations found in the present study were associated with biological processes such as oocyte maturation, embryo development, placenta development and response to reproductive hormones. CONCLUSIONS: The identification of heterozygous sires by positionally phasing SNP chip data and contrasting haplotype effects for previously detected QTL can be used for fine mapping to identify potential causal mutations and candidate genes. Genomic variants on genes MTERF2, RTBC, miRNA ENSBTAG00000037306, ENSBTAG00000040351, PRKDC, and RGS20, which are known to have influence on reproductive biological processes, were detected.

Accounting for population structure in selective cow genotyping strategies.

The objective of the present study was to investigate the impact of considering population structure in cow genotyping strategies over the accuracy and bias of genomic predictions. A small dairy cattle population was simulated to address these objectives. Based on four main traditional designs (random, top-yield, extreme-yield and top-accuracy cows), different numbers (1,000; 2,000 and 5,000) of cows were sampled and included in the reference population. Traditional designs were replicated considering or not population structure and compared among and with a reference population containing only bulls. The inclusion of cows increased accuracy in all scenarios compared with using only bulls. Scenarios accounting for population structure when choosing cows to the reference population slightly outperformed their traditional versions by yielding higher accuracy and lower bias in genomic predictions. Building a cow-based reference population from groups of related individuals considering the frequency of individuals from those same groups in the validation population yielded promising results with applications on selection for expensive- or difficult-to-measure traits. Methods here presented may be easily implemented in both new or already established breeding programs, as they improved prediction and reduced bias in genomic evaluations while demanding no additional costs.