TY - JOUR AU - Vandenberg, Laura N AB - Abstract Unconventional oil and gas (UOG) operations, which combine hydraulic fracturing (fracking) and directional drilling, involve the use of hundreds of chemicals, including many with endocrine-disrupting properties. Two previous studies examined mice exposed during early development to a 23-chemical mixture of UOG compounds (UOG-MIX) commonly used or produced in the process. Both male and female offspring exposed prenatally to one or more doses of UOG-MIX displayed alterations to endocrine organ function and serum hormone concentrations. We hypothesized that prenatal UOG-MIX exposure would similarly disrupt development of the mouse mammary gland. Female C57Bl/6 mice were exposed to ~3, ~30, ~ 300, or ~3000 μg/kg/d UOG-MIX from gestational day 11 to birth. Although no effects were observed on the mammary glands of these females before puberty, in early adulthood, females exposed to 300 or 3000 μg/kg/d UOG-MIX developed more dense mammary epithelial ducts; females exposed to 3 μg/kg/d UOG-MIX had an altered ratio of apoptosis to proliferation in the mammary epithelium. Furthermore, adult females from all UOG-MIX–treated groups developed intraductal hyperplasia that resembled terminal end buds (i.e., highly proliferative structures typically seen at puberty). These results suggest that the mammary gland is sensitive to mixtures of chemicals used in UOG production at exposure levels that are environmentally relevant. The effect of these findings on the long-term health of the mammary gland, including its lactational capacity and its risk of cancer, should be evaluated in future studies. Unconventional oil and gas (UOG) operations combine hydraulic fracturing (fracking) and directional drilling. These techniques were developed to collect deposits of oil and natural gas found in deep underground shale beds in low-permeability geologic formations (1). During the fracking process, a mixture of water and chemicals is pumped deep into the shale bed under high pressure, fracturing the reservoir rock, and releasing deposits of gas and/or oil, which can then be recovered at the surface. More than 1000 different chemicals are reportedly used during UOG operations for a range of purposes, including compounds that act as bactericides, stabilizers for the clay in the ground, chemicals that alter friction and fluid viscosity, and others, although each individual site typically uses only 12 to 24 of these compounds. During UOG operations, up to several million gallons of water are injected per well, and a mixture of injected fluids and target formation water are collected throughout the life of the producing well. With >17 million Americans living within 1 mile of an oil and gas well (2), concerns have been raised about the possibility of contamination of surface and groundwater by the released oil and gas, the numerous inorganic compounds that are liberated from target geologic layers (e.g., trace metals, radioactive isotopes, minerals), and the chemicals used in well injection (3, 4). More than 1000 chemicals have been identified in hydraulic fracturing fluids and waste water and/or have been reported to be used by the industry (5, 6). Many of these chemicals are known developmental and reproductive toxicants (7). Furthermore, recent evaluations found that >100 of these chemicals are known or suspected endocrine-disrupting chemicals (EDCs) (4, 8–10) (i.e., compounds that interfere with hormone action) (11). Water samples collected in drilling-dense or UOG wastewater-affected areas of the United States have exhibited disruption of the estrogen, androgen, progesterone, glucocorticoid, and thyroid receptors (4, 8, 12). In 2015, Kassotis et al. (13) evaluated 24 chemicals that were reported by industry as commonly used and/or produced by UOG operations to determine whether they displayed endocrine-disrupting properties. With cell-based reporter gene assays, their study revealed antiandrogenic, antiestrogenic, antiprogestogenic, antithyroidogenic, and antiglucocorticogenic activities for many of these compounds. When evaluated as mixtures, additive and, in some cases, synergistic antagonism of these receptors was also observed. Kassotis et al. (13) then evaluated the effects of a 23-chemical mixture of UOG chemicals (UOG-MIX) on male mice. This mixture included the 24 chemicals originally assessed, absent bisphenol A, a well-characterized EDC that has been evaluated at length previously (14, 15) but that is not directly used in UOG extraction, as reported by the industry (16). Only one of the chemicals included in this list (Table 1), benzene, was evaluated in a review of 216 chemicals for carcinogenic effects in the mammary gland, highlighting the general lack of knowledge on these chemicals (17). Male mice exposed to environmentally relevant doses of UOG-MIX during prenatal development displayed an increased testicular weight before puberty and, in adulthood, decreased sperm counts, increased serum testosterone concentrations, and alterations to the weight of additional organs, including the heart and thymus (13). A second study revealed the effects of developmental exposures to UOG-MIX on the female siblings (18). Exposed females had alterations to the number and developmental stages of ovarian follicles measured before puberty and in adulthood. The weight of several organs, including the uterus, ovary, and heart, was also affected. Furthermore, the serum concentrations of several hormones were disrupted in these exposed females. These results suggest the possibility that other hormone-sensitive organs might be disrupted by exposure to UOG chemical mixtures. Table 1. Chemicals in UOG Mixture and Hormone Receptor Antagonist Activity, as Listed in (13) Chemical Name  Receptor Antagonist Activity  ER  AR  PR  TR  GR  1,2,4-Trimethylbenzene            2-(2-Methoxyethoxy) ethanol  X  X      X  2-Ethylhexanol  X  X  X    X  2-Methyl-4-isothiazolin-3-one  X  X      X  Acrylamide    X        Benzene  X  X        Bronopol  X  X  X    X  Cumene  X  X  X      Diethanolamine  X  X  X      Ethoxylated nonylphenol  X  X  X  X  X  Ethoxylated octylphenol  X  X  X  X  X  Ethylbenzene  X  X        Ethylene glycol  X  X  X  X  X  Ethylene glycol monobutyl ether  X  X    X    Naphthalene  X  X  X  X  X  N,N-dimethylformamide  X  X  X      Phenol  X  X  X      Propylene glycol  X          Sodium tetraborate decahydrate  X  X        Styrene  X      X  X  Toluene  X  X        Triethylene glycol  X  X        Chemical Name  Receptor Antagonist Activity  ER  AR  PR  TR  GR  1,2,4-Trimethylbenzene            2-(2-Methoxyethoxy) ethanol  X  X      X  2-Ethylhexanol  X  X  X    X  2-Methyl-4-isothiazolin-3-one  X  X      X  Acrylamide    X        Benzene  X  X        Bronopol  X  X  X    X  Cumene  X  X  X      Diethanolamine  X  X  X      Ethoxylated nonylphenol  X  X  X  X  X  Ethoxylated octylphenol  X  X  X  X  X  Ethylbenzene  X  X        Ethylene glycol  X  X  X  X  X  Ethylene glycol monobutyl ether  X  X    X    Naphthalene  X  X  X  X  X  N,N-dimethylformamide  X  X  X      Phenol  X  X  X      Propylene glycol  X          Sodium tetraborate decahydrate  X  X        Styrene  X      X  X  Toluene  X  X        Triethylene glycol  X  X        X indicates the presence of receptor antagonist activity as measured via transiently transfected reporter gene assays in human cells. Abbreviations: GR, growth receptor; PR, progesterone receptor; TR, thyroid receptor. View Large The mouse mammary gland has proved to be an excellent model to study the effects of EDCs (19–21). The development of the mammary gland is dependent on estrogen, progesterone, prolactin, testosterone, and growth hormone, making it an integrated biological endpoint that is sensitive to the agonists and antagonists of different hormone receptors (22). The estrogen receptor (ER) is expressed in the mesenchymal compartment of the fetal gland; however, its expression shifts to the epithelial compartment in the adult gland (23). Although ERα knockout mice have mammary glands that are indistinguishable from those of wild-type controls before puberty, they are visibly stunted compared with controls once puberty begins, highlighting the importance of estrogen for growth of the gland (24, 25). The androgen receptor (AR) is also expressed in the mesenchyme of the fetal mammary gland, and the testosterone produced by the testes in male fetuses causes the mesenchyme to condense around the epithelium, detaching the epithelium from the skin. Thus, the male mice will typically not have nipples (26, 27). Antiandrogenic compounds can lead to nipple retention in exposed males (28). Testosterone is likely to have a physiological role during postnatal development of the female gland because the mammary glands of female AR knockout mice show impaired ductal growth in postnatal life, indicating a role in ductal elongation (29). A wide range of EDCs, including several pharmaceutical compounds, naturally occurring EDCs, and industrial compounds with varied modes of action, have been demonstrated to alter the development of the mammary gland, with effects that typically manifest at puberty and in adulthood, when endogenous hormones induce its growth (20, 30). The characteristics that are commonly evaluated to determine the level of development and the disruption of development by EDCs include a number of branching points, extension of the ductal epithelium into the stroma, the number and size of terminal end buds (TEBs) (e.g., highly proliferative epithelial structures present during puberty), and the presence of alveolar buds and lobuloalveolar units (e.g., structures that will produce milk in lactating females), among others (31). Based on the hormone receptor antagonism of the 23-chemical mixture and the effects this mixture induced in the endpoints relevant to the hypothalamic-pituitary-gonadal axis in exposed mice, we hypothesized that prenatal exposure to the UOG-MIX would disrupt development of the mouse mammary gland. We evaluated the mammary glands in female mice after gestational exposure to one of four doses of the UOG-MIX. Consistent with our hypothesis, we characterized the substantial effects of this mixture on mammary gland morphology in adulthood and documented the presence of hyperplastic lesions. Materials and Methods Animal husbandry and chemical administration C57BL/6J mice were housed in sterile polysulfone cages under temperature- and light-controlled (12-hour light, 12-hour dark) conditions in a barrier animal facility. The mice were fed LabDiet 5053 and provided acidified water ad libitum from glass bottles. Ten-week-old mice were mated, and the day of vaginal plug formation was denoted as gestational day 0. On gestational day 11, the dams were randomized to the treatment groups and provided with experimental treatment in their drinking water. The test concentrations included a 0.2% ethanol vehicle; flutamide, a known antiandrogenic pharmaceutical agent (mechanistic antiandrogen control, 167 μg/mL); and four concentrations of the 23-chemical mixture (UOG-MIX), with each individual chemical present at 0.01, 0.10, 1.0, and 10 μg/mL. The composition of the chemical mixture and hormone receptor antagonist activity are listed in Table 1. Water intake was monitored by weighing the drinking water bottle every day of the experiment. Based on intake, the chemical exposure was estimated at 3, 30, 300, or 3000 μg/kg body weight/d for the mixtures. These treatment groups are referred to as MIX-3, MIX-30, MIX-300, and MIX-3000, respectively. Intake of the antiandrogenic control was estimated at 50 mg flutamide/kg body weight/d. Water intake did not differ in experimental groups relative to the vehicle control (data not shown). Experimental treatment was provided until birth; the dams were then reverted to standard acidified water when the pups were first detected. Litters with <2 males and females were removed from the analyses owing to concerns of gestational hormone exposure. Litters with >2 males and females were left unaltered, because culling within litters has been shown to alter the feeding, behavior, and physiology of the remaining pups (32). All experimental procedures were performed according to an approved University of Missouri Animal Care and Use Committee protocol and were in accordance with the National Research Council's Guide for the Care and Use of Laboratory Animals. Tissue collection All mice were euthanized by carbon dioxide asphyxiation and cardiac puncture, and the mammary tissues were excised. One randomly selected female pup from each litter underwent necropsy on either postnatal day (PND) 21 or PND85. At both ages, one fourth of the inguinal mammary gland was whole mounted on a glass slide, fixed in neutral-buffered formalin overnight, and processed for carmine staining using protocols described previously (33). Whole mounts were preserved in K-pax sealed bags with methyl salicylate. The contralateral fourth inguinal mammary gland was formalin-fixed for 24 hours, washed in phosphate-buffered saline, and stored in 70% ethanol for histological evaluation. Morphological analysis of whole mount mammary glands Whole mount specimens were imaged using a Zeiss dissecting microscope with AxioCam HRc digital camera. For analyses of the prepubertal mammary gland (PND21), images were captured of the full ductal tree and its position relative to the central lymph node. The following measurements were taken using methods developed previously (33): total ductal area, measured by quantifying the area subtended by the ducts; ductal extension, quantified as the furthest growth of the ductal tree measured from the center of the lymph node; and the number of branching points, counted throughout the entire gland. No mice had visible TEBs; thus, these structures were not quantified. For analyses of the adult gland, two photographs were taken, one of the entire mammary gland (at ×3 magnification) and one anterior to the central lymph node (at ×13.5 magnification); the former was used to evaluate the presence of TEB-like structures and the latter for unbiased stereological evaluations using ZEN imaging software (Zeiss). In brief, to quantify the volume fraction of epithelial structures, a 130-point grid was superimposed over each anterior photograph. The structure present on each crosshair was counted; the individual epithelial structures that were evaluated included ducts (when the crosshair hit the middle of a duct), terminal ducts/terminal ends (when the crosshair hit a blunt end of a duct), and alveolar buds (34). The volume fraction of each type of structure was calculated by counting the number of crosshairs that hit each structure divided by the total number of crosshairs hitting the mammary tissue (typically, 130). Volume fraction of all epithelium was calculated by summing all epithelial structures (i.e., ducts, terminal ends, and alveolar buds). Excision of mammary tissue from whole mounts For whole mounts that displayed unusual or TEB-like structures, small areas of the gland, including these abnormal structures, were excised using a scalpel and dissecting microscope. The remainder of the whole mount was then rebagged using methyl salicylate. The excised parts of the whole mounts were washed, processed through a series of alcohols, embedded in paraffin, and sectioned using the methods described in the following sections. Additional areas from the whole mount specimens with a normal appearance were also excised for use as control tissues. All excised tissues were coded such that additional analyses could be conducted by experimenters who were unaware of their origin (e.g., unusual vs normal). Histological evaluation One mammary gland from each sample was processed through a series of dehydrating alcohol washes and embedded in paraffin under vacuum; 5-µm sections were produced using a rotary microtome and placed on Superfrost positively charged slides (Fisher Scientific). This process was also used to evaluate excised tissue from whole mount mammary glands. For the histological evaluations, the sections were deparaffinized in xylene and rehydrated in a series of alcohol washes. They were then stained using Harris hematoxylin and eosin (Fisher Scientific), dehydrated through an alcohol series, washed in xylene, and mounted using a permanent mounting medium (Fisher Scientific). The slides were examined using a Zeiss Observer Z1 inverted light microscope at ×40 magnification. Images were captured using an AxioCam HRc digital camera and evaluated with ZEN imaging software (Zeiss) for the presence of hyperplastic ducts. Immunohistochemistry Immunohistochemical analysis for ERα and Ki67, a marker of proliferation, was performed as previously described (34). Primary antibodies were used at 1:1000 (ERα; catalog no. 06-935; Millipore; Ki67; catalog no. 9106-S1; Thermo Fisher Scientific; Table 2). For feasibility, the samples were examined in the control, MIX-3, and MIX-3000 groups only; these groups were selected to evaluate a wide range of doses. The samples were visualized using a Zeiss Observer Z1 inverted light microscope at ×40 magnification, and the images were captured using an AxioCam HRc digital camera. The images were analyzed using ZEN imaging software (Zeiss). For each sample, two or three fields of view were selected arbitrarily (depending on the size of the epithelial ducts) and imaged at ×40; the expression of each marker was quantified in all ducts within these images. At least 200 epithelial cells were assessed for each antigen; each cell was counted as either positively expressing the marker of interest (brown owing to diaminobenzidine, the colorimetric reaction used to visualize immunohistochemical reactions), or no expression (blue; hematoxylin counterstain). Table 2. Antibody Table Peptide/ Protein Target  Antigen Sequence (if Known)  Name of Antibody  Manufacturer; Catalog No.  Species Raised in; Monoclonal or Polyclonal  RRID  Dilution Used  ERα    Anti-ERα (C1355)  Millipore; 06-935  Rabbit; polyclonal  AB_310305  1:1000  Ki67    Ki67  Thermo Fisher Scientific; RM-9106-S1  Rabbit; monoclonal  AB_149792  1:1000  Secondary    Biotinylated goat anti-rabbit IgG  Abcam; ab64256  Goat; polyclonal  AB_2661852  Ready to use (5 μg/mL)  Peptide/ Protein Target  Antigen Sequence (if Known)  Name of Antibody  Manufacturer; Catalog No.  Species Raised in; Monoclonal or Polyclonal  RRID  Dilution Used  ERα    Anti-ERα (C1355)  Millipore; 06-935  Rabbit; polyclonal  AB_310305  1:1000  Ki67    Ki67  Thermo Fisher Scientific; RM-9106-S1  Rabbit; monoclonal  AB_149792  1:1000  Secondary    Biotinylated goat anti-rabbit IgG  Abcam; ab64256  Goat; polyclonal  AB_2661852  Ready to use (5 μg/mL)  Abbreviations: IgG, immunoglobulin G; RRID, Research Resource Identifier. View Large Terminal deoxynucleotidyltransferase-mediated dUTP nick end labeling assay The Trevigen TACS 2 TdT-DAB in situ apoptosis detection kit was used for detection of apoptotic cells in the mammary tissue sections. All samples were counterstained with Harris hematoxylin, dehydrated, mounted with a permanent mounting medium, and imaged with a Zeiss Observer Z1 inverted light microscope at ×40 magnification. Quantification of terminal deoxynucleotidyltransferase-mediated dUTP nick end labeling (TUNEL) was completed in two or three fields of view selected arbitrarily; all ducts within these images were evaluated. At least 200 epithelial cells were assessed; each cell was counted as either TUNEL-positive (brown owing to diaminobenzidine, the colorimetric reaction used to visualize the TUNEL reactions) or no expression (blue; hematoxylin counterstain). Statistical analysis The SPSS statistical software package, version 22 (IBM Corp.), was used for all statistical analyses. Analysis of variance followed by Fisher post hoc tests, was used to assess for differences between the control and UOG-MIX treatment groups for each age (PND21 and PND85). Independent samples t tests were used to compare the control and flutamide-treated groups (35). A χ2 test was performed to compare the incidence of hyperplasia in the mammary gland epithelium of control and UOG-MIX–exposed mice. To account for litter effects, only 1 mouse was selected from each litter for each age evaluated. For all statistical tests, the results were considered statistically significant at P < 0.05. All results are presented as the mean ± standard error of the mean and were collected and analyzed by experimenters who were unaware of the treatment. Results Prepubertal mammary gland morphology was not altered by developmental exposure to UOG-MIX The mice were exposed to vehicle or one of four doses of a 23-chemical mixture (UOG-MIX) during gestation. To determine the effects of the UOG chemicals on the female prepubertal mammary gland, morphological evaluations were first conducted at PND21. TEBs were not present in any glands from any treatment group, consistent with the prepubertal stage of mammary gland development (Fig. 1A; data not shown). Morphometric evaluations revealed an inverse association between treatment and total ductal area (e.g., smaller epithelial trees in UOG-MIX–treated females), although no statistically significant differences were found in this growth parameter between controls and the UOG-MIX treatments (Fig. 1B). Flutamide-treated females had significantly smaller ductal trees compared with the controls (P < 0.05; independent samples t test; Fig. 1B). Ductal extension and the total number of branching points were not different in any treatment group, including the flutamide-treated females (Table 3). Collectively, these data suggest that prenatal exposure to the 23-chemical mixture did not alter the morphological features of the female mammary glands before puberty. Figure 1. View largeDownload slide No substantial effects from the developmental exposure to fracking mixtures were observed on the morphology of the prepubertal mammary gland. (A) Examples of whole mount mammary glands from vehicle, MIX-3000, and flutamide (FLUT)-treated female mice. Scale bar = 0.5 mm. (B) Quantification of ductal area suggested an inverse relationship between UOG-MIX treatment and ductal area, although this difference was not statistically significant. The flutamide-treated female mice had significantly smaller ductal trees. ^P < 0.05, independent samples t test, comparing control and flutamide-treated female mice. LN, lymph node. Figure 1. View largeDownload slide No substantial effects from the developmental exposure to fracking mixtures were observed on the morphology of the prepubertal mammary gland. (A) Examples of whole mount mammary glands from vehicle, MIX-3000, and flutamide (FLUT)-treated female mice. Scale bar = 0.5 mm. (B) Quantification of ductal area suggested an inverse relationship between UOG-MIX treatment and ductal area, although this difference was not statistically significant. The flutamide-treated female mice had significantly smaller ductal trees. ^P < 0.05, independent samples t test, comparing control and flutamide-treated female mice. LN, lymph node. Table 3. UOG-MIX Treatment Did not Alter Growth Parameters in Prepubertal Mammary Gland Variable  Ductal Extension (mm) a  Number of Branching Points  Control (n = 6)  −13.64 ± 1.07  25.0 ± 2.2  MIX-3 (n = 8)  −14.19 ± 0.96  24.9 ± 2.1  MIX-30 (n = 5)  −12.96 ± 0.43  20.8 ± 2.3  MIX-300 (n = 6)  −15.74 ± 1.46  25.7 ± 2.5  MIX-3000 (n = 6)  −13.88 ± 0.34  26.2 ± 2.4  Flutamide (n = 9)  −14.03 ± 0.84  20.9 ± 1.4  Variable  Ductal Extension (mm) a  Number of Branching Points  Control (n = 6)  −13.64 ± 1.07  25.0 ± 2.2  MIX-3 (n = 8)  −14.19 ± 0.96  24.9 ± 2.1  MIX-30 (n = 5)  −12.96 ± 0.43  20.8 ± 2.3  MIX-300 (n = 6)  −15.74 ± 1.46  25.7 ± 2.5  MIX-3000 (n = 6)  −13.88 ± 0.34  26.2 ± 2.4  Flutamide (n = 9)  −14.03 ± 0.84  20.9 ± 1.4  a Negative values for ductal extension indicate ductal trees that have not yet grown past the central lymph node. View Large Adult mammary gland morphology was altered by developmental exposure to UOG chemical mixtures Increased volume of ducts and epithelial compartment Unbiased stereological methods revealed treatment-related effects on mammary gland morphology in early adulthood (PND85; Fig. 2A). Females exposed to MIX-300 had significantly more ducts compared with vehicle-treated controls (P < 0.05, Fisher post hoc test; Fig. 2A and 2B), and a trend for an increase was also seen in females exposed to MIX-30 and MIX-3000 (P < 0.1, Fisher post hoc test; Fig. 2A and 2B). The volume fraction of total mammary epithelium was significantly increased in females exposed to MIX-300 and MIX-3000 (P < 0.05, Fisher post hoc test; Fig. 2A and 2C). Neither of these parameters were significantly altered by flutamide treatment. Alveolar buds were not observed in females of any treatment group (Fig. 2A; data not shown). Figure 2. View largeDownload slide Prenatal UOG-MIX treatment induced increased epithelial density in female mammary glands at adulthood. (A) Example of whole mount mammary glands (magnification ×13) from control, MIX-300, and MIX-3000 treatment groups on PND85. Vehicle-treated females had the least dense mammary epithelium. Note that alveolar buds were not observed in any treatment group. (B) Volume fraction of ducts and (C) volume fraction of all epithelium was increased in UOG-MIX–treated groups. *P < 0.05, Fisher post hoc test; aP < 0.1, Fisher post hoc test. Figure 2. View largeDownload slide Prenatal UOG-MIX treatment induced increased epithelial density in female mammary glands at adulthood. (A) Example of whole mount mammary glands (magnification ×13) from control, MIX-300, and MIX-3000 treatment groups on PND85. Vehicle-treated females had the least dense mammary epithelium. Note that alveolar buds were not observed in any treatment group. (B) Volume fraction of ducts and (C) volume fraction of all epithelium was increased in UOG-MIX–treated groups. *P < 0.05, Fisher post hoc test; aP < 0.1, Fisher post hoc test. Increased proliferation/apoptosis ratio The growth of the mammary epithelium is dependent on a balance of proliferation (to extend ductal structures into the mammary fat pad) and apoptosis (to produce hollow ducts capable of transporting milk) (36). To quantify apoptosis in the epithelium of the adult female mammary glands, TUNEL staining was used to compare the controls, MIX-3, MIX-3000, and flutamide-treated mice. No statistically significant changes in TUNEL incorporation were seen in MIX-3 and MIX-3000 females (Fig. 3A and 3B). A borderline statistically significant decrease was seen in TUNEL-positive cells in the flutamide-treated group (P = 0.052, independent samples t test; Fig. 3B). Figure 3. View largeDownload slide Apoptosis, proliferation, and proliferation/apoptosis ratios were altered in UOG-MIX– and flutamide-treated females. (A) Examples of TUNEL and immunohistochemistry testing for Ki67 and ERα in mammary gland sections from control and MIX-3–, MIX-3000–, and flutamide-treated females. Positive cells indicated by red arrows. Scale bar = 20 μm. (B) Quantification of TUNEL incorporation in epithelium from vehicle-, MIX-3–, MIX-3000–, and flutamide-treated females revealed a dose-dependent decrease in the percentage of epithelial cells undergoing apoptosis. (C) Quantification of Ki67, a marker of proliferation, in epithelial cells. (D) Substantial alterations to the ratio of cell proliferation and cell death were observed in both MIX-3 and MIX-3000 treatment groups and the flutamide-treated group. (E) Quantification of the percentage of mammary epithelial cells expressing ERα. *P < 0.05, Fisher post hoc test; δP < 0.1, independent samples t test, comparing control and flutamide-treated females. Figure 3. View largeDownload slide Apoptosis, proliferation, and proliferation/apoptosis ratios were altered in UOG-MIX– and flutamide-treated females. (A) Examples of TUNEL and immunohistochemistry testing for Ki67 and ERα in mammary gland sections from control and MIX-3–, MIX-3000–, and flutamide-treated females. Positive cells indicated by red arrows. Scale bar = 20 μm. (B) Quantification of TUNEL incorporation in epithelium from vehicle-, MIX-3–, MIX-3000–, and flutamide-treated females revealed a dose-dependent decrease in the percentage of epithelial cells undergoing apoptosis. (C) Quantification of Ki67, a marker of proliferation, in epithelial cells. (D) Substantial alterations to the ratio of cell proliferation and cell death were observed in both MIX-3 and MIX-3000 treatment groups and the flutamide-treated group. (E) Quantification of the percentage of mammary epithelial cells expressing ERα. *P < 0.05, Fisher post hoc test; δP < 0.1, independent samples t test, comparing control and flutamide-treated females. We next evaluated proliferation in the mammary epithelium using antibodies for Ki67, a marker of proliferation. Although the proliferation levels were low, as expected for adult mammary glands, we observed statistically significant effects of the 23-chemical mixture on the number of cells expressing Ki67 in the MIX-3 group (P < 0.05; Fisher post hoc test; Fig. 3A and 3C). Females from the MIX-3 group had 427% more Ki67-positive cells compared with the controls. The MIX-3000 treatment group had 54% more Ki67-positive cells compared with the controls, although these differences were not statistically significant. Ki67 expression was not altered by the use of flutamide (Fig. 3C). We evaluated the ratio of proliferation/apoptosis in the control, MIX-3, MIX-3000, and flutamide-treated females. We found a striking 397% increase in the proliferation/apoptosis ratio in the MIX-3 treatment group (P < 0.05; Fisher post hoc test) and nonsignificant increases in the MIX-3000 group (131%) compared with the controls (Fig. 3D). No effect on the proliferation/apoptosis ratio was observed in the flutamide-treated females (Fig. 3D). ERα expression tended to be associated with developmental UOG-MIX treatment EDCs have been shown to not only bind to hormone receptors, but also to alter expression of hormone receptors in a dose-, age-, and tissue-specific manner (37). We next asked whether developmental exposure to fracking chemicals would alter expression of ERα in the mammary epithelium. We observed a general increase in ERα expression with an increasing UOG-MIX dose, although these differences were not statistically significant (Fig. 3A and 3E). Flutamide also did not affect the percentage of epithelial cells expressing ERα (Fig. 3E). TEB-like intraductal hyperplasia in mammary glands after developmental mix treatment A striking observation made during the morphological assessment of the adult (PND85) mammary glands was the appearance of structures that resembled TEBs (Fig. 4). These TEB-like structures were not observed in the control females (Table 4). The unusual structures were excised from whole mount mammary glands and further characterized using histological and immunohistochemistry tools (Fig. 5). Hematoxylin and eosin staining revealed ducts with excessive layers of epithelial cells, consistent with intraductal hyperplasia (Fig. 5A). These TEB-like lesions were highly proliferative: >40% of epithelial cells in lesions were Ki67-positive compared with <4% of epithelial cells in normal ducts collected from the same mice or ducts excised from control females (Fig. 5A and 5B). ERα expression was also greater in many, although not all, of the excised lesions (Fig. 5A and 5C). Collectively, these results suggest that the retained TEB-like structures are highly proliferative intraductal hyperplasias that are likely to be estrogen-responsive. Figure 4. View largeDownload slide TEB-like structures observed in whole mount mammary glands from prenatally UOG-MIX–treated females on PND85. (A) Micrograph of whole mount mammary gland from a control mouse with typical terminal ducts. These blunt ends (arrowheads) are the normal structures observed at the ends of ducts. (B) These examples illustrate the TEB-like structures excised from whole mounts from UOG-MIX– and flutamide-treated females. Arrows indicate abnormal structures. Treatment groups for the individual female mice are indicated. LN, lymph node. Figure 4. View largeDownload slide TEB-like structures observed in whole mount mammary glands from prenatally UOG-MIX–treated females on PND85. (A) Micrograph of whole mount mammary gland from a control mouse with typical terminal ducts. These blunt ends (arrowheads) are the normal structures observed at the ends of ducts. (B) These examples illustrate the TEB-like structures excised from whole mounts from UOG-MIX– and flutamide-treated females. Arrows indicate abnormal structures. Treatment groups for the individual female mice are indicated. LN, lymph node. Table 4. Presence of TEB-Like Structures in Control or UOG-MIX–Treated Females Variable  TEB-Like Structures, %  Control  0 (n = 9)  MIX-3  29 (n = 7)  MIX-30  17 (n = 6)  MIX-300  25 (n = 8)  MIX-3000  40 (n = 5)  Flutamide  22 (n = 9)  Variable  TEB-Like Structures, %  Control  0 (n = 9)  MIX-3  29 (n = 7)  MIX-30  17 (n = 6)  MIX-300  25 (n = 8)  MIX-3000  40 (n = 5)  Flutamide  22 (n = 9)  View Large Figure 5. View largeDownload slide Histological and immunohistochemical evaluation of TEB-like structures excised from UOG-MIX–treated females. (A) Representative lesions from three MIX groups (lesion 1, MIX-30; lesion 2, MIX-3; lesion 3, MIX-300) were evaluated using hematoxylin and eosin (H&E) staining and immunohistochemistry for Ki67 (a marker of proliferation) and ERα. (B) Quantification of Ki67 in excised TEB-like structures compared with Ki67 expression in other regions of the same glands. (C) Quantification of ERα in excised TEB-like structures compared with ERα expression in other regions of the same glands. Figure 5. View largeDownload slide Histological and immunohistochemical evaluation of TEB-like structures excised from UOG-MIX–treated females. (A) Representative lesions from three MIX groups (lesion 1, MIX-30; lesion 2, MIX-3; lesion 3, MIX-300) were evaluated using hematoxylin and eosin (H&E) staining and immunohistochemistry for Ki67 (a marker of proliferation) and ERα. (B) Quantification of Ki67 in excised TEB-like structures compared with Ki67 expression in other regions of the same glands. (C) Quantification of ERα in excised TEB-like structures compared with ERα expression in other regions of the same glands. Discussion We examined the effects of exposure to a mixture of 23 chemicals used in unconventional oil and gas extraction, most of which were previously shown to exhibit antagonistic properties on one or more hormone receptors [(13) and Table 1], including the ER, progesterone receptor, and AR. Also, 21, 20, and 11 of these chemicals were previously demonstrated to antagonize human ER, AR, and PR, respectively, in a reporter gene assay in human endometrial cells (summarized in Table 1). To the best of our knowledge, we have shown for the first time that the mouse mammary gland is sensitive to developmental UOG-MIX exposure, with dose-specific effects on tissue morphology (after exposure to MIX-300 and MIX-3000), cell proliferation (after exposure to MIX-3), and the induction of unique intraductal hyperplasias (observed in all MIX groups). The effects observed in the present study occurred at low doses; the two lowest dose groups (MIX-3 and MIX-30) are equivalent to the concentrations measured in drinking water in regions experiencing drilling and the highest dose group (MIX-3000) is equivalent to the concentrations of many UOG-MIX components measured in industry wastewater (4, 13, 18, 38). In addition, the concentrations for several of the 23 chemicals in the UOG-MIX used in the present study have not yet been determined in either drinking water or wastewater. The mammary gland is a hormone-sensitive organ that is responsive to multiple endocrine inputs during early development. Testosterone has a unique role in establishing the sexually dimorphic development of the mouse mammary gland (26, 27) but is not thought to play a role in the female or postnatally. AR expression remains high in the mammary stroma until birth (26, 39), suggesting that UOG-MIX exposure could affect mammary development via actions at the AR. Although some endpoints were similarly affected (e.g., retention of TEBs), other effects of UOG-MIX exposure we observed in the present study were distinct from the effects of flutamide, suggesting that the UOG mixture might not working through an antiandrogenic mechanism as we originally hypothesized. To date, few studies have evaluated the effects of prenatal exposures to antiandrogenic chemicals on the female mammary gland; most studies of antiandrogens have examined their effects on male offspring, including nipple retention and other disruptions to mammary gland morphology [reviewed in (40–42)]. These studies suggested that evaluations of mammary glands in male mice exposed to UOG mixtures are a priority. A better mechanistic understanding of the effects of AR antagonists on the female mammary gland, including a wider range of exposure to flutamide, is also needed. Previous studies of other EDCs have shown few effects on the female gland before puberty (33, 43), suggesting that the most obvious effects of EDCs manifest visibly only after the onset of ovarian hormone production. In support of this, we observed striking effects of developmental exposure to UOG chemical mixtures in the adult female mammary gland (at PND85). Not only were mammary glands more developed in the UOG-MIX–treated adult females, as indicated by the increased volume fraction of epithelium (Fig. 2), but UOG-MIX treatment also altered the ratios of proliferation and apoptosis (Fig. 3), cell parameters important for dictating the growth and function of the mammary gland (33, 44, 45). Normal growth of the mammary epithelium is dependent on a balance of proliferation (to extend ductal structures into the mammary fat pad) and apoptosis (to produce hollow ducts capable of transporting milk) (36); thus, disruptions to these cellular features could predispose animals to mammary gland diseases (e.g., cancer) or abnormal function (e.g., disruptions to lactation). Additional studies are needed to evaluate these outcomes in UOG-MIX–treated females. The substantial alterations to the apoptosis/proliferation ratio suggest that the mammary glands from UOG-MIX–treated females might continue to manifest complex hyperplastic and preneoplastic lesions in later adulthood. Additional studies are needed to evaluate this possibility. We were also surprised to see TEB-like hyperplastic lesions in the mammary glands collected from mice in the UOG-MIX–treatment groups, as well as in the flutamide-treated mice (Fig. 4). Other EDCs have also been shown to induce development of intraductal hyperplasias, although these typically will have a “beaded duct,” rather than a TEB-like, appearance (34, 46, 47). TEBs are a characteristic structure of glands undergoing puberty (45, 48); these highly proliferative structures drive the growth of the mammary epithelium into the surrounding fat pad. Once the epithelial tree is fully formed, the TEB structures recede and are not seen in adult glands (19, 49). One possibility is that UOG-MIX exposure during early life delays the appearance of TEBs, and thus their presence in the glands of adult mice is indicative of a shift in the timing of puberty. Although the timing of vaginal opening and the age of first vaginal estrus were not affected by UOG-MIX exposure (18), the timing of pubertal growth in the mammary gland involves distinct events (50). Alternatively, it is possible that the timing of mammary puberty is unaffected by UOG-MIX exposure and that the retention of TEBs is indicative of failure to progress to blunt ductal ends. Thus, the presence of TEBs could be interpreted as diminished development of the mammary gland in UOG-MIX–treated groups, in contrast to the effects of these chemical mixtures on epithelial density, which are more consistent with advanced development of the gland (Fig. 2). EDCs have previously been shown to produce competing effects on growth parameters in the mammary gland, including some that advance one aspect of development and seem to delay other developmental landmarks (51). Studies of bisphenol A, for example, have shown that developmental exposure can both decrease ductal extension (growth) and increase the size of TEBs (33, 44, 52–54). These results are consistent with the two competing roles that hormones such as estrogen can have in the developing gland, including its ability to promote proliferation of some mammary epithelial cells and inducing apoptosis of other mammary epithelial cells inside the duct, allowing the lumen to form (55, 56). TEBs are one of the sites where cancers are thought to arise; thus, delays in TEB recession could increase the gland’s sensitivity to carcinogens (51, 57). Future studies are needed to determine whether these TEBs are retained into later adulthood and whether prenatal UOG-MIX exposure increases the sensitivity of the gland to carcinogens. The mammary gland provides an in vivo tool to evaluate the effects of EDCs with different modes of action, including chemicals administered as mixtures. To date, in vivo studies of EDC chemical mixtures have been limited, with even fewer studies examining mammals [examples include (58–60)]. The results from in vitro studies have suggested that compounds with a similar mode of action can have additive effects, including some xenoestrogen mixtures, which have been described as producing “something from nothing” (61–63). In the present study, the assessment of the mixture of 23 chemicals in human endometrial cells demonstrated synergistic responses for ER and thyroid receptor antagonism and less than additive effects for AR and glucocorticoid receptor antagonism (13). Other groups have found that the toxic effects of pesticides are compounded when examined as formulations compared with studying only the so-called active ingredient (64), suggesting that novel insights can be gained from examining chemical mixtures that would not be detected, or even anticipated, from examining the mixture’s components individually. Concerns have been raised that any effects observed after exposure to mixtures will be difficult to evaluate mechanistically because they cannot be attributed to any single component of the chemical mixture unless all individual components are also evaluated. However, it is worth noting that humans are exposed to chemical mixtures, rather than single compounds; thus, the study of mixtures in laboratory animals might provide better understanding of the human condition (65). To date, >1000 different chemicals have been reported for use during unconventional oil and gas extraction (5). Previous studies have shown that prenatal exposure to this same 23-chemical mixture induced adverse health outcomes in the male offspring (the brothers of the female mice examined in the present study) (13). Furthermore, other endocrine organs were affected in the female mice examined in the present study, in addition to alterations to serum hormone concentrations, including decreased levels of the pituitary hormones luteinizing hormone, prolactin, and follicle-stimulating hormone, among others (18). Decreased pituitary hormone concentrations could influence mammary gland health, including the development of the gland during pregnancy, an endpoint that deserves future attention. Determining whether mixtures of fracking chemicals affect human populations is an important goal, especially as the number of drilling sites continues to increase (66). A recent systematic review evaluating the strength of the data for the association between conventional oil and gas and UOG operations and human reproductive outcomes found moderate evidence for an increased risk of preterm birth, miscarriage, birth defects, decreased semen quality, and prostate cancer (67). The evidence for an association between UOG operations and breast cancer, or other diseases of the breast, remains inadequate. The results from our study suggest that longitudinal studies evaluating women exposed to UOG chemical mixtures during early life are needed to address this data gap. Conclusion Future studies are needed to evaluate the many additional chemicals used in, and produced by, UOG processes to better quantify the concentrations of these and other contaminants in environmental samples and to assess the effects of exposure during other sensitive windows of development, including pregnancy and lactation, puberty, and the aging female. Future studies should also evaluate whether developmental UOG-MIX treatments can sensitize animals to hormones or carcinogens, as would be expected in females with retained TEBs. For mechanistic insights, additional examination of some of the individual components in the fracking chemical mixture might be warranted. Abbreviations: AR androgen receptor EDC endocrine-disrupting chemical ER estrogen receptor PND postnatal day TEB terminal end bud TUNEL terminal deoxynucleotidyltransferase-mediated dUTP nick end labeling UOG unconventional oil and gas UOG-MIX 23-chemical mixture of unconventional oil and gas compounds. Acknowledgments The authors thank the other members of the Vandenberg and Nagel Laboratories for feedback on our study. Financial Support: This work was supported by the National Institute of Environmental Health Sciences of the National Institutes of Health (Award K22-ES025811 to L.N.V. and Grants R21-ES026395 and R01-ES021394-04S1 to S.C.N.). The content of this manuscript is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health. Disclosure Summary: L.N.V. has received travel reimbursement from universities, governments, nongovernmental organizations, and industry to speak about endocrine-disrupting chemicals. The remaining authors have nothing to disclose. References 1. Wilson DG. Energy supplies and future engines for land, sea, and air. J Air Waste Manag Assoc . 2012; 62( 6): 607– 624. Google Scholar CrossRef Search ADS PubMed  2. Vidic RD, Brantley SL, Vandenbossche JM, Yoxtheimer D, Abad JD. Impact of shale gas development on regional water quality. Science . 2013; 340( 6134): 1235009. Google Scholar CrossRef Search ADS PubMed  3. Vengosh A, Jackson RB, Warner N, Darrah TH, Kondash A. A critical review of the risks to water resources from unconventional shale gas development and hydraulic fracturing in the United States. Environ Sci Technol . 2014; 48( 15): 8334– 8348. Google Scholar CrossRef Search ADS PubMed  4. Kassotis CD, Iwanowicz LR, Akob DM, Cozzarelli IM, Mumford AC, Orem WH, Nagel SC. Endocrine disrupting activities of surface water associated with a West Virginia oil and gas industry wastewater disposal site. Sci Total Environ . 2016; 557-558: 901– 910. Google Scholar CrossRef Search ADS PubMed  5. Elliott EG, Ettinger AS, Leaderer BP, Bracken MB, Deziel NC. A systematic evaluation of chemicals in hydraulic-fracturing fluids and wastewater for reproductive and developmental toxicity. J Expo Sci Environ Epidemiol . 2017; 27( 1): 90– 99. Google Scholar CrossRef Search ADS PubMed  6. Wylie S, Schultz K, Thomas D, Kassotis C, Nagel S. Inspiring collaboration: the legacy of Theo Colborn’s transdisciplinary research on fracking [published online ahead of print September 13, 2016]. New Solut . doi:1048291116666037. 7. Webb E, Bushkin-Bedient S, Cheng A, Kassotis CD, Balise V, Nagel SC. Developmental and reproductive effects of chemicals associated with unconventional oil and natural gas operations. Rev Environ Health . 2014; 29( 4): 307– 318. Google Scholar CrossRef Search ADS PubMed  8. Kassotis CD, Tillitt DE, Davis JW, Hormann AM, Nagel SC. Estrogen and androgen receptor activities of hydraulic fracturing chemicals and surface and ground water in a drilling-dense region. Endocrinology . 2014; 155( 3): 897– 907. Google Scholar CrossRef Search ADS PubMed  9. Bolden AL, Kwiatkowski CF, Colborn T. New look at BTEX: are ambient levels a problem? Environ Sci Technol . 2015; 49( 9): 5261– 5276. Google Scholar CrossRef Search ADS PubMed  10. Colborn T, Kwiatkowski C, Schultz K, Bachran M. Natural gas operations from a public health perspective. Hum Ecol Risk Assess . 2011; 17( 5): 1039– 1056. Google Scholar CrossRef Search ADS   11. Zoeller RT, Brown TR, Doan LL, Gore AC, Skakkebaek NE, Soto AM, Woodruff TJ, Vom Saal FS. Endocrine-disrupting chemicals and public health protection: a statement of principles from The Endocrine Society. Endocrinology . 2012; 153( 9): 4097– 4110. Google Scholar CrossRef Search ADS PubMed  12. Cozzarelli IM, Skalak KJ, Kent DB, Engle MA, Benthem A, Mumford AC, Haase K, Farag A, Harper D, Nagel SC, Iwanowicz LR, Orem WH, Akob DM, Jaeschke JB, Galloway J, Kohler M, Stoliker DL, Jolly GD. Environmental signatures and effects of an oil and gas wastewater spill in the Williston Basin, North Dakota. Sci Total Environ . 2017; 579: 1781– 1793. Google Scholar CrossRef Search ADS PubMed  13. Kassotis CD, Klemp KC, Vu DC, Lin CH, Meng CX, Besch-Williford CL, Pinatti L, Zoeller RT, Drobnis EZ, Balise VD, Isiguzo CJ, Williams MA, Tillitt DE, Nagel SC. Endocrine-disrupting activity of hydraulic fracturing chemicals and adverse health outcomes after prenatal exposure in male mice. Endocrinology . 2015; 156( 12): 4458– 4473. Google Scholar CrossRef Search ADS PubMed  14. Vandenberg LN, Ehrlich S, Belcher SM, Ben-Jonathan N, Dolinoy DC, Hugo ER, Hunt PA, Newbold RR, Rubin BS, Saili KS, Soto AM, Wang HS, Vom Saal FS. Low dose effects of bisphenol A: an integrated review of in vitro, laboratory animal and epidemiology studies. Endocr Disruptors (Austin) . 2013; 1: e25078. Google Scholar CrossRef Search ADS   15. Rochester JR. Bisphenol A and human health: a review of the literature. Reprod Toxicol . 2013; 42: 132– 155. Google Scholar CrossRef Search ADS PubMed  16. Environmental Protection Agency. Assessment of the potential impacts of hydraulic fracturing for oil and gas on drinking water resources. Available at: http://www.epa.gov/hfstudy2015. Accessed 7 September 2017. 17. Rudel RA, Attfield KR, Schifano JN, Brody JG. Chemicals causing mammary gland tumors in animals signal new directions for epidemiology, chemicals testing, and risk assessment for breast cancer prevention. Cancer . 2007; 109( 12 suppl) 2635– 2666. Google Scholar CrossRef Search ADS PubMed  18. Kassotis CD, Bromfield JJ, Klemp KC, Meng CX, Wolfe A, Zoeller RT, Balise VD, Isiguzo CJ, Tillitt DE, Nagel SC. Adverse reproductive and developmental health outcomes following prenatal exposure to a hydraulic fracturing chemical mixture in female C57Bl/6 mice. Endocrinology . 2016; 157( 9): 3469– 3481. Google Scholar CrossRef Search ADS PubMed  19. Rudel RA, Fenton SE, Ackerman JM, Euling SY, Makris SL. Environmental exposures and mammary gland development: state of the science, public health implications, and research recommendations. Environ Health Perspect . 2011; 119( 8): 1053– 1061. Google Scholar CrossRef Search ADS PubMed  20. Fenton SE. The mammary gland: a tissue sensitive to environmental exposures. Rev Environ Health . 2009; 24( 4): 319– 325. Google Scholar CrossRef Search ADS PubMed  21. Soto AM, Vandenberg LN, Maffini MV, Sonnenschein C. Does breast cancer start in the womb? Basic Clin Pharmacol Toxicol . 2008; 102( 2): 125– 133. Google Scholar CrossRef Search ADS PubMed  22. Brisken C, Ataca D. Endocrine hormones and local signals during the development of the mouse mammary gland. Wiley Interdiscip Rev Dev Biol . 2015; 4( 3): 181– 195. Google Scholar CrossRef Search ADS PubMed  23. Saji S, Jensen EV, Nilsson S, Rylander T, Warner M, Gustafsson JA. Estrogen receptors alpha and beta in the rodent mammary gland. Proc Natl Acad Sci USA . 2000; 97( 1): 337– 342. Google Scholar CrossRef Search ADS PubMed  24. Bocchinfuso WP, Korach KS. Mammary gland development and tumorigenesis in estrogen receptor knockout mice. J Mammary Gland Biol Neoplasia . 1997; 2( 4): 323– 334. Google Scholar CrossRef Search ADS PubMed  25. Mallepell S, Krust A, Chambon P, Brisken C. Paracrine signaling through the epithelial estrogen receptor alpha is required for proliferation and morphogenesis in the mammary gland. Proc Natl Acad Sci USA . 2006; 103( 7): 2196– 2201. Google Scholar CrossRef Search ADS PubMed  26. Kratochwil K. Development and loss of androgen responsiveness in the embryonic rudiment of the mouse mammary gland. Dev Biol . 1977; 61( 2): 358– 365. Google Scholar CrossRef Search ADS PubMed  27. Kratochwil K, Schwartz P. Tissue interaction in androgen response of embryonic mammary rudiment of mouse: identification of target tissue for testosterone. Proc Natl Acad Sci USA . 1976; 73( 11): 4041– 4044. Google Scholar CrossRef Search ADS PubMed  28. Axelstad M, Christiansen S, Boberg J, Scholze M, Jacobsen PR, Isling LK, Kortenkamp A, Hass U. Mixtures of endocrine-disrupting contaminants induce adverse developmental effects in preweaning rats. Reproduction . 2014; 147( 4): 489– 501. Google Scholar CrossRef Search ADS PubMed  29. Yeh S, Hu Y-C, Wang P-H, Xie C, Xu Q, Tsai M-Y, Dong Z, Wang R-S, Lee T-H, Chang C. Abnormal mammary gland development and growth retardation in female mice and MCF7 breast cancer cells lacking androgen receptor. J Exp Med . 2003; 198( 12): 1899– 1908. Google Scholar CrossRef Search ADS PubMed  30. Soto AM, Brisken C, Schaeberle C, Sonnenschein C. Does cancer start in the womb? Altered mammary gland development and predisposition to breast cancer due to in utero exposure to endocrine disruptors. J Mammary Gland Biol Neoplasia . 2013; 18( 2): 199– 208. Google Scholar CrossRef Search ADS PubMed  31. Plante I, Stewart MK, Laird DW. Evaluation of mammary gland development and function in mouse models. J Vis Exp . 2011; ( 53): 2828. 32. Suvorov A, Vandenberg LN. To cull or not to cull? Considerations for studies of endocrine-disrupting chemicals. Endocrinology . 2016; 157( 7): 2586– 2594. Google Scholar CrossRef Search ADS PubMed  33. Markey CM, Luque EH, Munoz De Toro M, Sonnenschein C, Soto AM. In utero exposure to bisphenol A alters the development and tissue organization of the mouse mammary gland. Biol Reprod . 2001; 65( 4): 1215– 1223. Google Scholar CrossRef Search ADS PubMed  34. Vandenberg LN, Maffini MV, Schaeberle CM, Ucci AA, Sonnenschein C, Rubin BS, Soto AM. Perinatal exposure to the xenoestrogen bisphenol-A induces mammary intraductal hyperplasias in adult CD-1 mice. Reprod Toxicol . 2008; 26( 3-4): 210– 219. Google Scholar CrossRef Search ADS PubMed  35. Haseman JK, Bailer AJ, Kodell RL, Morris R, Portier K. Statistical issues in the analysis of low-dose endocrine disruptor data. Toxicol Sci . 2001; 61( 2): 201– 210. Google Scholar CrossRef Search ADS PubMed  36. Soto AM, Sonnenschein C. The two faces of Janus: sex steroids as mediators of both cell proliferation and cell death. J Natl Cancer Inst . 2001; 93( 22): 1673– 1675. Google Scholar CrossRef Search ADS PubMed  37. Yin W, Maguire SM, Pham B, Garcia AN, Dang NV, Liang J, Wolfe A, Hofmann HA, Gore AC. Testing the critical window hypothesis of timing and duration of estradiol treatment on hypothalamic gene networks in reproductively mature and aging female rats. Endocrinology . 2015; 156( 8): 2918– 2933. Google Scholar CrossRef Search ADS PubMed  38. Kassotis CD, Tillitt DE, Lin CH, McElroy JA, Nagel SC. Endocrine-disrupting chemicals and oil and natural gas operations: potential environmental contamination and recommendations to assess complex environmental mixtures. Environ Health Perspect . 2016; 124( 3): 256– 264. Google Scholar PubMed  39. Wasner G, Hennermann I, Kratochwil K. Ontogeny of mesenchymal androgen receptors in the embryonic mouse mammary gland. Endocrinology . 1983; 113( 5): 1771– 1780. Google Scholar CrossRef Search ADS PubMed  40. Skarda J. Bioassay of steroid hormone agonist and antagonist activities of anti-androgens on mammary gland, seminal vesicles and spleen of male mice. J Vet Med A Physiol Pathol Clin Med . 2003; 50( 4): 204– 212. Google Scholar CrossRef Search ADS PubMed  41. Boberg J, Christiansen S, Axelstad M, Kledal TS, Vinggaard AM, Dalgaard M, Nellemann C, Hass U. Reproductive and behavioral effects of diisononyl phthalate (DINP) in perinatally exposed rats. Reprod Toxicol . 2011; 31( 2): 200– 209. Google Scholar CrossRef Search ADS PubMed  42. Christiansen S, Boberg J, Axelstad M, Dalgaard M, Vinggaard AM, Metzdorff SB, Hass U. Low-dose perinatal exposure to di(2-ethylhexyl) phthalate induces anti-androgenic effects in male rats. Reprod Toxicol . 2010; 30( 2): 313– 321. Google Scholar CrossRef Search ADS PubMed  43. Mandrup K, Boberg J, Isling LK, Christiansen S, Hass U. Low-dose effects of bisphenol A on mammary gland development in rats. Andrology . 2016; 4( 4): 673– 683. Google Scholar CrossRef Search ADS PubMed  44. Markey CM, Coombs MA, Sonnenschein C, Soto AM. Mammalian development in a changing environment: exposure to endocrine disruptors reveals the developmental plasticity of steroid-hormone target organs. Evol Dev . 2003; 5( 1): 67– 75. Google Scholar CrossRef Search ADS PubMed  45. McNally S, Martin F. Molecular regulators of pubertal mammary gland development. Ann Med . 2011; 43( 3): 212– 234. Google Scholar CrossRef Search ADS PubMed  46. Murray TJ, Maffini MV, Ucci AA, Sonnenschein C, Soto AM. Induction of mammary gland ductal hyperplasias and carcinoma in situ following fetal bisphenol A exposure. Reprod Toxicol . 2007; 23( 3): 383– 390. Google Scholar CrossRef Search ADS PubMed  47. Durando M, Kass L, Piva J, Sonnenschein C, Soto AM, Luque EH, Muñoz-de-Toro M. Prenatal bisphenol A exposure induces preneoplastic lesions in the mammary gland in Wistar rats. Environ Health Perspect . 2007; 115( 1): 80– 86. Google Scholar CrossRef Search ADS PubMed  48. Paine IS, Lewis MT. The terminal end bud: the little engine that could. J Mammary Gland Biol Neoplasia . 2017; 22( 2): 93– 108. Google Scholar CrossRef Search ADS PubMed  49. Hinck L, Silberstein GB. Key stages in mammary gland development: the mammary end bud as a motile organ. Breast Cancer Res . 2005; 7( 6): 245– 251. Google Scholar CrossRef Search ADS PubMed  50. Hermes GL, McClintock MK. Isolation and the timing of mammary gland development, gonadarche, and ovarian senescence: implications for mammary tumor burden. Dev Psychobiol . 2008; 50( 4): 353– 360. Google Scholar CrossRef Search ADS PubMed  51. Macon MB, Fenton SE. Endocrine disruptors and the breast: early life effects and later life disease. J Mammary Gland Biol Neoplasia . 2013; 18( 1): 43– 61. Google Scholar CrossRef Search ADS PubMed  52. Muñoz-de-Toro M, Markey CM, Wadia PR, Luque EH, Rubin BS, Sonnenschein C, Soto AM. Perinatal exposure to bisphenol-A alters peripubertal mammary gland development in mice. Endocrinology . 2005; 146( 9): 4138– 4147. Google Scholar CrossRef Search ADS PubMed  53. Wadia PR, Vandenberg LN, Schaeberle CM, Rubin BS, Sonnenschein C, Soto AM. Perinatal bisphenol A exposure increases estrogen sensitivity of the mammary gland in diverse mouse strains. Environ Health Perspect . 2007; 115( 4): 592– 598. Google Scholar CrossRef Search ADS PubMed  54. Moral R, Wang R, Russo IH, Lamartiniere CA, Pereira J, Russo J. Effect of prenatal exposure to the endocrine disruptor bisphenol A on mammary gland morphology and gene expression signature. J Endocrinol . 2008; 196( 1): 101– 112. Google Scholar CrossRef Search ADS PubMed  55. Bresciani F. Effect of ovarian hormones on duration of DNA synthesis in cells of the C3H mouse mammary gland. Exp Cell Res . 1965; 38( 1): 13– 32. Google Scholar CrossRef Search ADS PubMed  56. Williams JM, Daniel CW. Mammary ductal elongation: differentiation of myoepithelium and basal lamina during branching morphogenesis. Dev Biol . 1983; 97( 2): 274– 290. Google Scholar CrossRef Search ADS PubMed  57. Tucker DK, Macon MB, Strynar MJ, Dagnino S, Andersen E, Fenton SE. The mammary gland is a sensitive pubertal target in CD-1 and C57Bl/6 mice following perinatal perfluorooctanoic acid (PFOA) exposure. Reprod Toxicol . 2015; 54: 26– 36. Google Scholar CrossRef Search ADS PubMed  58. Hass U, Boberg J, Christiansen S, Jacobsen PR, Vinggaard AM, Taxvig C, Poulsen ME, Herrmann SS, Jensen BH, Petersen A, Clemmensen LH, Axelstad M. Adverse effects on sexual development in rat offspring after low dose exposure to a mixture of endocrine disrupting pesticides. Reprod Toxicol . 2012; 34( 2): 261– 274. Google Scholar CrossRef Search ADS PubMed  59. Buñay J, Larriba E, Moreno RD, Del Mazo J. Chronic low-dose exposure to a mixture of environmental endocrine disruptors induces microRNAs/isomiRs deregulation in mouse concomitant with intratesticular estradiol reduction. Sci Rep . 2017; 7( 1): 3373. Google Scholar CrossRef Search ADS PubMed  60. Sobolewski M, Conrad K, Allen JL, Weston H, Martin K, Lawrence BP, Cory-Slechta DA. Sex-specific enhanced behavioral toxicity induced by maternal exposure to a mixture of low dose endocrine-disrupting chemicals. Neurotoxicology . 2014; 45: 121– 130. Google Scholar CrossRef Search ADS PubMed  61. Rajapakse N, Silva E, Kortenkamp A. Combining xenoestrogens at levels below individual no-observed-effect concentrations dramatically enhances steroid hormone action. Environ Health Perspect . 2002; 110( 9): 917– 921. Google Scholar CrossRef Search ADS PubMed  62. Rajapakse N, Silva E, Scholze M, Kortenkamp A. Deviation from additivity with estrogenic mixtures containing 4-nonylphenol and 4-tert-octylphenol detected in the E-SCREEN assay. Environ Sci Technol . 2004; 38( 23): 6343– 6352. Google Scholar CrossRef Search ADS PubMed  63. Silva E, Rajapakse N, Kortenkamp A. Something from “nothing”—eight weak estrogenic chemicals combined at concentrations below NOECs produce significant mixture effects. Environ Sci Technol . 2002; 36( 8): 1751– 1756. Google Scholar CrossRef Search ADS PubMed  64. Defarge N, Takács E, Lozano VL, Mesnage R, Spiroux de Vendômois J, Séralini GE, Székács A. Co-formulants in glyphosate-based herbicides disrupt aromatase activity in human cells below toxic levels. Int J Environ Res Public Health . 2016; 13( 3): 13. 65. Kortenkamp A. Low dose mixture effects of endocrine disrupters: implications for risk assessment and epidemiology. Int J Androl . 2008; 31( 2): 233– 240. Google Scholar CrossRef Search ADS PubMed  66. Czolowski ED, Santoro RL, Srebotnjak T, Shonkoff SBC. Toward consistent methodology to quantify populations in proximity to oil and gas development: a national spatial analysis and review. Environ Health Perspect . 2017; 125( 8): 086004. Google Scholar CrossRef Search ADS PubMed  67. Balise VD, Meng CX, Cornelius-Green JN, Kassotis CD, Kennedy R, Nagel SC. Systematic review of the association between oil and natural gas extraction processes and human reproduction. Fertil Steril . 2016; 106( 4): 795– 819. Google Scholar CrossRef Search ADS PubMed  Copyright © 2018 Endocrine Society TI - Prenatal Exposure to Unconventional Oil and Gas Operation Chemical Mixtures Altered Mammary Gland Development in Adult Female Mice JO - Endocrinology DO - 10.1210/en.2017-00866 DA - 2018-03-01 UR - https://www.deepdyve.com/lp/oxford-university-press/prenatal-exposure-to-unconventional-oil-and-gas-operation-chemical-jXjDXiTudS SP - 1277 EP - 1289 VL - 159 IS - 3 DP - DeepDyve ER -