|Ahead of print publication
Proteomic analysis reveals dysregulated cell signaling in ejaculated spermatozoa from infertile men
Luna Samanta1,2, Rakesh Sharma1, Zhihong Cui1,3, Ashok Agarwal1
1 American Center for Reproductive Medicine, Cleveland Clinic, Cleveland, OH 44195, USA
2 Department of Zoology, School of Life Sciences, Ravenshaw University, Cuttack, Odisha 753003, India
3 Institute of Toxicology, The Third Military Medical University, Chongqing 400038, China
|Date of Submission||20-Dec-2017|
|Date of Acceptance||13-Jun-2018|
|Date of Web Publication||26-Oct-2018|
American Center for Reproductive Medicine, Cleveland Clinic, Cleveland, OH 44195, USA
Source of Support: None, Conflict of Interest: None
Dysfunctional sperm maturation is the primary reason for the poor sperm motility and morphology in infertile men. Spermatozoa from infertile men were fractioned on three-layer density gradient (80%, 60%, and 40%). Fraction 1 (F1) refers to the least mature stage having the lowest density, whereas the fraction 4 (F4) includes the most dense and morphologically mature motile spermatozoa. Fraction 2 (F2) and fraction 3 (F3) represent the intermediate stages. Proteins were extracted and separated by 1-dimensional gel. Bands were digested with trypsin and analyzed on a LTQ-Orbitrap Elite hybrid mass spectrometer system. Functional annotations of proteins were obtained using bioinformatics tools and pathway databases. A total of 1585 proteins were detected in the four fractions of spermatozoa. A dysregulated protein turnover and protein folding may lead to accumulation of defective proteins or proteins that otherwise would have been eliminated during the process of maturation, resulting in the impairment of sperm function. Aberrant chaperone expression may be a major contributing factor to the defective sperm function. Androgen receptor was predicted as a transcription regulator in one of the networks and the affected pathways were chaperone-mediated stress response, proteosomal pathway, and sperm function. The downregulation of key pathways and proteins which compromises the fertilizing potential of spermatozoa may provide insight into the mechanisms that lead to male infertility.
Keywords: androgen receptor; chaperone; immature sperm; infertile men; proteasome; spermatozoa
Article in PDF
|How to cite this URL:|
Samanta L, Sharma R, Cui Z, Agarwal A. Proteomic analysis reveals dysregulated cell signaling in ejaculated spermatozoa from infertile men. Asian J Androl [Epub ahead of print] [cited 2019 Jan 23]. Available from: http://www.ajandrology.com/preprintarticle.asp?id=244246
| Introduction|| |
Male infertility is a multifactorial condition and there is no identifiable cause in 50% of the cases. The human testis produces spermatozoa at a rate of 1000 cells per second and these cells are highly differentiated and unique. Spermatozoa originate from the complex process of spermatogenesis in three major steps as follows: (1) proliferation and differentiation of spermatogonia; (2) divisions during the spermatocyte stage; and (3) spermiogenesis. Spermiogenesis involves major morphological and molecular changes, including the removal of cytoplasm, formation of the acrosome and flagella, mitochondrial rearrangement, and nuclear remodeling. During mid-spermiogenesis, the nucleus of the round spermatid changes from spherical to a unique elongated and flattened shape. This reshaping protects the male genome during sperm transport and also facilitates the penetration of spermatozoa into ovum. Thus, spermatozoa are terminally differentiated and possess specialized organelles. However, they undergo maturation during epididymal transit to acquire the ability to fertilize.
Spermatozoa have different pathologies from those of somatic cells, which result in different sperm phenotypes in the ejaculated semen. Seven sperm phenotypes have been detected in human semen from electron microscopy, which include spermatozoa with dysplasia of the fibrous sheath, nonspecific flagellar defects, immotile cilia, acrosomal hypoplasia, defective chromatin condensation and compaction, pin head, and even sperm cells without heads. These conditions cannot be identified by routine semen analysis or functional tests since the deficiencies demonstrated by these methods are secondary manifestations of an underlying pathology.
The generation of high-quality spermatozoa is governed by a number of selective mechanisms within the testes and epididymis. The dramatic changes that occur during spermiogenesis, sperm maturation, and capacitation involve loss and gain of specific proteins.,, Recently, we have reported that distinct proteomic signatures distinguish high-quality spermatozoa from their low-quality counterparts in fertile donors. Improper spermatogenesis produces abnormal spermatozoa that are generally earmarked for elimination by apoptosis and appear in the ejaculate when they escape apoptosis. Furthermore, differential localization of Fas, a membrane receptor of the tumor necrosis factor family that initiates apoptosis, also segregates the spermatozoa into different subsets., Thus, spermatozoa marked for apoptosis are of lower reproductive potential than their unmarked counterparts. Defects in epididymal maturation lead to increased morphological abnormalities in the spermatozoa and poor sperm motility.,, In addition, immature spermatozoa exhibit metabolic alterations, presence of excess cytoplasm in the ejaculate, increased production of reactive oxygen species (ROS), lipid peroxidation, and DNA fragmentation.,
Sperm preparation methods such as density gradient separation are routinely used to obtain highly motile and morphologically normal sperm for assisted reproductive technology (ART). Spermatozoa separated on a three-layer density gradient (40%, 60%, and 80%) demonstrate cell-to-cell variation in both fertile and infertile men. Spermatozoa were separated into four fractions on the basis of their density and maturity. The lowest level of ROS production and DNA damage correlates with morphologically normal, motile spermatozoa obtained in the mature subset from fertile men compared with abnormal spermatozoa from infertile men.,, Different subsets of spermatozoa obtained from the ejaculate of fertile men after separation on a three-layer density gradient differ in their proteome profile. We have demonstrated an increasing trend in proteins involved in key biological processes during sperm maturation such as reproductive cellular process, gamete production, motility, oxidative phosphorylation, and energy metabolism. A decreasing trend was seen in the expression of proteins that were involved in key biological processes such as protein synthesis, protein transport, and response to oxidative stress.
Division of spermatozoa into phenotypes or subsets is of importance in the evaluation of its true fertility potential, particularly when using testicular spermatozoa for intracytoplasmic sperm injection (ICSI). A recent review by Esteves et al. documented that infertile couples may benefit from ICSI with testicular spermatozoa instead of ejaculated spermatozoa if the male partners exhibit high sperm DNA fragmentation (SDF) in the ejaculate. We have also demonstrated that proteins critical for sperm maturation, motility, and fertilization are involved in biological processes that are activated or suppressed in different subsets of spermatozoa from fertile men. However, the underlying pathways and the distribution of proteins in immature and mature sperm from infertile men have not been explored utilizing a proteomic approach. The present study is a continuation of our previous report on fertile donors deciphering the proteomic signatures in the spermatozoa of infertile patients to understand the underlying mechanism(s) of defective sperm maturation.
| Patients and Methods|| |
Following the approval of the study by the Institutional Review Board (IRB) of Cleveland Clinic (Cleveland, OH, USA), semen samples were collected from 11 infertile men. Men with leukocytospermia (Endtz test positive) and female factor infertility were excluded. All enrolled patients provided written consent to participate in the study.
Sample collection and semen analysis
Semen samples were examined according to 2010 World Health Organization (WHO) criteria. All specimens were collected by masturbation after sexual abstinence of 48–72 h and were allowed to liquefy for 20 min at 37°C before further processing. Following liquefaction, manual semen analysis was performed, including evaluation of presence of round cells, presence of white blood cells by the peroxidase test (Endtz test), viability, and morphology as previously described.
Separation of sperm phenotypes and proteomic analysis
For separating immature and mature spermatozoa, a three-layer density gradient was used as described earlier. It consisted of 2 ml of 40%, 60%, and 80% of the upper layer, intermediate layer, and lower layer, respectively. The three layers were reconstituted from the stock (100%) solution of the gradient with the SpermRinse medium (Vitrolife, San Diego, CA, USA). The stock gradient was an antibiotic-free bicarbonate and HEPES-buffered medium containing silane-coated, colloid silica particles. The SpermRinse medium was a bicarbonate and HEPES-buffered medium containing human serum albumin and gentamycin as an antibiotic (Vitrolife, San Diego, CA, USA). The three-layer gradient is a slight modification of the 2-layer density gradient method routinely used for preparing sperm for ART techniques, especially intrauterine insemination. Briefly, 1–2 ml of liquefied semen sample was carefully loaded on the 40% gradient and centrifuged at 300 g for 20 min. The resulting interfaces between the seminal plasma and 40% (fraction 1), 40% and 60% (fraction 2), 60% and 80% (fraction 3), and the 80% pellet (fraction 4, mature fraction) were carefully aspirated, resuspended in human tubal fluid media (HTF, Irvine Scientific, Santa Ana, CA, USA), and centrifuged at 300 g for 7 min. The pellets of each fraction were resuspended in 0.5–1 ml HTF, and the total sperm count, motility, and morphology were assessed again.
Sperm proteins from two individual samples and one pooled (from four individuals after normalization for spermatozoa number and protein content) sample were dissolved in RIPA lysis buffer (Sigma-Aldrich, St. Louis, MO, USA) supplemented with proteinase inhibitor cocktail (Roche, Indianapolis, IN, USA). Protein concentration was determined by the Bicinchoninic acid (BCA) kit (Thermo, Rockford, IL, USA). Proteins were extracted and separated by 1-dimensional gel electrophoresis. Bands were digested with trypsin and analyzed on a LTQ-Orbitrap Elite hybrid mass spectrometer (Thermo) system as described. Each sample was run in triplicate and the average was taken.
Database searching and protein identification
Tandem mass spectra were extracted by Proteome Discoverer (version 188.8.131.528, Thermo Fisher Scientific, San Jose, CA, USA). Twelve tandem mass spectrometry or MS/MS samples (3 runs per sample) were analyzed by using Mascot (version 2.3.02, Matrix Science, London, UK), Sequest (version 184.108.40.2068, Thermo Fisher Scientific), and X! Tandem (TheGPM, thegpm.org; version CYCLONE 2010.12.01.1). Mascot, Sequest, and X! Tandem were set up to search the human database (33 292 entries) assuming the digestion enzyme trypsin. To validate MS/MS-based peptide and protein identifications, Scaffold (version Scaffold 220.127.116.11, Proteome Software Inc., Portland, OR, USA) was used. Peptide identifications were accepted if they could be established at >95.0% probability by the PeptideProphet™ algorithm with Scaffold delta-mass correction.
Quantitation of the relative abundance of protein in spermatozoa
The relative abundance of sperm proteins was determined by comparing the number of spectra, termed spectral counts (SC), used to identify each protein. The numerical values used in the quantitation correspond to the normalized spectral abundance factor (NSAF, SC/[ΣSC] × protein length). NSAF approach was applied before relative protein quantification. Differentially expressed proteins (DEPs) were obtained by applying different constraints for significance tests or fold change cutoffs from the average SC of the protein from multiple runs.
The categorization of overall abundance along with the filtering criteria used for differential expression analysis is summarized below:
- Very low abundance: spectral count range 1.7–7; P ≤ 0.001; and NSAF ratio ≥2.5 for overexpressed and ≤0.4 for underexpressed proteins
- Low abundance: spectral count range 8–19; P ≤ 0.01; and NSAF ratio ≥2.5 for overexpressed and ≤0.4 for underexpressed proteins
- Medium abundance: spectral count range 20–79; P ≤ 0.05; and NSAF ratio ≥2.0 for overexpressed and ≤0.5 for underexpressed proteins
- High abundance: spectral counts >80; P ≤ 0.05; and NSAF ratio ≥1.5 for overexpressed and ≤0.67 for underexpressed proteins.
Functional bioinformatics analyses were done with publicly available tools such as Gene Ontology (GO) annotations from GO Term Finder (http://search.cpan.org/dist/GO-TermFinder/), GO Term Mapper (http://go.princeton.edu/cgi-bin/GOTermMapper), UNIPROT (The UniProt Consortium; http://www.uniprot.org/), Database for Annotation, Visualization and Integrated Discovery (DAVID; http://david.niaid.nih.gov), and proprietary software packages such as Ingenuity Pathway Analysis (IPA from Ingenuity® Systems; https://www.qiagenbioinformatics.com/products/ingenuity-pathway-analysis/) and MetaCore™ (GeneGo Inc., Encinitas, CA, USA) to identify the cellular distribution of proteins and differentially affected processes and pathways.
The results were expressed as mean ± standard deviation (s.d.). To compare the differences between different fractions of the ejaculate, we used Jonckheere–Terpstra test (or Jonckheere's trend test). It is similar to the Kruskal–Wallis test where the null hypothesis is that several independent samples are from the same population. However, there is no priori ordering of the populations from which the samples are drawn. When there is a priori ordering, the Jonckheere–Terpstra test has more power than the Kruskal–Wallis test. In this test, there is no issue of normality and does not require log transformation of the data. The statistical analysis was performed using the MedCalc (version 17.9.7, MedCalc Software bvba, Ostend, Belgium). For IPA and MetaCore™, the right-tailed Fisher's exact testwas used. Differences were considered statistically significant for P < 0.05.
| Results|| |
Total sperm count (TSC), presence of round cells, motility, and total motile sperm count (TMS) before and after density gradient centrifugation are shown in [Table 1] and [Figure 1]. The TSC (×106, mean ± s.d.) recovered in F3 and F4 was similar (32.51 ± 25.78 and 34.32 ± 32.53, respectively) and significantly different (P = 0.004665) from the TSC in the prewash sample (79.98 ± 84.58) [Figure 1]a. There was a significant decrease (P < 0.000001) in the number of round cells (median, 25th and 75th percentiles) from the prewash sample (1.90 [1.07, 3.75]) to that of F3 (0 [0, 0]) and F4 (0 [0, 0]) [Figure 1]b. Spermatozoa recovered from F4 displayed the highest average motility (63.0 (53.5, 83.1)) compared with F1, F2, and F3 followed by F3 (36.6 [32.3, 43.2]) compared with F1, F2, and F4 (P = 0.000001) [Figure 1]c. A higher recovery of TMS (×106 spermatozoa) was observed in F4 (23.46 ± 23.39) than F3 (12.68 ± 11.90) [Figure 1]d.
|Table 1: Semen parameters before and after separation on a 3-layer density gradient|
Click here to view
|Figure 1: Semen parameters in four fractions compared to PW sample. (a) Total sperm count (×106); (b) round cell count (×106); (c) motility (%); and (d) total motile count (×106) in F1–F4. PW: prewash; F1: least mature stage having the lowest density; F2, F3: intermediate stages; F4: includes the most dense and morphologically mature motile spermatozoa.|
Click here to view
Overall protein abundance
The relative abundance of the identified proteins in the sperm samples was quantified. For higher accuracy of proteomic results, the protein abundance should be similar between the analyzed samples. The range of total SC was 31188–51131 which showed that the abundance of proteins in these samples was similar.
Distribution of proteins in different fractions
All fractions, i.e. F1, F2, and F3 were compared with mature fraction F4. A total of 1585 proteins were identified in the four fractions together. Among those, 1202, 1140, 1129, and 890 proteins were detected in F1, F2, F3, and F4, respectively. By comparing the F1 and F4 proteomes, 136 proteins were overexpressed in F1, 177 were underexpressed in F1, 158 were unique to F1, and 51 were unique to F4. When F2 was compared with F4, 113 proteins were overexpressed in F2, 111 were underexpressed in F2, 114 were unique to F2, and 24 were unique to F4. Comparison of F3 versus F4 showed 89 proteins overexpressed in F3, 53 underexpressed in F3, 38 unique to F3, and 8 unique to F4. All the four fractions were compared to identify the differentially expressed proteins (DEPs). Of the 656 DEPs, 75 proteins showed an increasing trend from fraction F1 to F4 [Supplementary Table 1 [Additional file 1]], while 279 showed a decreasing trend [Supplementary Table 2 [Additional file 2]].
Identification of functional proteins
The Reactome pathway analysis showed that the majority of DEPs were involved in metabolism, particularly protein metabolism, disease processes, gene expression, and signal transduction [Figure 2]a. The percentage of proteins involved in various pathways that were either underexpressed or showed a decreasing trend or were overexpressed and showed an increasing trend from F1 to F4 is shown in [Figure 2]b.
|Figure 2: DEPs in fraction 1 (F1) through fraction 4 (F4). (a) Reactome pathway showing the number of DEPs involved in various pathways in the four fractions. F1 being the most immature and F4 the most mature. (b) Percentage of proteins in the data set that were increasing or decreasing for the different pathways. (c) IPA network showing molecular chaperons that were involved in developmental disorder, posttranslational modifications, and protein folding. (d) IPA network showing proteins involved in molecular transport, protein trafficking, and cell cycle and (e) IPA network showing molecular chaperons involved in cell death and survival, posttranslational modifications, and protein folding. IPA: ingenuity pathway analysis; DEPs: differentially expressed proteins; F1: least mature stage having the lowest density; F2, F3: intermediate stages; F4: includes the most dense and morphologically mature motile spermatozoa. Full names of abbreviated proteins are presented in Supplementary Information.|
Click here to view
The functional annotations from use of the DAVID tool and their enrichment analysis with the number of key proteins involved are presented in [Supplementary Table 3]. F1 was characterized by proteins involved in translation, elongation, protein transport, oxidoreductase activity, reproductive processes, spermatid development/differentiation, and regulation of apoptotic pathways. The proteins identified in F2 are involved in cell differentiation, protein–protein interactions, protein transport and localization, oxidoreductase activity, spermatogenesis, gonad development, and proteolytic pathways. The key processes of F3 proteins include generation of precursor metabolites and integration of energy metabolism, oxidative phosphorylation, protein catabolic process, and protein ubiquitination. Finally, the proteins identified in F4 were largely involved in reproductive functions [Table 2]. The proteins involved in intracellular transport, oxidation–reduction, cellular amino acid catabolic process, and alternative splicing showed an increasing trend from F1 to F4 [Supplementary Table 2]. On the other hand, the proteins involved in spermatogenesis, protein metabolism, cell cycle, integration of energy metabolism, regulation of apoptosis, cell redox homeostasis, translational elongation, and response to protein folding showed a decreasing trend from F1 to F4 [Supplementary Table 3 [Additional file 3]].
|Table 2: Key proteins showing an increasing trend of expression from fraction 1 (F1: most immature) through faction 4 (F4: most mature) involved in structural assembly of spermatozoa, spermatogenesis, reproduction, sperm motility, energy metabolism, and oxidation–reduction processes|
Click here to view
DEPs involved in various networks
The DEPs were further subjected to network analysis by using IPA. A total of 161 pathways were linked to 279 proteins showing a decreasing trend from F1 to F4. Among these proteins, 58 were involved in cell death and survival, posttranslational modification, and protein folding [Figure 2]c. In addition, there were 58 proteins, mostly molecular chaperones, which were involved in developmental disorder, posttranslational modification, and protein folding [Figure 2]d. A key network was identified with 58 focal proteins involved in molecular transport, protein trafficking, and cell cycle [Figure 2]e, in which 18 proteins were involved in cell signaling, cancer, and cellular development processes. The proteins showing an increasing trend from immature to mature fraction (F1 to F4) were linked to 66 networks. The top networks were lipid metabolism, small molecule biochemistry, and molecular transport that included 25 proteins. Of these, 22 proteins participated in cell-to-cell signaling and interaction, cellular function and maintenance, and inflammatory response and 12 proteins were involved in cell death and survival, cellular compromise, and cancer [Table 2].
Proteins involved in top pathways and upstream regulators
On the basis of the common genes, 15 top pathways showed an overlap with each other and their translated protein product was shown to be gradually underexpressed from immature to mature fractions, while 8 connecting signaling pathways were predicted [Figure 3]. The proteins identified in our data set were analyzed using MetaCore™ to predict the upstream regulator(s) and the pathways involved [Figure 4]a,[Figure 4]b,[Figure 4]c. The androgen receptor was predicted as one of the transcription regulators and the pathways affected were chaperone-mediated stress response, proteosomal pathway, and sperm function [Figure 4]d.
|Figure 3: Fifteen overlapping pathways and eight connecting signaling pathways were predicted based on the identified common genes. Their translated protein products were underexpressed starting from F1 to F4. F1: least mature stage having the lowest density; F2, F3: intermediate stages; F4: includes the most dense and morphologically mature motile spermatozoa. Full names of abbreviated proteins are presented in Supplementary Information.|
Click here to view
|Figure 4: The biochemical process networks of the differentially expressed proteins generated by MetaCore™ software using a transcription regulation algorithm. The biochemical process regulation networks were ranked by a value and interpreted in terms of GO. Gene network illustrates proteins and interactions of the differentially expressed proteins. MetaCore™ pathway analysis predicting androgen receptor as the major regulator and pathways involved were; (a) chaperone-mediated stress response; (b) sperm function; (c) proteosomal pathway; and (d) predicted impaired androgen receptor signaling in spermatozoa of infertile patients. GO: gene ontology; F1: least mature stage having the lowest density; F2, F3: intermediate stages; F4: includes the most dense and morphologically mature motile spermatozoa. Full names of abbreviated proteins are presented in Supplementary Information.|
Click here to view
| Discussion|| |
We have previously reported the proteome profile changes in these four subsets of immature and mature spermatozoa from fertile men separated by three layer density gradient. Our results demonstrated that proteins critical for sperm maturation are altered in stage specific manner in fertile men. Of these, 4 proteins, namely, heat shock protein HSP701A, clusterin, tektin 2, and tektin 3 were validated by western blot to corroborate the proteomic findings. In this study, the goal was to compare the same four subsets of immature and mature spermatozoa in infertile men to understand why spermatozoa obtained in the mature fraction (F4), which separate the good quality spermatozoa in terms of motility, TMC, and even TSC, are still dysfunctional in achieving pregnancy. Therefore, in this study, a high-throughput shotgun proteomic approach followed by pathway analysis was used to understand the impaired molecular mechanism involved in sperm dysfunction.
Despite the absence of effective transcription and translation, spermatozoa undergo important functional transformation during epididymal transit. These modifications are essential to produce functionally efficient spermatozoa and depend on loss, modification, or remodeling of existing sperm proteins in response to the signals conveyed by the male reproductive tract. These signals are carefully regulated by an array of gene products in a stage-specific manner in which molecular chaperones play a key role., Molecular chaperones are structurally diverse proteins that are expressed nearly in all cells in response to cellular stress. In the present study, we observed that two of the top three networks identified by IPA were associated with underexpression of molecular chaperones. These networks were related to posttranslational modification, protein folding, and molecular transport. Furthermore, these networks were involved in protein trafficking related to cellular processes such as cell cycle, cell survival, cell death, and developmental disorder.
Proposed pathway 1
In the first network, eight proteins of the chaperonin-containing T-complex/TCP1-ring complex (CCT2, CCT3, CCT4, CCT5, CCT6A, CCT7, CCT8, and TCP1), four HSP 70 family chaperones (HSPA2, HSPA1L, HSP90AA1, and HSP8), and three members of the co-chaperonin DNAJ (DNAJA1, DNAJB1, and DNAJC3) were identified as key proteins [Figure 2]c. Vesicular trafficking is conducted by gene products of Bardet–Biedl syndrome (BBS) and McKusick–Kaufman syndrome (MKKS)., Disruptions of BBS or MKKS genes result in male infertility owing to the failure of flagellum formation in spermatozoa.,,, The chaperonin-containing T-complex/TCP1-ring complex as part of the BBS/CCT complex may play a role in the assembly of BBSome, a complex involved in ciliogenesis regulating transport vesicles to the cilia. These proteins also have a role in the folding of actin and tubulin in vitro, and TUBB4B protein was identified in this network. Our liquid chromatography (LC)-MS/MS data also showed a decrease in tubulin (TUBB4B) as well as in both dynein heavy and light chains (DYNC1H1 and DYNC1LI1), further corroborating the concept of flagellar disassembly. Owing to the defects in flagellar disassembly, we did not find any change in Tektin (structural components of outer doublet microtubules) proteins as reported by us in spermatozoa of fertile men.
The ribosomal subunits that include two translation elongation factors (EEF1E1 and EEF1D) and two tRNA ligases (DARS and VARS) identified in the network were also underexpressed. These are responsible for linking aspartate and valine to their cognate tRNAs, thus inhibiting the incorporation of these amino acids into the protein. In fact, Talevi et al. have reported that treatment of spermatozoa with aspartate in vitro enhances sperm motility in oligozoospermic samples along with inhibition of DNA damage and membrane lipid peroxidation.
The HSP family proteins identified in the network were inducible by stress (e.g. HSPA1), constitutively expressed, or both (e.g. HSPH1, HSPA8, and HSP90AA1). Expression of some HSPs is developmentally regulated or restricted to specific cells. HSPA2 is expressed on the sperm surface, and it is essential for sperm membrane remodeling during sperm–oocyte interaction. It may be used as a biochemical marker of human sperm function and male fertility. Levels of human HSPA2 expression have been correlated with sperm maturity and male fertility., Therefore, a decline in HSPA2 expression may be responsible for improper maturation, which was not so in our previous study in fertile men.
We found HSPA1A to be downregulated from F1 to F4 in this study in infertile men as well as in our previous study in fertile donors. HSPA1A plays a pivotal role in the protein quality control system, ensuring the correct folding of proteins, re-folding of misfolded proteins, and controlling the targeting of proteins for subsequent degradation. It maintains protein homeostasis during cellular stress through two opposing mechanisms: protein refolding and degradation.
The other HSPs and DNAJ family proteins are reported to affect spermatogenesis or sperm function. Proteins responsible for mitochondria-mediated cell death (FAM162) and the negative regulator of mitochondrial respiration (TRAP1) were also downregulated [Figure 2]c. Therefore, this network indicates that during the process of sperm production and maturation in infertile men, there is an impairment of the stress response protein that is demonstrated by the underexpression of chaperones in the network. This may lead to defects in the formation of the flagella and abortive apoptosis. Therefore, though the spermatozoa obtained in the most mature fraction (F4) appear normal, they are incapable of proper functioning after ejaculation.
The above findings were corroborated by the second predicted network [Figure 2]d showing deregulated posttranslational modification and protein folding, leading to developmental disorder.
Proposed pathway 2
The most important proteins were thioredoxins and peroxiredoxins which are reported to be abundantly expressed in developing spermatocytes. Thioredoxin domain-containing protein 2 (TXNDC2) is a spermatid-specific thioredoxin-1 (also known as Sptrx-1) which is found exclusively in the tail of elongating spermatids and spermatozoa. TXNDC2 transiently associates with the longitudinal columns of the fibrous sheath during sperm tail assembly but does not remain as a permanent component of the fibrous sheath in the mature sperm cell. Although the expression of this protein showed a gradual decline from F1 to F4, still its presence in the most mature fraction implies defects in sperm flagella, which corroborates the findings in pathway 1 and our previous study regarding the expression of tektins. Similarly, peroxiredoxin 4 (PRDX4) [Figure 2]d is present as a membrane-bound form only in testes, and the unprocessed form may be involved in acrosome formation during spermiogenesis. Peroxiredoxins are involved in hydrogen peroxide-mediated signaling in spermatozoa and prevent oxidative stress during capacitation. PRDX4 knockout mice showed elevated spermatogenic cell death via oxidative stress. We observed a decrease in the expression of voltage-dependent mitochondrial outer membrane ion channels (VDAC1, VDAC2, and VDAC3) along with mitochondrial ATP synthase complex subunits (ATP5I, ATP5F1, and ATP5J). These findings suggest mitochondrial dysfunction, which may lead to oxidative stress and cell death as mentioned above. In addition, underexpression of isocitrate dehydrogenase (IHD1), lipoprotein lipase (LPL), and apolipoprotein A1 (APOA1) may be responsible for impaired lipid metabolism. This may result in altered lipid profile of the sperm membrane that in turn may affect fertilization. Furthermore, underexpression of dynactin (DCTN2), contractin (actin-related protein 1 homolog A [ACTR1A]), and F-actin capping protein (capping protein Z line alpha [CAPZA]) observed in our study may lead to defects in microtubule packaging and centrosome formation. Since the spermatozoon donates its centriole during fertilization and is essential for cleavage, these spermatozoa from infertile patients may not be able to initiate cleavage in embryos after fertilization. One of the major groups from the network that we identified were the endoplasmic reticulum-specific protein suppressor enhancer lin-12 1 like (SEL1L) and endoplasmic reticulum protein 29 (ERP29) along with the Ras-related proteins RBA6A and RAB7A. The presence of these proteins suggests impaired membrane traffic from the Golgi apparatus toward the endoplasmic reticulum. In addition, proteins involved in posttranslational modification and folding of proteins such as chaperones (HSP90B1, HSPA5, calnexin [CANX], calreticulin [CALR]), protein disulfide isomerase (prolyl 4-hydroxylase, beta polypeptide [P4HB], and protein disulfide isomerase-associated 3 [PDIA3]), peptidyl prolyl-cis-trans isomerase (FK506-binding protein [FKBP]), and serine-threonine phosphatase (peptidylpropyl isomerase B [PPIB]) were also identified in this network. Thus, an endoplasmic reticulum stress may be responsible for improper protein formation leading to developmental disorders in the proteins obtained in F4 fraction in infertile men.
Proposed pathway 3
The third network demonstrated a deregulated molecular transport, protein trafficking, and cell cycle. In the present study, we identified nuclear transport receptor IOP5 or importin-5 and alpha karyopherins such as KPNA2, KPNA3, KPNA4, and KNPNB. These karyopherins are expressed in specific stages during spermatogenesis in the spermatogonium (KNNA1, KPNA2, and KPNA3), spermatocyte (KPNA2, KPNA3, and KPNA4), round spermatid (KPNA3 and KPNA4), and elongating spermatids (KPNA2 and IPO5). These importins were all underexpressed in ejaculated sperm in our study. This finding suggests that together with RNA GTPase they may be responsible for impaired import of testis specific H1 histone (HIST1H1T). HIST1H1T is required for less compaction of nuclear DNA during meiosis and is subsequently replaced by protamines. Thus, these importins may play a role in nuclear abnormalities that may eventually render the spermatozoa incapable of fertilization or successful completion of postfertilization steps.
The other two important proteins involved in exosome-mediated cargo delivery were ANXA1 and ANXA2. These were underexpressed in this network. It is hypothesized that sperm maturation in the epididymis and regulation of sperm function are mediated by regulatory RNAs that are delivered by exosomes., In fact, Yang et al. reported the presence of annexins in the seminal exosomes. Thus, an improper delivery of regulatory factors to the spermatozoa via exosomes might result in sperm dysfunction.
Proteasomes present in mammalian sperm play a pivotal role during fertilization. The enzymatic activity of the proteasome is modulated by protein kinase A and is involved in the progesterone-induced acrosome reaction. Furthermore, after capacitation, the acrosomal proteasomes facilitate the degradation of zona pellucida glycoproteins leading up to fertilization. In this network, we identified eight subunits of 26S proteasome (PSMC3, PSMC4, PSMC6, PSMD5, PSMD7, PSMD8, PSMD11, and PSMD13). One of the nodal proteins, matrix-metalloproteinase (MME) or Neprilysin, was also identified in the network. Pinto et al. reported that tachykinins present in human spermatozoa participate in the regulation of sperm motility, and their activity is regulated by Neprilysin. Another nodal protein involved in this network was 14-3-3 protein theta (YWHAQ) adapter protein which has been implicated in the regulation of a large spectrum of both general and specialized signaling pathways, thus implying a signaling failure in F4 fraction in our study. We identified 15 overlapping networks on the basis of mutual genes that were affected during the process of maturation in these infertile patients. These networks may be responsible for poor sperm quality.
MetaCore™ pathway analysis
While IPA reports protein–protein interactions, MetaCore™ is an integrated software suite for functional analysis of experimental data of human protein–protein interactions, protein–DNA interactions, transcription factors, and signaling and metabolic pathways including disease and toxicity. The top pathway affected during maturation was identified to be mediated through TXN, PRDX1, signaling adaptor molecules (YWHAQ and YWHAB), and the molecular chaperons (HSP90AA1, HSP90B1, HSC70, HSPA5, DNAJA1, DNAJ2, DNAJB1, DNAJB11, DNAJC3, and STIP1). Taken together, all the networks and pathways described above suggest that the spermatozoa at different stages of maturation suffer from endoplasmic reticulum stress (ER stress) due to accumulation of unfolded or misfolded proteins as a result of underexpressed chaperone and proteasome activities, which will lead to induction of oxidative stress and abortive apoptosis. This is corroborated by the fact that the Clusterin is not differentially expressed in the spermatozoa of infertile men like that reported for fertile men. Clusterin isoforms are responsible for differential regulation of apoptosis where the nuclear form promotes the process and the mitochondrial form opposes the process (Uniprot). In fertile men, its proper expression prevented apoptosis of mature spermatozoa where the mitochondria are intact and eliminated the spermatozoa with nuclear defects. Therefore, it is suggested that the process was impaired in infertile men leading to abortive apoptosis. Taken together, the cells that may appear morphologically mature may carry structural anomalies in the cytoskeleton, microtubules, and flagellum, thereby affecting the function of the spermatozoa. Androgen receptor is a transcription regulator. The presence of the AR in human spermatozoa has been demonstrated by Western blot and immunofluorescence assay. AR is also localized in the head region. In addition to stimulating cell growth, androgens or the AR plays an important role in apoptosis involving both intrinsic and extrinsic pathways. Short exposure of ejaculated spermatozoa to androgens produces an increase in AR phosphorylation, especially on the 110 kDa band which is the less expressed isoform in sperm cells.
Taken together, all the networks and pathways described in the present study show that the spermatozoa at different stages of maturation suffer from ER stress due to accumulation of unfolded or misfolded proteins as a result of underexpressed chaperone and proteasome activities. This may lead to induction of oxidative stress and abortive apoptosis. Thus, although the spermatozoa in F4 fraction may appear morphologically mature and exhibit good motility and morphology following separation on density gradient as is the case during intrauterine insemination or in vitro fertilization, these sperm can potentially carry structural anomalies in the cytoskeleton, microtubule, and flagellum, thereby affecting the functional capabilities of the spermatozoa. A limitation of our study was that we could not perform the western blot validation of the MS-spectrometry data as the quantity of protein from the fractions, particularly F3 and F4, was insufficient.
In conclusion, the results of the present study show that a defective signaling cascade is responsible for the defective sperm function observed in infertile patients, particularly a decline in mitochondrial function and oxidative phosphorylation in the most mature fraction (F4), implying a state of energy deprivation. Furthermore, a dysregulated protein turnover and protein folding may lead to accumulation of defective proteins or proteins that otherwise would have been eliminated during the process of maturation and thus result in the impairment of sperm function. The present study provides mounting evidence that aberrant chaperone expression may be a major contributing factor to the defective sperm function seen in many cases of male infertility. Our study also identified the involvement of AR as the core player in the downregulation of the signaling leading to production of defective spermatozoa incapable of fertilization. Future studies conducted in larger study population are necessary to examine the distribution of proteins in infertile men with different clinical diagnosis such as varicocele and validate the important proteins to further help in understanding the underlying pathology of male infertility.
| Author Contributions|| |
AA and RS planned the experiments. ZC selected the samples according to the clinical histories. ZC and LS analyzed the collected data. LS was involved in the manuscript preparation. RS, AA, and LS critically reviewed the manuscript. All authors read and approved the final manuscript.
| Competing Interests|| |
All authors declared no competing interests.
| Acknowledgments|| |
The authors are grateful to the Andrology Center technologists for scheduling the study subjects; Cleveland Clinic Arts Department for creating the four figures; Belinda Willard, Director, Proteomic Core Lab, Lerner Research Institute for providing assistance with proteomic analysis and Banu Gopalan for the bioinformatics analysis of the results. The Orbitrap Elite mass spectrometer used in this study was purchased with funds from an NIH shared instrument grant 1S10RR031537-01 to Belinda Willard. Financial support was provided by the American Center for Reproductive Medicine, Cleveland Clinic. Dr. Zhihong Cui's visit was supported by a fellowship from the Chinese Government and the Institute of Toxicology, College of Preventive Medicine, the Third Military Medical University, Chongqing, China.
[Supplementary information [Additional file 4]] is linked to the online version of the paper on the Asian Journal of Andrology website.
| References|| |
Jones RE, Kristen L. The male reproductive system. In: Jones R, Kristen L, editors. Human Reproductive Biology. Cambridge: Academic Press, Elsevier; 2013.
Sharma R, Agarwal A. Spermatogenesis: An overview. In: Armand Zini A, editor. Sperm chromatin: biological and clinical applications in male infertility and assisted reproduction. New York: Springer; 2011. p
Kishigami S, Wakayama S, Nguyen VT, Wakayama T. Similar time restriction for intracytoplasmic sperm injection and round spermatid injection into activated oocytes for efficient offspring production. Biol Reprod
2004; 70: 1863–9.
Chemes HE. Phenotypes of sperm pathology: genetic and acquired forms in infertile men. J Androl
2000; 21: 799–808.
Gil-Guzman E, Ollero M, Lopez MC, Sharma RK, Alvarez JG, et al
. Differential production of reactive oxygen species by subsets of human spermatozoa at different stages of maturation. Hum Reprod
2001; 16: 1922–30.
Gaucher J, Reynoird N, Montellier E, Boussouar F, Rousseaux S, et al
. From meiosis to postmeiotic events: the secrets of histone disappearance. FEBS J
2010; 277: 599–604.
Guo W, Qu F, Xia L, Guo Q, Ying X, et al
. Identification and characterization of ERp29 in rat spermatozoa during epididymal transit. Reproduction
2007; 133: 575–84.
Govin J, Gaucher J, Ferro M, Debernardi A, Garin J, et al
. Proteomic strategy for the identification of critical actors in reorganization of the post-meiotic male genome. Mol Hum Reprod
2012; 18: 1–13.
Cui Z, Sharma R, Agarwal A. Proteomic analysis of mature and immature ejaculated spermatozoa from fertile men. Asian J Androl
2016; 18: 735–46.
Sakkas D, Mariethoz E, St John JC. Abnormal sperm parameters in humans are indicative of an abortive apoptotic mechanism linked to the Fas-mediated pathway. Exp Cell Res
1999; 251: 350–5.
Sakkas D, Mariethoz E, Manicardi G, Bizzaro D, Bianchi PG, et al
. Origin of DNA damage in ejaculated human spermatozoa. Rev Reprod
1999; 4: 31–7.
Dun MD, Aitken RJ, Nixon B. The role of molecular chaperones in spermatogenesis and the post-testicular maturation of mammalian spermatozoa. Hum Reprod Update
2012; 18: 420–35.
Lin C, Tholen E, Jennen D, Ponsuksili S, Schellander K, et al
. Evidence for effects of testis and epididymis expressed genes on sperm quality and boar fertility traits. Reprod Domest Anim
2006; 41: 538–43.
Ollero M, Gil-Guzman E, Lopez MC, Sharma RK, Agarwal A, et al
. Characterization of subsets of human spermatozoa at different stages of maturation: implications in the diagnosis and treatment of male infertility. Hum Reprod
2001; 16: 1912–21.
Brahem S, Mehdi M, Elghezal H, Saad A. Semen processing by density gradient centrifugation is useful in selecting sperm with higher double-strand DNA integrity. Andrologia
2011; 43: 196–202.
Castillo J, Simon L, de Mateo S, Lewis S, Oliva R. Protamine/DNA ratios and DNA damage in native and density gradient centrifuged sperm from infertile patients. J Androl
2011; 32: 324–32.
Esteves SC, Roque M, Bradley CK, Garrido N. Reproductive outcomes of testicular versus ejaculated sperm for intracytoplasmic sperm injection among men with high levels of DNA fragmentation in semen: systematic review and meta-analysis. Fertil Steril
2017; 108: 456–67.e1.
World Health Organization. WHO laboratory manual for the examination and processing of human semen. 5th
ed. Geneva: WHO Press; 2010.
Keller A, Nesvizhskii AI, Kolker E, Aebersold R. Empirical statistical model to estimate the accuracy of peptide identifications made by MS/MS and database search. Anal Chem
2002; 74: 5383–92.
Zybailov B, Coleman MK, Florens L, Washburn MP. Correlation of relative abundance ratios derived from peptide ion chromatograms and spectrum counting for quantitative proteomic analysis using stable isotope labeling. Anal Chem
2005; 77: 6218–24.
Boyle EI, Weng S, Gollub J, Jin H, Botstein D, et al
. GO: termfinder – open source software for accessing gene ontology information and finding significantly enriched gene ontology terms associated with a list of genes. Bioinformatics
2004; 20: 3710–5.
Hermo L, Pelletier RM, Cyr DG, Smith CE. Surfing the wave, cycle, life history, and genes/proteins expressed by testicular germ cells. Part 1: background to spermatogenesis, spermatogonia, and spermatocytes. Microsc Res Tech
2010; 73: 241–78.
Hermo L, Pelletier RM, Cyr DG, Smith CE. Surfing the wave, cycle, life history, and genes/proteins expressed by testicular germ cells. Part 5: intercellular junctions and contacts between germs cells and sertoli cells and their regulatory interactions, testicular cholesterol, and genes/proteins associated with more than one germ cell generation. Microsc Res Tech
2010; 73: 409–94.
Ellis RJ. Discovery of molecular chaperones. Cell Stress Chaperones
1996; 1: 155–60.
Nachury MV, Loktev AV, Zhang Q, Westlake CJ, Peranen J, et al
. A core complex of BBS proteins cooperates with the GTPase Rab8 to promote ciliary membrane biogenesis. Cell
2007; 129: 1201–13.
Lechtreck KF, Johnson EC, Sakai T, Cochran D, Ballif BA, et al
. The chlamydomonas reinhardtii BBSome is an IFT cargo required for export of specific signaling proteins from flagella. J Cell Biol
2009; 187: 1117–32.
Mykytyn K, Mullins RF, Andrews M, Chiang AP, Swiderski RE, et al
. Bardet-Biedl syndrome type 4 (BBS4)-null mice implicate Bbs4 in flagella formation but not global cilia assembly. Proc Natl Acad Sci U S A
2004; 101: 8664–9.
Nishimura DY, Swiderski RE, Searby CC, Berg EM, Ferguson AL, et al
. Comparative genomics and gene expression analysis identifies BBS9, a new Bardet-Biedl syndrome gene. Am J Hum Genet
2005; 77: 1021–33.
Davis RE, Swiderski RE, Rahmouni K, Nishimura DY, Mullins RF, et al
. A knockin mouse model of the Bardet-Biedl syndrome 1 M390R mutation has cilia defects, ventriculomegaly, retinopathy, and obesity. Proc Natl Acad Sci U S A
2007; 104: 19422–7.
Talevi R, Barbato V, Fiorentino I, Braun S, Longobardi S, et al
. Protective effects of in vitro
treatment with zinc, d-aspartate and coenzyme q10 on human sperm motility, lipid peroxidation and DNA fragmentation. Rep Biol Endocrinol
2013; 11: 81.
Rupik W, Jasik K, Bembenek J, Widlak W. The expression patterns of heat shock genes and proteins and their role during vertebrate's development. Comp Biochem Physiol A Mol Integr Physiol
2011; 159: 349–66.
Naaby-Hansen S, Diekman A, Shetty J, Flickinger CJ, Westbrook A, et al
. Identification of calcium-binding proteins associated with the human sperm plasma membrane. Reprod Biol Endocrinol
2010; 8: 6.
Mamelak D, Lingwood C. Expression and sulfogalactolipid binding specificity of the recombinant testis-specific cognate heat shock protein 70. Glycoconj J
1997; 14: 715–22.
Seo JH, Park JH, Lee EJ, Vo TT, Choi H, et al
. ARD1-mediated Hsp70 acetylation balances stress-induced protein refolding and degradation. Nat Commun
2016; 7: 12882.
Hanschmann EM, Godoy JR, Berndt C, Hudemann C, Lillig CH. Thioredoxins, glutaredoxins, and peroxiredoxins – molecular mechanisms and health significance: from cofactors to antioxidants to redox signaling. Antioxid Redox Signal
2013; 19: 1539–605.
Yu Y, Oko R, Miranda-Vizuete A. Developmental expression of spermatid-specific thioredoxin-1 protein: transient association to the longitudinal columns of the fibrous sheath during sperm tail formation. Biol Reprod
2002; 67: 1546–54.
Ranney MK, Ahmed IS, Potts KR, Craven RJ. Multiple pathways regulating the anti-apoptotic protein clusterin in breast cancer. Biochim Biophys Acta
2007; 1772: 1103–11.
Matsuki S, Sasagawa I, Iuchi Y, Fujii J. Impaired expression of peroxiredoxin 4 in damaged testes by artificial cryptorchidism. Redox Report
2002; 7: 276–8.
O'Flaherty C, de Souza AR. Hydrogen peroxide modifies human sperm peroxiredoxins in a dose-dependent manner. Biol Reprod
2011; 84: 238–47.
Lee D, Moawad AR, Morielli T, Fernandez MC, O'Flaherty C. Peroxiredoxins prevent oxidative stress during human sperm capacitation. Mol Hum Reprod
2017; 23: 106–15.
Iuchi Y, Okada F, Tsunoda S, Kibe N, Shirasawa N, et al
. Peroxiredoxin 4 knockout results in elevated spermatogenic cell death via oxidative stress. Biochem J
2009; 419: 149–58.
Koh HJ, Lee SM, Son BG, Lee SH, Ryoo ZY, et al
. Cytosolic NADP+
-dependent isocitrate dehydrogenase plays a key role in lipid metabolism. J Biol Chem
2004; 279: 39968–74.
Major AT, Whiley PA, Loveland KL. Expression of nucleocytoplasmic transport machinery: clues to regulation of spermatogenic development. Biochim Biophys Acta
2011; 1813: 1668–88.
Leoni G, Neumann PA, Kamaly N, Quiros M, Nishio H, et al
. Annexin A1-containing extracellular vesicles and polymeric nanoparticles promote epithelial wound repair. J Clin Invest
2015; 125: 1215–27.
Vojtech L, Woo S, Hughes S, Levy C, Ballweber L, et al
. Exosomes in human semen carry a distinctive repertoire of small non-coding RNAs with potential regulatory functions. Nucleic Acids Res
2014; 42: 7290–304.
Yang C, Guo WB, Zhang WS, Bian J, Yang JK, et al
. Comprehensive proteomics analysis of exosomes derived from human seminal plasma. Andrology
2017; 5: 1007–15.
Kerns K, Morales P, Sutovsky P. Regulation of sperm capacitation by the 26S proteasome: an emerging new paradigm in spermatology. Biol Reprod
2016; 94: 117.
Pinto FM, Ravina CG, Subiran N, Cejudo-Roman A, Fernandez-Sanchez M, et al
. Autocrine regulation of human sperm motility by tachykinins. Reprod Biol Endocrinol
2010; 8: 104.
Olsen JV, Blagoev B, Gnad F, Macek B, Kumar C, et al
. Global, in vivo
, and site-specific phosphorylation dynamics in signaling networks. Cell
2006; 127: 635–48.
Solakidi S, Psarra AM, Nikolaropoulos S, Sekeris CE. Estrogen receptors alpha and beta (ERalpha and ERbeta) and androgen receptor (AR) in human sperm: localization of ERbeta and AR in mitochondria of the midpiece. Hum Reprod
2005; 20: 3481–7.
Aquila S, Middea E, Catalano S, Marsico S, Lanzino M, et al
. Human sperm express a functional androgen receptor: effects on PI3K/AKT pathway. Hum Reprod
2007; 10: 2594–605.
Nakagawa T, Zhu H, Morishima N, Li E, Xu J, et al
. Caspase-12 mediates endoplasmic-reticulum-specific apoptosis and cytotoxicity by amyloid-beta. Nature
2000; 403: 98–103.
[Figure 1], [Figure 2], [Figure 3], [Figure 4]
[Table 1], [Table 2]