ELIXIR and Toxicology: a community in development

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.

Omics technologies to understand drug toxicity mechanisms

Drug side effects are an important study subject in pharmacology. Recent omics technologies provide a range of omics data and help to understand the biological mechanisms involved in drug effects. These modern technologies provide significant support to all biological disciplines, including drug toxicology. In this review, we provide an overview the use of omics applications to understand drug side effects at the molecular level. We discuss by available omics technologies, their possible uses, as well as their advantages and limitations.

Manganese-Induced Neurotoxicity through Impairment of Cross-Talk Pathways in Human Neuroblastoma Cell Line SH-SY5Y Differentiated with Retinoic Acid

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, MnCl₂ 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 MnCl₂ 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.

A Model-Based Workflow to Benchmark the Clinical Cholestasis Risk of Drugs

We present a generic workflow combining physiology-based computational modeling and in vitro data to assess the clinical cholestatic risk of different drugs systematically. Changes in expression levels of genes involved in the enterohepatic circulation of bile acids were obtained from an in vitro assay mimicking 14 days of repeated drug administration for 10 marketed drugs. These changes in gene expression over time were contextualized in a physiology-based bile acid model of glycochenodeoxycholic acid. The simulated drug-induced response in bile acid concentrations was then scaled with the applied drug doses to calculate the cholestatic potential for each compound. A ranking of the cholestatic potential correlated very well with the clinical cholestasis risk obtained from medical literature. The proposed workflow allows benchmarking the cholestatic risk of novel drug candidates. We expect the application of our workflow to significantly contribute to the stratification of the cholestatic potential of new drugs and to support animal-free testing in future drug development.

Identifying novel transcript biomarkers for hepatocellular carcinoma (HCC) using RNA-Seq datasets and machine learning

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.

Contextualization of causal regulatory networks from toxicogenomics data applied to drug-induced liver injury

In recent years, network-based methods have become an attractive analytical approach for toxicogenomics studies. They can capture not only the global changes of regulatory gene networks but also the relationships between their components. Among them, a causal reasoning approach depicts the mechanisms of regulation that connect upstream regulators in signaling networks to their downstream gene targets. In this work, we applied CARNIVAL, a causal network contextualisation tool, to infer upstream signaling networks deregulated in drug-induced liver injury (DILI) from gene expression microarray data from the TG-GATEs database. We focussed on six compounds that induce observable histopathologies linked to DILI from repeated dosing experiments in rats. We compared responses in vitro and in vivo to identify potential cross-platform concordances in rats as well as network preservations between rat and human. Our results showed similarities of enriched pathways and network motifs between compounds. These pathways and motifs induced the same pathology in rats but not in humans. In particular, the causal interactions "LCK activates SOCS3, which in turn inhibits TFDP1" was commonly identified as a regulatory path among the fibrosis-inducing compounds. This potential pathology-inducing regulation illustrates the value of our approach to generate hypotheses that can be further validated experimentally.

Translational Proteomics Analysis of Anthracycline-Induced Cardiotoxicity From Cardiac Microtissues to Human Heart Biopsies

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.

FuSe: a tool to move RNA-Seq analyses from chromosomal/gene loci to functional grouping of mRNA transcripts

Summary

Typical RNA sequencing (RNA-Seq) analyses are performed either at the gene level by summing all reads from the same locus, assuming that all transcripts from a gene make a protein or at the transcript level, assuming that each transcript displays unique function. However, these assumptions are flawed, as a gene can code for different types of transcripts and different transcripts are capable of synthesizing similar, different or no protein. As a consequence, functional changes are not well illustrated by either gene or transcript analyses. We propose to improve RNA-Seq analyses by grouping the transcripts based on their similar functions. We developed FuSe to predict functional similarities using the primary and secondary structure of proteins. To estimate the likelihood of proteins with similar functions, FuSe computes two confidence scores: knowledge (KS) and discovery (DS) for protein pairs. Overlapping protein pairs exhibiting high confidence are grouped to form ‘similar function protein groups’ and expression is calculated for each functional group. The impact of using FuSe is demonstrated on in vitro cells exposed to paracetamol, which highlight genes responsible for cell adhesion and glycogen regulation which were earlier shown to be not differentially expressed with traditional analysis methods.

Availability and implementation

The source code is available at https://github.com/rajinder4489/FuSe. Data for APAP exposure are available in the BioStudies database (http://www.ebi.ac.uk/biostudies) under accession numbers S-HECA143, S-HECA(158) and S-HECA139.

Supplementary information

Supplementary data are available at Bioinformatics online.

Altered sphingolipid function in Alzheimer's disease; a gene regulatory network approach

Sphingolipids (SLs) are bioactive lipids involved in various important physiological functions. The SL pathway has been shown to be affected in several brain-related disorders, including Alzheimer's disease (AD). Recent evidence suggests that epigenetic dysregulation plays an important role in the pathogenesis of AD as well. Here, we use an integrative approach to better understand the relationship between epigenetic and transcriptomic processes in regulating SL function in the middle temporal gyrus of AD patients. Transcriptomic analysis of 252 SL-related genes, selected based on GO term annotations, from 46 AD patients and 32 healthy age-matched controls, revealed 103 differentially expressed SL-related genes in AD patients. Additionally, methylomic analysis of the same subjects revealed parallel hydroxymethylation changes in PTGIS, GBA, and ITGB2 in AD. Subsequent gene regulatory network-based analysis identified 3 candidate genes, that is, SELPLG, SPHK1 and CAV1 whose alteration holds the potential to revert the gene expression program from a diseased towards a healthy state. Together, this epigenomic and transcriptomic approach highlights the importance of SL-related genes in AD, and may provide novel biomarkers and therapeutic alternatives to traditionally investigated biological pathways in AD.

Comparing in vitro human liver models to in vivo human liver using RNA-Seq

The liver plays an important role in xenobiotic metabolism and represents a primary target for toxic substances. Many different in vitro cell models have been developed in the past decades. In this study, we used RNA-sequencing (RNA-Seq) to analyze the following human in vitro liver cell models in comparison to human liver tissue: cancer-derived cell lines (HepG2, HepaRG 3D), induced pluripotent stem cell-derived hepatocyte-like cells (iPSC-HLCs), cancerous human liver-derived assays (hPCLiS, human precision cut liver slices), non-cancerous human liver-derived assays (PHH, primary human hepatocytes) and 3D liver microtissues. First, using CellNet, we analyzed whether these liver in vitro cell models were indeed classified as liver, based on their baseline expression profile and gene regulatory networks (GRN). More comprehensive analyses using non-differentially expressed genes (non-DEGs) and differential transcript usage (DTU) were applied to assess the coverage for important liver pathways. Through different analyses, we noticed that 3D liver microtissues exhibited a high similarity with in vivo liver, in terms of CellNet (C/T score: 0.98), non-DEGs (10,363) and pathway coverage (highest for 19 out of 20 liver specific pathways shown) at the beginning of the incubation period (0 h) followed by a decrease during long-term incubation for 168 and 336 h. PHH also showed a high degree of similarity with human liver tissue and allowed stable conditions for a short-term cultivation period of 24 h. Using the same metrics, HepG2 cells illustrated the lowest similarity (C/T: 0.51, non-DEGs: 5623, and pathways coverage: least for 7 out of 20) with human liver tissue. The HepG2 are widely used in hepatotoxicity studies, however, due to their lower similarity, they should be used with caution. HepaRG models, iPSC-HLCs, and hPCLiS ranged clearly behind microtissues and PHH but showed higher similarity to human liver tissue than HepG2 cells. In conclusion, this study offers a resource of RNA-Seq data of several biological replicates of human liver cell models in vitro compared to human liver tissue.

Valproic acid promotes mitochondrial dysfunction in primary human hepatocytes in vitro; impact of C/EBPα-controlled gene expression

Valproic acid (VPA) is a frequently prescribed anti-epileptic drug which is known to cause liver toxicity and steatosis through mitochondrial dysfunction. Nevertheless the mechanisms underlying these adverse effects are incompletely understood. In this study, we determined the effect of relatively short (3 h) or prolonged (72 h) exposure to VPA on mitochondrial function in primary human hepatocytes (PHHs). While 3 h VPA exposure did not affect oxygen consumption rates (OCRs) in PHHs, prolonged exposure (24–72 h) significantly reduced basal and maximal OCRs. Given that in particular prolonged VPA exposure is required to cause mitochondrial dysfunction, we investigated gene expression data after VPA exposure for 24, 48, 72 h and 72 h VPA followed by a 72 h washout period. We were able to reduce the comprehensive gene expression changes into a more comprehensible set of 18 TFs that were predicted to be persistently activated after 72 h of VPA exposure. Lentiviral knock-down of one of the candidate TFs, C/EBPα, partly rescued VPA-induced mitochondrial dysfunction. Furthermore, RNA-Seq analysis of shC/EBPα and shGFP control PHHs identified 24 genuine C/EBPα target genes that are regulated in response to prolonged VPA exposure in PHHs. Altogether this provides new insights on the involvement of C/EBPα in driving VPA-induced mitochondrial dysfunction in human liver cells. This hub gene, with its downstream regulators involved in this deregulation, thus represent potential new biomarkers for VPA-induced mitochondrial dysfunction.

MicroRNA regulation of persistent stress-enhanced memory

Disruption of persistent, stress-associated memories is relevant for treating posttraumatic stress disorder (PTSD) and related syndromes, which develop in a subset of individuals following a traumatic event. We previously developed a stress-enhanced fear learning (SEFL) paradigm in inbred mice that produces PTSD-like characteristics in a subset of mice, including persistently enhanced memory and heightened cFos in the basolateral amygdala complex (BLC) with retrieval of the remote (30-day-old) stress memory. Here, the contribution of BLC microRNAs (miRNAs) to stress-enhanced memory was investigated because of the molecular complexity they achieve through their ability to regulate multiple targets simultaneously. We performed small-RNA sequencing (smRNA-Seq) and quantitative proteomics on BLC tissue collected from mice 1 month after SEFL and identified persistently changed microRNAs, including mir-135b-5p, and proteins associated with PTSD-like heightened fear expression. Viral-mediated overexpression of mir-135b-5p in the BLC of stress-resilient animals enhanced remote fear memory expression and promoted spontaneous renewal 14 days after extinction. Conversely, inhibition of BLC mir-135b-5p in stress-susceptible animals had the opposite effect, promoting a resilient-like phenotype. mir-135b-5p is highly conserved across mammals and was detected in post mortem human amygdala, as well as human serum samples. The mir-135b passenger strand, mir-135b-3p, was significantly elevated in serum from PTSD military veterans, relative to combat-exposed control subjects. Thus, miR-135b-5p may be an important therapeutic target for dampening persistent, stress-enhanced memory and its passenger strand a potential biomarker for responsivity to a mir-135-based therapeutic.

The MicroRNA-based Liquid Biopsy Improves Early Assessment of Lethal Acetaminophen Poisoning: A Case Report

BACKGROUND Acetaminophen overdose is the most common cause of acute liver failure. Nevertheless, new biomarker approaches enabling early prediction of the outcome of the acetaminophen overdose are needed. Recently, using next-generation sequencing analysis of serum from human study participants we uncovered injury-specific signatures of circulating microRNAs (miRNAs) that represented underlying molecular mechanisms of toxicity. This case study is first to show the application of miRNA profiling to assess prognosis of acetaminophen poisoning. CASE REPORT The patient was admitted to the hospital following supra therapeutic acetaminophen ingestion. The patient showed elevated levels of biomarkers of hepatocellular injury alanine aminotransferase, aspartate transaminase, and glutamate dehydrogenase. Even though treatment with N-acetyl cysteine was initiated 24 hours post-ingestion, levels of alanine-aminotransferase and aspartate transaminase peaked at about 40 hours post ingestion of acetaminophen. We analyzed global circulating miRNA levels from 24 consecutive serum samples from this study participant covering the period from admission to time of death. CONCLUSIONS The resulting global miRNA profiles were compared with profiles from study participants with non-lethal acetaminophen poisoning and healthy controls. At the admission, the miRNA profiles of both lethal and non-lethal acetaminophen poisoning showed induction of cellular stress and oxidative damage. Later, the miRNA profiles of the lethal poisoning featured fibrosis and coagulation pathways while profiles of non-lethal cases resembled those of healthy study participants. Although additional confirmatory studies are needed, our case study is first to indicate that global miRNA profiles to be used as liquid biopsies have potential to facilitate the assessment of acetaminophen poisoning.

Circulating Serum MicroRNAs as Potential Diagnostic Biomarkers of Posttraumatic Stress Disorder: A Pilot Study

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.

Balance of Active, Passive, and Anatomical Cardiac Properties in Doxorubicin-Induced Heart Failure

Late-onset heart failure (HF) is a known side effect of doxorubicin chemotherapy. Typically, patients are diagnosed when already at an irreversible stage of HF, which allows few or no treatment options. Identifying the causes of compromised cardiac function in this patient group may improve early patient diagnosis and support treatment selection. To link doxorubicin-induced changes in cardiac cellular and tissue mechanical properties to overall cardiac function, we apply a multiscale biophysical biomechanics model of the heart to measure the plausibility of changes in model parameters representing the passive, active, or anatomical properties of the left ventricle for reproducing measured patient phenotypes. We create representative models of healthy controls (N = 10) and patients with HF induced by (N = 22) or unrelated to (N = 25) doxorubicin therapy. The model predicts that HF in the absence of doxorubicin is characterized by a 2- to 3-fold stiffness increase, decreased tension (0–20%), and ventricular dilation (of order 10–30%). HF due to doxorubicin was similar but showed stronger bias toward reduced active contraction (10–30%) and less dilation (0–20%). We find that changes in active, passive, and anatomical properties all play a role in doxorubicin-induced cardiotoxicity phenotypes. Differences in parameter changes between patient groups are consistent with doxorubicin cardiotoxicity having a greater dependence on reduced cellular contraction and less anatomical remodeling than HF not caused by doxorubicin.

Genes associated with Parkinson's disease respond to increasing polychlorinated biphenyl levels in the blood of healthy females

Polychlorinated biphenyls (PCBs) are a class of widespread environmental pollutants, commonly found in human blood, that have been suggested to be linked to the occurrence of sporadic Parkinson's disease (PD). It has been reported that some non-coplanar PCBs accumulate in the brains of female PD patients. To improve our understanding of the association between PCB exposure and PD risk we have applied whole transcriptome gene expression analysis in blood cells from 594 PCB-exposed subjects (369 female, 225 male). Interestingly, we observe that in females, blood levels of non-coplanar PCBs appear to be associated with expression levels of PD-specific genes. However, no such association was detected in males. Among the 131 PD-specific genes affected, 39 have been shown to display similar changes in expression levels in the substantia nigra of deceased PD patients. Especially among the down-regulated genes, transcripts of genes involved in neurotransmitter vesicle-related functions were predominant.

Embracing the Dark Side: Computational Approaches to Unveil the Functionality of Genes Lacking Biological Annotation in Drug-Induced Liver Injury

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.

Moving forward in carcinogenicity assessment: Report of an EURL ECVAM/ESTIV workshop

There is an increased need to develop novel alternative approaches to the two-year rodent bioassay for the carcinogenicity assessment of substances where the rodent bioassay is still a basic requirement, as well as for those substances where animal use is banned or limited or where information gaps are identified within legislation. The current progress in this area was addressed in a EURL ECVAM- ESTIV workshop held in October 2016, in Juan les Pins. A number of initiatives were presented and discussed, including data-driven, technology-driven and pathway-driven approaches. Despite a seemingly diverse range of strategic developments, commonalities are emerging. For example, providing insight into carcinogenicity mechanisms is becoming an increasingly appreciated aspect of hazard assessment and is suggested to be the best strategy to drive new developments. Thus, now more than ever, there is a need to combine and focus efforts towards the integration of available information between sectors. Such cross-sectorial harmonisation will aid in building confidence in new approach methods leading to increased implementation and thus a decreased necessity for the two-year rodent bioassay.

DNA methylation and exposure to ambient air pollution in two prospective cohorts

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.

Transcriptome-wide analyses indicate mitochondrial responses to particulate air pollution exposure

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 PM₁₀ 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) PM₁₀ 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 <0.05) related to mitochondrial genome maintenance and apoptosis for short-term exposure and to the electron transport chain (ETC) for medium-term exposure in women. For men, medium-term PM₁₀ exposure was associated with the Tri Carbonic Acid cycle. In an independent study population, we validated several ETC genes, including UQCRH and COX7C (q-value <0.05), and some genes crucial for the maintenance of the mitochondrial genome, including LONP1 (q-value: 0.07) and POLG (q-value: 0.04) in women.

CONCLUSIONS

In this exploratory study, we identified mitochondrial genes and pathways associated with particulate air pollution indicating upregulation of energy producing pathways as a potential mechanism to compensate for PM-induced mitochondrial damage.

Newborn sex-specific transcriptome signatures and gestational exposure to fine particles: findings from the ENVIRONAGE birth cohort

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 < 2.5 μm (PM_2.5) in order to elucidate potential underlying mechanisms of action.

METHODS

Whole genome gene expression was investigated in cord blood of 142 mother-newborn pairs that were enrolled in the ENVIRONAGE birth cohort. Daily PM_2.5 exposure levels were calculated for each mother's home address using a spatial-temporal interpolation model in combination with a dispersion model to estimate both long- (annual average before delivery) and short- (last month of pregnancy) term exposure. We explored the association between gene expression levels and PM_2.5 exposure, and identified modulated pathways by overrepresentation analysis and gene set enrichment analysis.

RESULTS

Some processes were altered in both sexes for long- (e.g. DNA damage) or short-term exposure (e.g. olfactory signaling). For long-term exposure in boys neurodevelopment and RhoA pathways were modulated, while in girls defensin expression was down-regulated. For short-term exposure we identified pathways related to synaptic transmission and mitochondrial function (boys) and immune response (girls).

CONCLUSIONS

This is the first whole genome gene expression study in cord blood to identify sex-specific pathways altered by PM_2.5. The identified transcriptome pathways could provide new molecular insights as to the interaction pattern of early life PM_2.5 exposure with the biological development of the fetus.

Sex-Specific Associations between Particulate Matter Exposure and Gene Expression in Independent Discovery and Validation Cohorts of Middle-Aged Men and Women

BACKGROUND

Particulate matter (PM) exposure leads to premature death, mainly due to respiratory and cardiovascular diseases.

OBJECTIVES

Identification of transcriptomic biomarkers of air pollution exposure and effect in a healthy adult population.

METHODS

Microarray analyses were performed in 98 healthy volunteers (48 men, 50 women). The expression of eight sex-specific candidate biomarker genes (significantly associated with PM₁₀ in the discovery cohort and with a reported link to air pollution-related disease) was measured with qPCR in an independent validation cohort (75 men, 94 women). Pathway analysis was performed using Gene Set Enrichment Analysis. Average daily PM_2.5 and PM₁₀ exposures over 2-years were estimated for each participant's residential address using spatiotemporal interpolation in combination with a dispersion model.

RESULTS

Average long-term PM₁₀ was 25.9 (± 5.4) and 23.7 (± 2.3) μg/m³ in the discovery and validation cohorts, respectively. In discovery analysis, associations between PM₁₀ and the expression of individual genes differed by sex. In the validation cohort, long-term PM₁₀ was associated with the expression of DNAJB5 and EAPP in men and ARHGAP4 (p = 0.053) in women. AKAP6 and LIMK1 were significantly associated with PM₁₀ in women, although associations differed in direction between the discovery and validation cohorts. Expression of the eight candidate genes in the discovery cohort differentiated between validation cohort participants with high versus low PM₁₀ exposure (area under the receiver operating curve = 0.92; 95% CI: 0.85, 1.00; p = 0.0002 in men, 0.86; 95% CI: 0.76, 0.96; p = 0.004 in women).

CONCLUSIONS

Expression of the sex-specific candidate genes identified in the discovery population predicted PM₁₀ exposure in an independent cohort of adults from the same area. Confirmation in other populations may further support this as a new approach for exposure assessment, and may contribute to the discovery of molecular mechanisms for PM-induced health effects.

The exposome in practice: Design of the EXPOsOMICS project

EXPOsOMICS is a European Union funded project that aims to develop a novel approach to the assessment of exposure to high priority environmental pollutants, by characterizing the external and the internal components of the exposome. It focuses on air and water contaminants during critical periods of life. To this end, the project centres on 1) exposure assessment at the personal and population levels within existing European short and long-term population studies, exploiting available tools and methods which have been developed for personal exposure monitoring (PEM); and 2) multiple "omic" technologies for the analysis of biological samples (internal markers of external exposures). The search for the relationships between external exposures and global profiles of molecular features in the same individuals constitutes a novel advancement towards the development of "next generation exposure assessment" for environmental chemicals and their mixtures. The linkage with disease risks opens the way to what are defined here as 'exposome-wide association studies' (EWAS).

From the exposome to mechanistic understanding of chemical-induced adverse effects

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.

Biological marks of early-life socioeconomic experience is detected in the adult inflammatory transcriptome

Consistent evidence is accumulating to link lower socioeconomic position (SEP) and poorer health, and the inflammatory system stands out as a potential pathway through which socioeconomic environment is biologically embedded. Using bloodderived genome-wide transcriptional profiles from 268 Italian participants of the European Prospective Investigation into Cancer and Nutrition (EPIC) cohort, we evaluated the association between early life, young and later adulthood SEP and the expression of 845 genes involved in human inflammatory responses. These were examined individually and jointly using several inflammatory scores. Our results consistently show that participants whose father had a manual (as compared to nonmanual) occupation exhibit, later in life, a higher inflammatory score, hence indicating an overall increased level of expression for the selected inflammatory-related genes. Adopting a life course approach, these associations remained statistically significant upon adjustment for later-in-life socioeconomic experiences. Sensitivity analyses indicated that our findings were not affected by the way the inflammatory score was calculated, and were replicated in an independent study. Our study provides additional evidence that childhood SEP is associated with a sustainable upregulation of the inflammatory transcriptome, independently of subsequent socioeconomic experiences. Our results support the hypothesis that early social inequalities impacts adult physiology.

Persistent transcriptional responses show the involvement of feed-forward control in a repeated dose toxicity study

Chemical carcinogenesis, albeit complex, often relies on modulation of transcription through activation or repression of key transcription factors. While analyzing extensive networks may hinder the biological interpretation, one may focus on dynamic network motifs, among which persistent feed-forward loops (FFLs) are known to chronically influence transcriptional programming. Here, to investigate the relevance a FFL-oriented approach in depth, we have focused on aflatoxin B1-induced transcriptomic alterations during distinct states of exposure (daily administration during 5days followed by a non-exposed period) of human hepatocytes, by exploring known interactions in human transcription. Several TF-coding genes were persistently deregulated after washout of AFB1. Oncogene MYC was identified as the prominent regulator and driver of many FFLs, among which a FFL comprising MYC/HIF1A was the most recurrent. The MYC/HIF1A FFL was also identified and validated in an independent set as the master regulator of metabolic alterations linked to initiation and progression of carcinogenesis, i.e. the Warburg effect, possibly as result of persistent intracellular alterations arising from AFB1 exposure (nuclear and mitochondrial DNA damage, oxidative stress, transcriptional activation by secondary messengers). In summary, our analysis shows the involvement of FFLs as modulators of gene expression suggestive of a carcinogenic potential even after termination of exposure.

Vitamin D insufficient levels during pregnancy and micronuclei frequency in peripheral blood T lymphocytes mothers and newborns (Rhea cohort, Crete)

BACKGROUND & AIMS

Vitamin D deficiency is common among pregnant women and may be associated with several adverse health outcomes including cancer. Micronuclei frequency is a biomarker of early genetic effects and has been used to examine the association between genotoxic exposures and cancer. We examined maternal vitamin D levels during pregnancy in associations with micronuclei frequency in maternal blood and in cord blood.

METHODS

173 mothers and 171 newborns born between 2007 and 2008 in Heraklion (Crete, Greece) were included in the study. Between 14th and 18th weeks of gestation we collected information on maternal diet using food frequency questionnaires (FFQs). We measured maternal serum concentrations of 25-hydroxyvitamin D [25(OH)D] between the first and second trimester of pregnancy. We estimated dietary vitamin D intake using information from FFQ. After delivery we collected cord blood and maternal peripheral blood. We used the cytokinesis-block micronucleus (CBMN) assay to assess the frequencies of micronucleated cells in binucleated T lymphocytes (MNBN).

RESULTS

Maternal insufficient serum levels of 25(OH)D (<50 nmol/L) during pregnancy were associated with increased MNBN frequency in cord blood [IRR = 1.32 (95%CI: 1.00, 1.72)]. This increase was higher for newborns with birth weight above the third quartile [≥3.500 kg; IRR = 2.21 (1.26, 3.89)]. Similarly, low levels of dietary vitamin D were associated with increased MNBN frequency in cord blood [middle tertile IRR = 1.08 (0.78, 1.47), lower tertile IRR = 1.51 (1.06, 2.14)]. Insufficient levels of vitamin D were not associated with MNBN in mothers.

CONCLUSION

Our results suggest that vitamin D deficiency during pregnancy increases genotoxic risks in newborns. The prevalence of vitamin D deficiency globally is high and it is important to further investigate whether vitamin D supplementation or similar interventions during pregnancy could prevent DNA damage at early stages of life.

New insights into BaP-induced toxicity: role of major metabolites in transcriptomics and contribution to hepatocarcinogenesis

Benzo(a)pyrene (BaP) is a ubiquitous carcinogen resulting from incomplete combustion of organic compounds and also present at high levels in cigarette smoke. A wide range of biological effects has been attributed to BaP and its genotoxic metabolite BPDE, but the contribution to BaP toxicity of intermediary metabolites generated along the detoxification path remains unknown. Here, we report for the first time how 3-OH-BaP, 9,10-diol and BPDE, three major BaP metabolites, temporally relate to BaP-induced transcriptomic alterations in HepG2 cells. Since BaP is also known to induce AhR activation, we additionally evaluated TCDD to source the expression of non-genotoxic AhR-mediated patterns. 9,10-Diol was shown to activate several transcription factor networks related to BaP metabolism (AhR), oxidative stress (Nrf2) and cell proliferation (HIF-1α, AP-1) in particular at early time points, while BPDE influenced expression of genes involved in cell energetics, DNA repair and apoptotic pathways. Also, in order to grasp the role of BaP and its metabolites in chemical hepatocarcinogenesis, we compared expression patterns from BaP(-metabolites) and TCDD to a signature set of approximately nine thousand gene expressions derived from hepatocellular carcinoma (HCC) patients. While transcriptome modulation by TCDD appeared not significantly related to HCC, BaP and BPDE were shown to deregulate metastatic markers via non-genotoxic and genotoxic mechanisms and activate inflammatory pathways (NF-κβ signaling, cytokine-cytokine receptor interaction). BaP also showed strong repression of genes involved in cholesterol and fatty acid biosynthesis. Altogether, this study provides new insights into BaP-induced toxicity and sheds new light onto its mechanism of action as a hepatocarcinogen.

Inter-laboratory study of human in vitro toxicogenomics-based tests as alternative methods for evaluating chemical carcinogenicity: a bioinformatics perspective

The assessment of the carcinogenic potential of chemicals with alternative, human-based in vitro systems has become a major goal of toxicogenomics. The central read-out of these assays is the transcriptome, and while many studies exist that explored the gene expression responses of such systems, reports on robustness and reproducibility, when testing them independently in different laboratories, are still uncommon. Furthermore, there is limited knowledge about variability induced by the data analysis protocols. We have conducted an inter-laboratory study for testing chemical carcinogenicity evaluating two human in vitro assays: hepatoma-derived cells and hTERT-immortalized renal proximal tubule epithelial cells, representing liver and kidney as major target organs. Cellular systems were initially challenged with thirty compounds, genome-wide gene expression was measured with microarrays, and hazard classifiers were built from this training set. Subsequently, each system was independently established in three different laboratories, and gene expression measurements were conducted using anonymized compounds. Data analysis was performed independently by two separate groups applying different protocols for the assessment of inter-laboratory reproducibility and for the prediction of carcinogenic hazard. As a result, both workflows came to very similar conclusions with respect to (1) identification of experimental outliers, (2) overall assessment of robustness and inter-laboratory reproducibility and (3) re-classification of the unknown compounds to the respective toxicity classes. In summary, the developed bioinformatics workflows deliver accurate measures for inter-laboratory comparison studies, and the study can be used as guidance for validation of future carcinogenicity assays in order to implement testing of human in vitro alternatives to animal testing.

Outdoor air pollution exposures and micronuclei frequencies in lymphocytes from pregnant women and newborns in Crete, Greece (Rhea cohort)

BACKGROUND

Micronuclei (MN) are biomarkers of early genetic effects that have been used to investigate the association between environmental exposures and cancer. However, few studies have examined the association between environmental exposures during pregnancy and MN in mothers and newborns.

OBJECTIVES

We examined MN frequency in maternal blood and in cord blood, in relation to maternal air pollution exposure, and the potential interaction with maternal vitamin C intake and maternal smoking.

METHODS

We used the cytokinesis-block micronucleus assay to assess MN frequency per 1000 bi-nucleated T-lymphocytes from 181 mothers and 183 newborns born in 2007-2008 in Heraklion (Crete, Greece). The ESCAPE land-use regression methods were used to estimate annual mean exposure to outdoor air pollution [particulate matter (PM), black carbon, nitrogen dioxide (NO2) and nitrogen oxides (NOx)] at maternal home addresses. Food frequency questionnaires were used to estimate maternal dietary vitamin C intake during pregnancy. Smoking habits were self-reported using questionnaires which were checked by measuring maternal urinary cotinine levels.

RESULTS

Exposure to PM2.5 was associated with increased MN frequencies in pregnant women [rate ratio [RR (95%CI)] per 5 µg/m(3)=1.53 (1.02, 2.29)]. This increase was considerably higher among women who did not fulfill the recommended vitamin C dietary allowances [RR=9.35 (2.77, 31.61); n=20]. Exposure to PM2.5-10, PM10, NO2 and NOx were also associated with a higher incidence of MN frequencies in smoker women (n=56). No associations were found for newborns.

CONCLUSIONS

We found an association between air pollution, particularly PM2.5, and MN frequency in mothers but not in newborns. This association was more pronounced among women with a lower dietary intake of vitamin C during pregnancy and among women who smoked during pregnancy. While results are clear in mothers, the association between maternal carcinogenic exposures during pregnancy and biomarkers of early biologic effect in the newborn remains poorly understood.

Transcriptomic and metabolomic data integration

Many studies now produce parallel data sets from different omics technologies; however, the task of interpreting the acquired data in an integrated fashion is not trivial. This review covers those methods that have been used over the past decade to statistically integrate and interpret metabolomics and transcriptomic data sets. It defines four categories of approaches, correlation-based integration, concatenation-based integration, multivariate-based integration and pathway-based integration, into which all existing statistical methods fit. It also explores the choices in study design for generating samples for analysis by these omics technologies and the impact that these technical decisions have on the subsequent data analysis options.

High-throughput data integration of RNA-miRNA-circRNA reveals novel insights into mechanisms of benzo[a]pyrene-induced carcinogenicity

The chain of events leading from a toxic compound exposure to carcinogenicity is still barely understood. With the emergence of high-throughput sequencing, it is now possible to discover many different biological components simultaneously. Using two different RNA libraries, we sequenced the complete transcriptome of human HepG2 liver cells exposed to benzo[a]pyrene, a potent human carcinogen, across six time points. Data were integrated in order to reveal novel complex chemical–gene interactions. Notably, we hypothesized that the inhibition of MGMT, a DNA damage response enzyme, by the over-expressed miR-181a-1_3p induced by BaP, may lead to liver cancer over time.

Free radical scavenging and formation by multi-walled carbon nanotubes in cell free conditions and in human bronchial epithelial cells

Background

Certain multi-walled carbon nanotubes (MWCNTs) have been shown to elicit asbestos-like toxicological effects. To reduce needs for risk assessment it has been suggested that the physicochemical characteristics or reactivity of nanomaterials could be used to predict their hazard. Fibre-shape and ability to generate reactive oxygen species (ROS) are important indicators of high hazard materials. Asbestos is a known ROS generator, while MWCNTs may either produce or scavenge ROS. However, certain biomolecules, such as albumin – used as dispersants in nanomaterial preparation for toxicological testing in vivo and in vitro - may reduce the surface reactivity of nanomaterials.

Methods

Here, we investigated the effect of bovine serum albumin (BSA) and cell culture medium with and without BEAS 2B cells on radical formation/scavenging by five MWCNTs, Printex 90 carbon black, crocidolite asbestos, and glass wool, using electron spin resonance (ESR) spectroscopy and linked this to cytotoxic effects measured by trypan blue exclusion assay. In addition, the materials were characterized in the exposure medium (e.g. for hydrodynamic size-distribution and sedimentation rate).

Results

The test materials induced highly variable cytotoxic effects which could generally be related to the abundance and characteristics of agglomerates/aggregates and to the rate of sedimentation. All carbon nanomaterials were found to scavenge hydroxyl radicals (^•OH) in at least one of the solutions tested. The effect of BSA was different among the materials. Two types of long, needle-like MWCNTs (average diameter >74 and 64.2 nm, average length 5.7 and 4.0 μm, respectively) induced, in addition to a scavenging effect, a dose-dependent formation of a unique, yet unidentified radical in both absence and presence of cells, which also coincided with cytotoxicity.

Conclusions

Culture medium and BSA affects scavenging/production of ^•OH by MWCNTs, Printex 90 carbon black, asbestos and glass-wool. An unidentified radical is generated by two long, needle-like MWCNTs and these two CNTs were more cytotoxic than the other CNTs tested, suggesting that this radical could be related to the adverse effects of MWCNTs.

Characterization of the peripheral blood transcriptome in a repeated measures design using a panel of healthy individuals

A repeated measures microarray design with 22 healthy, non-smoking volunteers (aging 32±5years) was set up to study transcriptome profiles in whole blood samples. The results indicate that repeatable data can be obtained with high within-subject correlation. Probes that could discriminate between individuals are associated with immune and inflammatory functions. When investigating possible time trends in the microarray data, we have found no differential expression within a sampling period (within-season effect). Differential expression was observed between sampling seasons and the data suggest a weak response of genes related to immune system functioning. Finally, a high number of probes showed significant season-specific expression variability within subjects. Expression variability increased in springtime and there was an association of the probe list with immune system functioning. Our study suggests that the blood transcriptome of healthy individuals is reproducible over a time period of several months.