Front. Mol. Biosci.Frontiers in Molecular BiosciencesFront. Mol. Biosci.2296-889XFrontiers Media S.A.122682910.3389/fmolb.2023.1226829Molecular BiosciencesOriginal ResearchTranscriptomic signatures for human male infertilityHodžić et al.10.3389/fmolb.2023.1226829HodžićAlenka1MaverAleš1ZornBranko2PetrovičDaniel3KunejTanja4PeterlinBorut1*1Clinical Institute of Genomic Medicine, University Medical Centre Ljubljana, Ljubljana, Slovenia2Andrology Unit, Reproductive Unit, Department of Obstetrics and Gynecology, University Medical Centre Ljubljana, Ljubljana, Slovenia3Institute of Histology and Embryology, Faculty of Medicine, University of Ljubljana, Ljubljana, Slovenia4Biotechnical Faculty, University of Ljubljana, Ljubljana, Slovenia
Edited by:Zhiyun Wei, Shanghai First Maternity and Infant Hospital, China
Reviewed by:Shubo Jin, Chinese Academy of Fishery Sciences (CAFS), China
Hemakumar M. Reddy, Brown University, United States
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
Introduction: Male infertility is a common, complex disorder. A better understanding of pathogenesis and etiology is needed for timely diagnosis and treatment. The aim of this study, therefore, was to identify genes involved in the pathogenesis of idiopathic male infertility based on data from transcriptomic level supported with data from genomic level.
Materials and methods: First, we performed whole gene expression analysis in 20 testis biopsy samples of patients with severely impaired (10) and normal spermatogenesis (10). Further, we have performed systematic review of comparable male infertility studies and overlapped the most significantly expressed genes identified in our study with the most differentially expressed genes from selected studies. Gene Ontology analysis and KEGG functional enrichment have been performed with Enrichr analysis tool. Additionally, we have overlapped these genes with the genes where rare variants have been identified previously.
Results: In 10 patients with severely impaired spermatogenesis and 10 controls, we identified more than 1,800 differentially expressed genes (p < 0.001). With the systematic review of three previously performed microarray studies that have met inclusion criteria we identified 257 overlapped differentialy expressed genes (144 downregulated and 113 upregulated). Intersection of genes from transcriptomic studies with genes with identified rare variants revealed a total of 7 genes linked with male infertility phenotype (CYP11A1, CYP17A1, RSPH3, TSGA10, AKAP4, CCIN, NDNF).
Conclusion: Our comprehensive study highlighted the role of four genes in pathogenesis of male infertility and provided supporting evidence for three promising candidate genes which dysfunction may result in a male infertility disorder.
idiopathic male infertilitygene expresiontranscriptometestisspermatogenesissection-at-acceptanceMolecular Diagnostics and Therapeutics1 Introduction
Male infertility is a worldwide health problem affecting about 7% of adult human males (Krausz and Riera-Escamilla, 2018). Despite the considerable research effort, in a large proportion of infertile males the cause of their infertility is unknown (idiopathic), whereas genetic factors have been considered as a major contributing cause for them (Hwang et al., 2010).
Studies of animal models have identified hundreds of genes that influence reproduction and have uncovered various pathogenic pathways that could affect infertility in humans. Still, only a few mutations in human genes have been associated with male infertility (Jamsai and O’Bryan, 2010; Massart et al., 2011). One probable explanation is that for many years genetic research on male infertility involved the targeted sequencing of individual genes in cohorts of infertile men and fertile controls with a limited number of whole-exome sequencing studies performed in the last few years (Fakhro et al., 2018; Hodžić et al., 2020).
Human spermatogenesis is a complex process; therefore, these targeted approaches proved ineffective in identifying genetic variants responsible for infertility. Global gene-expression profiling has been used extensively in studying complex diseases and have already demonstrated potential for identifying genomic features associated with various diseases (Cookson et al., 2009). Contrary to rodent models and other pathologies in humans, attempts to perform similar investigation in the field of male infertility have been limited due to various restrictions (e.g., restrictions to human testis tissue).
A few microarray studies have been performed previously and have provided valuable insights into the complexity of spermatogenesis by identifying numerous testis expressed and testis-specific genes in humans (Feig et al., 2006; Ellis et al., 2007; Spiess et al., 2007; Okada et al., 2008). However, individual studies are usually not reproducible, and comparing studies with different designs may conceal the specific mechanisms underlying infertility. Kui et al. (2019) integrated results from different microarray data sets and focused on differentially expressed genes enriched in four divided pathological phenotypes (Kui et al., 2019). Furthermore, we hypothesized that integrating the data extracted from the transcriptomic level with the data originating from the genomic level could led us to identify potential biomarkers for male infertility. Namely, like all “omic” approaches, findings of rare, potentially pathogenic variants in infertile men need further evidence to distinguish the real signal from noise.
For this purpose, we performed global gene expression profiling on human testis samples in patients with severely impaired and normal spermatogenesis. In addition to our study, we systematically reviewed overlapping genes from comparable previous studies. We identified genes demonstrating differential expression and rare potentially pathogenic variants as a final step.
2 Materials and methods2.1 Ethics statement
The study was approved by National medical ethics committee (reference number: 73/05/12). All patients gave informed written consent to participate in the study.
2.2 Testicular biopsy samples
All biopsy samples were obtained from patients attending the outpatient infertility clinic of the Andrology Centre, Department of Obstetrics and Gynecology in Ljubljana. The study population comprised ten patients with non-obstructive azoospermia and ten patients with obstructive azoospermia (controls).
Testicular biopsies were performed under local anaesthesia. Following unilateral hemiscrototomy, a small testicular incision was made and at least two samples of testicular tissue were taken from each testis. The first sample was fixed in Bouin’s solution, routinely embedded in paraffin, and cut at a section thickness of 5 μm. The sections were stained with HE, PAS and van Gieson–Weigert. A systematic histological evaluation was performed under light microscopy. More than 100 seminiferous tubules were scored for each patient. The results were expressed as a relative number of tubules showing Sertoli cells, spermatogonia, spermatocytes, round and elongated spermatids, and spermatozoa. All examinations were made by the same observer (J.S.). The diagnoses were as follows: 10 patients with normal spermatogenesis, 2 with early and late maturation arrest and 8 with Sertoli-cell-only syndrome (SCOS). Before sperm recovery, medical history was established, testicular volume measured and serum FSH level was assessed. The clinical characteristics of patients are presented in Table 1.
Clinical characteristics of patients included in the transcriptome profiling study.
Patient
Clinical findingsa
Testicular volume (mL)
FSH level (IU/l)
Sperm counta
Sperm motilityb
Histologyb
Johnsen’s score
1
bilateral epididymal enlargement
12/13
3,7
3
2
normal spermatogenesis
10
2
bilateral epididymal enlargement
25/28
4,75
3
2
normal spermatogenesis
9
3
bilateral epididymal enlargement
13/16
3,2
3
2
normal spermatogenesis
9
4
bilateral epididymal enlargement
22/25
7,02
3
2
normal spermatogenesis
9
5
CABVD
14/14
3,63
3
2
normal spermatogenesis
9
6
CABVD
28/28
1,68
2
0
normal spermatogenesis
9
7
enlarged right epididymis
15/15
8,3
3
2
Hypospermatogenesis, mononuclear infiltration
9
8
enlarged right epididymis 8 years after
15/15
2,2
3
2
hypospermatogenesis
9
9
no enlargement
25/25
1,79
2
1
hypospermatogenesis
9
10
enlarged right epididymis
30/30
1,27
2
2
hypospermatogenesis, EMA
9
11
NA
NA
NA
0
0
LMA
6
12
no
15/15
12,5
0
0
SCOS, EMA, hypospermatogenesis
5
13
no
20/20
14,2
0
0
SCOS, EMA
5
14
no
12/15
12,8
0
0
SCOS, EMA
5
15
no
8/10
13
0
0
SCOS, EMA
5
16
hypogonadism
1/1
82
0
0
SCOS, EMA
2
17
hypogonadism
6/4
59,8
0
0
SCOS, LCH
2
18
no
15/15
44
0
0
SCOS, EMA
2
19
8/10
28,7
0
0
SCOS, EMA
1
20
hypogonadism
2/2
36,4
0
0
SCOS, LCH
1
CABVD, congenital bilateral absence of vasa deferentia.
EMA, early maturation arrest; LMA, late maturation arrest; LCH, Leydig cell hyperplasia; NA, not available.
Total RNA was isolated from frozen testis biopsy samples by using the TissueLyserLT (Qiagen) apparatus and RNeasy Plus Micro kit (Qiagen) according to the manufacturer’s instructions. About 3.5–4.5 mg weight of tissue was disrupted and homogenized in 2 mL centrifuge tubes containing 5-mm stainless steel beads and lysis buffer (Buffer RLT plus) at a frequency of 50 Hz for 10 min in a TissueLyser device. Homogenate was subsequently centrifuged for 3 min, collected and eluted in a gDNA eliminator spin column. Ethanol was added to the flow-through to provide appropriate binding conditions for RNA and samples were then applied to an RNeasy MinElute spin columns, where total RNA was bound to the membrane and contaminants were efficiently whashed away. RNA was then eluted in 14 µL of water. RNA concentrations were determined using NanoDrop 2000C Spectrophotometer (Thermo Scientific) and RNA quality was verified by Agilent Bioanalyzer 2,100 (Agilent Technologies).
2.4 Microarray experiment and data analysis
The measure of expression was performed on isolated RNA from biopsy samples of patients with obstructive and non-obstructive azoospermia using Agilent Whole Human Genome 4 × 44 microarrays (Agilent design id: 14850). Gene expression platform contains 264 × 10 biological probes, targeting altogether 19.596 unique mRNA sequences according to NCBI Reference Human Genome Build version 33. Sample preparing, labeling and amplification RNA (Low Input Quick Amp Labeling Kit, two color, Agilent Technologies), hybridization, washing and scanning were performed according to the manufacturer’s recommendations (Agilent Technologies). Differential gene expression was measured in relation to the Agilent’s Universal Human Reference RNA and each experimental sample was hybridized against this common reference sample.
After hybridization, microarray slides were scanned using Agilent High Resolution Microarray Scanner System, using the manufacturer’s recommended scanning settings. Subsequently, microarray features were extracted using Agilent Feature Extraction software v10.7.3.1. Post-processing steps included intra-array lowess and inter-array quantile normalization to correct for potential bias resulting from differential stability of cyanine dyes. Fluorescent values were offset by 100 units to reduce the anomalous dispersion of fold change (FC) values at lower signal intensities. MA and multidimensional scaling (MDS) plots were inspected for each array to detect any systemic error resulting from the preceding steps.
The raw data are deposited in NCBI’s Gene Expression Omnibus (Edgar et al., 2002) and are accessible through GEO Series accession number GSE145467 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE145467).
Results were statistically analyzed using a linear model fit in limma package for Bioconductor in an R statistical environment. To account for multiple testing, obtained significance values were corrected using the Benjamini-Hochberg method, and the adjusted significance threshold was set at <0.05.
2.5 Systematic review
To identify overlapping genes from studies that could be included in the systematic review, a PubMed search (www.ncbi.nlm.nih.gov/pubmed/) until June 2022 was performed using the search strings: ‘male infertility’ AND ‘human spermatogenesis’ AND ‘gene expression’ AND ‘testis’. The systematic review was conducted according to the Preferred Reporting Items for Systematic Reviews and Meta-analyses (PRISMA) statement (Moher et al., 2009). The following criteria were used: (a) comparable disease states as in our experiment, (b) the same tissue samples, (c) common reference design approach towards measuring differential expression in the microarray experiment and (d) availability of raw datasets at either GEO or ArrayExpress repositories.
All the steps described in the following sections were performed in using R statistical language version 2.9.2, using Bioconductor environment.
Due to the differences in methodology and statistical analyses among the 3 studies, FDR values were not used, and instead, differentially expressed transcripts with p < 0.001, were used for the comparisons. Shared differentially expressed genes were visualized using an online VENN tool, available from the VIB/UGent Bioinformatics & Systems Biology (http://bioinformatics.psb.ugent.be/cgi-bin/liste/Venn/calculate_venn.htpl).
To gain more insight into key processes that may possibly explain functional differences among the testicular samples from severely impaired and normal spermatogenesis types, we carried out functional annotation analysis of the top overlapping differentially expressed genes, attaining adjusted p-values below 0.001 and using hypergeometric test, based on genes’ relation to GeneOntology (GO Biological Process 2021; https://maayanlab.cloud/Enrichr/) and KEGG terms (KEGG 2021 Human, https://maayanlab.cloud/Enrichr/).
2.6 Common genes identified in both transcriptomic and sequencing studies searching for monogenic etiology of male infertility
To provide additional support for identified genes from transcriptomic studies, we searched for overlapped differentially expressed genes with evidence to be putative candidates for the monogenic aetiology of male infertility. For this purpose, we have used the selected genes from a systematic review of next-generation sequencing studies, identifying rare, potentially pathogenic variants published recently (Houston et al., 2021). We have included genes moderately, strongly, or definitively linked to male infertility phenotypes and genes with limited evidence. We have excluded genes scored as “no evidence”.
3 Results3.1 Global gene-expression profiling
Using Agilent Whole Human Genome 4 × 44 microarrays (Agilent design id: 14850), we analyzed global gene expression changes in testis biopsy samples from 10 subjects with severely impaired spermatogenesis and 10 with normal spermatogenesis. Excluding one sample identified as an outlier left 19 samples for analysis. With a strict statistical criterion (adj. p < 0.001), we identified more than 1801 differentially expressed genes. Of these 1801, 781 were downregulated and 1,020 were up-regulated.
3.2 Systematic review and gene set enrichment analyses
Among the identified 628 articles, transcriptomic profiling of the human testis was identified in 10 articles. Of these 10 articles, only three previously performed transcriptomic studies have met inclusion criteria. Therefore GEO datasets with accession numbers GSE9210, GSE4797 and GSE45885, in addition to our own microarray results, were finally included in the analysis. Details on the three studies included in the systematic review are given in Table 2.
Description of datasets included in the meta-analysis.
Dataset accession
Number of samples—azoospermia
Number of samples—obstructive forms of infertility
Array platform
Number of probes on the array
GSE145467 (our study)
10
10
Agilent-014850 Whole Human Genome Microarray 4 × 44K G4112F
∼19.596 probes
GSE9210
47
11
Agilent-012097 Human 1A Microarray (V2) G4110B
∼22000 probes
GSE4797
5
18
GE Healthcare/Amersham Biosciences CodeLink™ UniSet Human 20K I Bioarray
∼23000 probes
GSE45885
27
4
Affymetrix Human Gene 1.0 ST Array
∼28132 probes
Aiming to identify the specific mechanisms underlying infertility, we focused on overlapping genes among four analyzed studies. With a strict statistical criterion (p < 0.001), we successfully identified a signature comprised of 257 shared genes (144 down-regulated and 113 up-regulated) with changed expression in infertile men (Figures 1A, B). However, one study did not successfully pass through this sorting criterion. Considering good overlap between the other three studies, the possible reason for this deviation could be in the small set of control samples.
Venn diagram of the most significant downregulated (A) and upregulated genes (B) genes among all studies (p value less than 0.0001).
Top biological processes according to GO enriched with common genes across all studies were genes related to collagen fibril organization and spindle and mitotic assembly checkpoint signaling, while according to KEGG pathways, common genes across all studies were genes related to ovarian steroidogenesis and p53 signaling pathway (Bar Chart 1 and Bar Chart 2).
Gene set enrichment analysis of common down-regulated and up-regulated genes across all four studies using GENE ONTOLOGY (biological process) [(A) up-regulated, (B) down-regulated].
Gene set enrichment analysis of common down-regulated and up-regulated genes across all four studies using KEGG pathway functional annotation [(A) up-regulated, (B) down-regulated].
3.3 Common genes identified in transcriptomic and next-generation sequencing studies
The intersection of genes from transcriptomic studies and those reported as candidates for monogenic aetiology revealed seven genes linked with male infertility phenotype. CYP11A1, CYP17A1, RSPH3, and TSGA10 genes were classified as confidently associated with the phenotype (moderate, strong, or definitive), whereas AKAP4, CCIN, NDNF present genes classified as limited evidence.
4 Disscusion
In the present study, we performed global gene expression profiling on human testis samples in patients with severely impaired spermatogenesis and systematically reviewed previous transcriptomic studies to identify genes consistently differently expressed across the studies.
To provide additional support, we searched for differentially expressed genes with evidence to be putative candidates for the monogenic aetiology of male infertility. The integratomic approach revealed four genes classified as confidently causally associated with phenotype, whereas another three present candidate genes whose dysfunction may result in a male infertility disorder.
CYP11A1, CYP17A1 and RSPH3 were upregulated in patients with normal spermatogenesis when compared to the patients with impaired spermatogenesis.
Considering that CYP17A1 and CYP11A1 are required for proper synthesis of androgens, noticed upregulation in patients with impaired spermatogenesis could result from organism efforts to maintain stability. Biallelic variants in CYP11A1 and CYP17A1 gene can impair sexual differentiation caused by a complete or partial loss of steroid hormone production. The broad clinical spectrum has been reported in affected 46XY males and may vary from normal genitalia and surgically repairable defects, including cryptorchidism and hypospadias, to complete feminization of external gonads, accompanied by adrenal dysfunction (Sherbet et al., 2003; Rubtsov et al., 2009; Kok et al., 2010; Parajes et al., 2012; Lara-Velazquez et al., 2017; Kolli et al., 2019; Kallali et al., 2020). The broad phenotypic spectrum has been reported for numerous genes, and in this manner, establishing a specific genotype-phenotype relationship is not necessarily straightforward. Variability in genotype-phenotype correlation have been previously recognized in reproductive system syndromes with endocrine disorders. One such example is the NR5A1 gene, where pathogenic variants are associated with a phenotypic spectrum ranging from disorders of sex development to oligo/azoospermia (Domenice et al., 2016).
RSPH3 biallelic variants represent a recurrent genetic cause of Primary ciliary dyskinesia (PCD) in human (Jeanson et al., 2015). Although the RSPH3 gene is abundantly expressed in testicular tissue, the phenotype of PCD-associated male infertility caused by the defect in this ciliopathy-related gene has been rarely described (Wu et al., 2020). Mutations in RSPH3 have been reported in a sterile PCD patient with severe asthenoteratospermia characterized by multiple flagellar malformations (Wu et al., 2020). However, the mutation in this gene in association with non-obstructive azoospermia has not been identified so far.
Based on our analyses, the expression levels of the TSGA10 gene were downregulated across studies. It has been suggested that decreased expression of TSGA10 is associated with the reduction of autophagy process and increased ROS levels, which lead to aberrant spermatid differentiation and maturation (Asgari et al., 2021). The TSGA10 gene is involved in active cell division, differentiation and cell migration and is highly conserved among different species (Behnam et al., 2009). Functional analysis of the Tsga10 gene in knockout mice demonstrated that TSGA10 contribute to the correct arrangement of a mitochondrial sheath in spermatozoa, and its deficiency leads to male infertility (Luo et al., 2020). Previous studies have reported a homozygous mutation in this gene in patients with acephalic spermatozoa (Sha et al., 2018; Ye et al., 2020).
According to the literature, there has been limited evidence for the association of AKAP4, CCIN and NDNF genes with male infertility. Two recently published studies identified a homozygous missense mutation and a compound heterozygous mutation of the CCIN gene in patients with teratozoospermia, while a missense variant in the AKAP4 gene was associated with morphological abnormalities of the sperm flagella phenotype (He et al., 2023; Fan et al., 2022; Zhang et al., 2021). Both genes were significantly downregulated in patients with impaired spermatogenesis across included studies, making them interesting candidates for future research of male infertility.
5 Conclusion
In the first step, we identified 257 differentially expressed genes, which are according to KEGG, mainly involved in p53 signaling and steroidogenesis pathways. In the second step we identified seven genes demonstrating both differential expression and evidence of rare, putative pathologic genetic variants (CYP11A1,CYP17A1,RSPH3,TSGA10,AKAP4,CCIN, NDNF). With the integratomic approach, we contributed to building evidence for genetic etiology and mechanisms in male infertility.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Ethics statement
The studies involving humans were approved by National medical ethics committee (reference number: 73/05/12). The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
AH contributed to research design, performed analyses and wrote the manuscript. AM performed analyses. DP performed analyses. BZ conducted the clinical andrological examination. TK revised the manuscript. BP designed and supervised the study and revised the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This study was supported by grant P3-0326 from the Slovenian Research Agency.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
ReferencesAsgariR.BakhtiariM.RezazadehD.YaraniR.EsmaeiliF.MansouriK. (2021). TSGA10 as a potential key factor in the process of spermatid differentiation/maturation: deciphering its association with autophagy pathway. Reprod. Sci.28, 3228–3240. 10.1007/s43032-021-00648-6BehnamB.ChahlaviA.PattisapuJ.WolfeJ. (2009). TSGA10 is specifically expressed in astrocyte and over-expressed in brain tumors. Avicenna J. Med. Biotechnol.1 (3), 161–6.CooksonW.LiangL.AbecasisG.MoffattM.LathropM. (2009). Mapping complex disease traits with global gene expression. Nat. Rev. Genet.10 (3), 184–194. 10.1038/nrg2537DomeniceS.MachadoA. Z.FerreiraF. M.Ferraz-de-SouzaB.LerarioA. M.LinL. (2016). Wide spectrum of NR5A1-related phenotypes in 46,XY and 46,XX individuals. Birth Defects Res. Part C, Embryo Today Rev. [online]108 (4), 309–320. 10.1002/bdrc.21145EdgarR.DomrachevM.LashA. E. (2002). Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res.30 (1), 207–10.EllisP. J. I.FurlongR. A.ConnerS. J.Kirkman-BrownJ.AfnanM.BarrattC. (2007). Coordinated transcriptional regulation patterns associated with infertility phenotypes in men. J. Med. Genet.44 (8), 498–508. 10.1136/jmg.2007.049650FakhroK. A.ElbardisiH.ArafaM.RobayA.Rodriguez-FloresJ. L.Al-ShakakiA. (2018). Point-of-care whole-exome sequencing of idiopathic male infertility. Genet. Med. Official J. Am. Coll. Med. Genet.20, 1365–1373. 10.1038/gim.2018.10FanY.HuangC.ChenJ.ChenY.WangY.YanZ. (2022). Mutations in CCIN cause teratozoospermia and male infertility. Sci. Bull.67 (20), 2112–2123. 10.1016/j.scib.2022.09.026FeigC.KirchhoffC.IvellR.NaetherO.SchulzeW.SpiessA.-N. (2006). A new paradigm for profiling testicular gene expression during normal and disturbed human spermatogenesis. MHR Basic Sci. reproductive Med.13 (1), 33–43. 10.1093/molehr/gal097HeJ.LiuQ.WangW.SuL.MengL.TanC. (2023). Novel homozygous variant of CCIN causes male infertility owing to the abnormal sperm head with a nuclear subsidence phenotype. Clin. Genet.103 (4), 495–497. 10.1111/cge.14281HodžićA.MaverA.Plaseska‐KaranfilskaD.RistanovićM.NoveskiP.ZornB. (2020). De novo mutations in idiopathic male infertility—a pilot study. Andrology9 (1), 212–220. 10.1111/andr.12897HoustonB. J.Riera-EscamillaA.WyrwollM. J.Salas-HuetosA.XavierM. J.NagirnajaL. (2021). A systematic review of the validated monogenic causes of human male infertility: 2020 update and a discussion of emerging gene-disease relationships. Hum. Reprod. Update28 (1), 15–29. 10.1093/humupd/dmab030HwangK.YatsenkoA. N.JorgezC. J.MukherjeeS.NalamR. L.MatzukM. M. (2010). Mendelian genetics of male infertility. Ann. N. Y. Acad. Sci.1214 (1), E1–E17. 10.1111/j.1749-6632.2010.05917.xJamsaiD.O’BryanM. K. (2010). Mouse models in male fertility research. Asian J. Androl.13 (1), 139–151. 10.1038/aja.2010.101JeansonL.CopinB.PaponJ.-F.Dastot-Le MoalF.DuquesnoyP.MontantinG. (2015). RSPH3 mutations cause primary ciliary dyskinesia with central-complex defects and a near absence of radial spokes. Am. J. Hum. Genet.97 (1), 153–162. 10.1016/j.ajhg.2015.05.004KallaliW.GrayE.MehdiM. Z.LindsayR.MetherellL. A.BuonocoreF. (2020). Long-term outcome of partial P450 side-chain cleavage enzyme deficiency in three brothers: the importance of early diagnosis. Eur. J. Endocrinol.182, K15–K24. 10.1530/EJE-19-0696KokR. C.TimmermanM. A.WolffenbuttelK. P.DropS. L. S.de JongF. H. (2010). Isolated 17,20-lyase deficiency due to the cytochrome b5 mutation W27X. J. Clin. Endocrinol. Metabolism95 (3), 994–999. 10.1210/jc.2008-1745KolliV.KimH.TorkyA.LaoQ.TatsiC.MallappaA. (2019). Characterization of the CYP11A1 nonsynonymous variant p.E314K in children presenting with adrenal insufficiency. J. Clin. Endocrinol. Metabolism104 (2), 269–276. 10.1210/jc.2018-01661KrauszC.Riera-EscamillaA. (2018). Genetics of male infertility. Nat. Rev. Urol.15 (6), 369–384. 10.1038/s41585-018-0003-3KuiF.YeH.ChenX.ZhangJ. (2019). Microarray meta‐analysis identifies candidate genes for human spermatogenic arrest. Andrologia51 (8), e13301. 10.1111/and.13301Lara-VelazquezM.Perdomo-PantojaA.BlackburnP. R.GassJ. M.CaulfieldT. R.AtwalP. S. (2017). A novel splice site variant inCYP11A1intranswith the p.E314K variant in a male patient with congenital adrenal insufficiency. Mol. Genet. Genomic Med.5 (6), 781–787. 10.1002/mgg3.322LuoG.HouM.WangB.LiuZ.LiuW.HanT. (2020). Tsga10 is essential for arrangement of mitochondrial sheath and male fertility in mice. Andrology9, 368–375. 10.1111/andr.12889MassartA.LissensW.TournayeH.StouffsK. (2011). Genetic causes of spermatogenic failure. Asian J. Androl.14 (1), 40–48. 10.1038/aja.2011.67MoherD.LiberatiA.TetzlaffJ.AltmanD. G.The PRISMA Group (2009). Preferred reporting Items for systematic reviews and meta-analyses: the PRISMA statement. PLoS Med.6 (7), e1000097. 10.1371/journal.pmed.1000097OkadaH.TajimaA.ShichiriK.TanakaA.TanakaK.InoueI. (2008). Genome-wide expression of azoospermia testes demonstrates a specific profile and implicates ART3 in genetic susceptibility. PLoS Genet.4 (2), e26. 10.1371/journal.pgen.0040026ParajesS.ChanA. O.ButW. M.RoseI. T.TaylorA. E.DhirV. (2012). Delayed diagnosis of adrenal insufficiency in a patient with severe penoscrotal hypospadias due to two novel P450 side-change cleavage enzyme (CYP11A1) mutations (p.R360W; p.R405X). Eur. J. Endocrinol.167 (6), 881–885. 10.1530/eje-12-0450RubtsovP.KarmanovM.SverdlovaP.SpirinP.TiulpakovA. (2009). A novel homozygous mutation inCYP11A1Gene is associated with late-onset adrenal insufficiency and hypospadias in a 46,XY patient. J. Clin. Endocrinol. Metabolism94 (3), 936–939. 10.1210/jc.2008-1118ShaY.-W.ShaY.-K.JiZ.-Y.MeiL.-B.DingL.ZhangQ. (2018). TSGA10is a novel candidate gene associated with acephalic spermatozoa. Clin. Genet.93 (4), 776–783. 10.1111/cge.13140SherbetD. P.TiosanoD.KwistK. M.HochbergZ.AuchusR. J. (2003). CYP17 mutation E305G causes isolated 17,20-lyase deficiency by selectively altering substrate binding. J. Biol. Chem.278 (49), 48563–48569. 10.1074/jbc.m307586200SpiessA.-N.FeigC.SchulzeW.ChalmelF.Cappallo-ObermannH.PrimigM. (2007). Cross-platform gene expression signature of human spermatogenic failure reveals inflammatory-like response. Hum. Reprod. Oxf. Engl.22, 2936–2946. 10.1093/humrep/dem292WuH.WangJ.ChengH.GaoY.LiuW.ZhangZ. (2020). Patients with severe asthenoteratospermia carrying SPAG6 or RSPH3 mutations have a positive pregnancy outcome following intracytoplasmic sperm injection. J. Assisted Reproduction Genet.37 (4), 829–840. 10.1007/s10815-020-01721-wYeY.WeiX.ShaY.LiN.YanX.ChengL. (2020). Loss‐of‐function mutation in TSGA10 causes acephalic spermatozoa phenotype in human. Mol. Genet. Genomic Med.8 (7), e1284. 10.1002/mgg3.1284ZhangG.LiD.TuC.MengL.TanY.JiZ. (2021). Loss-of-function missense variant of AKAP4 induced male infertility through reduced interaction with QRICH2 during sperm flagella development. Hum. Mol. Genet.31 (2), 219–231. 10.1093/hmg/ddab234