|INVITED ORIGINAL ARTICLE
|Ahead of print publication
From azoospermia to macrozoospermia, a phenotypic continuum due to mutations in the ZMYND15 gene
Zine-Eddine Kherraf1,2, Caroline Cazin1,2, Florence Lestrade3, Jana Muronova1, Charles Coutton1,4, Christophe Arnoult1, Nicolas Thierry-Mieg5, Pierre F Ray1,2
1 Univ. Grenoble Alpes, INSERM U1209, CNRS UMR 5309, Institute for Advanced Biosciences, Team Genetics Epigenetics and Therapies of Infertility, Grenoble 38000, France
2 CHU Grenoble Alpes, UM GI-DPI, Grenoble 38000, France
3 CHR Metz Thionville Hospital, 1 allée Château, METZ 57000, France
4 CHU Grenoble Alpes, UM of Chromosomic Genetics, Grenoble 38000, France
5 Univ. Grenoble Alpes, CNRS, TIMC-IMAG, Grenoble F-38000, France
|Date of Submission||26-Jul-2021|
|Date of Acceptance||18-Oct-2021|
|Date of Web Publication||11-Jan-2022|
Pierre F Ray,
Univ. Grenoble Alpes, INSERM U1209, CNRS UMR 5309, Institute for Advanced Biosciences, Team Genetics Epigenetics and Therapies of Infertility, Grenoble 38000; CHU Grenoble Alpes, UM GI-DPI, Grenoble 38000
Source of Support: None, Conflict of Interest: None
Thanks to tremendous advances in sequencing technologies and in particular to whole exome sequencing (WES), many genes have now been linked to severe sperm defects. A precise genetic diagnosis is obtained for a minority of patients and only for the most severe defects like azoospermia or macrozoospermia which is very often due to defects in the aurora kinase C (AURKC gene. Here, we studied a subject with a severe oligozoospermia and a phenotypic diagnosis of macrozoospermia. AURKC analysis did not reveal any deleterious variant. WES was then initiated which permitted to identify a homozygous loss of function variant in the zinc finger MYND-type containing 15 (ZMYND15 gene. ZMYND15 has been described to serve as a switch for haploid gene expression, and mice devoid of ZMYND15 were shown to be sterile due to nonobstructive azoospermia (NOA). In man, ZMYND15 has been associated with NOA and severe oligozoospermia. We confirm here that the presence of a bi-allelic ZMYND15 variant induces a severe oligozoospermia. In addition, we show that severe oligozoospermia can be associated macrozoospermia, and that a phenotypic misdiagnosis is possible, potentially delaying the genetic diagnosis. In conclusion, genetic defects in ZMYND15 can induce complete NOA or severe oligozoospermia associated with a very severe teratozoospermia. In our experience, severe oligozoospermia is often associated with severe teratozoospermia and can sometimes be misinterpreted as macrozoospermia or globozoospermia. In these instances, specific AURKC or dpy-19 like 2 (DPY19L2) diagnosis is usually negative and we recommend the direct use of a pan-genomic techniques such as WES.
Keywords: azoospermia; gene defect; macrozoospermia; male infertility
Article in PDF
|How to cite this URL:|
Kherraf ZE, Cazin C, Lestrade F, Muronova J, Coutton C, Arnoult C, Thierry-Mieg N, Ray PF. From azoospermia to macrozoospermia, a phenotypic continuum due to mutations in the ZMYND15 gene. Asian J Androl [Epub ahead of print] [cited 2022 Jan 25]. Available from: https://www.ajandrology.com/preprintarticle.asp?id=335651
| Introduction|| |
For many years, the only genetic tests carried out to investigate the genetic causes of male infertility were, for patients with nonobstructive azoospermia (NOA), the realization of a karyotype and the search for microdeletions of the azoospermia factor (AZF) regions in the Y chromosome, and for subject with obstructive azoospermia, the search for mutations in the cystic fibrosis transmembrane conductance regulator (CFTR) gene. Then, genetic defects were identified in critical genes responsible for specific and severe forms of teratozoospemia like aurora kinase C (AURKC) in macrozoospermia, dpy-19 like 2 (DPY19L2) in globozoospermia or Sad1 and UNC84 domain containing 5 (SUN5) in acephalic spermatozoa syndrome. Macrozoospermia, also known as large-headed spermatozoa is a rare cause of male infertility characterized by the presence of a majority of large-headed multi-flagellar spermatozoa in the ejaculate., Since its first description in 1977, numerous studies were conducted to better describe its pathophysiology. In 1996, Yurov et al. studied an infertile man affected by macrozoospermia using fluorescent in situ hybridization (FISH) and observed a high level of large headed spermatozoa presenting a diploid content in contrast to normal-sized sperm cells that had an apparent haploid content. Since this publication, other teams performed similar works and confirmed that large-headed spermatozoa had an abnormal chromosomal content and were aneuploid.,,,, In 2007, Dieterich et al. investigated by genome-wide microsatellite mapping a cohort of 14 North-African patients with macrozoospermia and identified a homozygous deletion of a single base pair in the AURKC gene (c.144delC; p.Leu49TrpfsTer23) in all studied subjects. Further studies from the same group demonstrated that all spermatozoa from patients with an AURKC gene defects were tetraploid, indicating that both meiotic divisions I and II were affected by incomplete segregation of homologous chromosomes during meiosis I and sister chromatids during meiosis II and associated with failure of nuclear cleavage and cytokinesis. AURKC is by far the most frequent gene found to be associated with macrozoospermia, as 72 out of 87 men (83%) with macrozoospermia were found to carry bi-allelic AURKC variants, and no other gene has so far been clearly associated with this phenotype in human. AURCK genetic diagnosis is one of the rare genetic diagnoses for male infertility with a clear-cut clinical relevance as it was shown that all spermatozoa are tetraploid and no pregnancy has ever been reported following intracytoplasmic sperm injection (ICSI) for men with AURKC bi-allelic loss of function variant. Following a positive AURKC diagnosis, patients should thus be orientated towards sperm donation or adoption. For subjects with a pure phenotype of macrozoospermia, with a fair sperm concentration (>5 × 106 ml−1), and with 100% abnormal spermatozoa with much enlarged heads and a high frequency of multiple flagella (often 4 flagella), we noticed that AURKC diagnosis efficiency was very high. On the contrary, diagnosis efficiency dropped for patients with severe (<1 × 106 ml−1) oligozoospermia (SO) and mixed morphological anomalies.
Here, the studied patient had some large-headed spermatozoa and was addressed to Grenoble molecular genetic laboratory for AURKC diagnosis. However, he did not present a pure phenotype of macrozoospermia as he also showed a severe oligozoospermia, with mixed morphological anomalies and a very high rate of necrospermia. Not surprisingly, AURKC diagnosis was negative and we pursued the genetic analysis by realizing whole exome sequencing (WES).
| Participant and Methods|| |
Subject and initial investigations
A single subject of Turkish origin was addressed to Grenoble Hospital, France from the Metz Infertility Center, France, for the genetic investigation of macrozoospermia. Semen samples were obtained by masturbation after 2–7 days of sexual abstinence. Semen samples were incubated for liquefaction at 37°C for 30 min; ejaculate volume and pH together with sperm concentration, morphology and motility were evaluated according to the World Health Organization (WHO) guidelines. Sperm vitality was assessed by eosin staining. Informed consent was obtained according to local protocols and the principles of the Declaration of Helsinki. The study was approved by Grenoble Ethics Committee, and samples were then stored in the Fertithèque collection declared to the French Ministry of Health (DC-2015-2580) and the French Data Protection Authority (DR-2016-392).
The whole AURKC gene was sequenced as described previously, but no deleterious variants were identified. Whole exome sequence was subsequently realized.
WES and variants filtering
Genomic DNA was isolated from blood samples and sequenced by Novogene (Cambridge, UK). Coding regions and intron/exon boundaries were sequenced after enrichment using SureSelect Human All Exon V6 (Agilent).
An alignment-ready GRCh38 reference genome (including ALT, decoy and HLA) was produced using “run-gen-ref hs38DH” from Heng Li's bwakit package (https://github.com/lh3/bwa). The exomes were analyzed using a bioinformatics pipeline developed in-house. The pipeline consists of two modules, both distributed under the GNU General Public License version 3.0 and available on github.
The first module (https://github.com/ntm/grexome-TIMC-Primary) takes FASTQ files as input and produces a single merged GVCF file, as follows. Adaptors are trimmed and low-quality reads filtered with fastp 0.20.0, reads are aligned with BWA-MEM 0.7.17, duplicates are marked using samblaster 0.1.24, and BAM files are sorted and indexed with samtools 1.9. Single nucleotide variations (SNVs) and short indels are called from each BAM file using strelka 2.9.1016 to produce individual GVCF files. These are finally merged with mergeGVCFs_strelka.pl to obtain a single multi-sample GVCF, which combines all exomes available in our laboratory.
The second module (https://github.com/ntm/grexome-TIMC-Secondary) takes each merged GVCF as input and produces annotated analysis-ready TSV files. This is achieved by performing up to 15 streamlined tasks, including the following. Low-quality variant calls (DP<10, GQ<20, or <15% of reads supporting the ALT allele) are discarded. Variant Effect Predictor is used to annotate the variants and predict their impact, allowing to filter low-impact (MODIFIER) variants and/or prioritize high-impact ones (e.g., stop-gain or frameshift variants). Variants with a minor allele frequency greater than 1% in gnomAD v2.0 or 3% in 1000 Genomes Project phase 3 are filtered. Additional information can be found in the study of Arafah et al.
Sanger verification of the variant
The zinc finger MYND-type containing 15 (ZMYND15) variant identified by WES was subjected to a Sanger verification using an ABI 3500XL (Applied Biosystems). Analyses were performed using SeqScape software (Applied Biosystems). PCR amplification was achieved using the following primers: ZMYND15_F_1, TACCCGCTGACCGTGTACTA; ZMYND15_R_1, AAATGCTGCTCGTCGCTTTC.
| Results|| |
Couple infertility investigation and treatment
The subject and his spouse had a long history of infertility. The couple suffered from primary infertility since 2000. An initial sperm investigation was carried out in 2002. Sperm volume was normal (3.2 ml) with a normal pH (7.6) and viscosity. Sperm analysis highlighted a severe oligozoospermia (0.6 × 106 ml−1) with 100% teratozoospermia, with mainly large-headed spermatozoa, associated with 0% sperm motility correlated with a very high necrospermia (<10% vitality). A high number of round germ cells were detected (0.4 × 106 ml−1) with few leukocytes (100 ml−1). The first ICSI attempt was realized. The smaller spermatozoa (as normal as possible) were selected for injection but with no sign of motility. A total of 10 mature oocytes were injected but no embryos could be obtained. Sperm activation with Pentoxifylline did not permit to initiate any sperm movement. A testicular sperm extraction (TESE) was then carried out with the hope of obtaining viable, motile spermatozoa to achieve a successful ICSI. Unfortunately, no motile spermatozoa could be identified and the TESE was unsuccessful.
In 2010, the couple came back to initiate some genetic investigations. Due to high proportion of large-headed spermatozoa, a specific diagnosis of the AURKC gene was initially requested.
Identification of a homozygous frameshift variant in ZMYND15
All of AURKC coding exons were sequenced using the Sanger technique as previously described. No deleterious variants were identified.
Exome sequencing was performed on Patient 1 (P1) using an Agilent V6 capture kit. Overall, 83.0% and 96.0% of all nucleotides from the targeted genes were sequenced with a sequencing depth of at least twenty and ten times, respectively. As most of the gene defects responsible for male infertility identified so far have a recessive transmission, we focused our search on hemizygous or bi-allelic variants. Furthermore, we eliminated the variants scored as “low” (synonymous variant and intronic variants not predicted to impact splicing). Overall, a total of 19 genes remained as potential candidates after default filtering (Participants and Methods) and exclusion of heterozygous and “low” variants [Table 1]. A total of 7 genes had a homozygous variant and 13 carried two or three variants. For these last genes, as no familial study had been undertaken, it was not possible to know if the identified variants were allelic or bi-allelic [Table 1]. Only one gene, ZMYND15, carried a homozygous variant characterized as “high,” predicted to code for a truncated and thus nonfunctional protein. All the other identified variants had a low probability of being deleterious, mainly because they were predicted to have only a moderate impact and/or because the concerned genes are not strongly expressed in the testis and/or because of the possibility that the identified variants within a same gene are allelic.
The presence of the homozygous ZMYND15 variant (c.1520_1523delAACA; p.Lys507SerfsTer3; rs587777432) was confirmed by Sanger sequencing [Figure 1]. The identified variant is a frameshift deletion of 4 nucleotides creating a stop codon in exon 9 of 14 of the ZMYND15 gene [Figure 2]. As the same variant has already been found in three brothers with azoospermia and other ZMYND15 homozygous loss of function variants have been identified in subjects with SO, we are confident that the identified variant is responsible for P1's infertility.
|Figure 1: Sanger sequencing highlighting the deletion of nucleotides c.1520_1523 in ZMYND15 exon 9. Electropherograms from Sanger sequencing indicating the homozygous state of the identified variant c.1520_1523del. The deleted nucleotides are highlighted in gray. ZMYND15: zinc finger MYND-type containing 15.|
Click here to view
|Figure 2: Location and impact of identified ZMYND15 variant. Location of the variant within ZMYND15 gene (transcript ENST00000269289.10) and protein ENSP00000269289.6. Coiled-coil domain (109–129 AA) is indicated by blue box; zinc finger, MYND-type domain (313–359 AA) is indicated in purple. ZMYND15: zinc finger MYND-type containing 15; AA: amino acid.|
Click here to view
| Discussion|| |
Here, we studied a subject (P1) with long history of infertility presenting with SO, macrozoospermia and a complete absence of sperm motility (OATS). Most patients with a typical macrozoospermia with 100% sperm anomalies and multiple flagella carry AURKC variants. For these subjects, ICSI should be ruled out as all spermatozoa are tetraploid. For P1, AURKC analysis was negative, ruling out a typical diagnosis of macrozoospermia. Exome sequencing revealed a homozygous loss of function variant in ZMYND15, a gene already described in relation with azoospermia and SO. Zmynd15 has been described in mice to code for a histone deacetylase-dependent transcriptional repressor essential for spermiogenesis and male fertility and was shown to control the temporal expression of haploid cell genes during spermiogenesis. In mice, inactivation of Zmynd15 results in early activation of transcription of important haploid genes including Prm1, Tnp1, Spem1, and Catsper3 resulting in late spermatids loss and complete azoospermia. The presence of ZMYND15 loss of function homozygous variants has also been reported to induce male infertility in human. In fact, the same homozygous variant (c.1520_1523delAACA) has already been described in three azoospermic brothers of Turkish origin. As the variant is extremely rare and P1 is also of Turkish origin, it is likely that these individuals have a common ancestor. Recently, a large cohort of 414 infertile Chinese patients composed of 219 subjects with SO and 195 subjects with NOA was screened by WES. Homozygous loss of function variants were identified in three unrelated patients, all three presenting with SO and not NOA. The three identified variants are located in exons 6 (one subject) and 10 (two subjects), relatively close to the variant identified here in P1 present in exon 9 [Figure 2].
This is therefore the third report of homozygous loss of function variants in ZMYND15 in infertile men. Surprisingly, in each publication, the described phenotype was different: the three Turkish brothers had azoospermia, whereas our patient with the same mutation had SO with macrozoospermia, and the three Chinese subjects with different variants were all described to have SO but with no precise description of their sperm count or the presence of morphological defects. Also, in addition to SO and macrozoospermia, none of the spermatozoa from P1 showed any motility due to a very high frequency of necrospermia. Testicular sperm were extracted which showed the same defects. In consequence, ICSI was unsuccessful and no pregnancy could be achieved. Again this is in contradiction with what was achieved for one of the Chinese subject, who had a successful ICSI and an ongoing pregnancy. We have frequently observed that SO is often associated with nonspecific teratozoospermia. We can conclude that a severely altered spermatogenesis will often produce few and morphologically abnormal spermatozoa. This is what was observed here, similarly to what we had observed in a patient with a defect in the serine peptidase inhibitor Kazal type 2 (SPINK2) gene.
These observations highlight that even an identical genetic event can lead to a different phenotype (NOA, SO, and oligoasthenoteratozoospermia [OATS]). This confirms that there is a phenotypic continuum between these sperm defects. This is especially true for postmeiotic genes such as ZMYND15 which can lead to NOA or SO, potentially depending on the patient's genetic background and/or environmental factors, including living habits. As expressed before, the phenotype could also become more severe with age and SO/OATS could become NOA in an aging subject.
At present, treatment counseling for patients with ZMYND15 defects cannot be formally established. Here, we observed a very high rate of necrospermia which was not compatible with a successful ICSI. This was not reported in the Chinese patients with SO, one of whom was reported to have a successful ICSI. Further cases should be reported to assess if ICSI or TESE followed by ICSI should be proposed to subjects harboring ZMYND15 defects.
| Author Contributions|| |
Z-EK, C Cazin, CA, and PFR analyzed the data and wrote the manuscript. Z-EK, C Cazin, C Coutton, JM and NT-M performed and analyzed the genetic data. FL provided clinical samples and data. Z-EK and PFR designed the study, supervised all molecular laboratory work, had full access to all of the data in the study, and took responsibility for the integrity of the data and its accuracy. All authors read and approved the final manuscript.
| Competing Interests|| |
All authors declare no competing interests.
| Acknowledgments|| |
This work was financed in part by the French Research Agency: grant to PFR (FLAGEL-OME: ANR-19-CE17-0014).
| References|| |
Nistal M, Paniagua R, Herruzo A. Multi-tailed spermatozoa in a case with asthenospermia and teratospermia. Virchows Arch B Cell Pathol
1977; 26: 111–8.
Ray PF, Toure A, Metzler-Guillemain C, Mitchell MJ, Arnoult C, et al
. Genetic abnormalities leading to qualitative defects of sperm morphology or function. Clin Gene
2017; 91: 217–32.
Yurov YB, Saias MJ, Vorsanova SG, Erny R, Soloviev IV, et al
. Rapid chromosomal analysis of germ-line cells by FISH: an investigation of an infertile male with large-headed spermatozoa. Mol Hum Reprod
1996; 2: 665–8.
Benzacken B, Gavelle FM, Martin-Pont B, Dupuy O, Lièvre N, et al
. Familial sperm polyploidy induced by genetic spermatogenesis failure: case report. Hum Reprod
2001; 16: 2646–51.
Devillard F, Metzler-Guillemain C, Pelletier R, DeRobertis C, Bergues U, et al.
Polyploidy in large-headed sperm: FISH study of three cases. Hum Reprod
2002; 17: 1292–8.
Dieterich K, Rifo RS, Faure AK, Hennebicq S, Amar BB, et al
. Homozygous mutation of AURKC
yields large-headed polyploid spermatozoa and causes male infertility. Nat Genet
2007; 39: 661–5.
Mateu E, Rodrigo L, Prados N, Gil-Salom M, Remohí J, et al
. High incidence of chromosomal abnormalities in large-headed and multiple-tailed spermatozoa. J Androl
2006; 27: 6–10.
Perrin A, Morel F, Moy L, Colleu D, Amice V, et al
. Study of aneuploidy in large-headed, multiple-tailed spermatozoa: case report and review of the literature. Fertil Steril
2008; 90: e13–7.
Dieterich K, Zouari R, Harbuz R, Vialard F, Martinez D, et al
. The aurora kinase C c.144delC mutation causes meiosis I arrest in men and is frequent in the North African population. Hum Mol Genet
2009; 18: 1301–9.
Khelifa MB, Coutton C, Blum MG, Abada F, Harbuz R, et al
. Identification of a new recurrent aurora kinase C mutation in both European and African men with macrozoospermia. Hum Reprod
2012; 27: 3337–46.
Cooper TG, Noonan E, von Eckardstein S, Auger J, Baker HW, et al.
World Health Organization reference values for human semen characteristics. Hum Reprod Update
2010; 16: 231–45.
Chen S, Zhou Y, Chen Y, Gu J. Fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics
2018; 34: i884–90.
Li H. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. arXiv:1303.3997 [q-bio]
Faust GG, Hall IM. SAMBLASTER: fast duplicate marking and structural variant read extraction. Bioinformatics
2014; 30: 2503–5.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, et al
. The sequence alignment/map format and SAMtools. Bioinformatics
2009; 25: 2078–9.
Kim S, Scheffler K, Halpern AL, Bekritsky MA, Noh E, et al
. Strelka2: fast and accurate calling of germline and somatic variants. Nat Methods
2018; 15: 591–4.
McLaren W, Gil L, Hunt SE, Riat HS, Ritchie GR, et al
. The ensembl variant effect predictor. Genome Biol
2016; 17: 122.
Arafah K, Lopez F, Cazin C, Kherraf ZE, Tassistro V, et al
. Defect in the nuclear pore membrane glycoprotein 210-like gene is associated with extreme uncondensed sperm nuclear chromatin and male infertility: a case report. Hum Reprod
2021; 36: 693–701.
Ayhan Ö, Balkan M, Guven A, Hazan R, Atar M, et al
. Truncating mutations in TAF4B
causing recessive azoospermia. J Med Genet
2014; 51: 239–44.
Hu TY, Zhang H, Meng LL, Yuan SM, Tu CF, et al
. Novel homozygous truncating variants in ZMYND15
causing severe oligozoospermia and their implications for male infertility. Hum Mutat
2021; 42: 31–6.
Yan W, Si Y, Slaymaker S, Li J, Zheng H, et al
encodes a histone deacetylase-dependent transcriptional repressor essential for spermiogenesis and male fertility. J Biol Chem
2010; 285: 31418–26.
Kherraf ZE, Christou-Kent M, Karaouzene T, Amiri-Yekta A, Martinez G, et al
. SPINK2 deficiency causes infertility by inducing sperm defects in heterozygotes and azoospermia in homozygotes. EMBO Mol Med
2017; 9: 1132–49.2
[Figure 1], [Figure 2]