Increasing confidence in new approach methodologies for inhalation risk assessment with multiple end point assays using 5-day repeated exposure to 1,3-dichloropropene

New Approach Methodologies (NAMs) are being widely used to reduce, refine, and replace, animal use in studying toxicology. For respiratory toxicology, this includes both in silico and in vitro alternatives to replace traditional in vivo inhalation studies. 1,3-Dichloropropene (1,3-DCP) is a volatile organic compound that is widely used in agriculture as a pre-planting fumigant. Short-term exposure of humans to 1,3-DCP can result in mucous membrane irritation, chest pain, headache, and dizziness. In our previous work, we exposed differentiated cells representing different parts of the respiratory epithelium to 1,3-DCP vapor, measured cytotoxicity, and did In Vitro to In Vivo Extrapolation (IVIVE). We have extended our previous study with 1,3-DCP vapors by conducting transcriptomics on acutely exposed nasal cultures and have implemented a separate 5-day repeated exposure with multiple endpoints to gain further molecular insight into our model. MucilAir™ Nasal cell culture models, representing the nasal epithelium, were exposed to six sub-cytotoxic concentrations of 1,3-DCP vapor at the air-liquid interface, and the nasal cultures were analyzed by different methodologies, including histology, transcriptomics, and glutathione (GSH) -depletion assays. We observed the dose-dependent effect of 1,3-DCP in terms of differential gene expression, change in cellular morphology from pseudostratified columnar epithelium to squamous epithelium, and depletion of GSH in MucilAir™ nasal cultures. The MucilAir™ nasal cultures were also exposed to 3 concentrations of 1,3-DCP using repeated exposure 4 h per day for 5 days and the histological analyses indicated changes in cellular morphology and a decrease in ciliated bodies and an increase in apoptotic bodies, with increasing concentrations of 1,3-DCP. Altogether, our results suggest that sub-cytotoxic exposures to 1,3-DCP lead to several molecular and cellular perturbations, providing significant insight into the mode-of-action (MoA) of 1,3-DCP using an innovative NAM model.

Adopting duplex sequencing technology for genetic toxicity testing: A proof-of-concept mutagenesis experiment with N-ethyl-N-nitrosourea (ENU)-exposed rats

Duplex sequencing (DS) is an error-corrected next-generation sequencing method in which molecular barcodes informatically link PCR-copies back to their source DNA strands, enabling computational removal of errors in consensus sequences. The resulting background of less than one artifactual mutation per 107 nucleotides allows for direct detection of somatic mutations. TwinStrand Biosciences, Inc. has developed a DS-based mutagenesis assay to sample the rat genome, which can be applied to genetic toxicity testing. To evaluate this assay for early detection of mutagenesis, a time-course study was conducted using male Hsd:Sprague Dawley SD rats (3 per group) administered a single dose of 40 mg/kg N-ethyl-N-nitrosourea (ENU) via gavage, with mutation frequency (MF) and spectrum analyzed in stomach, bone marrow, blood, and liver tissues at 3 h, 24 h, 7 d, and 28 d post-exposure. Significant increases in MF were observed in ENU-exposed rats as early as 24 h for stomach (site of contact) and bone marrow (a highly proliferative tissue) and at 7 d for liver and blood. The canonical, mutational signature of ENU was established by 7 d post-exposure in all four tissues. Interlaboratory analysis of a subset of samples from different tissues and time points demonstrated remarkable reproducibility for both MF and spectrum. These results demonstrate that MF and spectrum can be evaluated successfully by directly sequencing targeted regions of DNA obtained from various tissues⁠, a considerable advancement compared to currently used in vivo gene mutation assays.

Genotoxicity assessment in HepaRG™ cells as a new approach methodology follow up to a positive response in the human TK6 cell micronucleus assay: Naphthalene case study

We are evaluating the use of metabolically competent HepaRG™ cells combined with CometChip® for DNA damage and the micronucleus (MN) assay as a New Approach Methodology (NAM) alternative to animals for follow up genotoxicity assessment to in vitro positive genotoxic response. Naphthalene is genotoxic in human TK6 cells inducing a nonlinear dose-response for the induction of micronuclei in the presence of rat liver S9. of naphthalene. In HepaRG™ cells, naphthalene genotoxicity was assessed using either 6 (CometChip™) or 12 concentrations of naphthalene (MN assay) with the top dose used for assessment of genotoxicity for the Comet and MN assay was 1.25 and 1.74 mM respectively, corresponding to approximately 45% cell survival. In contrast to human TK6 cell with S9, naphthalene was not genotoxic in either the HepaRG™ MN assay or the Comet assay using CometChip®. The lack of genotoxicity in both the MN and comet assays in HepaRG™ cells is likely due to Phase II enzymes removing phenols preventing further bioactivation to quinones and efficient detoxication of naphthalene quinones or epoxides by glutathione conjugation. In contrast to CYP450 mediated metabolism, these Phase II enzymes are inactive in rat liver S9 due to lack of appropriate cofactors causing a positive genotoxic response. Rat liver S9-derived BMD10 over-predicts naphthalene genotoxicity when compared to the negative genotoxic response observed in HepaRG™ cells. Metabolically competent hepatocyte models like HepaRG™ cells should be considered as human-relevant NAMs for use genotoxicity assessments to reduce reliance on rodents.

Early microRNA responses in rodent liver mediated by furan exposure establish dose thresholds for later adverse outcomes

To reduce reliance on long-term in vivo studies, short-term data linking early molecular-based measurements to later adverse health effects is needed. Although transcriptional-based benchmark dose (BMDT) modeling has been used to estimate potencies and stratify chemicals based on potential to induce later-life effects, dose-responsive epigenetic alterations have not been routinely considered. Here, we evaluated the utility of microRNA (miRNA) profiling in mouse liver and blood, as well as in mouse primary hepatocytes in vitro, to indicate mechanisms of liver perturbation due to short-term exposure of the known rodent liver hepatotoxicant and carcinogen, furan. Benchmark dose modeling of miRNA measurements (BMDmiR) were compared to the referent transcriptional (BMDT) and apical (BMDA) estimates. These analyses indicate a robust dose response for 34 miRNAs to furan and involvement of p53-linked pathways in furan-mediated hepatotoxicity, supporting mRNA and apical measurements. Liver-sourced miRNAs were also altered in the blood and primary hepatocytes. Overall, these results indicate mechanistic involvement of miRNA in furan carcinogenicity and provide evidence of their potential utility as accessible biomarkers of exposure and disease.

Error-corrected duplex sequencing enables direct detection and quantification of mutations in human TK6 cells with strong inter-laboratory consistency

Error-corrected duplex sequencing (DS) enables direct quantification of low-frequency mutations and offers tremendous potential for chemical mutagenicity assessment. We investigated the utility of DS to quantify induced mutation frequency (MF) and spectrum in human lymphoblastoid TK6 cells exposed to a prototypical DNA alkylating agent, N-ethyl-N-nitrosourea (ENU). Furthermore, we explored appropriate experimental parameters for this application, and assessed inter-laboratory reproducibility. In two independent experiments in two laboratories, TK6 cells were exposed to ENU (25-200 µM) and DNA was sequenced 48, 72, and 96 h post-exposure. A DS mutagenicity panel targeting twenty 2.4-kb regions distributed across the genome was used to sample diverse, genome-representative sequence contexts. A significant increase in MF that was unaffected by time was observed in both laboratories. Concentration-response in the MF from the two laboratories was strongly positively correlated (r = 0.97). C:G>T:A, T:A>C:G, T:A>A:T, and T:A>G:C mutations increased in consistent, concentration-dependent manners in both laboratories, with high proportions of C:G>T:A at all time points. The consistent results across the three time points suggest that 48 h may be sufficient for mutation analysis post-exposure. The target sites responded similarly between the two laboratories and revealed a higher average MF in intergenic regions. These results, demonstrating remarkable reproducibility across time and laboratory for both MF and spectrum, support the high value of DS for characterizing chemical mutagenicity in both research and regulatory evaluation.

Adopting Duplex Sequencing™ Technology for Genetic Toxicity Testing: A Proof-of-Concept Mutagenesis Experiment with N-Ethyl-N-Nitrosourea (ENU)-Exposed Rats

Duplex sequencing (DuplexSeq) is an error-corrected next-generation sequencing (ecNGS) method in which molecular barcodes informatically link PCR-copies back to their source DNA strands, enabling computational removal of errors by comparing grouped strand sequencing reads. The resulting background of less than one artifactual mutation per 10 7 nucleotides allows for direct detection of somatic mutations. TwinStrand Biosciences, Inc. has developed a DuplexSeq-based mutagenesis assay to sample the rat genome, which can be applied to genetic toxicity testing. To evaluate this assay for early detection of mutagenesis, a time-course study was conducted using male Hsd:Sprague Dawley SD rats (3 per group) administered a single dose of 40 mg/kg N-ethyl-N-nitrosourea (ENU) via gavage, with mutation frequency (MF) and spectrum analyzed in stomach, bone marrow, blood, and liver tissues at 3 h, 24 h, 7 d, and 28 d post-exposure. Significant increases in MF were observed in ENU-exposed rats as early as 24 h for stomach (site of contact) and bone marrow (a highly proliferative tissue) and at 7 d for liver and blood. The canonical, mutational signature of ENU was established by 7 d post-exposure in all four tissues. Interlaboratory analysis of a subset of samples from different tissues and time points demonstrated remarkable reproducibility for both MF and spectrum. These results demonstrate that MF and spectrum can be evaluated successfully by directly sequencing targeted regions of DNA obtained from various tissues, a considerable advancement compared to currently used in vivo gene mutation assays.

HIGHLIGHTS

DuplexSeq is an ultra-accurate NGS technology that directly quantifies mutationsENU-dependent mutagenesis was detected 24 h post-exposure in proliferative tissuesMultiple tissues exhibited the canonical ENU mutation spectrum 7 d after exposureResults obtained with DuplexSeq were highly concordant between laboratoriesThe Rat-50 Mutagenesis Assay is promising for applications in genetic toxicology.

Comparison of Rat Hepatocyte 2D-Monocultures and Hepatocytes Non-Parenchymal Cell Co-Cultures for Assessing Chemical Toxicity

Liver responses are the most common endpoints used as the basis for setting exposure standards. Liver hepatocytes play a vital role in biotransformation of xenobiotics, but non-parenchymal cells (NPCs) in the liver are also involved in certain liver responses. Development of in vitro systems that more faithfully capture liver responses to reduce reliance on animals is a major focus of New Approach Methodology (NAMs). Since rodent regulatory studies are frequently the sole source safety assessment data, mode-of-action data, and used for risk assessments, in vitro rodent models that reflect in vivo responses need to be developed to reduce reliance on animal models. In the work presented in this paper, we developed a 2-D hepatocyte monoculture and 2-D liver cell co-culture system using rat liver cells. These models were assessed for conditions for short-term stability of the cultures and phenotypic and transcriptomic responses of 2 prototypic hepatotoxicants compounds - acetaminophen and phenobarbital. The optimized multi-cellular 2-D culture required use of freshly prepared hepatocytes and NPCs from a single rat, a 3:1 ratio of hepatocytes to NPCs and growth medium using 50% Complete Williams E medium (WEM) and 50% Endothelial Cell Medium (ECM). The transcriptomic responses of the 2 model systems to PB were compared to previous studies from TG-Gates on the gene expression changes in intact rats and the co-culture model responses were more representative of the in vivo responses. Transcriptomic read-outs promise to move beyond conventional phenotypic evaluations with these in vitro NAMs and provide insights about modes of action.

Key Challenges and Recommendations for In Vitro Testing of Tobacco Products for Regulatory Applications: Consideration of Test Materials and Exposure Parameters

The Institute for In Vitro Sciences (IIVS) is sponsoring a series of workshops to identify, discuss and develop recommendations for optimal scientific and technical approaches for conducting in vitro assays, to assess potential toxicity within and across tobacco and various next generation nicotine and tobacco products (NGPs), including heated tobacco products (HTPs) and electronic nicotine delivery systems (ENDS). The third workshop (24-26 February 2020) summarised the key challenges and made recommendations concerning appropriate methods of test article generation and cell exposure from combustible cigarettes, HTPs and ENDS. Expert speakers provided their research, perspectives and recommendations for the three basic types of tobacco-related test articles: i) pad-collected material (PCM); ii) gas vapour phase (GVP); and iii) whole smoke/aerosol. These three types of samples can be tested individually, or the PCM and GVP can be combined. Whole smoke/aerosol can be bubbled through media or applied directly to cells at the air-liquid interface. Summaries of the speaker presentations and the recommendations developed by the workgroup are presented. Following discussion, the workshop concluded the following: that there needs to be greater standardisation in aerosol generation and collection processes; that methods for testing the NGPs need to be developed and/or optimised, since simply mirroring cigarette smoke testing approaches may be insufficient; that understanding and quantitating the applied dose is fundamental to the interpretation of data and conclusions from each study; and that whole smoke/aerosol approaches must be contextualised with regard to key information, including appropriate experimental controls, environmental conditioning, analytical monitoring, verification and performance criteria.

The common indoor air pollutant α-pinene is metabolized to a genotoxic metabolite α-pinene oxide

1. α-Pinene caused a concentration-responsive increase in bladder hyperplasia and decrease in sperm counts in rodents following inhalation exposure. Additionally, it formed a prospective reactive metabolite, α-pinene oxide.2. To provide human relevant context for data generated in animal models and explore potential mechanism, we undertook studies to investigate the metabolism of α-pinene to α-pinene oxide and mutagenicity of α-pinene and α-pinene oxide.3. α-Pinene oxide was formed in rat and human microsomes and hepatocytes with some species differences. Based on area under the concentration versus time curves, the formation of α-pinene oxide was up to 4-fold higher in rats than in humans.4. While rat microsomes cleared α-pinene oxide faster than human microsomes, the clearance of α-pinene oxide in hepatocytes was similar between species.5. α-Pinene was not mutagenic with or without induced rat liver S9 in Salmonella typhimurium or Escherichia coli when tested up to 10,000 μg/plate while α-pinene oxide was mutagenic at ≥25 μg/plate.6. α-Pinene was metabolized to α-pinene oxide under the conditions of the bacterial mutation assay although the concentration was approximately 3-fold lower than the lowest α-pinene oxide concentration that was positive in the assay, potentially explaining the lack of mutagenicity observed with α-pinene.

A Modern Genotoxicity Testing Paradigm: Integration of the High-Throughput CometChip® and the TGx-DDI Transcriptomic Biomarker in Human HepaRG™ Cell Cultures

Higher-throughput, mode-of-action-based assays provide a valuable approach to expedite chemical evaluation for human health risk assessment. In this study, we combined the high-throughput alkaline DNA damage-sensing CometChip^® assay with the TGx-DDI transcriptomic biomarker (DDI = DNA damage-inducing) using high-throughput TempO-Seq^®, as an integrated genotoxicity testing approach. We used metabolically competent differentiated human HepaRG™ cell cultures to enable the identification of chemicals that require bioactivation to cause genotoxicity. We studied 12 chemicals (nine DDI, three non-DDI) in increasing concentrations to measure and classify chemicals based on their ability to damage DNA. The CometChip^® classified 10/12 test chemicals correctly, missing a positive DDI call for aflatoxin B1 and propyl gallate. The poor detection of aflatoxin B1 adducts is consistent with the insensitivity of the standard alkaline comet assay to bulky lesions (a shortcoming that can be overcome by trapping repair intermediates). The TGx-DDI biomarker accurately classified 10/12 agents. TGx-DDI correctly identified aflatoxin B1 as DDI, demonstrating efficacy for combined used of these complementary methodologies. Zidovudine, a known DDI chemical, was misclassified as it inhibits transcription, which prevents measurable changes in gene expression. Eugenol, a non-DDI chemical known to render misleading positive results at high concentrations, was classified as DDI at the highest concentration tested. When combined, the CometChip^® assay and the TGx-DDI biomarker were 100% accurate in identifying chemicals that induce DNA damage. Quantitative benchmark concentration (BMC) modeling was applied to evaluate chemical potencies for both assays. The BMCs for the CometChip^® assay and the TGx-DDI biomarker were highly concordant (within 4-fold) and resulted in identical potency rankings. These results demonstrate that these two assays can be integrated for efficient identification and potency ranking of DNA damaging agents in HepaRG™ cell cultures.

Recommendations for conducting the rodent erythrocyte Pig-a assay: A report from the HESI GTTC Pig-a Workgroup

The rodent Pig-a assay is a flow cytometric, phenotype-based method used to measure in vivo somatic cell mutation. An Organization for Economic Co-operation and Development (OECD) test guideline is currently being developed to support routine use of the assay for regulatory purposes (OECD project number 4.93). This article provides advice on best practices for designing and conducting rodent Pig-a studies in support of evaluating test substance safety, with a focus on the rat model. Various aspects of assay conduct, including laboratory proficiency, minimum number of animals per dose group, preferred treatment and blood sampling schedule, and statistical analysis are described.

Transcriptomic pathway and benchmark dose analysis of Bisphenol A, Bisphenol S, Bisphenol F, and 3,3',5,5'-Tetrabromobisphenol A in H9 human embryonic stem cells

Bisphenol A (BPA) is a chemical used in the manufacturing of plastics to which human exposure is ubiquitous. Numerous studies have linked BPA exposure to many adverse health outcomes prompting the replacement of BPA with various analogues including bisphenol-F (BPF) and bisphenol S (BPS). Other bisphenols are used in various consumer applications, such as 3,3',5,5'-Tetrabromobisphenol A (TBBPA), which is used as a flame retardant. Few studies to date have examined the effects of BPA and its analogues in stem cells to explore potential developmental impacts. Here we used transcriptomics to investigate similarities and differences of BPA and three of its analogues in the estrogen receptor negative, human embryonic stem cell line H9 (WA09). H9 cells were exposed to increasing concentrations of the bisphenols and analyzed using RNA-sequencing. Our data indicate that BPA, BPF, and BPS have similar potencies in inducing transcriptional changes and perturb many of the same pathways. TBBPA, the least structurally similar bisphenol of the group, exhibited much lower potency. All bisphenols robustly impacted gene expression in these cells, albeit at concentrations well above those observed in estrogen-positive cells. Overall, we provide a foundational data set against which to explore the transcriptional similarities of other bisphenols in embryonic stem cells, which may be used to assess the suitability of chemical grouping for read-across and for preliminary potency evaluation.

Identification of compound causing yellow bone discoloration following alpha-glycosyl isoquercitrin exposure in Sprague-Dawley rats

Previous rat toxicity studies of alpha-glycosyl isoquercitrin (AGIQ), a water-soluble flavonol glycoside derived from rutin, revealed systemic yellow bone discoloration. This investigative study was conducted to determine the AGIQ metabolite(s) responsible for the discoloration. Female Sprague-Dawley rats were administered dietary AGIQ at doses of 0%, 1.5%, 3.0%, or 5.0% (0, 1735.0, 3480.8, and 5873.7 mg/kg/day, respectively) for 14 days, followed by a 14- or 28-day recovery period. Measurements of quercetin in urine and quercetin, quercetin 3-O-glucuronide, kaempferol, and 3-o-methylquercetin metabolites of AGIQ in bone (femur), white and brown fat, and cerebrum samples were conducted following the exposure period and each recovery period. Gross examination of the femur revealed yellow discoloration that increased in intensity with dose and was still present in a dose-related manner following both recovery periods. Quercetin, at levels correlating with AGIQ dose, was measured in the urine following the 14-day exposure period and, at lower concentrations, 14 or 28 days following cessation of AGIQ exposure. All four metabolites were present in a dose-dependent manner in the femur following 14 days of dietary exposure; only quercetin, quercetin 3-O-glucuronide, and 3-o-methylquercetin were present during the recovery periods. Quercetin, quercetin 3-O-glucuronide, and 3-o-methylquercetin were detected in white fat (along with kaempferol), brown fat (excluding quercetin due to analytical interference), and cerebrum samples, indicating systemic availability of the metabolites. Collectively, these data implicate quercetin, quercetin 3-O-glucuronide, or 3-o-methylquercetin (or a combination thereof) as the most likely metabolite of AGIQ causing the yellow discoloration of bone in rats administered dietary AGIQ.

Use of Frozen Tissue in the Comet Assay for the Evaluation of DNA Damage

The comet assay is gaining popularity as a means to assess DNA damage in cultured cells and tissues, particularly following exposure to chemicals or other environmental stressors. Use of the comet assay in regulatory testing for genotoxic potential in rodents has been driven by adoption of an Organisation for Economic Co-operation and Development (OECD) test guideline in 2014. Comet assay slides are typically prepared from fresh tissue at the time of necropsy; however, freezing tissue samples can avoid logistical challenges associated with simultaneous preparation of slides from multiple organs per animal and from many animals per study. Freezing also enables shipping samples from the exposure facility to a different laboratory for analysis, and storage of frozen tissue facilitates deferring a decision to generate DNA damage data for a given organ. The alkaline comet assay is useful for detecting exposure-related DNA double- and single-strand breaks, alkali-labile lesions, and strand breaks associated with incomplete DNA excision repair. However, DNA damage can also result from mechanical shearing or improper sample processing procedures, confounding the results of the assay. Reproducibility in collection and processing of tissue samples during necropsies may be difficult to control due to fluctuating laboratory personnel with varying levels of experience in harvesting tissues for the comet assay. Enhancing consistency through refresher training or deployment of mobile units staffed with experienced laboratory personnel is costly and may not always be feasible. To optimize consistent generation of high quality samples for comet assay analysis, a method for homogenizing flash frozen cubes of tissue using a customized tissue mincing device was evaluated. Samples prepared for the comet assay by this method compared favorably in quality to fresh and frozen tissue samples prepared by mincing during necropsy. Moreover, low baseline DNA damage was measured in cells from frozen cubes of tissue following prolonged storage.

Sensitive CometChip assay for screening potentially carcinogenic DNA adducts by trapping DNA repair intermediates

Genotoxicity testing is critical for predicting adverse effects of pharmaceutical, industrial, and environmental chemicals. The alkaline comet assay is an established method for detecting DNA strand breaks, however, the assay does not detect potentially carcinogenic bulky adducts that can arise when metabolic enzymes convert pro-carcinogens into a highly DNA reactive products. To overcome this, we use DNA synthesis inhibitors (hydroxyurea and 1-β-d-arabinofuranosyl cytosine) to trap single strand breaks that are formed during nucleotide excision repair, which primarily removes bulky lesions. In this way, comet-undetectable bulky lesions are converted into comet-detectable single strand breaks. Moreover, we use HepaRG™ cells to recapitulate in vivo metabolic capacity, and leverage the CometChip platform (a higher throughput more sensitive comet assay) to create the 'HepaCometChip', enabling the detection of bulky genotoxic lesions that are missed by current genotoxicity screens. The HepaCometChip thus provides a broadly effective approach for detection of bulky DNA adducts.

Flow cytometric micronucleus assay and TGx-DDI transcriptomic biomarker analysis of ten genotoxic and non-genotoxic chemicals in human HepaRG™ cells

Background

Modern testing paradigms seek to apply human-relevant cell culture models and integrate data from multiple test systems to accurately inform potential hazards and modes of action for chemical toxicology. In genetic toxicology, the use of metabolically competent human hepatocyte cell culture models provides clear advantages over other more commonly used cell lines that require the use of external metabolic activation systems, such as rat liver S9. HepaRG™ cells are metabolically competent cells that express Phase I and II metabolic enzymes and differentiate into mature hepatocyte-like cells, making them ideal for toxicity testing. We assessed the performance of the flow cytometry in vitro micronucleus (MN) test and the TGx-DDI transcriptomic biomarker to detect DNA damage-inducing (DDI) chemicals in human HepaRG™ cells after a 3-day repeat exposure. The biomarker, developed for use in human TK6 cells, is a panel of 64 genes that accurately classifies chemicals as DDI or non-DDI. Herein, the TGx-DDI biomarker was analyzed by Ion AmpliSeq whole transcriptome sequencing to assess its classification accuracy using this more modern gene expression technology as a secondary objective.

Methods

HepaRG™ cells were exposed to increasing concentrations of 10 test chemicals (six genotoxic chemicals, including one aneugen, and four non-genotoxic chemicals). Cytotoxicity and genotoxicity were measured using the In Vitro MicroFlow® kit, which was run in parallel with the TGx-DDI biomarker.

Results

A concentration-related decrease in relative survival and a concomitant increase in MN frequency were observed for genotoxic chemicals in HepaRG™ cells. All five DDI and five non-DDI agents were correctly classified (as genotoxic/non-genotoxic and DDI/non-DDI) by pairing the test methods. The aneugenic agent (colchicine) yielded the expected positive result in the MN test and negative (non-DDI) result by TGx-DDI.

Conclusions

This next generation genotoxicity testing strategy is aligned with the paradigm shift occurring in the field of genetic toxicology. It provides mechanistic insight in a human-relevant cell-model, paired with measurement of a conventional endpoint, to inform the potential for adverse health effects. This work provides support for combining these assays in an integrated test strategy for accurate, higher throughput genetic toxicology testing in this metabolically competent human progenitor cell line.

Carcinogenic activity of pentabrominated diphenyl ether mixture (DE-71) in rats and mice

•

Pentabrominated diphenyl ether (PBDE) mixture was a multispecies carcinogen causing liver tumors in male and female rats and mice.

•

Hras or Ctnnb1 mutations characterized the PBDE-induced liver tumors.

•

PBDE-induced liver tumors increased with increasing PBDE exposure.

Pentabrominated diphenyl ether (PBDE) flame retardants have been phased out in Europe and in the United States, but these lipid soluble chemicals persist in the environment and are found human and animal tissues. PBDEs have limited genotoxic activity. However, in a 2-year cancer study of a PBDE mixture (DE-71) (0, 3, 15, or 50 mg/kg (rats); 0, 3, 30, or 100 mg/kg (mice)) there were treatment-related liver tumors in male and female Wistar Han rats [Crl:WI(Han) after in utero/postnatal/adult exposure, and in male and female B6C3F1 mice, after adult exposure. In addition, there was evidence for a treatment-related carcinogenic effect in the thyroid and pituitary gland tumor in male rats, and in the uterus (stromal polyps/stromal sarcomas) in female rats. The treatment-related liver tumors in female rats were unrelated to the AhR genotype status, and occurred in animals with wild, mutant, or heterozygous Ah receptor. The liver tumors in rats and mice had treatment-related Hras and Ctnnb mutations, respectively. The PBDE carcinogenic activity could be related to oxidative damage, disruption of hormone homeostasis, and molecular and epigenetic changes in target tissue. Further work is needed to compare the PBDE toxic effects in rodents and humans.

Comprehensive evaluation of the flavonol anti-oxidants, alpha-glycosyl isoquercitrin and isoquercitrin, for genotoxic potential

Quercetin and its glycosides possess potential benefits to human health. Several flavonols are available to consumers as dietary supplements, promoted as anti-oxidants; however, incorporation of natural quercetin glycosides into food and beverage products has been limited by poor miscibility in water. Enzymatic conjugation of multiple glucose moieties to isoquercitrin to produce alpha-glycosyl isoquercitrin (AGIQ) enhances solubility and bioavailability. AGIQ is used in Japan as a food additive and has been granted generally recognized as safe (GRAS) status. However, although substantial genotoxicity data exist for quercetin, there is very little available data for AGIQ and isoquercitrin. To support expanded global marketing of food products containing AGIQ, comprehensive testing of genotoxic potential of AGIQ and isoquercitrin was conducted according to current regulatory test guidelines. Both chemicals tested positive in bacterial reverse mutation assays, and exposure to isoquercitrin resulted in chromosomal aberrations in CHO-WBL cells. All other in vitro mammalian micronucleus and chromosomal aberration assays, micronucleus and comet assays in male and female B6C3F1 mice and Sprague Dawley rats, and Muta™ Mouse mutation assays evaluating multiple potential target tissues, were negative for both chemicals. These results supplement existing toxicity data to further support the safe use of AGIQ in food and beverage products.

Impact of Acrylamide on Calcium Signaling and Cytoskeletal Filaments in Testes From F344 Rat

Acrylamide (AA) at high exposure levels is neurotoxic, induces testicular toxicity, and increases dominant lethal mutations in rats. RNA-sequencing in testes was used to identify differentially expressed genes (DEG), explore AA-induced pathway perturbations that could contribute to AA-induced testicular toxicity and then used to derive a benchmark dose (BMD). Male F344/DuCrl rats were administered 0.0, 0.5, 1.5, 3.0, 6.0, or 12.0 mg AA/kg bw/d in drinking water for 5, 15, or 31 days. The experimental design used exposure levels that spanned and exceeded the exposure levels used in the rat dominant lethal, 2-generation reproductive toxicology, and cancer bioassays. The time of sample collection was based on previous studies that developed gene expression-based BMD. At 12.0 mg/kg, there were 38, 33, and 65 DEG ( P value <.005; fold change >1.5) in the testes after 5, 15, or 31 days of exposure, respectively. At 31 days, there was a dose-dependent increase in the number of DEG, and at 12.0 mg/kg/d the top three functional clusters affected by AA exposure were actin filament organization, response to calcium ion, and regulation of cell proliferation. The BMD lower 95% confidence limit using DEG ranged from 1.8 to 6.8 mg/kg compared to a no-observed-adverse-effect-level of 2.0 mg/kg/d for male reproductive toxicity. These results are consistent with the known effects of AA on calcium signaling and cytoskeletal actin filaments leading to neurotoxicity and suggest that AA can cause rat dominant lethal mutations by these same mechanisms leading to impaired chromosome segregation during cell division.

Application of the TGx-28.65 transcriptomic biomarker to classify genotoxic and non-genotoxic chemicals in human TK6 cells in the presence of rat liver S9

In vitro transcriptional signatures that predict toxicities can facilitate chemical screening. We previously developed a transcriptomic biomarker (known as TGx-28.65) for classifying agents as genotoxic (DNA damaging) and non-genotoxic in human lymphoblastoid TK6 cells. Because TK6 cells do not express cytochrome P450s, we confirmed accurate classification by the biomarker in cells co-exposed to 1% 5,6 benzoflavone/phenobarbital-induced rat liver S9 for metabolic activation. However, chemicals may require different types of S9 for activation. Here we investigated the response of TK6 cells to higher percentages of Aroclor-, benzoflavone/phenobarbital-, or ethanol-induced rat liver S9 to expand TGx-28.65 biomarker applicability. Transcriptional profiles were derived 3 to 4 hr following a 4 hr co-exposure of TK6 cells to test chemicals and S9. Preliminary studies established that 10% Aroclor- and 5% ethanol-induced S9 alone did not induce the TGx-28.65 biomarker genes. Seven genotoxic and two non-genotoxic chemicals (and concurrent solvent and positive controls) were then tested with one of the S9s (selected based on cell survival and micronucleus induction). Relative survival and micronucleus frequency was assessed by flow cytometry in cells 20 hr post-exposure. Genotoxic/non-genotoxic chemicals were accurately classified using the different S9s. One technical replicate of cells co-treated with dexamethasone and 10% Aroclor-induced S9 was falsely classified as genotoxic, suggesting caution in using high S9 concentrations. Even low concentrations of genotoxic chemicals (those not causing cytotoxicity) were correctly classified, demonstrating that TGx-28.65 is a sensitive biomarker of genotoxicity. A meta-analysis of datasets from 13 chemicals supports that different S9s can be used in TK6 cells, without impairing classification using the TGx-28.65 biomarker.

Comparison of toxicogenomics and traditional approaches to inform mode of action and points of departure in human health risk assessment of benzo[a]pyrene in drinking water

Toxicogenomics is proposed to be a useful tool in human health risk assessment. However, a systematic comparison of traditional risk assessment approaches with those applying toxicogenomics has never been done. We conducted a case study to evaluate the utility of toxicogenomics in the risk assessment of benzo[a]pyrene (BaP), a well-studied carcinogen, for drinking water exposures. Our study was intended to compare methodologies, not to evaluate drinking water safety. We compared traditional (RA1), genomics-informed (RA2) and genomics-only (RA3) approaches. RA2 and RA3 applied toxicogenomics data from human cell cultures and mice exposed to BaP to determine if these data could provide insight into BaP's mode of action (MOA) and derive tissue-specific points of departure (POD). Our global gene expression analysis supported that BaP is genotoxic in mice and allowed the development of a detailed MOA. Toxicogenomics analysis in human lymphoblastoid TK6 cells demonstrated a high degree of consistency in perturbed pathways with animal tissues. Quantitatively, the PODs for traditional and transcriptional approaches were similar (liver 1.2 vs. 1.0 mg/kg-bw/day; lungs 0.8 vs. 3.7 mg/kg-bw/day; forestomach 0.5 vs. 7.4 mg/kg-bw/day). RA3, which applied toxicogenomics in the absence of apical toxicology data, demonstrates that this approach provides useful information in data-poor situations. Overall, our study supports the use of toxicogenomics as a relatively fast and cost-effective tool for hazard identification, preliminary evaluation of potential carcinogens, and carcinogenic potency, in addition to identifying current limitations and practical questions for future work.

Mining the Archives: A Cross-Platform Analysis of Gene Expression Profiles in Archival Formalin-Fixed Paraffin-Embedded Tissues

Formalin-fixed paraffin-embedded (FFPE) tissue samples represent a potentially invaluable resource for transcriptomic research. However, use of FFPE samples in genomic studies has been limited by technical challenges resulting from nucleic acid degradation. Here we evaluated gene expression profiles derived from fresh-frozen (FRO) and FFPE mouse liver tissues preserved in formalin for different amounts of time using 2 DNA microarray protocols and 2 whole-transcriptome sequencing (RNA-seq) library preparation methodologies. The ribo-depletion protocol outperformed the other methods by having the highest correlations of differentially expressed genes (DEGs), and best overlap of pathways, between FRO and FFPE groups. The effect of sample time in formalin (18 h or 3 weeks) on gene expression profiles indicated that test article treatment, not preservation method, was the main driver of gene expression profiles. Meta- and pathway analyses indicated that biological responses were generally consistent for 18 h and 3 week FFPE samples compared with FRO samples. However, clear erosion of signal intensity with time in formalin was evident, and DEG numbers differed by platform and preservation method. Lastly, we investigated the effect of time in paraffin on genomic profiles. Ribo-depletion RNA-seq analysis of 8-, 19-, and 26-year-old control blocks resulted in comparable quality metrics, including expected distributions of mapped reads to exonic, untranslated region, intronic, and ribosomal fractions of the transcriptome. Overall, our results indicate that FFPE samples are appropriate for use in genomic studies in which frozen samples are not available, and that ribo-depletion RNA-seq is the preferred method for this type of analysis in archival and long-aged FFPE samples.

Impact of Genomics Platform and Statistical Filtering on Transcriptional Benchmark Doses (BMD) and Multiple Approaches for Selection of Chemical Point of Departure (PoD)

Many regulatory agencies are exploring ways to integrate toxicogenomic data into their chemical risk assessments. The major challenge lies in determining how to distill the complex data produced by high-content, multi-dose gene expression studies into quantitative information. It has been proposed that benchmark dose (BMD) values derived from toxicogenomics data be used as point of departure (PoD) values in chemical risk assessments. However, there is limited information regarding which genomics platforms are most suitable and how to select appropriate PoD values. In this study, we compared BMD values modeled from RNA sequencing-, microarray-, and qPCR-derived gene expression data from a single study, and explored multiple approaches for selecting a single PoD from these data. The strategies evaluated include several that do not require prior mechanistic knowledge of the compound for selection of the PoD, thus providing approaches for assessing data-poor chemicals. We used RNA extracted from the livers of female mice exposed to non-carcinogenic (0, 2 mg/kg/day, mkd) and carcinogenic (4, 8 mkd) doses of furan for 21 days. We show that transcriptional BMD values were consistent across technologies and highly predictive of the two-year cancer bioassay-based PoD. We also demonstrate that filtering data based on statistically significant changes in gene expression prior to BMD modeling creates more conservative BMD values. Taken together, this case study on mice exposed to furan demonstrates that high-content toxicogenomics studies produce robust data for BMD modelling that are minimally affected by inter-technology variability and highly predictive of cancer-based PoD doses.

Genotoxicity evaluation of the flavonoid, myricitrin, and its aglycone, myricetin

Myricitrin, a flavonoid extracted from the fruit, leaves, and bark of Chinese bayberry (Myrica rubra SIEBOLD), is currently used as a flavor modifier in snack foods, dairy products, and beverages in Japan. Myricitrin is converted to myricetin by intestinal microflora; myricetin also occurs ubiquitously in plants and is consumed in fruits, vegetables, and beverages. The genotoxic potential of myricitrin and myricetin was evaluated in anticipation of worldwide marketing of food products containing myricitrin. In a bacterial reverse mutation assay, myricetin tested positive for frameshift mutations under metabolic activation conditions whereas myricitrin tested negative for mutagenic potential. Both myricitrin and myricetin induced micronuclei formation in human TK6 lymphoblastoid cells under conditions lacking metabolic activation; however, the negative response observed in the presence of metabolic activation suggests that rat liver S9 homogenate may detoxify reactive metabolites of these chemicals in mammalian cells. In 3-day combined micronucleus/Comet assays using male and female B6C3F1 mice, no induction of micronuclei was observed in peripheral blood, or conclusive evidence of damage detected in the liver, glandular stomach, or duodenum following exposure to myricitrin or myricetin. Our studies did not reveal evidence of genotoxic potential of myricitrin in vivo, supporting its safe use in food and beverages.

Integration of metabolic activation with a predictive toxicogenomics signature to classify genotoxic versus nongenotoxic chemicals in human TK6 cells

The use of integrated approaches in genetic toxicology, including the incorporation of gene expression data to determine the molecular pathways involved in the response, is becoming more common. In a companion article, a genomic biomarker was developed in human TK6 cells to classify chemicals as genotoxic or nongenotoxic. Because TK6 cells are not metabolically competent, we set out to broaden the utility of the biomarker for use with chemicals requiring metabolic activation. Specifically, chemical exposures were conducted in the presence of rat liver S9. The ability of the biomarker to classify genotoxic (benzo[a]pyrene, BaP; aflatoxin B1, AFB1) and nongenotoxic (dexamethasone, DEX; phenobarbital, PB) agents correctly was evaluated. Cells were exposed to increasing chemical concentrations for 4 hr and collected 0 hr, 4 hr, and 20 hr postexposure. Relative survival, apoptosis, and micronucleus frequency were measured at 24 hr. Transcriptome profiles were measured with Agilent microarrays. Statistical modeling and bioinformatics tools were applied to classify each chemical using the genomic biomarker. BaP and AFB1 were correctly classified as genotoxic at the mid- and high concentrations at all three time points, whereas DEX was correctly classified as nongenotoxic at all concentrations and time points. The high concentration of PB was misclassified at 24 hr, suggesting that cytotoxicity at later time points may cause misclassification. The data suggest that the use of S9 does not impair the ability of the biomarker to classify genotoxicity in TK6 cells. Finally, we demonstrate that the biomarker is also able to accurately classify genotoxicity using a publicly available dataset derived from human HepaRG cells.

Integration of metabolic activation with a predictive toxicogenomics signature to classify genotoxic versus nongenotoxic chemicals in human TK6 cells

The use of integrated approaches in genetic toxicology, including the incorporation of gene expression data to determine the molecular pathways involved in the response, is becoming more common. In a companion article, a genomic biomarker was developed in human TK6 cells to classify chemicals as genotoxic or nongenotoxic. Because TK6 cells are not metabolically competent, we set out to broaden the utility of the biomarker for use with chemicals requiring metabolic activation. Specifically, chemical exposures were conducted in the presence of rat liver S9. The ability of the biomarker to classify genotoxic (benzo[a]pyrene, BaP; aflatoxin B1, AFB1) and nongenotoxic (dexamethasone, DEX; phenobarbital, PB) agents correctly was evaluated. Cells were exposed to increasing chemical concentrations for 4 hr and collected 0 hr, 4 hr, and 20 hr postexposure. Relative survival, apoptosis, and micronucleus frequency were measured at 24 hr. Transcriptome profiles were measured with Agilent microarrays. Statistical modeling and bioinformatics tools were applied to classify each chemical using the genomic biomarker. BaP and AFB1 were correctly classified as genotoxic at the mid- and high concentrations at all three time points, whereas DEX was correctly classified as nongenotoxic at all concentrations and time points. The high concentration of PB was misclassified at 24 hr, suggesting that cytotoxicity at later time points may cause misclassification. The data suggest that the use of S9 does not impair the ability of the biomarker to classify genotoxicity in TK6 cells. Finally, we demonstrate that the biomarker is also able to accurately classify genotoxicity using a publicly available dataset derived from human HepaRG cells.

Comet assay evaluation of six chemicals of known genotoxic potential in rats

As a part of an international validation of the in vivo rat alkaline comet assay (comet assay) initiated by the Japanese Center for the Validation of Alternative Methods (JaCVAM) we examined six chemicals for potential to induce DNA damage: 2-acetylaminofluorene (2-AAF), N-nitrosodimethylamine (DMN), o-anisidine, 1,2-dimethylhydrazine dihydrochloride (1,2-DMH), sodium chloride, and sodium arsenite. DNA damage was evaluated in the liver and stomach of 7- to 9-week-old male Sprague Dawley rats. Of the five genotoxic carcinogens tested in our laboratory, DMN and 1,2-DMH were positive in the liver and negative in the stomach, 2-AAF and o-anisidine produced an equivocal result in liver and negative results in stomach, and sodium arsenite was negative in both liver and stomach. 1,2-DMH and DMN induced dose-related increases in hedgehogs in the same tissue (liver) that exhibited increased DNA migration. However, no cytotoxicity was indicated by the neutral diffusion assay (assessment of highly fragmented DNA) or histopathology in response to treatment with any of the tested chemicals. Therefore, the increased DNA damage resulting from exposure to DMN and 1,2-DMH was considered to represent a genotoxic response. Sodium chloride, a non-genotoxic non-carcinogen, was negative in both tissues as would be predicted. Although only two (1,2-DMH and DMN) out of five genotoxic carcinogens produced clearly positive results in the comet assay, the results obtained for o-anisidine and sodium arsenite in liver and stomach cells are consistent with the known mode of genotoxicity and tissue specificity exhibited by these carcinogens. In contrast, given the known genotoxic mode-of-action and target organ carcinogenicity of 2-AAF, it is unclear why this chemical failed to convincingly increase DNA migration in the liver. Thus, the results of the comet assay validation studies conducted in our laboratory were considered appropriate for five out of the six test chemicals.

Gene expression analysis of livers from female B6C3F1 mice exposed to carcinogenic and non-carcinogenic doses of furan, with or without bromodeoxyuridine (BrdU) treatment

Standard methodology for identifying chemical carcinogens is both time-consuming and resource intensive. Researchers are actively investigating how new technologies can be used to identify chemical carcinogens in a more rapid and cost-effective manner. Here we performed a toxicogenomic case study of the liver carcinogen furan. Full study and mode of action details were previously published in the Journal of Toxicology and Applied Pharmacology. Female B6C3F1 mice were sub-chronically treated with two non-carcinogenic (1 and 2 mg/kg bw) and two carcinogenic (4 and 8 mg/kg bw) doses of furan for 21 days. Half of the mice in each dose group were also treated with 0.02% bromodeoxyuridine (BrdU) for five days prior to sacrifice [13]. Agilent gene expression microarrays were used to measure changes in liver gene and long non-coding RNA expression (published in Toxicological Sciences). Here we describe the experimental and quality control details for the microarray data. We also provide the R code used to analyze the raw data files, produce fold change and false discovery rate (FDR) adjusted p values for each gene, and construct hierarchical clustering between datasets.

Bromodeoxyuridine (BrdU) treatment to measure hepatocellular proliferation does not mask furan-induced gene expression changes in mouse liver

Bromodeoxyuridine (BrdU) is a synthetic nucleoside used to detect cellular proliferation. BrdU incorporates in the place of thymine but pairs with guanine, thereby increasing the risk of transition mutations in dividing cells. Given its mutagenicity, standard practice is to use a second cohort of animals for parallel toxicogenomics studies; however, the impact of BrdU on global gene expression is unknown. To test this, we performed a case study to determine whether the molecular mode of action of furan, a liver carcinogen, could be detected in BrdU-treated samples. We measure global hepatic gene expression using Agilent DNA microarrays in female B6C3F1 mice that were sub-chronically exposed to 0, 1, 4, or 8mg/kg bodyweight (bw) per day furan either in the presence (+BrdU) or absence (-BrdU) of BrdU. Exposure to 0.02% BrdU in drinking water for five days resulted in minimal gene expression changes. A comparison of +BrdU versus -BrdU control mice revealed only 11 probes with fold change≥1.5 and false discovery rate (FDR) corrected p≤0.05. The same comparison in the high dose group yielded only 3 differentially expressed probes. Differentially expressed gene lists generated for furan-treated versus control mice and were compared for the -BrdU and +BrdU groups. The high dose of furan had 452 shared probes and 27 and 90 unique probes for -BrdU and +BrdU groups, respectively. These differences did not impact hierarchical clustering. Further, they did not impair detection of the previously reported furan mode of action, which was well represented in the BrdU-treated samples. Taken together, we demonstrate that BrdU treatment does not mask important furan-induced transcriptional changes. We suggest that BrdU-treated mice could be used for toxicogenomic analysis, which would generally halve the number of rodents required for toxicogenomics studies. However, we also recommend that this type of case study be repeated for other chemicals before the use of BrdU-treated animals in omics studies becomes common practice.

Relative potency for altered humoral immunity induced by polybrominated and polychlorinated dioxins/furans in female B6C3F1/N mice

The use of brominated flame retardants and incineration of bromine-containing materials has lead to an increase in polybrominated dibenzo-p-dioxins and dibenzofurans (PBDD/Fs) in the environment. Measurable amounts of PBDD/Fs have been detected in soil, seafood, and human breast milk and serum. Studies indicate that the relative potencies of some PBDD/Fs based on enzyme induction are equivalent to those of some polychlorinated dibenzo-p-dioxins and dibenzofurans. To assess the humoral immunity relative potencies of PBDD/Fs and compare them to their chlorinated analogs, female B6C3F1/N mice received a single oral exposure to 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), 2,3,7,8-tetrabromodibenzofuran (TBDF), 2,3,7,8-tetrachlorodibenzofuran (TCDF), 1,2,3,7,8-pentabromodibenzofuran (1PeBDF), 1,2,3,7,8-pentachlorodibenzofuran (1PeCDF), 2,3,4,7,8-pentabromodibenzofuran (4PeBDF), 2,3,4,7,8-pentachlorodibenzofuran (4PeCDF), 2,3-dibromo-7,8-dichlorodibenzo-p-dioxin (DBDCDD), or 2,3,7-tribromodibenzo-p-dioxin (TriBDD). Inhibition of the immunoglobulin M (IgM) antibody forming cell response was measured 4 days following immunization with sheep red blood cells. The data were fit to a Hill model to estimate the ED50 for inhibition. Expression of xenobiotic metabolizing enzyme (XME) and thyroxine transport protein (Ttr) genes in liver was measured by PCR to assess aryl hydrocarbon-mediated responses. TCDD, TBDF, TCDF, 1PeBDF, 4PeBDF, 4PeCDF, and DBDCDD suppressed the IgM antibody response and Ttr gene expression, and upregulated phase I XME genes. 1PeCDF suppressed the IgM antibody response but only upregulated phase I XME genes; TriBDD had no effect on antibody response. The rank order of potency (ED50) for these chemicals was TCDD>TBDF>4PeBDF>TCDF/4PeCDF/1PeBDF>1PeCDF. Whereas TCDD was the most potent compound tested, the brominated analogs were more potent than their chlorinated analogs, suggesting that these compounds should be considered in toxic equivalency factor evaluation and risk assessment.

Cytotoxicity of naphthalene toward cells from target and non-target organs in vitro

Chronic inhalation exposure to high concentrations of naphthalene produced nasal tumors in rats and lung tumors in female mice. Naphthalene bioactivation is required for target organ toxicity and cytotoxicity in target organs may be involved in tumor development. The present studies characterized the dose-response relationships for naphthalene-induced glutathione (GSH) depletion, effects on cellular ATP, and cytotoxicity in cells from both target (lung, nasal epithelium) and non-target (liver) organs in vitro using cells from F-344 rats, B6C3F1 mice and humans. The cells were incubated with various concentrations of naphthalene in sealed glass flasks for 3h, then placed in monolayer culture in fresh media for 24h to examine the repair or progression of damage. Naphthalene was a low potency cytotoxicant in vitro, with 500 μM frequently observed as a no-observed adverse effect concentration or lowest observed adverse effect concentration. Naphthalene exposure produced dose-dependent decreases in cellular GSH, ATP and viability in rat, mouse and human hepatocytes at concentrations >500 μM. Human nasal respiratory epithelial cells exhibited greater naphthalene cytotoxicity than rat or mouse nasal respiratory epithelial cell preparations. Rat nasal respiratory epithelial cell preparations metabolized naphthalene through pathways leading to the preferential formation of 1,2-naphthoquinone GSH conjugates rather than 1,4-naphthoquinone GSH conjugates observed in rat hepatocytes or mouse nasal respiratory epithelial cells, consistent with the suggestion that this bioactivation pathway may be involved in rat nasal tumor development. Naphthalene exposures of ≥500 μM decreased cellular GSH and ATP in rat, mouse and human lung cell preparations. The variability of the responses of the human lung cell preparations was consistent with the known variability of CYP activities in human lung tissue. The results of these studies can be used as the basis for future studies of the mechanisms involved in naphthalene-induced cytotoxicity and the relevance of the bioactivation pathways for human exposure to naphthalene.

Nonlinear responses for chromosome and gene level effects induced by vinyl acetate monomer and its metabolite, acetaldehyde in TK6 cells

Vinyl acetate monomer (VAM) produced rat nasal tumors at concentrations in the hundreds of parts per million. However, VAM is weakly genotoxic in vitro and shows no genotoxicity in vivo. A European Union Risk Assessment concluded that VAM's hydrolysis to acetaldehyde (AA), via carboxylesterase, is a critical key event in VAM's carcinogenic potential. In the following study, we observed increases in micronuclei (MN) and thymidine kinase (Tk) mutants that were dependent on the ability of TK6 cell culture conditions to rapidly hydrolyze VAM to AA. Heat-inactivated horse serum demonstrated a high capacity to hydrolyze VAM to AA; this activity was highly correlated with a concomitant increase in MN. In contrast, heat-inactivated fetal bovine serum (FBS) did not hydrolyze VAM and no increase in MN was observed. AA's ability to induce MN was not impacted by either serum since it directly forms Schiff bases with DNA and proteins. Increased mutant frequency at the Tk locus was similarly mitigated when AA formation was not sufficiently rapid, such as incubating VAM in the presence of FBS for 4 hr. Interestingly, neither VAM nor AA induced mutations at the HPRT locus. Finally, cytotoxicity paralleled genotoxicity demonstrating that a small degree of cytotoxicity occurred prior to increases in MN. These results established 0.25 mM as a consistent concentration where genotoxicity first occurred for both VAM and AA provided VAM is hydrolyzed to AA. This information further informs significant key events related to the mode of action of VAM-induced nasal mucosal tumors in rats.

ACB-PCR measurement of spontaneous and furan-induced H-ras codon 61 CAA to CTA and CAA to AAA mutation in B6C3F1 mouse liver

Furan is a rodent liver carcinogen, but the mode of action for furan hepatocarcinogenicity is unclear. H-ras codon 61 mutations have been detected in spontaneous liver tumors of B6C3F1 mice, and the fraction of liver tumors carrying H-ras codon 61 CAA to AAA mutation increased in furan-treated mice. Allele-specific competitive blocker PCR (ACB-PCR) has been used previously to quantify early, carcinogen-induced increases in tumor-associated mutations. The present pilot study investigated whether furan drives clonal expansion of pre-existing H-ras mutant cells in B6C3F1 mouse liver. H-ras codon 61 CAA to CTA and CAA to AAA mutations were measured in DNA isolated from liver tissue of female mice treated with 0, 1, 2, 4, or 8 mg furan/kg body weight, five days per week for three weeks, using five mice per treatment group. Spontaneous levels of mutation were low, with two of five control mice having an H-ras codon 61 CTA or AAA mutant fraction (MF) greater than 10(-5) . Several furan-treated mice had H-ras codon 61 AAA or CTA MFs greater than those measured in control mice and lower bound estimates of induced MF were calculated. However, no statistically-significant differences were observed between treatment groups. Therefore, while sustained exposure to furan is carcinogenic, at the early stage of carcinogenesis examined in this study (three weeks), there was not a significant expansion of H-ras mutant cells.

Differential expression of long noncoding RNAs in the livers of female B6C3F1 mice exposed to the carcinogen furan

The mammalian genome is transcribed into mRNAs that code for protein and a broad spectrum of other noncoding (nc) RNA products. Long ncRNAs (lncRNA), defined as ncRNA species > 200 nucleotides long, are emerging as important epigenetic regulators of gene expression that are involved in a spectrum of biological processes of relevance to toxicology. We conducted a gene expression profiling study in the livers of female B6C3F1 mice exposed to the carcinogen furan at 0.0, 1.0, and 2.0mg/kg (noncarcinogenic doses) and at 4.0 and 8.0mg/kg (carcinogenic doses) for 3 weeks. LncRNA differential expression showed a nonlinear dose response with none differentially expressed at 1.0 or 2.0mg/kg, 2 lncRNAs at 4.0mg/kg furan, and 83 at 8mg/kg, representing 13.3% (83/632) of the total number of differentially expressed transcripts. Among the lncRNAs observed, two lncRNAs examined showed transcriptional clustering with nearby protein-coding genes. LincRNA-p21 is an antisense transcript that is 15kb downstream from Cdkn1a locus and appears to be cotranscribed with the protein coding gene Cdkn1a at 8.0mg/kg furan. In a separate independent study, RNA samples from the livers of mice administered benzo(a)pyrene also demonstrated increased levels of Cdkn1a and the antisense lincRNA-p21 transcript. These data demonstrate that lncRNAs are transcriptional targets of furan exposures associated with levels of furan that are cytotoxic and induce cell proliferation. In addition, certain lncRNA transcripts are associated with the expression of nearby coding protein genes. We hypothesize that lncRNAs have potential as epigenetic biomarkers of carcinogenic exposures.

Biomarkers of exposure and effect in human lymphoblastoid TK6 cells following [13C2]-acetaldehyde exposure

The dose-response relationship for biomarkers of exposure (N(2)-ethylidene-dG adducts) and effect (cell survival and micronucleus formation) was determined across 4.5 orders of magnitude (50nM-2mM) using [(13)C2]-acetaldehyde exposures to human lymphoblastoid TK6 cells for 12h. There was a clear increase in exogenous N (2)-ethylidene-dG formation at exposure concentrations ≥ 1µM, whereas the endogenous adducts remained nearly constant across all exposure concentrations, with an average of 3.0 adducts/10(7) dG. Exogenous adducts were lower than endogenous adducts at concentrations ≤ 10µM and were greater than endogenous adducts at concentrations ≥ 250µM. When the endogenous and exogenous adducts were summed together, statistically significant increases in total adduct formation over the endogenous background occurred at 50µM. Cell survival and micronucleus formation were monitored across the exposure range and statistically significant decreases in cell survival and increases in micronucleus formation occurred at ≥ 1000µM. This research supports the hypothesis that endogenously produced reactive species, including acetaldehyde, are always present and constitute the majority of the observed biological effects following very low exposures to exogenous acetaldehyde. These data can replace default assumptions of linear extrapolation to very low doses of exogenous acetaldehyde for risk prediction.

Preparation of archival formalin-fixed paraffin-embedded mouse liver samples for use with the Agilent gene expression microarray platform

INTRODUCTION

Tissue samples are routinely formalin-fixed and paraffin-embedded (FFPE) for long term preservation. Gene expression analysis of archival FFPE tissues may advance knowledge of the molecular perturbations contributing to disease. However, formalin causes extensive degradation of RNA.

METHODS

We compared RNA quality/yield from FFPE samples using six commercial FFPE RNA extraction kits. In addition we compared four DNA microarray protocols for the Agilent 8×60K platform using 16year old FFPE mouse liver samples treated with phenobarbital or vehicle.

RESULTS

Despite low quality RNA, archival phenobarbital samples exhibited strong induction of the positive control genes Cyp2b9 and Cyp2b10 by quantitative real-time PCR (qPCR). We tested one- and two-color microarray designs and evaluated the effects of increasing the amount of hybridized cDNA. Canonical gene responders to phenobarbital were measurably induced under each experimental condition. Increasing the amount of labeled cDNA did not improve the overall signal intensity. One-color experiments yielded larger fold changes than two-color and the number of differentially expressed genes varied between protocols. Gene expression changes were validated by qPCR and literature searches. Individual protocols exhibited high rates of false positives; however, pathway analysis revealed that nine of the top ten canonical pathways were consistent across experiments. Genes that were differentially expressed in more than one experiment were more likely to be validated. Thus, we recommend that experiments on FFPE samples be done in duplicate to reduce false positives.

DISCUSSION

In this analysis of archival FFPE samples we were able to identify pathways that are consistent with phenobarbital's mechanism of action. Therefore, we conclude that FFPE samples can be used for meaningful microarray gene expression analyses.

Genotoxicity of styrene-acrylonitrile trimer in brain, liver, and blood cells of weanling F344 rats

Styrene-acrylonitrile Trimer (SAN Trimer), a by-product in production of acrylonitrile styrene plastics, was identified at a Superfund site in Dover Township, NJ, where childhood cancer incidence rates were elevated for a period of several years. SAN Trimer was therefore tested by the National Toxicology Program in a 2-year perinatal carcinogenicity study in F344/N rats and a bacterial mutagenicity assay; both studies gave negative results. To further characterize its genotoxicity, SAN Trimer was subsequently evaluated in a combined micronucleus (MN)/Comet assay in juvenile male and female F344 rats. SAN Trimer (37.5, 75, 150, or 300 mg/kg/day) was administered by gavage once daily for 4 days. Micronucleated reticulocyte (MN-RET) frequencies in blood were determined by flow cytometry, and DNA damage in blood, liver, and brain cells was assessed using the Comet assay. Highly significant dose-related increases (P < 0.0001) in MN-RET were measured in both male and female rats administered SAN Trimer. The RET population was reduced in high dose male rats, suggesting chemical-related bone marrow toxicity. Results of the Comet assay showed significant, dose-related increases in DNA damage in brain cells of male (P < 0.0074) and female (P < 0.0001) rats; increased levels of DNA damage were also measured in liver cells and leukocytes of treated rats. Chemical-related cytotoxicity was not indicated in any of the tissues examined for DNA damage. The results of this subacute MN/Comet assay indicate induction of significant genetic damage in multiple tissues of weanling F344 male and female rats after oral exposure to SAN Trimer.

Comparison of Comet assay dose-response for ethyl methanesulfonate using freshly prepared versus cryopreserved tissues

The National Toxicology Program (NTP) is using the Comet assay to evaluate genotoxic potential, and is investigating the integration of this assay into repeat-dose toxicity studies. To reduce sample-to-sample variability, address logistical concerns associated with evaluating multiple tissues from many animals, and accommodate sample collection at geographically distant testing facilities, tissue samples collected for Comet analysis by the NTP are routinely flash-frozen in liquid nitrogen and stored in a -80°C freezer until evaluation. To compare data obtained from frozen tissues to data from freshly isolated tissues, we conducted a dose-response study in male Sprague Dawley rats. Rats (5 per treatment group) were administered ethyl methanesulfonate (EMS; 0, 25, 50, 100, or 200 mg/kg) by gavage twice at an interval of 21 hr; blood, liver, stomach, and colon tissues were harvested 3 hr after the second treatment. Single-cell preparations from each of the four tissues were put into Hank's balanced salt solution with 10% fresh dimethyl sulfoxide. One aliquot of each tissue preparation was used for immediate analysis, while additional aliquots were flash-frozen in liquid nitrogen and stored in a -80°C freezer for 1 or 8 weeks. One set of 8-week frozen samples was shipped roundtrip via air courier from Research Triangle Park, NC to Rochester, NY prior to analysis. For all four tissues, results from frozen, nontransported samples showed a similar dose-response pattern for EMS-induced genotoxicity. We also demonstrated that for three tissues (blood, liver, stomach), air transport did not alter the sensitivity of the Comet assay for detecting DNA damage.

Dose-response assessment of naphthalene-induced genotoxicity and glutathione detoxication in human TK6 lymphoblasts

The dose-response relationship for the induction of micronuclei (MN) and the impact of glutathione (GSH) detoxication on naphthalene-induced cytotoxicity and genotoxicity were investigated in human TK6 cells. TK6 cells were exposed to 10 concentrations ranging from 0.0625 to 30μM naphthalene in the presence of β-naphthoflavone- and phenobarbital (βNP/PB)-induced rat liver S9 with a nicotinamide adenine dinucleotide phosphate-generating system. Three approaches were used to identify a no-observed-effect level (NOEL) for naphthalene-induced genotoxicity: (1) laboratory criteria of ≥ twofold increase over the concurrent solvent controls (NOEL = 10μM), (2) ANOVA with Bonferroni correction (NOEL = 2.5μM), and (3) the benchmark dose approach (BMCL(10) = 3.35μM). The NOEL and point of departure micronucleus frequency for naphthalene-induced MN are between the tested naphthalene concentrations of 2.5-10.0μM in this experimental system. Supplementation of the exposure system with physiological relevant concentrations of 5mM GSH eliminated naphthalene-induced cytotoxicity and genotoxicity; no increased cytotoxicity or genotoxicity was observed at concentrations of up to 500μM naphthalene in the presence of GSH compared with 2.5-10.0μM in the absence of GSH. Naphthalene bioactivation by βNP/PB-induced rat liver S9 exhibits a nonlinear dose-response for the induction of MN in TK6 cells with a NOEL of 2.5-10μM that in the presence of GSH is shifted upward greater than 50- to 200-fold. These data demonstrate a nonlinear dose-response for naphthalene-induced genotoxicity that is eliminated by GSH, and both observations should be considered when assessing human risk from naphthalene exposures.

No increases in biomarkers of genetic damage or pathological changes in heart and brain tissues in male rats administered methylphenidate hydrochloride (Ritalin) for 28 days

Following a 2005 report of chromosomal damage in children with attention deficit/hyperactivity disorder (ADHD) who were treated with the commonly prescribed medication methylphenidate (MPH), numerous studies have been conducted to clarify the risk for MPH-induced genetic damage. Although most of these studies reported no changes in genetic damage endpoints associated with exposure to MPH, one recent study (Andreazza et al. [2007]: Prog Neuropsychopharmacol Biol Psychiatry 31:1282-1288) reported an increase in DNA damage detected by the Comet assay in blood and brain cells of Wistar rats treated by intraperitoneal injection with 1, 2, or 10 mg/kg MPH; no increases in micronucleated lymphocyte frequencies were observed in these rats. To clarify these findings, we treated adult male Wistar Han rats with 0, 2, 10, or 25 mg/kg MPH by gavage once daily for 28 consecutive days and measured micronucleated reticulocyte (MN-RET) frequencies in blood, and DNA damage in blood, brain, and liver cells 4 hr after final dosing. Flow cytometric evaluation of blood revealed no significant increases in MN-RET. Comet assay evaluations of blood leukocytes and cells of the liver, as well as of the striatum, hippocampus, and frontal cortex of the brain showed no increases in DNA damage in MPH-treated rats in any of the three treatment groups. Thus, the previously reported observations of DNA damage in blood and brain tissue of rats exposed to MPH for 28 days were not confirmed in this study. Additionally, no histopathological changes in brain or heart, or elevated serum biomarkers of cardiac injury were observed in these MPH-exposed rats.