Evaluation of 12 mouse marker genes in rat toxicogenomics public data, Open TG-GATEs: Discrimination of genotoxic from non-genotoxic hepatocarcinogens

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.

Quinoline is more genotoxic than 4-methylquinoline in hiHeps cells and rodent liver

Environmental pollutants, such as quinoline (QN) and 4-methylquinoline (4-MeQ), may be genotoxic and carcinogenic. Earlier studies, including in vitro genotoxicity tests, indicated that 4-MeQ is more mutagenic than QN. However, we hypothesized that the methyl group of 4-MeQ favors detoxication over bioactivation, and this factor may be overlooked in in vitro tests that do not incorporate supplementation with cofactors for enzymes that catalyze conjugation reactions. We used human induced hepatocyte cells (hiHeps), which express such enzymes, and compared the genotoxicity of 4-MeQ and QN. We also carried out an in vivo micronucleus (MN) test in rat liver, since 4-MeQ is not genotoxic in rodent bone marrow. In the Ames test and the Tk gene mutation assay, with rat S9 activation, 4-MeQ was more mutagenic than QN. However, QN induced significantly higher MN frequencies in hiHeps and rat liver than did 4-MeQ. Furthermore, QN upregulated genotoxicity marker genes much more than did 4-MeQ. We also investigated the roles of two important detoxication enzymes, UDP-glucuronosyltransferases (UGTs) and cytosolic sulfotransferases (SULTs). When hiHeps were preincubated with hesperetin (UGT inhibitor) and 2,6-dichloro-4-nitrophenol (SULT inhibitor), MN frequencies were elevated approximately 1.5-fold for 4-MeQ, whereas no significant effects were seen for QN. This study shows that QN is more genotoxic than 4-MeQ, when the roles of SULTs and UGTs in detoxication are considered and our results may improve understanding the structure-activity relationships of quinoline derivatives.

Integration of Epigenetic Mechanisms into Non-Genotoxic Carcinogenicity Hazard Assessment: Focus on DNA Methylation and Histone Modifications

Epigenetics involves a series of mechanisms that entail histone and DNA covalent modifications and non-coding RNAs, and that collectively contribute to programing cell functions and differentiation. Epigenetic anomalies and DNA mutations are co-drivers of cellular dysfunctions, including carcinogenesis. Alterations of the epigenetic system occur in cancers whether the initial carcinogenic events are from genotoxic (GTxC) or non-genotoxic (NGTxC) carcinogens. NGTxC are not inherently DNA reactive, they do not have a unifying mode of action and as yet there are no regulatory test guidelines addressing mechanisms of NGTxC. To fil this gap, the Test Guideline Programme of the Organisation for Economic Cooperation and Development is developing a framework for an integrated approach for the testing and assessment (IATA) of NGTxC and is considering assays that address key events of cancer hallmarks. Here, with the intent of better understanding the applicability of epigenetic assays in chemical carcinogenicity assessment, we focus on DNA methylation and histone modifications and review: (1) epigenetic mechanisms contributing to carcinogenesis, (2) epigenetic mechanisms altered following exposure to arsenic, nickel, or phenobarbital in order to identify common carcinogen-specific mechanisms, (3) characteristics of a series of epigenetic assay types, and (4) epigenetic assay validation needs in the context of chemical hazard assessment. As a key component of numerous NGTxC mechanisms of action, epigenetic assays included in IATA assay combinations can contribute to improved chemical carcinogen identification for the better protection of public health.

Identification of Time-Invariant Biomarkers for Non-Genotoxic Hepatocarcinogen Assessment

Non-genotoxic hepatocarcinogens (NGHCs) can only be confirmed by 2-year rodent studies. Toxicogenomics (TGx) approaches using gene expression profiles from short-term animal studies could enable early assessment of NGHCs. However, high variance in the modulation of the genes had been noted among exposure styles and datasets. Expanding from our previous strategy in identifying consensus biomarkers in multiple experiments, we aimed to identify time-invariant biomarkers for NGHCs in short-term exposure styles and validate their applicability to long-term exposure styles. In this study, nine time-invariant biomarkers, namely A2m, Akr7a3, Aqp7, Ca3, Cdc2a, Cdkn3, Cyp2c11, Ntf3, and Sds, were identified from four large-scale microarray datasets. Machine learning techniques were subsequently employed to assess the prediction performance of the biomarkers. The biomarker set along with the Random Forest models gave the highest median area under the receiver operating characteristic curve (AUC) of 0.824 and a low interquartile range (IQR) variance of 0.036 based on a leave-one-out cross-validation. The application of the models to the external validation datasets achieved high AUC values of greater than or equal to 0.857. Enrichment analysis of the biomarkers inferred the involvement of chronic inflammatory diseases such as liver cirrhosis, fibrosis, and hepatocellular carcinoma in NGHCs. The time-invariant biomarkers provided a robust alternative for NGHC prediction.

Using FFPE RNA-Seq with 12 marker genes to evaluate genotoxic and non-genotoxic rat hepatocarcinogens

Introduction

Various challenges have been overcome with regard to applying 'omics technologies for chemical risk assessments. Previously we published results detailing targeted mRNA sequencing (RNA-Seq) on a next generation sequencer using intact RNA derived from freshly frozen rat liver tissues. We successfully discriminated genotoxic hepatocarcinogens (GTHCs) from non-genotoxic hepatocarcinogens (NGTHCs) using 11 selected marker genes. Based on this, we next attempted to use formalin-fixed paraffin-embedded (FFPE) pathology specimens for RNA-Seq analyses.

Findings

In this study we performed FFPE RNA-Seq to compare a typical GTHC, 2-acetylaminofluorene (AAF) to genotoxicity equivocal p-cresidine (CRE). CRE is used as a synthetic chemical intermediate, and this compound is classified as an IARC 2B carcinogen and is mutagenic in S. typhimurium, which is non-genotoxic to rat livers as assessed by single strand DNA damage analysis. RNA-Seq was used to examine liver FFPE samples obtained from groups of five 10-week-old male F344 rats that were fed with chemicals (AAF: 0.025% and CRE: 1% in food) for 4 weeks or from controls that were fed a basal diet. We extracted RNAs from FFPE samples and RNA-Seq was performed on a MiniSeq (Illumina) using the TruSeq custom RNA panel. AAF induced remarkable differences in the expression of eight genes (Aen, Bax, Btg2, Ccng1, Gdf15, Mbd1, Phlda3 and Tubb4b) from that in the control group, while CRE only induced expression changes in Gdf15, as shown using Tukey's test. Gene expression profiles for nine genes (Aen, Bax, Btg2, Ccng1, Cdkn1a, Gdf15, Mbd1, Phlda3, and Plk2) differed.between samples treated with AAF and CRE. Finally, principal component analysis (PCA) of 12 genes (Aen, Bax, Btg2, Ccnf, Ccng1, Cdkn1a, Gdf15, Lrp1, Mbd1, Phlda3, Plk2, and Tubb4b) using our previous Open TG-GATE data plus FFPE-AAF and FFPE-CRE successfully differentiated FFPE-AAF, as GTHC, from FFPE-CRE, as NGHTC.

Conclusion

Our results suggest that FFPE RNA-Seq and PCA are useful for evaluating typical rat GTHCs and NGTHCs.