Meta-analysis of transcriptomic datasets using benchmark dose modeling shows value in supporting radiation risk assessment

Signature analysis of high-throughput transcriptomics screening data for mechanistic inference and chemical grouping

High-throughput transcriptomics (HTTr) uses gene expression profiling to characterize the biological activity of chemicals in in vitro cell-based test systems. As an extension of a previous study testing 44 chemicals, HTTr was used to screen an additional 1,751 unique chemicals from the EPA's ToxCast collection in MCF7 cells using 8 concentrations and an exposure duration of 6 h. We hypothesized that concentration-response modeling of signature scores could be used to identify putative molecular targets and cluster chemicals with similar bioactivity. Clustering and enrichment analyses were conducted based on signature catalog annotations and ToxPrint chemotypes to facilitate molecular target prediction and grouping of chemicals with similar bioactivity profiles. Enrichment analysis based on signature catalog annotation identified known mechanisms of action (MeOAs) associated with well-studied chemicals and generated putative MeOAs for other active chemicals. Chemicals with predicted MeOAs included those targeting estrogen receptor (ER), glucocorticoid receptor (GR), retinoic acid receptor (RAR), the NRF2/KEAP/ARE pathway, AP-1 activation, and others. Using reference chemicals for ER modulation, the study demonstrated that HTTr in MCF7 cells was able to stratify chemicals in terms of agonist potency, distinguish ER agonists from antagonists, and cluster chemicals with similar activities as predicted by the ToxCast ER Pathway model. Uniform manifold approximation and projection (UMAP) embedding of signature-level results identified novel ER modulators with no ToxCast ER Pathway model predictions. Finally, UMAP combined with ToxPrint chemotype enrichment was used to explore the biological activity of structurally related chemicals. The study demonstrates that HTTr can be used to inform chemical risk assessment by determining in vitro points of departure, predicting chemicals' MeOA and grouping chemicals with similar bioactivity profiles.

Considerations for application of benchmark dose modeling in radiation research: workshop highlights

BACKGROUND

Exposure to different forms of ionizing radiation occurs in diverse occupational, medical, and environmental settings. Improving the accuracy of the estimated health risks associated with exposure is therefore, essential for protecting the public, particularly as it relates to chronic low dose exposures. A key aspect to understanding health risks is precise and accurate modeling of the dose-response relationship. Toward this vision, benchmark dose (BMD) modeling may be a suitable approach for consideration in the radiation field. BMD modeling is already extensively used for chemical hazard assessments and is considered statistically preferable to identifying low and no observed adverse effects levels. BMD modeling involves fitting mathematical models to dose-response data for a relevant biological endpoint and identifying a point of departure (the BMD, or its lower bound). Recent examples in chemical toxicology show that when applied to molecular endpoints (e.g. genotoxic and transcriptional endpoints), BMDs correlate to points of departure for more apical endpoints such as phenotypic changes (e.g. adverse effects) of interest to regulatory decisions. This use of BMD modeling may be valuable to explore in the radiation field, specifically in combination with adverse outcome pathways, and may facilitate better interpretation of relevant in vivo and in vitro dose-response data. To advance this application, a workshop was organized on June 3rd, 2022, in Ottawa, Ontario that brought together BMD experts in chemical toxicology and the radiation scientific community of researchers, regulators, and policy-makers. The workshop's objective was to introduce radiation scientists to BMD modeling and its practical application using case examples from the chemical toxicity field and demonstrate the BMDExpress software using a radiation dataset. Discussions focused on the BMD approach, the importance of experimental design, regulatory applications, its use in supporting the development of adverse outcome pathways, and specific radiation-relevant examples.

CONCLUSIONS

Although further deliberations are needed to advance the use of BMD modeling in the radiation field, these initial discussions and partnerships highlight some key steps to guide future undertakings related to new experimental work.

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.

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.

COHERE - strengthening cooperation within the Canadian government on radiation research

PURPOSE

Canadian Organization on Health Effects from Radiation Exposure (COHERE) is a government initiative to better understand biological and human health risks from ionizing radiation exposures relevant to occupational and environmental settings (<100 mGy, <6 mGy/h). It is currently a partnership between two federal agencies, Health Canada (HC) and the Canadian Nuclear Safety Commission (CNSC). COHERE's vision is to contribute knowledge to reduce scientific uncertainties from low dose and dose-rate exposures. COHERE will advance our understanding by bridging the knowledge gap between human health risks and linkages to molecular- and cellular-level responses to radiation. Research focuses on identifying sensitive, early, and key molecular events of relevance to risk assessment.

CONCLUSIONS

The initiative will address questions of relevance to better apprize Canadians, including radiation workers and members of the public and Indigenous peoples, on health risks from low dose radiation exposure and inform radiation protection frameworks at a national and international level. Furthermore, it will support global efforts to conduct collaborative undertakings and better coordinate research. Here, we describe a historical overview of the research conducted, the strategic research agenda that outlines the scientific framework, stakeholders, opportunities to harmonize internationally, and how research outcomes will better inform communication of risk to Canadians.