Integrative bioinformatics approaches for identifying potential biomarkers and pathways involved in non-obstructive azoospermia

Cancer risk and male Infertility: Unravelling predictive biomarkers and prognostic indicators

In recent years, there has been a global increase in cases of male infertility. There are about 30 million cases of male infertility worldwide and male reproductive health is showing rapid decline in last few decades. It is now recognized as a potential risk factor for developing certain types of cancer, particularly genitourinary malignancies like testicular and prostate cancer. Male infertility is considered a potential indicator of overall health and an early biomarker for cancer. Cases of unexplained male factor infertility have high levels of oxidative stress and oxidative DNA damage and this induces both denovo germ line mutations and epimutations due to build up of 8-hydroxy 2 deoxygunaosine abase which is highly mutagenic and also induces hypomethylation and genomic instability. Consequently, there is growing evidence to explore the various factors contributing to an increased cancer risk. Currently, the available prognostic and predictive biomarkers associated with semen characteristics and cancer risk are limited but gaining significant attention in clinical research for the diagnosis and treatment of elevated cancer risk in the individual and in offspring. The male germ cell being transcriptionally and translationally inert has a highly truncated repair mechanism and has minimal antioxidants and thus most vulnerable to oxidative injury due to environmental factors and unhealthy lifestyle and social habits. Therefore, advancing our understanding requires a thorough evaluation of the pathophysiologic mechanisms at the DNA, RNA, protein, and metabolite levels to identify key biomarkers that may underlie the pathogenesis of male infertility and associated cancer. Advanced methodologies such as genomics, epigenetics, proteomics, transcriptomics, and metabolomics stand at the forefront of cutting-edge approaches for discovering novel biomarkers, spanning from infertility to associated cancer types. Henceforth, in this review, we aim to assess the role and potential of recently identified predictive and prognostic biomarkers, offering insights into the success of assisted reproductive technologies, causes of azoospermia and idiopathic infertility, the impact of integrated holistic approach and lifestyle modifications, and the monitoring of cancer susceptibility, initiation and progression. Comprehending these biomarkers is crucial for providing comprehensive counselling to infertile men and cancer patients, along with their families.

Common as well as unique methylation-sensitive DNA regulatory elements in three mammalian SLC9C1 genes

The SLC9C1 gene (which encodes the NHE10 protein) is essential for male fertility in both mice and humans, however the epigenetic mechanisms regulating its testis/sperm-specific gene expression have yet to be studied. Here we identify and characterize DNA regulatory elements of the SLC9C1 gene across three mammalian species: mouse, rat, and human. First, in silico analysis of these mammalian SLC9C1 genes identified a CpG island located upstream of the transcription start site in the same relative position in all three genes. Further analysis reveals that this CpG island behaves differently, with respect to gene regulatory activity, in the mouse SLC9C1 gene than it does in the rat and human SLC9C1 gene. The mouse SLC9C1 CpG island displays strong promoter activity by itself and seems to have a stronger gene regulatory effect than either the rat or human SLC9C1 CpG islands. While the function of the upstream SLC9C1 CpG island may be divergent across the three studied species, it appears that the promoters of these three mammalian SLC9C1 genes share similar DNA methylation-sensitive regulatory mechanisms. All three SLC9C1 promoter regions are differentially methylated in lung and testis, being more hypermethylated in lung relative to the testis, and DNA sequence alignments provide strong evidence of primary sequence conservation. Luciferase assays reveal that in vitro methylation of constructs containing different elements of the three SLC9C1 genes largely exhibit methylation-sensitive promoter activity (reduced promoter activity when methylated) in both HEK 293 and GC-1spg cells. In total, our data suggest that the DNA methylation-sensitive elements of the mouse, rat, and human SLC9C1 promoters are largely conserved, while the upstream SLC9C1 CpG island common to all three species seems to perform a different function in mouse than it does in rat and human. This work provides evidence that while homologous genes can all be regulated by DNA methylation-dependent epigenetic mechanisms, the location of the specific cis-regulatory elements responsible for this regulation can differ across species.

The human sperm proteome-Toward a panel for male fertility testing

BACKGROUND

Although male factor accounts for 40%-50% of unintended childlessness, we are far from fully understanding the detailed causes. Usually, affected men cannot even be provided with a molecular diagnosis.

OBJECTIVES

We aimed at a higher resolution of the human sperm proteome for better understanding of the molecular causes of male infertility. We were particularly interested in why reduced sperm count decreases fertility despite many normal-looking spermatozoa and which proteins might be involved.

MATERIAL AND METHODS

Applying mass spectrometry analysis, we qualitatively and quantitatively examined the proteomic profiles of spermatozoa from 76 men differing in fertility. Infertile men had abnormal semen parameters and were involuntarily childless. Fertile subjects exhibited normozoospermia and had fathered children without medical assistance.

RESULTS

We discovered proteins from about 7000 coding genes in the human sperm proteome. These were mainly known for involvements in cellular motility, response to stimuli, adhesion, and reproduction. Numbers of sperm proteins showing at least threefold deviating abundances increased from oligozoospermia (N = 153) and oligoasthenozoospermia (N = 154) to oligoasthenoteratozoospermia (N = 368). Deregulated sperm proteins primarily engaged in flagellar assembly and sperm motility, fertilization, and male gametogenesis. Most of these participated in a larger network of male infertility genes and proteins.

DISCUSSION

We expose 31 sperm proteins displaying deviant abundances under infertility, which already were known before to have fertility relevance, including ACTL9, CCIN, CFAP47, CFAP65, CFAP251 (WDR66), DNAH1, and SPEM1. We propose 18 additional sperm proteins with at least eightfold differential abundance for further testing of their diagnostic potential, such as C2orf16, CYLC1, SPATA31E1, SPATA31D1, SPATA48, EFHB (CFAP21), and FAM161A.

CONCLUSION

Our results shed light on the molecular background of the dysfunctionality of the fewer spermatozoa produced in oligozoospermia and syndromes including it. The male infertility network presented may prove useful in further elucidating the molecular mechanism of male infertility.

Common as well as unique methylation-sensitive DNA regulatory elements in three mammalian SLC9C1 genes

The SLC9C1 gene (which encodes the NHE10 protein) is essential for male fertility in both mice and humans, however the epigenetic mechanisms regulating its testis/sperm-specific gene expression have yet to be studied. Here we identify and characterize DNA regulatory elements of the SLC9C1 gene across three mammalian species: mouse, rat, and human. First, in silico analysis of these mammalian SLC9C1 genes identified a CpG island located upstream of the transcription start site in the same relative position in all three genes. Further analysis reveals that this CpG island behaves differently, with respect to gene regulatory activity, in the mouse SLC9C1 gene than it does in the rat and human SLC9C1 gene. The mouse SLC9C1 CpG island displays strong promoter activity by itself and seems to have a stronger gene regulatory effect than either the rat or human SLC9C1 CpG islands. While the function of the upstream SLC9C1 CpG island may be divergent across the three studied species, it appears that the promoters of these three mammalian SLC9C1 genes share similar DNA methylation-sensitive regulatory mechanisms. All three SLC9C1 promoter regions are differentially methylated in lung and testis, being more hypermethylated in lung relative to the testis, and DNA sequence alignments provide strong evidence of primary sequence conservation. Luciferase assays reveal that in vitro methylation of constructs containing different elements of the three SLC9C1 genes largely exhibit methylation-sensitive promoter activity (reduced promoter activity when methylated) in both HEK 293 and GC-1spg cells. In total, our data suggest that the DNA methylation-sensitive elements of the mouse, rat, and human SLC9C1 promoters are largely conserved, while the upstream SLC9C1 CpG island common to all three species seems to perform a different function in mouse than it does in rat and human. This work provides evidence that while homologous genes can all be regulated by DNA methylation-dependent epigenetic mechanisms, the location of the specific cis-regulatory elements responsible for this regulation can differ across species.

Constructing a seventeen-gene signature model for non-obstructive azoospermia based on integrated transcriptome analyses and WGCNA

BACKGROUND

Non-obstructive azoospermia (NOA) affects approximately 1% of the male population worldwide. The underlying mechanism and gene transcription remain unclear. This study aims to explore the potential pathogenesis for the detection and management of NOA.

METHODS

Based on four microarray datasets from the Gene Expression Omnibus database, integrated analysis and weighted correlation network analysis (WGCNA) were used to obtain the intersected common differentially expressed genes (DESs). Differential signaling pathways were identified via GO and GSVA-KEGG analyses. We constructed a seventeen-gene signature model using least absolute shrinkage and selection operation (LASSO) regression, and validated its efficacy in another two GEO datasets. Three patients with NOA and three patients with obstructive azoospermia were recruited. The mRNA levels of seven key genes were measured in testicular samples, and the gene expression profile was evaluated in the Human Protein Atlas (HPA) database.

RESULTS

In total, 388 upregulated and 795 downregulated common DEGs were identified between the NOA and control groups. ATPase activity, tubulin binding, microtubule binding, and metabolism- and immune-associated signaling pathways were significantly enriched. A seventeen-gene signature predictive model was constructed, and receiver operating characteristic (ROC) analysis showed that the area under the curve (AUC) values were 1.000 (training group), 0.901 (testing group), and 0.940 (validation set). The AUCs of seven key genes (REC8, CPS1, DHX57, RRS1, GSTA4, SI, and COX7B) were all > 0.8 in both the testing group and the validation set. The qRT-PCR results showed that consistent with the sequencing data, the mRNA levels of RRS1, GSTA4, and COX7B were upregulated, while CPS1, DHX57, and SI were downregulated in NOA. Four genes (CPS1, DHX57, RRS1, and SI) showed significant differences. Expression data from the HPA database showed the localization characteristics and trajectories of seven key genes in spermatogenic cells, Sertoli cells, and Leydig cells.

CONCLUSIONS

Our findings suggest a novel seventeen-gene signature model with a favorable predictive power, and identify seven key genes with potential as NOA-associated marker genes. Our study provides a new perspective for exploring the underlying pathological mechanism in male infertility.

Transcriptomic Analysis Revealed Candidate Genes Involved in Pseudomale Sperm Abnormalities in Chinese Tongue Sole (Cynoglossus semilaevis)

Chinese tongue sole (Cynoglossus semilaevis) has a ZZ/ZW sex determination system, but the genotypic female (ZW) can be sex-reversed into phenotypic males, namely, pseudomales. Pseudomale fish can produce only Z-type sperm but not W sperm. However, the molecular mechanism is unclear. To screen the key genes involved in pseudomale sperm abnormalities, we analysed the transcriptomic profiles of pseudomale and male sperm. In comparison to male sperm, 592 differentially expressed genes (DEGs) were identified in pseudomale sperm, including 499 upregulated and 93 downregulated genes. KEGG analysis indicated that the FoxO signalling pathway, especially the foxo3a and foxo6-like genes, may play an important role in spermatogenesis. The DEGs were mainly distributed on sex chromosomes, with 158 downregulated genes on the Z chromosome and 41 upregulated genes on the W chromosome. A specific area (14-15 M) on the Z chromosome was identified, which enriched eight DEGs inside the ~1 M region. In addition, there were five gene alleles on the sex chromosomes, which showed the opposite transcription pattern (upregulated for the W allele, downregulated for the Z allele). This study has provided valuable data for screening candidate genes involved in the pseudomale sperm abnormality.

A Comparative Cross-Platform Analysis to Identify Potential Biomarker Genes for Evaluation of Teratozoospermia and Azoospermia

Male infertility is a global public health concern. Teratozoospermia is a qualitative anomaly of spermatozoa morphology, contributing significantly to male infertility, whereas azoospermia is the complete absence of spermatozoa in the ejaculate. Thus, there is a serious need for unveiling the common origin and/or connection between both of these diseases, if any. This study aims to identify common potential biomarker genes of these two diseases via an in silico approach using a meta-analysis of microarray data. In this study, a differential expression analysis of genes was performed on four publicly available RNA microarray datasets, two each from teratozoospermia (GSE6872 and GSE6967) and azoospermia (GSE145467 and GSE25518). From the analysis, 118 DEGs were found to be common to teratozoospermia and azoospermia, and, interestingly, sperm autoantigenic protein 17 (SPA17) was found to possess the highest fold change value among all the DEGs (9.471), while coiled-coil domain-containing 90B (CCDC90B) and coiled-coil domain-containing 91 (CCDC91) genes were found to be common among three of analyses, i.e., Network Analyst, ExAtlas, and GEO2R. This observation indicates that SPA17, CCDC90B, and CCDC91 genes might have significant roles to play as potential biomarkers for teratozoospermia and azoospermia. Thus, our study opens a new window of research in this area and can provide an important theoretical basis for the diagnosis and treatment of both these diseases.

MicroRNA-targeting in male infertility: Sperm microRNA-19a/b-3p and its spermatogenesis related transcripts content in men with oligoasthenozoospermia

Objective: To elucidate and validate the potential regulatory function of miR-19a/b-3p and its spermatogenesis-related transcripts content in sperm samples collected from men with oligoasthenozoospermia.

Methods: Men presenting at an infertility clinic were enrolled. MicroRNA (miRNA) and target genes evaluation were carried out using in silico prediction analysis, Reverse transcription-quantitative PCR (RT-qPCR) validation, and Western blot confirmation.

Results: The expression levels of miRNA-19a/b-3p were significantly up-regulated and 51 target genes were significantly down-regulated in oligoasthenozoospermic men compared with age-matched normozoospermic men as determined by RT-qPCR. Correlation analysis highlighted that sperm count, motility, and morphology were negatively correlated with miRNA-19a/b-3p and positively correlated with the lower expression level of 51 significantly identified target genes. Furthermore, an inverse correlation between higher expression levels of miRNA-19a/b-3p and lower expression levels of 51 target genes was observed. Consistent with the results of the RT-qPCR, reduced expression levels of STK33 and DNAI1 protein levels were identified in an independent cohort of sperm samples collected from men with oligoasthenozoospermia.

Conclusion: Findings suggest that the higher expression of miRNA-19a/b-3p or the lower expression of target genes are associated with oligoasthenozoospermia and male infertility, probably through influencing basic semen parameters. This study lay the groundwork for future studies focused on investigating therapies for male infertility.

Transcriptome-Wide m6A Analysis Provides Novel Insights Into Testicular Development and Spermatogenesis in Xia-Nan Cattle

Testis is the primary organ of the male reproductive tract in mammals that plays a substantial role in spermatogenesis. Improvement of our knowledge regarding the molecular mechanisms in testicular development and spermatogenesis will be reflected in producing spermatozoa of superior fertility. Evidence showed that N6-Methyladenosine (m6A) plays a dynamic role in post-transcription gene expression regulation and is strongly associated with production traits. However, the role of m6A in bovine testis has not been investigated yet. In this study, we conducted MeRIP-Seq analysis to explore the expression profiles of the m6A and its potential mechanism underlying spermatogenesis in nine bovine testes at three developmental stages (prepuberty, puberty and postpuberty). The experimental animals with triplicate in each stage were chosen based on their semen volume and sperm motility except for the prepuberty bulls and used for testes collection. By applying MeRIP-Seq analysis, a total of 8,774 m6A peaks and 6,206 m6A genes among the studied groups were identified. All the detected peaks were found to be mainly enriched in the coding region and 3'- untranslated regions. The cross-analysis of m6A and mRNA expression exhibited 502 genes with concomitant changes in the mRNA expression and m6A modification. Notably, 30 candidate genes were located in the largest network of protein-protein interactions. Interestingly, four key node genes (PLK4, PTEN, EGR1, and PSME4) were associated with the regulation of mammal testis development and spermatogenesis. This study is the first to present a map of RNA m6A modification in bovine testes at distinct ages, and provides new insights into m6A topology and related molecular mechanisms underlying bovine spermatogenesis, and establishes a basis for further studies on spermatogenesis in mammals.

Altered Gene Expression in the Testis of Infertile Patients with Nonobstructive Azoospermia

Background

The molecular mechanism of nonobstructive azoospermia (NOA) remains unclear. The aim of this study was to identify gene expression changes in NOA patients and to explore potential biomarkers and therapeutic targets.

Methods

The gene expression profiles of GSE45885 and GSE145467 were collected from the Gene Expression Omnibus (GEO) database, and the differences between NOA and normal spermatogenesis were analyzed. Enrichment analysis was performed to explore biological functions for common differentially expressed genes (DEGs) in GSE45885 and GSE145467. Coexpression analysis of DEGs in GSE45885 was performed, and two modules with the highest correlation with NOA were screened. Key genes were then screened from the intersection genes of the two modules and common DEGs and PPI network. The expression of key genes was validated by quantitative real-time polymerase chain reaction (qRT-PCR) experiments. Finally, through miRTarBase, miRDB, and RAID, the miRNAs were predicted to regulate key genes, respectively.

Results

A total of 345 common DEGs were identified and they were mainly related to spermatogenesis, insulin signaling pathway. Coexpression analysis of DEGs in GSE45885 yielded eight modules; MEblack and MEturquoise had the highest correlation with NOA. Six genes in MEturquoise and RNF141 in MEblack were identified as key genes. qRT-PCR experiments validated the differential expression of key genes between NOA and control. Furthermore, RNF141 was regulated by the largest number of miRNAs.

Conclusion

Our findings suggest that the significant change expression of key genes may be potential markers and therapeutic targets of NOA and may have some impact on the development of NOA.