Search

Toxicology has been an active research field for many decades, with academic, industrial and government involvement. Modern omics and computational approaches are changing the field, from merely disease-specific observational models into target-specific predictive models. Traditionally, toxicology has strong links with other fields such as biology, chemistry, pharmacology, and medicine. With the rise of synthetic and new engineered materials, alongside ongoing prioritisation needs in chemical risk assessment for existing chemicals, early predictive evaluations are becoming of utmost importance to both scientific and regulatory purposes. ELIXIR is an intergovernmental organisation that brings together life science resources from across Europe. To coordinate the linkage of various life science efforts around modern predictive toxicology, the establishment of a new ELIXIR Community is seen as instrumental. In the past few years, joint efforts, building on incidental overlap, have been piloted in the context of ELIXIR. For example, the EU-ToxRisk, diXa, HeCaToS, transQST, and the nanotoxicology community have worked with the ELIXIR TeSS, Bioschemas, and Compute Platforms and activities. In 2018, a core group of interested parties wrote a proposal, outlining a sketch of what this new ELIXIR Toxicology Community would look like. A recent workshop (held September 30th to October 1st, 2020) extended this into an ELIXIR Toxicology roadmap and a shortlist of limited investment-high gain collaborations to give body to this new community. This Whitepaper outlines the results of these efforts and defines our vision of the ELIXIR Toxicology Community and how it complements other ELIXIR activities.

Abstract The data currently described was generated within the EU/FP7 HeCaToS project (Hepatic and Cardiac Toxicity Systems modeling). The project aimed to develop an in silico prediction system to contribute to drug safety assessment for humans. For this purpose, multi-omics data of repeated dose toxicity were obtained for 10 hepatotoxic and 10 cardiotoxic compounds. Most data were gained from in vitro experiments in which 3D microtissues (either hepatic or cardiac) were exposed to a therapeutic (physiologically relevant concentrations calculated through PBPK-modeling) or a toxic dosing profile (IC20 after 7 days). Exposures lasted for 14 days and samples were obtained at 7 time points (therapeutic doses: 2-8-24-72-168-240-336 h; toxic doses 0-2-8-24-72-168-240 h). Transcriptomics (RNA sequencing & microRNA sequencing), proteomics (LC-MS), epigenomics (MeDIP sequencing) and metabolomics (LC-MS & NMR) data were obtained from these samples. Furthermore, functional endpoints (ATP content, Caspase3/7 and O2 consumption) were measured in exposed microtissues. Additionally, multi-omics data from human biopsies from patients are available. This data is now being released to the scientific community through the BioStudies data repository ( https://www.ebi.ac.uk/biostudies/ ).

Aims: Anthracyclines (ANTs) are essential chemotherapeutic agents; however, their adverse effects can lead to heart failure in cancer survivors. While long non-coding RNAs (lncRNAs) have become new players in cellular processes, there is limited knowledge on lncRNA expression related to anthracyclines-induced cardiotoxicity. This study investigates the lncRNA profiles in human cardiac microtissues exposed to 3 popular ANTs, namely doxorubicin, epirubicin, and idarubicin, as well as in heart biopsies from ANT-treated patients. Methods and results: The in vitro microtissues were exposed to each ANT at 2 doses over 2 weeks; the transcriptome data was collected at 7 time points. The human biopsies were collected from heart failure patients who underwent ANT treatment and control subjects. Over 100 lncRNAs were differentially expressed in each in vitro ANT treatment condition compared to control samples; 16 of them were differentially expressed across all ANT-treated conditions. The lncRNA databases and literature revealed insight on how these lncRNAs relate to heart failure and cellular functions. For instance, H19 and RMRP are involved in heart failure progression, while BDNF-AS is a cardiomyocyte damage-associated gene; SNHG7 is a cardiac hypertrophy regulator. PCAT19 can promote the miR‐182/PDK4 axis and modulate p53 expression, whereas SNHG29 can regulate the Wnt/β-catenin signaling pathway via the miR-223–3p/CTNND1 axis. Other lncRNAs, which were only differentially expressed in particular ANT-treated conditions, are also involved in cardiomyocyte damage and heart failure disease. The alterations of these lncRNA expressions in the in vitro cardiac tissue were also affirmed by similar changes in the human biopsies. Conclusion: This study revealed several lncRNAs that can be potential biomarkers or targets for further ANT-induced cardiotoxicity investigation, according to the transcriptome in both human cardiac microtissues expose to ANTs as well as in heart biopies form ANT-treated patients. Especially, H19 lncRNA showed its contribution to on-target toxicity, in which it is involved in both chemoresistance and cardiotoxic mechanism.

Abstract Background Hepatocellular carcinoma (HCC) is one of the leading causes of cancer death in the world owing to limitations in its prognosis. The current prognosis approaches include radiological examination and detection of serum biomarkers, however, both have limited efficiency and are ineffective in early prognosis. Due to such limitations, we propose to use RNA-Seq data for evaluating putative higher accuracy biomarkers at the transcript level that could help in early prognosis. Methods To identify such potential transcript biomarkers, RNA-Seq data for healthy liver and various HCC cell models were subjected to five different machine learning algorithms: random forest, K-nearest neighbor, Naïve Bayes, support vector machine, and neural networks. Various metrics, namely sensitivity, specificity, MCC, informedness, and AUC-ROC (except for support vector machine) were evaluated. The algorithms that produced the highest values for all metrics were chosen to extract the top features that were subjected to recursive feature elimination. Through recursive feature elimination, the least number of features were obtained to differentiate between the healthy and HCC cell models. Results From the metrics used, it is demonstrated that the efficiency of the known protein biomarkers for HCC is comparatively lower than complete transcriptomics data. Among the different machine learning algorithms, random forest and support vector machine demonstrated the best performance. Using recursive feature elimination on top features of random forest and support vector machine three transcripts were selected that had an accuracy of 0.97 and kappa of 0.93. Of the three transcripts, two were protein coding (PARP2–202 and SPON2–203) and one was a non-coding transcript (CYREN-211). Lastly, we demonstrated that these three selected transcripts outperformed randomly taken three transcripts (15,000 combinations), hence were not chance findings, and could then be an interesting candidate for new HCC biomarker development. Conclusion Using RNA-Seq data combined with machine learning approaches can aid in finding novel transcript biomarkers. The three biomarkers identified: PARP2–202, SPON2–203, and CYREN-211, presented the highest accuracy among all other transcripts in differentiating the healthy and HCC cell models. The machine learning pipeline developed in this study can be used for any RNA-Seq dataset to find novel transcript biomarkers. Code: www.github.com/rajinder4489/ML_biomarkers

Anthracycline-induced cardiotoxicity is the most severe collateral effect of chemotherapy originated by an excess of oxidative stress in cardiomyocytes that leads to cardiac dysfunction. We assessed clinical data from patients with breast cancer receiving anthracyclines and searched for discriminating microRNAs between patients that developed cardiotoxicity (cases) and those that did not (controls), using RNA sequencing and regression analysis. Serum levels of 25 microRNAs were differentially expressed in cases versus controls within the first year after anthracycline treatment, as assessed by three different regression models (elastic net, Robinson and Smyth exact negative binomial test and random forest). MiR-4732-3p was the only microRNA identified in all regression models and was downregulated in patients that experienced cardiotoxicity. MiR-4732-3p was also present in neonatal rat cardiomyocytes and cardiac fibroblasts and was modulated by anthracycline treatment. A miR-4732-3p mimic was cardioprotective in cardiac and fibroblast cultures, following doxorubicin challenge, in terms of cell viability and ROS levels. Notably, administration of the miR-4732-3p mimic in doxorubicin-treated rats preserved cardiac function, normalized weight loss, induced angiogenesis, and decreased apoptosis, interstitial fibrosis and cardiac myofibroblasts. At the molecular level, miR-4732-3p regulated genes of TGFβ and Hippo signaling pathways. Overall, the results indicate that miR-4732-3p is a novel biomarker of cardiotoxicity that has therapeutic potential against anthracycline-induced heart damage.

Anthracyclines, including doxorubicin, idarubicin, and epirubicin, are common antitumor drugs as well as well-known cardiotoxic agents. This study analyzed the proteomics alteration in cardiac tissues caused by these 3 anthracyclines analogs. The in vitro human cardiac microtissues were exposed to drugs in 2 weeks; the proteomic data were measured at 7 time points. The heart biopsy data were collected from heart failure patients, in which some patients underwent anthracycline treatment. The anthracyclines-affected proteins were separately identified in the in vitro and in vivo dataset using the WGCNA method. These proteins engage in different cellular pathways including translation, metabolism, mitochondrial function, muscle contraction, and signaling pathways. From proteins detected in 2 datasets, a protein-protein network was established with 4 hub proteins, and 7 weighted proteins from both cardiac microtissue and human biopsies data. These 11 proteins, which involve in mitochondrial functions and the NF-κB signaling pathway, could provide insights into the anthracycline toxic mechanism. Some of them, such as HSPA5, BAG3, and SH3BGRL, are cardiac therapy targets or cardiotoxicity biomarkers. Other proteins, such as ATP5F1B and EEF1D, showed similar responses in both the in vitro and in vivo data. This suggests that the in vitro outcomes could link to clinical phenomena in proteomic analysis.

Background: Epirubicin (EPI) is an important anticancer drug that is well-known for its cardiotoxic side effect. Studying epigenetic modification such as DNA methylation can help to understand the EPI-related toxic mechanisms in cardiac tissue. In this study, we analyzed the DNA methylation profile in a relevant human cell model and inspected the expression of differentially methylated genes at the transcriptome level to understand how changes in DNA methylation could affect gene expression in relation to EPI-induced cardiotoxicity. Methods: Human cardiac microtissues were exposed to either therapeutic or toxic (IC20) EPI doses during 2 weeks. The DNA and RNA were collected from microtissues in triplicates at 2, 8, 24, 72, 168, 240, and 336 hours of exposure. Methylated DNA immunoprecipitation-sequencing (MeDIP-seq) analysis was used to detect DNA methylation levels in EPI-treated and control samples. The MeDIP-seq data were analyzed and processed using the QSEA package with a recently published workflow. RNA sequencing (RNA-seq) was used to measure global gene expression in the same samples. Results: After processing the MeDIP-seq data, we detected 35, 37, 15 candidate genes which show strong methylated alterations between all EPI-treated, EPI therapeutic and EPI toxic dose-treated samples compared to control, respectively. For several genes, gene expressions changed compatibly reflecting the DNA methylation regulation. Conclusions: The observed DNA methylation modifications provide further insights into the EPI-induced cardiotoxicity. Multiple differentially methylated genes under EPI treatment, such as SMARCA4, PKN1, RGS12, DPP9, NCOR2, SDHA, POLR2A, and AGPAT3, have been implicated in different cardiac dysfunction mechanisms. Together with other differentially methylated genes, these genes can be candidates for further investigations of EPI-related toxic mechanisms. Data Repository: The data has been generated by the HeCaToS project (http://www.ebi.ac.uk/biostudies) under accession numbers S-HECA433 and S-HECA434 for the MeDIP-seq data and S-HECA11 for the RNA-seq data. The R code is available on Github (https://github.com/NhanNguyen000/MeDIP).

Posttraumatic stress disorder (PTSD) is a psychiatric disorder that can develop upon exposure to a traumatic event. While most people are able to recover promptly, others are at increased risk of developing PTSD. However, the exact underlying biological mechanisms of differential susceptibility are unknown. Identifying biomarkers of PTSD could assist in its diagnosis and facilitate treatment planning. Here, we identified serum microRNAs (miRNAs) of subjects that underwent a traumatic event and aimed to assess their potential to serve as diagnostic biomarkers of PTSD. Next-generation sequencing was performed to examine circulating miRNA profiles of 24 members belonging to the Dutch military cohort Prospective Research in Stress-Related Military Operations (PRISMO). Three groups were selected: “susceptible” subjects who developed PTSD after combat exposure, “resilient” subjects without PTSD, and nonexposed control subjects (N = 8 per group). Differential expression analysis revealed 22 differentially expressed miRNAs in PTSD subjects compared to controls and 1 in PTSD subjects compared to resilient individuals (after multiple testing correction and a log2 fold-change cutoff of ≥|1|). Weighted Gene Coexpression Network Analysis (WGCNA) identified a module of coexpressed miRNAs which could distinguish between the three groups. In addition, receiver operating characteristic curve analyses suggest that the miRNAs with the highest module memberships could have a strong diagnostic accuracy as reflected by high areas under the curves. Overall, the results of our pilot study suggest that serum miRNAs could potentially serve as diagnostic biomarkers of PTSD, both individually or grouped within a cluster of coexpressed miRNAs. Larger studies are now needed to validate and build upon these preliminary findings.

Manganese (Mn) is an important element; yet acute and/or chronic exposure to this metal has been linked to neurotoxicity and neurodegenerative illnesses such as Parkinson’s disease and others via an unknown mechanism. To better understand it, we exposed a human neuroblastoma cell model (SH-SY5Y) to two Mn chemical species, MnCl2 and Citrate of Mn(II) (0–2000 µM), followed by a cell viability assay, transcriptomics, and bioinformatics. Even though these cells have been chemically and genetically modified, which may limit the significance of our findings, we discovered that by using RA-differentiated cells instead of undifferentiated SH-SY5Y cell line, both chemical species induce a similar toxicity, potentially governed by disruption of protein metabolism, with some differences. The MnCl2 altered amino acid metabolism, which affects RNA metabolism and protein synthesis. Citrate of Mn(II), however, inhibited the E3 ubiquitin ligases–target protein degradation pathway, which can lead to the buildup of damaged/unfolded proteins, consistent with histone modification. Finally, we discovered that Mn(II)-induced cytotoxicity in RA-SH-SY5Y cells shared 84 percent of the pathways involved in neurodegenerative diseases.

Long-term exposure to air pollution has been associated with several adverse health effects including cardiovascular, respiratory diseases and cancers. However, underlying molecular alterations remain to be further investigated. The aim of this study is to investigate the effects of long-term exposure to air pollutants on (a) average DNA methylation at functional regions and, (b) individual differentially methylated CpG sites. An assumption is that omic measurements, including the methylome, are more sensitive to low doses than hard health outcomes.This study included blood-derived DNA methylation (Illumina-HM450 methylation) for 454 Italian and 159 Dutch participants from the European Prospective Investigation into Cancer and Nutrition (EPIC). Long-term air pollution exposure levels, including NO2, NOx, PM2.5, PMcoarse, PM10, PM2.5 absorbance (soot) were estimated using models developed within the ESCAPE project, and back-extrapolated to the time of sampling when possible. We meta-analysed the associations between the air pollutants and global DNA methylation, methylation in functional regions and epigenome-wide methylation. CpG sites found differentially methylated with air pollution were further investigated for functional interpretation in an independent population (EnviroGenoMarkers project), where (N=613) participants had both methylation and gene expression data available.Exposure to NO2 was associated with a significant global somatic hypomethylation (p-value=0.014). Hypomethylation of CpG island's shores and shelves and gene bodies was significantly associated with higher exposures to NO2 and NOx. Meta-analysing the epigenome-wide findings of the 2 cohorts did not show genome-wide significant associations at single CpG site level. However, several significant CpG were found if the analyses were separated by countries. By regressing gene expression levels against methylation levels of the exposure-related CpG sites, we identified several significant CpG-transcript pairs and highlighted 5 enriched pathways for NO2 and 9 for NOx mainly related to the immune system and its regulation.Our findings support results on global hypomethylation associated with air pollution, and suggest that the shores and shelves of CpG islands and gene bodies are mostly affected by higher exposure to NO2 and NOx. Functional differences in the immune system were suggested by transcriptome analyses. Keywords: Air pollution, Epigenome-wide DNA methylation, Illumina 450k human methylation array, Particulate matter, NOx, EPIC

Abstract Background Due to their lack of repair capacity mitochondria are critical targets for environmental toxicants. We studied genes and pathways reflecting mitochondrial responses to short- and medium-term PM10 exposure. Methods Whole genome gene expression was measured in peripheral blood of 98 adults (49% women). We performed linear regression analyses stratified by sex and adjusted for individual and temporal characteristics to investigate alterations in gene expression induced by short-term (week before blood sampling) and medium-term (month before blood sampling) PM10 exposure. Overrepresentation analyses (ConsensusPathDB) were performed to identify enriched mitochondrial associated pathways and gene ontology sets. Thirteen Human MitoCarta genes were measured by means of quantitative real-time polymerase chain reaction (qPCR) along with mitochondrial DNA (mtDNA) content in an independent validation cohort (n = 169, 55.6% women). Results Overrepresentation analyses revealed significant pathways (p-value

Abstract Background Air pollution exposure during pregnancy has been associated with adverse birth outcomes and health problems later in life. We investigated sex-specific transcriptomic responses to gestational long- and short-term exposure to particulate matter with a diameter

The exposome encompasses an individual's exposure to exogenous chemicals, as well as endogenous chemicals that are produced or altered in response to external stressors. While the exposome concept has been established for human health, its principles can be extended to include broader ecological issues. The assessment of exposure is tightly interlinked with hazard assessment. Here, we explore if mechanistic understanding of the causal links between exposure and adverse effects on human health and the environment can be improved by integrating the exposome approach with the adverse outcome pathway (AOP) concept that structures and organizes the sequence of biological events from an initial molecular interaction of a chemical with a biological target to an adverse outcome. Complementing exposome research with the AOP concept may facilitate a mechanistic understanding of stress-induced adverse effects, examine the relative contributions from various components of the exposome, determine the primary risk drivers in complex mixtures, and promote an integrative assessment of chemical risks for both human and environmental health. Keywords: Exposome, AOP, Systems toxicology, Systems biology, Systems chemistry, Risk assessment

In toxicogenomics, functional annotation is an important step to gain additional insights into genes with aberrant expression that drive pathophysiological mechanisms. Nevertheless, there exists a gap on annotation of these genes which often hampers the interpretation of results and limits their applicability in translational medicine. In this study, we evaluated the coverage of functional annotations of differentially expressed genes (DEGs) induced by 10 selected compounds from the TG-GATEs database identified as high- or no-risk in causing drug-induced liver injury (most-DILI or no-DILI, respectively) using in vitro human data. Functional roles of DEGs not present in the most common biological annotation databases – termed “dark genes” – were unveiled via literature mining and via the identification of shared regulatory transcription factors or signaling pathways. Our results demonstrated that there were approximately 13% of dark genes induced by these compounds in vitro and we were able to obtain additional relevant information for up to 76% of those. Using interactome data from several sources, we have uncovered genes such as LRBA, and WDR26 as highly connected in the protein network that play roles in drug response. Genes such as MALAT1, H19, and MIR29C – whose links to hepatotoxicity have been confirmed – were identified as markers for the most-DILI group and appeared as top hits across all literature-based mining methods. Furthermore, we investigated the potential impact of dark genes on liver toxicity by identifying their rat orthologs in combination with their correlation to drug-induced liver pathologies observed in vivo following chemical exposure. We identified a set of important regulatory transcription factors of dark genes for all most-DILI compounds including E2F1 and JUND with supporting evidences in literature and we found Magee1 correlated with chemically induced bile duct hyperplasia and adverse responses at 29 days in rats in vivo. In conclusion, in this study we show the potential role of these poorly annotated genes in mechanisms underlying hepatotoxicity and offer a number of computational approaches that may help to minimize current gaps in gene annotation and highlight their values as potential biomarkers in toxicological studies.

When studying time courses of biological measurements and comparing these to other measurements eg. gene expression and phenotypic endpoints, the analysis is complicated by the fact that although the associated elements may show the same patterns of behaviour, the changes do not occur simultaneously. In these cases standard correlation-based measures of similarity will fail to find significant associations. Dynamic time warping (DTW) is a technique which can be used in these situations to find the optimal match between two time courses, which may then be assessed for its significance. We implement DTW4Omics, a tool for performing DTW in R. This tool extends existing R scripts for DTW making them applicable for "omics" datasets where thousands entities may need to be compared with a range of markers and endpoints. It includes facilities to estimate the significance of the matches between the supplied data, and provides a set of plots to enable the user to easily visualise the output. We illustrate the utility of this approach using a dataset linking the exposure of the colon carcinoma Caco-2 cell line to oxidative stress by hydrogen peroxide (H2O2) and menadione across 9 timepoints and show that on average 85% of the genes found are not obtained from a standard correlation analysis between the genes and the measured phenotypic endpoints. We then show that when we analyse the genes identified by DTW4Omics as significantly associated with a marker for oxidative DNA damage (8-oxodG), through over-representation, an Oxidative Stress pathway is identified as the most over-represented pathway demonstrating that the genes found by DTW4Omics are biologically relevant. In contrast, when the positively correlated genes were similarly analysed, no pathways were found. The tool is implemented as an R Package and is available, along with a user guide from http://web.tgx.unimaas.nl/svn/public/dtw/.

Toxicology has been an active research field for many decades, with academic, industrial and government involvement. Modern omics and computational approaches are changing the field, from merely disease-specific observational models into target-specific predictive models. Traditionally, toxicology has strong links with other fields such as biology, chemistry, pharmacology and medicine. With the rise of synthetic and new engineered materials, alongside ongoing prioritisation needs in chemical risk assessment for existing chemicals, early predictive evaluations are becoming of utmost importance to both scientific and regulatory purposes. ELIXIR is an intergovernmental organisation that brings together life science resources from across Europe. To coordinate the linkage of various life science efforts around modern predictive toxicology, the establishment of a new ELIXIR Community is seen as instrumental. In the past few years, joint efforts, building on incidental overlap, have been piloted in the context of ELIXIR. For example, the EU-ToxRisk, diXa, HeCaToS, transQST, and the nanotoxicology community have worked with the ELIXIR TeSS, Bioschemas, and Compute Platforms and activities. In 2018, a core group of interested parties wrote a proposal, outlining a sketch of what this new ELIXIR Toxicology Community would look like. A recent workshop (held September 30th to October 1st, 2020) extended this into an ELIXIR Toxicology roadmap and a shortlist of limited investment-high gain collaborations to give body to this new community. This Whitepaper outlines the results of these efforts and defines our vision of the ELIXIR Toxicology Community and how it complements other ELIXIR activities.

Toxicology has been an active research field for many decades, with academic, industrial and government involvement. Modern omics and computational approaches are changing the field, from merely disease-specific observational models into target-specific predictive models. Traditionally, toxicology has strong links with other fields such as biology, chemistry, pharmacology and medicine. With the rise of synthetic and new engineered materials, alongside ongoing prioritisation needs in chemical risk assessment for existing chemicals, early predictive evaluations are becoming of utmost importance to both scientific and regulatory purposes. ELIXIR is an intergovernmental organisation that brings together life science resources from across Europe. To coordinate the linkage of various life science efforts around modern predictive toxicology, the establishment of a new ELIXIR Community is seen as instrumental. In the past few years, joint efforts, building on incidental overlap, have been piloted in the context of ELIXIR. For example, the EU-ToxRisk, diXa, HeCaToS, transQST, and the nanotoxicology community have worked with the ELIXIR TeSS, Bioschemas, and Compute Platforms and activities. In 2018, a core group of interested parties wrote a proposal, outlining a sketch of what this new ELIXIR Toxicology Community would look like. A recent workshop (held September 30th to October 1st, 2020) extended this into an ELIXIR Toxicology roadmap and a shortlist of limited investment-high gain collaborations to give body to this new community. This Whitepaper outlines the results of these efforts and defines our vision of the ELIXIR Toxicology Community and how it complements other ELIXIR activities.