Modeling PAH Mixture Interactions in a Human In Vitro Organotypic Respiratory Model.

One of the most significant challenges in human health risk assessment is to evaluate hazards from exposure to environmental chemical mixtures. Polycyclic aromatic hydrocarbons (PAHs) are a class of ubiquitous contaminants typically found as mixtures in gaseous and particulate phases in ambient air pollution associated with petrochemicals from Superfund sites and the burning of fossil fuels. However, little is understood about how PAHs in mixtures contribute to toxicity in lung cells. To investigate mixture interactions and component additivity from environmentally relevant PAHs, two synthetic mixtures were created from PAHs identified in passive air samplers at a legacy creosote site impacted by wildfires. The primary human bronchial epithelial cells differentiated at the air-liquid interface were treated with PAH mixtures at environmentally relevant proportions and evaluated for the differential expression of transcriptional biomarkers related to xenobiotic metabolism, oxidative stress response, barrier integrity, and DNA damage response. Component additivity was evaluated across all endpoints using two independent action (IA) models with and without the scaling of components by toxic equivalence factors. Both IA models exhibited trends that were unlike the observed mixture response and generally underestimated the toxicity across dose suggesting the potential for non-additive interactions of components. Overall, this study provides an example of the usefulness of mixture toxicity assessment with the currently available methods while demonstrating the need for more complex yet interpretable mixture response evaluation methods for environmental samples.

Benzo[a]pyrene toxicokinetics in humans following dietary supplementation with 3,3'-diindolylmethane (DIM) or Brussels sprouts.

Utilizing the atto-zeptomole sensitivity of UPLC-accelerator mass spectrometry (UPLC-AMS), we previously demonstrated significant first-pass metabolism following escalating (25-250 ng) oral micro-dosing in humans of [14C]-benzo[a]pyrene ([14C]-BaP). The present study examines the potential for supplementation with Brussels sprouts (BS) or 3,3'-diindolylmethane (DIM) to alter plasma levels of [14C]-BaP and metabolites over a 48-h period following micro-dosing with 50 ng (5.4 nCi) [14C]-BaP. Volunteers were dosed with [14C]-BaP following fourteen days on a cruciferous vegetable restricted diet, or the same diet supplemented for seven days with 50 g of BS or 300 mg of BR-DIM® prior to dosing. BS or DIM reduced total [14C] recovered from plasma by 56-67% relative to non-intervention. Dietary supplementation with DIM markedly increased Tmax and reduced Cmax for [14C]-BaP indicative of slower absorption. Both dietary treatments significantly reduced Cmax values of four downstream BaP metabolites, consistent with delaying BaP absorption. Dietary treatments also appeared to reduce the T1/2 and the plasma AUC(0,∞) for Unknown Metabolite C, indicating some effect in accelerating clearance of this metabolite. Toxicokinetic constants for other metabolites followed the pattern for [14C]-BaP (metabolite profiles remained relatively consistent) and non-compartmental analysis did not indicate other significant alterations. Significant amounts of metabolites in plasma were at the bay region of [14C]-BaP irrespective of treatment. Although the number of subjects and large interindividual variation are limitations of this study, it represents the first human trial showing dietary intervention altering toxicokinetics of a defined dose of a known human carcinogen.

3,3'-Diindolylmethane Exhibits Significant Metabolism after Oral Dosing in Humans.

3,3'-Diindolylmethane (DIM), a major phytochemical derived from ingestion of cruciferous vegetables, is also a dietary supplement. In preclinical models, DIM is an effective cancer chemopreventive agent and has been studied in a number of clinical trials. Previous pharmacokinetic studies in preclinical and clinical models have not reported DIM metabolites in plasma or urine after oral dosing, and the pharmacological actions of DIM on target tissues is assumed to be solely via the parent compound. Seven subjects (6 males and 1 female) ranging from 26-65 years of age, on a cruciferous vegetable-restricted diet prior to and during the study, took 2 BioResponse DIM 150-mg capsules (45.3 mg DIM/capsule) every evening for one week with a final dose the morning of the first blood draw. A complete time course was performed with plasma and urine collected over 48 hours and analyzed by UPLC-MS/MS. In addition to parent DIM, two monohydroxylated metabolites and 1 dihydroxylated metabolite, along with their sulfate and glucuronide conjugates, were present in both plasma and urine. Results reported here are indicative of significant phase 1 and phase 2 metabolism and differ from previous pharmacokinetic studies in rodents and humans, which reported only parent DIM present after oral administration. 3-((1H-indole-3-yl)methyl)indolin-2-one, identified as one of the monohydroxylated products, exhibited greater potency and efficacy as an aryl hydrocarbon receptor agonist when tested in a xenobiotic response element-luciferase reporter assay using Hepa1 cells. In addition to competitive phytochemical-drug adverse reactions, additional metabolites may exhibit pharmacological activity highlighting the importance of further characterization of DIM metabolism in humans. SIGNIFICANCE STATEMENT: 3,3'-Diindolylmethane (DIM), derived from indole-3-carbinol in cruciferous vegetables, is an effective cancer chemopreventive agent in preclinical models and a popular dietary supplement currently in clinical trials. Pharmacokinetic studies to date have found little or no metabolites of DIM in plasma or urine. In marked contrast, we demonstrate rapid appearance of mono- and dihydroxylated metabolites in human plasma and urine as well as their sulfate and glucuronide conjugates. The 3-((1H-indole-3-yl)methyl)indolin-2-one metabolite exhibited significant aryl hydrocarbon receptor agonist activity, emphasizing the need for further characterization of the pharmacological properties of DIM metabolites.

Linking Coregulated Gene Modules with Polycyclic Aromatic Hydrocarbon-Related Cancer Risk in the 3D Human Bronchial Epithelium.

Exposure to polycyclic aromatic hydrocarbons (PAHs) often occurs as complex chemical mixtures, which are linked to numerous adverse health outcomes in humans, with cancer as the greatest concern. The cancer risk associated with PAH exposures is commonly evaluated using the relative potency factor (RPF) approach, which estimates PAH mixture carcinogenic potential based on the sum of relative potency estimates of individual PAHs, compared to benzo[a]pyrene (BAP), a reference carcinogen. The present study evaluates molecular mechanisms related to PAH cancer risk through integration of transcriptomic and bioinformatic approaches in a 3D human bronchial epithelial cell model. Genes with significant differential expression from human bronchial epithelium exposed to PAHs were analyzed using a weighted gene coexpression network analysis (WGCNA) two-tiered approach: first to identify gene sets comodulated to RPF and second to link genes to a more comprehensive list of regulatory values, including inhalation-specific risk values. Over 3000 genes associated with processes of cell cycle regulation, inflammation, DNA damage, and cell adhesion processes were found to be comodulated with increasing RPF with pathways for cell cycle S phase and cytoskeleton actin identified as the most significantly enriched biological networks correlated to RPF. In addition, comodulated genes were linked to additional cancer-relevant risk values, including inhalation unit risks, oral cancer slope factors, and cancer hazard classifications from the World Health Organization's International Agency for Research on Cancer (IARC). These gene sets represent potential biomarkers that could be used to evaluate cancer risk associated with PAH mixtures. Among the values tested, RPF values and IARC categorizations shared the most similar responses in positively and negatively correlated gene modules. Together, we demonstrated a novel manner of integrating gene sets with chemical toxicity equivalence estimates through WGCNA to understand potential mechanisms.

Bioinformatics Resource Manager: a systems biology web tool for microRNA and omics data integration.

BACKGROUND: The Bioinformatics Resource Manager (BRM) is a web-based tool developed to facilitate identifier conversion and data integration for Homo sapiens (human), Mus musculus (mouse), Rattus norvegicus (rat), Danio rerio (zebrafish), and Macaca mulatta (macaque), as well as perform orthologous conversions among the supported species. In addition to providing a robust means of identifier conversion, BRM also incorporates a suite of microRNA (miRNA)-target databases upon which to query target genes or to perform reverse target lookups using gene identifiers. RESULTS: BRM has the capability to perform cross-species identifier lookups across common identifier types, directly integrate datasets across platform or species by performing identifier retrievals in the background, and retrieve miRNA targets from multiple databases simultaneously and integrate the resulting gene targets with experimental mRNA data. Here we use workflows provided in BRM to integrate RNA sequencing data across species to identify common biomarkers of exposure after treatment of human lung cells and zebrafish to benzo[a]pyrene (BAP). We further use the miRNA Target workflow to experimentally determine the role of miRNAs as regulators of BAP toxicity and identify the predicted functional consequences of miRNA-target regulation in our system. The output from BRM can easily and directly be uploaded to freely available visualization tools for further analysis. From these examples, we were able to identify an important role for several miRNAs as potential regulators of BAP toxicity in human lung cells associated with cell migration, cell communication, cell junction assembly and regulation of cell death. CONCLUSIONS: Overall, BRM provides bioinformatics tools to assist biologists having minimal programming skills with analysis and integration of high-content omics' data from various transcriptomic and proteomic platforms. BRM workflows were developed in Java and other open-source technologies and are served publicly using Apache Tomcat at https://cbb.pnnl.gov/brm/ .

Toxicokinetics of benzo[a]pyrene in humans: Extensive metabolism as determined by UPLC-accelerator mass spectrometry following oral micro-dosing.

Benzo[a]pyrene (BaP), is a known human carcinogen (International Agency for Research on Cancer (IARC) class 1). The remarkable sensitivity (zepto-attomole 14C in biological samples) of accelerator mass spectrometry (AMS) makes possible, with de minimus risk, pharmacokinetic (PK) analysis following [14C]-BaP micro-dosing of humans. A 46 ng (5 nCi) dose was given thrice to 5 volunteers with minimum 2 weeks between dosing and plasma collected over 72 h. [14C]-BaPeq PK analysis gave plasma Tmax and Cmax values of 1.25 h and 29-82 fg/mL, respectively. PK parameters were assessed by non- compartment and compartment models. Intervals between dosing ranged from 20 to 420 days and had little impact on intra-individual variation. DNA, extracted from peripheral blood mononuclear cells (PBMCs) of 4 volunteers, showed measurable levels (LOD ~ 0.5 adducts/1011 nucleotides) in two individuals 2-3 h post-dose, approximately three orders of magnitude lower than smokers or occupationally-exposed individuals. Little or no DNA binding was detectable at 48-72 h. In volunteers the allelic variants CYP1B1*1/*⁎1, *1/*3 or *3/*3 and GSTM1*0/0 or *1 had no impact on [14C]-BaPeq PK or DNA adduction with this very limited sample. Plasma metabolites over 72 h from two individuals (one CYP1B1*1/*1 and one CYP1B1*3/*3) were analyzed by UPLC-AMS. In both individuals, parent [14C]-BaP was a minor constituent even at the earliest time points and metabolite profiles markedly distinct. AMS, coupled with UPLC, could be used in humans to enhance the accuracy of pharmacokinetics, toxicokinetics and risk assessment of environmental carcinogens.

Comparative mechanisms of PAH toxicity by benzo[a]pyrene and dibenzo[def,p]chrysene in primary human bronchial epithelial cells cultured at air-liquid interface.

Current assumption for assessing carcinogenic risk of polycyclic aromatic hydrocarbons (PAHs) is that they function through a common mechanism of action; however, recent studies demonstrate that PAHs can act through unique mechanisms potentially contributing to cancer outcomes in a non-additive manner. Using a primary human 3D bronchial epithelial culture (HBEC) model, we assessed potential differences in mechanism of toxicity for two PAHs, benzo[a]pyrene (BAP) and dibenzo[def,p]chrysene (DBC), compared to a complex PAH mixture based on short-term biosignatures identified from transcriptional profiling. Differentiated bronchial epithelial cells were treated with BAP (100-500 µg/ml), DBC (10 µg/ml), and coal tar extract (CTE 500-1500 µg/ml, SRM1597a) for 48 h and gene expression was measured by RNA sequencing or quantitative PCR. Comparison of BAP and DBC gene signatures showed that the majority of genes (~60%) were uniquely regulated by treatment, including signaling pathways for inflammation and DNA damage by DBC and processes for cell cycle, hypoxia and oxidative stress by BAP. Specifically, BAP upregulated targets of AhR, NRF2, and KLF4, while DBC downregulated these same targets, suggesting a chemical-specific pattern in transcriptional regulation involved in antioxidant response, potentially contributing to differences in PAH potency. Other processes were regulated in common by all PAH treatments, BAP, DBC and CTE, including downregulation of genes involved in cell adhesion and reduced functional measurements of barrier integrity. This work supports prior in vivo studies and demonstrates the utility of profiling short-term biosignatures in an organotypic 3D model to identify mechanisms linked to carcinogenic risk of PAHs in humans.

Transcriptional changes in innate immunity genes in head kidneys from Aeromonas salmonicida-challenged rainbow trout fed a mixture of polycyclic aromatic hydrocarbons.

We previously observed that exposure to a complex mixture of high molecular weight polycyclic aromatic hydrocarbons (PAHs) increased sensitivity of rainbow trout (Oncorhynchus mykiss) to subsequent challenge with Aeromonas salmonicida, the causative agent of furunculosis. In this study, we evaluate potential mechanisms associated with disease susceptibility from combined environmental factors of dietary PAH exposure and pathogen challenge. Rainbow trout were fed a mixture of ten high molecular weight PAHs at an environmentally relevant concentration (7.82µg PAH mixture/g fish/day) or control diet for 50 days. After 50 days of PAH exposure, fish were challenged with either Aeromonas salmonicida at a lethal concentration 30 (LC30) or growth media without the pathogen (mock challenge). Head kidneys were collected 2, 4, 10 and 20 days after challenge and gene expression (q<0.05) was evaluated among treatments. In animals fed the PAH contaminated diet, we observed down-regulation of expression for innate immune system genes in pathways (p<0.05) for the terminal steps of the complement cascade (complement component C6) and other bacteriolytic processes (lysozyme type II) potentially underlying increased disease susceptibility after pathogen challenge. Increased expression of genes associated with hemorrhage/tissue remodeling/inflammation pathways (p<0.05) was likely related to more severe head kidney damage due to infection in PAH-fed compared to control-fed fish. This study is the first to evaluate transcriptional signatures associated with the impact of chronic exposure to an environmentally relevant mixture of PAHs in disease susceptibility and immunity.

Muscle Segment Homeobox Genes Direct Embryonic Diapause by Limiting Inflammation in the Uterus.

Embryonic diapause is a reproductive strategy widespread in the animal kingdom. This phenomenon is defined by a temporary arrest in blastocyst growth and metabolic activity within a quiescent uterus without implantation until the environmental and maternal milieu become favorable for pregnancy to progress. We found that uterine Msx expression persists during diapause across species; their inactivation in the mouse uterus results in termination of diapause with the development of implantation-like responses ("pseudoimplantation") that ultimately succumbed to resorption. To understand the cause of this failure, we compared proteome profiles between floxed and Msx-deleted uteri. In deleted uteri, several functional networks, including transcription/translation, ubiquitin-proteasome, inflammation, and endoplasmic reticulum stress, were dysregulated. Computational modeling predicted intersection of these pathways on an enhanced inflammatory signature. Further studies showed that this signature was reflected in increased phosphorylated IκB levels and nuclear NFκB in deleted uteri. This was associated with enhanced proteasome activity and endoplasmic reticulum stress. Interestingly, treatment with anti-inflammatory glucocorticoid (dexamethasone) reduced the inflammatory signature with improvement of the diapause phenotype. These findings highlight an unexpected role of uterine Msx in limiting aberrant inflammatory responses to maintain embryonic diapause.

Data integration reveals key homeostatic mechanisms following low dose radiation exposure.

The goal of this study was to define pathways regulated by low dose radiation to understand how biological systems respond to subtle perturbations in their environment and prioritize pathways for human health assessment. Using an in vitro 3-D human full thickness skin model, we have examined the temporal response of dermal and epidermal layers to 10 cGy X-ray using transcriptomic, proteomic, phosphoproteomic and metabolomic platforms. Bioinformatics analysis of each dataset independently revealed potential signaling mechanisms affected by low dose radiation, and integrating data shed additional insight into the mechanisms regulating low dose responses in human tissue. We examined direct interactions among datasets (top down approach) and defined several hubs as significant regulators, including transcription factors (YY1, MYC and CREB1), kinases (CDK2, PLK1) and a protease (MMP2). These data indicate a shift in response across time - with an increase in DNA repair, tissue remodeling and repression of cell proliferation acutely (24-72h). Pathway-based integration (bottom up approach) identified common molecular and pathway responses to low dose radiation, including oxidative stress, nitric oxide signaling and transcriptional regulation through the SP1 factor that would not have been identified by the individual data sets. Significant regulation of key downstream metabolites of nitrative stress was measured within these pathways. Among the features identified in our study, the regulation of MMP2 and SP1 was experimentally validated. Our results demonstrate the advantage of data integration to broadly define the pathways and networks that represent the mechanisms by which complex biological systems respond to perturbation.

Cytochrome P450 1b1 in polycyclic aromatic hydrocarbon (PAH)-induced skin carcinogenesis: Tumorigenicity of individual PAHs and coal-tar extract, DNA adduction and expression of select genes in the Cyp1b1 knockout mouse.

FVB/N mice wild-type, heterozygous or null for Cyp 1b1 were used in a two-stage skin tumor study comparing PAH, benzo[a]pyrene (BaP), dibenzo[def,p]chrysene (DBC), and coal tar extract (CTE, SRM 1597a). Following 20 weeks of promotion with TPA the Cyp 1b1 null mice, initiated with DBC, exhibited reductions in incidence, multiplicity, and progression. None of these effects were observed with BaP or CTE. The mechanism of Cyp 1b1-dependent alteration of DBC skin carcinogenesis was further investigated by determining expression of select genes in skin from DBC-treated mice 2, 4 and 8h post-initiation. A significant reduction in levels of Cyp 1a1, Nqo1 at 8h and Akr 1c14 mRNA was observed in Cyp 1b1 null (but not wt or het) mice, whereas no impact was observed in Gst a1, Nqo 1 at 2 and 4h or Akr 1c19 at any time point. Cyp 1b1 mRNA was not elevated by DBC. The major covalent DNA adducts, dibenzo[def,p]chrysene-(±)-11,12-dihydrodiol-cis and trans-13,14-epoxide-deoxyadenosine (DBCDE-dA) were quantified by UHPLC-MS/MS 8h post-initiation. Loss of Cyp1 b1 expression reduced DBCDE-dA adducts in the skin but not to a statistically significant degree. The ratio of cis- to trans-DBCDE-dA adducts was higher in the skin than other target tissues such as the spleen, lung and liver (oral dosing). These results document that Cyp 1b1 plays a significant role in bioactivation and carcinogenesis of DBC in a two-stage mouse skin tumor model and that loss of Cyp 1b1 has little impact on tumor response with BaP or CTE as initiators.

Release of severe acute respiratory syndrome coronavirus nuclear import block enhances host transcription in human lung cells.

The severe acute respiratory syndrome coronavirus accessory protein ORF6 antagonizes interferon signaling by blocking karyopherin-mediated nuclear import processes. Viral nuclear import antagonists, expressed by several highly pathogenic RNA viruses, likely mediate pleiotropic effects on host gene expression, presumably interfering with transcription factors, cytokines, hormones, and/or signaling cascades that occur in response to infection. By bioinformatic and systems biology approaches, we evaluated the impact of nuclear import antagonism on host expression networks by using human lung epithelial cells infected with either wild-type virus or a mutant that does not express ORF6 protein. Microarray analysis revealed significant changes in differential gene expression, with approximately twice as many upregulated genes in the mutant virus samples by 48 h postinfection, despite identical viral titers. Our data demonstrated that ORF6 protein expression attenuates the activity of numerous karyopherin-dependent host transcription factors (VDR, CREB1, SMAD4, p53, EpasI, and Oct3/4) that are critical for establishing antiviral responses and regulating key host responses during virus infection. Results were confirmed by proteomic and chromatin immunoprecipitation assay analyses and in parallel microarray studies using infected primary human airway epithelial cell cultures. The data strongly support the hypothesis that viral antagonists of nuclear import actively manipulate host responses in specific hierarchical patterns, contributing to the viral pathogenic potential in vivo. Importantly, these studies and modeling approaches not only provide templates for evaluating virus antagonism of nuclear import processes but also can reveal candidate cellular genes and pathways that may significantly influence disease outcomes following severe acute respiratory syndrome coronavirus infection in vivo.

Retinoic acid-dependent regulation of miR-19 expression elicits vertebrate axis defects.

Retinoic acid (RA) is involved in multifarious and complex functions necessary for vertebrate development. RA signaling is reliant on strict enzymatic regulation of RA synthesis and metabolism. Improper spatiotemporal expression of RA during development can result in vertebrate axis defects. microRNAs (miRNAs) are also pivotal in orchestrating developmental processes. While mechanistic links between miRNAs and axial development are established, the role of miRNAs in regulating metabolic enzymes responsible for RA abundance during axis formation has yet to be elucidated. Our results uncovered a role of miR-19 family members in controlling RA metabolism through the regulation of CYP26A1 during vertebrate axis formation. Global miRNA expression profiling showed that developmental RA exposure suppressed the expression of miR-19 family members during zebrafish somitogenesis. A reporter assay confirmed that cyp26a1 is a bona fide target of miR-19 in vivo. Transient knockdown of miR-19 phenocopied axis defects caused by RA exposure. Exogenous miR-19 rescued the axis defects induced by RA exposure. Taken together, these results indicate that the teratogenic effects of RA exposure result, in part, from repression of miR-19 expression and subsequent misregulation of cyp26a1. This highlights a previously unidentified role of miR-19 in facilitating vertebrate axis development via regulation of RA signaling.

Comparative iron oxide nanoparticle cellular dosimetry and response in mice by the inhalation and liquid cell culture exposure routes.

BACKGROUND: Toxicity testing the rapidly growing number of nanomaterials requires large scale use of in vitro systems under the presumption that these systems are sufficiently predictive or descriptive of responses in in vivo systems for effective use in hazard ranking. We hypothesized that improved relationships between in vitro and in vivo models of experimental toxicology for nanomaterials would result from placing response data in vitro and in vivo on the same dose scale, the amount of material associated with cells. METHODS: Balb/c mice were exposed nose-only to an aerosol (68.6 nm CMD, 19.9 mg/m(3), 4 hours) generated from of 12.8 nm superparamagnetic iron oxide particles (SPIO). Target cell doses were calculated, histological evaluations conducted, and biomarkers of response were identified by global transcriptomics. Representative murine epithelial and macrophage cell types were exposed in vitro to the same material in liquid suspension for four hours and levels of nanoparticle regulated cytokine transcripts identified in vivo were quantified as a function of measured nanoparticle cellular dose. RESULTS: Target tissue doses of 0.009-0.4 µg SPIO/cm(2) in lung led to an inflammatory response in the alveolar region characterized by interstitial inflammation and macrophage infiltration. In vitro, higher target tissue doses of ~1.2-4 µg SPIO/ cm(2) of cells were required to induce transcriptional regulation of markers of inflammation, CXCL2 & CCL3, in C10 lung epithelial cells. Estimated in vivo macrophage SPIO nanoparticle doses ranged from 1-100 pg/cell, and induction of inflammatory markers was observed in vitro in macrophages at doses of 8-35 pg/cell. CONCLUSIONS: Application of target tissue dosimetry revealed good correspondence between target cell doses triggering inflammatory processes in vitro and in vivo in the alveolar macrophage population, but not in the epithelial cells of the alveolar region. These findings demonstrate the potential for target tissue dosimetry to enable the more quantitative comparison of in vitro and in vivo systems and advance their use for hazard assessment and extrapolation to humans. The mildly inflammogentic cellular doses experienced by mice were similar to those calculated for humans exposed to the same material at the existing permissible exposure limit of 10 mg/m(3) iron oxide (as Fe).

Assessing susceptibility for polycyclic aromatic hydrocarbon toxicity in an in vitro 3D respiratory model for asthma.

There is increased emphasis on understanding cumulative risk from the combined effects of chemical and non-chemical stressors as it relates to public health. Recent animal studies have identified pulmonary inflammation as a possible modifier and risk factor for chemical toxicity in the lung after exposure to inhaled pollutants; however, little is known about specific interactions and potential mechanisms of action. In this study, primary human bronchial epithelial cells (HBEC) cultured in 3D at the air-liquid interface (ALI) are utilized as a physiologically relevant model to evaluate the effects of inflammation on toxicity of polycyclic aromatic hydrocarbons (PAHs), a class of contaminants generated from incomplete combustion of fossil fuels. Normal HBEC were differentiated in the presence of IL-13 for 14 days to induce a profibrotic phenotype similar to asthma. Fully differentiated normal and IL-13 phenotype HBEC were treated with benzo[a]pyrene (BAP; 1-40 µg/mL) or 1% DMSO/PBS vehicle at the ALI for 48 h. Cells were evaluated for cytotoxicity, barrier integrity, and transcriptional biomarkers of chemical metabolism and inflammation by quantitative PCR. Cells with the IL-13 phenotype treated with BAP result in significantly (p < 0.05) decreased barrier integrity, less than 50% compared to normal cells. The effect of BAP in the IL-13 phenotype was more apparent when evaluating transcriptional biomarkers of barrier integrity in addition to markers of mucus production, goblet cell hyperplasia, type 2 asthmatic inflammation and chemical metabolism, which all resulted in dose-dependent changes (p < 0.05) in the presence of BAP. Additionally, RNA sequencing data showed that the HBEC with the IL-13 phenotype may have increased potential for uncontrolled proliferation and decreased capacity for immune response after BAP exposure compared to normal phenotype HBEC. These data are the first to evaluate the role of combined environmental factors associated with inflammation from pre-existing disease and PAH exposure on pulmonary toxicity in a physiologically relevant human in vitro model.

Application of a fuzzy neural network model in predicting polycyclic aromatic hydrocarbon-mediated perturbations of the Cyp1b1 transcriptional regulatory network in mouse skin.

Polycyclic aromatic hydrocarbons (PAHs) are present in the environment as complex mixtures with components that have diverse carcinogenic potencies and mostly unknown interactive effects. Non-additive PAH interactions have been observed in regulation of cytochrome P450 (CYP) gene expression in the CYP1 family. To better understand and predict biological effects of complex mixtures, such as environmental PAHs, an 11 gene input-1 gene output fuzzy neural network (FNN) was developed for predicting PAH-mediated perturbations of dermal Cyp1b1 transcription in mice. Input values were generalized using fuzzy logic into low, medium, and high fuzzy subsets, and sorted using k-means clustering to create Mamdani logic functions for predicting Cyp1b1 mRNA expression. Model testing was performed with data from microarray analysis of skin samples from FVB/N mice treated with toluene (vehicle control), dibenzo[def,p]chrysene (DBC), benzo[a]pyrene (BaP), or 1 of 3 combinations of diesel particulate extract (DPE), coal tar extract (CTE) and cigarette smoke condensate (CSC) using leave-one-out cross-validation. Predictions were within 1 log(2) fold change unit of microarray data, with the exception of the DBC treatment group, where the unexpected down-regulation of Cyp1b1 expression was predicted but did not reach statistical significance on the microarrays. Adding CTE to DPE was predicted to increase Cyp1b1 expression, whereas adding CSC to CTE and DPE was predicted to have no effect, in agreement with microarray results. The aryl hydrocarbon receptor repressor (Ahrr) was determined to be the most significant input variable for model predictions using back-propagation and normalization of FNN weights.

Structurally distinct polycyclic aromatic hydrocarbons induce differential transcriptional responses in developing zebrafish.

Polycyclic aromatic hydrocarbons (PAHs) are ubiquitous in the environment as components of fossil fuels and by-products of combustion. These multi-ring chemicals differentially activate the aryl hydrocarbon receptor (AHR) in a structurally dependent manner, and induce toxicity via both AHR-dependent and -independent mechanisms. PAH exposure is known to induce developmental malformations in zebrafish embryos, and recent studies have shown cardiac toxicity induced by compounds with low AHR affinity. Unraveling the potentially diverse molecular mechanisms of PAH toxicity is essential for understanding the hazard posed by complex PAH mixtures present in the environment. We analyzed transcriptional responses to PAH exposure in zebrafish embryos exposed to benz(a)anthracene (BAA), dibenzothiophene (DBT) and pyrene (PYR) at concentrations that induced developmental malformations by 120 h post-fertilization (hpf). Whole genome microarray analysis of mRNA expression at 24 and 48 hpf identified genes that were differentially regulated over time and in response to the three PAH structures. PAH body burdens were analyzed at both time points using GC-MS, and demonstrated differences in PAH uptake into the embryos. This was important for discerning dose-related differences from those that represented unique molecular mechanisms. While BAA misregulated the least number of transcripts, it caused strong induction of cyp1a and other genes known to be downstream of the AHR, which were not induced by the other two PAHs. Analysis of functional roles of misregulated genes and their predicted regulatory transcription factors also distinguished the BAA response from regulatory networks disrupted by DBT and PYR exposure. These results indicate that systems approaches can be used to classify the toxicity of PAHs based on the networks perturbed following exposure, and may provide a path for unraveling the toxicity of complex PAH mixtures.

Diet-induced obesity reprograms the inflammatory response of the murine lung to inhaled endotoxin.

The co-occurrence of environmental factors is common in complex human diseases and, as such, understanding the molecular responses involved is essential to determine risk and susceptibility to disease. We have investigated the key biological pathways that define susceptibility for pulmonary infection during obesity in diet-induced obese (DIO) and regular weight (RW) C57BL/6 mice exposed to inhaled lipopolysaccharide (LPS). LPS induced a strong inflammatory response in all mice as indicated by elevated cell counts of macrophages and neutrophils and levels of proinflammatory cytokines (MDC, MIP-1γ, IL-12, RANTES) in the bronchoalveolar lavage fluid. Additionally, DIO mice exhibited 50% greater macrophage cell counts, but decreased levels of the cytokines, IL-6, TARC, TNF-α, and VEGF relative to RW mice. Microarray analysis of lung tissue showed over half of the LPS-induced expression in DIO mice consisted of genes unique for obese mice, suggesting that obesity reprograms how the lung responds to subsequent insult. In particular, we found that obese animals exposed to LPS have gene signatures showing increased inflammatory and oxidative stress response and decreased antioxidant capacity compared with RW. Because signaling pathways for these responses can be common to various sources of environmentally induced lung damage, we further identified biomarkers that are indicative of specific toxicant exposure by comparing gene signatures after LPS exposure to those from a parallel study with cigarette smoke. These data show obesity may increase sensitivity to further insult and that co-occurrence of environmental stressors result in complex biosignatures that are not predicted from analysis of individual exposures.

Impaired transcriptional response of the murine heart to cigarette smoke in the setting of high fat diet and obesity.

Smoking and obesity are each well-established risk factors for cardiovascular heart disease, which together impose earlier onset and greater severity of disease. To identify early signaling events in the response of the heart to cigarette smoke exposure within the setting of obesity, we exposed normal weight and high fat diet-induced obese (DIO) C57BL/6 mice to repeated inhaled doses of mainstream (MS) or sidestream (SS) cigarette smoke administered over a two week period, monitoring effects on both cardiac and pulmonary transcriptomes. MS smoke (250 µg wet total particulate matter (WTPM)/L, 5 h/day) exposures elicited robust cellular and molecular inflammatory responses in the lung with 1466 differentially expressed pulmonary genes (p < 0.01) in normal weight animals and a much-attenuated response (463 genes) in the hearts of the same animals. In contrast, exposures to SS smoke (85 µg WTPM/L) with a CO concentration equivalent to that of MS smoke (~250 CO ppm) induced a weak pulmonary response (328 genes) but an extensive cardiac response (1590 genes). SS smoke and to a lesser extent MS smoke preferentially elicited hypoxia- and stress-responsive genes as well as genes predicting early changes of vascular smooth muscle and endothelium, precursors of cardiovascular disease. The most sensitive smoke-induced cardiac transcriptional changes of normal weight mice were largely absent in DIO mice after smoke exposure, while genes involved in fatty acid utilization were unaffected. At the same time, smoke exposure suppressed multiple proteome maintenance genes induced in the hearts of DIO mice. Together, these results underscore the sensitivity of the heart to SS smoke and reveal adaptive responses in healthy individuals that are absent in the setting of high fat diet and obesity.

MicroRNAs control neurobehavioral development and function in zebrafish.

microRNAs (miRNAs) have emerged as regulators of a broad spectrum of neurodevelopmental processes, including brain morphogenesis, neuronal differentiation, and survival. While the role of miRNAs in establishing and maintaining the developing nervous system is widely appreciated, the developmental neurobehavioral role of miRNAs has yet to be defined. Here we show that transient disruption of brain morphogenesis by ethanol exposure results in behavioral hyperactivity in larval zebrafish challenged with changes in lighting conditions. Aberrations in swimming activity persist in juveniles that were developmentally exposed to ethanol. During early neurogenesis, multiple gene expression profiling studies revealed widespread changes in mRNA and miRNA abundance in ethanol-exposed embryos. Consistent with a role for miRNAs in neurobehavioral development, target prediction analyses identified multiple miRNAs misexpressed in the ethanol-exposed cohorts that were also predicted to target inversely expressed transcripts known to influence brain morphogenesis. In vivo knockdown of miR-9/9* or miR-153c persistently phenocopied the effect of ethanol on larval and juvenile swimming behavior. Structural analyses performed on adults showed that repression of miR-153c during development impacts craniofacial skeletal development. Together, these data support an integral role for miRNAs in the establishment of vertebrate neurobehavioral and skeletal systems.