- Research article
- Open Access
- Open Peer Review
CYP19A1 promoter methylation in saliva associated with milestones of pubertal timing in urban girls
BMC Pediatricsvolume 14, Article number: 78 (2014)
Childhood obesity and early puberty are intermediate risk factors for later metabolic and reproductive disorders including diabetes, polycystic ovarian syndrome (PCOS), and breast cancer. Atypical methylation patterns in genes related to hormone and adipose metabolism, such as CYP19A1 (aromatase) and PPARG (peroxisome proliferator-activated receptor gamma), are associated with alterations in gene expression which may contribute to pathogenesis of these diseases. If present in early life, it is conceivable similar methylation aberrations may result in hormone perturbations that alter pubertal timing.
We used Cox proportional hazard models to investigate whether promoter methylation of CYP19A1 and PPARG, independently or in concert with body weight, was associated with age at breast (B2) or pubic hair development (PH2) when assayed in saliva DNA collected from a cohort of New York City, Black and Hispanic girls (N = 130) enrolled in a study of pubertal timing between 6–8 years of age.
An inverse association between CYP19A1 methylation and risk of early PH2 was suggested (HR = 0.95, 95% CI = 0.90-1.00, p = 0.05). CYP19A1 methylation also appeared to modify risk of early B2 associated with body weight. Specifically, compared to normal weight girls with ‘high’ CYP19A1 methylation, significantly increased risk of early B2 was observed in overweight girls with ‘low’ but not ‘high’ CYP19A1 methylation (HR = 2.15; 95% CI = 1.23- 3.76). However, in formal tests for effect modification, the interaction between body weight and methylation did not reach statistical significance (p for interaction = 0.085). PPARG methylation was not significantly associated with PH2 or B2.
Though limited by sample size, our findings suggest methylation of CYP19A1, a critical gene in estrogen biosynthesis, may influence timing of breast development in overweight girls. Consistent with emerging reports, these data support the notion that epigenetic marks in surrogate tissues may improve risk prediction when added to standard plasma and anthropometric indicators, and warrant further study.
Emerging evidence indicates modifiable lifestyle factors and time-stable epigenetic determinants influence childhood obesity [1, 2] and timing of puberty . Childhood obesity and early puberty are risk factors for several metabolic and reproductive disorders, including: early menarche , adult obesity , diabetes , polycystic ovarian syndrome (PCOS) , and breast cancer . PCOS and its common childhood antecedent, premature pubarche, are associated with hyperinsulinemia and androgen excess in overweight females [7, 9]. Similarly, plasma leptin  and local adiposal estrogen  are associated with breast tumor growth in obese postmenopausal cases, and are hypothesized to drive early breast development and menarche in overweight girls [12–15]. Given that the ovaries are largely quiescent pre-menarche and post-menopause, it is conceivable biomarkers that link adiposity to pubertal timing may inform risk of androgen and estrogen-related disease development throughout the life course. In this exploratory pilot study we asked whether promoter methylation marks in PPARG or CYP19A1, two genes that connect energy balance to lifetime estrogen exposure, are associated with pubertal development in girls.
The transcription factor PPARγ is a master regulator of adipose differentiation and endocrine function. Human and animal data link PPARG hypermethylation to reduced PPARγ expression that is observed in cases of diabetes , breast cancer , and hyperandrogenic PCOS . Because PCOS and early pubarche share hyperinsulemia and androgen excess as common clinical features, it is conceivable PPARG methylation aberrations may alter timing of pubarche if present in early development .
Aromatase (product of the CYP19A1 gene) catalyzes all bodily estrogen biosynthesis via aromatization of androgen precursors, and inhibitors of this enzyme are highly effective therapies for estrogen-sensitive cancers of the breast . Aromatase expression varies across tissues and individuals owing to differential activation and repression of several tissue-specific gene promoters. Each of these promoters regulates a unique untranslated ‘first exon’ (‘exons I’) that is spliced to the common coding exons II-X and can be mapped back to its promoters for purposes of deducing what factors drive transcription in a particular tissue . In healthy adipose, CYP19A1 is expressed in fibroblasts primarily from activation of a distal glucocorticoid-regulated promoter termed ‘pI.4’, and is lost in the path of adipocyte differentiation driven by PPARγ. In malignant breast biopsies, CYP19A1 is overexpressed 3–4 fold in tumor and proximal adipose tissue from several tissue-specific promoters , but primarily from two largely ‘gonad-specific’ cAMP-responsive promoters termed “pII/I.3”  that are activated by cancer-associated transcription factors PPARγ suppresses in healthy breast [17, 23]. Though this phenomenon of increased tissue-specific promoter usage with preference for gonadal promoters (termed ‘promoter switching’) is well-documented in the breast cancer literature [19, 24], its timing and extent of tissue distribution related to intermediate risk indicators is less understood. Demura and Bulun  recently described hypomethylation of a CpG dinucleotide in the ‘cAMP-response element-like sequence’ (CLS) of pII/I.3, which they detected in aromatase-overexpressing fibroblasts derived from a skin-punch biopsy of a healthy patient. In light of the aggregate of findings regarding CYP19A1 misexpression from gonadal promoters in breast cancer cases, Demura and Bulun postulated CYP19A1 pII/I.3 hypomethyaltion may contribute to the phenomenon of ‘promoter switching’ and inter-individual variability in lifetime estrogen exposure.
In the present study we sought to determine whether methylation of this CYP19A1 pII/I.3 locus  or the average of five CpG dinucleotides in a differentially methylated region of the PPARG promoter [9, 16, 26] was associated with timing of pubic hair (PH2) or breast development (B2) in a cohort of New York City, Black and Hispanic girls who were enrolled in a study of pubertal timing between 6–8 years of age.
Growing up healthy
Prospective cohort study and part of the Puberty Studies of the Breast Cancer and Environment Research Program (BCERP) . The overarching goal of this longitudinal investigation is to identify genetic and environmental risk factors related to altered timing of puberty onset in girls. Girls, 6 to 8 years of age, from East Harlem schools, community health centers, and the Mount Sinai Pediatric clinic were recruited for this study between 2004–2007, as previously reported . Consent was obtained from parents or guardians and child assent was independently verified; the study was approved by the Institution Review Board at Mount Sinai. Eligibility included age, female sex, no underlying endocrine medical conditions, and self-identified Black or Hispanic race/ethnicity. A total of 416 girls were enrolled in the study at baseline; we restricted our analysis to the 130 who had whole saliva collected. The distribution of major demographic and physiological variables, including race, caregiver’s education, BMI, breast stage and pubic stage, showed no significant difference between those who donated saliva samples and those who did not (data not shown).
Demographic and anthropometric data collection
Uniformly trained interviewers conducted annual in-person interviews and standardized anthropometric measurements. Annual pubertal stage assessments were performed by physicians or nurse practitioners according to BCERP consortium-standardized protocols; the principal endpoints were age at first pubic hair and breast development as described in detail previously . A structured questionnaire was administered to the girl’s parent or guardian in either English or Spanish. Information ascertained through the questionnaire included medical history and demographic variables. Body mass index (BMI) was calculated as weight (in kg) divided by height (in cm)-squared. We classified girls as ‘overweight’ according to Centers for Disease Control (CDC) and Prevention criteria, where overweight girls had a BMI at or above the 85th percentile of their age and sex-specific BMI distribution. Age at B2 was defined as the midpoint between the age at the last visit where the girl was staged B1 with no prior staging greater than B1 and the age at the first visit where the girl was staged B2 with no subsequent staging less than B2. Girls who entered the study at B2 were assigned age at B2 as six months prior, and girls with a breast stage less than B2 at their last visit were right censored. Age at PH2 was assigned in the same manner.
Saliva DNA collection and processing
Interviewers instructed study participants to deposit saliva in pre-barcoded 2 ml Oragene DNA Self-Collection Kit tubes (DNA Genotek, Canada; REF OG-100) according to the manufacturer’s instructions. Barcoded vials were logged in our database by scanning upon receipt in the laboratory. DNA was extracted from whole saliva collected in Oragene tubes (DNA Genotek) with the ITprep kit (DNA Genotek) according to the manufacturer’s instructions.
Methylation assessment by pyrosequencing
Genomic DNA was bisulfite-converted with the Epitect DNA kit (Qiagen, Carslbad, CA). Pyrosequencing was performed on CYP19A1 and PPARG PCR products amplified from bisulfite-treated DNA with the Pyromark PCR kit (Qiagen) and the primers listed in Additional file 1. Cycling conditions were 95°C for 15 minutes followed by: 45 cycles of 94°C for 30 seconds; annealing temperature for 30 seconds; 72°C for 30 seconds; with a final extension step of 72°C for 10 minutes. PCR products were sequenced using the Pyromark Q24 system and kit (Qiagen). Percent methylation for each region of interest was quantified using Pyromark Q24 software version 1.0.1 (Qiagen). Genomic coordinates for the promoter regions amplified are included in Additional file 1; coordinates were obtained from the UCSC Genome Browser (http://genome.ucsc.edu), human assembly: February 2009, GRCh37/hg19). Laboratory personnel performing DNA methylation analysis were blinded to subject information.
We examined relationships among methylation and study characteristics with parametric (t-tests) and non-parametric (Mann-Whitey U tests) statistics and multivariate linear regression. Cox proportional hazard models were used to identify associations between DNA methylation and age at PH2 or B2. Interaction was examined by including a group variable that was constructed by pairing the dichotomized methylation (below vs. greater than or equal to their medians) and dichotomized body size (normal weight vs. overweight). All models were adjusted for Hispanic ethnicity, Black race, and caregiver education level. All analyses were performed using SAS (version 9.1.3 for PC; SAS Institute Inc., Cary, NC).
Study population demographics according to CYP19A1 and PPARGmethylation
Study subjects (N = 130) were Black or Hispanic girls living in the East Harlem neighborhood of New York City. Girls were recruited in local clinics and community centers between 2004–2007, and were 6 (72 months) to 8.9 (107 months) years old with a mean age of 7.5 years (89.5 months) at time of enrollment. Based on CDC criteria, 39.2% of our study subjects were considered overweight (≥ 85th percentile of their age and sex-specific BMI distribution) and 25.4% were considered obese (≥ 95th percentile). Of the study subjects’ primary caregivers, 59% had completed high school. Among the 130 whole saliva samples collected, 5 failed the pyrosequencing assay for CYP19A1 and 1 for PPARG, leaving 125 and 129 samples, respectively, with methylation data. CYP19A1 methylation values ranged from 77% to 95% (mean: 88.8% ± 3.5% SD). PPARG methylation ranged from 5.6% to 19% (mean: 10.4% ± 2.1% SD). Associations between methylation levels and key demographic variables are summarized in Table 1. No significant differences were observed with respect to race, ethnicity, BMI percentile, or caregiver’s education level.
Gene methylation related to milestones of pubertal development
We investigated whether methylation of CYP19A1 or PPARG was related to age at B2 or PH2 using Cox Proportion Hazards Models (Table 2). For PH2, we observed an inverse association with CYP19A1 methylation in unadjusted models; for a one percent increase in CYP19A1 methylation, girls were 5% more likely to be older at PH2 (HR = 0.95; 95% CI = 0.90-0.99). This association was attenuated in models adjusted for ethnicity, BMI percentile, and caregiver’s education (HR = 0.95, 95% CI = 0.90-1.00). Conversely, no significant associations between age at B2 and CYP19A1 methylation were observed. In addition, no significant associations among PPARG methylation and PH2 or B2 were observed.
Effect of body size modified by gene methylation
Obesity is one of the strongest predictors of pubertal onset . Therefore we next sought to determine whether gene methylation modifies the relationship between BMI and age at PH2 and B2. We created ‘normal weight’ and ‘overweight’ categories of body size (below and greater than or equal to the 85th percentile for age and sex-specific BMI distribution), and ‘high’ and ‘low’ methylation (above and below the median). As shown in Table 3, compared to normal weight girls with high CYP19A1 methylation (low risk referent), risk of earlier breast development was greater among overweight girls with low CYP19A1 methylation (HR = 2.15, 95% CI = 1.23 - 3.76). This BMI-methylation interaction reached borderline significance in formal tests for effect modification (p for interaction = 0.085). A similar effect was observed for CYP19A1 methylation and age at PH2, although the interaction did not reach statistical significance (p for interaction = 0.21). Lastly, no significant interactions between BMI and PPARG methylation in relation to PH2 or B2 were detected.
Early breast and pubic hair development have been associated with disordered leptin, insulin, and IGF-1 profiles in overweight girls in numerous studies . Perturbations in estrogens and androgens, critical drivers of breast and pubic hair development, remain clinically more challenging to detect . Given national trends, there is great motivation to identify biomarkers that add value to current plasma and anthropometric measures used in predicting puberty onset . In this exploratory study we aimed to ascertain whether salivary methylation of the CYP19A1 and PPARG promoters was related to age at breast or pubic hair development in girls, both independently and in concert with body size. In light of the current literature, we anticipated overweight girls with CYP19A1 hypomethylation and PPARG hypermethylation might be predisposed to early breast development [33–35], and those with PPARG hypermethylation to early pubic hair development [9, 16].
Our main observations were that relative hypomethylation of a CpG in the gonadal CYP19A1 promoter termed “pII “ was associated with earlier age at B2 among overweight girls only (Table 3), and with earlier age at PH2 independent of body size (Table 2). While only correlative and based on a relatively small number of samples, our B2 findings are supported by a case report authored by Demura and Bulun , which describes hypomethylation of pI.3/II in CYP19A1-overexpressing fibroblasts relative to CYP19A1-quiescent fibroblasts derived from punch biopsies of four healthy subjects. In their report, CYP19A1 activity was robustly induced in the former upon cAMP stimulation, while fibroblasts from the other three subjects were cAMP-refractory. Further investigation revealed CpG dinuleotides within and proximal to the CLS (CRE-like sequence) of gonadal pI.3/II were relatively hypomethylated in cAMP-responsive CYP19A1-overexpressing fibroblasts, and were relatively hypermethylated in non-responsive fibroblasts. These results support the hypothesis that CYP19A1 hypomethylation may be an early ‘permissive’ event, which renders one susceptible to subsequent intrinsic/extrinsic transcriptional activators of CYP19A1, and concomitant local or systemic estrogen excess. Such a ‘two-hit’ mechanism of derepression (hypomethylation) and activation (e.g., obesity-related cytokines) may also explain why CYP19A1 hypomethylation was associated with early B2 in overweight, but not normal weight girls in the present study.
Aromatase catalyzes estrogen biosynthesis from androgen precursors. Elevated androgen, insulin, and IGF-1 signaling are widely accepted co-determinants of early pubarche in overweight girls [7, 30, 36]. Thus, our finding that CYP19A1 hypomethylation (theoretical increase in expression) was related to earlier age at PH2, independent of BMI (Table 2), was unanticipated. While intriguing, the statistical significance of this association was attenuated after adjustment for covariates, and we can only speculate as to its implications without further study.
We assessed methylation of the −383 to −281 bp region of the PPARG promoter as PPARγ suppresses CYP19A1 expression in breast tissues in culture [20, 33], and relative hypermethylation of this region has been associated with reduced PPARγ expression in hyperandrogenic PCOS  and diabetes models . Though we detected no statistically significant effects related to PPARG methylation in the present study, puberty-associated methylation patterns may exist in genes for PPARγ co-factors, effectors, or downstream targets in salivary or other surrogate tissue DNA. Indeed, methylation biomarkers of childhood adiposity and maternal BMI have been described in RXRA and PPARGC1A when assayed in umbilical tissue [2, 37].
This exploratory investigation has several limitations regarding generalizability, including but not limited to: small sample size, lack of perceived stress assessments, use of candidate genes, and DNA derived from whole saliva samples collected only from Black and Hispanic girls. We describe salivary CYP19A1 hypomethylation not as a ‘causal’ event, but merely as a ‘surrogate biomarker’ that with further study may have utility in predicting risk of premature breast development in overweight girls. Specifically, the CpG we describe is contained in a critical transcription factor binding site (CLS), located in a strong CYP19A1 gonadal promoter termed ‘pII’, which is activated by the ubiquitous pleiotropic second messenger cAMP in the follicular phase of the menstrual cycle . DNA methylation is highly tissue-specific, and CYP19A1 is not likely expressed in buccal epithelial cells from gonadal pII to any significant degree. However, weighing these considerations together supports the notion that hypomethylation of CYP19A1 pII in whole saliva- a tissue in which expression from gonadal pII is likely silenced, may in fact represent a methylation aberration, possibly established early in life. Such a mark could be deemed a ‘surrogate’ aberration if it portends risk generalizable to more functional tissues with niche transcriptional machinery requisite to affect CYP19A1 expression changes that promote disease, as has been extensively reported for pII and estrogen-related disorders [22, 24].
Our findings are only suggestive and only extend to saliva samples we collected from Black and Hispanic girls. Procurement of effector/target tissues (e.g. adipose/breast) to investigate the validity of salivary pII methylation as a risk surrogate in a statistically robust manner in more diverse pediatric populations is precluded by ethical, logistical, and economic considerations. However, animal and cell line co-culture models designed to capture adipose and developing breast tissue interactions are emerging [39, 40], and it will be interesting to follow developments that functionally characterize the complex biological and environmental interactions that orchestrate epigenetic factors related to thelarche and pubarche onset.
Consistent with emerging human and animal studies, our findings suggest methylation of CYP19A1 may influence timing of breast development in overweight girls. These data warrant further investigation, and support the notion that epigenetic biomarkers may one day add value to current plasma and anthropometric measures used in predicting timing of puberty onset.
Breast development Tanner stage 2
Breast Cancer and Environment Research Center
Body mass index
Cyclic adenosine monophosphate
CRE binding protein
CYP19A1 gene (aromatase)
Insulin-like growth factor
National Cancer Institute
National Center for Research Resources
National Institutes of Environmental Health Sciences
Polycystic ovarian syndrome
Pubic hair development Tanner stage 2 stages 2 & higher collapsed in this study
Gene encoding peroxisome proliferator-activated receptor gamma
Trascription start site.
Relton CL, Groom A, St Pourcain B, Sayers AE, Swan DC, Embleton ND, Pearce MS, Ring SM, Northstone K, Tobias JH, Trakalo J, Ness AR, Shaheen SO, Davey Smith G: DNA methylation patterns in cord blood DNA and body size in childhood. PLoS One. 2012, 7 (3): e31821-10.1371/journal.pone.0031821.
Godfrey KM, Sheppard A, Gluckman PD, Lillycrop KA, Burdge GC, McLean C, Rodford J, Slater-Jefferies JL, Garratt E, Crozier SR, Emerald BS, Gale CR, Inskip HM, Cooper C, Hanson MA: Epigenetic gene promoter methylation at birth is associated with child’s later adiposity. Diabetes. 2011, 60 (5): 1528-1534. 10.2337/db10-0979.
Lomniczi A, Loche A, Castellano JM, Ronnekleiv OK, Bosch M, Kaidar G, Knoll JG, Wright H, Pfeifer GP, Ojeda SR: Epigenetic control of female puberty. Nat Neurosci. 2013, 16 (3): 281-289. 10.1038/nn.3319.
Ibanez L, Jimenez R, de Zegher F: Early puberty-menarche after precocious pubarche: relation to prenatal growth. Pediatrics. 2006, 117 (1): 117-121. 10.1542/peds.2005-0664.
Nader PR, O’Brien M, Houts R, Bradley R, Belsky J, Crosnoe R, Friedman S, Mei Z, Susman EJ: National Institute of Child Health and Human Development Early Child Care Research Network: identifying risk for obesity in early childhood. Pediatrics. 2006, 118 (3): e594-e601. 10.1542/peds.2005-2801.
American Diabetes Association: Type 2 diabetes in children and adolescents. Pediatrics. 2000, 105 (3 Pt 1): 671-680.
Ibanez L, Valls C, Marcos MV, Ong K, Dunger DB, De Zegher F: Insulin sensitization for girls with precocious pubarche and with risk for polycystic ovary syndrome: effects of prepubertal initiation and postpubertal discontinuation of metformin treatment. J Clin Endocrinol Metab. 2004, 89 (9): 4331-4337. 10.1210/jc.2004-0463.
Must A, Phillips SM, Naumova EN: Occurrence and timing of childhood overweight and mortality: findings from the Third Harvard Growth Study. J Pediatr. 2012, 160 (5): 743-750. 10.1016/j.jpeds.2011.10.037.
Qu F, Wang FF, Yin R, Ding GL, El-Prince M, Gao Q, Shi BW, Pan HH, Huang YT, Jin M, Leung PC, Sheng JZ, Huang HF: A molecular mechanism underlying ovarian dysfunction of polycystic ovary syndrome: hyperandrogenism induces epigenetic alterations in the granulosa cells. J Mol Med (Berl). 2012, 90 (8): 911-923. 10.1007/s00109-012-0881-4.
Tessitore L, Vizio B, Pesola D, Cecchini F, Mussa A, Argiles JM, Benedetto C: Adipocyte expression and circulating levels of leptin increase in both gynaecological and breast cancer patients. Int J Oncol. 2004, 24 (6): 1529-1535.
Geisler J: Breast cancer tissue estrogens and their manipulation with aromatase inhibitors and inactivators. J Steroid Biochem Mol Biol. 2003, 86 (3–5): 245-253.
Smith JT: Sex steroid regulation of kisspeptin circuits. Adv Exp Med Biol. 2013, 784: 275-295. 10.1007/978-1-4614-6199-9_13.
Jasik CB, Lustig RH: Adolescent obesity and puberty: the “perfect storm”. Ann N Y Acad Sci. 2008, 1135: 265-279. 10.1196/annals.1429.009.
Biro FM, Greenspan LC, Galvez MP: Puberty in girls of the 21st century. J Pediatr Adolesc Gynecol. 2012, 25 (5): 289-294. 10.1016/j.jpag.2012.05.009.
Hilakivi-Clarke L, Clarke R, Onojafe I, Raygada M, Cho E, Lippman M: A maternal diet high in n - 6 polyunsaturated fats alters mammary gland development, puberty onset, and breast cancer risk among female rat offspring. Proc Natl Acad Sci U S A. 1997, 94 (17): 9372-9377. 10.1073/pnas.94.17.9372.
Fujiki K, Kano F, Shiota K, Murata M: Expression of the peroxisome proliferator activated receptor gamma gene is repressed by DNA methylation in visceral adipose tissue of mouse models of diabetes. BMC Biol. 2009, 38: 7007-7-38
Meng L, Zhou J, Sasano H, Suzuki T, Zeitoun KM, Bulun SE: Tumor necrosis factor alpha and interleukin 11 secreted by malignant breast epithelial cells inhibit adipocyte differentiation by selectively down-regulating CCAAT/enhancer binding protein alpha and peroxisome proliferator-activated receptor gamma: mechanism of desmoplastic reaction. Cancer Res. 2001, 61 (5): 2250-2255.
Santen RJ, Brodie H, Simpson ER, Siiteri PK, Brodie A: History of aromatase: saga of an important biological mediator and therapeutic target. Endocr Rev. 2009, 30 (4): 343-375. 10.1210/er.2008-0016.
Bulun SE, Chen D, Lu M, Zhao H, Cheng Y, Demura M, Yilmaz B, Martin R, Utsunomiya H, Thung S, Su E, Marsh E, Hakim A, Yin P, Ishikawa H, Amin S, Imir G, Gurates B, Attar E, Reierstad S, Innes J, Lin Z: Aromatase excess in cancers of breast, endometrium and ovary. J Steroid Biochem Mol Biol. 2007, 106 (1–5): 81-96.
Rubin GL, Duong JH, Clyne CD, Speed CJ, Murata Y, Gong C, Simpson ER: Ligands for the peroxisomal proliferator-activated receptor gamma and the retinoid X receptor inhibit aromatase cytochrome P450 (CYP19) expression mediated by promoter II in human breast adipose. Endocrinology. 2002, 143 (8): 2863-2871.
Bulun SE, Price TM, Aitken J, Mahendroo MS, Simpson ER: A link between breast cancer and local estrogen biosynthesis suggested by quantification of breast adipose tissue aromatase cytochrome P450 transcripts using competitive polymerase chain reaction after reverse transcription. J Clin Endocrinol Metab. 1993, 77 (6): 1622-1628.
Agarwal VR, Bulun SE, Leitch M, Rohrich R, Simpson ER: Use of alternative promoters to express the aromatase cytochrome P450 (CYP19) gene in breast adipose tissues of cancer-free and breast cancer patients. J Clin Endocrinol Metab. 1996, 81 (11): 3843-3849.
Clyne CD, Speed CJ, Zhou J, Simpson ER: Liver receptor homologue-1 (LRH-1) regulates expression of aromatase in preadipocytes. J Biol Chem. 2002, 277 (23): 20591-20597. 10.1074/jbc.M201117200.
Bulun SE, Lin Z, Imir G, Amin S, Demura M, Yilmaz B, Martin R, Utsunomiya H, Thung S, Gurates B, Tamura M, Langoi D, Deb S: Regulation of aromatase expression in estrogen-responsive breast and uterine disease: from bench to treatment. Pharmacol Rev. 2005, 57 (3): 359-383. 10.1124/pr.57.3.6.
Demura M, Bulun SE: CpG dinucleotide methylation of the CYP19 I.3/II promoter modulates cAMP-stimulated aromatase activity. Mol Cell Endocrinol. 2008, 283 (1–2): 127-132.
Pancione M, Sabatino L, Fucci A, Carafa V, Nebbioso A, Forte N, Febbraro A, Parente D, Ambrosino C, Normanno N, Altucci L, Colantuoni V: Epigenetic silencing of peroxisome proliferator-activated receptor gamma is a biomarker for colorectal cancer progression and adverse patients’ outcome. PLoS One. 2010, 5 (12): e14229-10.1371/journal.pone.0014229.
Biro FM, Galvez MP, Greenspan LC, Succop PA, Vangeepuram N, Pinney SM, Teitelbaum S, Windham GC, Kushi LH, Wolff MS: Pubertal assessment method and baseline characteristics in a mixed longitudinal study of girls. Pediatrics. 2010, 126 (3): e583-e590. 10.1542/peds.2009-3079.
Wolff MS, Teitelbaum SL, Pinney SM, Windham G, Liao L, Biro F, Kushi LH, Erdmann C, Hiatt RA, Rybak ME, Calafat AM: Investigation of Relationships between Urinary Biomarkers of Phytoestrogens, Phthalates, and Phenols and Pubertal Stages in Girls. Environ Health Perspect. 2010, 118 (7): 1039-1046. 10.1289/ehp.0901690.
Davison KK, Susman EJ, Birch LL: Percent body fat at age 5 predicts earlier pubertal development among girls at age 9. Pediatrics. 2003, 111 (4 Pt 1): 815-821.
Ibanez L, Lopez-Bermejo A, Diaz M, Marcos MV, de Zegher F: Metformin treatment for four years to reduce total and visceral fat in low birth weight girls with precocious pubarche. J Clin Endocrinol Metab. 2008, 93 (5): 1841-1845. 10.1210/jc.2008-0013.
Janfaza M, Sherman TI, Larmore KA, Brown-Dawson J, Klein KO: Estradiol levels and secretory dynamics in normal girls and boys as determined by an ultrasensitive bioassay: a 10 year experience. J Pediatr Endocrinol Metab. 2006, 19 (7): 901-909.
Landrigan P, Garg A, Droller DB: Assessing the effects of endocrine disruptors in the National Children’s Study. Environ Health Perspect. 2003, 111 (13): 1678-1682. 10.1289/ehp.5799.
Subbaramaiah K, Howe LR, Zhou XK, Yang P, Hudis CA, Kopelovich L, Dannenberg AJ: Pioglitazone, a PPARgamma agonist, suppresses CYP19 transcription: evidence for involvement of 15-hydroxyprostaglandin dehydrogenase and BRCA1. Cancer Prev Res (Phila). 2012, 5 (10): 1183-1194. 10.1158/1940-6207.CAPR-12-0201.
Subbaramaiah K, Morris PG, Zhou XK, Morrow M, Du B, Giri D, Kopelovich L, Hudis CA, Dannenberg AJ: Increased levels of COX-2 and prostaglandin E2 contribute to elevated aromatase expression in inflamed breast tissue of obese women. Cancer Discov. 2012, 2 (4): 356-365. 10.1158/2159-8290.CD-11-0241.
Morris PG, Hudis CA, Giri D, Morrow M, Falcone DJ, Zhou XK, Du B, Brogi E, Crawford CB, Kopelovich L, Subbaramaiah K, Dannenberg AJ: Inflammation and increased aromatase expression occur in the breast tissue of obese women with breast cancer. Cancer Prev Res (Phila). 2011, 4 (7): 1021-1029. 10.1158/1940-6207.CAPR-11-0110.
Bremer AA: Polycystic ovary syndrome in the pediatric population. Metab Syndr Relat Disord. 2010, 8 (5): 375-394. 10.1089/met.2010.0039.
Gemma C, Sookoian S, Alvarinas J, Garcia SI, Quintana L, Kanevsky D, Gonzalez CD, Pirola CJ: Maternal pregestational BMI is associated with methylation of the PPARGC1A promoter in newborns. Obesity (Silver Spring). 2009, 17 (5): 1032-1039. 10.1038/oby.2008.605.
Michael MD, Michael LF, Simpson ER: A CRE-like sequence that binds CREB and contributes to cAMP-dependent regulation of the proximal promoter of the human aromatase P450 (CYP19) gene. Mol Cell Endocrinol. 1997, 134 (2): 147-156. 10.1016/S0303-7207(97)00178-0.
Liu E, Samad F, Mueller BM: Local adipocytes enable estrogen-dependent breast cancer growth: Role of leptin and aromatase. Adipocyte. 2013, 2 (3): 0-4.
Bulun SE, Chen D, Moy I, Brooks DC, Zhao H: Aromatase, breast cancer and obesity: a complex interaction. Trends Endocrinol Metab. 2012, 23 (2): 83-89. 10.1016/j.tem.2011.10.003.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2431/14/78/prepub
This project was made possible by the Breast Cancer and the Environment Research Program (BCERP) award numbers U01ES012771, U01ES019454 and U01ES019457 from the National Institute of Environmental Health Sciences (NIEHS) and the National Cancer Institute (NCI), and CSTA-UL1RR029887 from the National Center for Research Resources (NCRR). Its contents are solely the responsibility of the authors and do not necessarily represent the official views of the NIEHS or NCI, the National Institutes of Health. TRS was supported by a training grant from the National Institute of Environmental Health Sciences (NIEHS, 5T32ES013678-07) while writing the manuscript.
All authors declare that they have no competing interests.
MSW and SLT are the co-PIs of the parent study. TRS and JC contributed to the design of this sub-study. TRS performed pyrosequencing analysis, all related laboratory procedures, and assumes the primary role of composing the manuscript with the supervision of JC. SLT and AP performed statistical analysis. All authors were responsible for critical revision of the manuscript and approved the final copy.