Email updates

Keep up to date with the latest news and content from BMC Genomics and BioMed Central.

Open Access Highly Accessed Research article

The plasticizer butyl benzyl phthalate induces genomic changes in rat mammary gland after neonatal/prepubertal exposure

Raquel Moral1, Richard Wang1, Irma H Russo1, Daniel A Mailo1, Coral A Lamartiniere2 and Jose Russo1*

Author Affiliations

1 Breast Cancer Research Laboratory, Fox Chase Cancer Center, Philadelphia, PA 19111, USA

2 Department of Pharmacology and Toxicology, University of Alabama at Birmingham, Birmingham, AL 35294, USA

For all author emails, please log on.

BMC Genomics 2007, 8:453  doi:10.1186/1471-2164-8-453

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2164/8/453


Received:24 April 2007
Accepted:6 December 2007
Published:6 December 2007

© 2007 Moral et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Phthalate esters like n-butyl benzyl phthalate (BBP) are widely used plasticizers. BBP has shown endocrine-disrupting properties, thus having a potential effect on hormone-sensitive tissues. The aim of this study is to determine the effect of neonatal/prepubertal exposure (post-natal days 2–20) to BBP on maturation parameters and on the morphology, proliferative index and genomic signature of the rat mammary gland at different ages of development (21, 35, 50 and 100 days).

Results

Here we show that exposure to BBP increased the uterine weight/body weight ratio at 21 days and decreased the body weight at time of vaginal opening. BBP did not induce significant changes on the morphology of the mammary gland, but increased proliferative index in terminal end buds at 35 days and in lobules 1 at several ages. Moreover, BBP had an effect on the genomic profile of the mammary gland mainly at the end of the exposure (21 days), becoming less prominent thereafter. By this age a significant number of genes related to proliferation and differentiation, communication and signal transduction were up-regulated in the glands of the exposed animals.

Conclusion

These results suggest that BBP has an effect in the gene expression profile of the mammary gland.

Background

Butyl benzyl phthalate (BBP) is a plasticizer commonly used in pipes, vinyl floor tiles, vinyl foams, and carpet backing, and to a minor extent, in cellulose plastics and polyurethane [1]. This compound has shown to be weakly estrogenic [2], and to induce adverse effects on the development of male reproductive system [3-5], thus acting as an endocrine disruptor. Adult BBP intake has been estimated at 2 μg/kg body weight/day and it has been stated that exposure to infants and children could be up to three-fold higher [1]. Environmental endocrine disruptors are of growing concern due to the potential impact on human health. Actually, BBP is one of the phthalates selected by the National Toxicology Program (NTP) Center for the Evaluation of Risks to Human Reproduction (CERHR), because of its high production volume, potential for human exposure from their widespread use and occurrence within the environment, concern to the public, and published evidence of reproductive or developmental toxicity [6]. The evaluation concluded that there was negligible concern for adverse reproductive effects in exposed men, but data was insufficient to reach conclusions on possible reproductive effects in exposed women [6]. Nevertheless, the European Commission has banned BBP from toys and child care items at concentrations of greater than 0,1% by mass of the plasticized material [7] and new data is providing demonstration of possible developmental effects in humans [8]. Urinary concentrations of phthalate metabolites in pregnant mothers, including the major BBP metabolite monobenzyl phthalate, have been inversely correlated to anogenital distance among male infants [8]. Moreover an etiological association of serum BBP concentration with endometriosis has been suggested, as severity of endometriosis was strongly correlated with BBP concentration in blood [9]. Although more data is needed regarding the possible effects of this compound on human health, experimental studies have given evidence of an antiandrogenic effect of BBP and monobenzyl phthalate. Prenatal and prepubertal exposure of rats to such compounds induced significant alterations on reproductive system of males, like undescended testes, decrease in the anogenital distance and other malformations of the external genitalia, sex accessory glands, epididymis and testes [3-5]. It has also been reported a decrease in the weight of the ovaries in females exposed during adulthood, as well as changes in the offspring of exposed rats, including a decrease of the body weight of male and female at birth and an increase in the anogenital distance in females [10]. All these changes are compatible with an endocrine disrupting action of BBP.

Endocrine disruptors have also been suggested to contribute to the development of hormone-dependent cancers [11]. Although the risk for cancer is multifactorial, there is substantial contribution of environmental factors, including diet and environmental chemicals. Breast cancer is an estrogen-dependent malignancy with a significant mortality in women worldwide [12]. Changes in the hormonal environment during critical stages of development can modify the architecture and biological characteristics of the mammary gland and thus affect the future susceptibility to develop breast cancer [13,14]. The effects of environmental endocrine disruptors seem to be dependent on the time when the exposure occurs. Diethylstilbestrol has been associated to increased or decreased susceptibility for chemically-induced mammary cancer in rats depending on the period of exposure [15,16].

The present study is designed to determine the effects of neonatal and prepubertal exposure to BBP on rat mammary gland at different stages of development using morphological parameters of gland differentiation [13,14], proliferative index and the genomic expression profile. Here we report that BBP has an effect on the gene expression profile of the mammary gland.

Results

Maturation parameters

Neonatal/prepubertal exposure to BBP did not affect body weight, uterine weight or day at vaginal opening. The uterine to body weight ratios were significantly increased in the exposed animals at 21 days of age (Table 1). BBP exposure did not modify the day at vaginal opening (control group: 30.46 ± 0.13 days; BBP group: 30.44 ± 0.11 days), but decreased significantly (p < 0.05) the body weight at that day (control group: 101.96 ± 1.09 g; BBP group: 97.26 ± 1.49 g).

Table 1. Parameters for post-natal maturation and development

Mammary gland architecture

Histoarchitectural analysis of the mammary gland was performed by determining the number of specific epithelial structures at 21, 35, 50 and 100 days. The number of undifferentiated TEB decreased from 21 to 100 days in both control and BBP treated animals, while the TDs and Lob1 increased over time. Neonatal/prepubertal exposure to BBP did not modify significantly the number of these epithelial structures when compared to control groups (Table 1).

Proliferative index of the mammary gland

BrdU incorporation into the DNA was used as an index of cell proliferation (Figure 1). TEB from BBP exposed rats had significantly (p < 0.05) higher proliferation index than controls by 35 days, whereas in TD and AB the differences were not significant. The proliferative index of Lob1 from controls increased from 35 to 50 days and decreased at 100 days, while in BBP exposed rats this structure showed maximum index of proliferation at 35 days and decreased thereafter, showing significant differences when compared to controls at all these ages (Figure 1).

thumbnailFigure 1. Proliferative index analyses of the mammary gland. A: Immunohistochemical detection of BrdU incorporation in proliferating cells (brown cells) in the epithelial structures TEB, TD and ducts, AB and Lob1. Olympus BX40 microscope with 40× objective. B: Proliferative index (mean ± SEM) in each epithelial structure at 21, 35, 50 and 100 days of age. *: significantly different compared to control group (p < 0.05).

Gene expression analysis by microarrays

Neonatal/prepubertal exposure to BBP induced changes in the gene expression pattern of the mammary gland mainly by 21 days, just after the end of the treatment, and the number of modulated genes was low thereafter (Figure 2). In 21 days-old rats, BBP exposure resulted in the up-regulation of 515 genes, 141 of which were known, and included genes related to morphogenesis and cell differentiation, transcription factors, cell proliferation, response to stress, signal transduction, metabolism, transport and cell organization (Table 2). There was only one down-regulated gene (gad1). At 35 days of age the number of up-regulated genes decreased significantly to four unknown; two genes were down-regulated (one known -gad1-). By 50 days of age, 25 genes were up-regulated (10 known), and 14 genes were down-regulated (eight known). In the 100 day-old group, two known and one unknown genes were up-modulated (Table 2).

thumbnailFigure 2. Genomic changes induced by BBP. Number of differentially expressed genes in the mammary glands of rats exposed neonatally/prepubertally to BBP when compared to controls at different ages. Genes were obtained by confident analysis of microarrays data at p < 0.05.

Table 2. Microarrays results

Functional analysis of the genes found up-regulated in 21 days-old rats showed some categories that were significantly (p < 0.05) over-represented: cell proliferation and differentiation, signal transduction (TGF-beta signaling pathway, metabotropic glutamate receptor group I pathway, PI3 kinase pathway, endothelin signaling pathway and interleukin signaling pathway) and cell communication, (neuronal activities, anterior/posterior patterning). We also found over-representation of genes coding for cytoskeletal proteins among the down-modulated by 50 days.

Validation of microarrays by real time RT-PCR

To validate the results from microarrays we have chosen several genes that were differentially expressed in the mammary glands from BBP treated rats: transcription factors related to proliferation and differentiation (Foxg1, Hfh1, Nfyc, Ahr, Wt1, Tsn), other genes related to differentiation (Fabp3), adhesion (Cdh8) and communication (Gad1). We analyzed the gene expression of those genes by real time RT-PCR. Results of the validation are shown in Figure 3. In the BBP exposed animals Foxg1, Hfh1, Nfyc and Ahr were between 1.6- and 2-fold up-regulated in comparison to their gene expression in untreated animals. These differences, although discrete, were statistically significant for Hfh1 and Ahr, and close to significance for the other genes. We also obtained statistically significant increases in the expression of Fabp3, Wt1 and Cadh8. By 35 and 50 days, BBP exposure induced a significant decrease in the expression of Gad1. The decrease of this gene found by microarrays at 21 days as well as the increase in Tsn observed by 100 days of age was not confirmed by real time PCR.

thumbnailFigure 3. Comparative expression analyses of selected genes. Gene expression analyses by real time RT-PCR in comparison with results from microarrays. Grey bars represent results obtained in the microarrays, black bars represent results obtained in the real time RT-PCR. Results show the mean fold expression of the treated group related to control group. The functional categories significantly over-represented are indicated (21 days). *: significant differences compared to control group (p < 0.05).

Discussion

The objective of the present study was to investigate the effects of neonatal/prepubertal exposure to BBP in the rat mammary gland at different times of development. We have characterized the genomic signature of the mammary gland, observing modifications in the gene expression profile mainly at the end of the exposure and fewer changes thereafter. We also studied maturation parameters, finding also alterations compatible with a transitory effect of this environmental compound.

Although it is difficult to determine doses for experimental studies, the dose chosen for this study was based on the work of Singletary et al. [17], where they reported that treatment of female rats with BBP for 7 days prior to dosing with DMBA decreased mammary tumor incidence by 37%. The number of palpable tumors per rat at week 15 was significantly inhibited by 58% and 71% for animals administered BBP at 250 and 500 mg/kg, respectively, compared to controls. Our goal was to determine if neonatal and prepubertal only exposure to BBP would have a direct effect at day 21 and if this treatment would leave permanent alterations on mammary gland maturation, cell proliferation and gene expression at later time points. The offspring were exposed to the BBP via the mother's milk. Considering that the lactating dams were treated with 500 mg BBP/kg BW and this is dissipated in the mother before it makes its way via the milk to the 10 offspring, hence it is estimated that they receive 1/100 to 1/1000 of the mother's dose, i.e. 0,5–5 mg BBP/kg BW/day. This approximation places the rat BBP exposure near the EPA safe dose for humans of 0.2 mg/kg/day [18] and similar to a 2 μg/kg BW/day for adults, and lower than estimated exposure to infants and children (6 μg/kg BW/day) [1].

Neonatal/prepubertal exposure to BBP increased the relative uterine weight at 21 days but not thereafter, and decreased the body weight at day of vaginal opening (approximately day 30). In the literature there are conflicting results regarding the effect of BBP because it showed estrogenic effects in vitro but did not modify the estrogen-modulated gene calbindin-D9k in vivo [19]. However, BBP has also shown to increase the expression of progesterone receptor mRNA, a gene regulated by estrogens, in the preoptic area of adult ovariectomized rats [20]. Exposure to BBP in other phases of development induced alterations in the onset of puberty. Ashby et al. reported advanced vaginal opening by 1.1 days without modifications on uterine weight in rats by effect of prenatal plus prepubertal exposure to BBP, although the authors considered that this result could be related to the variability of the observation or, more interestingly, to an increase in the weight of the BBP treated animals [21]. In two generational studies, Aso et al. reported changes in anogenital distance of female offspring of rats treated with more than 100 mg/kg [22].

Little is known on the effects of BBP on mammary gland, another target of sexual maturity, and the present work indicates that BBP did not induce significant changes in the morphology of that gland, but changed the proliferative index of TEB by 35 days and in Lob1 at 35, 50 and 100 days. Although these modifications are subtle, we cannot rule out that they may have an effect if the mammary gland susceptibility to carcinogenesis, since TEB are the most susceptible epithelial structures to malignant transformation [23] and prevention of chemically-induced mammary carcinomas is accompanied by inhibition of cell proliferation [13,14]. Changes in the proliferative index of Lob1 in BBP exposed animals in comparison to control group can be associated with subtle alterations in mammary development. The proliferation of Lob1 is related to the differentiation of these structures to lobules type 2 [23]. Control group had the highest proliferative index in Lob1 by 50 days, concomitantly with the highest number of Lob1, at the age of profound changes that the mammary gland faces as a result of an intense hormonal activity [13,14]. BBP exposed animals had the maximum mitotic index in Lob1 by 35 days, decreasing by 50 days. In addition, the reduction of proliferation in Lob1 of BBP group by 100 days was not as important as in control group. The significant changes observed are compatible with slight changes in the timing of the breast tissue development by effect of BBP exposure.

We are the first to report that neonatal/prepubertal exposure to BBP induced modifications in the gene expression of the mammary tissue. Whereas the full relevance of these findings requires tumorigenesis studies, the characterization of the genomic signature of the mammary gland could be used as a predictor of the susceptibility to carcinogenesis. This statement is based on previous works in which changes in susceptibility to DMBA-induced carcinogenesis are accompanied by the induction of a specific genomic profile [13,14,24]. Neonatal/prepubertal exposure to BBP induced changes in the genomic profile of the mammary gland, mainly at the end of treatment (21 days) with a focus on genes involved in cell proliferation/differentiation, communication and signal transduction. Some of those genes related to differentiation were validated by real time RT-PCR, as Fabp3 and Wt1. Fabp3 belongs to the family of the fatty acid binding proteins, which have a role in fatty acid transport or metabolism. This member, however, is considered an identical homologue of the mammary-derived growth inhibitor (MDGI), with a role in the differentiation of the mammary gland [25]. MDGI is a suppressor of tumor growth and an inhibitor of tumor cell proliferation [26]. We have previously observed that this gene is up-regulated in the most differentiated structures of the mammary gland (lobules type 4) and silenced in breast cancer progression [27]. Another gene found significantly up-regulated by 21 days is Wt1, which has also been related to development of the mammary gland although its specific function is not well known [28,29]. This gene is essential for the development of the urogenital system but its role in the breast tissue is controversial. Wt1 presents a developmental pattern of expression in normal mammary gland, and changes observed in cell lines and tumors may be indicative of breast-tumor-related perturbations of Wt1 expression [29]. Cadherin 8 (cdh8), also found up-regulated in the exposed animals, is an adhesion molecule with an important role in brain and kidney development and differentiation [30,31]. We have also observed transitional changes in the expression of these three genes by effect of prepubertal exposure to other endocrine disruptor, bisphenol A [32]. Moreover, we also observed changes in the expression of Fabp3 when this endocrine disruptor was administered in utero. Thus, the exposure during the prenatal period to bisphenol A induced up-regulation or down-regulation of Fabp3 depending on the age tested, suggesting a "shift" in the normal development of the mammary gland, while prepubertal exposure to this endocrine disruptor induced a transient up-regulation just at the end of the treatment [32]. Fabp3 mRNA has also been reported as sensitive to the action of other hormone interfering agents, since the antiestrogens tamoxifen and ICI up-regulate its expression by stabilization of the mRNA [33]. These evidences suggest that genes related to differentiation of the mammary gland can be very sensitive to the effects of several compounds that may interfere with hormonal action, suggesting that, even if the effects are not strong enough to modify mammary histoarchitecture, they may change the molecular milieu of the gland.

By 21 days of age we also found a moderate increase in the expression of the transcription factors Foxg1, Hfh1 and Ahr. Foxg1 and Hfh1 have a role in the development and differentiation of several tissues, but to date their function in the mammary gland is unknown [34,35]. The aryl hydrocarbon receptor (Ahr) belongs to the PAS (per-arnt-sim) family of transcription factors. PAS proteins control important physiological processes, as toxicant metabolism, circadian rhythms and hormone signaling [36]. We have previously reported a role of Ahr in mammary differentiation, as an increased expression in primary cultures of Lob3 from parous women compared to primary cultures of Lob1 from nulliparous women breast was observed [37]. Ahr is known to interact with signaling pathways that are mediated by estrogen receptor and other hormone receptors. The mechanisms of Ahr-estrogen receptor crosstalk seem to be complex, and recent studies suggest direct activation of estrogen-regulated target genes by Ahr [38]. Thus, our results suggest that the increase in Ahr could be one of the mechanisms by which BBP treatment exerts the observed transitory effects at 21 days.

We also observed and validated down-regulation of Gad1 by 35 and 50 days of age. This gene codes for the protein Gad67, which catalyzes the production of the neurotransmitter GABA. Gad1 can be regulated by estrogens in rat brain [39,40], and expression changes in discrete regions of the rostral preoptic area during estrous cycle and with age has been described [41]. In addition, Gad1 has also a role in the control of puberty in nonhuman primates [42]. Modifications in the expression of Gad1 has been investigated in relation to some brain disorders, since Gad1 mRNA decreased in prefrontal cortex and other brain areas of schizophrenia and bipolar disorder patients with psychosis [43]. Experiments in rats have found down-modulation of Gad1 accompanied by hypermethylation of Gad1 promoter very likely mediated by the overexpression of DNA methyltransferase 1 in cortical GABAergic interneurons [44]. Thus, we can consider the possibility that the observed changes in Gad1 are due to an effect of BBP on the methylation of Gad1 promoter. BBP has also shown an effect on the methylation status of other genes in vitro, as the reported demethylation of ERα promoter-associated CpG islands in MCF-7 cells [45]. The GABA-ergic system has been involved in hormonal regulation and pathogenesis of breast cancer, since human and mouse tumors have shown higher GABA level and GAD activity than normal tissue [46]. Hence, although prepubertal exposure to BBP seems to have mainly a transitory effect on the genomic profile of the mammary gland, we cannot discard that subtle changes in the gland may have an effect later in life. Finally, very few genes were regulated by BBP exposure in the oldest age tested (100 days).

Conclusion

In summary, neonatal/prepubertal exposure to BBP had a transitory effect by increasing the relative uterine weight only at the end of exposure. This compound also changed the gene expression profile of the mammary gland mainly at the same age, although long-term subtle alterations can not be ruled out.

Methods

Animals and experimental design

Animals were housed in 12 hour light/dark cycle and controlled temperature in the University of Alabama at Birmingham animal facility. Pregnant Sprague Dawley CD rats (Charles River, Raleigh NC) were bred and maintained on phytoestrogen-free AIN-93G diet (Harlan Teklad, Madison WI). After parturition, litters were adjusted to 10 offspring per mother, and the lactating dams (10 dams per group) were gavaged with 500 mg BBP/kg body weight or an equivalent volume of sesame oil on days 2–20 (Mondays-Fridays only). We used a stock solution of 200 mg BBP/ml sesame oil, and dams were gastric intubated with 2,5 μl BBP/g BW (0.5 ml for a dam of 200 g) each day, preferable each late afternoon. The offspring were weaned at day 21 and continued on AIN-93G diet until day 70 where they were switched to AIN-93M diet. Vaginal opening was assessed as an index of maturation and body weight was recorded from a minimum of 27 female offsprings/group.

Female offspring were processed at 21, 35 ± 1, 50 ± 1 and 100 ± 2 days. For the cycling animals (the latter 3 ages), all females were killed in the estrus phase of the estrous cycle. The day of euthanasia the animals were anesthetized with ketamine (10 mg/100 g)/xylazine (1.5 mg/100 g), the fourth mammary glands were rapidly dissected and frozen for microarrays analysis (1 female offspring from each litter at each time point). Transsection of the aorta was performed after the tissues were properly collected. From another set of females from each litter, mammary glands were dissected for whole mount preparation and gland differentiation analysis, and the contralateral mammary glands were paraffin-blocked for cell proliferations studies. Ten glands representing all litters per group at each time point were used for each of the molecular, morphological and proliferation analyses. Uterine weight was also recorded at day of sacrifice and uterine weight/body weight ratios were calculated for each animal. Differences among groups were analyzed using Tukey's Test (p < 0.05). All animal studies were conducted in accordance with the University of Alabama at Birmingham Guidelines for Animal Use and Care.

Mammary whole mounts analysis

The excised mammary glands were fixed in 10% neutral-buffered formalin followed by de-fatting in acetone. The tissues were stained in alum carmine overnight and dehydrated in a series of graded ethanol. The glands were cleared in xylene and coverslipped with mounting media. Analysis of mammary epithelial structures (terminal end buds -TEBs-, terminal ducts -TDs-, alveolar buds -ABs- and lobules type 1 -Lob1-) was made through visual evaluation and computer-assisted image analysis using MetaMorph Imaging System (Universal Imaging Corporation, Downingtown PA). The outer 5-mm margin of the mammary whole mount was examined by light microscopy using the criteria of Russo and Russo [23]. This area represents the location of most actively proliferating TEB structures of the mammary gland for a young virgin rat [47], The total number of each epithelial structure was determined under an Olympus microscope using a 40× magnification objective. The collected data were analyzed using two-tailed unpaired t-tests.

Proliferative index analysis

Two hours before euthanasia animals were weighed and injected i.p. with 100 mg bromodeoxyuridine (BrdU, Sigma, St. Louis, MO) per kg body weight. The fourth abdominal mammary glands were excised and fixed in 10% neutral-buffered formalin, embedded in paraffin and sectioned at 4 μm thickness. Tissue sections were mounted on positively charged slides and immunocytochemically reacted with anti-BrdU monoclonal antibody (BioGenex, San Ramon, CA) using an automatic slide stainer (BioGenex). The incorporation of BrdU in the mammary glands was visualized using the streptavidin-biotin labeling system with 3, 3'-diaminobenzidine (DAB) as a color reaction substrate (BioGenex). The proliferative index was determined quantitatively using an Olympus BX40 microscope (60× magnification objective) as the percentage of DAB positive cells within each specific epithelial structure, i.e., TEBs, TD and ducts, ABs and Lob1. All cells in each structure were counted blinded to the groups. A minimum of 1,000 cells per tissue sample was examined. Data from different groups were analyzed by unpaired t-tests.

Characterization of the genomic profile

Gene expression profile was determined from total RNA using Agilent 60-mer oligo-microarrays slides containing 22,000 sequences. RNAs from the mammary glands were extracted with the RNeasy Mini Kit (Qiagen Inc., Valencia, CA). The quantity and purity of each RNA sample was determined from the 260- and 280-nm absorbance, and its integrity was determined (2100 BioAnalyzer; Agilent Technologies Inc., Palo Alto, CA). Only high quality RNA, with absorbance 260/280 greater than 1.8 and RNA integrity number (RIN) greater than 8 were considered for further analysis. 250 ng of pooled RNA from two-three different mammary glands were fluorescent labeled with Agilent Low Input RNA Fluorescent Linear Amplification Kit (Agilent Technologies) in the presence of Cy3-dCTP or Cy5-dCTP (Perkin Elmer, Wellesley, MA), and purified using RNeasy Mini Kit (Qiagen Inc., Valencia, CA).

Four separate experiments, each including two to three mammary gland cRNA pool were performed at each time point studied. 750 ng of Cy3 labeled control cRNA and 750 ng of Cy5 labeled BBP group cRNA were used to hybridize each slide at 60°C for 18 hours. Slides were washed and fluorescence images were captured with a laser scanner (Agilent Technologies). The obtained images were analyzed with Feature extraction software (Agilent Technologies) to verify the quality of the hybridization by the report of the outlier data.

Feature extraction software automatically flagged non-uniform features and population outliers, considering acceptable outlier numbers lower than 5% of total feature number on the array. Intensity of each fluorescent signal was quantified by ImaGene 5.6 software (BioDiscovery Inc., El Segundo, CA) and statistically analyzed using GeneSight 4.1.6 software (BioDiscovery) with Lowess normalization considering the arrays per group as replicates. We determined the genes with ≥ 1.4-fold differences by confidence analyses at p < 0.05 following manufacturer's recommendations. We performed functional analysis of the genes obtained as modulated by BBP exposure using PANTHER (Protein ANalysis THrough Evolutionary Relationships). Genes were annotated according to known biological process, molecular function and biological pathways. The number of genes in each category was statistically compared to NCBI: R. norvegicus genes database to look for under- and over-represented functional categories using the binomial test [48] and the Bonferroni correction for multiple testing. Genes with known function were also classified using Gene Ontology database.

Gene expression analysis by real-time RT-PCR

Several genes found differentially expressed in the mammary glands of rats exposed to BBP were chosen for confirmation by real-time reverse transcription PCR (RT-PCR). Beta-actin was used as a control. All RT-PCR reactions were performed on the ABI Prism 7700 Sequence Detection System (Applied Biosystems, Foster City, CA) using the fluorescent TaqMan methodology (TaqMan One Step RT-PCR Master Mix Reagents, Applied Biosytems) following manufacturer's protocol. 100 ng of total RNA in 50 μl of final volume was reverse-transcribed (30 min at 48°C) and amplified following 1 cycle of 10 min at 95°C, and 40 cycles of 15 s denaturation at 95°C and 60 s annealing at 60°C. Each gene was normalized to beta-actin and data was analyzed using two-tailed unpaired t-test.

Abbreviations

AB: alveolar bud;

BBP: butyl benzyl phthalate;

Lob1: lobule type 1;

PCR: polymerase chain reaction;

RT: reverse transcription;

TD: terminal duct;

TEB: terminal end bud.

Authors' contributions

RM performed the molecular experiments and contributed in the writing of the manuscript. RW carried out the morphological and proliferative index analyses. IHR participated in the design and coordination of the study, and helped in the writing of the manuscript. DAM and GAB participated in the analysis of the microarrays. CAL participated in the coordination of the study and performed the animal assays. JR conceived of the study, coordinated the work and helped in the writing of the manuscript. All authors have read and approved the final manuscript.

Acknowledgements

This research work was supported by the Breast Cancer and the Environment Research Centers grant number U01 ES/CA 12771 from the National Institute of Environmental Health Sciences (NIEHS), and the National Cancer Institute (NCI), NIH, DHHS. Its contents are solely the responsibility of the authors and do not necessarily represent the official views of the NIEHS or NCI, NIH.

References

  1. IPCS: Concise international chemical assessment document 17 - Butyl benzyl phthalate. Edited by WHO . Geneva, Switzerland ; 1999.

  2. Okubo T, Suzuki T, Yokoyama Y, Kano K, Kano I: Estimation of estrogenic and anti-estrogenic activities of some phthalate diesters and monoesters by MCF-7 cell proliferation assay in vitro.

    Biol Pharm Bull 2003, 26(8):1219-1224. PubMed Abstract | Publisher Full Text OpenURL

  3. Ema M, Miyawaki E: Effects on development of the reproductive system in male offspring of rats given butyl benzyl phthalate during late pregnancy.

    Reprod Toxicol 2002, 16(1):71-76. PubMed Abstract | Publisher Full Text OpenURL

  4. Ema M, Miyawaki E, Hirose A, Kamata E: Decreased anogenital distance and increased incidence of undescended testes in fetuses of rats given monobenzyl phthalate, a major metabolite of butyl benzyl phthalate.

    Reprod Toxicol 2003, 17(4):407-412. PubMed Abstract | Publisher Full Text OpenURL

  5. Gray LE Jr., Ostby J, Furr J, Price M, Veeramachaneni DN, Parks L: Perinatal exposure to the phthalates DEHP, BBP, and DINP, but not DEP, DMP, or DOTP, alters sexual differentiation of the male rat.

    Toxicol Sci 2000, 58(2):350-365. PubMed Abstract | Publisher Full Text OpenURL

  6. NTP: NTP-CERHR Monograph on the Potential Human Reproductive and Developmental Effects of Butyl Benzyl Phthalate (BBP).

    NTP CERHR MON 2003, (5):i-III90. OpenURL

  7. European-Commission: Directive 2005/84/EC of the European Parliament and of the Council of 14 December 2005.

    Official Journal of the European Union 2005. OpenURL

  8. Swan SH, Main KM, Liu F, Stewart SL, Kruse RL, Calafat AM, Mao CS, Redmon JB, Ternand CL, Sullivan S, Teague JL: Decrease in anogenital distance among male infants with prenatal phthalate exposure.

    Environ Health Perspect 2005, 113(8):1056-1061. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  9. Reddy BS, Rozati R, Reddy BV, Raman NV: Association of phthalate esters with endometriosis in Indian women.

    Bjog 2006, 113(5):515-520. PubMed Abstract | Publisher Full Text OpenURL

  10. Nagao T, Ohta R, Marumo H, Shindo T, Yoshimura S, Ono H: Effect of butyl benzyl phthalate in Sprague-Dawley rats after gavage administration: a two-generation reproductive study.

    Reprod Toxicol 2000, 14(6):513-532. PubMed Abstract | Publisher Full Text OpenURL

  11. Birnbaum LS, Fenton SE: Cancer and developmental exposure to endocrine disruptors.

    Environ Health Perspect 2003, 111(4):389-394. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  12. Althuis MD, Dozier JM, Anderson WF, Devesa SS, Brinton LA: Global trends in breast cancer incidence and mortality 1973-1997.

    Int J Epidemiol 2005, 34(2):405-412. PubMed Abstract | Publisher Full Text OpenURL

  13. Russo IH, Russo J: Role of hormones in cancer initiation and progression.

    J Mam Gland Biol Neoplasia 1998, 3:49-61. Publisher Full Text OpenURL

  14. Russo J, Russo IH: Role of differentiation in pathogenesis and prevention of breast cancer.

    Endrocr Related Cancer 1997, 4:7-21. OpenURL

  15. Boylan ES, Calhoon RE: Mammary tumorigenesis in the rat following prenatal exposure to diethylstilbestrol and postnatal treatment with 7,12-dimethylbenz[a]anthracene.

    J Toxicol Environ Health 1979, 5(6):1059-1071. PubMed Abstract OpenURL

  16. Lamartiniere CA, Holland MB: Neonatal diethylstilbestrol prevents spontaneously developing mammary tumors. In Proceed First Intern Symp Hormonal Carcinogenesis. Edited by Li JJ, Nandi S, Li SA. New York , Springer-Verlag; 1992 :305-308. [] OpenURL

  17. Singletary K, MacDonald C, Wallig M: The plasticizer benzyl butyl phthalate (BBP) inhibits 7,12-dimethylbenz[a]anthracene (DMBA)-induced rat mammary DNA adduct formation and tumorigenesis.

    Carcinogenesis 1997, 18(8):1669-1673. PubMed Abstract | Publisher Full Text OpenURL

  18. US-EPA: Butyl benzyl phthalate (CASRN 85-68-7). Washington, DC , U.S. Environmental Protection Agency; 2003:http://www.epa.gov/iris/subst/0293.htm.

  19. Hong EJ, Ji YK, Choi KC, Manabe N, Jeung EB: Conflict of estrogenic activity by various phthalates between in vitro and in vivo models related to the expression of Calbindin-D9k.

    J Reprod Dev 2005, 51(2):253-263. PubMed Abstract | Publisher Full Text OpenURL

  20. Funabashi T, Kawaguchi M, Kimura F: The endocrine disrupters butyl benzyl phthalate and bisphenol A increase the expression of progesterone receptor messenger ribonucleic acid in the preoptic area of adult ovariectomized rats.

    Neuroendocrinology 2001, 74(2):77-81. PubMed Abstract | Publisher Full Text OpenURL

  21. Ashby J, Tinwell H, Lefevre PA, Odum J, Paton D, Millward SW, Tittensor S, Brooks AN: Normal sexual development of rats exposed to butyl benzyl phthalate from conception to weaning.

    Regul Toxicol Pharmacol 1997, 26(1 Pt 1):102-118. Publisher Full Text OpenURL

  22. Aso S, Ehara H, Miyata K, Hosyuyama S, Shiraishi K, Umano T, Minobe Y: A two-generation reproductive toxicity study of butyl benzyl phthalate in rats.

    J Toxicol Sci 2005, 30 Spec No.:39-58. PubMed Abstract | Publisher Full Text OpenURL

  23. Russo J, Russo IH: DNA labeling index and structure of the rat mammary gland as determinants of its susceptibility to carcinogenesis.

    J Natl Cancer Inst 1978, 61(6):1451-1459. PubMed Abstract OpenURL

  24. Balogh GA, Heulings R, Mailo DA, Russo PA, Sheriff F, Russo IH, Moral R, Russo J: Genomic signature induced by pregnancy in the human breast.

    Int J Oncol 2006, 28(2):399-410. PubMed Abstract | Publisher Full Text OpenURL

  25. Yang Y, Spitzer E, Kenney N, Zschiesche W, Li M, Kromminga A, Muller T, Spener F, Lezius A, Veerkamp JH, et al.: Members of the fatty acid binding protein family are differentiation factors for the mammary gland.

    J Cell Biol 1994, 127(4):1097-1109. PubMed Abstract | Publisher Full Text OpenURL

  26. Huynh HT, Larsson C, Narod S, Pollak M: Tumor suppressor activity of the gene encoding mammary-derived growth inhibitor.

    Cancer Res 1995, 55(11):2225-2231. PubMed Abstract | Publisher Full Text OpenURL

  27. Hu YF, Russo IH, Ao X, Russo J: Mammary derived growth inhibitor (MDGI) cloned from human breast epithelial cells is expressed in fully differentiated lobular structures.

    Int J Oncol 1997, 11:5-11. OpenURL

  28. Loeb DM, Evron E, Patel CB, Sharma PM, Niranjan B, Buluwela L, Weitzman SA, Korz D, Sukumar S: Wilms' tumor suppressor gene (WT1) is expressed in primary breast tumors despite tumor-specific promoter methylation.

    Cancer Res 2001, 61(3):921-925. PubMed Abstract | Publisher Full Text OpenURL

  29. Silberstein GB, Van Horn K, Strickland P, Roberts CT Jr., Daniel CW: Altered expression of the WT1 wilms tumor suppressor gene in human breast cancer.

    Proc Natl Acad Sci U S A 1997, 94(15):8132-8137. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  30. Blaschke S, Mueller CA, Markovic-Lipkovski J, Puch S, Miosge N, Becker V, Mueller GA, Klein G: Expression of cadherin-8 in renal cell carcinoma and fetal kidney.

    Int J Cancer 2002, 101(4):327-334. PubMed Abstract | Publisher Full Text OpenURL

  31. Korematsu K, Nishi T, Okamura A, Goto S, Morioka M, Hamada J, Ushio Y: Cadherin-8 protein expression in gray matter structures and nerve fibers of the neonatal and adult mouse brain.

    Neuroscience 1998, 87(1):303-315. PubMed Abstract | Publisher Full Text OpenURL

  32. Moral R, Balogh GA, Mailo D, Russo IH, Lamartiniere C, Russo J: Effect of the contaminant bisphenol A (BPA) on the genomic signature of the rat mammary gland.

    Proc Amer Assoc Cancer Res 2005, 46:491-492. OpenURL

  33. Huynh H, Pollak M: Stabilization of mammary-derived growth inhibitor messenger RNA by antiestrogens.

    Clin Cancer Res 1997, 3(11):2151-2156. PubMed Abstract | Publisher Full Text OpenURL

  34. Hong HK, Noveroske JK, Headon DJ, Liu T, Sy MS, Justice MJ, Chakravarti A: The winged helix/forkhead transcription factor Foxq1 regulates differentiation of hair in satin mice.

    Genesis 2001, 29(4):163-171. PubMed Abstract | Publisher Full Text OpenURL

  35. Martynoga B, Morrison H, Price DJ, Mason JO: Foxg1 is required for specification of ventral telencephalon and region-specific regulation of dorsal telencephalic precursor proliferation and apoptosis.

    Dev Biol 2005, 283(1):113-127. PubMed Abstract | Publisher Full Text OpenURL

  36. Gu YZ, Hogenesch JB, Bradfield CA: The PAS superfamily: sensors of environmental and developmental signals.

    Annu Rev Pharmacol Toxicol 2000, 40:519-561. PubMed Abstract | Publisher Full Text OpenURL

  37. Guo S, Russo IH, Russo J: Difference in gene expression profile in breast epithelial cells from women with different reproductive history.

    Int J Oncol 2003, 23(4):933-941. PubMed Abstract OpenURL

  38. Wiseman A: Oestrogen-receptors (ER) are likely to be promiscuous: wider role for oestrogens and mimics.

    Med Hypotheses 2005, 65(4):760-765. PubMed Abstract | Publisher Full Text OpenURL

  39. Nakamura NH, Rosell DR, Akama KT, McEwen BS: Estrogen and ovariectomy regulate mRNA and protein of glutamic acid decarboxylases and cation-chloride cotransporters in the adult rat hippocampus.

    Neuroendocrinology 2004, 80(5):308-323. PubMed Abstract | Publisher Full Text OpenURL

  40. Wagner EJ, Ronnekleiv OK, Bosch MA, Kelly MJ: Estrogen biphasically modifies hypothalamic GABAergic function concomitantly with negative and positive control of luteinizing hormone release.

    J Neurosci 2001, 21(6):2085-2093. PubMed Abstract | Publisher Full Text OpenURL

  41. Cashion AB, Smith MJ, Wise PM: Glutamic acid decarboxylase 67 (GAD67) gene expression in discrete regions of the rostral preoptic area change during the oestrous cycle and with age.

    J Neuroendocrinol 2004, 16(8):711-716. PubMed Abstract | Publisher Full Text OpenURL

  42. Kasuya E, Nyberg CL, Mogi K, Terasawa E: A role of gamma-amino butyric acid (GABA) and glutamate in control of puberty in female rhesus monkeys: effect of an antisense oligodeoxynucleotide for GAD67 messenger ribonucleic acid and MK801 on luteinizing hormone-releasing hormone release.

    Endocrinology 1999, 140(2):705-712. PubMed Abstract | Publisher Full Text OpenURL

  43. Guidotti A, Auta J, Davis JM, Di-Giorgi-Gerevini V, Dwivedi Y, Grayson DR, Impagnatiello F, Pandey G, Pesold C, Sharma R, Uzunov D, Costa E: Decrease in reelin and glutamic acid decarboxylase67 (GAD67) expression in schizophrenia and bipolar disorder: a postmortem brain study.

    Arch Gen Psychiatry 2000, 57(11):1061-1069. PubMed Abstract | Publisher Full Text OpenURL

  44. Dong E, Agis-Balboa RC, Simonini MV, Grayson DR, Costa E, Guidotti A: Reelin and glutamic acid decarboxylase67 promoter remodeling in an epigenetic methionine-induced mouse model of schizophrenia.

    Proc Natl Acad Sci U S A 2005, 102(35):12578-12583. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  45. Kang SC, Lee BM: DNA methylation of estrogen receptor alpha gene by phthalates.

    J Toxicol Environ Health A 2005, 68(23-24):1995-2003. PubMed Abstract | Publisher Full Text OpenURL

  46. Opolski A, Mazurkiewicz M, Wietrzyk J, Kleinrok Z, Radzikowski C: The role of GABA-ergic system in human mammary gland pathology and in growth of transplantable murine mammary cancer.

    J Exp Clin Cancer Res 2000, 19(3):383-390. PubMed Abstract OpenURL

  47. Russo IH, Russo J: Mammary gland neoplasia in long-term rodent studies.

    Environ Health Perspect 1996, 104(9):938-967. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  48. Cho RJ, Campbell MJ: Transcription, genomes, function.

    Trends Genet 2000, 16(9):409-415. PubMed Abstract | Publisher Full Text OpenURL