Next Article in Journal
Exposure to Dichlorodiphenyldichloroethylene (DDE) and Metallothionein Levels in Rats Fed with Normocaloric or High-Fat Diet: A Review
Previous Article in Journal
Primary Human Chondrocytes Affected by Cigarette Smoke—Therapeutic Challenges
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 21
Issue 5
10.3390/ijms21051902
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Gene Polymorphisms in Boar Spermatozoa and Their Associations with Post-Thaw Semen Quality

by
Anna Mańkowska
1,
Paweł Brym
2,
Łukasz Paukszto
3,
Jan P. Jastrzębski
3 and
Leyland Fraser
1,*
1
Department of Animal Biochemistry and Biotechnology, Faculty of Animal Bioengineering, University of Warmia and Mazury in Olsztyn, 10-719 Olsztyn, Poland
2
Department of Animal Genetics, Faculty of Animal Bioengineering, University of Warmia and Mazury in Olsztyn, 10-719 Olsztyn, Poland
3
Department of Plant Physiology and Biotechnology, Faculty of Biology and Biotechnology, University of Warmia and Mazury in Olsztyn, 10-719 Olsztyn, Poland
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2020, 21(5), 1902; https://doi.org/10.3390/ijms21051902
Submission received: 1 February 2020 / Revised: 4 March 2020 / Accepted: 6 March 2020 / Published: 10 March 2020
(This article belongs to the Section Biochemistry)
Browse Figures
Versions Notes

Abstract

:
Genetic markers have been used to assess the freezability of semen. With the advancement in molecular genetic techniques, it is possible to assess the relationships between sperm functions and gene polymorphisms. In this study, variant calling analysis of RNA-Seq datasets was used to identify single nucleotide polymorphisms (SNPs) in boar spermatozoa and to explore the associations between SNPs and post-thaw semen quality. Assessment of post-thaw sperm quality characteristics showed that 21 boars were considered as having good semen freezability (GSF), while 19 boars were classified as having poor semen freezability (PSF). Variant calling demonstrated that most of the polymorphisms (67%) detected in boar spermatozoa were at the 3’-untranslated regions (3’-UTRs). Analysis of SNP abundance in various functional gene categories showed that gene ontology (GO) terms were related to response to stress, motility, metabolism, reproduction, and embryo development. Genomic DNA was isolated from sperm samples of 40 boars. Forty SNPs were selected and genotyped, and several SNPs were significantly associated with motility and membrane integrity of frozen-thawed (FT) spermatozoa. Polymorphism in SCLT1 gene was associated with significantly higher motility and plasma membrane integrity of FT spermatozoa from boars of the GSF group compared with those of the PSF group. Likewise, polymorphisms in MAP3K20, MS4A2, and ROBO1 genes were significantly associated with reduced cryo-induced lipid peroxidation and DNA damage of FT spermatozoa from boars of the GSF group. Candidate genes with significant SNP associations, including APPL1, PLBD1, FBXO16, EML5, RAB3C, OXSR1, PRICKLE1, and MAP3K20 genes, represent potential markers for post-thaw semen quality, and they might be relevant for future improvement in the selection procedure of boars for cryopreservation. The findings of this study provide evidence indicating that polymorphisms in genes expressed in spermatozoa could be considered as factors associated with post-thaw semen quality.

1. Introduction

Cryopreservation of semen allows for prolonged storage of genetically important reproductive traits through the use of assisted reproductive techniques (ARTs), such as artificial insemination (AI) [1,2]. Due to poor post-thaw semen quality and reduced fertility, the widespread application of frozen-thawed (FT) boar semen in the AI practices is limited compared to liquid-stored semen [2]. Selection of high proportions of viable FT spermatozoa is the major challenge in the pig AI industry [3]. Moreover, the selection of a high number of functionally viable FT boar spermatozoa is required to minimize economic loss and increase the worldwide application of FT semen in AI technologies [1,2].
Analysis of different sperm attributes has confirmed that cryo-induced damage to spermatozoa differs among individual boars, suggesting varying responses to the freezing–thawing process [4,5]. Evidence has indicated that the freezability of boar semen is influenced by genetic factors [4,6]. Recent developments in high-throughput sequencing techniques, such as transcriptome sequencing (RNA-Seq), have enabled a thorough analysis of gene expression and genetic variations in the pig reproductive tract [7]. Our recent study has confirmed that there are differences in the transcriptome profiles of boar spermatozoa [8]. Furthermore, it has been confirmed that boars with poor freezability ejaculates are characterized by an overexpression of differentially expressed genes (DEGs) that are mainly associated with inflammation and apoptosis, which increase the sperm susceptibility to cryo-induced damage [8]. Previous studies confirmed that genes enriched in the cytokine–cytokine receptor interaction and inflammatory response-related pathways predisposed boar sperm to cryo-induced damage [9,10]. Hence, transcriptome alterations in spermatozoa have been considered to be one of the main factors affecting the cryo-survival of boar semen [8,9,10]. It should be emphasized that variant calling procedure based on RNA-Seq data has been suggested to be an important screening tool to identify polymorphisms in differentially expressed genes [7] and potential genetic markers associated with production traits in the pig industry [11].
Several authors have reported an increasing number of markers related to reproduction traits in pigs [12,13,14]. Advances in molecular genetics have resulted in the identification of several single nucleotide polymorphisms (SNPs) that are associated with the quality parameters of boar spermatozoa (including motility and morphology), and there has been increasing interest to understand the molecular processes that affect the sperm phenotype traits [12,15]. Evidence has shown that polymorphisms in boar sperm could be used as markers for semen quality [12,15,16]. However, to date, there is limited information about the associations of SNPs with post-thaw quality of boar semen. To our knowledge, no study has yet explored the associations between SNP markers and post-thaw quality semen of the Polish large white (PLW) boars. The objectives of this study were to i) identify polymorphisms in candidate genes (CGs) of boar spermatozoa using variant calling analysis of RNA-Seq dataset, ii) identify the functional categories of SNP abundance in the CGs by gene ontology (GO) enrichment analysis, iii) validate the selected SNPs using the Kompetitive Allele Specific Polymorphism (KASP) genotyping assay, and iv) explore the associations of gene polymorphisms with post-thaw (PT) semen quality.

2. Results

2.1. Assessment of Semen Quality

Repeated measures analysis of variance (ANOVA) did not show significant differences (p > 0.05) among the boars, with respect to fresh, pre-freeze semen quality. However, ANOVA results demonstrated that boar was a significant (p < 0.001) factor affecting post-thaw motility, mitochondrial membrane potential (MMP), plasma membrane integrity (PMI), normal apical ridge (NAR) acrosome integrity, DNA fragmentation, and lipid peroxidation (LPO). The quality of fresh, pre-freeze, and post-thaw semen for 40 boars is shown in Table 1. Boars showing more than 30% (>30%) total motility were considered as having good semen freezability (GSF), whereas boars with motility less than 30% (<30%) were considered as having poor semen freezability (PSF). Based on post-thaw sperm analysis, the semen of 21 boars was classified as having GSF, whereas the semen of 19 boars was classified as having PSF (Table 1). Besides post-thaw motility, membrane integrity of FT spermatozoa was significantly higher (p < 0.05) in boars of the GSF group compared with those of the PSF group (Table 1).

2.2. RNA-Seq Mapping and Variant Calling

Clean reads that were mapped to the reference genome ranged from 38,651,853 to 58,239,476, corresponding to about 72.8% uniquely mapped reads (range, 69.8% to 76.7%). The results of the mapping procedure were merged into a single binary alignment map (BAM) file for each boar. The fastq sequence dataset of each library is accessible in the NCBI-SRA database (Bioproject: PRJNA415904; accession number: SRP121647).
Using rMATS-DVR software, in conjunction with the genome analysis tool kit (GATK), a total of 1,389,568 variants were identified in the sperm RNA-Seq datasets of the six boars. Analysis of the filtered/trimmed data base single nucleotide polymorphisms (dbSNPs) resulted in an average of 599,428 raw SNPs, corresponding to approximately 88,340 SNPs (≥10% reads ≥90%) and 82,356 SNPs (≥10% reads = 100%), as shown in Table 2. Following stringent filtering analysis of the dbSNPs, 1371 SNP variants were annotated (Supplementary Table S1), of which about 67% putative polymorphisms (919/1371) were located at the 3’-untranslated regions (3’-UTRs) (Table 3). Furthermore, variant calling detected approximately 5.4% missense polymorphism (74/1371) in genes expressed in boar spermatozoa (Table 3).

2.3. SNP Abundance in Functional Categories of Candidate Genes (CGs)

A total of approximately 180 CG genes were enriched in at least 27 pathways, and were mainly distributed in the inflammatory-related pathways followed by the glutamatergic synapse and signaling pathways (Supplementary Table S2). Furthermore, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway showed a high number of the CG genes were enriched in the PI3K-Akt signaling and metabolic pathways (Supplementary Table S2).
Enrichment analysis showed the gene ontology (GO) terms of the CGs (Figure 1, Supplementary Table S3). Generally, the CGs were associated with various functional categories, such as biological process categories (Figure 1A), molecular function categories (Figure 1B), and cellular component categories (Figure 1C). The more prominent enriched terms of the GO biological process categories included “cellular response to stimulus” (APPL1, TXNIP, and OXSR1) and “cell communication” (MYO3B, EXOC4, and KIF1B), as shown in (Figure 1A). The GO molecular function categories were mainly associated with “hydrolase activity” (PARG, GBA3, and GM2A) and “catalytic activity, acting on a protein” (HARS2, PLBD1, and SARS) (Figure 1B). In addition, the GO cellular component categories included enriched terms associated with “cell” (CFAP52, ANKRD50, and SKAP2) and “intracellular” (EVA1A, TRIM9, and PRICKLE1) (Figure 1C).
The selected CGs for KASP genotyping markers are highly associated with various reproductive processes and sperm functions (Supplementary Table S4). Analysis of the reproduction quantitative traits (QTLs) of the selected CGs showed that age at puberty, corpus luteum number, litter size, and number of stillborns were commonly detected among the pig reproductive traits, according to the Animal Genome pigQTLdb (https://www.animalgenome.org/cgi-bin/QTLdb/SS/index).

2.4. KASP Genotyping and Validation

Following KASP genotyping, the allele call rate was 92%, and the genotyping of three SNP markers—AHI1, GLMN, and IFNAR2—was unsuccessful, probably due to poor design of the primers. It was found that MY3OB was monomorphic, whereas three other SNP markers—CDK17, CLNK, and CRISP2—showed an excess of homozygotes or heterozygotes in either freezability group. These SNPs were removed from further analysis. Among the SNPs analyzed in this study, the χ2-test and probability estimates showed that most of the SNP loci were within in the Hardy–Weinberg equilibrium (HWE) expectations (Table 4). Noticeable significant departures (p < 0.05) from HWE were common in either the GSF group or PSF group, or in both freezability groups, such as at A2M g.62485457G > A (rs339026428) and COMMD2 g.189987893G > A (rs318435440) loci (Table 4).
There were variations in the frequency distributions of the alleles of the CGs (Table 4). Irrespective of the freezability group, higher frequency of allele G was observed in CYP7B1 (0.90) gene locus (Table 4A), whereas allele A was more frequent in CFAP52 (0.72) gene locus (Table 4B). Higher frequency of allele C was observed in EML6 (0.72) gene locus (Table 4C), whereas higher frequency distributions of allele T were observed in ROBO1 (0.87) gene locus (Table 4F).
Genotyping cluster plots of a few genes with polymorphism used for KASP genotyping assay are shown in Figure 2A–F (APPL1, OXSR1, FBXO16, RAB3C, PLBD1, and SARS). Samples shown in pink dots were not assigned to the cluster.

2.5. Association Analysis of SNPs with Post-Thaw Semen Quality

ANOVA analysis showed that polymorphisms in APPL1 (p < 0.017), CYP7B1 (p < 0.031), RAB3C (p < 0.043), OXSR1 (p < 0.047), FBXO16 (p < 0.045), and PLBD1 (p < 0.046) genes had significant effects on post-thaw sperm motility. Polymorphisms in APPL1 and PLBD1 genes markedly affected post-thaw PMI (p < 0.048 and p < 0.025, respectively) and NAR acrosome integrity (p < 0.016 and p < 0.018, respectively). Also, polymorphism in PLBD1 gene affected post-thaw MMP and LPO (p < 0.018 and p < 0.043, respectively). It was observed that polymorphisms in SARS gene had marked effects on post-thaw MMP (p < 0.046) and NAR acrosome integrity (p < 0.001), while EML5 polymorphism affected post-thaw DNA integrity (p < 0.015), LPO (p < 0.035), and NAR acrosome integrity (p < 0.048). Furthermore, polymorphisms in RAB3C and PRICKLE1 genes had marked effects on post-thaw NAR acrosome integrity (p < 0.041 and p < 0.040, respectively) and LPO (p < 0.037 and p < 0.013, respectively). In addition, polymorphisms in FBXO16 gene affected post-thaw MMP (p < 0.048).
Association analysis of post-thaw semen quality showed that motility was lower (p < 0.05) for AG genotype of APPL1 gene (a newly predicted SNP) (Figure 3A). Likewise, post-thaw motility was lower (p < 0.05) for AG genotype (AG vs. GG) at rs338842672 (CYP7B1 gene) and rs81210636 (RAB3C gene) (Figure 3B,C, respectively). Furthermore, post-thaw motility was lower (p < 0.05) for the heterozygous genotype at rs339379734 (OXSR1 gene) compared with AA genotype (Figure 3D), and for GG genotype at rs341614458 (FBXO16 gene), as shown in Figure 3E. As regards PLBD1 gene polymorphism, post-thaw motility was the lowest for TT genotype at rs321497623 (Figure 3F).
Significantly lower (p < 0.05) post-thaw MMP (Figure 4A), PMI (Figure 4B), and acrosome integrity (Figure 4C) were concurrent with higher post-thaw LPO for TT genotype at rs321497623 (Figure 4D).
Boars with AG genotype of APPL1 gene had significantly lower (p < 0.05) post-thaw PMI and NAR acrosome integrity (Figure 5A,B, respectively), while lower post-thaw NAR acrosome integrity and higher LPO were observed for AA genotype at rs345056502 (EML5 gene) (Figure 5C,D, respectively). Furthermore, post-thaw DNA damage was highest (p < 0.05) for AA genotype at rs345056502 (Figure 5E), while lower post-thaw MMP was observed for GG genotype at rs341614458 (Figure 5F).
Post-thaw analysis showed that boars with AG genotype had lower (p < 0.05) MMP and acrosome integrity at rs344846507 (SARS gene) (Figure 6A,B, respectively). It was found that lower (p < 0.05) post-thaw NAR acrosome integrity was concomitant with higher LPO for genotype AG (AG vs. GG) at rs81210636 (Figure 6C,D, respectively). Higher post-thaw acrosome damage and LPO were associated with CT or TT genotype at rs694366781 (PRICKLE1 gene) (Figure 6E,F, respectively).
ANOVA analysis revealed that freezability (F) and SNP genotype (G) and their interactions (F × G) significantly affected (p < 0.05) post-thaw semen quality (Table 5). Significant F × G interactions with post-thaw motility and PMI were observed for rs337913978 (SCLT1 gene, Table 5A). Also, freezability (p < 0.001) and F × G interaction (p < 0.010) with post-thaw NAR acrosome integrity were observed for rs340643892 (MAP3K20 gene, Table 5B). Post-thaw LPO and DNA fragmentation showed significant SNP effects and F × G interaction for rs339836492 (MS4A2 gene) and rs331568674 (ROBO1 gene), respectively (Table 5C,D, respectively).
Association analysis of post-thaw sperm quality showed that PSF boars with CT and TT genotypes at rs337913978 (SCLT1 gene) exhibited lower (p < 0.05) motility and PMI (Figure 7A,B, respectively). Significantly (p < 0.05) lower proportions of FT spermatozoa with acrosome integrity were concomitant with higher LPO for AA genotype at rs340643892 (MAP3K20 gene) in boars from the PSF group (Figure 7C,D, respectively). Post-thaw LPO was higher (p < 0.05) at rs339836492 (MS4A2 gene) (Figure 7E), while post-thaw DNA damage was greater (p < 0.05) for TT genotype at rs331568674 (ROBO1 gene) in boars of the PSF group (Figure 7F).

3. Discussion

3.1. SNP Functional Classes and Validation

In the present study, among the RNA-Seq SNPs (Supplementary Table S1), approximately 82% of the putative polymorphisms in boar spermatozoa had a corresponding dbSNP entry. In a previously published paper, about 88–91% of the detected SNPs in the testis tissue of boars had a dbSNP entry, and a large number of SNPs were at the 3’-UTRs, the untranslated region of noncoding mRNAs [7]. In this study, we found that about 67% of putative polymorphisms were at the 3’-UTR (Table 3), and the impact of most of these polymorphisms on sperm function is not fully understood. Studies have reported that the 3’-UTR plays an important role in the translation efficiency and stability of mRNA [17,18]. It has been reported that UTRs have strong impacts on post-transcriptional regulation of gene expression, suggesting that dysfunction in the UTRs might have a significant effect on gene expression and the associated cellular viability, growth, and development [17,18]. Among the UTR variants at the 3’-UTR, two SNP polymorphisms, rs339379734 (OXSR1) and rs337913978 (SCLT1), showed significant associations with post-thaw semen quality. Evidence has indicated that SNPs at the 3’-UTR of the targeted mRNA have been the principal elements of microRNA adhesion, and they are associated with semen quality [19]. We suggest that the 3’-UTR SNPs, detected in this study, might play a significant role at the transcription level. Furthermore, most of the coding sequences detected in spermatozoa were synonymous variants, which do not affect protein sequence, but could be in linkage disequilibrium to other causative mutations [7,12] and could have functional effects on mRNA stability and, ultimately, the phenotype traits [20]. In the present study, eight SNP markers that showed significant associations with post-thaw semen quality were synonymous variants (in CYP7B1, FBXO16, RAB3C, SARS, MAP3K20, PLBD1, EML5, and PROCKLE1 genes), which could be linked to causative mutations with functional effects on the reproductive traits in the pig. Furthermore, two missense polymorphisms in APPL1 and MS4A2 genes showed significant associations with post-thaw semen quality. Missense variants affect amino acid translation, protein structure and function, and several missense polymorphisms have been shown to affect enzyme activity [21]. It is reasonable to deduce that the missense polymorphisms in APPL1 (a novel predicted SNP) and MS4A2 (rs339836492) genes could have a marked effect on their protein expression levels, which could consequently compromise the sperm response to the freezing–thawing conditions.
Quality assessment of a panel of SNP markers was performed by KASP genotyping assay. It is noteworthy that KASP is a novel competitive allele specific PCR, which is mainly based on the amplification of DNA with a thermal cycler using allele specific primers [22]. Due to its low cost and high efficiency, the KASP genotyping assay has been widely used in SNP genotyping studies related to molecular marker-assisted selection breeding [23]. In this study, the genotype frequencies were consistent with the HWE assumptions for most of the analyzed loci. Screening of SNPs for HWE departure is frequently used when performing association studies [24]. There is little consensus on the appropriate p-value threshold for the identification of SNPs that violate the HWE assumption in association studies [25]. In the current study, it is unclear why HWE departure was more common in gene loci of individuals from the PSF group. It seems likely that the marked departure from the HWE could be related to long-term artificial selection and breeding of pigs. Several authors suggested that long-term selective breeding for production traits, not associated with reproductive traits, could be attributed to the marked differences in the genotype frequencies [14,26]. Moreover, it has been confirmed that differences in the allele frequencies in AI boars highly selected for a number of traits that are responsible for HWE departure [26]. It is likely that HWE departure is not necessarily caused by the selection of the analyzed traits, but rather other traits under the control of the selected genes or closely related genes [26]. We suggest that further studies are needed to find out whether such a phenomenon is a common feature of polymorphisms in genes expressed in poor freezability ejaculates.

3.2. SNP Associations with Post-Thaw Semen Quality

Evidence has shown that the evaluation of various sperm traits rather than a single trait analysis provides a better fertility prediction of FT semen [27,28,29]. Several sperm attributes have been used to assess the post-thaw quality of boar semen, and among these sperm motility is widely used in semen quality assessment because it correlates with fertility outcomes [30,31]. Sperm motility is assessed subjectively or with the computer-assisted semen analysis (CASA) system; however, subjective motility assessment has been shown to provide reliable estimates in association studies [26,32,33]. In the current study, several sperm attributes have been used to monitor the associations of gene polymorphisms with post-thaw quality of boar semen.
The distribution of genotype frequencies differs among the genotype groups, with respect to polymorphisms in the genes related to GO term, “response to stress”, such as APPL1, OXSR1 (rs339379734), MAP3K20 (rs340643892), and MS4A2 (rs339836492). Since stress plays an important role in semen cryopreservation [34], the function and localization of these genes are critical in understanding their impact on post-thaw semen quality. In general, the heterogeneous genotype of the APPL1 gene polymorphism has been associated with reduced post-thaw semen quality. It is noteworthy that the APPL1 polymorphism is a newly predicted SNP variant detected in boar spermatozoa, and it is very difficult to make comparisons with other studies. The APPL1 gene is related to several GO terms, such as “signal transduction”, “cell motility”, and “transmembrane transport”, suggesting its significant relevance in sperm function. Accumulating evidence has shown that APPL1 regulates the function of the adiponectin receptor (AdipoR1), which is implicated in various biological processes in the reproductive system [35,36]. Adiponectin is abundantly found in the tail region of bull spermatozoa, whereas AdipoR1 has been detected mainly in the equatorial and acrosome regions of sperm [35]. Furthermore, a direct role of the AdipoR1 receptor has been suggested in sperm capacitation [36], and adiponectin and its receptors have been associated with motility parameters of ram spermatozoa [37]. These findings reaffirm the important modulating function of the APPL1 gene on the action of the AdipoR1 receptor in regulating sperm motility. Consistent with this notion, it seems that polymorphisms in the APPL1 gene might be a contributing factor to compromised function of FT spermatozoa.
The OXSR1 gene regulates downstream kinases in response to environmental stress, and it plays a role in the Na+-K+-2Cl− cotransporter (NKCC1) in the testis [38]. Moreover, OXSR1 is one of the upstream phosphorylators of the NKCC1, which is activated through a phosphorylation-dependent mechanism and is implicated in sperm–oocyte fertilization events [38]. It seems, therefore, that an OXSR1 gene polymorphism could reduce the abundance of phosphorylated-NKCC1, resulting in reduced sperm motility and compromised IVF potential [38]. Association analysis shows that heterogeneous genotype of OXSR1 gene significantly increases the sperm susceptibility to reduced post-thaw motility. Hence, it seems likely that the OXSR1 gene polymorphism compromises the sperm’s ability to confer protection against oxidative stress, particularly during the freezing–thawing process. Cryo-induced stress increases the sperm susceptibility to reactive oxygen species (ROS)-related damage, such as mitochondrial dysfunction, peroxidation of sperm membrane lipids, and DNA damage, resulting in compromised fertility of FT spermatozoa [5,29,39]. Furthermore, it seems that the AA genotype at rs340643892 (MAP3K20 gene) might predispose spermatozoa to increased LPO and acrosome damage after freezing–thawing. It should be noted that MAPK is a serine/threonine kinase, which is implicated in the activation of gene transcription and expression [40]. Moreover, the MAPK cascade elements have been detected in the flagellum of human spermatozoa and are involved in the regulation of capacitation and the acrosome reaction processes [40]. Analysis shows that GC and CC genotypes at rs339836492 (MS4A2 gene) were associated with increased post-thaw LPO in poor freezability ejaculates, suggesting that spermatozoa from these ejaculates were more susceptible to increased damage incurred by ROS action during the freezing–thawing process. Although the functions of many of the MS4A proteins are currently not well defined, it was reported that the MS4A family member could form oligomers in sperm membranes, which might be involved in the interaction with the zona pellucida or cumulus during fertilization [41]. It is noteworthy that the MS4A2 transcript has been detected in the sperm head [42], suggesting its relevance in sperm–oocyte interaction mechanisms [41]. Our results show that polymorphisms in MAP3K20 and MS4A2 genes are associated with reduced post-thaw semen quality, which was manifested mainly in increased susceptibility of FT spermatozoa to LPO, as observed in poor freezability ejaculates. Moreover, an increase in LPO of FT spermatozoa results in a substantial loss of membrane integrity and motility, and ultimately leads to reduced fertility [43]. Additional studies are required to determine the biological relevance of these gene polymorphisms in sperm functions.
Interestingly, it has been confirmed that heterozygous AG (AG vs. GG) of CYP7B1 gene (rs338842672) is a contributing factor to reduced post-thaw motility. The CYP7B gene is enriched in the primary bile acid biosynthesis pathway (ssc00120) and 14 GO terms including “reproduction”, “oxidoreductase activity”, “cell motility”, and “signal transduction”, suggesting the important role of the gene in sperm functions. However, we would like to emphasize that the consequences of CYP7B1 polymorphism on post-thaw quality of boar semen are not yet known. Similar to the CYP7B1 gene, only two genotypes (AG vs. GG) of the RAB3C gene were considered in the association studies with post-thaw semen quality. In support of this notion, analysis showed that heterozygous AG at rs81210636 (RAB3C gene) is also a contributing factor to reduced post-thaw semen quality compared with the GG genotype. The impact of RAB3C gene on post-thaw semen quality is unclear. Not much information is available about the role of the RAB3C gene in sperm function, but earlier research reported that the RAB3 protein could interact with a series of proteins and could be associated with the cAMP/PKA (cAMP/protein kinase A) messenger system, in conjunction with phospholipase A2 (PLA2), to modulate exocytosis of the sperm acrosome [44]. Moreover, the RAB3 gene is related to 7 GO enriched terms, such as “GTPase activity”, “signal transduction”, and “vesicle-mediated transport”, suggesting the important role of the gene in sperm biological functions. It seems that the higher genotype frequency of heterozygous AG (60%) might predispose FT spermatozoa to compromise fertilization-related events. However, it is unclear how the RAB3C protein interacts with boar sperm membrane, and its role in cryo-survival has not yet been elucidated.
In our study, polymorphisms in FBXO16 (rs341614458) and PLBD1 (rs321497623) genes were associated with motility and membrane integrity of FT spermatozoa. It should be emphasized that there is not enough information in the literature about the specific role of these genes in the physiological functions of the sperm. It seems that the GG genotype of the FBXO16 gene is a contributing factor for reduced post-thaw semen quality. The effect of the FBXO16 polymorphism on sperm function is not fully known; however, members of the F-box family bind to phosphorylated proteins to promote their ubiquitination and degradation, which are required to protect and maintain sperm quality [45,46]. Evidence has shown that ubiquitin-mediated proteolysis of FBXO proteins is indispensable for the stability of sperm organelles [45]. Therefore, dysfunction in the ubiquitin-related pathway might have relevance in the physiological functions of the sperm. Furthermore, the effects of PLBD1 gene polymorphism on the biological function of the sperm is still unclear, even though it has been reported that PLBD1 is implicated in the fertilization processes [47]. Moreover, the PLBDI gene is related to 12 GO enriched terms, such as “lipid metabolic process”, “signal transduction”, and “kinase activity”, suggesting the important role of the gene in various sperm biological functions. It can be suggested that an understanding of the functional significance of FBXO16 and PLBD1 genes is necessary to elucidate their biological effects on sperm functions following cryopreservation.
Analysis showed that polymorphisms in the genes associated with cytoskeleton proteins, namely EML5 (rs345056502) and PRICKLE1 (rs694366781), were significantly correlated with reduced motility and acrosome integrity, and increased LPO of FT spermatozoa. It is noteworthy that the cytoskeleton is implicated in various cellular functions including mitosis, membrane translocations, and cellular motility [48]. The impact of polymorphisms of these genes on sperm function is not fully understood. However, our results provide evidence indicating that boars with the homologous AA genotype at rs345056502 (EML5 gene) were more susceptible to post-thaw acrosome damage and LPO, while the homologous TT genotype of the PRICKLE1 gene contributed to greater cryo-induced DNA damage. Moreover, microtubule-associated proteins (MAPs) in sperm cells have various functions including modulation of actin cytoskeletal function during spermatogenesis [49]. Collectively, EML5 and PRICKLE1 are implicated in microtubule (MT) dynamics in the testis, and mutations in their genotypes could affect the physiological and structural function of the sperm following cryopreservation. Interestingly, no association of either EML5 or PRICKLE1 genotypes with post-thaw motility was observed this study. This is surprising due to the critical role of MT or MAPs in the motility apparatus of spermatozoa. A possible explanation could be that the microscope assessment used in this study gives only one measured parameter of sperm motility, while the CASA system provides a variety of different motility parameters [29,50], which might be useful in association studies.
Of the gene associated with t-aminoacylation is SARS, polymorphism (rs344846507) in this gene was associated with mitochondrial function and acrosome integrity of FT spermatozoa. It is noteworthy that SARS is a nuclei acid binding gene related to tRNA-aminoacylation biosynthesis (ssc00970) and mitochondrial disorders. Earlier observation showed that increased abundance of SARS proteins in spermatozoa resulted in poor blastocyst development [51]. The finding of this study raises the possibility that the SARS polymorphism could have an important effect on male fertility, and this necessitates further studies. The SCLT1 gene is also of importance, in which the frequency distributions of either CT or TT genotype at rs337913978 are associated with reduced post-thaw semen quality. Enrichment analysis showed that the SCLT1 gene is associated with 9 GO terms, such as “microtubule organizing center” and “cilium”, suggesting the critical role of the gene in sperm function. It has been reported that the SCLT1 gene polymorphism causes dysfunction in axoneme assembly [52], which could compromise sperm motility. Presently, no significant associations between the SCLT1 gene and post-thaw semen quality have been reported. Another polymorphism (rs331568674) of interest was detected in the ROBO1 gene. Our findings show that the ROBO1 gene is enriched in the axon guidance pathway (ssc04360) and 12 GO terms involved in biological functions, such as “cell motility”, “signal transduction”, and “protein modification process”. In a recent study, it has been reported that the axon guidance pathway has a potential biological function in boar fertility [53], and its pathway has a significant role in DNA methylation alterations in sperm cell development [54]. In the current study, it seems that the TT genotype of the ROBO1 gene polymorphism might increase the sperm susceptibility to cryo-induced DNA fragmentation, as shown in the poor freezability ejaculates. It should be stressed that the homologous genotype was not detected in SCLT1 (CT-TT), MAP3K20 (GA-AA), MS4A2 (GC-CC), and ROBO1 (AT-TT) genes. Previous studies reported that the absence of a homologous genotype might indicate its removal through artificial selection and breeding [32,33]. Artificial selection and breeding is constantly practiced in the pig AI industry, and it is possible that such phenomenon might explain the absence of the homologous genotype in the above-mentioned gene polymorphisms. Furthermore, we are unable to explain the presence of low frequency of alleles in the homologous genotype of genes, such as CYP7B1 and RAB3C. Besides intensive selective breeding, it has been suggested that some genotypes might be eliminated by adaption of individuals to environmental stress [33]. However, it is unclear at this point whether the gene polymorphisms that were associated with post-thaw semen quality present a direct functional effect or if they are in linkage disequilibrium with other functional SNPs. Previous studies have reported that the effect on genotype or SNP variant might be influenced by the other SNPs, suggesting the interactions of multiple SNPs on the semen quality traits [32,33]. Accordingly, polymorphisms in boar spermatozoa could be used as markers associated with causative mutations within the gene [12,26].
In this study, the variant calling procedure detects polymorphisms associated with fertilization-related traits based on their relevance in sperm motility, metabolism, reproduction, and embryo development. We described here the associations of sperm-related genetic polymorphisms with post-thaw quality of boar semen. Candidate genes with significant SNP associations, including APPL1, PLBD1, FBXO16, EML5, RAB3C, OXSR1, PRICKLE1, and MAP3K20 genes, are promising markers for post-thaw semen quality, and that they might be relevant for future improvement in the selection procedure of boars for cryopreservation. We suggest that polymorphisms in genes expressed in spermatozoa could be considered as factors associated with post-thaw semen quality. However, further well-designed studies, with a larger animal population, are required to investigate the effect of gene polymorphisms in boar spermatozoa on semen freezability.

4. Materials and Methods

4.1. Animals and Semen Collections

A total of forty PLW boars were used in this study. The six boars used for RNA-Seq were stationed at the Cryopreservation laboratory, Faculty of Animal Bioengineering, University of Warmia and Mazury in Olsztyn [8]. For the association studies, most of the ejaculates were collected from boars, which were stationed at the Cryopreservation laboratory, and a total of 51 ejaculates were collected from boars stationed at three AI centers [5]. A total of approximately 296 ejaculates, at least three ejaculates per boar, were collected from the 40 boars (during the autumn–winter period), using the gloved-hand technique. The animals were fed with a commercial porcine ration and were kept in individual pens throughout the experimental period. Water was available ad libitum. Only sperm samples with more than 70% motility and 85% normal morphology were used for cryopreservation. Animal experiments were carried out in accordance with the guidelines set out by the Local Ethics Committee, Olsztyn (Poland). Experiments on boars (semen collection procedures) do not require the approval of the Local Ethics Committee (Olsztyn) from 15/01/2015.

4.2. Cryopreservation Procedure

Semen was frozen according to a cryopreservation protocol using lyophilized lipoprotein fractions of ostrich egg yolk (LPFo), as described in previous studies [39,50,55]. The LPFo-extended semen was cooled to 5 °C for 2h and diluted (2:1) with another freezing extender (89.5 mL lactose-LPFo extender, 9 mL glycerol, and 1.5 mL Orvus Es Paste). All samples (500 × 106 spermatozoa/mL) were frozen in a programmable computer freezer (Ice Cube 1810, SY-LAB, Austria), using an appropriate cooling rate [55], prior to storage in liquid nitrogen. Frozen samples were thawed in a water bath for 60 sec at 50 °C for post-thaw sperm analysis (motility and membrane integrity). Following post-thaw, the samples (50 × 106 spermatozoa/mL) were held in a water bath for 10 min at 37 °C prior to semen quality assessment.

4.3. Semen Quality Assessment

Quality assessment was performed in fresh, pre-freeze, and post-thaw semen.

4.3.1. Sperm Motility

The percentages of motile spermatozoa were assessed by the same technician throughout the study. Briefly, semen samples (6 μL) were placed on pre-warmed slide and assessed under a light microscope at 200 × magnification (Olympus BX 40, Tokyo, Japan) equipped with an attached heated stage (37 °C). Sperm motility was evaluated randomly in at least five fields per sample [5].

4.3.2. Mitochondrial Membrane Potential (MMP)

Sperm MMP was assessed with the fluorescent lipophilic cation JC-1 and propidium iodide (PI) fluorescent dyes [5,56].

4.3.3. Plasma Membrane Integrity (PMI)

Sperm PMI was assessed with the SYBR-14 and PI fluorescent probes, using the Live/Dead Sperm Viability Kit (Molecular Probes, Eugene, OR, USA) [57].

4.3.4. Normal Apical Ridge (NAR) Acrosome Integrity

A staining protocol was used to assess the sperm normal apical ridge (NAR) acrosome integrity [34,50].

4.3.5. DNA Fragmentation

The procedure used to assess sperm DNA fragmentation has been described in previous studies [39,55]. Briefly, agarose-embedded sperm samples on microscopic slides were stained with ethidium bromide and were assessed at 400× magnification under a fluorescence microscope (Olympus BX 41, Tokyo, Japan).

4.3.6. Lipid Peroxidation (LPO)

Sperm LPO was determined spectrophotometrically by malondialdehyde (MDA) production [5]. The LPO was defined as the production of nM MDA by 1 × 108 spermatozoa following 1 h incubation at 37 °C (nM MDA/108 spz/h).

4.4. Genomic DNA Isolation

Genomic DNA was isolated from washed sperm cells (200 × 106 spermatozoa/mL) of 40 boars (21 boars with GSF and 19 boars with PSF) using the protocol of the Sherlock AX Purification kit (A&A Biotechnology, Gdynia, Poland), according to the manufacturer’s instructions. Lysing solution containing proteinase K and dithiothreitol was added to sperm pellets. The mixture was incubated for 60 min at 50 °C and was subjected to a column filtration and purification procedure. Genomic DNA was treated with a precipitation enhancer isopropanol, and the DNA pellets were washed 2× in 70% ethanol and air-dried (10 min) at room temperature. The isolated DNA samples were dissolved in Tris-EDTA buffer and stored at −20 °C, until further analysis.

4.5. Variant Calling and Filtering Analysis

We used the RNA-Seq datasets that were submitted to the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA) database (Bioproject: PRJNA415904; accession number SRP121647) (2018/05/27). The datasets represent RNA-Seq from three PLW boars each of the GSF and PSF groups [8].
For the RNA-Seq, the clean reads were mapped to the Sus scrofa reference genome from Ensembl (genome-build Scrofa 11.1.91) using the two-pass mapping strategy in the spliced transcripts alignment to a reference (STAR) software [58]. The BAM alignment files were processed, and single nucleotide variants (SNVs) were called by the Picard tool (http://broadinstitute.github.io/picard) and GATK [59]. The data were filtered based on read depth greater than 10 (≥10) and quality score greater than 90 (≥ 90), as described in a previous study [23]. Comparison of frequency variations of SNVs was performed using the rMATS-DVR software [60]. Each potentially heterozygous genomic position was statistically examined using the False Discovery Rate (FDR< 0.05) to identify variations in the alternate allele frequencies between boars of the freezability groups. The parameters used to identify each SNV in the freezability groups were as follows: gene position, position of Sus scrofa chromosome (SSC), accession number in dbSNP database (version 11.91), differences in genetic coverage of both reference and alternate allele counts, allele fraction differences, number of counts for the reference and allele fractions, and the FDR value. Validation of SNPs was performed by matching putative polymorphisms to known pig dbSNP entries using standalone BLAT v.36 [61]. The SNPs that were not present in the dbSNP were considered as unknown (novel predicted SNP). As an additional variant filtering procedure, we retained SNPs that were on the coding sequences (CDS), 3′-untranslated regions (3′-UTRs), 5′untranslated regions (5′-UTRs), putative promotor regions (approximately 200 bp from the beginning of the transcript), and non-coding sequences (ncRNAs). After multi-filtering analysis, the identified variants were annotated by the SnpEff software v.4.1 [62] and Variant Effect Predictor (VEP) Ensembl [63] to retrieve the significant SNPs (Supplementary Table S1). We searched the literature to identify the functions of the CGs in the reproductive processes, sperm physiology or QTLs, according to the Animal Genome pigQTLdb. The CGs containing variants were selected according to their polymorphism coverage in the dataset (Supplementary Table S1) or their role in reproductive processes, sperm functions, or reproduction traits based on QTLs. Thus, by considering these criteria, forty CGs with polymorphisms (Supplementary Table S4) were selected for validation using the KASP genotyping assay (LGC Genomics Ltd., Trident Industrial Estate, Hoddesdon, Hertfordshire, UK).

4.6. KEGG Pathways and GO Enrichment Analysis

The Kobas functional annotation tool (v3.0) was used to identify in the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways statistical enrichment (p < 0.05) of the CG genes with polymorphisms [64]. Blast2Go Pro software, v.5.2.5 [65] was used to perform functional annotations, according to gene ontology (GO) categories (biological process, molecular function, and cellular components). Sus scrofa Ensembl database was downloaded from Ensembl BioMart Martview application to perform blastx, Blast2GO mapping, and the GO enrichment analysis. The GO significance levels (p < 0.05) were computed for multiple testing in the Blast2GO software program [65].

4.7. SNP KASP Genotyping Assay

The accuracy of the genotype calls was validated using the KASP genotyping assay (LGC Genomics Ltd., Trident Industrial Estate, Hoddesdon, Hertfordshire, UK). Forty SNPs were genotyped using 50 up and downstream flanking regions (Supplementary Table S5). Genomic DNA isolated from spermatozoa of forty PLW boars was shipped to the LGC Genomics Lab (http://www.lgcgroup.com) to perform the SNP genotyping assay based on KBioscience’s Kompetitive allele-specific PCR amplification. Primers designed for KASP™ genotyping procedure were performed by the LGC Genomics laboratory using the KBioscience PrimerPicker software [22]. The data were analyzed with the Kluster-caller software and SNPViewer (LGC Genomics, Ltd.; http://results.lgcgenomics.com/software/snpviewer,) to identify SNP genotypes.

4.8. Statistical Analysis

Data were statistically analyzed with the IBM SPSS Statistics 25 software package (IBM SPSS Statistics for Windows, version 25.0, IBM Corp., Armonk, NY, USA) and Statistica software package, version 12.5 (StatSoft Incorporation, Tulsa, OK, USA). The normality of the data distribution was analyzed with ANOVA assumption (Shapiro Wilk test), and the Levene’s test was used to check for homogeneity of variance. The general linear modeling (GLM) procedure was used for ANOVA analysis. The effect of boar on fresh, pre-freeze, and post-thaw semen quality was analyzed using Model (1).
Yij = µ +Bi + eij (Model 1)
where Y is the measured semen quality traits; µ is the overall mean; Bi is the fixed effect of boar; eij is the random residual effect.
The effect of SNP (G) on post-thaw semen quality was analyzed using Model (2).
Yij =µ +Gi +eij (Model 2)
where Y is the measured semen quality traits; µ is the overall mean; Gi is the fixed effect of the genotype; eij is the random residual effect.
In the boar population not all possible genotypes were detected for SCLT1, MAP3K20, MS4A2, and ROBO1. The association analysis of these gene polymorphisms with post-thaw semen quality was performed using Model (3).
Yijk =µ + Fi + Gj + eijk (Model 3)
where Y is the measured semen quality traits; µ is the overall mean; Fi is the fixed effect of freezability; Gj is the fixed effect of the genotype; eijk is the random residual effect.
Pairwise comparisons were analyzed with an independent two-tailed T-test, while multiple comparisons were performed with Tukey’s honest significant difference (HSD) test as appropriate. Descriptive variables are presented as the mean ± SEM. Significant differences were considered at p < 0.05.
The GENEPOP v.4.7.2 software package [66] (https://kimura.univ-montp2.fr/~rousset/Genepop.htm) was used to a) calculate the allele and genotype frequencies (Fisher’s exact test; b) determine Hardy–Weinberg equilibrium (HWE) for all loci using both the Fisher’s exact test and chi-square test (χ2-test), and c) examine the population differences between the two freezability groups.

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/1422-0067/21/5/1902/s1. Supplementary Table S1. Identification of single-nucleotide variants (SNVs) in spermatozoa boars with good and poor semen freezability (GSF and PSF, respectively). Supplementary Table S2. KEGG pathway of enriched candidate genes (CGs) of boar spermatozoa. Supplementary Table S3. Gene Ontology (GO) of enriched terms associated with selected candidate genes (CGs) of boar spermatozoa. Table S4. Reproductive processes, sperm functions and quantitative trait loci (QTLs) for pig reproductive traits of selected candidate genes (CGs) with single nucleotide polymorphisms (SNPs). Supplementary Table S5. Single nucleotide polymorphism (SNP) markers and primer sequences used for KASP genotyping assay. References [67,68,69,70,71,72,73,74,75,76,77,78,79,80,81,82,83,84,85,86,87,88,89,90,91,92] are cited in the supplementary materials.

Author Contributions

Conceptualization, L.F., A.M., and P.B.; methodology, A.M., and L.F.; software, Ł.P., J.P.J., and L.F.; validation, A.M., L.F., and P.B..; formal analysis, A.M., P.B., Ł.P., J.P.J., and L.F.; investigation, A.M., and L.F.; data curation, A.M., P.B., and L.F.; writing—original draft preparation, A.M., and L.F.; writing—review and editing, L.F., A.M., P.B., Ł.P., and J.P.J.; visualization, A.M., P.B., and LF; supervision, L.F.; project administration, L.F; funding acquisition, L.F. All authors have agreed to the final version of the manuscript.

Funding

This study was supported by a research project from the National Science Centre, Poland (2015/19/B/NZ9/01333).

Acknowledgments

A.M. is a recipient of a scholarship from the Program Interdisciplinary Doctoral Studies in Bioeconomy (POWR.03.02.00-00-I034/16-00) financed by the European Social Funds.

Conflicts of Interest

None of the authors have any conflicts of interest to declare.

Availability of data and materials

Raw sequencing data in the fastq format were submitted to the National Center for Biotechnology Information (NCBI) SRA database (https://www.ncbi.nlm.nih.gov/Traces/sra), with accession number SRP121647.

References

  1. Didion, B.A.; Braun, G.D.; Duggan, M.V. Field fertility of frozen boar semen: A retrospective report comprising over 2600 AI services spanning a four year period. Anim. Reprod. Sci. 2013, 137, 189–196. [Google Scholar] [CrossRef]
  2. Yeste, M.; Joan, E.; Rodríguez-Gil, J.E.; Bonet, S. Artificial insemination frozen-thawed boar sperm. Mol. Reprod. Dev. 2017, 84, 802–813. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Roca, J.; Broekhuijse, M.L.; Parrilla, I.; Rodríguez-Martínez, H.; Martinez, E.A.; Bolarin, A. Boar differences in artificial insemination outcomes: Can they be minimized? Reprod. Domest. Anim. 2015, 50 (Suppl. 2), 48–55. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Holt, W.V.; Medrano, A.; Thurston, L.M.; Watson, P.F. The significance of cooling rates and animal variability for boar sperm cryopreservation: Insights from the cryomicroscope. Theriogenology 2005, 63, 370–382. [Google Scholar] [CrossRef] [PubMed]
  5. Fraser, L.; Strzeżek, J.; Kordan, W. Post-thaw sperm characteristics following long-term storage of boar semen in liquid nitrogen. Anim. Reprod. Sci. 2014, 147, 119–127. [Google Scholar] [CrossRef]
  6. Thurston, L.M.; Siggins, K.; Mileham, A.J.; Watson, P.F.; Holt, W.V. Identification of amplified restriction fragment length polymorphism markers linked to genes controlling boar sperm viability following cryopreservation. Biol. Reprod. 2002, 66, 545–554. [Google Scholar] [CrossRef]
  7. Van Son, M.; Tremoen, N.H.; Gaustad, A.H.; Myromslien, F.D.; Våge, D.I.; Stenseth, E.-B.; Zeremichael, T.T.; Grindflek, E. RNA sequencing reveals candidate genes and polymorphisms related to sperm DNA integrity in testis tissue from boars. BMC Vet. Res. 2007, 13, 362. [Google Scholar] [CrossRef] [Green Version]
  8. Fraser, L.; Brym, P.; Pareek, C.S.; Mogielnicka-Brzozowska, M.; Paukszto, Ł.; Jastrzębski, J.P.; Wasilewska-Sakowska, K.; Mańkowska, A.; Sobiech, P.; Żukowski, K. Transcriptome analysis of boar spermatozoa with different freezability using RNA-Seq. Theriogenology 2020, 142, 400–413. [Google Scholar] [CrossRef]
  9. Dai, D.-H.; Qazi, I.H.; Ran, M.-X.; Liang, K.; Zhang, Y.; Zhang, M.; Zhou, G.-B.; Angel, C.; Zeng, C.-J. Exploration of miRNA and mRNA profiles in fresh and frozen-thawed boar sperm by transcriptome and small RNA sequencing. Int. J. Mol. Sci. 2019, 20, 802. [Google Scholar] [CrossRef] [Green Version]
  10. Ran, M.X.; Zhou, Y.M.; Liang, K.; Wang, W.C.; Zhang, Y.; Zhang, M.; Yang, J.D.; Zhou, G.B.; Wu, K.; Wang, C.D.; et al. Comparative analysis of microRNA and mRNA profiles of sperm with different freeze tolerance capacities in boar (Sus scrofa) and giant panda (Ailuropoda melanoleuca). Biomolecules 2019, 9, 432. [Google Scholar] [CrossRef] [Green Version]
  11. Ropka-Molik, K.; Bereta, A.; Żukowski, K.; Tyra, M.; Piórkowska, K.; Żak, G.; Oczkowicz, M. Screening for candidate genes related with histological microstructure, meat quality and carcass characteristics in pig based on RNA-seq data. Asian Aust. J. Anim. Sci. 2018, 31, 1565–1574. [Google Scholar] [CrossRef] [PubMed]
  12. Wimmers, K.; Lin, C.L.; Tholen, E.; Jennen, G.J.; Schellander, K.; Ponsuksili, S. Polymorphisms in candidate genes as markers for sperm quality and boar fertility. Anim. Genet. 2005, 36, 152–155. [Google Scholar] [CrossRef] [PubMed]
  13. Diniz, D.B.; Lopes, M.S.; Broekhuijse, M.L.W.J.; Lopes, P.S.; Harlizius, B.; Guimaraes, S.E.F.; Duijvesteijn, N.; Knol, E.F.; Silva, F.F. A genome-wide association study reveals a novel candidate gene for sperm motility in pigs. Anim. Reprod. Sci. 2014, 151, 201–207. [Google Scholar] [CrossRef] [PubMed]
  14. Mucha, A.; Piórkowska, K.; Ropka-Molik, K.; Szyndler-Nędza, M. New polymorphic changes in the WNT7A gene and their effect on reproductive traits in pigs. Ann. Anim. Sci. 2018, 18, 375–385. [Google Scholar] [CrossRef] [Green Version]
  15. Marques, D.B.D.; Bastiaansen, J.W.M.; Broekhuijse, M.L.W.J.; Lopes, M.S.; Knol, E.F.; Harlizius, B.; Guimarães, S.E.F.; Silva, F.F.; Lopes, P.S. Weighted single-step GWAS and gene network analysis reveal new candidate genes for semen traits in pigs. Genet. Sel. Evol. 2018, 50, 40. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Xing, Y.; Ren, J.; Ren, D.; Guo, Y.; Wu, Y.; Yang, G.; Mao, H.; Brenig, B.; Huang, L. A whole genome scanning for quantitative trait loci on traits related to sperm quality and ejaculation in pigs. Anim. Reprod. Sci. 2008, 114, 210. [Google Scholar] [CrossRef]
  17. Matoulkova, E.; Michalova, E.; Borivoj Vojtesek, B.; Hrstka, R. The role of the 3′ untranslated region in post-transcriptional regulation of protein expression in mammalian cells. RNA Biol. 2012, 9, 563–576. [Google Scholar] [CrossRef] [Green Version]
  18. Mayr, C. Evolution and biological roles of alternative 3′UTRs. Trends Cell Biol. 2016, 26, 227–237. [Google Scholar] [CrossRef] [Green Version]
  19. Liu, J.; Sun, Y.; Yang, C.; Zhang, Y.; Jiang, Q.; Huang, J.; Ju, Z.; Wang, X.; Zhong, J.; Wang, C. Functional SNPs of INCENP affect semen quality by alternative splicing mode and binding affinity with the target bta-mir-378 in Chinese Holstein bulls. PLoS ONE 2016, 11, e0162730. [Google Scholar] [CrossRef]
  20. Liu, L.; Yu, S.; Chen, R.; Lv, X.; Pan, C. A novel synonymous SNP (A47A) of the TMEM95 gene is significant;y associated with the reproductive tracts related to testis in male piglets. Arch. Anim. Breed. 2017, 60, 235–241. [Google Scholar] [CrossRef] [Green Version]
  21. Zhao, D.; Wu, W.; Xu, B.; Niu, X.; Cui, H.; Zhang, Y.; Wang, Z.; Wang, X. Variants in the SRD5A2 gene are associated with quality of semen. Mol. Med. Rep. 2012, 6, 639–644. [Google Scholar] [CrossRef] [PubMed]
  22. KBiosciences. PrimerPicker Lite for KASPar v.0.26; KBiosciences Ltd.: Hoddesdon, UK, 2009. [Google Scholar]
  23. Pareek, C.S.; Smoczyński, R.; Kadarmideen, H.N.; Dziuba, P.; Blaszczyk, P.; Sikora, M.; Walendzik, P.; Grzybowski, T.; Pierzchała, M.; Horbańczuk, J.; et al. Single nucleotide polymorphism discovery in bovine pituitary gland using RNA-Seq technology. PLoS ONE 2016, 11, e0161370. [Google Scholar] [CrossRef] [PubMed]
  24. Fardo, D.W.; Becker, K.D.; Bertram, L.; Tanzi, R.E.; Lange, C. Recovering unused information in gemone-wide association studies: The benefit of analyzing SNPs out of Hardy-Weinberg equilibrium. Europ. J. Hum. Genet. 2009, 17, 1676–1682. [Google Scholar] [CrossRef] [PubMed]
  25. Minelli, C.; Hompson, J.R.; Abrams, K.R.; Thakkinstian, A.; Attia, J. How should we use information about HWE in the meta-analyses of genetic association studeis? Int. J. Epidemiol. 2007, 37, 136–146. [Google Scholar] [CrossRef] [PubMed]
  26. Lin, C.L.; Ponsuksili, S.; Tholen, E.; Jennen, D.G.J.; Schellander, K.; Wimmers, K. Candidate gene markers for sperm quality and fertility of boars. Anim. Reprod. Sci. 2006, 349–363. [Google Scholar] [CrossRef]
  27. Rodríguez-Martínez, H. Laboratory semen assessment and predictions of fertility: Still Utopia. Reprod. Domest. Anim. 2003, 38, 312–318. [Google Scholar] [CrossRef] [PubMed]
  28. Gadea, J. Sperm factors related to in vitro and in vivo porcine fertility. Theriogenology 2005, 63, 431–444. [Google Scholar] [CrossRef] [PubMed]
  29. Yeste, M. Sperm cryopreservation update: Cryodamage, markers, and factors affecting the sperm freezability in pigs. Theriogenology 2016, 85, 47–64. [Google Scholar] [CrossRef] [PubMed]
  30. Daigneault, B.W.; McNamara, K.A.; Purdy, P.H.; Krisher, R.L.; Knox, R.V.; Rodriguez-Zas, S.; Miller, D.J. Enhanced fertility prediction of cryopreserved boar spermatozoa using novel sperm function assessment. Andrology 2015, 3, 558–568. [Google Scholar] [CrossRef] [Green Version]
  31. Gadea, J.; Sellés, E.; Marco, M. The predictive value of porcine seminal parameters on fertility outcome under commercial conditions. Reprod. Domest. Anim. 2004, 239, 303–308. [Google Scholar] [CrossRef] [Green Version]
  32. Dai, L.; Zhao, Z.; Zhao, R.; Xiao, S.; Jiang, H.; Yue, X.; Li, X.; Gao, Y.; Liu, J.; Zhang, J. Effects of novel single nucleotide polymorphisms of the FSH beta-subunit gene on semen quality and fertility in bulls. Anim. Reprod. Sci. 2009, 114, 14–22. [Google Scholar] [CrossRef] [PubMed]
  33. Liu, X.; Ju, Z.; Wang, L.; Zhang, Y.; Huang, J.; Li, Q.; Zhong, J.; An, L.; Wang, C. Six novel single-nucleotide polymorphisms in SPAG11 gene and their association with sperm quality traits in Chinese Holstein bull. Anim. Reprod. Sci. 2011, 129, 14–21. [Google Scholar] [CrossRef] [PubMed]
  34. Fraser, L.; Dziekońska, A.; Strzeżek, R.; Strzeżek, J. Dialysis of boar semen prior to freezing thawing: Its effects on post-thaw sperm characteristics. Theriogenology 2007, 67, 994–1003. [Google Scholar] [CrossRef] [PubMed]
  35. Kasimanickam, V.R.; Kasimanickam, R.K.; Kastelic, J.P.; Stevenson, J.S. Associations of adiponectin and fertility estimates in Holstein bulls. Theriogenology 2013, 79, 766–777. [Google Scholar] [CrossRef]
  36. Thomas, S.; Kratzsch, D.; Schaab, M.; Scholz, M.; Grunewald, S.; Thiery, J. Seminal plasma adipokine levels are correlated with functional characteristics of spermatozoa. Fertil. Steril. 2013, 99, 1256–1263. [Google Scholar] [CrossRef] [PubMed]
  37. Kadivar, A.; Heidari Khoei, H.; Hassanpour, H.; Golestanfar, A.; Ghanaei, H. Correlation of adiponectin mRNA abundance and its receptors with quantitative parameters of sperm motility in rams. Int. J. Fertil. Steril. 2016, 10, 127–135. [Google Scholar]
  38. Yung-Liang, L.; Sung-Sen, Y.; Shyi-Jou, C.; Yu-Chun, L.; Chin-Chen, C.; Hsin-Hui, H.; Chang, F.W.; Yu, M.H.; Lin, S.H.; Wu, G.J.; et al. OSR1 and SPAK cooperatively modulate Sertoli cell support of mouse spermatogenesis. Sci. Rep. 2016, 6, 37205. [Google Scholar]
  39. Fraser, L.; Strzeżek, J. Effect of different procedures of ejaculate collection, extenders and packages on DNA integrity of boar spermatozoa following freezing-thawing. Anim. Reprod. Sci. 2007, 99, 317–329. [Google Scholar] [CrossRef]
  40. Li, M.W.M.; Mruk, D.D.; Cheng, C.Y. Mitogen-activated protein kinases in male reproductive function. Trends Mol. Med. 2009, 15, 159–168. [Google Scholar] [CrossRef] [Green Version]
  41. Kaneko, T.; Toshimori, K.; Lida, H. Subcellular localization of MS4A13 isoform 2 in mouse spermatozoa. Reproduction 2017, 154, 843–857. [Google Scholar] [CrossRef]
  42. Schuster, A.; Tang, C.; Xie, Y.; Ortogero, N.; Yua, S.; Yan, W. SpermBase—A database for sperm-borne RNA contents. Biol. Reprod. 2016, 95, 1–12. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Neild, D.M.; Brouwers, J.F.H.M.; Colenbrander, B.; Agüero, A.; Gadella, B.M. Lipid peroxide formation in relation to membrane stability of fresh and frozen thawed stallion spermatozoa. Mol. Reprod. Dev. 2005, 72, 230–238. [Google Scholar] [CrossRef] [PubMed]
  44. Garde, J.; Roldan, E.R.S. rab3-Peptide stimulates exocytosis of the ram sperm acrosome via interaction with cyclic AMP and phospholipase A2 metabolites. FEBS Lett. 1996, 391, 263–268. [Google Scholar] [CrossRef] [Green Version]
  45. Sutovsky, P.; Moreno, R.; Ramalho-Santos, J.; Dominko, T.; Thompson, W.E.; Schatten, G.A. Putative, ubiquitin-dependent mechanism for the recognition and elimination of defective spermatozoa in the mammalian epididymis. J. Cell Sci. 2001, 114, 1665–1675. [Google Scholar] [PubMed]
  46. Suryawanshi, A.R.; Khan, S.A.; Gajbhiye, R.K.; Gurav, M.Y.; Khole, V.V. Differential proteomics leads to identification of domain-specific epididymal sperm proteins. J. Androl. 2011, 32, 240–259. [Google Scholar] [CrossRef]
  47. Torabi, F.; Bogle, O.A.; Estanyol, J.M.; Oliva, R.; Miller, D. Zona pellucida-binding protein 2 (ZPBP2) and several proteins containing BX7B motifs in human sperm may have hyaluronic acid binding or recognition properties. Mol. Human Reprod. 2017, 23, 803–816. [Google Scholar] [CrossRef]
  48. Morris, R.L.; Hoffman, M.P.; Obar, R.A.; McCafferty, S.S.; Gibbons, I.R.; Leone, A.D.; Cool, J.; Allgood, E.L.; Musante, A.M.; Judkins, K.M.; et al. Analysis of cytoskeletal and motility proteins in the sea urchin genome assembly. Dev. Biol. 2006, 300, 219–237. [Google Scholar] [CrossRef] [Green Version]
  49. Chen, H.; Cheng, C.Y. Planar cell polarity (PCP) proteins and spermatogenesis. Semin. Cell Dev. Biol. 2016, 59, 99–109. [Google Scholar] [CrossRef] [Green Version]
  50. Wasilewska, K.; Fraser, L. Boar variability in sperm cryo-tolerance after cooling of semen in different long term extenders at various temperatures. Anim. Reprod. Sci. 2017, 185, 161–173. [Google Scholar] [CrossRef]
  51. McReynolds, S.; Dzieciatkowska, M.; Stevens, J.; Hansen, K.C.; Schoolcraft, W.B.; Katz-Jaffe, M.G. Toward the identification of a subset of unexplained infertility: A sperm proteomic approach. Fertil. Steril. 2014, 102, 692–699. [Google Scholar] [CrossRef]
  52. Yuan, S.; Liu, Y.; Peng, H.; Tang, C.; Hennig, G.W.; Wang, Z.; Wang, L.; Yu, T.; Klukovich, R.; Zhang, Y.; et al. Motile cilia of the male reproductive system require miR-34/miR-449 for development and function to generate luminal turbulence. Proc. Natl. Acad. Sci. USA 2019, 116, 3584–3593. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Zhao, Y.; Gao, N.; Li, X.; El-Ashram, S.; Wang, Z.; Zhu, L.; Jiang, W.; Peng, X.; Zhang, C.; Chen, Y.; et al. Identifying candidate genes associated with sperm morphology abnormalities using weighted single-step GWAS in a Duroc boar population. Theriogenology 2020, 141, 9–15. [Google Scholar] [CrossRef] [PubMed]
  54. Beck, D.; Sadler-Riggleman, I.; Skinner, M.K. Generational comparisons (F1 versus F3) of vinclozolin induced epigenetic transgenerational inheritance of sperm differential DNA methylation regions (epimutations) using MeDIP-Seq. Environ. Epigenetics 2017, 3, dvx016. [Google Scholar] [CrossRef] [PubMed]
  55. Fraser, L.; Parda, A.; Filipowicz, K.; Strzeżek, J. Comparison of post-thaw DNA integrity of boar spermatozoa assessed with the neutral Comet assay and sperm-Sus Halomax test kit. Reprod. Domest. Anim. 2010, 45, e155–e160. [Google Scholar] [CrossRef] [PubMed]
  56. Dziekońska, A.; Fraser, L.; Strzeżek, J. Effect of different storage temperatures on the metabolic activity of spermatozoa following liquid storage of boar semen. J. Anim. Feed Sci. 2009, 18, 638–649. [Google Scholar] [CrossRef]
  57. Garner, D.L.; Johnson, L.A. Viability assessment of mammalian sperm using SYBR-14 and propidium iodide. Biol. Reprod. 1995, 53, 276–284. [Google Scholar] [CrossRef]
  58. Dobin, A.; Davis, C.A.; Schlesinger, F.; Drenkow, J.; Zaleski, C.; Jha, S.; Batut, P.; Chaisson, M.; Gingeras, T.R. STAR: Ultrafast universal RNA-seq aligner. Bioinformatics 2013, 29, 15–21. [Google Scholar] [CrossRef]
  59. McKenna, A.; Hanna, M.; Banks, E.; Sivachenko, A.; Cibulskis, K.; Kernytsky, A.; Garimella, K.; Altshuler, D.; Gabriel, S.; Daly, M.; et al. The Genome Analysis Toolkit: A MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010, 20, 1297–1303. [Google Scholar] [CrossRef] [Green Version]
  60. Wang, J.; Pan, Y.; Shen, S.; Lin, L.; Xing, Y. rMATS-DVR: rMATS discovery of differential variants in RNA. Bioinformatics 2017, 33, 2216–2217. [Google Scholar] [CrossRef]
  61. Kent, W.J.; Sugnet, C.W.; Furey, T.S.; Roskin, K.M.; Pringle, T.H.; Zahler, A.M.; Haussler, D. The human genome browser at UCSC. Genome Res. 2002, 12, 996–1006. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Cingolani, P.; Platts, A.; Wang, L.L.; Coon, M.; Nguyen, T.; Wang, L.; Land, S.J.; Lu, X.; Ruden, D.M. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNP in the genome Drosophilla melanogaster w1118; iso-2; iso-3. Fly 2012, 6, 80–92. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. McLaren, W.; Gil, L.; Hunt, S.E.; Riat, H.S.; Ritchie, G.R.; Thormann, A.; Flicek, P.; Cunningham, F. The Ensembl variant effect predictor. Genome Biol. 2016, 17, 122. [Google Scholar] [CrossRef] [Green Version]
  64. Mao, X.; Cai, T.; Olyarchuk, J.G.; Wei, L. Automated genome annotation and pathway identification using the KEGG Orthology (KO) as a controlled vocabulary. Bioinformatics 2005, 21, 3787–3793. [Google Scholar] [CrossRef]
  65. Götz, S.; García-Gómez, J.M.; Terol, J.; Williams, T.D.; Nagaraj, S.H.; Nueda, M.J.; Robles, M.; Talon, M.; Dopazo, J.; Conesa, A. High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res. 2008, 36, 3420–3435. [Google Scholar] [CrossRef] [PubMed]
  66. Rousset, F. geneopop’007: A complete re-implementation of the genepop software for Windows and Linux. Mol. Ecol. Resour. 2008, 8, 103–106. [Google Scholar] [CrossRef] [PubMed]
  67. Vitku, J.; Kolatorova, L.; Ricco, C.; Ferroud, C.; Hennebert, O.; Skodova, T.; Heracek, J.; Starka, L. The quantitation of 7 beta-hydroxy-epiandrosterone in the plasma and seminal plasma of men with different degrees of fertility. Physiol. Res. 2018, 67 (Suppl. 3), S511–S519. [Google Scholar] [CrossRef] [PubMed]
  68. Hilz, S.; Fogarty, E.A.; Modzelewska, A.J.; Cohen, P.E.; Grimson, A. Transcriptome profiling of the developing male germ line identifies the miR-29 family as a global regulator during meiosis. RNA Biol. 2017, 14, 219–235. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Choi, J.W.; Herr, D.R.; Noguchi, K.; Yung, Y.C.; Lee, C.W.; Mutoh, T.; Lin, M.E.; Teo, S.T.; Park, K.E.; Mosley, A.N.; et al. LPA receptors: Subtypes and biological actions. Annu. Rev. Pharmacol. Toxicol. 2010, 50, 157–186. [Google Scholar] [CrossRef] [Green Version]
  70. O’Flaherty, C. Redox regulation of mammalian sperm capacitation. Asian J. Androl. 2015, 17, 583–590. [Google Scholar] [CrossRef]
  71. Jamsai, D.; Bianco, D.M.; Smith, S.J.; Merriner, D.J.; Ly-Huynh, J.D.; Herlihy, A.; Niranjan, B.; Gibbs, G.M.; O’Bryan, M.K. Characterization of gametogenetin 1 (GGN1) and its potential role in male fertility through the interaction with the ion channel regulator, cysteine-rich secretory protein 2 (CRISP2) in the sperm tail. Reproduction 2008, 135, 751–759. [Google Scholar] [CrossRef] [Green Version]
  72. Naz, R.K.; Sellamuthu, R. Receptors in spermatozoa: Are they real? J. Androl. 2006, 27, 627–636. [Google Scholar] [CrossRef] [PubMed]
  73. Zhou, W.; Song, P. Molecular cloning of a novel gene ZAhi-1 and its expression analysis during zebrafish gametogenesis. Mol. Biol. Rep. 2006, 33, 111–116. [Google Scholar] [CrossRef] [PubMed]
  74. Lobo, V.; Parte, P. Membrane bound glucose regulated protein 78 interacts with alpha-2-macroglobulin to promote actin reorganization in sperm during epididymal maturation. Mol. Hum. Reprod. 2019, 25, 137–155. [Google Scholar] [CrossRef] [PubMed]
  75. Silva, J.V.; Korrodi-Gregorio, L.; Luers, G.; Cardosa, M.J.; Patricio, A.; Maia, N.; da E Silva Cruz, E.F.; Fardilha, M. Characterisation of several ankyrin repeat protein variant 2, a phosphoprotein phosphatase 1-interacting protein, in testis and spermatozoa. Reprod. Fertil. Dev. 2016, 7, 1009–1019. [Google Scholar] [CrossRef]
  76. Owa, M.; Uchihashi, T.; Yanagisaw, H.A.; Yamano, T.; Iguchi, H.; Fukuzawa, H.; Wakabayashi, K.I.; Ando, T.; Kikkawa, M. Inner lumen proteins stabilize doublet microtubules in cilia and flagella. Nat. Commun. 2019, 1143, 1–10. [Google Scholar] [CrossRef] [Green Version]
  77. Du, X.; Shen, J.; Kugan, N.; Furth, E.E.; Lombard, D.B.; Cheung, C.; Pak, S.; Luo, G.; Pignolo, R.J.; DePinho, R.A.; et al. Telomere shortening exposes functions for the mouse Werner and Bloom Syndrome genes. Mol. Cell. Biol. 2004, 24, 8437–8446. [Google Scholar] [CrossRef] [Green Version]
  78. Sun, Y.; Cheng, L.; Gu, Y.; Xin, Y.; Wu, B.; Zhou, S.; Guo, S.; Liu, Y.; Diao, H.; Shi, H.; et al. A human lectin microarray for sperm surface glycosylation analysis. Mol. Cell. Proteom. 2016, 15, 2839–2851. [Google Scholar] [CrossRef] [Green Version]
  79. Kropp, J.; Carrillo, J.A.; Namous, H.; Daniels, A.; Salih, S.M.; Song, J.; Khatib, H. Male fertility status is associated with DNA methylation signatures in sperm and transcriptomic profiles of bovine preimplantation embryos. BMC Genom. 2017, 18, 280. [Google Scholar] [CrossRef] [Green Version]
  80. Augustine, L.M.; Markelewicz, R.J.; Boekelheide, K.; Cherrington, N.J. Xenobiotic and endobiotic transporter mRNA expression in the blood–testis barrier. Drug Metab. Dispos. 2005, 33, 182–189. [Google Scholar] [CrossRef] [Green Version]
  81. Pegg, A.E. The function of spermine. IUBMB Life 2014, 66, 8–18. [Google Scholar] [CrossRef]
  82. Mikolcevic, P.; Sigl, R.; Rauch, V.; Hess, M.W.; Pfaller, K.; Barisic, M.; Pelliniemi, L.J.; Boesl, M.; Geley, S. Cyclin-dependent kinase 16/PCTAIRE kinase 1 is activated by cyclin Y and is essential for spermatogenesis. Mol. Cell. Biol. 2012, 32, 868–879. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  83. Smith, M.E.; Saraceno, G.E.; Capani, F.; Castilla, R. Long-chain acyl-CoA synthetase 4 is regulated by phosphorylation. Biochem. Biophys. Res. Commun. 2013, 430, 272–277. [Google Scholar]
  84. Moscatelli, N.; Lunetti, P.; Braccia, C.; Armirotti, A.; Pisanello, F.; De Vittorio, M.; Zara, V.; Ferramosca, A. Comparative proteomic analysis of proteins involved in bioenergetics pathways associated with human sperm motility. Int. J. Mol. Sci. 2019, 20, 3000. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  85. Spinaci, M.; Volpe, S.; Bernardini, C.; De Ambrogi, M.; Tamanini, C.; Seren, E.; Galeati, G. Immunolocalization of heat shock protein 70 (Hsp 70) in boar spermatozoa and its role during fertilization. Mol. Reprod. Dev. 2005, 72, 534–544. [Google Scholar] [CrossRef] [PubMed]
  86. Braas, K.M.; Harakall, S.A.; Ouafik, L.; Eipper, B.A.; May, V. Expression of peptidylglycine alpha-amidating monooxygenase: An in situ hybridization and immunocytochemical study. Endocrinology 1992, 130, 2778–2788. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  87. Zhuang, Q.; Feng, T.; Coleman, M.L. Modifying the maker: Oxygenases target ribosome biology. Translation 2015, 3, e1009331. [Google Scholar] [CrossRef] [Green Version]
  88. Bragado, M.J.; Gil, M.C.; Martin-Hidalgo, D.; Hurtado de Llera, A.; Bravo, N.; Moreno, A.D.; Garcia-Marin, L.L. Src family tyrosine kinase regulates acrosome reaction but not motility in porcine spermatozoa. Reproduction 2012, 144, 67–75. [Google Scholar] [CrossRef] [Green Version]
  89. Modarressi, M.H.; Behnam, B.; Cheng, M.; Taylor, K.E.; Wolfe, J.; van der Hoorn, F.A. Tsga10 Encodes a 65-kilodalton protein that is processed to the 27-kilodalton fibrous sheath protein. Biol. Reprod. 2004, 70, 608–615. [Google Scholar] [CrossRef] [Green Version]
  90. Vicens, A.; Luke, L.; Roldan, E.R. Proteins involved in motility and sperm-egg interaction evolve more rapidly in mouse spermatozoa. PLoS ONE 2014, 9, e91302. [Google Scholar] [CrossRef] [Green Version]
  91. Dickinson, R.E.; Duncan, W.C. The SLIT–ROBO pathway: A regulator of cell function with implications for the reproductive system. Reproduction 2010, 139, 697–704. [Google Scholar] [CrossRef] [Green Version]
  92. Demain, L.A.; Conway, G.S.; Newman, W.G. Genetics of mitochondrial dysfunction and infertility. Clin. Genet. 2017, 91, 199–207. [Google Scholar] [CrossRef] [PubMed]
Ijms 21 01902 g001 550
Figure 1. Gene ontology (GO) of enriched terms associated with candidate genes (CGs) in boar spermatozoa. A. Biological process categories; B. molecular function categories; C. cellular component categories. The values in the parenthesis indicate the number of inputs and the percentage of CGs enriched to the category.
Figure 1. Gene ontology (GO) of enriched terms associated with candidate genes (CGs) in boar spermatozoa. A. Biological process categories; B. molecular function categories; C. cellular component categories. The values in the parenthesis indicate the number of inputs and the percentage of CGs enriched to the category.
Ijms 21 01902 g001
Ijms 21 01902 g002 550
Figure 2. Genotyping cluster plots of A. APPL1, B. OXSR1, C. FBXO16, D. RAB3C, E. PLBD1, and F. SARS genes using the Kompetitive Allele Specific PCR (KASP) assay. Samples shown in pink were considered questionable genotypes as they were not assigned to the cluster. NTC, non-template control. Visualization of the genotype data was performed with the SNPViewer software (http://results.lgcgenomics.com/software/snpviewer).
Figure 2. Genotyping cluster plots of A. APPL1, B. OXSR1, C. FBXO16, D. RAB3C, E. PLBD1, and F. SARS genes using the Kompetitive Allele Specific PCR (KASP) assay. Samples shown in pink were considered questionable genotypes as they were not assigned to the cluster. NTC, non-template control. Visualization of the genotype data was performed with the SNPViewer software (http://results.lgcgenomics.com/software/snpviewer).
Ijms 21 01902 g002
Ijms 21 01902 g003 550
Figure 3. Effect of polymorphisms in A. APPL1, B. CYP7B1, C. RAB3C, D. OXSR1, E. FBXO16, and F. PLBD1 genes on post-thaw motility of boar spermatozoa. Data are presented as the mean ± SEM; a,b Differences between genotypes were significantly differed at p < 0.05.
Figure 3. Effect of polymorphisms in A. APPL1, B. CYP7B1, C. RAB3C, D. OXSR1, E. FBXO16, and F. PLBD1 genes on post-thaw motility of boar spermatozoa. Data are presented as the mean ± SEM; a,b Differences between genotypes were significantly differed at p < 0.05.
Ijms 21 01902 g003
Ijms 21 01902 g004 550
Figure 4. Effect of PLBD1 gene polymorphism on post-thaw A. MMP, B. PMI, C. NAR acrosome integrity, and D. LPO of boar spermatozoa. Data are presented as the mean ± SEM. a,b Differences between genotypes were significantly differed at p < 0.05.
Figure 4. Effect of PLBD1 gene polymorphism on post-thaw A. MMP, B. PMI, C. NAR acrosome integrity, and D. LPO of boar spermatozoa. Data are presented as the mean ± SEM. a,b Differences between genotypes were significantly differed at p < 0.05.
Ijms 21 01902 g004
Ijms 21 01902 g005 550
Figure 5. Effect of polymorphisms in A–B. APPL1, C–E. EML5, and F. FBXO16 genes on post-thaw membrane integrity of boar spermatozoa. Data are presented as the mean ± SEM. a,b Differences between genotypes were significantly different at p < 0.05.
Figure 5. Effect of polymorphisms in A–B. APPL1, C–E. EML5, and F. FBXO16 genes on post-thaw membrane integrity of boar spermatozoa. Data are presented as the mean ± SEM. a,b Differences between genotypes were significantly different at p < 0.05.
Ijms 21 01902 g005
Ijms 21 01902 g006 550
Figure 6. Effect of polymorphisms in A–B. SARS, C–D. RAB3C, and E–F. PRICKLE1 genes on post-thaw membrane integrity of boar spermatozoa. Data are presented as the mean ± SEM. a,b Differences between genotypes were significantly different at p < 0.05.
Figure 6. Effect of polymorphisms in A–B. SARS, C–D. RAB3C, and E–F. PRICKLE1 genes on post-thaw membrane integrity of boar spermatozoa. Data are presented as the mean ± SEM. a,b Differences between genotypes were significantly different at p < 0.05.
Ijms 21 01902 g006
Ijms 21 01902 g007 550
Figure 7. Effect of freezability and gene polymorphisms in A–B. SCLTL1, C–D. MAP3K20, E. MS4A2, and F. ROBO1 genes on post-thaw quality of boar semen. Data are presented as the mean ± SEM. a,b Differences between genotypes were significantly different at p < 0.05; GSF—good semen freezability; PSF—poor semen freezability.
Figure 7. Effect of freezability and gene polymorphisms in A–B. SCLTL1, C–D. MAP3K20, E. MS4A2, and F. ROBO1 genes on post-thaw quality of boar semen. Data are presented as the mean ± SEM. a,b Differences between genotypes were significantly different at p < 0.05; GSF—good semen freezability; PSF—poor semen freezability.
Ijms 21 01902 g007
Table
Table 1. Fresh, pre-freeze, and post-thaw quality of boar semen. Data are presented as the mean ± SEM of 21 and 19 boars with good and poor semen freezability (GSF and PSF, respectively).
Table 1. Fresh, pre-freeze, and post-thaw quality of boar semen. Data are presented as the mean ± SEM of 21 and 19 boars with good and poor semen freezability (GSF and PSF, respectively).
Sperm ParametersFresh,
Pre-Freeze Semen		Post-Thaw SemenTotal
GSFPSFGSF + PSF
Motility (%)75.78 ± 0.2842.96 ± 1.63 a26.85 ± 1.28 b35.31 ± 1.65
MMP (%)87.83 ± 2.6448.82 ± 1.04 a31.55 ± 1.55 b40.62 ± 1.64
PMI (%)87.43 ± 1.2249.91 ± 0.84 a34.55 ± 1.62 b42.64 ± 1.51
NAR acrosome integrity (%)91.62 ± 1.0753.93 ± 0.60 a40.85 ± 0.86 b47.72 ± 1.38
DNA fragmentation (%)2.97 ± 0.137.85 ± 0.45 a10.56 ± 0.83 b9.16 ± 0.50
LPO (nMol MDA/109spz/h)21.34 ± 0.6827.37 ± 0.47 a31.55 ± 0.48 b29.36 ± 0.45
Values with different letters (a and b) are significantly differed (p < 0.05); MMP—mitochondrial membrane potential; PMI—plasma membrane integrity; NAR—normal apical ridge; LPO—lipid peroxidation; MDA—malondialdehyde.
Table
Table 2. Analysis of filtered and trimmed dbSNPs in boar spermatozoa using variant calling of RNA-Seq datasets.
Table 2. Analysis of filtered and trimmed dbSNPs in boar spermatozoa using variant calling of RNA-Seq datasets.
Sample
ID		Raw SNV≥10% Reads ≥90%% Raw Reads≥10 Reads
= 100%		Raw dbSNPs≥10% Reads ≥90%% Raw Reads≥10 Reads
= 100%
G01941,35580,7148.6%66,625632,89368,57710.8%57,182
G09873,32778,0848.9%73,323583,61066,13111.3%62,796
G17941,859144,41515.3%135,223618,586122,91719.9%116,302
P30866,451100,80611.6%94,399575,54984,14814.6%79,695
P38890,984110,57412.4%103,347591,52794,16015.9%88,990
P39898,293110,75312.3%103,664594,40294,10715.8%89,169
dbSNP—data base single nucleotide polymorphism; SNV—single nucleotide variant.
Table
Table 3. Putative SNPs detected in boar spermatozoa.
Table 3. Putative SNPs detected in boar spermatozoa.
SNP EffectSNP CountPercent (%)
3 prime UTR variant91967.031
5 prime UTR variant1208.752
Synonymous variant1007.294
Missense variant745.398
Non coding transcript exon variant574.158
Splice region variant40.292
Intron variant, non-coding transcript variant30.219
Stop gained20.145
Others40.292
Unknown886.419
Table
Table 4. Genic differentiation, frequencies of single nucleotide polymorphism (SNP) alleles, and Hardy–Weinberg equilibrium (HWE) with respective p-value, chi-square (χ2), and probability. GSF—good semen freezability; PSF—poor semen freezability.
Table 4. Genic differentiation, frequencies of single nucleotide polymorphism (SNP) alleles, and Hardy–Weinberg equilibrium (HWE) with respective p-value, chi-square (χ2), and probability. GSF—good semen freezability; PSF—poor semen freezability.
Table 4A
SNP IDLocusFreezability
Groups		Allele CountsTotalAllele Frequenciesp-valueHWE
AGAGp-valueχ2Probability
APPL1GSF1923420.4520.5481.000001.00005.86590.2094
unknownPSF1519340.4410.5590.0532
Total3442760.4470.553
CYP7B1GSF535400.1250.8750.711800.01059.10440.0585
rs338842672PSF335380.0790.9211.0000
Total870780.1030.897
EML5GSF1230420.2860.7140.062210.12836.16400.1872
rs345056502PSF1717340.5000.5000.3575
Total2947760.3820.618
LPAR1GSF2614400.6500.3500.628450.63472.02280.7316
unknownPSF2410340.7060.2940.5731
Total5024740.6760.324
RAB3CGSF1131420.2620.7380.100381.00007.57320.1085
rs81210636PSF1721380.4470.5530.0227
Total3852800.3500.650
SARSGSF1921400.4750.5250.501101.00008.84840.0650
rs344846507PSF1523380.3950.6050.0120
Total3444780.4360.564
TXNIPGSF1230420.2860.7140.106310.130114.09700.0070
rs340075321PSF1820380.4740.5260.0067
Total3050800.3750.625
Table 4B
SNP IDLocusFreezability
Groups		Allele CountsTotalAllele Frequenciesp-valueHWE
GAGAp-valueχ2Probability
OXSR1GSF1824420.4290.5710.481380.37758.41710.0774
rs339379734PSF1715320.5310.4690.0394
Total3539740.4730.527
A2MGSF834420.1900.8100.366580.000126.94320.0000
rs339026428PSF236380.0530.9470.0271
Total1070800.1250.875
ANKRD42GSF1626420.3810.6190.498790.641318.66500.0009
rs81210697PSF2018380.5260.4740.0015
Total3644800.4500.550
CCDC149GSF1725420.4050.5950.488040.36817.59010.1078
rs332902509PSF1717340.5000.5000.0611
Total3442760.4470.553
CFAP52GSF2812400.7000.3000.798200.61783.52450.4742
unknownPSF259340.7350.2650.2779
Total5321740.7160.284
COMMD2GSF1921400.4750.5250.645180.005019.75780.0006
rs318435440PSF1420340.4120.5880.0102
Total3341740.4460.554
FBXO16GSF1624400.4000.6000.166681.00000.94490.9180
rs341614458PSF2115360.5830.4170.6235
Total3739760.4870.513
MAP3K20GSF1127380.2890.7111.000000.25315.33600.2545
rs340643892PSF925340.2650.7350.2741
Total2052720.2780.722
WRNGSF929380.2370.7630.764701.00000.00001.0000
rs319208708PSF521260.1920.8081.0000
Total1450640.2120.788
Table 4C
SNP IDLocusFreezability
Groups		Allele CountsTotalAllele Frequenciesp-valueHWE
ACACp-valueχ2Probability
CLEC7AGSF1527420.3570.6431.000000.66473.21220.5230
rs325939188PSF1321340.3820.6180.3019
Total2848760.3680.632
EML6GSF1030400.2500.7500.605881.00004.43700.3501
rs322659685PSF1123340.3240.6760.3501
Total2153730.2840.716
CA CA
ABCB11GSF1624400.4000.6000.628970.000515.28340.0041
rs324930519PSF1123340.3240.6761.0000
Total2747740.6350.365
SMSGSF1226380.3160.6840.333520.001513.81620.0079
rs343122214PSF1620360.4440.5560.6577
Total2846740.3780.622
Table 4D
SNP IDLocusFreezability GroupAllele CountsAllele FrequenciesTotalp-valueHWE
CGCGp-valueχ2Probability
TMEM177GSF2117380.5530.4470.816660.000123.06610.0001
rs80954196PSF1818360.5000.5000.1750
Total3935740.5270.473
GC GC
MS4A2GSF1030400.2500.7500.323550.27849.03530.0602
rs339836492PSF1323360.3610.6390.0392
Total2353760.3030.697
Table 4E
SNP IDLocusFreezability GroupsAllele CountsTotalAllele Frequenciesp-valueHWE
CTCTp-valueχ2Probability
PRICKLE1GSF1527420.3570.6431.000001.00002.29180.6823
rs694366781PSF1220320.3750.6250.3179
Total2747740.3650.635
SCLT1GSF933420.2140.7860.774160.53221.26140.8679
rs337913978PSF630360.1670.8331.0000
Total1563780.1920.808
TC TC
ACSL4GSF537420.1190.8810.360371.00000.00001.0000
rs334625232PSF828360.2220.7781.0000
Total1365780.1670.833
ATP5F1AGSF2220420.5240.4760.066970.0000>43.1295<0.0000
rs328079913PSF2810380.7370.2630.0036
Total5030800.6250.375
HSPA13GSF1119300.3670.6330.577560.29603.79250.4348
rs335938037PSF822300.2670.7330.5072
Total1941600.3170.683
PAMGSF1129400.2750.7250.778660.25552.72920.6041
rs81217594PSF620260.2310.7691.0000
Total1749660.2580.742
RIOX2GSF929380.2370.7630.574281.00000.00001.0000
unknownPSF818260.3080.6921.0000
Total1747640.2660.734
SKAP2GSF1626420.3810.6190.216180.64222.19350.7002
rs336351767PSF826340.2350.7650.5200
Total2452760.3160.684
Table 4F
SNP IDLocusFreezability groupAllele countsTotalAllele frequenciesp-valueHWE
ATATp-valueχ2Probability
ROBO1GSF537420.1190.8810.745971.00000.00001.0000
rs331568674PSF529340.1470.8531.0000
Total1066760.1320.868
TA TA
PLBD1GSF1428420.3330.6670.166540.04876.04320.1959
rs321497623PSF1818360.5000.5001.0000
Total3246780.4100.590
Table 4G
SNP IDLocusFreezability groupAllele countsTotalAllele frequenciesp-valueHWE
GTGTp-valueχ2Probability
HARS2GSF832400.2000.8000.784010.54822.51520.6419
rs336003721PSF826340.2350.7650.5187
Total1658740.2160.784
Table
Table 5. ANOVA analysis showing the effects of freezability (F) and gene polymorphism (G) on post-thaw semen quality.
Table 5. ANOVA analysis showing the effects of freezability (F) and gene polymorphism (G) on post-thaw semen quality.
A. rs337913978 (SCLT1)Motility Plasma Membrane Integrity (PMI)
FGF × GFGF × G
p-valuep-valuep-valuep-valuep-valuep-value
<0.001>0.305<0.043<0.001>0.063<0.028
B. rs340643892 (MAP3K20)NAR Acrosome IntegrityLipid Peroxidation (LPO)
FGF × GFGF × G
p-valuep-valuep-valuep-valuep-valuep-value
<0.001>0.146<0.010<0.001<0.022>0.981
C. rs339836492 (MS4A2)Lipid Peroxidation (LPO)
FGF × G
p-valuep-valuep-value
<0.001<0.039>0.302
D. rs331568674 (ROBO1)DNA Fragmentation
FGF × G
p-valuep-valuep-value
<0.030>0.800<0.047
Significantly differed at p < 0.05; NAR—normal apical ridge.

Share and Cite

MDPI and ACS Style

Mańkowska, A.; Brym, P.; Paukszto, Ł.; Jastrzębski, J.P.; Fraser, L. Gene Polymorphisms in Boar Spermatozoa and Their Associations with Post-Thaw Semen Quality. Int. J. Mol. Sci. 2020, 21, 1902. https://doi.org/10.3390/ijms21051902

AMA Style

Mańkowska A, Brym P, Paukszto Ł, Jastrzębski JP, Fraser L. Gene Polymorphisms in Boar Spermatozoa and Their Associations with Post-Thaw Semen Quality. International Journal of Molecular Sciences. 2020; 21(5):1902. https://doi.org/10.3390/ijms21051902

Chicago/Turabian Style

Mańkowska, Anna, Paweł Brym, Łukasz Paukszto, Jan P. Jastrzębski, and Leyland Fraser. 2020. "Gene Polymorphisms in Boar Spermatozoa and Their Associations with Post-Thaw Semen Quality" International Journal of Molecular Sciences 21, no. 5: 1902. https://doi.org/10.3390/ijms21051902

APA Style

Mańkowska, A., Brym, P., Paukszto, Ł., Jastrzębski, J. P., & Fraser, L. (2020). Gene Polymorphisms in Boar Spermatozoa and Their Associations with Post-Thaw Semen Quality. International Journal of Molecular Sciences, 21(5), 1902. https://doi.org/10.3390/ijms21051902

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop