AOP report: Development of an adverse outcome pathway for deposition of energy leading to cataracts

Cataracts are one of the leading causes of blindness, with an estimated 95 million people affected worldwide. A hallmark of cataract development is lens opacification, typically associated not only with aging but also radiation exposure as encountered by interventional radiologists and astronauts during the long-term space mission. To better understand radiation-induced cataracts, the adverse outcome pathway (AOP) framework was used to structure and evaluate knowledge across biological levels of organization (e.g., macromolecular, cell, tissue, organ, organism and population). AOPs identify a sequence of key events (KEs) causally connected by key event relationships (KERs) beginning with a molecular initiating event to an adverse outcome (AO) of relevance to regulatory decision-making. To construct the cataract AO and retrieve evidence to support it, a scoping review methodology was used to filter, screen, and review studies based on the modified Bradford Hill criteria. Eight KEs were identified that were moderately supported by empirical evidence (e.g., dose-, time-, incidence-concordance) across the adjacent (directly linked) relationships using well-established endpoints. Over half of the evidence to justify the KER linkages was derived from the evidence stream of biological plausibility. Early KEs of oxidative stress and protein modifications had strong linkages to downstream KEs and could be the focus of countermeasure development. Several identified knowledge gaps and inconsistencies related to the quantitative understanding of KERs which could be the basis of future research, most notably directed to experiments in the range of low or moderate doses and dose-rates, relevant to radiation workers and other occupational exposures.

Empirical Characterization of False Discovery Rates of Differentially Expressed Genes and Transcriptomic Benchmark Concentrations in Zebrafish Embryos

High-throughput transcriptomics (HTTr) is increasingly applied to zebrafish embryos to survey the toxicological effects of environmental chemicals. Before the adoption of this approach in regulatory testing, it is essential to characterize background noise in order to guide experimental designs. We thus empirically quantified the HTTr false discovery rate (FDR) across different embryo pool sizes, sample sizes, and concentration groups for toxicology studies. We exposed zebrafish embryos to 0.1% dimethyl sulfoxide (DMSO) for 5 days. Pools of 1, 5, 10, and 20 embryos were created (n = 24 samples for each pool size). Samples were sequenced on the TempO-Seq platform and then randomly assigned to mock treatment groups before differentially expressed gene (DEG), pathway, and benchmark concentration (BMC) analyses. Given that all samples were treated with DMSO, any significant DEGs, pathways, or BMCs are false positives. As expected, we found decreasing FDRs for DEG and pathway analyses with increasing pool and sample sizes. Similarly, FDRs for BMC analyses decreased with increasing pool size and concentration groups, with more stringent BMC premodel filtering reducing BMC FDRs. Our study provides foundational data for determining appropriate experiment designs for regulatory toxicity testing with HTTr in zebrafish embryos.

Next Generation Sequencing Workshop at the Royal Society of Medicine (London, May 2022): how genomics is on the path to modernizing genetic toxicology

The use of error-corrected Next Generation Sequencing (ecNG) to determine mutagenicity has been a subject of growing interest and potentially a disruptive technology that could supplement, and in time, replace current testing paradigms in preclinical safety assessment. Considering this, a Next Generation Sequencing Workshop was held at the Royal Society of Medicine in London in May 2022, supported by the United Kingdom Environmental Mutagen Society (UKEMS) and TwinStrand Biosciences (WA, USA), to discuss progress and future applications of this technology. In this meeting report, the invited speakers provide an overview of the Workshop topics covered and identify future directions for research. In the area of somatic mutagenesis, several speakers reviewed recent progress made with correlating ecNGS to classic in vivo transgenic rodent mutation assays as well as exploring the use of this technology directly in humans and animals, and in complex organoid models. Additionally, ecNGS has been used for detecting off-target effects of gene editing tools and emerging data suggest ecNGS potential to measure clonal expansion of cells carrying mutations in cancer driver genes as an early marker of carcinogenic potential and for direct human biomonitoring. As such, the workshop demonstrated the importance of raising awareness and support for advancing the science of ecNGS for mutagenesis, gene editing, and carcinogenesis research. Furthermore, the potential of this new technology to contribute to advances in drug and product development and improve safety assessment was extensively explored.

Applying genomics in regulatory toxicology: a report of the ECETOC workshop on omics threshold on non-adversity

In a joint effort involving scientists from academia, industry and regulatory agencies, ECETOC's activities in Omics have led to conceptual proposals for: (1) A framework that assures data quality for reporting and inclusion of Omics data in regulatory assessments; and (2) an approach to robustly quantify these data, prior to interpretation for regulatory use. In continuation of these activities this workshop explored and identified areas of need to facilitate robust interpretation of such data in the context of deriving points of departure (POD) for risk assessment and determining an adverse change from normal variation. ECETOC was amongst the first to systematically explore the application of Omics methods, now incorporated into the group of methods known as New Approach Methodologies (NAMs), to regulatory toxicology. This support has been in the form of both projects (primarily with CEFIC/LRI) and workshops. Outputs have led to projects included in the workplan of the Extended Advisory Group on Molecular Screening and Toxicogenomics (EAGMST) group of the Organisation for Economic Co-operation and Development (OECD) and to the drafting of OECD Guidance Documents for Omics data reporting, with potentially more to follow on data transformation and interpretation. The current workshop was the last in a series of technical methods development workshops, with a sub-focus on the derivation of a POD from Omics data. Workshop presentations demonstrated that Omics data developed within robust frameworks for both scientific data generation and analysis can be used to derive a POD. The issue of noise in the data was discussed as an important consideration for identifying robust Omics changes and deriving a POD. Such variability or "noise" can comprise technical or biological variation within a dataset and should clearly be distinguished from homeostatic responses. Adverse outcome pathways (AOPs) were considered a useful framework on which to assemble Omics methods, and a number of case examples were presented in illustration of this point. What is apparent is that high dimension data will always be subject to varying processing pipelines and hence interpretation, depending on the context they are used in. Yet, they can provide valuable input for regulatory toxicology, with the pre-condition being robust methods for the collection and processing of data together with a comprehensive description how the data were interpreted, and conclusions reached.

Per- and Polyfluoroalkyl Substances (PFAS) in mixtures show additive effects on transcriptomic points of departure in human liver spheroids

Per- and polyfluoroalkyl substances (PFAS) are a wide range of chemicals that are used in a variety of consumer and industrial products leading to direct human exposure. Many PFAS are chemically non-reactive and persistent in the environment, resulting in additional exposure from water, soil, and dietary intake. While some PFAS have documented negative health effects, data on simultaneous exposures to multiple PFAS (PFAS mixtures) are inadequate for making informed decisions for risk assessment. The current study leverages data from previous work in our group using Templated Oligo-Sequencing (TempO-Seq™) for high-throughput transcriptomic analysis of PFAS-exposed primary human liver cell spheroids; herein, we determine the transcriptomic potency of PFAS in mixtures. Gene expression data from single PFAS and mixture exposures of liver cell spheroids were subject to benchmark concentration (BMC) analysis. We used the 25th lowest gene BMC as the point of departure to compare the potencies of single PFAS to PFAS mixtures of varying complexity and composition. Specifically, the empirical potency of eight PFAS mixtures were compared to predicted mixture potencies calculated using the principal of concentration addition (i.e., dose addition) in which mixture component potencies are summed by proportion to predict mixture potency. In this study, for most mixtures, empirical mixture potencies were comparable to potencies calculated through concentration addition. This work supports that the effects of PFAS mixtures on gene expression largely follow the concentration addition predicted response and suggests that effects of these individual PFAS in mixtures are not strongly synergistic or antagonistic.

A Transformative Vision for an Omics-Based Regulatory Chemical Testing Paradigm

Use of molecular data in human and ecological health risk assessments of industrial chemicals and agrochemicals has been anticipated by the scientific community for many years; however, these data are rarely used for risk assessment. Here, a logic framework is proposed to explore the feasibility and future development of transcriptomic methods to refine and replace the current apical endpoint-based regulatory toxicity testing paradigm. Four foundational principles are outlined and discussed that would need to be accepted by stakeholders prior to this transformative vision being realized. Well-supported by current knowledge, the first principal is that transcriptomics is a reliable tool for detecting alterations in gene expression that result from endogenous or exogenous influences on the test organism. The second principle states that alterations in gene expression are indicators of adverse or adaptive biological responses to stressors in an organism. Principle three is that transcriptomics can be employed to establish a benchmark dose-based point of departure (POD) from short-term, in vivo studies at a dose level below which a concerted molecular change (CMC) is not expected. Finally, principle four states that the use of a transcriptomic POD (set at the CMC dose level) will support a human health-protective risk assessment. If all four principles are substantiated, this vision is expected to transform aspects of the industrial chemical and agrochemical risk assessment process that are focused on establishing safe exposure levels for mammals across numerous toxicological contexts resulting in a significant reduction in animal use while providing equal or greater protection of human health. Importantly, these principals and approach are also generally applicable for ecological safety assessment.

Evaluating the Influences of Confounding Variables on Benchmark Dose using a Case Study in the Field of Ionizing Radiation

Purpose A vast amount of data regarding the effects of radiation stressors on transcriptional changes has been produced over the past few decades. These data have shown remarkable consistency across platforms and experimental design, enabling increased understanding of early molecular effects of radiation exposure. However, the value of transcriptomic data in the context of risk assessment is not clear and represents a gap that is worthy of further consideration. Recently, benchmark dose (BMD) modeling has shown promise in correlating a transcriptional point of departure (POD) to that derived using phenotypic outcomes relevant to human health risk assessment. Although frequently applied in chemical toxicity evaluation, our group has recently demonstrated application within the field of radiation research. This approach allows the possibility to quantitatively compare radiation-induced gene and pathway alterations across various datasets using BMD values and derive meaningful biological effects. However, before BMD modeling can confidently be used, an understanding of the impact of confounding variables on BMD outputs is needed. Methods: To this end, BMD modeling was applied to a publicly available microarray dataset (Gene Expression Omnibus #GSE23515) that used peripheral blood ex-vivo gamma-irradiated at 0.82 Gy/min, at doses of 0, 0.1, 0.5 or 2 Gy, and assessed 6 hours post-exposure. The dataset comprised six female smokers (F-S), six female non-smokers (F-NS), six male smokers (M-S), and six male non-smokers (M-NS). Results: A combined total of 412 genes were fit to models and the BMD distribution was noted to be bi-modal across the four groups. A total of 74, 41, 62 and 62 genes were unique to the F-NS, M-NS, F-S and M-S groups. Sixty-two BMD modeled genes and nine pathways were common across all four groups. There were no differential sensitivity of responses in the robust common genes and pathways. Conclusion: For radiation-responsive genes and pathways common across the study groups, the BMD distribution of transcriptional activity was unaltered by sex and smoking status. Although further validation of the data is needed, these initial findings suggest BMD values for radiation relevant genes and pathways are robust and could be explored further in future studies.

Deploying elements of scoping review methods for Adverse Outcome Pathway development: A space travel case example

Purpose Health protection agencies require scientific information for evidence-based decision-making and guideline development. However, vetting and collating large quantities of published research to identify relevant high-quality studies is a challenge. One approach to address this issue is the use of Adverse Outcome Pathways (AOPs) that provide a framework to assemble toxicological knowledge into causally linked chains of key events across levels of biological organization to culminate in an adverse health outcome of significance. Traditionally, AOPs have been constructed using a narrative review approach where the collection of evidence that supports each pathway is based on prior knowledge of influential studies that can also be supplemented by individually selecting and reviewing relevant references. Objectives: We aimed to create a protocol for AOP weight of evidence gathering that harnesses elements of both scoping review methods and artificial intelligence tools to increase transparency while reducing bias and workload of human screeners. Methods: To develop this protocol, an existing space-health AOP in the workplan of the Organisation for Economic Co-operation and Development (OECD) AOP program was used as a case example. To balance the benefits of both scoping review tools and narrative approaches, a study protocol outlining a screening and search strategy was developed, and three reference collection workflows were tested to identify the most efficient method to inform weight of evidence. The workflows differed in their literature search strategies, and combinations of software tools used. Results: Across the three tested workflows, over 59 literature searches were completed, retrieving over 34000 references of which over 3300 were human reviewed. The most effective of the three methods used a search strategy with searches across each component of the AOP network, SWIFT Review as a pre-filtering software, and DistillerSR to create structured screening and data extraction forms. This methodology effectively retrieved relevant studies while balancing efficiency in data retrieval without compromising transparency, leading to a well-synthesized evidence base to support the AOP. Conclusions: The workflow is still exploratory in the context of AOP development, and we anticipate adaptations to the protocol with further experience. To further the systematicity, future iterations of the workflow could include structured quality assessment and risk of bias analysis. Overall, the workflow provides a transparent and unbiased approach to support AOP development, which in turn will support the need for rigorous methods to identify relevant scientific evidence while being practical to allow uptake by the broader community.

Benchmark dose modeling of transcriptional data: a systematic approach to identify best practices for study designs used in radiation research

PURPOSE

Benchmark dose (BMD) modeling is a method commonly used in chemical toxicology to identify the point of departure (POD) from a dose-response curve linked to a health-related outcome. Recently, it is being explored on transcriptional data and in adverse outcome pathways (AOPs). As AOPs are informed by diverse data types, it is important to understand the impact of study parameters such as dose selection, number of replicates and dose range on BMD outputs for radiation induced genes and pathways.

MATERIALS AND METHODS

Data were selected from the Gene Expression Omnibus (GSE52403) that featured gene expression profiles of peripheral blood samples from C57BL/6 mice 6 hours post-exposure to ¹³⁷Cs gamma-radiation at 0, 1, 2, 3, 4.5, 6, 8 and 10.5 Gy. The dataset comprised a broad dose-range over multiple dose-points with consistent dose spacing and multiple biological replicates. This dataset was ideal for systematically transforming across three categories: (1) dose-range, (2) dose-spacing and (3) number of controls/replicates. Across these categories, 29 transformed datasets were compared to the original dataset to determine the impact of each transformation on the BMD outputs.

RESULTS

Most of the experimental changes did not impact the BMD outputs. The transformed datasets were largely consistent with the original dataset in terms of number of reproduced genes modeled and absolute BMD values for genes and pathways. Variations in dose selection identified the importance of the absolute value of the lowest and second dose. It was determined that dose selection should include at least two doses <1 Gy and two >5 Gy to achieve meaningful BMD outputs. Changes to the number of biological replicates in the control and non-zero dose groups impacted the overall accuracy and precision of the BMD outputs as well as the ability to fit dose-response models consistent with the original dataset.

CONCLUSION

Successful application of transcriptomic BMD modeling for radiation datasets requires considerations of the exposure dose and the number of biological replicates. Most important is the selection of the lowest doses and dose spacing. Reflections on these parameters in experimental design will provide meaningful BMD outputs that could correlate well to apical endpoints of relevance to radiation exposure assessment.

Reactive Oxygen Species in the Adverse Outcome Pathway Framework: Toward Creation of Harmonized Consensus Key Events

Reactive oxygen species (ROS) and reactive nitrogen species (RNS) are formed as a result of natural cellular processes, intracellular signaling, or as adverse responses associated with diseases or exposure to oxidizing chemical and non-chemical stressors. The action of ROS and RNS, collectively referred to as reactive oxygen and nitrogen species (RONS), has recently become highly relevant in a number of adverse outcome pathways (AOPs) that capture, organize, evaluate and portray causal relationships pertinent to adversity or disease progression. RONS can potentially act as a key event (KE) in the cascade of responses leading to an adverse outcome (AO) within such AOPs, but are also known to modulate responses of events along the AOP continuum without being an AOP event itself. A substantial discussion has therefore been undertaken in a series of workshops named "Mystery or ROS" to elucidate the role of RONS in disease and adverse effects associated with exposure to stressors such as nanoparticles, chemical, and ionizing and non-ionizing radiation. This review introduces the background for RONS production, reflects on the direct and indirect effects of RONS, addresses the diversity of terminology used in different fields of research, and provides guidance for developing a harmonized approach for defining a common event terminology within the AOP developer community.

Brief Developmental Exposure to Fluoxetine Causes Life-Long Alteration of the Brain Transcriptome in Zebrafish

Fluoxetine (FLX) and other selective serotonin reuptake inhibitors are widely used to treat depressive disorders during pregnancy. Early-life exposure to FLX is known to disrupt the normal function of the stress axis in humans, rodents, and teleosts. We used a zebrafish line with a cortisol-inducible fluorescent transgene to study the effects of developmental daily exposure to FLX (54 µg/L) on the transcriptomic profile of brain tissues in exposed larvae and later as 6-month-old adults. High throughput RNA sequencing was conducted on brain tissues in unstressed and stressed conditions. Long-lasting effects of FLX were observed in telencephalon (Tel) and hypothalamus (Hyp) of adult zebrafish with 1927 and 5055 genes significantly (≥1.2 fold-change, false-discovery p-value < 0.05) dysregulated in unstressed condition, respectively. Similar findings were observed in Hyp with 1245 and 723 genes being significantly dysregulated in stressed adults, respectively. Differentially expressed genes converted to Homo sapiens orthologues were used for Ingenuity Pathway Analysis. The results showed alteration of pathways involved in neuroendocrine signaling, cholesterol metabolism and synaptogenesis. Enriched networks included lipid metabolism, molecular transport, and nervous system development. Analysis of putative upstream transcription regulators showed potential dysregulation of clocka and nr3c1 which control circadian rhythm, stress response, cholesterol metabolism and histone modifications. Several genes involved in epigenetic regulation were also affected by FLX, including dnmt3a, adarb1, adarb2, hdac4, hdac5, hdac8, and atf2. We report life-long disruptive effects of FLX on pathways associated with neuroendocrine signaling, stress response and the circadian rhythm, and all of which are implicated in the development of depressive disorders in humans. Our results raise concern for the persistent endocrine-disrupting potential of brief antidepressant exposure during embryonic development.

High-Throughput Transcriptomic Analysis of Human Primary Hepatocyte Spheroids Exposed to Per- and Polyfluoroalkyl Substances as a Platform for Relative Potency Characterization

Per- and poly-fluoroalkyl substances (PFAS) are widely found in the environment because of their extensive use and persistence. Although several PFAS are well studied, most lack toxicity data to inform human health hazard and risk assessment. This study focused on 4 model PFAS: perfluorooctanoic acid (PFOA; 8 carbon), perfluorobutane sulfonate (PFBS; 4 carbon), perfluorooctane sulfonate (PFOS; 8 carbon), and perfluorodecane sulfonate (PFDS; 10 carbon). Human primary liver cell spheroids (pooled from 10 donors) were exposed to 10 concentrations of each PFAS and analyzed at 4 time points. The approach aimed to: (1) identify gene expression changes mediated by the PFAS, (2) identify similarities in biological responses, (3) compare PFAS potency through benchmark concentration analysis, and (4) derive bioactivity exposure ratios (ratio of the concentration at which biological responses occur, relative to daily human exposure). All PFAS induced transcriptional changes in cholesterol biosynthesis and lipid metabolism pathways, and predicted PPARα activation. PFOS exhibited the most transcriptional activity and had a highly similar gene expression profile to PFDS. PFBS induced the least transcriptional changes and the highest benchmark concentration (ie, was the least potent). The data indicate that these PFAS may have common molecular targets and toxicities, but that PFOS and PFDS are the most similar. The transcriptomic bioactivity exposure ratios derived here for PFOA and PFOS were comparable to those derived using rodent apical endpoints in risk assessments. These data provide a baseline level of toxicity for comparison with other known PFAS using this testing strategy.

Bringing together scientific disciplines for collaborative undertakings: a vision for advancing the adverse outcome pathway framework

BACKGROUND

Decades of research to understand the impacts of various types of environmental occupational and medical stressors on human health have produced a vast amount of data across many scientific disciplines. Organizing these data in a meaningful way to support risk assessment has been a significant challenge. To address this and other challenges in modernizing chemical health risk assessment, the Organisation for Economic Cooperation and Development (OECD) formalized the adverse outcome pathway (AOP) framework, an approach to consolidate knowledge into measurable key events (KEs) at various levels of biological organisation causally linked to disease based on the weight of scientific evidence (http://oe.cd/aops). Currently, AOPs have been considered predominantly in chemical safety but are relevant to radiation. In this context, the Nuclear Energy Agency's (NEA's) High-Level Group on Low Dose Research (HLG-LDR) is working to improve research co-ordination, including radiological research with chemical research, identify synergies between the fields and to avoid duplication of efforts and resource investments. To this end, a virtual workshop was held on 7 and 8 October 2020 with experts from the OECD AOP Programme together with the radiation and chemical research/regulation communities. The workshop was a coordinated effort of Health Canada, the Electric Power Research Institute (EPRI), and the Nuclear Energy Agency (NEA). The AOP approach was discussed including key issues to fully embrace its value and catalyze implementation in areas of radiation risk assessment.

CONCLUSIONS

A joint chemical and radiological expert group was proposed as a means to encourage cooperation between risk assessors and an initial vision was discussed on a path forward. A global survey was suggested as a way to identify priority health outcomes of regulatory interest for AOP development. Multidisciplinary teams are needed to address the challenge of producing the appropriate data for risk assessments. Data management and machine learning tools were highlighted as a way to progress from weight of evidence to computational causal inference.

Transcriptomic pathway and benchmark dose analysis of Bisphenol A, Bisphenol S, Bisphenol F, and 3,3',5,5'-Tetrabromobisphenol A in H9 human embryonic stem cells

Bisphenol A (BPA) is a chemical used in the manufacturing of plastics to which human exposure is ubiquitous. Numerous studies have linked BPA exposure to many adverse health outcomes prompting the replacement of BPA with various analogues including bisphenol-F (BPF) and bisphenol S (BPS). Other bisphenols are used in various consumer applications, such as 3,3',5,5'-Tetrabromobisphenol A (TBBPA), which is used as a flame retardant. Few studies to date have examined the effects of BPA and its analogues in stem cells to explore potential developmental impacts. Here we used transcriptomics to investigate similarities and differences of BPA and three of its analogues in the estrogen receptor negative, human embryonic stem cell line H9 (WA09). H9 cells were exposed to increasing concentrations of the bisphenols and analyzed using RNA-sequencing. Our data indicate that BPA, BPF, and BPS have similar potencies in inducing transcriptional changes and perturb many of the same pathways. TBBPA, the least structurally similar bisphenol of the group, exhibited much lower potency. All bisphenols robustly impacted gene expression in these cells, albeit at concentrations well above those observed in estrogen-positive cells. Overall, we provide a foundational data set against which to explore the transcriptional similarities of other bisphenols in embryonic stem cells, which may be used to assess the suitability of chemical grouping for read-across and for preliminary potency evaluation.

Demodifying RNA for Transcriptomic Analyses of Archival Formalin-Fixed Paraffin-Embedded Samples

Archival formalin-fixed paraffin-embedded (FFPE) tissue samples offer a vast but largely untapped resource for genomic research. The primary technical issues limiting use of FFPE samples are RNA yield and quality. In this study, we evaluated methods to demodify RNA highly fragmented and crosslinked by formalin fixation. Primary endpoints were RNA recovery, RNA-sequencing quality metrics, and transcriptional responses to a reference chemical (phenobarbital, PB). Frozen mouse liver samples from control and PB groups (n = 6/group) were divided and preserved for 3 months as follows: frozen (FR); 70% ethanol (OH); 10% buffered formalin for 18 h followed by ethanol (18F); or 10% buffered formalin (3F). Samples from OH, 18F, and 3F groups were processed to FFPE blocks and sectioned for RNA isolation. Additional sections from 3F received the following demodification protocols to mitigate RNA damage: short heated incubation with Tris-Acetate-EDTA buffer; overnight heated incubation with an organocatalyst using 2 different isolation kits; or overnight heated incubation without organocatalyst. Ribo-depleted, stranded, total RNA libraries were built and sequenced using the Illumina HiSeq 2500 platform. Overnight incubation (± organocatalyst) increased RNA yield >3-fold and RNA integrity numbers and fragment analysis values by > 1.5- and >3.0-fold, respectively, versus 3F. Postsequencing metrics also showed reduced bias in gene coverage and deletion rates for overnight incubation groups. All demodification groups had increased overlap for differentially expressed genes (77%-84%) and enriched pathways (91%-97%) with FR, with the highest overlap in the organocatalyst groups. These results demonstrate simple changes in RNA isolation methods that can enhance genomic analyses of FFPE samples.

Is there a role for the adverse outcome pathway framework to support radiation protection?

PURPOSE

In 2012, the Organization for Economic Cooperation and Development (OECD) formally launched the Adverse Outcome Pathway (AOP) Programme. The AOP framework has the potential for predictive utility in identifying early biological endpoints linked to adverse effects. It uses the weight of correlative evidence to identify a minimal set of measurable key events that link molecular initiating events to an adverse outcome. AOPs have the capability to identify knowledge gaps and priority areas for future research based on relevance to an adverse outcome. In addition, AOPs can identify pathways that are common among multiple stressors, thereby allowing for the possibility of refined risk assessments based on co-exposure considerations. The AOP framework is increasingly being used in chemical and ecological risk assessment; however, its use in the development of radiation-specific pathways has yet to be fully explored. To bring awareness of the AOP framework to the Canadian radiation community, a workshop was held in Canada in June 2018 that brought together radiation experts from Health Canada, the Canadian Nuclear Laboratories, and the Canadian Nuclear Safety Commission.

METHODS

The purpose of the workshop was to share knowledge on the AOP framework, specifically (1) to introduce the concept of the AOP framework and its possible utility to Canadian radiation experts; (2) to provide examples on how it has advanced risk assessment; (3) to discuss an illustrative example specific to ionizing radiation; and lastly (4) to identify the broad benefits and challenges of the AOP framework to the radiation community.

RESULTS

The participants showed interest in the framework, case examples were described and areas of challenge were identified. Herein, we summarize the outcomes of the workshop.

CONCLUSIONS

Overall, participants agreed that by building AOPs in the radiation field, a network of data-sharing initiatives will enhance our interpretation of existing knowledge where current scientific evidence is minimal. They would provide new avenues to understand effects at low-dose and dose-rates and help to quantify the combined effect of multiple stressors on shared mechanistic pathways.

Framework for the quality assurance of 'omics technologies considering GLP requirements

'Omics technologies are gaining importance to support regulatory toxicity studies. Prerequisites for performing 'omics studies considering GLP principles were discussed at the European Centre for Ecotoxicology and Toxicology of Chemicals (ECETOC) Workshop Applying 'omics technologies in Chemical Risk Assessment. A GLP environment comprises a standard operating procedure system, proper pre-planning and documentation, and inspections of independent quality assurance staff. To prevent uncontrolled data changes, the raw data obtained in the respective 'omics data recording systems have to be specifically defined. Further requirements include transparent and reproducible data processing steps, and safe data storage and archiving procedures. The software for data recording and processing should be validated, and data changes should be traceable or disabled. GLP-compliant quality assurance of 'omics technologies appears feasible for many GLP requirements. However, challenges include (i) defining, storing, and archiving the raw data; (ii) transparent descriptions of data processing steps; (iii) software validation; and (iv) ensuring complete reproducibility of final results with respect to raw data. Nevertheless, 'omics studies can be supported by quality measures (e.g., GLP principles) to ensure quality control, reproducibility and traceability of experiments. This enables regulators to use 'omics data in a fit-for-purpose context, which enhances their applicability for risk assessment.

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.

The Next Generation of Risk Assessment Multi-Year Study-Highlights of Findings, Applications to Risk Assessment, and Future Directions

BACKGROUND

The Next Generation (NexGen) of Risk Assessment effort is a multi-year collaboration among several organizations evaluating new, potentially more efficient molecular, computational, and systems biology approaches to risk assessment. This article summarizes our findings, suggests applications to risk assessment, and identifies strategic research directions.

OBJECTIVE

Our specific objectives were to test whether advanced biological data and methods could better inform our understanding of public health risks posed by environmental exposures.

METHODS

New data and methods were applied and evaluated for use in hazard identification and dose-response assessment. Biomarkers of exposure and effect, and risk characterization were also examined. Consideration was given to various decision contexts with increasing regulatory and public health impacts. Data types included transcriptomics, genomics, and proteomics. Methods included molecular epidemiology and clinical studies, bioinformatic knowledge mining, pathway and network analyses, short-duration in vivo and in vitro bioassays, and quantitative structure activity relationship modeling.

DISCUSSION

NexGen has advanced our ability to apply new science by more rapidly identifying chemicals and exposures of potential concern, helping characterize mechanisms of action that influence conclusions about causality, exposure-response relationships, susceptibility and cumulative risk, and by elucidating new biomarkers of exposure and effects. Additionally, NexGen has fostered extensive discussion among risk scientists and managers and improved confidence in interpreting and applying new data streams.

CONCLUSIONS

While considerable uncertainties remain, thoughtful application of new knowledge to risk assessment appears reasonable for augmenting major scope assessments, forming the basis for or augmenting limited scope assessments, and for prioritization and screening of very data limited chemicals. Citation: Cote I, Andersen ME, Ankley GT, Barone S, Birnbaum LS, Boekelheide K, Bois FY, Burgoon LD, Chiu WA, Crawford-Brown D, Crofton KM, DeVito M, Devlin RB, Edwards SW, Guyton KZ, Hattis D, Judson RS, Knight D, Krewski D, Lambert J, Maull EA, Mendrick D, Paoli GM, Patel CJ, Perkins EJ, Poje G, Portier CJ, Rusyn I, Schulte PA, Simeonov A, Smith MT, Thayer KA, Thomas RS, Thomas R, Tice RR, Vandenberg JJ, Villeneuve DL, Wesselkamper S, Whelan M, Whittaker C, White R, Xia M, Yauk C, Zeise L, Zhao J, DeWoskin RS. 2016. The Next Generation of Risk Assessment multiyear study-highlights of findings, applications to risk assessment, and future directions. Environ Health Perspect 124:1671-1682; http://dx.doi.org/10.1289/EHP233.

Abstract 5226: Liver and blood miRNA alterations as putative biomarkers of hepatotoxic response to short-term furan exposure in mice

Epigenetic miRNA-based changes measured in accessible matrices may serve as useful biomarkers of environmental exposures and human health effects. We investigated miRNA profiles following short-term exposure to a known cytotoxic hepatocarcinogen, furan. We measured global liver and blood miRNA changes in female B3C6F1 mice exposed to furan for 3 weeks p.o. at daily doses of 0, 1, 2, 4, and 8 mg/kg. Small RNA was extracted from frozen samples using Exiqon miRCURY RNA Isolation kit and quantified by Qubit. RNA quality was checked by Nanodrop and Bioanalyzer, and liver miRNA were measured by miRNA-seq on the Illumina HiScanSQ platform. Blood miRNA were measured from pooled RNA samples by Affymetrix GeneChip v3 microarrays. In the liver, the 1, 2, 4, and 8 mg/kg doses resulted in 0, 10, 17 and 14 differentially altered miRNAs, respectively (&gt;1.5-fold; B-H corrected p-value&lt;0.05). The majority of miRNA changes in the low-dose exposure group (2mg/kg) were not seen with other treatments (7 of 10 miRNA). Conversely, most altered miRNA observed with high-dose (4 or 8 mg/kg) exposures were shared (10 miRNA), indicating different miRNA-mediated pathways for non-carcinogenic and carcinogenic furan exposures. Using Ingenuity Pathway Analysis (IPA), liver miRNAs altered by the carcinogenic doses of furan (4 and 8 mg/kg) significantly enriched (B-H corrected p-value&lt;0.05) functions related to cellular development, growth and proliferation, movement, cell cycle, death and survival; other categories identified by IPA were related to hepatotoxicity, liver inflammation, and cancer. Using predicted miRNA-mRNA interactions, 68 mRNAs were expressed in a contrasting pattern to miRNA expression (upregulated miRNA/downregulated mRNAs and vice versa). Associated mRNAs showed enrichment of pathways similar to miRNAs. Measurements in whole blood samples indicated that 6 out of the 21 miRNAs altered by the carcinogenic doses of furan were also altered in the blood (4 in the same direction). Blood miRNAs included mmu-miR-34a, -146b, -183, -5099 (upregulated), and -10a, -99b (downregulated). Results demonstrate distinctive miRNA profiles in rodent liver for carcinogenic doses of furan with corresponding changes for a subset of liver miRNAs in blood. Our findings support ongoing efforts to identify novel miRNA biomarkers in accessible matrices related to environmental health effects. This abstract does not necessarily reflect the policy of the US EPA.

Citation Format: Brian Norris Chorley, Gail Nelson, Gleta Carswell, Holly Mortensen, James Crooks, William Ward, Charles Wood, Anna F. Jackson, Carole Yauk, Les Recio, Susan Hester. Liver and blood miRNA alterations as putative biomarkers of hepatotoxic response to short-term furan exposure in mice. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 5226. doi:10.1158/1538-7445.AM2014-5226

Expression of FUS-CHOP fusion protein in immortalized/transformed human mesenchymal stem cells drives mixoid liposarcoma formation

Increasing evidence supports that mesenchymal stromal/stem cells (MSCs) may represent the target cell for sarcoma development. Although different sarcomas have been modeled in mice upon expression of fusion oncogenes in MSCs, sarcomagenesis has not been successfully modeled in human MSCs (hMSCs). We report that FUS-CHOP, a hallmark fusion gene in mixoid liposarcoma (MLS), has an instructive role in lineage commitment, and its expression in hMSC sequentially immortalized/transformed with up to five oncogenic hits (p53 and Rb deficiency, hTERT over-expression, c-myc stabilization, and H-RAS(v12) mutation) drives the formation of serially transplantable MLS. This is the first model of sarcoma based on the expression of a sarcoma-associated fusion protein in hMSC, and allowed us to unravel the differentiation processes and signaling pathways altered in the MLS-initiating cells. This study will contribute to test novel therapeutic approaches and constitutes a proof-of-concept to use hMSCs as target cell for modeling other fusion gene-associated human sarcomas.

Functional analysis of the TRIB1 associated locus linked to plasma triglycerides and coronary artery disease

Background

The TRIB1 locus has been linked to hepatic triglyceride metabolism in mice and to plasma triglycerides and coronary artery disease in humans. The lipid‐associated single nucleotide polymorphisms (SNPs), identified by genome‐wide association studies, are located ≈30 kb downstream from TRIB1, suggesting complex regulatory effects on genes or pathways relevant to hepatic triglyceride metabolism. The goal of this study was to investigate the functional relationship between common SNPs at the TRIB1 locus and plasma lipid traits.

Methods and Results

Characterization of the risk locus reveals that it encompasses a gene, TRIB1‐associated locus (TRIBAL), composed of a well‐conserved promoter region and an alternatively spliced transcript. Bioinformatic analysis and resequencing identified a single SNP, rs2001844, within the promoter region that associates with increased plasma triglycerides and reduced high‐density lipoprotein cholesterol and coronary artery disease risk. Further, correction for triglycerides as a covariate indicated that the genome‐wide association studies association is largely dependent on triglycerides. In addition, we show that rs2001844 is an expression trait locus (eQTL) for TRIB1 expression in blood and alters TRIBAL promoter activity in a reporter assay model. The TRIBAL transcript has features typical of long noncoding RNAs, including poor sequence conservation. Modulation of TRIBAL expression had limited impact on either TRIB1 or lipid regulatory genes mRNA levels in human hepatocyte models. In contrast, TRIB1 knockdown markedly increased TRIBAL expression in HepG2 cells and primary human hepatocytes.

Conclusions

These studies demonstrate an interplay between a novel locus, TRIBAL, and TRIB1. TRIBAL is located in the genome‐wide association studies identified risk locus, responds to altered expression of TRIB1, harbors a risk SNP that is an eQTL for TRIB1 expression, and associates with plasma triglyceride concentrations.

Impact of cigarette smoke on the human and mouse lungs: a gene-expression comparison study

Cigarette smoke is well known for its adverse effects on human health, especially on the lungs. Basic research is essential to identify the mechanisms involved in the development of cigarette smoke-related diseases, but translation of new findings from pre-clinical models to the clinic remains difficult. In the present study, we aimed at comparing the gene expression signature between the lungs of human smokers and mice exposed to cigarette smoke to identify the similarities and differences. Using human and mouse whole-genome gene expression arrays, changes in gene expression, signaling pathways and biological functions were assessed. We found that genes significantly modulated by cigarette smoke in humans were enriched for genes modulated by cigarette smoke in mice, suggesting a similar response of both species. Sixteen smoking-induced genes were in common between humans and mice including six newly reported to be modulated by cigarette smoke. In addition, we identified a new conserved pulmonary response to cigarette smoke in the induction of phospholipid metabolism/degradation pathways. Finally, the majority of biological functions modulated by cigarette smoke in humans were also affected in mice. Altogether, the present study provides information on similarities and differences in lung gene expression response to cigarette smoke that exist between human and mouse. Our results foster the idea that animal models should be used to study the involvement of pathways rather than single genes in human diseases.

Paternal lifestyle as a potential source of germline mutations transmitted to offspring

Paternal exposure to high levels of radioactivity causes heritable germline minisatellite mutations. However, the effect of more general paternal exposures, such as cigarette smoking, on germline mutations remains unexplored. We analyzed two of the most commonly used minisatellite loci (CEB1 and B6.7) to identify germline mutations in blood samples of complete mother-father-child triads from the Norwegian Mother and Child Cohort Study (MoBa). The presence of mutations was subsequently related to general lifestyle factors, including paternal smoking before the partner became pregnant. Paternally derived mutations at the B6.7 locus (mutation frequency 0.07) were not affected by lifestyle. In contrast, high gross yearly income as a general measure of a healthy lifestyle coincided with low-mutation frequencies at the CEB1 locus (P=0.047). Income was inversely related to smoking behavior, and paternally derived CEB1 mutations were dose dependently increased when the father smoked in the 6 mo before pregnancy, 0.21 vs. 0.05 in smoking and nonsmoking fathers, respectively (P=0.061). These results suggest that paternal lifestyle can affect the chance of heritable mutations in unstable repetitive DNA sequences. To our knowledge, this is the first study reporting an effect of lifestyle on germline minisatellite mutation frequencies in a human population with moderate paternal exposures.

Subchronic oral exposure to benzo(a)pyrene leads to distinct transcriptomic changes in the lungs that are related to carcinogenesis

We have previously shown that acute oral exposure to the environmental carcinogen benzo(a)pyrene (BaP) elicits comparable levels of DNA adducts, but distinct transcriptomic changes, in mouse lungs and livers, the two main BaP bioactivating organs. Oral BaP exposure is predominantly associated with lung cancer and not hepatic cancer in some animal models, suggesting that gene expression differences may provide insight into the drivers of tissue-specific carcinogenesis. In the present study, we examine pulmonary DNA adduct formation, lacZ mutant frequency, and mRNA profiles in adult male MutaMouse following subchronic (28 day) oral exposure to BaP (0, 25, 50, and 75 mg/kg/day) and sacrificed 3 days postexposure. The results are compared with those obtained from livers of the same mice (previously published). Although there was a 1.8- to 3.3-fold increase in the levels of DNA adducts in lung compared with liver, the lacZ transgene mutant frequency was similar in both tissues. At the transcriptomic level, a transition from activation of the DNA damage response p53 pathway at the low dose to the induction of genes involved in angiogenesis, evasion of apoptosis and growth signals at the high doses was evident only in the lungs. These results suggest that tissue DNA adducts and mutant frequency are sensitive markers of target tissue exposure and mode of action, whereas early changes in gene expression may provide a better indication of the likelihood of carcinogenic transformation in selected tissues. Moreover, the study provides new information on the underlying mechanisms that contribute to tissue-specific responses to BaP.

Abstract 2550: Expanded simple tandem repeat (ESTR) mutations in offspring of benzo(a)pyrene exposed DNA repair deficient male mice

Germ line mutations can be transmitted to the offspring and may thus affect an infinite number of generations thereafter. To study germ line mutations, expanded simple tandem repeats (ESTR) can be used as a sensitive tool, because they are known to be highly unstable and therefore exhibit high mutation rates. This facilitates the assessment of germ line mutations in a relatively small number of offspring following exposures that may pose a genetic risk to the germ line. We used two ESTR loci Ms6-hm and Hm-2 in families of benzo(a)pyrene (BaP) exposed and unexposed male mice, proficient or deficient for global genome repair (Xpc-/-). Male mice (C57BL/6) were exposed to BaP during 6 weeks (13 mg/kg bw, 3 times per week) and were crossed with female mice (Balb/C) 6 weeks after the last exposure to obtain offspring and to ensure the analysis of mutations originating from spermatogonial stem cells. The overall ESTR mutation rate at both ESTR loci, irrespective of the DNA repair capacity, was higher in offspring of benzo(a)pyrene exposed mice than in control mice (9.4% vs. 4.5%, P=0.054). This effect was predominantly seen in BaP exposed DNA repair deficient Xpc-/- mice; mutation frequencies were 12.5% in the exposed group (12 out of 96 offspring carried a mutant allele) and 3.4% in the unexposed group (3 out of 87), respectively (P=0.024). In DNA repair proficient mice, BaP was unable to increase the mutation frequency transmitted to the offspring (5.8% vs. 5.5% for offspring of exposed and unexposed fathers). These results indicate that analysis of tandem repetitive sequences is a sensitive method for the detection of germline mutations transmitted to the offspring and that the environmental mutagen BaP is able to induce germ line mutations in DNA repair deficient mice.

Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 103rd Annual Meeting of the American Association for Cancer Research; 2012 Mar 31-Apr 4; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2012;72(8 Suppl):Abstract nr 2550. doi:1538-7445.AM2012-2550

Gene expression analysis to identify molecular correlates of pre- and post-conditioning derived neuroprotection

Mild ischaemic exposures before or after severe injurious ischaemia that elicit neuroprotective responses are referred to as preconditioning and post-conditioning. The corresponding molecular mechanisms of neuroprotection are not completely understood. Identification of the genes and associated pathways of corresponding neuroprotection would provide insight into neuronal survival, potential therapeutic approaches and assessments of therapies for stroke. The objectives of this study were to use global gene expression approach to infer the molecular mechanisms in pre- and post-conditioning-derived neuroprotection in cortical neurons following oxygen and glucose deprivation (OGD) in vitro and then to apply these findings to predict corresponding functional pathways. To this end, microarray analysis was applied to rat cortical neurons with or without the pre- and post-conditioning treatments at 3-h post-reperfusion, and differentially expressed transcripts were subjected to statistical, hierarchical clustering and pathway analyses. The expression patterns of 3,431 genes altered under all conditions of ischaemia (with and without pre- or post-conditioning). We identified 1,595 genes that were commonly regulated within both the pre- and post-conditioning treatments. Cluster analysis revealed that transcription profiles clustered tightly within controls, non-conditioned OGD and neuroprotected groups. Two clusters defining neuroprotective conditions associated with up- and downregulated genes were evident. The five most upregulated genes within the neuroprotective clusters were Tagln, Nes, Ptrf, Vim and Adamts9, and the five most downregulated genes were Slc7a3, Bex1, Brunol4, Nrxn3 and Cpne4. Pathway analysis revealed that the intracellular and second messenger signalling pathways in addition to cell death were predominantly associated with downregulated pre- and post-conditioning associated genes, suggesting that modulation of cell death and signal transduction pathways plays a role in the neuroprotection. A high degree of similarity in the pathways associated with the differentially expressed genes in the pre- and post-conditioning treatments suggests that similar molecular mechanisms may mediate their neuroprotective effects.

Maternal inhalation of surface-coated nanosized titanium dioxide (UV-Titan) in C57BL/6 mice: effects in prenatally exposed offspring on hepatic DNA damage and gene expression

We investigated effects of maternal pulmonary exposure to titanium dioxide (UV-Titan) on prenatally exposed offspring. Time-mated mice (C57BL/6BomTac) were inhalation exposed (1 h/day to 42 mg UV-Titan/m(3) aerosolised powder or filtered air) during gestation days (GDs) 8-18. We evaluated DNA strand breaks using the comet assay in bronchoalveolar lavage (BAL) cells and livers of the time-mated mice (5 and 26-27 days after inhalation exposure), and in livers of the offspring (post-natal days (PND) 2 and 22). We also analysed hepatic gene expression in newborns using DNA microarrays. UV-Titan exposure did not induce DNA strand breaks in time-mated mice or their offspring. Transcriptional profiling of newborn livers revealed changes in the gene expression related to the retinoic acid signalling pathway in the females, while gene expression in male offspring was unaffected. Changes may be a secondary response to maternal inflammation although no direct link was evident through gene expression analysis.

Gene expression profiling of the fathead minnow (Pimephales promelas) neuroendocrine brain in response to pulp and paper mill effluents

The toxicity of pulp and paper mill effluents (PPMEs) has been greatly decreased, yet some continue to negatively affect fish reproduction. We hypothesized that PPMEs are affecting the brain resulting in decreased reproductive performance. Our goal was to use gene expression profiling to test the hypothesis that PPMEs are having an effect on neural systems in the fathead minnow (FHM; Pimephales promelas) in vivo. Sexually mature male and female FHM were exposed to 100% final biotreated PPMEs from 5 different sources for 5 days. Using an oligo-array (15K genes) we examined the effect of PPMEs on gene expression in the hypothalamus of female fish. We validated selected genes (cholecystokinin, RevErbbeta2, and urotensin I) that were identified by microarray analysis using real-time reverse-transcriptase polymerase chain reaction (RT-PCR). We compared the FHM microarray dataset to multiple microarray datasets from experiments conducted with goldfish injected with different dopaminergic pharmaceuticals to examine whether PPMEs could be affecting the dopamine system. Exposure of FHM to PPMEs resulted in varying degrees of spawning inhibition. Microarray analysis revealed surprisingly few genes in the brain that were commonly affected by the different PPMEs. Real-time RT-PCR confirmed the changes in expression for cholecystokinin, RevErbbeta2, and urotensin I. Comparison of the FHM and goldfish microarray datasets suggest that some PPMEs may be acting on the dopamine system. We show that PPMEs are neuroactive in fish and may be acting through some of the pathways in a manner similar to dopamine.

Thyroid hormone may regulate mRNA abundance in liver by acting on microRNAs

MicroRNAs (miRNAs) are extensively involved in diverse biological processes. However, very little is known about the role of miRNAs in mediating the action of thyroid hormones (TH). Appropriate TH levels are known to be critically important for development, differentiation and maintenance of metabolic balance in mammals. We induced transient hypothyroidism in juvenile mice by short-term exposure to methimazole and perchlorate from post natal day (PND) 12 to 15. The expression of miRNAs in the liver was analyzed using Taqman Low Density Arrays (containing up to 600 rodent miRNAs). We found the expression of 40 miRNAs was significantly altered in the livers of hypothyroid mice compared to euthyroid controls. Among the miRNAs, miRs-1, 206, 133a and 133b exhibited a massive increase in expression (50- to 500-fold). The regulation of TH on the expression of miRs-1, 206, 133a and 133b was confirmed in various mouse models including: chronic hypothyroid, short-term hyperthyroid and short-term hypothyroid followed by TH supplementation. TH regulation of these miRNAs was also confirmed in mouse hepatocyte AML 12 cells. The expression of precursors of miRs-1, 206, 133a and 133b were examined in AML 12 cells and shown to decrease after TH treatment, only pre-mir-206 and pre-mir-133b reached statistical significance. To identify the targets of these miRNAs, DNA microarrays were used to examine hepatic mRNA levels in the short-term hypothyroid mouse model relative to controls. We found transcripts from 92 known genes were significantly altered in these hypothyroid mice. Web-based target predication software (TargetScan and Microcosm) identified 14 of these transcripts as targets of miRs-1, 206, 133a and 133b. The vast majority of these mRNA targets were significantly down-regulated in hypothyroid mice, corresponding with the up-regulation of miRs-1, 206, 133a and 133b in hypothyroid mouse liver. To further investigate target genes, miR-206 was over-expressed in AML 12 cells. TH treatment of cells over-expressing miR-206 resulted in decreased miR-206 expression, and a significant increase in two predicted target genes, Mup1 and Gpd2. The results suggest that TH regulation of these genes may occur secondarily via miR-206. These studies provide new insight into the role of miRNAs in mediating TH regulation of gene expression.

Testing for mean and correlation changes in microarray experiments: an application for pathway analysis

Background

Microarray experiments examine the change in transcript levels of tens of thousands of genes simultaneously. To derive meaningful data, biologists investigate the response of genes within specific pathways. Pathways are comprised of genes that interact to carry out a particular biological function. Existing methods for analyzing pathways focus on detecting changes in the mean or over-representation of the number of differentially expressed genes relative to the total of genes within the pathway. The issue of how to incorporate the influence of correlation among the genes is not generally addressed.

Results

In this paper, we propose a non-parametric rank test for analyzing pathways that takes into account the correlation among the genes and compared two existing methods, Global and Gene Set Enrichment Analysis (GSEA), using two publicly available data sets. A simulation study was conducted to demonstrate the advantage of the rank test method.

Conclusions

The data indicate the advantages of the rank test. The method can distinguish significant changes in pathways due to either correlations or changes in the mean or both. From the simulation study the rank test out performed Global and GSEA. The greatest gain in performance was for the sample size case which makes the application of the rank test ideal for microarray experiments.

Global transcriptional characterization of a mouse pulmonary epithelial cell line for use in genetic toxicology

Prior to its application for in vitro toxicological assays, thorough characterization of a cell line is essential. The present study uses global transcriptional profiling to characterize a lung epithelial cell line (FE1) derived from MutaMouse [White, P.A., Douglas, G.R., Gingerich, J., Parfett, C., Shwed, P., Seligy, V., Soper, L., Berndt, L., Bayley, J., Wagner, S., Pound, K., Blakey, D., 2003. Development and characterization of a stable epithelial cell line from Muta Mouse lung. Environmental and Molecular Mutagenesis 42, 166-184]. Results presented here demonstrate the origin of the FE1 lung cell line as epithelial, presenting both type I and type II alveolar phenotype. An assessment of toxicologically-relevant genes, including those involved in the response to stress and stimuli, DNA repair, cellular metabolism, and programmed cell death, revealed changes in expression of 22-27% of genes in one or more culture type (proliferating and static FE1 cultures, primary epithelial cultures) compared with whole lung isolates. Gene expression analysis at 4 and 24h following benzo(a)pyrene exposure revealed the induction of cyp1a1, cyp1a2, and cyp1b1 in FE1 cells and lung isolates. The use of DNA microarrays for gene expression profiling allows an improved understanding of global, coordinated cellular events arising in cells under different physiological conditions. Taken together, these data indicate that the FE1 cell line is derived from a cell type relevant to toxic responses in vivo, and shows some similarity in response to chemical insult as the original tissue.

Lack of acute phase response in the livers of mice exposed to diesel exhaust particles or carbon black by inhalation

BACKGROUND

Epidemiologic and animal studies have shown that particulate air pollution is associated with increased risk of lung and cardiovascular diseases. Although the exact mechanisms by which particles induce cardiovascular diseases are not known, studies suggest involvement of systemic acute phase responses, including C-reactive protein (CRP) and serum amyloid A (SAA) in humans. In this study we test the hypothesis that diesel exhaust particles (DEP) - or carbon black (CB)-induced lung inflammation initiates an acute phase response in the liver.

RESULTS

Mice were exposed to filtered air, 20 mg/m3 DEP or CB by inhalation for 90 minutes/day for four consecutive days; we have previously shown that these mice exhibit pulmonary inflammation (Saber AT, Bornholdt J, Dybdahl M, Sharma AK, Loft S, Vogel U, Wallin H. Tumor necrosis factor is not required for particle-induced genotoxicity and pulmonary inflammation., Arch. Toxicol. 79 (2005) 177-182). As a positive control for the induction of an acute phase response, mice were exposed to 12.5 mg/kg of lipopolysaccharide (LPS) intraperitoneally. Quantitative real time RT-PCR was used to examine the hepatic mRNA expression of acute phase proteins, serum amyloid P (Sap) (the murine homologue of Crp) and Saa1 and Saa3. While significant increases in the hepatic expression of Sap, Saa1 and Saa3 were observed in response to LPS, their levels did not change in response to DEP or CB. In a comprehensive search for markers of an acute phase response, we analyzed liver tissue from these mice using high density DNA microarrays. Globally, 28 genes were found to be significantly differentially expressed in response to DEP or CB. The mRNA expression of three of the genes (serine (or cysteine) proteinase inhibitor, clade A, member 3C, apolipoprotein E and transmembrane emp24 domain containing 3) responded to both exposures. However, these changes were very subtle and were not confirmed by real time RT-PCR.

CONCLUSION

Our findings collectively suggest that Sap, Saa1 and Saa3 are not induced in livers of mice exposed to DEP or CB. Despite pulmonary inflammation in these mice, global transcriptional profiling of liver did not reveal any hepatic response following exposure by inhalation.

Germ-line mutations, DNA damage, and global hypermethylation in mice exposed to particulate air pollution in an urban/industrial location

Particulate air pollution is widespread, yet we have little understanding of the long-term health implications associated with exposure. We investigated DNA damage, mutation, and methylation in gametes of male mice exposed to particulate air pollution in an industrial/urban environment. C57BL/CBA mice were exposed in situ to ambient air near two integrated steel mills and a major highway, alongside control mice breathing high-efficiency air particulate (HEPA) filtered ambient air. PCR analysis of an expanded simple tandem repeat (ESTR) locus revealed a 1.6-fold increase in sperm mutation frequency in mice exposed to ambient air for 10 wks, followed by a 6-wk break, compared with HEPA-filtered air, indicating that mutations were induced in spermatogonial stem cells. DNA collected after 3 or 10 wks of exposure did not exhibit increased mutation frequency. Bulky DNA adducts were below the detection threshold in testes samples, suggesting that DNA reactive chemicals do not reach the germ line and cause ESTR mutation. In contrast, DNA strand breaks were elevated at 3 and 10 wks, possibly resulting from oxidative stress arising from exposure to particles and associated airborne pollutants. Sperm DNA was hypermethylated in mice breathing ambient relative to HEPA-filtered air and this change persisted following removal from the environmental exposure. Increased germ-line DNA mutation frequencies may cause population-level changes in genetic composition and disease. Changes in methylation can have widespread repercussions for chromatin structure, gene expression and genome stability. Potential health effects warrant extensive further investigation.

Consensus recommendations to promote and advance predictive systems toxicology and toxicogenomics

The number of high throughput -omics technologies continues to grow. Toxicogenomic application of these technologies is poised to greatly influence current regulatory toxicology. However, many changes are needed before a systems biology level approach can be effectively incorporated into the regulatory toxicology framework. A workshop was held at the Annual Environmental Mutagen Society meeting in Vancouver, British Columbia, on advances in -omics applications. A number of recommendations emerged from the workshop discussion (beyond what activities are currently ongoing) aimed at advancing the ultimate goal of predictive systems toxicology from the present formative state of toxicogenomics. Recommendations include: (1) encouraging investigators to embrace open-access data sharing, (2) increasing current database and curation capacity, (3) establishment of large collaborative projects investigating multiple -omics endpoints within the same groups of animals, (4) mechanisms to encourage collaborative science including increasing the value of junior authorship on multi-authored papers and changes in the promotion process, (5) further development of standardized protocols, and (6) investment from the funding agencies and toxicology community to build the required infrastructure.

A single-molecule PCR approach to the measurement of induced expanded simple tandem repeat instability in vitro

Sensitive and precise models are needed to identify potential genotoxicity at environmentally relevant doses of mutagens. The size length alterations in expanded simple tandem repeat (ESTR) loci have been used as a biomarker of genetic instability caused by a variety of agents in the mouse germline. The mechanisms operating in both spontaneous and induced instability are poorly understood. We have developed a single-molecule polymerase chain reaction (SM-PCR) method to investigate mutation at the mouse ESTR locus Ms6-hm in the murine C3H/10T1/2 embryonic cell line. Growth of cells to levels of high cell density induced increased ESTR instability, with mutation frequencies 5.1-fold (+/-2.8) over sub-confluent cultures. Accordingly, cell cultures were maintained at sub-confluent levels for further investigations of the induction of ESTR mutation by genotoxic agents. Treatment with the DNA alkylating agent N-nitroso-N-ethylurea (ENU) resulted in a 1.94-fold (+/-1.1) increase in mutation frequency, similar to responses measured previously in the germline in vivo. Therefore, mutagen exposure can also affect somatic (non-meiotic) rapidly dividing mouse cells. This SM-PCR approach eliminates the requirement of sub-cloning individual treated cells, thereby, reducing the time needed to screen for ESTR mutation, and will be a very useful tool for future investigations into the mechanisms involved in ESTR mutation.

Automation of cDNA microarray hybridization and washing yields improved data quality

Microarray technology allows the analysis of whole-genome transcription within a single hybridization, and has become a standard research tool. It is extremely important to minimize variation in order to obtain high quality microarray data that can be compared among experiments and laboratories. The majority of facilities implement manual hybridization approaches for microarray studies. We developed an automated method for cDNA microarray hybridization that uses equivalent pre-hybridization, hybridization and washing conditions to the suggested manual protocol. The automated method significantly decreased variability across microarray slides compared to manual hybridization. Although normalized signal intensities for buffer-only spots across the chips were identical, significantly reduced variation and inter-quartile ranges were obtained using the automated workstation. This decreased variation led to improved correlation among technical replicates across slides in both the Cy3 and Cy5 channels.

Molecular insight into the effects of hypothyroidism on the developing cerebellum

Despite the recognized importance of thyroid hormones for normal brain development, little is known about the critical molecular events underlying this role. We investigated the molecular basis of thyroid hormone action on the developing brain by comparing genome-wide gene expression patterns in the cerebellum between euthyroid and hypothyroid juvenile mice using microarrays. Pregnant dams were treated with 0.1% or 0.04% 6-propyl-2-thiouracil (PTU) in drinking water continuously from day 13 post conception until weaning to produce hypothyroid pups. Cerebella were collected from vehicle and 0.1% PTU treated pups at post-natal day (PND) 15, and mRNA from these was subjected to microarray analysis using Agilent high-density oligonucleotide chips. Statistical analysis (MAANOVA) revealed significant differential expression in 2940 genes including 1357 up- and 1583 down-regulated genes. Further analysis (combined MAANOVA and ANOVA) identified 204 significantly altered genes. Hypothyroidism had a greater effect on gene expression in male than in female pups. Transcriptional changes in several genes [Syt12 (Synaptotagmin 12), Rcor (RE1-silencing transcription factor co-repressor), Bag3 (Bcl-associated athanogene 3), p21, cyclin D, Bax (Bcl2-associated X protein), and Pcp2 (Purkinje cell protein 2)] were confirmed using real-time (RT) PCR analysis. Significantly altered expression of Bag3 in cerebella from PND 15 and PND 60 pups exposed to PTU suggests permanent functional alterations in the hypothyroid brain. The thyroid hormone negative regulation of Rcor expression was confirmed in vitro using HepG2 cells. In addition to Rcor, expression of several other genes that code for critical components of the REST (RE1-silencing transcription factor) pathway was shown to be altered in hypothyroid animals. These results suggest that modification of this pathway may have a significant role in causing impaired development in the hypothyroid brain.

Monitoring for induced heritable mutations in natural populations: application of minisatellite DNA screening

The need to understand the role that anthropogenic chemicals play in generating germline mutations is critical, both from an ecological and a human health perspective. Exposure to complex mixtures of urban and industrial chemicals is widespread and we have little understanding of the long-term implications to populations and gene pools. It has recently been suggested that minisatellite DNA mutations may be sensitive biomarkers for induced heritable mutations in populations exposed to radioactive and non-radioactive contamination in their environments. Minisatellite loci are attractive targets for mutational analyses because they undergo a rate of mutation much greater than unique sequence DNA and with DNA fingerprinting many loci can be scanned simultaneously. As a result, the technique is statistically powerful requiring relatively small sample sizes (compared to other in situ mutation assays) and is reasonably cost and time efficient. This paper will review the application of minisatellite mutation screening to the field of genetic toxicology.