Novel mutations of the SRD5A2 and AR genes in Thai patients with 46, XY disorders of sex development

Clinical, Hormonal, and Genetic Spectrum of 46 XY Disorders of Sexual Development (DSD) Patients

OBJECTIVES

To evaluate the clinical, hormonal and genetic characteristics of 46XY disorders of sexual development (DSD) patients from South India.

METHODS

46XY DSD patients with a provisional diagnosis of 17β-hydroxysteroid dehydrogenase 3 (17BHSD3) deficiency, 5 alpha-reductase type 2 deficiency (5ARD2) or partial androgen insensitivity syndrome (PAIS) based on clinical and hormonal analysis were included in this study. All the patients underwent detailed clinical and hormonal evaluations. Targeted next-generation sequencing for all three genes (AR, HSD17B3, and SRD5A2) in parallel was carried out for all the included patients and their parents.

RESULTS

Based upon the clinical and hormonal analysis, among the 37 children with 46XY DSD in the present study, 21 children were diagnosed with 5ARD2, 10 with PAIS, and six with 17BHSD3 deficiency. However, genetic analysis revealed pathogenic mutations in nine patients - six in the AR gene, two in the SRD5A2 gene, and one in the HSD17B3 gene. The concordance rate between provisional hormonal and genetic diagnosis was only 22.2%. Two out of six subjects with AR gene variants were positive for somatic mosaicism.

CONCLUSIONS

In the present study, a positive genetic diagnosis was detected in nine patients (24%), including five novel variants. In this study, mutations in the AR gene was the most reported. The authors did not find the testosterone: dihydrotestosterone (T: DHT) ratio to be an accurate hormonal diagnostic tool.

46,XY disorders of sex development: the use of NGS for prevalent variants

46,XY disorders of sex development (DSD) present with diverse phenotypes and complicated genetic causes. Precise genetic diagnosis contributes to accurate management, and targeted next-generation sequencing (NGS) and whole-exome sequencing are powerful tools for investigating DSD. However, the prevalent variants resulting in 46,XY DSD remain unclear, especially those associated with mild forms, such as isolated hypospadias, inguinal cryptorchidism, and micropenis. From 2019 to 2021, 74 patients with 46,XY DSD (48 typical and 26 mild) from the First Affiliated Hospital of Sun Yat-sen University were enrolled in our cohort study for targeted NGS or whole-exome sequencing. Our targeted 46,XY DSD panel included 108 genes involved in disorders of gonadal development and differentiation, steroid hormone synthesis and activation, persistent Müllerian duct syndrome, idiopathic hypogonadotropic hypogonadism, syndromic disorder, and others. Variants were classified as pathogenic, likely pathogenic, variant of uncertain significance, likely benign, or benign following the American College of Medical Genetics guidelines. As a result, 28 of 74 (37.8%) patients with pathogenic and/or likely pathogenic variants acquired genetic diagnoses. The Mild DSD patients acquired a diagnosis rate of 30.7%. We detected 44 variants in 28 DSD genes from 31 patients, including 33 novel and 11 reported variants. Heterozygous (65%) and missense (70.5%) variants were the most common. Variants associated with steroid hormone synthesis and activation were the main genetic causes of 46,XY DSD. In conclusion, 46,XY DSD manifests as a series of complicated polygenetic diseases. NGS reveals prevalent variants and improves the genetic diagnoses of 46,XY DSD, regardless of severity.

Clinical, Hormonal, and Genetic Characteristics of 5α-Reductase Type 2 Deficiency in 103 Chinese Patients

OBJECTIVE

5α-Reductase type 2 (5α-RD2) deficiency causes variable degrees of undervirilization in patients. The correlation between its genotype and phenotype is unclear.

METHODS

We retrospectively evaluated 103 patients with 46,XY disorders of sex development who were diagnosed with 5α-RD2 deficiency.

RESULTS

The prevalence of female sex assignment (P = .008) and the incidences of cryptorchidism (P = .0003) and bifid scrotum (P = .0002) in the non-p.R227Q variant group were higher, but there were no significant differences in the incidences of hypospadias and isolated microphallus. The external masculinization score in the non-p.R227Q variant group was lower than that in the homozygous p.R227Q variant (P = .019) and compound heterozygous p.R227Q variant groups (P = .013). The level of anti-Mullerian hormone in the non-p.R227Q variant group was lower than that in the homozygous p.R227Q variant (P < .001) and compound heterozygous p.R227Q variant groups (P = .006). The testosterone-to-dihydrotestosterone ratio of the homozygous p.R227Q variant group was higher than that of the non-p.R227Q variant (P = .018) and compound heterozygous p.R227Q variant groups (P = .029). Twenty-three reportedly pathogenic variants and 11 novel steroid 5α-reductase 2 (SRD5A2) variants were identified.

CONCLUSION

Compared with patients without p.R227Q, patients with p.R227Q exhibited higher external masculinization scores and anti-Mullerian hormone expression, a lower prevalence of female sex assignment, and lower incidences of cryptorchidism and bifid scrotum. We identified 23 reportedly pathogenic SRD5A2 variants and 11 novel SRD5A2 variants that led to 5α-RD2 deficiency. We established a genotype-phenotype correlation, and patients with p.R227Q showed a relatively mild phenotype.

DETECTION OF DISEASE‐CAUSING MUTATIONS IN PROSTATE CANCER BY NGS SEQUENCING

Gene mutations may affect the fate of many tumors including prostate cancer (PCa); therefore, the research of specific mutations associated with tumor outcomes might help the urologist to identify the best therapy for PCa patients such as surgical resection, adjuvant therapy or active surveillance. Genomic DNA (gDNA) was extracted from 48 paraffin‐embedded PCa samples and normal paired tissues. Next, gDNA was amplified and analyzed by next‐generation sequencing (NGS) using a specific gene panel for PCa. Raw data were refined to exclude false‐positive mutations; thus, variants with coverage and frequency lower than 100× and 5%, respectively were removed. Mutation significance was processed by Genomic Evolutionary Rate Profiling, ClinVar, and Varsome tools. Most of 3000 mutations (80%) were single nucleotide variants and the remaining 20% indels. After raw data elaboration, 312 variants were selected. Most mutated genes were KMT2D (26.45%), FOXA1 (16.13%), ATM (15.81%), ZFHX3 (9.35%), TP53 (8.06%), and APC (5.48%). Hot spot mutations in FOXA1, ATM, ZFHX3, SPOP, and MED12 were also found. Truncating mutations of ATM, lesions lying in hot spot regions of SPOP and FOXA1 as well as mutations of TP53 correlated with poor prognosis. Importantly, we have also found some germline mutations associated with hereditary cancer‐predisposing syndrome. gDNA sequencing of 48 cancer tissues by NGS allowed to detect new tumor variants as well as confirmed lesions in genes linked to prostate cancer. Overall, somatic and germline mutations linked to good/poor prognosis could represent new prognostic tools to improve the management of PCa patients.

Clinical, Biochemical, and Molecular Characterization of Indian Children with Clinically Suspected Androgen Insensitivity Syndrome

This study describes the clinical, biochemical, and molecular characteristics of Indian children with 46,XY DSD and suspected androgen insensitivity syndrome (AIS). Fifty children (median age 3.0 years, range 0-16.5 years) with 46,XY DSD and a suspected diagnosis of AIS were enrolled. Sanger sequencing was performed to identify pathogenic variants in the androgen receptor (AR) gene and to study genotype-phenotype correlations. All 5 (100%) patients with CAIS and 14/45 (31%) patients with PAIS had pathogenic/likely pathogenic variants in the AR gene (overall, 14 different variants in 19 patients; 38.8%). There was no significant difference in clinical (cryptorchidism, hypospadias, or external masculinizing score) or biochemical parameters (gonadotropins and testosterone) between patients with or without pathogenic variants. However, patients with AIS were more likely to have a positive family history, be assigned female gender at birth, and present with gynaecomastia at puberty. Three novel pathogenic/likely pathogenic variants, including one splice donor site variant c.2318+1G>A, one frameshift variant p.H790Lfs*40, and one missense variant p.G821E, were identified in 3 patients with CAIS. The missense variant p.G821E was predicted as deleterious, damaging, disease-causing, and likely functionally inactive by in silico analysis and protein modelling study. Two previously not reported pathogenic/likely pathogenic variants, including p.R386H and p.G396R, were identified in patients with PAIS. This study contributes in expanding the spectrum of pathogenic variants in the AR gene in patients with AIS. Only 31% patients with a provisional diagnosis of PAIS had pathogenic variants in the AR gene, suggesting other possible mechanisms or candidate genes may be responsible for such a phenotypic presentation.

Accuracy of Urinary Etiocholanolone/Androsterone Ratio as Alternative to Serum Testosterone/Dihydrotestosterone Ratio for Diagnosis of 5 Alpha-reductase Type 2 Deficiency Patients and Carriers in Indonesia

BACKGROUND

The 5 Alpha-reductase type 2 deficiency (5ARD2) is an inherited condition, which clinically presents as variable degrees of under virilization in affected 46,XY individuals. In the diagnostic pathway of 5ARD2, the testosterone/dihydrotestosterone (T/DHT) ratio is broadly employed before molecular analysis of the SRD5A2 gene. However, due to cost-benefit considerations, the DHT test in our country is routinely lacking in clinical settings; therefore, we considered applying the urinary etiocholanolone/androsterone (Et/An) ratio as an alternative test.

OBJECTIVES

We aimed to determine the diagnostic value of the urinary Et/An ratio versus the T/DHT ratio in diagnosing 5ARD2 patients and carriers.

METHODS

Sixty-six suspected 5ARD2 46,XY disorders of sex development (DSD) individuals and 95 family members were recruited in the study. Their clinical manifestations, T/DHT and urinary Et/An ratios, and SRD5A2 genes were analyzed. Using molecular analysis of the SRD5A2 gene as the gold standard, we compared the accuracy of both ratios in diagnosing 5ARD2 patients and carriers with receiver operating characteristic (ROC) curve analysis.

RESULTS

Thirty-seven patients were confirmed molecularly to have 5ARD2, and the rest (n = 29) were assessed as normal controls, while in the carrier group, 53 were molecularly confirmed as carriers and 42 as controls. The AUCs (areas under the curve) of the T/DHT and urinary Et/An ratios were 57.7% (95% CI 43.0 - 72.4%, P > 0.05) and 79.7% (95% CI 69.0 - 90.4%, P < 0.001), respectively, in diagnosing 5ARD2 patients and 54.1% (95% CI 42.4 - 65.8%, P > 0.05) and 75.1% (95% CI 65.1 - 85.1%, P < 0.001), respectively, in diagnosing carriers. The cutoff value of the urinary Et/An ratio was set at ≥ 0.95 for detecting 5ARD2 patients and ≥ 0.99 for detecting carriers.

CONCLUSIONS

The testosterone/DHT ratio was inaccurate in diagnosing 5ARD2 patients. When molecular analysis for the SRD5A2 gene is lacking, the urinary Et/An ratio may be a useful test to diagnose 5ARD2 patients and carriers.

Discrete Hedgehog Factor Expression and Action in the Developing Phallus

Hypospadias is a failure of urethral closure within the penis occurring in 1 in 125 boys at birth and is increasing in frequency. While paracrine hedgehog signalling is implicated in the process of urethral closure, how these factors act on a tissue level to execute closure itself is unknown. This study aimed to understand the role of different hedgehog signalling members in urethral closure. The tammar wallaby (Macropus eugenii) provides a unique system to understand urethral closure as it allows direct treatment of developing offspring because mothers give birth to young before urethral closure begins. Wallaby pouch young were treated with vehicle or oestradiol (known to induce hypospadias in males) and samples subjected to RNAseq for differential expression and gene ontology analyses. Localisation of Sonic Hedgehog (SHH) and Indian Hedgehog (IHH), as well as the transcription factor SOX9, were assessed in normal phallus tissue using immunofluorescence. Normal tissue culture explants were treated with SHH or IHH and analysed for AR, ESR1, PTCH1, GLI2, SOX9, IHH and SHH expression by qPCR. Gene ontology analysis showed enrichment for bone differentiation terms in male samples compared with either female samples or males treated with oestradiol. Expression of SHH and IHH localised to specific tissue areas during development, akin to their compartmentalised expression in developing bone. Treatment of phallus explants with SHH or IHH induced factor-specific expression of genes associated with bone differentiation. This reveals a potential developmental interaction involved in urethral closure that mimics bone differentiation and incorporates discrete hedgehog activity within the developing phallus and phallic urethra.

Growth Pattern in Chinese Children With 5α-Reductase Type 2 Deficiency: A Retrospective Multicenter Study

BACKGROUND

5α-reductase type 2 deficiency (5αRD) is an autosomal recessive hereditary disease of the group of 46, XY disorders of sex development (DSD).

OBJECTIVE

To study the growth pattern in Chinese pediatric patients with 5αRD.

SUBJECTS

Data were obtained from 141 patients with 5αRD (age: 0-16 years old) who visited eight pediatric endocrine centers from January 2010 to December 2017.

METHODS

In this retrospective cohort study, height, weight, and other relevant data were collected from the multicenter hospital registration database. Baseline luteinizing hormone (LH), follicle stimulating hormone (FSH), testosterone (T), and dihydrotestosterone (DHT) after human chorionic gonadotropin (HCG) stimulation test were measured by enzyme enhanced chemiluminescence assay. Bone age (BA) was assessed using the Greulich-Pyle (G-P) atlas. Growth curve was constructed based on λ-median-coefficient of variation method (LMS).

RESULTS

The height standard deviation scores (HtSDS) and weight standard deviation scores (WtSDS) in 5αRD children were in the normal range as compared to normal boys. Significantly higher HtSDS was observed in patients with 5αRD who were <1 year old (t = 3.658, 2.103, P = 0.002, 0.048, respectively), and higher WtSDS in those <6 months old (t = 2.756, P = 0.012). Then HtSDS and WtSDS decreased gradually and fluctuated near the median of the same age until 13 years. WtSDS in 5αRD children from northern China were significantly higher than those from the south (Z = -2.670, P = 0.008). The variation tendency of HtSDS in Chinese 5αRDs was consistent with the trend of stimulating T. HtSDS and stimulating T in the external masculinization score (EMS) <7 group were slightly higher than those in EMS ≥ 7 group without significant difference. Additionally, the ratio of BA over chronological age (BA/CA) was significantly <1 in children with 5αRD.

CONCLUSION

Children with 5αRD had a special growth pattern that was affected by high levels of T, while DHT played a very small role in it. Their growth accelerated at age <1 year, followed by slowing growth and fluctuating height near normal median boys' height. The BA was delayed in 5αRD children. Androgen treatment, which may be considered anyway for male 5αRD patients with a micropenis, may also be beneficial for growth.

Characterising SRD5A2 Gene Variants in 37 Indonesian Patients with 5-Alpha-Reductase Type 2 Deficiency

The 5-alpha-reductase type 2 deficiency (5ARD2) is an autosomal recessive condition associated with impairment in the conversion of testosterone to dihydrotestosterone. This condition leads to undervirilisation in 46,XY individuals. To date, there have been more than 100 variations identified in the gene responsible for 5ARD2 development (steroid 5-alpha-reductase 2, SRD5A2). However, few studies have examined the molecular characterisation of Indonesian 5ARD2 cases. In the current study, we analysed 37 subjects diagnosed with 46,XY DSD (disorders of sex development) with confirmed variations in the SRD5A2 gene. We examined results from testosterone/dihydrotestosterone (T/DHT) and urinary etiocholanolone/androsterone (Et/An) ratios, as well as from molecular and clinical analyses. Twelve variants in the SRD5A2 gene were identified, and 6 of which were novel, namely, c.34-38delGinsCCAGC, p.Arg50His, p.Tyr136 ^∗ , p.Gly191Arg, p.Phe194Ile, and p.Ile253Val variants. Moreover, we determined that 20 individuals contained harmful mutations, while the remaining 17 variants were benign. Those containing harmful mutations exhibited more severe phenotypes with median external genitalia masculinisation scores (EMS) of 3 (1.5-9) and were more likely to be diagnosed at a later age, reared as female, and virilised at pubertal age. In addition, the respective sensitivities for detecting severe 5ARD2 cases using T/DHT (cutoff: 10) and urinary Et/An ratios (cutoff: 0.95) were 85% and 90%, whereas mild cases were only identified with 64.7% and 47.1% sensitivity, respectively. Although we were unable to identify clear correlations between genotypic and phenotypic characteristics in this study, we clearly showed that individuals who were homozygous or compound heterozygous for any of the harmful mutations were more likely to exhibit classic 5ARD2 phenotypes, lower EMS, female assignment at birth, and virilisation during puberty. These results serve to inform the development of improved clinical and molecular 5ARD2 diagnostic approaches, specifically in Indonesian patients.