Passivating the Background of Living Microbes with a Zwitterionic Peptide for Therapies.

Living microbial therapies have been proposed as a course of action for a variety of diseases. However, problematic interactions between the host immune system and the microbial organism present significant clinical concerns. Previously, we developed a genetically encoded superhydrophilic zwitterionic peptide, termed EKP, to mimic low-immunogenic zwitterionic materials, which have been used for the chemical modification of biologics such as protein and nucleic acid drugs to increase their in vivo circulation time and reduce their immunogenicity. Herein, we demonstrate the protective effects of the EKP polypeptide genetically cloaking the surface of Saccharomyces cerevisiae as a model microbe in both in vitro and in vivo systems. First, we show that EKP peptide cloaking suppresses the interactions between yeast cells and their specific antibodies, thereby illustrating its cloaking behavior. Then, we examine the in vitro interactions between EKP peptide surface cloaked yeast cells and murine macrophage cells, which exhibit phagocytotic behavior in the presence of foreign microbes. Our results indicate that EKP cloaking suppresses macrophage interactions and thus reduces phagocytosis. Furthermore, EKP cloaked yeast cells demonstrate a prolonged circulation time in mice in vivo.

Ocular manifestations of CHARGE syndrome in a pediatric cohort with genotype/phenotype analysis.

CHARGE syndrome is a rare multi-system condition associated with CHD7 variants. However, ocular manifestations and particularly ophthalmic genotype-phenotype associations, are not well-studied. This study evaluated ocular manifestations and genotype-phenotype associations in pediatric patients with CHARGE syndrome. A retrospective chart review included pediatric patients under 20 years-old with clinical diagnosis of CHARGE syndrome and documented ophthalmic examination. Demographics, genetic testing, and ocular findings were collected. Comprehensive literature review enhanced the genotype-phenotype analysis. Forty-two patients (20 male) underwent eye examination at an average age of 9.45 ± 6.52 years-old. Thirty-nine (93%) had ophthalmic manifestations in at least one eye. Optic nerve/chorioretinal colobomas were most common (38 patients), followed by microphthalmia (13), cataract (6), and iris colobomas (4). Extraocular findings included strabismus (32 patients), nasolacrimal duct obstructions (11, 5 with punctal agenesis), and cranial nerve VII palsy (10). Genotype-phenotype analyses (27 patients) showed variability in ocular phenotypes without association to location or variant types. Splicing (10 patients) and frameshift (10) variants were most prevalent. Patients with CHARGE syndrome may present with a myriad of ophthalmic manifestations. There is limited data regarding genotype-phenotype correlations and additional studies are needed.

Comparison of Rat Hepatocyte 2D-Monocultures and Hepatocytes Non-Parenchymal Cell Co-Cultures for Assessing Chemical Toxicity.

Liver responses are the most common endpoints used as the basis for setting exposure standards. Liver hepatocytes play a vital role in biotransformation of xenobiotics, but non-parenchymal cells (NPCs) in the liver are also involved in certain liver responses. Development of in vitro systems that more faithfully capture liver responses to reduce reliance on animals is a major focus of New Approach Methodology (NAMs). Since rodent regulatory studies are frequently the sole source safety assessment data, mode-of-action data, and used for risk assessments, in vitro rodent models that reflect in vivo responses need to be developed to reduce reliance on animal models. In the work presented in this paper, we developed a 2-D hepatocyte monoculture and 2-D liver cell co-culture system using rat liver cells. These models were assessed for conditions for short-term stability of the cultures and phenotypic and transcriptomic responses of 2 prototypic hepatotoxicants compounds - acetaminophen and phenobarbital. The optimized multi-cellular 2-D culture required use of freshly prepared hepatocytes and NPCs from a single rat, a 3:1 ratio of hepatocytes to NPCs and growth medium using 50% Complete Williams E medium (WEM) and 50% Endothelial Cell Medium (ECM). The transcriptomic responses of the 2 model systems to PB were compared to previous studies from TG-Gates on the gene expression changes in intact rats and the co-culture model responses were more representative of the in vivo responses. Transcriptomic read-outs promise to move beyond conventional phenotypic evaluations with these in vitro NAMs and provide insights about modes of action.

Comparative transcriptomics and eQTL mapping of response to Melampsora americana in selected Salix purpurea F2 progeny.

BACKGROUND: Melampsora spp. rusts are the greatest pathogen threat to shrub willow (Salix spp.) bioenergy crops. Genetic resistance is key to limit the effects of these foliar diseases on host response and biomass yield, however, the genetic basis of host resistance has not been characterized. The addition of new genomic resources for Salix provides greater power to investigate the interaction between S. purpurea and M. americana, species commonly found in the Northeast US. Here, we utilize 3' RNA-seq to investigate host-pathogen interactions following controlled inoculations of M. americana on resistant and susceptible F2 S. purpurea genotypes identified in a recent QTL mapping study. Differential gene expression, network analysis, and eQTL mapping were used to contrast the response to inoculation and to identify associated candidate genes. RESULTS: Controlled inoculation in a replicated greenhouse study identified 19 and 105 differentially expressed genes between resistant and susceptible genotypes at 42 and 66 HPI, respectively. Defense response gene networks were activated in both resistant and susceptible genotypes and enriched for many of the same defense response genes, yet the hub genes of these common response modules showed greater mean expression among the resistant plants. Further, eight and six eQTL hotspots were identified at 42 and 66 HPI, respectively. The combined results of three analyses highlight 124 candidate genes in the host for further analysis while analysis of pathogen RNA showed differential expression of 22 genes, two of which are candidate pathogen effectors. CONCLUSIONS: We identified two differentially expressed M. americana transcripts and 124 S. purpurea genes that are good candidates for future studies to confirm their role in conferring resistance.

Impacts of a Zwitterionic Peptide on its Fusion Protein.

Therapeutic proteins frequently need to be modified with high-molecular-weight molecules to improve their pharmacokinetic properties. The genetic linkage of therapeutic proteins to a high-molecular-weight zwitterionic peptide, termed EKP, offers a promising approach. As with any protein modification, EKP could impact the structural behavior and receptor binding properties of the linked therapeutic protein, thereby altering its bioactivity. To evaluate the effects of EKP on therapeutic proteins, we study the receptor binding properties of high-molecular-weight EKP linked to the growth colony-stimulating factor (GCSF) using the genetically based yeast display platform. We find that yeast-displayed EKP-GCSF and GCSF exhibits similar binding to its receptor GCSF-R, suggesting that EKP does not hinder receptor binding. Furthermore, yeast-displayed EKP-GCSF demonstrates protection against thermal denaturation compared to GCSF. Similarly, to study the structural effects of EKP on GCSF, we employ in silico modeling using alphaFold2 in conjunction with molecular dynamics (MD) simulations. Likewise, in silico modeling reveals that EKP does not alter the structural behavior of GCSF. Finally, we demonstrate the functional benefits of EKP, by which the EKP-GCSF fusion protein produced in Escherichia coli exhibits improved pharmacokinetics and prolonged bioactivity in vivo.

The Melampsora americana Population on Salix purpurea in the Great Lakes Region Is Highly Diverse with a Contributory Influence of Clonality.

Shrub willows (Salix spp.) are emerging as a viable lignocellulosic, second-generation bioenergy crop with many growth characteristics favorable for marginal lands in New York State and surrounding areas. Willow rust, caused by members of the genus Melampsora, is the most limiting disease of shrub willow in this region and remains extremely understudied. In this study, genetic diversity, genetic structure, and pathogen clonality were examined in Melampsora americana over two growing seasons via genotyping-by-sequencing to identify single-nucleotide polymorphism markers. In conjunction with this project, a reference genome of rust isolate R15-033-03 was generated to aid in variant discovery. Sampling between years allowed regional and site-specific investigation into population dynamics, in the context of both wild and cultivated hosts within high-density plantings. This work revealed that this pathogen is largely panmictic over the sampled areas, with few sites showing moderate genetic differentiation. These data support the hypothesis of sexual recombination between growing seasons because no genotype persisted across the two years of sampling. Additionally, clonality was determined as a driver of pathogen populations within cultivated fields and single shrubs; however, there is also evidence of high genetic diversity of rust isolates in all settings. This work provides a framework for M. americana population structure in the Great Lakes region, providing crucial information that can aid in future resistance breeding efforts.

A systematic approach to evaluate plausible modes of actions for mouse lung tumors in mice exposed to 4-methylimidozole.

The National Toxicology Program (NTP) reported that chronic dietary exposure to 4-methylimidazole (4-MeI) increased the incidence of lung adenomas/carcinomas beyond the normally high spontaneous rate in B6C3F1 mice. To examine plausible modes of action (MoAs) for mouse lung tumors (MLTs) upon exposure to high levels of 4-MeI, and their relevance in assessing human risk, a systematic approach was used to identify and evaluate mechanistic data (in vitro and in vivo) in the primary and secondary literature, along with high-throughput screening assay data. Study quality, relevance, and activity of mechanistic data identified across the evidence-base were organized according to key characteristics of carcinogens (KCCs) to identify potential key events in known or novel MLT MoAs. Integration of these evidence streams provided confirmation that 4-MeI lacks genotoxic and cytotoxic activity with some evidence to support a lack of mitogenic activity. Further evaluation of contextual and chemical-specific characteristics of 4-MeI was consequently undertaken. Due to lack of genotoxicity, along with transcriptomic and histopathological lung changes up to 28 and 90 days of exposure, the collective evidence suggests MLTs observed following exposure to high levels of 4-MeI develop at a late stage in the mouse chronic bioassay, albeit the exact MoA remains unclear.

Zwitterionic Polymer Conjugated Glucagon-like Peptide-1 for Prolonged Glycemic Control.

Glucagon-like peptide-1 (GLP-1) is of particular interest for treating type 2 diabetes mellitus (T2DM), as it induces insulin secretion in a glucose-dependent fashion and has the potential to facilitate weight control. However, native GLP-1 is a short incretin peptide that is susceptible to fast proteolytic inactivation and rapid clearance from the circulation. Various GLP-1 analogs and bioconjugation of GLP-1 analogs have been developed to counter these issues, but these modifications are frequently accompanied by the sacrifice of potency and the induction of immunogenicity. Here, we demonstrated that with the conjugation of a zwitterionic polymer, poly(carboxybetaine) (pCB), the pharmacokinetic properties of native GLP-1 were greatly enhanced without serious negative effects on its potency and secondary structure. The pCB conjugated GLP-1 further provided glycemic control for up to 6 days in a mouse study. These results illustrate that the conjugation of pCB could realize the potential of using native GLP-1 for prolonged glycemic control in treating T2DM.

A toxicogenomic approach for the risk assessment of the food contaminant acetamide.

Acetamide (CAS 60-35-5) is detected in common foods. Chronic rodent bioassays led to its classification as a group 2B possible human carcinogen due to the induction of liver tumors in rats. We used a toxicogenomics approach in Wistar rats gavaged daily for 7 or 28 days at doses of 300 to 1500 mg/kg/day (mkd) to determine a point of departure (POD) and investigate its mode of action (MoA). Ki67 labeling was increased at doses ≥750 mkd up to 3.3-fold representing the most sensitive apical endpoint. Differential gene expression analysis by RNA-Seq identified 1110 and 1814 differentially expressed genes in male and female rats, respectively, following 28 days of treatment. Down-regulated genes were associated with lipid metabolism while up-regulated genes included cell signaling, immune response, and cell cycle functions. Benchmark dose (BMD) modeling of the Ki67 labeling index determined the BMD10 lower confidence limit (BMDL10) as 190 mkd. Transcriptional BMD modeling revealed excellent concordance between transcriptional POD and apical endpoints. Collectively, these results indicate that acetamide is most likely acting through a mitogenic MoA, though specific key initiating molecular events could not be elucidated. A POD value of 190 mkd determined for cell proliferation is suggested for risk assessment purposes.

The role of fit-for-purpose assays within tiered testing approaches: A case study evaluating prioritized estrogen-active compounds in an in vitro human uterotrophic assay.

Chemical risk assessment relies on toxicity tests that require significant numbers of animals, time and costs. For the >30,000 chemicals in commerce, the current scale of animal testing is insufficient to address chemical safety concerns as regulatory and product stewardship considerations evolve to require more comprehensive understanding of potential biological effects, conditions of use, and associated exposures. We demonstrate the use of a multi-level new approach methodology (NAMs) strategy for hazard- and risk-based prioritization to reduce animal testing. A Level 1/2 chemical prioritization based on estrogen receptor (ER) activity and metabolic activation using ToxCast data was used to select 112 chemicals for testing in a Level 3 human uterine cell estrogen response assay (IKA assay). The Level 3 data were coupled with quantitative in vitro to in vivo extrapolation (Q-IVIVE) to support bioactivity determination (as a surrogate for hazard) in a tissue-specific context. Assay AC50s and Q-IVIVE were used to estimate human equivalent doses (HEDs), and HEDs were compared to rodent uterotrophic assay in vivo-derived points of departure (PODs). For substances active both in vitro and in vivo, IKA assay-derived HEDs were lower or equivalent to in vivo PODs for 19/23 compounds (83%). Activity exposure relationships were calculated, and the IKA assay was as or more protective of human health than the rodent uterotrophic assay for all IKA-positive compounds. This study demonstrates the utility of biologically relevant fit-for-purpose assays and supports the use of a multi-level strategy for chemical risk assessment.

Elucidating Molecular Design Principles for Charge-Alternating Peptides.

The therapeutic potential of protein drugs has been hindered by difficulties with long-term stability and rapid clearance from the body. Recombinant fusion proteins provide a scalable platform for engineered biologics, whereby a polypeptide domain is appended to alter the physical characteristics of a therapeutic protein and enhance its pharmaceutical viability. Two simple design principles for recombinant fusion proteins, based on the physical properties of the polypeptide domain, have been separately applied to address issues with the stability and delivery of biologics. "Conformationally disordered" peptides, exemplified by the homo amino acid peptide polyG, have been shown to increase the circulation half-life and bioactivity of protein therapeutics in vivo. Superhydrophilic peptides, exemplified by the alternating-charge peptide poly(EK), have been shown to increase the thermostability of proteins in vitro. The combination of superhydrophilicity and conformational disorder in a single fusion peptide could simultaneously address concerns regarding the stability and therapeutic lifetime of biologics. In the current work, we use enhanced sampling molecular dynamics (MD) simulations to investigate the conformational ensemble of poly(EK) and glycine-substituted poly(EK) variants and validate our structural predictions with circular dichroism (CD). We find the (EK)15 peptide exhibits a high propensity for forming antiparallel ß-strand secondary structures, which are stabilized by extensive salt bridging of the positive and negative side chains. MD simulations predict that limited glycine substitutions effectively disrupt the secondary structure and promote disordered conformations at physiologically relevant temperatures. We conclude that the conformational disorder of alternating-charge peptides should be taken into account to improve their suitability for drug delivery applications. We also contribute a computational approach to quantify conformational disorder in polypeptides, which should facilitate the de novo design of effective fusion proteins.

Differential Susceptibility of Diverse Salix spp. to Melampsora americana and Melampsora paradoxa.

Melampsora spp. willow rust is the most serious disease of shrub willow bioenergy production in the northeastern United States. Recent phylogenetic studies have identified several Melampsora spp. present on willow in the Northeast; however, in-depth understanding of Melampsora spp. host susceptibility remain unresolved. In this study, a panel of 82 rust isolates collected from the northeastern United States were genotyped via ribosomal DNA sequencing and a subset of these isolates were assayed for host susceptibility. This work revealed that Melampsora americana is the most prevalent species in the sampled geographic region and that there is potential for rust resistance breeding using the Salix spp. taxa assayed. Additionally, leaf morphology traits of these Salix spp. hosts were quantified for correlation analysis, revealing that trichome density and stomata density are possible contributors to resistance. This work provides foundational rust pathology information, which is crucial for M. americana resistance breeding.

Zwitterionic Peptide Cloak Mimics Protein Surfaces for Protein Protection.

Inspired by the amino acid composition of natural protein surfaces, we developed a zwitterionic cloak containing multi-layers of short alternating glutamic acid and lysine (EK) peptides as a facile, highly effective and low-immunogenicity approach for the protection and delivery of biotherapeutics. Each EK layer grafted to proteins provides multiple times of new lysine reaction sites for the growth of subsequent EK layers. This unique design allows EK peptides to achieve high coating density on proteins, overcoming the limitation of traditional conjugation strategies that rely on the number of innate lysine groups. A triple-layer EK cloak manifests to successfully eliminate the specific and non-specific interactions of protected asparaginase with biological media while prolong the drug circulation time and significantly mitigate its immunogenicity in vivo, suggesting an EK peptide cloak as a promising approach to improve the safety and efficacy of biotherapeutics.

An in vitro approach for prioritization and evaluation of chemical effects on glucocorticoid receptor mediated adipogenesis.

Rising obesity rates worldwide have socio-economic ramifications. While genetics, diet, and lack of exercise are major contributors to obesity, environmental factors may enhance susceptibility through disruption of hormone homeostasis and metabolic processes. The obesogen hypothesis contends that chemical exposure early in development may enhance adipocyte differentiation, thereby increasing the number of adipocytes and predisposing for obesity and metabolic disease. We previously developed a primary human adipose stem cell (hASC) assay to evaluate the effect of environmental chemicals on PPARG-dependent adipogenesis. Here, the assay was modified to determine the effects of chemicals on the glucocorticoid receptor (GR) pathway. In differentiation cocktail lacking the glucocorticoid agonist dexamethasone (DEX), hASCs do not differentiate into adipocytes. In the presence of GR agonists, adipocyte maturation was observed using phenotypic makers for lipid accumulation, adipokine secretion, and expression of key genes. To evaluate the role of environmental compounds on adipocyte differentiation, progenitor cells were treated with 19 prioritized compounds previously identified by ToxPi as having GR-dependent bioactivity, and multiplexed assays were used to confirm a GR-dependent mode of action. Five chemicals were found to be strong agonists. The assay was also modified to evaluate GR-antagonists, and 8/10 of the hypothesized antagonists inhibited adipogenesis. The in vitro bioactivity data was put into context with extrapolated human steady state concentrations (Css) and clinical exposure data (Cmax). These data support using a human adipose-derived stem cell differentiation assay to test the potential of chemicals to alter human GR-dependent adipogenesis.

Strain-related differences in mouse lung gene expression over a two-year period of inhalation exposure to styrene: Relevance to human risk assessment.

Both CD-1 and C57BL/6 wildtype (C57BL/6-WT) mice show equivalent short-term lung toxicity from exposures to styrene, while long-term tumor responses are greater in CD-1 mice. We analyzed lung gene expression from styrene exposures lasting from 1-day to 2-years in male mice from these two strains, including a Cyp2f2(-/-) knockout (C57BL/6-KO) and a Cyp2F1/2A13/2B6 transgenic mouse (C57BL/6-TG). With short term exposures (1-day to 1-week), CD-1 and C57BL/6-WT mice had thousands of differentially expressed genes (DEGs), consistent with changes in pathways for cell proliferation, cellular lipid metabolism, DNA-replication and inflammation. C57BL/6-WT mice responded within a single day; CD-1 mice required several days of exposure. The numbers of exposure related DEGs were greatly reduced at longer times (4-weeks to 2-years) with enrichment only for biological oxidations in C57BL/6-WT and metabolism of lipids and lipoproteins in CD-1. Gene expression results indicate a non-genotoxic, mouse specific mode of action for short-term styrene responses related to activation of nuclear receptor signaling and cell proliferation. Greater tumor susceptibility in CD-1 mice correlated with the presence of the Pas1 loci, differential Cytochrome P450 gene expression, down-regulation of Nr4a, and greater inflammatory pathway activation. Very few exposure-related responses occurred at any time in C57BL/6-KO or -TG mice indicating that neither the short term nor long term responses of styrene in mice are relevant endpoints for assessing human risks.

Assessing molecular initiating events (MIEs), key events (KEs) and modulating factors (MFs) for styrene responses in mouse lungs using whole genome gene expression profiling following 1-day and multi-week exposures.

Styrene increased lung tumors in mice at chronic inhalation exposures of 20ppm and greater. MIEs, KEs and MFs were examined using gene expression in three strains of male mice (the parental C57BL/6 strain, a CYP2F2(-/-) knock out and a CYP2F2(-/-) transgenic containing human CYP2F1, 2A13 and 2B6). Exposures were for 1-day and 1, 4 and 26weeks. After 1-day exposures at 1, 5, 10, 20, 40 and 120ppm significant increases in differentially expressed genes (DEGs) occurred only in parental strain lungs where there was already an increase in DEGs at 5ppm and then many thousands of DEGs by 120ppm. Enrichment for 1-day and 1-week exposures included cell cycle, mitotic M-M/G1 phases, DNA-synthesis and metabolism of lipids and lipoproteins pathways. The numbers of DEGs decreased steadily over time with no DEGs meeting both statistical significance and fold-change criteria at 26weeks. At 4 and 26weeks, some key transcription factors (TFs) - Nr1d1, Nr1d2, Dbp, Tef, Hlf, Per3, Per2 and Bhlhe40 - were upregulated (|FC|>1.5), while others - Npas, Arntl, Nfil3, Nr4a1, Nr4a2, and Nr4a3 - were down-regulated. At all times, consistent changes in gene expression only occurred in the parental strain. Our results support a MIE for styrene of direct mitogenicity from mouse-specific CYP2F2-mediated metabolites activating Nr4a signaling. Longer-term MFs include down-regulation of Nr4a genes and shifts in both circadian clock TFs and other TFs, linking circadian clock to cellular metabolism. We found no gene expression changes indicative of cytotoxicity or activation of p53-mediated DNA-damage pathways.

Combining transcriptomics and PBPK modeling indicates a primary role of hypoxia and altered circadian signaling in dichloromethane carcinogenicity in mouse lung and liver.

Dichloromethane (DCM) is a lung and liver carcinogen in mice at inhalation exposures≥2000ppm. The modes of action (MOA) of these responses have been attributed to formation of genotoxic, reactive metabolite(s). Here, we examined gene expression in lung and liver from female B6C3F1 mice exposed to 0, 100, 500, 2000, 3000 and 4000ppm DCM for 90days. We also simulated dose measures - rates of DCM oxidation to carbon monoxide (CO) in lung and liver and expected blood carboxyhemoglobin (HbCO) time courses with a PBPK model inclusive of both conjugation and oxidation pathways. Expression of large numbers of genes was altered at 100ppm with maximal changes in the numbers occurring by 500 or 2000ppm. Most changes in genes common to the two tissues were related to cellular metabolism and circadian clock. At the lower concentrations, the changes in metabolism-related genes were discordant - up in liver and down in lung. These processes included organelle biogenesis, TCA cycle, and respiratory electron transport. Changes in circadian cycle genes - primarily transcription factors - showed strong concentration-related response at higher concentrations (Arntl, Npas2, and Clock were down-regulated; Cry2, Wee1, Bhlhe40, Per3, Nr1d1, Nr1d2 and Dbp) were up-regulated with similar directionality in both tissues. Overall, persistently elevated HbCO from DCM oxidation appears to cause extended periods of hypoxia, leading to altered circadian coupling to cellular metabolism. The dose response for altered circadian processes correlates with the cancer outcome. We found no evidence of changes in genes indicative of responses to cytotoxic, DNA-reactive metabolites.

Multiple receptors shape the estrogen response pathway and are critical considerations for the future of in vitro-based risk assessment efforts.

Current in life toxicity testing paradigms are being challenged as the future of risk assessment moves towards more comprehensive mode of action/adverse outcome pathway based approaches. In particular, endocrine disruption screening is now a global activity and key initiatives in the United States focus on the use of high throughput in vitro assays to prioritize compounds for further testing of estrogen, androgen or thyroid disruption. Of these pathways, much of the emphasis to date has been on high-throughput methods for estrogenic activity primarily using ligand binding and trans-activation assays. However, as the knowledge regarding estrogen receptor signaling pathways continues to evolve, it is clear that the assumption of a simple one-receptor pathway underlying current in vitro screening assays is out of date. To develop more accurate models for estrogen-initiated pathways useful for quantitative safety assessments, we must design assays that account for the key signaling processes driving cellular dose response based on up-to-date understanding of the biological network. In this review, we summarize the state of the science for the estrogen receptor signaling network, particularly with regard to proliferative effects, and highlight gaps in current high throughput approaches. From the sum of this literature, we propose a model for the estrogen-signaling pathway that should serve as a starting point for development of new in vitro methods fit for the purpose of predicting dose response for estrogenic chemicals in the human.

Information-dependent enrichment analysis reveals time-dependent transcriptional regulation of the estrogen pathway of toxicity.

The twenty-first century vision for toxicology involves a transition away from high-dose animal studies to in vitro and computational models (NRC in Toxicity testing in the 21st century: a vision and a strategy, The National Academies Press, Washington, DC, 2007). This transition requires mapping pathways of toxicity by understanding how in vitro systems respond to chemical perturbation. Uncovering transcription factors/signaling networks responsible for gene expression patterns is essential for defining pathways of toxicity, and ultimately, for determining the chemical modes of action through which a toxicant acts. Traditionally, transcription factor identification is achieved via chromatin immunoprecipitation studies and summarized by calculating which transcription factors are statistically associated with up- and downregulated genes. These lists are commonly determined via statistical or fold-change cutoffs, a procedure that is sensitive to statistical power and may not be as useful for determining transcription factor associations. To move away from an arbitrary statistical or fold-change-based cutoff, we developed, in the context of the Mapping the Human Toxome project, an enrichment paradigm called information-dependent enrichment analysis (IDEA) to guide identification of the transcription factor network. We used a test case of activation in MCF-7 cells by 17ß estradiol (E2). Using this new approach, we established a time course for transcriptional and functional responses to E2. ERα and ERß were associated with short-term transcriptional changes in response to E2. Sustained exposure led to recruitment of additional transcription factors and alteration of cell cycle machinery. TFAP2C and SOX2 were the transcription factors most highly correlated with dose. E2F7, E2F1, and Foxm1, which are involved in cell proliferation, were enriched only at 24 h. IDEA should be useful for identifying candidate pathways of toxicity. IDEA outperforms gene set enrichment analysis (GSEA) and provides similar results to weighted gene correlation network analysis, a platform that helps to identify genes not annotated to pathways.

Applying 'omics technologies in chemicals risk assessment: Report of an ECETOC workshop.

Prevailing knowledge gaps in linking specific molecular changes to apical outcomes and methodological uncertainties in the generation, storage, processing, and interpretation of 'omics data limit the application of 'omics technologies in regulatory toxicology. Against this background, the European Centre for Ecotoxicology and Toxicology of Chemicals (ECETOC) convened a workshop Applying 'omics technologies in chemicals risk assessment that is reported herein. Ahead of the workshop, multi-expert teams drafted frameworks on best practices for (i) a Good-Laboratory Practice-like context for collecting, storing and curating 'omics data; (ii) the processing of 'omics data; and (iii) weight-of-evidence approaches for integrating 'omics data. The workshop participants confirmed the relevance of these Frameworks to facilitate the regulatory applicability and use of 'omics data, and the workshop discussions provided input for their further elaboration. Additionally, the key objective (iv) to establish approaches to connect 'omics perturbations to phenotypic alterations was addressed. Generally, it was considered promising to strive to link gene expression changes and pathway perturbations to the phenotype by mapping them to specific adverse outcome pathways. While further work is necessary before gene expression changes can be used to establish safe levels of substance exposure, the ECETOC workshop provided important incentives towards achieving this goal.