Potential biomarkers of nonobstructive azoospermia identified in microarray gene expression analysis. Fertil Steril. 2013;100:1686-94.e1-7. [PMID: 24012201 DOI: 10.1016/j.fertnstert.2013.07.1999]

Genomic characterization of dryland adaptation in endangered Anxi cattle in China

Local species exhibit distinctive indigenous characteristics while showing unique productive and phenotypic traits. However, the advent of commercialization has posed a substantial threat to the survival of indigenous species. Anxi cattle, an endangered native breed in China, have evolved unique growth and reproductive characteristics in extreme desert and semidesert ecosystems. In this study, we conducted a genomic comparison of 10 Anxi cattle genomes with those of five other global populations/breeds to assess genetic diversity and identify candidate genomic regions in Anxi cattle. Population structure and genetic diversity analyses revealed that Anxi cattle are part of the East Asian cattle clade, exhibiting higher genetic diversity than commercial breeds. Through selective sweep analysis, we identified specific genetic variations linked to the environmental adaptability of Anxi cattle. Notably, we identified several candidate genes, including CERS3 involved in regulating skin permeability and antimicrobial functions, RBFOX2 associated with cardiac development, SLC16A7 participated in the regulation of pancreatic endocrine function, and SPATA3 related to reproduction. Our findings revealed the distinctive genomic features of Anxi cattle in dryland environments, provided invaluable insights for further research and breed preservation, and had important significance for enriching the domestic cattle breeding gene bank.

Transcriptional Regulation Analysis Provides Insight into the Function of GSK3β Gene in Diannan Small-Ear Pig Spermatogenesis

Glycogen synthase kinase-3β (GSK3β) not only plays a crucial role in regulating sperm maturation but also is pivotal in orchestrating the acrosome reaction. Here, we integrated single-molecule long-read and short-read sequencing to comprehensively examine GSK3β expression patterns in adult Diannan small-ear pig (DSE) testes. We identified the most important transcript ENSSSCT00000039364 of GSK3β, obtaining its full-length coding sequence (CDS) spanning 1263 bp. Gene structure analysis located GSK3β on pig chromosome 13 with 12 exons. Protein structure analysis reflected that GSK3β consisted of 420 amino acids containing PKc-like conserved domains. Phylogenetic analysis underscored the evolutionary conservation and homology of GSK3β across different mammalian species. The evaluation of the protein interaction network, KEGG, and GO pathways implied that GSK3β interacted with 50 proteins, predominantly involved in the Wnt signaling pathway, papillomavirus infection, hippo signaling pathway, hepatocellular carcinoma, gastric cancer, colorectal cancer, breast cancer, endometrial cancer, basal cell carcinoma, and Alzheimer's disease. Functional annotation identified that GSK3β was involved in thirteen GOs, including six molecular functions and seven biological processes. ceRNA network analysis suggested that DSE GSK3β was regulated by 11 miRNA targets. Furthermore, qPCR expression analysis across 15 tissues highlighted that GSK3β was highly expressed in the testis. Subcellular localization analysis indicated that the majority of the GSK3β protein was located in the cytoplasm of ST (swine testis) cells, with a small amount detected in the nucleus. Overall, our findings shed new light on GSK3β's role in DSE reproduction, providing a foundation for further functional studies of GSK3β function.

Dysregulation of MTFR2, ATP5IF1 and BAK1 in Sertoli cells relates to idiopathic non-obstructive azoospermia via inhibiting mitochondrial fission and inducing mitochondrial dysfunction†

Non-obstructive azoospermia affects more than 10% of infertile men with over 70% patients are idiopathic with uncharacterized molecular mechanisms, which is referred as idiopathic non-obstructive azoospermia. In this study, we checked the morphology of Sertoli cell mitochondria in testis biopsies from patients with idiopathic non-obstructive azoospermia and patients with obstructive azoospermia who have normal spermiogenesis. The expression of 104 genes controlling mitochondria fission and fusion were analyzed in three gene expression datasets including a total of 60 patients with non-obstructive azoospermia. The levels of 7 candidate genes were detected in testis biopsies from 38 patients with idiopathic non-obstructive azoospermia and 24 patients with obstructive azoospermia who have normal spermatogenesis by RT-qPCR. Cell viability, apoptosis, mitochondria membrane potential, adenosine triphosphate production, oxygen consumption, and mitochondria morphology were examined in primary human Sertoli cells. Mouse spermatogonial stem cells were used to detect the cell supporting capacity of Sertoli cells. We observed that patients with idiopathic non-obstructive azoospermia had elongated mitochondria. MTFR2 and ATP5IF1 were downregulated, whereas BAK1 was upregulated in idiopathic non-obstructive azoospermia testis and Sertoli cells. Sertoli cells from patients with idiopathic non-obstructive azoospermia had reduced viability, mitochondria membrane potential, adenosine triphosphate production, oxygen consumption rate, glycolysis and increased apoptosis. Knockdown MTFR2 in Sertoli cells increased the mitochondria size. Knockdown ATP5IF1 did not change mitochondrial morphology but increased adenosine triphosphate hydrolysis. Overexpression of BAK1 reduced membrane potential and upregulated cell apoptosis. The dysregulation of all these three genes contributed to the dysfunction of Sertoli cells, which provides a clue for idiopathic non-obstructive azoospermia treatment.

TMEM225 Is Essential for Sperm Maturation and Male Fertility by Modifying Protein Distribution of Sperm in Mice

Nonobstructive azoospermia is the leading cause of male infertility. Abnormal levels of transmembrane protein 225 (TMEM225), a testis-specific protein, have been found in patients with nonobstructive azoospermia, suggesting that TMEM225 plays an essential role in male fertility. Here, we generated a Tmem225 KO mouse model to explore the function and mechanism of TMEM225 in male reproduction. Male Tmem225 KO mice were infertile. Surprisingly, Tmem225 deletion did not affect spermatogenesis, but TMEM225-null sperm exhibited abnormalities during epididymal maturation, resulting in reduced sperm motility and an abnormal hairpin-loop configuration. Furthermore, proteomics analyses of cauda sperm revealed that signaling pathways related to mitochondrial function, the glycolytic pathway, and sperm flagellar morphology were abnormal in Tmem225 KO sperm, and spermatozoa lacking TMEM225 exhibited high reactive oxygen species levels, reduced motility, and flagellar folding, leading to typical asthenospermia. These findings suggest that testicular TMEM225 may control the sperm maturation process by regulating the expression of proteins related to mitochondrial function, glycolysis, and sperm flagellar morphology in epididymal spermatozoa.

Comparative genomic survey and functional analysis of DKKL1 during spermatogenesis in the Chinese soft-shelled turtle (Pelodiscus sinensis)

A feature of the Chinese soft-shelled turtle (Pelodiscus sinensis) is seasonal spermatogenesis; however, the underlying molecular mechanism is not well clarified. Here, we firstly cloned and characterized P. sinensis DKKL1, and then performed comparative genomic studies, expression analysis, and functional validation. P. sinensis DKKL1 had 2 putative N-glycosylation sites and 16 phosphorylation sites. DKKL1 also had classic transmembrane structures that were extracellularly localized. DKKL1's genetic distance was close to turtles, followed by amphibians and mammals, but its genetic distance was far from fishes. DKKL1 genes from different species shared distinct genomic characteristics. Meanwhile, they were also relatively conserved among themselves, at least from the perspective of classes. Notably, the transcription factors associated with spermatogenesis were also identified, containing CTCF, EWSR1, and FOXL2. DKKL1 exhibited sexually dimorphic expression only in adult gonads, which was significantly higher than that in other somatic tissues (P < 0.001), and was barely expressed in embryonic gonads. DKKL1 transcripts showed a strong signal in sperm, while faint signals were detected in other male germ cells. DKKL1 in adult testes progressively increased per month (P < 0.05), displaying a seasonal expression trait. DKKL1 was significantly downregulated in testes cells after the sex hormones (17β-estradiol and 17α-methyltestosterone) and Wnt/β-catenin inhibitor treatment (P < 0.05). Likewise, the Wnt/β-catenin inhibitor treatment dramatically repressed CTCF, EWSR1, and FOXL2 expression. Conversely, they were markedly upregulated after the 17β-estradiol and 17α-methyltestosterone treatment, suggesting that the three transcription factors might bind to different promoter regions, thereby negatively regulating DKKL1 transcription in response to the changes in the estrogen and androgen pathways, and positively controlling DKKL1 transcription in answer to the alterations in the Wnt/β-catenin pathway. Knockdown of DKKL1 significantly reduced the relative expression of HMGB2 and SPATS1 (P < 0.01), suggesting that it may be involved in seasonal spermatogenesis of P. sinensis through a positive regulatory interaction with these two genes. Overall, our findings provide novel insights into the genome evolution and potential functions of seasonal spermatogenesis of P. sinensis DKKL1.

Identification of ferroptotic genes and phenotypes in idiopathic nonobstructive azoospermia

Effective treatments for nonobstructive azoospermia (NOA), which affects 1% of all men globally, are limited by undefined pathogenic mechanisms, especially in idiopathic NOA (iNOA). Here, we tried to identify the functional ferroptosis-related genes and phenotypes involved in iNOA. Differentially expressed ferroptotic genes were identified from iNOA mRNA microarray datasets by bioinformatic analyses, and these ferroptotic genes were subsequently filtered by various algorithms. Then, receiver operating characteristic (ROC) curves were generated to evaluate the diagnostic ability of the abovementioned genes for iNOA. Generally, 11 differentially expressed ferroptotic genes were downregulated, and five genes were upregulated in iNOA samples. Four genes, including DUSP1, GPX4, HSD17B11, and SLC2A8, were technically selected and determined to be potential biomarkers for iNOA. Subsequently, similar expression levels were validated at both the RNA and protein levels in the iNOA specimens. Finally, morphologic and biochemical assays were applied to define the ferroptotic phenotypes in testes. The ferroptotic features, like shrunken mitochondria with electron-dense membranes and a reduction in cristae were observed across various cell types within iNOA patients, accompanied by the overload of ferrous ions and increased lipid peroxidation production. Our findings demonstrated that these ferroptosis genes could be involved in the underlying pathogenesis mechanisms of iNOA by regulating ferroptosis and serve as potential diagnostic biomarkers. Also, the ferroptotic phenotypes were identified in iNOA patients.

Acquired Male Hypogonadism in the Post-Genomic Era-A Narrative Review

Although precision medicine took its first steps from genomic medicine, it has gone far beyond genomics, considering the full complexity of cellular physiology. Therefore, the present time can be considered as the "post-genomic era". In detail, proteomics captures the overall protein profile of an analyzed sample, whilst metabolomics has the purpose of studying the molecular aspects of a known medical condition through the measurement of metabolites with low molecular weight in biological specimens. In this review, the role of post-genomic platforms, namely proteomics and metabolomics, is evaluated with a specific interest in their application for the identification of novel biomarkers in male hypogonadism and in the identification of new perspectives of knowledge on the pathophysiological function of testosterone. Post-genomic platforms, including MS-based proteomics and metabolomics based on ultra-high-performance liquid chromatography-HRMS, have been applied to find solutions to clinical questions related to the diagnosis and treatment of male hypogonadism. In detail, seminal proteomics helped us in identifying novel non-invasive markers of androgen activity to be translated into clinical practice, sperm proteomics revealed the role of testosterone in spermatogenesis, while serum metabolomics helped identify the different metabolic pathways associated with testosterone deficiency and replacement treatment, both in patients with insulin sensitivity and patients with insulin resistance.

Identification of a novel pyroptosis-related gene signature in human spermatogenic dysfunction

PURPOSE

To reveal the underlying roles that pyroptosis-related genes (PRGs) played in human spermatogenic dysfunction.

METHODS

One discovery set and three validation sets were employed to inspect the previously reported 33 PRGs in the human testis with different status of spermatogenesis. PRGs that differentially expressed in all sets were considered as key differentially expressed pyroptosis-related genes (PR-DEGs). The relationships between key PR-DEGs and samples' clinicopathological, therapeutic, and immune patterns were respectively studied. Single-cell RNA sequencing (scRNS-seq) analyses were conducted to show the expression changes and related mechanisms of key PR-DEGs at a single-cell resolution.

RESULTS

CASP4 and GPX4 were identified as two key PR-DEGs. These two genes were significantly dysregulated in spermatogenic dysfunctional samples, but with opposite tendency. CASP4 was negatively correlated with Johnsen scores but positively correlated with follicle-stimulating hormone (FSH) levels (all p < 0.05), while GPX4 exhibited significant positive correlations with Johnsen scores and negative relevance with FSH. For treatments, both molecules showed a prospective value of being predictors for sperm retrieval surgeries. Moreover, CASP4 and GPX4 were potential immunoregulators in the testicular immune microenvironment and showed significant correlations to testicular macrophages and mast cell infiltration. In scRNA-seq analyses, GPX4 was highly expressed in germ cells, which therefore suffered a sharp reduction with the loss of germ cells in spermatogenic dysfunction. On the other hand, CASP4 were basically somatic cell-derived, and the proportion of CASP4-positive Leydig cells significantly increased in disease testes (p = 0.0001).

CONCLUSION

In all, we revealed two key PRGs of human testes that might be functional in spermatogenic dysfunction.

Integrated molecular-network analysis reveals infertility-associated key genes and transcription factors in the non-obstructive azoospermia

BACKGROUND

Male infertility is a multifactorial reproductive health problem with complex causes. Non-obstructive azoospermia (NOA) is characterized by failure of spermatogenesis, leading to the absence of spermatozoa in ejaculates. The molecular mechanism underlying the NOA is still not well understood.

OBJECTIVES

This study aims to identify the key genes involved in male infertility that could be a potential biomarker in the diagnosis and prognosis of azoospermia.

STUDY DESIGN

The microarray expression profiles dataset GSE45885 and GSE45887 were downloaded from the NCBI's Gene Expression Omnibus (GEO) database and analyzed for male infertility-associated differentially expressed genes (DEGs) using the GEO2R tool. The common DEGs between the two datasets were combined and their protein-protein interaction (PPI) network was constructed using Cytoscape to reveal the hub genes by topology and module analysis. In addition, transcription factors (TFs) and protein kinases regulating the hub genes were identified using the X2K tool. Then, the expression of the hub genes was validated by analyzing the GSE190752 microarray dataset. Further, the PPI network was screened for biological roles and enriched pathways using DAVID software.

RESULTS

About 256 DEGs associated with NOA were identified and constructed the PPI network to find the infertility-associated proteins. The biological processes linked with these proteins were spermatogenesis, cell differentiation, flagellated sperm motility, and spermatid development. The topology and module analysis of the infertility-associated protein network identified the hub genes TEX38, FAM71F, PRR30, FAM166A, LYZL6, TPPP2, ARMC12, SPACA4, and FAM205A, which were found to be upregulated in the non-obstructive azoospermia. In addition, a total of 23 transcription factors and 3 protein kinases that are regulating these key hub genes were identified. Further these hub genes expression was validated using the microarray data and found that their expression was increased in the testicular biopsies obtained from NOA subjects, compared to healthy individuals.

CONCLUSION

The identified key genes and its associated transcription factors are known to regulate the infertility-related processes in the non-obstructive azoospermia. Also, the clinical sample-based microarray data validation for the expression of these key hub genes indicates their potentiality to develop them as diagnostic or prognostic biomarkers for NOA.

Spermatogenic cell-specific SPACA4 is essential for efficient sperm-zona pellucida binding in vitro

Fertilization is a complex and highly regulated process that involves a series of molecular interactions between sperm and oocytes. However, the mechanisms of proteins involved in human fertilization, such as that of testis-specific SPACA4, remain poorly understood. Here we demonstrated that SPACA4 is a spermatogenic cell-specific protein. SPACA4 is expressed during spermatogenesis, upregulated in early-stage spermatids, and downregulated in elongating spermatids. SPACA4 is an intracellular protein that locates in the acrosome and is lost during the acrosome reaction. Incubation with antibodies against SPACA4 inhibited the binding of spermatozoa to zona pellucida. SPACA4 protein expression levels across different semen parameters were similar but varied significantly among patients. A prospective clinical study found no association between SPACA4 protein levels and fertilization or cleavage rates. Thus, the study suggests a novel function for SPACA4 in human fertilization in a non-dose-dependent manner. However, a larger clinical trial is required to evaluate the potential use of sperm SPACA4 protein levels to predict fertilization potential.

Identification and validation of CCL2 as a potential biomarker relevant to mast cell infiltration in the testicular immune microenvironment of spermatogenic dysfunction

BACKGROUND

Spermatogenic dysfunction is an important cause of azoospermia. Numerous studies have focused on germ-cell-related genes that lead to spermatogenic impairment. However, based on the immune-privileged characteristics of the testis, the relationship of immune genes, immune cells or immune microenvironment with spermatogenic dysfunction has rarely been reported.

RESULTS

Using integrated methods including single-cell RNA-seq, microarray data, clinical data analyses and histological/pathological staining, we found that testicular mast cell infiltration levels were significantly negatively related to spermatogenic function. We next identified a functional testicular immune biomarker, CCL2, and externally validated that testicular CCL2 was significantly upregulated in spermatogenic dysfunctional testes and was negatively correlated with Johnsen scores (JS) and testicular volumes. We also demonstrated that CCL2 levels showed a significant positive correlation with testicular mast cell infiltration levels. Moreover, we showed myoid cells and Leydig cells were two of the important sources of testicular CCL2 in spermatogenic dysfunction. Mechanistically, we drew a potential "myoid/Leydig cells-CCL2-ACKR1-endothelial cells-SELE-CD44-mast cells" network of somatic cell-cell communications in the testicular microenvironment, which might play roles in spermatogenic dysfunction.

CONCLUSIONS

The present study revealed CCL2-relevant changes in the testicular immune microenvironment in spermatogenic dysfunction, providing new evidence for the role of immunological factors in azoospermia.

Machine learning and integrative analysis identify the common pathogenesis of azoospermia complicated with COVID-19

Background

Although more recent evidence has indicated COVID-19 is prone to azoospermia, the common molecular mechanism of its occurrence remains to be elucidated. The aim of the present study is to further investigate the mechanism of this complication.

Methods

To discover the common differentially expressed genes (DEGs) and pathways of azoospermia and COVID-19, integrated weighted co-expression network (WGCNA), multiple machine learning analyses, and single-cell RNA-sequencing (scRNA-seq) were performed.

Results

Therefore, we screened two key network modules in the obstructive azoospermia (OA) and non-obstructive azoospermia (NOA) samples. The differentially expressed genes were mainly related to the immune system and infectious virus diseases. We then used multiple machine learning methods to detect biomarkers that differentiated OA from NOA. Enrichment analysis showed that azoospermia patients and COVID-19 patients shared a common IL-17 signaling pathway. In addition, GLO1, GPR135, DYNLL2, and EPB41L3 were identified as significant hub genes in these two diseases. Screening of two different molecular subtypes revealed that azoospermia-related genes were associated with clinicopathological characteristics of age, hospital-free-days, ventilator-free-days, charlson score, and d-dimer of patients with COVID-19 (P < 0.05). Finally, we used the Xsum method to predict potential drugs and single-cell sequencing data to further characterize whether azoospermia-related genes could validate the biological patterns of impaired spermatogenesis in cryptozoospermia patients.

Conclusion

Our study performs a comprehensive and integrated bioinformatics analysis of azoospermia and COVID-19. These hub genes and common pathways may provide new insights for further mechanism research.

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.

Transcriptomic analysis of the Non-Obstructive Azoospermia (NOA) to address gene expression regulation in human testis

There is a need to understand the molecular basis of testes under Non-Obstructive Azoospermia (NOA), a state of failed spermatogenesis. There has been a lack of attention to the transcriptome at the level of alternatively spliced mRNAs (iso-mRNAs) and the mechanism of gene expression regulation. Hence, we aimed to establish a reliable iso-mRNA profile of NOA-testes, and explore molecular mechanisms - especially those related to gene expression regulation. We sequenced mRNAs from testicular samples of donors with complete spermatogenesis (control samples) and a failure of spermatogenesis (NOA samples). We identified differentially expressed genes and their iso-mRNAs via standard NGS data analyses. We then listed these iso-mRNAs hierarchically based on the extent of consistency of differential quantities across samples and groups, and validated the lists via RT-qPCRs (for 80 iso-mRNAs). In addition, we performed extensive bioinformatic analysis of the splicing features, domains, interactions, and functions of differentially expressed genes and iso-mRNAs. Many top-ranking down-regulated genes and iso-mRNAs, i.e., those down-regulated more consistently across the NOA samples, are associated with mitosis, replication, meiosis, cilium, RNA regulation, and post-translational modifications such as ubiquitination and phosphorylation. Most down-regulated iso-mRNAs correspond to full-length proteins that include all expected domains. The predominance of alternative promoters and termination sites in these iso-mRNAs indicate their gene expression regulation via promoters and UTRs. We compiled a new, comprehensive list of human transcription factors (TFs) and used it to identify TF-'TF gene' interactions with potential significance in down-regulating genes under the NOA condition. The results indicate that RAD51 suppression by HSF4 prevents SP1-activation, and SP1, in turn, could regulate multiple TF genes. This potential regulatory axis and other TF interactions identified in this study could explain the down-regulation of multiple genes in NOA-testes. Such molecular interactions may also have key regulatory roles during normal human spermatogenesis.

DDX58 expression promotes inflammation and growth arrest in Sertoli cells by stabilizing p65 mRNA in patients with Sertoli cell-only syndrome

Sertoli cell -only syndrome (SCOS) is a type of testicular pathological failure that causes male infertility and no effective treatment strategy, is available for this condition. Moreover, the molecular mechanism underlying its development remains unknown. We identified DExD/H-Box helicase 58 (DDX58) as a key gene in SCOS based on four datasets of testicular tissue samples obtained from the Gene Expression Synthesis database. DDX58 was significantly upregulated in SCOS testicular Sertoli cells. Moreover, high expression of DDX58 was positively correlated with the expression of several testicular inflammatory factors, such as IL -1β, IL-18, and IL-6. Interestingly, DDX58 could be induced in the D-galactose (D-gal)-stimulated TM4 cell injury model. Whereas silencing of DDX58 inhibited D-gal -mediated p65 expression, inflammatory cytokine release, and growth arrest. Mechanistically, we found that DDX58 acts as an RNA-binding protein, which enhances p65 expression by promoting mRNA stability. Furthermore, p65 gene silencing decreased the expression of inflammatory cytokines and inhibition of cell growth in D-gal-induced cells. In conclusion, our findings demonstrate that DDX58 promotes inflammatory responses and growth arrest in SCOS Sertoli cells by stabilizing p65 mRNA. Accordingly, the DDX58/p65 regulatory axis might be a therapeutic target for SCOS.

Whole-genome sequencing identifies new candidate genes for nonobstructive azoospermia

BACKGROUND

Genetic causes that lead to spermatogenetic failure in patients with nonobstructive azoospermia (NOA) have not been yet completely established.

OBJECTIVE

To identify low-frequency NOA-associated single nucleotide variants (SNVs) using whole-genome sequencing (WGS).

MATERIALS AND METHODS

Men with various types of NOA (n = 39), including samples that had been previously tested with whole-exome sequencing (WES; n = 6) and did not result in diagnostic conclusions. Variants were annotated using the Ensembl Variant Effect Predictor, utilizing frequencies from GnomAD and other databases to provide clinically relevant information (ClinVar), conservation scores (phyloP), and effect predictions (i.e., MutationTaster). Structural protein modeling was also performed.

RESULTS

Using WGS, we revealed potential NOA-associated SNVs, such as: TKTL1, IGSF1, ZFPM2, VCX3A (novel disease causing variants), ESX1, TEX13A, TEX14, DNAH1, FANCM, QRICH2, FSIP2, USP9Y, PMFBP1, MEI1, PIWIL1, WDR66, ZFX, KCND1, KIAA1210, DHRSX, ZMYM3, FAM47C, FANCB, FAM50B (genes previously known to be associated with infertility) and ALG13, BEND2, BRWD3, DDX53, TAF4, FAM47B, FAM9B, FAM9C, MAGEB6, MAP3K15, RBMXL3, SSX3 and FMR1NB genes, which may be involved in spermatogenesis.

DISCUSSION AND CONCLUSION

In this study, we identified novel potential candidate NOA-associated genes in 29 individuals out of 39 azoospermic males. Note that in 5 out of 6 patients subjected previously to WES analysis, which did not disclose potentially causative variants, the WGS analysis was successful with NOA-associated gene findings.

Chromosome 17 translocation affects sperm morphology: Two case studies and literature review

We present two cases of infertile males with teratozoospermia stemming from chromosome 17 translocation. The patients present karyotypes that have not been previously reported. Genes located on breakpoints (17p11.2, 9q31, and 11p15) were analysed to find the probable mechanism affecting sperm morphology. Our results suggest that ALKBH5, TOP3A, and LLGL1 interactions may be an underlying cause of abnormal sperm head morphology. Translocation of chromosome 17 occurred in conjunction with chromosome 9 and chromosome 11 translocation in the two cases, resulting in oligozoospermia and asthenozoospermia, respectively. These abnormal phenotypes may involve meiosis- and motility-related genes such as LDHC, DNHD1, UBQLN3, and NUP98. Translocation is thus a risk factor for sperm morphological abnormalities and motility deficiency. The interaction network of 22 genes on breakpoints suggests that they contribute to spermatogenesis as a group. In conclusion, this study highlighted the importance of investigating genes linked to sperm morphology, together with chromosome 17 translocation and reproductive risks. For patients interested in screening before a future pregnancy, we recommend preimplantation genetic diagnosis to reduce the risk of karyotypically unbalanced foetuses and birth defects.

Integrated analysis of miRNA and mRNA expression profiles in testes of Landrace and Hezuo boars

Precocious puberty is closely related to testicular development and spermatogenesis, and there is increasing evidence that miRNAs are involved in regulation of testicular development and spermatogenesis. However, little is known about the regulation of microRNAs (miRNAs) during precocious maturation in Hezuo (HZ) boars. In this study, serum Testosterone (T), Estradiol (E₂), Follicle-stimulating hormone (FSH) and Luteinizing hormone (LH) levels were detected in HZ and Landrace (LC) boars in the postnatal period at 30, 90, 120, 180, and 240 days, and the testes of HZ and LC boars at 30 and 120 days were used for histological observation. In addition, we performed small RNA-Seq to identify miRNA at sexual immaturity (30-days-old) and maturity (120-days-old) of HZ boar testis (using LC boar as control) to reveal the key miRNA in regulation of precocious puberty. Hormone assay results showed that high levels of T, E₂, FSH, and LH may be related to precocious sexual maturity of HZ boars, and that FSH may play an important function before sexual maturity. Histological observation showed that HZ boars developed earlier than LC boars and had reached sexual maturity at 120 days. Small RNA-Seq yielded a total of 359 exist miRNAs, 767 known miRNAs and 322 novel miRNAs in 12 samples; 549, 468, 133, and 247 differentially expressed (DE) miRNAs were identified between Ha vs. Hb, La vs. Lb, Ha vs. La, and Hb vs. Lb (log₂ fold change >1 and p < 0.05). Enrichment analysis showed that target genes of these DE miRNAs were enriched in many gene ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) signaling pathways (such as PI3K-Akt, Hippo and Rap1 signaling pathways) were related to testicular development and spermatogenesis. Further screening, some miRNAs (such as ssc-miR-29b, ssc-miR-199b, ssc-miR-383, ssc-miR-149, ssc-miR-615, and ssc-miR-370) were possibly associated with precocious puberty. These results provide new light on miRNA regulatory mechanisms involved in precocious puberty.

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.

Aberrant Gene Expression Profiling in Men With Sertoli Cell-Only Syndrome

Sertoli cell-only syndrome (SCOS) is the most severe and common pathological type of non-obstructive azoospermia. The etiology of SCOS remains largely unknown to date despite a handful of studies reported in this area. According to the gene expression of testicular tissue samples in six datasets from the Gene Expression Omnibus, we detected 1441 differentially expressed genes (DEGs) between SCOS and obstructive azoospermia (OA) testicular tissue samples. Enriched GO terms and KEGG pathways for the downregulated genes included various terms and pathways related to cell cycle and reproduction, while the enrichment for the upregulated genes yielded many inflammation-related terms and pathways. In accordance with the protein-protein interaction (PPI) network, all genes in the most critical module belonged to the downregulated DEGs, and we obtained nine hub genes, including CCNB1, AURKA, CCNA2, BIRC5, TYMS, UBE2C, CDC20, TOP2A, and OIP5. Among these hub genes, six were also found in the most significant SCOS-specific module obtained from consensus module analysis. In addition, most of SCOS-specific modules did not have a consensus counterpart. Based on the downregulated genes, transcription factors (TFs) and kinases within the upstream regulatory network were predicted. Then, we compared the difference in infiltrating levels of immune cells between OA and SCOS samples and found a significantly higher degree of infiltration for most immune cells in SCOS than OA samples. Moreover, CD56bright natural killer cell was significantly associated with six hub genes. Enriched hallmark pathways in SCOS had remarkably more upregulated pathways than the downregulated ones. Collectively, we detected DEGs, significant modules, hub genes, upstream TFs and kinases, enriched downstream pathways, and infiltrated immune cells that might be specifically implicated in the pathogenesis of SCOS. These findings provide new insights into the pathogenesis of SCOS and fuel future advances in its theranostics.

The Role of Testosterone in Spermatogenesis: Lessons From Proteome Profiling of Human Spermatozoa in Testosterone Deficiency

Testosterone is essential to maintain qualitative spermatogenesis. Nonetheless, no studies have been yet performed in humans to analyze the testosterone-mediated expression of sperm proteins and their importance in reproduction. Thus, this study aimed to identify sperm protein alterations in male hypogonadism using proteomic profiling. We have performed a comparative proteomic analysis comparing sperm from fertile controls (a pool of 5 normogonadic normozoospermic fertile men) versus sperm from patients with secondary hypogonadism (a pool of 5 oligozoospermic hypogonadic patients due to isolated LH deficiency). Sperm protein composition was analyzed, after peptide labelling with Isobaric Tags, via liquid chromatography followed by tandem mass spectrometry (LC-MS/MS) on an LTQ Velos-Orbitrap mass spectrometer. LC-MS/MS data were analyzed using Proteome Discoverer. Criteria used to accept protein identification included a false discovery rate (FDR) of 1% and at least 1 peptide match per protein. Up to 986 proteins were identified and, of those, 43 proteins were differentially expressed: 32 proteins were under-expressed and 11 were over-expressed in the pool of hypogonadic patients compared to the controls. Bioinformatic analyses were performed using UniProt Knowledgebase, and the Gene Ontology Consortium database based on PANTHER. Notably, 13 of these 43 differentially expressed proteins have been previously reported to be related to sperm function and spermatogenesis. Western blot analyses for A-Kinase Anchoring Protein 3 (AKAP3) and the Prolactin Inducible Protein (PIP) were used to confirm the proteomics data. In summary, a high-resolution mass spectrometry-based proteomic approach was used for the first time to describe alterations of the sperm proteome in secondary male hypogonadism. Some of the differential sperm proteins described in this study, which include Prosaposin, SMOC-1, SERPINA5, SPANXB1, GSG1, ELSPBP1, fibronectin, 5-oxoprolinase, AKAP3, AKAP4, HYDIN, ROPN1B, ß-Microseminoprotein and Protein S100-A8, could represent new targets for the design of infertility treatments due to androgen deficiency.

Testis-specific fascin component FSCN3 is dispensable for mouse spermatogenesis and fertility

BACKGROUND

Fascins belong to a family of actin-bundling proteins that are involved in a wide range of biological functions. FSCN3, a newly identified testis-specific actin-bundling protein, is specifically expressed in elongated spermatids. However, its in vivo function in mouse spermiogenesis remains unknown.

METHODS AND RESULTS

We generated Fscn3 knockout mice through CRISPR/Cas9 gene-editing technology. Fscn3^-/- mice displayed normal testis morphology and testis to bodyweight ratio, and sperm concentrations did not differ significantly between Fscn3^+/+ and Fscn3^-/- mice. Fertility assays consistently revealed that Fscn3^-/- mice are completely fertile and their reproductive status does not differ from that of wild-type. Moreover, hematoxylin and eosin staining of the testis sections of Fscn3^-/- mice detected various germ cells, ranging from spermatogonia to mature spermatozoa. Furthermore, the swimming velocity of the sperm of Fscn3^-/- mice was comparable to that of their wild-type littermates. Both Fscn3^+/+ and Fscn3^-/-mice had normal sperm morphology, indicating that the disruption of Fscn3 does not affect sperm morphology. The analysis of meiotic prophase I progression demonstrated normal prophase-I phases (leptonema to diplonema) in both Fscn3^+/+ and Fscn3^-/- mice, suggesting that Fscn3 is not essential for meiosis I.

CONCLUSION

Our study provides the first evidence that FSCN3 is a testis-specific actin-bundling protein that is not required for mouse spermatogenesis. Our results will help reproductive biologists focus their efforts on genes that are crucial for fertility and avoid research duplication.

Y-chromosome genes associated with sertoli cell-only syndrome identified by array comparative genome hybridization

BACKGROUND

The precise contribution of each chromosome gene or gene family in achieving male fertility is still the subject of debate. Most studies have examined male populations with heterogeneous causes of infertility, and have therefore reached controversial or uncertain conclusions. This study uses Y-chromosome array-based comparative genomic hybridization (aCGH) to examine a population of males with a uniform sertoli cell-only syndrome (SCOS) infertility phenotype.

METHODS

Initial analysis of gene copy number variations in 8 SCOS patients, with determination of the log-ratio of probe signal intensity against a DNA reference, was performed using the Y-chromosome NimbleGen aCGH. To confirm the role of candidate genes, real-time quantitative RT-PCR was used to compare 19 patients who had SCOS non-obstructive azoospermia with 15 patients who had obstructive azoospermia but normal spermatogenesis.

RESULTS

Our initial aCGH experiments identified CDY1a and CDY1b double deletions in all 8 patients who had total germ cell depletion. However, 5 patients had DAZ1/2 and DAZ3/4 deletions, 1 patient had a DAZ2 and DAZ3/4 deletion, and 2 patients had no DAZ1/2 or DAZ3/4 deletions. Examination of testicular mRNA expression in another 19 patients with SCOS indicated all patients had no detectable levels of CDY1.

CONCLUSIONS

Our findings indicate that CDY1 deletion in SCOS patients, and analysis of the expression of DAZ and CDY1 genes using aCGH and quantitative RT-PCR, may be useful to predict the presence of mature spermatozoa.

Identification and characterization of circular RNAs in the testicular tissue of patients with non-obstructive azoospermia

Circular RNAs (circRNAs) are highly conserved and ubiquitously expressed noncoding RNAs that participate in multiple reproduction-related diseases. However, the expression pattern and potential functions of circRNAs in the testes of patients with non-obstructive azoospermia (NOA) remain elusive. In this study, according to a circRNA array, a total of 37 881 circRNAs were identified that were differentially expressed in the testes of NOA patients compared with normal controls, including 19 874 upregulated circRNAs and 18 007 downregulated circRNAs. Using quantitative real-time polymerase chain reaction (qRT-PCR) analysis, we confirmed that the change tendency of some specific circRNAs, including hsa_circ_0137890, hsa_circ_0136298, and hsa_circ_0007273, was consistent with the microarray data in another larger sample. The structures and characteristics of these circRNAs were confirmed by Sanger sequencing, and fluorescence in situ hybridization revealed that these circRNAs were primarily expressed in the cytoplasm. Bioinformatics analysis was used to construct the competing endogenous RNA (ceRNA) network, and numerous miRNAs that could be paired with circRNAs validated in this study were reported to be vital for spermatogenesis regulation. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses indicated that genes involved in axoneme assembly, microtubule-based processes, and cell proliferation were significantly enriched. Our data suggest that there are aberrantly expressed circRNA profiles in patients with NOA and that these circRNAs may help identify key diagnostic and therapeutic molecular biomarkers for NOA patients.

Using publicly available transcriptomic data to identify mechanistic and diagnostic biomarkers in azoospermia and overall male infertility

Azoospermia, which is the absence of spermatozoa in an ejaculate occurring due to defects in sperm production, or the obstruction of the reproductive tract, affects about 1% of all men and is prevalent in up to 10–15% of infertile males. Conventional semen analysis remains the gold standard for diagnosing and treating male infertility; however, advances in molecular biology and bioinformatics now highlight the insufficiency thereof. Hence, the need to widen the scope of investigating the aetiology of male infertility stands pertinent. The current study aimed to identify common differentially expressed genes (DEGs) that might serve as potential biomarkers for non-obstructive azoospermia (NOA) and overall male infertility. DEGs across different datasets of transcriptomic profiling of testis from human patients with different causes of infertility/ impaired spermatogenesis and/or azoospermia were explored using the gene expression omnibus (GEO) database. Following the search using the GEOquery, 30 datasets were available, with 5 meeting the inclusion criteria. The DEGs for datasets were identified using limma R packages through the GEO2R tool. The annotated genes of the probes in each dataset were intersected with DEGs from all other datasets. Enriched Ontology Clustering for the identified genes was performed using Metascape to explore the possible connection or interaction between the genes. Twenty-five DEGs were shared between most of the datasets, which might indicate their role in the pathogenesis of male infertility. Of the 25 DEGs, eight genes (THEG, SPATA20, ROPN1L, GSTF1, TSSK1B, CABS1, ADAD1, RIMBP3) are either involved in the overall spermatogenic processes or at specific phases of spermatogenesis. We hypothesize that alteration in the expression of these genes leads to impaired spermatogenesis and, ultimately, male infertility. Thus, these genes can be used as potential biomarkers for the early detection of NOA.

Evaluation of seminal plasma HSPA2 protein as a biomarker of human spermatogenesis status

In mammals, testicular Heat shock-related 70 kDa protein 2 (HSPA2) is a chaperon strictly linked to spermatogenesis status, whereas its presence in spermatozoa ensures successful oocyte fertilization. However, there is little information on this protein in seminal plasma in infertile males. Based on our previous two independent studies, we have selected HSPA2 to evaluate this seminal plasma protein is a potential biomarker of correct spermatogenesis. Using immunoblotting and mass spectrometry (MS) we have screened human seminal plasma samples for the presence of HSPA2. Samples were obtained from individuals with normozoospermia, cryptozoospermia, non-obstructive and obstructive azoospermia. Our results showed a lack of HSPA2 in seminal plasma in all azoospermic males however, in cryptozoospermia the results were extremely diversified. Additionally, the application of 2-dimensional gel electrophoresis (2-DE) indicated the presence of additional protein isoforms suggesting possible mechanisms underlying the male infertility. Our findings suggest seminal plasma HSPA2 protein as a possible biomarker not only of spermatogenesis status, especially in cryptozoospermic males, but also as a biomarker predicting the success of reproductive treatment including assisted reproductive techniques (ART).

In Silico Analysis of Sperm From Ejaculates with Different Semen Characteristics

Background:

Male infertility is associated with altered characteristics of the sperm within the ejaculate. It is possible to find molecular explanations for the observed phenotypes and their consequences. This study aimed to analyze, using a specialized software, a gene set of transcriptomic data from different types of ejaculates.

Methods:

Data from ejaculate samples categorized as normal, oligospermia, and teratozoospermia were obtained from Gene Expression Omnibus (GEO). After normalization, the data average for each sample category was calculated and analyzed independently using Ingenuity Pathway Analysis (IPA).

Results:

Five important canonical pathways are involved in normal and altered semen samples (Oligospermia and teratozoospermia) except sirtuin signaling and mitochondrial dysfunction pathways. The five most important biological processes are identified in all semen phenotypes, but the only difference is the genes connected with initiation of RNA transcription in oligospermic and asthenospermic samples.

Conclusion:

Surprisingly, different types of ejaculates share many pathways and biological processes; sperm proteomics as a new global approach gives clues for the development of strategies to explain the reason for observed phenotypes of ejaculated spermatozoa, their possible effect on fertility, and for implementing research strategies in the context of infertility diagnosis and treatment.

FAM71F1 binds to RAB2A and RAB2B and is essential for acrosome formation and male fertility in mice

The acrosome is a cap-shaped, Golgi-derived membranous organelle that is located over the anterior of the sperm nucleus and highly conserved throughout evolution. Although morphological changes during acrosome biogenesis in spermatogenesis have been well described, the molecular mechanism underlying this process is still largely unknown. Family with sequence similarity 71, member F1 and F2 (FAM71F1 and FAM71F2) are testis-enriched proteins that contain a RAB2B-binding domain, a small GTPase involved in vesicle transport and membrane trafficking. Here, by generating mutant mice for each gene, we found that Fam71f1 is essential for male fertility. In Fam71f1-mutant mice, the acrosome was abnormally expanded at the round spermatid stage, likely because of enhanced vesicle trafficking. Mass spectrometry analysis after immunoprecipitation indicated that, in testes, FAM71F1 binds not only RAB2B, but also RAB2A. Further study suggested that FAM71F1 binds to the GTP-bound active form of RAB2A/B, but not the inactive form. These results indicate that a complex of FAM71F1 and active RAB2A/B suppresses excessive vesicle trafficking during acrosome formation.

Downregulation of KIF2C and TEKT2 is associated with male infertility and testicular carcinoma

Background: Genetic factors are important in spermatogenesis and fertility maintenance, and are potentially significant biomarkers for the early detection of infertility. However, further understanding of these biological processes is required.

Methods: In the present study, we sought to identify associated genes by reanalyzing separate studies from Gene Expression Omnibus datasets (GSE45885, GSE45887 and GSE9210) and validation datasets (GSE4797, 145467, 108886, 6872). The differential genes were used the limma package in R language. Gene ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses were performed by the clusterprofier package. The protein-protein interaction network was constructed by the STRING database. The interaction between mRNA and TF was predicted by miRWalk web. At last, The Cancer Genome Atlas data were used to identify hub gene expression levels in GEPIA web.

Results: The results showed that 27 shared genes associated with spermatogenesis. We effectively screen out two genes (KIF2C and TEKT2) and both validated by GSE4797, 145467, 108886 and 6872. Among 27 shared genes, KIF2C and TEKT2 both down-regulated in spermatogenesis. The network of TF-miRNA-target gene was established, we found KIF2C-miRNAs (has-miR-3154, 6075, 6760-5p, 1251-5p, 186-sp)-TFs (EP300, SP1) might work in spermatogenesis.

Conclusions: Our study might help to improve our understanding of the mechanisms in spermatogenesis and provide diagnostic biomarkers and therapeutics targets.

Identification and Characterization of lncRNA and mRNA in Testes of Landrace and Hezuo Boars

Simple Summary

Precocious puberty is an excellent reproductive trait in domestic animals, which can generate higher breeding benefits in livestock production. However, regulators associated with this sexual maturation process remain largely unknown. Chinese Hezuo (HZ) boars are known for their early sexual maturity. In this work, the characteristics of precocious puberty in HZ pigs were confirmed by histological analysis, and some important long noncoding RNA (lncRNA) and mRNA were identified in the testes of immature (30-day-old) and mature (120-day-old) HZ boars, which could play a key role in precocious puberty. These results will provide a theoretical basis for further research on the regulatory mechanism of precocious puberty, which is important for accelerating the breeding process of highly fertile animals.

Abstract

Chinese HZ boars are typical plateau miniature boars characterized by precocious puberty, which is closely related to testicular development and spermatogenesis. Accumulating evidence indicates that lncRNA is involved in the testicular development and regulation of spermatogenesis. However, little is known about the lncRNA precocious regulation in testicular development and spermatogenesis on early sexual maturity of HZ boars. Thus, we investigated the expression and characterization of lncRNA and mRNA in 30-day-old and 120-day-old HZ boar testes using transcriptome to explore precocious puberty. Landrace (LC) boar was treated as the control. Histological analyses indicated that HZ boar underwent puberty development at an earlier stage than LC boar and had achieved sexual maturity at 120 days old. RNA-Seq yielded a total of 187 lncRNAs and 984 mRNAs; these molecules were identified as possible candidates for precocious puberty. GO terms and KEGG pathways enrichment analyses revealed that the differentially expressed lncRNA and their targeted genes were involved in metabolic pathways regulating testis development and spermatogenesis, such as the PI3K-Akt, TGF-beta and Wnt pathways. Further screening, some lncRNA (such as LOC102166140, LOC110259451, and MSTRG.15011.2), and mRNA (such as PDCL2, HSD17B4, SHCBP1L, CYP21A2, and SPATA3) were found to be possibly associated with precocious puberty, which would add to our understanding of the molecular regulatory mechanisms of precocious puberty. This study provided valuable information for further study of the role of lncRNA and mRNA in the process of precocious puberty.

Evaluation of immune status in testis and macrophage polarization associated with testicular damage in patients with nonobstructive azoospermia

OBJECTIVE

Immune cells residing in the testicular interstitial space form the immunological microenvironment of the testis. They are assumed to play a role in maintaining testicular homeostasis and immune privilege. However, the immune status and related cell polarization in patients with nonobstructive azoospermia (NOA) remains poorly characterized. System evaluation of the testis immunological microenvironment in NOA patients may help to reveal the mechanisms of idiopathic azoospermia.

STUDY DESIGN

The gene expression patterns of immune cells in normal human testes were systematically analyzed by single-cell RNA sequencing (scRNA-seq) and preliminarily verification by the human protein atlas (HPA) online database. The immune cell infiltration profiles and immune status of patients with NOA was analyzed by single-sample gene set enrichment analysis (ssGSEA) and gene set variation analysis (GSVA) based on four independent public microarray datasets (GSE45885, GSE45887, GSE9210, and GSE145467), obtained from Gene Expression Omnibus (GEO) online database. The relationship between immune cells and spermatogenesis score was further analyzed by Spearman correlation analysis. Finally, immunohistochemistry (IHC) staining was performed to identify the main immune cell types and their polarization status in patients with NOA.

RESULTS

Both scRNA-seq and HPA analysis showed that testicular macrophages represent the largest pool of immune cells in the normal testis, and also exhibit an attenuated inflammatory response by expressing high levels of tolerance proteins (CD163, IL-10, TGF-β, and VEGF) and reduced expression of TLR signaling pathway-related genes. Correlation analysis revealed that the testicular immune score and macrophages including M1 and M2 macrophages were significantly negatively correlated with spermatogenesis score in patients with NOA (GSE45885 and GSE45887). In addition, the number of M1 and M2 macrophages was significantly higher in patients with NOA (GSE9210 and GSE145467) than in normal testis. GSVA analysis indicated that the immunological microenvironment in NOA tissues was manifested by activated immune system and pro-inflammatory status. IHC staining results showed that the number of M1 and M2 macrophages was significantly higher in NOA tissues than in normal testis and negatively correlated with the Johnson score.

CONCLUSION

Testicular macrophage polarization may play a vital role in NOA development and is a promising potential therapeutic target.

Whole-Exome Sequencing Analysis of Human Semen Quality in Russian Multiethnic Population

The global trend toward the reduction of human spermatogenic function observed in many countries, including Russia, raised the problem of extensive screening and monitoring of male fertility and elucidation of its genetic and ethnic mechanisms. Recently, whole-exome sequencing (WES) was developed as a powerful tool for genetic analysis of complex traits. We present here the first Russian WES study for identification of new genes associated with semen quality. The experimental 3 × 2 design of the WES study was based on the analysis of 157 samples including three ethnic groups—Slavs (59), Buryats (n = 49), and Yakuts (n = 49), and two different semen quality groups—pathozoospermia (n = 95) and normospermia (n = 62). Additionally, our WES study group was negative for complete AZF microdeletions of the Y-chromosome. The normospermia group included men with normal sperm parameters in accordance with the WHO-recommended reference limit. The pathozoospermia group included men with impaired semen quality, namely, with any combined parameters of sperm concentration <15 × 10⁶/ml, and/or progressive motility <32%, and/or normal morphology <4%. The WES was performed for all 157 samples. Subsequent calling and filtering of variants were carried out according to the GATK Best Practices recommendations. On the genotyping stage, the samples were combined into four cohorts: three sets corresponded to three ethnic groups, and the fourth set contained all the 157 whole-exome samples. Association of the obtained polymorphisms with semen quality parameters was investigated using the χ2 test. To prioritize the obtained variants associated with pathozoospermia, their effects were determined using Ensembl Variant Effect Predictor. Moreover, polymorphisms located in genes expressed in the testis were revealed based on the genomic annotation. As a result, the nine potential SNP markers rs6971091, rs557806, rs610308, rs556052, rs1289658, rs278981, rs1129172, rs12268007, and rs17228441 were selected for subsequent verification on our previously collected population sample (about 1,500 males). The selected variants located in seven genes FAM71F1, PPP1R15A, TRIM45, PRAME, RBM47, WDFY4, and FSIP2 that are expressed in the testis and play an important role in cell proliferation, meiosis, and apoptosis.

Gene and histomorphology alteration analysis in spermatogenesis arrest mouse model: a probable novel approach for infertility

Introduction

Approximately 15% of couples in the reproductive age are struggling with infertility which, in nearly half of them, is caused by male factors.

Material and methods

The present study comprised of two groups of sixteen C57BL/6 mice; each mouse received either an intraperitoneal injection of 30 mg/kg of an alkylating agent or the same amount of distilled water. Testes were harvested 30 days following the injection. Morphometric analysis of hematoxylin and eosin (H&E) stained slides including mean tubular area, diameter and intratubular particles were performed. Spermatogenesis rate was assessed by spermatogonial markers including promyelocytic leukemia zinc finger protein (PLZF) and neurogenin-3 (NGN3). Moreover, the expression rate of Wilms Tumor-1 (WT-1), A-Kinase Anchoring Protein 4 (AKAP4) and adenosine deaminase domain containing 1 (ADAD1) genes were evaluated via real-time polymerase chain reaction (RT-PCR).

Results

The body weight gradually increased in both groups after a period of 30 days, however, the increase was significantly (p-value = 0.023) lower in the chemically treated group. All the morphometric parameters were considerably decreased in the azoospermic mice. Also, promyelocytic leukemia zinc finger protein and neurogenin-3 expression dramatically declined (p-value <0.001 for both markers). In comparison with the negative control group, the expression rates of A-Kinase Anchoring Protein 4 and adenosine deaminase domain containing 1, two genes participating in the sperm structure, were remarkably reduced in the intervention group (p-value <0.001); however, our investigations demonstrated that the azoospermia model could induce a 5-fold upregulation in Wilms Tumor-1 gene expression.

Conclusions

Development of an azoospermia model can upregulate Wilms Tumor-1 gene expression in a higher rate after 30 days; however, expression of the testis-specific genes, A-Kinase Anchoring Protein 4 and adenosine deaminase domain containing 1, decreased after the intervention. To the best of our knowledge, this upregulation could be related to spermatogenesis recovery after the follow-up period.

Comparative Transcriptomic Analysis of Spermatozoa From High- and Low-Fertile Crossbred Bulls: Implications for Fertility Prediction

Crossbred bulls produced by crossing Bos taurus and Bos indicus suffer with high incidence of infertility/subfertility problems; however, the etiology remains poorly understood. The uncertain predictability and the inability of semen evaluation techniques to maintain constant correlation with fertility demand for alternate methods for bull fertility prediction. Therefore, in this study, the global differential gene expression between high- and low-fertile crossbred bull sperm was assessed using a high-throughput RNA sequencing technique with the aim to identify transcripts associated with crossbred bull fertility. Crossbred bull sperm contained transcripts for 13,563 genes, in which 2,093 were unique to high-fertile and 5,454 were unique to low-fertile bulls. After normalization of data, a total of 776 transcripts were detected, in which 84 and 168 transcripts were unique to high-fertile and low-fertile bulls, respectively. A total of 176 transcripts were upregulated (fold change > 1) and 209 were downregulated (<1) in low-fertile bulls. Gene ontology analysis identified that the sperm transcripts involved in the oxidative phosphorylation pathway and biological process such as multicellular organism development, spermatogenesis, and in utero embryonic development were downregulated in low-fertile crossbred bull sperm. Sperm transcripts upregulated and unique to low-fertile bulls were majorly involved in translation (biological process) and ribosomal pathway. With the use of RT-qPCR, selected sperm transcripts (n = 12) were validated in crossbred bulls (n = 12) with different fertility ratings and found that the transcriptional abundance of ZNF706, CRISP2, TNP2, and TNP1 genes was significantly (p < 0.05) lower in low-fertile bulls than high-fertile bulls and was positively (p < 0.05) correlated with conception rate. It is inferred that impaired oxidative phosphorylation could be the predominant reason for low fertility in crossbred bulls and that transcriptional abundance of ZNF706, CRISP2, TNP2, and TNP1 genes could serve as potential biomarkers for fertility in crossbred bulls.

The RNA content of human sperm reflects prior events in spermatogenesis and potential post-fertilization effects

Transcriptome analyses using high-throughput methodologies allow a deeper understanding of biological functions in different cell types/tissues. The present study provides an RNA-seq profiling of human sperm mRNAs and lncRNAs (messenger and long non-coding RNAs) in a well-characterized population of fertile individuals. Sperm RNA was extracted from twelve ejaculate samples under strict quality controls. Poly(A)-transcripts were sequenced and aligned to the human genome. mRNAs and lncRNAs were classified according to their mean expression values (FPKM: Fragments Per Kilobase of transcript per Million mapped reads) and integrity. Gene Ontology analysis of the Expressed and Highly Expressed mRNAs showed an involvement in diverse reproduction processes, while the Ubiquitously Expressed and Highly Stable mRNAs were mainly involved in spermatogenesis. Transcription factor enrichment analyses revealed that the Highly Expressed and Ubiquitously Expressed sperm mRNAs were primarily regulated by zinc-fingers and spermatogenesis-related proteins. Regarding the Expressed lncRNAs, only one-third of their potential targets corresponded to Expressed mRNAs and were enriched in cell-cycle regulation processes. The remaining two-thirds were absent in sperm and were enriched in embryogenesis-related processes. A significant amount of post-testicular sperm mRNAs and lncRNAs was also detected. Even though our study is solely directed to the poly-A fraction of sperm transcripts, results indicate that both sperm mRNAs and lncRNAs constitute a footprint of previous spermatogenesis events and are configured to affect the first stages of embryo development.

SPATS1 (spermatogenesis-associated, serine-rich 1) is not essential for spermatogenesis and fertility in mouse

SPATS1 (spermatogenesis-associated, serine-rich 1) is an evolutionarily conserved, testis-specific protein that is differentially expressed during rat male meiotic prophase. Some reports have suggested a link between SPATS1 underexpression/mutation and human pathologies such as male infertility and testicular cancer. Given the absence of functional studies, we generated a Spats1 loss-of-function mouse model using CRISPR/Cas9 technology. The phenotypic analysis showed no overt phenotype in Spats1^-/- mice, with both males and females being fertile. Flow cytometry and histological analyses did not show differences in the testicular content and histology between WT and knockout mice. Moreover, no significant differences in sperm concentration, motility, and morphology, were observed between WT and KO mice. These results were obtained both for young adults and for aged animals. Besides, although an involvement of SPATS1 in the Wnt signaling pathway has been suggested, we did not detect changes in the expression levels of typical Wnt pathway-target genes in mutant individuals. Thus, albeit Spats1 alteration might be a risk factor for male testicular health, we hereby show that this gene is not individually essential for male fertility and spermatogenesis in mouse.

Impact of platelet-rich plasma versus selenium in ameliorating induced toxicity in rat testis: histological, immunohistochemical, and molecular study

This study was conducted on forty adult rats divided into four groups: Group I (control) that is divided into subgroups A, B, and C and Group II (methotrexate (MTX)-treated); the rats were injected intraperitoneally with MTX at a dose of 1 mg/kg/week, for 8 weeks. Group III (MTX-Se co-treated) was injected with MTX like Group II plus an oral administration of selenium at a dose of 10 μg/kg b.w/day, for 8 weeks. Group IV (MTX-PRP co-treated), rats were injected intraperitoneally with MTX like Group II plus platelet-rich plasma (PRP) injection under the scrotum, three times with 2-week intervals (volume-0.1 ml per injection) and euthanized after 8 weeks. Histological, immunohistochemical, and genetic expression using qPCR and western blotting technique were conducted. There was improvement in histological structure of testes in most specimens of Group IV. The latter group revealed a significant decrease in Bax and an increase in Bcl-2. The regeneration of testicular tissue was more observed in Group IV as measured by an increase in mean number of PCNA. Moreover, Group IV revealed an increased genetic level of FSCN3, GCNF, UBQLN3, and DAZL. Both MTX-Se and MTX-PRP have an anti-inflammatory effect as measured by a reduction in NF-κb. The anti-oxidative effect of selenium and PRP was noticed by a decrease in the level of the iNos and an increase in eNos protein and the autophagy marker LC3. PRP has ameliorative effects on induced rat testicular toxicity as evaluated by morphological changes and confirmed by immunohistochemical reactions, genetic expression, and western blotting analyses including oxidative and anti- oxidative markers.

Weighted Correlation Gene Network Analysis Reveals New Potential Mechanisms and Biomarkers in Non-obstructive Azoospermia

Non-obstructive azoospermia (NOA) denotes a severe form of male infertility, whose etiology is still poorly understood. This is mainly due to limited knowledge on the molecular mechanisms that lead to spermatogenesis failure. In this study, we acquired microarray data from GEO DataSets and identified differentially expressed genes using the limma package in R. We identified 1,261 differentially expressed genes between non-obstructive and obstructive azoospermia. Analysis of their possible biological functions and related signaling pathways using the cluster profiler package revealed an enrichment of genes involved in germ cell development, cilium organization, and oocyte meiosis. Immune infiltration analysis indicated that macrophages were the most significant immune component of NOA, cooperating with mast cells and natural killer cells. The weighted gene coexpression network analysis algorithm generated three related functional modules, which correlated closely with clinical parameters derived from histopathological subtypes of NOA. The resulting data enabled the construction of a protein–protein interaction network of these three modules, with CDK1, CDC20, CCNB1, CCNB2, and MAD2L1 identified as hub genes. This study provides the basis for further investigation of the molecular mechanism underlying NOA, as well as indications about potential biomarkers and therapeutic targets of NOA. Finally, using tissues containing different tissue types for differential expression analysis can reflect the expression differences in different tissues to a certain extent. But this difference in expression is only related and not causal. The specific causality needs to be verified later.

MicroRNA-targeting in spermatogenesis: Over-expressions of microRNA-23a/b-3p and its affected targeting of the genes ODF2 and UBQLN3 in spermatozoa of patients with oligoasthenozoospermia

BACKGROUND

Male infertility is a multifactorial syndrome with diverse phenotypic representations. MicroRNAs (miRNAs) are small, non-coding RNAs that are involved in the post-transcriptional regulation of gene expression. Altered abundance levels of ODF2 and UBQLN3 have been reported in patients with different spermatogenic impairments. However, the transcriptional regulation of these two genes by miR-23a/b-3p is still unclear.

OBJECTIVES

To investigate experimentally whether miR-23a/b-3p targets the genes ODF2 and UBQLN3 and whether this targeting impacts abundance levels of ODF2 and UBQLN3 in patients with oligoasthenozoospermia.

MATERIALS AND METHODS

A total of 92 men attending a fertility clinic were included in the study, including 46 oligoasthenozoospermic men and 46 age-matched normozoospermic volunteers who served as controls. Reverse transcription-quantitative PCR (RT-qPCR), Western blot, and dual-luciferase (Firefly-Renilla) assays were used to validate the miRNAs and their target genes.

RESULTS

RT-qPCR revealed that miR-23a/b-3p was more abundant and ODF2 and UBQLN3 targets were less abundant in men with impaired spermatogenesis. Besides, Western blot shows that ODF2 and UBQLN3 protein levels were reduced in men with impaired spermatogenesis. In silico prediction and dual-luciferase assays revealed that potential links exist between the higher abundance level of miR-23a/b-3p and the lower abundance level of ODF2 and UBQLN3 targets. Mutations in the miR-23a/b-3p-binding site within the 3'UTRs (3'untranslated regions) of ODF2 and UBQLN3 genes resulted in abrogated responsiveness to miR-23a/b-3p. Correlation analysis showed that sperm count, motility, and morphology were negatively correlated with miR-23a/b-3p and positively correlated with the lower abundance level of UBQLN3, while ODF lower abundance level was positively correlated with sperm motility.

CONCLUSION

Findings indicate that the higher abundance level of miR-23a/b-3p and the lower abundance level of ODF2 and UBQLN3 targets are associated with oligoasthenozoospermia and male subfertility.

The Roles of MicroRNAs in Male Infertility

MicroRNAs applications were vastly studied throughout the years, spanning from potential cancer biomarkers to targeted therapies for various diseases. Out of these utilizations, this paper focuses on their role in male infertility. Approximately 10–15% of worldwide couples are affected by infertility. Out of these, 50% are due to male determinants. The majority of cases still have an undetermined cause. Previous studies have found that the aberrant expression of microRNAs could be linked to certain reproductive dysfunctions in males. Further on, this study looked into the most recent literature published on this subject in order to assess the connection between the up-/down-regulation of various microRNAs and the roles they play in male infertility. MicroRNAs were found to be abundant and stable in the seminal liquid, which led to a facile identification using regular RNA detection methods. It was observed that the concentration of microRNAs in semen was modified in the case of patients suffering from asthenozoospermia and azoospermia. Moreover, idiopathic male infertility was associated with a single nucleotide polymorphism of the microRNA binding site. Future studies should focus their attention on discovering future treatments against male infertility targeting specific microRNAs and also on developing new and improved contraceptive methods.

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

Non-obstructive azoospermia (NOA) is one of the most important causes of male infertility. It is mainly characterized by the absence of sperm in semen repeatedly or the number of sperm is small and not fully developed. At present, its pathogenesis remains largely unknown. The goal of this study is to identify hub genes that might affect biomarkers related to spermatogenesis. Using the clinically significant transcriptome and single-cell sequencing data sets on the Gene Expression Omnibus (GEO) database, we identified candidate hub genes related to spermatogenesis. Based on them, we performed Gene Ontology (GO) functional enrichment analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment pathway analyses, protein-protein interaction (PPI) network analysis, principal component analysis (PCA), cell cluster analysis, and pseudo-chronological analysis. We identified a total of 430 differentially expressed genes, of which three have not been reported related to spermatogenesis (C22orf23, TSACC, and TTC25), and the expression of these three hub genes was different in each type of sperm cells. The results of the pseudo-chronological analysis of the three hub genes indicated that TTC25 was in a low expression state during the whole process of sperm development, while the expression of C22orf23 had two fluctuations in the differentiating spermatogonia and late primary spermatocyte stages, and TSACC showed an upward trend from the spermatogonial stem cell stage to the spermatogenesis stage. Our research found that the three hub genes were different in the trajectory of sperm development, indicating that they might play important roles in different sperm cells. This result is of great significance for revealing the pathogenic mechanism of NOA and further research.

Infertility network and hub genes for nonobstructive azoospermia utilizing integrative analysis

Non-obstructive azoospermia (NOA) is the most severe form of male infertility owing to the absence of sperm during ejaculation as a result of failed spermatogenesis. The molecular mechanisms of NOA have not been well studied. Here, we revealed the dysregulated differentially expressed genes in NOA and related signaling pathways or biological processes. Cluster features of biological processes include spermatogenesis, fertilization, cilium movement, penetration of zona pellucida, sperm chromatin condensation, and being significantly enriched metabolic pathways in proximal tubule bicarbonate reclamation, aldosterone synthesis and secretion, glycolysis and glycogenesis pathways in NOA using Gene Ontology analysis and pathway enrichment analysis. The NOA gene co-expression network was constructed by weighted gene co-expression network analysis to identify the hub genes (CHD5 and SPTBN2). In addition, we used another Gene Expression Omnibus dataset (GSE45887) to validate these hub genes. Furthermore, we used the Seurat package to classify testicular tissue cells from NOA patients and to characterize the differential expression of hub genes in different cell types from different adult males based on the scRNA-seq dataset (GSE106487). These results provide new insights into the pathogenesis of NOA. Of particular note, CHD5 and SPTBN2 may be potential biomarkers for the diagnosis and treatment of NOA.

Deletion of the Spata3 Gene Induces Sperm Alterations and In Vitro Hypofertility in Mice

Thanks to the analysis of an Interspecific Recombinant Congenic Strain (IRCS), we previously defined the Mafq1 quantitative trait locus as an interval on mouse Chromosome 1 associated with male hypofertility and ultrastructural abnormalities. We identified the Spermatogenesis associated protein 3 gene (Spata3 or Tsarg1) as a pertinent candidate within the Mafq1 locus and performed the CRISPR-Cas9 mediated complete deletion of the gene to investigate its function. Male mice deleted for Spata3 were normally fertile in vivo but exhibited a drastic reduction of efficiency in in vitro fertilization assays. Mobility parameters were normal but ultrastructural analyses revealed acrosome defects and an overabundance of lipids droplets in cytoplasmic remnants. The deletion of the Spata3 gene reproduces therefore partially the phenotype of the hypofertile IRCS strain.

Decreased expression of AKAP4 and TyrPho proteins in testis, epididymis, and spermatozoa with low sexual performance of mice induced by modified CUMS

The molecular mechanism of chronic stress especially reduced motility, a major cause of male infertility, has not been proved. It is known that A-kinase anchor protein 4 (AKAP4) and tyrosine-phosphorylated (TyrPho) proteins are involved in progressive motility. This study aimed to investigate the effect of chronic unpredictable mild stress (CUMS) on sexual behaviours, sperm quality, and expressions of AKAP4 and TyrPho proteins in testis, epididymis, and spermatozoa. Sixteen male mice were divided into control and CUMS groups (n = 8/group). Animals were induced by a stressor from twelve stressors for 36 days. Sexual behaviours, corticosterone and testosterone, sperm parameters, and histopathology were observed. The expressions of AKAP4 and TyrPho proteins in testis, epididymis, and spermatozoa were examined. Results showed that CUMS significantly increased corticosterone while serum testosterone level was decreased. Sexual behaviours and sperm parameter quality were significantly decreased. CUMS mice showed vacuolisation and pyknotic cells in seminiferous epithelium and less sperm mass was observed within epididymal lumen. CUMS decreased expressions of AKAP4 and TyrPho proteins in testis, epididymis, and spermatozoa. In conclusion, the decreased expression of AKAP4 and TyrPho proteins may be a mechanism associated with low semen qualities particularly decrease of sperm motility in CUMS.

Testicular STAC3 regulates Leydig cell steroidogenesis through potentiating mitochondrial membrane potential and StAR processing

SH3 and cysteine-rich protein 3 (STAC3), a small adapter protein originally identified as a core component of excitation–contraction coupling machinery, regulates the voltage-induced Ca²⁺ release in skeletal muscle. However, the possibility of additional, as yet unknown, non-muscle effects of STAC3 cannot be ruled out. Herein, we provide the evidence for the expression and functional involvement of STAC3 in spermatogenesis. STAC3 expression was localized in the testicular interstitium of rodent and human testes. By using the cytotoxic drug ethylene dimethane sulfonate (EDS), STAC3 expression was observed to be decreased sharply in rat testis after selective withdrawal of Leydig cells (LCs), and reappeared immediately after LCs repopulation, indicating that testicular expression of STAC3 mainly stems from LCs. From a functional standpoint, in vivo lentiviral vector–mediated suppression of STAC3 resulted in a significant decrease in testosterone production, and thereafter caused impairment of male fertility by inducing oligozoospermia and asthenospermia. The indispensible involvement of STAC3 in testicular steroidogenesis was validated using the in vivo knockdown model with isolated primary LCs as well as in vitro experiments with primary LCs. By generating the TM3^Stac3−/− cells, we further revealed that STAC3 depletion attenuated mitochondrial membrane potential and StAR processing in db-cAMP-stimulated LCs. Thus, the inhibitory effect of STAC3 deficiency on testicular steroidogenesis may be ascribed to a disturbed mitochondrial homeostasis. Collectively, the present results strongly suggest that STAC3 may function as a novel regulator linking mitochondrial homeostasis and testicular steroidogenesis in LCs. Our data underscore an unexpected reproductive facet of this muscle-derived factor.

Downregulation of gene expression and the outcome of ICSI in severe oligozoospermic patients: A preliminary study

Preimplantation embryo development might be influenced by a specific set of transcripts that are delivered to the oocyte by the sperm. The aim of the study was to determine the relationship between the level of selected transcripts in spermatozoa and preimplantation development of the embryos in couples with severe oligozoospermia undergoing intracytoplasmic sperm injection (ICSI) procedure. Therefore, we assessed messenger RNA (mRNA) levels of genes involved in fertilization events, oocyte activation, chromatin remodeling, and DNA repair in severe oligozoospermic compared with normozoospermic men as well as morphokinetic parameters of embryos using the time-lapse imaging system. mRNA profiling (44 genes), in mature sperm, was carried out with custom-designed 384-well TLDA Cards. The morphokinetic parameters of zygotes and embryos were recorded by using a time-lapse imaging system. The transcript levels of 21 genes were significantly decreased in the severe oligozoospermic group. Most were associated with fertilization events, oocyte activation and embryonic genome activation. Among them, mRNA of AKAP4 and PTK7 was greatly reduced, moreover, the transcripts of PLCζ and POU5F1, essential for OA and EGA, were not detected at all in patients with severe oligozoospermia. Moreover, the reduced expression of genes important for spermatogenesis, chromatin remodeling and DNA repair was also observed in this group. Time-lapse analysis revealed that fertilization failure occurred in 14% of retrieved oocytes and 90% of all degenerated embryos did not reach morula stage. This study provides preliminary results indicating a significant decrease in transcripts of genes important for spermatogenesis and early preimplantation development in the mature sperm of men with severe oligozoospermia.

Altered gene expression profiles of testicular tissues from azoospermic patients with maturation arrest

Maturation arrest is a common cause of male infertility which has caused worldwide concern, and its pathophysiological process remains further elucidation. Our study aimed to identify genetic characteristics of maturation arrest by comparing gene expression between maturation arrest and normal samples using microarray technology. A total of 6,373 genes were identified differentially expressed (p < .05, fold change > 2.0 or <-2.0) and 1,594 genes were selected as statistically significant after Bonferroni correction, including 419 up-regulated and 1,175 down-regulated genes. Microarray data were validated by quantitative reverse transcriptase-polymerase chain reaction. Bioinformation analysis was performed to explore genetic function of statistically significant genes. Gene Ontology results showed the statistically significant genes enriched in sexual reproduction, spermatogenesis and male gamete generation. Reactome pathway analysis highlighted the olfactory signalling pathway, fertilisation, developmental biology, etc. One module and eight hub genes were found to be involved in ubiquitin-mediated proteolysis and may affect as indicators of spermatogenic process through protein-protein interaction analysis. Our study provided a comprehensive genetic characteristic of differential expressed genes in testicular tissues from maturation arrest patients and speculated several genes as potential indicators of disease.

Levels of endocrine-disrupting chemicals are associated with changes in the peri-pubertal epigenome

Puberty marks a transition period, which leads to the attainment of adult sexual maturity. Timing of puberty is a strongly heritable trait. However, large genetic association studies can only explain a fraction of the observed variability and striking secular trends suggest that lifestyle and/or environmental factors are important. Using liquid-chromatography tandem-mass-spectrometry, we measured endocrine-disrupting chemicals (EDCs; triclosan, bisphenol A, benzophenone-3, 2,4-dichlorophenol, 11 metabolites from 5 phthalates) in longitudinal urine samples obtained biannually from peri-pubertal children included in the COPENHAGEN puberty cohort. EDC levels were associated with blood DNA methylation profiles from 31 boys and 20 girls measured both pre- and post-pubertally. We found little evidence of single methylation sites that on their own showed association with urinary excretion levels of EDCs obtained either the same-day or measured as the yearly mean of dichotomized EDC levels. In contrast, methylation of several promoter regions was found to be associated with two or more EDCs, overlap with known gene-chemical interactions, and form a core network with genes known to be important for puberty. Furthermore, children with the highest yearly mean of dichotomized urinary phthalate metabolite levels were associated with higher promoter methylation of the thyroid hormone receptor interactor 6 gene (TRIP6), which again was mirrored by lower circulating TRIP6 protein levels. In general, the mean TRIP6 promoter methylation was mirrored by circulating TRIP6 protein levels. Our results provide a potential molecular mode of action of how exposure to environmental chemicals may modify pubertal development.

RNA profiling of human testicular cells identifies syntenic lncRNAs associated with spermatogenesis

STUDY QUESTION

Is the noncoding transcriptional landscape during spermatogenesis conserved between human and rodents?

SUMMARY ANSWER

We identified a core group of 113 long noncoding RNAs (lncRNAs) and 20 novel genes dynamically and syntenically transcribed during spermatogenesis.

WHAT IS KNOWN ALREADY

Spermatogenesis is a complex differentiation process driven by a tightly regulated and highly specific gene expression program. Recently, several studies in various species have established that a large proportion of known lncRNAs are preferentially expressed during meiosis and spermiogenesis in a testis-specific manner.

STUDY DESIGN, SIZE, DURATION

To further investigate lncRNA expression in human spermatogenesis, we carried out a cross-species RNA profiling study using isolated testicular cells.

PARTICIPANTS/MATERIALS, SETTING, METHODS

Human testes were obtained from post-mortem donors (N = 8, 51 years old on average) or from prostate cancer patients with no hormonal treatment (N = 9, 80 years old on average) and only patients with full spermatogenesis were used to prepare enriched populations of spermatocytes, spermatids, Leydig cells, peritubular cells and Sertoli cells. To minimize potential biases linked to inter-patient variations, RNAs from two or three donors were pooled prior to RNA-sequencing (paired-end, strand-specific). Resulting reads were mapped to the human genome, allowing for assembly and quantification of corresponding transcripts.

MAIN RESULTS AND THE ROLE OF CHANCE

Our RNA-sequencing analysis of pools of isolated human testicular cells enabled us to reconstruct over 25 000 transcripts. Among them we identified thousands of lncRNAs, as well as many previously unidentified genes (novel unannotated transcripts) that share many properties of lncRNAs. Of note is that although noncoding genes showed much lower synteny than protein-coding ones, a significant fraction of syntenic lncRNAs displayed conserved expression during spermatogenesis.

LARGE SCALE DATA

Raw data files (fastq) and a searchable table (.xlss) containing information on genomic features and expression data for all refined transcripts have been submitted to the NCBI Gene Expression Omnibus under accession number GSE74896.

LIMITATIONS, REASONS FOR CAUTION

Isolation procedures may alter the physiological state of testicular cells, especially for somatic cells, leading to substantial changes at the transcriptome level. We therefore cross-validated our findings with three previously published transcriptomic analyses of human spermatogenesis. Despite the use of stringent filtration criteria, i.e. expression cut-off of at least three fragments per kilobase of exon model per million reads mapped, fold-change of at least three and false discovery rate adjusted P-values of less than <1%, the possibility of assembly artifacts and false-positive transcripts cannot be fully ruled out.

WIDER IMPLICATIONS OF THE FINDINGS

For the first time, this study has led to the identification of a large number of conserved germline-associated lncRNAs that are potentially important for spermatogenesis and sexual reproduction. In addition to further substantiating the basis of the human testicular physiology, our study provides new candidate genes for male infertility of genetic origin. This is likely to be relevant for identifying interesting diagnostic and prognostic biomarkers and also potential novel therapeutic targets for male contraception.

STUDY FUNDING/COMPETING INTEREST(S)

This work was supported by l'Institut national de la santé et de la recherche médicale (Inserm); l'Université de Rennes 1; l'Ecole des hautes études en santé publique (EHESP); INERIS-STORM to B.J. [N 10028NN]; Rennes Métropole 'Défis scientifiques émergents' to F.C (2011) and A.D.R (2013). The authors have no competing financial interests.

Identification of new protein biomarkers associated with the boar fertility using iTRAQ-based quantitative proteomic analysis

In this study, we performed the isobaric tags for relative and absolute quantitation (iTRAQ) proteomic analysis in the spermatozoa of Landrace boars with different fertility potentials and investigated the ability of sperm acrosome associated 4 (SPACA4) and IZUMO family member 2 (IZUMO2) to predict the reproductive perform of boars. The iTRAQ results revealed that 202 proteins were up-regulated and 43 proteins were down-regulated in the spermatozoa from high fertility boars. SPACA4 and IZUMO2 protein levels were significantly up-regulated in the spermatozoa from high fertility boars. SPACA4 and IZUMO2 expression were specifically detected in the adult boar testis. SPACA4 levels were positively correlated with Sow's farrowing rate and reproductive efficiency, but not litter size. IZUMO2 were positively correlated with litter size, Sow's farrowing rate and reproductive efficiency. Treating the boar semen with SPACA4 or IZUMO2 antibodies for 30 min and 60 min failed to affect the sperm motility; while treating the semen with SPACA4 antibody significantly reduced the fertilization and cleavage rates. Similar results for fertilization and cleavage rates were found in IZUMO2 antibody-treated semen. Collectively, our results indicated that protein levels of SPACA4 and IZUMO2 in the spermatozoa were positively related to the reproductive performance of Landrace boars.