|Ahead of print publication
Gene transcription profiling of astheno- and normo-zoospermic sperm subpopulations
Pedro Caballero-Campo1,2, Saúl Lira-Albarrán1, David Barrera1, Elizabeth Borja-Cacho3, Héctor S Godoy-Morales3, Claudia Rangel-Escareño4, Fernando Larrea1, Mayel Chirinos1
1 Department of Reproduction Biology, National Institute of Medical Sciences and Nutrition Salvador Zubirán, Mexico City 14080, Mexico
2 Tambre Foundation, Madrid 28002, Spain
3 Reproductive Medicine Unit, Angeles del Pedregal Hospital, Mexico City 10700, Mexico
4 Computational Genomic and Integrative Biology Laboratory, National Institute of Genomic Medicine, Mexico City 14610, Mexico
|Date of Submission||27-Mar-2019|
|Date of Acceptance||11-Nov-2019|
|Date of Web Publication||14-Feb-2020|
Department of Reproduction Biology, National Institute of Medical Sciences and Nutrition Salvador Zubirán, Mexico City 14080
Department of Reproduction Biology, National Institute of Medical Sciences and Nutrition Salvador Zubirán, Mexico City 14080
Source of Support: None, Conflict of Interest: None
Spermatozoa contain a repertoire of RNAs considered to be potential functional fertility biomarkers. In this study, the gene expression of human sperm subpopulations with high (F1) and low (F2) motility from healthy normozoospermic (N) and asthenozoospermic (A) individuals was evaluated using RNA microarray followed by functional genomic analysis of differentially expressed genes. Results from A–F1 versus N–F1, A–F2 versus N–F2, N–F1 versus N–F2, and A–F1 versus A–F2 comparisons showed a considerably larger set of downregulated genes in tests versus controls. Gene ontology (GO) analysis of A–F1 versus N–F1 identified 507 overrepresented biological processes (BPs), several of which are associated with sperm physiology. In addition, gene set enrichment analysis of the same contrast showed 110 BPs, 36 cellular components, and 31 molecular functions, several of which are involved in sperm motility. A leading-edge analysis of selected GO terms resulted in several downregulated genes encoding to dyneins and kinesins, both related to sperm physiology. Furthermore, the predicted activation state of asthenozoospermia was increased, while fertility, cell movement of sperm, and gametogenesis were decreased. Interestingly, several downregulated genes characteristic of the canonical pathway protein ubiquitination were involved in asthenozoospermia activation. Conversely, GO analysis of A–F2 versus N–F2 did not identify overrepresented BPs, although the gene set enrichment analysis detected six enriched BPs, one cellular component, and two molecular functions. Overall, the results show differences in gene transcription between sperm subpopulations from asthenozoospermic and normozoospermic semen samples and allowed the identification of gene sets relevant to sperm physiology and reproduction.
Keywords: asthenozoospermia; male infertility; microarray; sperm; transcriptome
Article in PDF
|How to cite this URL:|
Caballero-Campo P, Lira-Albarrán S, Barrera D, Borja-Cacho E, Godoy-Morales HS, Rangel-Escareño C, Larrea F, Chirinos M. Gene transcription profiling of astheno- and normo-zoospermic sperm subpopulations. Asian J Androl [Epub ahead of print] [cited 2020 Mar 28]. Available from: http://www.ajandrology.com/preprintarticle.asp?id=278424
Pedro Caballero-Campo, Saúl Lira-Albarrán
These authors contributed equally to this work.
| Introduction|| |
Spermatogenesis and sperm maturation are multifactorial processes that result in the production of differentiated spermatozoa with the ability to become capacitated and fertilize. In andrology laboratories, standard semen evaluation is the primary approach to assess sperm-fertilizing ability, by detecting anomalies in sperm concentration, motility, and morphology. However, it has a limited diagnostic power, and identification of novel markers of sperm function is increasingly demanded., To this end, several studies have characterized the presence of RNAs in mature spermatozoa,, including comparisons between infertile and fertile men. The potential utility of sperm RNA as molecular fertility markers, mainly those consistently present in samples from fertile men, has been previously evaluated., Sperm RNA is a complex mixture of different types of RNA, some of which may represent a particular RNA population with functional roles after delivery into the oocyte., Besides this, some RNAs may also have physiological roles during the epigenetic reprogramming of sperm chromatin and sperm maturation and may therefore be involved in their ability to fertilize., However, recent studies indicate that a substantial portion of sperm mRNA is fragmented, though some transcripts are maintained intact with potential roles during sperm transit through the female reproductive tract, fertilization, and early embryogenesis.,,
Asthenozoospermia, defined as a semen sample with a proportion of progressively motile spermatozoa below the lower WHO reference limit, is a major cause of male infertility, and several genes and molecular markers have been linked to this condition. This is relevant because assisted reproductive technologies require the isolation of progressively motile spermatozoa from those that are immotile and incapable of fertilizing. Indeed, only a small percentage of spermatozoa can achieve fertilization, and their isolation and full characterization represents a significant challenge in the study of male infertility. In this study, we compared the gene expression profiling of high (F1) and low (F2) sperm motility fractions isolated from the semen samples of asthenozoospermic and normozoospermic men by microarray technology. The aim was to identify transcripts that differ between sperm subpopulations and to determine their functional relevance in sperm physiology and reproduction.
| Materials and Methods|| |
Semen samples and sperm analysis
This study was approved by the Human Ethical and Scientific Review Committees of the National Institute of Medical Sciences and Nutrition Salvador Zubirán, Mexico City, Mexico (reference number 1516, from 04/27/2015), and all participants gave written informed consent for the use of their semen in research. Human semen samples obtained by masturbation from normozoospermic (N) and asthenozoospermic (A) volunteers were collected in sterile containers and evaluated within 1 h after ejaculation according to the WHO reference guidelines. To isolate sperm subpopulations, liquefied semen samples were centrifuged at 800 g for 30 min through a 90%/60% discontinuous density gradient (Isolate; Irvine Scientific, Irvine, CA, USA). For each sample, two sperm subpopulations were recovered and named as A–F1/N–F1 (sperm pellets) and A–F2/N–F2 (90%/60% interphases). Sperm fractions were washed with phosphate-buffered saline; re-assessed for cell concentration, motility, and purity (visual absence of somatic cell contamination); and processed for RNA extractions. For analysis by microarray technology, F1 and F2 samples recovered from three A and three N semen donors were grouped into four study groups (three A–F1, three A–F2, three N–F1, and three N–F2), for a total of 12 microarrays.
RNA isolation and microarray hybridization
Total RNA from sperm subpopulations was extracted as previously described., Briefly, sperm samples were suspended in 1 ml of TRIzol® reagent (Invitrogen, Carlsbad, CA, USA), passed through 20G hypodermic needles, and separated by the addition of 0.2-ml chloroform. The RNA concentration was determined by using an ND-1000 Spectrophotometer (NanoDrop Technologies, Wilmington, DE, USA), and its quality was checked using an Agilent Bioanalyzer 2100 (Agilent Technologies Inc., Palo Alto, CA, USA) by loading samples onto RNA nano chips included in the RNA 6000 Nano kit (Agilent Technologies) for analysis with the Eukaryote Total RNA Nano Assay (version 2.6) and following sperm RNA quality criteria previously described. To confirm the absence of contamination by somatic cell RNA, reverse transcription polymerase chain reaction (RT-PCR) for CD45, c-kit, and E-cadherin was performed with primers indicated in [Supplementary Table 1 [Additional file 1]], and the resulting PCR amplicons were analyzed by gel electrophoresis.,
For microarray experiments, target cDNAs were prepared according to the Whole-Transcript PLUS (WT) Sense Target Labeling Protocol (Affymetrix, Inc., Santa Clara, CA, USA) as previously described., The cRNA products used as templates for second cycles of cDNA synthesis in the presence of dUTPs were fragmented by uracil-DNA glycosidase and apurin apirymidin endonuclease. The fragments (40–70 mer) were then biotin labeled with deoxynucleotide terminal addition reaction. Labeled single-stranded cDNAs were hybridized during 17 h at 45°C onto a GeneChip® Human Gene 2.0 ST Array (Affymetrix), containing 40 716 total RefSeq transcripts for 30 654 NM RefSeq well-established annotations for the measurement of protein coding and long intergenic noncoding RNA transcripts. Washing and staining with streptavidin-phycoerythrin were performed with the Affymetrix Fluidics Station 450 using the Affymetrix Staining Kit. Finally, the microarrays were scanned for fluorescence signals in a GeneChip Scanner 3000 7G (Affymetrix).
Microarray raw data were background corrected by the Robust Multiarray Average method and normalized by quantile normalization. The analyses involved four contrasts (A–F1 vs N–F1, A–F2 vs N–F2, N–F1 vs N–F2, and A–F1 vs A–F2), and the differentially expressed genes (DEGs) were determined by linear statistical models with arbitrary coefficients from the Bioconductor library limma., Multiple hypotheses were corrected by applying a false discovery rate (FDR). The up- and down-DEGs were selected on the bases of a fold-change |FC| >2 and P < 0.01. All analyses were performed by R software (version 3.2.1, 2015; The R Foundation for Statistical Computing, Vienna, Austria; https://www.r-project.org/).
Gene ontology analysis
Gene ontology (GO) was analyzed by Babelomics 5.0's FatiGO (http://babelomics.bioinfo.cipf.es) from GO terms generated by the list of DEGs, as previously described., The nomenclature of the enriched biological processes (BPs) included terms of the Gene Ontology Consortium. The adjusted P-values were calculated by an exact Fisher's test that evaluated the significant overrepresentation of functional terms in the list of DEGs with respect to the rest of the human genome after correcting for multiple tests (multiple hypotheses, one for each functional term) by FDR. A REVIGO (reduce + visualize Gene Ontology; http://revigo.irb.hr) analysis was done to identify and avoid redundancy of the most representative BPs, as well as for the construction of TreeMaps. Moreover, gene set enrichment analysis (GSEA; http://www.broadinstitute.org/gsea) was performed with an FDR <0.05 and normalized enrichment scores (NES) ≥2.0, and a leading-edge analysis was done to identify common gene subsets within the selected GO terms.
Ingenuity pathway analysis
Selected bio-functions associated with sperm physiology and reproduction were analyzed by using ingenuity pathway analysis (IPA; https://www.qiagenbioinformatics.com/products/ingenuity-pathway-analysis/). Bio-functions were considered relevant when an absolute z-score >2 and P < 0.01 were fulfilled. Furthermore, from a hypothesis-driven approach to the DEGs in A–F1 versus N–F1 comparison on sperm physiology bio-functions, a molecule activity prediction (MAP) by IPA was analyzed as previously described.,,
Quantitative real-time PCR
For microarray validation of DEGs, 1-μg total RNA from each sperm sample was reversely transcribed by using the Transcriptor First Strand cDNA Synthesis Kit (Roche, Mannheim, Germany) according to the manufacturer's instructions. Quantitative real-time PCR (qPCR) was performed on a Light Cycler® 2.0 Detection System with a TaqMan PCR Master Mix (Applied Biosystems, Carlsbad, CA, USA). Primers and probes for PCR amplifications were designed with the online software from Universal Probe Library Assay Design Center of Roche (https://lifescience.roche.com/en_mx/brands/universal-probe-library.html#assay-design-center), and each pair of primers was submitted to the National Center for Biotechnology Information BLAST search to ensure specificity for the target mRNA. The oligonucleotide sequences used for gene amplification are shown in [Supplementary Table 1]. The comparative CT method was used to quantify the expression of target genes after normalizing against ribosomal protein L32 (RPL32)as the housekeeping gene.,
NCSS software (NCSS, Kaysville, UT, USA) was used for data analysis. From the homogeneity and distribution of the data, either Student's t- test or nonparametric Mann–Whitney U test was performed to compare differences between the A and N control samples. P ≤ 0.05 was considered statistically significant.
| Results|| |
The analysis of A semen samples showed significantly lower progressive motility than N samples. Moreover, the comparison of sperm subpopulations showed a lower percentage of progressive motility in A–F1 than N–F1, but there was no statistical difference between A–F2 and N–F2 [Supplementary Table 2 [Additional file 2]].
RNA sample quality
Sperm RNA samples' integrity and purity were evaluated by microelectropherograms, with controls of RNA from peripheral blood mononuclear cells (PBMC). RNA isolated from PBMC showed two peaks corresponding to 18S and 28S rRNA that were absent from sperm RNA [Supplementary Figure 1 [Additional file 3]]a. Moreover, no amplicons were observed when the RNA from sperm subpopulations was analyzed for specific somatic RNA markers [Supplementary Figure 1]b.
Differential gene expression in sperm subpopulations and GO analysis
The unsupervised hierarchical clustering of A–F1/N–F1 [Figure 1]a and A–F2/N–F2 [Figure 1]b revealed two sample groups in each contrast. Analysis of A–F1 versus N–F1 resulted in 1863 DEGs, comprising 116 (6.2%) up- and 1747 (93.8%) down-regulated genes. Conversely, the A–F2 versus N–F2 contrast identified 207 DEGs, including 17 (8.2%) up- and 190 (91.8%) down-regulated genes. However, intragroup subpopulation comparisons analyzed as N–F1 versus N–F2 and A–F1 versus A–F2 resulted in 61 (7 up- and 54 down-regulated genes) and 908 DEGs (42 up- and 866 down-regulated genes), respectively. The corresponding heatmaps are shown in [Supplementary Figure 2 [Additional file 4]]. Moreover, the top up- and down-regulated genes from each of the four contrasts are summarized in [Supplementary Table 3 [Additional file 5]], [Supplementary Table 4 [Additional file 6]], [Supplementary Table 5 [Additional file 7]], [Supplementary Table 6 [Additional file 8]]. The directionality of three down- and one up-regulated DEG from A–F1 versus N–F1 was confirmed by qPCR [Supplementary Table 7 [Additional file 9]].
|Figure 1: Gene clustering of the microarray data showing pair-wise comparisons. (a) A–F1 versus N–F1; (b) A–F2 versus N–F2. Heatmaps show one sample for each column and one gene for each horizontal line. Color indicates gene expression value intensities (z-score), where red signifies upregulation, green downregulation, and black unchanged. A–F1: asthenozoospermic sperm subpopulation with high motility; N–F1: normozoospermic sperm subpopulation with high motility; A-F2: asthenozoospermic sperm subpopulation with low motility; N-F2: normozoospermic sperm subpopulation with low motility.|
Click here to view
GO analysis identified 507 BPs overrepresented in A–F1 versus N–F1. The ten BPs with the highest statistical significance (log10P > -9.9) are shown in [Table 1]. BPs involved in reproduction, spermatid differentiation, and cell movement were significantly overrepresented. Most genes within each BP were downregulated. A REVIGO TreeMap was constructed to summarize the findings of the GO analysis [Supplementary Figure 3 [Additional file 10]]. On the other hand, no BP was overrepresented in the A–F2 versus N–F2 comparison.
|Table 1: Top ten biological processes identified from 507 gene ontology terms summarized by REVIGO in A–F1 versus N–F1|
Click here to view
Gene set enrichment analysis
Leading-edge analysis of the 110 BPs identified in the A-F1 versusN-F1 comparison showed a set of genes, including dynein axonemal intermediate chain 1 (DNAI1), dynein axonemal heavy chain 1 (DNAH1) and 7 (DNAH7), leucine-rich repeat containing 6 (LRRC6), tetratricopeptide repeat domain 26 (TTC26), sperm-associated antigen 16 (SPAG16), CF transmembrane conductance regulator (CFTR), ubiquitin conjugating enzyme E2 B (UBE2B), and zona pellucida binding protein 2 (ZPBP2), that overlapped with several physiologically relevant BPs (with NES ≥2 and FDR <0.05) [Figure 2]a. Most of the enriched BPs identified were involved in spermatogenesis and sperm physiology, including fertilization, sperm motility, axonemal dynein complex assembly, male gamete generation, cilium morphogenesis, and microtubule-based movement, with the last two showing the highest NES [Supplementary Table 8 [Additional file 11]]. Furthermore, 36 cellular components (CC) were identified in the same contrast [Supplementary Table 9 [Additional file 12]], which highlights a subset of genes comprising the kinesin family members 2A (KIF2A), 2B (KIF2B), and 2C (KIF2C); intraflagellar transport 81 (IFT81); NME/NM23 family member 8 (TXNDC3); DNAH8, DNAH1, and DNAH7;A-kinase anchoring protein 3 (AKAP3) and 4 (AKAP4) involved in motile cilium, dynein complex, axonemal part, microtubule-associated complex, sperm flagellum, and kinesin complex [Figure 2]b. Finally, 31 molecular functions (MFs) were also identified [Supplementary Table 10 [Additional file 13]], where the kinesin family members 3B (KIF3B), C3 (KIFC3), 16 (KIF16), and 20 (KIF20), KIF2A, KIF2C,KIF2B, DNAH1, DNAH7, DNAH8, DNAH6, and DNAH9 genes overlapped with the MFs' ATPase activity signal and microtubule motor activity. Moreover, the genes KIF16, KIF3B, KIFC3, and KIF20 also overlapped with the MF ATP-dependent microtubule motor activity [Figure 2]c.
|Figure 2: Leading-edge analysis of the A–F1 versus N–F1. (a) Biological processes; (b) cellular components; (c) molecular functions. Columns: genes within the core enrichment of the shown GO terms; rows: selected GO terms from the gene set enrichment analysis (Supplementary Table 8–10). Gene expression values are represented in colors, and the range of colors (from light to dark blue) shows the level of expression values (from low to lowest) (http://www.broadinstitute.org/gsea). GO: gene ontology; A–F1: asthenozoospermic sperm subpopulation with high motility; N–F1: normozoospermic sperm subpopulation with high motility.|
Click here to view
The core enrichment of microtubule-based movement showed that genes KIF2A, KIF2B, KIF2C, DNAH1, DNAH7, and DNAH8 [Figure 3] were also identified in the main CCs and MFs [Figure 2]b and [Figure 2]c. Conversely, A–F2 versus N–F2 showed the enrichment of only six BPs (highlighting cilium morphogenesis and cilium organization; [Supplementary Table 11 [Additional file 14]], one CC (cilium), and two MFs (ubiquitin-like protein-specific protease activity and RAN GTPase binding).
|Figure 3: Heatmap of the core enrichment of the biological process microtubule-based movement. Columns correspond to sperm samples analyzed, whereas lines represent genes within the core enrichment. Arrows indicate genes downregulated in A–F1 that were also identified in the leading-edge analysis of biological processes, cellular components, and molecular functions, as shown in Figure 2. Color indicates gene expression value intensities (z-score), where red signifies upregulation, green downregulation, and black unchanged. A–F1: asthenozoospermic sperm subpopulation with high motility; N–F1: normozoospermic sperm subpopulation with high motility.|
Click here to view
Predicted activation state of bio-functions
IPA core analysis of the list of DEGs from the A–F1 versus N–F1 comparison showed that the predicted activation state of several bio-functions involved in reproduction, either increased or decreased [Table 2]. Sperm disorders, asthenozoospermia, oligozoospermia, and infertility were predicted to be increased. In addition, fertility was predicted to be decreased, along with microtubule dynamics, organization of cytoskeleton, formation of cilia, and cell movement of sperm.
|Table 2: Bio-functions predicted activation state in the A–F1 versus N–F1 comparison by ingenuity pathway analysis core analysis|
Click here to view
Molecule activity prediction of bio-functions
The MAP analysis by IPA of DEGs from the A–F1 versus N–F1 comparison predicted bio-function fertility to be inhibited, on the basis of downregulation of 13 genes where Angiotensin-converting enzyme (ACE) and ADAM Metallopeptidase Domain 2 (ADAM2) are involved in the canonical pathway (CP) axonal guidance signaling [Figure 4]a. Similarly, MAP analyses of the bio-function capacitation, acrosome reaction, and binding of sperm were also predicted to be inhibited [Supplementary Figure 4 [Additional file 15]]. Conversely, asthenozoospermia was predicted to be activated on the basis of the downregulation of 21 genes including Ubiquitin Conjugating Enzyme E2 J1 (UBE2J1), Heat Shock Protein Family A 4 Like (HSPA4L), and DnaJ Heat Shock Protein Family (Hsp40) A1 (DNAJA1), which areall involved in the CP protein ubiquitination pathway [Figure 4]b.
|Figure 4: Molecular activity prediction by IPA from the DEGs in A–F1 versus N–F1 contrast. (a) Fertility (inhibited); (b) asthenozoospermia (activated). Colors indicate predicted relationships of gene expression levels and bio-functions, and color intensities reflect the degree of gene expression or bio-function activity. IPA: ingenuity pathway analysis (https://www.qiagenbioinformatics. com/products/ingenuity-pathway-analysis/); DEGs: differentially expressed genes; A–F1: asthenozoospermic sperm subpopulation with high motility; N–F1: normozoospermic sperm subpopulation with high motility; CP: canonical pathway.|
Click here to view
| Discussion|| |
Approximately 13% of couples worldwide seek medical advice for fertility purposes. Unexplained infertility is commonly treated with expensive and possibly unneeded high-complexity technologies. The recurrent use of in-vitro fertilization and intracytoplasmic sperm injection in assisted reproductive technologies is a consequence of the inadequacy of current male fertility diagnostic methods., Therefore, new approaches to studying spermatozoa at the cellular and molecular level are still required. In this regard, the characterization of sperm RNA profiling has become a promising strategy to reveal the mechanisms leading to altered sperm parameter values. The study of sperm physiology using high-throughput measurement technology such as microarray analysis has previously revealed several biomarkers of human sperm function., This approach may help to increase our understanding of the basic molecular mechanisms underlying the pathogenesis of male infertility.
In natural conception, a small fraction of the spermatozoa in an ejaculate reaches the fertilization site. Sperm selection takes places in the female genital tract, where sperm motility is a key feature to reach the oocyte. Hence, analysis of sperm subpopulations based on differential motility may represent a valuable model for studying sperm physiology and reproduction. In 2012, Jodar et al. performed the first study of human sperm transcriptome of asthenozoospermic semen samples using microarray analysis and showed a differential expression pattern between the F1 of patients and controls. We also found a transcriptomic profile in A–F1 different from that of N–F1 samples, but decided further to compare sperm subpopulations with different motilities obtained from density gradients. As expected, F2 subpopulations presented lower motility values than the corresponding F1 from A and N samples, but there were considerably fewer DEGs in A–F2 versus N–F2 than A–F1 versus N–F1, which indicates that F2 from A and N samples had similar expression profiles. Moreover, these data indicate that genomic transcriptional differences between the A and N samples reside mainly in the highly motile sperm subpopulation. However, when the internal contrasts were analyzed, N–F1 versus N–F2 showed considerably fewer differences in DEGs than A–F1 versus A–F2, which suggests that normozoospermic sperm subpopulations are more homogeneous than asthenozoospermic. Furthermore, overlapping of DEGs from N–F1 versus N–F2 and A–F1 versus A–F2 evidenced the downregulation of nicotinamide phosphoribosyltransferase, serglycin, and lysosomal protein transmembrane 5, but these transcripts did not show differential expression in A–F1 versus N–F1 or A–F2 versus N–F2, which suggests that their overexpression in F2 fractions may be associated with sperm's reduced motility. Similarly, both A–F1 versus N–F1 and A–F2 versus N–F2 contrasts identified the downregulation in A fractions of apolipoprotein B mRNA editing enzyme catalytic polypeptide-like 4, glyceronephosphate O-acyltransferase, and long intergenic nonprotein coding RNA 1364, indicating that the underexpression of these three transcripts could be related to asthenozoospermia and therefore may be used as potential candidates for future investigations.
As the A–F1 versus N–F1 contrast resulted in a higher number of DEGs, a deeper analysis showed that among the more upregulated RNAs in A–F1, several of the mitochondrial origin (including mitochondrially encoded nicotinamide adenine dinucleotide reduced [NADH]:ubiquinone oxidoreductase core subunits) and microRNAs were highlighted. The differential expression of microRNAs in asthenozoospermia has been described previously,, and upregulation of microRNA 27b seems to be specifically involved in reduced cysteine-rich secretory protein 2 expression in asthenozoospermia. All of these indicate that microRNAs may play an important role in the manifestations of this pathology.
As expected, microtubule-based movement was among the most overrepresented BP, and a core enrichment analysis of this BP revealed several genes involved in sperm motility that were also identified by the leading-edge GSEA. This analysis also identified KIF2A, KIF2B, KIF2C, DNAH1, DNAH7, and DNAH8 as consistently downregulated in A–F1, all of them being members of the axonemal dynein and kinesin families, which play fundamental roles in the mammalian sperm flagellum. Previous investigations indicate that other members of the dynein and kinesin gene families are also differentially expressed in asthenozoospermia.DNAH1 mutations have formerly been associated with sperm defects and infertility.,, Furthermore, naturally occurring mutations in some other genes identified in the present investigation, such as CFTR and LRRC6, have also been associated with several pathologies that involve morphological anomalies of the flagellum, resulting in asthenozoospermia and male infertility.,
Protein polyubiquitination was identified as a top BP, and the canonical pathway protein ubiquitination was also associated with asthenozoospermia, which suggests that proper regulation of protein degradation in mature spermatozoa is fundamental for obtaining healthy normozoospermic spermatozoa. Indeed, previous investigations have demonstrated that increased sperm ubiquitination is inversely associated with sperm motility and that there is a positive correlation between ubiquitinated spermatozoa and the percentage of spermatozoa with abnormal chromatin.
In addition, data mining analysis of the A–F1 versus N–F1 comparison identified specific GO terms associated with activation of the bio-function asthenozoospermia and with the inhibition of fertility, acrosome reaction, capacitation, and sperm binding. Previously, GO analyses of different pathological sperm profiles have shown several GO terms consistent with alterations in morphology, motility, and sperm count,,, that were also detected in the present study, such as spermatid differentiation and development and repair of DNA. This suggests that asthenozoospermia has common mechanisms with some sperm bio-functions essential for fertilization, and should therefore be seen as a complex pathology not only involving a compromised cell motility, but also other cellular malfunctions, jeopardizing the proper delivery of DNA/RNA cargo to the oocyte. This hypothesis may have implications for the management of asthenozoospermic patients in reproductive assisted facilities. However, whether the RNAs herein identified are of relevance in sperm capacitation, cell fusion, or early stages of embryo development requires further investigation.
Finally, owing to the high level of downregulated mRNAs and upregulated microRNAs observed in spermatozoa from asthenozoospermic patients, we suggest that an environmental origin via epigenetic genomic markers could be considered among the causes of semen abnormalities. In conclusion, our analysis of sperm RNA fractions from asthenozoospermic patients identified a number of gene sets significantly associated with male infertility due to low sperm motility. These results provide new ideas for further studies on male infertility diagnosis and treatment and for the development of strategies for research and development in male contraception.
| Author Contributions|| |
PCC, SLA, FL, and MC designed the experimental study. PCC, EBC, and HSGM selected and processed the sperm samples. PCC, MC, and DB acquired and interpreted experimental results. SLA and CRE performed the bioinformatics analyses. PCC, SLA, CRE, FL, and MC analyzed the data and wrote the manuscript. All authors read and approved the final manuscript.
| Competing Interests|| |
All authors declare no competing interests.
| Acknowledgments|| |
We thank Elisabeth Raab MD, MPH (University of Southern California, Los Angeles, CA, USA), and Paolo Rinaudo MD, PhD (University of California San Francisco, San Francisco, CA, USA), for the critical revision of the manuscript. Financial support was provided by a grant from Cathedra Salvador Zubirán (National Autonomous University of Mexico and National Institute of Medical Sciences and Nutrition Salvador Zubirán, Mexico City, Mexico) and funds from the Department of Reproduction Biology of the National Institute of Medical Sciences and Nutrition Salvador Zubirán, Patronage of the National Institute of Medical Sciences and Nutrition Salvador Zubirán, FunSaEd A.C. (Mexico), and Tambre Foundation (Madrid, Spain).
Supplementary Information is linked to the online version of the paper on the Asian Journal of Andrology website.
| References|| |
Sakkas D, Ramalingam M, Garrido N, Barratt CL. Sperm selection in natural conception: what can we learn from mother nature to improve assisted reproduction outcomes? Hum Reprod Update
2015; 21: 711–26.
Sousa AP, Amaral A, Baptista M, Tavares R, Caballero Campo P, et al.
Not all sperm are equal: functional mitochondria characterize a subpopulation of human sperm with better fertilization potential. PLoS One
2011; 6: e18112.
Jodar M, Selvaraju S, Sendler E, Diamond MP, Krawetz SA, et al.
The presence, role and clinical use of spermatozoal RNAs. Hum Reprod Update
2013; 19: 604–24.
Garrido N, Garcia-Herrero S, Meseguer M. Assessment of sperm using mRNA microarray technology. Fertil Steril
2013; 99: 1008–22.
Bansal SK, Gupta N, Sankhwar SN, Rajender S. Differential genes expression between fertile and infertile spermatozoa revealed by transcriptome analysis. PLoS One
2015; 10: e0127007.
Bonache S, Mata A, Ramos MD, Bassas L, Larriba S. Sperm gene expression profile is related to pregnancy rate after insemination and is predictive of low fecundity in normozoospermic men. Hum Reprod
2012; 27: 1556–67.
Ostermeier GC, Dix DJ, Miller D, Khatri P, Krawetz SA. Spermatozoal RNA profiles of normal fertile men. Lancet
2002; 360: 772–7.
Boerke A, Dieleman SJ, Gadella BM. A possible role for sperm RNA in early embryo development. Theriogenology
2007; 68: S147–55.
Hamatani T. Human spermatozoal RNAs. Fertil Steril
2012; 97: 275–81.
Li Y, Li RH, Ran MX, Zhang Y, Liang K, et al.
High throughput small RNA and transcriptome sequencing reveal capacitation-related microRNAs and mRNA in boar sperm. BMC Genomics
2018; 19: 736.
Johnson GD, Mackie P, Jodar M, Moskovtsev S, Krawetz SA. Chromatin and extracellular vesicle associated sperm RNAs. Nucleic Acids Res
2015; 43: 6847–59.
Selvaraju S, Parthipan S, Somashekar L, Kolte AP, Binsila BK, et al.
Occurrence and functional significance of the transcriptome in bovine (Bos taurus
) spermatozoa. Sci Rep
2017; 7: 42392.
Sendler E, Johnson GD, Mao S, Goodrich RJ, Diamond MP, et al.
Stability, delivery and functions of human sperm RNAs at fertilization. Nucleic Acids Res
2013; 41: 4104–17.
World Health Organization. WHO Laboratory Manual for the Examination and Processing of Human Semen. Geneva: World Health Organization Press; 2010.
Jodar M, Kalko S, Castillo J, Ballesca JL, Oliva R. Differential RNAs in the sperm cells of asthenozoospermic patients. Hum Reprod
2012; 27: 1431–8.
Henkel R. Sperm preparation: state-of-the-art-physiological aspects and application of advanced sperm preparation methods. Asian J Androl
2012; 14: 260–9.
Caballero-Campo P, Lin W, Simbulan R, Liu X, Feuer S, et al.
Advanced paternal age affects sperm count and anogenital distance in mouse offspring. Reprod Sci
2018; 25: 515–22.
Lambard S, Galeraud-Denis I, Martin G, Levy R, Chocat A, et al.
Analysis and significance of mRNA in human ejaculated sperm from normozoospermic donors: relationship to sperm motility and capacitation. Mol Hum Reprod
2004; 10: 535–41.
Lira-Albarran S, Durand M, Larrea-Schiavon MF, Gonzalez L, Barrera D, et al.
Ulipristal acetate administration at mid-cycle changes gene expression profiling of endometrial biopsies taken during the receptive period of the human menstrual cycle. Mol Cell Endocrinol
2017; 447: 1–11.
Lira-Albarran S, Larrea-Schiavon MF, Gonzalez L, Durand M, Rangel C, et al.
The effects of levonorgestrel on FSH-stimulated primary rat granulosa cell cultures through gene expression profiling are associated to hormone and folliculogenesis processes. Mol Cell Endocrinol
2017; 439: 337–45.
Irizarry RA, Bolstad BM, Collin F, Cope LM, Hobbs B, et al.
Summaries of Affymetrix GeneChip probe level data. Nucleic Acids Res
2003; 31: e15.
Bolstad BM, Irizarry RA, Astrand M, Speed TP. A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics
2003; 19: 185–93.
Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, et al.
Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res
2015; 43: e47.
Smyth GK. Linear models and empirical Bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol
2004; 3: 1–25.
Al-Shahrour F, Diaz-Uriarte R, Dopazo J. FatiGO: a web tool for finding significant associations of gene ontology terms with groups of genes. Bioinformatics
2004; 20: 578–80.
Alonso R, Salavert F, Garcia-Garcia F, Carbonell-Caballero J, Bleda M, et al.
Babelomics 5.0: functional interpretation for new generations of genomic data. Nucleic Acids Res
2015; 43: W117–21.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, et al.
Gene ontology: tool for the unification of biology. Nat Genet
2000; 25: 25–9.
Colquhoun D. An investigation of the false discovery rate and the misinterpretation of p-values. R Soc Open Sci
2014; 1: 140216.
Supek F, Bosnjak M, Skunca N, Smuc T. REVIGO summarizes and visualizes long lists of gene ontology terms. PLoS One
2011; 6: e21800.
Lira-Albarran S, Durand M, Barrera D, Vega C, Becerra RG, et al.
A single preovulatory administration of ulipristal acetate affects the decidualization process of the human endometrium during the receptive period of the menstrual cycle. Mol Cell Endocrinol
2018; 476: 70–8.
Lira-Albarran S, Vega CC, Durand M, Rangel C, Larrea F. Functional genomic analysis of the human receptive endometrium transcriptome upon administration of mifepristone at the time of follicle rupture. Mol Cell Endocrinol
2019; 485: 88–96.
Schefe JH, Lehmann KE, Buschmann IR, Unger T, Funke-Kaiser H. Quantitative real-time RT-PCR data analysis: current concepts and the novel “gene expression's CT difference” formula. J Mol Med (Berl)
2006; 84: 901–10.
Barragan M, Martinez A, Llonch S, Pujol A, Vernaeve V, et al.
Effect of ribonucleic acid (RNA) isolation methods on putative reference genes messenger RNA abundance in human spermatozoa. Andrology
2015; 3: 797–804.
Perez-Rico A, Crespo F, Sanmartin ML, De Santiago A, Vega-Pla JL. Determining ACTB
as optimal reference genes for quantitative RT-PCR studies of cryopreserved stallion semen. Anim Reprod Sci
2014; 149: 204–11.
Feuer S, Rinaudo P. From embryos to adults: a DOHaD perspective onin vitro
fertilization and other assisted reproductive technologies. Healthcare (Basel)
2016; 4: 51.
Dyer S, Chambers GM, de Mouzon J, Nygren KG, Zegers-Hochschild F, et al.
International committee for monitoring assisted reproductive technologies world report: assisted reproductive technology 2008, 2009 and 2010. Hum Reprod
2016; 31: 1588–609.
Zhang T, Wu J, Liao C, Ni Z, Zheng J, et al.
System analysis of teratozoospermia mRNA profile based on integrated bioinformatics tools. Mol Med Rep
2018; 18: 1297–304.
Zhang X, Zhang P, Song D, Xiong S, Zhang H, et al.
Expression profiles and characteristics of human lncRNA in normal and asthenozoospermia sperm. Biol Reprod
2018; 100: 982–93.
Abu-Halima M, Hammadeh M, Schmitt J, Leidinger P, Keller A, et al.
Altered microRNA expression profiles of human spermatozoa in patients with different spermatogenic impairments. Fertil Steril
2013; 99: 1249–55.e16.
Salas-Huetos A, Blanco J, Vidal F, Godo A, Grossmann M, et al.
Spermatozoa from patients with seminal alterations exhibit a differential micro-ribonucleic acid profile. Fertil Steril
2015; 104: 591–601.
Zhou JH, Zhou QZ, Lyu XM, Zhu T, Chen ZJ, et al.
The expression of cysteine-rich secretory protein 2 (CRISP2) and its specific regulator miR-27b in the spermatozoa of patients with asthenozoospermia. Biol Reprod
2015; 92: 28.
Lindemann CB, Lesich KA. Functional anatomy of the mammalian sperm flagellum. Cytoskeleton
2016; 73: 652–69.
Amiri-Yekta A, Coutton C, Kherraf ZE, Karaouzene T, Le Tanno P, et al.
Whole-exome sequencing of familial cases of multiple morphological abnormalities of the sperm flagella (MMAF) reveals new DNAH1
mutations. Hum Reprod
2016; 31: 2872–80.
Sha Y, Yang X, Mei L, Ji Z, Wang X, et al.DNAH1
gene mutations and their potential association with dysplasia of the sperm fibrous sheath and infertility in the Han Chinese population. Fertil Steril
2017; 107: 1312–8.e2.
Ben Khelifa M, Coutton C, Zouari R, Karaouzene T, Rendu J, et al.
Mutations in DNAH1
, which encodes an inner arm heavy chain dynein, lead to male infertility from multiple morphological abnormalities of the sperm flagella. Am J Hum Genet
2014; 94: 95–104.
Li CY, Jiang LY, Chen WY, Li K, Sheng HQ, et al.
CFTR is essential for sperm fertilizing capacity and is correlated with sperm quality in humans. Hum Reprod
2010; 25: 317–27.
Kott E, Duquesnoy P, Copin B, Legendre M, Dastot-Le Moal F, et al.
Loss-of-function mutations in LRRC6
, a gene essential for proper axonemal assembly of inner and outer dynein arms, cause primary ciliary dyskinesia. Am J Hum Genet
2012; 91: 958–64.
Sutovsky P, Hauser R, Sutovsky M. Increased levels of sperm ubiquitin correlate with semen quality in men from an andrology laboratory clinic population. Hum Reprod
2004; 19: 628–38.
Hodjat M, Akhondi MA, Al-Hasani S, Mobaraki M, Sadeghi MR. Increased sperm ubiquitination correlates with abnormal chromatin integrity. Reprod Biomed Online
2008; 17: 324–30.
Fu G, Wei Y, Wang X, Yu L. Identification of candidate causal genes and their associated pathogenic mechanisms underlying teratozoospermia based on the spermatozoa transcript profiles. Andrologia
2016; 48: 576–83.
Liu XX, Cai L, Liu FJ. An in silico
analysis of human sperm genes associated with asthenozoospermia and its implication in male infertility. Medicine (Baltimore)
2018; 97: e13338.
Montjean D, De La Grange P, Gentien D, Rapinat A, Belloc S, et al.
Sperm transcriptome profiling in oligozoospermia. J Assist Reprod Genet
2012; 29: 3–10.
[Figure 1], [Figure 2], [Figure 3], [Figure 4]
[Table 1], [Table 2]