Genetic dissection of bull fertility in US Jersey dairy cattle

Genomic insights into key genes and QTLs involved in cattle reproduction

From an economic standpoint, reproductive characteristics are fundamental for sustainable production, particularly for monotocous livestock like cattle. A longer inter-calving interval is indicative of low reproductive capacity. This issue changes the dynamics of current and future lactations since it necessitates more inseminations, veterinary care, and hormone interventions. Various reproductive phenotypes, including ovulation, mating, fertility, pregnancy, embryonic growth, and calving-related traits, are observed in dairy cattle, and these traits have been associated with several QTLs. Calving ease, age at puberty, scrotal circumference, and inseminations per conception have been associated with 4437, 10623, 10498, and 2476 Quantitative Trait Loci (QTLs), respectively. This data offers valuable insights into enhancing and comprehending reproductive traits in livestock breeding. Studying QTLs associated with reproductive traits has far-reaching implications across various fields, from agriculture and animal husbandry to human health, evolutionary biology, and conservation. It provides the foundation for informed breeding practices, advances in biotechnology, and a deeper understanding of the genetic underpinnings of reproduction.

Genome-wide association study considering genotype-by-environment interaction for productive and reproductive traits using whole-genome sequencing in Nellore cattle

BACKGROUND

The genotype-by-environment interaction (GxE) in beef cattle can be investigated using reaction norm models to assess environmental sensitivity and, combined with genome-wide association studies (GWAS), to map genomic regions related to animal adaptation. Including genetic markers from whole-genome sequencing in reaction norm (RN) models allows us to identify high-resolution candidate genes across environmental gradients through GWAS. Hence, we performed a GWAS via the RN approach using whole-genome sequencing data, focusing on mapping candidate genes associated with the expression of reproductive and growth traits in Nellore cattle. For this purpose, we used phenotypic data for age at first calving (AFC), scrotal circumference (SC), post-weaning weight gain (PWG), and yearling weight (YW). A total of 20,000 males and 7,159 females genotyped with 770k were imputed to the whole sequence (29 M). After quality control and linkage disequilibrium (LD) pruning, there remained ∼ 2.41 M SNPs for SC, PWG, and YW and ∼ 5.06 M SNPs for AFC.

RESULTS

Significant SNPs were identified on Bos taurus autosomes (BTA) 10, 11, 14, 18, 19, 20, 21, 24, 25 and 27 for AFC and on BTA 4, 5 and 8 for SC. For growth traits, significant SNP markers were identified on BTA 3, 5 and 20 for YW and PWG. A total of 56 positional candidate genes were identified for AFC, 9 for SC, 3 for PWG, and 24 for YW. The significant SNPs detected for the reaction norm coefficients in Nellore cattle were found to be associated with growth, adaptative, and reproductive traits. These candidate genes are involved in biological mechanisms related to lipid metabolism, immune response, mitogen-activated protein kinase (MAPK) signaling pathway, and energy and phosphate metabolism.

CONCLUSIONS

GWAS results highlighted differences in the physiological processes linked to lipid metabolism, immune response, MAPK signaling pathway, and energy and phosphate metabolism, providing insights into how different environmental conditions interact with specific genes affecting animal adaptation, productivity, and reproductive performance. The shared genomic regions between the intercept and slope are directly implicated in the regulation of growth and reproductive traits in Nellore cattle raised under different environmental conditions.

CFAP58 is involved in the sperm head shaping and flagellogenesis of cattle and mice

CFAP58 is a testis-enriched gene that plays an important role in the sperm flagellogenesis of humans and mice. However, the effect of CFAP58 on bull semen quality and the underlying molecular mechanisms involved in spermatogenesis remain unknown. Here, we identified two single-nucleotide polymorphisms (rs110610797, A>G and rs133760846, G>T) and one indel (g.-1811_ g.-1810 ins147bp) in the promoter of CFAP58 that were significantly associated with semen quality of bulls, including sperm deformity rate and ejaculate volume. Moreover, by generating gene knockout mice, we found for the first time that the loss of Cfap58 not only causes severe defects in the sperm tail, but also affects the manchette structure, resulting in abnormal sperm head shaping. Cfap58 deficiency causes an increase in spermatozoa apoptosis. Further experiments confirmed that CFAP58 interacts with IFT88 and CCDC42. Moreover, it may be a transported cargo protein that plays a role in stabilizing other cargo proteins, such as CCDC42, in the intra-manchette transport/intra-flagellar transport pathway. Collectively, our findings reveal that CFAP58 is required for spermatogenesis and provide genetic markers for evaluating semen quality in cattle.

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.

Molecular quantitative trait loci in reproductive tissues impact male fertility in cattle

Breeding bulls are well suited to investigate inherited variation in male fertility because they are genotyped and their reproductive success is monitored through semen analyses and thousands of artificial inseminations. However, functional data from relevant tissues are lacking in cattle, which prevents fine-mapping fertility-associated genomic regions. Here, we characterize gene expression and splicing variation in testis, epididymis, and vas deferens transcriptomes of 118 mature bulls and conduct association tests between 414,667 molecular phenotypes and 21,501,032 genome-wide variants to identify 41,156 regulatory loci. We show broad consensus in tissue-specific and tissue-enriched gene expression between the three bovine tissues and their human and murine counterparts. Expression- and splicing-mediating variants are more than three times as frequent in testis than epididymis and vas deferens, highlighting the transcriptional complexity of testis. Finally, we identify genes (WDR19, SPATA16, KCTD19, ZDHHC1) and molecular phenotypes that are associated with quantitative variation in male fertility through transcriptome-wide association and colocalization analyses.

A Single-Step Genome-Wide Association Study for Semen Traits of Egyptian Buffalo Bulls

The present study aimed to contribute to the limited research on buffalo (Bubalus bubalis) semen traits by incorporating genomic data. A total of 8465 ejaculates were collected. The genotyping procedure was conducted using the Axiom® Buffalo Genotyping 90 K array designed by the Affymetrix Expert Design Program. After conducting a quality assessment, we utilized 67,282 SNPs genotyped in 192 animals. We identified several genomic loci explaining high genetic variance by employing single-step genomic evaluation. The aforementioned regions were located on buffalo chromosomes no. 3, 4, 6, 7, 14, 16, 20, 22, and the X-chromosome. The X-chromosome exhibited substantial influence, accounting for 4.18, 4.59, 5.16, 5.19, and 4.31% of the genomic variance for ejaculate volume, mass motility, livability, abnormality, and concentration, respectively. In the examined genomic regions, we identified five novel candidate genes linked to male fertility and spermatogenesis, four in the X-chromosome and one in chromosome no. 16. Additional extensive research with larger sample sizes and datasets is imperative to validate these findings and evaluate their applicability for genomic selection.

Genome-wide scan for runs of homozygosity in South American Camelids

BACKGROUND

Alpaca (Vicugna pacos), llama (Lama glama), vicugna (Vicugna vicugna) and guanaco (Lama guanicoe), are the camelid species distributed over the Andean high-altitude grasslands, the Altiplano, and the Patagonian arid steppes. Despite the wide interest on these animals, most of the loci under selection are still unknown. Using whole-genome sequencing (WGS) data we investigated the occurrence and the distribution of Runs Of Homozygosity (ROHs) across the South American Camelids (SACs) genome to identify the genetic relationship between the four species and the potential signatures of selection.

RESULTS

A total of 37 WGS samples covering the four species was included in the final analysis. The multi-dimensional scaling approach showed a clear separation between the four species; however, admixture analysis suggested a strong genetic introgression from vicugna and llama to alpaca. Conversely, very low genetic admixture of the guanaco with the other SACs was found. The four species did not show significant differences in the number, length of ROHs (100-500 kb) and genomic inbreeding values. Longer ROHs (> 500 kb) were found almost exclusively in alpaca. Seven overlapping ROHs were shared by alpacas, encompassing nine loci (FGF5, LOC107034918, PRDM8, ANTXR2, LOC102534792, BSN, LOC116284892, DAG1 and RIC8B) while nine overlapping ROHs were found in llama with twenty-five loci annotated (ERC2, FZD9, BAZ1B, BCL7B, LOC116284208, TBL2, MLXIPL, PHF20, TRNAD-AUC, LOC116284365, RBM39, ARFGEF2, DCAF5, EXD2, HSPB11, LRRC42, LDLRAD1, TMEM59, LOC107033213, TCEANC2, LOC102545169, LOC116278408, SMIM15, NDUFAF2 and RCOR1). Four overlapping ROHs, with three annotated loci (DLG1, KAT6B and PDE4D) and three overlapping ROHs, with seven annotated genes (ATP6V1E1, BCL2L13, LOC116276952, BID, KAT6B, LOC116282667 and LOC107034552), were detected for vicugna and guanaco, respectively.

CONCLUSIONS

The signatures of selection revealed genomic areas potentially selected for production traits as well as for natural adaptation to harsh environment. Alpaca and llama hint a selection driven by environment as well as by farming purpose while vicugna and guanaco showed selection signals for adaptation to harsh environment. Interesting, signatures of selection on KAT6B gene were identified for both vicugna and guanaco, suggesting a positive effect on wild populations fitness. Such information may be of interest to further ecological and animal production studies.

Review: Genetic mutations affecting bull fertility

Cattle are a well-suited "model organism" to study the genetic underpinnings of variation in male reproductive performance. The adoption of artificial insemination and genomic prediction in many cattle breeds provide access to microarray-derived genotypes and repeated measurements for semen quality and insemination success in several thousand bulls. Similar-sized mapping cohorts with phenotypes for male fertility are not available for most other species precluding powerful association testing. The repeated measurements of the artificial insemination bulls' semen quality enable the differentiation between transient and biologically relevant trait fluctuations, and thus, are an ideal source of phenotypes for variance components estimation and genome-wide association testing. Genome-wide case-control association testing involving bulls with either aberrant sperm quality or low insemination success revealed several causal recessive loss-of-function alleles underpinning monogenic reproductive disorders. These variants are routinely monitored with customised genotyping arrays in the male selection candidates to avoid the use of subfertile or infertile bulls for artificial insemination and natural service. Genome-wide association studies with quantitative measurements of semen quality and insemination success revealed quantitative trait loci for male fertility, but the underlying causal variants remain largely unknown. Moreover, these loci explain only a small part of the heritability of male fertility. Integrating genome-wide association studies with gene expression and other omics data from male reproductive tissues is required for the fine-mapping of candidate causal variants underlying variation in male reproductive performance in cattle.

A potent allele marker related to low bull conception rate in Japanese Black bulls

Over the years, there has been considerable variation in the bull conception rate (BCR) of Japanese Black cattle; moreover, several Japanese Black bulls with a low BCR of ≤10% have been identified. However, the alleles responsible for the low BCR are not determined yet. Therefore, in this study, we aimed to identify single-nucleotide polymorphisms (SNPs) for predicting low BCR. To this end, the genome of Japanese Black bulls was comprehensively examined by a genome-wide association study with whole-exome sequencing (WES), and the effect of the identified marker regions on BCR was determined. The WES analysis of six sub-fertile bulls with a BCR of ≤10% and 73 normal bulls with a BCR of ≥40% identified a homozygous genotype for low BCR in Bos taurus autosome 5 in the region between 116.2 and 117.9 Mb. The g.116408653G > A SNP in this region had the most significant effect on the BCR (P-value = 1.0 × 10^-23), and the GG (55.4 ± 11.2%) and AG (54.4 ± 9.4%) genotypes in the SNP had a higher phenotype than the AA (9.5 ± 6.1%) genotype for the BCR. The mixed model analysis revealed that g.116408653G > A was related to approximately 43% of the total genetic variance. In conclusion, the AA genotype of g.116408653G > A is a useful index for identifying sub-fertile Japanese Black bulls. Some positive and negative effects of SNP on the BCR were presumed to identify the causative mutations, which can help evaluate bull fertility.

Haplotype-based association study of production and reproduction traits in multigenerational Vrindavani population

Haplotype-based association analysis promises to reveal important information regarding the effect of genetic variants on economic traits of interest. The present study aimed to evaluate the haplotype structure of Vrindavani cattle and explore the association of haplotypes with (re)production traits of economic interest. Genotyping array data of medium density (Bovine50KSNP BeadChip) on 96 randomly selected Vrindavani cows was used in the present study. Genotypes were called in GenomeStudio program while quality control was undertaken in PLINK using standard thresholds. The phenotypic traits used in the present study included age at first calving, dry days, lactation length, peak yield, total lactation milk yield, inter-calving period and service period. The haplotype structure of Vrindavani population was assessed, using a sliding window of 20 SNP with a shift of 5 SNPs at a time, in terms of the size of haplotype blocks regarding their length (in Kb) and frequency in chromosome-wise fashion. Haplotype blocks were assessed for possible association with important production and reproduction traits across three lactation cycles in Vrindavani cattle population. The first ten principal components were included in the model for haplotype-based association analysis to correct for stratification effects of assessed individuals. Multiple haplotypes were found to be associated with age at first calving, total lactation milk yield, peak yield, dry days, inter-calving period and service period. Various candidate genes were found to overlap haplotypes that were significantly associated with age at first calving (CDH18, MARCHF11, MYO10, FBXL7), total lactation milk yield (TGF, PDE1A, and COL8A1), peak yield (PPARGC1A, RCAN1, KCNE1, SMIM34 and MRPS6), dry days (CPNE4, ACAD11 and MRAS), inter-calving period (ABCG5, ABCG8 and COX7A2L) and service period (FOXL2 and PIK3CB). The putative candidate genes overlapping the significantly associated haplotypes revealed important pathways affecting the production and reproduction performance of animals. The identified genes and pathways may serve as good candidate markers to select animals for improved production and reproduction performance in future generations.

Identification of genomic regions associated with total and progressive sperm motility in Italian Holstein bulls

Sperm motility is directly related to the ability of sperm to move through the female reproductive tract to reach the ovum. Sperm motility is a complex trait that is influenced by environmental and genetic factors and is associated with male fertility, oocyte penetration rate, and reproductive success of cattle. In this study we carried out a GWAS in Italian Holstein bulls to identify candidate regions and genes associated with variations in progressive and total motility (PM and TM, respectively). After quality control, the final data set consisted of 5,960 records from 949 bulls having semen collected in 10 artificial insemination stations and genotyped at 412,737 SNPs (call rate >95%; minor allele frequency >5%). (Co)variance components were estimated using single trait mixed models, and associations between SNPs and phenotypes were assessed using a genomic BLUP approach. Ten windows that explained the greatest percentage of genetic variance were located on Bos taurus autosomes 1, 2, 4, 6, 7, 23, and 26 for TM and Bos taurus autosomes 1, 2, 4, 6, 8, 16, 23, and 26 for PM. A total of 150 genes for TM and 72 genes for PM were identified within these genomic regions. Gene Ontology enrichment analyses identified significant Gene Ontology terms involved with energy homeostasis, membrane functions, sperm-egg interactions, protection against oxidative stress, olfactory receptors, and immune system. There was significant enrichment of quantitative trait loci for fertility, calving ease, immune response, feed intake, and carcass weight within the candidate windows. These results contribute to understanding the architecture of the genetic control of sperm motility and may aid in the development of strategies to identify subfertile bulls and improve reproductive success.

Favored single nucleotide variants identified using whole genome Re-sequencing of Austrian and Chinese cattle breeds

Cattle have been essential for the development of human civilization since their first domestication few thousand years ago. Since then, they have spread across vast geographic areas following human activities. Throughout generations, the cattle genome has been shaped with detectable signals induced by various evolutionary processes, such as natural and human selection processes and demographic events. Identifying such signals, called selection signatures, is one of the primary goals of population genetics. Previous studies used various selection signature methods and normalized the outputs score using specific windows, in kbp or based on the number of SNPs, to identify the candidate regions. The recent method of iSAFE claimed for high accuracy in pinpointing the candidate SNPs. In this study, we analyzed whole-genome resequencing (WGS) data of ten individuals from Austrian Fleckvieh (Bos taurus) and fifty individuals from 14 Chinese indigenous breeds (Bos taurus, Bos taurus indicus, and admixed). Individual WGS reads were aligned to the cattle reference genome of ARS. UCD1.2 and subsequently undergone single nucleotide variants (SNVs) calling pipeline using GATK. Using these SNVs, we examined the population structure using principal component and admixture analysis. Then we refined selection signature candidates using the iSAFE program and compared it with the classical iHS approach. Additionally, we run Fst population differentiation from these two cattle groups. We found gradual changes of taurine in north China to admixed and indicine to the south. Based on the population structure and the number of individuals, we grouped samples to Fleckvieh, three Chinese taurines (Kazakh, Mongolian, Yanbian), admixed individuals (CHBI_Med), indicine individuals (CHBI_Low), and a combination of admixed and indicine (CHBI) for performing iSAFE and iHS tests. There were more significant SNVs identified using iSAFE than the iHS for the candidate of positive selection and more detectable signals in taurine than in indicine individuals. However, combining admixed and indicine individuals decreased the iSAFE signals. From both within-population tests, significant SNVs are linked to the olfactory receptors, production, reproduction, and temperament traits in taurine cattle, while heat and parasites tolerance in the admixed individuals. Fst test suggests similar patterns of population differentiation between Fleckvieh and three Chinese taurine breeds against CHBI. Nevertheless, there are genes shared only among the Chinese taurine, such as PAX5, affecting coat color, which might drive the differences between these yellowish coated breeds, and those in the greater Far East region.

Along the Bos taurus genome, uncover candidate imprinting control regions

Background

In mammals, Imprinting Control Regions (ICRs) regulate a subset of genes in a parent-of-origin-specific manner. In both human and mouse, previous studies identified a set of CpG-rich motifs occurring as clusters in ICRs and germline Differentially Methylated Regions (gDMRs). These motifs consist of the ZFP57 binding site (ZFBS) overlapping a subset of MLL binding units known as MLL morphemes. MLL or MLL1 (Mixed Lineage Leukemia 1) is a relatively large multidomain protein that plays a central role in the regulation of transcription. The structures of both MLL1 and MLL2 include a domain (MT) that binds CpG-rich DNA and a conserved domain (SET) that methylates lysine 4 in histone H3 producing H3K4me3 marks in chromatin.

Results

Since genomic imprinting impacts many developmental and key physiological processes, we followed a previous bioinformatics strategy to pinpoint ICR positions in the Bos taurus genome. Initial genome-wide analyses involved finding the positions of ZFP57 binding sites, and the CpG-rich motifs (ZFBS-morph overlaps) along cattle chromosomal DNA. By creating plots displaying the density of ZFBS-morph overlaps, we removed background noise and thus improved signal detection. With the density-plots, we could view the positions of peaks locating known and candidate ICRs in cattle DNA. Our evaluations revealed the correspondence of peaks in plots to reported known and inferred ICRs/DMRs in cattle. Beside peaks pinpointing such ICRs, the density-plots also revealed additional peaks. Since evaluations validated the robustness of our approach, we inferred that the additional peaks may correspond to candidate ICRs for imprinted gene expression.

Conclusion

Our bioinformatics strategy offers the first genome-wide approach for systematically localizing candidate ICRs. Furthermore, we have tailored our datasets for upload onto the UCSC genome browser so that researchers could find known and candidate ICRs with respect to a wide variety of annotations at all scales: from the positions of Single Nucleotide Polymorphisms (SNPs), to positions of genes, transcripts, and repeated DNA elements. Furthermore, the UCSC genome browser offers tools to produce enlarged views: to uncover the genes in the vicinity of candidate ICRs and thus discover potential imprinted genes for experimental validations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08694-3.

Deciphering the genetic basis of male fertility in Italian Brown Swiss dairy cattle

Improving reproductive performance remains a major goal in dairy cattle worldwide. Service sire has been recognized as an important factor affecting herd fertility. The main objective of this study was to reveal the genetic basis of male fertility in Italian Brown Swiss dairy cattle. Dataset included 1102 Italian Brown Swiss bulls with sire conception rate records genotyped with 454k single nucleotide polymorphisms. The analysis included whole-genome scans and gene-set analyses to identify genomic regions, individual genes and genetic mechanisms affecting Brown Swiss bull fertility. One genomic region on BTA1 showed significant additive effects. This region harbors gene RABL3 which is implicated cell proliferation and motility. Two genomic regions, located on BTA6 and BTA26, showed marked non-additive effects. These regions harbor genes, such as WDR19 and ADGRA1, that are directly involved in male fertility, including sperm motility, acrosome reaction, and embryonic development. The gene-set analysis revealed functional terms related to cell adhesion, cellular signaling, cellular transport, immune system, and embryonic development. Remarkably, a gene-set analysis also including Holstein and Jersey data, revealed significant processes that are common to the three dairy breeds, including cell migration, cell–cell interaction, GTPase activity, and the immune function. Overall, this comprehensive study contributes to a better understanding of the genetic basis of male fertility in cattle. In addition, our findings may guide the development of novel genomic strategies for improving service sire fertility in Brown Swiss cattle.

Genome-Wide Selection Signatures and Human-Mediated Introgression Events in Bos taurus indicus-influenced Composite Beef Cattle

Genetic introgression from interbreeding hybridization of European Bos taurus taurus (EBT) and Indian Bos taurus indicus (IBI) cattle breeds have been widely used to combine the climatic resilience of the IBI cattle and the higher productivity of EBT when forming new composite beef cattle (CB) populations. The subsequent breeding strategies have shifted their initial genomic compositions. To uncover population structure, signatures of selection, and potential introgression events in CB populations, high-density genotypes [containing 492,954 single nucleotide polymorphisms (SNPs) after the quality control] of 486 individuals from 15 cattle breeds, including EBT, IBI, and CB populations, along with two Bos grunniens genotypes as outgroup were used in this study. Then, in-depth population genetics analyses were performed for three CB breeds of Beefmaster, Brangus, and Santa Gertrudis. Neighbor-joining, principal components, and admixture analyses confirmed the historical introgression of EBT and IBI haplotypes into CB breeds. The f_dM statistics revealed that only 12.9% of CB populations' genetic components are of IBI origin. The results of signatures of selection analysis indicated different patterns of selection signals in the three CB breeds with primary pressure on pathways involved in protein processing and stress response in Beefmaster, cell proliferation regulation and immune response in Brangus, and amino acids and glucose metabolisms in Santa Gertrudis. An average of >90% of genomic regions underlying selection signatures were of EBT origin in the studied CB populations. Investigating the CB breeds' genome allows the estimation of EBT and IBI ancestral proportions and the locations within the genome where either taurine or indicine origin alleles are under selective pressure. Such findings highlight various opportunities to control the selection process more efficiently and explore complementarity at the genomic level in CB populations.

A Random Forest-Based Genome-Wide Scan Reveals Fertility-Related Candidate Genes and Potential Inter-Chromosomal Epistatic Regions Associated With Age at First Calving in Nellore Cattle

This study aimed to perform a genome-wide association analysis (GWAS) using the Random Forest (RF) approach for scanning candidate genes for age at first calving (AFC) in Nellore cattle. Additionally, potential epistatic effects were investigated using linear mixed models with pairwise interactions between all markers with high importance scores within the tree ensemble non-linear structure. Data from Nellore cattle were used, including records of animals born between 1984 and 2015 and raised in commercial herds located in different regions of Brazil. The estimated breeding values (EBV) were computed and used as the response variable in the genomic analyses. After quality control, the remaining number of animals and SNPs considered were 3,174 and 360,130, respectively. Five independent RF analyses were carried out, considering different initialization seeds. The importance score of each SNP was averaged across the independent RF analyses to rank the markers according to their predictive relevance. A total of 117 SNPs associated with AFC were identified, which spanned 10 autosomes (2, 3, 5, 10, 11, 17, 18, 21, 24, and 25). In total, 23 non-overlapping genomic regions embedded 262 candidate genes for AFC. Enrichment analysis and previous evidence in the literature revealed that many candidate genes annotated close to the lead SNPs have key roles in fertility, including embryo pre-implantation and development, embryonic viability, male germinal cell maturation, and pheromone recognition. Furthermore, some genomic regions previously associated with fertility and growth traits in Nellore cattle were also detected in the present study, reinforcing the effectiveness of RF for pre-screening candidate regions associated with complex traits. Complementary analyses revealed that many SNPs top-ranked in the RF-based GWAS did not present a strong marginal linear effect but are potentially involved in epistatic hotspots between genomic regions in different autosomes, remarkably in the BTAs 3, 5, 11, and 21. The reported results are expected to enhance the understanding of genetic mechanisms involved in the biological regulation of AFC in this cattle breed.

Candidate Genes in Bull Semen Production Traits: An Information Approach Review

Semen quality plays a crucial role in the successful implementation of breeding programs, especially where artificial insemination (AI) is practiced. Bulls with good semen traits have good fertility and can produce a volume of high semen per ejaculation. The aim of this review is to use an information approach to highlight candidate genes and their relation to bull semen production traits. The use of genome-wide association studies (GWAS) has been demonstrated to be successful in identifying genomic regions and individual variations associated with production traits. Studies have reported over 40 genes associated with semen traits using Illumina BeadChip single-nucleotide polymorphism (SNPs).

Bull fertility and semen quality are not correlated with dairy and production traits in Brown Swiss cattle

Undisturbed reproduction is key for successful breeding of beef and dairy cattle. Improving reproductive ability can be difficult because of antagonistic relationships with other economically relevant traits. In cattle, thorough investigation of female fertility revealed unfavorable genetic correlations with various production phenotypes. However, the correlation between male reproductive ability and production traits remains poorly understood. Here, we investigated the genetic relationships among and between male fertility characteristics and economically relevant traits in a population of Brown Swiss cattle. We performed GWAS with imputed genotypes at nearly 12 million sequence variants for semen quality (sperm head and tail anomalies, motility, concentration, and volume), male fertility, and 57 production phenotypes. Allele substitution effects were then correlated on a trait-by-trait basis to estimate genetic correlations. Correlations between male reproductive characteristics and traits of economic value were small and ranged from -0.0681 to 0.0787. Among the semen quality parameters, sperm motility was negatively correlated with anomalies (head: r = -0.7083 ± 0.0002; tail: r = -0.7739 ± 0.0002) and volume (r = -0.1266 ± 0.0003), whereas volume was negatively correlated with concentration (r = -0.3503 ± 0.0002). Sire nonreturn rate was negatively correlated with sperm anomalies (head: r = -0.1640 ± 0.0002; tail: r = -0.1580 ± 0.0002) and positively correlated with motility (r = 0.1598 ± 0.0002). A meta-analysis of male reproductive traits identified 2 quantitative trait loci: a previously described region on chromosome 6 showed pleiotropic effects and a novel region on chromosome 11 was associated with sperm head anomalies. In conclusion, our results suggest that selection for economically important dairy and production phenotypes has little impact on semen quality and fertility of Brown Swiss bulls.

Genome-wide association study of trypanosome prevalence and morphometric traits in purebred and crossbred Baoulé cattle of Burkina Faso

In this study, single-SNP GWAS analyses were conducted to find regions affecting tolerance against trypanosomosis and morphometrics traits in purebred and crossbred Baoulé cattle of Burkina Faso. The trypanosomosis status (positive and negative) and a wide set of morphological traits were recorded for purebred Baoulé and crossbred Zebu x Baoulé cattle, and genotyped with the Illumina Bovine SNP50 BeadChip. After quality control, 36,203 SNPs and 619 animals including 343 purebred Baoulé and 279 crossbreds were used for the GWAS analyses. Several important genes were found that can influence morphological parameters. Although there were no genes identified with a reported strong connection to size traits, many of them were previously identified in various growth-related studies. A re-occurring theme for the genes residing in the regions identified by the most significant SNPs was pleiotropic effect on growth of the body and the cardiovascular system. Regarding trypanosomosis tolerance, two potentially important regions were identified in purebred Baoulé on chromosomes 16 and 24, containing the CFH, CRBN, TRNT1 and, IL5RA genes, and one additional genomic region in Baoulé, x Zebu crossbreds on chromosome 5, containing MGAT4C and NTS. Almost all of these regions and genes were previously related to the trait of interest, while the CRBN gene was to our knowledge presented in the context of trypanosomiasis tolerance for the first time.

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.

Targeted sequencing reveals candidate causal variants for dairy bull subfertility

Bull fertility is a key factor for successful reproductive performance in dairy cattle. Since the semen from a single bull can be used to inseminate hundreds of cows, one subfertile bull could have a major impact on herd reproductive efficiency. We have previously identified five genomic regions, located on BTA8 (72.2 Mb), BTA9 (43.7 Mb), BTA13 (60.2 Mb), BTA17 (63.3 Mb), and BTA27 (34.7 Mb), that show large dominance effects on bull fertility. Each of these regions explained about 5–8% of the observed differences in sire conception rate between Holstein bulls. Here, we aimed to identify candidate causal variants responsible for this variation using targeted sequencing (10 Mb per region). For each genomic region, two DNA pools were constructed from n≈20 high‐fertility and n≈20 low‐fertility Holstein bulls. The DNA‐sequencing analysis included reads quality control (using FastQC), genome alignment (using BWA and ARS‐UCD1.2), variant calling (using GATK) and variant annotation (using Ensembl). The sequencing depth per pool varied from 39× to 51×. We identified a set of nonsense mutations, missense mutations, and frameshift variants carried by low‐fertility bulls. Notably, some of these variants were classified as strong candidate causal variants, i.e., mutations with deleterious effects located on genes exclusively/highly expressed in testis. Genes affected by these candidate causal variants include AK9, TTLL9, TCHP, and FOXN4. These results could aid in the development of novel genomic tools that allow early detection and culling of subfertile bull calves.

Autosomal recessive loci contribute significantly to quantitative variation of male fertility in a dairy cattle population

BACKGROUND

Cattle are ideally suited to investigate the genetics of male fertility. Semen from individual bulls is used for thousands of artificial inseminations for which the fertilization success is monitored. Results from the breeding soundness examination and repeated observations of semen quality complement the fertility evaluation for each bull.

RESULTS

In a cohort of 3881 Brown Swiss bulls that had genotypes at 683,609 SNPs, we reveal four novel recessive QTL for male fertility on BTA1, 18, 25, and 26 using haplotype-based association testing. A QTL for bull fertility on BTA1 is also associated with sperm head shape anomalies. All other QTL are not associated with any of the semen quality traits investigated. We perform complementary fine-mapping approaches using publicly available transcriptomes as well as whole-genome sequencing data of 125 Brown Swiss bulls to reveal candidate causal variants. We show that missense or nonsense variants in SPATA16, VWA3A, ENSBTAG00000006717 and ENSBTAG00000019919 are in linkage disequilibrium with the QTL. Using whole-genome sequence data, we detect strong association (P = 4.83 × 10^- 12) of a missense variant (p.Ile193Met) in SPATA16 with male fertility. However, non-coding variants exhibit stronger association at all QTL suggesting that variants in regulatory regions contribute to variation in bull fertility.

CONCLUSION

Our findings in a dairy cattle population provide evidence that recessive variants may contribute substantially to quantitative variation in male fertility in mammals. Detecting causal variants that underpin variation in male fertility remains difficult because the most strongly associated variants reside in poorly annotated non-coding regions.

Regulatory Potential of Long Non-Coding RNAs (lncRNAs) in Boar Spermatozoa with Good and Poor Freezability

Long non-coding RNAs (lncRNAs) are suggested to play an important role in the sperm biological processes. We performed de novo transcriptome assembly to characterize lncRNAs in spermatozoa, and to investigate the role of the potential target genes of the differentially expressed lncRNAs (DElncRNAs) in sperm freezability. We detected approximately 4007 DElncRNAs, which were differentially expressed in spermatozoa from boars classified as having good and poor semen freezability (GSF and PSF, respectively). Most of the DElncRNAs were upregulated in boars of the PSF group and appeared to significantly affect the sperm's response to the cryopreservation conditions. Furthermore, we predicted that the potential target genes were regulated by DElncRNAs in cis or trans. It was found that DElncRNAs of both freezability groups had potential cis- and trans-regulatory effects on different protein-coding genes, such as COX7A2L, TXNDC8 and SOX-7. Gene Ontology (GO) enrichment revealed that the DElncRNA target genes are associated with numerous biological processes, including signal transduction, response to stress, cell death (apoptosis), motility and embryo development. Significant differences in the de novo assembled transcriptome expression profiles of the DElncRNAs between the freezability groups were confirmed by quantitative real-time PCR analysis. This study reveals the potential effects of protein-coding genes of DElncRNAs on sperm functions, which could contribute to further research on their relevance in semen freezability.

Elucidating the processes and pathways enriched in buffalo sperm proteome in regulating semen quality

Sperm carries a reservoir of proteins regulating the molecular functions to attain functional competence. Semen samples collected from buffalo bulls were assessed for sperm functional attributes (n = 11) and proteome profiling (n = 6). Sperm proteins were extracted and profiled by employing LC-MS/MS. Overall, the buffalo sperm contained 1365 proteins, of which 458 were common between the groups. The unique proteins were 477 and 430 in good and poor quality semen, respectively. In the whole proteome of buffalo sperm, sexual reproduction with phosphatidylethanolamine-binding protein1 (PEBP1), fetuin-B (FETUB) and acrosin (ACR) was the most enriched (p = 8.44E-19) biological process, also with thermogenesis (p = 0.003), oocyte meiosis (p = 0.007) and vascular smooth muscle contraction (p = 0.009) apart from metabolic pathways. In good quality semen, mesenchyme migration (p = 1.24E-07) and morphogenesis (p = 0.001) were abundant biological processes. In good quality semen, the fluid shear stress (p = 0.01) and, in poor quality semen, valine, leucine and isoleucine degradation (p = 3.8E-05) pathways were enriched. In good quality semen, 7 proteins were significantly (p < 0.05) upregulated and 33 proteins were significantly (p < 0.05) downregulated. On validating the abundantly expressed sperm proteins, serine protease inhibitor Kazal-type 2-like (SPINK2; 2.17-fold) and neddylin (NEDD8; 1.13-fold) were upregulated and YBX2 was downregulated (0.41-fold) in good quality semen as compared with poor quality semen (1-fold). The present findings revealed the importance of sperm proteins in oocyte maturation, fertilization process and early embryonic development. The variations in the proteomic composition can be used as potential markers for the selection of breeding bulls.

Genome-Wide Association Study and Pathway Analysis for Female Fertility Traits in Iranian Holstein Cattle

Female fertility is an important trait that contributes to cow’s profitability and it can be improved by genomic information. The objective of this study was to detect genomic regions and variants affecting fertility traits in Iranian Holstein cattle. A data set comprised of female fertility records and 3,452,730 pedigree information from Iranian Holstein cattle were used to predict the breeding values, which were then employed to estimate the de-regressed proofs (DRP) of genotyped animals. A total of 878 animals with DRP records and 54k SNP markers were utilized in the genome-wide association study (GWAS). The GWAS was performed using a linear regression model with SNP genotype as a linear covariate. The results showed that an SNP on BTA19, ARS-BFGL-NGS-33473, was the most significant SNP associated with days from calving to first service. In total, [69] significant SNPs were located within 27 candidate genes. Novel potential candidate genes include OSTN, DPP6, EphA5, CADPS2, Rfc1, ADGRB3, Myo3a, C10H14orf93, KIAA1217, RBPJL, SLC18A2, GARNL3, NCALD, ASPH, ASIC2, OR3A1, CHRNB4, CACNA2D2, DLGAP1, GRIN2A and ME3. These genes are involved in different pathways relevant to female fertility and other characteristics in mammals. Gene set enrichment analysis showed that thirteen GO terms had significant overrepresentation of genes statistically associated with female fertility traits. The results of network analysis identified CCNB1 gene as a hub gene in the progesterone-mediated oocyte maturation pathway, significantly associated with age at first calving. The candidate genes identified in this study can be utilized in genomic tests to improve reproductive performance in Holstein cattle.

Post-Translational Protein Deimination Signatures in Serum and Serum-Extracellular Vesicles of Bos taurus Reveal Immune, Anti-Pathogenic, Anti-Viral, Metabolic and Cancer-Related Pathways for Deimination

The bovine immune system is known for its unusual traits relating to immunoglobulin and antiviral responses. Peptidylarginine deiminases (PADs) are phylogenetically conserved enzymes that cause post-translational deimination, contributing to protein moonlighting in health and disease. PADs also regulate extracellular vesicle (EV) release, forming a critical part of cellular communication. As PAD-mediated mechanisms in bovine immunology and physiology remain to be investigated, this study profiled deimination signatures in serum and serum-EVs in Bos taurus. Bos EVs were poly-dispersed in a 70-500 nm size range and showed differences in deiminated protein cargo, compared with whole sera. Key immune, metabolic and gene regulatory proteins were identified to be post-translationally deiminated with some overlapping hits in sera and EVs (e.g., immunoglobulins), while some were unique to either serum or serum-EVs (e.g., histones). Protein-protein interaction network analysis of deiminated proteins revealed KEGG pathways common for serum and serum-EVs, including complement and coagulation cascades, viral infection (enveloped viruses), viral myocarditis, bacterial and parasitic infections, autoimmune disease, immunodeficiency intestinal IgA production, B-cell receptor signalling, natural killer cell mediated cytotoxicity, platelet activation and hematopoiesis, alongside metabolic pathways including ferroptosis, vitamin digestion and absorption, cholesterol metabolism and mineral absorption. KEGG pathways specific to EVs related to HIF-1 signalling, oestrogen signalling and biosynthesis of amino acids. KEGG pathways specific for serum only, related to Epstein-Barr virus infection, transcription mis-regulation in cancer, bladder cancer, Rap1 signalling pathway, calcium signalling pathway and ECM-receptor interaction. This indicates differences in physiological and pathological pathways for deiminated proteins in serum-EVs, compared with serum. Our findings may shed light on pathways underlying a number of pathological and anti-pathogenic (viral, bacterial, parasitic) pathways, with putative translatable value to human pathologies, zoonotic diseases and development of therapies for infections, including anti-viral therapies.

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.

Structural Insights into the Intracellular Region of the Human Magnesium Transport Mediator CNNM4

The four member family of "Cyclin and Cystathionine β-synthase (CBS) domain divalent metal cation transport mediators", CNNMs, are the least-studied mammalian magnesium transport mediators. CNNM4 is abundant in the brain and the intestinal tract, and its abnormal activity causes Jalili Syndrome. Recent findings show that suppression of CNNM4 in mice promotes malignant progression of intestinal polyps and is linked to infertility. The association of CNNM4 with phosphatases of the regenerating liver, PRLs, abrogates its Mg2+-efflux capacity, thus resulting in an increased intracellular Mg2+ concentration that favors tumor growth. Here we present the crystal structures of the two independent intracellular domains of human CNNM4, i.e., the Bateman module and the cyclic nucleotide binding-like domain (cNMP). We also derive a model structure for the full intracellular region in the absence and presence of MgATP and the oncogenic interacting partner, PRL-1. We find that only the Bateman module interacts with ATP and Mg2+, at non-overlapping sites facilitating their positive cooperativity. Furthermore, both domains dimerize autonomously, where the cNMP domain dimer forms a rigid cleft to restrict the Mg2+ induced sliding of the inserting CBS1 motives of the Bateman module, from a twisted to a flat disk shaped dimer.

Selection for bull fertility: a review

Fertility is a critically important factor in cattle production because it directly relates to the ability to produce the offspring necessary to offset costs in production systems. Female fertility has received much attention and has been enhanced through assisted reproductive technologies, as well as genetic selection; however, improving bull fertility has been largely ignored. Improvements in bull reproductive performance are necessary to optimize the efficiency of cattle production. Selection and management to improve bull fertility not only have the potential to increase conception rates but also have the capacity to improve other economically relevant production traits. Bull fertility has reportedly been genetically correlated with traits such as average daily gain, heifer pregnancy, and calving interval. Published studies show that bull fertility traits are low to moderately heritable, indicating that improvements in bull fertility can be realized through selection. Although female fertility has continued to progress according to increasing conception rates, the reported correlation between male and female fertility is low, indicating that male fertility cannot be improved by selection for female fertility. Correlations between several bull fertility traits, such as concentration, number of spermatozoa, motility, and number of spermatozoa abnormalities, vary among studies. Using male fertility traits in selection indices would provide producers with more advanced selection tools. The objective of this review was to discuss current beef bull fertility measurements and to discuss the future of genetic evaluation of beef bull fertility and potential genetic improvement strategies.

Across-breed validation study confirms and identifies new loci associated with sexual precocity in Brahman and Nellore cattle

The aim of this study was to identify candidate regions associated with sexual precocity in Bos indicus. Nellore and Brahman were set as validation and discovery populations, respectively. SNP selected in Brahman to validate in Nellore were from gene regions affecting reproductive traits (G1) and significant SNP (p ≤ 10^-3 ) from a meta-analysis (G2). In the validation population, early pregnancy (EP) and scrotal circumference (SC) were evaluated. To perform GWAS in validation population, we used regression and Bayes C. SNP with p ≤ 10^-3 in regression and Bayes factor ≥3 in Bayes C were deemed significant. Significant SNP (for EP or SC) or SNP in their ±250 Kb vicinity region, which were in at least one discovery set (G1 or G2), were considered validated. SNP identified in both G1 and G2 were considered candidate. For EP, 145 SNP were validated in G1 and 41 in G2, and for SC, these numbers were 14 and 2. For EP, 21 candidate SNP were detected (G1 and G2). For SC, no candidate SNP were identified. Validated SNP and their vicinity region were located close to quantitative trait loci or genes related to reproductive traits and were enriched in gene ontology terms related to reproductive success. These are therefore strong candidate regions for sexual precocity in Nellore and Brahman.

Genomic prediction of bull fertility in US Jersey dairy cattle

Service sire has a major effect on reproductive success in dairy cattle. Recent studies have reported accurate predictions for Holstein bull fertility using genomic data. The objective of this study was to assess the feasibility of genomic prediction of sire conception rate (SCR) in US Jersey cattle using alternative predictive models. Data set consisted of 1.5k Jersey bulls with SCR records and 95k SNP covering the entire genome. The analyses included the use of linear and Gaussian kernel-based models fitting either all the SNP or subsets of markers with presumed functional roles, such as SNP significantly associated with SCR or SNP located within or close to annotated genes. Model predictive ability was evaluated using 5-fold cross-validation with 10 replicates. The entire SNP set exhibited predictive correlations around 0.30. Interestingly, either SNP marginally associated with SCR or genic SNP achieved higher predictive abilities than their counterparts using random sets of SNP. Among alternative SNP subsets, Gaussian kernel models fitting significant SNP achieved the best performance with increases in predictive correlation up to 7% compared with the standard whole-genome approach. Notably, the use of a multi-breed reference population including the entire US Holstein SCR data set (11.5k bulls) allowed us to achieve predictive correlations up to 0.315, gaining 8% in accuracy compared with the standard model fitting a pure Jersey reference set. Overall, our findings indicate that genomic prediction of Jersey bull fertility is feasible. The use of Gaussian kernels fitting markers with relevant roles and the inclusion of Holstein records in the training set seem to be promising alternatives to the standard whole-genome approach. These results have the potential to help the dairy industry improve US Jersey sire fertility through accurate genome-guided decisions.