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.

Within- and between-Breed Selection Signatures in the Original and Improved Valachian Sheep

This study explored the genomic diversity and selection signatures in two Slovakian national breeds, the Original Valachian and the Improved Valachian sheep. As they are an important animal genetic resource within the country, but with decreasing population size, our aim is to identify potentially valuable genomic regions. A total of 97 sheep (18 male and 79 female) from the Original Valachian, and 69 sheep (25 male and 44 female) from the Improved Valachian populations were genotyped using the GeneSeek GGP Ovine 50 K chip. The inbreeding levels were assessed with runs of homozygosity (ROH). The selection signatures within breeds were identified based on the top 1% of most homozygous regions within the breed, the so-called ROH islands. The selection signatures between breeds were assessed based on variance in linkage disequilibrium. Overall, we have identified selection signatures with quantitative trait loci (QTL) and genes pointing towards all three production purposes of the Valachian sheep, milk, meat, and wool, including their quality characteristics. Another group with apparent large importance was the various traits related to health and resistance to parasites, which is well in line with the sturdy nature of this breed.

Characteristics of miRNAs Present in Bovine Sperm and Associations With Differences in Fertility

Small non-coding RNAs have been linked to different phenotypes in bovine sperm, however attempts to identify sperm-borne molecular biomarkers of male fertility have thus far failed to identify a robust profile of expressed miRNAs related to fertility. We hypothesized that some differences in bull fertility may be reflected in the levels of different miRNAs in sperm. To explore such differences in fertility that are not due to differences in visible metrics of sperm quality, we employed Next Generation Sequencing to compare the miRNA populations in Bos taurus sperm from bulls with comparable motility and morphology but varying Sire Conception Rates. We identified the most abundant miRNAs in both populations (miRs -34b-3p; -100-5p; -191-5p; -30d-4p; -21-5p) and evaluated differences in the overall levels and specific patterns of isomiR expression. We also explored correlations between specific pairs of miRNAs in each population and identified 10 distinct pairs of miRNAs that were positively correlated in bulls with higher fertility and negatively correlated in comparatively less fertile individuals. Furthermore, 8 additional miRNA pairs demonstrated the opposite trend; negatively correlated in high fertility animals and positively correlated in less fertile bulls. Finally, we performed pathway analysis to identify potential roles of miRNAs present in bull sperm in the regulation of specific genes that impact spermatogenesis and embryo development. Together, these results present a comprehensive picture of the bovine sperm miRNAome that suggests multiple potential roles in fertility.

Taurine and Indicine Haplotype Representation in Advanced Generation Individuals From Three American Breeds

Development of the American Breeds of beef cattle began in the 1920s as breeders and U. S. Experiment Station researchers began to create Bos taurus taurus × Bos taurus indicus hybrids using Brahman as the B. t. indicus source. By 1954, U.S. Breed Associations had been formed for Brangus (5/8 Angus × 3/8 Brahman), Beefmaster (½ Brahman × ¼ Shorthorn × ¼ Hereford), and Santa Gertrudis (5/8 Shorthorn × 3/8 Brahman). While these breeds were developed using mating designs expected to create base generation animals with the required genome contributions from progenitor breeds, each association has now registered advanced generation animals in which selection or drift may have caused the realized genome compositions to differ from initial expected proportions. The availability of high-density SNP genotypes for 9,161 Brangus, 3,762 Beefmaster, and 1,942 Santa Gertrudis animals allowed us to compare the realized genomic architectures of breed members to the base generation expectations. We used RFMix to estimate local ancestry and identify genomic regions in which the proportion of Brahman ancestry differed significantly from a priori expectations. For all three breeds, lower than expected levels of Brahman composition were found genome-wide, particularly in early-generation animals where we demonstrate that selection on beef production traits was likely responsible for the taurine enrichment. Using a proxy for generation number, we also contrasted the genomes of early- and advanced-generation animals and found that the indicine composition of the genome has increased with generation number likely due to selection on adaptive traits. Many of the most-highly differentiated genomic regions were breed specific, suggesting that differences in breeding objectives and selection intensities exist between the breeds. Global ancestry estimation is commonly performed in admixed animals to control for stratification in association studies. However, local ancestry estimation provides the opportunity to investigate the evolution of specific chromosomal segments and estimate haplotype effects on trait variation in admixed individuals. Investigating the genomic architecture of the American Breeds not only allows the estimation of indicine and taurine genome proportions genome-wide, but also the locations within the genome where either taurine or indicine alleles confer a selective advantage.

Genetic and Genomic Analyses of Service Sire Effect on Female Reproductive Traits in Holstein Cattle

Fertility and reproductive performance are key drivers of dairy farm profitability. Hence, reproduction traits have been included in a large majority of worldwide dairy cattle selection indexes. The reproductive traits are lowly heritable but can be improved through direct genetic selection. However, most scientific studies and dairy cattle breeding programs have focused solely on the genetic effects of the dam (GED) on reproductive performance and, therefore, ignored the contribution of the service sire in the phenotypic outcomes. This study aimed to investigate the service sire effects on female reproductive traits in Holstein cattle from a genomic perspective. Genetic parameter estimation and genome-wide association studies (GWAS) were performed for the genetic effect of service sire (GESS) on conception rate (CR), 56-day non-return rate (NRR56), calving ease (CE), stillbirth (SB), and gestation length (GL). Our findings indicate that the additive genetic effects of both sire and dam contribute to the phenotypic variance of reproductive traits measured in females (0.0196 vs. 0.0109, 0.0237 vs. 0.0133, 0.0040 vs. 0.0289, 0.0782 vs. 0.0083, and 0.1024 vs. 0.1020 for GESS and GED heritability estimates for CR, NRR56, CE, SB, and GL, respectively), and these two genetic effects are positively correlated for SB (0.1394) and GL (0.7871). Interestingly, the breeding values for GESS on insemination success traits (CR and NRR56) are unfavorably and significantly correlated with some production, health, and type breeding values (ranging from -0.449 to 0.274), while the GESS values on calving traits (CE, SB, and GL) are usually favorably associated with those traits (ranging from -0.493 to 0.313). One hundred sixty-two significant single-nucleotide polymorphisms (SNPs) and their surrounding protein-coding genes were identified as significantly associated with GESS and GED, respectively. Six genes overlapped between GESS and GED for calving traits and 10 genes overlapped between GESS for success traits and calving traits. Our findings indicate the importance of considering the GESS when genetically evaluating the female reproductive traits in Holstein cattle.

Genome-wide association and genomic prediction for scrotal circumference in Hereford and Braford bulls

Scrotal circumference (SC) is widely used as a selection criterion for bulls in breeding programs, since it is easily assessed and correlated with several desirable reproductive traits. The objectives of this study were: to perform a genome-wide association study (GWAS) to identify genomic regions associated with SC adjusted for age (SCa) and for both age and weight (SCaw); to select Tag SNPs from GWAS to construct low-density panel for genomic prediction; and to compare the prediction accuracy of the SC through different methods for Braford and Hereford bulls from the same genetic breeding program. Data of SC from 18,172 bulls (30.4 ± 3.7 cm) and of genotypes from 131 sires and 3,545 animals were used. From GWAS, the top 1% of 1-Mb windows were observed on chromosome (BTA) 2, 20, 7, 8, 15, 3, 16, 27, 6 and 8 for SCa and on BTA 8, 15, 16, 21, 19, 2, 6, 5 and 10 for SCaw, representing 17.4% and 18.8% of the additive genetic variance of SCa and SCaw, respectively. The MeSH analysis was able to translate genomic information providing biological meanings of more specific gene functions related to the SCa and SCaw. The genomic enhancement methods, especially single step GBLUP, that combined phenotype and pedigree data with direct genomic values generated gains in accuracy in relation to pedigree BLUP, suggesting that genomic predictions should be applied to improve genetic gain and to narrow the generation interval compared to traditional methods. The proposed Tag-SNP panels may be useful for lower-cost commercial genomic prediction applications in the future, when the number of bulls in the reference population increases for SC in Hereford and Braford breeds.

Development of a RAPD marker-based classification criterion for quality semen production in Holstein crossbred bulls

In cattle production systems, an intense selection pressure for production traits has resulted in the decline of fertility traits. To optimize an efficient reproduction system, the inclusion of both male and female fertility traits in the selection process is very much essential. RAPD (Random Amplified Polymorphic DNA) was developed as a molecular biology tool and has been extensively used, to study intra- and interspecific genetic diversity. The present study was undertaken to utilize RAPD primers to investigate the association between DNA markers and semen quality traits viz. Sperm concentration, total sperm count ejaculate and initial sperm motility and thereby to identify good/poor semen producers. DNA isolated from the blood samples of healthy bulls was subjected to RAPD-PCR. The multiple regression analysis followed by independent t test was carried out to identify suitable markers. Based on the results, only 12 bands were identified as marker suitable for any of the quality trait. This includes, OPA2 ~ 760, OPA2 ~ 700, OPA6 ~ 1,200, OPA9 ~ 400, OPA9 ~ 380, OPA12 ~ 970, OPA14 ~ 715, OPA14 ~ 605, OPA16 ~ 485, OPA17 ~ 860 and OPA18 ~ 480. Multiple regression analysis selected, OPA2 ~ 760 and OPA2 ~ 1,750 for sperm concentration and OPA2 ~ 760, OPA2 ~ 700, OPA9 ~ 620, OPA4 ~ 670 and OPA18 ~ 1,015 for total sperm count/ejaculate. But the t test revealed a significant association between OPA2 ~ 760 and total sperm count. Further, discriminant function analysis also identified this marker in the first step itself. The results of the present study can be exploited as a low-cost alternative strategy for identification of good /poor semen producers in crossbred bulls at an early age.

Genetic polymorphism in MAP1B gene associated with conception rate in Holstein Friesian crossbred breeding bulls

The microtubule associated protein 1 B (MAP1B) gene is one of the candidate genes which was found to be strongly associated with the sire conception rate. In the present study, identification of single nucleotide polymorphisms (SNPs) in the targeted region of MAP1B gene and their effect on conception rate (CR) in Karan Fries (KF) breeding bulls were sought. A tetra-primer ARMS PCR technique was adapted to genotype an SNP, g. 3,066A>G located on intron 1 of MAP1B gene. The relationship between the CR of KF breeding bulls and polymorphic locus of MAP1B gene was assessed. Genotyping of the targeted region revealed three kinds of genotype (AA, AG and GG) with respect to g. 3,066A>G SNP locus. The identified SNP had highly significant (P=0.0009) influence on CR. The KF bulls with AA and AG genotype had comparatively higher CR than the bulls with GG genotype. The identified SNP was significantly associated with CR, thus it may be utilized as a molecular marker for developing marker-assisted selection strategy to enhance the CR in KF cattle population.

Association of α/β-Hydrolase D16B with Bovine Conception Rate and Sperm Plasma Membrane Lipid Composition

We have identified a Holstein sire named Tarantino who had been approved for artificial insemination that is based on normal semen characteristics (i.e., morphology, thermoresistance, motility, sperm concentration), but had no progeny after 412 first inseminations, resulting in a non-return rate (NRdev) of -29. Using whole genome association analysis and next generation sequencing, an associated nonsense variant in the α/β-hydrolase domain-containing 16B gene (ABHD16B) on bovine chromosome 13 was identified. The frequency of the mutant allele in the German Holstein population was determined to be 0.0018 in 222,645 investigated cattle specimens. The mutant allele was traced back to Whirlhill Kingpin (bornFeb. 13th, 1959) as potential founder. The expression of ABHD16B was detected by Western blotting and immunohistochemistry in testis and epididymis of control bulls. A lipidome comparison of the plasma membrane of fresh semen from carriers and controls showed significant differences in the concentration of phosphatidylcholine (PC), diacylglycerol (DAG), ceramide (Cer), sphingomyelin (SM), and phosphatidylcholine (-ether) (PC O-), indicating that ABHD16B plays a role in lipid biosynthesis. The altered lipid contents may explain the reduced fertilization ability of mutated sperms.

Whole-genome scan reveals significant non-additive effects for sire conception rate in Holstein cattle

Background

Service sire has a considerable impact on reproductive success in dairy cattle. Most gene mapping studies for bull fertility have focused on additive effects, while non-additive effects have been largely ignored. The main goal of this study was to assess the relevance of non-additive effects on Sire Conception Rate (SCR) in Holstein dairy cattle. The analysis included 7.5 k Holstein bulls with both SCR records and 57.8 k single nucleotide polymorphism (SNP) markers spanning the entire genome.

Results

The importance of non-additive effects was evaluated using an efficient two-step mixed model-based approach. Four genomic regions located on chromosomes BTA8, BTA9, BTA13 and BTA17 showed marked dominance and/or recessive effects. Most of these regions harbor genes, such as ADAM28, DNAJA1, TBC1D20, SPO11, PIWIL3 and TMEM119, that are directly implicated in testis development, male germ line maintenance, and sperm maturation.

Conclusions

This study provides further evidence for the relevance of non-additive effects in fitness-related traits, such as male fertility. In addition, these findings may point out new strategies for improving service sire fertility in dairy cattle via marker-assisted selection.

A dual targeted β-defensin and exome sequencing approach to identify, validate and functionally characterise genes associated with bull fertility

Bovine fertility remains a critical issue underpinning the sustainability of the agricultural sector. Phenotypic records collected on >7,000 bulls used in artificial insemination (AI) were used to identify 160 reliable and divergently fertile bulls for a dual strategy of targeted sequencing (TS) of fertility-related β-defensin genes and whole exome sequencing (WES). A haplotype spanning multiple β-defensin genes and containing 94 SNPs was significantly associated with fertility and functional analysis confirmed that sperm from bulls possessing the haplotype showed significantly enhanced binding to oviductal epithelium. WES of all exons in the genome in 24 bulls of high and low fertility identified 484 additional SNPs significantly associated with fertility. After validation, the most significantly associated SNP was located in the FOXJ3 gene, a transcription factor which regulates sperm function in mice. This study represents the first comprehensive characterisation of genetic variation in bovine β-defensin genes and functional analysis supports a role for β-defensins in regulating bull sperm function. This first application of WES in AI bulls with divergent fertility phenotypes has identified a novel role for the transcription factor FOXJ3 in the regulation of bull fertility. Validated genetic variants associated with bull fertility could prove useful for improving reproductive outcomes in cattle.

Identification of single nucleotide polymorphism in protein phosphatase 1 regulatory subunit 11 gene in Murrah bulls

Aim:

This study was conducted with the objective to identify single nucleotide polymorphism (SNP) in protein phosphatase 1 regulatory subunit 11 (PPP1R11) gene in Murrah bulls.

Materials and Methods:

Genomic DNA was isolated by phenol–chloroform extraction method from the frozen semen samples of 65 Murrah bulls maintained at Artificial Breeding Research Centre, ICAR-National Dairy Research Institute, Karnal. The quality and concentration of DNA was checked by spectrophotometer reading and agarose gel electrophoresis. The target region of PPP1R11 gene was amplified using four sets of primer designed based on Bos taurus reference sequence. The amplified products were sequenced and aligned using Clustal Omega for identification of SNPs. Animals were genotyped by polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) using EcoNI restriction enzyme.

Results:

The sequences in the NCBI accession number NW_005785016.1 for Bubalus bubalis were compared and aligned with the edited sequences of Murrah bulls with Clustal Omega software. A total of 10 SNPs were found, out of which 1 at 5’UTR, 3 at intron 1, and 6 at intron 2 region. PCR-RFLP using restriction enzyme EcoNI revealed only AA genotype indicating monomorphism in PPP1R11 gene of all Murrah animals included in the study.

Conclusion:

A total of 10 SNPs were found. PCR-RFLP revealed only AA genotype indicating monomorphism in PPP1R11 gene of all Murrah animals included in the study, due to which association analysis with conception rate was not feasible.

Unravelling the genomic architecture of bull fertility in Holstein cattle

BACKGROUND

Fertility is considered an important economic trait in dairy cattle. Most studies have investigated cow fertility while bull fertility has received much less consideration. The main objective of this study was to perform a comprehensive genomic analysis in order to unravel the genomic architecture underlying sire fertility in Holstein dairy cattle. The analysis included the application of alternative genome-wide association mapping approaches and the subsequent use of diverse gene set enrichment tools.

RESULTS

The association analyses identified at least eight genomic regions strongly associated with bull fertility. Most of these regions harbor genes, such as KAT8, CKB, TDRD9 and IGF1R, with functions related to sperm biology, including sperm development, motility and sperm-egg interaction. Moreover, the gene set analyses revealed many significant functional terms, including fertilization, sperm motility, calcium channel regulation, and SNARE proteins. Most of these terms are directly implicated in sperm physiology and male fertility.

CONCLUSIONS

This study contributes to the identification of genetic variants and biological processes underlying sire fertility. These findings can provide opportunities for improving bull fertility via marker-assisted selection.

Inferring quantitative trait pathways associated with bull fertility from a genome-wide association study

Whole-genome association studies typically focus on genetic markers with the strongest evidence of association. However, single markers often explain only a small component of the genetic variance and hence offer a limited understanding of the trait under study. As such, the objective of this study was to perform a pathway-based association analysis in Holstein dairy cattle in order to identify relevant pathways involved in bull fertility. The results of a single-marker association analysis, using 1,755 bulls with sire conception rate data and genotypes for 38,650 single nucleotide polymorphisms (SNPs), were used in this study. A total of 16,819 annotated genes, including 2,767 significantly associated with bull fertility, were used to interrogate a total of 662 Gene Ontology (GO) terms and 248 InterPro (IP) entries using a test of proportions based on the cumulative hypergeometric distribution. After multiple-testing correction, 20 GO categories and one IP entry showed significant overrepresentation of genes statistically associated with bull fertility. Several of these functional categories such as small GTPases mediated signal transduction, neurogenesis, calcium ion binding, and cytoskeleton are known to be involved in biological processes closely related to male fertility. These results could provide insight into the genetic architecture of this complex trait in dairy cattle. In addition, this study shows that quantitative trait pathways inferred from single-marker analyses could enhance our interpretations of the results of genome-wide association studies.