Skip to main content

Whole exome sequencing analysis of 167 men with primary infertility

Abstract

Background

Spermatogenic failure is one of the leading causes of male infertility and its genetic etiology has not yet been fully understood.

Methods

The study screened a cohort of patients (n = 167) with primary male infertility in contrast to 210 normally fertile men using whole exome sequencing (WES). The expression analysis of the candidate genes based on public single cell sequencing data was performed using the R language Seurat package.

Results

No pathogenic copy number variations (CNVs) related to male infertility were identified using the the GATK-gCNV tool. Accordingly, variants of 17 known causative (five X-linked and twelve autosomal) genes, including ACTRT1, ADAD2, AR, BCORL1, CFAP47, CFAP54, DNAH17, DNAH6, DNAH7, DNAH8, DNAH9, FSIP2, MSH4, SLC9C1, TDRD9, TTC21A, and WNK3, were identified in 23 patients. Variants of 12 candidate (seven X-linked and five autosomal) genes were identified, among which CHTF18, DDB1, DNAH12, FANCB, GALNT3, OPHN1, SCML2, UPF3A, and ZMYM3 had altered fertility and semen characteristics in previously described knockout mouse models, whereas MAGEC1,RBMXL3, and ZNF185 were recurrently detected in patients with male factor infertility. The human testis single cell-sequencing database reveals that CHTF18, DDB1 and MAGEC1 are preferentially expressed in spermatogonial stem cells. DNAH12 and GALNT3 are found primarily in spermatocytes and early spermatids. UPF3A is present at a high level throughout spermatogenesis except in elongating spermatids. The testicular expression profiles of these candidate genes underlie their potential roles in spermatogenesis and the pathogenesis of male infertility.

Conclusion

WES is an effective tool in the genetic diagnosis of primary male infertility. Our findings provide useful information on precise treatment, genetic counseling, and birth defect prevention for male factor infertility.

Peer Review reports

Introduction

Worldwide, infertility affects 8–12% of fertile-aged couples, and male factor infertility accounts for approximately 50% of the condition [1, 2]. Primary spermatogenic failure is a major cause of male infertility [3]. Taking advantage of next-generation sequencing technologies (e.g., whole exome sequencing, whole genome sequencing), a growing number of monogenic variants have been associated with human spermatogenic failure [4, 5].

Programmed differentiation of spermatogonial stem cells (SSCs) involves over 2,000 genes enriched in human testis [6]. However, only a small proportion of these genes have been linked to spermatogenic failure/defects [7]. The genetic causes of male factor infertility remain largely unknown.

In the present study, we retrospectively analyzed the blood samples of 167 men with primary infertility using whole exome sequencing (WES) combined with gene expression analysis based on single cell sequencing data [8]. The results are expected to increase our knowledge on monogenic variants associated with spermatogenic failure/ defects, as well as to provide new information on gene diagnosis, genetic counseling, and individualized treatment for primary male infertility.

Materials and methods

Clinical samples

Patients with primary male infertility and altered semen parameters (e.g., azoospermia, asthenozoospermia, oligozoospermia, oligoasthenoteratozoospermia) according to the WHO guidelines (2010) were enrolled for the study at the cytogenetics laboratory of a tertiary medical center. The eligible patients had a normal karyotype and a negative result of AZF deletion test. Patients with a medical history of ejaculatory duct obstructions, cryptorchidism, epididymis and vas deferens, orchitis, hypogonadism, chemotherapy or radiation therapy were excluded for the study.The blood samples collected from the veins were stored at -80 C until use. .

WES

According to the manufacturer’s instructions, genomic DNA was isolated from peripheral blood leukocytes of all participants using a QIAamp mini DNA kit (Qiagen, Hilden, Germany) for library construction. WES was performed at the BGI Clinical Testing Center (Shenzhen, Guangdong, China). Exons were captured using the BGI-Exome kit V4 and sequenced by BGI-seq 500 with 100 bp paired end reads at an average sequencing depth of 100× and average coverage of 99.5%. Low-quality reads were removed by SOAPnuke, and then the reads were mapped to the human genome reference (GRCh38/hg38) by the Burrows-Wheeler aligner (BWA-MEM, version 0.7.10). Variants Calling was performed using the Genome Analysis Tool Kit (GATK, version 3.3). The ANNOVAR tool was used to annotate and classify all variants.

Variants filtering and prioritization

All variants were filtered based on their frequency in public databases including gnomAD (http://gnomad.broadinstitute.org/), the 1000 Genome Project (http://browser.1000genomes.org), dbSNP (http://www.ncbi.nlm.nih.gov/snp) and BGI internal database; and variants with minor allele frequency (MAF) < 0.01 were retained. Autosomal homozygous or double heterozygous variants, and X-Chromosomal variants were considered in this study. Prioritization of candidate variants was based on a recent review on 657 gene-disease relationships (GDR) for monogenic causes of human male infertility [4], and the genes associated with male infertility in the Mouse Genome Informatics database (MGI: http://www.informatics.jax.org/). Online tools including SIFT (http://provean.jcvi.org/), PolyPhen2 (http://genetics.bwh.harvard.edu/pph2/), MutationTaster (http://www.mutationtaster.org/), CADD (https://cadd.gs.washington.edu/snv), and GERP (http://mendel.stanford.edu/SidowLab/downloads/gerp/) were used to predict the damage of candidate variants.The exome data of 210 normally fertile men served as control to further screen the candidate variants.

CNV analysis using WES data

The GATK-gCNV pipeline [9] was used to analyze the germline CNVs from WES data of 167 patients in the study. CNVs with a quality score (QS) > 1000 were considered for annotation using DECIPHER database (https://decipher.sanger.ac.uk/). In order to identify potential male infertility-related genes within the pathogenic CNVs, coding genes covered by pathogenic CNVs were filtered based on the human candidate genes for male infertility [4].

Sanger sequencing

Forward and reverse primers for PCR amplification were designed using the Primer3 program based on the human genome reference sequence GRCh38/hg38. The purified PCR product was loaded for capillary electrophoresis on an ABI 3730XL DNA analyzer (ABI PRISM, Foster City, CA, USA).

Differential expression analysis of candidate genes

Single cell sequencing data (https://ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE109037) were downloaded and subjected to quality control and filtering following the scripts provided by the submitter [8]. Gene expression analysis was performed using the R language (Version 4.2.3) Seurat package (V4.0).

Results

Clinical features

The blood samples of a total of 167 patients, with a mean age ± standard deviation of 30.4 ± 4.8 yrs and a range of 22 to 49 years, were retrospectively analyzed. Based on medical records, semen characteristics of these patients include non-obstructive azoospermia (NOA) (n = 58, 34.7%), asthenozoospermia (n = 27, 16.2%), oligoasthenzoospermia (n = 63, 37.7%), oligoasthenoteratozoospermia (n = 6, 3.6%), oligozoospermia (n = 10, 6.0%), teratozoospermia (n = 3, 1.8%). All patients had a normal karyotype and no Y-chromosome infertility.

Genetic findings

In this study, variants of the genes that have been substantially demonstrated as causative genes of human male infertility were selected in priority for diagnostic purpose. Based on a recent review on 657 gene-disease relationships (GDR) for monogenic causes of human male infertility [4] and a complementary literature search, we identified deleterious variants of 17 known causative (five X-linked and twelve autosomal) genes, including ACTRT1, ADAD2, AR, BCORL1, CFAP47, CFAP54, DNAH17, DNAH6, DNAH7, DNAH8 DNAH9, FSIP2, MSH4, SLC9C1, TDRD9, TTC21A, and WNK3 in 23 patients (Table 1). The mutation spectrum of these known causative genes in patients included missense (64.1%), splice site (17.9%), frameshift (5.1%), nonsense (5.1%), start codon (2.6%), and in-frame deletion (2.6%) mutations, with a predominance of missense mutations resulting in amino acid substitution. We used the exome data from 210 normally fertile men as control. X-linked variants that were present in both male factor infertility patients and normal fertile controls were discarded. With respect to autosomal recessive variants DNAH17 c.12488T > G p.(Phe4163Cys), c.5278G > T p.(Asp1760Tyr), and c.3601 C > T p.(Arg1201Trp), DNAH8 c.6565 A > G p.(Met2189Val), and FSIP c.3365 C > G p.(Ser1122Cys), though one or two heterozygous carriers of these variants in the control group, they were still selected as causal variants for patients with male factor infertility in this study (Table 1).

Table 1 Variants of known causative genes for male infertility detected in the study

Based on the genes that have been linked to male infertility in mouse models or recurrently affected in patients, we selected 12 candidate genes for further investigation, involving seven X-linked and five autosomal genes, namely CHTF18, DDB1, DNAH12, FANCB, GALNT3, MAGEC1, OPHN1, RBMXL3, SCML2, UPF3A, ZMYM3, ZNF185 (Table 2). The mutation spectrum of candidate genes in patients included missense (50%), nonsense (15.4%), frameshift (11.5%), splice site (19.2%) and intron/exon boundary deletion (3.8%) mutations. These variants were not found in the 210 fertile men, except two variants DNAH12 c.6286 C > T p.(Gln2096X) and: c.6021–3 A > G, which had one heterozygous carrier, respectively. Considering the autosomal recessive inheritance pattern of the DNAH12 gene, the two variants were still considered as candidate causal variants for male factor infertility.

Table 2 Variants of candidate genes for male infertility

The results of the Sanger sequencing validation for the identified variants are shown in the Supplementary Data Sheet 1.

A total of 29 CNVs were identified using the GATK-gCNV tool, and none of them were directly related to male infertility in DECIPHER (accessed on August 27th, 2024; Supplementary Data Sheet 2). In one patient (M1241), a heterozygous deletion of chromosome 10q26.3 resulted in loss of a part of the spermatogenic failure gene SYCE1. This CNV was not considered to be the cause of infertility in the patient based on the recessive inheritance pattern of SYCE1.

Differential expression of candidate genes

We explored the differential expression of candidate genes at the developmental stages of spermatogenesis by referencing the human testis single cell-sequencing database [8]. Of the 12 candidate genes, CHTF18, DDB1 and MAGEC1 were preferentially expressed in spermatogonial stem cells (Fig. 1B, C and G). DNAH12 and GALNT3 were found primarily in spermatocytes and early spermatids (Fig. 1D and F). Obviously, UPF3A was present at a high level throughout spermatogenesis except in elongating spermatids (Fig. 1K). The testicular expression profiles of the genes mentioned above underlie their potential roles in spermatogenesis and the pathophysiology of male infertility. Other candidate genes were found to be expressed only in a small part of the cells during spermatogenesis and remain to be further investigated (Fig. 1).

Fig. 1
figure 1

(A) Uniform Manifold Approximation and Projection (UMAP) plot of testicular cells. The cells are colored based on cell types. SPG, SPC and SPT stand for spermatogonia, spermatocytes and spermatids, respectively. (B-M) Expression patterns of the candidate genes on the UMAP space

Discussion

Several previous studies have explored the genetic etiology of spermatogenic failure using WES, with a diagnostic yield ranging from 3.2% (10 out of 314 patients) to 47.6% (10 out of 21 patients) [31, 37,38,39,40,41]. In a more recent study that screened likely pathogenic and pathogenic (LP/P) variants of a male infertility gene panel in 521 patients with idiopathic primary spermatogenic failure, a molecular diagnosis rate of 12% (n = 64) was achieved [42]. In the present study, we identified 23 out of 167 patients with primary male infertility to carry candidate variants of known causative genes for spermatogenic failure, giving a diagnostic rate of 13.8%. We reasoned that performance of WES for detection of causal mutations of spermatogenic failure might be affected by the genetic background and inclusion criteria (e.g. pathological type, severity) of patients. Definitely, accumulating evidence on causative genes for spermatogenic failure would improve the diagnostic performance of WES.

Among the 17 known causative genes for male infertility identified in this study, AR, MSH4, and TDRD9 have been reported to be associated with azoospermia [12, 23, 43, 44]. The human androgen receptor (AR) gene, located on X chromosome q11-12, encodes a ligand-activated transcription factor that plays a key role in spermatogenesis [45]. Recently, a novel missense mutation [NM_000044.6: c.2051G > C p.(Gly684Ala)] within the ligand-binding domain of AR in azoospermic individuals from a Chinese family has been identified using WES [44]. In this study, we identified a start-loss variant of AR in a patient with NOA. MSH4 (MutShhomologue 4) is a member of the mammalian mismatch repair (MMR) gene family, and plays a vital role in meiotic chromosome synapsis [46]. Homozygous missense [NM_002440.4:c.1913 C > T p. (Pro638Leu); c.2261 C > T p. (Ser754Leu)], homozygous deletion frameshift [c.805_812del p.(Val269Glnfs*15); c.2220_2223del p.(Lys741Argfs*2)], and compound heterozygous [c.1950G > A p.(Trp650X) and c.2179delG p.(Asp727Metfs*11); c.244G > A p.(Gly82Ser) and c.670delT p.(Leu224Cysfs*3)] variants of MSH4 have been reported in patients with NOA [37, 43]. Here, we have identified a homozygous splice-site variant in a patient with NOA. The mouse Tdrd9 gene encodes an ATPase/DExH-type helicase, which has been linked to male sterility and meiotic failure [47]. A previous study has identified a homozygous 4 bp deletion frameshift mutation in TDRD9 [c.720_723 del TAGT p.(Ser241Profs*4)] as the causative mutation in five azoospermic infertile men of a large consanguineous Bedouin family [23]. In our study, we identified a homozygous missense variant of TDRD9 in a NOA patient. Together, our findings expand the mutation spectrum of these azoopsermia genes. Interestingly, it has been shown that biallelic deleterious mutations of CFAP54 can induce severe MMAF and NOA in humans [26].With respect to the recently identified MMAF causative gene CFAP47 [14], one patient (M1294) hemizygous for a missense variant of the gene had NOA from our study. This finding may expand the phenotype of male infertility caused by CFAP47. However, further evidence is required to confirm the relationship between CFAP47 and NOA.

Five of the 17 known causative genes, namely DNAH17, DNAH6, DNAH7, DNAH8, and DNAH9, belong to the DNAH family, members of which encode the axonemal dynein heavy chains [48]. Dyneins are components of the inner and outer dynein arms attached to the doublets of the axonemal microtubules in the motile cilia and sperm flagella. Mutations in the DNAH genes are known to be associated with cilia-related disorders, such as immotile cilia syndrome [49].Using WES, a recent study analyzed 90 cases of Chinese patients with MMAF and identified two patients with biallelic heterozygous missense variants of DNAH8 [c.11771 C > T p.(Thr3924Met) plus c.6689 A > G p.(Lys2230Arg); c.9427 C > T p.(Arg3143Cys) plus c.12721G > A p.(Ala4241Thr)]. Meanwhile, reanalysis of the exome sequencing data of 167 MMAF patients from France, Iran, and North Africa identified a patient of homozygous frameshift mutation [c.6962_6968del p.(His2321Profs*4)] in DNAH8 [18]. In this study, we identified two patients (M1237 and M1761) with asthenozoopsermia carrying double heterozygous variants of DNAH8. In the study by Whitfield et al., five patients with asthenozoospermia were found to carry variants of DNAH17 [15]. In our study, three patients with azoospermia were found to carry double heterozygous missense mutations in DNAH17. In line with our finding, a most recent study also detected multiple-heterozygous damaging variants in DNAH17 in Japanese men with isolated non-obstructive azoospermia [50]. Regarding DNAH9, we noticed the difference between the NOA phenotype of the patient M1470 in this study and the asthenozoospermia phenotype of the patients described previously [19]. Further investigation is needed to clarify this issue.

In the present study, we have identified 12 candidate genes for human spermatogenic failure/defects. Based on testicular single cell sequencing data, six of the 12 candidate genes (CHTF18, DDB1, DNAH12, GALNT3, MAGEC1, and UPF3A) are deferentially expressed at spermatogenic stages, suggestive of distinct roles of these genes in spermatogenesis. Variants of DNAH12 were previously identified in patients with MMAF or sperm motility disorder [29, 31]. In accordance with these findings, a homozygous nonsense variant of DNAH12 was identified in a patient (M842) with asthenozoospermia in our study, adding evidence to DNAH12 being a causative gene for male infertility. However, the patient M867 who carried double heterozygous splice-site variants in DNAH12 had azoospermia by routine semen analysis, seeming to be inconsistent with the asthenozoospermia and MMAF phenotype of patients with homozygous or compound heterozygous variants in DNAH12 [30]. We reason that the possibility of DNAH12 leading to azoospermia could not be excluded, considering that single cell sequencing data revealed its expression in spermatocytes and early spermatids, and that Dnah12−/− male mice had a dramatically decreased sperm number compared with the wild type [30]. In addition, the MMAF causative genes DNAH6 and DNAH17 have also been associated with azoospermia [50, 51]. To confirm whether DNAH12 is a causative for NOA, more azoopsermic patients with pathogenic variants in DNAH12 are anticipated to be identified in future studies.

Interestingly, loss of Chtf18 in mice leads to oligospermia due to defective meiotic recombination [26]; but the patient (M1286) double heterozygous for missense variants in CHTF18 from our study was diagnosed as NOA. The Upf3a knockout mice were described to have significant reduction in sperm count and a defect in spermatocyte progression [35], whereas the patient M847 carrying a homozygous frameshift variant of UPF3A from our study had a semen characteristic of asthenozoospermia. The phenotypic differences may be attributed to genetic differences between humans and mouse strains. Though our findings provide important information on the roles of CHTF18 and UPF3A in human spermatogenesis, whether they are casual genes for spermatogenic failure needs to be further verified.

There are some limitations to the study. First, the sample size of the current study was relatively small. Secondly, this study was a retrospective study and parental samples were not available. Therefore, the maternal and paternal origins of the candidate variants identified in patients are unclear. Furthermore, the patients did not undergo testicular puncture or biopsy, so it is impossible to assess the pathological types in the testes and detect the expression of related genes.

Taken together, our study shows that WES is an effective tool in the genetic diagnosis of primary male infertility. Our study increases a comprehensive and in-depth understanding of the genetic factors involved in spermatogenesis and has meaning for the diagnosis, clinical treatment and genetic counseling of male infertility.

Data availability

The data that support the findings of this study have been deposited into CNGB Sequence Archive (CNSA) of China National GeneBank DataBase (CNGBdb) with accession number CNP0005865.

References

  1. Agarwal A, Baskaran S, Parekh N, Cho CL, Henkel R, Vij S, et al. Male infertility. Lancet. 2021;397(10271):319–33.

    Article  PubMed  Google Scholar 

  2. Minhas S, Bettocchi C, Boeri L, Capogrosso P, Carvalho J, Cilesiz NC, et al. European Association of Urology Guidelines on male sexual and Reproductive Health: 2021 update on male infertility. Eur Urol. 2021;80(5):603–20.

    Article  PubMed  Google Scholar 

  3. Cham TC, Chen X, Honaramooz A. Current progress, challenges, and future prospects of testis organoidsdagger. Biol Reprod. 2021;104(5):942–61.

    Article  PubMed  Google Scholar 

  4. Houston BJ, Riera-Escamilla A, Wyrwoll MJ, Salas-Huetos A, Xavier MJ, Nagirnaja L, et al. A systematic review of the validated monogenic causes of human male infertility: 2020 update and a discussion of emerging gene-disease relationships. Hum Reprod Update. 2021;28(1):15–29.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Nagirnaja L, Lopes AM, Charng WL, Miller B, Stakaitis R, Golubickaite I, et al. Diverse monogenic subforms of human spermatogenic failure. Nat Commun. 2022;13(1):7953.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Krausz C, Riera-Escamilla A. Genetics of male infertility. Nat Rev Urol. 2018;15(6):369–84.

    Article  CAS  PubMed  Google Scholar 

  7. Malcher A, Stokowy T, Berman A, Olszewska M, Jedrzejczak P, Sielski D, et al. Whole-genome sequencing identifies new candidate genes for nonobstructive azoospermia. Andrology. 2022;10(8):1605–24.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Hermann BP, Cheng K, Singh A, Roa-De La Cruz L, Mutoji KN, Chen IC, et al. The mammalian spermatogenesis single-cell transcriptome, from Spermatogonial Stem cells to Spermatids. Cell Rep. 2018;25(6):1650–e16671658.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Babadi M, Fu JM, Lee SK, Smirnov AN, Gauthier LD, Walker M, et al. GATK-gCNV enables the discovery of rare copy number variants from exome sequencing data. Nat Genet. 2023;55(9):1589–97.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Sha Y, Liu W, Li L, Serafimovski M, Isachenko V, Li Y, et al. Pathogenic variants in ACTRT1 cause Acephalic Spermatozoa Syndrome. Front Cell Dev Biol. 2021;9:676246.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Dai S, Liu M, Liu M, Jiang C, Yang Y, Han H, et al. Population-based genetic analysis in infertile men reveals novel mutations of ADAD family members in patients with impaired spermatogenesis. Hum Mol Genet. 2023;32(11):1814–25.

    Article  CAS  PubMed  Google Scholar 

  12. Giwercman YL, Nikoshkov A, Bystrom B, Pousette A, Arver S, Wedell A. A novel mutation (N233K) in the transactivating domain and the N756S mutation in the ligand binding domain of the androgen receptor gene are associated with male infertility. Clin Endocrinol (Oxf). 2001;54(6):827–34.

    Article  CAS  PubMed  Google Scholar 

  13. Lu C, Zhang Y, Qin Y, Xu Q, Zhou R, Cui Y, et al. Human X chromosome exome sequencing identifies BCORL1 as contributor to spermatogenesis. J Med Genet. 2021;58(1):56–65.

    Article  CAS  PubMed  Google Scholar 

  14. Liu C, Tu C, Wang L, Wu H, Houston BJ, Mastrorosa FK, et al. Deleterious variants in X-linked CFAP47 induce asthenoteratozoospermia and primary male infertility. Am J Hum Genet. 2021;108(2):309–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Whitfield M, Thomas L, Bequignon E, Schmitt A, Stouvenel L, Montantin G, et al. Mutations in DNAH17, encoding a sperm-specific axonemal outer Dynein Arm Heavy Chain, cause isolated male infertility due to Asthenozoospermia. Am J Hum Genet. 2019;105(1):198–212.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Tu C, Nie H, Meng L, Yuan S, He W, Luo A, et al. Identification of DNAH6 mutations in infertile men with multiple morphological abnormalities of the sperm flagella. Sci Rep. 2019;9(1):15864.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Wei X, Sha Y, Wei Z, Zhu X, He F, Zhang X, et al. Bi-allelic mutations in DNAH7 cause asthenozoospermia by impairing the integrality of axoneme structure. Acta Biochim Biophys Sin (Shanghai). 2021;53(10):1300–9.

    Article  CAS  PubMed  Google Scholar 

  18. Liu C, Miyata H, Gao Y, Sha Y, Tang S, Xu Z, et al. Bi-allelic DNAH8 variants lead to multiple morphological abnormalities of the sperm flagella and primary male infertility. Am J Hum Genet. 2020;107(2):330–41.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Tang D, Sha Y, Gao Y, Zhang J, Cheng H, Zhang J, et al. Novel variants in DNAH9 lead to nonsyndromic severe asthenozoospermia. Reprod Biol Endocrinol. 2021;19(1):27.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Martinez G, Kherraf ZE, Zouari R, Fourati Ben Mustapha S, Saut A, Pernet-Gallay K, et al. Whole-exome sequencing identifies mutations in FSIP2 as a recurrent cause of multiple morphological abnormalities of the sperm flagella. Hum Reprod. 2018;33(10):1973–84.

    Article  CAS  PubMed  Google Scholar 

  21. Wyrwoll MJ, van Walree ES, Hamer G, Rotte N, Motazacker MM, Meijers-Heijboer H, et al. Bi-allelic variants in DNA mismatch repair proteins MutS Homolog MSH4 and MSH5 cause infertility in both sexes. Hum Reprod. 2021;37(1):178–89.

    Article  CAS  PubMed  Google Scholar 

  22. Cavarocchi E, Whitfield M, Chargui A, Stouvenel L, Lores P, Coutton C, et al. The sodium/proton exchanger SLC9C1 (sNHE) is essential for human sperm motility and fertility. Clin Genet. 2021;99(5):684–93.

    Article  CAS  PubMed  Google Scholar 

  23. Arafat M, Har-Vardi I, Harlev A, Levitas E, Zeadna A, Abofoul-Azab M, et al. Mutation in TDRD9 causes non-obstructive azoospermia in infertile men. J Med Genet. 2017;54(9):633–9.

    Article  CAS  PubMed  Google Scholar 

  24. Liu W, He X, Yang S, Zouari R, Wang J, Wu H, et al. Bi-allelic mutations in TTC21A induce asthenoteratospermia in humans and mice. Am J Hum Genet. 2019;104(4):738–48.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Fakhro KA, Elbardisi H, Arafa M, Robay A, Rodriguez-Flores JL, Al-Shakaki A, et al. Point-of-care whole-exome sequencing of idiopathic male infertility. Genet Med. 2018;20(11):1365–73.

    Article  CAS  PubMed  Google Scholar 

  26. Tian S, Tu C, He X, Meng L, Wang J, Tang S, et al. Biallelic mutations in CFAP54 cause male infertility with severe MMAF and NOA. J Med Genet. 2023;60(8):827–34.

    Article  CAS  PubMed  Google Scholar 

  27. Berkowitz KM, Sowash AR, Koenig LR, Urcuyo D, Khan F, Yang F, et al. Disruption of CHTF18 causes defective meiotic recombination in male mice. PLoS Genet. 2012;8(11):e1002996.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Zheng W, Nazish J, Wahab F, Khan R, Jiang X, Shi Q. DDB1 regulates sertoli cell proliferation and testis cord remodeling by TGFbeta Pathway. Genes (Basel) 2019, 10(12).

  29. Li Y, Wang Y, Wen Y, Zhang T, Wang X, Jiang C, et al. Whole-exome sequencing of a cohort of infertile men reveals novel causative genes in teratozoospermia that are chiefly related to sperm head defects. Hum Reprod. 2021;37(1):152–77.

    Article  CAS  PubMed  Google Scholar 

  30. Geng H, Wang K, Liang D, Ni X, Yu H, Tang D, et al. Further evidence from DNAH12 supports favorable fertility outcomes of infertile males with dynein axonemal heavy chain gene family variants. iScience. 2024;27(7):110366.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Oud MS, Houston BJ, Volozonoka L, Mastrorosa FK, Holt GS, Alobaidi BKS, et al. Exome sequencing reveals variants in known and novel candidate genes for severe sperm motility disorders. Hum Reprod. 2021;36(9):2597–611.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Kato Y, Alavattam KG, Sin HS, Meetei AR, Pang Q, Andreassen PR, et al. FANCB is essential in the male germline and regulates H3K9 methylation on the sex chromosomes during meiosis. Hum Mol Genet. 2015;24(18):5234–49.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Miyazaki T, Mori M, Yoshida CA, Ito C, Yamatoya K, Moriishi T, et al. Galnt3 deficiency disrupts acrosome formation and leads to oligoasthenoteratozoospermia. Histochem Cell Biol. 2013;139(2):339–54.

    Article  CAS  PubMed  Google Scholar 

  34. Sakashita A, Ooga M, Otsuka K, Maezawa S, Takeuchi C, Wakayama S, et al. Polycomb protein SCML2 mediates paternal epigenetic inheritance through sperm chromatin. Nucleic Acids Res. 2023;51(13):6668–83.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Shum EY, Jones SH, Shao A, Chousal JN, Krause MD, Chan WK, et al. The antagonistic gene paralogs Upf3a and Upf3b govern nonsense-mediated RNA decay. Cell. 2016;165(2):382–95.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Hu X, Shen B, Liao S, Ning Y, Ma L, Chen J, et al. Gene knockout of Zmym3 in mice arrests spermatogenesis at meiotic metaphase with defects in spindle assembly checkpoint. Cell Death Dis. 2017;8(6):e2910.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Krausz C, Riera-Escamilla A, Moreno-Mendoza D, Holleman K, Cioppi F, Algaba F, et al. Genetic dissection of spermatogenic arrest through exome analysis: clinical implications for the management of azoospermic men. Genet Med. 2020;22(12):1956–66.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Chen S, Wang G, Zheng X, Ge S, Dai Y, Ping P, et al. Whole-exome sequencing of a large Chinese azoospermia and severe oligospermia cohort identifies novel infertility causative variants and genes. Hum Mol Genet. 2020;29(14):2451–9.

    Article  CAS  PubMed  Google Scholar 

  39. Tang D, Li K, Geng H, Xu C, Lv M, Gao Y, et al. Identification of deleterious variants in patients with male infertility due to idiopathic non-obstructive azoospermia. Reprod Biol Endocrinol. 2022;20(1):63.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Zhang H, Li W, Jiang Y, Li J, Chen M, Wang R, et al. Whole exome sequencing identifies genes Associated with Non-obstructive Azoospermia. Front Genet. 2022;13:872179.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Kherraf ZE, Cazin C, Bouker A, Fourati Ben Mustapha S, Hennebicq S, Septier A, et al. Whole-exome sequencing improves the diagnosis and care of men with non-obstructive azoospermia. Am J Hum Genet. 2022;109(3):508–17.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Lillepea K, Juchnewitsch AG, Kasak L, Valkna A, Dutta A, Pomm K, et al. Toward clinical exomes in diagnostics and management of male infertility. Am J Hum Genet. 2024;111(5):877–95.

    Article  CAS  PubMed  Google Scholar 

  43. Li P, Ji Z, Zhi E, Zhang Y, Han S, Zhao L, et al. Novel bi-allelic MSH4 variants causes meiotic arrest and non-obstructive azoospermia. Reprod Biol Endocrinol. 2022;20(1):21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Yuan Y, Xu WQ, Chen Y, Luo T, Chen HY. A Gly684Ala substitution in the androgen receptor is the cause for azoospermia in a Chinese family with mild androgen insensitivity syndrome and normal hormone levels. Front Genet. 2022;13:988202.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Smith LB, Walker WH. The regulation of spermatogenesis by androgens. Semin Cell Dev Biol. 2014;30:2–13.

    Article  CAS  PubMed  Google Scholar 

  46. Neyton S, Lespinasse F, Moens PB, Paul R, Gaudray P, Paquis-Flucklinger V, et al. Association between MSH4 (MutS homologue 4) and the DNA strand-exchange RAD51 and DMC1 proteins during mammalian meiosis. Mol Hum Reprod. 2004;10(12):917–24.

    Article  CAS  PubMed  Google Scholar 

  47. Shoji M, Tanaka T, Hosokawa M, Reuter M, Stark A, Kato Y, et al. The TDRD9-MIWI2 complex is essential for piRNA-mediated retrotransposon silencing in the mouse male germline. Dev Cell. 2009;17(6):775–87.

    Article  CAS  PubMed  Google Scholar 

  48. Silvanovich A, Li MG, Serr M, Mische S, Hays TS. The third P-loop domain in cytoplasmic dynein heavy chain is essential for dynein motor function and ATP-sensitive microtubule binding. Mol Biol Cell. 2003;14(4):1355–65.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Inaba K, Mizuno K. Sperm dysfunction and ciliopathy. Reprod Med Biol. 2016;15(2):77–94.

    Article  CAS  PubMed  Google Scholar 

  50. Muranishi Y, Kobori Y, Katoh-Fukui Y, Tamaoka S, Hattori A, Osaka A, et al. Systematic molecular analyses for 115 karyotypically normal men with isolated non-obstructive azoospermia. Hum Reprod. 2024;39(5):1131–40.

    Article  PubMed  Google Scholar 

  51. Gershoni M, Hauser R, Yogev L, Lehavi O, Azem F, Yavetz H, et al. A familial study of azoospermic men identifies three novel causative mutations in three new human azoospermia genes. Genet Med. 2017;19(9):998–1006.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We thank the patients for their cooperation and continued support for this study.

Funding

This research was supported by Natural Science Foundation of Hunan Province (2022JJ30322), and Hunan Provincial Hospital of Maternal and Child Health Care’s High-Level Talent Development Scheme.

Author information

Authors and Affiliations

Authors

Contributions

HZ collected and prepared the samples for exome sequencing, and wrote the first version of the manuscript. ZY interpreted the WES data. BN and JL collected the clinical information of the patients. SL performed the gene expression analysis and wrote the manuscript. WX designed the study and wrote the manuscript and take responsibility for the integrity of the work. Each author believes that the manuscript represents honest work. All authors read and approved the final version of manuscript.

Corresponding authors

Correspondence to Shuwei Luo or Wanqin Xie.

Ethics declarations

Ethics approval and consent to participate

The present study was approved by the Medical Ethics Committee of Hunan Provincial Maternal and Child Health Care Hospital (No. K202031), and written informed consent was obtained from all study participants.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhou, H., Yin, Z., Ni, B. et al. Whole exome sequencing analysis of 167 men with primary infertility. BMC Med Genomics 17, 230 (2024). https://doi.org/10.1186/s12920-024-02005-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12920-024-02005-3

Keywords