|INVITED ORIGINAL ARTICLE
|Ahead of print publication
Investigation of the genetic etiology in male infertility with apparently balanced chromosomal structural rearrangements by genome sequencing
Matthew Hoi Kin Chau1,2,3, Ying Li1,2, Peng Dai4, Mengmeng Shi1,2, Xiaofan Zhu4, Jacqueline Pui Wah Chung1, Yvonne K Kwok1,2, Kwong Wai Choy1,2,3,5, Xiangdong Kong4, Zirui Dong1,2,3
1 Department of Obstetrics and Gynaecology, The Chinese University of Hong Kong, Hong Kong, China
2 Shenzhen Research Institute, The Chinese University of Hong Kong, Shenzhen 518057, China
3 Hong Kong Hub of Paediatric Excellence, The Chinese University of Hong Kong, Hong Kong, China
4 Genetics and Prenatal Diagnosis Center, The First Affiliated Hospital of Zhengzhou University, Zhengzhou 450052, China
5 The Chinese University of Hong Kong-Baylor College of Medicine Joint Center for Medical Genetics, Hong Kong, China
|Date of Submission||28-Jul-2021|
|Date of Acceptance||27-Oct-2021|
|Date of Web Publication||07-Jan-2022|
Department of Obstetrics and Gynaecology, The Chinese University of Hong Kong, Hong Kong; Shenzhen Research Institute, The Chinese University of Hong Kong, Shenzhen 518057; Hong Kong Hub of Paediatric Excellence, The Chinese University of Hong Kong, Hong Kong
Source of Support: None, Conflict of Interest: None
Apparently balanced chromosomal structural rearrangements are known to cause male infertility and account for approximately 1% of azoospermia or severe oligospermia. However, the underlying mechanisms of pathogenesis and etiologies are still largely unknown. Herein, we investigated apparently balanced interchromosomal structural rearrangements in six cases with azoospermia/severe oligospermia to comprehensively identify and delineate cryptic structural rearrangements and the related copy number variants. In addition, high read-depth genome sequencing (GS) (30-fold) was performed to investigate point mutations causative of male infertility. Mate-pair GS (4-fold) revealed additional structural rearrangements and/or copy number changes in 5 of 6 cases and detected a total of 48 rearrangements. Overall, the breakpoints caused truncations of 30 RefSeq genes, five of which were associated with spermatogenesis. Furthermore, the breakpoints disrupted 43 topological-associated domains. Direct disruptions or potential dysregulations of genes, which play potential roles in male germ cell development, apoptosis, and spermatogenesis, were found in all cases (n = 6). In addition, high read-depth GS detected dual molecular findings in case MI6, involving a complex rearrangement and two point mutations in the gene DNAH1. Overall, our study provided the molecular characteristics of apparently balanced interchromosomal structural rearrangements in patients with male infertility. We demonstrated the complexity of chromosomal structural rearrangements, potential gene disruptions/dysregulation and single-gene mutations could be the contributing mechanisms underlie male infertility.
Keywords: azoospermia; balanced structural rearrangements; genome sequencing; male infertility; severe oligospermia
Article in PDF
|How to cite this URL:|
Kin Chau MH, Li Y, Dai P, Shi M, Zhu X, Wah Chung JP, Kwok YK, Choy KW, Kong X, Dong Z. Investigation of the genetic etiology in male infertility with apparently balanced chromosomal structural rearrangements by genome sequencing. Asian J Androl [Epub ahead of print] [cited 2022 Jan 25]. Available from: https://www.ajandrology.com/preprintarticle.asp?id=335265
Matthew Hoi Kin Chau, Ying Li
These authors contributed equally to this work.
| Introduction|| |
Infertility is defined as the failure to achieve pregnancy after 12 months or more of regular unprotected intercourse. Approximately 15% of couples are affected by infertility, and among them, male factor contributes to 50%. Male factor infertility is often defined by abnormal semen parameters, which could be associated with other medical conditions, developmental factors, lifestyle factors, and genetics. However, the exact nature of these associations is still largely unclear. The most common etiology of male infertility is primary testicular failure, resulting in quantitative impairment of spermatogenesis, accounting for 75% of male factor infertility. However, the etiology of primary testicular failure is unknown in about 40% of cases, in which genetic factors are thought to contribute to a significant proportion of these cases.
Genetic factors account for at least 15% of male infertility, and they have been known to contribute to all four major etiological categories, including (1) spermatogenic quantitative defects, (2) ductal obstructions or dysfunction, (3) hypothalamic–pituitary axis disturbances, and (4) spermatogenic qualitative defects. The genetic landscape of male infertility is highly complex and heterogeneous, with over 2000 genes involved in spermatogenesis and testicular function. Diagnosing the genetic causes of male infertility has clinical implications for prognosis of testicular sperm retrieval and personalizing therapy, improving the reproductive health and general health of the couple.
Men with azoospermia have the highest risk of being affected by genetic factors (25%). There are two well-known genetic causes of nonobstructive azoospermia. The most common causative genetic defects include chromosome Yq11 microdeletions of the azoospermia factor region (AZF), which can be detected in 13% of men with azoospermia. In patients with complete AZFb or AZFb + AZFc deletions, testicular biopsy usually reveals Sertoli cells only or diffuse early maturation arrest; hence, these patients must use donor sperm or adoption. In contrast, men with partial AZFa deletions often have hypospermatogenesis, with much better surgical sperm retrieval rates.
Chromosomal abnormalities including sex chromosomal abnormalities (SCA) also cause male infertility. Among them, balanced chromosomal abnormalities (BCA; such as translocation and insertion), defined as exchange of chromosomal segments between two or more chromosomes without apparent gains or losses detectable by karyotyping, account for 0.5%–1.0% of patients with severe oligospermia or azoospermia. Male infertility caused by BCA may be explained by the failure of pairing of homologous elements on derivative chromosomes during meiosis I. There is also evidence that failure of pairing of homologous elements promotes association between the quadrivalent and the X-Y bivalent, which affects sperm count. However, conception can still be achieved in a significant proportion of male BCA carriers. In addition, semen analysis showed no significant differences in BCA carriers and males with normal karyotype. In contrast, GS-based investigations of BCAs have suggested that gene disruption or dysregulation caused by structural rearrangements and cryptic complexities could also contribute to male infertility., In complex rearrangements, the number of chromosomes and breakpoints involved, the location of breakpoints, and their relative sizes are presumed to affect fertility.
Although GS has been applied to investigate the molecular breakpoints of chromosomal structural rearrangements in patients with abnormal phenotypes,,,, standard GS (with small-insert DNA) is not well suited for identifying breakpoint junctions, particularly those that are mediated by repetitive elements. To overcome this limitation, we developed an in-house mate-pair genome sequencing method which utilizes large-insert size DNA (3–8 kb) libraries, sequenced at a low read-depth. Our previous studies demonstrated that this method increased the sensitivity in the detection of BCAs and demonstrated the detection of additional cryptic and complex rearrangements in karyotypically simple chromosomal abnormalities. In this study, we aim to investigate BCA breakpoint junctions and the related copy number deletions/duplications by mate-pair GS and identify other genomic variants that potentially contribute to male infertility.
| Patients And Methods|| |
Ethics approval and case recruitment
This study was approved by The Chinese University of Hong Kong-New Territories East Cluster Clinical Research Ethics Committee (The Joint CUHK-NTEC CREC; approval No. 2020.046). All participants in this study have provided written informed consent. We studied six idiopathic nonobstructive azoospermic or severe oligospermic males with BCAs recruited from the Prince of Wales Hospital, 30–32 Ngan Shing Street, Shatin, New Territories, Hong Kong, China. BCAs were previously ascertained in these subjects by karyotyping [Table 1], and chromosome Y microdeletions have been excluded. One patient carried a simple two-way reciprocal translocation, one carried double two-way translocations, two had three-way translocations, one case had a two-way chromosomal insertion (described in our previous study), and one case had both a complex insertion and a translocation.
|Table 1: Case summary of the azoospermic/oligospermic cohort with apparently balanced structural rearrangements identified|
Click here to view
Mate-pair genome sequencing
Genomic DNA (gDNA) was extracted with the DNeasy Blood & Tissue Kit (category No. 69506; Qiagen, Hilden, Germany) from each sample. gDNA was then quantified with a Qubit dsDNA HS Assay Kit (Invitrogen, Carlsbad, CA, USA). After quality control (QC), a 1-μg aliquot of gDNA was sheared to 3–8 kb for mate-pair library construction. The libraries were pooled and sequenced to paired-end 100 bp on an MGISEQ-2000 platform (MGI Technology Co., Ltd., Shenzhen, China). A minimal of 60 million read-pairs were generated for each case, equating to an average 4 × read-depth and 60 × to 100 × physical coverage. In addition, for each sample except MI5, 400 ng of gDNA was also subjected for standard high read-depth GS. Each library was sequenced in each lane (paired-end 150 bp) on an MGISEQ-2000 platform (MGI Technology Co., Ltd.) for a minimal of 30-fold read-depth.
Variant identification and interpretation
Genomic variants were detected using our in-house bioinformatics pipelines. (1) Structural rearrangements were detected at a resolution of 10 kb. In brief, chimeric read-pairs that aligned to different chromosomes or mapped to the same chromosome with a genomic distance of >10 kb, were clustered and filtered against a control dataset of >2500 GS data (to remove systematic errors due to short-read alignment or variants with high frequencies as polymorphisms). (2) Copy number variants (CNVs) were detected using our reported increment-ratio-of-coverage method at a resolution of 50 kb.,, Rare CNVs with allele frequencies of <1% in our in-house Chinese subjects (n > 2000) were selected for review. (3) Absence of heterozygosity (AOH) was detected at a resolution of 5 Mb as previously described to investigate any uniparental disomy and parental consanguinity.
The high read-depth GS data analysis involved the following steps: fastq data QC and alignment (BWA-MEM), detection of single-nucleotide variants (SNVs), and small insertions/deletions (InDels) by HaplotypeCaller version 3.4 from the Genome Analysis Toolkit (GATK; Broad Institute, MA, USA). The variants were subsequently annotated using ANNOVAR and InterVAR with in-house and public databases.
The following criteria were used for prioritization of SNVs/InDels: (1) variants reported by ClinVAR or human gene mutation database (HGMD); (2) variants with a minor allele frequency ≤5% in the ExAC (http://exac.broadinstitute.org) and gnomAD (https://gnomad.broadinstitute.org) databases; (3) variants that are located in coding regions or exon–intron junctions; (4) variants with damaging/intolerant or splicing-change effects predicted by multiple computational algorithms (Revel, SIFT, Polyphen-2, MutationTaster, Human Splicing Finder, and MaxEntScan); and (5) variants that contain or are located in OMIM disease-causing genes. For known mutations, clinical correlation will be performed (for male infertility). Further classification of novel mutations will be performed for variants (1) located in autosomal dominant or X-linked dominant genes or (2) being homozygous or compound heterozygous in autosomal recessive or X-linked recessive genes. Potential disease-causing mutations will be selected for validation and parental confirmation when available.
Genomic variants (CNVs, AOHs, SNVs, and InDels) involving genes associated with male infertility, oligospermia or azoospermia, spermatogenesis deficiency, and ciliopathy are interpreted based on the guidelines of the American College of Medical Genetics and Genomics (ACMG) using the 5-tier classification.
Variant verification and junction annotation
Rearrangement breakpoints were resolved to the single-nucleotide level by junction-specific PCR followed by Sanger sequencing, which enabled accurate junction annotation for interpretation. Sanger sequencing results (fa format) were mapped to the reference genome (GRCh37) using UCSC BLAT. To analyze and assess the functional and phenotypic association of structural rearrangements, the breakpoint junctions were annotated for (1) direct disruption of genes, (2) disruption of regulatory elements, (3) disruption of topologically associating domain (TAD) defined for hESC cell line (3D Genome browser: http://3dgenome.fsm.northwestern.edu/index.html), and (4) genes or regulatory elements within the same TAD as the breakpoint. Furthermore, the breakpoint junction sequence features were annotated to investigate the mechanisms involved in the rearrangements.
Orthogonal validation of CNVs and AOHs was performed using the 8 × 60K Fetal DNA Chip version 2.0 (Agilent Technologies, Santa Clara, CA, USA) microarray., Sanger sequencing was performed to validate SNVs and InDels.
| Results|| |
We performed mate-pair GS on six male cases with azoospermia/severe oligospermia [Table 1] and apparently balanced chromosomal structural rearrangements previously ascertained by G-banded chromosome analysis, who were negative for chromosome Y microdeletions. Karyotyping revealed reciprocal-balanced translocation in one case (sample ID: MI1), two independent balanced translocations in one case (MI2), three-way balanced translocations in two cases (MI3 and MI4), a two-way chromosomal insertion in one case (MI5), and a complex insertion and translocation in one case (MI6; [Table 1]). Overall, mate-pair GS detected 48 structural rearrangements. Cryptic deletions were identified in two cases [Supplementary Table 1 [Additional file 1]]. Thirty-six rearrangements were precisely mapped to single-nucleotide resolution [Supplementary Table 2 [Additional file 2]] by Sanger sequencing, while the remaining 12 rearrangements were verified by gap-PCR but could not be resolved to nucleotide level due to the presence of repetitive elements.
In the case with the simplest rearrangement, two breakpoints were detected in a reciprocal translocation in case MI1. In the most complex case (MI6), 19 rearrangements were detected with deletions of four segments. Of the 48 rearrangements, 39 disrupted genes. Overall, a total of 30 unique RefSeq genes [Supplementary Table 3 [Additional file 3]] and 43 unique TADs were disrupted [Supplementary Table 4 [Additional file 4]]. Potential gene disruption, dysregulation, and/or point mutations were investigated in each case to reveal the underlying mechanisms of azoospermia/severe oligospermia.
Simple-BCA: case MI1
The 39-year-old patient suffered from severe oligospermia with a karyotype of 46,XY,t(3;19)(p21.3;q13.3). Mate-pair GS revealed a simple reciprocal balanced translocation without copy number changes at the breakpoint junctions. Mate-pair GS revised the affected bands on both chromosomes to seq[GRCh37] t(3;19)(p14.3q13.2). Validated by PCR and Sanger sequencing, the breakpoint junctions were found to directly truncate the gene ERC2 at intron 12. The expression of ERC2 has been reported to be heat-sensitive in rat spermatocytes and round spermatid. Expression of genes that are heat-sensitive may be crucial because spermatogenesis occurs at approximately 3°C lower than body temperature. Mate-pair GS revealed no additional clinically significant CNVs and AOH. In addition, high read-depth GS not only confirmed the findings from mate-pair GS but also indicated that no SNVs/InDels associated with azoospermia/oligospermia were detected. As ERC2 is likely haploinsufficient (gnomAD: pLoF=0.99), its disruption could be associated with the phenotype in the patient.
Two independent sets of BCA: case MI2
This patient has primary infertility and azoospermia. The patient has two independent balanced translocations 46,XY,t(5;9)(p13.3;p22),t(7;21)(p13;q22.1) reported by G-banded chromosome analysis. Mate-pair GS revealed two balanced translocations t(7;21)(p13;q22.11) and t(5;9)(p14.1;p23) with a 15.4-kb deletion seq[GRCh37] del(9)(p23) chr9:g.10081969_10097410del at one of the breakpoint junctions (intron 3 of gene PTPRD). Sanger sequencing enabled fine mapping of three rearrangements [Figure 1]. Two breakpoints had 12-bp and 15-bp sequence insertions, and one had a 3-bp microhomology [Supplementary Table 5 [Additional file 5]]. These features suggested that the structural rearrangements were formed by DNA replication-based mechanisms.
|Figure 1: Two reciprocal translocations of case MI2 with a karyotype of 46,XY,t(5;9)(p13.3;p22),t(7;21)(p13;q22.1). The upper panel shows the normal chromosomes 5, 7, 9, and 21, and the lower panel shows the respective derivative chromosomes. Bars in red, green, orange, and purple indicate chromosomal segments on chromosomes 5, 7, 9, and 21, respectively. Each dotted line indicates a breakpoint junction on each chromosome. A small white bar indicates a copy number loss of 15.4 kb involving the intron 3 of PTPRD in chromosome 9. Disruption of candidate genes associated with male infertility is shown above the chromosomes. PTPRD: protein tyrosine phosphatase receptor type D; chr: chromosome; der: derivative chromosome.|
Click here to view
The genes COA1 and PTPRD were disrupted. Among them, PTPRD has been reported as a candidate gene for nonobstructive azoospermia., PTPRD is also predicted to be haploinsufficient (pLoF=1; gnomAD). Although four TADs were disrupted involving a total of 41 genes, none of the genes located within the disrupted TADs was associated with azoospermia/severe oligospermia. Mate-pair GS and high read-depth GS revealed no additional SNVs/InDels, CNVs, and AOH associated with male infertility in this case. Thus, the disruption of PTPRD may potentially explain the infertility.
Three-way BCA: case MI3
The patient has severe oligospermia and karyotype results indicated a three-way balanced translocation 46,XY,t(4;11;6)(q22;q21;q16). However, mate-pair GS detected additional cryptic complexities underlying the translocations which resulted in an overall of seven rearrangements without CNVs detected at the breakpoint junctions [Figure 2]. Among six rearrangements fine-mapped by Sanger sequencing, the features of the breakpoint junctions included blunt ends in four rearrangements and 1 bp microhomology in the other two rearrangements, suggesting end-joining as the rearrangement formation mechanism.
|Figure 2: Three-way translocation of case MI3 with a karyotype of 46,XY,t(4;11;6)(q22;q21;q16). The upper panel shows the normal chromosomes, and the lower panel shows the derivative chromosomes. Seven rearrangements were detected by mate-pair genome sequencing, revealing more complex rearrangements than a three-way translocation. Bars in red, green, and orange indicate chromosomal segments on chromosomes 4, 6, and 11, respectively. Segments originating from chromosome 4 were found on both der(6) and der(11). Each dotted line indicates the breakpoint on each chromosome. Bars in white represent cryptic chromosomal segments on 4q21.21q21.23, while boxes in light green indicate a chromosomal segment of 6q14.2q16.2. Segments that have an inverse genomic orientation are shown by reversed numbers. Disruption of candidate genes associated with male infertility is shown above the chromosomes. CFAP299: cilia and flagella associated protein 299; chr: chromosome; der: derivative chromosome.|
Click here to view
Genes CEP162 and CFAP299 were disrupted by the breakpoint junctions. Among them, the role of CFAP299 in spermatogenesis has been implicated in animal studies. Messenger RNA of CFAP299 is predominantly expressed in mouse testis at increased levels from 2–8-week-old testes in the cytoplasm of spermatogonia and primary spermatocytes. Silencing of CFAP299 resulted in an increase in the number of apoptotic cells and arrested cells at the G2/M phase. As such, CFAP299 plays a potential role in spermatogenesis in regulating cell apoptosis. In addition, five TADs were disrupted, involving 16 genes. However, none of the genes were known to be associated with spermatogenesis. Mate-pair GS and high read-depth GS confirmed that there were no additional clinically significant SNVs/InDels, CNVs, and AOH associated with azoospermia/oligospermia detected in this case. Therefore, disruption of CFAP299 could be a candidate cause of the infertility.
Three-way BCA: case MI4
The azoospermic patient carried a three-way reciprocal translocation: 46,XY,t(8;12;10)(q24.1;p13;q22), as shown in [Figure 3]. Mate-pair GS revealed additional cryptic complexities, with a total of seven rearrangements. There were no CNVs involved at the breakpoint junctions. All breakpoints were precisely mapped with Sanger sequencing. The sequence features showed blunt-ends in four rearrangements and 1 bp microhomology in the other three rearrangements. The formation of these rearrangements was likely mediated by end-joining mechanism.
|Figure 3: Three-way translocation of case MI4 with a karyotype of 46,XY,t(8;12;10)(q24.1;p13;q22). (a) The upper panel shows the normal chromosomes, and the lower panel shows the derivative chromosomes. Seven breakpoint junctions were detected by mate-pair genome sequencing. Bars in red, green and orange indicate chromosomal segments in chromosomes 8, 10 and 12, respectively. Each dotted line indicates the breakpoint locus on each chromosome. Bars in white represent chromosomal segments on 12p12.3. Characters in a reverse direction indicate a reverse orientation of the chromosomal segments in derivative chromosomes. Disruption of candidate genes associated with male infertility is shown above the chromosomes. (b) Sequence resolved breakpoint junction of a cryptic rearrangement der(12) reveals blunt end as the sequence characteristics. (c) The 3D genome map (http://3dgenome.fsm.northwestern.edu/) visualizing the 3D interaction and topologically associating domains. The black triangle shows a topologically associating domain disrupted by a breakpoint (shown by the vertical black line). Genes encompassed under each topologically associating domain are listed in the lower panel. EPS8: epidermal growth factor receptor pathway substrate 8; DHSs: DNase I hypersensitive sites; chr: chromosome; der: derivative chromosome.|
Click here to view
Three genes including EPS8, SOX5, and CTNNA3 were disrupted. EPS8 is known to be sensitive to haploinsufficiency (gnomAD: pLoF=0.97) and has been implicated to play a role in the induction of spermatid release from seminiferous epithelium during spermiation. Loss of EPS8 at the apical Sertoli cell–spermatid interface induced premature spermatid release and infertility. In addition, seven TADs were disrupted by the breakpoint junctions, which contained a total of 49 genes. Among them, PDE3A is a candidate gene associated with hypogonadism and Kallmann syndrome. Translocation breakpoints near the PDE3A gene have been reported in translocation carriers with Kallmann syndrome, potentially due to positional effects. In addition, female mice with PDE3A homozygous deletions were infertile, and ovulated oocytes were arrested in the germinal stage and could not be fertilized. Mate-pair GS and high read-depth GS found no additional SNVs/InDels, CNVs, and AOH associated with male infertility in this case. The disruption of EPS8 and the potential dysregulation of PDE3A are the possible explanations for infertility.
Two-way insertion: case MI5
This patient suffered from severe oligoasthenospermia and has been described in our previous study. An apparently balanced chromosomal insertion was identified by karyotyping analysis: 46,XY,ins(6;2)(q23;p13p22). Mate-pair GS detected nine rearrangements and indicated that the acceptor chromosome was likely involved in a chromothripsis event. Among them, six rearrangements were fine mapped by Sanger sequencing, including 1–2 bp microhomology (n = 3), small insertion (1 bp and 9 bp; n = 2), and blunt end (n = 1).
Seven genes were directly truncated by the breakpoints; however, none of which have been implicated in azoospermia/severe oligospermia. Eight TADs involving a total of 65 genes were disrupted. Although no direct evidence has shown associations of SLC17A5 (directly disrupted by the breakpoint junction) with male infertility, the gene EEF1A1 whose super enhancer was highly correlated to SLC17A5 was haploinsufficient (gnomAD: pLoF=0.98) and essential for spermatogenesis. Therefore, male infertility might be explained by the dysregulation of EEF1A1 expression through disruption of SLC17A5, leading to spermatogenesis failure. In addition, mutations in MCM9 are known to result in ovarian dysgenesis. MCM9-knockout male mice were sterile, resulted by arrest of spermatocytes in meiotic prophase I, while MCM9-knockout female mice were also sterile as their ovaries only contained arrested primary follicles and frequently developed tumors. There was insufficient DNA to pursue high read-depth GS in this case. Dysregulation of gene EEF1A1 and/or MCM9 by the structural rearrangement might explain the infertility in this case.
Complex insertion and translocation: case MI6
This patient had a complex karyotype: 46,XY,t(4;20)(q28;q12),der(20)ins(20;3)(q12;q28q13.3)inv(3)?(q13.3q25.3) involving three chromosomes. Mate-pair GS detected multiple complex rearrangements with 19 rearrangements [Figure 4]. Four copy number losses were found on chromosome 3, and the largest was a 1.15-Mb deletion involving the gene CFAP91. Mutations in CFAP91 are known to cause spermatogenic failure 51 (OMIM#619177) in an autosomal recessive manner (OMIM*609910). However, no additional variants were detected in CFAP91. Among the 19 rearrangements, 12 were successfully fine mapped by Sanger sequencing. The results revealed blunt ends at five junctions, 1–2 bp microhomology at six junctions, and 1 bp insertion at one junction. The breakpoint junction features along with frequent rearrangements suggested that this complex rearrangement was likely formed by DNA replicative mechanisms [Supplementary Table 5].
|Figure 4: Complex insertion and translocation of case MI6 detected by mate-pair genome sequencing. Bars in red, green, and orange indicate chromosomal segments in chromosomes 3, 4, and 20, respectively. Each dotted line indicates the breakpoint locus on each chromosome. Bars in white present cryptic chromosomal segment in 3q13.32q28 indicating there are 22 segments after chromosome scatting. Characters in a reverse direction indicate a reverse orientation of the chromosomal segments in derivative chromosomes. Disruption of candidate genes associated with male infertility is shown above the chromosomes. Copy number losses of four chromosome segments originating on chromosome 3 are indicated in a dotted frame. TP63: tumor protein p63; chr: chromosome; der: derivative chromosome.|
Click here to view
Overall, 14 genes were directly truncated by the rearrangements including TP63. TP63 is sensitive to haploinsufficiency (gnomAD: pLoF=1) and plays a role in mediating male germ cell development, apoptosis, and regulation of spermatogenesis., In addition, 15 TADs, involving a total of 128 genes, were disrupted. Among them, three genes were implicated in male infertility, including GMNC, SOX2, and CCDC39. GMNC (also known as Gemc1 in mice) plays a role in later stages of spermatogenesis. Gemc1−/− mice had a significantly reduced number of round spermatids and elongating spermatids. SOX2 is also a candidate gene for hypogonadotropic hypogonadism and resultant testicular failure. At last, mutations in CCDC39 are known to cause primary ciliary dyskinesia and Kartagener syndrome in an autosomal recessive manner, of which the common features include situs inversus, infertility, and oligoasthenospermia.
There were no additional SNVs/InDels in CFAP91, which was involved in a heterozygous deletion as described above. However, high read-depth GS detected two heterozygous SNVs in the gene DNAH1. Mutations in DNAH1 cause spermatogenic failure through multiple morphological abnormalities of the sperm flagella and primary ciliary dyskinesia in an autosomal recessive manner., However, the inheritance pattern of these two SNVs (in cis- or trans-) could not be determined as parental samples were not available. Nonetheless, male infertility in this case may be contributed by dual diagnoses of a complex rearrangement and point mutations.
| Discussion|| |
In this study, we investigated the genetic etiology of six cases with male infertility and apparently balanced interchromosomal structural rearrangements by mate-pair GS in combination with high read-depth GS. Overall, our study revealed that disruption or potential dysregulation of candidate genes implicated in male germ cell development, apoptosis, and spermatogenesis could be the underlying mechanisms of male infertility in these cases.
Balanced chromosomal structural rearrangements are thought to cause male infertility through failure of pairing of homologous elements. Recently, studies have shown male infertility in patients with BCAs may be caused by gene disruption or dysregulation by breakpoints or other genomic variants, for example, SNVs and InDels., While standard high read-depth GS cannot detect chromosomal rearrangements in low complexity or highly repetitive regions of the genome,, our large-insert size DNA mate-pair GS approach had increased sensitivity in detecting structural rearrangements. We identified the composition of each derivative chromosome in each case and revealed additional rearrangements cryptic to conventional G-banded chromosome analysis. We showed direct disruption or potential dysregulation of candidate genes associated with male infertility. In addition, high read-depth GS not only reported a dual molecular findings in case MI6, including a complex rearrangement and two point mutations in an autosomal recessive gene DNAH1, but also excluded known causative variants and noncoding regions (such as deep intronic region).
Our study on interchromosomal structural rearrangements revealed the complexities of BCAs that are largely underappreciated by G-banded chromosome analysis. Except MI1, additional cryptic rearrangements and/or copy number changes were detected in the remaining cases. The simplest was a cryptic deletion (15.4 kb) in case MI2 with two independent reciprocal translocations. In comparison, the rearrangements were most complex in MI6, which had 19 rearrangements and four deletions. The proportion of cases (5/6, 83.3%) with additional cryptic complex rearrangement identified by mate-pair GS over karyotyping was significantly higher than that reported in our previous study of males carrying BCAs in recurrent miscarriage couples (5/33, 15.2%, Fisher's exact test: P = 0.0023), who were able to conceive. It indicated that the complexities of structural rearrangements may be correlated with the phenotypic presentation (azoospermia or severe oligospermia).
Interestingly, among these cases, direct gene disruptions and/or potential dysregulation due to TAD disruptions are the potential underlying disease-causing mechanisms. In case MI2, the cryptic deletion (15.4 kb) was located in the intron of gene PTPRD. In case MI6, there were four copy number losses including gene CFAP91; a gene known to cause spermatogenic failure 51 (OMIM#619177) in an autosomal recessive manner (OMIM*609910). We excluded the possibility of compound heterozygous disease-causing alleles in autosomal recessive genes by identifying any causative hemizygous point mutations using high read-depth GS.
The limitations of our study include (1) limited sample size as there were six cases enrolled in this study, partly due to the rarity of apparently balanced chromosomal structural rearrangements and negative AZF region deletions in patients with severe oligospermia or azoospermia (0.5%–1.0%); (2) no parental samples available to study the inheritance mode of the detected rearrangements variants; for the two point mutations identified in DNAH1, long-read sequencing may resolve the haplotype spanning the two variants that are approximately 33 kb apart; however, no additional sample was available for this experiment; and (3) no additional samples were available for transcriptome analysis to confirm the potential consequences of the genomic variants. Nonetheless, our study provided a comprehensive composition of genomic structural rearrangements in each case and showed that the disruption or potential dysregulation of candidate genes could be the underlying causes of male infertility.
Overall, our study provided molecular characteristics of apparently balanced interchromosomal structural rearrangements in patients with male infertility. We showed the complexity of chromosomal structural rearrangements involved in male infertility and revealed that disruption or potential dysregulation of genes potentially implicated in male gametes development, apoptosis, and spermatogenesis could be a cause for infertility.
| Author Contributions|| |
MHKC, YL, KWC, XK, and ZD designed the study. JPWC collected the samples and followed up the clinical outcomes. MHKC, YL, MS, YKK, and ZD performed the analysis and data interpretation. MHKC, YL, PD and XZ conducted the validation. MHKC, YL, KWC, and ZD wrote the manuscript. All authors read and approved the final manuscript.
| Competing Interests|| |
All authors declare no competing interests.
| Acknowledgments|| |
This project is supported by the National Natural Science Foundation of China (No. 31801042), the Health and Medical Research Fund (No. 04152666 and No. 07180576), General Research Fund (No. 14115418), and Direct Grant (No. 2020.052).
Supplementary Information is linked to the online version of the paper on the Asian Journal of Andrology website.
| References|| |
Zegers-Hochschild F, Adamson GD, de Mouzon J, Ishihara O, Mansour R, et al
. International Committee for Monitoring Assisted Reproductive Technology (ICMART) and the World Health Organization (WHO) revised glossary of ART terminology, 2009. Fertil Steril
2009; 92: 1520–4.
Pfeifer S, Butts S, Dumesic D, Fossum G, Gracia C, et al
. Diagnostic evaluation of the infertile male: a committee opinion. Fertil Steril
2015; 103: e18–25.
Choy JT, Eisenberg ML. Male infertility as a window to health. Fertil Steri
2018; 110: 810–4.
Krausz C, Consultant S. Endocrinology & metabolism male infertility: pathogenesis and clinical diagnosis. Best Pract Res Clin Endocrinol Metab
2011; 25: 271–85.
Krausz C, Riera-Escamilla A. Genetics of male infertility. Nat Rev Urol
2018; 15: 369–84.
Kasak L, Laan M. Monogenic causes of non-obstructive azoospermia: challenges, established knowledge, limitations and perspectives. Hum Genet
2021; 140: 135–54.
Hamada A, Esteves S, Agarwal A. A comprehensive review of genetics and genetic testing in azoospermia. Clinics
2013; 68: 39–60.
Neto FT, Bach PV, Najari BB, Li PS, Goldstein M. Spermatogenesis in humans and its affecting factors. Semin Cell Dev Biol
2016; 59: 10–26.
Paoloni-Giacobino A, Kern I, Rumpler Y, Djlelati R, Morris M, et al
. Familial t(6;21)(p21.1;p13) translocation associated with male-only sterility. Clin Genet
2000; 58: 324–8.
Pastuszek E, Kiewisz J, Kulwikowska PM, Lukaszuk M, Lukaszuk K. Sperm parameters and DNA fragmentation of balanced chromosomal rearrangements carriers. Folia Histochem Cytobiol
2016; 53: 314–21.
Schilit SL, Menon S, Friedrich C, Kammin T, Wilch E, et al
. SYCP2 translocation-mediated dysregulation and frameshift variants cause human male infertility. Am J Hum Genet
2020; 106: 41–57.
Yammine T, Reynaud N, Lejeune H, Diguet F, Rollat-Farnier P, et al
. Deciphering balanced translocations in infertile males by next generation sequencing to identify candidate genes for spermatogenesis disorders. Mol Hum Reprod
2021; 27: gaab034.
Nguyen MH, Morel F, Pennamen P, Parent P, Douet-Guilbert N, et al
. Balanced complex chromosome rearrangement in male infertility: case report and literature review. Andrologia
2015; 47: 178–85.
Redin C, Brand H, Collins RL, Kammin T, Mitchell E, et al
. The genomic landscape of balanced cytogenetic abnormalities associated with human congenital anomalies. Nat Genet
2017; 49: 36–45.
David D, Freixo JP, Fino J, Carvalho I, Marques M, et al
. Comprehensive clinically oriented workflow for nucleotide level resolution and interpretation in prenatal diagnosis of de novo
apparently balanced chromosomal translocations in their genomic landscape. Hum Genet
2020; 139: 531–43.
Dong Z, Yan J, Xu F, Yuan J, Jiang H, et al
. Genome sequencing explores complexity of chromosomal abnormalities in recurrent miscarriage. Am J Hum Genet
2019; 105: 1102–11.
Dong Z, Chau MH, Zhang Y, Dai P, Zhu X, et al
. Deciphering the complexity of simple chromosomal insertions by genome sequencing. Hum Genet
2021; 140: 361–80.
Dong Z, Zhao X, Li Q, Yang Z, Xi Y, et al
. Development of coupling controlled polymerizations by adapter-ligation in mate-pair sequencing for detection of various genomic variants in one single assay. DNA Res
2019; 26: 313–25.
Dong Z, Wang H, Chen H, Jiang H, Yuan J, et al
. Identification of balanced chromosomal rearrangements previously unknown among participants in the 1000 Genomes Project: implications for interpretation of structural variation in genomes and the future of clinical cytogenetics. Genet Med
2018; 20: 697–707.
Dong Z, Jiang L, Yang C, Hu H, Wang X, et al
. A robust approach for blind detection of balanced chromosomal rearrangements with whole-genome low-coverage sequencing. Hum Mutat
2014; 35: 625–36.
Dong Z, Zhang J, Hu P, Chen H, Xu J, et al
. Low-pass whole-genome sequencing in clinical cytogenetics: a validated approach. Genet Med
2016; 18: 940–8.
Wang H, Dong Z, Zhang R, Chau MH, Yang Z, et al
. Low-pass genome sequencing versus chromosomal microarray analysis: implementation in prenatal diagnosis. Genet Med
2020; 22: 500–10.
Chau MH, Wang H, Lai Y, Zhang Y, Xu F, et al
. Low-pass genome sequencing: a validated method in clinical cytogenetics. Hum Genet
2020; 139: 1403–15.
Dong Z, Chau MH, Zhang Y, Yang Z, Shi M, et al
. Low-pass genome sequencing-based detection of absence of heterozygosity: validation in clinical cytogenetics. Genet Med
2021; 23: 1225–33.
McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, et al
. The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res
2010; 20: 1297–303.
Wang K, Li M, Hakonarson H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res
2010; 38: e164.
Li Q, Wang K. InterVar: clinical interpretation of genetic variants by the 2015 ACMG-AMP Guidelines. Am J Hum Genet
2017; 100: 267–80.
Richards S, Aziz N, Bale S, Bick D, Das S, et al
. Standards and guidelines for the interpretation of sequence variants: a joint consensus recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology. Genet Med
2015; 17: 405–24.
Nalbandian K, Piña-Aguilar RE, Morton CC. Resolving breakpoints of chromosomal rearrangements at the nucleotide level using Sanger sequencing. Curr Protoc Hum Genet
2020; 108: e107.
Dixon JR, Selvaraj S, Yue F, Kim A, Li Y, et al
. Topological domains in mammalian genomes identified by analysis of chromatin interactions. Nature
2012; 485: 376–80.
Leung TY, Vogel I, Lau TK, Chong W, Hyett JA, et al
. Identification of submicroscopic chromosomal aberrations in fetuses with increased nuchal translucency and apparently normal karyotype. Ultrasound Obstet Gynecol
2011; 38: 314–9.
Huang J, Poon LC, Akolekar R, Choy KW, Leung TY, et al
. Is high fetal nuchal translucency associated with submicroscopic chromosomal abnormalities on array CGH? Ultrasound Obstet Gynecol
2014; 43: 620–4.
Yadav SK, Pandey A, Kumar L, Devi A, Kushwaha B, et al
. The thermo-sensitive gene expression signatures of spermatogenesis. Reprod Biol Endocrinol
2018; 16: 56.
Yu J, Wu H, Wen Y, Liu Y, Zhou T, et al
. Identification of seven genes essential for male fertility through a genome-wide association study of non-obstructive azoospermia and RNA interference-mediated large-scale functional screening in Drosophila
. Hum Mol Genet
2015; 24: 1493–503.
Barban N, Jansen R, de Vlaming R, Vaez A, Mandemakers JJ, et al
. Genome-wide analysis identifies 12 loci influencing human reproductive behavior. Nat Genet
2016; 48: 1462–72.
Li H, Dai Y, Luo Z, Nie D. Cloning of a new testis-enriched gene C4orf22 and its role in cell cycle and apoptosis in mouse spermatogenic cells. Mol Biol Rep
2019; 46: 2029–38.
Cheng CY, Mruk DD. Regulation of spermiogenesis, spermiation and blood-testis barrier dynamics: novel insights from studies on Eps8 and Arp3. Biochem J
2011; 435: 553–62.
Quaynor SD, Bosley ME, Duckworth CG, Porter KR, Kim S, et al
. Molecular and cellular endocrinology targeted next generation sequencing approach identifies eighteen new candidate genes in normosmic hypogonadotropic hypogonadism and Kallmann syndrome. Mol Cell Endocrinol
2016; 437: 86–96.
Bhagavath B, Podolsky RH, Ozata M, Bolu E, Bick DP, et al
. Clinical and molecular characterization of a large sample of patients with hypogonadotropic hypogonadism. Fertil Steril
2006; 85: 706–13.
Masciarelli S, Horner K, Liu C, Park SH, Hinckley M, et al
. Cyclic nucleotide phosphodiesterase 3A-deficient mice as a model of female infertility. J Clin Invest
2004; 114: 196–205.
Chen J, Jiang D, Tan D, Fan Z, Wei Y, et al
. Heterozygous mutation of eEF1A1b
resulted in spermatogenesis arrest and infertility in male tilapia, Oreochromis niloticus
. Sci Rep
2017; 7: 43733.
Lutzmann M, Grey C, Traver S, Ganier O, Maya-Mendoza A, et al
. MCM8- and MCM9-deficient mice reveal gametogenesis defects and genome instability due to impaired homologous recombination. Mol Cell
2012; 47: 523–34.
Wang H, Yuan Q, Niu M, Zhang W, Wen L, et al
. Transcriptional regulation of P63 on the apoptosis of male germ cells and three stages of spermatogenesis in mice. Cell Death Dis
2018; 9: 76.
Petre-Lazar B, Livera G, Moreno SG, Trautmann E, Duquenne C, et al
. The role of p63 in germ cell apoptosis in the developing testis. J Cell Physiol
2007; 210: 87–98.
Terré B, Lewis M, Gil-Gómez G, Han Z, Hao L, et al
. Defects in efferent duct multiciliogenesis underlie male infertility in GEMC1, MCIDAS or CCNO deficient mice. Development
2019; 146: dev162628.
Song SH, Chiba K, Ramasamy R, Lamb D. Recent advances in the genetics of testicular failure. Asian J Androl
2016; 18: 350–5.
Merveille AC, Davis EE, Becker-Heck A, Legendre M, Amirav I, et al
. CCDC39 is required for assembly of inner dynein arms and the dynein regulatory complex and for normal ciliary motility in humans and dogs. Nat Genet
2011; 43: 72–8.
Amiri-Yekta A, Coutton C, Kherraf ZE, Karaouzène T, Le Tanno P, et al
. Whole-exome sequencing of familial cases of multiple morphological abnormalities of the sperm flagella (MMAF) reveals new DNAH1
mutations. Hum Reprod
2016; 31: 2872–80.
Ben Khelifa M, Coutton C, Zouari R, Karaouzène T, Rendu J, et al
. Mutations in DNAH1
, which encodes an inner arm heavy chain dynein, lead to male infertility from multiple morphological abnormalities of the sperm flagella. Am J Hum Genet
2014; 94: 95–104.
Mahmoud M, Gobet N, Cruz-Dávalos DI, Mounier N, Dessimoz C, et al
. Structural variant calling: the long and the short of it. Genome Biol
2019; 20: 246.
Amarasinghe SL, Su S, Dong X, Zappia L, Ritchie ME, et al
. Opportunities and challenges in long-read sequencing data analysis. Genome Biol
2020; 21: 30.
[Figure 1], [Figure 2], [Figure 3], [Figure 4]