Abstract
Mutations in bone morphogenetic protein receptor type 2 (BMPR2) cause familial pulmonary arterial hypertension (FPAH), but the penetrance is reduced and females are significantly overrepresented. In addition, gene expression data implicating the oestrogen-metabolising enzyme CYP1B1 suggests a detrimental role of oestrogens or oestrogen metabolites. We examined genetic and metabolic markers of altered oestrogen metabolism in subjects with a BMPR2 mutation.
Genotypes for CYP1B1 Asn453Ser (N453S) were determined for 140 BMPR2 mutation carriers (86 females and 54 males). Nested from those subjects, a case–control study of urinary oestrogen metabolite levels (2-hydroxyoestrogen (2-OHE) and 16α-hydroxyoestrone (16α-OHE1)) was conducted in females (five affected mutation carriers versus six unaffected mutation carriers).
Among females, there was four-fold higher penetrance among subjects homozygous for the wild-type genotype (N/N) than those with N/S or S/S genotypes (p = 0.005). Consistent with this finding, the 2-OHE/16α-OHE1 ratio was 2.3-fold lower in affected mutation carriers compared to unaffected mutation carriers (p = 0.006).
Our findings suggest that variations in oestrogens and oestrogen metabolism modify FPAH risk. Further investigation of the role of oestrogens in this disease with profound sex bias may yield new insights and, perhaps, therapeutic interventions.
Familial pulmonary arterial hypertension (FPAH) is caused, in 80% of families, by autosomal dominant mutations in the gene encoding bone morphogenetic protein receptor type 2 (BMPR2). However, it is unclear why BMPR2 mutations have reduced penetrance, and why, after puberty, females with mutations are about 2.5-fold more likely to develop FPAH than males 1. The excess of adult females suggests oestrogen as a possible disease modifier, as does a lower prevalence among pre-pubertal females compared with males 2, 3. A detrimental role for oestrogens in FPAH is supported by studies that show oestrogen is a potent mitogen of pulmonary vascular smooth muscle cells 4.
Recently, our studies of BMPR2 mutation carriers implicated altered oestrogen metabolism as a key factor in the penetrance of FPAH in females. Using genome expression arrays confirmed by RT-PCR, we found significantly decreased transcript levels of cytochrome P450 1B1 (CYP1B1; a cytochrome p450 family enzyme critical to oestrogen metabolism) in affected female compared with unaffected female BMPR2 mutation carriers 5. CYP1B1 is highly expressed in the lung, and it is a modifier for several cancers, including lung and breast cancers 6–14. CYP1B1 catalyses the oxidation of oestrogens to 2-hydroxy (2-OHE) and 4-hydroxy (4-OHE) oestrogens, and metabolises environmental toxins, including tobacco smoke 15.
Oxidation of oestrogens also occurs by hydroxylation at the C-16 position by other P450 enzymes, predominantly resulting in 16α-hydroxyoestrone (16α-OHE1) in extrahepatic tissues 16, 17. Data suggest that, unlike the weak mitogen 2-OHE, 16α-OHE1 stimulates cellular proliferation by constitutively activating the oestrogen receptor. In addition to being more mitogenic than 2-OHE, 16α-OHE1 may also be more genotoxic. Thus, individuals who metabolise a larger proportion of oestrogen to 16α-OHE1 may be at increased risk of diseases that result from both the mitogenic and genotoxic effects of oestrogens, such as breast and prostate cancers 18.
We hypothesised that oestrogens may contribute to FPAH, and that the risk may be altered by variations in CYP1B1. To test this, we studied the relationship between the CYP1B1 Asn453Ser (N453S) genetic polymorphism and FPAH penetrance in BMPR2 mutation carriers 19. We also determined the urinary 2-OHE/16α-OHE1 ratio, one biological measure of CYP1B1 activity, among a group of females nested from our polymorphism study. Our genetic and biochemical results suggest that imbalanced oestrogen metabolism may contribute to the pathogenesis of FPAH, and that further analysis of this pathway may lead to new biological markers and therapeutic treatments for all forms of PAH.
MATERIALS AND METHODS
Study population
The Vanderbilt pulmonary arterial hypertension (PAH) research registry provided the study base from which the subjects in this study were recruited. The registry contains participants with idiopathic PAH, as well as participants affected with FPAH, bloodline relatives of those affected and unrelated family members of all ages. Study subjects for the registry, the majority of whom are now deceased, were recruited via the Vanderbilt Pulmonary Hypertension Center (Nashville, TN, USA), the Pulmonary Hypertension Association (Silver Spring, MD, USA) and the National Institutes of Health (Bethesda, MD, USA) clinical trials website (http://clinicaltrials.gov). BMPR2 mutations were detected in a large proportion of families with FPAH followed by the registry. These mutations were heterogeneous and included nonsense, missense and frameshift alterations, as well as insertion–deletion mutations that lead to splicing errors. In addition the detected BMPR2 mutations were heterogeneous in location across the gene.
For this study, subjects diagnosed with FPAH who had a BMPR2 mutation were considered affected BMPR2 mutation carriers (AMCs). Subjects who have never been diagnosed with FPAH but have a BMPR2 mutation were considered unaffected BMPR2 mutation carriers (UMCs). Clinical information on AMCs, including cardiac catheterisation data and therapeutic regimens, was available from medical history and records. Specialist physicians in their geographic regions identified subjects as affected with FPAH, and our investigators reviewed all medical records for accuracy of diagnosis. We defined FPAH diagnostically, either by autopsy results showing plexogenic pulmonary arteriopathy in the absence of alternative causes, such as congenital heart disease, or by clinical and cardiac catheterisation criteria. These criteria included a mean pulmonary arterial pressure >25 mmHg with a pulmonary capillary or left atrial pressure <15 mmHg, and exclusion of other causes of pulmonary hypertension, in accordance with accepted international standards of diagnostic criteria 20. All subjects designated as UMCs were seen by a physician within 1 yr of sample acquisition, and without symptoms or signs of PAH at time of sample acquisition.
The Vanderbilt University Medical Center (Nashville, TN, USA) institutional review board approved all study protocols. All participants gave informed written consent to participate in genetic and clinical studies. Samples were obtained at the time of clinic visit or hospitalisation, or by post via a kit for collection of whole blood and DNA.
For our CYP1B1 N453S polymorphism study, 140 BMPR2 mutation carriers were recruited, including 86 females and 54 males. Relevant clinical characteristics for these subjects are listed in table 1⇓. Subsequently, a nested case–control cohort study of urinary oestrogen metabolites was performed in 11 female BMPR2 mutation carriers. This cohort was composed of five female patients previously diagnosed with FPAH (AMCs) and six females who were UMCs. Relevant clinical characteristics for the AMCs included in this nested cohort are listed in table 2⇓. The following criteria had to be met for inclusion in the nested case–control cohort study: 1) living; 2) female; 3) not currently pregnant; 4) not pregnant within the past 1 yr; 5) not hysterectomised or with an oophorectomy; and 6) not currently taking exogenous steroid hormones. UMCs were not eligible if taking prescribed or over-the-counter medications or supplements.
Overnight 12-h urine void collections were obtained from all participants following previously established protocols. Urine samples from pre-menopausal females were obtained in the mid-luteal phase of the menstrual cycle (4–6 days prior to the estimated start of the next menses) and subsequent menses confirmed. Urine was stored at 4°C on ice or in a refrigerator until processing for storage. To prevent oxidation of metabolites during storage, 50 mg of ascorbic acid was added to 50 mL of urine (1 mg·mL−1) and the sample mixed thoroughly. Samples were aliquoted into conical tubes and stored at -80°C within 24 h of collection.
Laboratory analyses
Genomic DNA was isolated from whole blood using Puregene® DNA Purification Kit (Gentra, Minneapolis, MN, USA) according to the manufacturer's protocol for genotyping. As previously described, established lymphoblastoid cell lines of participants were grown in RPMI-1640 containing 15% fetal bovine serum, penicillin 100 μg·mL−1 and streptomycin 100 μg·mL−1 in T-25 flasks to provide RNA. Cells were grown to ∼1×106 cells per flask and incubated in the presence or absence of puromycin (100 μg·mL−1) for 16 h before harvesting. Cells were harvested using RNAeasy Mini Kit including the optional DNase treatment (Qiagen Inc., Valencia, CA, USA) 21. We performed BMPR2 gene mutation detection by sequencing exons and exon intron boundaries of genomic DNA and RT-PCR analysis as previously described 21, 22. All detected BMPR2 mutations in this study have been previously reported, and are included in a recent summary of detectable BMPR2 mutations 23.
Genotyping of the samples for the CYP1B1 N453S genetic polymorphism (rs1800440) was performed using the TaqMan SNP Assay, C_11642651_30 from Applied Biosystems (Foster City, CA, USA). The samples were amplified with Universal Master Mix (Applied Biosystems) and run and analysed on the 7500 Real Time PCR System (Applied Biosystems) using allelic discrimination cycling conditions. Sample genotyping was performed in triplicate by staff blinded to the subjects' phenotype; nontemplate controls were also included. The accuracy of genotyping data was tested by direct DNA sequencing of ∼15% of samples analysed, and all results were in 100% concordance.
All urinary oestrogen metabolite study samples were run together in a blinded manner. The samples were thawed at room temperature, centrifuged at 2000×g for 10 min, and 10 mL of clear urine was extracted for each assay. Urinary oestrogen metabolites were quantitated using a monoclonal antibody-based enzyme assay (Estramet 2/16TM; Immuna Care, Blue Bell, PA, USA). This assay was previously validated by gas chromatography–mass spectroscopy methods and found to reliably measure 2-OHE and 16α-OHE1 levels 24–26. Urine samples from pre-menopausal females were pre-diluted 1:4 using the manufacturer's diluent and assayed according to the manufacturer's protocol. All urine samples were run in sextuplet, with controls, and standards in triplicate. Samples, controls and standards were pipetted into individual microtubes in sextuplet, with 190 μL hydrolysis buffer added to each tube followed by a 2-h covered incubation at room temperature. Neutralisation buffer (200 μL per well) was added to the enzyme digest to bring the sample to a pH of 7.0, and 75 μL of each sample was added to two 96-well microtitre plates (one for the determination of 2-OHE and one for 16α-OHE1) coated with polyclonal rabbit anti-mouse immunoglobulin. Immediately prior to use, the plates were activated by washing with buffer using an automatic plate washer (Bio-Rad model 1575 Immunowash; Bio-Rad Laboratories Inc., Hercules, CA, USA) according to the manufacturer's package insert. To each respective assay plate was added 75 μL of a conjugate buffer solution containing metabolite-specific 2-OHE:alkaline phosphatase or 16α-OHE1:alkaline phosphatase. The plates were sealed and incubated for 3 h at room temperature, then washed with an automatic plate washer according to the manufacturer's insert. Stabilised paranitrophenyl phosphate (100 μL) was added to each well, and the enzyme product (a colour dye) was inversely proportional to the concentration of free antigen (oestrogen metabolite).
Absorbance was read kinetically by end-point analysis with a Bio-Rad iMark Microplate Absorbance Reader (Bio-Rad Laboratories Inc.) at 405 nm 60 min after adding the paranitrophenyl phosphate solution. The data were reduced using Microplate Manager 6 Application Software (Bio-Rad Laboratories Inc.). A four-parameter fit curve of the log of the concentration versus absorbance at 405 nm yielded a sigmoidal standard curve with a linear range of 1–15 ng·mL−1 of urine. Urinary creatinine was quantitated in duplicate using a commercially available colorimetric assay (Ortho-Clinical Diagnostics, Rochester, NY, USA) and was used to normalise oestrogen metabolite levels. Results for the concentration of each metabolite were reported as nanogrammes per millilitre per milligram of creatinine. Samples from each patient were assayed together in a single assay to eliminate inter-assay variation.
Statistical analysis
We used the Chi-squared test to determine if the CYP1B1 genotypes were in Hardy–Weinberg equilibrium and to determine the p-values for differences observed in genotype frequencies between AMCs and UMCs. Associations between CYP1B1 genotypes and FPAH penetrance were examined by conditional logistic regression to calculate odds ratios and 95% confidence intervals. Because of sample size considerations, we combined data from heterozygotes (N/S) and homozygotes (S/S) for statistical analyses. Detailed pedigree analysis showed that subjects were too distantly related to justify the use of statistical methods to evaluate association between genetic markers and disease status using family-based data.
All urinary oestrogen metabolite laboratory experiments were performed in sextuplet, and the mean value recorded for each subject. To account for any non-normal data distributions, and for consistency of analysis, statistical significance was evaluated using nonparametric statistical techniques with two-sided testing, with a p-value of <0.05 considered statistically significant. Urinary oestrogen metabolite levels were compared between the AMC and UMC groups using the nonparametric Mann–Whitney U-test (Wilcoxon rank-sum test). Oestrogen metabolite values are expressed as mean±sd or median (interquartile range).
All p-values are two-sided. Statistical analysis was performed using the statistical package SPSS for Windows (Version 16.0; SPSS Inc., Chicago, IL, USA).
RESULTS
The distributions of the CYP1B1 N453S genotypes among the 140 subjects tested were compared between AMCs (n = 92) and UMCs (n = 48) according to sex. Among females (n = 86), the genotype frequencies from AMCs (n = 62) were significantly different from the UMCs (n = 24) (Chi-squared 8.06, degrees of freedom (df) 1, p = 0.005), with a greater penetrance among subjects homozygous for the wild-type genotype (N/N) (table 3⇓). Compared with the N/N genotype, the unadjusted odds ratio for the group with N/S or S/S genotypes was 0.248 (95% CI 0.092–0.669), indicating a four-fold higher penetrance among female N/N homozygotes compared with those with N/S or S/S genotypes. There was no difference among males (n = 54) (Chi-squared 0.69, df 1, p = 0.407) between AMCs and UMCs.
The nested case–control study evaluated urinary metabolites in 11 nonsmoking females from the genotype study. The five AMCs (age 52.4±15.3 yrs) had a mean body mass index (BMI) of 23.7 kg·m−2, and three out of five were post-menopausal. The six UMCs (age 57.3±10.8 yrs) were healthy nonsmokers with no known diseases and taking no medications, with a mean BMI of 25.4 kg·m−2, and five out of six were post-menopausal.
Cases (AMCs) and controls (UMCs) were not statistically different in terms of age (p = 0.463), weight and height (BMI; p = 0.268), age at menarche, number of pregnancies, parity or menstrual status. Urine creatinine levels were not statistically different between cases and controls (0.84 versus 1.16 mg·mL−1; p = 0.361)
There were no statistically significant differences in oestrogen metabolite concentrations, although the sum of urinary metabolites (2-OHE and 16α-OHE1) measured was higher in the AMCs than UMCs, as was the measured level of 16α-OHE1. Importantly, the 2-OHE/16α-OHE1 ratio was 2.3-fold lower in AMCs compared with UMCs (0.65 versus 1.48 ng per mg creatinine per mL; p = 0.006), with a distinct separation in the ranges of values between the two groups (table 4⇓ and fig. 1⇓). These findings were unchanged by the stratification into post-menopausal subjects (0.57 (n = 3) versus 1.47 (n = 4) ng per mg creatinine per mL; p = 0.034) and pre-menopausal subjects (0.76 (n = 2) versus 1.83 (n = 1) ng per mg creatinine per mL; p = 0.221).
DISCUSSION
The 2.5-fold female to male disparity in penetrance, and our recent discovery that transcript levels of a major oestrogen-metabolising enzyme (CYP1B1) may modify disease risk in females, suggests that oestrogens or oestrogen metabolites may influence PAH risk in BMPR2 mutation carriers 5. To pursue this further, we determined the genotype of a CYP1B1 polymorphism previously implicated in breast, endometrial and prostate cancers, and also measured oestrogen metabolites in AMCs and UMCs. We found that AMCs are more likely to be CYP1B1 N/N than UMCs, and to have lower 2-OHE/16α-OHE1 ratios in urinary oestrogen metabolites consistent with diminished CYP1B1 activity. Our results suggest that oestrogen or oestrogen metabolites participate in the pathogenesis of FPAH, and provide an important clue toward understanding the mechanisms behind sex differences in disease penetrance.
While FPAH is an autosomal dominant disease, the female predominance (∼2.5 females to one male), reduced penetrance and variable age of diagnosis strongly suggest that combinations of modifying exposures may interact to promote its pathogenesis 27. However, sex remains the most definitive modifier of disease, suggesting that oestrogens may be important. Additional features that also implicate oestrogens include the association of pharmacological oestrogen exposures with PAH onset and the onset of PAH around pregnancy 2. Furthermore, data from our research registry also demonstrates that fewer of the female diagnoses and deaths occur during childhood compared with males (13% of females versus 21% of males, p = 0.031 for difference prior to age 18 yrs; fig. 2⇓), but the reverse for adults. The smaller difference between adult female and male groups likely reflects the importance of additional factors in the pathophysiology of disease as a BMPR2 mutation carrier ages.
Experimental data concerning the role of oestrogens in pulmonary vascular disease implicate oestrogens as promoters of inappropriate vascular proliferation and cell damage, but also as inhibitors of pulmonary vasoreactivity. In vitro, oestrogen is mitogenic and promotes proliferation of pulmonary vascular smooth muscle cells, while the anti-oestrogen drug tamoxifen blocks this effect 4. In chronic hypoxia models oestrogen attenuates pulmonary artery vasoreactivity and acute hypoxic pulmonary vasoconstriction 28. The implication of oestrogens in the pathogenesis of BMPR2-associated FPAH is consistent with hypotheses that PAH is primarily driven by pro-proliferative and anti-apoptotic forces within the pulmonary vascular bed 29.
The CYP multigene family of enzymes metabolise many endogenous and exogenous compounds 30. CYP1B1 is a key enzyme in oestrogen metabolism associated with several hormone-regulated cancers, including breast, endometrial and prostate. It is an active extrahepatic enzyme that is highly expressed in oestrogen-regulated tissues, such as the breast, but also in the lung, suggesting that CYP1B1 is important in the localised metabolic control of oestrogen balance 7. Interestingly, CYP1B1 is transcriptionally activated by oestrogens and oestrogen metabolites, via the oestrogen receptor, to hydroxylate oestrogens into the 2-OHE and 4-OHE metabolites 31. Thus, oestrogens are not only substrates, but activators of CYP1B1 activity by direct interaction between the oestrogen receptor and oestrogen response elements on the CYP1B1 gene.
Oestrogens and most oestrogen metabolites are pro-proliferative. However, not all oestrogen metabolites are equal in their effects. In particular, 4-OHE and 16α-OHE1 exert far greater mitogenic and genotoxic effects than 2-OHE, making the ratio of 2-OHE to 4-OHE and/or 16α-OHE1 important in the genesis of hormonal-based cancers 32–36. For example, it has been hypothesised that metabolism favouring 2-OHE over the 16α-OHE1 pathway may be protective of breast cancer in at-risk subjects 18, 37. The substantially lower ratio of 2-OHE/16α-OHE1 we observed in AMCs compared with UMCs is consistent with the expectation that a lower ratio is detrimental to those subjects already at increased genetic risk of developing PAH.
Consistent with this finding, 2-methoxyoestradiol and its synthetic analogue 2-ethoxyoestradiol possess potent anti-mitogenic properties. These compounds, similar to 2-OHE, are metabolites of oestrogen without oestrogenic effects that may antagonise oestrogen. In fact, in a rat model of monocrotaline-induced pulmonary hypertension, treatment with both 2-methoxyoestradiol and 2-ethoxyoestradiol prevented pulmonary hypertension exacerbation and reduced mortality via attenuation of the inflammatory response and vascular remodelling 38, 39. The findings in our study would support the conclusion than certain non-oestrogenic metabolites of oestrogen may be protective, and perhaps of therapeutic benefit, in pulmonary hypertension.
While the present study is the first, to our knowledge, to suggest the possible contribution of CYP1B1 N453S genotypes and oestrogen metabolites to PAH, our results are consistent with those of oncology researchers who have implicated this polymorphism in several oestrogen-regulated cancers. For example, studying individuals at risk of familial breast cancer, Greenlee et al. 40 found 40% higher 2-OHE/16α-OHE1 ratios in subjects heterozygous (N/S) or homozygous (S/S) compared with those with the homozygous (N/N) genotype. Our findings of increased penetrance (higher percentage of AMCs) among N/N homozygotes, as well as lower 2-OHE/16α-OHE1 ratios in AMCs compared with UMCs, suggest that variations in oestrogen metabolism due to CYP1B1 genetic variation may influence disease in FPAH as well.
Our study has several strengths, including the availability of specimens from one of the largest registries of FPAH in the world and standardised methods of sample collection. However, several limitations in our study warrant caution for interpreting its findings, and suggest the importance of further investigation. First, within the limits of exclusion by family history, the majority of subjects included in our analysis are unrelated. While a small percentage are distantly related to one another, this population of subjects did not meet recommended criteria for the use of family-based tests of genetic analysis, which may have improved power to detect differences by genotypes 41. Secondly, our hypothesis that oestrogens and oestrogen metabolites are important in PAH makes it likely that additional variations in CYP1B1 and/or other enzymes involved in oestrogen metabolism may also be important. This is particularly relevant because the formation and metabolism of oestrogen is a complicated process, involving multiple enzymes in addition to CYP1B1, including aromatase, CYP1A1, CYP3A4, and catechol-O-methyltransferase. In fact, pulmonary hypertension associated with liver disease, which also predominantly occurs in females, has recently been linked to alterations in oestrogen activity associated with polymorphisms in the gene that encodes the aromatase enzyme 42. Thirdly, although we attempted to recruit all AMC and UMC members of our PAH research registry who met study criteria, and the number of subjects available for this study represents a significant fraction of all known FPAH cases in North America, the number remains small. Although the large 2-OHE/16α-OHE1 ratio difference we observed between AMCs and UMCs was statistically significant despite the small number of subjects, we need to expand our study to include more subjects. Finally, we studied oestrogen metabolites previously reported to be important in other diseases with sex disparity, but did not measure other potentially important oestrogen metabolites which may impact risk.
In conclusion, our results support the hypothesis that variations in oestrogen metabolism may play a detrimental role in the development of FPAH, and that the risk of FPAH is influenced by variations in CYP1B1. Female BMPR2 mutation carriers with the N/N genotype for CYP1B1 N453S have greater penetrance of FPAH, consistent with prior evidence that this genotype is associated with a lower 2-OHE/16α-OHE1 ratio and, thus, more mitogenic and genotoxic activity. As expected, female BMPR2 AMCs have statistically lower 2-OH/16α-OHE1 ratios than female BMPR2 UMCs. In aggregate, our findings support the possibility that oestrogens play a direct and detrimental role in the development of FPAH, and that the risk of FPAH is affected by variations in CYP1B1.
Support statement
Sources of support: National Institutes of Health, National Heart, Lung, and Blood Institute P01 HL72058, K12 RR1 7697 and P01 HL072058 (Bethesda, MD, USA); and the Vanderbilt General Clinical Research Center grant M01 RR 000095 NCRR/NIH (Nashville, TN, USA).
Statement of interest
None declared.
Acknowledgments
The authors thank the many patients and families who graciously contributed to this work.
- Received January 21, 2009.
- Accepted March 23, 2009.
- © ERS Journals Ltd