Next Article in Journal
Transglutaminase 3 Reduces the Severity of Psoriasis in Imiquimod-Treated Mouse Skin
Next Article in Special Issue
Proteomic Analyses of Human Sperm Cells: Understanding the Role of Proteins and Molecular Pathways Affecting Male Reproductive Health
Previous Article in Journal
Periodic Membrane Potential and Ca2+ Oscillations in T Cells Forming an Immune Synapse
Previous Article in Special Issue
Molecular Mechanisms Involved in the Impairment of Boar Sperm Motility by Peroxynitrite-Induced Nitrosative Stress
 
 
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/ijms21051572
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

The Transcriptome of Pig Spermatozoa, and Its Role in Fertility

by
Manuel Alvarez-Rodriguez
1,*,
Cristina Martinez
1,
Dominic Wright
2,
Isabel Barranco
3,
Jordi Roca
4 and
Heriberto Rodriguez-Martinez
1
1
Department of Biomedical & Clinical Sciences (BKV), BKH/Obstetrics & Gynaecology, Faculty of Medicine and Health Sciences, Linköping University, SE-58185 Linköping, Sweden
2
Department of Physics, Chemistry and Biology, Faculty of Science and Engineering, Linköping University, SE-58183 Linköping, Sweden
3
Biotechnology of Animal and Human Reproduction (TechnoSperm), Department of Biology, Institute of Food and Agricultural Technology, University of Girona, 17003 Girona, Spain
4
Department of Medicine and Animal Surgery, Faculty of Veterinary Medicine, International Campus for Higher Education and Research “Campus Mare Nostrum”, University of Murcia, 30100 Murcia, Spain
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2020, 21(5), 1572; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms21051572
Submission received: 6 February 2020 / Revised: 19 February 2020 / Accepted: 24 February 2020 / Published: 25 February 2020
(This article belongs to the Special Issue Advances in Molecular Regulation of Spermatozoa Function)
Browse Figures
Versions Notes

Abstract

:
In the study presented here we identified transcriptomic markers for fertility in the cargo of pig ejaculated spermatozoa using porcine-specific micro-arrays (GeneChip® miRNA 4.0 and GeneChip® Porcine Gene 1.0 ST). We report (i) the relative abundance of the ssc-miR-1285, miR-16, miR-4332, miR-92a, miR-671-5p, miR-4334-5p, miR-425-5p, miR-191, miR-92b-5p and miR-15b miRNAs, and (ii) the presence of 347 up-regulated and 174 down-regulated RNA transcripts in high-fertility breeding boars, based on differences of farrowing rate (FS) and litter size (LS), relative to low-fertility boars in the (Artificial Insemination) AI program. An overrepresentation analysis of the protein class (PANTHER) identified significant fold-increases for C-C chemokine binding (GO:0019957): CCR7, which activates B- and T-lymphocytes, 8-fold increase), XCR1 and CXCR4 (with ubiquitin as a natural ligand, 1.24-fold increase), cytokine receptor activity (GO:0005126): IL23R receptor of the IL23 protein, associated to JAK2 and STAT3, 3.4-fold increase), the TGF-receptor (PC00035) genes ACVR1C and ACVR2B (12-fold increase). Moreover, two micro-RNAs (miR-221 and mir-621) were down- and up-regulated, respectively, in high-fertility males. In conclusion, boars with different fertility performance possess a wide variety of differentially expressed RNA present in spermatozoa that would be attractive targets as non-invasive molecular markers for predicting fertility.

1. Introduction

Just minutes after natural mating or artificial insemination (AI) boar spermatozoa reach the sow oviduct [1], and increasing evidence indicates an interaction occurs between the semen and the lining epithelium of the female genital tract that initiates changes aiding the establishment and maintenance of pregnancy [2]. Transcriptomic changes in the pig uterus and oviduct in the presence of semen have been previously described [3,4]. Semen starts an inflammatory response by the female [5], and further induces changes in the expression of immune-related genes [4], apparently leading to the establishment of an immune tolerance status to paternal antigens by the female immune system [6]. This cross-talk between semen and the female genital tract seems initiated by signals originating from components of the seminal plasma (proteins, cytokines, extracellular vesicles) and spermatozoa (proteins, associated exosomes, RNAs) [4].
In mammals, spermiogenesis implies crucial changes of cell differentiation in round spermatids that leads, among other specific modifications, to the loss of most of their cytoplasm before becoming a testicular spermatozoon that soon moves from the seminiferous epithelium towards the epididymis. Further maturational changes occur in the epididymis, including compaction of the nuclear genome [7], that leads to the general belief that the ejaculated spermatozoon is unlikely to perform transcription and translation until its nucleus decondenses following fertilization [8]. However, the presence of RNA, maybe only representing remnants of spermatogenesis, has been described in boar spermatozoa [9], including non-coding RNA and miRNAs [10] (short non-coding RNA molecules that regulate gene expression via post-transcriptional gene silencing, chromatin-dependent gene silencing or RNA activation [11]). The presence of different miRNAs has been reported in both epidydimal and ejaculated mammalian spermatozoa [12,13]. In boars, miRNA profiling has been previously reported in testis, epididymis and ejaculated spermatozoa [14,15]. An increased expression of certain miRNAs in boar spermatozoa (let-7a, -7d, -7e, and miR-22) has been correlated with low sperm motility [16]. However, we still lack conclusive evidence of the origin [9], and the precise roles of these miRNAs in spermatozoa, including their impact on fertilization and the subsequent development of the zygote, ultimately affecting fertility.
In this relation, decreased levels of miR-10b, and miR-135b in asthenospermic human sperm samples compared to healthy controls [17] and other miRNAs has been associated with asthenozoospermic, teratozoospermic, and oligozoospermic populations [17,18]. However, to the best of our knowledge, similar studies of differential expression of miRNAs in relation to fertility performance are not available in other species including pigs, a species increasingly used as an animal-model for human fertility due to its similarities [19]. The porcine model allows us to not only analyse the farrowing rate and litter size but to correlate it with a wide range of sperm quality parameters or even differential RNA sperm load. Moreover, in pigs, AI has, in most countries, replaced natural mating and breeding boars are continuously renewed after a certain number of AIs. Thus finding suitable fertility markers in semen is an ongoing strategic goal, since it would reduce the numbers of less fertile boars being used and allow for enhanced selection of the most fertile candidates [20].
The same reasoning would also apply to markers based on gene expression, even when spermatozoa are supposed to be transcriptionally silent (with the exception of their mitochondria [21,22]). In fact, the total RNA cargo of boar spermatozoa present in the ejaculate and selected for robustness has been analyzed via RNA-seq [23], and indicates the presence of readily available transcripts.
The present study aimed therefore to determine the relative load (in terms of expression levels) of miRNAs and mRNAs in ejaculated boar spermatozoa, including the relation of any differential RNA cargo with fertility characteristics post AI. Microarray-based technology was used to study either (i) relative abundance of miRNAs on boar spermatozoa and (ii) total RNA differential expression between high fertile and low fertile boars. The hypothesis tested was that the available RNA-cargo in ejaculated spermatozoa includes specific transcripts whose level of expression is enriched in high-fertility breeding boars and thus could be valuable biomarkers for fertility.

2. Results

2.1. Experiment 1: Screening of Sperm miRNAs in the Sperm-rich Fraction (SRF) of the Boar Ejaculate

Total RNA was extracted from spermatozoa retrieved from the SRF of the ejaculate of healthy mature boars (n = 3), and then analyzed with the GeneChip™ miRNA 4.0 Array. Of the 326 pig (Sus scrofa)-specific miRNAs included in the array, 10 miRNAs had an expression of more than 10 average expression units (log 2 average expression = >1031 times expression from the basal-normalized signal), with these being miR-1285, miR-16, miR-4332, miR-92a, miR-671-5p, miR-4334-5p, miR-425-5p, miR-191, miR-92b-5p and miR-15b (Table 1). A total of 71 miRNAs were expressed between 5 and 10 average expression units and 245 miRNAs showed only <5 average expression units (Supplementary Table S1).

2.2. Experiment 2. Differential Cargo of RNAs in Ejaculated Spermatozoa from Breeding Boars with High- or Low-fertility (as Farrowing Rate and Litter Size) after Artificial Insemination (AI)

The Bioinformatic analysis of the normalized data collected from the gene expression analysis using the GeneChip® Porcine Gene 1.0 ST Array, using a Fold-change threshold value of ± 1 (p-value < 0.05) identified up to 521 differentially expressed genes (DEGs), among which 374 were up-regulated in high-fertility (n = 4) as compared to low-fertility (n = 3) males (Figure 1 and Supplementary Table S2).
When we performed a more stringent analysis of these DEGs (Table 2), 67 out of 347 up-regulated genes in high-fertility boars had a fold-change higher than 2, whereas 10 out of 174 down-regulated genes in high-fertility boars had a fold-change lower than −2.
Figure 2 depicts a PANTHER analysis of the biological functions of the differentially expressed genes (DEGs) (Figure 2A,B). The majority of the functions attributed for both up- and down-regulated genes were: cellular processes (GO:0009987), metabolic processes (GO:0008152) and biological regulation (GO:0065007). Interestingly, genes related to biological adhesion (GO:0022610; CDH10, CDSN, ITGB8, ANGPTL1, CTNNA3) and cell proliferation (GO:0008283; CCNJ) were only found in down-regulated genes in high-fertility boars whereas related genes to signaling (GO:0023052; FAM109A) and cellular component organization or biogenesis (GO:0071840; KRI1, ZNHIT6) were only identified in up-regulated genes in high-fertility boars.
In terms of molecular function (Figure 2C,D), both dysregulated groups of genes presented and overrepresentation of genes mainly involved in binding (GO:0005488) and catalytic activity (GO:0005488). Genes related to cargo receptor activity (GO:0038024; CXCL16) and translation regulator activity (GO:0045182; CPEB3) were, on the other hand, only found in up-regulated DEGs in high-fertility boars.
Figure 3 summarizes the analysis of the pathways in both groups of DEGs (up- (Figure 3A) and down-regulated (Figure 3B) genes) in high-fertility boars, detailed pie-chart include a large variety of functions common to up- and down-regulated genes in high-fertile boars.
An over-representation analysis (PANTHER) (Table 3) allowed us to restrict the number of pathways according to their representation significance (p < 0.05) relative to the expected genes present in the genome of Sus scrofa.
Finally, the analysis of a plausible interconnection between target genes of those miRNAs detected in experiment 1 (SRF-derived spermatozoa, with more than 10 average expression units) and the DEGs in the PORGENE-array in experiment 2 is presented in Table 4. Only a few genes, mostly up-regulated in high-fertility males, appear to be potentially shared. In this regard, results from DEGs in terms of miRNA dysregulation found only miR-615 up- and miR-221 to be down-regulated in high-fertility boars relative to the low-fertility boars (Supplementary Table S2).

3. Discussion

The microarray platform analysis used here and, in particular the multispecies miRNA array, have allowed us to define the relative abundance or cargo for large numbers of non-coding small RNAs present in spermatozoa retrieved from the SRF of boar ejaculates. An alternative would have been to use qPCR, but the latter is highly reliant on the number and selection of housekeeping genes [24]. The microarray platform is more robust and stable than a single qPCR because it is normalized against the total expression of the entire transcriptome that is measured in the array. The information provided by the microarray has been widely studied in terms both of the spermatogenesis process, in which the cytoplasm is reduced to a minimal portion, and the function of this restricted load plays on further endometrium receptivity and even on embryo development.
Only a few miRNAs out of the detected 326 pig-specific miRNAs were abundant. Of those in abundance, the miR-1285 has been related to sperm production, and considered inhibited by 17β-estradiol, the latter (estradiol) also promoting AMPK phosphorylation [25]. Moreover, this miR-1285 has been previously implicated in oxidative stress, with a response seen in human retinal pigment epithelial using a hydrogen peroxide challenge as an in vitro model for the induction of oxidative stress [26]. Another highly abundant miRNA, miR-15/miR-16 seems to suppress the TGF-β Signaling Pathway (a key superfamily cytokines that promotes apoptosis), mainly through the inhibition of the expression of endogenous Smad3 and ACVR2A proteins [27]. This is likely related to our previous finding which demonstrated the differential expression of immune regulatory genes in the porcine female genital tract in the presence of spermatozoa and seminal plasma [4]. The miRNA miR-4332 is a key regulator of lipid deposition found in pig muscle tissue [28] but has not previously been identified as a candidate for any known function directly related to reproduction. Increased levels of expression of miR-92 relate to decreases in the estrogen receptor β1 [29] and also suppresses proliferation and induces apoptosis by targeting EP4/Notch1 Axis, regulated by the NF-kB pathway [30]. Both of these miRNAs might be coupled to sperm survival and the triggering of certain responses in the female through miR-92 release from the spermatozoon. The miRNA miR-671-5p is an important inhibitor of cell proliferation and inductor of apoptosis by targeting URGCP (upregulator of cell proliferation) [31] and also the suppression of FGFR2 [32] (gene up-regulated in high-fertility boars in our results). Although this gene FGFR2 has been defined as non-essential for spermatogenesis and fertility in mice [33], it could play a different and more essential role in boar spermatozoa, with the potential of being used as a fertility biomarker. MiR-4334 is found in milk exosomes reduced lipopolysaccharide (LPS)-induced apoptosis via the p53 pathway; LPS is a bacterial endotoxin that induces inflammation [34]. Its presence could indicate a further regulation of this process in the female genital tract, induced via a sperm mediator. Mir-425-5p negatively regulates E2F6, a key regulator of the cell cycle and apoptosis [35] and, together with miR-4332, is linked to a negative regulation of the differentiation and proliferation of intramuscular pre-adipocytes [36]. Thus, the relative abundance of these two regulators should be studied to link the anti-apoptotic signaling from the spermatozoa to the female genital tract. Moreover, miR-425 and IL-23 have been correlated [37] and its receptor (IL-23R) was down-regulated in our comparison of high- vs. low-fertility boars. IL-23 play critical roles in the induction of Th17 differentiation and governs local inflammation [38], so it might be possible that a reduction in its receptor expression could translate in a reduction of the local immune reaction in the female to the presence of semen. Finally, among these overexpressed miRNAs, miR-191, in combination with miR-425, increases cell proliferation, survival and migration [39], potentially related to the travel of the spermatozoon through the entire female genital tract until the site of fertilization in the mid-ampulla of the oviduct or via repressive signals to the female genital tract during this sperm progression. It has been also described as key activator of the NF-κB signaling pathway [40] as well as a regulator of the estrogen receptor [41] and an inhibitor of the Wnt pathway [42]. Interestingly, some studies in humans have reported the reduced levels of miR-191 in the semen of infertile human patients [17,18].
The main focus of the second experiment was to identify, via microarray analysis, which transcripts were differentially expressed in ejaculated spermatozoa from breeding boars, highly selected for fertility after artificial insemination in commercial programs. The boars selected had fertility data based on more than 100 artificial inseminations (up to 540) per boar. The analyses of fertility resulted in two clear subgroups according to their deviation in farrowing rate (FR) and litter size (LS) from the average direct boar effect recorded by the boar population of the same genetic line: high-fertility boars (FR > 0.45 and a LS > 0.15); and low-fertility boars (FR < −0.12 and LS < −0.18). These fertility results show the true differences in the fertile potential of boars, since they were specifically recorded to show the direct effect of the boar on fertility, which was possible by using the statistical model developed by Broekhuijse et al. [43], and represent a very strong selection of individuals who would scape conventional screening of semen quality, which is continually done in breeding AI-programs. The microarrays used allow for the identification of more than 25,000 transcripts from total RNA, including isoforms, of Sus scrofa annotated genes and also of non-coding RNAs. The overall results showed differential gene expression with 347 up-regulated genes and 174 down-regulated genes in high-fertility boars relative to low-fertility boars. These findings are relevant based on the fact that the boars with different specific fertility included in the present study were breeding boars, highly selected for sperm quality, a pre-requisite for inclusion in programs of artificial insemination. All boars had normal semen quality, in terms of sperm numbers, sperm motility, sperm morphology and absence of inflammatory; contrasting with other species as human, or even other animal classes [44] that had been explored for the specific presence of miRNAs in individuals with clear differences in sperm quality and the influence of high- or low-sperm quality on progeny, often caused by specific spermatogenic impairment [45]. Of interest, these boars, deviating in true fertility but having normal semen quality, had differences in gene expression that ought to be relevant as future biomarkers. Two miRNAs: miR-615 appeared up-regulated while the miR-221 was down-regulated in high-fertility boars. The miR-615 is considered an inhibitor of apoptosis targeting EGFR [31], and it has been studied as being active in the presence of epidermal growth factor receptor (EGFR) in spermatozoa, as well as having a role in the transactivation by direct activation of the G-protein coupled receptor (GPCR), as well as being involved in sperm capacitation and in the acrosome reaction [46]. Thus, it could be a possible link between capacitation and the miR-615 expression. On the other hand, The miR-221 has been related to Wnt2, BDNF and CREB-related genes [47] and ties in to a previous result by our group describing the downregulation of CREB3L2 (p<0.05) in the uterus and utero-tubal junction in the porcine female genital tract in response to both mating and artificial insemination. Finally, it is also related to PI3K-Akt and the estrogen signaling pathway [4].
Overall (See Table 2 and Supplementary Table S2, DEGs), we have found several transcripts related to membrane channels that were up-regulated in high-fertility boars: KCNA3 (potassium channel, voltage gated shaker related subfamily A, member 3), KCNIP3 (Kv channel interacting protein 3, calsenilin), KCNH4; LOC100620256 (potassium channel, voltage gated related subfamily H, member 4; potassium voltage-gated channel subfamily H member 4-like), LOC102159765 (voltage-dependent P/Q-type calcium channel subunit alpha-1A-like), KCTD9 (potassium channel tetramerization domain containing 9), CNGA4 (cyclic nucleotide gated channel alpha 4). Among these channels, only KCNA3 has been characterized in rat testis as well as identified not only in the cytoplasm of primary spermatocytes and post-meiotic elongating spermatids, but also in the head and tails of epididymal spermatozoa [48]. Whether these channels are directly related to capacitation is still unknown, but the role of this gene should be examined in future studies in this context, in particular at the protein level. Finally, CATSPERG (catsper channel auxiliary subunit gamma) is up-regulated whereas the CATSPERB (catsper channel auxiliary subunit beta) was down-regulated in high-fertility boars. These two last transcripts are subunits of catsper, a channel involved in boar sperm motility during in vitro capacitation [49]. The CATSPER protein complex contains CATSPERB and CATSPERG, the latter associated with CATSPER1 [50] and the inhibition of CATSPER is linked to an increase on a reactive oxygen species (ROS) generation [51] and it is involved in key processes of mammalian fertilization by Calcium signaling [52].
Moreover, a representative number of transcripts included in the list of our DEGs were related to Zinc finger proteins (ZNF), one of the most abundant group of proteins with a wide variety of functions, including transcriptional regulation, ubiquitin-mediated protein degradation, signal transduction, actin targeting, DNA repair and cell migration [53]. The majority of the ZNF-related transcripts were up-regulated in high-fertility boars: LOC100739821 (zinc finger MIZ domain-containing protein 2), ZNRF4 (zinc and ring finger 4), LOC100512674 (zinc finger protein 367), LOC100511053; LOC100624021 (zinc finger protein PLAGL2), ZFAND3 (zinc finger, AN1-type domain 3), ZNF25 (zinc finger protein 25), ZDHHC7 (zinc finger, DHHC-type containing 7), ZNF268 (zinc finger protein 268), ZNHIT6 (zinc finger, HIT-type containing 6), ZNF568 (zinc finger protein 568), LOC102163983 (zinc finger protein 341). Among these, ZDHHC7 is directly involved in signaling via membrane-localized estrogen receptors (ER) [54] and the ER interacts directly with the Phosphatidylinositol-3-OH Kinase/Akt Pathway [55]. We conjecture that this up-regulation in the boar spermatozoa from high-fertility males could interact with the ER in the female genital tract. It is also relevant for the role of ZNF568 in a later event, being expressed in 3-5 week old embryos and potentially is related to embryogenesis [56]. In contrast, only a few of these genes were found to be down-regulated in high-fertility boars: ZNF283 (zinc finger protein 283), FEZF2; LOC100622205 (FEZ family zinc finger 2; fez family zinc finger protein 2), GLI1; LOC100621250 (GLI family zinc finger 1; zinc finger protein GLI1-like). Interestingly, down-regulation in high fertility males is in accordance with a previous study on GLI1 up-regulation in transgenic mice leading to a blockage of spermatogenesis in pachytene [57]. This gene is directly involved in the sonic hedgehog (Shh) signaling pathway, with important roles in pattern orientation in the developing embryo [58].
Transcription factors are also represented among our results of DEGs. Not only are they up-regulated in high-fertility boars: LOC100152461; LOC100627430; NFYA (nuclear transcription factor Y subunit alpha; nuclear transcription factor Y, alpha), TCF21 (transcription factor 21), LOC100524199; LOC100620736 (transcription factor 25), MBTPS2 (membrane-bound transcription factor peptidase, site 2), MTF2 (metal response element binding transcription factor 2), but some are also down-regulated in high-fertility boars: UBTFL1 (upstream binding transcription factor, RNA polymerase I-like 1).
The matrix metalloproteinases (ADAM) are involved in several processes such as cell adhesion and immune related signaling, including the activation of TNF-alfa and generation of the active forms of Epidermal Growth Factor Receptor (EGFR) [59]. In fact, the relevance of this ADAM gene has been previously described by our group, mainly in terms of biological adhesion but also in terms of pH regulation of the utero tubal junction in sows [60,61]. In the present study, we identify down-regulation in high-fertility boars of ADAM7 (ADAM metallopeptidase domain 7) and ADAM29 (ADAM metallopeptidase domain 29). ADAM7 in mouse spermatozoa is required for normal fertility [62], and here we its presence and the decreased abundancy in high-fertility spermaozoa. Both ADAMs have been linked to the increase of adhesion capacities of extracellular matrix proteins [63]. Whether this increase in adhesion is related to the progression of spermatozoa through the female genital tract is yet to be determined. Other relevant immune-related transcripts that appeared down-regulated in high-fertility boars were LOC102162067 (interleukin-6 receptor subunit alpha-like), IL23R (interleukin 23 receptor), IFN-DELTA-4; LOC100736862 (interferon-delta-4; interferon tau-11-like), IFNK (interferon, kappa), LOC100736882 (interferon lambda receptor 1-like) and IFIT3 (interferon-induced protein with tetratricopeptide repeats 3), however similar genes were also found to be up-regulated in high-fertility boars: IRF3 (interferon regulatory factor 3) and ISG20L2 (interferon stimulated exonuclease gene 20kDa-like 2).
In order to categorize the overall DEGs resulted from our study, we performed a biological function analysis (PHANTER) of the DEGs which resulted in a differential function depending on the fertility performance of the boars. Thus, biological adhesion-related genes (GO:0022610; CDH10, CDSN, ITGB8, ANGPTL1, CTNNA3) were found only to be down-regulated in high-fertility boars. In contrast, genes related to signaling (GO:0023052; FAM109A) and cellular component organization or biogenesis (GO:0071840; KRI1, ZNHIT6) were only identified as up-regulated in high-fertility boars. Cadherin-10 (CDH10), the type II classical cadherin, a protein binding molecule in terms of molecular function and is used as a biomarker in the prognosis of certain types of cancer [64], whilst CDSN is involved in skin immune-related disease [65]. Whether the down-regulation of both these genes in our study might be related to a decrease in the attachment of the spermatozoa to certain sections of the female tract remains to be studied. ITGB8 is involved in the endometrial receptivity through VAV-RAC1 signaling axis via FAK to facilitate the endometrial epithelial cell receptivity for the attachment of blastocyst [66] and, surprisingly, is less abundant in high-fertility boars. ANGPTL1, Angiopoietin-like proteins is involved in several processes from lipid and glucose metabolism to inflammation [67]. The last of this biological adhesion group is CTNNA3 (alpha-T-catenin), whose gene dosage compensation of CTNNA3 and p57KIP2 in the placenta shares a conserved regulatory mechanism that correlates with an early step in trophoblast determination [68]. Regarding cell proliferation, the cyclin J (CCNJ), down-regulated in high-fertile boars, is involved in cell cycle regulation, but without any apparent described function in spermatogenesis or fertility performance. Regarding up-regulated genes in high-fertility boars, ZRI1 permits DNA damage-induced apoptosis of cells in the germline by an unknown cell non-autonomous mechanism [69] and ZNHIT6 is required for the maintenance of correct protein folding under heat stress conditions [70].
If we categorize according to molecular function, one of the main findings was the cytoplasmic polyadenylation element-binding protein 3 (CPEB3), a translational regulator involved in long-term memory storage maintenance in synaptic connections in neurons [71], which has also been involved in EGFR axis regulation [72] and also in SUMOylation [73]. SUMOylation is a process that controls the activity of pre-existing proteins through a reversible post-transcriptional modification [74]. In addition, decidualization renders the endometrium transiently receptive to an implanting blastocyst although the underlying mechanisms remain incompletely understood, but CXCL16 appears to be an initiator for the molecular crosstalk that enhances decidualization of the endometrial stromal cells [75], despite this process as such has not been described in pig. Nevertheless, this fact is certainly in agreement with our results, observing up-regulation of this gene in high-fertility boars, but whether a “remodeling” or decidualization-like process is present in pig remains to be studied.
A further analysis of the overrepresentation analysis in terms of molecular function reveals an important component on chemokine binding (GO:0019957). Regulatory T (Treg) cells facilitate maternal immune tolerance of the hemi-allogeneic conceptus in early pregnancy. In this regard, seminal fluid induced expression of mRNA encoding the Treg chemokine CCL19 (MIP3beta), which acts through the CCR7 receptor to regulate Treg cell recruitment and retention in peripheral tissues [76]. Our results showed a decrease in abundancy of the CCR7 receptor in high-fertile boars. In addition, CXCR4, is the main G-coupled receptor of CXCL12, involved in several diseases in the female reproductive tract, mainly through inflammatory processes [77]. In agreement with this, the lower abundancy of this receptor on high-fertile boars, perhaps in relation to a reduced inflammatory response of the female genital tract in response to the entrance of spermatozoa with a lower load of this transcript.
We have also observed in our analysis of protein class overrepresentation, that the ACVR1C and ACVR2B were down-regulated in high-fertility males. Activin affects many aspects of cellular development, including those essential for reproductive fitness and, in fact, activin receptor transcript levels also change, with Acvr1 (encoding ALK2) and Acvr2b (ActRIIB) significantly higher and lower, respectively in the proliferation of mouse testis [78].
Finally, in an attempt to address a prospective functional approach to the high-expressed miRNAs, we contrasted the target genes of the most abundant miRNA (with a target match of more than 0.9. miRbase and miRDB [79]), and found several genes of interest that were mostly up-regulated in the high-fertility boars. In Table 4, as detailed in the results section, we found target genes differentially expressed in relation with miR-16, miR-92a and miR-15b. Interestingly, up-regulation of certain genes (UBQLNL, CPEB3, DESI1 and ATG9A) and down-regulation in high- relative with low-fertile boars were found both link to miR-16 and miR15b. In common with miR-92b we identified the transcript CPEB3, and the specific genes related solely to this miRNA: TCF21, PDZD8. Finally, in relation with miR-16, the gene BAG4 appeared up-regulated in the DEGs analysis performed with microarrays. The UBQLNL is a testis-specific gene, not indispensable due to a compensatory mechanism in mice, but involved in the elongation of spermatids and found up-regulated in high-fertile boar samples [80]. It is a member of the ubiquilin gene family, highly conserved in eukaryotes and with high expression levels in testis, and mainly in the “spermatid” stage [81]. In relation with the already mentioned SUMOylation process, our results of DEGs included the up-regulation in high-fertility boars of a SUMO1 activating enzyme subunit 1 (SAE1), essential for the maintenance of hematopoietic stem/progenitor cell in zebrafish [82] and also DESI1 (desumoylating isopeptidase 1), a second class od SUMO-related proteases [83]. The gene ATG9A (autophagy related 9A) is essential for autophagy, being trafficked to sites of autophagosome formation [84]. Finally, we identified a down-regulation of the LRRN3 gene in high-fertility boars, with LRRN3 being involved in Ras-MAPK signaling by facilitating internalization of EGF in clathrin-coated vesicles [85]. TCF21 (transcription factor 21) is a key inactivation factor of the PI3K/AKT signaling pathway and also of matrix metalloproteinases [86]. This is in agreement with previous proteomic results in bull spermatozoa in relation with fertility [87], with PDZD8 being up-regulated in high-fertility boars. BAG4 (BCL2-associated athanogene 4) was found to be up-regulated in high-fertility boars, forms a complex with annexin 7 and is positively correlated with levels of Hsp70 and Bcl-2, these two last inhibiting the Cytochrome C dependent activation of Caspase-3, and also involved in mitochondrial apoptosis [88].

4. Materials and Methods

4.1. Ethics Statement

Animal husbandry and experimental handling were performed in compliance with the European Community (Directive 2010/63/EU) and current Swedish legislation (SJVFS 2017:40). All procedures were performed according to the European Directive 2010/63/EU, 22/09/2010 for animal experiments and current Swedish legislation (SJVFS 2017:40. Date: 06/2016) and approved by the Bioethics Committee of Murcia University (research code: 639/2012) and the “Regional Committee for Ethical Approval of Animal Experiments” (Linköpings Djurförsöksetiska nämnd) in Linköping, Sweden (permits no. 75-12 and no. ID1400. Date: (12/2017)).
The chemicals used in the experiments were of analytical grade. Unless otherwise stated, all reagents were acquired from Sigma-Aldrich (St. Louis, MO, USA).

4.2. Ejaculated Spermatozoa

For experiment 1, young mature boars (9–11 mo, n = 3) of proven sperm quality (concentration, motility and morphology) of Swedish Landrace breed were recruited from a controlled breeding farm. The animals were individually kept in separate pens at the Translational Medicine Centre (TMC/CBR-3) of Linköping University under controlled temperature and light regimes (12 h:12 h light/dark cycle). For experiment 2, healthy, mature (1–2 years old) boars of Landrace and Large White breeds with proven fertility were housed in a commercial artificial insemination (AI) enterprise (Topigs Norsvin España, Calasparra, Murcia, Spain). All boars were fed with commercial feedstuff well-adjusted to nutritional requirements of adult AI boars, provided with water ad libitum and receiving the same management. Throughout all experiments, animals were handled carefully and in such a way as to avoid any unnecessary stress. The sperm-rich fraction (SRF) was collected using the gloved-hand method yielding, in experiment 1 a total of thirty ejaculates (ten ejaculates per boar) and, in experiment 2, 28 ejaculates from 7 different boars (one ejaculate/boar/month). All ejaculates used in the experiments fulfilled the standards of quantity and sperm quality requirements for AI-doses (˃ 200 × 106 spermatozoa/mL, 70% motile spermatozoa, and 75% of morphologically normal sperm). In all cases, spermatozoa were separated from seminal plasma immediately after ejaculate collection by double centrifugation (1500×g at room temperature for 10 min, Rotofix 32A, Hettich Zentrifugen, Tuttlingen, Germany) and the resulting sperm pellets were stored at −80°C (Ultra Low Freezer; Haier Inc., Qingdao, China).

4.3. Experimental Design

4.3.1. Experiment 1. Screening of Sperm miRNAs in the Sperm-rich Fraction (SRF) of the Boar Ejaculate

Each collected SRF was evaluated for sperm concentration and motility (velocity and forward progressive motility) using a light microscope (Zeiss, Stockholm Sweden) equipped with a thermal plate (38 °C), positive phase contrast optics (10× objective), a Charge Coupled Device (CCD) camera (UI-1540LE-M-HQ, Ueye, IDS Imaging Development Systems GmbH, Ubersulm, Germany), and the Qualisperm® Software (Biophos SA, Lausanne, Switzerland). Ejaculates with at least 70% motile and 75% morphologically normal-looking spermatozoa immediately after collection were used for the experiments. Ten boar ejaculates (SRF) from each male (n = 3) were used for analyze the miRNA expression profile, after RNA isolation.

4.3.2. Experiment 2: Assessment of the Cargo of RNAs in Ejaculated Spermatozoa from Breeding Boars with High- or Low-fertility (as Farrowing Rate and Litter Size) after AI.

Ejaculated spermatozoa came from 7 AI boars routinely used in a commercial AI program. Noteworthy, all boars presented normal semen quality, as expected for boars selected for breeding via commercial AI. The ejaculates (n = 28) used in the experiment were collected over a period of four months (one ejaculation per boar and month), while in the same period two other ejaculates were collected per week from the same boars and used to inseminate a total of 1722 sows, whose fertility data in terms of farrowing rate and litter size were recorded. The raw fertility data were statistically corrected for parameters related to farm and sow in order to isolate the contribution of the individual boars (direct boar effect [44]). The direct boar effect on fertility was shown as the deviation of both the farrowing rate and litter size of codified selected boars including a correction from the boar population of the same genetic line (Figure 4).
Thus, the 7 boars were classified as having high (n = 4) or low (n = 3) fertility based in more than 100 artificial inseminations (up to 540) per boar, according to their farrowing rate (FR) and litter size (LS): high-fertility boars with a FR > 0.45 and a LS > 0.15; and low-fertility boars with FR < −0.12 and LS < −0.18. Four ejaculates from each boar were used for analyzing the differential gene expression profile, after RNA isolation.

4.4. RNA Isolation

For both experiments, the RNA was obtained using a commercially available kit designed to isolate cell-free RNAs for low quantity samples (miRNeasy kit, Qiagen, Hilden, Germany), following the manufacturer protocol. In brief, each sample (200 μL) was thawed on ice, homogenized and RNA was extracted using 200 μL of QIAzol Lysis Reagent (Qiagen, Hilden, Germany). Samples were mixed by pipetting and incubated (5 min, RT) followed by addition of 200 μL of chloroform. The mixture was shaken vigorously (15 s, RT) and incubated (3 min, RT). Phases were separated by centrifugation (12,000× g, 15 min, 4 °C) and the upper aqueous phase (400 μL) was transferred into a gDNA eliminator column (Quiagen), to avoid DNA contamination in the purified RNA. After centrifugation (8000× g, 30 s, 4 °C), the resulting phase, free of genomic DNA, were mixed with 600 μL of 100% ethanol by gentle pipetting. A total of 700 μL of sample were transferred into a RNeasy MinElute spin column (Qiagen, Hilden, Germany) in a 2 mL collection tube and centrifuged (8000× g, 15 s, RT) in order to allow that RNA bound to the membrane. The flow-through was removed and this step was repeated with the rest of the sample (300 μL) using the same column. Later, the column was sequentially washed with 700 μL of RWT buffer (8000× g, 15 s, RT), 500 μL of RPE buffer (8000× g, 15 s, RT) and 500 μL of 80% ethanol (8000× g, 2 min, RT), discarding the flow-through in each washing. Finally, the column was dried (15,000× g, 6 min, RT) and the flow-through was discarded. RNA was eluted by placing 14 μL nuclease-free water in the spin column membrane and centrifuging (15,000× g, 1 min, RT). Total RNA content and its quality was determined by NanoDrop® 1000 (Thermo Fisher Scientific, Waltham, MA, USA). Only samples with high RNA quality were snap-frozen in liquid nitrogen and stored at −80 °C for its further analysis.

4.5. MiRNA Microarray Protocol (Experiment 1)

MiRNA expression profiling was carried out using Affymetrix® GeneChip miRNA 4.0 Array (Thermo Fisher Scientific, Göteborg, Sweden), following the manufacturer guidelines. Biotin-labeled RNA was synthesized using Affymetrix® FlashTagTM Biotin HSR RNA Labelling Kit (Thermo Fisher Scientific, Göteborg, Sweden). Briefly, 100 ng of total RNA were subjected to poly-A tail incorporation to 3′-end and, subsequently, a biotin-labeled 3DNA molecule was joined to the 3′-end using DNA ligase.

4.6. Microarray Protocol (Experiment 2)

Total RNA (75 ng) from each sample were used to make cDNA using GeneChip® Whole Transcript Plus reagent kit (Thermo Fisher Scientific, Göteborg, Sweden) following the manufacturer protocol [4].

4.7. miRNA Arrays and Microarrays Hybridization and Scanning

Finally, 21.5 μL of biotin-labeled sample (experiment 1: miRNA arrays) or 41 μL of fragmented and labelled single stranded cDNA (experiment 2: microarrays) were mixed with 110.5 μL or 109 μL, respectively, of hybridization master mix to make a cocktail hybridization mix for a single reaction. The hybridization cocktail was then incubated first at 99 °C for 5 min, followed by a descent to 45 °C until loaded on the array chip (GeneChip miRNA 4.0 for miRNA and the Affymetrix GeneChip® Porcine Gene 1.0 ST Array for gene expression; Thermo Fisher Scientific, Gothenburg, Sweden). A total of 130 μL of the cocktail hybridization mix was loaded into each array chip and was incubated at 45 °C under 60 rotations per min, for 18 h (16 h for the microarrays in experiment 2). The hybridized cartridge array chip was then unloaded and subjected to washing and staining using a GeneChip® Fluidics Station 450 (Thermo Fisher Scientific, Göteborg, Sweden), to be finally scanned using the Affymetrix GeneChip® scanner GCS3000 (Thermo Fisher Scientific, Göteborg, Sweden).

4.8. Bioinformatics and Statistics

The intensity data of each array chip was processed using the robust multi-array average (RMA) normalization, computing average intensity values by background adjustment, quantile normalization among arrays and finally log2 transformation for extracting the expression values of each transcript in the probe set, as implemented in the official Transcriptome Analysis Console (TAC; version 4.0) from Affymetrix. Molecular and biological functions as well as of overrepresentation analysis were identified using PANTHER GO [89].
Statistical analysis was performed by Affymetrix® Transcriptome Analysis Console software (version 4.0, Thermo Fisher Scientific, Göteborg, Sweden). Experiment 1: a pool of 10 ejaculates for each male (n = 3) were used for miRNA analysis. The average expression in the arrays (Mean ± Standard Deviation of the log2 average expression) were used for describing the differential expression intensities analyzed by the miRNAs array. Experiment 2: one comparation were performed between the high (n = 4, and four ejaculates per each boar) and low (n = 3, and four ejaculates per each boar) fertility boars. The differential expressed genes (DEGs) that fulfilled the criteria (ANOVA p-value < 0.05) were extracted as significantly up- or down-regulated transcripts. Finally, the analysis of the target genes, with a score of more than 0.9. in the miRDB [80] was used for finding the relation between miRNAs (experiment 1) and transcripts (experiment 2).

5. Conclusions

An overrepresentation analysis of the protein class (PANTHER) identified significant fold-increases for C-C chemokine binding (GO:0019957): CCR7 which activates B- and T-lymphocytes, 8-fold increase), XCR1 and CXCR4 (with ubiquitin as natural ligand, 1.24-fold increase), cytokine receptor activity (GO:0005126): IL23R receptor of the IL23 protein, associated to the JAK2 and STAT3, 3.4-fold increase), the TGF-receptor (PC00035) genes ACVR1C and ACVR2B (12-fold increase). Moreover, two micro-RNAs (miR-221 and mir-621) were down- and up-regulated, respectively, in high-fertile males. In conclusion, boars with different fertility performance depict a differential RNA cargo that would be an attractive battery of non-invasive molecular markers predicting fertility.

Supplementary Materials

Supplementary materials can be found at https://0-www-mdpi-com.brum.beds.ac.uk/1422-0067/21/5/1572/s1.

Author Contributions

Conceptualization, H.R.-M., D.W., J.R. and M.A.-R.; methodology, C.M., M.A.-R., I.B. and H.R.-M.; software, C.M. and M.A.R.; validation, M.A.-R.; data curation, M.A.-R..; writing—original draft preparation, M.A.R. and C.M.; writing—review and editing, D.W., J.R. and H.R.-M.; visualization, M.A.-R.; supervision, J.R. and H.R.-M.; project administration, H.R.-M.; funding acquisition, J.R., H.R.-M. and M.A.-R. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Research Council FORMAS, Stockholm (Project 2017-00946 and Project 2019-00288), the Swedish Research Council (Vetenskapsrådet, VR; project 2015-05919), the MINECO and FEDER (AGL2015-69738-R) Madrid (Spain) and the Seneca Foundation (20780/PD/18) Murcia (Spain).

Acknowledgments

To Topigs Norsvin España for supplying the boar ejaculates. To Annette Molbaek and Åsa Schippert, from the genomics Core Facility at LiU for expert assistance when running the microarrays. To Linköping University for supporting the open access funding.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Baker, R.D.; Degen, A.A. Transport of live and dead boar spermatozoa within the reproductive tract of gilts. J. Reprod. Fertil. 1972, 28, 369–377. [Google Scholar] [CrossRef] [Green Version]
  2. Langendijk, P.; Soede, N.M.; Kemp, B. Uterine activity, sperm transport, and the role of boar stimuli around insemination in sows. Theriogenology 2005, 63, 500–513. [Google Scholar] [CrossRef] [PubMed]
  3. Alminana, C.; Caballero, I.; Heath, P.R.; Maleki-Dizaji, S.; Parrilla, I.; Cuello, C.; Gil, M.A.; Vazquez, J.L.; Vazquez, J.M.; Roca, J.; et al. The battle of the sexes starts in the oviduct: Modulation of oviductal transcriptome by X and Y-bearing spermatozoa. BMC Genom. 2014, 15, 293. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Alvarez-Rodriguez, M.; Atikuzzaman, M.; Venhoranta, H.; Wright, D.; Rodriguez-Martinez, H. Expression of Immune Regulatory Genes in the Porcine Internal Genital Tract Is Differentially Triggered by Spermatozoa and Seminal Plasma. Int. J. Mol. Sci. 2019, 20, 513. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Rodriguez-Martinez, H.; Kvist, U.; Saravia, F.; Wallgren, M.; Johannisson, A.; Sanz, L.; Pena, F.J.; Martinez, E.A.; Roca, J.; Vazquez, J.M.; et al. The physiological roles of the boar ejaculate. Soc. Reprod. Fertil. Suppl. 2009, 66, 1–21. [Google Scholar] [PubMed]
  6. Robertson, S.A.; Sharkey, D.J. The role of semen in induction of maternal immune tolerance to pregnancy. Semin. Immunol. 2001, 13, 243–254. [Google Scholar] [CrossRef] [PubMed]
  7. Rodriguez-Martinez, H.; Courtens, J.L.; Kvist, U.; Ploen, L. Immunocytochemical localization of nuclear protamine in boar spermatozoa during epididymal transit. J. Reprod. Fertil. 1990, 89, 591–595. [Google Scholar] [CrossRef] [Green Version]
  8. Grunewald, S.; Paasch, U.; Glander, H.J.; Anderegg, U. Mature human spermatozoa do not transcribe novel RNA. Andrologia 2005, 37, 69–71. [Google Scholar] [CrossRef]
  9. Yang, C.C.; Lin, Y.S.; Hsu, C.C.; Wu, S.C.; Lin, E.C.; Cheng, W.T.K. Identification and sequencing of remnant messenger RNAs found in domestic swine (Sus scrofa) fresh ejaculated spermatozoa. Anim. Reprod. Sci. 2009, 113, 143–155. [Google Scholar] [CrossRef]
  10. 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]
  11. Zhang, C. Novel functions for small RNA molecules. Curr. Opin. Mol. Ther. 2009, 11, 641–651. [Google Scholar] [PubMed]
  12. Miller, D. Analysis and Significance of Messenger RNA in Human Ejaculated Spermatozoa. Mol. Reprod. Dev. 2000, 56, 259–264. [Google Scholar] [CrossRef]
  13. Brunet, A.; Sweeney, L.B.; Sturgill, J.F.; Chua, K.F.; Greer, P.L.; Lin, Y.; Tran, H.; Ross, S.E.; Mostoslavsky, R.; Cohen, H.Y.; et al. Stress-dependent regulation of FOXO transcription factors by the SIRT1 deacetylase. Science 2004, 303, 2011–2015. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Luo, Z.; Liu, Y.; Chen, L.; Ellis, M.; Li, M.; Wang, J.; Zhang, Y.; Fu, P.; Wang, K.; Li, X.; et al. microRNA profiling in three main stages during porcine spermatogenesis. J. Assist. Reprod. Genet. 2015, 32, 451–460. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Luo, L.; Ye, L.; Liu, G.; Shao, G.; Zheng, R.; Ren, Z.; Zuo, B.; Xu, D.; Lei, M.; Jiang, S.; et al. Microarray-Based approach identifies differentially expressed MicroRNAs in porcine sexually immature and mature testes. PLoS ONE 2010, 5, e11744. [Google Scholar] [CrossRef] [PubMed]
  16. Curry, E.; Safranski, T.J.; Pratt, S.L. Differential expression of porcine sperm microRNAs and their association with sperm morphology and motility. Theriogenology 2011, 76, 1532–1539. [Google Scholar] [CrossRef]
  17. Tian, H.; Li, Z.; Peng, D.; Bai, X.; Liang, W. Expression difference of miR-10b and miR-135b between the fertile and infertile semen samples (p). Forensic Sci. Int. Genet. Suppl. Ser. 2017, 6, e257–e259. [Google Scholar] [CrossRef] [Green Version]
  18. Salas-Huetos, A.; Blanco, J.; Vidal, F.; Godo, A.; Grossmann, M.; Pons, M.C.; F-Fernández, S.; Garrido, N.; Anton, E. Spermatozoa from patients with seminal alterations exhibit a differential micro-ribonucleic acid profile. Fertil. Steril. 2015, 104, 591–601. [Google Scholar] [CrossRef] [Green Version]
  19. Swindle, M.M.; Makin, A.; Herron, A.J.; Clubb, F.J.; Frazier, K.S. Swine as Models in Biomedical Research and Toxicology Testing. Vet. Pathol. 2012, 49, 344–356. [Google Scholar] [CrossRef]
  20. Roca, J.; Parrilla, I.; Bolarin, A.; Martinez, E.A.; Rodriguez-Martinez, H. Will AI in pigs become more efficient? Theriogenology 2016, 86, 187–193. [Google Scholar] [CrossRef] [Green Version]
  21. Alcivar, A.A.; Hake, L.E.; Millette, C.F.; Trasler, J.M.; Hecht, N.B. Mitochondrial gene expression in male germ cells of the mouse. Dev. Biol. 1989, 135, 263–271. [Google Scholar] [CrossRef]
  22. Gur, Y.; Breitbart, H. Mammalian sperm translate nuclear-encoded proteins by mitochondrial-type ribosomes. Genes Dev. 2006, 20, 411–416. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Gòdia, M.; Estill, M.; Castelló, A.; Balasch, S.; Rodríguez-Gil, J.E.; Krawetz, S.A.; Sánchez, A.; Clop, A. A RNA-seq analysis to describe the boar sperm transcriptome and its seasonal changes. Front. Genet. 2019, 10, 299. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Bustin, S.A.; Beaulieu, J.F.; Huggett, J.; Jaggi, R.; Kibenge, F.S.B.; Olsvik, P.A.; Penning, L.C.; Toegel, S. MIQE précis: Practical implementation of minimum standard guidelines for fluorescence-based quantitative real-time PCR experiments. BMC Mol. Biol. 2010, 11, 1–5. [Google Scholar] [CrossRef] [Green Version]
  25. Jiao, Z.J.; Yi, W.; Rong, Y.W.; Kee, J.D.; Zhong, W.X. MicroRNA-1285 regulates 17β-estradiol-inhibited immature boar sertoli cell proliferation via adenosine monophosphate-Activated protein kinase activation. Endocrinology 2015, 156, 4059–4070. [Google Scholar] [CrossRef] [Green Version]
  26. Ayaz, L.; Dinç, E. Evaluation of microRNA responses in ARPE-19 cells against the oxidative stress. Cutan. Ocul. Toxicol. 2018, 37, 121–126. [Google Scholar] [CrossRef]
  27. Jin, W.; Chen, F.; Wang, K.; Song, Y.; Fei, X.; Wu, B. miR-15a/miR-16 cluster inhibits invasion of prostate cancer cells by suppressing TGF-β signaling pathway. Biomed. Pharmacother. 2018, 104, 637–644. [Google Scholar] [CrossRef]
  28. Wang, Z.; Li, Q.; Chamba, Y.; Zhang, B.; Shang, P.; Zhang, H.; Wu, C. Identification of Genes Related to Growth and Lipid Deposition from Transcriptome Profiles of Pig Muscle Tissue. PLoS ONE 2015, 10, e0141138. [Google Scholar] [CrossRef] [Green Version]
  29. He, K.; Niu, G.; Gao, J.; Liu, J.X.; Qu, H. MicroRNA-92 expression may be associated with reduced estrogen receptor β1 mRNA levels in cervical portion of uterosacral ligaments in women with pelvic organ prolapse. Eur. J. Obstet. Gynecol. Reprod. Biol. 2016, 198, 94–99. [Google Scholar] [CrossRef]
  30. Shin, V.Y.; Siu, M.-T.; Liu, X.; Ng, E.K.O.; Kwong, A.; Chu, K.-M. MiR-92 suppresses proliferation and induces apoptosis by targeting EP4/Notch1 axis in gastric cancer. Oncotarget 2018, 9, 24209. [Google Scholar] [CrossRef] [Green Version]
  31. Qiu, T.; Wang, K.; Li, X.; Jin, J. MiR-671-5p inhibits gastric cancer cell proliferatioand promotes cell apoptosis by targeting URGCP. Exp. Ther. Med. 2018, 16, 4753–4758. [Google Scholar] [PubMed]
  32. Li, X.; Nie, C.; Tian, B.; Tan, X.; Han, W.; Wang, J.; Jin, Y.; Li, Y.; Guan, X.; Hong, A.; et al. MiR-671-5p blocks the progression of human esophageal squamous cell carcinoma by suppressing FGFR2. Int. J. Biol. Sci. 2019, 15, 1892–1904. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Li, S.; Lan, Z.-J.; Li, X.; Lin, J.; Lei, Z. Role of postnatal expression of fgfr1 and fgfr2 in testicular germ cells on spermatogenesis and fertility in mice. J. Reprod. Infertil. 2014, 15, 122–133. [Google Scholar] [PubMed]
  34. Xie, M.Y.; Hou, L.J.; Sun, J.J.; Zeng, B.; Xi, Q.Y.; Luo, J.Y.; Chen, T.; Zhang, Y.L. Porcine Milk Exosome MiRNAs Attenuate LPS-Induced Apoptosis through Inhibiting TLR4/NF-κB and p53 Pathways in Intestinal Epithelial Cells. J. Agric. Food Chem. 2019, 67, 9477–9491. [Google Scholar] [CrossRef]
  35. Cai, Q.; Zhao, A.; Ren, L.G.; Chen, J.; Liao, K.S.; Wang, Z.S.; Zhang, W. MiR-425 Involves in the Development and Progression of Renal Cell Carcinoma by Inhibiting E2F6. Eur. Rev. Med. Pharmacol. Sci. 2018, 22, 6300–6307. [Google Scholar]
  36. Chen, F.F.; Xiong, Y.; Peng, Y.; Gao, Y.; Qin, J.; Chu, G.Y.; Pang, W.J.; Yang, G.S. MiR-425-5p inhibits differentiation and proliferation in porcine intramuscular preadipocytes. Int. J. Mol. Sci. 2017, 18, 2101. [Google Scholar] [CrossRef] [Green Version]
  37. Lu, Y.; Wu, X.; Wang, J. Correlation of miR-425-5p and IL-23 with pancreatic cancer. Oncol. Lett. 2019, 17, 4595–4599. [Google Scholar] [CrossRef]
  38. Kastelein, R.A.; Hunter, C.A.; Cua, D.J. Discovery and Biology of IL-23 and IL-27: Related but Functionally Distinct Regulators of Inflammation. Annu. Rev. Immunol. 2007, 25, 221–242. [Google Scholar] [CrossRef] [Green Version]
  39. Zhang, X.; Wu, M.; Chong, Q.Y.; Zhang, W.; Qian, P.; Yan, H.; Qian, W.; Zhang, M.; Lobie, P.E.; Zhu, T. Amplification of hsa-miR-191/425 locus promotes breast cancer proliferation and metastasis by targeting DICER1. Carcinogenesis 2018, 39, 1506–1516. [Google Scholar] [CrossRef]
  40. Gu, Y.; Ampofo, E.; Menger, M.D.; Laschke, M.W. MIR-191 suppresses angiogenesis by activation of NF-kB signaling. FASEB J. 2017, 31, 3321–3333. [Google Scholar] [CrossRef] [Green Version]
  41. Tian, X.; Zhang, Z. miR-191/DAB2 axis regulates the tumorigenicity of estrogen receptor-positive breast cancer. IUBMB Life 2018, 70, 71–80. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Zhou, L.; Zhang, F.; Tong, J.; Liu, F. MiR-191-5p inhibits lung adenocarcinoma by repressing SATB1 to inhibit Wnt pathway. Mol. Genet. Genomic Med. 2020, 8, e1043. [Google Scholar] [CrossRef] [PubMed]
  43. Broekhuijse, M.; Sostaric, E.; Feitsma, H.; Gadella, B. The value of microscopic semen motility assessment at collection for a commercial artificial insemination center, a retrospective study on factors explaining variation in pig fertility. Theriogenology 2012, 77, 1466–1479. [Google Scholar] [CrossRef] [PubMed]
  44. Riesco, M.F.; Valcarce, D.G.; Martínez-Vázquez, J.M.; Robles, V. Effect of low sperm quality on progeny: A study on zebrafish as model species. Sci. Rep. 2019, 9, 11192. [Google Scholar] [CrossRef] [Green Version]
  45. Abu-Halima, M.; Hammadeh, M.; Schmitt, J.; Leidinger, P.; Keller, A.; Meese, E.; Backes, C. Altered microRNA expression profiles of human spermatozoa in patients with different spermatogenic impairments. Fertil. Steril. 2013, 99, 1249–1255. [Google Scholar] [CrossRef]
  46. Michailov, Y.; Ickowicz, D.; Breitbart, H. Zn2+-stimulation of sperm capacitation and of the acrosome reaction is mediated by EGFR activation. Dev. Biol. 2014, 396, 246–255. [Google Scholar] [CrossRef] [Green Version]
  47. Lian, N.; Niu, Q.; Lei, Y.; Li, X.; Li, Y.; Song, X. MiR-221 is involved in depression by regulating Wnt2/CREB/BDNF axis in hippocampal neurons. Cell Cycle 2018, 17, 2745–2755. [Google Scholar] [CrossRef] [Green Version]
  48. Jacob, A. Molecular characterization of a voltage-gated potassium channel expressed in rat testis. Mol. Hum. Reprod. 2000, 6, 303–313. [Google Scholar] [CrossRef] [Green Version]
  49. Vicente-Carrillo, A.; Álvarez-Rodríguez, M.; Rodríguez-Martínez, H. The CatSper channel modulates boar sperm motility during capacitation. Reprod. Biol. 2017, 17, 69–78. [Google Scholar] [CrossRef] [Green Version]
  50. Wang, H.; Liu, J.; Cho, K.-H.; Ren, D. A Novel, Single, Transmembrane Protein CATSPERG Is Associated with CATSPER1 Channel Protein1. Biol. Reprod. 2009, 81, 539–544. [Google Scholar] [CrossRef] [Green Version]
  51. Ghanbari, H.; Keshtgar, S.; Zare, H.R.; Gharesi-Fard, B. Inhibition of CatSper and Hv1 Channels and NOX5 Enzyme Affect Progesterone-Induced Increase of Intracellular Calcium Concentration and ROS Generation in Human Sperm. Iran. J. Med. Sci. 2019, 44, 127–134. [Google Scholar] [PubMed]
  52. Ren, D.; Xia, J. Calcium signaling through CatSper channels in mammalian fertilization. Physiology 2010, 25, 165–175. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Cassandri, M.; Smirnov, A.; Novelli, F.; Pitolli, C.; Agostini, M.; Malewicz, M.; Melino, G.; Raschellà, G. Zinc-finger proteins in health and disease. Cell Death Discov. 2017, 3, 1–2. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Tonn Eisinger, K.R.; Woolfrey, K.M.; Swanson, S.P.; Schnell, S.A.; Meitzen, J.; Dell’Acqua, M.; Mermelstein, P.G. Palmitoylation of caveolin-1 is regulated by the same DHHC acyltransferases that modify steroid hormone receptors. J. Biol. Chem. 2018, 293, 15901–15911. [Google Scholar] [CrossRef] [Green Version]
  55. Aquila, S.; Sisci, D.; Gentile, M.; Middea, E.; Catalano, S.; Carpino, A.; Rago, V.; Andò, S. Estrogen Receptor (ER)α and ERβ Are Both Expressed in Human Ejaculated Spermatozoa: Evidence of Their Direct Interaction with Phosphatidylinositol-3-OH Kinase/Akt Pathway. J. Clin. Endocrinol. Metab. 2004, 89, 1443–1451. [Google Scholar] [CrossRef] [Green Version]
  56. Gou, D.M.; Sun, Y.; Gao, L.; Chow, L.M.C.; Huang, J.; Feng, Y.D.; Jiang, D.H.; Li, W.X. Cloning and characterization of a novel Krüppel-like zinc finger gene, ZNF268, expressed in early human embryo. Biochim. Biophys. Acta-Gene Struct. Expr. 2001, 1518, 306–310. [Google Scholar] [CrossRef]
  57. Kroft, T.L.; Patterson, J.; Won Yoon, J.; Doglio, L.; Walterhouse, D.O.; Iannaccone, P.M.; Goldberg, E. GLI1 Localization in the Germinal Epithelial Cells Alternates Between Cytoplasm and Nucleus: Upregulation in Transgenic Mice Blocks Spermatogenesis in Pachytene1. Biol. Reprod. 2001, 65, 1663–1671. [Google Scholar] [CrossRef] [Green Version]
  58. Turner, T.T.; Bang, H.J.; Attipoe, S.A.; Johnston, D.S.; Tomsig, J.L. Sonic hedgehog pathway inhibition alters epididymal function as assessed by the development of sperm motility. J. Androl. 2006, 27, 225–232. [Google Scholar] [CrossRef]
  59. Edwards, D.R.; Handsley, M.M.; Pennington, C.J. The ADAM metalloproteinases. Mol. Asp. Med. 2009, 29, 258–289. [Google Scholar] [CrossRef]
  60. Atikuzzaman, M.; Alvarez-Rodriguez, M.; Vicente-Carrillo, A.; Johnsson, M.; Wright, D.; Rodriguez-Martinez, H. Conserved gene expression in sperm reservoirs between birds and mammals in response to mating. BMC Genom. 2017, 18, 98. [Google Scholar] [CrossRef] [Green Version]
  61. Atikuzzaman, M.; Hoglund, A.; Johnsson, M.; Alvarez-Rodriguez, M.; Mehta, R.; Wright, D.; Rodriguez-Martinez, H. Spermatozoa rather than seminal fluid drives immune-modulation in the sperm storage tubuli of Red Jungle Fowl hens. Reprod. Domest. Anim. 2015, 50, 44. [Google Scholar]
  62. Choi, H.; Han, C.; Jin, S.; Kwon, J.T.; Kim, J.; Jeong, J.; Kim, J.; Ham, S.; Jeon, S.; Yoo, Y.J.; et al. Reduced Fertility and Altered Epididymal and Sperm Integrity in Mice Lacking ADAM71. Biol. Reprod. 2015, 93, 1–11. [Google Scholar] [CrossRef]
  63. Wei, X.; Moncada-Pazos, A.; Cal, S.; Soria-Valles, C.; Gartner, J.; Rudloff, U.; Lin, J.C.; Rosenberg, S.A.; López-Otín, C.; Samuels, Y. Analysis of the disintegrin-metalloproteinases family reveals ADAM29 and ADAM7 are often mutated in melanoma. Hum. Mutat. 2011, 32, E2148–E2175. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Jinawath, N.; Shiao, M.S.; Norris, A.; Murphy, K.; Klein, A.P.; Yonescu, R.; Iacobuzio-Donahue, C.; Meeker, A.; Jinawath, A.; Yeo, C.J.; et al. Alterations of type II classical cadherin, cadherin-10 (CDH10), is associated with pancreatic ductal adenocarcinomas. Genes Chromosom. Cancer 2017, 56, 427–435. [Google Scholar] [CrossRef] [PubMed]
  65. Wiśniewski, A.; Matusiak, Ł.; Szczerkowska-Dobosz, A.; Nowak, I.; Kuśnierczyk, P. HLA-C*06:02-independent, gender-related association of PSORS1C3 and PSORS1C1/CDSN single-nucleotide polymorphisms with risk and severity of psoriasis. Mol. Genet. Genom. 2018, 293, 957–966. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Kumar, V.; Soni, U.K.; Maurya, V.K.; Singh, K.; Jha, R.K. Integrin beta8 (ITGB8) activates VAV-RAC1 signaling via FAK in the acquisition of endometrial epithelial cell receptivity for blastocyst implantation. Sci. Rep. 2017, 7, 1–18. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Santulli, G. Angiopoietin-like proteins: A comprehensive look. Front. Endocrinol. 2014, 5, 4. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Van Dijk, M.; Mulders, J.; Könst, A.; Janssens, B.; van Roy, F.; Blankenstein, M.; Oudejans, C. Differential downregulation of αT-catenin expression in placenta: Trophoblast cell type-dependent imprinting of the CTNNA3 gene. Gene Expr. Patterns 2004, 5, 61–65. [Google Scholar] [CrossRef]
  69. Chapman, E.M.; Lant, B.; Ohashi, Y.; Yu, B.; Schertzberg, M.; Go, C.; Dogra, D.; Koskimäki, J.; Girard, R.; Li, Y.; et al. A conserved CCM complex promotes apoptosis non-autonomously by regulating zinc homeostasis. Nat. Commun. 2019, 10, 1–15. [Google Scholar] [CrossRef] [Green Version]
  70. Bragantini, B.; Tiotiu, D.; Rothé, B.; Saliou, J.M.; Marty, H.; Cianférani, S.; Charpentier, B.; Quinternet, M.; Manival, X. Functional and Structural Insights of the Zinc-Finger HIT protein family members Involved in Box C/D snoRNP Biogenesis. J. Mol. Biol. 2016, 428, 2488–2506. [Google Scholar] [CrossRef]
  71. Stephan, J.S.; Fioriti, L.; Lamba, N.; Colnaghi, L.; Karl, K.; Derkatch, I.L.; Kandel, E.R. The CPEB3 Protein Is a Functional Prion that Interacts with the Actin Cytoskeleton. Cell Rep. 2015, 11, 1772–1785. [Google Scholar] [CrossRef] [PubMed]
  72. Drisaldi, B.; Colnaghi, L.; Fioriti, L.; Rao, N.; Myers, C.; Snyder, A.M.; Metzger, D.J.; Tarasoff, J.; Konstantinov, E.; Fraser, P.E.; et al. SUMOylation Is an Inhibitory Constraint that Regulates the Prion-like Aggregation and Activity of CPEB3. Cell Rep. 2015, 11, 1694–1702. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Liu, F.; Zhang, G.; Lv, S.; Wen, X.; Liu, P. miRNA-301b-3p accelerates migration and invasion of high-grade ovarian serous tumor via targeting CPEB3/EGFR axis. J. Cell. Biochem. 2019, 120, 12618–12627. [Google Scholar] [CrossRef]
  74. Wilson, V.G. Introduction to sumoylation. In Advances in Experimental Medicine and Biology; Springer: New York, NY, USA, 2017; Volume 963, pp. 1–12. [Google Scholar]
  75. Mei, J.; Yan, Y.; Li, S.Y.; Zhou, W.J.; Zhang, Q.; Li, M.Q.; Sun, H.X. CXCL16/CXCR6 interaction promotes endometrial decidualization via the PI3K/AKT pathway. Reproduction 2019, 157, 273–282. [Google Scholar] [CrossRef] [Green Version]
  76. Guerin, L.R.; Moldenhauer, L.M.; Prins, J.R.; Bromfield, J.J.; Hayball, J.D.; Robertson, S.A. Seminal Fluid Regulates Accumulation of FOXP3+ Regulatory T Cells in the Preimplantation Mouse Uterus Through Expanding the FOXP3+ Cell Pool and CCL19-Mediated Recruitment1. Biol. Reprod. 2011, 85, 397–408. [Google Scholar] [CrossRef] [Green Version]
  77. Krikun, G. The CXL12/CXCR4/CXCR7 axis in female reproductive tract disease: Review. Am. J. Reprod. Immunol. 2018, 80, e13028. [Google Scholar] [CrossRef]
  78. Mendis, S.H.S.; Meachem, S.J.; Sarraj, M.A.; Loveland, K.L. Activin A Balances Sertoli and Germ Cell Proliferation in the Fetal Mouse Testis. Biol. Reprod. 2011, 84, 379–391. [Google Scholar] [CrossRef] [Green Version]
  79. Kozomara, A.; Birgaoanu, M.; Griffiths-Jones, S. miRBase: From microRNA sequences to function. Nucleic Acids Res. 2019, 47, D155–D162. [Google Scholar] [CrossRef]
  80. Yuan, S.; Swiggin, H.M.; Zheng, H.; Yan, W. A testis-specific gene, Ubqlnl, is dispensable for mouse embryonic development and spermatogenesis. Mol. Reprod. Dev. 2015, 82, 408–409. [Google Scholar] [CrossRef] [Green Version]
  81. Marín, I. The ubiquilin gene family: Evolutionary patterns and functional insights. BMC Evol. Biol. 2014, 14, 63. [Google Scholar] [CrossRef] [Green Version]
  82. Li, X.; Lan, Y.; Xu, J.; Zhang, W.; Wen, Z. SUMO1-activating enzyme subunit 1 is essential for the survival of hematopoietic stem/progenitor cells in zebrafish. Development 2012, 139, 4321–4329. [Google Scholar] [CrossRef] [Green Version]
  83. Shin, E.J.; Shin, H.M.; Nam, E.; Kim, W.S.; Kim, J.H.; Oh, B.H.; Yun, Y. DeSUMOylating isopeptidase: A second class of SUMO protease. EMBO Rep. 2012, 13, 339–346. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  84. Judith, D.; Jefferies, H.B.J.; Boeing, S.; Frith, D.; Snijders, A.P.; Tooze, S.A. ATG9A shapes the forming autophagosome through Arfaptin 2 and phosphatidylinositol 4-kinase IIIβ. J. Cell Biol. 2019, 218, 1634–1652. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  85. Fukamachi, K.; Matsuoka, Y.; Ohno, H.; Hamaguchi, T.; Tsuda, H. Neuronal leucine-rich repeat protein-3 amplifies MAPK activation by epidermal growth factor through a carboxyl-terminal region containing endocytosis motifs. J. Biol. Chem. 2002, 277, 43549–43552. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  86. Dai, Y.; Duan, H.; Duan, C.; Zhu, H.; Zhou, R.; Pei, H.; Shen, L. TCF21 functions as a tumor suppressor in colorectal cancer through inactivation of PI3K/AKT signaling. Onco Targets Ther. 2017, 10, 1603–1611. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  87. Muhammad Aslam, M.K.; Kumaresan, A.; Yadav, S.; Mohanty, T.K.; Datta, T.K. Comparative proteomic analysis of high- and low-fertile buffalo bull spermatozoa for identification of fertility-associated proteins. Reprod. Domest. Anim. 2019, 54, 786–794. [Google Scholar]
  88. Du, Y.; Huang, Y.; Gao, Y.; Song, B.; Mao, J.; Chen, L.; Bai, L.; Tang, J. Annexin A7 modulates BAG4 and BAG4-binding proteins in mitochondrial apoptosis. Biomed. Pharmacother. 2015, 74, 30–34. [Google Scholar] [CrossRef]
  89. Thomas, P.D.; Campbell, M.J.; Kejariwal, A.; Mi, H.; Karlak, B.; Daverman, R.; Diemer, K.; Muruganujan, A.; Narechania, A. PANTHER: A library of protein families and subfamilies indexed by function. Genome Res. 2003, 13, 2129–2141. [Google Scholar] [CrossRef] [Green Version]
Ijms 21 01572 g001 550
Figure 1. (A) Scatter plot and (B) Vulcano plot of differential gene expression using a microarrays platform (GeneChip® Porcine Gene 1.0 ST Array) of 521 transcripts (347 transcripts up-regulated and 174 transcripts down-regulated) in high-fertile relative with low-fertile males. Annotation of the full gene name for the transcript and up- (red color-positive fold change) and down-regulation (green color-negative fold change in high-fertile males in comparison to low-fertile males. The Porcine Genome Array provides comprehensive coverage of the S. scrofa transcriptome. The array contains 23,937 probe sets that interrogate approximately 23.256 transcripts from 20.201 Sus scrofa genes. The sequence information for this array was selected from public data sources including UniGene GenBank® mRNAs, and GenBank porcine mitochondrial and rRNA sequences.
Figure 1. (A) Scatter plot and (B) Vulcano plot of differential gene expression using a microarrays platform (GeneChip® Porcine Gene 1.0 ST Array) of 521 transcripts (347 transcripts up-regulated and 174 transcripts down-regulated) in high-fertile relative with low-fertile males. Annotation of the full gene name for the transcript and up- (red color-positive fold change) and down-regulation (green color-negative fold change in high-fertile males in comparison to low-fertile males. The Porcine Genome Array provides comprehensive coverage of the S. scrofa transcriptome. The array contains 23,937 probe sets that interrogate approximately 23.256 transcripts from 20.201 Sus scrofa genes. The sequence information for this array was selected from public data sources including UniGene GenBank® mRNAs, and GenBank porcine mitochondrial and rRNA sequences.
Ijms 21 01572 g001
Ijms 21 01572 g002 550
Figure 2. Differentially expressed genes (DEGs) classification (PANTHER) according to 1) biological process (A). up-regulated genes in high-fertile boars and (B). down-regulated genes in high-fertile boars; and 2) molecular function (C). up-regulated genes in high-fertile boars and (D). down-regulated genes in high-fertile boars.
Figure 2. Differentially expressed genes (DEGs) classification (PANTHER) according to 1) biological process (A). up-regulated genes in high-fertile boars and (B). down-regulated genes in high-fertile boars; and 2) molecular function (C). up-regulated genes in high-fertile boars and (D). down-regulated genes in high-fertile boars.
Ijms 21 01572 g002
Ijms 21 01572 g003 550
Figure 3. Differential expressed genes (DEGs) classification (PANTHER) according to pathway (A). up-regulated genes in high-fertile boars and (B). down-regulated genes in high-fertile boars.
Figure 3. Differential expressed genes (DEGs) classification (PANTHER) according to pathway (A). up-regulated genes in high-fertile boars and (B). down-regulated genes in high-fertile boars.
Ijms 21 01572 g003
Ijms 21 01572 g004 550
Figure 4. Bar chat of the deviation normalized from the average for each genetic line in terms of boar effect in farrowing rate and litter size.
Figure 4. Bar chat of the deviation normalized from the average for each genetic line in terms of boar effect in farrowing rate and litter size.
Ijms 21 01572 g004
Table
Table 1. Expression (Mean ± Standard Deviation of the log2 average expression or transcript detection) of Sus scrofa specific miRNAs detected by the GeneChip™ miRNA 4.0 Array (ThermoFisher Scientific) in spermatozoa retrieved from the sperm-rich fraction (SRF) of the ejaculate of healthy mature boars (n = 3). The miRNA is designed to interrogate all mature miRNA sequences in miRBase v20. The array includes 30.424 mature miRNA (all organisms) and we select specifically 10 of the 326 Sus scrofa-specific miRNAs included in the array with expression value > 10.
Table 1. Expression (Mean ± Standard Deviation of the log2 average expression or transcript detection) of Sus scrofa specific miRNAs detected by the GeneChip™ miRNA 4.0 Array (ThermoFisher Scientific) in spermatozoa retrieved from the sperm-rich fraction (SRF) of the ejaculate of healthy mature boars (n = 3). The miRNA is designed to interrogate all mature miRNA sequences in miRBase v20. The array includes 30.424 mature miRNA (all organisms) and we select specifically 10 of the 326 Sus scrofa-specific miRNAs included in the array with expression value > 10.
Transcript IDExpression Mean ± Standard Deviation (log2) AccessionSequence LengthSequence
ssc-miR-128513.98 ± 0.13MIMAT001395424CUGGGCAACAUAGCGAGACCCCGU
ssc-miR-1612.6 ± 0.74MIMAT000775422UAGCAGCACGUAAAUAUUGGCG
ssc-miR-433212.32 ± 0.29MIMAT001796220CACGGCCGCCGCCGGGCGCC
ssc-miR-92a12.06 ± 0.09MIMAT001390822UAUUGCACUUGUCCCGGCCUGU
ssc-miR-671-5p11.73 ± 0.54MIMAT002538124AGGAAGCCCUGGAGGGGCUGGAGG
ssc-miR-4334-5p11.31 ± 0.05MIMAT001796619CCCUGGAGUGACGGGGGUG
ssc-miR-425-5p10.99 ± 0.15MIMAT001391723AAUGACACGAUCACUCCCGUUGA
ssc-miR-19110.57 ± 0.22MIMAT001387623CAACGGAAUCCCAAAAGCAGCUG
ssc-miR-92b-5p10.53 ± 0.18MIMAT001737724AGGGACGGGACGCGGUGCAGUGUU
ssc-miR-15b10.01 ± 0.9MIMAT000212522UAGCAGCACAUCAUGGUUUACA
Table
Table 2. Differential gene expression using a microarrays platform (GeneChip® Porcine Gene 1.0 ST Array) of 77 transcripts (67 transcripts up-regulated and 10 transcripts down-regulated) in high-fertile relative with low-fertile males (Fold change >2 & <–2). Annotation of the full gene name for the transcript and up- (positive fold change) and down-regulation (negative fold change in high-fertile males in comparison to low-fertile males. The Porcine Genome Array provides comprehensive coverage of the Sus scrofa transcriptome. The array contains 23,937 probe sets that interrogate approximately 23,256 transcripts from 20,201 Sus scrofa genes. The sequence information for this array was selected from public data sources including UniGene GenBank® mRNAs, and GenBank porcine mitochondrial and rRNA sequences.
Table 2. Differential gene expression using a microarrays platform (GeneChip® Porcine Gene 1.0 ST Array) of 77 transcripts (67 transcripts up-regulated and 10 transcripts down-regulated) in high-fertile relative with low-fertile males (Fold change >2 & <–2). Annotation of the full gene name for the transcript and up- (positive fold change) and down-regulation (negative fold change in high-fertile males in comparison to low-fertile males. The Porcine Genome Array provides comprehensive coverage of the Sus scrofa transcriptome. The array contains 23,937 probe sets that interrogate approximately 23,256 transcripts from 20,201 Sus scrofa genes. The sequence information for this array was selected from public data sources including UniGene GenBank® mRNAs, and GenBank porcine mitochondrial and rRNA sequences.
UP-REGULATED GENES in High-Fertile Boars vs. Low-Fertile BoarsDESCRIPTIONFold ChangeDOWN-REGULATED GENES in High-Fertile Boars vs. Low-Fertile BoarsDESCRIPTIONFold Change
LOC100739568replication factor C subunit 1-like5.72IFN-DELTA-4; LOC100736862interferon-delta-4; interferon tau-11-like−2
AGPAT21-acylglycerol-3-phosphate O-acyltransferase 23.53FOXP2; LOC100620234forkhead box P2; forkhead box protein P2-like−2.03
FBXW5F-box and WD repeat domain containing 53.52LOC100155405olfactory receptor 4K1−2.03
PPP1R16Aprotein phosphatase 1, regulatory subunit 16A3.34LOC100516093olfactory receptor 52N4−2.27
PIM1Pim-1 proto-oncogene, serine/threonine kinase3.26LOC100517289E3 ubiquitin-protein ligase RNF220−2.41
RAB7ARAB7A, member RAS oncogene family3.26LOC100738965; LOC100739445olfactory receptor 1J4-like−2.49
C12H17orf97chromosome 12 open reading frame, human C17orf973.1LOC100525437ATP synthase subunit d, mitochondrial-like−3.06
TMEM239; C17H20orf141transmembrane protein 239; chromosome 17 open reading frame, human C20orf1413.1LOC100157815olfactory receptor 4K15−3.9
SAE1SUMO1 activating enzyme subunit 13.01LOC100522664; LOC100622735olfactory receptor 2AJ1-like−4.97
TSSK6testis-specific serine kinase 62.91LOC100157785olfactory receptor 8H1-like−5.22
WNT3wingless-type MMTV integration site family, member 32.91
CARHSP1calcium regulated heat stable protein 1, 24kDa2.88
UBQLNLubiquilin-like2.79
TSSK2testis-specific serine kinase 22.77
MAP1LC3Amicrotubule-associated protein 1 light chain 3 alpha2.71
RABAC1Rab acceptor 1 (prenylated)2.68
PHKG2phosphorylase kinase, gamma 2 (testis)2.66
MTHFSDmethenyltetrahydrofolate synthetase domain containing2.61
CSNK1G2casein kinase 1, gamma 22.56
C9H11orf71chromosome 9 open reading frame, human C11orf712.54
AGBL5ATP/GTP binding protein-like 52.51
SGSM2small G protein signaling modulator 22.47
TIMP2TIMP metallopeptidase inhibitor 22.47
DUSP18dual specificity phosphatase 182.45
SARNPSAP domain containing ribonucleoprotein2.45
TNP1transition protein 1 (during histone to protamine replacement)2.45
FAM57Afamily with sequence similarity 57, member A2.44
EHD1EH-domain containing 12.43
NSUN4NOP2/Sun domain family, member 42.4
C7H6orf125uncharacterized protein C6orf125 homolog2.39
EIF4A3; LOC100739660eukaryotic translation initiation factor 4A3; eukaryotic initiation factor 4A-III-like2.38
NDUFA10NADH dehydrogenase (ubiquinone) 1 alpha subcomplex, 10, 42kDa2.36
KHDRBS3; LOC100625959KH domain containing, RNA binding, signal transduction associated 3; KH domain-containing, RNA-binding, signal transduction-associated protein 3-like2.33
CHCHD5coiled-coil-helix-coiled-coil-helix domain containing 52.32
OSBP2oxysterol binding protein 22.3
C1H9orf16; LCN2chromosome 1 open reading frame, human C9orf16; lipocalin 22.29
GRK4G protein-coupled receptor kinase 42.28
PSMF1; LOC100736938proteasome inhibitor subunit 1; proteasome inhibitor PI31 subunit2.28
CATSPERGcatsper channel auxiliary subunit gamma2.27
MFAP3Lmicrofibrillar-associated protein 3-like2.26
TSACCTSSK6 activating co-chaperone2.25
LOC100155534verprolin-like2.24
LOC102166949E3 ubiquitin-protein ligase RNF1692.24
NEURL1neuralized E3 ubiquitin protein ligase 12.24
FAM46Cfamily with sequence similarity 46, member C2.22
ISG20L2interferon stimulated exonuclease gene 20kDa-like 22.19
MIEF1mitochondrial elongation factor 12.19
HTR4; LOC1007378355-hydroxytryptamine (serotonin) receptor 4, G protein-coupled; 5-hydroxytryptamine receptor 42.18
VRK3vaccinia related kinase 32.16
CARM1coactivator-associated arginine methyltransferase 12.15
LOC100152461; LOC100627430; NFYAnuclear transcription factor Y subunit alpha; nuclear transcription factor Y, alpha2.15
LOC100157936; LOC102162592E3 ubiquitin-protein ligase NRDP1; E3 ubiquitin-protein ligase NRDP1-like2.13
LOC102162328uncharacterized LOC1021623282.13
SYCE2synaptonemal complex central element protein 22.12
UBE2D2ubiquitin-conjugating enzyme E2D 22.12
LOC100625908; LOC102162228uncharacterized LOC100625908; uncharacterized LOC1021622282.1
ACTR3BARP3 actin-related protein 3 homolog B (yeast)2.09
CCNYcyclin Y2.09
CPEB3cytoplasmic polyadenylation element binding protein 32.09
LOC100511504; LOC100515480maestro heat-like repeat-containing protein family member 12.09
LOC100737397; C1H9orf9uncharacterized protein C9orf9 homolog; chromosome 1 open reading frame, human C9orf92.07
LOC102158209uncharacterized LOC1021582092.04
REEP2receptor accessory protein 22.04
TTC7Atetratricopeptide repeat domain 7A2.04
GNA12guanine nucleotide binding protein (G protein) alpha 122.02
PIM3; BRD1Pim-3 proto-oncogene, serine/threonine kinase; bromodomain containing 12.02
DESI1desumoylating isopeptidase 12
Table
Table 3. Overrepresentation analysis (PANTHER) of the DEGs between high-fertile and low-fertile boars, according to A. PANTHER molecular function class and B. PANTHER protein class. In both overrepresentation analysis, the expected number of genes in the Sus scrofa reference genome is compared with the observed number of genes in the category, resulting in a fold enrichment is indicated in the table (p < 0.05 as significant level). The column “Genes IDs” depicts the up-regulated (bold) and down-regulated (italics) genes between high-fertile boars and low-fertile boars. (n = 7).
Table 3. Overrepresentation analysis (PANTHER) of the DEGs between high-fertile and low-fertile boars, according to A. PANTHER molecular function class and B. PANTHER protein class. In both overrepresentation analysis, the expected number of genes in the Sus scrofa reference genome is compared with the observed number of genes in the category, resulting in a fold enrichment is indicated in the table (p < 0.05 as significant level). The column “Genes IDs” depicts the up-regulated (bold) and down-regulated (italics) genes between high-fertile boars and low-fertile boars. (n = 7).
A. PANTHER Molecular Function ClassSus Scrofa RefExpectedFold EnrichmentP ValueGene IDs
C-C chemokine binding (GO:0019957)240.387.936.76E-03CCR7, XCR1, CXCR4
chemokine binding (GO:0019956)240.387.936.76E-03CCR7, XCR1, CXCR4
cytokine binding (GO:0019955)881.394.332.98E-03CCR7, XCR1, CXCR4, ACVR2B, IL23R
protein binding (GO:0005515)282444.491.514.32E-04WDR34, MTF2, ELOF1, KIF15, KAT7, FGFR2, GPAT3, ATP6V0A2, AGBL5, CCNY,BAG6, TRAF3, EHD1, WDTC1, PPP1R16A, CXCL16, MICU1, AKAP11, MAP1LC3A, SGSM2, DVL1, UBQLNL, ERC1, MOBKL3, SNX14, WNT3, BAG6, NRG4, SHC3, GNA12, UBE2D2, STYX, FIBCD1, TSACC, ATN1, CDH10, STXBP6, TMEFF2, CCR7, SERPINE2, FGD4, AKAP7, ACVR1C, CDSN, NUF2, ACVR2B, PPARGC1B, SYT1, GDF9, ANKRD55, GNAT3, ANGPTL1, WFDC9, SYT10, IL23R, SERPINI2, TFPI2, CTNNA3, XCR1, TEAD3, CXCR4, BRCC3, CDSN, CCNJ
cytokine receptor activity (GO:0005126)751.183.393.20E-02CCR7, XCR1, CXCR4, IL23R
B. PANTHER Protein ClassSus Scrofa RefExpectedFold EnrichmentP ValueGene IDs
TGF-beta receptor (PC00035)150.2412.691.83E-03ACVR1C, ACVR2B
serine/threonine protein kinase receptor (PC00205)200.329.524.10E-03ACVR1C, ACVR2B
protein kinase (PC00193)2664.192.391.05E-02TSSK2, TSSK6, RPS6KC1, PRKCQ, VRK3, CSNK1G2, GRK4, ACVR1C, ACVR2B
kinase (PC00137)3866.081.972.10E-02TSSK2, TSSK6, NDUFA10, RPS6KC1, PRKCQ, VRK3, CSNK1G2, PI4K2A, GRK4, ACVR1C, ACVR2B
transferase (PC00220)90714.291.822.79E-03NTMT1, CHST13, NSUN4, TSSK2, TSSK6, KAT7, GPAT3, METTL2A, GALNT18, NDUFA10, RPS6KC1, PRKCQ, VRK3, FAM177A1, CSNK1G2, CHPT1, PI4K2A, RNGTT, GRK4, ACCSL, ME1, ACVR1C, ACVR2B, TGM6
Table
Table 4. Target genes of the 10 miRNA Sus scrofa-specific miRNAs detected by the GeneChip™ miRNA 4.0 Array (ThermoFisher Scientific) expressed more than 10 times average expression units (log2 average expression or transcript detection average expression units; experiment 1) in relation to the DEGs using a microarrays platform (GeneChip® Porcine Gene 1.0 ST Array) of 521 transcripts between high- and low-fertile boars (experiment 2). Up-regulated in bold and down-regulated in italics.
Table 4. Target genes of the 10 miRNA Sus scrofa-specific miRNAs detected by the GeneChip™ miRNA 4.0 Array (ThermoFisher Scientific) expressed more than 10 times average expression units (log2 average expression or transcript detection average expression units; experiment 1) in relation to the DEGs using a microarrays platform (GeneChip® Porcine Gene 1.0 ST Array) of 521 transcripts between high- and low-fertile boars (experiment 2). Up-regulated in bold and down-regulated in italics.
miRNA IDGene TargetGO – Molecular FunctionGO – Biological Process
miR-16UBQLNL
(Ubiquilin-like protein)		GO:0031593
polyubiquitin modification-dependent protein binding		GO:0006511
ubiquitin-dependent protein catabolic process
CPEB3
(Cytoplasmic polyadenylation element-binding protein 3)		GO:0000900
translation repressor activity, mRNA regulatory element binding		GO:0071230
cellular response to amino acid stimulus
GO:0007616
long-term memory
DESI1
(Desumoylating isopeptidase 1)		GO:0042802
identical protein binding
ATG9A
(Autophagy-related protein 9A)		 GO:0000045
autophagosome assembly
BAG4
(BAG family molecular chaperone regulator 4)		GO:0031625
ubiquitin protein ligase binding		GO:0071364
cellular response to epidermal growth factor stimulus
GO:0045785
positive regulation of cell adhesion
LRRN3
(Leucine-rich repeat neuronal protein 3)		 GO:0051965
positive regulation of synapse assembly
miR-92aCPEB3
(Cytoplasmic polyadenylation element-binding protein 3)		GO:0000900
translation repressor activity, mRNA regulatory element binding		GO:0071230
cellular response to amino acid stimulus
GO:0007616
long-term memory
TCF21
(ranscription factor 21)		GO:0050681
androgen receptor binding		GO:0060766
negative regulation of androgen receptor signaling pathway
PDZD8
(PDZ domain-containing protein 8)		GO:0008289
lipid binding		GO:0035556
intracellular signal transduction
miR-15bUBQLNL
(Ubiquilin-like protein)		GO:0031593
polyubiquitin modification-dependent protein binding		GO:0006511
ubiquitin-dependent protein catabolic process
CPEB3
(Cytoplasmic polyadenylation element-binding protein 3)		GO:0000900
translation repressor activity, mRNA regulatory element binding		GO:0071230
cellular response to amino acid stimulus
GO:0007616
long-term memory
DESI1
(Desumoylating isopeptidase 1)		GO:0042802
identical protein binding
ATG9A
(Autophagy-related protein 9A)		 GO:0000045
autophagosome assembly
LRRN3
(Leucine-rich repeat neuronal protein 3)		 GO:0051965
positive regulation of synapse assembly

Share and Cite

MDPI and ACS Style

Alvarez-Rodriguez, M.; Martinez, C.; Wright, D.; Barranco, I.; Roca, J.; Rodriguez-Martinez, H. The Transcriptome of Pig Spermatozoa, and Its Role in Fertility. Int. J. Mol. Sci. 2020, 21, 1572. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms21051572

AMA Style

Alvarez-Rodriguez M, Martinez C, Wright D, Barranco I, Roca J, Rodriguez-Martinez H. The Transcriptome of Pig Spermatozoa, and Its Role in Fertility. International Journal of Molecular Sciences. 2020; 21(5):1572. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms21051572

Chicago/Turabian Style

Alvarez-Rodriguez, Manuel, Cristina Martinez, Dominic Wright, Isabel Barranco, Jordi Roca, and Heriberto Rodriguez-Martinez. 2020. "The Transcriptome of Pig Spermatozoa, and Its Role in Fertility" International Journal of Molecular Sciences 21, no. 5: 1572. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms21051572

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