Scholarly article on topic 'Exposure to an environmentally relevant mixture of brominated flame retardants affects fetal development in Sprague-Dawley rats'

Exposure to an environmentally relevant mixture of brominated flame retardants affects fetal development in Sprague-Dawley rats Academic research paper on "Environmental engineering"

CC BY-NC-ND
0
0
Share paper
Academic journal
Toxicology
OECD Field of science
Keywords
{"Environmental contaminant" / "House dust" / "Polybrominated diphenyl ethers" / Mixture / "In utero" / "Skeletal development"}

Abstract of research paper on Environmental engineering, author of scientific article — Robert G. Berger, Pavine L.C. Lefèvre, Sheila R. Ernest, Michael G. Wade, Yi-Qian Ma, et al.

Abstract Brominated flame retardants are incorporated into a wide variety of consumer products and are known to enter into the surrounding environment, leading to human exposure. There is accumulating evidence that these compounds have adverse effects on reproduction and development in humans and animal models. Animal studies have generally characterized the outcome of exposure to a single technical mixture or congener. Here, we determined the impact of exposure of rats prior to mating and during gestation to a mixture representative of congener levels found in North American household dust. Adult female Sprague-Dawley rats were fed a diet containing 0, 0.75, 250 or 750mg/kg of a mixture of flame retardants (polybrominated diphenyl ethers, hexabromocyclododecane) from two weeks prior to mating to gestation day 20. This formulation delivered nominal doses of 0, 0.06, 20 and 60mg/kg body weight/day. The lowest dose approximates high human exposures based on house dust levels and the dust ingestion rates of toddlers. Litter size and resorption sites were counted and fetal development evaluated. No effects on maternal health, litter size, fetal viability, weights, crown rump lengths or sex ratios were detected. The proportion of litters with fetuses with anomalies of the digits (soft tissue syndactyly or malposition of the distal phalanges) was increased significantly in the low (0.06mg/kg/day) dose group. Skeletal analysis revealed a decreased ossification of the sixth sternebra at all exposure levels. Thus, exposure to an environmentally relevant mixture of brominated flame retardants results in developmental abnormalities in the absence of apparent maternal toxicity. The relevance of these findings for predicting human risk is yet to be determined.

Academic research paper on topic "Exposure to an environmentally relevant mixture of brominated flame retardants affects fetal development in Sprague-Dawley rats"

lllll^^^B AMIN E IN PRESS

Toxicology xxx (2014) xxx-xxx

Contents lists available at ScienceDirect

ELSEVIER

Toxicology

journal homepage www.elsevier.com/locate/toxicol

Exposure to an environmentally relevant mixture of brominated flame retardants affects fetal development in Sprague-Dawley rats

Robert G. Berger3, Pavine L.C. Lefèvrea, Sheila R. Ernest3, Michael G. Wadec, Yi-Qian Maa, Dorothea F.K. Rawnd, Dean W. Gaertnerd, Bernard Robaireab, Barbara F. Halesa *

a Department of Pharmacology and Therapeutics, McGill University, Montreal, QC, Canada H3G 1Y6 b Department of Obstetrics and Gynecology, McGill University, Montreal, QC, Canada H3G 1Y6

c Environmental Health Science and Research Bureau, Healthy Environments and Consumer Safety Branch, Health Canada, Ottawa, ON, Canada K1A 0K9 d Food Research Division, Bureau of Chemical Safety, Health Products and Food Branch, Health Canada, Ottawa, ON, Canada K1A OL2

ABSTRACT

Brominated flame retardants are incorporated into a wide variety of consumer products and are known to enter into the surrounding environment, leading to human exposure. There is accumulating evidence that these compounds have adverse effects on reproduction and development in humans and animal models. Animal studies have generally characterized the outcome of exposure to a single technical mixture or congener. Here, we determined the impact of exposure of rats prior to mating and during gestation to a mixture representative of congener levels found in North American household dust. Adult female Sprague-Dawley rats were fed a diet containing 0,0.75,250 or 750 mg/kg of a mixture of flame retardants (polybrominated diphenyl ethers, hexabromocyclododecane) from two weeks prior to mating to gestation day 20. This formulation delivered nominal doses of 0,0.06, 20 and 60 mg/kg body weight/day. The lowest dose approximates high human exposures based on house dust levels and the dust ingestion rates of toddlers. Litter size and resorption sites were counted and fetal development evaluated. No effects on maternal health, litter size, fetal viability, weights, crown rump lengths or sex ratios were detected. The proportion of litters with fetuses with anomalies of the digits (soft tissue syndactyly or malposition of the distal phalanges) was increased significantly in the low (0.06mg/kg/day) dose group. Skeletal analysis revealed a decreased ossification of the sixth sternebra at all exposure levels. Thus, exposure to an environmentally relevant mixture of brominated flame retardants results in developmental abnormalities in the absence of apparent maternal toxicity. The relevance of these findings for predicting human risk is yet to be determined.

© 2014 Published by Elsevier Ireland Ltd.

ARTICLE INFO

Article history:

Received 7 August 2013

Received in revised form 12 March 2014

Accepted 14 March 2014

Available online xxx

Keywords:

Environmental contaminant House dust

Polybrominated diphenyl ethers

Mixture

In utero

Skeletal development

1. Introduction

Brominated flame retardants (BFRs) are incorporated into a wide variety of household and commercial goods to reduce flame propagation and ignition rates (Camino et al., 1991). Many BFRs do not form covalent bonds with the polymer matrices in products to which they are added and leach out into the environment. Recently, a major class of BFRs, the polybrominated diphenyl ethers (PBDEs), were withdrawn from commerce in North America and Europe

Abbreviations: BFRs, brominated flame retardants; PBDEs, polybrominated diphenyl ethers; HBCD, hexabromocyclododecane.

* Corresponding author at: Department of Pharmacology and Therapeutics, McGill University, 3655 Promenade Sir William Osler, Montreal, QC, Canada H3G 1Y6. Tel.: +1 514 398 3610; fax: +1 514 398 7120.

E-mail address: barbara.hales@mcgill.ca (B.F. Hales).

http://dx.doi.org/10.1016/j.tox.2014.03.005 0300-483X/© 2014 Published by Elsevier Ireland Ltd.

due to their environmental persistence and bioaccumulative properties (EC, 2011; Tullo, 2003). Despite this, products containing PBDEs, such as furniture or electronic devices, continue to be used (Stapleton et al., 2012). In spite of their being banned in North America, analyses of house dust show that BFRs are present in the majority of North American homes; PBDEs and hexabromocyclododecane (HBCD) are prominent among the BFRs that are detected (Allen et al., 2008; Stapleton et al., 2008) and have not declined appreciably since being removed from commerce (Dodson et al., 2012); recent studies demonstrate that BFRs are still found in the serum of pregnant women (Zota et al., 2013).

Brominated flame retardants are detected in various human tissues, with the highest concentrations detected in North Americans (Hites, 2004; Sjödin et al., 2004; Toms et al., 2009). Primary routes of exposure have been identified as inhalation and ingestion of dust (Johnson et al., 2010; Stapleton et al., 2005; Wu et al., 2007) and dietary consumption (Fraser et al., 2009; Wu et al., 2007).

lllll^^^B AIIIIII E IN PRESS

2 R.G. Berger et al. / Toxicology xxx (2014) xxx-xxx

Recent investigations have shown that BFRs are detected in serum (Johnson et al., 2010; Meeker et al., 2009), hair (Aleksa et al., 2012), adipose tissue (Johnson-Restrepo et al., 2005), breast milk and ovarian follicular fluid (Johnson et al., 2012; Petro et al., 2012; Marchitti et al., 2013). Children have the highest serum BFR concentrations, with peak measurements up to 2-10-fold higher than in adults (Rose et al., 2010). Increased BFR levels in infants and toddlers may be a result of additional exposure through hand to mouth activities (Stapleton et al., 2008, 2012) and ingestion of contaminated breast milk (Carrizo et al., 2007; Siddique et al., 2012). Prenatal development is also vulnerable to BFR exposure through placental transfer (Mazdai et al., 2003). In the Montreal area, increasing levels of PBDEs were found in human fetal livers from 1998 through 2006, mirroring trends in adults and children (Doucet et al., 2009).

Exposure to BFRs during in utero development is of significant concern. Indeed, maternal PBDE levels during pregnancy have been associated with deficits in children's attention, fine motor coordination, and cognitive functioning (Eskenazi et al., 2013). PBDEs are structurally similar to thyroid hormone and may disrupt thyroid homeostasis (Darnerud et al., 2001; Ernest et al., 2012; Fowles et al., 1994; Hallgren et al., 2001; Blake et al., 2011). Alternatively, BFRs may act as endocrine disruptors by altering estrogen (Ceccatelli et al., 2006; Meerts et al., 2001) or androgen signaling (Stoker et al., 2005), or perturbing insulin-like growth factor 1 (IGF1) signaling (Haave et al., 2011; Suvorov et al., 2009; Suvorov andTakser, 2010). In addition to their effects on the brain and developing nervous system (Dreiem et al., 2010; Saegusa et al., 2009; van der Ven et al., 2008; Ta et al., 2011), exposure to BFRs has been associated with detrimental effects on the male and female reproductive systems (Talsness et al., 2007; van der Ven et al., 2008; Ema et al., 2008; Lilienthal et al., 2006; Zhang et al., 2013), bone (van der Ven et al., 2009; Blanco et al., 2012), the immune system (van der Ven et al., 2009; Liu et al., 2012; Bondy et al., 2013), and the liver (Blanco et al., 2012; Ema et al., 2008). The consequences of exposure to single BFR congeners, such as PBDE-47 (Talsness et al., 2008; Ta et al., 2011) or PBDE-99 (Talsness et al., 2005), or to BFR technical mixtures, such as DE-71 (Blake et al., 2011) have been reported. For example, in utero exposure to PBDE-99 delayed bone ossification and enlarged the liver (Blanco et al., 2012), decreased fertility in females (Talsness et al., 2005), and altered steroid hormone production and the onset of puberty (Lilienthal et al., 2006). Bone development was altered in animals exposed to HBCD from premating until weaning, with a decrease in trabecular bone mineral density of the tibia in females (van der Ven et al., 2009).

Although humans are exposed to a complex mixture of BFRs, animal studies have focused on the effects of a single congener or technical mixture. Here, we determined the effects of an environmentally relevant BFR mixture, typically found in household dust, on female fertility, the establishment and maintenance of pregnancy, and fetal development in Sprague-Dawley rats.

2. Materials and methods

2.1. BFR mixture formulation

Formulation of the BFR mixture was described previously (Ernest et al., 2012). Briefly, three technical PBDE mixtures (DE-71, DE-79 and BDE-209) and one HBCD mixture were combined to yield a ratio of PBDE congeners and HBCD comparable to the median levels observed in Boston house dust (Allen et al., 2008; Stapleton et al., 2008). This BFR mixture was incorporated into an isoflavone-free diet (Teklad Global 2019 diet; Harlan Laboratories, Madison, WI) with 4.3 g/kg corn oil. Diets were formulated to contain 0, 0.75, 250 or 750 mg of BFR mixture/kg. Diet samples from each experimental condition were collected for later analysis of BFR content via gas

chromatography/mass spectrometry (GC/MS), as described below. The diet formulations were intended to deliver nominal doses of

0, 0.06, 20 and 60 mg/kg body weight/day, assuming a daily food consumption of 80 g/kg body weight/day. The lowest dose was estimated to be a close approximation of maximum human exposure, based on a dust ingestion rate of 100mg/day in children (16.5 kg body weight) and the scaling of dose from humans to rodents (1:6.9, human to rat body surface area ratio) (Allen et al., 2008; Stapleton et al., 2008).

2.2. GC/MS analysis of dietary BFR content

Feed samples were thawed and ground to a fine powder, thoroughly mixed and weighed (~0.1 g) into Erlenmeyer flasks. Each sample was fortified with surrogate standards (twelve 13C12 labeled PBDE congeners: 15, 28, 47, 77, 99,100,126,138,153,154, 183, 209 and 13C12 labeled a-, P- and 7-HBCD) (Cambridge Isotope Laboratories, Andover, MA and Wellington Laboratories, Guelph, ON, Canada). Distilled water (20 mL) was added to each sample and allowed to soak for approximately 30 min, followed by the addition of 50 mL acetone:hexane (2:1, v/v). Samples were homogenized for approximately 30 s using a Silverson homogenizer. Solid material was removed from the extract by transferring the extract to a sep-aratory funnel through a funnel packed with glass wool. Aqueous sodium chloride (50 mL) was added to each separatory funnel, and the funnel was gently shaken and vented as necessary. Following phase separation, the aqueous layer was drained from the funnel, discarded and the extract was transferred to a round bottomed flask (rbf) after passing through a bed of anhydrous sodium chloride. Each sample was concentrated to approximately 1 mL using rotary evaporation and quantitatively transferred to a clean pre-weighed vial which was placed, uncapped, in the fume hood until a stable weight was achieved, to allow for lipid content determination.

Samples were then re-dissolved in hexane (1-2 mL) with gentle swirling and added to a column containing 2% deactivated Florisil (6g), topped with approximately 1 cm of each of glass beads and anhydrous sodium sulfate (Na2SO4). Once the extract was washed onto the column and drained to within a few millimeters of the top of the Na2SO4, enough hexane was added to wet the entire column (~5 mL). PBDEs were eluted with 70 mL hexane and collected in a 250 mL rbf; the flask was then removed and retained. A second rbf was placed under the Florisil column and HBCD isomers were eluted using 70 mL dichloromethane (DCM): hexane (30:70, v/v). Both fractions were reduced in volume to ~1 mL using rotary evaporation. The PBDE extract was transferred to a v-notch vial and evaporated to dryness using a gentle stream of nitrogen, re-diluted in isooctane and mixed. This final isooctane extract was transferred to a chromatographic vial for analysis. The HBCD fraction was similarly transferred to a v-notch vial, but taken only to near dryness (~50 |L) using a gentle stream of nitrogen prior to the addition of C122H18Br6 analogs of a-, P- and 7-HBCD. The samples were then placed in a fume hood where they remained until dryness was attained. To each dry HBCD extract, 100 |L of methanol:water (80:20, v/v) was added, the vial was mixed and the final extract was transferred to a chromatography vial. Dilution of feed sample extracts, to which higher levels of PBDEs and HBCD had been added (e.g., 20 ppm, 60 ppm nominal concentrations), was necessary to ensure that accurate concentration measurements could be achieved.

PBDE analysis was performed using an Agilent 6890 gas chromatograph coupled to a 5973N mass selective detector (Agilent, Mississauga, ON, Canada). A 15 m J&W DB5 MS (0.25 mm

1.d. x 0.10 |im film thickness) column was used for separation of the PBDEs, following introduction of 2 |L sample volumes into the cool on-column injection system which was set to track the oven temperature. The initial temperature was set to 80 ° C and held for 1 min,

III AIIIIII E IN PRESS

R.G. Berger et al. / Toxicology xxx (2014) xxx-xxx 3

increased to 225 °C at a rate of 32 °C min-1, followed by an increase at a rate of 3 ° C min-1 to 230 ° C and held for 1 min. The final increase was set at a rate of 40 °C min-1 to 315 °C and held for 7 min. Helium was used as the carrier gas set at constant flow (1.2 mL min-1) for all PBDE analyses and the head pressure initially was 27 kPa. The mass spectrometer was operating in electron impact ionization mode and the source temperature was 250°C.

HBCD analysis was performed using an Agilent 1100 (Agilent Technologies, Mississauga, ON, Canada) high pressure liquid chromatograph coupled to a Quattro Ultima triple quadrupole MS/MS (Waters Corporation, Milford, MA) with electrospray ionization in the negative ion detection mode using a Phenomenex Kinetex C18 column (2.1 mm x 150 mm, 2.6 |im) (Phenomenex, Torrance, CA). Water (mobile phase A) and acetonitrile:methanol (2:1) (mobile phase B) were used to separate the HBCD isomers and the gradient was as follows: 40% mobile phase B for 1 min, 40-80% phase B for 4 min, 80-90% phase B for 13 min, then returned to 40% phase B over 0.5 min where it remained until 18 min. The flow rate was maintained at 0.22 mLmin-1 and the column temperature at 30 °C, to completely resolve the deuterated HBCD analogs from the 13C-labeled and native analogs. The capillary and cone voltage were -3.0 kV and 35 V, respectively. The source temperature and des-olvation temperature were 140 °C and 350°C, respectively. Cone gas flow and desolvation gas flow were 77 Lh-1 and 723 Lh-1, respectively. Argon was the collision gas set at 5 x 10-3 mbar and resolution was established at 90% valley at base for both quadrupole analyzers. Dwell times were set to 5 ms.

All reported sample concentrations were corrected for recoveries using isotope dilution. Two reagent blanks, an aliquot of corn oil known to be free of PBDEs and HBCD fortified with known concentrations of these analytes and a reference material of fish certified to contain known amounts of PBDEs were included with each set of samples prepared for analysis. Laboratory background was accounted for by removing the average of the PBDE and HBCD concentrations observed in the two reagent blank samples for each set of samples analyzed.

Average PBDE surrogate recoveries ranged from 71% to 81% (PBDE 183 and PBDE 47, respectively) and 13C HBCD recoveries ranged from 84% to 97% (a- and P-HBCD, respectively) in the diet samples. PBDE concentrations were within the expected concentration range in the reference material tested and the average oil spike recovery was 100% for PBDEs (86-109%) and 101% for HBCD (100-102%) in the fortified oil.

2.3. Animals and treatment

All animal studies were conducted in accordance with the procedures and principles outlined in the Guide to the Care and Use of Experimental Animals prepared by the Canadian Council on Animal Care (McGill Animal Research Centre protocol 5862). Virgin female Sprague-Dawley rats (200-250 g) were obtained from Charles River (St-Constant, QC, Canada). Animals were individually housed at the Animal Resources Centre of McGill University in rooms maintained at 20 °C on a 12h:12h light:dark cycle. Food and water were provided ad libitum. Following one week of acclimation to the control diet, female rats were randomly assigned to one of four experimental conditions (n = 35-38/group) and fed a diet supplemented with 0, 0.75, 250 or 750 mg BFR mixture/kg for two to four weeks prior to mating. During this exposure period, estrus cyclicity was assessed daily by analyzing vaginal cytology, as previously described (Goldman et al., 2007). After this premat-ing exposure, females in proestrus were caged with proven breeder male Sprague-Dawley rats overnight. The next morning was designated as gestation day (GD) 0 if spermatozoa were detected in the vaginal smear. Following insemination, females continued on their

respective diets until euthanasia on GD 20. Throughout the pre-mating period, dam physical examinations were undertaken and body weight and food intake were monitored once a week. Following mating, body weight and food intake were assessed on GDs 3, 8,13,16 and 20.

2.4. Tissue collection

Dams were euthanized by CO2 asphyxiation followed by exsan-guination via cardiac puncture on GD 20. The sequence in which dams from each treatment group were euthanized was random and investigators were blinded as to the treatment group. Liver, kidneys, spleen, thymus, heart, lung, uterus and ovaries were removed and weighed. Uterine horns were inspected for implantation and resorption sites. The trachea, including the thyroid gland, was removed gently and fixed in neutral buffered formalin for morphological analysis. Whole blood was transferred to a Vacutainer SST (BD Biosciences Canada, Mississauga, ON, Canada), allowed to clot for 30 min at room temperature, then held on ice for no longer than 4 h until centrifugation at 1300 x g for 20 min. Serum was aliquoted and stored at -80 °C.

2.5. Thyroid gland histology

Fixed thyroid glands, in situ on the trachea, were embedded in paraffin, sectioned (5 |im) and stained with periodic acid Schiff. The height of the thyroid gland epithelium was measured using digital analysis in sections known to be in the interior of the gland. Epithelial height measurements were taken from at least 5 separate follicles per image, selected using an overlaid, numbered grid and 5 randomly generated numbers; at least 10 images per animal were analyzed. All measurements were made by a single observer blinded to the animal treatment (n = 10/group).

2.6. Serum measures

Calcium, magnesium, phosphorus, blood urea nitrogen, uric acid, albumin, total protein, creatinine kinase, creatinine, triglycerides, cholesterol, glucose, alkaline phosphatase (ALP) and alanine aminotransferase (ALT) were measured in serum with an ABX Pentra 400 clinical chemistry analyzer (Horiba ABX, Montpellier, France). Total serum thyroxine (T4) was assayed using commercial RIA kits following the manufacturer's protocol (MP Biomedicals, Lachine, QC, Canada).

2.7. Fetal examinations

The position of fetuses in the uterine horns was recorded prior to their excision from the yolk sac and separation from the placenta. Placental weights and diameters were recorded. All live fetuses were evaluated for gender, body weight, crown-rump length and anogenital distance index (AGDI), calculated by dividing the anogenital distance (mm) by the cubed root of the body mass (Gallavan et al., 1999). All fetuses were examined for external physical malformations (n = 255, 365, 285 and 283 for the control, 0.06, 20 and 60 mg/kg/day BFR doses, respectively); an average of two fetuses/litter (n = 18, 41, 26, and 38 for the control, 0.06, 20 and 60 mg/kg/day BFR doses, respectively) were prepared for the evaluation of skeletal development. Accordingly, treatment-related effects on external morphology were analyzed on the basis of the proportion of litters with affected fetuses whereas skeletal evaluations were based on the proportion of individual fetuses affected.

Fetuses were euthanized by hypothermia prior to fixation and storage in 95% ethanol. The preparation and staining methods were adapted from Whitaker and Dix (1979) with minor alterations. Briefly, at the time of processing, fixed fetuses were immersed in

III AIIIIII E IN PRESS

4 R.G. Berger et al. / Toxicology xxx (2014) xxx-xxx

Table 1

Dietary BFR content.

Congener (^g/g) Nominal BFR mixture dose (mg/kg/day)

Control (0) 0.06 20 60

BDE-15 NDa ND ND ND

BDE-17 ND ND ND ND

BDE-28 0.0001 ± 0.00001 0.0009 ±0.00003 0.195 ±0.0281 0.445 ± 0.0109

BDE-37 ND ND ND ND

BDE-47 0.001 ± 0.0002 0.098 ± 0.0013 32.99 ± 3.86 100.35 ±1.53

BDE-66 ND 0.001 ± 0.00004 0.357 ± 0.0391 1.135 ±0.027

BDE-71 ND ND ND ND

BDE-75 ND ND ND ND

BDE-77 ND ND ND ND

BDE-85 ND 0.006 ±0.0001 1.331 ±0.124 3.360 ±0.102

BDE-99 0.002 ±0.0003 0.124 ±0.0014 44.65 ±4.95 136.19 ±3.37

BDE-100 0.001 ± 0.0001 0.031 ± 0.0003 11.19± 1.26 32.94 ±0.86

BDE-119 ND ND ND ND

BDE-126 ND ND ND ND

BDE-138 ND 0.002 ± 0.0004 0.475 ±0.051 1.662 ±0.048

BDE-153 0.003 ±0.0015 0.011 ± 0.0002 3.62 ±0.388 10.66 ±0.16

BDE-154 0.001 ± 0.00001 0.011 ± 0.0002 3.78 ±0.416 11.36 ±0.20

BDE-160 ND ND ND ND

BDE-181 ND ND ND ND

BDE-183 ND 0.001 ± 0.00004 0.506 ±0.0798 1.230 ±0.051

BDE-190 ND ND ND ND

BDE-205 ND ND ND ND

BDE-209 ND 0.362 ± 0.0141 145.77 ±11.47 361.90 ±6.27

HBCD-a 0.0009 ±0.0003 0.018 ± 0.0003 0.304 ±0.778 19.43 ±0.32

HBCD-p 0.0001 ±4 x10-6 0.001 ± 0.00005 0.381 ±0.252 1.194 ±0.013

HBCD-7 0.0003 ± 0.0006 0.0012 ±0.00007 0.468 ±0.0519 1.486 ±0.03

J>BDE(|xg/g) 0.015 ± 0.001 0.647 ± 0.033 245.68 ±12.14 661.23 ±33.19

VHBCD^g/g) 0.001 ± 0.0002 0.020 ±0.006 7.15 ±1.96 22.11 ±6.03

^2^PBDE nominal dose NA 88.87 101.13 91.11

Number of samples 24 18 18 18

Values are expressed as mean ± SEM.

a PBDE limits of detection ranged from 0.03 ngg-1 to 6.5 ngg-1 (PBDE15 and 209, respectively) and HBCD limits of detection ranged from 8.9 to 9.1 ngg-1 (p- and a-HBCD, respectively).

a water bath (50 °C) for 5 min, skinned, and cervical and dorsal muscles were gently removed. Fetuses were placed in the alcian blue/alizarin red stain solution for 24 h at 37 °C, followed by 1% potassium hydroxide for 48 h, with fresh 1% potassium hydroxide added at 24 h, and then 0.5% potassium hydroxide for an additional 24 h.The potassium hydroxide was then replaced with a 2:2:1 solution consisting of two parts 70% ethanol, two parts glycerine, and one part benzyl alcohol. After 24 h, stained skeletons were placed in a 1:1 solution (70% ethanol:glycerine) for evaluation and storage. The skull, sternebrae, vertebrae, ribs, pectoral and pelvic girdles, fore and hind limbs and paws were examined.

2.8. Statistical analyses

Data were analyzed by one-way ANOVA followed by post hoc Dunnett's test to compare means to control. When tests for assumptions of homogeneity of variance and normality failed, data were log transformed and retested. When homoscedasticity and normality were still not satisfied after this transformation, data were retested using Kruskal-Wallis (K-W) ANOVA on ranks. Analyses of external (n = litter) and skeletal (n = fetus) malformations were conducted using a 1-tail Chi-Square test to compare the proportions of affected fetuses. The level of significance was p <0.05.

3. Results

3.1. BFR contents of the diets

Low levels of BDE-28, -47, -99, -100, -153 and -154 and the HBCD isomers were found in the control diet (Table 1). The relative levels

of BDE-47, -99, -100, -153, -154, and -209 did not differ across the experimental groups; they were determined to be approximately 15%, 19%, 5%, 1.6% 1.7% and 57%, respectively, of the total mixture and were consistent throughout the experimental period. These congeners accounted for 99% of the PBDEs measured in the dietary mixture. Values of total PBDE were 88-101% of the nominal doses (Table 1).

3.2. Diet consumption during treatment

Ingestion of the BFR mixture during the exposure period did not alter dam food consumption or weight gain (Table 2). Since the BFR dietary mixture remained constant, there was a decrease in the average BFR dose ingested with increased weight gain during gestation (Table 2). Thus, BFR ingestion was calculated to be on average 0.05,17.4 and 47.5 mg BFR/kg/day or 74-86% of target doses for the three experimental conditions.

3.3. Effects of BFR exposure on dams

Dam survival and appearance were not affected by ingestion of the BFR mixture. BFR exposure did not affect absolute tissue weights (Supplementary data Table S1); however, it did result in a significant increase in relative liver and kidney weights at the highest dose compared to controls (Table 3). Analysis of variance showed a significant effect of BFR exposure on serum T4 levels (p = 0.047) although post hoc tests did not indicate a significant difference from controls (Table 4). No significant differences in the inner thyroid epithelial follicle height were observed among the treatment groups (control: 7.9 ±0.4 |im; low dose: 7.1 ±0.2 |im;

I'l'llB^^^M AIIIHI.E IN PRESS

R.G. Berger et al. / Toxicology xxx (2014) xxx-xxx

Table 2

Consumption of the BFR diet during gestation.

Parameter

Gestational day

Nominal BFR mixture dose (mg/kg/day)

Control (0) 0.06a 20 60

0-3 24.25 ±0.81 (13) 22.91 ± 0.84 (19) 25.32 ± 0.96 (14) 23.32 ± 0.69 (17)

Food consumption (g/day) 3-8 27.63 ±0.96(13) 25.53 ± 1.04 (18) 27.95 ± 1.22 (14) 28.69 ± 1.26 (16)

8-13 27.70 ±1.16(14) 27.76 ± 0.84 (22) 28.23 ± 0.84 (17) 25.82 ± 1.18 (17)

13-16 26.27 ±0.97(13) 25.77 ± 1.00 (17) 24.29 ± 0.60 (13) 25.81 ± 0.89 (13)

16-20 25.48 ±1.11 (16) 25.09 ± 1.05 (22) 24.03 ± 0.76 (15) 24.33 ± 0.88 (17)

3 332.31 ±5.10(14) 330.37 ± 5.83 (19) 333.95 ± 7.57 (14) 320.73 ± 6.38 (16)

8 358.50 ±5.46(14) 350.13 ± 5.87 (21) 354.06 ± 7.77 (16) 341.12 ± 6.89 (17)

Gestational weight (g) 13 390.70 ±5.88(13) 383.55 ± 5.95 (19) 382.30 ± 8.33 (15) 370.55 ± 9.17 (15)

16 413.77 ±5.10(16) 402.80 ± 6.40 (20) 399.49 ± 8.85 (12) 391.19 ± 8.14 (16)

20 462.41 ±7.51 (16) 447.43 ± 7.16 (21) 446.74 ± 7.49 (17) 441.81 ± 8.72 (17)

3 0 0.046 19.17 49.68

Approximate BFR 8 0 0.049 19.96 57.47

dose (mg BFR/kg 13 0 0.048 18.67 47.62

BW) 16 0 0.043 15.37 45.08

20 0 0.037 13.60 37.64

Values are expressed as mean ± SEM of (n) observations per treatment.

a The maximum measured content of house dust PBDEs is 192 |ig per g dust in the studies that we chose as a basis for our BFR mixture (Allen et al., 2008). This represents a dose of 1.16 |ig PBDE per kg/day if a toddler (16.5 kg) ingests 100 mg dust. When we apply allometric scaling (x6.9) to adjust for the difference in body weight to surface area between humans and rats, our calculated rat equivalent to human dose is 8.0 |ig/kg/day. Thus, our low dose group represents an exposure that is approximately 5 times the estimated high exposure for toddlers.

Table 3

Effects of BFR exposure on dam organ weights.

Nominal BFR mixture treatment dose (mg/kg/day)

Control (0) 0.06 20 60

n 15 21 17 17

Liver 34.2 ± 0.9 33.5 ± 0.7 35.3 ±0.8 37.9 ±0.7"

Kidney 2.15 ±0.04 2.23 ±0.04 2.30 ±0.05 2.34 ±0.04*

Spleen 1.41 ±0.06 1.44 ±0.05 1.47 ±0.08 1.43 ±0.04

Thymus 0.83 ±0.05 0.72 ±0.03 0.73 ± 0.04 0.79 ±0.05

Heart 2.71 ±0.06 2.69 ±0.03 2.64 ±0.04 2.66 ±0.04

Lung 2.71 ±0.08 2.66 ±0.04 2.82 ±0.07 2.77 ±0.10

Values are expressed as mean tissue weight (g) normalized to body weight (kg) ± SEM.

* p <0.05 compared with control.

** p <0.01 compared with control.

Table 4

Effects of BFR exposure on dam serum biochemistry.

Parameter Nominal BFR mixture dose (mg/kg/day)

Control (0) 0.06 20 60

n 14 21 17 17

Thyroxin (T4) (|xg/dL) 1.52 ±0.12 1.59 ±0.09 1.68 ±0.13 1.30 ±0.08

Calcium (mg/dL) 12.38 ±0.24 12.06 ±0.21 11.92 ±0.17 12.07 ±0.22

Magnesium (mg/dL) 3.23 ±0.07 3.19 ±0.07 3.25 ±0.06 3.07 ± 0.07

Phosphorus (mg/dL) 10.40 ±0.45 10.21 ± 0.63 8.84 ±0.26 9.04 ± 0.20

Blood urea N (mg/dL) 15.3 ±0.7 14.6 ±0.5 13.5 ±0.6 14.1 ±0.4

Uric acid (mg/dL) 3.49 ± 0.27 2.91 ± 0.21 3.18 ±0.15 2.98 ± 0.20

Albumin (g/dL) 2.78 ±0.04 2.72 ±0.04 2.77 ±0.05 2.72 ±0.07

Total protein (g/dL) 6.52 ±0.11 6.49 ±0.19 6.21 ±0.11 5.97 ±0.15*

Creatinine kinase (U/L) 311 ± 33 319 ±42 354 ±35 308 ± 31

Creatinine (mg/dL) 0.53 ±0.02 0.50 ±0.02 0.52 ±0.01 0.52 ± 0.01

Triglycerides (mg/dL) 446 ±52 463 ±55 329 ±34 292 ±17

Cholesterol (mg/dL) 102± 4 92 ± 3 95 ± 3 109 ±4

Glucose (mg/Dl) 125 ±6 119± 6 123 ±3 114± 5

ALP (U/L) 89.2 ±14.1 83.7 ±9.1 111.4± 18.1 94.8 ± 10.5

ALT (U/L) 66.0 ± 3.2 61.3 ±3.0 56.0 ±2.5 63.5 ± 3.3

n = no. of dams. Values are expressed as mean ± SEM.

* p <0.05 compared with control.

349 middle dose: 8.1 ±0.3 |im; high dose: 9.0 ±0.4 |im). BFR expo-

350 sure also had minimal effects on the other serum parameters

351 assessed in the dams (Table 4). In the high dose BFR treatment

352 group total protein levels were significantly reduced compared

353 to controls. Although ANOVA indicated an influence of BFR expo-

354 sure on phosphorus (p = 0.029), total cholesterol (p = 0.012), and

triglycerides (p = 0.017), post hoc tests failed to reveal significant 355

differences from controls. Neither estrous cycle length nor fer- 356

tility measures were affected by BFR exposure (Table 5). Finally, 357

BFR treatment did not affect the mean numbers of implanta- 358

tion and resorption sites, postimplantation loss, or the litter 359

size. 360

IMllB^^^B ARIIII E IN PRESS

6 R.G. Berger et al. / Toxicology xxx (2014) xxx-xxx

Table 5

Effects of BFR exposures on measures of dam fertility.

Parameter Nominal BFR mixture dose (mg/kg/day)

Control (0) 0.06 20 60

Estrous cycle length - days 4.8 ±0.4 (31) 4.6 ±0.3(29) 4.7 ± 0.2 (29) 5.6 ±0.6(29)

Mating indexa 70.4 (27) 85.2 (27) 75.0 (24) 70.8 (24)

Fecundity indexb 79.0(15) 91.3(21) 94.4(17) 100.0(17)

No. ofimplants per litter 17.9 ±0.4(15) 18.4 ±0.4(21) 18.4 ±0.9(17) 17.7 ±0.7(17)

No. of resorptions per litter 0.9 ±0.4(15) 1.0 ±0.3(21) 1.5 ±0.4(17) 1.1 ±0.3(17)

Postimplantation lossesc 4.8 ±1.8(15) 5.4 ±1.5(21) 8.2 ±1.6(17) 6.8 ± 3.0(17)

No. of dead fetuses per litter 0(15) 0.05 ±0.05(21) 0.06 ±0.06(17) 0(17)

No. of live fetuses per litter 17.0 ±0.3(15) 17.4 ±0.4(21) 16.8 ±0.80(17) 16.7 ±0.9(17)

Values are expressed as mean ± SEM with the total number of observations. a Mating index = (no. of sperm positive females/no. of mated females) x 100. b Fecundity index = (no. of pregnant females/no. of sperm positive females) x 100. c Postimplantation losses = [(no. of implants - no. of viable fetuses)/no. of implants] x 10.

3.4. Fetal developmental measures

BFR exposure had no effect on placental weights (Table 6). Fetal body weights were not affected, although a significant decrease in crown rump length in female fetuses was observed in the 20 mg/kg/day BFR treatment group. Neither sex ratios nor the anogenital distance indices in male and female fetuses were affected.

External fetal examinations revealed malpositioned, fused and spread digits, crooked paws and single observations of paw edema (control fetus) and a shortened tail (60 mg/kg/day fetus) in control and BFR exposed fetuses (Fig. 1A-C). Fetuses with malpositioned or fused digits were observed in 25-52% of the litters (Fig. 1D). Analysis ofthe proportion of litters with fetuses with malpositioned or fused digits revealed a significant increase above control only in the low dose treatment group (Fig. 1D). Among BFR exposed fetuses, 80-84% of digit defects were found in the hind paws, compared to 67% in control fetuses (data not shown).

A number of skeletal variations, including offset ossification of sternebrae, unossified sternebrae and vertebrae, and incompletely fused ossified centrums, were observed (Fig. 2A-D). A significant increase in the incidence of fetuses with any skeletal variation was found in both the middle (p = 0.018) and high (p = 0.024) doses (Fig. 2E). The mean numbers of affected fetuses per litter did not differ between control and treated groups; however, affected fetuses were found in multiple litters (Supplementary data, Fig. S1A). A decrease or lack of ossification in the sternebrae was the most sensitive skeletal variation observed after BFR exposure; the incidence of fetuses with unossified sternebrae was significantly increased in all BFR treatment groups (Table 7, Supplementary data, Fig. S1B). There were no significant changes in ossification or cartilage development in other areas of the skeleton, including the skull, phalanges, carpals or metacarpals. Skeletal analysis of the paws revealed no significant changes in ossification or cartilage development (Supplementary data, Fig. S2).

Supplementary Figs. SI and SII related to this article found, in the online version, at http://dx.doi.org/10.1016Zj.tox.2014.03.005.

4. Discussion

The current study demonstrated that exposure during pregnancy to environmentally relevant complex BFR mixture, with relative congener levels similar to those found in household dust, significantly affected fetal development. This is the first investigation to report that exposure to BFRs increased fetal digit malformations. In addition, these data are the first to demonstrate an increase in the incidence of fetal skeletal variations in the absence of any indications of maternal toxicity following exposure to an environmentally relevant BFR dose (0.06 mg/kg/day).

Effects on any parameter of maternal health were observed only in dams exposed to the high BFR dose (60 mg/kg/day); these dams had significant liver and kidney enlargement and a reduction in serum protein compared to controls. Increases in liver and kidney weights are indicative of toxicity; however, there were no changes in serum biomarkers that are typically associated with liver and kidney damage. Our findings are consistent with a previous investigation in which increased liver weight and drug metabolizing enzyme activity were observed in GD 20 dams following exposure to DE-71 (30 mg/kg/day) (Zhou et al., 2002). Adult males fed the same complex BFR mixture used here showed an increase in both liver and kidney weights following exposure to BFRs at 20 mg/kg/day (Ernest et al., 2012). Others have also reported that exposures to PBDE congeners, as well as HBCD, are associated with increased kidney (van der Ven et al., 2006) and liver (Ema et al., 2008; Stoker et al., 2004; van der Ven et al., 2006, 2008) weights.

We report here that exposure of female rats prior to conception and during gestation (for approximately 35 days) to 60 mg/kg/day of our BFR mixture did not affect either serum thyroxin (T4) concentrations or thyroid epithelial cell height in the dams. In contrast, exposure of adult males to this BFR mixture (20 mg/kg/day for 70 days) induced a significant decrease in serumT4 levels (Ernest et al., 2012). Serum T4 concentrations in the male rats, independently of exposure to BFRs, were approximately five times higher than those in the pregnant female rats. This observation is consistent with previous studies showing a significant reduction of thyroid hormone levels in pregnant rats near term (Calvo et al., 1990). A number of in vivo rodent studies have reported a decrease in serum T4 concentrations after exposure to BFRs, with results varying depending on the specific BFRs, dose, length of exposure and gender or strain of animals studied (reviewed by Legler, 2008). With respect to the effects of BFR exposures during pregnancy on thyroid hormone status, the exposure of Sprague-Dawley rats to DE-71 (60 or 120 mg/kg/day) from gestation day 6 to 19 lowered plasma T4 concentrations on gestation day 20 (Ellis-Hutchings et al., 2009), while exposure to DE-71, at a dose of 60 |g/kg/day, from gestation day 1.5 to lactation day 20 had no effect on serum T4 or T3 in the dams (Blake et al., 2011). Positive associations between TSH and the serum concentrations of lower-brominated PBDE congeners and OH-PBDEs were reported in second trimester pregnant women from California; however, the associations with free and total T4 were generally weak, inconsistent, and not statistically significant in this study (Zota et al., 2011). The lack of an effect of BFR exposure on serum T4 concentrations or thyroid gland histomorphology in this study suggests that the effects we observed are not ascribed to this mechanism.

The majority of the anomalies observed in BFR-exposed fetuses consisted of malpositioned or fused digits in the hind paws. To the best of our knowledge, no previous investigation has observed a

AIIIHI.E IN PRESS

R.G. Berger et al. / Toxicology xxx (2014) xxx-xxx 7

Table 6

Effects of BFR exposure on fetal developmental measures on gestation day 20.

Parameter Sex BFR mixture dose (mg/kg/day)

0 0.06 20 60

Total number of live pups (litters) 255(15) 365(21) 285(17) 283(17)

Placental weight (mg) 477.0 ±10.5 470.4 ±9.57 464.6 ± 14.2 494.4 ±21.7

Sex ratio (M/F) 1.07 ±0.14 1.22 ±0.11 1.12 ±0.17 1.18 ±0.12

Body weight (g) M 3.64 ±0.07 3.50 ±0.05 3.45 ± 0.08 3.57 ± 0.05

F 3.39 ±0.07 3.32 ±0.05 3.20 ±0.07 3.35 ± 0.05

Crown rump length M 36.7 ± 0.2 36.3 ±0.2 36.0 ± 0.3 36.5 ± 0.2

(mm) F 35.7 ± 0.2 35.2 ±0.1 34.8 ±0.2' 35.4 ±0.2

Anogenital M 139.9 ±1.5 140.0 ±0.9 139.1 ±1.1 136.6 ±1.4

distance index3 F 62.1 ± 0.9 60.4 ±0.6 61.1 ± 0.7 60.6 ± 0.8

Values are expressed as mean ± SEM using the litter as the data point. a Anogenital distance index = (anogenital distance (mm)/cubed body weight (g)) x 100. * p <0.05 compared with control.

similar phenotype. One previous investigation reported a significant increase in tail malformations in the F2 generation following a single administration of BDE-99 on GD 6 (Talsness et al., 2005). The digit defects we observed were not associated with changes in ossi-

fication or cartilage development in the paws of exposed fetuses, implying that only soft tissues are affected. The proportion of litters with fetuses with malpositioned or fused digits was significantly increased in the low (0.06 mg/kg/day) dose treatment groups. Digit

Fig. 1. Effects of gestational BFR exposure on fetal external morphology. Representative pictures of fetal hindpaws on gestation day 20. (A) Normal paw (from the 60 mg BFR/kg/body weight/day treatment group). (B) Paw with fused digits (60 mg BFR/kg/body weight/day). (C) Paw with malpositioned digits (0.06 mg BFR/kg body weight/day). (D) Scatterplot depicting the numbers of fetuses with malpositioned or fused digits in individual litters; each data point represents a single litter (n =16, 21,17 and 18) for the control (0), 0.06, 20 and 60 mg/kg/day BFR groups, respectively. The proportion of litters with affected fetuses was significantly increased only in the 0.06 mg/kg/day BFR treatment group (0.06 mg/kg/day, p = 0.05; 20 mg/kg/day: p = 0.09; 60 mg/kg/day: p = 0.21).

fill"И^^^И AIIIIILE IN PRESS

R.G. Berger et al. / Toxicology xxx (2014) xxx-xxx

Fig. 2. Effects of gestational BFR exposure on fetal skeletal development. Representative pictures of fetal skeletal ossification on gestation day 20 with all six sternebrae ossified (A, 0.06 mg BFR/kg body weight), two sternebrae sites unossified (B, 60 mg BFR/kg body weight), offset ossification of the sternebrae 1-4 and incomplete ossification Q7 of sternebrae 5 and 6 (C, 0.06 mg BFR/kg body weight), and an incompletely fused ossification vertebral centrum (D, 0.06 mg BFR/kg body weight). White arrows show sites of typical ossification, yellow circles mark areas with decreased or no ossification, and red arrows show markedly offset sternebrae (C) or an incompletely fused ossification vertebral centrum (D). (E) Proportion of fetuses with any type of skeletal variation, including lack of ossification of the sternebrae, markedly offset fusion of sternebrae and incomplete or lack of ossification of the fused vertebral centrum. (F) Proportion of fetuses with no ossification of the sixth sternebra. Bars represent the proportion of affected fetuses relative to the total number of fetuses examined (n = 18, 41, 26 and 38 for the control (0), 0.06, 20 and 60 doses, respectively). *p <0.05, **p< 0.01 compared with control. (For interpretation of the references to color in this figure legend, the reader is referred to the web version of the article.)

Table 7

Effects of BFR exposure on fetal skeletal development.

Variation BFR mixture dose (mg/kg/day)

0 0.06 20 60

Number of fetuses examined 18 41 26 38

Two or more unossified sternebrae 1(5.6) 11 (26.8)' 8 (30.8)' 12(31.6)'

Incomplete vertebral centrum fusion 2(11.1) 8(19.5) 6 (23.1) 8(21.1)

AH spinal skeletal variations 2(11.1) 9(22.0) 6(23.1) 11(29.0)

Values are the number of observed variations per group with the percentage of the total number of examined fetuses in brackets. * p <0.05 compared with control.

defects were also observed in the 20 mg/kg/day and 60 mg/kg/day treatment groups but statistical significance was not attained.

An increase in the incidence of fetuses with skeletal variations was observed following exposure to the BFR mixture. The most

common variation was a decrease in ossification in the sternebrae, particularly the sixth sternebra. Ossification in the fetus is sensitive to maternal toxicity (Beck, 1989); however, no indications of toxicity were observed in dams exposed to either the low or

ARTICLE IN PRESS

R.G. Berger et al. / Toxicology xxx (2014) xxx-xxx

472 middle dose BFR regimens. In utero exposures to exogenous

473 chemicals such as ethanol or tobacco smoke alter fetal skeletal

474 ossification; however, these insults are most often associated with

475 endpoints that were not observed here, such as decreased fetal

476 survival and growth retardation (Chernoff, 1977; Seller and Bnait,

477 1995).

478 Ossification of the rat skeleton typically commences around GDs

479 15-16 and is considered an indicator of maturity. At birth typ-

480 ically all (95-100%) of sternebrae are ossified. Variations in the

481 ossification of the sternum of rodents are relatively common, with

482 incomplete ossification one day before birth ranging from 1-7%

483 for sternebrae 1-4 and 6, and up to 35% for sternebra 5 (Fritz

484 and Hess, 1970). In the current investigation, the sixth sterne-

485 bra was unossified in 35-40% of GD20 fetuses exposed to any

486 BFR dose, compared to 5.5% in controls. These data suggest that

487 exposure to BFRs may induce a delay in ossification. Delays in

488 ossification may also be induced by fetal hypothyroidism (Capelo

489 et al., 2008). Previously, single BFR congeners have been reported

490 to induce an increase in the incidence of variations in ossifica-

491 tion following in utero exposure. The administration of PBDE-99

492 (2mg/kg/day) to Sprague-Dawley rats from GDs 6-19 increased

493 skeletal defects, including incomplete ossification in the skull, cau-

494 dal vertebrae and ribs, but did not affect the ossification of the

495 sternebrae (Blanco et al., 2012). An earlier investigation in rabbits

496 reported the delayed ossification of sternebrae following exposure

497 to polybromodiphenyl oxide at a high dose (15 mg/kg/day) on GDs

498 7-20 (Breslin et al., 1989). This report, to the best of our knowledge,

499 is the first to show a significant impact on fetal bone ossification

500 following exposure to an environmentally relevant BFR mixture.

501 Since reduced ossification may represent a delay in developmen-

502 tal timing it would be useful to examine the skeletons of postnatal

503 pups to determine whether these effects are persistent (Carney and

504 Kimmel, 2007).

505 The current investigation shows that BFR exposure induced

506 developmental alterations at environmentally relevant doses that

507 have not been reported previously. This is also the first investigation

508 to determine the effects of exposure to an environmentally relevant

509 complex mixture of BFRs in an animal model. Previous experi-

510 ments have shown that mixtures of endocrine disrupting chemicals

511 may produce effects that are not observed when the components

512 are administered individually (reviewed by Kortenkamp, 2007).

513 For instance, exposure to a mixture of thyroid hormone disrupt-

514 ing chemicals, including dioxins and polychlorinated biphenyls,

515 produced effects suggesting synergistic responses and the acti-

516 vation of different mechanisms than those observed when the

517 individual components were administered (Crofton et al., 2005).

518 Since there are low concentrations of polybrominated dibenzo-

519 p-dioxins (PBDDs) and dibenzofurans (PBDFs) in the commercial

520 BFR mixtures (Hanari et al., 2006), these compounds may also

521 contribute to the effects observed. Measurements of BFRs in

522 serum from pregnant women (Buttke et al., 2013; Woodruff et al.,

523 2011), umbilical cord blood (Frederiksen et al., 2010), and fetal

524 liver and placental tissues (Doucet et al., 2009) demonstrate that

525 humans are exposed to multiple BFR congeners and technical mix-

526 tures throughout gestation. Therefore, elucidation of the effects

527 and potential dangers of exposure to complex, environmentally

528 relevant, mixtures is critical since the impact of exposure to

529 environmentally relevant BFR mixtures may differ from that of

530 individual congeners.

531 In conclusion, we have shown evidence that exposure to the

532 mixture of BFRs commonly found in North American household

533 dust alters fetal development. The phenotypes observed, includ-

534 ing malpositioned or fused digits and delayed ossification in the

535 sternebrae, have not been reported previously in investigations

536 of the effects of individual BFR congeners or technical mixtures.

537 The absence of concurrent effects on measures of toxicity and

pregnancy maintenance in the dams suggests that exposure to a 538

wide range of BFR doses may have a direct impact on fetal devel- 539

opment. 540

Supplementary data 541

Supplementary data Table S1 presents the absolute values for 542

dam organ weights. Scatterplots of the incidence of skeletal vari- 543

ations are provided in Supplementary data Fig. S1. Images of the 544

skeletal staining of paws from control and BFR exposed fetuses are 545

shown in Supplementary data Fig. S2. 546

Uncited reference Q4 547

Costa et al. (2008). 548

Conflict of interest 549

The authors declare that there are no conflicts of interest. Q5 550

Transparency document 551

The Transparency document associated with this article can be 552

found in the online version. 553

Acknowledgements 554

The authors acknowledge Mr. Donald Demers, Animal Resources 555

Division Health Canada, for food pelleting and Dr. Chunwei Huang 556

and Ms. Lydia Goff, McGill University, for their assistance. This 557 research was supported by grant RHF100625 from the Canadian Q6 558

Institutes of Health Research (CIHR) Institute for Human Develop- 559

ment, Child and Youth Health. RGB is the recipient of an NSERC 560

Create award from the Réseau Québécois en Reproduction and PLCL 56i

of a Fonds de la Recherche Santé du Québec fellowship. BR and BFH 562

are James McGill Professors. 563

References 564

Aleksa, K., Carnevale, A., Goodyer, C., Koren, G., 2012. Detection of polybrominated 565

biphenyl ethers (PBDEs) in pediatric hair as a tool for determining in utero 566

exposure. Forensic Sci. Int. 218,37-43. 567

Allen, J.G., McClean, M.D., Stapleton, H.M., Webster, T.F., 2008. Critical factors in 568

assessing exposure to PBDEs via house dust. Environ. Int. 34,1085-1091. 569

Beck, S.L., 1989. Prenatal ossification as an indicator of exposure to toxic agents. 570

Teratology 40,365-374. 571

Blake, C.A., McCoy, G.L., Hui, Y.Y., LaVoie, H.A., 2011. Perinatal exposure to low-dose 572

DE-71 increases serum thyroid hormones and gonadal osteopontin gene expres- 573

sion. Exp. Biol. Med. (Maywood) 236 (4), 445-455. 574

Blanco, J., Mulero, M., Domingo, J.L., Sanchez, D.J., 2012. Gestational exposure to 575

BDE-99 produces toxicity through upregulation of CYP isoforms and ROS pro- 576

duction in the fetal rat liver. Toxicol. Sci. 127, 296-302. 577

Breslin, W.J., Kirk, H.D., Zimmer, M.A., 1989. Teratogenic evaluation of a poly- 578

bromodiphenyl oxide mixture in New Zealand white rabbits following oral 579

exposure. Fundam. Appl. Toxicol. 12,151-157. 580

Buttke, D.E., Wolkin, A., Stapleton, H.M., Miranda, M.L., 2013. Associations between 58i

serum levels of polybrominated diphenyl ether (PBDE) flame retardants and 582

environmental and behavioral factors in pregnant women. J. Expo. Sci. Environ. 583

Epidemiol. 23,176-182. 584

Calvo, R., Obregon, M.J., De Ona, C.R., Ferreiro, B., Escobar Del Rey, F., Morreale De 585

Escobar, G., 1990. Thyroid hormone economy in pregnant rats near term: a 586

physiological animal model of nonthyroidal illness? Endocrinology 127,10-16. 587

Camino, G., Costa, L., Luda di Cortemiglia, M.P., 1991. Overview of fire retardant 588

mehcanisms. Polym. Degrad. Stab. 98,1662-1670. 589

Capelo, L.P., Beber, E.H., Huang, S.A., Zorn, T.M., Bianco, A.C., Gouveia, C.H., 2008. 590

Deiodinase-mediated thyroid hormone inactivation minimizes thyroid hor- 591

mone signaling in the early development of fetal skeleton. Bone 43,921-930. 592

Carney, E.W., Kimmel, C.A., 2007. Interpretation of skeletal variations for human 593

risk assessment: delayed ossification and wavy ribs. Birth Defects Res. B: Dev. 594

Reprod. Toxicol. 80, 473-496. 595 Carrizo, D., Grimalt, J.O., Ribas-Fito, N., Sunyer, J., Torrent, M., 2007. Influence of 596

breastfeeding in the accumulation of polybromodiphenyl ethers during the first 597

years of child growth. Environ. Sci. Technol. 41, 4907-4912. 598

li'lB^^^M AIIIIII.E IN PRESS

R.G. Berger et al. / Toxicology xxx (2014) xxx-xxx

Ceccatelli, R., Faass, O., Schlumpf, M., Lichtensteiger, W., 2006. Gene expression and estrogen sensitivity in rat uterus after developmental exposure to the poly-brominated diphenylether PBDE 99 and PCB. Toxicology 220,104-116.

Chernoff, G.F., 1977. The fetal alcohol syndrome in mice: an animal model. Teratology 15, 223-229.

Costa, L.G., Giordano, G., Tagliaferri, S., Caglieri, A., Mutti, A., 2008. Polybrominated diphenyl ether (PBDE) flame retardants: environmental contamination, human body burden and potential adverse health effects. Acta Biomed. 79,172-183.

Crofton, K.M., Craft, E.S., Hedge, J.M., Gennings, C., Simmons, J.E., Carchman, R.A., CarterJr., W.H., DeVito, M.J., 2005. Thyroid-hormone-disrupting chemicals: evidence for dose-dependent additivity or synergism. Environ. Health Perspect. 113,1549-1554.

Darnerud, P.O., Eriksen, G.S., Johannesson, T., Larsen, P.B., Viluksela, M., 2001. Polybrominated diphenyl ethers: occurrence, dietary exposure, and toxicology. Environ. Health Perspect. 109 (Suppl. 1), 49-68.

Dodson, R.E., Perovich, L.J., Covaci, A., Van den Eede, N., Ionas, A.C., Dirtu Brody, J.G., Rudel, R.A., 2012. After the PBDE phase-out: a broad suite of flame retardants in repeat house dust samples from California. Environ. Sci. Technol. 46, 13056-13066.

Doucet, J., Tague, B., Arnold, D.L., Cooke, G.M., Hayward, S., Goodyer, C.G., 2009. Persistent organic pollutant residues in human fetal liver and placenta from Greater Montreal, Quebec: a longitudinal study from 1998 through 2006. Environ. Health Perspect. 117,605-610.

Ellis-Hutchings, R.G., Cherr, G.N., Hanna, L.A., Keen, C.L., 2009. The effects of marginal maternal vitamin A status on penta-brominated diphenyl ether mixture-induced alterations in maternal and conceptal vitamin A and fetal development in the Sprague Dawley rat. Birth Defects Res. B: Dev. Reprod. Toxicol. 86, 48-57.

Environment Canada. 2011. Available: www.ec.gc.ca/ceparegistry/the_act/schedules.

Ema, M., Fujii, S., Hirata-Koizumi, M., Matsumoto, M., 2008. Two-generation reproductive toxicity study of the flame retardant hexabromocyclododecane in rats. Reprod. Toxicol. 25,335-351.

Ernest, S.R., Wade, M.G., Lalancette, C., Ma, Y.Q., Berger, R.G., Robaire, B., Hales, B.F., 2012. Effects of chronic exposure to an environmentally relevant mixture of brominated flame retardants on the reproductive and thyroid system in adult male rats. Toxicol. Sci. 127,496-507.

Eskenazi, B., Chevrier, J., Rauch, S.A., Kogut, K., Harley, K.G., Johnson, C., Trujillo, C., Sjodin, A., Bradman, A., 2013. In utero and childhood polybrominated diphenyl ether (PBDE) exposures and neurodevelopment in the CHAMACOS study. Environ. Health Perspect. 121 (2), 257-262.

Fowles, J.R., Fairbrother, A., Baecher-Steppan, L., Kerkvliet, N.I., 1994. Immunologic and endocrine effects ofthe flame-retardant pentabromodiphenyl ether (DE-71) in C57BL/6J mice. Toxicology 86,49-61.

Fraser, A.J., Webster, T.F., McClean, M.D., 2009. Diet contributes significantly to the body burden of PBDEs in the general US population. Environ. Health Perspect. 117,1520-1525.

Frederiksen, M., Thomsen, C., Frashaug, M., Vorkamp, K., Thomsen, M., Becher, G., Knudsen, L.E., 2010. Polybrominated diphenyl ethers in paired samples of maternal and umbilical cord blood plasmaand associations with house dust ina Danish cohort. Int. J. Hyg. Environ. Health 213, 233-242.

Fritz, H., Hess, R., 1970. Ossification of the rat and mouse skeleton in the perinatal period. Teratology 3,331-337.

GallavanJr., R.H., Holson,J.F., Stump, D.G., Knapp, J.F., Reynolds, V.L., 1999. Interpreting the toxicologic significance of alterations in anogenital distance: potential for confounding effects of progeny body weights. Reprod. Toxicol. 13,383-390.

Goldman, J.M., Murr, A.S., Cooper, R.L., 2007. The rodent estrous cycle: characterization of vaginal cytology and its utility in toxicological studies. Birth Defects Res. B: Dev. Reprod. Toxicol. 80, 84-97.

Haave, M., Folven, K.I., Carrol, T., Glover, C., Heegaard, E., Brattelid, T., Hogstrand, C., Lundebye, A.K., 2011. Cerebral gene expression and neurobehavioural development after perinatal exposure to an environmentally relevant polybrominated diphenylether (BDE47). Cell Biol. Toxicol. 27,343-361.

Hallgren, S., Sinjari, T., Hakansson, H., Darnerud, P.O., 2001. Effects of polybrominated diphenyl ethers (PBDEs) and polychlorinated biphenyls (PCBs) on thyroid hormone and vitamin a levels in rats and mice. Arch. Toxicol. 75, 200-208.

Hanari, N., Kannan, K., Miyake, Y., Okazawa, T., Kodavanti, P.R.S., Aldous, K.M., Yamashita, N., 2006. Occurrence of polybrominated dibenzo-rho-dioxins, and polybrominated dibenzofurans as impurities in commercial polybrominated biphenyl ether mixtures. Environ. Sci. Technol. 40,4400-4405.

Hites, R.A., 2004. Polybrominated diphenyl ethers in the environment and in people: a meta-analysis of concentrations. Environ. Sci. Technol. 38, 945-956.

Johnson-Restrepo, B., Kannan, K., Rapaport, D.P., Rodan, B.D., 2005. Polybrominated diphenyl ethers and polychlorinated biphenyls in human adipose tissue from New York. Environ. Sci. Technol. 39, 5177-5182.

Johnson, P.I., Stapleton, H.M., Sjodin, A., Meeker, J.D., 2010. Relationships between polybrominated diphenyl ether concentrations in house dust and serum. Environ. Sci. Technol. 44, 5627-5632.

Johnson, P.I., Altshul, L., Cramer, D.W., Missmer, S.A., Hauser, R., Meeker, J.D., 2012. Serum and follicular fluid concentrations of polybrominated diphenyl ethers and in-vitro fertilization outcome. Environ. Int. 45,9-14.

Kortenkamp, A., 2007. Ten years of mixing cocktails: a review of combination effects of endocrine-disrupting chemicals. Environ. Health Perspect. 115, 98-105.

Legler, J., 2008. New insights into the endocrine disrupting effects of brominated flame retardants. Chemosphere 73, 216-222.

Lilienthal, H., Hack, A., Roth-Harer, A., Grande, S.W., Talsness, C.E., 2006. Effects of developmental exposure to 2,2',4,4',5-pentabromodiphenyl ether (PBDE-99)

on sex steroids, sexual development, and sexually dimorphic behavior in rats. Environ. Health Perspect. 114,194-201.

Marchitti, S.A., LaKind, J.S., Naiman, D.Q., Berlin, C.M., Kenneke, F.J., 2013. Improving infant exposure and health risk estimates: using serum data to predict polybrominated diphenyl ether concentrations in breast milk. Envrion. Sci. Technol. 47, 4787-4795.

Mazdai, A., Dodder, N.G., Abernathy, M.P., Hites, R.A., Bigsby, R.M., 2003. Polybromi-nated diphenyl ethers in maternal and fetal blood samples. Environ. Health Perspect. 111,1249-1252.

Meeker, J.D., Johnson, P.I., Camann, D., Hauser, R., 2009. Polybrominated diphenyl ether(PBDE) concentrations in house dust are related to hormone levels in men. Sci. Total Environ. 407,3425-3429.

Meerts, I., Letcher, R.J., Hoving, S., Marsh, G., Bergman, A., Lemmen,J.G., van derBurg, B., Brouwer, A., 2001. In vitro estrogenicity of polybrominated diphenyl ethers, hydroxylated PDBES, and polybrominated bisphenol A compounds. Environ. Health Perspect. 109,399-407.

Petro, E.M., Leroy, J.L., Covaci, A., Fransen, E., De Neubourg, D., Dirtu, A.C., De Pauw, I., Bols, P.E., 2012. Endocrine-disrupting chemicals in human follicular fluid impair in vitro oocyte developmental competence. Hum. Reprod. 27, 1025-1033.

Rose, M., Bennett, D.H., Bergamn, A., Fangstrom, B., Pessah, I.N., Hertz-Picciotto, I., 2010. PBDEs in 2-5 year old children from California and associations with diet and indoor environment. Environ. Sci. Technol. 44, 2648-2653.

Seller, M.J., Bnait, K.S., 1995. Effects of tobacco smoke inhalation on the developing mouse embryo and fetus. Reprod. Toxicol. 9, 449-459.

Siddique, S., Xian, Q., Abdelouahab, N., Takser, L., Phillips, S.P., Feng, Y.L., Wang, B., Zhu, J., 2012. Levels of dechlorane plus and polybrominated diphenylethers in human milk in two Canadian cities. Environ. Int. 39, 50-55.

1.cfnSjodin, A., Jones, R.S., Focant, J.F., Lapeza, C., Wang, R.Y., McGahee, E.E., 2004. Retrospective time-trend study of polybrominated diphenyl ether and poly-brominated and polychlorinated biphenyl levels in human serum from the United States. Environ. Health Perspect. 112,654-658.

Stapleton, H.M., Dodder, N.G., Offenberg, J.H., Schantz, M.M., Wise, S.A., 2005. Poly-brominated diphenyl ethers in house dust and clothes dryer lint. Environ. Sci. Technol. 39, 925-931.

Stapleton, H.M., Allen, J.G., Kelly, S.M., Konstantinov, A., Klosterhaus, S., Watkins, D., McClean, M.D., Webster, T.F., 2008. Alternate and new brominated flame retardants detected in U.S. house dust. Environ. Sci. Technol. 42,6910-6916.

Stapleton, H.M., Eagle, S., Sjodin, A., Webster, T.F., 2012. Serum PBDEs in a North Carolina toddler cohort: associations with handwipes, house dust, and socioeconomic variables. Environ. Health Perspect. 120,1049-1054.

Stoker, T.E., Laws, S.C., Crofton, K.M., Hedge, J.M., Ferrell, J.M., Cooper, R.L., 2004. Assessment of DE-71, a commercial polybrominated diphenyl ether (PBDE) mixture, in the EDSP male and female pubertal protocols. Toxicol. Sci. 78, 144-155.

Stoker, T.E., Cooper, R.L., Lambright, C.S., Wilson, V.S., Furr,J., Gray, L.E., 2005. In vivo and in vitro anti-androgenic effects of DE-71, a commercial polybrominated diphenyl ether (PBDE) mixture. Toxicol. Appl. Pharmacol. 207, 78-88.

Suvorov, A., Battista, M.C., Takser, L., 2009. Perinatal exposure to low-dose 2,2',4,4'-tetrabromodiphenyl ether affects growth in rat offspring: what is the role of IGF-1? Toxicology 260,126-131.

Suvorov, A., Takser, L., 2010. Global gene expression analysis in the livers of rat offspring perinatally exposed to low doses of 2,2',4,4'-tetrabomodiphenyl ether. Environ. Health Perspect. 118,97-102.

Talsness, C.E., Shakibaei, M., Kuriyama, S.N., Grande, S.W., Sterner-Kock, A., Schnitker, P., de Souza, C., Grote, K., Chahoud, I., 2005. Ultrastructural changes observed in rat ovaries following in utero and lactational exposure to low doses of a polybrominated flame retardant. Toxicol. Lett. 157,189-202.

Toms, L.-M.L., Hearn, L., Kennedy, K., Harden, F., Bartkow, M., Temme, C., Jochen, F., Mueller, J.F., 2009. Concentrations ofpolybrominated diphenyl ethers (PBDEs) in matched samples of human milk, dust and indoor air. Environ. Int. 35, 864-869.

Tullo, A., 2003. Great lakes to phase out two flame retardants. Chem. Eng. News 81, 13.

van der Ven, L.T., Verhoef, A., van de Kuil, T., Slob, W., Leonards, P.E., Visser, T.J., Hamers, T., Herlin, M., Hakansson, H., Olausson, H., Piersma, A.H., Vos, J.G., 2006. A 28-day oral dose toxicity study enhanced to detect endocrine effects of hexabromocyclododecane in Wistar rats. Toxicol. Sci. 94, 281-292.

van der Ven, L.T., van de Kuil, T., Leonards, P.E., Slob, W., Canton, R.F., Germer, S., Visser, T.J., Litens, S., Hakansson, H., Schrenk, D., van den Berg, M., Piersma, A.H., Vos, J.G., Opperhuizen, A., 2008. A 28-day oral dose toxicity study in Wistar rats enhanced to detect endocrine effects of decabromodiphenyl ether (decaBDE). Toxicol. Lett. 179, 6-14.

van der Ven, L.T., van de Kuil, T., Leonards, P.E., Slob, W., Lilienthal, H., Litens, S., Herlin, M., Hakansson, H., Canton, R.F., van den Berg, M., Visser, T.J., van Loveren, H., Vos, J.G., Piersma, A.H., 2009. Endocrine effects of hexabromocyclododecane (HBCD) in a one-generation reproduction study in Wistar rats. Toxicol. Lett. 185, 51-62.

Whitaker, J., Dix, K.M., 1979. Double staining technique for rat foetus skeletons in teratological studies. Lab. Anim. 13,309-310.

Woodruff, T.J., Zota, A.R., Schwartz, J.M., 2011. Environmental chemicals in pregnant women in the United States: NHANES 2003-2004. Environ. Health Perspect. 119, 878-885.

Wu, N., Herrmann, T., Paepke, O., Tickner, J., Hale, R., Harvey, L.E., La Guardia, M., McClean, M.D., Webster, T.F., 2007. Human exposure to PBDEs: associations of PBDE body burdens with food consumption and house dust concentrations. Environ. Sci. Technol. 41,1584-1589.

I'lllB^^^M AIIIHI.E IN PRESS

R.G. Berger et al. / Toxicology xxx (2014) xxx-xxx 11

771 Zhou, T., Taylor, M.M., DeVito, M.J., Crofton, K.M., 2002. Developmental exposure to San Francisco General Hospital, California. Environ. Sci. Technol. (Epub ahead 776

772 brominated diphenyl ethers results in thyroid hormone disruption. Toxicol. Sci. of print). 777

773 66,105-116. Zota, A.R., Park, J.S., Wang, Y., Petreas, M., Zoeller, R.T., Woodruff, T.J., 2011. Poly- 778

774 Zota, A.R., Linderholm, L., Park, J.S., Petreas, M., Guo, T., Privalsky, M.L., Zoeller, brominated diphenyl ethers, hydroxylated polybrominated diphenyl ethers, and 779

775 R.T., Woodruff, T.J., 2013. Temporal comparison of PBDEs, OH-PBDEs, PCBs, and measures of thyroid function in second trimester pregnant women in California. 780 OH-PCBs in the serum of second trimester pregnant women recruited from Environ. Sci. Technol. 45 (18), 7896-7905. 781