Concomitant mutations in the P450 oxidoreductase and androgen receptor genes presenting with 46,XY disordered sex development and androgenization at adrenarche

Exploring Novel Variants of the Cytochrome P450 Reductase Gene (POR) from the Genome Aggregation Database by Integrating Bioinformatic Tools and Functional Assays

Cytochrome P450 oxidoreductase (POR) is an essential redox partner for steroid and drug-metabolizing cytochromes P450 located in the endoplasmic reticulum. Mutations in POR lead to metabolic disorders, including congenital adrenal hyperplasia, and affect the metabolism of steroids, drugs, and xenobiotics. In this study, we examined approximately 450 missense variants of the POR gene listed in the Genome Aggregation Database (gnomAD) using eleven different in silico prediction tools. We found that 64 novel variants were consistently predicted to be disease-causing by most tools. To validate our findings, we conducted a population analysis and selected two variations in POR for further investigation. The human POR wild type and the R268W and L577P variants were expressed in bacteria and subjected to enzyme kinetic assays using a model substrate. We also examined the activities of several cytochrome P450 proteins in the presence of POR (WT or variants) by combining P450 and reductase proteins in liposomes. We observed a decrease in enzymatic activities (ranging from 35% to 85%) of key drug-metabolizing enzymes, supported by POR variants R288W and L577P compared to WT-POR. These results validate our approach of curating a vast amount of data from genome projects and provide an updated and reliable reference for diagnosing POR deficiency.

Congenital adrenal hyperplasia, disorders of sex development, and infertility in patients with POR gene pathogenic variants: a systematic review of the literature

BACKGROUND

P450 oxidoreductase (POR) deficiency (PORD) is characterized by congenital adrenal hyperplasia (CAH) and disorders of sex development (DSD) in both sexes. PORD can also associate with skeletal defects. However, the prevalence of these phenotypes is unknown.

AIM

To evaluate the prevalence of CAH, DSD, and infertility of patients with POR gene pathogenic variants by a systematic review of the literature.

METHODS

The literature search was performed through PubMed, MEDLINE, Cochrane, Academic One Files, Google Scholar, and Scopus databases. All studies reporting information on CAH, DSD, testicular adrenal rest tumor (TARTs), and fertility in patients with POR gene pathogenic variants were included. Finally, the prevalence of abnormal phenotypes was calculated.

RESULTS

Of the 246 articles initially retrieved, only 48 were included for a total of 119 (46 males and 73 females) patients with PORD. We also included the case of a male patient who consulted us for CAH and TARTs but without DSD. This patient, found to be a carrier of combined heterozygous POR mutation, reached fatherhood spontaneously. All the patients found had CAH. The presence of DSD was found in 65.2%, 82.1%, and 82.1% of patients with compound heterozygosity, homozygosity, or monoallelic heterozygous variants, respectively. The prevalence was significantly higher in females than in males. The prevalence of TARTs in patients with PORD is 2.7%. Only 5 women with PORD became pregnant after assisted reproductive techniques and delivered a healthy baby. Except for the recently reported proband, no other studies focused on male infertility in patients with POR gene variants.

CONCLUSION

This systematic review of the literature reports the prevalence of CAH, DSD, and TARTs in patients with PORD. The unknown prevalence of POR gene pathogenetic variants and the paucity of studies investigating fertility do not allow us to establish whether PORD is associated with infertility. Further studies on both women and men are needed to clarify this relationship.

In Silico Analysis of PORD Mutations on the 3D Structure of P450 Oxidoreductase

Cytochrome P450 oxidoreductase (POR) is a membrane-bound flavoprotein that helps in transferring electrons from its NADPH domain to all cytochrome P450 (CYP450) enzymes. Mutations in the POR gene could severely affect the metabolism of steroid hormones and the development of skeletal muscles, a condition known as Cytochrome P450 oxidoreductase deficiency (PORD). PORD is associated with clinical presentations of disorders of sex development, Antley and Bixler’s syndrome (ABS), as well as an abnormal steroid hormone profile. We have performed an in silico analysis of POR 3D X-ray protein crystal structure to study the effects of reported mutations on the POR enzyme structure. A total of 32 missense mutations were identified, from 170 PORD patients, and mapped on the 3D crystal structure of the POR enzyme. In addition, five of the missense mutations (R457H, A287P, D210G, Y181D and Y607C) were further selected for an in-depth in silico analysis to correlate the observed changes in POR protein structure with the clinical phenotypes observed in PORD patients. Overall, missense mutations found in the binding sites of POR cofactors could lead to a severe form of PORD, emphasizing the importance of POR cofactor binding domains in transferring electrons to the CYP450 enzyme family.

Normal and Premature Adrenarche

Adrenarche is the maturational increase in adrenal androgen production that normally begins in early childhood. It results from changes in the secretory response to adrenocorticotropin (ACTH) that are best indexed by dehydroepiandrosterone sulfate (DHEAS) rise. These changes are related to the development of the zona reticularis (ZR) and its unique gene/enzyme expression pattern of low 3ß-hydroxysteroid dehydrogenase type 2 with high cytochrome b5A, sulfotransferase 2A1, and 17ß-hydroxysteroid dehydrogenase type 5. Recently 11-ketotestosterone was identified as an important bioactive adrenarchal androgen. Birth weight, body growth, obesity, and prolactin are related to ZR development. Adrenarchal androgens normally contribute to the onset of sexual pubic hair (pubarche) and sebaceous and apocrine gland development. Premature adrenarche causes ≥90% of premature pubarche (PP). Its cause is unknown. Affected children have a significantly increased growth rate with proportionate bone age advancement that typically does not compromise growth potential. Serum DHEAS and testosterone levels increase to levels normal for early female puberty. It is associated with mildly increased risks for obesity, insulin resistance, and possibly mood disorder and polycystic ovary syndrome. Between 5% and 10% of PP is due to virilizing disorders, which are usually characterized by more rapid advancement of pubarche and compromise of adult height potential than premature adrenarche. Most cases are due to nonclassic congenital adrenal hyperplasia. Algorithms are presented for the differential diagnosis of PP. This review highlights recent advances in molecular genetic and developmental biologic understanding of ZR development and insights into adrenarche emanating from mass spectrometric steroid assays.

Oligogenic Origin of Differences of Sex Development in Humans

Sex development is a very complex biological event that requires the concerted collaboration of a large network of genes in a spatial and temporal correct fashion. In the past, much has been learned about human sex development from monogenic disorders/differences of sex development (DSD), but the broad spectrum of phenotypes in numerous DSD individuals remains a conundrum. Currently, the genetic cause of less than 50% of DSD individuals has been solved and oligogenic disease has been proposed. In recent years, multiple genetic hits have been found in individuals with DSD thanks to high throughput sequencing. Our group has been searching for additional genetic hits explaining the phenotypic variability over the past years in two cohorts of patients: 46,XY DSD patients carriers of NR5A1 variants and 46,XY DSD and 46,XX DSD with MAMLD1 variants. In both cohorts, our results suggest that the broad phenotypes may be explained by oligogenic origin, in which multiple hits may contribute to a DSD phenotype, unique to each individual. A search for an underlying network of the identified genes also revealed that a considerable number of these genes showed interactions, suggesting that genetic variations in these genes may affect sex development in concert.

P450 Oxidoreductase Deficiency: A Systematic Review and Meta-analysis of Genotypes, Phenotypes, and Their Relationships

CONTEXT

P450 oxidoreductase deficiency (PORD) is a rare genetic disorder that is associated with significant morbidity. However there has been limited analysis of reported PORD cases.

OBJECTIVE

To determine, based on the cohort of reported PORD cases, genotype-phenotype relationships for skeletal malformations, maternal virilisation in pregnancy, adrenal insufficiency, and disorders of sexual development (DSD).

DATA SOURCES

PubMed and Web of Science from January 2004 to February 2018.

STUDY SELECTION

Published case reports/series of patients with PORD. Eligible patients were unique, had biallelic mutations, and their clinical features were reported.

DATA EXTRACTION

Patient data were manually extracted from the text of case reports/series. A malformation score, representing the severity of skeletal malformations, was calculated for each patient.

DATA SYNTHESIS

Of the 211 patients published in the literature, 90 were eligible for inclusion. More than 60 unique mutations were identified in this cohort. Four groups of mutations were identified, through regression modeling, as having significantly different skeletal malformation scores. Maternal virilization in pregnancy, reported for 21% of patients, was most common for R457H mutations. Adrenal insufficiency occurred for the majority of patients (78%) and was typically mild, with homozygous R457H mutations being the least deficient. DSD affected most patients (72%), but were less common for males (46XY) with homozygous R457H mutations.

CONCLUSIONS

PORD is a complex disorder with many possible mutations affecting a large number of enzymes. By analyzing the cohort of reported PORD cases, this study identified clear relationships between genotype and several important phenotypic features.

Monogenic Disorders of Adrenal Steroidogenesis

Disorders of adrenal steroidogenesis comprise autosomal recessive conditions affecting steroidogenic enzymes of the adrenal cortex. Those are located within the 3 major branches of the steroidogenic machinery involved in the production of mineralocorticoids, glucocorticoids, and androgens. This mini review describes the principles of adrenal steroidogenesis, including the newly appreciated 11-oxygenated androgen pathway. This is followed by a description of pathophysiology, biochemistry, and clinical implications of steroidogenic disorders, including mutations affecting cholesterol import and steroid synthesis, the latter comprising both mutations affecting steroidogenic enzymes and co-factors required for efficient catalysis. A good understanding of adrenal steroidogenic pathways and their regulation is crucial as the basis for sound management of these disorders, which in the majority present in early childhood.

P450 Oxidoreductase deficiency: Analysis of mutations and polymorphisms

Cytochrome P450 oxidoreductase (POR) is required for metabolic reactions of steroid and drug metabolizing cytochrome P450 proteins located in endoplasmic reticulum. Mutations in POR cause a complex set of disorders resembling combined deficiencies of multiple steroid metabolizing enzymes. The P450 oxidoreductase deficiency (PORD) was first reported in patients with symptoms of defects in steroidogenic cytochrome P450 enzymes and ambiguous genitalia, and bone malformation features resembling Antley-Bixler syndrome. POR is now classified as a separate and rare form of congenital adrenal hyperplasia (CAH), which may cause disorder of sexual development (DSD). Since the initial description of PORD in 2004, a large number of POR mutations and polymorphisms have been described. In this report we have performed computational analysis of mutations and polymorphisms in POR linked to metabolism of steroids and xenobiotics and pathology of PORD from the reported cases. The mutations in POR that were identified in patients with disruption of steroidogenesis also have severe effects on cytochrome P450 proteins involved in metabolism of drugs. Different variations in POR show a range of diverse effects on different partner proteins that are often linked to the location of the particular variants. The variations in POR that cause defective binding of co-factors always have damaging effects on all partner proteins, while the mutations causing subtle structural changes may lead to altered interaction with partner proteins and the overall effect may be different for each individual partner. Computational analysis of available sequencing data and mutation analysis shows that Japanese (R457H), Caucasian (A287P) and Turkish (399-401) populations can be linked to unique founder mutations. Other mutations identified so far were identified as rare alleles or in single isolated reports. The common polymorphism of POR is the variant A503V which can be found in about 27% of alleles in general population but there are remarkable differences among different sub populations.

Identification of an AR Mutation-Negative Class of Androgen Insensitivity by Determining Endogenous AR Activity

CONTEXT

Only approximately 85% of patients with a clinical diagnosis complete androgen insensitivity syndrome and less than 30% with partial androgen insensitivity syndrome can be explained by inactivating mutations in the androgen receptor (AR) gene.

OBJECTIVE

The objective of the study was to clarify this discrepancy by in vitro determination of AR transcriptional activity in individuals with disorders of sex development (DSD) and male controls.

DESIGN

Quantification of DHT-dependent transcriptional induction of the AR target gene apolipoprotein D (APOD) in cultured genital fibroblasts (GFs) (APOD assay) and next-generation sequencing of the complete coding and noncoding AR locus.

SETTING

The study was conducted at a university hospital endocrine research laboratory.

PATIENTS

GFs from 169 individuals were studied encompassing control males (n = 68), molecular defined DSD other than androgen insensitivity syndrome (AIS; n = 18), AR mutation-positive AIS (n = 37), and previously undiagnosed DSD including patients with a clinical suspicion of AIS (n = 46).

INTERVENTION(S)

There were no interventions.

MAIN OUTCOME MEASURE(S)

DHT-dependent APOD expression in cultured GF and AR mutation status in 169 individuals was measured.

RESULTS

The APOD assay clearly separated control individuals (healthy males and molecular defined DSD patients other than AIS) from genetically proven AIS (cutoff < 2.3-fold APOD-induction; 100% sensitivity, 93.3% specificity, P < .0001). Of 46 DSD individuals with no AR mutation, 17 (37%) fell below the cutoff, indicating disrupted androgen signaling.

CONCLUSIONS

AR mutation-positive AIS can be reliably identified by the APOD assay. Its combination with next-generation sequencing of the AR locus uncovered an AR mutation-negative, new class of androgen resistance, which we propose to name AIS type II. Our data support the existence of cellular components outside the AR affecting androgen signaling during sexual differentiation with high clinical relevance.

Clinical review of 95 patients with 46,XX disorders of sex development based on the new Chicago classification

STUDY OBJECTIVE

The aim of our study was to determine the etiologic distribution of 46,XX disorder of sexual development (DSD) according to the new DSD classification system and to evaluate the clinical features of this DSD subgroup in our patient cohort.

PARTICIPANTS

The evaluation criteria and clinical findings of 95 46,XX patients were described by clinical presentation, gonadal morphology, genital anatomy, associated dysmorphic features, presence during prenatal period with/without postnatal virilization, hormonal characteristics, and presence or absence of steroidogenic defects among 319 patients with DSD.

RESULTS

Types and ratios of each presentation of our 95 patients with 46,XX DSD were as follows: 82 had androgen excess (86.3%): (74 had classical congenital adrenal hyperplasia, 2 had CAH variant possibility of P450-oxidoreductase gene defect), 6 had disorders of ovarian development (6.3%): (1 patient had gonadal dysgenesis with virilization at birth with bilateral streak gonad, 4 patients had complete gonadal dysgenesis, and 1 patient had ovotesticular DSD) and 7 had other 46,XX DSD. Two sisters, who had 46,XX complete gonadal dysgenesis,were diagnosed with Perrault Syndrome with ovarian failure due to streak gonads and associated with sensorineural deafness.

CONCLUSION

46,XX DSD are usually derived from intrauterine virilization and CAH is the most common cause of 46,XX DSD due to fetal androgen exposure.

Pharmacogenomics of human P450 oxidoreductase

Cytochrome P450 oxidoreductase (POR) supports reactions of microsomal cytochrome P450 which metabolize drugs and steroid hormones. Mutations in POR cause disorders of sexual development. P450 oxidoreductase deficiency (PORD) was initially identified in patients with Antley–Bixler syndrome (ABS) but now it has been established as a separate disorder of sexual development (DSD). Here we are summarizing the work on variations in POR related to metabolism of drugs and xenobiotics. We have compiled mutation data on reported cases of PORD from clinical studies. Mutations found in patients with defective steroid profiles impact metabolism of steroid hormones as well as drugs. Some trends are emerging that establish certain founder mutations in distinct populations, with Japanese (R457H), Caucasian (A287P), and Turkish (399–401) populations showing repeated findings of similar mutations. Most other mutations are found as single occurrences. A large number of different variants in POR gene with more than 130 amino acid changes are now listed in databases. Among the polymorphisms, the A503V is found in about 30% of all alleles but there are some differences across different population groups.

Prenatal diagnosis of congenital adrenal hyperplasia caused by P450 oxidoreductase deficiency

CONTEXT

Mutations in the electron donor enzyme P450 oxidoreductase (POR) result in congenital adrenal hyperplasia with apparent combined 17α-hydroxylase/17,20 lyase and 21-hydroxylase deficiencies, also termed P450 oxidoreductase deficiency (PORD). Major clinical features present in PORD are disordered sex development in affected individuals of both sexes, glucocorticoid deficiency, and multiple skeletal malformations.

OBJECTIVE

The objective of the study was to establish a noninvasive approach to prenatal diagnosis of PORD including assessment of malformation severity to facilitate optimized prenatal diagnosis and timely treatment.

DESIGN

We analyzed 20 pregnancies with children homozygous or compound heterozygous for disease-causing POR mutations and 1 pregnancy with a child carrying a heterozygous POR mutation by recording clinical and biochemical presentations and fetal ultrasound findings. In 4 of the pregnancies (3 homozygous and 1 heterozygous for disease-causing POR mutations), prenatal analysis of steroid metabolite excretion in maternal urine was carried out by gas chromatography/mass spectrometry during gestational weeks 11-23.

RESULTS

Pregnancy complications in our cohort included maternal virilization (6 of 20) with onset in the second trimester. Seven pregnant women presented with low unconjugated estriol at prenatal screening (triple or quadruple antenatal screening test). Overt dysmorphic features were noted in 19 of the 20 babies at birth but observed in only 5 by prenatal ultrasound. These 5 had the most severe malformation phenotypes and poor outcome, whereas the other babies showed normal development. Steroid profiling of maternal urine revealed significantly increased steroids of fetal origin, namely the pregnenolone metabolite epiallopregnanediol and the androgen metabolite androsterone, with concomitant low values for estriol. Diagnostic steroid ratios conclusively indicated PORD as early as gestational week 12. In the heterozygous pregnancy, steroid ratios were only slightly elevated and estriol excretion was normal.

CONCLUSION

Prenatal diagnosis in PORD is readily established via urinary steroid metabolite analysis of maternal urine. Visible malformations at prenatal ultrasound predict a severe malformation phenotype.

NADPH P450 oxidoreductase: structure, function, and pathology of diseases

Cytochrome P450 oxidoreductase (POR) is an enzyme that is essential for multiple metabolic processes, chiefly among them are reactions catalyzed by cytochrome P450 proteins for metabolism of steroid hormones, drugs and xenobiotics. Mutations in POR cause a complex set of disorders that often resemble defects in steroid metabolizing enzymes 17α-hydroxylase, 21-hydroxylase and aromatase. Since our initial reports of POR mutations in 2004, more than 200 different mutations and polymorphisms in POR gene have been identified. Several missense variations in POR have been tested for their effect on activities of multiple steroid and drug metabolizing P450 proteins. Mutations in POR may have variable effects on different P450 partner proteins depending on the location of the mutation. The POR mutations that disrupt the binding of co-factors have negative impact on all partner proteins, while mutations causing subtle structural changes may lead to altered interaction with specific partner proteins and the overall effect may be different for each partner. This review summarizes the recent discoveries related to mutations and polymorphisms in POR and discusses these mutations in the context of historical developments in the discovery and characterization of POR as an electron transfer protein. The review is focused on the structural, enzymatic and clinical implications of the mutations linked to newly identified disorders in humans, now categorized as POR deficiency.

Evidence for DNA-binding domain--ligand-binding domain communications in the androgen receptor

DNA binding as well as ligand binding by nuclear receptors has been studied extensively. Both binding functions are attributed to isolated domains of which the structure is known. The crystal structure of a complete receptor in complex with its ligand and DNA-response element, however, has been solved only for the peroxisome proliferator-activated receptor γ (PPARγ)-retinoid X receptor α (RXRα) heterodimer. This structure provided the first indication of direct interactions between the DNA-binding domain (DBD) and ligand-binding domain (LBD). In this study, we investigated whether there is a similar interface between the DNA- and ligand-binding domains for the androgen receptor (AR). Despite the structural differences between the AR- and PPARγ-LBD, a combination of in silico modeling and docking pointed out a putative interface between AR-DBD and AR-LBD. The surfaces were subjected to a point mutation analysis, which was inspired by known AR mutations described in androgen insensitivity syndromes and prostate cancer. Surprisingly, AR-LBD mutations D695N, R710A, F754S, and P766A induced a decrease in DNA binding but left ligand binding unaffected, while the DBD-residing mutations K590A, K592A, and E621A lowered the ligand-binding but not the DNA-binding affinity. We therefore propose that these residues are involved in allosteric communications between the AR-DBD and AR-LBD.

Genotype-phenotype analysis in congenital adrenal hyperplasia due to P450 oxidoreductase deficiency

CONTEXT

P450 oxidoreductase deficiency (PORD) is a unique congenital adrenal hyperplasia variant that manifests with glucocorticoid deficiency, disordered sex development (DSD), and skeletal malformations. No comprehensive data on genotype-phenotype correlations in Caucasian patients are available.

OBJECTIVE

The objective of the study was to establish genotype-phenotype correlations in a large PORD cohort.

DESIGN

The design of the study was the clinical, biochemical, and genetic assessment including multiplex ligation-dependent probe amplification (MLPA) in 30 PORD patients from 11 countries.

RESULTS

We identified 23 P450 oxidoreductase (POR) mutations (14 novel) including an exonic deletion and a partial duplication detected by MLPA. Only 22% of unrelated patients carried homozygous POR mutations. p.A287P was the most common mutation (43% of unrelated alleles); no other hot spot was identified. Urinary steroid profiling showed characteristic PORD metabolomes with variable impairment of 17α-hydroxylase and 21-hydroxylase. Short cosyntropin testing revealed adrenal insufficiency in 89%. DSD was present in 15 of 18 46,XX and seven of 12 46,XY individuals. Homozygosity for p.A287P was invariably associated with 46,XX DSD but normal genitalia in 46,XY individuals. The majority of patients with mild to moderate skeletal malformations, assessed by a novel scoring system, were compound heterozygous for missense mutations, whereas nearly all patients with severe malformations carried a major loss-of-function defect on one of the affected alleles.

CONCLUSIONS

We report clinical, biochemical, and genetic findings in a large PORD cohort and show that MLPA is a useful addition to POR mutation analysis. Homozygosity for the most frequent mutation in Caucasians, p.A287P, allows for prediction of genital phenotype and moderate malformations. Adrenal insufficiency is frequent, easily overlooked, but readily detected by cosyntropin testing.

Pubertal presentation in seven patients with congenital adrenal hyperplasia due to P450 oxidoreductase deficiency

CONTEXT

P450 oxidoreductase (POR) is a crucial electron donor to all microsomal P450 cytochrome (CYP) enzymes including 17α-hydroxylase (CYP17A1), 21-hydroxylase (CYP21A2) and P450 aromatase. Mutant POR causes congenital adrenal hyperplasia with combined glucocorticoid and sex steroid deficiency. P450 oxidoreductase deficiency (ORD) commonly presents neonatally, with disordered sex development in both sexes, skeletal malformations, and glucocorticoid deficiency.

OBJECTIVE

The aim of the study was to describe the clinical and biochemical characteristics of ORD during puberty.

DESIGN

Clinical, biochemical, and genetic assessment of seven ORD patients (five females, two males) presenting during puberty was conducted.

RESULTS

Predominant findings in females were incomplete pubertal development (four of five) and large ovarian cysts (five of five) prone to spontaneous rupture, in some only resolving after combined treatment with estrogen/progestin, GnRH superagonists, and glucocorticoids. Pubertal development in the two boys was more mildly affected, with some spontaneous progression. Urinary steroid profiling revealed combined CYP17A1 and CYP21A2 deficiencies indicative of ORD in all patients; all but one failed to mount an appropriate cortisol response to ACTH stimulation indicative of adrenal insufficiency. Diagnosis of ORD was confirmed by direct sequencing, demonstrating disease-causing POR mutations.

CONCLUSION

Delayed and disordered puberty can be the first sign leading to a diagnosis of ORD. Appropriate testosterone production during puberty in affected boys but manifest primary hypogonadism in girls with ORD may indicate that testicular steroidogenesis is less dependent on POR than adrenal and ovarian steroidogenesis. Ovarian cysts in pubertal girls may be driven not only by high gonadotropins but possibly also by impaired CYP51A1-mediated production of meiosis-activating sterols due to mutant POR.