Effects of chronic exposure to an environmentally relevant mixture of brominated flame retardants on the reproductive and thyroid system in adult male rats.

Brominated flame retardants (BFRs) are incorporated into a wide variety of consumer products, are readily released into home and work environments, and are present in house dust. Studies using animal models have revealed that exposure to polybrominated diphenyl ethers (PBDEs) may impair adult male reproductive function and thyroid hormone physiology. Such studies have generally characterized the outcome of acute or chronic exposure to a single BFR technical mixture or congener but not the impact of environmentally relevant BFR mixtures. We tested whether exposure to the BFRs found in house dust would have an adverse impact on the adult male rat reproductive system and thyroid function. Adult male Sprague Dawley rats were exposed to a complex BFR mixture composed of three commercial brominated diphenyl ethers (52.1% DE-71, 0.4% DE-79, and 44.2% decaBDE-209) and hexabromocyclododecane (3.3%), formulated to mimic the relative congener levels in house dust. BFRs were delivered in the diet at target doses of 0, 0.02, 0.2, 2, or 20 mg/kg/day for 70 days. Compared with controls, males exposed to the highest dose of BFRs displayed a significant increase in the weights of the kidneys and liver, which was accompanied by induction of CYP1A and CYP2B P450 hepatic drug-metabolizing enzymes. BFR exposure did not affect reproductive organ weights, serum testosterone levels, testicular function, or sperm DNA integrity. The highest dose caused thyroid toxicity as indicated by decreased serum thyroxine (T4) and hypertrophy of the thyroid gland epithelium. At lower doses, the thickness of the thyroid gland epithelium was reduced, but no changes in hormone levels (T4 and thyroid-stimulating hormone) were observed. Thus, exposure to BFRs affected liver and thyroid physiology but not male reproductive parameters.

Brominated flame retardants (BFRs) are frequently added to electronics, home furnishings, and a variety of other consumer items to reduce flammability. Polybrominated diphenyl ethers (PBDEs) and hexabromocyclododecane (HBCD) are among the most widely used additive flame retardants (Alaee ). There are three general commercial formulations of the PBDEs, penta, octa, and decabromodiphenyl ethers (Alaee ). As a result of evidence of persistence and bioaccumulation, the pentabromo and octabromo forms were banned in Europe in 2004, in California in 2008, and U.S. producers agreed to cease their production by 2004 (Tullo, 2003). In 2008, the Canadian government prohibited their import, production, or sale and in 2010, recommended that similar measures be applied to decabromodiphenyl ether (decaBDE) (Environment Canada, 2010b). The need to meet stringent flammability standards in many jurisdictions has resulted in producers replacing PBDEs with alternates such as HBCD. However, a recent assessment of HBCD by the Canadian government has proposed that HBCD be banned based on its persistence in the environment and ability to biomagnify through food chains (Environment Canada, 2010a).

As neither PBDEs nor HBCD are covalently linked with the polymers of the matrices to which they are added, they are released readily, leading to contamination of both home and work environments. It has been suggested that the largest source of human exposure to PBDEs is from ingestion of contaminated house dust (Jones-Otazo ), which is estimated to account for as much as 80–93% of daily exposure in toddlers (Wilford ). Studies of house dust in North America have revealed that the PBDEs and HBCD were the most prominent BFRs detected in the majority of houses sampled (Allen ; Stapleton ).

The potential cumulative effects of exposure to PBDEs and HBCD on human health are unknown although the study of the effects of PBDEs, as either technical mixtures or individual congeners, remains an active field of investigation. Several lines of evidence suggest that exposure to PBDEs may impair adult male reproductive function and thyroid hormone physiology. Exposure to house dust concentrations of penta-bromodiphenyl ether (pentaBDE) congeners was negatively correlated with luteinizing and follicle-stimulating hormone concentrations (Meeker ) and sperm numbers and testis size (Akutsu ). Although decaBDE has been detected in semen samples of Chinese workers who process electronic waste (Liu ), studies have not yet linked its exposure to effects on the male reproductive system in humans. In adult male rats, exposure to pentaBDE (van der Ven ) or decaBDE (Van der Ven ) reduced the weights of the epididymides and altered seminal vesicle weights.

There is evidence that PBDEs may interfere with thyroid hormone signaling in addition to affecting male reproductive parameters. PentaBDE levels in adult men exposed to house dust (Meeker ) and in Great Lakes sport fish consumers (Turyk ) were positively correlated with free thyroxine (T4); their associations with triiodothyronine (T3) or thyroid-stimulating hormone (TSH) were variable. In adult male rats, pentaBDE exposure decreased serum T4 levels without affecting T3 (van der Ven ).

Less information is available on the thyroid or male reproductive effects of HBCD exposure. The sole study to examine the correlation between HBCD exposure and measures of thyroid function in humans studied mother-infant pairs in Norway (Eggesbo ) and failed to observe any significant relationship between breast milk HBCD levels and neonatal TSH. In animal studies, HBCD did not significantly affect reproductive or thyroid parameters (van der Ven ).

Because humans are constantly exposed to complex mixtures of BFRs, and epidemiological studies have suggested that exposures to BFRs at environmentally relevant levels correlate with disruption in male reproductive and thyroid functions, our goal was to determine the effects of exposure to the BFRs found in house dust in an animal model. We reconstituted a complex BFR mixture based on the relative levels reported in a representative sample of house dust in North America (Allen ; Stapleton ) and evaluated the effects of chronic exposure across a broad range of exposure levels on adult male rats. Specific emphasis was placed on measures of reproductive- and thyroid-related effects.

MATERIALS AND METHODS

Mixture Formulation

The BFR mixture was composed by combining three technical commercial mixtures of PBDEs (DE-71, DE-79, and BDE-209) and one mixture of HBCD to yield a ratio of PBDE congeners and HBCD comparable to that observed in house dust in Boston, MA (Allen ; Stapleton ) (Table 1). The HBCD was composed of all three isomers (alpha, beta, and gamma) but was enriched for the more environmentally stable alpha congener to represent the composition that might be found in house dust (Roosens ). Technical DE-79 was purchased from Wellington Laboratories (Guelph, Canada), and BDE-209 (> 97% pure) was purchased from Sigma-Aldrich (St Louis, MO). DE-71 was generously provided by Doug Arnold (Health Canada), from stock provided by Chemtura (Lawrenceville, GA). The HBCD mixture was kindly provided by Ivan Curran (Health Canada). The BFR formulation was dissolved in corn oil and mixed into a powdered isoflavone-free rodent diet (Teklad Global 2019 diet; Harlan Laboratories, Madison, WI) with 15% vol/wt deionized water, pelleted, and allowed to dry at room temperature. Corn oil added to the diets amounted to 2.4 g/kg in addition to the 90 g/kg of soya oil present in this diet. Dried pellets were stored away from light at 4°C for no longer than 1 month prior to feeding. Diets contained 0, 0.25, 2.5, 25, and 250 mg of the BFR mixture/kg diet; assuming a daily food consumption of 80 g/kg body weight/day, this formulation delivered nominal doses of 0, 0.02, 0.2, 2.0, and 20 mg/kg body weight/day weight. The calculation of the lowest dose was based on assumptions of dust ingestions rates (100 mg/day) of children (16.5 kg body weight) and the scaling of dose from humans to rodents based on body surface area (1:6.3 human to rat) rather than weight. This dose estimates the highest levels of human exposure as indicated by the maximum BFR levels measured in house dust in the Boston studies (Allen ; Stapleton ).

TABLE 1

Diet Formulation

	House dust levelsa		Congeners in component technical mixtures (%)
Congener	Median (ng/g)	Percentage (%)	DE-71b (52.08%)e	DE-79c (0.36%)e	BDE-209d (44.18%)e	HBCD (3.38%)e	Congener levels in final mixture (%)
BDE-17	1.4	0.01	0.04				0.02
BDE-28	16.3	0.16	0.37				0.19
BDE-47	1865	17.93	33.93				17.19
BDE-49	29.6	0.28	0.77				0.40
BDE-66	17.2	0.17	1.02				0.53
BDE-85	124	1.19	3.18				1.66
BDE-99	2460	23.64	41.40	0.03			21.56
BDE-100	436.3	4.19	7.10				3.70
BDE-138	20.9	0.2	0.24	0.50			0.13
BDE-153	234.4	2.25	3.75	6.17			4.18
BDE-154	182.8	1.76	3.00	0.93			1.57
BDE-183	27.9	0.27	0.08	37.20			0.18
BDE-197	2.7	0.03		21.60			0.08
BDE-207	45.9	0.44		13.70			0.05
BDE-209	4502	43.27		1.81	98.00		43.30
HBCD		3.40				99.00	3.35

Values for PBDE congeners in dust are those reported by Allen , whereas total HBCD values are derived from the analyses of the same samples but reported by Stapleton

Congener levels in technical DE-71 are those reported by Konstantinov .

Congener levels in technical DE-79 are those reported in the certificate of analysis provided by supplier.

Congener levels in BDE-209 are those reported in the certificate of analysis provided by supplier.

Amount in final mixture.

Animals and Treatment

Adult Sprague Dawley male rats (350–375 g) were obtained from Charles River (St-Constant, Canada) and maintained on a 12-L:12-D cycle in the Animal Resources Centre of McGill University. Food and water were provided ad libitum. 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).

Animals were housed individually, and BFR mixture −72 consumption was monitored by measuring food weight on a weekly basis. Animals were acclimatized on the control diet for 1 week and randomly divided into the control and one of the BFR mixture four treatment groups (n = 15 per group) for a period of 70 days, a sufficient time for germ cells to progress from spermatogonia to mature spermatozoa. Animals were weighed and examined physically once a week.

Tissue Collection

At the end of the 70-day treatment period, 50 rats (n = 10 per group) were euthanized by CO2 asphyxiation followed by exsanguination via cardiac puncture. The following organs were dissected: liver, kidneys, spleen, heart, lung, thymus, epididymis, testes, ventral prostate, vas deferens, and empty seminal vesicles (with coagulating glands). All organs were examined, weighed, snap frozen in liquid nitrogen, and stored at −80°C. Mature spermatozoa from the right cauda epididymidis were collected, as previously described (Delbes ), and stored at −80°C until further analysis. Sperm from the contralateral cauda epididymidis were collected for direct motility analysis by computer-assisted sperm motility (CASA). The trachea, including the thyroid gland, was gently removed and fixed in neutral buffered formalin before being processed for morphological analysis. Whole blood was transferred to a Vacutainer SST (BD Biosciences Canada, Mississauga, ON), allowed to clot for 30 min at room temperature, then held on ice until centrifugation at 1300 × g for 20 min. Serum was aliquoted and stored at −80°C until assayed.

The remaining 25 rats (n = 5 per group) were anaesthetized by injection of a cocktail containing ketamine (50 mg/kg), xylazine (5 mg/kg), and acepromazine (1 mg/kg); the left testis was cleared with saline and perfused with Bouin’s fluid through the abdominal aorta. The tissue was then excised, postfixed for an additional 24 h in the same fixative, dehydrated, and embedded in paraffin (Bieber ).

Clinical Chemistry

Calcium, magnesium, phosphorus, blood urea nitrogen, uric acid, albumin, total protein, creatine kinase, creatinine, triglycerides, cholesterol, glucose, alkaline phosphatase (ALP), and alanine aminotransferase (ALT) were measured in the serum (n = 8–10 per group) with an ABX Pentra 400 clinical chemistry analyzer (Horiba ABX, Montpellier, France).

Liver RNA Extraction and Real-Time Quantitative Reverse Transcription-PCR

RNAs were extracted from liver samples (n = 9–10 per group) using TRIzol (Invitrogen, Burlington, Canada) following the manufacturer’s protocol. Briefly, frozen liver samples (50–100 mg) were homogenized in 1 ml of TRIzol using a steel ball in a Retch MM 400 mixer mill for 2 min at 25 Hz. After TRIzol extraction, the RNA was cleaned using an RNeasy Mini kit (Qiagen Inc., Mississauga, ON). RNA quantity and quality were assessed using a NanoDrop spectrophotometer (ThermoFisher Scientific, Ottawa, Canada).

RNA was reverse transcribed in 20 μl reaction mix containing 1 μg total RNA using the QuantiTech reverse transcription kit (Qiagen Inc.). Primer sequences for Cyp1a1 (cytochrome P450, family1, member 1A1) and putative genes responsive to thyroid hormone receptor (Me1, malic enzyme 1, NADP(+)-dependent, cytosolic; Dio1, deiodinase I; and Thrsp, thyroid hormone responsive) and estrogen receptor (Apoa4, apolipoprotein A-IV; Igbp1, insulin-like growth factor binding protein 1; Lpin1, lipin 1, and serpinb9, serpin peptidase inhibitor, clade B [ovalbumin], member 9) are shown in Supplementary table 1. Real-time PCR reactions using SYBR Green in a final volume of 25 μl were performed using a CFX96 (Bio-Rad Laboratories Ltd, Mississauga, ON). PCR efficiency was examined using standard curves for each gene. Complementary DNAs from all treatment groups were pooled, except for Cyp1a1 where pools from the highest dose animals were used. Relative gene expression was calculated using the CFX Manager Software, with RNA loading correction based on the expression of two housekeeping genes, beta actin (Actb) and hypoxanthine-guanine phosphoribosyltransferase (Hprt).

Testis RNA Extraction and Real-time Quantitative Reverse Transcription-PCR

Frozen left testes (n = 9–10 per group) were disrupted using a mortar and pestle followed by homogenization via a 20 gauge needle and syringe combined with QIAshredder (Qiagen Inc.). RNA extraction was carried out using an RNeasy Plus Mini Kit with on-column DNase digestion (Qiagen Inc.) according to the manufacturer’s protocol. RNA concentrations were determined using the NanoDrop 2000 (ThermoFisher Scientific), and integrity was assessed by conventional gel electrophoresis.

Primer sequences for genes involved in steroidogenesis (Star, steroidogenic acute regulatory protein; Cyp17a1, cytochrome P450 17-hydroxylase/lyase; Ar, androgen receptor; Srd5a1, steroid 5-reductase 1; Srd5a2, steroid 5-reductase 2; Cyp19a1, cytochrome P450, family 19, subfamily a, polypeptide 1; Esr1, estrogen receptor 1, and Esr2, estrogen receptor 2) are shown in Supplementary table 2. Real-time PCR reactions were run in duplicate using QuantiTect One-Step SYBR Green quantitative reverse transcription PCR (qRT-PCR) kits (Qiagen Inc.) in a volume of 20 μl with the Roche LightCycler system (Roche Diagnostics, Laval, QC). PCR efficiency was examined with standard curves for each gene utilizing RNA from control males. Relative gene expression was normalized to 18S ribosomal RNA.

Drug-Metabolizing Enzyme Activities

A 0.5 g piece of liver (n = 9–10 per group) was homogenized in 2.5 volumes (vol/wt) of cold 0.05M Tris-1.15% KCl buffer (pH 7.4). The 10,000 × g supernatant (S9 fraction) was retained and assayed for the activities of CYP1A1/1A2, 2B, 1A2/3A, and 2B/3A, using ethoxyresorufin-(EROD), pentoxyresorufin-(PROD), methoxyresorufin-(MROD), and benzyloxyresorufin-O-dealkylase (BROD) activities, respectively, according to the method of Burke , with modifications. Reactions were carried out at 37°C with triplicate samples in 96-well plates. The fluorescence of resorufin was measured with a SpectraMax M2 (Molecular Devices, Sunnyvale, CA), with excitation and emission wavelengths of 530 and 585 nm. Fluorescence values were converted to nanomoles using standard resorufin fluorescence. Protein content of the S9 fraction was measured using the Quick Start Bradford Protein Assay (Bio-Rad Laboratories, Hercules, CA) and a standard curve of bovine serum albumin.

Hormone Measurements

Serum testosterone concentrations were measured using a testosterone ELISA kit (IBL Immuno-Biological Laboratories, Hamburg, Germany), following the manufacturer’s protocol. Commercial radioimmunoassay kits were used to assay serum thyroxine (T4; MP Biomedicals, Lachine, QC) and TSH (MP Biomedicals). Sample size was 8–10 per group.

Histology

Testis.

Testicular sections (5 μm) were cut and mounted on charged slides and stained with hematoxylin and eosin (Sigma-Aldrich, Oakville, ON) for overall morphological evaluation (n = 5 per group). Images were captured using a Leica DM LB2 microscope (Leica Microsystems Inc., Concord, ON) using a 20× lens fitted with an Infinity-3 video camera (Lumenera Corp., Ottawa, ON).

Thyroid gland.

Fixed thyroid glands in situ on the trachea, were embedded in paraffin, sectioned (5 μm), and stained with periodic acid Schiff. The height of the thyroid gland epithelium was measured in sections known to be either in the interior or periphery of the gland using digital analysis of images. At least 10 images from each location (internal vs. periphery) per animal were analyzed, and measurements of epithelial height were taken from at least five separate follicles per image, selected using an overlaid, numbered grid, and five randomly generated numbers. All measurements were made blindly regarding the animal treatment (n = 9–10 per group) by a single observer.

Sperm Count and Motility

Portions of the frozen left testes (n = 10 per group) were thawed on ice and homogenized in 5 ml of 0.9% saline, 0.1% merthiolate, and 0.05% Triton X-100, for two intervals of 15 s, separated by a 30-s interval on ice. Sperm heads (condensed spermatids and spermatozoa) were counted using a hemocytometer (Robb ).

Spermatozoal motility was measured by CASA using our previously published method (Zubkova and Robaire, 2004). Briefly, spermatozoa from the cauda epididymidis were collected into 3 ml of motility buffer (Hanks Balanced Salt Solution, 0.35 mg/ml sodium bicarbonate, 0.1 mg/ml sodium pyruvate, 0.025 mg/ml soybean trypsin inhibitor, 2 mg/ml bovine serum albumin, pH 7.4 [Gibco Invitrogen, Burlington, ON], 4.2 mg/ml 4-(2-hydroxyethyl)-1-piperazineethanesulfonic acid, 0.9 mg/ml d-glucose [Sigma Chemical Co.], prewarmed to 37°C) and were allowed to disperse for 3 min. A volume of 20 μl of a 1:30-diluted suspension was loaded onto 2X Cel Sperm Analysis Chambers (Hamilton-Thorne Research, Beverly, MA), prewarmed to 37°C and analyzed using a Hamilton-Thorne IVOS Motility Analyzer, version 12 (Hamilton-Thorne Research); the following primary parameters were obtained: average path velocity (VAP), curvilinear velocity (VCL), straight line velocity (VSL), amplitude of lateral head displacement (ALH), as well as two derived parameters, linearity (LIN = BSL/VCL × 100) and straightness (STR = VSL/VAP × 100). Approximately 10,000–14,000 cauda epididymal spermatozoa were analyzed for each treatment group (n = 9–10 per group). Analysis was performed using the following settings: stage temperature, 37°C; frames acquired, 30; frame rate, 60 Hz; minimum contrast, 80; minimum cell size, 4 pixels; minimum static contrast, 15; cell intensity, 80; magnification, 0.82; static size limits, 0.29–8.82; and static intensity, 0.18–1.8. The chamber depth was 80 μm, where 20 μl were loaded into each chamber.

Sperm DNA Integrity

Comet assay.

DNA strand breaks in spermatozoa were evaluated using the alkaline comet assay, as previously described (Codrington ). Briefly, frozen sperm samples were thawed, diluted, and 50 μl of the cell suspension (1 × 105 cells/ml) were mixed with molten agarose, and samples were placed on a two-well CometSlide (Trevigen Inc., Gaithersburg, MD). Slides were immersed in lysis buffer containing 10% dimethyl sulfoxide, 1% Triton X-100, and 40mM dithiothreitol for 1 h on ice and then incubated for 3 h at 37°C in lysis buffer containing 0.1 mg/ml proteinase K. Slides were washed, kept at 4°C for 10 min to reset the agarose, and immersed in freshly prepared alkaline solution (1mM EDTA and 0.05M NaOH, pH 12.1) for 45 min. Electrophoresis was done at 14 V (0.7 V/cm) for 10 min in 1× Tris-borate-EDTA buffer, pH 8 (Mini-Sub Cell GT; Bio-Rad Laboratories). Slides were then fixed in ethanol and stored at room temperature. DNA was stained with SYBR Green (Trevigen Inc.) (1:10,000 in Tris-EDTA buffer, pH 7.5) and immediately analyzed using an Infinity 3-1 digital camera (Lumenera Corp.) attached to a Leica DM LB2 fluorescence microscope. Fifty cells per well were selected at random and analyzed, for a total of 100 cells per animal (n = 5 per group), and fluorescent images were scored for comet parameters. Tail length, percent tail DNA, and tail extent moment (tail length/fraction of tail DNA) were assessed using the KOMET 5.0 image analysis system (Kinetic Imaging Ltd, Liverpool, U.K.).

Acridine orange/sperm chromatin structure assay.

To assess the susceptibility of sperm nuclear DNA to acid-induced denaturation, the acridine orange (AO) assay, also known as the sperm chromatin structure assay (SCSA) (Evenson ), was applied, as previously described (Delbes ; Evenson ). Briefly, sperm samples (2.5 × 106 cells/ml in PBS) were thawed, sonicated on ice, and mixed with denaturation buffer (0.08 N HCl, 0.15M NaCl, and 0.1% Triton X-100, pH 1.4) for 30 s at 4°C to denature uncondensed sperm DNA. After 30 s, 1.2 ml of AO staining solution (0.126M Na2HPO4, 0.037M citric acid buffer, 1mM EDTA, 0.15M NaCl, pH 6.0 containing 6 μg/ml AO [Sigma Chemical Company]) was added. Exactly 3 min following the addition of the denaturation buffer, spermatozoa were analyzed using a FACScan flow cytometer (BD Biosciences) fitted with an argon ion laser (488-nm line excitation). A positive control was obtained by preincubating the spermatozoa with 100mM H2O2 at 56°C for 45 min. The raw data were processed using WinList software (Verity Software House, Topsham, ME). Results were expressed as mean DNA fragmentation index (DFI), standard deviation of DFI, percentage of DFI (% DFI, corresponding to the percentage of cells outside the main population), and as percentage of spermatozoa with high green fluorescence or DNA stainability (% HDS), as an indication of sperm DNA compaction (Evenson ). A total of 10,000 AO-labeled spermatozoa per sample (n = 5 per group) was analyzed and each sample was run in duplicate.

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. The level of significance was p < 0.05. For thyroid epithelium measures, data were analyzed by a two-way ANOVA followed by post hoc Duncan's multiple range test. Analyses were done using SigmaPlot version 12 (Systat Software, Chicago, IL) or JMP statistical software version 4.0.2 (SAS Institute, Cary, NC) except for K-W ANOVA on ranks and post hoc analyses, which were performed using SigmaStat version 3.5 (Systat Software).

RESULTS

The 70-day exposure to the BFR mixture did not affect survival; rats appeared normal, with no hair loss, impaired mobility, or hunched posture. Body weight gain (Fig. 1A), absolute body weight (Supplementary figure 1A), as well as weekly food consumption (Supplementary figure 1B), were unaffected. Based on an average food consumption of 30 g per day, daily exposure for each treatment group was estimated at 0.014, 0.14, 1.39, and 14.0 mg/kg/day (Supplementary table 3), resulting in 68–70% of target doses. Therefore, the lowest dose to which the rats were exposed was 0.007 mg/kg pentaBDE, 0.006 mg/kg decaBDE, 0.0005 mg/kg HBCD, and 0.00005 mg/kg octaBDE.

FIG. 1.

Effects of 70-day BFR exposure on body weights and liver and kidney weights in male rats. (A) Change in body weight during the BFR treatment. Male rats were randomly divided into five treatment groups: 0 (control), 0.02, 0.2, 2, and 20 mg/kg/day. Animals were weighed weekly for 70–72 days and the overall change in body weight during this period is shown. Bars represent the means ± SEM (n = 15 per group). At the end of BFR treatment, (B) liver and (C) kidneys were weighed and normalized to body weight (BW, g/kg BW). Values are expressed as mean ± SEM (n = 15 per group). *p < 0.05 compared with control.

The highest BFR dose significantly increased relative liver (Fig. 1B) and kidney (Fig. 1C) weights compared with control. There were no effects of exposure on spleen, thymus, lung, or heart weights (data not shown). To determine the effects of BFR on various organ systems, a variety of parameters of serum biochemistry were assessed (Table 2). Of these, only glucose was significantly reduced at the highest BFR exposure level (p = 0.03). Uric acid showed a dose-related decrease, but this effect was not statistically significant (p = 0.096; Table 2).

TABLE 2

Serum Biochemistry

	BFR mixture nominal doses (mg/kg/day)
Parameter	Control (0)	0.02	0.2	2	20	p
Calcium (mg/dl)	13.7 ± 0.2 (10)	13.6 ± 0.2 (8)	13.2 ± 0.3 (9)	13.6 ± 0.2 (10)	13.4 ± 0.2 (10)	0.51
Magnesium (mg/dl)	3.40 ± 0.10 (8)	3.65 ± 0.10 (7)	3.49 ± 0.09 (7)	3.52 ± 0.11 (9)	3.42 ± 0.10 (8)	0.51
Phosphorus (mg/dl)	9.11 ± 0.34 (8)	9.71 ± 0.30 (6)	9.65 ± 0.54 (7)	10.2 ± 0.42 (9)	9.38 ± 0.26 (10)	0.30
Blood urea N	17.1 ± 0.82 (9)	17.2 ± 0.53 (9)	17.1 ± 0.78 (9)	17.5 ± 0.55 (10)	17.5 ± 0.49 (10)	0.99
Uric acid (mg/dl)	3.18 ± 0.34 (10)	2.62 ± 0.25 (9)	2.89 ± 0.27 (9)	2.47 ± 0.31 (10)	2.09 ± 0.25 (10)	0.096
Albumin (g/dl)	3.40 ± 0.05 (8)	3.41 ± 0.10 (4)	3.42 ± 0.04 (7)	3.38 ± 0.05 (8)	3.44 ± 0.05 (6)	0.96
Total protein (g/dl)	7.77 ± 0.10 (10)	7.97 ± 0.12 (8)	7.81 ± 0.13 (9)	7.91 ± 0.14 (10)	7.95 ± 0.11 (10)	0.71
Creatinine kinase (U/l)	227 ± 22 (9)	179 ± 14 (9)	215 ± 29 (9)	275 ± 41 (10)	383 ± 87 (10)	0.20
Creatinine (mg/dl)	0.53 ± 0.02 (9)	0.49 ± 0.02 (9)	0.50 ± 0.01 (9)	0.49 ± 0.01 (10)	0.49 ± 0.01 (10)	0.37
Triglycerides (mg/dl)	184 ± 29 (8)	226 ± 29 (8)	177 ± 19 (7)	248 ± 49 (10)	173 ± 35 (8)	0.48
Cholesterol (mg/dl)	109 ± 5 (10)	117 ± 7 (9)	118 ± 8 (9)	118 ± 7 (10)	127 ± 8 (10)	0.48
Glucose (mg/dl)	296 ± 16 (9)	283 ± 9 (9)	289 ± 16 (8)	255 ± 15 (10)	238 ± 15* (9)	0.03
ALP (U/l)	184 ± 16 (8)	199 ± 21 (6)	202 ± 17 (8)	191 ± 18 (10)	187 ± 9 (8)	0.95
ALT (U/l)	67.9 ± 3.6 (10)	64.9 ± 3.6 (9)	75.4 ± 6.3 (8)	65.4 ± 2.1 (10)	69.5 ± 3.8 (10)	0.40

Notes. Values are expressed as mean ± SEM of (n) observations per treatment. Values in bold are significantly different from control.

*p < 0.05 compared with control.

To further evaluate the effects of BFR exposure, RNA levels and the catalytic activities of drug-metabolizing enzymes in the liver were measured. The highest BFR dose significantly increased Cyp1a1 messenger RNA (mRNA) levels by 50-fold (Fig. 2) and its EROD catalytic activity by 4.6-fold (Fig. 3A), in addition to increasing MROD (CYP2B activity, 3.2-fold) and PROD (CYP2B/3A activity, 5.6-fold) activities, compared with controls (Figs. 3B and C). The CYP1A2/3A-mediated BROD activity showed a dose-dependent increase in both the 2 (2-fold) and 20 (8.9-fold) mg/kg/day BFR exposure groups (Fig. 3D), suggesting that BROD may represent a more sensitive marker of response to our BFR mixture.

FIG. 2.

Effects of BFR exposure on liver Cyp1a1 mRNA levels. Relative levels of hepatic Cyp1a1 mRNA were determined by qRT-PCR and normalized to both Actb and Hprt, then compared with control levels. Values are expressed as mean ± SEM (n = 9–10 per group). *p < 0.01 compared with control.

FIG. 3.

Effects of BFR exposure on the hepatic activities of drug-metabolizing enzymes. The activities of ethoxy-(EROD) (A), methoxy-(MROD) (B), propyloxy-(PROD) (C), and benzyloxy-(BROD) (D)-resorufin-O-deethylation were measured in hepatic S9 fractions prepared from frozen liver tissues, as described in the text. Activities were corrected for S9 protein levels and are expressed as mean ± SEM (n = 9–10 per group). *p < 0.05 compared with control.

Various male reproductive parameters were evaluated following the 70-day BFR treatment. Androgen-responsive organ weights (Figs. 4A–E) and serum testosterone levels (Fig. 4F) were unaffected by BFR exposure. Although testosterone levels appeared to increase with BFR exposure, this did not reach statistical significance (p = 0.78); the variability observed may be due to the pulsatile nature of testosterone secretion (Bartke ). BFR exposure had no observable effect on testicular sperm count (Supplementary figure 2). Testicular histology did not reveal significant morphological changes in seminiferous tubular structure or in cellular composition between controls and animals in the highest BFR dose group (Supplementary figure 3).

FIG. 4.

Reproductive organ weights and serum testosterone levels in control and BFR-treated males. At the end of treatment, (A) epididymis, (B) testes, (C) empty seminal vesicles, (D) prostate, and (E) vas deferens were weighed. Blood was collected and (F) serum testosterone concentrations were measured by ELISA. Values are expressed as mean ± SEM (n = 8–10 per group).

Motility of spermatozoa isolated from the cauda epididymidis was also analyzed. There were no effects on sperm motility or motion parameters (Supplementary table 4). As toxicants may still affect sperm chromatin integrity without observable effects on sperm counts, motility, or testicular histology (Delbes ), DNA integrity in caudal sperm was also assessed. BFR treatment did not induce sperm DNA single or double strand breaks, as measured by the COMET assay (Supplementary figure 4), and did not affect the susceptibility of the spermatozoa to low pH denaturation, as assessed by the SCSA assay (Supplementary figure 5).

To further assess the impact of BFR treatment on testicular function, the RNA expression of key genes involved in the steroidogenesis pathway was examined by qRT-PCR (Supplementary figure 6). Although there was an overall treatment effect on Ar gene expression, the differences were not significant compared with control and were not dose related. No significant changes in the relative mRNA expression of any of the other genes were observed.

To assess BFR effects on thyroid physiology, serum T4 and TSH levels were measured (Fig. 5). In the group exposed to the highest BFR dose, T4 levels were significantly reduced to 73% of control levels (p < 0.0001; Fig. 5A), but TSH levels were unchanged (Fig. 5B). Although serum thyroid hormone levels indicate immediate effects of chemical exposure on the hypothalamic-pituitary-thyroid (HPT) axis, thyroid gland histomorphology provides a more integrated and persistent indicator of adverse chemical effects on thyroid hormone physiology. Qualitative evaluation of the inner follicles of thyroid glands did not reveal pathological changes induced by BFR exposure with the exception of increased vacuolation of the luminal apices of epithelial cells apparent in high dose animals (Fig. 6). Thyroid epithelial thickness (a reliable indicator of TSH stimulation) was measured in two general locations in the thyroid gland (inner and outer follicles) to determine if these sites, known to have different levels of physiological activity (Low, 1982), respond differently to BFRs. As expected, inner follicles had thicker epithelia than outer follicles (9.85 ± 1.2 vs. 4.74 ± 0.99 μm for inner vs. outer, respectively, p < 0.001) across all doses, reflecting the increased sensitivity of the inner follicles to TSH. Exposure to the BFR mixture caused a biphasic response in epithelial thickness in the inner follicles, with a reduction in thickness in response to lower doses but a significant increase at the highest dose (Fig. 7A). Follicles at the periphery showed no change in epithelial thickness in response to any dose of the flame retardant mixture (Fig. 7B). To determine if BFR exposure influenced thyroid hormone action at the level of a target tissue, as reported for structurally similar compounds such as polychlorinated biphenyls (PCBs) (Bansal and Zoeller, 2008), we investigated the expression of several thyroid hormone–regulated genes in the liver (Table 3). Gene expression of lipogenic enzyme Me1, known to be induced by thyroid hormone stimulation (Dozin ), was significantly increased in animals fed the highest BFR-supplemented diet, despite reduction in circulating T4 levels (Fig. 5). No dose of BFR significantly altered mRNA levels of either Thrsp or Dio1 (Table 3), although both are positively regulated by thyroid hormone action (Kohrle, 2002; Narayan ). Esr-responsive genes were also assessed because liver Me1 expression is induced by estrogenic substances (Schmutzler ); their mRNA levels were unaffected by the BFR treatment (Table 3).

FIG. 5.

Effects of BFR exposure on circulating (A) thyroxine (T4) and (B) TSH. Serum collected at termination was assayed for both hormones, as described in the text. Values are expressed as mean ± SEM (n = 9–10 per group). *p < 0.05 compared with control.

FIG. 6.

Light micrographs (×200, periodic acid Schiff staining) of the inner follicles of thyroid glands of control (A) and BFR-treated males (B) 0.02, (C) 0.2, (D) 2.0, or (E) 20 mg/kg/day. “c” indicates follicular colloid in the acellular lumen, whereas the arrow indicates the single layer of cuboidal epithelium that surrounds and encloses the follicle. Bar = 50 μm.

FIG. 7.

Effects of BFR exposure on height of the (A) inner and (B) outer thyroid follicular epithelium. Values are expressed as mean ± SEM (n = 9–10 per group). Bars with the same letter superscript are not significantly different.

TABLE 3

Hepatic Gene Expression

	BFR mixture nominal dose (mg/kg/day)
	Control (0)	0.02	0.2	2.0	20	p
Gene
Thyroid hormone receptor responsive
Me1	1.00 ± 0.19 (9)	1.02 ± 0.16 (10)	0.94 ± 0.13 (7)	1.70 ± 0.38 (9)	2.06 ± 0.28* (8)	0.006
Dio1	1.00 ± 0.09 (9)	0.86 ± 0.07 (10)	0.82 ± 0.09 (9)	1.09 ± 0.15 (10)	1.04 ± 0.08 (10)	0.27
Thrsp	1.00 ± 0.20 (9)	0.93 ± 0.11 (10)	0.92 ± 0.14 (9)	1.25 ± 0.39 (10)	1.05 ± 0.18 (10)	0.99
Estrogen receptor responsive
Apoa4	1.00 ± 0.12 (9)	0.89 ± 0.10 (10)	0.83 ± 0.04 (9)	1.12 ± 0.16 (10)	0.84 ± 0.05 (10)	0.56
Igfbp1	1.00 ± 0.21 (9)	1.35 ± 0.29 (10)	1.01 ± 0.19 (9)	1.52 ± 0.27 (10)	0.94 ± 0.15 (10)	0.32
Lpin1	1.00 ± 0.20 (9)	0.97 ± 0.24 (10)	0.81± 0.12 (9)	0.91 ± 0.24 (10)	1.19 ± 0.26 (10)	0.81
Serpinb9	1.00 ± 0.19 (9)	0.94 ± 0.15 (10)	0.93 ± 0.15 (9)	0.86 ± 0.14 (10)	1.01 ± 0.06 (10)	0.94

Notes. Values are expressed as mean levels of RNA, normalized to control levels, ± SEM of (n) observations per treatment. Values in bold are significantly different from control.

*p < 0.05 compared with control.

DISCUSSION

Adult male rats exposed for 70 days to a high dose of a complex BFR mixture mimicking the relative congener levels in house dust showed enlargement of the liver and kidney and altered thyroid hormone parameters; no impact on the reproductive system was noted. This BFR mixture was incorporated into experimental diets to deliver nominal exposures of 0.02, 0.2, 2, and 20 mg/kg. Based on the measured food consumption, rats ingested an average of 0.014, 0.14, 1.39, and 14.0 mg/kg/day or approximately 70% of our target doses.

Liver and kidney weights were increased by 18 and 12%, respectively, at the highest BFR exposure level; these increases are within the typical range for xenobiotic exposure (Maronpot ) and indicative of toxicity. However, these changes were not accompanied by increases in ALT and ALP, the usual serum biological markers indicative of liver damage, or in blood urea nitrogen, creatinine, and uric acid, indicative of kidney damage. These results are not unusual because circulating levels of these markers may not correlate well with the degree of organ injury, and significant increases are often associated with severe and irreversible injury (Ennulat ).

Liver enlargement was associated with increased activities in EROD/MROD and PROD activities, suggesting induction of CYP1A1/1A2 and CYP2B P450 drug-metabolizing enzymes, respectively. Induction of BROD, catalyzed by CYP1A, CYP2B, and CYP3A (Burke ), was observed in both the 2 and 20 mg/kg/day BFR dose groups. Although BROD does not provide a definitive measure of any single enzyme, it does indicate that phase I enzymes were induced at the second highest dose. Induction of EROD and PROD activities was also reported for male rats treated with pentaBDE (van der Ven ) or decaBDE (van der Ven ), whereas only PROD was induced in HBCD-treated males (Germer ). Induction of CYP1A and CYP2B suggests activation of the aryl hydrocarbon receptor (AhR) and the constitutive androstane receptor (CAR) pathways, respectively (Honkakoski and Negishi, 2000).

The significant decrease in serum glucose that was observed may be secondary to repression of gluconeogenesis through activation of the CAR pathway. This has also been observed with phenobarbital-like compounds (Rushmore and Kong, 2002). Alternatively, it may be secondary to thyroid hormone–induced alterations in transcription of the glucose transporter gene (Chidakel ). Our results are consistent with observations reported after exposure to pentaBDE (van der Ven ) and HBCD (van der Ven ).

Exposure to the BFR mixture throughout all stages of male germ cell development in adult rats had no effect on the reproductive system, as assessed by a wide range of parameters; these included the weights of reproductive organs, testicular function (testosterone levels, histology, sperm count and motility, and gene expression), and sperm DNA integrity (COMET assay and SCSA). Our results are consistent with observations reported after exposure of adult male rats to several technical BFR mixtures. No effects were seen on sperm counts (van der Ven , 2008a,b), and reproductive organ weights were affected after exposure to relatively high doses of pentaBDE (starting from 27.5 mg/kg/day) (van der Ven ). In contrast, low doses of decaBDE (starting from 0.2 mg/kg/day) were reported to increase seminal vesicle weights (van der Ven ). There is evidence that male rats exposed during development or weaning are likely to be more sensitive than adults to BFRs (Stoker ; van der Ven ). Indeed, maternal breast milk PBDE levels were directly correlated with the incidence of cryptorchidism in 3-month-old boys (Main ). Further studies will be necessary to establish the potential hazards of the BFR mixture to the developing male reproductive tract.

In contrast to the apparent insensitivity of the adult male reproductive system to BFRs, treatment with our BFR mixture reduced serum total T4; paradoxically, this was not accompanied by a compensatory increase in TSH. A variety of other substances known to activate the AhR or CAR pathways also have been reported to reduce circulating T4 levels with minimal or no compensatory increase in TSH; these include Aroclor 1254, a mixture of PCBs, phenobarbital, 3-methylcholanthrene (Barter and Klaassen, 1994), and 2,3,7,8-tetrachlorodibenzo-p-dioxin (Schuur ).

Histomorphological analysis of the thyroid gland epithelium revealed a biphasic effect of the BFR mixture on the HPT axis, suggesting a thyroid hormone–like action of the mixture, at least at the lowest dose. Epithelial thickness is largely a result of the trophic action of TSH (Capen and Martin, 1989) and provides an integrated measure reflecting the culmination of the recent history of TSH, whereas serum levels reflect a single point. Thyroid epithelial thickness was reduced at the lowest dose of the BFR mixture, suggesting a reduction in trophic stimulation in the gland, and was increased at the highest dose to significantly exceed control thickness. Although serum TSH levels did not change with any dose of the mixture, mechanistically relevant changes in the target organ of this hormone suggest that the cumulative action of TSH up until sacrifice was significantly altered by the mixture. The reduction in epithelial height at the lowest dose suggests the possibility that (a) component(s) of the BFR mixture was acting like thyroid hormone in suppressing TSH release. However, the lack of change in expression of any thyroid hormone–responsive genes (DioI and Thsrp) in the liver or genes in the ventricles (Myh6, myosin, heavy chain 6, cardiac muscle, alpha; Myh7, myosin, heavy chain 7, cardiac muscle, beta; Pln, phospholamban Atp2a2, ATPase, Ca++ transporting, cardiac muscle, and slow twitch 2; data not shown) at any dose does not support the argument that BFRs are being converted to thyroid hormone receptor agonists.

The sole thyroid hormone–responsive gene to show any response was Me1 and its transcript was elevated in the livers of animals consuming the highest BFR dose. The lack of effect on other hepatic thyroid hormone–responsive genes suggests that some of the Me1 response was not mediated through the activation of a thyroid receptor. Hepatic transcription of Me1 is also increased by treatment with estrogenic substances in vivo (Schmutzler ). Some PBDE metabolites have been shown to activate Esr (Kojima ), but the lack of change in hepatic estrogen–responsive genes in the current study argues against Esr activation by BFRs in the mixture. Specific induction of hepatic Me1 gene expression has also been reported in rats exposed to other organochlorines, such as hexachlorobenzene (HCB) (Loaiza-Perez ). The presence of a thyroid hormone response element in the promoter region was essential for the HCB-induced increase in Me1 transcription in hepatoma cells in vitro. This would suggest that the thyroid hormone receptor is required for this effect (Loaiza-Perez ), although HCB does not bind to the ligand-binding domain of the thyroid hormone receptor (Li ). Altogether, these results suggest that the expression of Me1 is influenced by BFRs and other organochlorines through an uncharacterized interaction with the thyroid hormone receptor and thyroid hormone–responsive elements.

In conclusion, we found no evidence that exposure to an environmentally relevant mixture of BFRs commonly found in North American residences presents a risk to adult male rat reproductive function. Further analyses are currently underway to determine whether this exposure produced serum levels similar to those found in humans. Nevertheless, the absence of effects on any of diverse measures of male reproductive function across a broad range of doses, including those that induced xenobiotic metabolism activity in the liver, provides some assurance that adult BFR exposure does not influence male fertility. These data provide a robust dataset for assessing the hazard of this exposure. At the highest dose, multiple end-points (reduced serum T4 and increased epithelial thickness and hepatic Me1 expression) of thyroid physiology were adversely affected. Thus, the weight of evidence argues that thyroid hormone physiology is influenced in adult male rats by high dose exposure to an environmentally relevant mixture of BFRs.

SUPPLEMENTARY DATA

Supplementary data are available online at http://toxsci.oxfordjournals.org/.

FUNDING

Canadian Institutes of Health Research (CIHR), Institute for Human Development, Child and Youth Health (RHF100625).