Genotype-phenotype correlation and identification of two novel SRD5A2 mutations in 33 Chinese patients with hypospadias

Novel Genomic Variants, Atypical Phenotypes and Evidence of a Digenic/Oligogenic Contribution to Disorders/Differences of Sex Development in a Large North African Cohort

In a majority of individuals with disorders/differences of sex development (DSD) a genetic etiology is often elusive. However, new genes causing DSD are routinely reported and using the unbiased genomic approaches, such as whole exome sequencing (WES) should result in an increased diagnostic yield. Here, we performed WES on a large cohort of 125 individuals all of Algerian origin, who presented with a wide range of DSD phenotypes. The study excluded individuals with congenital adrenal hypoplasia (CAH) or chromosomal DSD. Parental consanguinity was reported in 36% of individuals. The genetic etiology was established in 49.6% (62/125) individuals of the total cohort, which includes 42.2% (35/83) of 46, XY non-syndromic DSD and 69.2% (27/39) of 46, XY syndromic DSD. No pathogenic variants were identified in the 46, XX DSD cases (0/3). Variants in the AR, HSD17B3, NR5A1 and SRD5A2 genes were the most common causes of DSD. Other variants were identified in genes associated with congenital hypogonadotropic hypogonadism (CHH), including the CHD7 and PROKR2. Previously unreported pathogenic/likely pathogenic variants (n = 30) involving 25 different genes were identified in 22.4% of the cohort. Remarkably 11.5% of the 46, XY DSD group carried variants classified as pathogenic/likely pathogenic variant in more than one gene known to cause DSD. The data indicates that variants in PLXNA3, a candidate CHH gene, is unlikely to be involved in CHH. The data also suggest that NR2F2 variants may cause 46, XY DSD.

Whole exome sequencing applied to 42 Han Chinese patients with posterior hypospadias

Hypospadias, a malformation of male external genitalia, is characterized by an aberrant opening of the urethra on the ventral side of the penis. It is considered a complex disorder with both environmental and genetic factors involved in its pathogenesis. To identify the genetic abnormality involved in the pathogenesis of hypospadias, we performed whole exome sequencing (WES) analysis in 42 hypospadias patients with karyotype 46, XY in the Nanhai Meternity&Child Health Hospital of Foshan. All the likely pathogenic variants were confirmed by Sanger sequencing and assessed by Sorting Intolerant from Tolerant (SIFT), PROVEAN, PolyPhen2, ClinPred, LRT, Mutation Assessor, FATHMM, and GERP software. We discovered 27 gene mutations in 20 patients, including eight cases of the SRD5A2 gene, 4 cases of the AR gene, 3 cases of the CYP17A1 gene, 1 case of the WT1 gene, 1 case of the ANOS1 gene, 1 case of the NR5A1 gene, 1 case of the FGFR1 gene, and one case of the DHX37 gene. Our study is the first to describe six novel missense mutations, AR(c.302G > A, c.2593G > T, and c.1705G > T), CYP17A1(c.1298 T > C), FGFR1 (c.995C > T) and DHX37(c.923G > A). In summary, genetic defect detection was useful for early diagnosis of severe hypospadias in the Han Chinese population. Nevertheless, most cases remain unexplained, and the exact pathogenesis of hypospadias still needs further study.

Etiology of Hypospadias: A Comparative Review of Genetic Factors and Developmental Processes Between Human and Animal Models

Hypospadias is a congenital anomaly of the penis with an occurrence of approximately 1 in 200 boys, but the etiology of the majority of hypospadias has remained unknown. Numerous genes have been reported as having variants in hypospadias patients, and many studies on genetic deletion of key genes in mouse genital development have also been published. Until now, no comparative analysis in the genes related literature has been reported. The basic knowledge of penile development and hypospadias is mainly obtained from animal model studies. Understanding of the differences and similarities between human and animal models is crucial for studies of hypospadias. In this review, mutations and polymorphisms of hypospadias-related genes have been compared between humans and mice, and differential genotype–phenotype relationships of certain genes between humans and mice have been discussed using the data available in PubMed and MGI online databases, and our analysis only revealed mutations in seven out of 43 human hypospadias related genes which have been reported to show similar phenotypes in mutant mice. The differences and similarities in the processes of penile development and hypospadias malformation among human and commonly used animal models suggest that the guinea pig may be a good model to study the mechanism of human penile development and etiology of hypospadias.

Clinical and molecular characterization of 12 prenatal cases of Cri-du-chat syndrome

Background

This study aimed to define the molecular basis for 12 prenatal cases of Cri‐du‐chat syndrome (CdCS) and the potential genotyping‐phenotyping association.

Methods

Karyotyping and single nucleotide polymorphism array analyses for copy number variants were performed.

Results

Nine cases had 5p terminal deletions and three had 5p interstitial deletions, and these cases had variable deletion sizes with partial overlapping. Phenotypically, besides intrauterine growth restriction (IUGR) and brain as well as heart abnormalities, hypospadias, and lung dysplasia were observed. Potential genetic causes for specific phenotypes in these cases were identified.

Conclusion

This study defined the molecular bases for the patients of CdCS, which is important for genetic counseling for these families. The findings of present study expand the clinical features of CdCS in the fetal period, and provided important information for further refining the genotypic–phenotypic correlations for this syndrome.

Mutational analysis of SRD5A2: From gene to functional kinetics in individuals with steroid 5α-reductase 2 deficiency

Human steroid 5α-reductase 2 (SRD5A2) plays a determinative role in the masculinization of external genitalia. To date, approximately 114 different mutations of the SRD5A2 gene have been reported; however, little information is available about their impact on catalytic function or their three-dimensional (3D) structures. We determined the effect of point mutations on the testosterone-depend kinetic constants (K_m,app and V_max,app) and structural characteristics of SRD5A2 from Mexican patients with 46,XY-steroid 5α-reductase 2 deficiency. PCR-SSCP assays identified ten distinct gene variants and sequencing analysis identified missense mutations [p.V3I, p.S14R, p.A52T, p.F118L, p.R145W, p.R171S, p.L226P, p.F229S, p.S245Y, and p.A248V]. Mutations were re-created by site-directed mutagenesis and expressed in HEK293 cells. Functional studies demonstrated that 8 variants led to partial (K_m,app = 0.16-2.6 μM; V_max,app = 224-2640 pmol/mg P/min) or complete losses of activity compared to the wild-type enzyme (K_m,app = 0.7 μM; V_max,app = 4044 pmol/mg P/min). All the mutations were assessed using multiple software tools and the results predicted that all of the mutations were associated with disease or damage. Mapping mutations on the model of a 3D structure of SRD5A2 demonstrated alterations in contact sites with their proximal amino acids. Our data show that mutations affect the catalytic efficiency (V_max/K_m) or result in residual enzymatic activity, which could be due to erroneous interactions between amino acid residues, the substrate testosterone, or NADPH.

New frontiers on the molecular underpinnings of hypospadias according to severity

Hypospadias, which is characterised by the displacement of the urethral meatus from its typical anatomical location in males, shows various degrees of severity. In this systematic review, we surveyed our current understanding of the genetics of isolated hypospadias in humans according to the severity of the condition. We found that sequencing and genotyping approaches were the preferred methods of study and that single nucleotide polymorphisms were the most common finding associated with hypospadias. Most genes fell into four gene-pathway categories related to androgens, oestrogens, growth factors, or transcription factors. Few hypospadias studies classify their findings by severity. Taken together, we argue that it is advantageous to take into consideration the severity of the condition in search of novel candidates in the aetiology of hypospadias. Abbreviations: AR: androgen receptor; ATF3: activating transcription factor 3; BMP4: bone morphogenetic protein 4; BMP7: bone morphogenetic protein 7; CYP17: steroid 17-alpha-hydroxylase/17,20 lyase; CYP1A1: cytochrome P450 1A1; CYP3A4: cytochrome P450 3A4; CNVs: copy number variants; DGKK: diacylglycerol kinase kappa; ESR1: oestrogen receptor 1; ESR2: oestrogen receptor 2; FGF8: fibroblast growth factor 8; FGF10: fibroblast growth factor 10; FGFR2: fibroblast growth factor receptor 2; HOXA4: homeobox protein Hox-A4; HOXB6: homeobox protein Hox-B6; HSD17B3: hydroxysteroid 17-beta dehydrogenase 3; MAMLD1: mastermind-like domain-containing protein 1; SF-1: splicing factor 1; SHH: sonic hedgehog; SNPs: single nucleotide polymorphisms; SOX9: SRY-box 9; SRD5A2: steroid 5 alpha-reductase 2; SRY: sex-determining region Y protein; STAR: steroidogenic acute regulatory protein; STARD3: StAR-related lipid transfer protein 3; STS: steryl-sulfatase; WT1: Wilms tumour protein; ZEB1: zinc finger oestrogen-box binding homeobox 1.

Integrative and Analytical Review of the 5-Alpha-Reductase Type 2 Deficiency Worldwide

Introduction

The conversion of testosterone into dihydrotestosterone is catalyzed by the 5α-reductase type 2 enzyme which plays a crucial role in the external genitalia virilization. It is encoded by the SRD5A2 gene. Allelic variants in this gene cause a 46,XY DSD with no genotype-phenotype relationship. It was firstly reported in the early 70s from isolated clusters. Since then, several cases have been reported. Putting together, it will expand the knowledge on the molecular bases of androgen milieu.

Methods

We searched for SRD5A2 allelic variants (AV) in the literature (PubMed, Embase, MEDLINE) and websites (ensembl, HGMD, ClinVar). Only cases with AV in both alleles, either in homozygous or compound heterozygous were included. The included cases were analyzed according to ethnicity, exon, domain, aminoacid (aa) conservation, age at diagnosis, sex assignment, gender reassignment, external genitalia virilization and functional studies. External genitalia virilization was scored using Sinnecker scale. Conservation analysis was carried out using the CONSURF platform. For categorical variables, we used X2 test and Cramer's V. Continuous variables were analyzed by t test or ANOVA. Concordance was estimated by Kappa.

Results

We identified 434 cases of 5ARD2 deficiencies from 44 countries. Most came from Turkey (23%), China (17%), Italy (9%), and Brazil (7%). Sixty-nine percent were assigned as female. There were 70% of homozygous allelic variants and 30% compound heterozygous. Most were missense variants (76%). However, small indels (11%), splicing (5%) and large deletions (4%) were all reported. They were distributed along with all exons with exon 1 (33%) and exon 4 (25%) predominance. Allelic variants in the exon 4 (NADPH-binding domain) resulted in lower virilization (p<0.0001). The codons 55, 65, 196, 235 and 246 are hotspots making up 25% of all allelic variants. Most of them (76%) were located at conserved aa. However, allelic variants at non-conserved aa were more frequently indels (28% vs 6%; p<0.01). The overall rate of gender change from female to male ranged from 16% to 70%. The lowest rate of gender change from female to male occurred in Turkey and the highest in Brazil. External genitalia virilization was similar between those who changed and those who kept their assigned gender. The gender change rate was significantly different across the countries (V=0.44; p<0.001) even with similar virilization scores.

Conclusion

5ARD2 deficiency has a worldwide distribution. Allelic variants at the NADPH-ligand region cause lower virilization. Genitalia virilization influenced sex assignment but not gender change which was influenced by cultural aspects across the countries. Molecular diagnosis influenced on sex assignment, favoring male sex assignment in newborns with 5α-reductase type 2 deficiency.

Phenotypic and molecular characteristics of androgen insensitivity syndrome patients

Androgen insensitivity syndrome (AIS), an X-linked recessive genetic disorder of sex development, is caused by mutations in the androgen receptor (AR) gene, and is characterized by partial or complete inability of specific tissues to respond to androgens in individuals with the 46, XY karyotype. This study aimed to investigate AR gene mutations and to characterize genotype–phenotype correlations. Ten patients from unrelated families, aged 2–31 years, were recruited in the study. Based on karyotype, altered hormone profile, and clinical manifestations, nine patients were preliminarily diagnosed with complete AIS and one with partial AIS. Genetic analysis of AR gene revealed the existence of 10 different mutations, of which five were novel (c.2112 C>G[p.S704R], c.2290T>A[p.Y764N], c.2626C>T[p.Q876X], c.933dupC[p.K313Qfs*28], and c.1067delC[p.A356Efs*123]); the other five were previously reported (c.1789G>A[p.A597T], c.2566C>T[p.R856C], c.2668G>A[p.V890M], c.2679C>T[p.P893L], and c.1605C>G[p.Y535X]). Regarding the distribution of these mutations, 60.0% were clustered in the ligand-binding domain of AR gene. Exons 1 and 8 of AR gene each accounted for 30.0% (3/10) of all mutations. Most of the truncation mutations were in exon 1 and missense mutations were mainly located in exons 4–8. Our study expands the spectrum of AR gene mutations and confirms the usefulness of AR gene sequencing to support a diagnosis of AIS and to enable prenatal or antenatal screening.

Identification of gene variants in 130 Han Chinese patients with hypospadias by targeted next-generation sequencing

Background

Hypospadias is a common congenital malformation of male external genitalia, which mainly manifests as an abnormal urethral opening on the ventral side of the penis. The etiology and clinical phenotype of hypospadias is highly heterogeneous, and its clinical diagnosis is challenging. Currently, over 70% of patients have an unknown etiology. Here, we performed a targeted analysis of gene mutations in 130 patients with hypospadias of unknown etiology to find the precise genetic cause.

Methods

We developed a targeted next‐generation sequencing (NGS) panel, encompassing the exon coding regions of 105 genes involved in external genitalia and urogenital tract development and performed sequencing analysis on 130 children with hypospadias of unknown etiology.

Results

In total, 25 patients with hypospadias (19.2%) were found to have 20 mutations among the nine genes involved in external genitalia and urogenital tract development, including 16 reported and four novel mutation sites. Twenty‐two patients (16.9%) had diagnostic variants. Multiple genetic mutations were identified in three of the 25 patients. Hypospadias combined with micropenis was the most common phenotype (68%) in 25 patients.

Conclusions

Higher frequency mutations were identified in SRD5A2 (52%) and AR (24%) in our patient cohort. Middle or posterior hypospadias with micropenis may be significant indicators of genetic variations. Polygenic inheritance may be a rare genetic cause of hypospadias.

Polymorphisms of MAMLD1, SRD5A2, and AR Candidate Genes in Seven Dogs (78,XY; SRY-Positive) Affected by Hypospadias or Cryptorchidism

Knowledge of the molecular background of disorders of sex development (DSD) in dogs with normal sets of XY chromosomes (XY DSD) is very scarce. However, extensive studies have been carried out in humans, showing that polymorphisms and mutations of numerous genes, including SRY, MAMLD1, SRD5A2, and AR, are associated with or responsible for XY DSD. In this study, we analyzed the entire coding sequence of these genes in 7 dogs (78,XY) with ambiguous external genitalia (hypospadias, cryptorchidism, bifid scrotum, or rudimentary penis). The most common disorder was hypospadias (6 cases), followed by cryptorchidism (4 cases). The co-occurrence of both abnormalities was observed in 3 dogs. Polymorphisms were found in MAMLD1 (3 SNPs), SRD5A2 (5 SNPs), and AR (2 STRs and 1 SNP), while SRY was monomorphic. However, the distribution of the polymorphic variants in the DSD dogs and 11 control XY dogs did not differ significantly. Our study suggests that an association between the polymorphisms of the studied candidate genes and hypospadias or cryptorchidism is unlikely in dogs. We thus support the recent suggestion that hypospadias is not rare in this species, and moreover, we show that co-occurrence of hypospadias and cryptorchidism can be quite frequent.

Molecular diagnostics of disorders of sexual development: an Indian survey and systems biology perspective

We aimed to survey the monogenic causes of disorders of sex development (DSD) and thereby its prevalence in India. This study revealed mutations resulting in androgen insensitivity syndrome, 5α-reductase type 2 deficiency, and gonadal dysgenesis were commonly reported. Intriguingly, AR deficits were the most prevalent (32 mutations) and of 11/26 missense mutations were in exons 4-8 (encoding ligand binding domain). The unique features of SRD5A2 defects were p.R246Q (most prevalent) and p.G196S could be mutational hotspots, dual gene defects (p.A596T in AR and p.G196S in SRD5A2) in a patient with hypospadias and novel 8 nucleotide deletion (exon 1) found in a patient with perineal hypospadias. Deficits in SRY, WT1, DHH, NR5A1, and DMRT1 caused 46,XY gonadal dysgenesis. Notably, mutations in AR, SRD5A2, MAMLD1, WT1, and MAP3K1 have led to hypospadias and only one CYP19A1 mutation caused aromatase deficiency was reported to date. Data mining from various databases has not only reinforced the role of well-established genes (e.g., SRY, WT1, DHH, NR5A1, DMRT1, AR, SRD5A2, MAMLD1) involved in DSD but also provided us 12 more potential candidate genes (ACVR1, AMHR2, CTNNB1, CYP11A1, CYP19A1, FGFR2, FGF9, PRKACA, PRKACG, SMAD9, TERT, ZFPM2), which benefit from a close association with the well-established genes involved in DSD and might be useful to screen owing to their direct gene-phenotype relationship or through direct functional interaction. As more genes have been revealed in relation to DSD, we believe ultimately it holds a better scenario for therapeutic regimen. Despite the advances in translational medicine, hospitals are yet to adopt genetic testing and counseling facilities in India that shall have potential impact on clinical diagnosis. Abbreviations: 5α-RD2: 5α-Reductase type 2; AIS: androgen insensitivity syndrome; AMH: antimullerian hormone; AMHR: antimullerian hormone receptor; AR: androgen receptor gene; CAH: congenital adrenal hyperplasia; CAIS: complete AIS; CAH: congenital adrenal hyperplasia; CHH: congenital hypogonadotropic hypogonadism; CXORF6: chromosome X open reading frame 6 gene; CYP19A1: cytochrome P450 family 19 subfamily A member 1 gene; DHT: dihydrotestosterone; DMRT1: double sex and mab-3 related transcription factor 1 gene; DSD: disorders of sexual development; GD: gonadal dysgenesis; HGMD: human gene mutation database; IH: isolated hypospadias; MAMLD1: mastermind like domain containing 1 gene; MIS: mullerian inhibiting substance; NTD: N-terminal domain; OT DSD: ovotesticular DSD; PAIS: partial AIS; SOX9: SRY-related HMG-box 9 gene; SRY: sex-determining region Y gene; STAR: steroidogenic acute regulatory protein gene; SRD5A2: steroid 5 alpha-reductase 2 gene; T DSD: testicular DSD; T: testosterone; WNT4: Wnt family member 4 gene; WT1: Wilms tumor 1 gene; Δ₄: androstenedione.