Multi-breed genomic predictions and functional variants for fertility of tropical bulls.

Worldwide, most beef breeding herds are naturally mated. As such, the ability to identify and select fertile bulls is critically important for both productivity and genetic improvement. Here, we collected ten fertility-related phenotypes for 6,063 bulls from six tropically adapted breeds. Phenotypes were comprised of four bull conformation traits and six traits directly related to the quality of the bull's semen. We also generated high-density DNA genotypes for all the animals. In total, 680,758 single nucleotide polymorphism (SNP) genotypes were analyzed. The genomic correlation of the same trait observed in different breeds was positive for scrotal circumference and sheath score on most breed comparisons, but close to zero for the percentage of normal sperm, suggesting a divergent genetic background for this trait. We confirmed the importance of a breed being present in the reference population to the generation of accurate genomic estimated breeding values (GEBV) in an across-breed validation scenario. Average GEBV accuracies varied from 0.19 to 0.44 when the breed was not included in the reference population. The range improved to 0.28 to 0.59 when the breed was in the reference population. Variants associated with the gene HDAC4, six genes from the spermatogenesis-associated (SPATA) family of proteins, and 29 transcription factors were identified as candidate genes. Collectively these results enable very early in-life selection for bull fertility traits, supporting genetic improvement strategies currently taking place within tropical beef production systems. This study also improves our understanding of the molecular basis of male fertility in mammals.

Utility of multi-omics data to inform genomic prediction of heifer fertility traits.

Biologically informed single nucleotide polymorphisms (SNPs) impact genomic prediction accuracy of the target traits. Our previous genomics, proteomics, and transcriptomics work identified candidate genes related to puberty and fertility in Brahman heifers. We aimed to test this biological information for capturing heritability and predicting heifer fertility traits in another breed i.e., Tropical Composite. The SNP from the identified genes including 10 kilobases (kb) region on either side were selected as biologically informed SNP set. The SNP from the rest of the Bos taurus genes including 10-kb region on either side were selected as biologically uninformed SNP set. Bovine high-density (HD) complete SNP set (628,323 SNP) was used as a control. Two populations-Tropical Composites (N = 1331) and Brahman (N = 2310)-had records for three traits: pregnancy after first mating season (PREG1, binary), first conception score (FCS, score 1 to 3), and rebreeding score (REB, score 1 to 3.5). Using the best linear unbiased prediction method, effectiveness of each SNP set to predict the traits was tested in two scenarios: a 5-fold cross-validation within Tropical Composites using biological information from Brahman studies, and application of prediction equations from one breed to the other. The accuracy of prediction was calculated as the correlation between genomic estimated breeding values and adjusted phenotypes. Results show that biologically informed SNP set estimated heritabilities not significantly better than the control HD complete SNP set in Tropical Composites; however, it captured all the observed genetic variance in PREG1 and FCS when modeled together with the biologically uninformed SNP set. In 5-fold cross-validation within Tropical Composites, the biologically informed SNP set performed marginally better (statistically insignificant) in terms of prediction accuracies (PREG1: 0.20, FCS: 0.13, and REB: 0.12) as compared to HD complete SNP set (PREG1: 0.17, FCS: 0.10, and REB: 0.11), and biologically uninformed SNP set (PREG1: 0.16, FCS: 0.10, and REB: 0.11). Across-breed use of prediction equations still remained a challenge: accuracies by all SNP sets dropped to around zero for all traits. The performance of biologically informed SNP was not significantly better than other sets in Tropical Composites. However, results indicate that biological information obtained from Brahman was successful to predict the fertility traits in Tropical Composite population.

Prior biological information can be helpful in the genomic prediction of the traits. Previous multi-omics studies by our group identified genes relevant to puberty and fertility in Brahman cattle, a beef breed in northern Australia. We used this biological information in the genomic prediction of three heifer fertility traits, measured in another beef cattle breed: Tropical Composites. The three traits were: pregnancy status after the first mating season (PREG1), first conception score (FCS), and rebreeding score (REB). To test if prior biological information could capture genetic variation in the traits and improve genomic predictions, we compared the results obtained using three subsets of genetic information (i.e., subsets of DNA variants). The first subset contained only variants deemed biologically relevant (as per previous multi-omics studies). The second subset contained only variants considered biologically irrelevant. The third subset had all the variants contained in the commercial DNA assay known as the bovine high-density chip, intended as a practical control. The results indicate that multi-omics data was informative across breed scenario and can be useful in informing genomic predictions of traits of interest.

Meta-Analysis of Heifer Traits Identified Reproductive Pathways in Bos indicus Cattle.

Fertility traits measured early in life define the reproductive potential of heifers. Knowledge of genetics and biology can help devise genomic selection methods to improve heifer fertility. In this study, we used ~2400 Brahman cattle to perform GWAS and multi-trait meta-analysis to determine genomic regions associated with heifer fertility. Heifer traits measured were pregnancy at first mating opportunity (PREG1, a binary trait), first conception score (FCS, score 1 to 3) and rebreeding score (REB, score 1 to 3.5). The heritability estimates were 0.17 (0.03) for PREG1, 0.11 (0.05) for FCS and 0.28 (0.05) for REB. The three traits were highly genetically correlated (0.75-0.83) as expected. Meta-analysis was performed using SNP effects estimated for each of the three traits, adjusted for standard error. We identified 1359 significant SNPs (p-value < 9.9 × 10-6 at FDR < 0.0001) in the multi-trait meta-analysis. Genomic regions of 0.5 Mb around each significant SNP from the meta-analysis were annotated to create a list of 2560 positional candidate genes. The most significant SNP was in the vicinity of a genomic region on chromosome 8, encompassing the genes SLC44A1, FSD1L, FKTN, TAL2 and TMEM38B. The genomic region in humans that contains homologs of these genes is associated with age at puberty in girls. Top significant SNPs pointed to additional fertility-related genes, again within a 0.5 Mb region, including ESR2, ITPR1, GNG2, RGS9BP, ANKRD27, TDRD12, GRM1, MTHFD1, PTGDR and NTNG1. Functional pathway enrichment analysis resulted in many positional candidate genes relating to known fertility pathways, including GnRH signaling, estrogen signaling, progesterone mediated oocyte maturation, cAMP signaling, calcium signaling, glutamatergic signaling, focal adhesion, PI3K-AKT signaling and ovarian steroidogenesis pathway. The comparison of results from this study with previous transcriptomics and proteomics studies on puberty of the same cattle breed (Brahman) but in a different population identified 392 genes in common from which some genes-BRAF, GABRA2, GABR1B, GAD1, FSHR, CNGA3, PDE10A, SNAP25, ESR2, GRIA2, ORAI1, EGFR, CHRNA5, VDAC2, ACVR2B, ORAI3, CYP11A1, GRIN2A, ATP2B3, CAMK2A, PLA2G, CAMK2D and MAPK3-are also part of the above-mentioned pathways. The biological functions of the positional candidate genes and their annotation to known pathways allowed integrating the results into a bigger picture of molecular mechanisms related to puberty in the hypothalamus-pituitary-ovarian axis. A reasonable number of genes, common between previous puberty studies and this study on early reproductive traits, corroborates the proposed molecular mechanisms. This study identified the polymorphism associated with early reproductive traits, and candidate genes that provided a visualization of the proposed mechanisms, coordinating the hypothalamic, pituitary, and ovarian functions for reproductive performance in Brahman cattle.