Impact of Sequencing Depth and Library Preparation on Toxicological Interpretation of RNA-Seq Data in a "Three-Sample" Scenario

A comparison of the TempO-Seq and Affymetrix microarray platform using RTqPCR validation

Next-generation risk assessment relies on mechanistic data from new approach methods, including transcriptome data. Various technologies, such as high-throughput targeted sequencing methods and microarray technologies based on hybridization with complementary probes, are used to determine differentially expressed genes (DEGs). The integration of data from different technologies requires a good understanding of the differences arising from the use of various technologies.To better understand the differences between the TempO-Seq platform and Affymetrix chip technology, whole-genome data for the volatile compound dimethylamine were compared. Selected DEGs were also confirmed using RTqPCR validation. Although the overlap of DEGs between TempO-Seq and Affymetrix was no higher than 37%, a comparison of the gene regulation in terms of log2fold changes revealed a very high concordance. RTqPCR confirmed the majority of DEGs from either platform in the examined dataset. Only a few conflicts were found (11%), while 22% were not confirmed, and 3% were not detected.Despite the observed differences between the two platforms, both can be validated using RTqPCR. Here we highlight some of the differences between the two platforms and discuss their applications in toxicology.

Transcriptome analysis of a springtail, Folsomia candida, reveals energy constraint and oxidative stress during petroleum hydrocarbon exposure

Petroleum hydrocarbon (PHC) contamination in soil is ubiquitous and poses harmful consequences to many organisms. The toxicity of PHC-impacted soil is difficult to predict due to variations in mixture composition and the impacts of natural weathering processes. Hence, high-throughput methods to assess PHC-impacted soils is required to expedite land management decisions. Next-generation sequencing is a robust tool that allows researchers to investigate the effects of contaminants on the transcriptome of organisms and identify molecular biomarkers. In this study, the effects of PHCs on conventional endpoints (i.e., survival and reproduction) and gene expression rates of a model springtail species, Folsomia candida were investigated. Age-synchronized F. candida were exposed to ecologically-relevant concentrations of soils spiked with fresh crude oil to calculate the reproductive EC25 and EC50 values using conventional toxicity testing. Soils spiked to these concentrations were then used to evaluate effects on the F. candida transcriptome over a 7-day exposure period. RNA-seq analysis found 98 and 132 differentially expressed genes when compared to the control for the EC25 and EC50 treatment groups, respectively. The majority of up-regulated genes were related to xenobiotic biotransformation reactions and oxidative stress response, while down-regulated genes coded for carbohydrate and peptide metabolic processes. Promotion of the pentose phosphate pathway was also found. Results suggest that the decreased reproduction rates of F. candida exposed to PHCs is due to energy constraints caused by inhibition of carbohydrate metabolic processes and allocation of remaining energy to detoxify xenobiotics. These findings provide insights into the molecular effects in F. candida following exposure to crude oil for seven days and highlight their potential to be used as a high-throughput screening test for PHC-contaminated sites. Adverse molecular effects can be measured as early as 24 h following exposure, whereas conventional toxicity tests may require a minimum of four weeks.

Applying genomics in regulatory toxicology: a report of the ECETOC workshop on omics threshold on non-adversity

In a joint effort involving scientists from academia, industry and regulatory agencies, ECETOC's activities in Omics have led to conceptual proposals for: (1) A framework that assures data quality for reporting and inclusion of Omics data in regulatory assessments; and (2) an approach to robustly quantify these data, prior to interpretation for regulatory use. In continuation of these activities this workshop explored and identified areas of need to facilitate robust interpretation of such data in the context of deriving points of departure (POD) for risk assessment and determining an adverse change from normal variation. ECETOC was amongst the first to systematically explore the application of Omics methods, now incorporated into the group of methods known as New Approach Methodologies (NAMs), to regulatory toxicology. This support has been in the form of both projects (primarily with CEFIC/LRI) and workshops. Outputs have led to projects included in the workplan of the Extended Advisory Group on Molecular Screening and Toxicogenomics (EAGMST) group of the Organisation for Economic Co-operation and Development (OECD) and to the drafting of OECD Guidance Documents for Omics data reporting, with potentially more to follow on data transformation and interpretation. The current workshop was the last in a series of technical methods development workshops, with a sub-focus on the derivation of a POD from Omics data. Workshop presentations demonstrated that Omics data developed within robust frameworks for both scientific data generation and analysis can be used to derive a POD. The issue of noise in the data was discussed as an important consideration for identifying robust Omics changes and deriving a POD. Such variability or "noise" can comprise technical or biological variation within a dataset and should clearly be distinguished from homeostatic responses. Adverse outcome pathways (AOPs) were considered a useful framework on which to assemble Omics methods, and a number of case examples were presented in illustration of this point. What is apparent is that high dimension data will always be subject to varying processing pipelines and hence interpretation, depending on the context they are used in. Yet, they can provide valuable input for regulatory toxicology, with the pre-condition being robust methods for the collection and processing of data together with a comprehensive description how the data were interpreted, and conclusions reached.

Short-term in vivo testing to discriminate genotoxic carcinogens from non-genotoxic carcinogens and non-carcinogens using next-generation RNA sequencing, DNA microarray, and qPCR

Next-generation RNA sequencing (RNA-Seq) has identified more differentially expressed protein-coding genes (DEGs) and provided a wider quantitative range of expression level changes than conventional DNA microarrays. JEMS·MMS·Toxicogenomics group studied DEGs with targeted RNA-Seq on freshly frozen rat liver tissues and on formalin-fixed paraffin-embedded (FFPE) rat liver tissues after 28 days of treatment with chemicals and quantitative real-time PCR (qPCR) on rat and mouse liver tissues after 4 to 48 h treatment with chemicals and analyzed by principal component analysis (PCA) as statics. Analysis of rat public DNA microarray data (Open TG-GATEs) was also performed. In total, 35 chemicals were analyzed [15 genotoxic hepatocarcinogens (GTHCs), 9 non-genotoxic hepatocarcinogens (NGTHCs), and 11 non-genotoxic non-hepatocarcinogens (NGTNHCs)]. As a result, 12 marker genes (Aen, Bax, Btg2, Ccnf, Ccng1, Cdkn1a, Gdf15, Lrp1, Mbd1, Phlda3, Plk2, and Tubb4b) were proposed to discriminate GTHCs from NGTHCs and NGTNHCs. U.S. Environmental Protection Agency studied DEGs induced by 4 known GTHCs in rat liver using DNA microarray and proposed 7 biomarker genes, Bax, Bcmp1, Btg2, Ccng1, Cdkn1a, Cgr19, and Mgmt for GTHCs. Studies involving the use of whole-transcriptome RNA-Seq upon exposure to chemical carcinogens in vivo have also been performed in rodent liver, kidney, lung, colon, and other organs, although discrimination of GTHCs from NGTHCs was not examined. Candidate genes published using RNA-Seq, qPCR, and DNA microarray will be useful for the future development of short-term in vivo studies of environmental carcinogens using RNA-Seq.

Comprehensive microRNA-seq transcriptomic profiling across 11 organs, 4 ages, and 2 sexes of Fischer 344 rats

Rat is one of the most widely-used models in chemical safety evaluation and biomedical research. However, the knowledge about its microRNA (miRNA) expression patterns across multiple organs and various developmental stages is still limited. Here, we constructed a comprehensive rat miRNA expression BodyMap using a diverse collection of 320 RNA samples from 11 organs of both sexes of juvenile, adolescent, adult and aged Fischer 344 rats with four biological replicates per group. Following the Illumina TruSeq Small RNA protocol, an average of 5.1 million 50 bp single-end reads was generated per sample, yielding a total of 1.6 billion reads. The quality of the resulting miRNA-seq data was deemed to be high from raw sequences, mapped sequences, and biological reproducibility. Importantly, aliquots of the same RNA samples have previously been used to construct the mRNA BodyMap. The currently presented miRNA-seq dataset along with the existing mRNA-seq dataset from the same RNA samples provides a unique resource for studying the expression characteristics of existing and novel miRNAs, and for integrative analysis of miRNA-mRNA interactions, thereby facilitating better utilization of rats for biomarker discovery.

R-ODAF: Omics data analysis framework for regulatory application

Despite the widespread use of transcriptomics technologies in toxicology research, acceptance of the data by regulatory agencies to support the hazard assessment is still limited. Fundamental issues contributing to this are the lack of reproducibility in transcriptomics data analysis arising from variance in the methods used to generate data and differences in the data processing. While research applications are flexible in the way the data are generated and interpreted, this is not the case for regulatory applications where an unambiguous answer, possibly later subject to legal scrutiny, is required. A reference analysis framework would give greater credibility to the data and allow the practitioners to justify their use of an alternative bioinformatic process by referring to a standard. In this publication, we propose a method called omics data analysis framework for regulatory application (R-ODAF), which has been built as a user-friendly pipeline to analyze raw transcriptomics data from microarray and next-generation sequencing. In the R-ODAF, we also propose additional statistical steps to remove the number of false positives obtained from standard data analysis pipelines for RNA-sequencing. We illustrate the added value of R-ODAF, compared to a standard workflow, using a typical toxicogenomics dataset of hepatocytes exposed to paracetamol.