Stan Gaj

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 (R2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.

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.

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.

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.

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.

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.

Validating nutrient-related gene expression changes from microarrays using RT(2) PCR-arrays.

Microarray technology allows us to perform high-throughput screening of changes in gene expression. The outcome of microarray experiments largely depends on the applied analysis methods and cut-off values chosen. Results are often required to be verified using a more sensitive detection technique, such as quantitative real-time PCR (qPCR or RT-PCR). Throughout the years, this technique has become a de facto golden standard. Individual qPCRs are time-consuming, but the technology to perform high-throughput qPCR reactions has become available through PCR-arrays that allow up to 384 PCR reactions simultaneously. Our current aim was to investigate the usability of a RT2 Profiler™ PCR-array as validation in a nutritional intervention study, where the measured changes in gene expression were low. For some differentially expressed genes, the PCR-array confirmed the microarray prediction, though not for all. Furthermore, the PCR-array allowed picking up the expression of genes that were not measurable on the microarray platform but also vice versa. We conclude that both techniques have their own (dis)advantages and specificities, and for less pronounced changes using both technologies may be useful as complementation rather than validation.

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.

Linking microarray reporters with protein functions

BackgroundThe 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.ResultsThis 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.ConclusionOur 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/.