Genomic prediction of bull fertility in US Jersey dairy cattle

Integrating Genomic Selection and a Genome-Wide Association Study to Improve Days Open in Thai Dairy Holstein Cattle: A Comprehensive Genetic Analysis

Days open (DO) is a critical economic and reproductive trait that is commonly employed in genetic selection. Making improvements using conventional genetic techniques is exceedingly challenging. Therefore, new techniques are required to improve the accuracy of genetic selection using genomic data. This study examined the genetic approaches of traditional AIREML and single-step genomic AIREML (ssGAIREML) to assess genetic parameters and the accuracy of estimated breeding values while also investigating SNP regions associated with DO and identifying candidate genes through a genome-wide association study (GWAS). The dataset included 59415 DO records from 36368 Thai-Holstein crossbred cows and 882 genotyped animals. The cows were classified according to their Holstein genetic proportion (breed group, BG) as follows: BG1 (>93.7% Holstein genetics), BG2 (87.5% to 93.6% Holstein genetics), and BG3 (<87.5% Holstein genetics). AIREML was utilized to estimate genetic parameters and variance components. The results of this study reveal that the average DO values for BG1, BG2, and BG3 were 97.64, 97.25, and 96.23 days, respectively. The heritability values were estimated to be 0.02 and 0.03 for the traditional AIREML and ssGAIREML approaches, respectively. Depending on the dataset, the ssGAIREML method produced more accurate estimated breeding values than the traditional AIREML method, ranging from 40.5 to 45.6%. The highest values were found in the top 20% of the dam dataset. For the GWAS, we found 12 potential candidate genes (DYRK1A, CALCR, MIR489, MIR653, SLC36A1, GNA14, GNAQ, TRNAC-GCA, XYLB, ACVR2B, SLC22A14, and EXOC2) that are believed to have a significant influence on days open. In summary, the ssGAIREML method has the potential to enhance the accuracy and heritability of reproductive values compared to those obtained using conventional AIREML. Consequently, it is a viable alternative for transitioning from conventional methodologies to the ssGAIREML method in the breeding program for dairy cattle in Thailand. Moreover, the 12 identified potential candidate genes can be utilized in future studies to select markers for days open in regard to dairy cattle.

Genomic prediction of male fertility in Brown Swiss cattle

Bull fertility has been recognized as an important factor affecting dairy herd fertility. The objective of this study was to assess the feasibility of predicting male fertility in Brown Swiss cattle using genomic data. The dataset consisted of 1,102 Italian Brown Swiss bulls with sire conception rate (SCR) records and genotype data for roughly 480k SNP. The analyses included the use of linear kernel-based regression models fitting all SNPs or incorporating markers with large effect. Predictive performance was evaluated in 5-fold cross-validation using the correlation between observed and predicted SCR values and mean squared error of prediction. The entire SNP set exhibited predictive correlations around 0.19. Interestingly, the inclusion of 2 markers with large effect yielded predictive correlations around 0.32. Overall, using linear kernel-based models fitting markers with large effect is a promising approach. Our findings could help Brown Swiss breeders make enhanced genome-guided management and selection decisions on male fertility.

Genomics and Dairy Bull Fertility

Current evidence suggests that dairy bull fertility is influenced by genetic factors, and hence, it could be managed and improved by genetic means. There are major mutations that explain about 4% to 8% of the observed differences in conception rate between bulls segregating in most dairy breeds. Research has shown that genomic prediction of bull fertility is possible, and this could be used to make accurate genome-guided selection decisions, such as early culling of predicted subfertile bull calves. Inbreeding negatively influences bull fertility, and the increase in homozygosity seems an important risk factor for dairy bull subfertility.

Harnessing male germline epigenomics for the genetic improvement in cattle

Sperm is essential for successful artificial insemination in dairy cattle, and its quality can be influenced by both epigenetic modification and epigenetic inheritance. The bovine germline differentiation is characterized by epigenetic reprogramming, while intergenerational and transgenerational epigenetic inheritance can influence the offspring's development through the transmission of epigenetic features to the offspring via the germline. Therefore, the selection of bulls with superior sperm quality for the production and fertility traits requires a better understanding of the epigenetic mechanism and more accurate identifications of epigenetic biomarkers. We have comprehensively reviewed the current progress in the studies of bovine sperm epigenome in terms of both resources and biological discovery in order to provide perspectives on how to harness this valuable information for genetic improvement in the cattle breeding industry.

Genetic parameters for various semen production and quality traits and indicators of male and female reproductive performance in Nellore cattle

BACKGROUND

Given the economic relevance of fertility and reproductive traits for the beef cattle industry, investigating their genetic background and developing effective breeding strategies are paramount. Considering their late and sex-dependent phenotypic expression, genomic information can contribute to speed up the rates of genetic progress per year. In this context, the main objectives of this study were to estimate variance components and genetic parameters, including heritability and genetic correlations, for fertility, female precocity, and semen production and quality (andrological attributes) traits in Nellore cattle incorporating genomic information.

RESULTS

The heritability estimates of semen quality traits were low-to-moderate, while moderate-to-high estimates were observed for semen morphological traits. The heritability of semen defects ranged from low (0.04 for minor semen defects) to moderate (0.30 for total semen defects). For seminal aspect (SMN_ASPC) and bull reproductive fitness (BULL_FIT), low (0.19) and high (0.69) heritabilities were observed, respectively. The heritability estimates for female reproductive traits ranged from 0.16 to 0.39 for rebreeding of precocious females (REBA) and probability of pregnancy at 14 months (PP14), respectively. Semen quality traits were highly genetically correlated among themselves. Moderate-to-high genetic correlations were observed between the ability to remain productive in the herd until four years of age (stayability; STAY) and the other reproductive traits, indicating that selection for female reproductive performance will indirectly contribute to increasing fertility rates. High genetic correlations between BULL_FIT and female reproductive traits related to precocity (REBA and PP14) and STAY were observed. The genetic correlations between semen quality and spermatic morphology with female reproductive traits ranged from -0.22 (REBA and scrotal circumference) to 0.48 (REBA and sperm vigor). In addition, the genetic correlations between REBA with semen quality traits ranged from -0.23 to 0.48, and with the spermatic morphology traits it ranged from -0.22 to 0.19.

CONCLUSIONS

All male and female fertility and reproduction traits evaluated are heritable and can be improved through direct genetic or genomic selection. Selection for better sperm quality will positively influence the fertility and precocity of Nellore females. The findings of this study will serve as background information for designing breeding programs for genetically improving semen production and quality and reproductive performance in Nellore cattle.

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.

Prediction ability for growth and maternal traits using SNP arrays based on different marker densities in Nellore cattle using the ssGBLUP

This study aimed to investigate the prediction ability for growth and maternal traits using different low-density customized SNP arrays selected by informativeness and distribution of markers across the genome employing single-step genomic BLUP (ssGBLUP). Phenotypic records for adjusted weight at 210 and 450 days of age were utilized. A total of 945 animals were genotyped with high-density chip, and 267 individuals born after 2008 were selected as validation population. We evaluated 11 scenarios using five customized density arrays (40 k, 20 k, 10 k, 5 k and 2 k) and the HD array was used as desirable scenario. The GEBV predictions and BIF (Beef Improvement Federation) accuracy were obtained with BLUPF90 family programs. Linear regression was used to evaluate the prediction ability, inflation, and bias of GEBV of each customized array. An overestimation of partial GEBVs in contrast with complete GEBVs and increase of BIF accuracy with the density arrays diminished were observed. For all traits, the prediction ability was higher as the array density increased and it was similar with customized arrays higher than 10 k SNPs. Level of inflation was lower as the density array increased of and was higher for MW210 effect. The bias was susceptible to overestimation of GEBVs when the density customized arrays decreased. These results revealed that the BIF accuracy is sensible to overestimation using low-density customized arrays while the prediction ability with least 10,000 informative SNPs obtained from the Illumina BovineHD BeadChip shows accurate and less biased predictions. Low-density customized arrays under ssGBLUP method could be feasible and cost-effective in genomic selection.

The future of assessing bull fertility: Can the 'omics fields identify usable biomarkers?

Breeding soundness examinations for bulls rely heavily on the subjective, visual assessment of sperm motility and morphology. Although these criteria have the potential to identify infertile males, they cannot be used to guarantee fertility or provide information about varying degrees of bull fertility. Male factor fertility is complex, and the success of the male gamete is not necessarily realized until well after the spermatozoon enters the oocyte. This paper reviews our existing knowledge of the bull’s contribution from a standpoint of the sperm’s cargo and the impact that this can have on fertilization and the development of the embryo. There has been a plethora of recent research characterizing the many molecular attributes that can affect the functional competence of a spermatozoon. A better understanding of the molecular factors influencing fertilization and embryo development in cattle will lead to the identification of biomarkers for the selection of bulls of superior fertility, which will have major implications for livestock production. To see this improvement in reproductive performance, we believe incorporation of modern technology into breeding soundness examinations will be necessary—although many of the discussed technologies are not ready for large-scale field application. Each of the ‘omics fields discussed in this review have shown promise for the identification of biomarkers of fertility, with certain families of biomarkers appearing to be better suited to different evaluations throughout a bull’s lifetime. Further research is needed for the proposed biomarkers to be of diagnostic or predictive value.

Evaluation of bull fertility in Italian Brown Swiss dairy cattle using cow field data

Dairy bull fertility is traditionally evaluated using semen production and quality traits; however, these attributes explain only part of the differences observed in fertility among bulls. Alternatively, bull fertility can be directly evaluated using cow field data. The main objective of this study was to investigate bull fertility in the Italian Brown Swiss dairy cattle population using confirmed pregnancy records. The data set included a total of 397,926 breeding records from 1,228 bulls and 129,858 lactating cows between first and fifth lactation from 2000 to 2019. We first evaluated cow pregnancy success, including factors related to the bull under evaluation, such as bull age, bull inbreeding, and AI organization, and factors associated with the cow that receives the dose of semen, including herd-year-season, cow age, parity, and milk yield. We then estimated sire conception rate using only factors related to the bull. Model predictive ability was evaluated using 10-fold cross-validation with 10 replicates. Interestingly, our analyses revealed that there is a substantial variation in conception rate among Brown Swiss bulls, with more than 20% conception rate difference between high-fertility and low-fertility bulls. We also showed that the prediction of bull fertility is feasible as our cross-validation analyses achieved predictive correlations equal to 0.30 for sire conception rate. Improving reproduction performance is one of the major challenges of the dairy industry worldwide, and for this, it is essential to have accurate predictions of service sire fertility. This study represents the foundation for the development of novel tools that will allow dairy producers, breeders, and artificial insemination companies to make enhanced management and selection decisions on Brown Swiss male fertility.

Marker-assisted selection vis-à-vis bull fertility: coming full circle-a review

Bull fertility is considered an indispensable trait, as far as farm economics is concerned since it is the successful conception in a cow that provides calf crop, along with the ensuing lactation. This ensures sustainability of a dairy farm. Traditionally, bull fertility did not receive much attention by the farm managers and breeding animals were solely evaluated based on phenotypic predictors, namely, sire conception rate and seminal parameters in bull. With the advent of the molecular era in animal breeding, attempts were made to unravel the genetic complexity of bull fertility by the identification of genetic markers related to the trait. Marker-Assisted Selection (MAS) is a methodology that aims at utilizing the genetic information at markers and selecting improved populations for important traits. Traditionally, MAS was pursued using a candidate gene approach for identifying markers related to genes that are already known to have a physiological function related to the trait but this approach had certain shortcomings like stringent criteria for significance testing. Now, with the availability of genome-wide data, the number of markers identified and variance explained in relation to bull fertility has gone up. So, this presents a unique opportunity to revisit MAS by selection based on the information of a large number of genome-wide markers and thus, improving the accuracy of selection.

Across-country genomic prediction of bull fertility in Jersey dairy cattle

The use of information across populations is an attractive approach to increase the accuracy of genomic predictions for numerically small breeds and traits that are time-consuming and difficult to measure, such as male fertility in cattle. This study was conducted to evaluate genomic prediction of Jersey bull fertility using an across-country reference population combining records from the United States and Australia. The data set consisted of 1,570 US Jersey bulls with sire conception rate (SCR) records, 603 Australian Jersey bulls with semen fertility value (SFV) records and SNP genotypes for roughly 90,000 loci. Both SCR and SFV are evaluations of service sire fertility based on cow field data, and both are intended as phenotypic evaluations because the estimates include genetic and nongenetic effects. Within- and across-country genomic predictions were evaluated using univariate and bivariate genomic best linear unbiased prediction models. Predictive ability was assessed in 5-fold cross-validation using the correlation between observed and predicted fertility values and mean squared error of prediction. Within-country genomic predictions exhibited predictive correlations of around 0.28 and 0.02 for the United States and Australia, respectively. The Australian Jersey population is genetically diverse and small in size, so careful selection of the reference population by including only closely related animals (e.g., excluding New Zealand bulls, which is a less-related population) increased the predictive correlations up to 0.20. Notably, the use of bivariate models fitting all US Jersey records and the optimized Australian population resulted in predictive correlations around of 0.24 for SFV values, which is a relative increase in predictive ability of 20%. Conversely, for predicting SCR values, the use of an across-country reference population did not outperform the standard approach using pure US Jersey reference data set. Our findings indicate that genomic prediction of male fertility in dairy cattle is feasible, and the use of an across-country reference population would be beneficial when local populations are small and genetically diverse.

Optimizing genomic prediction model given causal genes in a dairy cattle population

As genotypic data are moving from SNP chip toward whole-genome sequence, the accuracy of genomic prediction (GP) exhibits a marginal gain, although all genetic variation, including causal genes, are contained in whole-genome sequence data. Meanwhile, genetic analyses on complex traits, such as genome-wide association studies, have identified an increasing number of genomic regions, including potential causal genes, which would be reliable prior knowledge for GP. Many studies have tried to improve the performance of GP by modifying the prediction model to incorporate prior knowledge. Although several plausible results have been obtained from model modification or strategy optimization, most of them were validated in a specific empirical population with a limited variety of genetic architecture for complex traits. An alternative approach is to use simulated genetic architecture with known causal genes (e.g., simulated causative SNP) to evaluate different GP models with given causal genes. Our objectives were to (1) evaluate the performance of GP under a variety of genetic architectures with a subset of known causal genes and (2) compare different GP models modified by highlighting causal genes and different strategies to weight causal genes. In this study, we simulated pseudo-phenotypes under a variety of genetic architectures based on the real genotypes and phenotypes of a dairy cattle population. Besides classical genomic best linear unbiased prediction, we evaluated 3 modified GP models that highlight causal genes as follows: (1) by treating them as fixed effects, (2) by treating them as a separate random component, and (3) by combining them into the genomic relationship matrix as random effects. Our results showed that highlighting the known causal genes, which explained a considerable proportion of genetic variance in the GP models, increased the predictive accuracy. Combining all given causal genes into the genomic relationship matrix was the optimal strategy under all the scenarios validated, and treating causal genes as a separate random component is also recommended, when more than 20% of genetic variance was explained by known causal genes. Moreover, assigning differential weights to each causal gene further improved the predictive accuracy.

Activation of cryptic splicing in bovine WDR19 is associated with reduced semen quality and male fertility

Cattle are ideally suited to investigate the genetics of male reproduction, because semen quality and fertility are recorded for all ejaculates of artificial insemination bulls. We analysed 26,090 ejaculates of 794 Brown Swiss bulls to assess ejaculate volume, sperm concentration, sperm motility, sperm head and tail anomalies and insemination success. The heritability of the six semen traits was between 0 and 0.26. Genome-wide association testing on 607,511 SNPs revealed a QTL on bovine chromosome 6 that was associated with sperm motility (P = 2.5 x 10⁻²⁷), head (P = 2.0 x 10⁻⁴⁴) and tail anomalies (P = 7.2 x 10⁻⁴⁹) and insemination success (P = 9.9 x 10⁻¹³). The QTL harbors a recessive allele that compromises semen quality and male fertility. We replicated the effect of the QTL on fertility (P = 7.1 x 10⁻³²) in an independent cohort of 2481 Brown Swiss bulls. The analysis of whole-genome sequencing data revealed that a synonymous variant (BTA6:58373887C>T, rs474302732) in WDR19 encoding WD repeat-containing protein 19 was in linkage disequilibrium with the fertility-associated haplotype. WD repeat-containing protein 19 is a constituent of the intraflagellar transport complex that is essential for the physiological function of motile cilia and flagella. Bioinformatic and transcription analyses revealed that the BTA6:58373887 T-allele activates a cryptic exonic splice site that eliminates three evolutionarily conserved amino acids from WDR19. Western blot analysis demonstrated that the BTA6:58373887 T-allele decreases protein expression. We make the remarkable observation that, in spite of negative effects on semen quality and bull fertility, the BTA6:58373887 T-allele has a frequency of 24% in the Brown Swiss population. Our findings are the first to uncover a variant that is associated with quantitative variation in semen quality and male fertility in cattle.

In cattle farming, artificial insemination is the most common method of breeding. To ensure high fertilization rates, ejaculate quality and insemination success are closely monitored in artificial insemination bulls. We analyse semen quality, insemination success and microarray-called genotypes at more than 600,000 genome-wide SNP markers of 794 bulls to identify a recessive allele that compromises semen quality. We take advantage of whole-genome sequencing to pinpoint a variant in the coding sequence of WDR19 encoding WD repeat-containing protein 19 that activates a novel exonic splice site. Our results indicate that cryptic splicing in WDR19 is associated with reduced male reproductive performance. This is the first report of a variant that contributes to quantitative variation in bovine semen quality.

Integration of whole-genome DNA methylation data with RNA sequencing data to identify markers for bull fertility

Predicting bull fertility prior to breeding is a current challenge for the dairy industry. The use of molecular biomarkers has been previously assessed. However, the integration of this information has not been performed to extract biologically relevant markers. The goal of this study was to integrate DNA methylation data with previously published RNA-sequencing results in order to identify candidate markers for sire fertility. A total of 1765 differentially methylated cytosines were found between high- and low-fertility sires. Ten genes associated with 11 differentially methylated cytosines were found in a previous study of gene expression between high- and low-fertility sires. Additionally, two of these genes code for proteins found exclusively in bull seminal plasma. Collectively, our results reveal 10 genes that could be used in the future as a panel for predicting bull fertility.

Deep learning versus parametric and ensemble methods for genomic prediction of complex phenotypes

Background

Transforming large amounts of genomic data into valuable knowledge for predicting complex traits has been an important challenge for animal and plant breeders. Prediction of complex traits has not escaped the current excitement on machine-learning, including interest in deep learning algorithms such as multilayer perceptrons (MLP) and convolutional neural networks (CNN). The aim of this study was to compare the predictive performance of two deep learning methods (MLP and CNN), two ensemble learning methods [random forests (RF) and gradient boosting (GB)], and two parametric methods [genomic best linear unbiased prediction (GBLUP) and Bayes B] using real and simulated datasets.

Methods

The real dataset consisted of 11,790 Holstein bulls with sire conception rate (SCR) records and genotyped for 58k single nucleotide polymorphisms (SNPs). To support the evaluation of deep learning methods, various simulation studies were conducted using the observed genotype data as template, assuming a heritability of 0.30 with either additive or non-additive gene effects, and two different numbers of quantitative trait nucleotides (100 and 1000).

Results

In the bull dataset, the best predictive correlation was obtained with GB (0.36), followed by Bayes B (0.34), GBLUP (0.33), RF (0.32), CNN (0.29) and MLP (0.26). The same trend was observed when using mean squared error of prediction. The simulation indicated that when gene action was purely additive, parametric methods outperformed other methods. When the gene action was a combination of additive, dominance and of two-locus epistasis, the best predictive ability was obtained with gradient boosting, and the superiority of deep learning over the parametric methods depended on the number of loci controlling the trait and on sample size. In fact, with a large dataset including 80k individuals, the predictive performance of deep learning methods was similar or slightly better than that of parametric methods for traits with non-additive gene action.

Conclusions

For prediction of traits with non-additive gene action, gradient boosting was a robust method. Deep learning approaches were not better for genomic prediction unless non-additive variance was sizable.

Gene mapping and genomic prediction of bull fertility using sex chromosome markers

Service sire has been recognized as an important factor affecting dairy herd fertility. Our group has reported promising results on gene mapping and genomic prediction of dairy bull fertility using autosomal SNP markers. Little is known, however, about the genetic contribution of sex chromosomes, which are enriched in genes related to sexual development and reproduction. As such, the main goal of this study was to investigate the effect of SNP markers on X and Y chromosomes (BTAX and BTAY, respectively) on sire conception rate (SCR) in US Holstein bulls. The analysis included a total of 5,014 bulls with SCR records and genotypes for roughly 291k SNP located on the autosomes, 1.5k SNP located on the pseudoautosomal region (PAR), 13.7k BTAX-specific SNP, and 24 BTAY-specific SNP. We first performed genomic scans of the sex chromosomes, and then we evaluated the genomic prediction of SCR including BTAX SNP markers in the predictive models. Two markers located on PAR and 3 markers located on the X-specific region showed significant associations with sire fertility. Interestingly, these regions harbor genes, such as FAM9B, TBL1X, and PIH1D3, that are directly implicated in testosterone concentration, spermatogenesis, and sperm motility. On the other hand, BTAY showed very low genetic variability, and none of the segregating markers were associated with SCR. Notably, model predictive ability was largely improved by including BTAX markers. Indeed, the combination of autosomal with BTAX SNP delivered predictive correlations around 0.343, representing an increase in accuracy of about 7.5% compared with the standard whole autosomal genome approach. Overall, this study provides evidence of the importance of both PAR and X-specific regions in male fertility in dairy cattle. These findings may help to improve conception rates in dairy herds through accurate genome-guided decisions on bull fertility.

Enhancing Genetic Gain through Genomic Selection: From Livestock to Plants

Although long-term genetic gain has been achieved through increasing use of modern breeding methods and technologies, the rate of genetic gain needs to be accelerated to meet humanity's demand for agricultural products. In this regard, genomic selection (GS) has been considered most promising for genetic improvement of the complex traits controlled by many genes each with minor effects. Livestock scientists pioneered GS application largely due to livestock's significantly higher individual values and the greater reduction in generation interval that can be achieved in GS. Large-scale application of GS in plants can be achieved by refining field management to improve heritability estimation and prediction accuracy and developing optimum GS models with the consideration of genotype-by-environment interaction and non-additive effects, along with significant cost reduction. Moreover, it would be more effective to integrate GS with other breeding tools and platforms for accelerating the breeding process and thereby further enhancing genetic gain. In addition, establishing an open-source breeding network and developing transdisciplinary approaches would be essential in enhancing breeding efficiency for small- and medium-sized enterprises and agricultural research systems in developing countries. New strategies centered on GS for enhancing genetic gain need to be developed.

Divergence Analyses of Sperm DNA Methylomes between Monozygotic Twin AI Bulls

Semen quality is critical for fertility. However, it is easily influenced by environmental factors and can induce subfertility in the next generations. Here, we aimed to assess the impacts of differentially methylated regions and genes on semen quality and offspring fertility. A specific pair of monozygotic (MZ) twin artificial insemination (AI) Holstein bulls with moderately different sperm qualities (Bull1 > Bull2) was used in the study, and each twin bull had produced ~6000 recorded daughters nationwide in China. Using whole genome bisulfite sequencing, we profiled the landscape of the twin bulls' sperm methylomes, and we observed markedly higher sperm methylation levels in Bull1 than in Bull2. Furthermore, we found 528 differentially methylated regions (DMR) between the MZ twin bulls, which spanned or overlapped with 309 differentially methylated genes (DMG). These DMG were particularly associated with embryo development, organ development, reproduction, and the nervous system. Several DMG were also shown to be differentially expressed in the sperm cells. Moreover, the significant differences in DNA methylation on gene INSL3 between the MZ twin bulls were confirmed at three different age points. Our results provided new insights into the impacts of AI bull sperm methylomes on offspring fertility.

In silico genome-wide miRNA-QTL-SNPs analyses identify a functional SNP associated with mastitis in Holsteins

Background

Single-nucleotide polymorphisms (SNPs) in microRNAs (miRNAs) and their target binding sites affect miRNA function and are involved in biological processes and diseases, including bovine mastitis, a frequent inflammatory disease. Our previous study has shown that bta-miR-2899 is significantly upregulated in the mammary gland tissue of mastitis-infected cow than that of healthy cows.

Results

In the present study, we used a customized miRNAQTLsnp software and identified 5252 SNPs in 691 bovine pre-miRNAs, which are also located within the quantitative trait loci (QTLs) that are associated with mastitis and udder conformation-related traits. Using luciferase assay in the bovine mammary epithelial cells, we confirmed a candidate SNP (rs109462250, g. 42,198,087 G > A) in the seed region of bta-miR-2899 located in the somatic cell score (SCS)-related QTL (Chr.18: 33.9–43.9 Mbp), which affected the interaction of bta-miR-2899 and its putative target Spi-1 proto-oncogene (SPI1), a pivotal regulator in the innate and adaptive immune systems. Quantitative real-time polymerase chain reaction results showed that the relative expression of SPI1 in the mammary gland of AA genotype cows was significantly higher than that of GG genotype cows. The SNP genotypes were associated with SCS in Holstein cows.

Conclusions

Altogether, miRNA-related SNPs, which influence the susceptibility to mastitis, are one of the plausible mechanisms underlying mastitis via modulating the interaction of miRNAs and immune-related genes. These miRNA-QTL-SNPs, such as the SNP (rs109462250) of bta-miR-2899 may have implication for the mastitis resistance breeding program in Holstein cattle.

Electronic supplementary material

The online version of this article (10.1186/s12863-019-0749-5) contains supplementary material, which is available to authorized users.

Improving Genomic Selection With Quantitative Trait Loci and Nonadditive Effects Revealed by Empirical Evidence in Maize

Genomic selection (GS), a tool developed for molecular breeding, is used by plant breeders to improve breeding efficacy by shortening the breeding cycle and to facilitate the selection of candidate lines for creating hybrids without phenotyping in various environments. Association and linkage mapping have been widely used to explore and detect candidate genes in order to understand the genetic mechanisms of quantitative traits. In the current study, phenotypic and genotypic data from three experimental populations, including data on six agronomic traits (e.g., plant height, ear height, ear length, ear diameter, grain yield per plant, and hundred-kernel weight), were used to evaluate the effect of trait-relevant markers (TRMs) on prediction accuracy estimation. Integrating information from mapping into a statistical model can efficiently improve prediction performance compared with using stochastically selected markers to perform GS. The prediction accuracy can reach plateau when a total of 500-1,000 TRMs are utilized in GS. The prediction accuracy can be significantly enhanced by including nonadditive effects and TRMs in the GS model when genotypic data with high proportions of heterozygous alleles and complex agronomic traits with high proportion of nonadditive variancein phenotypic variance are used to perform GS. In addition, taking information on population structure into account can slightly improve prediction performance when the genetic relationship between the training and testing sets is influenced by population stratification due to different allele frequencies. In conclusion, GS is a useful approach for prescreening candidate lines, and the empirical evidence provided by the current study for TRMs and nonadditive effects can inform plant breeding and in turn contribute to the improvement of selection efficiency in practical GS-assisted breeding programs.