Search

BackgroundB-cell lymphomas are a diverse group of hematological neoplasms with differential etiology and clinical trajectories. Increased insights in the etiology and the discovery of prediagnostic markers have the potential to improve the clinical course of these neoplasms.MethodsWe investigated in a prospective study global gene expression in peripheral blood mononuclear cells of 263 incident B-cell lymphoma cases, diagnosed between 1 and 17 years after blood sample collection, and 439 controls, nested within two European cohorts.ResultsOur analyses identified only transcriptomic markers for specific lymphoma subtypes; few markers of multiple myeloma (N = 3), and 745 differentially expressed genes in relation to future risk of chronic lymphocytic leukemia (CLL). The strongest of these associations were consistently found in both cohorts and were related to (B-) cell signaling networks and immune system regulation pathways. CLL markers exhibited very high predictive abilities of disease onset even in cases diagnosed more than 10 years after blood collection.ConclusionsThis is the first investigation on blood cell global gene expression and future risk of B-cell lymphomas. We mainly identified genes in relation to future risk of CLL that are involved in biological pathways, which appear to be mechanistically involved in CLL pathogenesis. Many but not all of the top hits we identified have been reported previously in studies based on tumor tissues, therefore suggesting that a mixture of preclinical and early disease markers can be detected several years before CLL clinical diagnosis.

The 32P-post-labelling assay for DNA adduct quantification gives the opportunity to examine endogenous exposure to DNA reactive compounds. Most human biomonitoring studies applied white blood cells (WBC) or cells obtained by broncho-alveolar lavages (BAL) as source of DNA, but still it is not clear what cell type represents the most reliable indicator for exposure to cigarette smoke-associated genotoxins. At first, we examined DNA adduct levels by means of nuclease P1 (NP1) enriched 32P-post-labelling in separated WBC subpopulations after in vitro incubations for 18 h with 10 μM benzo[a]pyrene (B[a]P). DNA adduct levels were highest in monocytes (10.7±2.9 adducts/108 nucleotides, n=8), followed by lymphocytes (5.9±1.7, n=8), and granulocytes (0.5±0.2, n=8). Secondly, aromatic-DNA adduct levels were determined in BAL cells and WBC-subsets from (non-)smoking volunteers. In smoking individuals, adduct levels were in the ranking order: BAL cells (3.7±1.0, n=5) > monocytes (2.0±0.5, n=8)  lymphocytes (1.6±0.4, n=8) >granulocytes (0.8±0.2, n=8) by NP1-enrichment and monocytes (9.0±3.2, n=5)  lymphocytes (8.0±2.1, n=6) > granulocytes (2.1±0.3, n=7) by butanol-enriched 32P-post-labelling. Aromatic-DNA adduct levels were significantly higher in WBC-subsets of smokers as compared with non-smokers except for DNA adducts in granulocytes using butanol enrichment. Thirdly, dose-response relationships were investigated in mononuclear white blood cells (MNC, i.e. monocytes plus lymphocytes) and BAL-cells of a larger group of smoking individuals (n=78). Adduct levels in MNC were related to daily exposure to cigarette tar (r=0.31, P<0.01). Adduct levels in BAL cells seemed to be correlated with pack-years, but after correction for age this relationship was lost. Butanol extraction resulted in 5-6 fold higher DNA adduct levels in MNC, whereas butanol extraction of BAL-DNA of the same individuals yielded only 2-fold higher adduct levels. The two enrichment procedures of 32P-post-labelling were correlated in BAL cells (r=0.86, P<0.001, n=12). We conclude that particularly MNC are good surrogates for the detection of smoking-related DNA adducts. [ABSTRACT FROM PUBLISHER]

Capecitabine is a chemotherapeutic drug that is widely used as a monotherapy option in advanced cancer patients. After administration, it is converted into its active metabolite 5-fluorouracil (5-FU), a cytotoxic compound that may also induce adverse side effects in the gastrointestinal (GI) tract. Although these side effects can interfere with the continuation of the chemotherapy, diagnostic tools to detect early onset and prevention strategies are not available. In this explorative case study, we aim to identify differentially expressed genes (DEGs) that provide insight into the molecular mechanisms of toxicity induced by 5-FU in healthy colon tissue of breast cancer patients receiving capecitabine. Gene expression responses observed in patients were compared with those established in an in vitro model of healthy colon organoids. Colon biopsies from two patients with advanced breast cancer were collected before and after the treatment with capecitabine and used for RNA sequencing to determine transcriptomic responses. Differential expression analysis resulted in 31 affected genes, showing that the most affected pathways were transport of small molecules, cellular responses to stress, folate metabolism, NF-kB signalling pathway and immune system responses. The most biologically relevant genes were haemoglobin subunits encoding genes, involved in several processes; ATP12A, SLC26A3 and AQP8, involved in the transport of ions and water; TRIM31, a regulator of NF-kB signalling pathway; MST1P2 and MST1L, stimulators of macrophages. Comparison of human in vitro and in vivo responses showed that the gene expression of TRIM31 was similarly altered in the colon organoids exposed to 5-FU. Therefore, this gene constitutes a potential biomarker of colon toxicity that might be used in future in vitro drug safety design and screening., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2022 The Authors. Published by Elsevier B.V. All rights reserved.)

Usage of injectable dermal fillers applied for aesthetic purposes has extensively increased over the years. As such, the number of related adverse reactions has increased, including patients showing severe complications such as product migration, topical swelling and inflammatory reactions of the skin. In order to understand the underlying molecular events of these adverse reactions we performed a genome-wide gene expression study on the multi-cell type human Phenion® Full-Thickness Skin Model exposed to five experimental hyaluronic acid (HA) preparations with increasing cross-linking degree, four commercial fillers from Perfectha®, and non-resorbable filler Bio-Alcamid®. In addition, we evaluated whether cross-linking degree or particle size of the HA-based fillers could be associated with the occurrence of adverse effects. In all cases, exposure to different HA fillers resulted in a clearly elevated gene expression of cytokines and chemokines related to acute inflammation as part of the foreign body response. Furthermore, for one experimental filler genes of OXPHOS complexes I-V were significantly down-regulated (adjusted p-value < 0.05), resulting in mitochondrial dysfunction which can be linked to over-expression of pro-inflammatory cytokines TNFα and IL-1β and chemokine CCL2. Our hypothesis that cross-linking degree or particle size of the HA-based fillers is related to the biological responses induced by these fillers could only partially be confirmed for particle size. In conclusion, our innovative approach resulted in gene expression changes from a human 3D skin model exposed to dermal fillers that mechanistically substantiate aforementioned adverse reactions, and thereby adds to the weight of evidence that these fillers may induce inflammatory and fibrotic responses.

Resorbable tissue fillers for aesthetic purposes can induce severe complications including product migration, late swelling, and inflammatory reactions. The relation between product characteristics and adverse effects is not well understood. We hypothesized that the degree of cross-linking hyaluronic acid (HA) fillers was associated with the occurrence of adverse effects. Five experimental HA preparations similar to HA fillers were synthesized with an increasing degree of cross-linking. Furthermore, a series of commercial fillers (Perfectha ® ) was obtained that differ in degradation time based on the size of their particulate HA components. Cytotoxic responses and cytokine production by human THP-1-derived macrophages exposed to extracts of the evaluated resorbable HA fillers were absent to minimal. Gene expression analysis of the HA-exposed macrophages revealed the responses related to cell cycle control and immune reactivity. Our results could not confirm the hypothesis that the level of cross-linking in our experimental HA fillers or the particulate size of commercial HA fillers is related to the induced biological responses. However, the evaluation of cytokine induction and gene expression in macrophages after biomaterial exposure presents promising opportunities for the development of methods to identify cellular processes that may be predictive for biomaterial-induced responses in patients.

Gefitinib is a tyrosine kinase inhibitor (TKI) that selectively inhibits the epidermal growth factor receptor (EGFR), hampering cell growth and proliferation. Due to its action, gefitinib has been used in the treatment of cancers that present abnormally increased expression of EGFR. However, side effects from gefitinib therapy may occur, among which diarrhoea is most common, that can lead to interruption of the planned therapy in the more severe cases. The mechanisms underlying intestinal toxicity induced by gefitinib are not well understood. Therefore, this study aims at providing insight into these mechanisms based on transcriptomic responses induced in vitro. A 3D culture of healthy human colon and small intestine (SI) organoids was exposed to 0.1, 1, 10 and 30 µM of gefitinib, for a maximum of three days. These drug concentrations were selected using physiologically-based pharmacokinetic simulation considering patient dosing regimens. Samples were used for the analysis of viability and caspase 3/7 activation, image-based analysis of structural changes, as well as RNA isolation and sequencing via high-throughput techniques. Differential gene expression analysis showed that gefitinib perturbed signal transduction pathways, apoptosis, cell cycle, FOXO-mediated transcription, p53 signalling pathway, and metabolic pathways. Remarkably, opposite expression patterns of genes associated with metabolism of lipids and cholesterol biosynthesis were observed in colon versus SI organoids in response to gefitinib. These differences in the organoids' responses could be linked to increased activated protein kinase (AMPK) activity in colon, which can influence the sensitivity of the colon to the drug. Therefore, this study sheds light on how gefitinib induces toxicity in intestinal organoids and provides an avenue towards the development of a potential tool for drug screening and development.

Doxorubicin is widely used in the treatment of different cancers, and its side effects can be severe in many tissues, including the intestines. Symptoms such as diarrhoea and abdominal pain caused by intestinal inflammation lead to the interruption of chemotherapy. Nevertheless, the molecular mechanisms associated with doxorubicin intestinal toxicity have been poorly explored. This study aims to investigate such mechanisms by exposing 3D small intestine and colon organoids to doxorubicin and to evaluate transcriptomic responses in relation to viability and apoptosis as physiological endpoints. The in vitro concentrations and dosing regimens of doxorubicin were selected based on physiologically based pharmacokinetic model simulations of treatment regimens recommended for cancer patients. Cytotoxicity and cell morphology were evaluated as well as gene expression and biological pathways affected by doxorubicin. In both types of organoids, cell cycle, the p53 signalling pathway, and oxidative stress were the most affected pathways. However, significant differences between colon and SI organoids were evident, particularly in essential metabolic pathways. Short time-series expression miner was used to further explore temporal changes in gene profiles, which identified distinct tissue responses. Finally, in silico proteomics revealed important proteins involved in doxorubicin metabolism and cellular processes that were in line with the transcriptomic responses, including cell cycle and senescence, transport of molecules, and mitochondria impairment. This study provides new insight into doxorubicin-induced effects on the gene expression levels in the intestines. Currently, we are exploring the potential use of these data in establishing quantitative systems toxicology models for the prediction of drug-induced gastrointestinal toxicity.

Arsenic (As) contamination in groundwater is responsible for numerous adverse health outcomes among millions of people. Epigenetic alterations are among the most widely studied mechanisms of As toxicity. To understand how As exposure alters gene expression through epigenetic modifications, a systematic genome-wide study was designed to address the impact of multiple important single nucleotide polymorphisms (SNPs) related to As exposure on the methylome of drinking water As-exposed rural subjects from Pakistan. Urinary As levels were used to stratify subjects into low, medium and high exposure groups. Genome-wide DNA methylation was investigated using MeDIP in combination with NimbleGen 2.1 M Deluxe Promotor arrays. Transcriptome levels were measured using Agilent 8 × 60 K expression arrays. Genotyping of selected SNPs (As3MT, DNMT1a, ERCC2, EGFR and MTHFR) was measured and an integrated genetic risk factor for each respondent was calculated by assigning a specific value to the measured genotypes based on known risk allele numbers. To select a representative model related to As exposure we compared 9 linear mixed models comprising of model 1 (including the genetic risk factor), model 2 (without the genetic risk factor) and models with individual SNPs incorporated into the methylome data. Pathway analysis was performed using ConsensusPathDB. Model 1 comprising the integrated genetic risk factor disclosed biochemical pathways including muscle contraction, cardio-vascular diseases, ATR signaling, GPCR signaling, methionine metabolism and chromatin modification in association with hypo- and hyper-methylated gene targets. A unique pathway (direct P53 effector) was found associated with the individual DNMT1a polymorphism due to hyper-methylation of CSE1L and TRRAP. Most importantly, we provide here the first evidence of As-associated DNA methylation in relation with gene expression of ATR, ATF7IP, TPM3, UBE2J2. We report the first evidence that integrating SNPs data with methylome data generates a more representative epigenome profile and discloses a better insight in disease risks of As-exposed individuals., (Copyright © 2021 The Authors. Published by Elsevier Ltd.. All rights reserved.)

5-Fluorouracil (5-FU) is a widely used chemotherapeutical that induces acute toxicity in the small and large intestine of patients. Symptoms can be severe and lead to the interruption of cancer treatments. However, there is limited understanding of the molecular mechanisms underlying 5-FU-induced intestinal toxicity. In this study, well-established 3D organoid models of human colon and small intestine (SI) were used to characterize 5-FU transcriptomic and metabolomic responses. Clinically relevant 5-FU concentrations for in vitro testing in organoids were established using physiologically based pharmacokinetic simulation of dosing regimens recommended for cancer patients, resulting in exposures to 10, 100 and 1000 µM. After treatment, different measurements were performed: cell viability and apoptosis; image analysis of cell morphological changes; RNA sequencing; and metabolome analysis of supernatant from organoids cultures. Based on analysis of the differentially expressed genes, the most prominent molecular pathways affected by 5-FU included cell cycle, p53 signalling, mitochondrial ATP synthesis and apoptosis. Short time-series expression miner demonstrated tissue-specific mechanisms affected by 5-FU, namely biosynthesis and transport of small molecules, and mRNA translation for colon; cell signalling mediated by Rho GTPases and fork-head box transcription factors for SI. Metabolomic analysis showed that in addition to the effects on TCA cycle and oxidative stress in both organoids, tissue-specific metabolic alterations were also induced by 5-FU. Multi-omics integration identified transcription factor E2F1, a regulator of cell cycle and apoptosis, as the best key node across all samples. These results provide new insights into 5-FU toxicity mechanisms and underline the relevance of human organoid models in the safety assessment in drug development., (© 2021. The Author(s).)

The renal proximal tubule is responsible for re-absorption of the majority of the glomerular filtrate and its proper function is necessary for whole-body homeostasis. Aging, certain diseases and chemical-induced toxicity are factors that contribute to proximal tubule injury and chronic kidney disease progression. To better understand these processes, it would be advantageous to generate renal tissues from human induced pluripotent stem cells (iPSC). Here, we report the differentiation and characterization of iPSC lines into proximal tubular-like cells (PTL). The protocol is a step wise exposure of small molecules and growth factors, including the GSK3 inhibitor (CHIR99021), the retinoic acid receptor activator (TTNPB), FGF9 and EGF, to drive iPSC to PTL via cell stages representing characteristics of early stages of renal development. Genome-wide RNA sequencing showed that PTL clustered within a kidney phenotype. PTL expressed proximal tubular-specific markers, including megalin (LRP2), showed a polarized phenotype, and were responsive to parathyroid hormone. PTL could take up albumin and exhibited ABCB1 transport activity. The phenotype was stable for up to 7 days and was maintained after passaging. This protocol will form the basis of an optimized strategy for molecular investigations using iPSC derived PTL.

Purpose: Identify differentially expressed microRNAs (miRNAs) in aqueous humor (AH) and blood of primary open-angle glaucoma (POAG) patients by using small RNA sequencing. These may provide insight into POAG pathophysiology or serve as diagnostic biomarker., Methods: AH and plasma of nine POAG patients and 10 cataract control patients were small RNA sequenced on Illumina NovaSeq 6000. Identification of gene transcripts targeted by differentially expressed miRNAs was done with miRWalk and MirPath. These targets were used for pathway analysis and Gene Ontology enrichment. Diagnostic potential was evaluated by receiver operating characteristics analysis., Results: We identified 715 miRNAs in plasma and 62 miRNAs in AH. Plasma miRNA profile did not differ between POAG and control. In contrast, in AH, seven miRNAs were differentially expressed. Hsa-miR-30a-3p, hsa-miR-143-3p, hsa-miR-211-5p, and hsa-miR-221-3p were upregulated, whereas hsa-miR-92a-3p, hsa-miR-451a, and hsa-miR-486-5p were downregulated in POAG. Compared to previous studies, hsa-mir-143-3p, hsa-miR-211-5p, and hsa-miR-221-3p were reported previously, strengthening their involvement in POAG whereas hsa-miR-30a-3p, hsa-miR-92a-3p, and hsa-miR-486-5p are implicated in POAG for the first time. Identified gene transcripts were involved in several pathways, some implicated in glaucoma before (e.g., TGF-β and neurotrophin signaling), whereas others are new (e.g., prolactin and apelin signaling). In respect to diagnostics, AH concentration of hsa-mir-143-3p had an area under the curve (AUC) of 0.889. Combined with hsa-miR-221-3p, AUC improved to 0.96., Conclusions: Small RNA sequencing identified seven differentially expressed miRNAs in AH of POAG patients. The differentially expressed miRNAs may be useful as POAG biomarkers or could become targets for new therapeutic strategies.

Cancer is a complex disease where cancer cells express epigenetic and transcriptomic mechanisms to promote tumor initiation, progression, and survival. To extract relevant features from the 2019 Cancer Cell Line Encyclopedia (CCLE), a multi-layer nonnegative matrix factorization approach is used. We used relevant feature genes and DNA promoter regions to construct genomic interaction network to study gene-gene and gene-DNA promoter methylation relationships. Here, we identified a set of gene transcripts and methylated DNA promoter regions for different clusters, including one homogeneous lymphoid neoplasms cluster. In this cluster, we found different methylated transcription factors that affect transcriptional activation of EGFR and downstream interactions. Furthermore, the hippo-signaling pathway might not function properly because of DNA hypermethylation and low gene expression of both LATS2 and YAP1. Finally, we could identify a potential dysregulation of the CD28-CD86-CTLA4 axis. Characterizing the interaction of the epigenome and the transcriptome is vital for our understanding of cancer cell line behavior, not only for deepening insights into cancer-related processes but also for future disease treatment and drug development. Here we have identified potential candidates that characterize cancer cell lines, which give insight into the development and progression of cancers.

An increasing number of findings from epidemiological studies support associations between exposure to air pollution and the onset of several diseases, including pulmonary, cardiovascular and neurodegenerative diseases, and malignancies. However, intermediate, and potentially mediating, biological mechanisms associated with exposure to air pollutants are largely unknown. Previous studies on the human exposome have shown that the expression of certain circulating microRNAs (miRNAs), regulators of gene expression, are altered upon exposure to traffic-related air pollutants. In the present study, we investigated the relationship between particulate matter (PM) smaller than 2.5 μm (PM 2.5 ), PM 2.5 absorbance (as a proxy of black carbon and soot), and ultrafine-particles (UFP, smaller than 0.1 μm), measured in healthy volunteers by 24 h personal monitoring (PEM) sessions and global expression levels of peripheral blood miRNAs. The PEM sessions were conducted in four European countries, namely Switzerland (Basel), United Kingdom (Norwich), Italy (Turin), and The Netherlands (Utrecht). miRNAs expression levels were analysed using microarray technology on blood samples from 143 participants. Seven miRNAs, hsa-miR-24-3p, hsa-miR-4454, hsa-miR-4763-3p, hsa-miR-425-5p, hsa-let-7d-5p, hsa-miR-502-5p, and hsa-miR-505-3p were significantly (FDR corrected) expressed in association with PM 2.5 personal exposure, while no significant association was found between miRNA expression and the other pollutants. The results obtained from this investigation suggest that personal exposure to PM 2.5 is associated with miRNA expression levels, showing the potential for these circulating miRNAs as novel biomarkers for air pollution health risk assessment., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2020 Elsevier Ltd. All rights reserved.)

Circulating microRNAs (cimiRNAs) are a class of non-encoding RNAs found in bodily fluids such as blood, cerebrospinal fluid (CSF) and tears. CimiRNAs have been implicated as promising biomarkers for central nervous system (CNS) disorders because they are actively secreted as messengers and are profoundly involved in fine-tuning of developmental and differentiation processes. Furthermore, they are attractive biomarkers because they are extremely stable, tissue enriched and can be determined in a quantitative manner. This review aims to provide a comprehensive assessment on the current progress regarding the potential value of cimiRNAs as CNS biomarkers. Within this framework five CNS disorders are explored which share a common pathological hallmark namely cognitive impairment. The CNS disorders include Major depression disorder (MDD), Bipolar disorder (BD), Schizophrenia (SZ), Alzheimer's disease (AD) and Parkinson disease (PD). The similarities and differences between altered cimiRNAs in the different disorders are described. The miR-29 family, miR-34a-5p and miR-132-3p are discussed as common dysregulated cimiRNAs found in the CNS disorders. Furthermore, it is shown that the type of bodily fluid used for measuring cimiRNAs is important as inconsistencies in cimiRNAs expression directions are found when comparing CSF, blood cell-free and blood cell-bound samples., (Copyright © 2019 Elsevier Ltd. All rights reserved.)

Background: Groundwater Arsenic (As) contamination is a global public health concern responsible for various health implications and a neglected area of environmental health research in Pakistan. Because of interindividual differences in genetic predisposition, As-related health issues may not be equally distributed among the As-exposed population. However, till date, no studies have been conducted including multiple SNPs involved in As metabolism and disease risk using a linear mixed effect model approach to analyze peripheral blood transcriptomics results., Objectives: In order to detect early responses on the gene expression level and to evaluate the impact of selected SNPs inferring disease risks associated with As exposure, we designed a systematic study to investigate blood transcriptomics profiles of 57 differentially exposed rural subjects living in drinking water As-contaminated settings of Lahore and Kasur districts in Punjab Province in southeast Pakistan. Exposure among the subjects was correlated with individual transcriptome responses applying urinary As profiles as the main biomarker for risk stratification., Methods: We performed whole genome gene expression analysis in blood of subjects using microarrays. Linear effect mixed models were applied for evaluating the combined impact of SNPs hypothetically increasing the risk for As exposure-induced health effects (GSTM1, GSTT1, As3MT, DNMT1, MTHFR, ERCC2 and EGFR)., Results: Our findings confirmed important signaling, growth factor, cancer and other disease related pathways known to be associated with increased As exposure levels. In addition, upon implementing our integrative SNPs-based genetic risk factor, pathways associated with an increased risk of NAFLD and diabetes appeared significantly enhanced by down-regulation of genes NDUFV3, IKBKB, IL6R, ADIPOR1, PPARA, OGT and FOXO1., Conclusion: We report the first comprehensive study applying state-of-the-art bioinformatics approaches to address multiple SNP-based inter-individual variability in adverse molecular responses among subjects exposed to drinking water As contamination in Pakistan thereby providing strong evidence of various gene expression targets associated with development of known As-related diseases., (Copyright © 2019 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Four complementary approaches were used to investigate acetaminophen overdose as a risk factor for Parkinson's disease (PD). Circulating microRNAs (miRNAs) serum profiles from acetaminophen-overdosed patients were compared with patients with terminal PD, revealing four shared miRNAs. Similarities were found among molecular structures of dopamine (DA), acetaminophen, and two known PD inducers indicating affinity for dopaminergic transport. Potential interactions between acetaminophen and the human DA transporter were confirmed by molecular docking modeling and binding free energy calculations. Thus, it is plausible that acetaminophen is taken up by the dopaminergic transport system into the substantia nigra (SN). A ChEMBL query identified proteins that are similarly targeted by DA and acetaminophen. Here, we highlight CYP3A4, present in the SN, a predominant metabolizer of acetaminophen into its toxic metabolite N-acetyl-p-benzoquinone imine and shown to be regulated in PD. Overall, based on our results, we hypothesize that overdosing of acetaminophen is a potential risk factor for parkinsonism., (© 2019 Maastricht University. Clinical and Translational Science published by Wiley Periodicals Inc. on behalf of the American Society of Clinical Pharmacology & Therapeutics.)

Background: The development of high throughput sequencing techniques provides us with the possibilities to obtain large data sets, which capture the effect of dynamic perturbations on cellular processes. However, because of the dynamic nature of these processes, the analysis of the results is challenging. Therefore, there is a great need for bioinformatics tools that address this problem., Results: Here we present DynOVis, a network visualization tool that can capture dynamic dose-over-time effects in biological networks. DynOVis is an integrated work frame of R packages and JavaScript libraries and offers a force-directed graph network style, involving multiple network analysis methods such as degree threshold, but more importantly, it allows for node expression animations as well as a frame-by-frame view of the dynamic exposure. Valuable biological information can be highlighted on the nodes in the network, by the integration of various databases within DynOVis. This information includes pathway-to-gene associations from ConsensusPathDB, disease-to-gene associations from the Comparative Toxicogenomics databases, as well as Entrez gene ID, gene symbol, gene synonyms and gene type from the NCBI database., Conclusions: DynOVis could be a useful tool to analyse biological networks which have a dynamic nature. It can visualize the dynamic perturbations in biological networks and allows the user to investigate the changes over time. The integrated data from various online databases makes it easy to identify the biological relevance of nodes in the network. With DynOVis we offer a service that is easy to use and does not require any bioinformatics skills to visualize a network.

One of the major complications that patients experience during pharmacological treatment is the occurrence of adverse drug reactions (ADRs). The most affected organs are the liver, kidney, heart and the gastrointestinal-immune system. In comparison to the other organs, less progress has been made on human-relevant prediction of drug-induced intestinal toxicity, evidencing current large data gaps. The most widely used drugs that are associated with intestinal damage include chemotherapeutics, such as 5-Fluorouracil or Tyrosine Kinase Inhibitors (TKIs), as well as non-steroidal anti-inflammatory drugs (NSAIDs). Chemotherapeutics are regarded as inducers of acute intestinal toxicity whereas NSAIDs are associated with chronic inflammation of the intestine. In view of the fact that only a few studies have been dedicated to studying cellular and genomic responses in relation to drug-induced intestinal ADRs, little is known about how intestinal toxicity develops after exposure to such drugs or which molecular mechanisms are involved. Therefore, new models and experiments are required to establish transcriptomic responses and alterations of molecular markers induced by different medicines. This review summarizes the available information about transcriptomic responses and biomarkers of toxicity induced by 5-FU, NSAIDS or TKIs in different experimental models. Future investigation should address the challenges in predicting intestinal toxicity induced by drugs and unveil specific gene expression profiles that can be applied in the development of safer drugs., (Copyright © 2019 The Author(s). Published by Elsevier Ltd.. All rights reserved.)

Objectives: To characterize the impact of PCB exposure on DNA methylation in peripheral blood leucocytes and to evaluate the corresponding changes in relation to possible health effects, with a focus on B-cell lymphoma., Methods: We conducted an epigenome-wide association study on 611 adults free of diagnosed disease, living in Italy and Sweden, in whom we also measured plasma concentrations of 6 PCB congeners, DDE and hexachlorobenzene., Results: We identified 650 CpG sites whose methylation correlates strongly (FDR < 0.01) with plasma concentrations of at least one PCB congener. Stronger effects were observed in males and in Sweden. This epigenetic exposure profile shows extensive and highly statistically significant overlaps with published profiles associated with the risk of future B-cell chronic lymphocytic leukemia (CLL) as well as with clinical CLL (38 and 28 CpG sites, respectively). For all these sites, the methylation changes were in the same direction for increasing exposure and for higher disease risk or clinical disease status, suggesting an etiological link between exposure and CLL. Mediation analysis reinforced the suggestion of a causal link between exposure, changes in DNA methylation and disease. Disease connectivity analysis identified multiple additional diseases associated with differentially methylated genes, including melanoma for which an etiological link with PCB exposure is established, as well as developmental and neurological diseases for which there is corresponding epidemiological evidence. Differentially methylated genes include many homeobox genes, suggesting that PCBs target stem cells. Furthermore, numerous polycomb protein target genes were hypermethylated with increasing exposure, an effect known to constitute an early marker of carcinogenesis., Conclusions: This study provides mechanistic evidence in support of a link between exposure to PCBs and the etiology of CLL and underlines the utility of omic profiling in the evaluation of the potential toxicity of environmental chemicals., (Copyright © 2019 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Cellular adaptation is important to cope with various stresses induced by altered environmental conditions. By controlling mRNA translation rates cells may adapt to stress to promote survival. Phosphorylation of eIF2α at serine 51 is one of the pathways controlling mRNA translation. Here we investigated the role of phosphorylated eIF2α during exposure to the environmental carcinogen benzo(a)pyrene (BaP). For our study we used mouse embryonic fibroblasts with a wild type eIF2α (MEF WT) and mouse embryonic fibroblasts with an eIF2α S51A knock-in mutation that cannot be phosphorylated. Here, we show that eIF2α phosphorylation occurs in MEF WT cells but not in MEF S51A cells. Survival of MEF S51A cells is profoundly reduced compared to MEF WT controls after BaP exposure. No differences in DNA damage or ROS production were observed between MEF WT and S51A cells. Disruption of eIF2α phosphorylation caused increased levels of apoptosis in response to BaP. This work demonstrates that eIF2α phosphorylation is important for reducing apoptosis and promoting cell survival in order to adapt to BaP exposure., (Copyright © 2018 The Authors. Published by Elsevier Ltd.. All rights reserved.)

PCBs are classified as xenoestrogens and carcinogens and their health risks may be sex-specific. To identify potential sex-specific responses to PCB-exposure we established gene expression profiles in a population study subdivided into females and males. Gene expression profiles were determined in a study population consisting of 512 subjects from the EnviroGenomarkers project, 217 subjects who developed lymphoma and 295 controls were selected in later life. We ran linear mixed models in order to find associations between gene expression and exposure to PCBs, while correcting for confounders, in particular distribution of white blood cells (WBC), as well as random effects. The analysis was subdivided according to sex and development of lymphoma in later life. The changes in gene expression as a result of exposure to the six studied PCB congeners were sex- and WBC type specific. The relatively large number of genes that are significantly associated with PCB-exposure in the female subpopulation already indicates different biological response mechanisms to PCBs between the two sexes. The interaction analysis between different PCBs and WBCs provides only a small overlap between sexes. In males, cancer-related pathways and in females immune system-related pathways are identified in association with PCBs and WBCs. Future lymphoma cases and controls for both sexes show different responses to the interaction of PCBs with WBCs, suggesting a role of the immune system in PCB-related cancer development.

The discovery of the epigenetic regulation of transcription has provided a new source of mechanistic understanding to long lasting effects of chemicals. However, this information is still seldom exploited in a toxicological context and studies of chemical effect after washout remain rare. Here we studied the effects of two nephrocarcinogens on the human proximal tubule cell line RPTEC/TERT1 using high-content mRNA microarrays coupled with miRNA, histone acetylation (HA) and DNA methylation (DM) arrays and metabolomics during a 5-day repeat-dose exposure and 3 days after washout. The mycotoxin ochratoxin A (OTA) was chosen as a model compound for its known impact on HA and DM. The foremost effect observed was the modulation of thousands of mRNAs and histones by OTA during and after exposure. In comparison, the oxidant potassium bromate (KBrO 3 ) had a milder impact on gene expression and epigenetics. However, there was no strong correlation between epigenetic modifications and mRNA changes with OTA while with KBrO 3 the gene expression data correlated better with HA for both up- and down-regulated genes. Even when focusing on the genes with persistent epigenetic modifications after washout, only half were coupled to matching changes in gene expression induced by OTA, suggesting that while OTA causes a major effect on the two epigenetic mechanisms studied, these alone cannot explain its impact on gene expression. Mechanistic analysis confirmed the known activation of Nrf2 and p53 by KBrO 3 , while OTA inhibited most of the same genes, and genes involved in the unfolded protein response. A few miRNAs could be linked to these effects of OTA, albeit without clear contribution of epigenetics to the modulation of the pathways at large. Metabolomics revealed disturbances in amino acid balance, energy catabolism, nucleotide metabolism and polyamine metabolism with both chemicals. In conclusion, the large impact of OTA on transcription was confirmed at the mRNA level but also with two high-content epigenomic methodologies. Transcriptomic data confirmed the previously reported activation (by KBrO 3 ) and inhibition (by OTA) of protective pathways. However, the integration of omic datasets suggested that HA and DM were not driving forces in the gene expression changes induced by either chemical.

Background: Obesity is an established risk factor for several common chronic diseases such as breast and colorectal cancer, metabolic and cardiovascular diseases; however, the biological basis for these relationships is not fully understood. To explore the association of obesity with these conditions, we investigated peripheral blood leucocyte (PBL) DNA methylation markers for adiposity and their contribution to risk of incident breast and colorectal cancer and myocardial infarction., Methods: DNA methylation profiles (Illumina Infinium ® HumanMethylation450 BeadChip) from 1941 individuals from four population-based European cohorts were analysed in relation to body mass index, waist circumference, waist-hip and waist-height ratio within a meta-analytical framework. In a subset of these individuals, data on genome-wide gene expression level, biomarkers of glucose and lipid metabolism were also available. Validation of methylation markers associated with all adiposity measures was performed in 358 individuals. Finally, we investigated the association of obesity-related methylation marks with breast, colorectal cancer and myocardial infarction within relevant subsets of the discovery population., Results: We identified 40 CpG loci with methylation levels associated with at least one adiposity measure. Of these, one CpG locus (cg06500161) in ABCG1 was associated with all four adiposity measures (P = 9.07×10 - 8 to 3.27×10 -18 ) and lower transcriptional activity of the full-length isoform of ABCG1 (P = 6.00×10 -7 ), higher triglyceride levels (P = 5.37×10 - 9 ) and higher triglycerides-to-HDL cholesterol ratio (P = 1.03×10 -10 ). Of the 40 informative and obesity-related CpG loci, two (in IL2RB and FGF18) were significantly associated with colorectal cancer (inversely, P < 1.6×10 -3 ) and one intergenic locus on chromosome 1 was inversely associated with myocardial infarction (P < 1.25×10 -3 ), independently of obesity and established risk factors., Conclusion: Our results suggest that epigenetic changes, in particular altered DNA methylation patterns, may be an intermediate biomarker at the intersection of obesity and obesity-related diseases, and could offer clues as to underlying biological mechanisms.

Background: One of the potential mechanisms linking air pollution to health effects is through changes in DNA-methylation, which so far has mainly been analyzed globally or at candidate sites., Objective: We investigated the association of personal and ambient air pollution exposure measures with genome-wide DNA-methylation changes., Methods: We collected repeated 24-hour personal and ambient exposure measurements of particulate matter (PM 2.5 ), PM 2.5 absorbance, and ultrafine particles (UFP) and peripheral blood samples from a panel of 157 healthy non-smoking adults living in four European countries. We applied univariate mixed-effects models to investigate the association between air pollution and genome-wide DNA-methylation perturbations at single CpG (cytosine-guanine dinucleotide) sites and in Differentially Methylated Regions (DMRs). Subsequently, we explored the association of air pollution-induced methylation alterations with gene expression and serum immune marker levels measured in the same subjects., Results: Personal exposure to PM 2.5 was associated with methylation changes at 13 CpG sites and 69 DMRs. Two of the 13 identified CpG sites (mapped to genes KNDC1 and FAM50B) were located within these DMRs. In addition, 42 DMRs were associated with personal PM 2.5 absorbance exposure, 16 DMRs with personal exposure to UFP, 4 DMRs with ambient exposure to PM 2.5 , 16 DMRs with ambient PM 2.5 absorbance exposure, and 15 DMRs with ambient UFP exposure. Correlation between methylation levels at identified CpG sites and gene expression and immune markers was generally moderate., Conclusion: This study provides evidence for an association between 24-hour exposure to air pollution and DNA-methylation at single sites and regional clusters of CpGs. Analysis of differentially methylated regions provides a promising avenue to further explore the subtle impact of environmental exposures on DNA-methylation., (Copyright © 2018 Elsevier Ltd. All rights reserved.)

Diesel vehicle emissions are the major source of genotoxic compounds in ambient air from urban areas. These pollutants are linked to risks of cardiovascular diseases, lung cancer, respiratory infections and adverse neurological effects. Biological events associated with exposure to some air pollutants are widely unknown but applying omics techniques may help to identify the molecular processes that link exposure to disease risk. Most data on health risks are related to long-term exposure, so the aim of this study is to investigate the impact of short-term exposure (two hours) to air pollutants on the blood transcriptome and microRNA expression levels. We analyzed transcriptomics and microRNA expression using microarray technology on blood samples from volunteers participating in studies in London, the Oxford Street cohort, and, in Barcelona, the TAPAS cohort. Personal exposure levels measurements of particulate matter (PM 10 , PM 2.5 ), ultrafine particles (UFPC), nitrogen oxides (NO 2 , NO and NOx), black carbon (BC) and carbon oxides (CO and CO 2 ) were registered for each volunteer. Associations between air pollutant levels and gene/microRNA expression were evaluated using multivariate normal models (MVN). MVN-models identified compound-specific expression of blood cell genes and microRNAs associated with air pollution despite the low exposure levels, the short exposure periods and the relatively small-sized cohorts. Hsa-miR-197-3p, hsa-miR-29a-3p, hsa-miR-15a-5p, hsa-miR-16-5p and hsa-miR-92a-3p are found significantly expressed in association with exposures. These microRNAs target also relevant transcripts, indicating their potential relevance in the research of omics-biomarkers responding to air pollution. Furthermore, these microRNAs are also known to be associated with diseases previously linked to air pollution exposure including several cancers such lung cancer and Alzheimer's disease. In conclusion, we identified in this study promising compound-specific mRNA and microRNA biomarkers after two hours of exposure to low levels of air pollutants during two hours that suggest increased cancer risks., (Copyright © 2018. Published by Elsevier Ltd.)

Translational control is a cellular response mechanism which initiates adaptation during various stress situations. Here, we investigated the role of translational control after benzo[a]pyrene (BaP) exposure in primary mouse hepatocytes. Translated mRNAs were separated and captured based on the number of associated ribosomes using sucrose gradients and subjected to RNA sequencing (RNAseq) to investigate translational changes. Furthermore, unseparated RNA (total RNA) was used for RNAseq to determine the transcriptional alterations. We showed that, after 24 h of exposure to 10 μM BaP, the number of genes altered by changes in mRNA translation was substantially higher compared with the number of genes altered by transcription. Although part of the BaP regulated genes were regulated by both transcription and translation, we identified genes that were uniquely regulated by mRNA translation. These mRNA transcripts encode proteins that are involved in biological processes that are not affected by transcriptional regulation. Al together this work identified a new layer of gene expression regulation that might contribute to BaP-induced carcinogenesis., (Copyright © 2018 Elsevier B.V. All rights reserved.)

Recent prospective studies have shown that dysregulation of the immune system may precede the development of B-cell lymphomas (BCL) in immunocompetent individuals. However, to date, the studies were restricted to a few immune markers, which were considered separately. Using a nested case-control study within two European prospective cohorts, we measured plasma levels of 28 immune markers in samples collected a median of 6 years before diagnosis (range 2.01-15.97) in 268 incident cases of BCL (including multiple myeloma [MM]) and matched controls. Linear mixed models and partial least square analyses were used to analyze the association between levels of immune marker and the incidence of BCL and its main histological subtypes and to investigate potential biomarkers predictive of the time to diagnosis. Linear mixed model analyses identified associations linking lower levels of fibroblast growth factor-2 (FGF-2 p = 7.2 × 10 -4 ) and transforming growth factor alpha (TGF-α, p = 6.5 × 10 -5 ) and BCL incidence. Analyses stratified by histological subtypes identified inverse associations for MM subtype including FGF-2 (p = 7.8 × 10 -7 ), TGF-α (p = 4.08 × 10 -5 ), fractalkine (p = 1.12 × 10 -3 ), monocyte chemotactic protein-3 (p = 1.36 × 10 -4 ), macrophage inflammatory protein 1-alpha (p = 4.6 × 10 -4 ) and vascular endothelial growth factor (p = 4.23 × 10 -5 ). Our results also provided marginal support for already reported associations between chemokines and diffuse large BCL (DLBCL) and cytokines and chronic lymphocytic leukemia (CLL). Case-only analyses showed that Granulocyte-macrophage colony stimulating factor levels were consistently higher closer to diagnosis, which provides further evidence of its role in tumor progression. In conclusion, our study suggests a role of growth-factors in the incidence of MM and of chemokine and cytokine regulation in DLBCL and CLL., (© 2018 The Authors International Journal of Cancer published by John Wiley & Sons Ltd on behalf of UICC.)

Batch effects are technical sources of variation introduced by the necessity of conducting gene expression analyses on different dates due to the large number of biological samples in population-based studies. The aim of this study is to evaluate the performances of linear mixed models (LMM) and Combat in batch effect removal. We also assessed the utility of adding quality control samples in the study design as technical replicates. In order to do so, we simulated gene expression data by adding "treatment" and batch effects to a real gene expression dataset. The performances of LMM and Combat, with and without quality control samples, are assessed in terms of sensitivity and specificity while correcting for the batch effect using a wide range of effect sizes, statistical noise, sample sizes and level of balanced/unbalanced designs. The simulations showed small differences among LMM and Combat. LMM identifies stronger relationships between big effect sizes and gene expression than Combat, while Combat identifies in general more true and false positives than LMM. However, these small differences can still be relevant depending on the research goal. When any of these methods are applied, quality control samples did not reduce the batch effect, showing no added value for including them in the study design., Competing Interests: The authors have declared that no competing interests exist.

We performed a multiple 'omics study by integrating data on epigenomic, transcriptomic, and proteomic perturbations associated with mitochondrial dysfunction in primary human hepatocytes caused by the liver toxicant valproic acid (VPA), to deeper understand downstream events following epigenetic alterations in the mitochondrial genome. Furthermore, we investigated persistence of cross-omics changes after terminating drug treatment. Upon transient methylation changes of mitochondrial genes during VPA-treatment, increasing complexities of gene-interaction networks across time were demonstrated, which normalized during washout. Furthermore, co-expression between genes and their corresponding proteins increased across time. Additionally, in relation to persistently decreased ATP production, we observed decreased expression of mitochondrial complex I and III-V genes. Persistent transcripts and proteins were related to citric acid cycle and β-oxidation. In particular, we identified a potential novel mitochondrial-nuclear signaling axis, MT-CO2-FN1-MYC-CPT1. In summary, this cross-omics study revealed dynamic responses of the mitochondrial epigenome to an impulse toxicant challenge resulting in persistent mitochondrial dysfunctioning. Moreover, this approach allowed for discriminating between the toxic effect of VPA and adaptation., (Copyright © 2018 Elsevier B.V. All rights reserved.)

Background: The regulatory management of chemicals and toxicants in the EU addresses hundreds of different chemicals and health hazards individually, one by one. An issue is that, so far, the possible interactions among chemicals or hazards are not considered as such. Another issue is the anticipated delay of several decades before effective protection of public health by regulatory decisions due to a time consuming process. Prenatal and early postnatal life is highly vulnerable to environmental health hazards with lifelong consequences, and a priority period for reduction of exposure. There are some initiatives regarding recommendations for pregnant women aiming at protection against one or another category of health hazard, however not validated by intervention studies., Hypothesis: Here, we aim at strengthening the management of exposure to individual health hazards during pregnancy and lactation, with protective measures in a global strategy of Environmental Hygiene. We hypothesize that such a strategy could reduce both the individual effects of harmful agents in complex mixtures and the possible interactions among them. A panel of experts should develop and endorse implementable measures towards a protective behavior. Their application is meant to be preferably as a package of measures in order to maximize protection and minimize interactions in causing adverse effects. Testing our hypothesis requires biomonitoring studies and longitudinal evaluation of health endpoints in the offspring. Favorable effects would legitimate further action towards equal opportunity access to improved environmental health., Conclusion: Environmental Hygiene is proposed as a global strategy aiming at effective protection of pregnant women, unborn children and infants against lifelong consequences of exposure to combinations of adverse lifestyle factors.

Regulation of DNA methylation plays a crucial role in biological processes and carcinogenesis. The formation of 5-hydroxymethylcytosine (5hmC) by oxidation of 5-methylcytosine (5mC) has been proposed as an intermediate of active demethylation. However, whether and how active demethylation is regulated by oxidative stress-related processes is not well understood. Here we investigated whether free oxygen radicals are capable of directly forming 5hmC and how this enhanced whole genome gene expression. We applied LC-MS/MS technology for the analysis of 5mC, 5hmC, 5-formylcytosine (5fC) and 5-hydroxymethyluracyl (5hmU) in HepG2 cells exposed to hydroxyl- and methyl radicals, formed by tert-butyl hydroperoxide (TBH) at multiple time points. We observed that TBH is able to induce a significant increase in 5hmC. A detailed evaluation of the hydroxymethylome using a combination of 5hmC-immunoprecipitation and microarrays resulted in the identification of highly dynamic modifications that appear to increase during prolonged oxidant exposure. Analyses of temporal gene expression changes in combination with network analysis revealed different subnetworks containing differentially expressed genes (DEGs) with differentially hydroxyl-methylated regions (DhMRs) in different regulatory kinases enriched with serine-threonine kinases. These serine-threonine kinases compromises MAPK14, RPSK6KA1, RIPK1, and PLK3 and were all previously identified as key-regulators in hepatocarcinogenesis and subject of study for chemotherapeutic interventions., (Copyright © 2018 Elsevier Ltd. All rights reserved.)

Valproic acid (VPA) is a very potent anti-cancer and neuro-protective drug probably by its HDAC inhibiting properties, which may cause steatosis in the liver. The present study investigates the effect of repetitive VPA treatment of primary human hepatocytes (PHH) on whole genome gene expression-, DNA methylation-, and miRNA changes, using microarrays and integrated data analyses. PHH were exposed to a non-cytotoxic dose of VPA for 5days daily which induced lipid accumulation. Part of the PHH was left untreated for 3days for studying the persistence of 'omics' changes. VPA treatment appeared to inhibit the expression of the transcription factors HNF1A and ONECUT1. HNF1A interacted with 41 differentially expressed genes of which 12 were also differentially methylated. None of the genes present in this network were regulated by a DE-miR. The subnetwork of ONECUT1 consisted of 44 differentially expressed genes of which 15 were differentially methylated, and 3 were regulated by a DE-miR. A number of genes in the networks are involved in fatty acid metabolism, and may contribute to the development of steatosis by increasing oxidative stress thereby causing mitochondrial dysfunction, and by shifting metabolism of VPA towards β-oxidation due to reduced glucuronidation. Part of the changes remained persistent after washing out of VPA, like PMAIP1 which is associated with cellular stress in liver of patients with NASH. The MMP2 gene showed the highest number of interactions with other persistently expressed genes, among which LCN2 which is a key modulator of lipid homeostasis. Furthermore, VPA modulated the expression and DNA methylation level of nuclear receptors and their target genes involved in the adverse outcome pathway of steatosis, thereby expanding our current knowledge of the pathway. In particular, VPA modulated PPARγ, and PPARα, AHR and CD36 on both the gene expression and the DNA methylation level, thereby inhibiting β-oxidation and increasing uptake of fatty acid into the hepatocytes, respectively. Overall, our integrative data analyses identified novel genes modulated by VPA, which provide more insight into the mechanisms of repeated dose toxicity of VPA, leading to steatosis., (Copyright © 2017 Elsevier B.V. All rights reserved.)

Valproic acid (VPA) is one of the most widely prescribed antiepileptic drugs in the world. Despite its pharmacological importance, it may cause liver toxicity and steatosis through mitochondrial dysfunction. The aim of this study is to further investigate VPA-induced mechanisms of steatosis by analyzing changes in patterns of methylation in nuclear DNA (nDNA) and mitochondrial DNA (mtDNA). Therefore, primary human hepatocytes (PHHs) were exposed to an incubation concentration of VPA that was shown to cause steatosis without inducing overt cytotoxicity. VPA was administered daily for 5 days, and this was followed by a 3 day washout (WO). Methylated DNA regions (DMRs) were identified by using the methylated DNA immunoprecipitation-sequencing (MeDIP-seq) method. The nDNA DMRs after VPA treatment could indeed be classified into oxidative stress- and steatosis-related pathways. In particular, networks of the steatosis-related gene EP300 provided novel insight into the mechanisms of toxicity induced by VPA treatment. Furthermore, we suggest that VPA induces a crosstalk between nDNA hypermethylation and mtDNA hypomethylation that plays a role in oxidative stress and steatosis development. Although most VPA-induced methylation patterns appeared reversible upon terminating VPA treatment, 31 nDNA DMRs (including 5 zinc finger protein genes) remained persistent after the WO period. Overall, we have shown that MeDIP-seq analysis is highly informative in disclosing novel mechanisms of VPA-induced toxicity in PHHs. Our results thus provide a prototype for the novel generation of interesting methylation biomarkers for repeated dose liver toxicity in vitro.

Perturbation of biological networks is often observed during exposure to xenobiotics, and the identification of disturbed processes, their dynamic traits, and dose-response relationships are some of the current challenges for elucidating the mechanisms determining adverse outcomes. In this scenario, reverse engineering of gene regulatory networks (GRNs) from expression data may provide a system-level snapshot embedded within accurate molecular events. Here, we investigate the composition of GRNs inferred from groups of chemicals with two distinct outcomes, namely carcinogenicity [azathioprine (AZA) and cyclophosphamide (CYC)] and drug-induced liver injury (DILI; diclofenac, nitrofurantoin, and propylthiouracil), and a non-carcinogenic/non-DILI group (aspirin, diazepam, and omeprazole). For this, we analyzed publicly available exposed in vitro human data, taking into account dose and time dependencies. Dose-Time Network Identification (DTNI) was applied to gene sets from exposed primary human hepatocytes using four stress pathways, namely endoplasmic reticulum (ER), NF-κB, NRF2, and TP53. Inferred GRNs suggested case specificity, varying in interactions, starting nodes, and target genes across groups. DILI and carcinogenic compounds were shown to directly affect all pathway-based GRNs, while non-DILI/non-carcinogenic chemicals only affected NF-κB. NF-κB-based GRNs clearly illustrated group-specific disturbances, with the cancer-related casein kinase CSNK2A1 being a target gene only in the carcinogenic group, and opposite regulation of NF-κB subunits being observed in DILI and non-DILI/non-carcinogenic groups. Target genes in NRF2-based GRNs shared by DILI and carcinogenic compounds suggested markers of hepatotoxicity. Finally, we indicate several of these group-specific interactions as potentially novel. In summary, our reversed-engineered GRNs are capable of revealing dose dependent, chemical-specific mechanisms of action in stress-related biological networks.

Drug-induced liver injury remains the most common cause of acute liver failure and a frequently indicated reason for withdrawal of drugs. For the purpose of evaluating the relevance of liver cell models for assessing hepatotoxic risks in intact humans, we here aimed to benchmark 'omics-derived mechanistic data from three in vitro models for parenchymal liver function, intended for the investigation of drug-induced cholestasis, against 'omics data from cholestatic patients. Transcriptomic changes in HepG2 cells, primary mouse hepatocytes and primary human hepatocytes exposed to known cholestatic compounds were analyzed using microarrays. Some of the differentially expressed genes in HepG2 cells were also differentially expressed into the same direction in human cholestasis. The overlap between drug-induced transcriptomic responses in primary mouse hepatocytes and primary human hepatocytes appeared limited and no genes overlapping with in vivo cholestasis were found. Thereupon, a pathway for drug-induced cholestasis was used to map the drug-induced transcriptomic modifications involved in bile salt homeostasis. Indications of an adaptive response to prevent and reduce intracellular bile salt accumulation were observed in vivo as well as in the in vitro liver models. Furthermore, drug-specific changes were found, which may be indicative for their cholestatic properties. Furthermore, connectivity mapping was applied in order to investigate the predictive value of the in vitro models for in vivo cholestasis. This analysis resulted in a positive connection score for most compounds, which may indicate that for identifying cholestatic compounds the focus should be on gene expression signatures rather than on differentially expressed genes., (Copyright © 2017 The Authors. Published by Elsevier Ltd.. All rights reserved.)

The mechanisms of idiosyncratic drug-induced hepatotoxicity remain largely unclear. It has demonstrated that the drug idiosyncrasy is potentiated in the context of inflammation and intracellular ceramides may play a role in this process. To study the mechanisms, HepG2 cells were co-treated with high and low doses of three idiosyncratic (I) and three non-idiosyncratic (N) compounds, with (I+ and N+) or without (I- and N-) a cytokine mix. Microarray, lipidomics and flow cytometry were performed to investigate the genome-wide expression patterns, the intracellular ceramide levels and the induction of apoptosis. We found that all I+ treatments significantly influenced the immune response- and response to stimulus-associated gene ontology (GO) terms, but the induction of apoptotic pathways, which was confirmed by flow cytometry, only appeared to be induced after the high-dose treatment. The ceramide signaling-, ER stress-, NF-kB activation- and mitochondrial activity-related pathways were biologically involved in apoptosis induced by the high-dose I+. Additionally, genes participating in ceramide metabolism were significantly altered resulting in a measurable increase in ceramide levels. The increases in ceramide concentrations may induce ER stress and activate the JNK pathway by affecting the expression of the related genes, and eventually trigger the mitochondria-independent apoptosis in hepatocytes. Overall, our study provides a potential mechanism to explain the role of inflammation in idiosyncratic drug reactions., (Copyright © 2017 Elsevier Inc. All rights reserved.)

Background: B-cell chronic lymphocytic leukemia (CLL) is a common type of adult leukemia. It often follows an indolent course and is preceded by monoclonal B-cell lymphocytosis, an asymptomatic condition, however it is not known what causes subjects with this condition to progress to CLL. Hence the discovery of prediagnostic markers has the potential to improve the identification of subjects likely to develop CLL and may also provide insights into the pathogenesis of the disease of potential clinical relevance., Results: We employed peripheral blood buffy coats of 347 apparently healthy subjects, of whom 28 were diagnosed with CLL 2.0-15.7 years after enrollment, to derive for the first time genome-wide DNA methylation, as well as gene and miRNA expression, profiles associated with the risk of future disease. After adjustment for white blood cell composition, we identified 722 differentially methylated CpG sites and 15 differentially expressed genes (Bonferroni-corrected p < 0.05) as well as 2 miRNAs (FDR < 0.05) which were associated with the risk of future CLL. The majority of these signals have also been observed in clinical CLL, suggesting the presence in prediagnostic blood of CLL-like cells. Future CLL cases who, at enrollment, had a relatively low B-cell fraction (<10%), and were therefore less likely to have been suffering from undiagnosed CLL or a precursor condition, showed profiles involving smaller numbers of the same differential signals with intensities, after adjusting for B-cell content, generally smaller than those observed in the full set of cases. A similar picture was obtained when the differential profiles of cases with time-to-diagnosis above the overall median period of 7.4 years were compared with those with shorted time-to-disease. Differentially methylated genes of major functional significance include numerous genes that encode for transcription factors, especially members of the homeobox family, while differentially expressed genes include, among others, multiple genes related to WNT signaling as well as the miRNAs miR-150-5p and miR-155-5p., Conclusions: Our findings demonstrate the presence in prediagnostic blood of future CLL patients, more than 10 years before diagnosis, of CLL-like cells which evolve as preclinical disease progresses, and point to early molecular alterations with a pathogenetic potential.

DNA damage mediates widespread changes in transcription through activation or repression of transcription factors (TFs). However, the consequences of regulating specific TFs for the outcome of the DNA repair process remain incompletely understood. Here, we combined transcriptomics and TF binding prediction with functional genomics to identify TFs essential for adequate DNA repair in HepG2 liver cells after a non-cytotoxic dose of carcinogens benzo(a)pyrene (BaP) (2μM) and aflatoxin B1 (AFB1) (5μM). BaP and AFB1 induced a largely common transcriptional response, mediated by similar TFs. A lentiviral shRNA screen knocking down the top31 identified TFs, was performed to determine their effect on DNA repair by assessing phosphorylation of H2AX (γ-H2AX). In addition to the top candidate p53, we identified several other interesting TFs that modulated γ-H2AX after BaP and AFB1 treatment. Validation studies confirmed the role of p53 in reducing γ-H2AX formation and DNA breaks measured by COMET assay after BaP and AFB1 exposure. Expression of the cell cycle inhibitor p21 was profoundly impaired upon p53 knock-down. In addition, the expression of 2 genes involved in nucleotide exchange repair, DDB2 and XPC was significantly reduced in p53 knock-down cells. Although p63 knock-down affected DNA damage upon BaP treatment this was not associated with altered expression of DDB2 or XPC. Finally, knock-down of ARNT reduced γ-H2AX in response to BaP, which was associated with reduced CYP1A1 expression. Importantly, our results suggest a new role for ARNT and its dimerization partner AHR in the occurrence of H2AX phosphorylation after AFB1 treatment. These data show that modulation of TF activity impacts on the repair of BaP- and AFB1-induced DNA damage. Our study also demonstrates the potential of combining functional genomics with genome-wide expression analysis to identify yet unknown causal relationships, thereby aiding in the interpretation of complex biological systems., (Copyright © 2017 Elsevier B.V. All rights reserved.)

Transcription factors (TFs) are important modulators of the inducible portion of the transcriptome, and therefore relevant in the context of exposure to exogenous compounds. Current approaches to predict the activity of TFs in biological systems are usually restricted to a few entities at a time due to low-throughput techniques targeting a limited fraction of annotated human TFs. Therefore, high-throughput alternatives may help to identify new targets of mechanistic and predictive value in toxicological investigations. In this study, we inferred the activity multiple TFs using publicly available microarray data from primary human hepatocytes exposed to hundreds of chemicals and evaluated these molecular profiles using multiple correspondence analysis. Our results demonstrate that the lowest dose and latest exposure time (24h) in a subset of chemicals generates a signature indicative of carcinogenicity possibly due to DNA-damaging properties. Furthermore, profiles from the earliest exposure time (2h) and highest dose creates clusters of chemicals implicated in the development of diverse forms of drug-induced liver injury (DILI). Both approaches yielded a number of TFs with similar activity across groups of chemicals, including TFs known in toxicological responses such as AhR, NFE2L2 (Nrf2), NF-κB and PPARG. FOXM1, IRF1 and E2F4 were some of the TFs identified that may be relevant in genotoxic carcinogenesis. SMADs (SMAD1, SMAD2, SMAD5) and KLF5 were identified as some of potentially new TFs whose inferred activities were linked to acute and progressive outcomes in DILI. In conclusion this study offers a novel mechanistic approach targeting TF activity during chemical exposure., (Copyright © 2017 Elsevier B.V. All rights reserved.)

The inflammatory stress has been associated with an increase in susceptibility to idiosyncratic drug-induced liver injury (DILI). However, the molecular mechanisms of this inflammation-associated idiosyncratic drug hepatotoxicity remain unknown. We exposed HepG2 cells with high and low doses of three idiosyncratic (I) and three non-idiosyncratic (N) compounds, in the presence (I+ and N+) or absence (I- and N-) of a cytokine mix for 6, 12 and 24 h. To investigate the genome-wide expression patterns, microarray was performed using the Agilent 4×44K Whole Human Genome chips. The data presented in this DIB include the expression of genes participating in the ceramide metabolism, ER stress, apoptosis and cell survival pathways. The functions of these genes were illustrated in our associated article (Jiang et al., 2017) [1]. Raw and normalized gene expression data are available through NCBI GEO (accession number GSE102006).

Unravelling gene regulatory networks (GRNs) influenced by chemicals is a major challenge in systems toxicology. Because toxicant-induced GRNs evolve over time and dose, the analysis of global gene expression data measured at multiple time points and doses will provide insight in the adverse effects of compounds. Therefore, there is a need for mathematical methods for GRN identification from time-over-dose-dependent data. One of the current approaches for GRN inference is Time Series Network Identification (TSNI). TSNI is based on ordinary differential equations (ODE), describing the time evolution of the expression of each gene, which is assumed to be dependent on the expression of other genes and an external perturbation (i.e. chemical exposure). Here, we present Dose-Time Network Identification (DTNI), a method extending TSNI by including ODE describing how the expression of each gene evolves with dose, which is supposed to depend on the expression of other genes and the exposure time. We also adapted TSNI in order to enable inclusion of time-over-dose-dependent data from multiple compounds. Here, we show that DTNI outperforms TSNI in inferring a toxicant-induced GRN. Moreover, we show that DTNI is a suitable method to infer a GRN dose- and time-dependently induced by a group of compounds influencing a common biological process. Applying DTNI on experimental data from TG-GATEs, we demonstrate that DTNI provides in-depth information on the mode of action of compounds, in particular key events and potential molecular initiating events. Furthermore, DTNI also discloses several unknown interactions which have to be verified experimentally.