The Novel Key Genes of Non-obstructive Azoospermia Affect Spermatogenesis: Transcriptomic Analysis Based on RNA-Seq and scRNA-Seq Data

The Odad3 Gene Is Necessary for Spermatozoa Development and Male Fertility in Mice

Odad3 gene loss-of-function mutation leads to Primary Ciliary Dyskinesia (PCD), a disease caused by motile cilia dysfunction. Previously, we demonstrated that knockout of the Odad3 gene in mice replicates several features of PCD, such as hydrocephalus, defects in left-right body symmetry, and male infertility, with a complete absence of sperm in the reproductive tract. The majority of Odad3 knockout animals die before sexual maturation due to severe hydrocephalus and failure to thrive, which precludes fertility studies. Here, we performed the expression analysis of the Odad3 gene during gonad development and in adult testes. We showed that Odad3 starts its expression during the first wave of spermatogenesis, specifically at the meiotic stage, and that its expression is restricted to the germ cells in the adult testes, suggesting that Odad3 plays a role in spermatozoa formation. Subsequently, we conditionally deleted the Odad3 gene in adult males and demonstrated that even partial ablation of the Odad3 gene leads to asthenoteratozoospermia with multiple morphological abnormalities of sperm flagella (MMAF) in mice. The analysis of the seminiferous tubules in Odad3-deficient mice revealed defects in spermatogenesis with accumulation of seminiferous tubules at the spermiogenesis and spermiation phases. Furthermore, analysis of fertility in heterozygous Odad3+/- knockout mice revealed a reduction in sperm count and motility as well as abnormal sperm morphology. Additionally, Odad3+/- males exhibited a shorter fertile lifespan. Overall, these results suggest the important role of Odad3 and Odad3 gene dosage in male fertility. These findings may have an impact on the genetic and fertility counseling practice of PCD patients carrying Odad3 loss-of-function mutations.

Whole transcriptome analysis to identify non-coding RNA regulators and hub genes in sperm of non-obstructive azoospermia by microarray, single-cell RNA sequencing, weighted gene co-expression network analysis, and mRNA-miRNA-lncRNA interaction analysis

BACKGROUND

The issue of male fertility is becoming increasingly common due to genetic differences inherited over generations. Gene expression and evaluation of non-coding RNA (ncRNA), crucial for sperm development, are significant factors. This gene expression can affect sperm motility and, consequently, fertility. Understanding the intricate protein interactions that play essential roles in sperm differentiation and development is vital. This knowledge could lead to more effective treatments and interventions for male infertility.

MATERIALS AND METHODS

Our research aim to identify new and key genes and ncRNA involved in non-obstructive azoospermia (NOA), improving genetic diagnosis and offering more accurate estimates for successful sperm extraction based on an individual's genotype.

RESULTS

We analyzed the transcript of three NOA patients who tested negative for genetic sperm issues, employing comprehensive genome-wide analysis of approximately 50,000 transcript sequences using microarray technology. This compared gene expression profiles between NOA sperm and normal sperm. We found significant gene expression differences: 150 genes were up-regulated, and 78 genes were down-regulated, along with 24 ncRNAs up-regulated and 13 ncRNAs down-regulated compared to normal conditions. By cross-referencing our results with a single-cell genomics database, we identified overexpressed biological process terms in differentially expressed genes, such as "protein localization to endosomes" and "xenobiotic transport." Overrepresented molecular function terms in up-regulated genes included "voltage-gated calcium channel activity," "growth hormone-releasing hormone receptor activity," and "sialic acid transmembrane transporter activity." Analysis revealed nine hub genes associated with NOA sperm: RPL34, CYB5B, GOL6A6, LSM1, ARL4A, DHX57, STARD9, HSP90B1, and VPS36.

CONCLUSIONS

These genes and their interacting proteins may play a role in the pathophysiology of germ cell abnormalities and infertility.

Genetic profiling of azoospermic men to identify the etiology and predict reproductive potential

PURPOSE

To identify germline mutations related to azoospermia etiology and reproductive potential of surgically retrieved spermatozoa, and to investigate the feasibility of predicting seminiferous tubule function of nonobstructive azoospermic men by transcriptomic profiling of ejaculates.

MATERIALS AND METHODS

Sperm specimens were obtained from 30 men (38.4 ± 6 years) undergoing epididymal sperm aspiration for obstructive azoospermia (OA, n = 19) acquired by vasectomy, or testicular biopsy for nonobstructive azoospermia (NOA, n = 11). To evaluate for a correlation with azoospermia etiology, DNAseq was performed on surgically retrieved spermatozoa, and cell-free RNAseq on seminal fluid (n = 23) was performed to predict spermatogenesis in the seminiferous tubule.

RESULTS

Overall, surgically retrieved sperm aneuploidy rates were 1.7% and 1.8% among OA and NOA cohorts, respectively. OA men carried housekeeping-related gene mutations, while NOA men displayed mutations on genes involved in crucial spermiogenic functions (AP1S2, AP1G2, APOE). We categorized couples within each cohort according to ICSI clinical outcomes to investigate genetic causes that may affect reproductive potential. All OA-fertile men (n = 9) carried mutations in ZNF749 (sperm production), whereas OA-infertile men (n = 10) harbored mutations in PRB1, which is essential for DNA replication. NOA-fertile men (n = 8) carried mutations in MPIG6B (stem cell lineage differentiation), whereas NOA-infertile individuals (n = 3) harbored mutations in genes involved in spermato/spermio-genesis (ADAM29, SPATA31E1, MAK, POLG, IFT43, ATG9B) and early embryonic development (MBD5, CCAR1, PMEPA1, POLK, REC8, REPIN1, MAPRE3, ARL4C). Transcriptomic assessment of cell-free RNAs in seminal fluid from NOA men allowed the prediction of residual spermatogenic foci.

CONCLUSIONS

Sperm genome profiling provides invaluable information on azoospermia etiology and identifies gene-related mechanistic links to reproductive performance. Moreover, RNAseq assessment of seminal fluid from NOA men can help predict sperm retrieval during testicular biopsies.

Molecular mechanisms of cellular dysfunction in testes from men with non-obstructive azoospermia

Male factor infertility affects 50% of infertile couples worldwide; the most severe form, non-obstructive azoospermia (NOA), affects 10-15% of infertile males. Treatment for individuals with NOA is limited to microsurgical sperm extraction paired with in vitro fertilization intracytoplasmic sperm injection. Unfortunately, spermatozoa are only retrieved in ~50% of patients, resulting in live birth rates of 21-46%. Regenerative therapies could provide a solution; however, understanding the cell-type-specific mechanisms of cellular dysfunction is a fundamental necessity to develop precision medicine strategies that could overcome these abnormalities and promote regeneration of spermatogenesis. A number of mechanisms of cellular dysfunction have been elucidated in NOA testicular cells. These mechanisms include abnormalities in both somatic cells and germ cells in NOA testes, such as somatic cell immaturity, aberrant growth factor signalling, increased inflammation, increased apoptosis and abnormal extracellular matrix regulation. Future cell-type-specific investigations in identifying modulators of cellular transcription and translation will be key to understanding upstream dysregulation, and these studies will require development of in vitro models to functionally interrogate spermatogenic niche dysfunction in both somatic and germ cells.

Epididymal segment-specific miRNA and mRNA regulatory network at the single cell level

Spermatozoa released from the testis cannot fertilize an egg before becoming mature and motile in the epididymis. Based on three bulk and one single-cell RNA-seq (scRNA-seq) data series, we compared mRNA or miRNA expression between epididymal segment-specific samples and the other samples. Hereby, we identified 570 differentially expressed mRNAs (DE-mRNAs) and 23 differentially expressed miRNAs (DE-miRNAs) in the caput, 175 DE-mRNAs and 15 DE-miRNAs in the corpus, 946 DE-mRNAs and 12 DE-miRNAs in the cauda. In accordance with respective DE-miRNAs, we predicted upstream transcription factors (TFs) and downstream target genes. Subsequently, we intersected target genes of respective DE-miRNAs with corresponding DE-mRNAs, thereby obtaining 127 upregulated genes in the caput and 92 upregulated genes in cauda. Enriched upregulated pathways included cell motility-related pathways for the caput, smooth muscle-related pathways for the corpus, and immune-associated pathways for the cauda. Protein-protein interaction (PPI) network was constructed to extract key module for the caput and cauda, followed by identifying hub genes through cytohubba. Epididymis tissues from six mice were applied to validate hub genes expression using qRT-PCR, and 7 of the 10 genes displayed identical expression trends in mice caput/cauda. These hub genes were found to be predominantly distributed in spermatozoa using scRNA-seq data. In addition, target genes of DE-miRNAs were intersected with genes in the PPI network for each segment. Subsequently, the miRNA and mRNA regulatory networks for the caput and cauda were constructed. Conclusively, we uncover segment-specific miRNA-mRNA regulatory network, upstream TFs, and downstream pathways of the human epididymis, warranting further investigation into epididymal segment-specific functions.

Single-cell RNA sequencing technology in human spermatogenesis: Progresses and perspectives

Spermatogenesis, a key part of the spermiation process, is regulated by a combination of key cells, such as primordial germ cells, spermatogonial stem cells, and somatic cells, such as Sertoli cells. Abnormal spermatogenesis can lead to azoospermia, testicular tumors, and other diseases related to male infertility. The application of single-cell RNA sequencing (scRNA-seq) technology in male reproduction is gradually increasing with its unique insight into deep mining and analysis. The data cover different periods of neonatal, prepubertal, pubertal, and adult stages. Different types of male infertility diseases including obstructive and non-obstructive azoospermia (NOA), Klinefelter Syndrome (KS), Sertoli Cell Only Syndrome (SCOS), and testicular tumors are also covered. We briefly review the principles and application of scRNA-seq and summarize the research results and application directions in spermatogenesis in different periods and pathological states. Moreover, we discuss the challenges of applying this technology in male reproduction and the prospects of combining it with other technologies.

An artificial neural network model to diagnose non-obstructive azoospermia based on RNA-binding protein-related genes

Non-obstructive azoospermia (NOA) is a severe form of male infertility, but its pathological mechanisms and diagnostic biomarkers remain obscure. Since the dysregulation of RNA-binding proteins (RBPs) had nonnegligible effects on spermatogenesis, we aimed to investigate the functions and diagnosis values of RBPs in NOA. 58 testicular samples (control = 11, NOA = 47) from Gene Expression Omnibus (GEO) were set as the training cohort. Three public datasets, containing GSE45885 (control = 4, NOA = 27), GSE45887 (control = 4, NOA = 16), and GSE145467 (control = 10, NOA = 10), and 44 clinical samples from the local hospital (control = 27, NOA = 17) were used for validation. Through a series of bioinformatical analyses and machine learning algorithms, including genomic difference detection, protein-protein interaction network analysis, LASSO, SVM-RFE, and Boruta, DDX20 and NCBP2 were determined as significant predictors of NOA. Single-cell RNA sequencing of 432 testicular cell samples from NOA patients indicated that DDX20 and NCBP2 were associated with spermatogenesis (false discovery rate < 0.05). Based on the transcriptome expressions of DDX20 and NCBP2, we constructed multiple diagnosis models using logistic regression, random forest, and artificial neural network (ANN). The ANN model exhibited the most reliable predictive performance in the training cohort (AUC = 0.840), GSE45885 (AUC = 0.731), GSE45887 (AUC = 0.781), GSE145467 (AUC = 0.850), and local cohort (AUC = 0.623). Totally, an ANN diagnosis model based on RBP DDX20 and NCBP2 was developed and externally validated in NOA, functioning as a promising tool in clinical practice.

Deciphering the Molecular Characteristics of Human Idiopathic Nonobstructive Azoospermia from the Perspective of Germ Cells

Nonobstructive azoospermia (NOA) is one of the most important causes of male infertility, accounting for 10-15% of infertile men worldwide. Among these, more than 70% of cases are idiopathic NOA (iNOA), whose pathogenesis and molecular basis remain unknown. This work profiles 3696 human testicular single-cell transcriptomes from 17 iNOA patients, which are classified into four classes with different arrest periods and variable cell proportions based on the gene expression patterns and pathological features. Genes related to the cell cycle, energy production, and gamete generation show obvious abnormalities in iNOA germ cells. This work identifies several candidate causal genes for iNOA, including CD164, LELP1, and TEX38, which are significantly downregulated in iNOA germ cells. Notably, CD164 knockdown promotes apoptosis in spermatogonia. Cellular communications between spermatogonial stem cells and Sertoli cells are disturbed in iNOA patients. Moreover, BOD1L2, C1orf194, and KRTCAP2 are found to indicate testicular spermatogenic capacity in a variety of testicular diseases, such as Y-chromosome microdeletions and Klinefelter syndrome. In general, this study analyzes the pathogenesis of iNOA from the perspective of germ cell development, transcription factor (TF) regulatory networks, as well as germ cell and somatic cell interactions, which provides new ideas for clinical diagnosis.

Expression profiles of meiotic genes in male vs. female gonads and gametes: Insights into fertility issues

Gametes are specialized cells that, at fertilization, give rise to a totipotent zygote capable of generating an entire organism. Female and male germ cells undergo meiosis to produce mature gametes; however, sex-specific events of oogenesis and spermatogenesis contribute to specific roles of gametes in reproductive issues. We investigate the differential gene expression (DGE) of meiosis-related genes in human female and male gonads and gametes in normal and pathological conditions. The transcriptome data for the DGE analysis was obtained through the Gene Expression Omnibus repository, comprising human ovary and testicle samples of the prenatal period and adulthood, additionally to male (non-obstructive azoospermia (NOA) and teratozoospermia), and female (polycystic ovary syndrome (PCOS) and advanced maternal age) reproductive conditions. Gene ontology terms related to meiosis were associated with 678 genes, of which 17 genes in common were differentially expressed between the testicle and ovary during the prenatal period and adulthood. Except for SERPINA5 and SOX9, the 17 meiosis-related genes were downregulated in the testicle during the prenatal period and upregulated in adulthood compared to the ovary. No differences were observed in the oocytes of PCOS patients; however, meiosis-related genes were differentially expressed according to the patient’s age and maturity of the oocyte. In NOA and teratozoospermia, 145 meiosis-related genes were differentially expressed in comparison to the control, including OOEP; despite no recognized role in male reproduction, OOEP was co-expressed with genes related to male fertility. Taking together, these results shed light on potential genes that might be relevant to comprehend human fertility disorders.

Application of single-cell RNA sequencing on human testicular samples: a comprehensive review

So far there has been no comprehensive review using systematic literature search strategies to show the application of single-cell RNA sequencing (scRNA-seq) in the human testis of the whole life cycle (from embryos to aging males). Here, we summarized the application of scRNA-seq analyses on various human testicular biological samples. A systematic search was conducted in PubMed and Gene Expression Omnibus (GEO), focusing on English researches published after 2009. Articles related to GEO data-series were also retrieved in PubMed or BioRxiv. 81 full-length studies were finally included in the review. ScRNA-seq has been widely used on different human testicular samples with various library strategies, and new cell subtypes such as State 0 spermatogonial stem cells (SSC) and stage_a/b/c Sertoli cells (SC) were identified. For the development of normal testes, scRNA-seq-based evidence showed dynamic transcriptional changes of both germ cells and somatic cells from embryos to adults. And dysregulated metabolic signaling or hedgehog signaling were revealed by scRNA-seq in aged SC or Leydig cells (LC), respectively. For infertile males, scRNA-seq studies revealed profound changes of testes, such as the increased proportion of immature SC/LC of Klinefelter syndrome, the somatic immaturity and altered germline autophagy of patients with non-obstructive azoospermia, and the repressed differentiation of SSC in trans-females receiving testosterone inhibition therapy. Besides, the re-analyzing of public scRNA-seq data made further discoveries such as the potential vulnerability of testicular SARS-CoV-2 infection, and both evolutionary conservatism and divergence among species. ScRNA-seq analyses would unveil mechanisms of testes' development and changes so as to help developing novel treatments for male infertility.

Influence of the Postmortem/Storage Time of Human Corneas on the Properties of Cultured Limbal Epithelial Cells

Besides being a powerful model to study the mechanisms of corneal wound healing, tissue-engineered human corneas (hTECs) are sparking interest as suitable substitutes for grafting purposes. To ensure the histological and physiological integrity of hTECs, the primary cultures generated from human cornea (identified as human limbal epithelial cells (hLECs) that are used to produce them must be of the highest possible quality. The goal of the present study consisted in evaluating the impact of the postmortem/storage time (PM/ST) on their properties in culture. hLECs were isolated from the entire cornea comprising the limbus and central cornea. When grown as monolayers, short PM/ST hLECs displayed increased daily doublings and generated more colonies per seeded cells than long PM/ST hLECs. Moreover, hLECs with a short PM/ST exhibited a markedly faster wound closure kinetic both in scratch wound assays and hTECs. Collectively, these results suggest that short PM/ST hLECs have a greater number of highly proliferative stem cells, exhibit a faster and more efficient wound healing response in vitro, and produce hTECs of a higher quality, making them the best candidates to produce biomaterial substitutes for clinical studies.

Advances in studying human gametogenesis and embryonic development in China

Reproductive medicine in China has developed rapidly since 1988 due to the support from the government and scientific exploration. However, the success rate of assisted reproduction technology (ART) is around 30-40% and many unknown "black boxes" in gametogenesis and embryo development are still present. With the development of single-cell and low-input sequencing technologies, the network of transcriptome and epigenetic regulation (DNA methylation, chromatin accessibility, and histone modifications) during the development of human primordial germ cells (PGCs), gametes and embryos has been investigated in depth. Furthermore, pre-implantation genetic testing (PGT) has also rapidly developed. In this review, we summarize and analyze China's outstanding progress in these fields.

KIF2C affects sperm cell differentiation in patients with Klinefelter syndrome, as revealed by RNA-Seq and scRNA-Seq data

Klinefelter syndrome (KS) is a leading contributor to male infertility and is characterised by complex and diverse clinical features; however, genetic changes in the KS transcriptome remain largely unknown. We therefore used transcriptomic and single‐cell RNA sequencing (scRNA‐seq) datasets from KS versus normal populations through the Gene Expression Omnibus (GEO) database to identify gene biomarkers associated with the occurrence of KS. We identified a total of 700 differentially expressed genes (DEGs) and completed Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), enrichment pathway analysis and protein‐protein interaction (PPI) network analysis. A total of four unreported KS‐related hub genes (KIF2C, MRPS2, RPS15 and TSFM) were identified. Validation of the single‐cell sequencing dataset showed that only KIF2C and RPS15 were expressed in spermatocytes and that they were differentially expressed in sperm cells. Further construction of the developmental trajectories of these two genes in sperm cells showed that the KIF2C gene showed an upward trend throughout the differentiation and development of sperm cells. In conclusion, we report here that KIF2C may be closely related to the differentiation and development of sperm cells in KS patients, which is important for revealing the molecular mechanism of KS and conducting further studies.

Construction and external validation of a 5-gene random forest model to diagnose non-obstructive azoospermia based on the single-cell RNA sequencing of testicular tissue

Non-obstructive azoospermia (NOA) is among the most severe factors for male infertility, but our understandings of the latent biological mechanisms remain insufficient. The single-cell RNA sequencing (scRNA-seq) data of 432 testicular cells isolated from the patient with NOA was analyzed, and the cell samples were grouped into 5 cell clusters. A sum of 455 cell markers was identified and then included in the protein-protein interaction network. The Top 5 most critical genes in the network, including CCT8, CDC6, PSMD1, RPS4X, RPL36A, were selected for the diagnosis model construction through the random forest (RF). The RF model was a strong classifier for NOA and obstructive azoospermia (OA), which was validated in the training cohort (n = 58, AUC = 1) and external validation cohort (n = 20, AUC = 0.9). We collected the seminal plasma samples and testicular biopsy samples from 20 OA and 20 NOA cases from the local hospital, and the gene expression was detected via Real-Time quantitative Polymerase Chain Reaction (RT-qPCR) and Immunohistochemistry. The RF model also exhibited high accuracy (AUC = 0.725) in the local cohort. In summary, a novel gene signature was developed and externally validated based on scRNA-seq analysis, providing some new biomarkers to uncover the underlying mechanisms and a promising clinical tool for diagnosis in NOA.