Next Article in Journal
Estimated Dietary Intake of Trace Metals from Swordfish Consumption: A Human Health Problem
Previous Article in Journal
Mitochondrial Morphology and Function of the Pancreatic β-Cells INS-1 Model upon Chronic Exposure to Sub-Lethal Cadmium Doses
Previous Article in Special Issue
Species and Sex Differences in the Morphogenic Response of Primary Rodent Neurons to 3,3′-Dichlorobiphenyl (PCB 11)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Transcriptomic Analysis of Gonadal Adipose Tissue in Male Mice Exposed Perinatally to 2,2′,4,4′-Tetrabromodiphenyl Ether (BDE-47)

Department of Environmental Health Sciences, University of Massachusetts Amherst, 686 N. Pleasant St., 149A Goessmann, Amherst, MA 01003, USA
*
Author to whom correspondence should be addressed.
Toxics 2018, 6(2), 21; https://doi.org/10.3390/toxics6020021
Submission received: 7 March 2018 / Revised: 21 March 2018 / Accepted: 23 March 2018 / Published: 29 March 2018
(This article belongs to the Special Issue Polyhalogenated Aromatic Hydrocarbons)

Abstract

:
For the majority of lipophilic compounds, adipose tissue is traditionally considered as a storage depot and only rarely as a target organ. Meanwhile, abnormalities in adipose tissue physiology induced by chemical exposure may contribute to the current epidemic of obesity and metabolic diseases. Polybrominated diphenyl ethers (PBDEs) are a group of lipophilic flame retardants found in the majority of human samples in North America. Their ability to alter the physiology of adipose tissue is unknown. We exposed pregnant mice to 0.2 mg/kg body weight/day of BDE-47 perinatally. Transcriptomic changes in gonadal adipose tissue were analyzed in male offspring using the RNA-seq approach with subsequent bioinformatic analysis. The expression of genes of coagulation and complement cascade, de novo lipogenesis, and xenobiotic metabolism was altered in response to BDE-47 exposure. The affected molecular network included the following hubs: PPARα, HNF1A, and HNF4. These findings suggest that adipose tissue should be considered a target tissue for BDE-47, in addition to its role as a storage depot. This study also builds a background for a targeted search of sensitive phenotypic endpoints of BDE-47 exposure, including lipid profile parameters and coagulation factors in circulation. Additional studies are needed to investigate the role of PBDEs as an obesogen.

1. Introduction

Metabolic disease represents one of the biggest public health challenges today. Approximately two out of every five Americans will develop type 2 diabetes at some point during their adult lives [1]. More than 60% of children above the age of 10 either are obese or will become obese later in life. One in four overweight children suffers from impaired glucose tolerance [2]. A growing body of evidence links the current worldwide epidemics of obesity, diabetes, and related metabolic diseases with environmental exposure [3]. Many environmental xenobiotics, such as dioxins, PCBs, PBDEs, and others, are highly lipophilic and therefore accumulate in adipose tissue in humans, livestock, and wildlife. Adipose tissue is usually considered a storage depot of lipophilic compounds and is only rarely considered a target of toxicity in toxicological research [4], although it is known to be involved in several physiological functions, including metabolic regulation, energy storage, and endocrine functions [4]. For example, white adipose tissue is a major secretory cell of a diverse list of proteins collectively known as adipokines involved in the regulation of metabolic functions of the organism [5]. Changes in adipose physiology induced by environmental xenobiotics may be associated with obesity, diabetes, and other metabolic diseases [6,7].
Polybrominated diphenyl ethers (PBDEs) are a type of brominated flame retardants that were used in many household products such as paints, plastics, textiles, furniture padding, mattresses, and various electrical appliances until recently [8]. Although the commercial use of PBDEs has ceased in many industrial countries, PBDEs are ubiquitously found in the environment due to their resistance to degradation. The level of PBDEs in human blood, breast milk, and adipose tissue has increased exponentially over the last three decades and continues to rise in North America [9]. Among the various PBDE congers, 2,2′,4,4′-tetrabromodiphenyl ether (BDE-47) is the most prevalent in human samples, including maternal milk [10,11,12,13,14,15]. PBDEs are lipophilic and bioaccumulate in fatty tissues with a half-life in human tissues estimated at 1.8–6.5 years [16]. Although PBDEs accumulate predominantly in adipose tissue, the majority of studies of PBDE toxicity were focusing on neurodevelopment, endocrine, and reproductive endpoints [17,18,19] and target organs such as the liver, thyroid gland, kidney, and reproductive organs. The ability of PBDEs to affect the physiology of adipose tissue was not thoroughly addressed in previous studies.
In this study, we analyze global gene expression in gonadal adipose tissue of mice exposed perinatally to a low dose (0.2 mg/kg body weight) of BDE-47. Our data demonstrate significant changes in several functional networks including coagulation and complement cascade, lipid and lipoprotein synthesis, and xenobiotic metabolism pathways.

2. Materials and Methods

2.1. Animals and Treatment

Eight-week-old male (30–35 g) and female (27–30 g) CD-1 mice were obtained from Charles River Laboratories (Kingston, NY, USA) and housed in a temperature (23 ± 2 °C)- and humidity (40 ± 10%)-controlled environment, with a 12-h light/dark cycle, and food and water available ad libitum. After three days of acclimation, the animals were bred and the day of vaginal plug detection was considered pregnancy day 1. The dams were assigned to one of two treatment groups (n = 5 per group) based on weight match, and were exposed to tocopherol-stripped corn oil (MP Biomedicals, Solon, OH, USA) or a 0.2 mg/mL solution of BDE-47 (AccuStandard, Inc., New Haven, CT, USA; 100% purity) in tocopherol-stripped corn oil daily from pregnancy day 8 through postpartum day 21; the females were fed 1 µL/g BW from the tip of a pipette, resulting in exposure of 0.2 mg/kg BW/day. This method of exposure is routinely used in our laboratory as a substitution of oral gavage as the latter method induces significant stress response by the endocrine system, which may interfere with the analyzed health outcomes [20]. The dams were allowed to deliver naturally, and the litters were not culled to maintain consistency of nutrient distribution among the same number of fetuses/pups at pre- and postnatal periods, and to avoid catch-up growth [21]. The dams and pups were kept together until weaning on PND 21, when the male and female pups were separated. Using cervical dislocation, one randomly selected male pup per litter was euthanized on PND21. The other pups were used in a different study. All animals were euthanized between 9 a.m. and 11 a.m., following 2 h of fasting. Gonadal fat from the right side of each animal was collected, weighed, and snap-frozen in liquid nitrogen and stored at −80 °C. All procedures met the guidelines of the National Institutes of Health Guide for the Care and Use of Laboratory Animals, and this study was approved by the Institutional Animal Care and Use Committee at University of Massachusetts, Amherst, MA, USA.

2.2. RNA Extraction and Sequencing

Four samples of gonadal fat per exposure group were used for RNA extraction, collected from four individual pups from four different litters per exposure group. Total RNA was isolated using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) and quantified using a NanoDrop 1000 instrument (Thermo Fisher Scientific, Wilmington, DE, USA). RNA quality was assessed using Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA). Samples of RNA with integrity values > 9 were used for library preparation in accordance with the recommendations of the library preparation and sequencing kit manuals. The NEBNext Poly(A) mRNA Magnetic Isolation Module (E7490, New England BioLabs Inc., Ipswich, MA, USA) was then used to isolate intact poly(A) + RNA from 3 µg of total RNA. Libraries were constructed using NEBNext mRNA Library Prep Reagent Set for Illumina with multiplexing indexes from NEBNext Multiplex Oligos for Illumina (E6100 and E7335 respectively, New England BioLabs Inc., Ipswich, MA, USA). The quality and purity of the libraries were assessed by Agilent 2100 Bioanalyzer, and the concentration of the libraries was measured by real-time PCR with primers for P5 and P7 flow cell oligo sequences using a KAPA Library Quantification Kit (KR0405, Kapa Biosystems, Boston, MA, USA) in a 384-well plate on a CFX384 Touch Real-Time PCR Detection System (Bio-Rad, Hercules, CA, USA). High-throughput sequencing was performed on a NextSeq500 sequencing system (Illumina, San-Diego, CA, USA) in the Genomic Resource Laboratory of the University of Massachusetts, Amherst. cDNA libraries were single-end sequenced in 75 cycles using NextSeq 500 High Output Kit (FC-404-1005, Illumina, San-Diego, CA, USA) in a multiplex run. Four samples were sequenced per exposure group, representing eight different litters. Sequencing was completed with minimum 20 million reads per sample. One control sample was sequenced with smaller number of reads and was excluded from further analysis. All sequencing data were uploaded to the GEO public repository and GEO were assigned series accession number GSE112207.

2.3. RT-qPCR

The RNA-seq results were validated using RT-qPCR for genes selected randomly from the list of the most significantly differentially regulated genes in RNA-seq experiment (Table 1). Total RNA was purified of genomic DNA contamination using DNase (RQ1 RNAse-free DNAse, Cat. # M610A, Promega, Madison, WI, USA) and reverse-transcribed using the High Capacity cDNA Reverse Transcription Kit (Cat. # 4368814, Applied Biosystems, Vilnus, Lithuania). Using free online software Primer3Plus (primer3plus.com), forward and reverse primers were designed to anneal different exons spanning long intron (Table 1). The housekeeping gene (Polr3c) was selected from genes that were not regulated in our RNA-seq dataset, with the consideration of its ubiquitous presence in different cell types. Triplicate 5-μL real-time PCR reactions, each containing iTaq Universal SYBR Green Supermix (Cat. 172-5124, BioRad), primers, and cDNA template were loaded onto a 384-well plate and run through 40 cycles on a CFX384 real time cycler (Bio-Rad Laboratories, Inc.). The data were analyzed using the manufacturer’s CFX manager software (version 3.1, BioRad, Hercules, CA, USA, 2013). Relative quantification was determined using the ΔΔCq method [22].

2.4. Bioinformatic Analysis of RNA-seq Data

To analyze the data from our mRNA sequencing experiment, read filtering, trimming and de-multiplexing were performed via the BaseSpace cloud computing service supported by Illumina (https://basespace.illumina.com/home/index). The preprocessed reads were mapped to the reference mouse genome (MM10) using TopHat 2 aligner [23,24]. Aligned reads were then used for assembly of novel transcripts with Cufflinks 2.1.1 and differential expression of novel and reference transcripts with Cuffdiff 2.1.1 [25]. The list of differentially expressed transcripts was filtered to select only transcripts associated with gene names and transcripts with FPKM values > 1 in at least one of conditions (control or exposed). The resulting gene list included 14,289 genes.
DAVID bioinformatics resources version 6.7, Gene Set Enrichment Analysis (GSEA), and Ingenuity Pathway Analysis (IPA) were used for bioinformatics analysis of RNA-seq data. A short list of 2-fold significantly differentially expressed genes was imported to DAVID or IPA for functional analysis with default settings. GSEA analysis was conducted using a ranked list of 14,289 genes against the following collections of gene sets: hallmark genes sets (H), curated gene sets from online pathway databases, publications in PubMed, and knowledge of domain experts (C2) and gene ontology gene sets (C5) obtained from the molecular signatures database v5.2 [26] and PPARα target genes obtained from [27]. Gene sets enriched with FDR q-value < 0.02 were considered significant.

3. Results

We found no significant relationship between litter size and exposure to BDE-47, with the number of pups varying from 11 to 15 per litter. No weight differences were observed between the control and exposed dams and pups throughout the experiment. No changes in adipose tissue weight were found in male mice on PND21.
The resulting filtered list of differentially expressed genes included 14,289 genes. 758 genes were deferentially expressed 2-fold or higher (Figure 1) and 94 genes also had FDR q-value of ≤ 0.05 (Table S1). Expression of select genes measured by RT-PCR showed good concordance with RNA-seq data (Table 1).

3.1. DAVID Analysis

DAVID analysis revealed four enriched annotation clusters with an enrichment score >2.0: ‘secreted proteins’, ‘peptidase inhibitor activity’, ‘complement and coagulation cascades’ (Figure 2), ‘biosynthesis of amino acid’, ‘microsomal metabolism’, and ‘lipid metabolism’ (Table 2).

3.2. GSEA Analysis

No significantly positively enriched gene sets were found in GSEA analysis (see Table S2 for details). The following hallmark gene sets were negatively enriched with FDR q < 0.2: ‘coagulation’, ‘xenobiotic metabolism’, and ‘bile acid metabolism’ (Figure 3A). Analysis of C2 collection of gene sets revealed several negatively regulated liver specific gene sets, including genes targets of HNF1A and HNF4A (Figure 3B). Several positively enriched (although not significant) gene sets indicated upregulation of lipid biosynthesis, including ‘SREBF targets’ (Figure 3C), ‘obesity’ and ‘cholesterol biosynthesis’. Significantly negatively enriched gene ontology (GO) gene sets included ‘coagulation’, ‘serine type endopeptidase activity’, and ‘lipid catabolism’. Genes—targets of PPARα alpha were also significantly downregulated (normalized enrichment score = −1.72, nominal p-value < 0.0001, FDR q < 0.001).

3.3. IPA Analysis

The top IPA enriched canonical pathways were ‘FXR/RXR activation’, ‘LXR/RXR activation’, ‘acute phase response signaling’, ‘coagulation system’, and ‘extrinsic prothrombin activation pathway’ (Table 3). HNF1A (p-value of overlap = 2.98 × 10−37, predicted activation—inhibited), and PPARα (p-value of overlap = 1.41 × 10−18, predicted activation—n/a) were identified as upstream regulators with significantly altered activity (Figure 4).

4. Discussion

We found that perinatal exposure to 0.2 mg/kg body weight of BDE-47 results in changes in the expression of genes involved in the coagulation and complement cascades, xenobiotic metabolism and lipid metabolism in gonadal adipose tissue of mice, suggesting that adipose tissue can also be a target tissue for BDE-47, in addition to its role as a storage depot.

4.1. Relevance of the Dosing Paradigm

Results reported in this study were produced using environmentally-relevant exposure paradigm. In our previous study [28], exposure of pregnant rats to 0.2 mg/kg body weight of BDE-47 resulted in 234.3 ng BDE-47/g lipid in adipose tissue of dams and 1054.7 ng BDE-47/g lipid in adipose tissue of pups. These concentrations are comparable with those of the North American human population (mean concentration in adipose tissue of adult subjects = 399 ng/g lipids) [29]. Given that the rate of BDE-47 elimination is around 10 times higher in mice than in rats [30,31], the exposure used in this study is also relevant to that of the North American human population. In this study we dosed animals during the pre- or neonatal period since human PBDE exposure peaks during the perinatal period of life due to the active transport of PBDE via cord blood and breast milk [32,33,34], higher rates of dust ingestion [35], and higher food intake per kilogram of body weight in toddlers [33].

4.2. Downregulation of the Complement and Coagulation Cascades Genes

Coagulation is the primary defense against blood loss at a site of injury. It is characterized by a serine protease mediated conversion of zymogens to their active form to induce blood clotting [36]. The major players in the coagulation and fibrinolytic system are circulating zymogens (fibrinogen, clotting factors, protein C and its cofactor protein S, plasminogen, kininogen), activating and converting enzymes (plasminogen and thrombin), a group of serine protease inhibitors (SERPINs), and high density lipoproteins (HDL). Our results indicate that the expression of the major players in the coagulation pathway and fibrinolytic pathway, namely fibrinogen protein chains (FGG, FGA, FGB), clotting factors (F2, F7, F10, F11, and F12), serine protease inhibitors (SERPINC1, SERPIND1, SERPINF2), plasminogen (PLG), and protein C (PROC) were downregulated by more than 2-fold. We also observed more than 2-fold downregulation of the complement cascade genes such as C9, C8A, C8B, C4bp Cfi, HC, Mbl1, and MBl2. Adipose tissue has an important role in direct regulation of the coagulation and complement system by producing coagulation tissue factors and complement proteins [37,38]. Therefore, the observed changes in gene expression may have significant consequences for the whole complement and coagulation system in exposed organisms.

4.3. Downregulation of Xenobiotic Metabolism Genes

A big group of heterogeneous xenobiotic metabolism genes, including but not limited to cytochromes P450 (Cyp1a, Cyp2c, and Cyp3a subfamilies), was downregulated in exposed mice. This group of genes is usually induced in toxicological experiments, including experiments with BDE-47 [39,40,41,42,43], while in our experiment enrichment was negative, corresponding to a significant decrease in the expression of many genes. We speculate that this may be a result of overcompensation for sub-chronic developmental exposure to a low dose of BDE-47. Similarly, changes in gene expression, opposite in direction to predicted, were reported recently for liver tissue in BDE-47-exposed mice [44]. Altered expression of P450 cytochrome genes may affect the availability of active forms of vitamin D [45], which was recently shown to regulate the adipogenic fate of adipose-derived stem cells [46].

4.4. Upregulation of De Novo Lipogenesis Genes

Our bioinformatic analysis suggests that the metabolic profile of adipose tissue in exposed animals is slightly shifted in the direction of activation of de novo lipogenesis (DNL). For example positive enrichment of GSEA ‘SREBF targets’ dataset was due to upregulation of several essential DNL genes, including 1.6-fold induction of Acaca, encoding DNL rate-limiting enzyme acetyl-CoA carboxylase [47]. Genes of cholesterol biosynthesis were also upregulated, including 1.5-fold induction of Hmgcs1, encoding cytosolic 3-hydroxy-3-methylglutaryl-CoA (HMG-CoA) synthase, a rate-limiting enzyme of cholesterol biosynthesis [48]. At the same time, bile acid metabolism genes were downregulated, suggesting decreased cleavage of cholesterol, concordant with the previously reported increase in circulating cholesterol in rats perinatally exposed to BDE-47 [39]. However, it is difficult to interpret the observed changes in molecular pathways responsible for cholesterol turnover in the current study, as adipose tissue is not a major player in cholesterol biosynthesis and elimination. Importantly, changes in gene expression discussed in this paragraph may affect the major function of the adipose tissue—to store energy in the form of lipids. This concern is supported by human population data reporting a positive association between PBDE concentrations in blood, milk, and placenta and BMI in nursing women [34,49,50,51]. In several animal studies, developmental BDE-47 exposures also resulted in increased body weight [52,53,54].

4.5. Network Affected by BDE-47

The coagulation and complement system is known to be regulated among other transcription factors by HNF1, HNF4α and PPARα—mechanisms enriched in our bioinformatic analysis. For example the α and β promoters of fibrinogen gene contain tissue specific motif for HNF1 [55]. SERPINC1 (antithrombin) is transcriptionally activated by the binding of HNF1α and HNF4α to SERPINC1 gene (α1AT) promoter/enhancer [56]. Additionally, synthetic PPARα agonists, fibrates, have been shown to suppress fibrinogen gene expression and impair the induction of tissue factors in endothelial cells in rodents [57,58] and downregulate the complement and coagulation pathways in cynomolgus monkey liver [59]. Both transcription factor HNF4A and PPARα regulate expression of major protein component of the high density lipoproteins (HDL) [60,61]. In mice (contrary to humans), PPARα upregulates expression of APO5, APOA1, and APOA2 [62,63,64]. In our study, the expression of several HDL genes APOA1, APOA2, APOA5, APOB, APOC2, APOC3, APOC4, APOH, APOM was more than 2-fold downregulated. HDL modulates protease activity by directly binding to the majority of SERPINs [65]. Finally the expression of PPARα is regulated by HNF4 in mice, and PPARα is capable of binding to HNF4αRE [66].

4.6. Role of PPARα in the Modulation of BDE-47 Induced Response

Our analysis identifies PPARα as a hub of a molecular network affected in adipose tissue by BDE-47. No studies reported previously BDE-47 interaction with PPARα. In fact, a theoretical structure–activity assessment predicted that PBDEs would not have an interaction with PPAR receptors [67]. However, several studies have indicated the ability of BDE-47 and/or its metabolites to activate another PPAR isoform—PPARγ [68,69,70]. This evidence, taken together with the fact that the PPAR nuclear family is considered to have the largest and the most promiscuous binding pocket [71], suggests the possibility that most of the observed changes in gene expression in the current study are a result of direct interaction of BDE-47 and/or its metabolites with PPARα. An alternative hypothesis may explain the changes in PPARα activity by the ability of BDE-47 to interact with regulatory mechanisms upstream of the receptor. In our recent study we have shown that the same dosing paradigm used in this study results in activation of mTOR signaling in mouse livers on PND21 [44]. PPARα is a well-known target of mTOR complex 1 (mTORC1) [72,73,74]. Additional studies are needed to determine the interaction between BDE-47 and the PPAR-regulated molecular network.

5. Conclusions

We have demonstrated that low-dose perinatal exposure to BDE-47 produces significant changes in gene expression in the gonadal adipose tissue in mice, suggesting that adipose tissue should be considered a target tissue for BDE-47, in addition to its role as a storage depot. Our data imply that functional networks, including coagulation and complement cascades, de novo lipogenesis, and xenobiotic metabolism, are affected by BDE-47. The affected molecular network is connected via the following hubs: PPARα, HNF1A, and HNF4. Altered metabolism of vitamin D due to the changes in P450 gene expression may interfere with adipose tissue function. This study builds a background for a targeted search of sensitive phenotypic endpoints of BDE-47 exposure, including lipid profile parameters and coagulation factors in circulation. Additional research is needed to dissect molecular cascades affected in the adipose tissue by PBDEs, and the role of PBDEs as obesogens.

Supplementary Materials

The following are available online at https://www.mdpi.com/2305-6304/6/2/21/s1, Table S1: Significant 2-fold differentially expressed genes in gonadal adipose tissue of mice exposed perinatally to 0.2 mg/kg body weight BDE-47; Table S2: GSEA analysis: summary of top enriched gene sets in the gonadal adipose tissue of mice exposed perinatally to 0.2 mg/kg body weight BDE-47.

Acknowledgments

This study was supported by startup funding from the University of Massachusetts, School of Public Health and Health Sciences to A.S. Library preparation for this study was done by Supei Ma, MPH.

Author Contributions

A.S. conceived and designed the experiments, carried out all experimental work with animals, and coordinated laboratory procedures, bioinformatic analysis, and manuscript preparation. A.A. conducted bioinformatic analysis of RNA-seq data and drafted the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gregg, E.W.; Zhuo, X.; Cheng, Y.J.; Albright, A.L.; Narayan, K.M.; Thompson, T.J. Trends in lifetime risk and years of life lost due to diabetes in the USA, 1985–2011: A modelling study. Lancet Diabetes Endocrinol. 2014, 2, 867–874. [Google Scholar] [CrossRef]
  2. National Institutes of Health (NIH). Funding Opportunity Announcement: Role of Environmental Chemical Exposures in the Development of Obesity, Type 2 Diabetes and Metabolic Syndrome (R21); National Institutes of Health: Bethesda, MD, USA, 2012.
  3. Heindel, J.J.; Vom Saal, F.S.; Blumberg, B.; Bovolin, P.; Calamandrei, G.; Ceresini, G.; Cohn, B.A.; Fabbri, E.; Gioiosa, L.; Kassotis, C.; et al. Parma consensus statement on metabolic disruptors. Environ. Health 2015, 14. [Google Scholar] [CrossRef] [PubMed]
  4. La Merrill, M.; Emond, C.; Kim, M.J.; Antignac, J.P.; Le Bizec, B.; Clément, K.; Birnbaum, L.S.; Barouki, R. Toxicological function of adipose tissue: Focus on persistent organic pollutants. Environ. Health Perspect. 2013, 121, 162–169. [Google Scholar] [PubMed]
  5. Trayhurn, P. Endocrine and signalling role of adipose tissue: New perspectives on fat. Acta Physiol. Scand. 2005, 184, 285–293. [Google Scholar] [CrossRef] [PubMed]
  6. Kelishadi, R.; Poursafa, P.; Jamshidi, F. Role of environmental chemicals in obesity: A systematic review on the current evidence. J. Environ. Public Health 2013, 2013, 896789. [Google Scholar] [CrossRef] [PubMed]
  7. Regnier, S.M.; Sargis, R.M. Adipocytes under assault: Environmental disruption of adipose physiology. Biochim. Biophys. Acta 2014, 1842, 520–533. [Google Scholar] [CrossRef] [PubMed]
  8. Suvorov, A.; Takser, L. Facing the challenge of data transfer from animal models to humans: The case of persistent organohalogens. Environ. Health 2008. [Google Scholar] [CrossRef] [PubMed]
  9. Hurley, S.; Goldberg, D.; Nelson, D.O.; Guo, W.; Wang, Y.; Baek, H.G.; Park, J.S.; Petreas, M.; Bernstein, L.; Anton-Culver, H.; et al. Temporal Evaluation of Polybrominated Diphenyl Ether (PBDE) Serum Levels in Middle-Aged and Older California Women, 2011–2015. Environ. Sci. Technol. 2017, 51, 4697–4704. [Google Scholar] [CrossRef] [PubMed]
  10. Herbstman, J.B.; Sjödin, A.; Apelberg, B.J.; Witter, F.R.; Patterson, D.G., Jr.; Halden, R.U.; Jones, R.S.; Park, A.; Zhang, Y.; Heidler, J.; et al. Determinants of prenatal exposure to polychlorinated biphenyls (PCBs) and polybrominated diphenyl ethers (PBDEs) in an urban population. Environ. Health Perspect. 2007, 115, 1794–1800. [Google Scholar] [CrossRef] [PubMed]
  11. Doucet, J.; Tague, B.; Arnold, D.L.; Cooke, G.M.; Hayward, S.; Goodyer, C.G. Persistent organic pollutant residues in human fetal liver and placenta from Greater Montreal, Quebec: A longitudinal study from 1998 through 2006. Environ. Health Perspect. 2009, 117, 605–610. [Google Scholar] [CrossRef] [PubMed]
  12. Guvenius, D.M.; Aronsson, A.; Ekman-Ordeberg, G.; Bergman, A.; Noren, K. Human prenatal and postnatal exposure to polybrominated diphenyl ethers, polychlorinated biphenyls, polychlorobiphenylols, and pentachlorophenol. Environ. Health Perspect. 2003, 111, 1235–1241. [Google Scholar] [CrossRef] [PubMed]
  13. Marchitti, S.A.; Lakind, J.S.; Naiman, D.Q.; Berlin, C.M.; Kenneke, J.F. Improving infant exposure and health risk estimates: Using serum data to predict polybrominated diphenyl ether concentrations in breast milk. Environ. Sci. Technol. 2013, 47, 4787–4795. [Google Scholar] [CrossRef] [PubMed]
  14. Schecter, A.; Pavuk, M.; Papke, O.; Ryan, J.J.; Birnbaum, L.; Rosen, R. Polybrominated diphenyl ethers (PBDEs) in U.S. mothers’ milk. Environ. Health Perspect. 2003, 111, 1723–1729. [Google Scholar] [CrossRef] [PubMed]
  15. Horton, M.K.; Bousleiman, S.; Jones, R.; Sjodin, A.; Liu, X.; Whyatt, R.; Wapner, R.; Factor-Litvak, P. Predictors of serum concentrations of polybrominated flame retardants among healthy pregnant women in an urban environment: A cross-sectional study. Environ. Health 2013, 12. [Google Scholar] [CrossRef] [PubMed]
  16. Geyer, H.J.; Schramm, K.W.; Darnerud, P.O.; Aune, M.; Feicht, E.A.; Fried, K.W.; Henkelmann, B.; Lenoir, D.; Schmid, P.; McDonald, T.A. Terminal Elimination Half-lives of the Brominated Flame Retardants TBBPA, HBCD, and Lower Brominated PBDEs in Humans. Organohalogen Compd. 2004, 66, 3820–3825. [Google Scholar]
  17. Renwick, A.G.; Walton, K. The use of surrogate endpoints to assess potential toxicity in humans. Toxicol. Lett. 2001, 120, 97–110. [Google Scholar] [CrossRef]
  18. Wade, M.G.; Parent, S.; Finnson, K.W.; Foster, W.; Younglai, E.; McMahon, A.; Cyr, D.G.; Hughes, C. Thyroid toxicity due to subchronic exposure to a complex mixture of 16 organochlorines, lead, and cadmium. Toxicol. Sci. 2002, 67, 207–218. [Google Scholar] [CrossRef] [PubMed]
  19. Marsh, G.; Bergman, A.; Bladh, L.; Gillner, M.; Jakobsson, E. Synthesis of p-hydrobromobiphenyl ethers and binding to the thyroid hormone receptor. Organohalogen Compd. 1998, 37, 305–308. [Google Scholar]
  20. Vandenberg, L.N.; Welshons, W.V.; Vom Saal, F.S.; Toutain, P.L.; Myers, J.P. Should oral gavage be abandoned in toxicity testing of endocrine disruptors? Environ. Health 2014, 13. [Google Scholar] [CrossRef] [PubMed]
  21. Suvorov, A.; Vandenberg, L.N. To Cull or not to Cull? Considerations for studies of endocrine disrupting chemicals. Endocrinology 2016. [Google Scholar] [CrossRef] [PubMed]
  22. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef] [PubMed]
  23. Kim, D.; Pertea, G.; Trapnell, C.; Pimentel, H.; Kelley, R.; Salzberg, S.L. TopHat2: Accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013, 14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Trapnell, C.; Pachter, L.; Salzberg, S.L. TopHat: Discovering splice junctions with RNA-Seq. Bioinformatics 2009, 25, 1105–1111. [Google Scholar] [CrossRef] [PubMed]
  25. Trapnell, C.; Hendrickson, D.G.; Sauvageau, M.; Goff, L.; Rinn, J.L.; Pachter, L. Differential analysis of gene regulation at transcript resolution with RNA-seq. Nat. Biotechnol. 2013, 31, 46–53. [Google Scholar] [CrossRef] [PubMed]
  26. Subramanian, A.; Tamayo, P.; Mootha, V.K.; Mukherjee, S.; Ebert, B.L.; Gillette, M.A.; Paulovich, A.; Pomeroy, S.L.; Golub, T.R.; Lander, E.S.; et al. Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl. Acad. Sci. USA 2005, 102, 15545–15550. [Google Scholar] [CrossRef] [PubMed]
  27. Rakhshandehroo, M.; Knoch, B.; Muller, M.; Kersten, S. Peroxisome proliferator-activated receptor alpha target genes. PPAR Res. 2010, 2010. [Google Scholar] [CrossRef] [PubMed]
  28. Suvorov, A.; Girard, S.; Lachapelle, S.; Abdelouahab, N.; Sebire, G.; Takser, L. Perinatal exposure to low-dose BDE-47, an emergent environmental contaminant, causes hyperactivity in rat offspring. Neonatology 2009, 95, 203–209. [Google Scholar] [CrossRef] [PubMed]
  29. Johnson-Restrepo, B.; Kannan, K.; Rapaport, D.P.; Rodan, B.D. Polybrominated diphenyl ethers and polychlorinated biphenyls in human adipose tissue from New York. Environ. Sci. Technol. 2005, 39, 5177–5182. [Google Scholar] [CrossRef] [PubMed]
  30. Sanders, J.M.; Chen, L.J.; Lebetkin, E.H.; Burka, L.T. Metabolism and disposition of 2,2′,4,4′-tetrabromodiphenyl ether following administration of single or multiple doses to rats and mice. Xenobiotica 2006, 36, 103–117. [Google Scholar] [CrossRef] [PubMed]
  31. Staskal, D.F.; Diliberto, J.J.; DeVito, M.J.; Birnbaum, L.S. Toxicokinetics of BDE 47 in female mice: Effect of dose, route of exposure, and time. Toxicol. Sci. 2005, 83, 215–223. [Google Scholar] [CrossRef] [PubMed]
  32. Antignac, J.P.; Cariou, R.; Maume, D.; Marchand, P.; Monteau, F.; Zalko, D.; Berrebi, A.; Cravedi, J.P.; Andre, F.; Le Bizec, B. Exposure assessment of fetus and newborn to brominated flame retardants in France: Preliminary data. Mol. Nutr. Food Res. 2008, 52, 258–265. [Google Scholar] [CrossRef] [PubMed]
  33. Schecter, A.; Päpke, O.; Harris, T.R.; Tung, K.C.; Musumba, A.; Olson, J.; Birnbaum, L. Polybrominated diphenyl ether (PBDE) levels in an expanded market basket survey of U.S. food and estimated PBDE dietary intake by age and sex. Environ. Health Perspect. 2006, 114, 1515–1520. [Google Scholar] [CrossRef] [PubMed]
  34. Shi, Z.; Jiao, Y.; Hu, Y.; Sun, Z.; Zhou, X.; Feng, J.; Li, J.; Wu, Y. Levels of tetrabromobisphenol A, hexabromocyclododecanes and polybrominated diphenyl ethers in human milk from the general population in Beijing, China. Sci. Total Environ. 2013, 452–453, 10–18. [Google Scholar] [CrossRef] [PubMed]
  35. Wilford, B.H.; Shoeib, M.; Harner, T.; Zhu, J.; Jones, K.C. Polybrominated diphenyl ethers in indoor dust in Ottawa, Canada: Implications for sources and exposure. Environ. Sci. Technol. 2005, 39, 7027–7035. [Google Scholar] [CrossRef] [PubMed]
  36. Hertle, E.; Stehouwer, C.D.; van Greevenbroek, M.M. The complement system in human cardiometabolic disease. Mol. Immunol. 2014, 61, 135–148. [Google Scholar] [CrossRef] [PubMed]
  37. Samad, F.; Ruf, W. Inflammation, obesity, and thrombosis. Blood 2013, 122, 3415–3422. [Google Scholar] [CrossRef] [PubMed]
  38. Samad, F.; Pandey, M.; Loskutoff, D.J. Regulation of tissue factor gene expression in obesity. Blood 2001, 98, 3353–3358. [Google Scholar] [CrossRef] [PubMed]
  39. Suvorov, A.; Takser, L. Global Gene Expression Analysis in the Livers of Rat Offspring Perinatally Exposed to Low Doses of 2,2′,4,4′-tetrabromodiphenyl ether. Environ. Health Perspect. 2010, 118, 97–102. [Google Scholar] [CrossRef] [PubMed]
  40. Sanders, J.M.; Burka, L.T.; Smith, C.S.; Black, W.; James, R.; Cunningham, M.L. Differential expression of CYP1A, 2B, and 3A genes in the F344 rat following exposure to a polybrominated diphenyl ether mixture or individual components. Toxicol. Sci. 2005, 88, 127–133. [Google Scholar] [CrossRef] [PubMed]
  41. Wahl, M.; Lahni, B.; Guenther, R.; Kuch, B.; Yang, L.; Straehle, U.; Strack, S.; Weiss, C. A technical mixture of 2,2′,4,4′-tetrabromo diphenyl ether (BDE47) and brominated furans triggers aryl hydrocarbon receptor (AhR) mediated gene expression and toxicity. Chemosphere 2008, 73, 209–215. [Google Scholar] [CrossRef] [PubMed]
  42. Pacyniak, E.K.; Cheng, X.; Cunningham, M.L.; Crofton, K.; Klaassen, C.D.; Guo, G.L. The flame retardants, polybrominated diphenyl ethers, are pregnane X receptor activators. Toxicol. Sci. 2007, 97, 94–102. [Google Scholar] [CrossRef] [PubMed]
  43. Dunnick, J.K.; Brix, A.; Cunny, H.; Vallant, M.; Shockley, K.R. Characterization of polybrominated diphenyl ether toxicity in Wistar Han rats and use of liver microarray data for predicting disease susceptibilities. Toxicol. Pathol. 2012, 40, 93–106. [Google Scholar] [CrossRef] [PubMed]
  44. Khalil, A.; Parker, M.; Mpanga, R.; Cevik, E.C.; Thorburn, C.; Suvorov, A. Developmental Exposure to 2,2′,4,4′–Tetrabromodiphenyl Ether Induces Long-Lasting Changes in Liver Metabolism in Male Mice. J. Endocr. Soc. 2017, 1, 323–344. [Google Scholar] [CrossRef] [PubMed]
  45. Wang, Z.; Lin, Y.S.; Zheng, X.E.; Senn, T.; Hashizume, T.; Scian, M.; Dickmann, L.J.; Nelson, S.D.; Baillie, T.A.; Hebert, M.F.; et al. An inducible cytochrome P450 3A4-dependent vitamin D catabolic pathway. Mol. Pharmacol. 2012, 81, 498–509. [Google Scholar] [CrossRef] [PubMed]
  46. Basoli, V.; Santaniello, S.; Cruciani, S.; Ginesu, G.C.; Cossu, M.L.; Delitala, A.P.; Serra, P.A.; Ventura, C.; Maioli, M. Melatonin and Vitamin D Interfere with the Adipogenic Fate of Adipose-Derived Stem Cells. Int. J. Mol. Sci. 2017, 18, 981. [Google Scholar] [CrossRef] [PubMed]
  47. Mashek, D.G. Hepatic fatty acid trafficking: Multiple forks in the road. Adv. Nutr. 2013, 4, 697–710. [Google Scholar] [CrossRef] [PubMed]
  48. Nguyen, P.; Leray, V.; Diez, M.; Serisier, S.; Bloc’h, J.L.; Siliart, B.; Dumon, H. Liver lipid metabolism. J. Anim. Physiol. Anim. Nutr. 2008, 92, 272–283. [Google Scholar] [CrossRef] [PubMed]
  49. Daniels, J.L.; Pan, I.J.; Jones, R.; Anderson, S.; Patterson, D.G., Jr.; Needham, L.L.; Sjödin, A. Individual characteristics associated with PBDE levels in U.S. human milk samples. Environ. Health Perspect. 2010, 118, 155–160. [Google Scholar] [PubMed]
  50. Croes, K.; Colles, A.; Koppen, G.; Govarts, E.; Bruckers, L.; Van de Mieroop, E.; Nelen, V.; Covaci, A.; Dirtu, A.C.; Thomsen, C.; et al. Persistent organic pollutants (POPs) in human milk: A biomonitoring study in rural areas of Flanders (Belgium). Chemosphere 2012, 89, 988–994. [Google Scholar] [CrossRef] [PubMed]
  51. Frederiksen, M.; Thomsen, M.; Vorkamp, K.; Knudsen, L.E. Patterns and concentration levels of polybrominated diphenyl ethers (PBDEs) in placental tissue of women in Denmark. Chemosphere 2009, 76, 1464–1469. [Google Scholar] [CrossRef] [PubMed]
  52. Gee, J.R.; Moser, V.C. Acute postnatal exposure to brominated diphenylether 47 delays neuromotor ontogeny and alters motor activity in mice. Neurotoxicol. Teratol. 2008, 30, 79–87. [Google Scholar] [CrossRef] [PubMed]
  53. Bondy, G.S.; Lefebvre, D.E.; Aziz, S.; Cherry, W.; Coady, L.; MacLellan, E.; Armstrong, C.; Barker, M.; Cooke, G.; Gaertner, D.; et al. Toxicologic and immunologic effects of perinatal exposure to the brominated diphenyl ether (BDE) mixture DE-71 in the Sprague-Dawley rat. Environ. Toxicol. 2013, 28, 215–228. [Google Scholar] [CrossRef] [PubMed]
  54. Suvorov, A.; Battista, M.C.; Takser, L. Perinatal exposure to low-dose 2,2′,4,4′-tetrabromodiphenyl ether affects growth in rat offspring: What is the role of IGF-1? Toxicology 2009, 260, 126–131. [Google Scholar] [CrossRef] [PubMed]
  55. Fuller, G.M.; Zhang, Z. Transcriptional control mechanism of fibrinogen gene expression. Ann. N. Y. Acad. Sci. 2001, 936, 469–479. [Google Scholar] [CrossRef] [PubMed]
  56. Zhao, H.; Friedman, R.D.; Fournier, R.E. The locus control region activates serpin gene expression through recruitment of liver-specific transcription factors and RNA polymerase II. Mol. Cell. Biol. 2007, 27, 5286–5295. [Google Scholar] [CrossRef] [PubMed]
  57. Golledge, J.; Mangan, S.; Clancy, P. Effects of peroxisome proliferator-activated receptor ligands in modulating tissue factor and tissue factor pathway inhibitor in acutely symptomatic carotid atheromas. Stroke 2007, 38, 1501–1508. [Google Scholar] [CrossRef] [PubMed]
  58. Kockx, M.; Gervois, P.P.; Poulain, P.; Derudas, B.; Peters, J.M.; Gonzalez, F.J.; Princen, H.M.; Kooistra, T.; Staels, B. Fibrates suppress fibrinogen gene expression in rodents via activation of the peroxisome proliferator-activated receptor-alpha. Blood 1999, 93, 2991–2998. [Google Scholar] [PubMed]
  59. Cariello, N.F.; Romach, E.H.; Colton, H.M.; Ni, H.; Yoon, L.; Falls, J.G.; Casey, W.; Creech, D.; Anderson, S.P.; Benavides, G.R.; et al. Gene expression profiling of the PPAR-alpha agonist ciprofibrate in the cynomolgus monkey liver. Toxicol. Sci. 2005, 88, 250–264. [Google Scholar] [CrossRef] [PubMed]
  60. Mogilenko, D.A.; Dizhe, E.B.; Shavva, V.S.; Lapikov, I.A.; Orlov, S.V.; Perevozchikov, A.P. Role of the nuclear receptors HNF4 alpha, PPAR alpha, and LXRs in the TNF alpha-mediated inhibition of human apolipoprotein A-I gene expression in HepG2 cells. Biochemistry 2009, 48, 11950–11960. [Google Scholar] [CrossRef] [PubMed]
  61. Malik, S. Transcriptional regulation of the apolipoprotein AI gene. Front. Biosci. 2003, 8, d360–d368. [Google Scholar] [CrossRef] [PubMed]
  62. Berthou, L.; Duverger, N.; Emmanuel, F.; Langouët, S.; Auwerx, J.; Guillouzo, A.; Fruchart, J.C.; Rubin, E.; Denèfle, P.; Staels, B.; et al. Opposite regulation of human versus mouse apolipoprotein A-I by fibrates in human apolipoprotein A-I transgenic mice. J. Clin. Investig. 1996, 97, 2408–2416. [Google Scholar] [CrossRef] [PubMed]
  63. Prieur, X.; Lesnik, P.; Moreau, M.; Rodríguez, J.C.; Doucet, C.; Chapman, M.J.; Huby, T. Differential regulation of the human versus the mouse apolipoprotein AV gene by PPARalpha. Implications for the study of pharmaceutical modifiers of hypertriglyceridemia in mice. Biochim. Biophys. Acta 2009, 1791, 764–771. [Google Scholar] [CrossRef] [PubMed]
  64. Pawlak, M.; Lefebvre, P.; Staels, B. Molecular mechanism of PPARalpha action and its impact on lipid metabolism, inflammation and fibrosis in non-alcoholic fatty liver disease. J. Hepatol. 2015, 62, 720–733. [Google Scholar] [CrossRef] [PubMed]
  65. Gordon, S.M.; Remaley, A.T. High density lipoproteins are modulators of protease activity: Implications in inflammation, complement activation, and atherothrombosis. Atherosclerosis 2017, 259, 104–113. [Google Scholar] [CrossRef] [PubMed]
  66. Pineda Torra, I.; Jamshidi, Y.; Flavell, D.M.; Fruchart, J.C.; Staels, B. Characterization of the human PPARalpha promoter: Identification of a functional nuclear receptor response element. Mol. Endocrinol. 2002, 16, 1013–1028. [Google Scholar] [PubMed]
  67. Luthe, G.; Jacobus, J.A.; Robertson, L.W. Receptor interactions by polybrominated diphenyl ethers versus polychlorinated biphenyls: A theoretical Structure-activity assessment. Environ. Toxicol. Pharmacol. 2008, 25, 202–210. [Google Scholar] [CrossRef] [PubMed]
  68. Tung, E.W.; Boudreau, A.; Wade, M.G.; Atlas, E. Induction of adipocyte differentiation by polybrominated diphenyl ethers (PBDEs) in 3T3-L1 cells. PLoS ONE 2014, 9, e94583. [Google Scholar] [CrossRef] [PubMed]
  69. Kamstra, J.H.; Hruba, E.; Blumberg, B.; Janesick, A.; Mandrup, S.; Hamers, T.; Legler, J. Transcriptional and epigenetic mechanisms underlying enhanced in vitro adipocyte differentiation by the brominated flame retardant BDE-47. Environ. Sci. Technol. 2014, 48, 4110–4119. [Google Scholar] [CrossRef] [PubMed]
  70. Fang, M.; Webster, T.F.; Ferguson, P.L.; Stapleton, H.M. Characterizing the peroxisome proliferator-activated receptor (PPARgamma) ligand binding potential of several major flame retardants, their metabolites, and chemical mixtures in house dust. Environ. Health Perspect. 2015, 123, 166–172. [Google Scholar] [PubMed]
  71. Grun, F.; Blumberg, B. Endocrine disrupters as obesogens. Mol. Cell. Endocrinol. 2009, 304, 19–29. [Google Scholar] [CrossRef] [PubMed]
  72. Laplante, M.; Sabatini, D.M. Regulation of mTORC1 and its impact on gene expression at a glance. J. Cell Sci. 2013, 126, 1713–1719. [Google Scholar] [CrossRef] [PubMed]
  73. Laplante, M.; Sabatini, D.M. mTOR signaling in growth control and disease. Cell 2012, 149, 274–293. [Google Scholar] [CrossRef] [PubMed]
  74. Sengupta, S.; Peterson, T.R.; Laplante, M.; Oh, S.; Sabatini, D.M. mTORC1 controls fasting-induced ketogenesis and its modulation by ageing. Nature 2010, 468, 1100–1104. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Heat-map for 758 2-fold differentially regulated genes in mouse gonadal adipose tissue on PND21.
Figure 1. Heat-map for 758 2-fold differentially regulated genes in mouse gonadal adipose tissue on PND21.
Toxics 06 00021 g001
Figure 2. KEGG coagulation and complement cascades pathway. Differentially expressed genes are marked with red stars.
Figure 2. KEGG coagulation and complement cascades pathway. Differentially expressed genes are marked with red stars.
Toxics 06 00021 g002
Figure 3. GSEA enrichment plots for transcriptional changes in gonadal adipose tissue of mice exposed perinatally to BDE-47: (A) bile acid metabolism gene set; (B) genes targets of HNF4A; and (C) genes targets of SREBF. GSEA plot legend: 1—running enrichment score for the gene set; 2—vertical lines show where the members of the gene set appear in the ranked list of genes; 3 and 4—ranked list of differentially expressed genes from the most up-regulated (left of each plot) to the most downregulated (right of each plot). Regulation is shown in the heat-bar (3) and in the bar plot (4).
Figure 3. GSEA enrichment plots for transcriptional changes in gonadal adipose tissue of mice exposed perinatally to BDE-47: (A) bile acid metabolism gene set; (B) genes targets of HNF4A; and (C) genes targets of SREBF. GSEA plot legend: 1—running enrichment score for the gene set; 2—vertical lines show where the members of the gene set appear in the ranked list of genes; 3 and 4—ranked list of differentially expressed genes from the most up-regulated (left of each plot) to the most downregulated (right of each plot). Regulation is shown in the heat-bar (3) and in the bar plot (4).
Toxics 06 00021 g003
Figure 4. Molecular networks of differentially expressed genes regulated by (A) HNF1A and (B) PPARα.
Figure 4. Molecular networks of differentially expressed genes regulated by (A) HNF1A and (B) PPARα.
Toxics 06 00021 g004aToxics 06 00021 g004b
Table 1. RT-qPCR validation of RNA-seq data for select differentially expressed genes.
Table 1. RT-qPCR validation of RNA-seq data for select differentially expressed genes.
GeneqPCR PrimersFold Change
RNA-seqqPCR
Grem2GAGGAGAGGGACAGGGAGAC2.82.6
AGCGAGAGCTTCCAGAACAT
CishCGTTGTCTCTGGGACATGGT4.43.9
ACAGCCAGCAAAGGACAAGA
Sfrp2AGGACAACGACCTCTGCATC4.34.2
TTCTTGGTTTTGCAGGCTTC
FggAACCCACGAGACAAGCATTC−15.7−16.3
CAAGCTGGGCTACCTTCTGT
F2CCCTGGTTAGCCTGGTACAC−16.2−13.9
CAGCGATAGTGCTTGCTGAG
ApobAAAACTTAAGCTTCAGCTGTCCA−20.3−17.5
CTTTCGGAGGTGCTTGAATC
B2mCCGGCCTGTATGCTATCCAG--housekeeping
TGTTCGGCTTCCCATTCTCC
Table 2. DAVID functional annotation clustering for 94 2-fold differentially expressed genes (FDR q ≤ 0.05).
Table 2. DAVID functional annotation clustering for 94 2-fold differentially expressed genes (FDR q ≤ 0.05).
Functional Annotation ClusterEnrichment ScoreGenesp-ValueBenjamini
Secreted proteins11.17386.7 × 10−6–3.0 × 10−187.0 × 10−4–5.9 × 10−16
Peptidase inhibitor activity5.15102.7 × 10−2–1.9 × 10−96.2 × 10−1–9.2 × 10−8
Blood coagulation3.14222.7 × 10−2–1.9 × 10−96.2 × 10−1–9.2 × 10−8
Biosynthesis of amino acids2.8164.6 × 10−2–1.2 × 10−43.0 × 10−1–5.7 × 10−3
Microsomal metabolism2.4071.1 × 10−1–1.9 × 10−55.3 × 10−1–4.5 × 10−4
Lipid metabolism2.2753.8 × 10−2–7.8 × 10−43.9 × 10−1–3.0 × 10−2
Table 3. Ingenuity Pathway Analysis: top enriched canonical pathways.
Table 3. Ingenuity Pathway Analysis: top enriched canonical pathways.
Canonical Pathwayp-ValueOverlapCount
FXR/RXR Activation2.4 × 102325.4%32/126
LXR/RXR Activation2.99 × 101721.5%26/121
Acute Phase Response Signaling3.02 × 101717.8%30/169
Coagulation System1.71 × 101542.9%15/35
Extrinsic Prothrombin Activation Pathway3.68 × 101156.2%9/16

Share and Cite

MDPI and ACS Style

Abrha, A.; Suvorov, A. Transcriptomic Analysis of Gonadal Adipose Tissue in Male Mice Exposed Perinatally to 2,2′,4,4′-Tetrabromodiphenyl Ether (BDE-47). Toxics 2018, 6, 21. https://doi.org/10.3390/toxics6020021

AMA Style

Abrha A, Suvorov A. Transcriptomic Analysis of Gonadal Adipose Tissue in Male Mice Exposed Perinatally to 2,2′,4,4′-Tetrabromodiphenyl Ether (BDE-47). Toxics. 2018; 6(2):21. https://doi.org/10.3390/toxics6020021

Chicago/Turabian Style

Abrha, Aser, and Alexander Suvorov. 2018. "Transcriptomic Analysis of Gonadal Adipose Tissue in Male Mice Exposed Perinatally to 2,2′,4,4′-Tetrabromodiphenyl Ether (BDE-47)" Toxics 6, no. 2: 21. https://doi.org/10.3390/toxics6020021

APA Style

Abrha, A., & Suvorov, A. (2018). Transcriptomic Analysis of Gonadal Adipose Tissue in Male Mice Exposed Perinatally to 2,2′,4,4′-Tetrabromodiphenyl Ether (BDE-47). Toxics, 6(2), 21. https://doi.org/10.3390/toxics6020021

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop