Integrative omics data analyses of repeated dose toxicity of valproic acid in vitro reveal new mechanisms of steatosis induction.

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.

Integrative "-Omics" Analysis in Primary Human Hepatocytes Unravels Persistent Mechanisms of Cyclosporine A-Induced Cholestasis.

Cyclosporine A (CsA) is an undecapeptide with strong immunosuppressant activities and is used a lot after organ transplantation. Furthermore, it may induce cholestasis in the liver. In general, the drug-induced cholestasis (DIC) pathway includes genes involved in the uptake, synthesis, conjugation, and secretion of bile acids. However, whether CsA-induced changes in the cholestasis pathway in vitro are persistent for repeated dose toxicity has not yet been investigated. To explore this, primary human hepatocytes (PHH) were exposed to a subcytotoxic dose of 30 µM CsA daily for 3 and 5 days. To investigate the persistence of induced changes upon terminating CsA exposure after 5 days, a subset of PHH was subjected to a washout period (WO-period) of 3 days. Multiple -omics analyses, comprising whole genome analysis of DNA methylation, gene expression, and microRNA expression, were performed. The CsA-treatment resulted after 3 and 5 days, respectively, in 476 and 20 differentially methylated genes (DMGs), 1353 and 1481 differentially expressed genes (DEGs), and in 22 and 29 differentially expressed microRNAs (DE-miRs). Cholestasis-related pathways appeared induced during CsA-treatment. Interestingly, 828 persistent DEGs and 6 persistent DE-miRs but no persistent DMGs were found after the WO-period. These persistent DEGs and DE-miRs showed concordance for 22 genes. Furthermore, 29 persistent DEGs changed into the same direction as observed in livers from cholestasis patients. None of those 29 DEGs which among others relate to oxidative stress and lipid metabolism are yet present in the DIC pathway or cholestasis adverse outcome pathway (AOP) thus presenting novel findings. In summary, we have demonstrated for the first time a persistent impact of repeated dose administration of CsA on genes and microRNAs related to DIC in the gold standard human liver in vitro model with PHH.

Dynamic Interplay between the Transcriptome and Methylome in Response to Oxidative and Alkylating Stress.

In recent years, it has been shown that free radicals not only react directly with DNA but also regulate epigenetic processes such as DNA methylation, which may be relevant within the context of, for example, tumorigenesis. However, how these free radicals impact the epigenome remains unclear. We therefore investigated whether methyl and hydroxyl radicals, formed by tert-butyl hydroperoxide (TBH), change temporal DNA methylation patterns and how this interferes with genome-wide gene expression. At three time points, TBH-induced radicals in HepG2 cells were identified by electron spin resonance spectroscopy. Total 5-methylcytosine (5mC) levels were determined by liquid chromatography and tandem mass spectrometry and genome-wide changes in 5mC and gene expression by microarrays. Induced methylome changes rather represent an adaptive response to the oxidative stress-related reactions observed in the transcriptome. More specifically, we found that methyl radicals did not induce DNA methylation directly. An initial oxidative and alkylating stress-related response of the transcriptome during the early phase of TBH treatment was followed by an epigenetic response associated with cell survival signaling. Also, we identified genes of which the expression seems directly regulated by DNA methylation. This work suggests an important role of the methylome in counter-regulating primary oxidative and alkylating stress responses in the transcriptome to restore normal cell function. Altogether, the methylome may play an important role in counter-regulating primary oxidative and alkylating stress responses in the transcriptome presumably to restore normal cell function.

Aflatoxin B1 induces persistent epigenomic effects in primary human hepatocytes associated with hepatocellular carcinoma.

Chronic exposure to aflatoxin B1 (AFB1) has, in certain regions in the world, been strongly associated with hepatocellular carcinoma (HCC) development. AFB1 is a very potent hepatotoxic and carcinogenic mycotoxin which is frequently reported as a food contaminant. Epigenetic modifications provoked by environmental exposures, such as AFB1, may create a persistent epigenetic footprint. Deregulation of epigenetic mechanisms has actually been reported in HCC patients following AFB1 exposure; however, no attempts have yet been made to investigate early effects on the epigenome level which may be persistent on longer term, thereby possibly initiating carcinogenic events. In this study, we aim to identify methyl DNA-mRNA-interactions representative for a persistent epigenetic footprint associated with the early onset of AFB1-induced HCC. For this, primary human hepatocytes were exposed to 0.3µM of AFB1 for 5 days. Persistent epigenetic effects were measured 3 days after terminating the carcinogenic exposure. Whole genome DNA methylation changes and whole genome transcriptomic analysis were analyzed applying microarray technologies, and cross-omics interactions were evaluated. Upon combining transcriptomics data with results on DNA methylation, a range of persistent hyper- and hypo-methylated genes was identified which also appeared affected on the transcriptome level. For six of the hypo-methylated and up-regulated genes, namely TXNRD1, PCNA, CCNK, DIAPH3, RAB27A and HIST1H2BF, a clear role in carcinogenic events could be identified. This study is the first to report on a carcinogen-induced persistent impact on the epigenetic footprint in relation with the transcriptome which could be indicative for the early onset of AFB1-related development of HCC.

Epigenetic mechanisms underlying arsenic-associated lung carcinogenesis.

Arsenic is an established human carcinogen, but the mechanisms through which it contributes to for instance lung cancer development are still unclear. As arsenic is methylated during its metabolism, it may interfere with the DNA methylation process, and is therefore considered to be an epigenetic carcinogen. In the present study, we hypothesize that arsenic is able to induce DNA methylation changes, which lead to changes in specific gene expression, in pathways associated with lung cancer promotion and progression. A549 human adenocarcinoma lung cells were exposed to a low (0.08 µM), intermediate (0.4 µM) and high (2 µM) concentration of sodium arsenite for 1, 2 and 8 weeks. DNA was isolated for whole-genome DNA methylation analyses using NimbleGen 2.1 M deluxe promoter arrays. In addition, RNA was isolated for whole-genome transcriptomic analysis using Affymetrix microarrays. Arsenic modulated DNA methylation and expression levels of hundreds of genes in a dose-dependent and time-dependent manner. By combining whole-genome DNA methylation and gene expression data with possibly involved transcription factors, a large molecular interaction network was created based on transcription factor-target gene pairs, consisting of 216 genes. A tumor protein p53 (TP53) subnetwork was identified, showing the interactions of TP53 with other genes affected by arsenic. Furthermore, multiple other new genes were discovered showing altered DNA methylation and gene expression. In particular, arsenic modulated genes which function as transcription factor, thereby affecting target genes which are known to play a role in lung cancer promotion and progression.

Evaluation of database-derived pathway development for enabling biomarker discovery for hepatotoxicity.

Current testing models for predicting drug-induced liver injury are inadequate, as they basically under-report human health risks. We present here an approach towards developing pathways based on hepatotoxicity-associated gene groups derived from two types of publicly accessible hepatotoxicity databases, in order to develop drug-induced liver injury biomarker profiles. One human liver 'omics-based and four text-mining-based databases were explored for hepatotoxicity-associated gene lists. Over-representation analysis of these gene lists with a hepatotoxicant-exposed primary human hepatocytes data set showed that human liver 'omics gene lists performed better than text-mining gene lists and the results of the latter differed strongly between databases. However, both types of databases contained gene lists demonstrating biomarker potential. Visualizing those in pathway format may aid in interpreting the biomolecular background. We conclude that exploiting existing and openly accessible databases in a dedicated manner seems promising in providing venues for translational research in toxicology and biomarker development.