Genome-wide association study for semen volume and total number of sperm in Holstein-Friesian bulls

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.

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.

Effects of Heterologous Bovine Seminal Plasma-Supplemented Egg Yolk-Based Extender on Cryosurvivability of Pantja Buck Semen

The present study was conducted to observe the effects of removal of seminal plasma of Pantja buck semen and supplementation of bovine seminal plasma (BSP) in the extender before cryopreservation. In a preliminary experiment, different levels of BSP were supplemented (1, 3, 5, 7, and 9% v/v) in egg yolk (7.5% egg yolk)-tris (EYT) extender and used for cryopreservation of Pantja buck semen. Results in terms of motility, viability, plasma membrane integrity, acrosome integrity, and lipid peroxidation showed that 5% BSP was suitable for maintaining Pantja buck semen quality during cryopreservation. In the final experiment, pooled semen from four Pantja bucks was split into three aliquots (I, II, and III). Aliquot I was directly diluted in EYT extender and grouped as the control (C); aliquot II and III were washed separately with TALP solution and diluted as D1 (Washed semen with EYT extender) and D2 (Washed semen with EYT extender containing 5% BSP), respectively. Seminal attributes (sperm individual motility, viability, plasma membrane integrity, acrosome integrity, and total morphological abnormalities) were assessed at the postdilution, postequilibration, and post-thawing stages. Alanine aminotransferase (ALT), aspartate aminotransferase (AST), malondialdehyde (MDA) concentration, and glutathione peroxidase (GSH-Px) activity were measured at post-thaw. Washed semen significantly improved (p < 0.05) seminal parameters at post-thaw compared with unwashed semen (control). A significant difference (p < 0.05) was observed in seminal attributes between freezing stages and between dilution groups. Significantly higher (p < 0.05) post-thaw sperm motility, viability, plasma membrane integrity, acrosome integrity, and GSH-Px activity, and significantly lower (p < 0.05) MDA concentration and extracellular release of enzymes (ALT, AST) were observed in group D2 compared with control and D1. The results of the present study demonstrated that cryopreservation of washed Pantja buck semen diluted with 5% BSP-supplemented EYT extender can improve post-thaw semen quality.

QTLs and Candidate Genes Associated with Semen Traits in Merino Sheep

Ram semen traits play a significant role in conception outcomes, which in turn may influence reproductive efficiency and the overall productivity and profitability of sheep enterprises. Since hundreds of ewes may be inseminated from a single ejaculate, it is important to evaluate semen quality prior to use in sheep breeding programs. Given that semen traits have been found to be heritable, genetic variation likely contributes to the variability observed in these traits. Identifying such genetic variants could provide novel insights into the molecular mechanisms underlying variability in semen traits. Therefore, this study aimed to identify quantitative trait loci (QTLs) associated with semen traits in Merino sheep. A genome-wide association study (GWAS) was undertaken using 4506 semen collection records from 246 Merino rams collected between January 2002 and May 2021. The R package RepeatABEL was used to perform a GWAS for semen volume, gross motility, concentration, and percent post-thaw motility. A total of 35 QTLs, located on 16 Ovis aries autosomes (OARs), were significantly associated with either of the four semen traits in this study. A total of 89, 95, 33, and 73 candidate genes were identified, via modified Bonferroni, within the QTLs significantly associated with volume, gross motility, concentration, and percent post-thaw motility, respectively. Among the candidate genes identified, SORD, SH2B1, and NT5E have been previously described to significantly influence spermatogenesis, spermatozoal motility, and high percent post-thaw motility, respectively. Several candidate genes identified could potentially influence ram semen traits based on existing evidence in the literature. As such, validation of these putative candidates may offer the potential to develop future strategies to improve sheep reproductive efficiency. Furthermore, Merino ram semen traits are lowly heritable (0.071-0.139), and thus may be improved by selective breeding.

Candidate genes and proteins regulating bull semen quality: a review

Poor semen profile reflected by suboptimum fertility statistics is a concern in bulls reared for breeding purpose. A critical review of research on candidate genes and proteins associated with semen quality traits will be useful to understand the progress of molecular marker development for bull semen quality traits. Here, we have tabulated and classified candidate genes and proteins associated with bull semen quality based on a literature survey. A total of 175 candidate genes are associated with semen quality traits in various breeds of cattle. Several studies using candidate gene approach have identified 26 genes carrying a total of 44 single nucleotide polymorphisms. Furthermore, nine genome-wide association studies (GWASes) have identified 150 candidate genes using bovine single nucleotide polymorphisms (SNP) chips. Three genes, namely membrane-associated ring-CH-type finger 1 (MARCH1), platelet-derived growth factor receptor beta, and phosphodiesterase type 1, were identified commonly in two GWASes, which, especially MARCH1, are required to explore their regulatory roles in bull semen quality in in-depth studies. With the advancement of high-throughput-omic technologies, more candidate genes associated with bull semen quality may be identified in the future. Therefore, the functional significance of candidate genes and proteins need to be delved further into future investigations to augment bull semen quality.

Runs of Homozygosity and Quantitative Trait Locus/Association for Semen Parameters in Selected Chinese and South African Beef Cattle

In this study, runs of homozygosity (ROH) and quantitative trait locus/association (QTL) for semen parameters in selected Chinese and South African beef cattle breed were estimated. The computed results showed 7516 ROH were observed between classes 0−5 Mb with no ROH observed in classes >40 Mb. Distribution of ROH showed high level of genomic coverage for ANG, NGU, CSI, and BEL breeds. Approximately 13 genomic regions with QTL were controlling sperm motility, sperm concentration, semen volume, sperm count, sperm head abnormalities, sperm tail abnormalities, sperm integrity, and percentage of abnormal sperm traits. Nine candidate genes, CDF9, MARCH1, WDR19, SLOICI, ST7, DOP1B, CFAF9, INHBA, and ADAMTS1, were suggested to be associated with above mentioned QTL traits. The results for inbreeding coefficient showed moderate correlation between FROH vs FHOM at 0.603 and high correlation between FROH 0−5 Mb 0.929, and lowest correlation for 0−>40 Mb 0.400. This study suggested recent inbreeding in CSI, BEL, ANG, BON, SIM, and NGU breeds. Furthermore, it highlighted varied inbreeding levels and identified QTL for semen traits and genes of association. These results can assist in implementation of genetic improvement strategies for bulls and provide awareness and proper guidelines in developing breeding programs.

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).

Predicting sperm production of young dairy bulls using collection history and management factors

Selection of elite young dairy bulls by using genomic data shortened the generation interval and increased pressure to collect and market germplasm at an early age. The objectives of this study were (1) develop prediction models for daily, weekly, and monthly total sperm (TSp) production from collection history, health status, and management factors, and (2) assess the ability of these models to forecast future TSp production, as well as differences in prediction accuracy by seasonality or age of bull. Data consisted of 43,918 daily processing records from 1,037 Holstein and Jersey bulls between 10 and 28 mo of age at collection. Potential explanatory variables included year and season of collection, barn location, collection frequency, breed, scrotal circumference, TSp in previous months, health events, and age at arrival, first collection, and current collection. Linear regression, random forest (RF), Bayesian regularized neural network, model tree, multilayer perceptron neural network with multiple layers, and extreme learning machine were used to predict daily, weekly, and monthly TSp (R v3.5.1, https://www.r-project.org/). In the additive approach, all prior data were used for training; however, in the fixed-window approach, records from 3 previous months were used for age-based prediction, records from 4 previous months or 1 yr were used for the monthly date-based analyses, and records from 1 previous month or year were used for the weekly date-based analyses. Model performance was measured by root mean squared error (RMSE) and the correlation (r) between actual and predicted TSp in testing sets. In monthly analyses, RF with additive training performed best in age-based (RMSE = 13.6 billion cells, r = 0.93) and date-based (RMSE = 11.9, r = 0.94) prediction, compared with linear regression (age-based RMSE = 16.6, r = 0.89; date-based RMSE = 15.5, r = 0.90) and Bayesian regularized neural network (age-based RMSE = 14.1, r = 0.92). On average, RMSE was 0.93 or 0.14 billion cells greater with fixed 4-mo or 1-yr training windows, respectively, than in the additive analyses. The most important management variables affecting TSp were collection frequency, TSp in previous months, and age at collection. Results indicate RF models with additive training can predict TSp output of individual bulls with ≥85% accuracy up to 4 mo into the future. Spikes in accuracy were associated with sire summary times and company processing changes, and accuracy tended to stabilize when bulls reached 19 to 20 mo of age.

Genome-wide association studies for sperm traits in Assaf sheep breed

Sperm quality traits routinely collected by artificial insemination (AI) center for rams progeny test are related with the capacity to produce sperm doses for AI and, in more or less grade, with males' fertility. Low-quality ejaculates are unuseful to perform AI sperm doses, which suppose high economic loses for the AI center. Moreover, sperm quality traits have low heritability values which make traditional genetic selection little efficient to its improvement. In this work, a genome-wide association study (GWAS) was conducted by using sperm quality traits data and 50 K Affymetrix custom chip genotypes of 429 rams of Assaf breed from OVIGEN AI centre. Furthermore, 47 of these rams were also genotyped with the Illumina HD Ovine BeadChip, and therefore HD genotypes were imputed for all rams with phenotype data. Previous to the GWAS, a linear regression model was fitted including sperm traits as dependent variables; the flock of origin, date of sperm collection, and jump number as fixed effects; rams age at collection in months as covariate; and ram permanent effect as random. Pseudo-phenotypes obtained from this model were used as input for GWAS. Associations at the chromosome-wise level (FDR 10%) of 76 single-nucleotide polymorphisms (SNPs) in 4 chromosomes for ejaculate concentration (CON), 20 SNPs in 3 chromosomes for ejaculate volume (VOL), 32 SNPs in 1 chromosome for ejaculate number of spermatozoa (SPZ), and 23 SNPs for spermatozoa mass motility (MOT) in 17 chromosomes were found. Only SNPs associated with MOT overcame the genome-wide significance level. Some candidate genes for sperm traits variability were SLC9C1 (OAR1), TSN (OAR2), and FUT10 (OAR26) for MOT;. DOCK2, CPLANE1, SPEF2, and RAI14 (OAR16) for CON; SCAPER and PSMA4 (OAR18) for VOL; and PARM1 and LOC101110593 (OAR6) for SPZ. SNPs associated with sperm traits were not found to be correlated with milk production genetic variation; however, the high frequencies of some SNPs with negative effect over sperm traits found in animals at the top milk yield estimated breeding values (EBVs) ranking would allow to exert some selective presure to improve rams sperm performances. Effects and frequencies of some of the SNPs detected over sperm quality traits make these variants good candidates to be used in marker-assisted selection to improve sperm characteristics of Assaf rams and AI center efficiency to produce sperm doses.

A systems biology framework integrating GWAS and RNA-seq to shed light on the molecular basis of sperm quality in swine

BACKGROUND

Genetic pressure in animal breeding is sparking the interest of breeders for selecting elite boars with higher sperm quality to optimize ejaculate doses and fertility rates. However, the molecular basis of sperm quality is not yet fully understood. Our aim was to identify candidate genes, pathways and DNA variants associated to sperm quality in swine by analysing 25 sperm-related phenotypes and integrating genome-wide association studies (GWAS) and RNA-seq under a systems biology framework.

RESULTS

By GWAS, we identified 12 quantitative trait loci (QTL) associated to the percentage of head and neck abnormalities, abnormal acrosomes and motile spermatozoa. Candidate genes included CHD2, KATNAL2, SLC14A2 and ABCA1. By RNA-seq, we identified a wide repertoire of mRNAs (e.g. PRM1, OAZ3, DNAJB8, TPPP2 and TNP1) and miRNAs (e.g. ssc-miR-30d, ssc-miR-34c, ssc-miR-30c-5p, ssc-miR-191, members of the let-7 family and ssc-miR-425-5p) with functions related to sperm biology. We detected 6128 significant correlations (P-value ≤ 0.05) between sperm traits and mRNA abundances. By expression (e)GWAS, we identified three trans-expression QTL involving the genes IQCJ, ACTR2 and HARS. Using the GWAS and RNA-seq data, we built a gene interaction network. We considered that the genes and interactions that were present in both the GWAS and RNA-seq networks had a higher probability of being actually involved in sperm quality and used them to build a robust gene interaction network. In addition, in the final network we included genes with RNA abundances correlated with more than four semen traits and miRNAs interacting with the genes on the network. The final network was enriched for genes involved in gamete generation and development, meiotic cell cycle, DNA repair or embryo implantation. Finally, we designed a panel of 73 SNPs based on the GWAS, eGWAS and final network data, that explains between 5% (for sperm cell concentration) and 36% (for percentage of neck abnormalities) of the phenotypic variance of the sperm traits.

CONCLUSIONS

By applying a systems biology approach, we identified genes that potentially affect sperm quality and constructed a SNP panel that explains a substantial part of the phenotypic variance for semen quality in our study and that should be tested in other swine populations to evaluate its relevance for the pig breeding sector.

Bull seminal plasma stimulates in vitro production of TGF-β, IL-6 and IL-8 from bovine endometrial epithelial cells, depending on dose and bull fertility

Seminal plasma (SP) regulates immune responses in the female reproductive tract through specific cytokines. It is not known whether SP from high fertility bulls (H) differs from SP from low fertility bulls (L). In this study, the cytokine response of bovine endometrial epithelial cells (bEEC) in culture was investigated after challenge with SP from two bulls of below average (L) or three bulls of above average fertility (H). The bEECs were challenged with 1% or 4% SP from l- or H-fertility bulls (L1, L4, H1, H4, respectively) or 1% or 4% PBS as control (C1, C4) for 72 h. The culture media were analysed for concentrations (pg/million cells) of transforming growth factor beta (TGF-β1, TGF-β2 and TGF-β3) by Luminex, and Interleukin 6 and 8 (IL-6, IL-8) by ELISA. Challenge significantly affected production of TGF-ß1, TGF-ß2 and IL-8 compared to controls and was affected by bull fertility (p < 0.0001), SP concentration (p < 0.0001) and their interaction (p < 0.0001). A higher production of TGF-β1, TGF-β2 and IL-8 (p < 0.0001), and also IL-6 (p < 0.01), resulted from challenge with high doses of SP, being higher for L than H (p < 0.05). For TGF-β3, fertility of bull (p < 0.05). For TGF-B3, fertility of bull (p < 0.05) and the interaction between fertility and concentration of SP were significant (p < 0.01). In conclusion, 4% SP from L bulls stimulated more TGF-β1, TGF-β2, TGF-β3, IL-6 and IL-8 production than SP from H bulls, indicating that stimulation of the endometrium is relevant for fertility. Seminal plasma from high fertility bulls seems to affect cytokine production in utero positively in inseminated cows.

Genetic and genomic characterization of vulva size traits in Yorkshire and Landrace gilts

BACKGROUND

Reproductive performance is critical for efficient swine production. Recent results indicated that vulva size (VS) may be predictive of reproductive performance in sows. Study objectives were to estimate genetic parameters, identify genomic regions associated, and estimate genomic prediction accuracies (GPA) for VS traits.

RESULTS

Heritability estimates of VS traits, vulva area (VA), height (VH), and width (VW) measurements, were moderately to highly heritable in Yorkshire, with 0.46 ± 0.10, 0.55 ± 0.10, 0.31 ± 0.09, respectively, whereas these estimates were low to moderate in Landrace, with 0.16 ± 0.09, 0.24 ± 0.11, and 0.08 ± 0.06, respectively. Genetic correlations within VS traits were very high for both breeds, with the lowest of 0.67 ± 0.29 for VH and VW for Landrace. Genome-wide association studies (GWAS) for Landrace, reveled genomic region associated with VS traits on Sus scrofa chromosome (SSC) 2 (154-157 Mb), 7 (107-110 Mb), 8 (4-6 Mb), and 10 (8-19 Mb). For Yorkshire, genomic regions on SSC 1 (87-91 and 282-287 Mb) and 5 (67 Mb) were identified. All regions explained at least 3.4% of the genetic variance. Accuracies of genomic prediction were moderate in Landrace, ranging from 0.30 (VH) to 0.61 (VA), and lower for Yorkshire, with 0.07 (VW) to 0.11 (VH). Between-breed and multi-breed genomic prediction accuracies were low.

CONCLUSIONS

Our findings suggest that VS traits are heritable in Landrace and Yorkshire gilts. Genomic analyses show that major QTL control these traits, and they differ between breed. Genomic information can be used to increase genetic gains for these traits in gilts. Additional research must be done to validate the GWAS and genomic prediction results reported in our study.

RAPID COMMUNICATION: Multi-breed validation study unraveled genomic regions associated with puberty traits segregating across tropically adapted breeds1

An efficient strategy to improve QTL detection power is performing across-breed validation studies. Variants segregating across breeds are expected to be in high linkage disequilibrium (LD) with causal mutations affecting economically important traits. The aim of this study was to validate, in a Tropical Composite cattle (TC) population, QTL associations identified for sexual precocity traits in a Nellore and Brahman meta-analysis genome-wide association study. In total, 2,816 TC, 8,001 Nellore, and 2,210 Brahman animals were available for the analysis. For that, genomic regions significantly associated with puberty traits in the meta-analysis study were validated for the following sexual precocity traits in TC: age at first corpus luteum (AGECL), first postpartum anestrus interval (PPAI), and scrotal circumference at 18 months of age (SC). We considered validated QTL those underpinned by significant markers from the Nellore and Brahman meta-analysis (P ≤ 10-4) that were also significant for a TC trait, i.e., presenting a P-value of ≤10-3 for AGECL, PPAI, or SC. We also considered as validated QTL those regions where significant markers in the reference population were at ±250 kb from significant markers in the validation population. Using this criteria, 49 SNP were validated for AGECL, 4 for PPAI, and 14 for SC, from which 5 were in common with AGECL, totaling 62 validated SNP for these traits and 30 candidate genes surrounding them. Considering just candidate genes closest to the top SNP of each chromosome, for AGECL 8 candidate genes were identified: COL8A1, PENK, ENSBTAG00000047425, BPNT1, ADAMTS17, CCHCR1, SUFU, and ENSBTAG00000046374. For PPAI, 3 genes emerged as candidates (PCBP3, KCNK10, and MRPS5), and for SC 8 candidate genes were identified (SNORA70, TRAC, ASS1, BPNT1, LRRK1, PKHD1, PTPRM, and ENSBTAG00000045690). Several candidate regions presented here were previously associated with puberty traits in cattle. The majority of emerging candidate genes are related to biological processes involved in reproductive events, such as maintenance of gestation, and some are known to be expressed in reproductive tissues. Our results suggested that some QTL controlling early puberty seem to be segregating across cattle breeds adapted to tropical conditions.

Semen quality and diversity of morphological sperm abnormalities in bulls: breed and strain effects

At present great attention is paid to studying genetic regulation of farm animal adaptations to environmental conditions. This problem is very important due to a wide expansion of highly productive cattle breeds created in Europe and North America. However, until the present no investigation of changing semen quality in bulls of imported breeds during their adaptations to environmental conditions of Western Siberia has been conducted. The aim of this study was to investigate semen quality peculiarities and the diversity of morphological sperm abnormalities in bulls of imported and local breeds kept in the environmental conditions of the southern part of Western Siberia. We determined sperm concentration, sperm count, and rate of sperm with progressive motility and percentage of morphologically normal spermatozoa. The rate of sperm abnormalities according to Blome’s classifcation was determined too. It was found that the mean values of sperm concentration, sperm motility and percentage of morphologically normal spermatozoa in the bulls investigated were similar to those in bulls kept in European countries. Inter­breed differences in these parameters were not found. However, bulls of the Red Danish, Angler, and Simmental breeds had a higher percentage of misshapen sperm head and pyriform sperm head than bulls of the Black­White breed. An inter­strain difference in sperm motility in bulls of the Black­White breed was observed. It was found that bulls of Reflection Sovereign 198998 strain had lover sperm motility than bulls of Wis Burke Ideal 1013415 strain. No inter­strain differences in sperm production, percentage of morphologically normal spermatozoa and rate of main sperm abnormalities were found. Thus, it has been found that the environmental conditions of the southern part of Western Siberia do not seriously affect the sperm production, sperm motility or percentage of morphologically normal spermatozoa in bulls. However, the increased rate of misshapen and pyriform sperm heads in the bulls of the foreign breeds points to a need to study sperm DNA fragmentation.

Genetic mechanisms underlying spermatic and testicular traits within and among cattle breeds: systematic review and prioritization of GWAS results

Reduced bull fertility imposes economic losses in bovine herds. Specifically, testicular and spermatic traits are important indicators of reproductive efficiency. Several genome-wide association studies (GWAS) have identified genomic regions associated with these fertility traits. The aims of this study were as follows: 1) to perform a systematic review of GWAS results for spermatic and testicular traits in cattle and 2) to identify key functional candidate genes for these traits. The identification of functional candidate genes was performed using a systems biology approach, where genes shared between traits and studies were evaluated by a guilt by association gene prioritization (GUILDify and ToppGene software) in order to identify the best functional candidates. These candidate genes were integrated and analyzed in order to identify overlapping patterns among traits and breeds. Results showed that GWAS for testicular-related traits have been developed for beef breeds only, whereas the majority of GWAS for spermatic-related traits were conducted using dairy breeds. When comparing traits measured within the same study, the highest number of genes shared between different traits was observed, indicating a high impact of the population genetic structure and environmental effects. Several chromosomal regions were enriched for functional candidate genes associated with fertility traits. Moreover, multiple functional candidate genes were enriched for markers in a species-specific basis, taurine (Bos taurus) or indicine (Bos indicus). For the different candidate regions identified in the GWAS in the literature, functional candidate genes were detected as follows: B. Taurus chromosome X (BTX) (TEX11, IRAK, CDK16, ATP7A, ATRX, HDAC6, FMR1, L1CAM, MECP2, etc.), BTA17 (TRPV4 and DYNLL1), and BTA14 (MOS, FABP5, ZFPM2). These genes are responsible for regulating important metabolic pathways or biological processes associated with fertility, such as progression of spermatogenesis, control of ciliary activity, development of Sertoli cells, DNA integrity in spermatozoa, and homeostasis of testicular cells. This study represents the first systematic review on male fertility traits in cattle using a system biology approach to identify key candidate genes for these traits.

Effect of adding heterologous versus homologous bovine seminal plasma prior to cryopreservation on bull sperm quality after thawing

SummaryThe aim of this study was to investigate the effect of adding homologous or heterologous bovine seminal plasma (SP) to SP-free sperm samples before freezing on sperm quality after thawing. Ejaculates from bulls of known fertility were used as a source of SP. The SP was removed from further aliquots of the same ejaculates by colloid centrifugation to create SP-free sperm samples; the resuspended sperm pellets were treated with homologous or heterologous SP from high or low fertility bulls at 0%, 1% or 5% before freezing. After thawing, sperm quality was evaluated by computer-assisted sperm analysis and flow cytometry for membrane integrity, reactive oxygen species, chromatin structure, mitochondrial membrane potential and protein tyrosine phosphorylation. Data were analysed using Proc MIXED, SAS®. Post-hoc comparisons were adjusted for multiplicity using Tukey's method. The addition of SP resulted in significant differences in sperm quality, namely velocity class A, Velocity Straight Line (VSL), Velocity Average Path (VAP), Velocity Curved Line (VCL), Amplitude of Lateral Head Displacement (ALH), Hyperactive (HYP), reactive oxygen species (ROS) production and % DNA fragmentation index (DFI) (P<0.05 for each). Although adding 5% homologous SP from high fertility bulls was beneficial to sperm kinematics, 5% heterologous SP from high fertility bulls had a deleterious effect on chromatin integrity and on sperm velocity. In conclusion, adding SP may have either a beneficial effect or a deleterious effect depending on the individuals involved. It might be feasible to use this method to improve sperm quality in some circumstances.

Characterization of biological pathways associated with semen traits in the Thai multibreed dairy population

The objective of this research was to characterize biological pathways associated with semen volume (VOL), number of sperm (NS), and sperm motility (MOT) of dairy bulls in the Thai multibreed dairy population. Phenotypes for VOL (n = 13,535), NS (n = 12,773), and MOT (n = 12,660) came from 131 bulls of the Dairy Farming Promotion Organization of Thailand. Genotypic data consisted of 76,519 imputed and actual single nucleotide polymorphisms (SNP) from 72 animals. The SNP variances for VOL, NS, and MOT were estimated using a three-trait genomic-polygenic repeatability model. Fixed effects were contemporary group, ejaculate order, age of bull, ambient temperature, and heterosis. Random effects were animal additive genetic, permanent environmental, and residual. Individual SNP explaining at least 0.001% of the total genetic variance for each trait were selected to identify associated genes in the NCBI database (UMD Bos taurus 3.1 assembly) using the R package Map2NCBI. A set of 1,999 NCBI genes associated with all three semen traits was utilized for the pathway analysis conducted with the ClueGO plugin of Cytoscape using information from the Kyoto Encyclopedia of Genes and Genomes database. The pathway analysis revealed seven significant biological pathways involving 127 genes that explained 1.04% of the genetic variance for VOL, NS, and MOT. These genes were known to affect cell structure, motility, migration, proliferation, differentiation, survival, apoptosis, signal transduction, oxytocin release, calcium channel, neural development, and immune system functions related to sperm morphology and physiology during spermatogenesis.

Functional haplotypes of ARID4A affect promoter activity and semen quality of bulls

The AT-rich interaction domain 4 A (ARID4A) has an important role in regulating Sertoli cell function and male fertility. Its molecular mechanisms, however, remain largely unknown. In this study, two single nucleotide polymorphisms (SNPs) (g.53 G > T, ss 1966531596, and g.826 G > A, rs 210809648) were identified in the promoter region of ARID4A in 215 Chinese Holstein bulls using polymerase chain reaction (PCR)-restriction fragment length polymorphism and created restriction site-PCR. Results revealed that bulls with g.53 G > T-GG and g.826 G > A-G G genotype exhibited higher sperm deformity rate than those with g.53 G > T-TT and g.826 G > A-AA genotype (P <  0.01). Furthermore, three haplotypes (H1 (GG), H3 (TG), H4 (TA)) and six haplotype combinations (H1H1, H1H3, H1H4, H3H3, H3H4, H4H4) were obtained. The bulls with H4H4 exhibited lower sperm deformity rate than those with H1H1 and H1H3 (P <  0.05). In addition, results of bioinformatics analysis revealed that ARID4A has two promoters and that two SNPs of ARID4A are located in transcription factor binding sites. Compared with g.53 G > T-G and g.826 G > A-G allele, there was a greater fluorescence intensity in g.53 G > T-T and g.826 G > A-A allele by transient transfection in MLTC-1 cells and the luciferase report assay. qRT-PCR indicated the ARID4A expression was greater in bull spermatozoa with H4H4 haplotype combination than those with H1H1 haplotype combination (P < 0.05). Results of the present study indicate that g.53 G > T and g.826 G > A are functional mutations that are involved in regulation of ARID4A gene expression by affecting promoter activity and thus semen quality of Chinese Holstein bulls.

Identification of genome regions determining semen quality in Holstein-Friesian bulls using information theory

Use of information theory can be an alternative statistical approach to detect genome regions and candidate genes that are associated with livestock traits. The aim of this study was to verify the validity of the SNPs effects on some semen quality variables of bulls using entropy analysis. Records from 288 Holstein-Friesian bulls from one AI station were included. The following semen quality variables were analyzed: CASA kinematic variables of sperm (total motility, average path velocity, straight line velocity, curvilinear velocity, amplitude of lateral head displacement, beat cross frequency, straightness, linearity), sperm membrane integrity (plazmolema, mitochondrial function), sperm ATP content. Molecular data included 48,192 SNPs. After filtering (call rate = 0.95 and MAF = 0.05), 34,794 SNPs were included in the entropy analysis. The entropy and conditional entropy were estimated for each SNP. Conditional entropy quantifies the remaining uncertainty about values of the variable with the knowledge of SNP. The most informative SNPs for each variable were determined. The computations were performed using the R statistical package. A majority of the loci had relatively small contributions. The most informative SNPs for all variables were mainly located on chromosomes: 3, 4, 5 and 16. The results from the study indicate that important genome regions and candidate genes that determine semen quality variables in bulls are located on a number of chromosomes. Some detected clusters of SNPs were located in RNA (U6 and 5S_rRNA) for all the variables for which analysis occurred. Associations between PARK2 as well GALNT13 genes and some semen characteristics were also detected.

Seminal plasma and sperm proteome of ring-tailed coatis (Nasua nasua, Linnaeus, 1766)

Ring-tailed coati is listed as a species of least concern in the International Union for Conservation of Nature (IUCN) Red List, however, there has been a sharp decline in their population. The present study was conducted to evaluate the major proteins of both seminal plasma and sperm in ring-tailed coatis. Semen sample was collected from three adult coatis and evaluated for their morphological characteristics. Further, the sample was centrifuged to separate spermatozoa from seminal plasma, and then stored in liquid nitrogen. The seminal plasma and sperm proteins were subjected to one-dimensional (1-D) sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) and identified by mass spectrometry. Gene ontology and protein networks were analyzed using bioinformatics tools. Based on sperm concentration and average protein content of the semen, the concentration of protein/spermatozoon was found to be 104.69 ± 44.43 μg. The analysis of SDS-PAGE gels showed 20.3 ± 3.1 and 17 ± 2 protein bands/lane for seminal plasma and sperm, respectively. In-gel protein digestion and peptide analysis by mass spectrometry revealed 238 and 246 proteins in the seminal plasma and sperm, respectively. The gene ontology analysis revealed that the proteins of seminal plasma mainly participated in cellular (35%) and regulatory (21%) processes. According to their cellular localization, seminal plasma proteins were categorized as structural (18%), extracellular (17%), and nuclear (14%) proteins with molecular functions, such as catalytic activity (43%) and binding (43%). The sperm proteins were also involved in cellular (38%) and regulatory (23%) processes, and mainly categorized as extracellular (17%), nuclear (13%), and cytoplasmic (10%) proteins. The major molecular functions of the sperm proteins were catalytic activity (44%) and binding (42%). These results indicated that the seminal plasma of ring-tailed coati has an array of proteins that can potentially modulate several sperm functions, from sperm protection to oocyte binding. However, further studies are necessary to interpret the roles of these major seminal plasma proteins in coatis.

Genomewide association study for production and meat quality traits in Canchim beef cattle

The commercial value of the bovine carcass is determined by a set of traits, such as weight, yield, back fat thickness, and marbling; therefore, the genetic improvement of growth, meat, and carcass quality traits is an important tool to add value to the supply chain. Genomewide association studies (GWAS) enable the identification of loci that control phenotypic expression of quantitative traits (QTL). Therefore, the objective of this work was to perform a GWAS to identify genomic regions and genes associated with growth, carcass traits, and meat quality in Canchim beef cattle. These traits were yearling weight (YW), rib eye area (REA), back fat thickness (BFT), and marbling (MARB). To increase sample size and marker density, genotype imputation was performed, and only markers imputed with greater than 95% accuracy were used. Genomewide association study was performed using a Bayesian approach, by the Bayes B statistical method, incorporating genotypes and phenotypes from 614 animals from both the Canchim breed and the MA genetic group (offspring of Charolais bulls and one-half Canchim + one-half Zebu cows). This investigation identified 1 and 4 genomic regions explaining 0.23 and 7.35% of the genetic variance for REA and YW, respectively. These regions harbor a total of 19 genes, 7 of which were classified for biological functions by functional analysis. Significant associations were not observed for BFT and MARB. The identification of QTL that had been previously described in the literature reinforces associations found in this study.

Genetic effects of PDGFRB and MARCH1 identified in GWAS revealing strong associations with semen production traits in Chinese Holstein bulls

Background

Using a genome-wide association study strategy, our previous study discovered 19 significant single-nucleotide polymorphisms (SNPs) related to semen production traits in Chinese Holstein bulls. Among them, three SNPs were within or close to the phosphodiesterase 3A (PDE3A), membrane associated ring-CH-type finger 1 (MARCH1) and platelet derived growth factor receptor beta (PDGFRB) genes. The present study was designed with the objectives of identifying genetic polymorphism of the PDE3A, PDGFRB and MARCH1 genes and their effects on semen production traits in a Holstein bull population.

Results

A total of 20 SNPs were detected and genotyped in 730 bulls. Association analyses using de-regressed estimated breeding values of each semen production trait revealed four statistically significant SNPs for one or more semen production traits (P < 0.05): one SNP was located downstream of PDGFRB and three SNPs were located in the promoter of MARCH1. Interestingly, for MARCH1, haplotype-based analysis revealed significant associations of haplotypes with semen volume per ejaculate. Furthermore, high expression of the MARCH1 gene was observed in sperm cells. One SNP (rs43445726) in the regulatory region of MARCH1 had a significant effect on gene expression.

Conclusion

Our study demonstrated the significant associations of genetic variants of the PDGFRB and MARCH1 genes with semen production traits. The identified SNPs may serve as genetic markers to optimize breeding programs for semen production traits in Holstein bull populations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-017-0527-1) contains supplementary material, which is available to authorized users.

Gene structure of the pregnancy-associated glycoprotein-like (PAG-L) in the Eurasian beaver (Castor fiber L.)

The pregnancy-associated glycoprotein-like family (PAG-L) is a large group of chorionic products, expressed in the pre-placental trophoblast and later in the post-implantational chorionic epithelium, and are involved in proper placenta development and embryo-maternal interaction in eutherians. This study describes identification of the PAG-L family in the genome of the Eurasian beaver (Castor fiber L.), named CfPAG-L. We identified 7657 bp of the CfPAG-L gDNA sequence (Acc. No. KX377932), encompassing nine exons (1-9) and eight introns (A-H). The length of the CfPAG-L exons (59-200 bp) was equivalently similar to the only known counterparts of bPAG1, bPAG2, and pPAG2. The length of the CfPAG-L introns ranged 288-1937 bp and was completely different from previously known PAG introns. The exonic CfPAG-L regions revealed 50.3-72.9% homology with equivalent segments of bPAG1 and pPAG2 structure. The intronic CfPAG-L regions alignments revealed a lack of homology. Within the entire CfPAG-L gene, 31 potential single nucleotide variants (SNV: 7 transversions and 24 transitions) were predicted. The identified exonic polymorphic loci did not affect the amino acid sequence of the CfPAG-L polypeptide precursor. This is the first report describing the CfPAG-L gene sequence, structural organization, and SNVs in the Eurasian beaver, one of the largest rodents.

Functional SNPs of INCENP Affect Semen Quality by Alternative Splicing Mode and Binding Affinity with the Target Bta-miR-378 in Chinese Holstein Bulls

Inner centromere protein (INCENP) plays an important role in mitosis and meiosis as the main member of chromosomal passenger protein complex (CPC). To investigate the functional markers of the INCENP gene associated with semen quality, the single nucleotide polymorphisms (SNPs) g.19970 A>G and g.34078 T>G were identified and analyzed. The new splice variant INCENP-TV is characterized by the deletion of exon 12. The g.19970 A>G in the exonic splicing enhancer (ESE) motif region results in an aberrant splice variant by constructing two minigene expression vectors using the pSPL3 exon capturing vector and transfecting vectors into MLTC-1 cells. INCENP-TV was more highly expressed than INCENP-reference in adult bull testes. The g.34078 T>G located in the binding region of bta-miR-378 could affect the expression of INCENP, which was verified by luciferase assay. To analyze comprehensively the correlation of SNPs with sperm quality, haplotype combinations constructed by g.19970 A>G and g.34078 T>G, as well as g.-692 C>T and g.-556 G>T reported in our previous studies, were analyzed. The bulls with H1H12 and H2H2 exhibited a higher ejaculate volume than those with H2H10 and H9H12, respectively (P < 0.05). Bulls with H11H11 and H2H10 exhibited higher initial sperm motility than those with H2H2 (P < 0.05). The expression levels of INCENP in bulls with H1H12 and H11H11 were significantly higher than those in bulls with H9H12 (P < 0.05), as determined by qRT-PCR. Findings suggest that g.19970 A>G and g.34078 T>G in INCENP both of which appear to change the molecular and biological characteristics of the mRNA transcribed from the locus may serve as a biomarkers of male bovine fertility by affecting alternative splicing mode and binding affinity with the target bta-miR-378.

Genome-wide association study for sperm membrane integrity in frozen-thawed semen of Holstein-Friesian bulls

The aim of the study was to screen the entire bull genome to identify SNP markers and propose candidate genes potentially involved in the variation of sperm membrane integrity in Holstein-Friesian bulls. Two hundred eighty eight bulls kept in one AI center were included in the study. Each bull was genotyped for 54.001 Single Nucleotide Polymorpisms (SNP) by the Illumina BovineSNP50 BeadChip. Commercial straws of frozen-thawed semen were used for the evaluation of sperm plasma membrane integrity (SYBR-14/PI staining) and sperm mitochondrial function (JC1/PI staining). An additive model for Linear Regression Analysis was applied to estimate the effect of SNP marker for sperm membrane integrity (by the use of GoldenHelix SVS7 software). Five significant markers (encompassing 2,2 MB region located on chromosome 6) for SYBR-14/PI were found. Among them one marker-rs41570391 passed Bonferroni correction test. Within approximately 3 Mb genomic region including significant markers three candidate genes: SGMS2 (Sphingomyelin Synthase 2), TET2 (Methylcytosine dioxygenase 2) and GSTCD genes (Gluthatione S-transferase C terminal domain) were proposed as potentially involved in sperm membrane integrity in frozen-thawed semen of Holstein-Friesian bulls.

Variations in Antioxidant Genes and Male Infertility

Oxidative stress and reactive oxygen species (ROS) are generated from both endogenous and environmental resources, which in turn may cause defective spermatogenesis and male infertility. Antioxidant genes, which include catalase (CAT), glutathione peroxidase (GPX), glutathione S-transferase (GST), nitric oxide synthase (NOS), nuclear factor erythroid 2-related factor 2 (NRF2), and superoxide dismutase (SOD), play important roles in spermatogenesis and normal sperm function. In this review, we discuss the association between variations in major antioxidant genes and male infertility. Numerous studies have suggested that genetic disruption or functional polymorphisms in these antioxidant genes are associated with a higher risk for male infertility, which include low sperm quality, oligoasthenoteratozoospermia, oligozoospermia, and subfertility. The synergistic effects of environmental ROS and functional polymorphisms on antioxidant genes that result in male infertility have also been reported. Therefore, variants in antioxidant genes, which independently or synergistically occur with environmental ROS, affect spermatogenesis and contribute to the occurrence of male infertility. Large cohort and multiple center-based population studies to identify new antioxidant genetic variants that increase susceptibility to male infertility as well as validate its potential as genetic markers for diagnosis and risk assessment for male infertility for precise clinical approaches are warranted.