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.

The exposome concept in a human nutrigenomics study: evaluating the impact of exposure to a complex mixture of phytochemicals using transcriptomics signatures.

The application of transcriptome analyses in molecular epidemiology studies has become a promising tool in order to evaluate the impact of environmental exposures. These analyses have a great value in establishing the exposome, the totality of human exposures, both by identifying the chemical nature of the exposures and the induced molecular responses. Transcriptomic signatures can be regarded as biomarker of exposure as well as markers of effect which reflect the interaction between individual genetic background and exposure levels. However, the biological interpretation of modulated gene expression profiles is a challenging task and translating affected molecular pathways into risk assessment, for instance in terms of cancer promoting or disease preventing responses, is a far from standardised process. Here, we describe the in-depth analyses of the gene expression responses in a human dietary intervention in which the interaction between genotype and exposure to a blueberry-apple juice containing a complex mixture of phytochemicals is investigated. We also describe how data on differences in genetic background combined with different effect markers can provide a better understanding of gene-environment interactions. Pathway analyses of differentially expressed genes in combination with gene were used to identify complex but strong changes in several biological processes like immune response, cell adhesion, lipid metabolism and apoptosis. These observed changes may lead to upgraded growth control, induced immunity, reduced platelet aggregation and activation, diminished production of reactive oxidative species by platelets, blood glucose homeostasis, regulation of blood lipid levels and increased apoptosis. Our findings demonstrate that applying transcriptomics to well-controlled human dietary intervention studies can provide insight into mechanistic pathways involved in disease prevention by dietary factors.

Integrating multiple omics to unravel mechanisms of Cyclosporin A induced hepatotoxicity in vitro.

In order to improve attrition rates of candidate-drugs there is a need for a better understanding of the mechanisms underlying drug-induced hepatotoxicity. We aim to further unravel the toxicological response of hepatocytes to a prototypical cholestatic compound by integrating transcriptomic and metabonomic profiling of HepG2 cells exposed to Cyclosporin A. Cyclosporin A exposure induced intracellular cholesterol accumulation and diminished intracellular bile acid levels. Performing pathway analyses of significant mRNAs and metabolites separately and integrated, resulted in more relevant pathways for the latter. Integrated analyses showed pathways involved in cell cycle and cellular metabolism to be significantly changed. Moreover, pathways involved in protein processing of the endoplasmic reticulum, bile acid biosynthesis and cholesterol metabolism were significantly affected. Our findings indicate that an integrated approach combining metabonomics and transcriptomics data derived from representative in vitro models, with bioinformatics can improve our understanding of the mechanisms of action underlying drug-induced hepatotoxicity. Furthermore, we showed that integrating multiple omics and thereby analyzing genes, microRNAs and metabolites of the opposed model for drug-induced cholestasis can give valuable information about mechanisms of drug-induced cholestasis in vitro and therefore could be used in toxicity screening of new drug candidates at an early stage of drug discovery.

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.

The concordance between RNA-seq and microarray data depends on chemical treatment and transcript abundance.

The concordance of RNA-sequencing (RNA-seq) with microarrays for genome-wide analysis of differential gene expression has not been rigorously assessed using a range of chemical treatment conditions. Here we use a comprehensive study design to generate Illumina RNA-seq and Affymetrix microarray data from the same liver samples of rats exposed in triplicate to varying degrees of perturbation by 27 chemicals representing multiple modes of action (MOAs). The cross-platform concordance in terms of differentially expressed genes (DEGs) or enriched pathways is linearly correlated with treatment effect size (R(2)0.8). Furthermore, the concordance is also affected by transcript abundance and biological complexity of the MOA. RNA-seq outperforms microarray (93% versus 75%) in DEG verification as assessed by quantitative PCR, with the gain mainly due to its improved accuracy for low-abundance transcripts. Nonetheless, classifiers to predict MOAs perform similarly when developed using data from either platform. Therefore, the endpoint studied and its biological complexity, transcript abundance and the genomic application are important factors in transcriptomic research and for clinical and regulatory decision making.

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.

Can transcriptomics provide insight into the chemopreventive mechanisms of complex mixtures of phytochemicals in humans?

Blueberries contain relatively large amounts of different phytochemicals, which are suggested to have chemopreventive properties, but little information is available on the underlying molecular modes of action. This study investigates whole genome gene expression changes in lymphocytes of 143 humans after a 4-week blueberry-apple juice dietary intervention. Differentially expressed genes and genes correlating with the extent of antioxidant protection were identified in four subgroups. The magnitude of the preventive effect after the intervention differed between these four subgroups. Furthermore, subjects in two groups carried genetic polymorphisms that were previously found to influence the chemopreventive response. Pathway analysis of the identified genes showed strong but complex gene expression changes in pathways signaling for apoptosis, immune response, cell adhesion, and lipid metabolism. These pathways indicate increased apoptosis, upgraded growth control, induced immunity, reduced platelet aggregation and activation, blood glucose homeostasis, and regulation of fatty acid metabolism. Based on these observations, we hypothesize that combining transcriptomic data with phenotypic markers of oxidative stress may provide insight into the relevant cellular processes and genetic pathways, which contribute to the antioxidant response of complex mixtures of phytochemicals, such as found in blueberry-apple juice.

User-friendly solutions for microarray quality control and pre-processing on ArrayAnalysis.org.

Quality control (QC) is crucial for any scientific method producing data. Applying adequate QC introduces new challenges in the genomics field where large amounts of data are produced with complex technologies. For DNA microarrays, specific algorithms for QC and pre-processing including normalization have been developed by the scientific community, especially for expression chips of the Affymetrix platform. Many of these have been implemented in the statistical scripting language R and are available from the Bioconductor repository. However, application is hampered by lack of integrative tools that can be used by users of any experience level. To fill this gap, we developed a freely available tool for QC and pre-processing of Affymetrix gene expression results, extending, integrating and harmonizing functionality of Bioconductor packages. The tool can be easily accessed through a wizard-like web portal at http://www.arrayanalysis.org or downloaded for local use in R. The portal provides extensive documentation, including user guides, interpretation help with real output illustrations and detailed technical documentation. It assists newcomers to the field in performing state-of-the-art QC and pre-processing while offering data analysts an integral open-source package. Providing the scientific community with this easily accessible tool will allow improving data quality and reuse and adoption of standards.

RNA-Seq provides new insights in the transcriptome responses induced by the carcinogen benzo[a]pyrene.

Whole-genome transcriptome measurements are pivotal for characterizing molecular mechanisms of chemicals and predicting toxic classes, such as genotoxicity and carcinogenicity, from in vitro and in vivo assays. In recent years, deep sequencing technologies have been developed that hold the promise of measuring the transcriptome in a more complete and unbiased manner than DNA microarrays. Here, we applied this RNA-seq technology for the characterization of the transcriptomic responses in HepG2 cells upon exposure to benzo[a]pyrene (BaP), a well-known DNA damaging human carcinogen. Based on EnsEMBL genes, we demonstrate that RNA-seq detects ca 20% more genes than microarray-based technology but almost threefold more significantly differentially expressed genes. Functional enrichment analyses show that RNA-seq yields more insight into the biology and mechanisms related to the toxic effects caused by BaP, i.e., two- to fivefold more affected pathways and biological processes. Additionally, we demonstrate that RNA-seq allows detecting alternative isoform expression in many genes, including regulators of cell death and DNA repair such as TP53, BCL2 and XPA, which are relevant for genotoxic responses. Moreover, potentially novel isoforms were found, such as fragments of known transcripts, transcripts with additional exons, intron retention or exon-skipping events. The biological function(s) of these isoforms remain for the time being unknown. Finally, we demonstrate that RNA-seq enables the investigation of allele-specific gene expression, although no changes could be observed. Our results provide evidence that RNA-seq is a powerful tool for toxicology, which, compared with microarrays, is capable of generating novel and valuable information at the transcriptome level for characterizing deleterious effects caused by chemicals.

GO-Elite: a flexible solution for pathway and ontology over-representation.

UNLABELLED: We introduce GO-Elite, a flexible and powerful pathway analysis tool for a wide array of species, identifiers (IDs), pathways, ontologies and gene sets. In addition to the Gene Ontology (GO), GO-Elite allows the user to perform over-representation analysis on any structured ontology annotations, pathway database or biological IDs (e.g. gene, protein or metabolite). GO-Elite exploits the structured nature of biological ontologies to report a minimal set of non-overlapping terms. The results can be visualized on WikiPathways or as networks. Built-in support is provided for over 60 species and 50 ID systems, covering gene, disease and phenotype ontologies, multiple pathway databases, biomarkers, and transcription factor and microRNA targets. GO-Elite is available as a web interface, GenMAPP-CS plugin and as a cross-platform application. AVAILABILITY: http://www.genmapp.org/go_elite

Benzo[a]pyrene-induced changes in microRNA-mRNA networks.

Toxicological studies assessing the safety of compounds for humans frequently use in vitro systems to characterize toxic responses in combination with transcriptomic analyses. Thus far, changes have mostly been investigated at the mRNA level. Recently, microRNAs have attracted attention because they are powerful negative regulators of mRNA levels and, thus, may be responsible for the modulation of important mRNA networks implicated in toxicity. This study aimed to identify possible microRNA-mRNA networks as novel interactions on the gene expression level after a genotoxic insult. We used benzo[a]pyrene (BaP), a polycyclic aromatic hydrocarbon, as a model genotoxic/carcinogenic compound. We analyzed time-dependent effects on mRNA and microRNA profiles in HepG2 cells, a widely used human liver cell line that expresses active p53 and is competent for the biotransformation of BaP. Changes in microRNA expression in response to BaP, in combination with multiple alterations of mRNA levels, were observed. Many of these altered mRNAs are targets of altered microRNAs. Using pathway analysis, we evaluated the relevance of such microRNA deregulations to genotoxicity. This revealed eight microRNAs that appear to participate in specific BaP-responsive pathways relevant to genotoxicity, such as apoptotic signaling, cell cycle arrest, DNA damage response, and DNA damage repair. Our results particularly highlight the potential of microRNA-29b, microRNA-26a-1*, and microRNA-122* as novel players in the BaP response. Therefore, this study demonstrates the added value of an integrated microRNA-mRNA approach for identifying molecular mechanisms induced by BaP in an in vitro human model.

'Omics analysis of low dose acetaminophen intake demonstrates novel response pathways in humans.

Acetaminophen is the primary cause of acute liver toxicity in Europe/USA, which led the FDA to reconsider recommendations concerning safe acetaminophen dosage/use. Unfortunately, the current tests for liver toxicity are no ideal predictive markers for liver injury, i.e. they only measure acetaminophen exposure after profound liver toxicity has already occurred. Furthermore, these tests do not provide mechanistic information. Here, 'omics techniques (global analysis of metabolomic/gene-expression responses) may provide additional insight. To better understand acetaminophen-induced responses at low doses, we evaluated the effects of (sub-)therapeutic acetaminophen doses on metabolite formation and global gene-expression changes (including, for the first time, full-genome human miRNA expression changes) in blood/urine samples from healthy human volunteers. Many known and several new acetaminophen-metabolites were detected, in particular in relation to hepatotoxicity-linked, oxidative metabolism of acetaminophen. Transcriptomic changes indicated immune-modulating effects (2g dose) and oxidative stress responses (4g dose). For the first time, effects of acetaminophen on full-genome human miRNA expression have been considered and confirmed the findings on mRNA level. 'Omics techniques outperformed clinical chemistry tests and revealed novel response pathways to acetaminophen in humans. Although no definitive conclusion about potential immunotoxic effects of acetaminophen can be drawn from this study, there are clear indications that the immune system is triggered even after intake of low doses of acetaminophen. Also, oxidative stress-related gene responses, similar to those seen after high dose acetaminophen exposure, suggest the occurrence of possible pre-toxic effects of therapeutic acetaminophen doses. Possibly, these effects are related to dose-dependent increases in levels of hepatotoxicity-related metabolites.

Effect of Trichostatin A on miRNA expression in cultures of primary rat hepatocytes.

In the present study, the effect of Trichostatin A (TSA), a histone deacetylase inhibitor, was investigated on the microRNA (miR, miRNA) expression profile in cultured primary rat hepatocytes by means of microarray analysis. Simultaneously, albumin secretory capacity and morphological features of the hepatocytes were evaluated throughout the culture time. In total, 25 out of 348 miRNAs were found to be differentially expressed between freshly isolated hepatocytes and 7-day cultured cells. Nineteen of these miRNAs were connected with 'general metabolism'. miR-21 and miR-126 were shown to be the most up and down regulated miRs upon cultivation and could be linked to the proliferative response triggered in the hepatocytes upon their isolation from the liver. miR-379 and miR-143, on the other hand, were found to be the most up and down regulated miRs upon TSA treatment. Together with the higher expression of miR-122 observed in TSA-treated versus non-treated cultures, we hypothesize that the changes observed for miR-122, miR-143 and miR-379 could be related to the inhibitory effects of TSA on hepatocellular proliferation.

Biotransformation pathway maps in WikiPathways enable direct visualization of drug metabolism related expression changes.

In recent decades, our knowledge of the genetics and functional genomics of drug-metabolizing enzymes has increased and a wealth of data on drug-related 'omics' has become available. Despite the availability of large amounts of biological information on xenobiotic biotransformation, the number of available biotransformation pathway maps that can easily be used for visualization of multiple omics data is limited. Here, we created integrated biotransformation pathway maps suitable for multiple omics analysis using PathVisio. The ease of visualizing data on these maps was demonstrated by using published microarray data from human hepatocyte-like cell models, exemplifying - where a sufficient capacity for metabolizing chemicals is a prerequisite for a suited model - how the biotransformation pathway maps can be used for model selection.

Four selenoproteins, protein biosynthesis, and Wnt signalling are particularly sensitive to limited selenium intake in mouse colon.

Selenium is an essential micronutrient. Its recommended daily allowance is not attained by a significant proportion of the population in many countries and its intake has been suggested to affect colorectal carcinogenesis. Therefore, microarrays were used to determine how both selenoprotein and global gene expression patterns in the mouse colon were affected by marginal selenium deficiency comparable to variations in human dietary intakes. Two groups of 12 mice each were fed a selenium-deficient (0.086 mg Se/kg) or a selenium-adequate (0.15 mg Se/kg) diet. After 6 wk, plasma selenium level, liver, and colon glutathione peroxidase (GPx) activity in the deficient group was 12, 34, and 50%, respectively, of that of the adequate group. Differential gene expression was analysed with mouse 44K whole genome microarrays. Pathway analysis by GenMAPP identified the protein biosynthesis pathway as most significantly affected, followed by inflammation, Delta-Notch and Wnt pathways. Selected gene expression changes were confirmed by quantitative real-time PCR. GPx1 and the selenoproteins W, H, and M, responded significantly to selenium intake making them candidates as biomarkers for selenium status. Thus, feeding a marginal selenium-deficient diet resulted in distinct changes in global gene expression in the mouse colon. Modulation of cancer-related pathways may contribute to the higher susceptibility to colon carcinogenesis in low selenium status.

Reduction of colonic inflammation in HLA-B27 transgenic rats by feeding Marie Ménard apples, rich in polyphenols.

Inflammatory bowel diseases (IBD) are immunomediated ailments affecting millions of individuals. Although diet is regarded as an important factor influencing IBD, there are no accepted dietary recommendations presently available. We administered 7.6 % lyophilised apples obtained from two cultivars (Golden Delicious and Marie Ménard, low and high in polyphenols, respectively) to HLA-B27 transgenic rats which develop spontaneous IBD. After 3 months feeding, rats fed Marie Ménard apples had reduced myeloperoxidase activity (3.6 (sem 0.3) v. 2.2 (sem 0.2) U/g tissue; P < 0.05) and reduced cyclo-oxygenase-2 (P < 0.05) and inducible NO synthase gene expression (P < 0.01) in the colon mucosa and significantly less diarrhoea (P < 0.05), compared with control rats. Cell proliferation in the colon mucosa was reduced significantly by feeding Golden Delicious apples, with a borderline effect of Marie Ménard apples. Gene expression profiling of the colon mucosa, analysed using the Whole Rat Genome 4 x 44 K Agilent Arrays, revealed a down-regulation of the pathways of PG synthesis, mitogen-activated protein kinase (MAPK) signalling and TNFalpha-NF-kappaB in Marie Ménard-fed rats. In the stools of the animals of this group we also measured a significant reduction of bacteria of the Bacteriodes fragilis group. In conclusion, the administration of Marie Ménard apples, rich in polyphenols and used at present only in the manufacturing of cider, ameliorates colon inflammation in transgenic rats developing spontaneous intestinal inflammation, suggesting the possible use of these and other apple varieties to control inflammation in IBD patients.

Linking microarray reporters with protein functions.

BACKGROUND: The analysis of microarray experiments requires accurate and up-to-date functional annotation of the microarray reporters to optimize the interpretation of the biological processes involved. Pathway visualization tools are used to connect gene expression data with existing biological pathways by using specific database identifiers that link reporters with elements in the pathways. RESULTS: This paper proposes a novel method that aims to improve microarray reporter annotation by BLASTing the original reporter sequences against a species-specific EMBL subset, that was derived from and crosslinked back to the highly curated UniProt database. The resulting alignments were filtered using high quality alignment criteria and further compared with the outcome of a more traditional approach, where reporter sequences were BLASTed against EnsEMBL followed by locating the corresponding protein (UniProt) entry for the high quality hits. Combining the results of both methods resulted in successful annotation of > 58% of all reporter sequences with UniProt IDs on two commercial array platforms, increasing the amount of Incyte reporters that could be coupled to Gene Ontology terms from 32.7% to 58.3% and to a local GenMAPP pathway from 9.6% to 16.7%. For Agilent, 35.3% of the total reporters are now linked towards GO nodes and 7.1% on local pathways. CONCLUSION: Our methods increased the annotation quality of microarray reporter sequences and allowed us to visualize more reporters using pathway visualization tools. Even in cases where the original reporter annotation showed the correct description the new identifiers often allowed improved pathway and Gene Ontology linking. These methods are freely available at http://www.bigcat.unimaas.nl/public/publications/Gaj_Annotation/.