Sex- and Dose-Specific Effects of Maternal Bisphenol A Exposure on Pancreatic Islets of First- and Second-Generation Adult Mice Offspring.
Over the past decade, the incidence of diabetes has increased dramatically. The WHO estimates mortality due to diabetes to increase by 50% in the next 10 y, making it a major cause of death by 2030 and placing a substantial economic burden on the health care system globally (Mathers and Loncar 2006; American Diabetes Association 2013). Therefore, understanding the underlying factors of this global epidemic is critical. It has been well established that factors such as food intake and sedentary lifestyle significantly influence the incidence of diabetes. Although these are important factors, changes in our lifestyle have also increased our exposure to a variety of synthetic chemicals that are linked to various diseases such as diabetes and obesity in experimental models and human epidemiological studies (Alonso-Magdalena et al. 2015b; Jasarevic et al. 2011; Lang et al. 2008; Newbold et al. 2008).
Bisphenol A (BPA) is one such synthetic chemical to which we are ubiquitously exposed through food, drink, and skin contact (Stahlhut et al. 2009). Detectable amounts of BPA are found in urine samples of populations from around the world (Nahar et al. 2012; Stahlhut et al. 2009; Zhang et al. 2011). Exposure to BPA has been associated with prediabetes and type 2 diabetes in human cross sectional (Aekplakorn et al. 2015; Ahmadkhaniha et al. 2014; Sabanayagam et al. 2013; Shankar and Teppala 2011) and longitudinal studies (Sun et al. 2014). Such associations have been further confirmed using animal models of exposure during adulthood (Alonso-Magdalena et al. 2006; Batista et al. 2012).
The pancreas undergoes substantial remodeling during late gestation and early neonatal life. This involves formation of new [beta] cells by differentiation of precursor cells during late gestation (Hellerstrom et al. 1988) followed by substantial remodeling, including high rates of replication and cell death in the neonatal period (Finegood et al. 1995; Scaglia et al. 1995, 1997; Swenne 1983). The [beta]-cell replication rate gradually declines (Finegood et al. 1995; Swenne 1983), and post-weaning [beta]-cell turnover is regulated by [beta]-cell death and very slow rates of replication (Bouwens and Rooman 2005; Scaglia et al. 1995). An exposure to synthetic chemicals such as BPA during this vulnerable period is even more concerning because of the potential to not only alter normal fetal development, but also to impact later life health and disease. Furthermore, the fetus is particularly vulnerable to the effects of BPA because the degradation pathways are immature, leading to accumulation and high levels of BPA (Takahashi and Oishi 2000; Zalko et al. 2003). Consistent with this finding, our previous studies and studies by others have shown that perinatal exposure to BPA leads to metabolic defects in offspring (Alonso-Magdalena et al. 2010; Angle et al. 2013; Mackay et al. 2013; van Esterik et al. 2014; Wei et al. 2011) and in subsequent generation offspring (Li et al. 2014; Susiarjo et al. 2015). We recently demonstrated that maternal exposure via maternal transmission to relevant human exposure levels of BPA, 10 [micro]g/kg/d (LowerB) or 10 mg/kg/d (UpperB) prior to conception and through lactation is associated with higher body fat, impaired glucose tolerance, and reduced glucose-stimulated insulin secretion (GSIS) across two generations in adult male mice offspring (Susiarjo et al. 2015). Despite our expanding knowledge on the association of BPA with increased risk of diabetes, the mechanisms underlying this link remain to be elucidated.
In the current study, we tested the hypothesis that early developmental exposure to BPA alters key aspects of [beta]-cell function leading to impaired insulin secretion in the offspring and that these changes would persist across two generations.
Materials and Methods
The animal work was conducted with the approval of the University of Pennsylvania Institutional Animal Care and Use Committee. All animals were treated humanely and with regard for alleviation of suffering. The animals used in this study were generated using the same procedure as described previously (Susiarjo et al. 2015). (See Figure S1 for the experimental design.) C57BL/6 virgin female mice (F0) were randomly assigned to the following three diets from 2 wk prior to mating until weaning: 10 [micro]g/kg/d (LowerB), 10 mg/kg/d (UpperB) BPA, and 7% corn oil (Control). These doses were selected based on our previous study (Susiarjo et al. 2013), where we reported that the levels of unconjugated BPA in serum from pregnant mice were representative of human exposure levels (Susiarjo et al. 2013). The offspring after weaning, and the subsequent generations, were maintained on the control diet. All terminal tests were performed in 16- to 21-wk-old first- (F1) and second-generation (F2) adult offspring, other than cell viability assays that were performed in 8- to 10-wk-old offspring. The experiments were blinded, where possible (e.g., immunofluorescence staining and Luminex assay). To account for robustness of the studied outcomes and cohort consistency in our paradigm, the end points measured were from animals generated across multiple cohorts studied in two different animal care facilities, with a control group included in all cohorts (see Figure S2). For most assessments, n = 5-6 litters per group were used other than for extensive cytokine assays, RNA seq, islet perifusion ramps, and islet cytosolic calcium determination where mice from n = 3-4 litters per group were used. Metabolic studies including GSIS and [alpha]-ketoisocaproate ramp studies were performed in both male and female offspring. However, because no metabolic phenotype was observed in female offspring across two generations, all mechanistic studies were performed on males only. Females were only additionally assessed for transcriptomic differences via RNA seq to determine if there were additional pathways we did not anticipate that might be disrupted in the islet.
F1 and F2 adult mice pancreata were perfused with Hanks' balanced salt solution (Life Technologies) supplemented with 2.5% bovine serum albumin (BSA) (wt/vol; Sigma), 0.35 g/L NaHC[O.sub.3] (Sigma), and 2 mg/mL Collagenase P (Roche); excised; and incubated at 37[degrees]C for 10-15 min. After digestion, islets were washed and then purified using a Histopaque 10771 and Histopaque 11191 (Sigma) gradient. To provide sufficient starting material for Oroboros, RNA and DNA extraction, islets had to be pooled from 2-3 male mice (either male or female for RNA) from within the same litter. This helped minimize the within-litter variability, which is not uncommon in endocrine-disruptor models. For remaining evaluations, one mouse per litter was randomly selected as a litter representative. More than one islet-specific end point was determined for each litter, where possible (see Figure S2). A combination of studied outcomes from pooled as well as from an individual animal from each litter helped in assessing both nonrandom and stochastic biological changes.
Insulin Measurement and Islet Perifusion Study
Isolation of pancreatic islets, perifusion, and insulin assays were performed as previously described (Li et al. 2010). In brief, islets were isolated by collagenase digestion and cultured with 10 mM glucose in RPMI 1640 medium (Sigma) for 2 d. Islets were perifused with a Krebs-Ringer bicarbonate buffer containing 0.25% BSA at a flow rate of 1 mL/min. Ramps at increments of 0.5 mM/min for glucose, and 0.3 mM/min for [alpha]-ketoisocaproate were performed. Thirty millimolar potassium chloride was used to determine maximum insulin release. Insulin was measured in the perfusates following the ramp studies, as well as in the lysates from freshly isolated islets by homogeneous time-resolved fluorescence technology (Cisbio Kit).
Islet Mitochondrial Function by High-Resolution Respirometry
A high-resolution respirometry was performed in ~800 overnight-cultured intact islets using Oxygraph-2K (Oroboros Instruments, Innsbruck, Austria). The intact islets were diluted with RPMI and 1-1.5 million cells/chamber were added to each chamber of the Oxygraph-2K in 2 mL volume per chamber. For each experiment, two Oroboros machines were used at a time, and in each machine islets from one of the BPA-exposed groups and Control were run in parallel. The use of chambers for a particular treatment and control was switched between the experiments to avoid any bias due to chamber differences. The experiment was performed at 37[degrees]C with constant magnetic stirring set at 750 rpm. After the addition of intact islets into the oxygraph chambers, the basal respiration referred to as the ROUTINE state was measured when a stable oxygen flux was reached. Following this step, inhibitors for the different mitochondrial respiratory complexes were added and oxygen consumption of mitochondria was measured according to established protocol published previously (Zhang et al. 2013). The inhibitors were added in the following order: oligomycin (3 [micro]M) (Sigma), to inhibit complex V (to assess nonmitochondrial respiratory capacity or proton LEAK rate), carbonyl cyanide-p-trifluoro-methoxyphenylhydrazone (FCCP) (Sigma) uncoupler with stepwise titration in 0.5-to1.5-[micro]M increments [to assess maximal electron transport system respiratory capacity rate (ETS)]. Rotenone (Sigma) was added in a 0.5-[micro]M final concentration to inhibit complex I, and antimycin A (Sigma) in a 2.5-[micro]M final concentration was added last to inhibit complex III. Data was analyzed using DatLab4 and DatLab6 (Oroboros, Austria) software, and normalized to the total DNA concentration.
Islet Cytosolic Calcium Levels
In isolated islets, cytosolic calcium [[[Ca.sup.2+]].sub.i] responses were measured by dual-wavelength fluorescence microscopy using Fura-2AM (Life Technologies) as a calcium indicator, and calcium dynamics were measured with the Zeiss AxioVision system (Carl Zeiss Microscopy, Thornwood, NY) (Li et al. 2013).
Pancreata were excised, and fixed in 10% formalin (pH 7.0) for 48 h at room temperature, and embedded in paraffin (Tissue-Tek, 4,583). Six sections of 5 [micro]m thickness, 100 [micro]m apart were selected for immunohistochemical (IHC) and immunofluorescence (IF) staining per animal. Prior to any antigen-specific staining, subsets of sections from each animal were H & E stained, and assessed by an independent, blinded pathologist for section quality (see Figure S3). All slides were pretreated in an antigen retriever with either 0.01 M citrate buffer pH 7.6 (Vector Labs), or pH9 Retrieval Solution (DAKO), and blocked using 2% fetal bovine serum.
For immunohistochemical staining, slides were incubated in primary F4/80 antibody (Invitrogen) 1:50 overnight at 4[degrees]C, or CD3 antibody (Santa Cruz) 1:250 for 2 h at room temperature. Secondary anti-rat IgG (Vector Lab) and anti-goat IgG (Vector Lab) antibodies were used in a 1:200 dilution for 30 min at room temperature.
For immunofluorescence staining, slides were incubated with antibodies against insulin (Abcam) 1:50 for 16 h at 4[degrees]C, glucagon (Biogenix) 1:50 for 1 h at room temperature, and somatostatin (Biorad) 1:50 for 1 h at room temperature. Secondary antibodies Alexa Fluor[R] 594 anti-guinea pig (Invitrogen), Alexa Fluor[R] 647 anti-rabbit (Invitrogen), and Alexa Fluor[R] 488 anti-rat (Invitrogen) 1:200 were applied for 1 h at room temperature in the dark. The nuclei were stained using DAPI (Sigma). To quantify staining, slides were digitally scanned at 20 x magnification on an Aperio Scanscope CS-O (brightfield) and Scanscope IF (fluorescent) (Leica Biosystems). Nonfluorescent staining was analyzed by ImageScope v184.108.40.20615 (Leica Biosystems) using Aperio-Color Deconvolution v9.1 algorithm. The area of CD3 and F4/80 staining was calculated within islets, pancreatic lymph nodes, and non-islet tissue by multiplying percent positive staining with total stained area and the product was normalized to total analysis area. Fluorescent staining was analyzed by HALO v2.0.1018 (Indica Lab) using Indica Labs Area Quantification FL v1.0 algorithm. Islets were identified as insulin surrounded by glucagon and somatostatin staining. Beta-, alpha-, and delta-cell mass was calculated by multiplying the average insulin, glucagon, and somatostatin stained area, respectively, within each islet relative to the total area scanned with the pancreatic weight and then corrected for body weight measured at tissue harvest.
The cytokine and chemokine panels were measured in pancreatic lysates by mouse Cytokine/Chemokine Luminex 25plex Assay with a detection limit of 3.2pg/mL to 10,000pg/mL (EMD Millipore), as described previously (Jaeckle Santos et al. 2014). Interleukin-6 (IL6) and monocyte chemoattractant protein-1 (MCP1) levels were confirmed using mouse ELISA kits (Thermo Fisher). Cytokine levels were normalized to total protein concentration measured by BCA (bicinchoninic acid) assay (Pierce).
Total RNA was extracted from pooled islets using the Trizol and Qiagen RNeasy mini kit following the manufacturer's instructions. RNA integrity was determined using Agilent RNA 6000 Nano Kit and samples with an RNA integrity number of >7 were used to make RNA libraries. TruSeq Illumina RNA libraries were sequenced using the high-throughput whole transcriptome next generation sequencing to generate 60-80 million reads per sample. Fastq data was processed with HTSeq-count v0.6.1, and aligned using STAR v2.4.0 software. The edgeR v3.12.1 was used to find differentially expressed genes via [log.sub.2]FC and false discovery rate values. Functional analyses were conducted using Ingenuity Pathway Analysis (Ingenuity Systems; http://www. ingenuity.com).
Islet mRNA Expression Analysis
RNA extracted from pooled islets with an RNA integrity number of >7 were used for gene expression analysis. Using Biorad iScript kit, cDNA was subsequently generated. Real-time polymerase chain reaction (PCR) was performed using TaqMan Gene expression assay in Applied Biosystems 7900HT Fast Real-Time PCR System. F1 data was normalized to Rpl19 and Beta actin, and F2 data was normalized to Rpl19 and Cyclophilin A (Ppia). The primers and probes are listed in Table S1.
Islet Protein Expression Analysis
The levels of estrogen receptor alpha (ESR1) and uncoupling protein 2 (UCP2) proteins were measured in lysates by using mouse ESR1 (Aviva Systems Biology), and mouse UCP2-Ab (MyBioSource) ELISA kits, following the manufacturer's instructions. Target protein levels were normalized to total protein concentration measured by BCA assay (Pierce).
Islet DNA Methylation Analysis
DNA was extracted from pooled islets using the QIAmp DNA Mini Kit (Qiagen), and 1 [micro]g of DNA was bisulfite treated using the EpiTect Bisulfite Kit (Qiagen) following the manufacturer's protocol. Pyrosequencing was conducted to measure the methylation status of the promoter region of the Igf2 DMR1 (chr7: 149,851,180-149,851,655; NCBI37/MM9) and H19/Igf2 ICR (chr7: 149,767,599-149,767,819; NCBI37/MM9), and exon A and exon C of the Esr1 gene as described previously by Susiarjo et al. (2013) and Kundakovic et al. (2013), respectively. Briefly, 50 ng of purified bisulfite-treated DNA was amplified using a PyroMark PCR kit (Qiagen) and primers specific for the Igf2 DMR1, H19/Igf2 ICR, and exon A and exon C of the Esr1 gene (see Table S2). Ten microliters of the biotinylated PCR product was used for each sequencing assay and sequenced using a PyroMark Q96MD (Qiagen) pyrosequencer and PyroMark Gold 96 reagents kit (Qiagen) and specific pyrosequencing primers (see Table S2). Percent methylation levels of CpG sites were quantified using Qiagen's Pyro Q-CpG software.
Islet Viability Assays
Lysates from freshly isolated islets from 8- to 10-wk-old F1 and F2 mice were assessed for apoptosis using the caspase 3 activity fluorometric assay (Abcam), BCL2 levels using the mouse BCL2 ELISA kit (LS Bio), and phospho-AKT (pSer473) and total AKT levels using the AKT (pSer473) + total AKT ELISA kit (Abcam) following the manufacturer's instructions.
For comparison between LowerB with Control, and UpperB with Control, we used Dunnett's test, which accounts for multiple testing. For pyrosequencing data, using Dunnett's test, we statistically analyzed methylation levels at individual CpG sites as well as mean methylation across all CpG sites for our two BPA treatment groups relative to Control. In addition, we performed multivariate analyses of variance followed by repeated measures response design to determine a) methylation levels differ between the groups, b) methylation levels differ between CpG sites, and c) differences in methylation levels across CpG sites interact with study groups. Nonparametric data were log-transformed to approximate a standard distribution, where required. All values are presented as mean [+ or -] SEM. A p-value of <0.05 was considered significant. All data were analyzed using JMP and Prism analysis software.
Dose- and Sex-Specific Effects of Maternal BPA Exposure on Islet Mitochondrial Function in F1 and F2 Adult Offspring
Using our mouse model of BPA exposure, we previously observed impaired GSIS in the F1 and F2 generation adult male offspring of mothers (F0) exposed to 10 [micro]g/kg/d (LowerB) or 10 mg/kg/d (UpperB) relative to Controls (7% corn oil diet) (Susiarjo et al. 2015). Nutrient-induced insulin secretion requires ATP production from mitochondria (Antinozzi et al. 2002; Ortsater et al. 2002). Mitochondrial dysfunction has been associated with [beta]-cell failure in islets of intrauterine growth-restricted rats (Simmons et al. 2005). To assess whether impaired GSIS in F1 and F2 adult male offspring is due to altered mitochondrial driven insulin secretion, we measured insulin secretion in response to [alpha]-ketoisocaproate, a substrate metabolized in mitochondria. Surprisingly, [alpha]-ketoisocaproate-stimulated insulin secretion in the islets of F1 and F2 LowerB males was comparable to levels in Controls. However, it trended to be reduced in F1 and F2 UpperB males relative to Controls (Figure 1), thus reflecting an abnormal mitochondrial phenotype in islets of UpperB males. We saw no differences in glucose or [alpha]-ketoisocaproate-stimulated insulin secretion in LowerB or UpperB F1 or F2 females compared with Control females (Figure 2). Therefore, all further studies were conducted exclusively in adult male offspring.
To determine whether the phenotype in males was associated with mitochondrial dysfunction, we performed high-resolution respirometry to assess oxygen consumption in intact islet mitochondria using Oroboros. Consistent with the above results, we did not see any difference in oxygen consumption in the intact islets of F1 and F2 LowerB males, but did observe reduced basal and maximal oxygen consumption in the intact islets of F1 and F2 UpperB males relative to Controls (Figure 3). These results are consistent with abnormal mitochondrial function in the islets of the UpperB males, and the reduction in GSIS in F1 and F2 LowerB males may be due to a factor other than a mitochondrial defect.
Others have shown that BPA exposure alters expression of key mitochondrial genes such as oxoglutarate dehydrogenase (Ogdh) and uncoupling protein 2 (Ucp2) in rodent islets (Wei et al. 2011; Song et al. 2012). Because we observed mitochondrial dysfunction in the current study, we sought to determine whether these previously demonstrated mitochondrial targets of BPA are affected. OGDH is the rate-limiting enzyme in the tricarboxylic acid (TCA) cycle, which plays a key role in mitochondrial function. UCP2 is a member of the mitochondrial carrier family whose loss of function mutation can cause hyperinsulinism, whereas its overexpression is associated with the efflux of TCA cycle intermediates out of the mitochondria, thereby reducing acetyl-CoA-oxidizing substrates in the mitochondria (Gonzalez-Barroso et al. 2008; Vozza et al. 2014). Ogdh expression was only mildly increased in LowerB F2 males compared with Controls (Table 1). In contrast, Ucp2 mRNA levels were significantly increased in LowerB and UpperB F1 and F2 males compared with Controls (Table 1). The changes in mitochondrial gene expression are consistent with UpperB F1 and F2 mitochondrial impairment, although these changes did not translate to statistically significant differences in protein expression levels. F1 and F2 UpperB males (F1 4.0 x [10.sup.-4] [+ or -] 1.0 x [10.sup.-5] ng/[micro]g, p = 0.43; F2 6.3 x [10.sup.-4] [+ or -] 1.7 x [10.sup.-5] ng/[micro]g, p = 0.74; n = 5) and LowerB males (F1 4.3 x [10.sup.-4] [+ or -] 1.0 x [10.sup.-5] ng/[micro]g, p = 0.31; F2 6.1 x [10.sup.-4] [+ or -] 1.7 x [10.sup.-5] ng/[micro]g, p = 0.80; n = 5) had increased UCP2 protein levels relative to Controls (F1 2.3 x [10.sup.-4] [+ or -] 1.0 x [10.sup.-5] ng/[micro]g; F2 4.7 x [10.sup.-4] [+ or -] 1.7 x [10.sup.-5] ng/[micro]g; n = 5); however, these findings were not statistically significant.
Dose-Specific Effects of Maternal BPA Exposure on [beta]-cell Mass in F1 and F2 Adult Male Offspring
An alternative hypothesis for reduced insulin secretion is that BPA exposure is associated with a reduction in the number of insulin-secreting cells. To determine whether reduced GSIS in BPA-exposed F1 and F2 male offspring is linked to reduced [beta]-cell mass, we stained pancreatic sections with anti-insulin, glucagon and somatostatin antibodies to quantify [beta]-cell, alpha-, and delta-cell mass, respectively. F1 and F2 LowerB males had reduced [beta]-cell mass adjusted for body weight relative to Controls (Figure 4). F1 and F2 UpperB and Control males had comparable [beta]-cell mass. Alpha-cell mass was similar among the F1 males, but F2 LowerB males had a modest increase in alpha-cell mass compared with Controls (see Figure S4). F1 and F2 LowerB and UpperB had comparable delta mass as Controls (see Figure S4). Thus, LowerB exposure specifically reduced [beta]-cell mass across two generations. We also measured pancreatic insulin content in F1 and F2 males. UpperB had similar (F1 0.005 [+ or -] 0.001 ng/[micro]g, p = 0.97; F2 0.002 [+ or -] 0.001 ng/[micro]g, p = 0.93; n = 6), whereas LowerB tended to have reduced insulin content (F1 0.002 [+ or -] 0.001 ng/[micro]g, p = 0.10; F2 0.001 [+ or -] 0.0006 ng/[micro]g, p = 0.08; n = 6) relative to Controls (F1 0.004 [+ or -] 0.001 ng/[micro]g; F2 0.003 [+ or -] 0.001 ng/[micro]g; n = 6).
We then examined possible mechanisms underlying the reduction in [beta]-cell mass. Caspase 3 activity, a marker of cell death, was significantly increased in islets of F1 and F2 LowerB males but not in islets from UpperB males (Figure 5). BCL2, an anti-apoptotic protein, did not differ in F1or F2 LowerB compared with Controls, but levels were increased in F1 UpperB males relative to Controls (p = 0.05; Figure 5). These changes were associated with increased AKT phosphorylation in F1 and F2 UpperB males relative to Controls. F1 and F2 LowerB males had similar AKT phosphorylation as Controls (Figure 5). AKT phosphorylation mediates the cell survival response by inhibiting programed cell death (Reed and Paternostro 1999; Yao and Cooper 1995). Thus, the UpperB group had increased AKT phosphorylation and no difference in islet cell death and [beta]-cell mass, whereas the LowerB group has increased islet cell death and decreased [beta]-cell mass across two generations.
Effects of Maternal BPA Exposure on Cytokine/ Chemokine Levels/Expression in Pancreata of F1 and F2 Adult Male Offspring
Several studies indicate that immune cells such as macrophages, dendritic cells, and lymphocytes are present in islets and may play an important role in normal islet development (Criscimanna et al. 2014; Geutskens et al. 2005; Jansen et al. 1993, 1994). In fact, immune cell defects are linked to abnormal pancreatic development and glucose homeostasis observed in athymic BALB/c nude mice (Zeidler et al. 1991). To explore whether [beta]-cell-specific effects of BPA are associated with a perturbed immune response, we assessed the cytokine/chemokine profile in pancreatic lysates of F1 and F2 mice via a Luminex assay. Because different immune cell types produce and respond to different cytokines/chemokines, we have provided a list of these cell types in Table 2; see also Tables S3 and S4. In F1 and F2 BPA-exposed versus Control mice levels of the Th2 cytokine IL4 was reduced, but IL5 was elevated, and pro-inflammatory cytokines IL6 and MCP1 were also elevated (Table 2). Compared with Controls, F1 UpperB males also had increased levels of G-CSF (granulocyte colony stimulating factor) and IL9 (see Table S3). IL17 was increased in F2 LowerB and UpperB males compared with Controls (Table 2). The levels of other cytokines measured were comparable among the three groups (see Tables S3 and S4). We validated the changes in levels of key cytokines via ELISA and saw a similar trend of increased cytokine levels of IL6 as well as significantly increased MCP1 levels in F1 and F2 LowerB and UpperB males relative to Controls (Figure 6). This indicated that both doses of BPA increased levels of pro-inflammatory cytokines in islets.
To verify the increased inflammatory phenotype in the pancreas of F1 and F2 BPA-exposed male mice, we immunostained pancreatic sections with CD3, a marker for T lymphocytes, and with F4/80, a marker for macrophages. We consistently observed increased immunostaining of CD3 and F4/80 in the islets of BPA-exposed F1 and F2 male mice compared with Controls (Figure 7; Table 3). Thus, lower- and upper-dose BPA exposure increased islet inflammation across two generations in males. Taken together these findings indicate the recruitment of macrophages and [T.sub.h]1 cells and, possibly, the recruitment of neutrophils and activation of eosinophils in the pancreas.
Effects of Maternal BPA Exposure on Estrogen Receptor Alpha Gene Expression in Islets of F1 and F2 Adult Male Offspring
The endocrine-disrupting actions of BPA are mediated by both estrogen receptor dependent and independent pathways (Yoon et al. 2014). Estrogen receptors (ESR1, or Esr1) are critically important to [beta]-cell function and development (Tiano et al. 2011; Yuchi et al. 2015) and have anti-inflammatory effects in other organs (Miller et al. 2012). We observed a significant reduction in mRNA expression of Esr1 in F1 and F2 LowerB and UpperB males compared with Controls (Table 1). The Esr1 mRNA changes in the LowerB group correlated with its protein expression, such that F1 and F2 LowerB males (F1 0.0002 [+ or -] 0.001 ng/[micro]g, p = 0.09; F2 0.0001 [+ or -] 0.0003 ng/[micro]g, p = 0.09; n = 5) trended toward reduced ESR1 protein levels relative to Controls (F1 0.0010 [+ or -] 0.001 ng/[micro]g; F2 0.0011 [+ or -] 0.0003 ng/[micro]g; n = 5); however, F1 and F2 UpperB males (F1 0.0028 [+ or -] 0.001 ng/[micro]g, p = 0.45; F2 0.0006 [+ or -] 0.0003 ng/[micro]g, p = 0.48; n = 5) had ESR1 protein levels that were similar to those of Controls. BPA also acts via the estrogen receptor beta (ESR2, or Esr2) and alters stimulus-secretion coupling through its effects on [K.sub.ATP] channels and increasing glucose-stimulated [[[Ca.sup.2+]].sub.i] levels (Soriano et al. 2009; Soriano et al. 2012). We did not observe any difference in either mRNA levels of Esr2, [K.sub.ATP] channels subunits Kcnj11 and Abcc8, and Snap25, which mediate the calcium signaling necessary for insulin secretion (Tables 1, 4, and 5), or in cytosolic [Ca.sup.2+] levels in isolated islets of F1 and F2 LowerB and UpperB adult males relative to Controls (see Figure S5). These results demonstrate that the effects of BPA exposure are mediated preferentially by Esr1 and not Esr2 or targets of Esr2 in our study.
To determine whether BPA mediates change in islets Esr1 gene expression via DNA methylation as shown by others in mouse (Kundakovic et al. 2013) and rat (Chang et al. 2016b) brain, we assessed methylation status of 16 CpG sites across untranslated exon A and exon C of Esr1 gene in islets. We observed similar percent methylation at all five CpG sites of exon C among the three groups (see Figure S6). We observed minimal (4-5%), yet dose-specific, changes in the DNA methylation levels of exon A in islets (Figure 8). F1 LowerB males trended toward increased methylation levels at CpG sites 1 and 4, but F2 LowerB males had similar percent methylation level at all CpG sites examined relative to Controls (Figure 8). The mean methylation across all CpG sites was also similar in F1 and F2 LowerB males compared with Controls. In contrast, F2 UpperB males trend toward increased methylation at CpG sites 3, 4, 5, and 10 and mean methylation across all CpG sites relative to Controls (Figure 8). F1 UpperB males had similar methylation at individual CpG sites and mean methylation across all CpG sites relative to Controls (Figure 8). On repeated measures analysis, we observed significant differences in methylation levels between CpG sites (p < 0.01); however, we saw no statistical difference in methylation levels between the groups (p > 0.1) or groups*CpG site interaction (p > 0.1) in islets of F1 and F2 males. Overall these findings indicate that BPA had minimal effects on Esr1 methylation in islets. It is possible that in addition to these minimal differences in methylation levels, BPA may be mediating its effect on Esr1 via other regulatory mechanisms such as changes in chromatin acetylation /methylation levels. This requires further investigation.
Effects of Maternal BPA Exposure on Igf2 Expression and DNA Methylation at Igf2 DMR1 in Islets ofF1 and F2 Adult Male Offspring
Insulin-like growth factor 2 (Igf2) plays a critical role in [beta]-cell function, proliferation, and apoptosis. We and others have reported that Igf2 overexpression is associated with impaired glucose tolerance (Hoog et al. 1996; Susiarjo et al. 2015) and impaired insulin secretion (Hoog et al. 1996; Petrik et al. 1999). Overexpression of Igf2 has also been reported to lead to [beta]-cell dysfunction and marked morphological abnormalities of the islet (Casellas et al. 2015). In previous studies, we observed an increased expression of Igf2 in F1 and F2 upper dose BPA-exposed embryos (Susiarjo et al. 2013). Thus, to determine whether Igf2 expression is also increased postnatally in islets and whether these changes persist into the next generation, we measured mRNA levels in islets from F1 and F2 LowerB and UpperB groups by quantitative PCR (qPCR). We found that Igf2 mRNA levels were significantly increased in F1 and F2 LowerB and UpperB males compared with Controls (Table 1).
Previously, we determined that the increase in expression of Igf2 in F1 and F2 islets in UpperB males was associated with increased methylation at the Igf2 DMR1 (Susiarjo et al. 2013, 2015). Thus, we measured DNA methylation at H19/Igf2 imprinting control region (ICR) and at Igf2 DMR1 in islets. F1 and F2 LowerB males had DNA methylation levels at Igf2 DMR1 and H19/Igf2 ICR comparable with Controls (Figure 8; see also Figure S6). In contrast, F1 UpperB males had increased methylation at CpG site 1 and increased mean methylation, and F2 UpperB males had increased methylation at CpG sites 1 and 2 as well as increased mean methylation at Igf2 DMR1 relative to Controls (Figure 8). Interestingly, on repeated measures, we observed differences in methylation levels across CpG sites (F1 p < 0.01; F2 p <0.01) and between groups (F1 p = 0.04; F2 p = 0.02), but no statistical difference in groups*CpG site interaction (F1 p > 0.1; F2 p > 0.1) at Igf2 DMR1. This analysis indicated that BPA affects methylation level at CpG sites but does not account for the methylation differences between CpG sites at Igf2 DMR1. However, we saw no difference in methylation at H19/Igf2 ICR (see Figure S6), suggesting that the increase in Igf2 expression in the islet results from increased expression of the normally expressed paternal allele and not loss of imprinting.
Effects of Maternal BPA Exposure on the Transcriptome of F1 Adult Male and Female Offspring
To identify novel pathways and to determine whether changes in islet phenotypes in BPA-exposed mice are regulated at the gene expression level, we conducted whole transcriptome analysis via RNA sequencing in males and females. Consistent with our phenotypic data, we saw no differentially expressed genes in females (Tables 6 and 7). Because of interanimal variability and the small sample size, there were very few differentially expressed genes in males when the analysis was corrected for multiple testing (q < 0.05). Interestingly, however, of the few differentially expressed genes, there was minimal overlap between the LowerB and UpperB males (Figure 9), consistent with the observed dose-specific phenotypes. Moreover, qPCR of several key islet genes demonstrated robust differences between the groups despite higher q-values in the RNAseq analysis (Tables 1, 4, 5). Thus, we performed ingenuity pathway analysis on the RNAseq data using a p <0.05 and a fold change of 1.5. Multiple genes associated with inflammation in LowerB males, and both inflammation and mitochondrial function including branched-chain [alpha]-keto acid dehydrogenase complex in UpperB males were altered compared with Controls, which is consistent with the phenotypic data (Figure 9). Other pathways that are critical for normal pancreatic development and regulation of insulin secretion such as sonic hedgehog (Thomas et al. 2000) and gap-junction signaling (Bavamian et al. 2007; Klee et al. 2011) were altered in LowerB males, whereas pathways involved in inflammation and mitochondrial dysfunction such as G-protein couple receptor (Patial et al. 2009; Sorriento et al. 2013) and endoplasmic reticulum stress (Rocha et al. 2016) and signaling were altered in UpperB males. Unexpectedly, there were minimal changes in mRNA expression levels as determined by RNA Seq and qPCR of other [beta]-cell regulatory genes such as Pdx1, Igf1, and Hnf1a in BPA groups compared with Controls (Tables 1, 4, and 5).
We recently reported that early-life exposure to the endocrine-disrupting chemical BPA in the mouse is associated with sexspecific metabolic health defects across multiple generations (Susiarjo et al. 2015). In the current study, we provide further evidence that BPA has sex-specific effects, such that the islet function of only males is affected. Additionally, we observed dose-specific effects in male offspring of dams fed a lower-dose versus a higher-dose BPA diet. We show here that the insulin secretory defect in male offspring is associated with impaired mitochondrial function in UpperB mice and reduced [beta]-cell mass in the LowerB mice (Table 8). It is not uncommon to find dose-specific effects following exposure to endocrine-disrupting chemicals (Vandenberg et al. 2012). Various studies in the mouse have shown that BPA affects metabolism (Alonso-Magdalena et al. 2010), DNA methylation levels (Kim et al. 2014), and behavior (Jasarevic et al. 2013) in a dose-specific manner. The dose-specific findings of the current study are consistent with these previous studies. Uniquely, our current study shows that the dose-specific effects of BPA on the pancreas are not limited to the first, but are transmitted to the second generation. These persistent phenotypic changes into at least the second generation are not due to an altered metabolic milieu in the pregnant F1 female, and thus may be transmitted through epigenetic mechanisms or via mitochondria (Susiarjo et al. 2015). Our current findings of persistent changes in DNA methylation at Igf2 and altered mitochondria function in [beta] cells in F1 and F2 offspring suggest that both mechanisms are involved.
Mitochondrial impairment, and altered expression levels of key mitochondrial genes such as Ucp2 and Ogdh have been previously reported in the pancreas of first-generation adult male offspring of BPA-exposed rat dams (Wei et al. 2011), in isolated rat islets (Song et al. 2012), and in BPA-exposed INS1 cells (Lin et al. 2013). Importantly, others have shown that increased Ucp2 expression increases efflux of acetyl-CoA-oxidizing substrates such as [alpha]-ketoglutarate of the TCA cycle out of the mitochondria (Vozza et al. 2014). It is known that [alpha]-ketoisocaproate promotes glutamate oxidation to produce [alpha]-ketoglutarate (Zielke et al. 1997). Therefore, the increased Ucp2 expression with an efflux of TCA intermediates from the mitochondria as previously shown (Vozza et al. 2014) suggests altered [alpha]-ketoisocaproate metabolism. Although in the current study we observed increased Ucp2 mRNA expression levels in BPA-exposed groups, we did not find statistically increased UCP2 protein levels. However, our transcriptome analysis revealed differential expression of other genes such as Ucp2 that are also involved in the branched-chain [alpha]-ketoacid dehydrogenase complex, which catalyzes the metabolism of branched-chain amino acids such as leucine and its derivative, [alpha]-ketoisocaproate. This is consistent with a reduction in [alpha]-ketoisocaproate-stimulated insulin secretion in the UpperB group. Therefore, it may be that genes other than Ucp2 impair [alpha]-ketoisocaproate metabolism, consequently affecting [alpha]-ketoisocaproate-stimulated insulin secretion.
Reduced [alpha]-ketoisocaproate-stimulated insulin secretion and mitochondrial dysfunction persisting until the second generation in islets of the UpperB group indicate that the higher dose of BPA is associated with mitochondrial impairment. Mitochondrial impairment is associated with the production of reactive oxygen species (ROS) such as superoxide, hydrogen peroxide, and hydroxyl radicals as seen in different type 2 diabetes models and IUGR rats (Bindokas et al. 2003; Ihara et al. 1999; Simmons et al. 2005). ROS can lead to the activation of NF-[kappa]B, which in turn increases proinflammatory cytokine production (Kepp et al. 2011), eventually leading to increased [beta]-cell death and dysfunction. Our study demonstrates that UpperB exposure is associated with impaired mitochondrial function and increased pro-inflammatory cytokine production. Whether this is due to oxidative stress remains to be determined. Elevated pro-inflammatory cytokine levels have been shown to induce [beta]-cell death (Sauter et al. 2015), which in turn reduces [beta]-cell mass. AKT phosphorylation inhibits programed cell death (Reed and Paternostro 1999; Yao and Cooper 1995); in the current study, increased AKT phosphorylation levels and the lack of changes in cell death observed in UpperB islets are consistent with these previous findings. The mechanism underlying this change remains to be determined.
Despite the lack of changes in mitochondrial function in the LowerB group, the first- and second-generation offspring of this group had significantly reduced [beta]-cell function as defined by impaired insulin secretion, which suggested that a lower dose of BPA may act via an alternate pathway in [beta] cells, such as a reduction in [beta]-cell mass. Impaired insulin secretion is known to be secondary to reduced [beta]-cell mass in humans (Butler et al. 2003; Rahier et al. 2008; Yoon et al. 2003) and animals (Bansal et al. 2015; Kjems et al. 2001; Movassat et al. 1997; Simmons et al. 2001). Indeed, we observed reduced [beta]-cell mass adjusted for body weight across two generations in adult male offspring of LowerB exposed dams, and these changes were associated with increased cell death. These finding are consistent with previous studies that reported reduced [beta]-cell mass at birth in the offspring of rats exposed to 10 [micro]g BPA/kg body weight per day throughout gestation and lactation (Chang et al. 2016a) and increased cell death in BPA-exposed INS1 cells (Gong et al. 2013; Lin et al. 2013). [beta]-cell mass was also reduced and cell death was increased at 6 mo after delivery in female mice that were exposed to BPA during pregnancy (Alonso-Magdalena et al. 2015a). Interestingly, when C57BL/6 mice dams were subcutaneously injected from either day 6 of pregnancy (P6) until birth (PND0), PND0 until weaning (PND21), or P6-PND21 with 100 [micro]g BPA/kg body weight per day, [beta]-cell mass was increased in 3-mo-old male offspring (Liu et al. 2013). Differences in timing, dose, route of exposure, and methods of measurement of [beta]-cell mass that were not corrected for body weight may explain the different results in our study compared with those of Liu et al. (2013). In rodents, ([beta]-cell mass has a linear correlation with body weight such that [beta]-cell mass increases with body weight (Montanya et al. 2000). Therefore, we corrected [beta]-cell mass for body weight in our study, so that any changes due to difference in body weight per se would be accounted for. We chose oral BPA exposure via diet as it has been shown to be effective and a representative model of human exposures because humans are chronically exposed to BPA primarily through ingestion, and our doses are well within the human exposure limits (Susiarjo et al. 2013; Sieli et al. 2011). Importantly, end points measured in this study were from animals generated across multiple cohorts studied in two different animal care facilities with a control group included in each cohort. This reflects robustness of the studied outcomes and cohort consistency in our paradigm.
Our study is the first to show that developmental BPA exposure induces inflammation in the islet in both F1 and F2 generations. Although immune cells are present in the pancreas and play an important role in normal pancreatic development (Criscimanna et al. 2014; Geutskens et al. 2005; Jansen et al. 1993, 1994), several studies have reported that increased cytokine levels, and increased immune cell infiltration in islets, is associated with [beta]-cell dysfunction in type 2 diabetes (Boni-Schnetzler et al. 2008; Ehses et al. 2007; Homo-Delarche et al. 2006; Kim et al. 2005; Lin et al. 2012; Maedler et al. 2002; Richardson et al. 2009). Increased levels of pro-inflammatory cytokines such as IL6 and MCP1 reduce insulin secretion in human and mouse islets (Sauter et al. 2015). Inflammatory cytokines also impair mitochondria function in the [beta] cell, which in turn blunts insulin secretion (Gao et al. 2015). Thus, it is likely that BPA-induced islet inflammation plays a key role in the abnormal [beta]-cell phenotype observed in F1 and F2 mice.
BPA is an estrogen mimic, and therefore, it is plausible that BPA could alter Esr1 expression. In fact, maternal BPA exposure has been shown to down-regulate Esr1 mRNA and protein expression at P7 and P21 in the hippocampus of male rat offspring (Xu et al. 2014). Two injections of ICI 182,780, a brain-permeable endoplasmic reticulum (ER) antagonist, in rat pups were able to reverse this change, demonstrating that BPA is causal in reducing Esr1 expression (Xu et al. 2014). Consistent with these findings, we observed reduced Esr1 gene expression in islets of both BPA exposure groups, and reduced ESR1 protein levels in lower-dose exposure group. However, because BPA has weak binding affinity for estrogen receptors relative to estrogens, it is likely that BPA mediates its effect on estrogen receptor expression without directly binding to the receptor. Other mechanisms, such as a direct effect on DNA methylation, may instead be responsible for altering Esr1 gene expression as shown in the mouse (Kundakovic et al. 2013) and rat brain (Chang et al. 2016b). We assessed DNA methylation status at 16 CpG sites in exon A and exon C of Esr1 gene in islets, and found minimal differences. It is possible that other gene regulatory mechanisms such as a change in chromatin conformation may be involved. This remains to be determined.
Estrogen receptors are widely distributed in the thymus (Kawashima et al. 1992), bone marrow (Smithson et al. 1998), and spleen (Samy et al. 2003), suggesting that they play a role in regulating normal immune response. In fact, loss of Esr1 leads to perturbed development of the thymus and spleen (Erlandsson et al. 2001) and to altered immune responses (Douin-Echinard et al. 2008; Lambert et al. 2005; Maret et al. 2003). Altered Esr1 expression could therefore influence the effect of developmental exposure to BPA on the pancreatic immune response, leading to the increased inflammation that we observed in the islets of the BPA-exposed animals in this study.
It is remarkable that the changes we reported in Igf2 expression and DNA methylation in the embryo persist in the islet in adulthood and into the next generation. Further, in previous studies, the metabolic phenotypes in the F1 and F2 BPA male offspring were linked to Igf2 overexpression (Susiarjo et al. 2015). Although Igf2 is protective in the adult [beta] cell and increased expression prevents [beta]-cell death (Cornu et al. 2009; Hughes et al. 2013), overexpression of Igf2 in the fetus results in marked [beta]-cell abnormalities in fetal and adult life (Casellas et al. 2015; Hoog et al. 1996; Petrik et al. 1999). These [beta]-cell abnormalities may be secondary to increased islet inflammation because Igf2 overexpression has been reported to be associated with increased immune cell infiltration in the rodent pancreas (Casellas et al. 2015). Therefore, it is plausible to suggest that Igf2 in combination with reduced Esr1 mediates the increased inflammatory response that we observed in islets of the BPA-exposed animals.
Finally, ingenuity pathway analysis of the RNAseq data revealed multiple pathways that were altered and could contribute to the abnormal pancreatic phenotype in the offspring of BPA-exposed mice, including sonic hedgehog (Shh) and gap-junction signaling in the LowerB dose group and ER stress in the UpperB dose group. Shh is known to play a critical role in pancreatic development (Thomas et al. 2000), and gap junctions help establish cross-talk within islet cells thereby regulating normal insulin secretion and [beta]-cell mass (Bavamian et al. 2007; Klee et al. 2011). ER stress is associated with mitochondrial dysfunction (Rocha et al. 2016). Although we did not further interrogate these pathways in our study, they remain important targets for future research.
The current study elucidates the underlying mechanisms responsible for impaired insulin secretion in mice exposed to relevant human exposure levels of BPA during development. The persistent changes in insulin secretion, mitochondria function, and [beta]-cell mass across multiple generations are consistent with an epigenetically mediated mechanism. These findings have public health ramifications and identify [beta]-cell-specific targets of BPA that could be used to develop potential therapeutic strategies to prevent environment-induced diabetes.
Please Note: Illustration(s) are not available due to copyright restrictions.
The authors would like to acknowledge D. Condon, Z. (Paul) Wang, T. Kim, and M. Li for their bioinformatics and biostatistics input. We are also grateful for the services of the University of Pennsylvania Next Generation Sequencing Core and the Radioimmunoassay and Biomarkers Core (P30 DK19525), and the Children's Hospital of Philadelphia Pathology Core.
This work is supported by the National Institute of Environmental Health Sciences/National Institutes of Health (ES023284 and ES013508 [M.S.B., R.A.S.]); March of Dimes (M.S.B.); the Perelman School of Medicine Center of Excellence in Environmental Toxicology (CEET-ES-013508-05 [R.A.S.], RO1DK098517 [C.L.], and T32 ES019851 [F.X.]), the Bioenergetics Core of the Children's Hospital of Philadelphia and Development Disabilities Research Center (IDDRC; U54-HD086984), and the Islet Cell Biology Core of the Penn Diabetes Institute for Diabetes, Obesity & Metabolism (P30DK19525).
Aekplakorn W, Chailurkit LO, Ongphiphadhanakul B. 2015. Relationship of serum bisphenol A with diabetes in the Thai population, National Health Examination Survey IV, 2009. J Diabetes 7(2):240-249, PMID: 24720399, https://doi.org/10. 1111/1753-0407.12159.
Ahmadkhaniha R, Mansouri M, Yunesian M, Omidfar K, Jeddi MZ, Larijani B, et al. 2014. Association of urinary bisphenol A concentration with type-2 diabetes mellitus. J Environ Health Sci Eng 12(1):64, PMID: 24625016, https://doi.org/10. 1186/2052-336X-12-64.
Alonso-Magdalena P, Garcia-Arevalo M, Quesada I, Nadal A. 2015a. Bisphenol-A treatment during pregnancy in mice: a new window of susceptibility for the development of diabetes in mothers later in life. Endocrinology 156:1659-1670, PMID: 25830705, https://doi.org/10.1210/en.2014-1952.
Alonso-Magdalena P, Morimoto S, Ripoll C, Fuentes E, Nadal A. 2006. The estrogenic effect of bisphenol A disrupts pancreatic beta-cell function in vivo and induces insulin resistance. Environ Health Perspect 114(1):106-112, PMID: 16393666, https://doi.org/10.1289/ehp.8451.
Alonso-Magdalena P, Quesada I, Nadal A. 2015b. Prenatal exposure to BPA and offspring outcomes: the diabesogenic behavior of BPA. Dose Response 13(2): 1559325815590395, PMID: 26676280, https://doi.org/10.1177/1559325815590395.
Alonso-Magdalena P, Vieira E, Soriano S, Menes L, Burks D, Quesada I, et al. 2010. Bisphenol A exposure during pregnancy disrupts glucose homeostasis in mothers and adult male offspring. Environ Health Perspect 118(9):1243-1250, PMID: 20488778, https://doi.org/10.1289/ehp.1001993.
American Diabetes Association. 2013. Economic costs of diabetes in the U.S. in 2012. Diabetes Care 36:1033-1046, PMID: 23468086, https://doi.org/10.2337/ dc12-2625.
Angle BM, Do RP, Ponzi D, Stahlhut RW, Drury BE, Nagel SC, et al. 2013. Metabolic disruption in male mice due to fetal exposure to low but not high doses of bisphenol A (BPA): evidence for effects on body weight, food intake, adipocytes, leptin, adiponectin, insulin and glucose regulation. Reprod Toxicol 42:256-268, PMID: 23892310, https://doi.org/10.1016/j.reprotox.2013.07.017.
Antinozzi PA, Ishihara H, Newgard CB, Wollheim CB. 2002. Mitochondrial metabolism sets the maximal limit of fuel-stimulated insulin secretion in a model pancreatic beta cell: a survey of four fuel secretagogues. J Biol Chem 277(14):11746-11755, PMID: 11821387, https://doi.org/10.1074/jbc. M108462200.
Bansal A, Bloomfield FH, Connor KL, Dragunow M, Thorstensen EB, Oliver MH, et al. 2015. Glucocorticoid-induced preterm birth and neonatal hyperglycemia alter ovine [beta]-cell development. Endocrinology 156(10):3763-3776, PMID: 26204462, https://doi.org/10.1210/en.2015-1095.
Batista TM, Alonso-Magdalena P, Vieira E, Amaral ME, Cederroth CR, Nef S, et al. 2012. Short-term treatment with bisphenol-A leads to metabolic abnormalities in adult male mice. PLoS One 7(3):e33814, PMID: 22470480, https://doi.org/10. 31371/journal.pone.0033814.
Bavamian S, Klee P, Britan A, Populaire C, Caille D, Cancela J, et al. 2007. Islet-cell-to-cell communication as basis for normal insulin secretion. Diabetes Obes Metab 9(suppl 2):118-132, PMID: 17919186, https://doi.org/10.111Vj.14631326.2007.00780.x
Bindokas VP, Kuznetsov A, Sreenan S, Polonsky KS, Roe MW, Philipson LH. 2003. Visualizing superoxide production in normal and diabetic rat islets of Langerhans. J Biol Chem 278(11):9796-9801, PMID: 12514170, https://doi.org/10. 1074/jbc.M206913200.
Boni-Schnetzler M, Thorne J, Parnaud G, Marselli L, Ehses JA, Kerr-Conte J, et al. 2008. Increased interleukin (IL)-1 b messenger ribonucleic acid expression in [beta]-cells of individuals with type 2 diabetes and regulation of IL-1 b in human islets by glucose and autostimulation. J Clin Endocrinol Metab 93(10):4065-4074, PMID: 18664535, https://doi.org/10.1210/jc.2008-0396.
Bouwens L, Rooman I. 2005. Regulation of pancreatic beta-cell mass. Physiol Rev 85(4):1255-1270, PMID: 16183912, https://doi.org/10.1152/physrev.00025.2004.
Butler AE, Janson J, Bonner-Weir S, Ritzel R, Rizza RA, Butler PC. 2003. Beta- cell deficit and increased [beta]-cell apoptosis in humans with type 2 diabetes. Diabetes 52(1):102-110, PMID: 12502499.
Casellas A, Mallol C, Salavert A, Jimenez V, Garcia M, Agudo J, et al. 2015. Insulin-like growth factor 2 overexpression induces [beta]-cell dysfunction and increases beta-cell susceptibility to damage. J Biol Chem 290(27):16772-16785, PMID: 25971976, https://doi.org/10.1074/jbc.M115.642041.
Cavaillon JM. 2001. Pro- versus anti-inflammatory cytokines: myth or reality. Cell Mol Biol (Noisy-le-grand) 47(4):695-702, PMID: 11502077.
Chang H, Wang M, Xia W, Chen T, Huo W, Mao Z, et al. 2016b. Perinatal exposure to low-dose bisphenol A disrupts learning/memory and DNA methylation of estrogen receptor alpha in the hippocampus. Toxicol Res 5:828-835, https://doi.org/10. 1039/C5TX00449G.
Chang H, Wang D, Xia W, Pan X, Huo W, Xu S, et al. 2016a. Epigenetic disruption and glucose homeostasis changes following low-dose maternal bisphenol A exposure. Toxicol Res 5:1400-1409, https://doi.org/10.1039/C6TX00047A.
Commins SP, Borish L, Steinke JW. 2010. Immunologic messenger molecules: cytokines, interferons, and chemokines. J Allergy Clin Immunol 125(2 suppl 2):S53-S72, PMID: 19932918, https://doi.org/10.1016/jjaci.2009.07.008.
Cornu M, Yang JY, Jaccard E, Poussin C, Widmann C, Thorens B. 2009. Glucagon- like peptide-1 protects [beta]-cells against apoptosis by increasing the activity of an Igf-2/Igf-1 receptor autocrine loop. Diabetes 58(8):1816-1825, PMID: 19401425, https://doi.org/10.2337/db09-0063.
Criscimanna A, Coudriet GM, Gittes GK, Piganelli JD, Esni F. 2014. Activated macrophages create lineage-specific microenvironments for pancreatic acinar- and [beta]-cell regeneration in mice. Gastroenterology 147(5):1106-1118.e11, PMID: 25128759, https://doi.org/10.1053/j.gastro.2014.08.008.
Douin-Echinard V, Laffont S, Seillet C, Delpy L, Krust A, Chambon P, et al. 2008. Estrogen receptor a, but not b, is required for optimal dendritic cell differentiation and [corrected] CD40-induced cytokine production. J Immunol 180(6): 3661-3669, PMID: 18322171.
Ehses JA, Perren A, Eppler E, Ribaux P, Pospisilik JA, Maor-Cahn R, et al. 2007. Increased number of islet-associated macrophages in type 2 diabetes. Diabetes 56(9):2356-2370, PMID: 17579207, https://doi.org/10.2337/db06-1650.
Erlandsson MC, Ohlsson C, Gustafsson JA, Carlsten H. 2001. Role of oestrogen receptors a and b in immune organ development and in oestrogen-mediated effects on thymus. Immunology 103(1):17-25, PMID: 11380688, https://doi.org/ 10.1046/j.1365-2567.2001.01212.x.
Finegood DT, Scaglia L, Bonner-Weir S. 1995. Dynamics of [beta]-cell mass in the growing rat pancreas. Estimation with a simple mathematical model. Diabetes 44(3):249-256, PMID: 7883109, https://doi.org/10.2337/diabetes.44 3.249.
Gao J, Sang M, Zhang X, Zheng T, Pan J, Dai M, et al. 2015. Miro1-mediated mitochondrial dysfunction under high nutrient stress is linked to NOD-like receptor 3 (NLRP3)-dependent inflammatory responses in rat pancreatic beta cells. Free Radic Biol Med 89:322-332, PMID: 26427885, https://doi.org/10.1016/j.freeradbiomed.2015. 09.002.
Geutskens SB, Otonkoski T, Pulkkinen MA, Drexhage HA, Leenen PJ. 2005. Macrophages in the murine pancreas and their involvement in fetal endocrine development in vitro. J Leukoc Biol 78(4):845-852, PMID: 16037409, https://doi.org/10. 1189/jlb.1004624.
Gong H, Zhang X, Cheng B, Sun Y, Li C, Li T, et al. 2013. Bisphenol A accelerates toxic amyloid formation of human islet amyloid polypeptide: a possible link between bisphenol A exposure and type 2 diabetes. PLoS One 8(1):e54198, PMID: 23372685, https://doi.org/10.1371/journal.pone.0054198.
Gonzalez-Barroso MM, Giurgea I, Bouillaud F, Anedda A, Bellanne-Chantelot C, Hubert L, et al. 2008. Mutations in UCP2 in congenital hyperinsulinism reveal a role for regulation of insulin secretion. PLoS One 3(12):e3850, PMID: 19065272, https://doi.org/10.1371/journal.pone.0003850.
Hellerstrom C, Swenne I, Andersson A. 1988. Islet replication and diabetes. In: The Pathology of the Endocrine Pancreas in Diabetes. Lefebvre PJ, Pipeleers DG, eds. Berlin, Germany:Springer-Verlag, 141-170.
Homo-Delarche F, Calderari S, Irminger JC, Gangnerau MN, Coulaud J, Rickenbach K, et al. 2006. Islet inflammation and fibrosis in a spontaneous model of type 2 diabetes, the GK rat. Diabetes 55(6):1625-1633, PMID: 16731824, https://doi.org/10.2337/db05-1526.
Hoog A, Sandberg-Nordqvist AC, Abdel-Halim SM, Carlsson-Skwirut C, Guenifi A, Tally M, et al. 1996. Increased amounts of a high molecular weight insulin-like growth factor II (IGF-II) peptide and IGF-II messenger ribonucleic acid in pancreatic islets of diabetic Goto-Kakizaki rats. Endocrinology 137(6):2415-2423, PMID: 8641194, https://doi.org/10.1210/endo.137.6.8641194.
Hughes A, Mohanasundaram D, Kireta S, Jessup CF, Drogemuller CJ, Coates PT. 2013. Insulin-like growth factor-II (IGF-II) prevents proinflammatory cytokine- induced apoptosis and significantly improves islet survival after transplantation. Transplantation 95(5):671-678, PMID: 23364485, https://doi.org/10.1097/TP. 0b013e31827fa453.
Ihara Y, Toyokuni S, Uchida K, Odaka H, Tanaka T, Ikeda H, et al. 1999. Hyperglycemia causes oxidative stress in pancreatic beta-cells of GK rats, a model of type 2 diabetes. Diabetes 48(4):927-932, PMID: 10102716, https://doi.org/10.2337/diabetes.48.4.927.
Jaeckle Santos LJ, Li C, Doulias PT, Ischiropoulos H, Worthen GS, Simmons RA. 2014. Neutralizing Th2 inflammation in neonatal islets prevents [beta]-cell failure in adult IUGR rats. Diabetes 63(5):1672-1684, PMID: 24408314, https://doi.org/10. 2337/db13-1226.
Jansen A, Homo-Delarche F, Hooijkaas H, Leenen PJ, Dardenne M, Drexhage HA. 1994. Immunohistochemical characterization of monocytes-macrophages and dendritic cells involved in the initiation of the insulitis and [beta]-cell destruction in NOD mice. Diabetes 43(5):667-675, PMID: 8168644, https://doi.org/10.2337/diabetes.43.5.667.
Jansen A, Voorbij HA, Jeucken PH, Bruining GJ, Hooijkaas H, Drexhage HA. 1993. An immunohistochemical study on organized lymphoid cell infiltrates in fetal and neonatal pancreases. A comparison with similar infiltrates found in the pancreas of a diabetic infant. Autoimmunity 15(1):31-38, PMID: 8218828, https://doi.org/10.3109/08916939309004836.
Jasarevic E, Sieli PT, Twellman EE, Welsh TH Jr, Schachtman TR, Roberts RM, et al. 2011. Disruption of adult expression of sexually selected traits by developmental exposure to bisphenol A. Proc Natl Acad Sci U S A 108(28):11715-11720, PMID: 21709224, https://doi.org/10.1073/pnas.1107958108.
Jasarevic E, Williams SA, Vandas GM, Ellersieck MR, Liao C, Kannan K, et al. 2013. Sex and dose-dependent effects of developmental exposure to bisphenol A on anxiety and spatial learning in deer mice (Peromyscus maniculatus bairdii) offspring. Horm Behav 63(1):180-189, PMID: 23051835, https://doi.org/10.1016/j. yhbeh.2012.09.009.
Kawashima I, Seiki K, Sakabe K, Ihara S, Akatsuka A, Katsumata Y. 1992. Localization of estrogen receptors and estrogen receptor-mRNA in female mouse thymus. Thymus 20(2):115-121, PMID: 1519316.
Kepp O, Galluzzi L, Kroemer G. 2011. Mitochondrial control of the NLRP3 inflammasome. Nat Immunol 12(3):199-200, PMID: 21321591, https://doi.org/10.1038/ni0311-199.
Kim JH, Sartor MA, Rozek LS, Faulk C, Anderson OS, Jones TR, et al. 2014. Perinatal bisphenol A exposure promotes dose-dependent alterations of the mouse methylome. BMC Genomics 15(1):30, PMID: 24433282, https://doi.org/10. 1186/1471-2164-15-30.
Kim WH, Lee JW, Gao B, Jung MH. 2005. Synergistic activation of JNK/SAPK induced by TNF-[alpha] and IFN-[gamma]: apoptosis of pancreatic [beta]-cells via the p53 and ROS pathway. Cell Signal 17(12):1516-1532, PMID: 15908180, https://doi.org/10. 1016/j.cellsig.2005.03.020.
Kimura A, Kishimoto T. 2010. Il-6: Regulator of Treg/Th17 balance. Eur J Immunol 40(7):1830-1835, PMID: 20583029, https://doi.org/10.1002/eji.201040391.
Kjems LL, Kirby BM, Welsh EM, Veldhuis JD, Straume M, McIntyre SS, et al. 2001. Decrease in [beta]-cell mass leads to impaired pulsatile insulin secretion, reduced postprandial hepatic insulin clearance, and relative hyperglucagonemia in the minipig. Diabetes 50(9):2001-2012, PMID: 11522665.
Klee P, Lamprianou S, Charollais A, Caille D, Sarro R, Cederroth M, et al. 2011. Connexin implication in the control of the murine beta-cell mass. Pediatr Res 70(2):142-147, PMID: 21527868, https://doi.org/10.1203/PDR.0b013e318220f106.
Kourilsky P, Truffa-Bachi P. 2001. Cytokine fields and the polarization of the immune response. Trends Immunol 22(9):502-509, PMID: 11525941, https://doi.org/10. 1016/S1471 -4906(01)02012-9.
Kundakovic M, Gudsnuk K, Franks B, Madrid J, Miller RL, Perera FP, et al. 2013. Sex-specific epigenetic disruption and behavioral changes following low-dose in utero bisphenol A exposure. Proc Natl Acad Sci U S A 110(24):9956-9961, PMID: 23716699, https://doi.org/10.1073/pnas.1214056110.
Lambert KC, Curran EM, Judy BM, Milligan GN, Lubahn DB, Estes DM. 2005. Estrogen receptor alpha (ERa) deficiency in macrophages results in increased stimulation of CD4+ T cells while 17[beta]-estradiol acts through ERa to increase IL-4 and GATA-3 expression in [CD4.sup.+] T cells independent of antigen presentation. J Immunol 175(9):5716-5723, PMID: 16237062.
Lang IA, Galloway TS, Scarlett A, Henley WE, Depledge M, Wallace RB, et al. 2008. Association of urinary bisphenol A concentration with medical disorders and laboratory abnormalities in adults. JAMA 300(11):1303-1310, PMID: 18799442, https://doi.org/10.1001/jama.300.11.1303.
Li C, Chen P, Palladino A, Narayan S, Russell LK, Sayed S, et al. 2010. Mechanism of hyperinsulinism in short-chain 3-hydroxyacyl-CoA dehydrogenase deficiency involves activation of glutamate dehydrogenase. J Biol Chem 285(41):31806-31818, PMID: 20670938, https://doi.org/10.1074/jbc.M110.123638.
Li C, Liu C, Nissim I, Chen J, Chen P, Doliba N, et al. 2013. Regulation of glucagon secretion in normal and diabetic human islets by y-hydroxybutyrate and glycine. J Biol Chem 288(6):3938-3951, PMID: 23266825, https://doi.org/10.1074/jbc. M112.385682.
Li G, Chang H, Xia W, Mao Z, Li Y, Xu S. 2014. F0 maternal BPA exposure induced glucose intolerance of F2 generation through DNA methylation change in Gck. Toxicol Lett 228(3):192-199, PMID: 24793715, https://doi.org/10.1016/j.toxlet.2014. 04.012.
Lin CY, Ni CC, Yin MC, Lii CK. 2012. Flavonoids protect pancreatic beta-cells from cytokines mediated apoptosis through the activation of PI3-kinase pathway. Cytokine 59(1):65-71, PMID: 22579112, https://doi.org/10.1016/j.cyto.2012.04.011.
Lin Y, Sun X, Qiu L, Wei J, Huang Q, Fang C, et al. 2013. Exposure to bisphenol A induces dysfunction of insulin secretion and apoptosis through the damage of mitochondria in rat insulinoma (INS-1) cells. Cell Death Dis 4:e460, PMID: 23328667, https://doi.org/10.1038/cddis.2012.206.
Liu J, Yu P, Qian W, Li Y, Zhao J, Huan F, et al. 2013. Perinatal bisphenol A exposure and adult glucose homeostasis: identifying critical windows of exposure. PLoS One 8(5):e64143, PMID: 23675523, https://doi.org/10.1371/journal.pone.0064143.
Mackay H, Patterson ZR, Khazall R, Patel S, Tsirlin D, Abizaid A. 2013. Organizational effects of perinatal exposure to bisphenol-A and diethylstilbestrol on arcuate nucleus circuitry controlling food intake and energy expenditure in male and female CD-1 mice. Endocrinology 154(4):1465-1475, PMID: 23493373, https://doi.org/10.1210/en.2012-2044.
Maedler K, Sergeev P, Ris F, Oberholzer J, Joller-Jemelka HI, Spinas GA, et al. 2002. Glucose-induced [beta] cell production of IL-1 b contributes to glucotoxicity in human pancreatic islets. J Clin Invest 110(6):851-860, PMID: 12235117, https://doi.org/10.1172/JCI15318.
Maret A, Coudert JD, Garidou L, Foucras G, Gourdy P, Krust A, et al. 2003. Estradiol enhances primary antigen-specific CD4 T cell responses and Th1 development in vivo. Essential role of estrogen receptor alpha expression in hematopoietic cells. Eur J Immunol 33(2):512-521, PMID: 12645950, https://doi.org/10.1002/ immu.200310027.
Mathers CD, Loncar D. 2006. Projections of global mortality and burden of disease from 2002 to 2030. PLoS Med 3(11):e442, PMID: 17132052, https://doi.org/10. 1371/journal.pmed.0030442.
Miller CN, Brown LM, Rayalam S, Della-Fera MA, Baile CA. 2012. Estrogens, inflammation and obesity: an overview. Front Biol 7(1):40-47, PMID: 26811110, https://doi.org/10.1007/s11515-011-1174-y.
Montanya E, Nacher V, Biarnes M, Soler J. 2000. Linear correlation between betacell mass and body weight throughout the lifespan in Lewis rats: role of betacell hyperplasia and hypertrophy. Diabetes 49[C]:1341-1346, PMID: 10923635, https://doi.org/10.2337/diabetes.49.8.1341.
Movassat J, Saulnier C, Serradas P, Portha B. 1997. Impaired development of pancreatic beta-cell mass is a primary event during the progression to diabetes in the GK rat. Diabetologia 40(8):916-925, PMID: 9267986, https://doi.org/10.1007/ s001250050768.
Nahar MS, Soliman AS, Colacino JA, Calafat AM, Battige K, Hablas A, et al. 2012. Urinary bisphenol A concentrations in girls from rural and urban Egypt: a pilot study. Environ Health 11:20, PMID: 22472083, https://doi.org/10.1186/1476- 069X11-20.
Newbold RR, Padilla-Banks E, Jefferson WN, Heindel JJ. 2008. Effects of endocrine disruptors on obesity. Int J Androl 31(2):201-208, PMID: 18315718, https://doi.org/10.1111/j.1365-2605.2007.00858.x.
Ortsater H, Liss P, Akerman KE, Bergsten P. 2002. Contribution of glycolytic and mitochondrial pathways in glucose-induced changes in islet respiration and insulin secretion. Pflugers Arch 444(4):506-512, PMID: 12136270, https://doi.org/10. 1007/s00424-002-0842-9.
Patial S, Luo J, Porter KJ, Benovic JL, Parameswaran N. 2009. G-protein-coupled- receptor kinases mediate TNF[alpha]-induced NF[kappa]B signalling via direct interaction with and phosphorylation of I[kappa]B[alpha]. Biochem J 425(1):169-178, PMID: 19796012, https://doi.org/10.1042/BJ20090908.
Petrik J, Pell JM, Arany E, McDonald TJ, Dean WL, Reik W, et al. 1999. Overexpression of insulin-like growth factor-II in transgenic mice is associated with pancreatic islet cell hyperplasia. Endocrinology 140(5):2353-2363, PMID: 10218989, https://doi.org/10.1210/endo.140.5.6732.
Rahier J, Guiot Y, Goebbels RM, Sempoux C, Henquin JC. 2008. Pancreatic [beta]- cell mass in European subjects with type 2 diabetes. Diabetes Obes Metab 10:32-42, PMID: 18834431, https://doi.org/10.1111/j.1463-1326.2008.00969.x.
Reed JC, Paternostro G. 1999. Postmitochondrial regulation of apoptosis during heartfailure. Proc Natl Acad Sci U S A 96(14):7614-7616, PMID: 10393865.
Richardson SJ, Willcox A, Bone AJ, Foulis AK, Morgan NG. 2009. Islet-associated macrophages in type 2 diabetes. Diabetologia 52(8):1686-1688, PMID: 19504085, https://doi.org/10.1007/s00125-009-1410-z.
Rocha M, Diaz-Morales N, Rovira-Llopis S, Escribano-Lopez I, Banuls C, Hernandez-Mijares A, et al. 2016. Mitochondrial dysfunction and endoplasmic reticulum stress in diabetes. Curr Pharm Des 22(18):2640-2649, PMID: 26861650, https://doi.org/10.2174/1381612822666160209152033.
Sabanayagam C, Teppala S, Shankar A. 2013. Relationship between urinary bisphenol A levels and prediabetes among subjects free of diabetes. Acta Diabetol 50(4):625-631, PMID: 23636267, https://doi.org/10.1007/s00592-0130472- z.
Samy TS, Zheng R, Matsutani T, Rue LW III, Bland KI, Chaudry IH. 2003. Mechanism for normal splenic T lymphocyte functions in proestrus females after trauma: enhanced local synthesis of 17[beta]-estradiol. Am J Physiol Cell Physiol 285(1):C139-C149, PMID: 12660147, https://doi.org/10.1152/ajpcell.00058. 2003.
Sauter NS, Thienel C, Plutino Y, Kampe K, Dror E, Traub S, et al. 2015. Angiotensin II induces interleukin-1 b-mediated islet inflammation and [beta]-cell dysfunction independently of vasoconstrictive effects. Diabetes 64(4):1273-1283, PMID: 25352639, https://doi.org/10.2337/db14-1282.
Scaglia L, Cahill CJ, Finegood DT, Bonner-Weir S. 1997. Apoptosis participates in the remodeling of the endocrine pancreas in the neonatal rat. Endocrinology 138(4):1736-1741, PMID: 9075738, https://doi.org/10.1210/endo.138A5069.
Scaglia L, Smith FE, Bonner-Weir S. 1995. Apoptosis contributes to the involution of beta cell mass in the post partum rat pancreas. Endocrinology 136(12):5461- 5468, PMID: 7588296, https://doi.org/10.1210/endo.136.12.7588296.
Shankar A, Teppala S. 2011. Relationship between urinary bisphenol A levels and diabetes mellitus. J Clin Endocrinol Metab 96(12):3822-3826, PMID: 21956417, https://doi.org/10.1210/jc.2011-1682.
Sieli PT, Jasarevic E, Warzak DA, Mao J, Ellersieck MR, Liao C, et al. 2011. Comparison of serum bisphenol A concentrations in mice exposed to bisphenol A through the diet versus oral bolus exposure. Environ Health Perspect 119(9):1260-1265, PMID: 21642047, https://doi.org/10.1289/ehp.1003385.
Simmons RA, Suponitsky-Kroyter I, Selak MA. 2005. Progressive accumulation of mitochondrial DNA mutations and decline in mitochondrial function lead to [beta]-cell failure. J Biol Chem 280(31):28785-28791, PMID: 15946949, https://doi.org/ 10.1074/jbc.M505695200.
Simmons RA, Templeton LJ, Gertz SJ. 2001. Intrauterine growth retardation leads to the development of type 2 diabetes in the rat. Diabetes 50(10):2279-2286, PMID: 11574409, https://doi.org/10.2337/diabetes.50.10.2279.
Smithson G, Couse JF, Lubahn DB, Korach KS, Kincade PW. 1998. The role of estrogen receptors and androgen receptors in sex steroid regulation of B lymphopoiesis. J Immunol 161(1):27-34, PMID: 9647203.
Song L, Xia W, Zhou Z, Li Y, Lin Y, Wei J, et al. 2012. Low-level phenolic estrogen pollutants impair islet morphology and [beta]-cell function in isolated rat islets. J Endocrinol 215(2):303-311, PMID: 22946080, https://doi.org/10.1530/JOE-120219.
Soriano S, Alonso-Magdalena P, Garcia-Arevalo M, Novials A, Muhammed SJ, Salehi A, et al. 2012. Rapid insulinotropic action of low doses of bisphenol-A on mouse and human islets of Langerhans: role of estrogen receptor b. PLoS One 7(2):e31109, PMID: 22347437, https://doi.org/10.1371/journal.pone.0031109.
Soriano S, Ropero AB, Alonso-Magdalena P, Ripoll C, Quesada I, Gassner B, et al. 2009. Rapid regulation of K(ATP) channel activity by 17[beta]-estradiol in pancreatic [beta]-cells involves the estrogen receptor b and the atrial natriuretic peptide receptor. Mol Endocrinol 23(12):1973-1982, PMID: 19855088, https://doi.org/10. 1210/me.2009-0287.
Sorriento D, Fusco A, Ciccarelli M, Rungi A, Anastasio A, Carillo A, et al. 2013. Mitochondrial G protein coupled receptor kinase 2 regulates proinflammatory responses in macrophages. FEBS Lett 587(21):3487-3494, PMID: 24036448, https://doi.org/10.1016/j.febslet.2013.09.002.
Stahlhut RW, Welshons WV, Swan SH. 2009. Bisphenol A data in NHANES suggest longer than expected half-life, substantial nonfood exposure, or both. Environ Health Perspect 117(5):784-789, PMID: 19479022, https://doi.org/10.1289/ehp.0800376.
Sun Q, Cornelis MC, Townsend MK, Tobias DK, Eliassen AH, Franke AA, et al. 2014. Association of urinary concentrations of bisphenol A and phthalate metabolites with risk of type 2 diabetes: a prospective investigation in the Nurses' Health Study (NHS) and NHSII cohorts. Environ Health Perspect 122(6):616-623, PMID: 24633239, https://doi.org/10.1289/ehp.1307201.
Susiarjo M, Sasson I, Mesaros C, Bartolomei MS. 2013. Bisphenol A exposure disrupts genomic imprinting in the mouse. PLoS Genet 9(4):e1003401, PMID: 23593014, https://doi.org/10.1371/journal.pgen.1003401.
Susiarjo M, Xin F, Bansal A, Stefaniak M, Li C, Simmons RA, et al. 2015. Bisphenol A exposure disrupts metabolic health across multiple generations in the mouse. Endocrinology 156(6):2049-2058, PMID: 25807043, https://doi.org/10.1210/en.2014- 2027.
Swenne I. 1983. Effects of aging on the regenerative capacity of the pancreatic Bcell of the rat. Diabetes 32(1):14-19, PMID: 6336699, https://doi.org/10.2337/ diab.32.1.14.
Takahashi O, Oishi S. 2000. Disposition of orally administered 2,2-bis(4- hydroxyphenyl) propane (bisphenol A) in pregnant rats and the placental transfer to fetuses. Environ Health Perspect 108(10):931-935, PMID: 11049811, https://doi.org/10.1289/ ehp.00108931.
Thomas MK, Rastalsky N, Lee JH, Habener JF. 2000. Hedgehog signaling regulation of insulin production by pancreatic beta-cells. Diabetes 49(12):2039-2047, PMID: 11118005, https://doi.org/10.2337/diabetes.49.12.2039.
Tiano JP, Delghingaro-Augusto V, Le May C, Liu S, Kaw MK, Khuder SS, et al. 2011. Estrogen receptor activation reduces lipid synthesis in pancreatic islets and prevents [beta] cell failure in rodent models of type 2 diabetes. J Clin Invest 121(8):3331-3342, PMID: 21747171, https://doi.org/10.1172/JCI44564.
Turner MD, Nedjai B, Hurst T, Pennington DJ. 2014. Cytokines and chemokines: at the crossroads of cell signalling and inflammatory disease. Biochim Biophys Acta 1843(11):2563-2582, PMID: 24892271, https://doi.org/10.1016/j.bbamcr.2014.05.014.
van Esterik JC, Dolle ME, Lamoree MH, van Leeuwen SP, Hamers T, Legler J, et al. 2014. Programming of metabolic effects in C57BL/6JxFVB mice by exposure to bisphenol A during gestation and lactation. Toxicology 321:40-52, PMID: 24726836, https://doi.org/10.1016/j.tox.2014.04.001.
Vandenberg LN, Colborn T, Hayes TB, Heindel JJ, Jacobs DR Jr, Lee DH, et al. 2012. Hormones and endocrine-disrupting chemicals: low-dose effects and nonmonotonic dose responses. Endocr Rev 33(3):378-455, PMID: 22419778, https://doi.org/10.1210/er.2011-1050.
Vozza A, Parisi G, De Leonardis F, Lasorsa FM, Castegna A, Amorese D, et al. 2014. UCP2 transports C4 metabolites out of mitochondria, regulating glucose and glutamine oxidation. Proc Natl Acad Sci U S A 111(3):960-965, PMID: 24395786, https://doi.org/10.1073/pnas.1317400111.
Wei J, Lin Y, Li Y, Ying C, Chen J, Song L, et al. 2011. Perinatal exposure to bisphenol A at reference dose predisposes offspring to metabolic syndrome in adult rats on a high-fat diet. Endocrinology 152(8):3049-3061, PMID: 21586551, https://doi.org/10.1210/en.2011-0045.
Xu XB, He Y, Song C, Ke X, Fan SJ, Peng WJ, et al. 2014. Bisphenol A regulates the estrogen receptor alpha signaling in developing hippocampus of male rats through estrogen receptor. Hippocampus 24(12):1570-1580, PMID: 25074486, https://doi.org/10.1002/hipo.22336.
Yao R, Cooper GM. 1995. Requirement for phosphatidylinositol-3 kinase in the prevention of apoptosis by nerve growth factor. Science 267(5206):2003-2006, PMID: 7701324, https://doi.org/10.1126/science.7701324.
Yoon KH, Ko SH, Cho JH, Lee JM, Ahn YB, Song KH, et al. 2003. Selective [beta]- cell loss and a-cell expansion in patients with type 2 diabetes mellitus in Korea. J Clin Endocrinol Metab 88(5):2300-2308, PMID: 12727989, https://doi.org/10.1210/ jc.2002-020735.
Yoon K, Kwack SJ, Kim HS, Lee BM. 2014. Estrogenic endocrine-disrupting chemicals: molecular mechanisms of actions on putative human diseases. J Toxicol Environ Health B Crit Rev 17(3):127-174, PMID: 24749480, https://doi.org/10. 1080/10937404.2014.882194.
Yuchi Y, Cai Y, Legein B, De Groef S, Leuckx G, Coppens V, et al. 2015. Estrogen receptor a regulates [beta]-cell formation during pancreas development and following injury. Diabetes 64(9):3218-3228, PMID: 26015547, https://doi.org/10.2337/db14- 1798.
Zalko D, Soto AM, Dolo L, Dorio C, Rathahao E, Debrauwer L, et al. 2003. Biotransformations of bisphenol A in a mammalian model: answers and new questions raised by low-dose metabolic fate studies in pregnant CD1 mice. Environ Health Perspect 111(3):309-319, PMID: 12611660, https://doi.org/10. 1289/ehp.5603.
Zeidler A, Shargill NS, Chan TM. 1991. Peripheral insulin insensitivity in the hyperglycemic athymic nude mouse: similarity to noninsulin-dependent diabetes mellitus. Proc Soc Exp Biol Med 196(4):457-460, PMID: 2008443, https://doi.org/ 10.3181/00379727-196-43216.
Zhang Z, Alomirah H, Cho HS, Li YF, Liao C, Minh TB, et al. 2011. Urinary bisphenol A concentrations and their implications for human exposure in several Asian countries. Environ Sci Technol 45(16):7044-7050, PMID: 21732633, https://doi.org/ 10.1021/es200976k.
Zhang Z, Tsukikawa M, Peng M, Polyak E, Nakamaru-Ogiso E, Ostrovsky J, et al. 2013. Primary respiratory chain disease causes tissue-specific dysregulation of the global transcriptome and nutrient-sensing signaling network. PLoS One 8(7):e69282, PMID: 23894440, https://doi.org/10.1371/journal.pone.0069282.
Zielke HR, Huang Y, Baab PJ, Collins RM Jr, Zielke CL, Tildon JT. 1997. Effect of [alpha]-ketoisocaproate and leucine on the in vivo oxidation of glutamate and glutamine in the rat brain. Neurochem Res 22(9):1159-1164, PMID: 9251107.
Amita Bansal, [1,2,3] Cetewayo Rashid, [1,2,3] Frances Xin, [2,4] Changhong Li,  Erzsebet Polyak,  Anna Duemler, [1,7] Tom van der Meer, [1,8] Martha Stefaniak,  Sana Wajid, 9] Nicolai Doliba,  Marisa S. Bartolomei, [1,2,4] and Rebecca A. Simmons [1,2,3]
 Center for Research on Reproduction and Women's Health, Perelman School of Medicine, University of Pennsylvania, Philadelphia, Pennsylvania, USA
 Center of Excellence in Environmental Toxicology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, Pennsylvania, USA
 Division of Neonatology, Children's Hospital of Philadelphia, Philadelphia, Pennsylvania, USA
 Epigenetics Institute, Department of Cell and Developmental Biology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, Pennsylvania, USA
 Division of Endocrinology and Diabetes, Department of Pediatrics, Children's Hospital of Philadelphia, Philadelphia, Pennsylvania, USA
 Division of Human Genetics, Department of Pediatrics, Children's Hospital of Philadelphia, Philadelphia, Pennsylvania, USA
 Eberly College of Science, Pennsylvania State University, University Park, Pennsylvania, USA
 Department of Pediatrics, University of Groningen, Groningen, Netherlands
 Exposure Biology Informatics Core, Perelman School of Medicine, University of Pennsylvania, Philadelphia, Pennsylvania, USA
 Department of Biochemistry and Biophysics, Perelman School of Medicine, University of Pennsylvania, Philadelphia, Pennsylvania, USA
Please address correspondence to R.A. Simmons, Center for Research on Reproduction and Women' s Health, Perelman School of Medicine, University of Pennsylvania, 1308 BRB II/III, 421 Curie Blvd., Philadelphia, PA 19104 USA. Telephone: (215) 662-3269. Email: firstname.lastname@example.org; or M.S. Bartolomei, Center for Research on Reproduction and Women's Health, Perelman School of Medicine, University of Pennsylvania, 1308 BRB II/III, 421 Curie Blvd., Philadelphia, PA 19104 USA. Telephone: (215) 898-9063. Email: email@example.com
Supplemental Material is available online (https://doi.org/10.1289/EHP1674).
The authors declare they have no actual or potential competing financial interests.
Received 25 January 2017; Revised 26 July 2017; Accepted 26 July 2017; Published 27 September 2017.
Note to readers with disabilities: EHP strives to ensure that all journal content is accessible to all readers. However, some figures and Supplemental Material published in EHP articles may not conform to 508 standards due to the complexity of the information being presented. If you need assistance accessing journal content, please contact firstname.lastname@example.org. Our staff will work with you to assess and meet your accessibility needs within 3 working days.
Caption: Figure 1. Mitochondrial driven insulin secretion in islets of F1 and F2 adult male offspring: (A, C) F1 and (B, D) F2 adult male offspring. The maximal KIC-stimulated insulin release was calculated as the insulin area under the curve for 15 mM KIC. Data are individual litter data (one animal per litter), presented as mean + SEM, and analyzed using Dunnett's test. p-Values are relative to Control. Note: KCL, potassium chloride; KIC, [alpha]-ketoisocaproate.
Caption: Figure 2. Glucose-stimulated and mitochondrial associated insulin secretion in islets of F1 and F2 adult female offspring: (A-B, E-F) F1, and (C-D, G-H) F2 adult female offspring. The maximal glucose- and KIC-stimulated insulin release was calculated as the insulin area under the curve for 25 mM glucose or 15 mM KIC, respectively. Data are individual litter data (one animal per litter), presented as mean + SEM, and analyzed using Dunnett's test. p-Values are relative to Control. Note: KCL, potassium chloride; KIC, [alpha]-ketoisocaproate.
Caption: Figure 3. Basal and maximal mitochondrial oxygen consumption in islets of F1 and F2 adult male offspring: (A-B) F1, and (C-D) F2 adult male offspring. Data are individual litter data (islets pooled from 2-3 males from the same litter) with mean superimposed. Data were analyzed using Dunnett's test. p-values are relative to Control. Oxygen flow per cell at: ROUTINE, i.e., basal state of cell respiration; electron transport system [ETS; uncoupled; carbonyl cyanide-p-trifluoro-methoxyphenylhydrazone (FCCP)]; and residual oxygen consumption (rox; rotenone, antimycin A).
Caption: Figure 4. Beta cell mass of F1 and F2 adult male offspring. Representative photomicrographs of pancreatic immunofluorescent staining in (A-C) F1 and (D-F) F2 Control, LowerB, UpperB male mice. All images have insulin (red), glucagon (green), somatostatin (yellow), and DAPI (blue). Image magnified 20 x. Data are individual litter data (one animal per litter) with mean superimposed from (G) F1 and (H) F2 males, analyzed by Dunnett's test. p-Values are relative to Control.
Caption: Figure 5. Cell viability in islets of (A-C) F1 and (D-F) F2 male offspring. Data are individual litter data (one animal per litter) with mean superimposed. Each parameter is normalized to total protein concentration. Data were analyzed using Dunnett's test performed on log-transformed data, where required (BCL2). Note: RFU: relative fluorescence unit. p-Values are relative to Control.
Caption: Figure 6. IL6, MCP1 levels in (A-B) F1 and (C-D) F2 adult male offspring. Note: IL6, interleukin 6; MCPm: monocyte chemoattractant protein-1. Cytokine levels for each sample were normalized to total protein concentration. Data are individual litter data (one animal per litter) with mean superimposed. Data were analyzed by Dunnett's test performed on log-transformed data, where required (F1 IL6, MCP1; F2 IL6). p-Values are relative to Control.
Caption: Figure 7. Representative photomicrographs of pancreatic immunohistochemical staining in F1 and F2 adult male offspring: (A-H) F1 and (I-P) F2 Control, LowerB, UpperB mice, and Positive Ctrl (internal Control--lymph node). Image magnified 20 x. Note: CD3, cluster of differentiation 3 (marker for T lymphocytes); F4/80, marker for macrophages. F1, n = 6 litters per group (one animal per litter); F2, n = 5 litters Control and UpperB, n = 4 litters LowerB (one animal per litter).
Caption: Figure 8. Percent DNA methylation at Igf2 DMR1 and Esr1 Exon A in islets of F1 and F2 adult male offspring: (A-B) F1 and (C-D) F2 adult male offspring. Data are percent CpG methylation values from individual litter (islets pooled from 2-3 males from the same litter) and presented as mean + SEM. Data were analyzed using Dunnett's test; p-Values are relative to Control. ** p < 0.05, * p > 0.05 and p < 0.09.
Caption: Figure 9. RNA Seq data in islets of F1 male mice. (A) F1 Lower BPA male vs. Control male canonical pathways from Ingenuity Pathway Analysis (IPA); IPA parameters: 962 genes; p <0. 05; inflammatory pathways highlighted by the rectangular box. (B) F1 Upper BPA male vs. Control male canonical Pathways from IPA; IPA parameters: 1,275 genes; p < 0. 05; inflammatory and mitochondrial dysfunction pathways highlighted by the rectangular box. (C) Venn diagram illustrating number of differentially expressed genes when comparing F1 Lower BPA (LowerB), Upper BPA (UpperB), and Control males; parameters: false discovery rate (FDR) <10%, which is equivalent to p adjusted (also known as q- value) <0.1; n = 3 litters (islets pooled from 2-3 males from the same litter) per group.
Table 1. mRNA Levels in islets of F1 and F2 Lower BPA and Upper BPA male mice. Gene Lower BPA (95% CI) Upper BPA (95% CI) Igf2 F1 1.67 (1.20, 2.31) * 2.19 (1.71, 2.79) * F2 1.48 (1.03, 2.30) * 1.23 (1.13, 1.34) * Ucp2 F1 1.51 (1.03, 2.21) * 2.54 (1.52, 3.24) * F2 2.05 (1.78, 2.36) * 1.74 (1.05, 3.08) * Ogdh F1 0.97 (0.72, 1.29) 1.40 (0.96, 2.03) F2 1.20 (1.18, 1.22) * 1.32 (0.97, 1.80) Esr1 F1 0.70 (0.53, 0.93) * 0.88 (0.82, 0.94) * F2 0.38 (0.18, 0.84) * 0.48 (0.24, 0.97) * Pdx1 F1 1.32 (0.82, 2.11) 1.72 (1.11, 2.64) * F2 1.31 (0.86, 2.00) 1.21 (0.48, 3.03) Igf1 F1 1.05 (0.95, 1.15) 1.55 (1.10, 2.19) * F2 1.08 (0.94, 1.25) 1.23 (1.13, 1.34) * Hnfl[alpha] F1 1.06 (0.89, 1.27) 1.43 (1.30, 1.59) * F2 1.00 (0.85, 1.17) 1.21 (0.94, 1.57) Kcnj11 F1 1.54 (0.70, 3.41) 0.94 (0.45, 1.99) F2 1.56 (0.92, 2.64) 1.06 (0.36, 3.12) Abcc8 F1 1.48 (0.75, 2.92) 0.99 (0.55, 1.78) F2 1.94 (1.48, 2.54) * 1.46 (0.63, 3.37) Snap25 F1 1.02 (0.43, 2.38) 0.71 (0.30, 1.65) F2 1.46 (0.66, 3.26) 1.17 (0.38, 3.61) Note: Data were normalized to Rpl19 and Beta actin or Ppia mRNA levels and expressed as fold change relative to Controls (where Control is set at 1.0) with 95% confidence intervals (CI) in islets of F1 and F2 Lower BPA (10 [micro]g-kg-d) and Upper BPA (10 mg-kg-d) male mice. n = 5 litters per group; islets pooled from 2-3 male mice per litter. * p < 0.05 vs. Controls. Table 2. Cytokine/chemokine levels on a LUMINEX assay in pancreatic lysates of F1 and F2 adult male offspring. Cytokine/ Inflammatory Produced by Chemokine action IL4 pro/anti Macrophages, eosinophils, basophills IL5 pro/anti [T.sub.h]2 cells, mast Cells JL6 pro/anti T-cells, macrophages IL17 pro [T.sub.h] cells MCP1 pro Chemoattractant produced by injured tissue, or infection TNF[alpha] pro Macrophages [CD4.sup.+], natural killer and mast cells IFN[gamma] pro Natural killer cells, natural killer T-cells, macrophages, [CD4.sup.+], [CD8.sup.+] Cytokine/ Target cell/action Chemokine IL4 Macrophages, T cells, [T.sub.h]0 [right arrow] [T.sub.h]2 IL5 B-cell growth, [up arrow] eosinophills JL6 [up arrow] B-cells, naive T cells [right arrow] [T.sub.h]17, [down arrow] [T.sub.reg] IL17 Recruits macrophages and neutrophils MCP1 [up arrow] Macrophages, memory T cells, dendritic cells, basophils TNF[alpha] Macrophages, IL1- and IL6-producing cells IFN[gamma] Macrophages, [up arrow] MHCII Cytokine/ F1 Offspring Chemokine Control LowerB (pg/[micro]g) (pg/[micro]g) IL4 3E-4 [+ or -] 6E-5 2E-4 [+ or -] 7E-5 IL5 4E-4 [+ or -] 3E-5 3E-4 [+ or -] 3E-5 JL6 6E-4 [+ or -] 5E-4 7E-4 [+ or -] 5E-4 IL17 2E-4 [+ or -] 7E-5 4E-4 [+ or -] 8E-5 MCP1 2E-3 [+ or -] 4E-4 3E-3 [+ or -] 5E-4 TNF[alpha] 4E-4 [+ or -] 9E-5 5E-4 [+ or -] 1E-4 IFN[gamma] 2E-3 [+ or -] IE-3 4E-3 [+ or -] IE-3 Cytokine/ F1 Offspring Chemokine UpperB LowerB UpperB (pg/[micro]g) (p-value) (p-value) IL4 3E-4 [+ or -] 7E-5 0.60 0.75 IL5 7E-4 [+ or -] 4E-5 0.08 0.01 JL6 22E-4 [+ or -] 5E-4 0.89 0.08 IL17 4E-4 [+ or -] 8E-5 0.20 0.33 MCP1 4E-3 [+ or -] 5E-4 0.18 0.02 TNF[alpha] 5E-4 [+ or -] 1E-4 0.59 0.53 IFN[gamma] 4E-3 [+ or -] IE-3 0.45 0.59 Cytokine/ F2 Offspring Chemokine Control LowerB (pg/[micro]g) (pg/[micro]g) IL4 3E-5 [+ or -] 5E-5 lE-4 [+ or -] 5E-5 IL5 2E-4 [+ or -] 4E-4 4E-4 [+ or -] 4E-4 JL6 3E-4 [+ or -] 6E-4 12E-4 [+ or -] 6E-4 IL17 lE-4 [+ or -] 9E-5 2E-4 [+ or -] 9E-5 MCP1 2E-3 [+ or -] 21E-4 7E-3 [+ or -] 21E-4 TNF[alpha] 2E-4 [+ or -] 2E-4 3E-4 [+ or -] 2E-4 IFN[gamma] 2E-3 [+ or -] 6E-3 2E-3 [+ or -] 6E-3 Cytokine/ F2 Offspring Chemokine UpperB LowerB UpperB (pg/[micro]g) (p-value) (p-value) IL4 2E-4 [+ or -] 5E-5 0.23 0.04 IL5 13E-4 [+ or -] 4E-4 0.43 0.09 JL6 19E-4 [+ or -] 5E-4 0.15 0.04 IL17 4E-4 [+ or -] 8E-5 0.14 0.03 MCP1 8E-3 [+ or -] 17E-4 0.27 0.09 TNF[alpha] 4E-4 [+ or -] 1E-4 0.62 0.28 IFN[gamma] 9E-3 [+ or -] 5E-3 0.99 0.85 Note: Data are normalized to total protein concentration as pg of cytokine or chemokine per [micro]g of total protein, and presented as mean + SEM. Decimal places are represented by E notation, where E represents base 10, followed by the power of 10 (example, E-3 = [10.sup.-3], E-4 = [10.sup.-4], E-5 = [10.sup.-5]). Information in columns 2, 3, and 4 is derived from Cavaillon 2001; Commins et al. 2010; Kimura and Kishimoto 2010; Kourilsky and Truffa-Bachi 2001; Turner et al. 2014. Control, 7% corn o11 diet; LowerB, Lower BP A (10 [micro]g-kg-d); UpperB, Upper BP A (10 mg-kg-d). n = 3 to 4 litters per group, p-Values are from Dunnett's test performed on log-transformed data, where required (FI IL5, IL6, MCP1, and IFN[gamma]; F2 IL4, IL5 IL6, IL17, TNF[alpha], and IFN[gamma]). Table 3. Quantification of CD3 and F4/80 staining area in pancreatic sections of F1 and F2 adult male offspring. Staining area Control ([mm.sup.2]/[mm.sup.2]) Area of percent positive CD3 staining in islets/ total analysis area F1 14.88 [+ or -] 5.57 F2 5.02 [+ or -] 2.89 Area of percent positive CD3 staining in pancreas excluding islets and lymph nodes/total analysis area F1 20.14 [+ or -] 3.43 F2 2.28 [+ or -] 1.10 Area of percent positive CD3 staining in pancreatic lymph nodes/total analysis area F1 23.81 [+ or -] 6.16 F2 17.99 [+ or -] 6.19 Area of percent positive F4/80 staining in islets/ total analysis area F1 1.91 [+ or -] 1.01 F2 3.52 [+ or -] 1.10 Area of percent positive F4/80 staining in pancreas excluding islets and lymph nodes/total analysis area F1 2.67 [+ or -] 1.62 F2 7.99 [+ or -] 1.74 Area of percent positive F4/80 staining in pancreatic lymph nodes/total analysis area F1 7.74 [+ or -] 3.34 F2 14.45 [+ or -] 4.17 Staining area LowerB ([mm.sup.2] /[mm.sup.2]) Area of percent positive CD3 staining in islets/ total analysis area F1 29.33 [+ or -] 5.57 F2 5.90 [+ or -] 3.23 Area of percent positive CD3 staining in pancreas excluding islets and lymph nodes/total analysis area F1 24.67 [+ or -] 3.43 F2 6.40 [+ or -] 1.23 Area of percent positive CD3 staining in pancreatic lymph nodes/total analysis area F1 42.85 [+ or -] 6.16 F2 23.37 [+ or -] 6.92 Area of percent positive F4/80 staining in islets/ total analysis area F1 4.49 [+ or -] 1.01 F2 6.62 [+ or -] 1.23 Area of percent positive F4/80 staining in pancreas excluding islets and lymph nodes/total analysis area F1 7.51 [+ or -] 1.62 F2 12.87 [+ or -] 1.95 Area of percent positive F4/80 staining in pancreatic lymph nodes/total analysis area F1 11.81 [+ or -] 3.34 F2 12.23 [+ or -] 4.66 Staining area UpperB LowerB vs. Control ([mm.sup.2]/[mm.sup.2]) (p-value) Area of percent positive CD3 staining in islets/ total analysis area F1 35.46 [+ or -] 5.57 0.04 F2 12.21 [+ or -] 2.89 0.51 Area of percent positive CD3 staining in pancreas excluding islets and lymph nodes/total analysis area F1 32.52 [+ or -] 3.43 0.42 F2 4.69 [+ or -] 1.10 0.04 Area of percent positive CD3 staining in pancreatic lymph nodes/total analysis area F1 44.98 [+ or -] 6.16 0.08 F2 15.47 [+ or -] 6.19 0.78 Area of percent positive F4/80 staining in islets/ total analysis area F1 5.94 [+ or -] 1.01 0.12 F2 7.66 [+ or -] 1.10 0.15 Area of percent positive F4/80 staining in pancreas excluding islets and lymph nodes/total analysis area F1 10.48 [+ or -] 1.62 0.09 F2 15.22 [+ or -] 1.74 0.10 Area of percent positive F4/80 staining in pancreatic lymph nodes/total analysis area F1 17.46 [+ or -] 3.34 0.60 F2 11.21 [+ or -] 4.17 0.94 Staining area UpperB vs Control (p-value) Area of percent positive CD3 staining in islets/ total analysis area F1 0.02 F2 0.06 Area of percent positive CD3 staining in pancreas excluding islets and lymph nodes/total analysis area F1 0.04 F2 0.18 Area of percent positive CD3 staining in pancreatic lymph nodes/total analysis area F1 0.05 F2 0.84 Area of percent positive F4/80 staining in islets/ total analysis area F1 0.02 F2 0.04 Area of percent positive F4/80 staining in pancreas excluding islets and lymph nodes/total analysis area F1 0.01 F2 0.03 Area of percent positive F4/80 staining in pancreatic lymph nodes/total analysis area F1 0.10 F2 0.81 Note: Data are normalized to total analysis area and presented as mean [+ or -] SEM. Control, 7% corn oil diet; LowerB, Lower BPA (10 [micro]g/kg/d); UpperB, Upper BPA (10 mg/kg/d). n = 4-6 litters per group (one animal per litter). p/Values are from Dunnett's test performed on log-transformed data, where required (F1 CD3 islets, F1 CD3 excluding islets and lymph nodes, F2 CD3 islets, F2 F4/80 lymph nodes). Table 4. Genes with q <0.1 in F1 Lower BPA male vs. Control male from RNA seq analysis. The genes that were assessed via qPCR are included in the table (rows 10-20) for completeness. Serial no. Gene name logFC logCPM p-Value QValue.BH 1 Lars2 -2.68 9.43 1.49E-06 0.04 2 Cth -2.81 2.43 1.60E-06 0.04 3 Klk1b22 -12.35 -0.04 8.85E-06 0.08 4 Cyp2c70 -3.67 -2.71 1.32E-05 0.08 5 Cpa5 -10.75 5.70 1.41E-05 0.08 6 Klk1b24 -5.32 1.56 1.65E-05 0.08 7 Tinag -2.48 1.01 2.22E-05 0.09 8 Reg3a -4.02 3.38 2.83E-05 0.09 9 Fgf21 -3.78 0.34 2.88E-05 0.09 10 Igf2 4.27 1.26 0.02 0.80 11 Esr1 -2.77 1.56 0.18 0.80 12 Esr2 -2.76 0.87 0.18 0.80 13 Snap25 -1.63 1.58 0.31 0.98 14 Kcnj11 -0.60 6.35 0.31 0.98 15 Ogdh 0.31 2.77 0.70 1.00 16 Abcc8 0.18 5.18 0.77 1.00 17 Pdx1 0.25 2.62 0.78 1.00 18 Igf1 0.18 2.41 0.83 1.00 19 Ucp2 -0.34 1.15 0.84 1.00 20 Hnf1a -0.09 6.06 0.86 1.00 Note: CPM, counts per million; FC, fold change. n = 3 litters per group (islets pooled from 2-3 male mice per litter). Table 5. Genes with q <0.1 in F1 Upper BPA male vs. Control male from RNA seq analysis. The genes that were assessed via qPCR are included in the table (rows 42-52) for completeness. Serial no. Gene name logFC logCPM p-Value QValue.BH 1 Lars2 -4.65 9.43 1.99E-14 4.87E-10 2 Rn45s -2.84 13.57 7.84E-10 9.59E-06 3 Dbt -3.6 3.19 3.93E-07 0.003 4 Klk1b5 -4.64 2.25 5.21E-06 0.03 5 Matn4 -1.92 9.06 7.90E-06 0.004 6 Ctnnbip1 -1.83 6.67 1.75E-05 0.007 7 Sdhc -2.06 4.94 6.00E-05 0.02 8 Spryd4 -2.34 4.13 6.11E-05 0.02 9 Chst2 -2.15 5.92 8.48E-05 0.02 10 Nccrp1 -2.27 4.16 0.0002 0.03 11 2610203C22Rik -1.83 4.51 0.0003 0.04 12 Gm5148 -1.53 10.35 0.0002 0.04 13 Zfhx2os 1.45 6.56 0.0003 0.05 14 Fam69a -1.54 8.97 0.0003 0.05 15 0610039K10Rik -1.37 9.10 0.0006 0.06 16 9530027J09Rik -0.87 7.67 0.0008 0.06 17 Casc1 -1.79 4.96 0.0005 0.06 18 Ccdc176 -1.48 7.85 0.0007 0.06 19 Ifi27l2a -1.51 4.68 0.0006 0.06 20 Loxl3 -1.16 7.96 0.0007 0.06 21 Lrrc75b -1.77 4.31 0.0007 0.06 22 Pigf -0.86 8.57 0.0007 0.06 23 4930429F24Rik 1.29 5.68 0.0009 0.06 24 Nphs1os -1.72 7.69 0.0009 0.06 25 1700028J19Rik -1.70 4.77 0.001 0.07 26 A230070E04Rik 1.45 5.51 0.001 0.07 27 Adam17 -1.30 7.30 0.002 0.08 28 Arhgap26 -1.57 4.72 0.002 0.08 29 C330022C24Rik -1.67 5.73 0.002 0.08 30 Cand2 -1.52 11.15 0.001 0.08 31 Casp9 -2.02 18.49 0.002 0.08 32 Snx22 -1.57 12.90 0.002 0.08 33 Spg11 -0.89 6.16 0.001 0.08 34 Tap2 -1.68 5.53 0.002 0.08 35 Vmn2r9 -1.62 6.12 0.002 0.08 36 Gm5088 -1.55 10.41 0.002 0.08 37 Gtf2i -2.10 4.38 0.002 0.08 38 Txnip -1.25 5.21 0.002 0.08 39 Prss57 -0.92 6.75 0.002 0.09 40 Tspan4 -0.88 6.53 0.002 0.09 41 Slc12a4 -1.78 4.98 0.002 0.09 42 Esr1 1.62 1.56 0.20 0.94 43 Abcc8 0.78 5.18 0.19 0.94 44 Igf2 1.24 1.26 0.28 0.94 45 Esr2 -1.50 0.87 0.35 0.94 46 Pdx1 0.67 2.62 0.39 0.97 47 Igf1 -0.47 2.41 0.56 1 48 Hnf1a -0.19 6.06 0.71 1 49 Snap25 0.34 1.58 0.76 1 50 Ogdh 0.16 2.77 0.83 1 51 Kcnj11 0.07 6.35 0.90 1 52 Ucp2 0.07 1.15 0.97 1 Note: CPM, counts per million; FC, fold change. n = 3 litters per group (islets pooled from 2-3 male mice per litter). Table 6. Top 25 genes, by p-values, in F1 Lower BPA female vs. Control female from RNA seq analysis. Serial no. Gene name logFC logCPM p-Value QValue.BH 1 Creb3l3 -5.83 2.61 9.51E-05 1 2 4833419F23Rik 5.26 2.44 2.03E-04 1 3 Gm10012 -6.86 2.66 3.88E-04 1 4 Mpzl2 5.01 1.96 3.88E-04 1 5 Fam204a 5.13 2.23 7.35E-04 1 6 Krt19 -5.48 2.51 7.40E-04 1 7 Nedd9 -5.07 2.14 7.86E-04 1 8 4930470O06Rik -5.65 2.03 8.05E-04 1 9 Egr1 -4.95 2.23 9.15E-04 1 10 1700019G17Rik -5.28 2.32 0.001 1 11 Trmt5 4.89 2.24 0.001 1 12 Espl1 -5.20 2.04 0.001 1 13 St18 -5.18 2.65 0.001 1 14 Tspan2os -4.07 2.74 0.001 1 15 Drg1 -4.89 2.25 0.001 1 16 Ggta1 5.51 2.06 0.001 1 17 Fbxw7 -4.84 2.28 0.001 1 18 Serf2 -4.85 2.29 0.001 1 19 4930509J09Rik -5.05 2.40 0.001 1 20 Dhrs1 -5.02 2.21 0.001 1 21 8430437L04Rik -5.29 1.96 0.001 1 22 Rragb 4.71 2.38 0.002 1 23 Sh3bgr 4.66 1.80 0.002 1 24 Taf11 -4.86 2.12 0.002 1 25 4931431B13Rik -5.67 4.75 0.002 1 Note: CPM, counts per million; FC, fold change. n = 3 litters per group (islets pooled from 2-3 female mice per litter). Table 7. Top 25 genes, by p-values, in F1 Upper BPA female vs. Control female from RNA seq analysis. Serial no. Gene name logFC logCPM p-Value QValue.BH 1 Fabp1 -9.26 5.54 4.93E-05 0.83 2 Spr 5.41 2.14 8.67E-05 0.83 3 Gda -6.75 3.71 1.82E-04 0.83 4 BC051537 3.19 5.86 2.27E-04 0.83 5 St3gal4 -5.99 3.01 2.86E-04 0.83 6 Rragb 5.26 2.38 3.79E-04 0.83 7 Wdr74 4.98 2.01 4.44E-04 0.83 8 Myo1a -7.73 4.63 4.77E-04 0.83 9 Rprm 5.64 2.26 5.47E-04 0.83 10 Pcsk5 -5.61 2.50 0.0006 0.83 11 Psmg3 -5.78 3.18 0.0006 0.83 12 Smim24 -4.89 5.48 0.0007 0.83 13 Apoa1 -6.47 5.39 0.0007 0.83 14 Rab40b 5.35 2.17 0.0007 0.83 15 Rbp2 -6.40 4.70 0.0008 0.83 16 Ephx2 -6.12 3.12 0.0009 0.83 17 Ccs -5.09 2.77 0.0009 0.83 18 Ano3 -5.57 2.56 0.0009 0.83 19 Tctn1 2.22 8.44 0.0009 0.83 20 4833419F23Rik 4.78 2.44 0.001 0.83 21 Adipor2 -3.01 3.21 0.001 0.83 22 Enpp3 -5.20 2.33 0.001 0.83 23 Slc15a1 -6.30 3.33 0.001 0.83 24 Pkdrej -3.77 2.69 0.001 0.83 25 Sis -5.91 6.32 0.001 0.83 Note: CPM, counts per million; FC, fold change. n = 3 litters per group (islets pooled from 2-3 female mice per litter). Table 8. Dose-specific effects of bisphenol A on pancreatic islets of male offspring. Group Insulin Cytosolic content calcium F1 GENERATION LowerB [down arrow] NS insulin content (#) UpperB NS NS F2 GENERATION LowerB [down arrow] NS insulin content (#) UpperB NS NS Group Mitochondria KIC ramp Oroboros mRNA/protein levels F1 GENERATION LowerB NS NS [up arrow] Ucp2 [up arrow] UCP2 protein levels NS UpperB [down arrow] [down arrow] [up arrow] Ucp2 [up KIC-stimulated Basal arrow] UCP2 protein insulin and maximal levels NS release respiration F2 GENERATION LowerB NS NS [up arrow] Ucp2, [up arrow] Ogdh [up arrow] UCP2 protein levels NS UpperB [down arrow] [down arrow] [up arrow] Ucp2 [up KIC-stimulated Basal arrow] UCP2 protein insulin and maximal levels NS release respiration Group Beta cell mass, and islet viability Histology Viability assay F1 GENERATION LowerB [down arrow] [up arrow] Cell death [beta]-cell mass UpperB NS [up arrow] Akt-phos, [up arrow] Bcl2 (#) F2 GENERATION LowerB [down arrow] [up arrow] Cell death (#) [beta]-cell mass UpperB NS [up arrow] Akt-phos Group Gene expression and methylation mRNA/protein Igf2/Esr1 levels methylation F1 GENERATION LowerB [up arrow] Igf2, [up arrow] Esr1 (#) [down arrow] Esr1 [down arrow] ESR1 (#) UpperB [up arrow] Igf2, [up arrow] Igf2 [down arrow] Esr1 ESR1 NS F2 GENERATION LowerB [up arrow] Igf2, NS [down arrow] Esr1 [down arrow] ESR1 (#) UpperB [up arrow] Igf2, [up arrow] Igf2 [up [down arrow] Esr1 arrow] Esr1 (#) ESR1 NS Group Inflammation Luminex ELISA F1 GENERATION LowerB [up arrow] IL5 (#) [up arrow] IL6 (#) [up arrow] MCP1 NS UpperB [up arrow] IL5, [up [up arrow] IL6, NS arrow] IL6 (#), [up [up arrow] MCP1 arrow] MCP1, [up arrow] G-CSF, [up arrow] IL9 (#) F2 GENERATION LowerB [up arrow] IL17 [up arrow] IL6 (#) [up arrow] MCP1 UpperB [down arrow] IL4, [up arrow] IL6, NS [up arrow] IL5 (#), [up arrow] MCP1 [up arrow] IL6, [up arrow] IL17, [up arrow] MCP1 (#) Group Inflammation Histology F1 GENERATION LowerB [up arrow] CD3, [up arrow] F4/80 UpperB [up arrow] CD3, [up arrow] F4/80 F2 GENERATION LowerB [up arrow] CD3, [up arrow] F4/ 80 (#) UpperB [up arrow] CD3, [up arrow] F4/ 80 Note: Arrows indicate levels are statistically (p < 0.05) increased or decreased relative to Controls. LowerB, Lower BPA (10 [micro]g/ kg/d); NS, not significant (p > 0.1); UpperB, Upper BPA (10 mg/kg/ d). (#) p< 0.05 and p < 0.1
|Printer friendly Cite/link Email Feedback|
|Author:||Bansal, Amita; Rashid, Cetewayo; Xin, Frances; Li, Changhong; Polyak, Erzsebet; Duemler, Anna; van d|
|Publication:||Environmental Health Perspectives|
|Date:||Sep 1, 2017|
|Previous Article:||Children's Lead Exposure: A Multimedia Modeling Analysis to Guide Public Health Decision-Making.|
|Next Article:||Respiratory, Dermal, and Eye Irritation Symptoms Associated with Corexit[TM] EC9527A/EC9500A following the Deepwater Horizon Oil Spill: Findings from...|