Case study: Targeted RNA-sequencing of aged formalin-fixed paraffin-embedded samples for understanding chemical mode of action.

Formalin-fixed paraffin-embedded (FFPE) samples are the only remaining biological archive for many toxicological and clinical studies, yet their use in genomics has been limited due to nucleic acid damage from formalin fixation. Older FFPE samples with highly degraded RNA pose a particularly difficult technical challenge. Probe-based targeted sequencing technologies show promise in addressing this issue but have not been directly compared to standard whole-genome RNA-Sequencing (RNA-Seq) methods. In this study, we evaluated dose-dependent transcriptional changes from paired frozen (FROZ) and FFPE liver samples stored for over 20 years using targeted resequencing (TempO-Seq) and whole-genome RNA-Seq methods. Samples were originally collected from male mice exposed to a reference chemical (dichloroacetic acid, DCA) at 0, 198, 313, and 427 mg/kg-day (n = 6/dose) by drinking water for 6 days. TempO-Seq showed high overlap in differentially expressed genes (DEGs) between matched FFPE and FROZ samples and high concordance in fold-change values across the two highest dose levels of DCA vs. control (R2 ≥ 0.94). Similarly, high concordance in fold-change values was observed between TempO-Seq FFPE and RNA-Seq FROZ results (R2 ≥ 0.92). In contrast, RNA-Seq FFPE samples showed few overlapping DEGs compared to FROZ RNA-Seq (≤5 for all dose groups). Modeling of DCA-dependent changes in gene sets identified benchmark doses from TempO-Seq FROZ and FFPE samples within 1.4-fold of RNA-Seq FROZ samples (93.9 mg/kg-d), whereas RNA-Seq FFPE samples were 3.3-fold higher (310.3 mg/kg-d). This work demonstrates that targeted sequencing may provide a more robust method for quantifying gene expression profiles from aged archival FFPE samples.

Organocatalyst treatment improves variant calling and mutant detection in archival clinical samples.

Formalin fixation of biological specimens damages nucleic acids and limits their use in genomic analyses. Previously, we showed that RNA isolation with an organocatalyst (2-amino-5-methylphenyl phosphonic acid, used to speed up reversal of formalin-induced adducts) and extended heated incubation (ORGΔ) improved RNA-sequencing data from formalin-fixed paraffin-embedded (FFPE) tissue samples. The primary goal of this study was to evaluate whether ORGΔ treatment improves DNA-sequencing data from clinical FFPE samples. We isolated RNA and DNA ± ORGΔ from paired FFPE and frozen human renal and ovarian carcinoma specimens collected as part of the National Cancer Institute Biospecimen Pre-analytical Variables program. Tumor types were microscopically confirmed from adjacent tissue sections. Following extraction, DNA was fragmented and sequenced and differences were compared between frozen and FFPE sample pairs. Treatment with ORGΔ improved concurrent SNP calls in FFPE DNA compared to non-ORGΔ FFPE samples and enhanced confidence in SNP calls for all FFPE DNA samples, beyond that of matched frozen samples. In general, the concordant SNPs identified in paired frozen and FFPE DNA samples agreed for both genotype and homozygosity vs. heterozygosity of calls regardless of ORGΔ treatment. The increased confidence in ORGΔ FFPE DNA variant calls relative to the matched frozen DNA suggests a novel application of this method. With further optimization, this method may improve quality of DNA-sequencing data in FFPE as well as frozen tissue samples.

Genomic comparisons between hepatocarcinogenic and non-hepatocarcinogenic organophosphate insecticides in the mouse liver.

Short-term biomarkers of toxicity have an increasingly important role in the screening and prioritization of new chemicals. In this study, we examined early indicators of liver toxicity for three reference organophosphate (OP) chemicals, which are among the most widely used insecticides in the world. The OP methidathion was previously shown to increase the incidence of liver toxicity, including hepatocellular tumors, in male mice. To provide insights into the adverse outcome pathway (AOP) that underlies these tumors, effects of methidathion in the male mouse liver were examined after 7 and 28 day exposures and compared to those of two other OPs that either do not increase (fenthion) or possibly suppress liver cancer (parathion) in mice. None of the chemicals caused increases in liver weight/body weight or histopathological changes in the liver. Parathion decreased liver cell proliferation after 7 and 28 days while the other chemicals had no effects. There was no evidence for hepatotoxicity in any of the treatment groups. Full-genome microarray analysis of the livers from the 7 and 28 day treatments demonstrated that methidathion and fenthion regulated a large number of overlapping genes, while parathion regulated a unique set of genes. Examination of cytochrome P450 enzyme activities and use of predictive gene expression biomarkers found no consistent evidence for activation of AhR, CAR, PXR, or PPARα. Parathion suppressed the male-specific gene expression pattern through STAT5b, similar to genetic and dietary conditions that decrease liver tumor incidence in mice. Overall, these findings indicate that methidathion causes liver cancer by a mechanism that does not involve common mechanisms of liver cancer induction.

Direct formalin fixation induces widespread transcriptomic effects in archival tissue samples.

Sequencing technologies now provide unprecedented access to genomic information in archival formalin-fixed paraffin-embedded (FFPE) tissue samples. However, little is known about artifacts induced during formalin fixation, which could bias results. Here we evaluated global changes in RNA-sequencing profiles between matched frozen and FFPE samples. RNA-sequencing was performed on liver samples collected from mice treated with a reference chemical (phenobarbital) or vehicle control for 7 days. Each sample was divided into four parts: (1) fresh-frozen, (2) direct-fixed in formalin for 18 h, (3) frozen then formalin-fixed, and (4) frozen then ethanol-fixed and paraffin-embedded (n = 6/group/condition). Direct fixation resulted in 2,946 differentially expressed genes (DEGs) vs. fresh-frozen, 98% of which were down-regulated. Freezing prior to formalin fixation had ≥ 95% fewer DEGs vs. direct fixation, indicating that most formalin-derived transcriptional effects in the liver occurred during fixation. This finding was supported by retrospective studies of paired frozen and FFPE samples, which identified consistent enrichment in oxidative stress, mitochondrial dysfunction, and transcription initiation pathways with direct fixation. Notably, direct formalin fixation in the parent study did not significantly impact response profiles resulting from chemical exposure. These results advance our understanding of FFPE samples as a resource for genomic research.

Enhanced Quality Metrics for Assessing RNA Derived From Archival Formalin-Fixed Paraffin-Embedded Tissue Samples.

Formalin-fixed paraffin-embedded (FFPE) tissues provide an important resource for toxicogenomic research. However, variability in the integrity or quality of RNA obtained from archival FFPE specimens can lead to unreliable data and wasted resources, and standard protocols for measuring RNA integrity do not adequately assess the suitability of FFPE RNA. The main goal of this study was to identify improved methods for evaluating FFPE RNA quality for whole-genome sequencing. We examined RNA quality metrics conducted prior to RNA-sequencing in paired frozen and FFPE samples with varying levels of quality based on age in block and time in formalin. RNA quality was measured by the RNA integrity number (RIN), a modified RIN called the paraffin-embedded RNA metric, the percentage of RNA fragments >100-300 nucleotides in size (DV100-300), and 2 quantitative PCR-based methods. This information was correlated to sequencing read quality, mapping, and gene detection. Among fragmentation-based methods, DV and PCR-based metrics were more informative than RIN or paraffin-embedded RNA metric in determining sequencing success. Across low- and high-quality FFPE samples, a minimum of 80% of RNA fragments >100 nucleotides (DV100 > 80) provided the best indication of gene diversity and read counts upon sequencing. The PCR-based methods further showed quantitative reductions in amplifiable RNA of target genes related to sample age and time in formalin that inform input quantity of FFPE RNA for sequencing. These results should aid in screening and prioritizing archival FFPE samples for retrospective analyses of gene expression.

Authentic Inquiry through Modeling in Biology (AIM-Bio): An Introductory Laboratory Curriculum That Increases Undergraduates' Scientific Agency and Skills.

Providing opportunities for science, technology, engineering, and mathematics undergraduates to engage in authentic scientific practices is likely to influence their view of science and may impact their decision to persist through graduation. Laboratory courses provide a natural place to introduce students to scientific practices, but existing curricula often miss this opportunity by focusing on confirming science content rather than exploring authentic questions. Integrating authentic science within laboratory courses is particularly challenging at high-enrollment institutions and community colleges, where access to research-active faculty may be limiting. The Authentic Inquiry through Modeling in Biology (AIM-Bio) curriculum presented here engages students in authentic scientific practices through iterative cycles of model generation, testing, and revision. AIM-Bio university and community college students demonstrated their ability to propose diverse models for biological phenomena, formulate and address hypotheses by designing and conducting experiments, and collaborate with classmates to revise models based on experimental data. Assessments demonstrated that AIM-Bio students had an enhanced sense of project ownership and greater identification as scientists compared with students in existing laboratory courses. AIM-Bio students also experienced measurable gains in their nature of science understanding and skills for doing science. Our results suggest AIM-Bio as a potential alternative to more resource-intensive curricula with similar outcomes.

From the Cover: Genomic Effects of Androstenedione and Sex-Specific Liver Cancer Susceptibility in Mice.

Current strategies for predicting carcinogenic mode of action for nongenotoxic chemicals are based on identification of early key events in toxicity pathways. The goal of this study was to evaluate short-term key event indicators resulting from exposure to androstenedione (A4), an androgen receptor agonist and known liver carcinogen in mice. Liver cancer is more prevalent in men compared with women, but androgen-related pathways underlying this sex difference have not been clearly identified. Short-term hepatic effects of A4 were compared with reference agonists of the estrogen receptor (ethinyl estradiol, EE) and glucocorticoid receptor (prednisone, PRED). Male B6C3F1 mice were exposed for 7 or 28 days to A4, EE, or PRED. EE increased and PRED suppressed hepatocyte proliferation, while A4 had no detectable effects. In a microarray analysis, EE and PRED altered >3000 and >670 genes, respectively, in a dose-dependent manner, whereas A4 did not significantly alter any genes. Gene expression was subsequently examined in archival liver samples from male and female B6C3F1 mice exposed to A4 for 90 days. A4 altered more genes in females than males and did not alter expression of genes linked to activation of the mitogenic xenobiotic receptors AhR, CAR, and PPARα in either sex. A gene expression biomarker was used to show that in female mice, the high dose of A4 activated the growth hormone-regulated transcription factor STAT5b, which controls sexually dimorphic gene expression in the liver. These findings suggest that A4 induces subtle age-related effects on STAT5b signaling that may contribute to the higher risk of liver cancer in males compared with females.

Metabolic Disruption Early in Life is Associated With Latent Carcinogenic Activity of Dichloroacetic Acid in Mice.

Early-life environmental factors can influence later-life susceptibility to cancer. Recent evidence suggests that metabolic pathways may mediate this type of latency effect. Previously, we reported that short-term exposure to dichloroacetic acid (DCA) increased liver cancer in mice 84 weeks after exposure was stopped. Here, we evaluated time course dynamics for key events related to this effect. This study followed a stop-exposure design in which 28-day-old male B6C3F1 mice were given the following treatments in drinking water for up to 93 weeks: deionized water (dH2O, control); 3.5 g/l DCA continuously; or 3.5 g/l DCA for 4-52 weeks followed by dH2O. Effects were evaluated at eight interim time points. A short-term biomarker study was used to evaluate DCA effects at 6, 15, and 30 days. Liver tumor incidence was higher in all DCA treatment groups, including carcinomas in 82% of mice previously treated with DCA for only 4 weeks. Direct effects of DCA in the short-term study included decreased liver cell proliferation and marked mRNA changes related to mitochondrial dysfunction and altered cell metabolism. However, all observed short-term effects of DCA were ultimately reversible, and prior DCA treatment did not affect liver cell proliferation, apoptosis, necrosis, or DNA sequence variants with age. Key intermediate events resulting from transient DCA exposure do not fit classical cytotoxic, mitogenic, or genotoxic modes of action for carcinogenesis, suggesting a distinct mechanism associated with early-life metabolic disruption.

Dynamics of growth zone patterning in the milkweed bug Oncopeltus fasciatus.

We describe the dynamic process of abdominal segment generation in the milkweed bug Oncopeltus fasciatus We present detailed morphological measurements of the growing germband throughout segmentation. Our data are complemented by cell division profiles and expression patterns of key genes, including invected and even-skipped as markers for different stages of segment formation. We describe morphological and mechanistic changes in the growth zone and in nascent segments during the generation of individual segments and throughout segmentation, and examine the relative contribution of newly formed versus existing tissue to segment formation. Although abdominal segment addition is primarily generated through the rearrangement of a pool of undifferentiated cells, there is nonetheless proliferation in the posterior. By correlating proliferation with gene expression in the growth zone, we propose a model for growth zone dynamics during segmentation in which the growth zone is functionally subdivided into two distinct regions: a posterior region devoted to a slow rate of growth among undifferentiated cells, and an anterior region in which segmental differentiation is initiated and proliferation inhibited.

Editor's Highlight: Genetic Targets of Acute Toluene Inhalation in Drosophila melanogaster.

Interpretation and use of data from high-throughput assays for chemical toxicity require links between effects at molecular targets and adverse outcomes in whole animals. The well-characterized genome of Drosophila melanogaster provides a potential model system by which phenotypic responses to chemicals can be mapped to genes associated with those responses, which may in turn suggest adverse outcome pathways associated with those genes. To determine the utility of this approach, we used the Drosophila Genetics Reference Panel (DGRP), a collection of ∼200 homozygous lines of fruit flies whose genomes have been sequenced. We quantified toluene-induced suppression of motor activity in 123 lines of these flies during exposure to toluene, a volatile organic compound known to induce narcosis in mammals via its effects on neuronal ion channels. We then applied genome-wide association analyses on this effect of toluene using the DGRP web portal (http://dgrp2.gnets.ncsu.edu), which identified polymorphisms in candidate genes associated with the variation in response to toluene exposure. We tested ∼2 million variants and found 82 polymorphisms located in or near 66 candidate genes that were associated with phenotypic variation for sensitivity to toluene at P < 5 × 10-5, and human orthologs for 52 of these candidate Drosophila genes. None of these orthologs are known to be involved in canonical pathways for mammalian neuronal ion channels, including GABA, glutamate, dopamine, glycine, serotonin, and voltage sensitive calcium channels. Thus this analysis did not reveal a genetic signature consistent with processes previously shown to be involved in toluene-induced narcosis in mammals. The list of the human orthologs included Gene Ontology terms associated with signaling, nervous system development and embryonic morphogenesis; these orthologs may provide insight into potential new pathways that could mediate the narcotic effects of toluene.

Editor's Highlight: Dose-Response Analysis of RNA-Seq Profiles in Archival Formalin-Fixed Paraffin-Embedded Samples.

Use of archival resources has been limited to date by inconsistent methods for genomic profiling of degraded RNA from formalin-fixed paraffin-embedded (FFPE) samples. RNA-sequencing offers a promising way to address this problem. Here, we evaluated transcriptomic dose responses using RNA-sequencing in paired FFPE and frozen (FROZ) samples from 2 archival studies in mice, one <2 years old and the other >20 years old. Experimental treatments included 3 different doses of di(2-ethylhexyl)phthalate or dichloroacetic acid for the recently archived and older studies, respectively. Total RNA was ribo-depleted and sequenced using the Illumina HiSeq platform. In the recently archived study, FFPE samples had 35% lower total counts compared to FROZ samples but high concordance in fold-change values of differentially expressed genes (DEGs) (r2 = 0.99), highly enriched pathways (90% overlap with FROZ), and benchmark dose estimates for preselected target genes (<5% difference vs FROZ). In contrast, older FFPE samples had markedly lower total counts (3% of FROZ) and poor concordance in global DEGs and pathways. However, counts from FFPE and FROZ samples still positively correlated (r2 = 0.84 across all transcripts) and showed comparable dose responses for more highly expressed target genes. These findings highlight potential applications and issues in using RNA-sequencing data from FFPE samples. Recently archived FFPE samples were highly similar to FROZ samples in sequencing quality metrics, DEG profiles, and dose-response parameters, while further methods development is needed for older lower-quality FFPE samples. This work should help advance the use of archival resources in chemical safety and translational science.

Dose and Effect Thresholds for Early Key Events in a PPARα-Mediated Mode of Action.

Current strategies for predicting adverse health outcomes of environmental chemicals are centered on early key events in toxicity pathways. However, quantitative relationships between early molecular changes in a given pathway and later health effects are often poorly defined. The goal of this study was to evaluate short-term key event indicators using qualitative and quantitative methods in an established pathway of mouse liver tumorigenesis mediated by peroxisome proliferator-activated receptor alpha (PPARα). Male B6C3F1 mice were exposed for 7 days to di (2-ethylhexyl) phthalate (DEHP), di-n-octyl phthalate (DNOP), and n-butyl benzyl phthalate (BBP), which vary in PPARα activity and liver tumorigenicity. Each phthalate increased expression of select PPARα target genes at 7 days, while only DEHP significantly increased liver cell proliferation labeling index (LI). Transcriptional benchmark dose (BMDT) estimates for dose-related genomic markers stratified phthalates according to hypothetical tumorigenic potencies, unlike BMDs for non-genomic endpoints (relative liver weights or proliferation). The 7-day BMDT values for Acot1 as a surrogate measure for PPARα activation were 29, 370, and 676 mg/kg/day for DEHP, DNOP, and BBP, respectively, distinguishing DEHP (liver tumor BMD of 35 mg/kg/day) from non-tumorigenic DNOP and BBP. Effect thresholds were generated using linear regression of DEHP effects at 7 days and 2-year tumor incidence values to anchor early response molecular indicators and a later phenotypic outcome. Thresholds varied widely by marker, from 2-fold (Pdk4 and proliferation LI) to 30-fold (Acot1) induction to reach hypothetical tumorigenic expression levels. These findings highlight key issues in defining thresholds for biological adversity based on molecular changes.

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.

Latent carcinogenicity of early-life exposure to dichloroacetic acid in mice.

Environmental exposures occurring early in life may have an important influence on cancer risk later in life. Here, we investigated carryover effects of dichloroacetic acid (DCA), a small molecule analog of pyruvate with metabolic programming properties, on age-related incidence of liver cancer. The study followed a stop-exposure/promotion design in which 4-week-old male and female B6C3F1 mice received the following treatments: deionized water alone (dH2O, control); dH2O with 0.06% phenobarbital (PB), a mouse liver tumor promoter; or DCA (1.0, 2.0 or 3.5g/l) for 10 weeks followed by dH2O or PB (n = 20-30/group/sex). Pathology and molecular assessments were performed at 98 weeks of age. In the absence of PB, early-life exposure to DCA increased the incidence and number of hepatocellular tumors in male and female mice compared with controls. Significant dose trends were observed in both sexes. At the high dose level, 10 weeks of prior DCA treatment induced comparable effects (≥85% tumor incidence and number) to those seen after continuous lifetime exposure. Prior DCA treatment did not enhance or inhibit the carcinogenic effects of PB, induce persistent liver cytotoxicity or preneoplastic changes on histopathology or alter DNA sequence variant profiles within liver tumors compared with controls. Distinct changes in liver messenger RNA and micro RNA profiles associated with prior DCA treatment were not apparent at 98 weeks. Our findings demonstrate that early-life exposure to DCA may be as carcinogenic as life-long exposures, potentially via epigenetic-mediated effects related to cellular metabolism.

Dose-response modeling of early molecular and cellular key events in the CAR-mediated hepatocarcinogenesis pathway.

Low-dose extrapolation and dose-related transitions are paramount in the ongoing debate regarding the quantification of cancer risks for nongenotoxic carcinogens. Phenobarbital (PB) is a prototypical nongenotoxic carcinogen that activates the constitutive androstane receptor (CAR) resulting in rodent liver tumors. In this study, male and female CD-1 mice administered dietary PB at 0, 0.15, 1.5, 15, 75, or 150 mg/kg-day for 2 or 7 days to characterize multiple apical and molecular endpoints below, at (∼75 mg/kg-day), and above the carcinogenic dose level of PB and examine these responses using benchmark dose modeling. Linear toxicokinetics were observed for all doses. Increased liver weight, hepatocellular hypertrophy, and mitotic figures were seen at 75 and 150 mg/kg-day. CAR activation, based on Cyp2b qPCR and pentoxyresorufin dealkylase activity, occurred at doses ≥ 1.5 mg/kg-day. The no-observable transcriptional effect level for global gene expression was 15 mg/kg-day. At 2 days, several xenobiotic metabolism and cell protective pathways were activated at lower doses and to a greater degree in females. However, hepatocellular proliferation, quantified by bromodeoxyuridine immunohistochemistry, was the most sensitive indicator of PB exposure with female mice more sensitive than males, contrary to sex-specific differences in sensitivity to hepatocarcinogenesis. Taken together, the identification of low-dose cellular and molecular transitions in the subtumorigenic dose range aids the understanding of early key events in CAR-mediated hepatocarcinogenesis.

Concordance of transcriptional and apical benchmark dose levels for conazole-induced liver effects in mice.

The ability to anchor chemical class-based gene expression changes to phenotypic lesions and to describe these changes as a function of dose and time informs mode-of-action determinations and improves quantitative risk assessments. Previous global expression profiling identified a 330-probe cluster differentially expressed and commonly responsive to 3 hepatotumorigenic conazoles (cyproconazole, epoxiconazole, and propiconazole) at 30 days. Extended to 2 more conazoles (triadimefon and myclobutanil), the present assessment encompasses 4 tumorigenic and 1 nontumorigenic conazole. Transcriptional benchmark dose levels (BMDL(T)) were estimated for a subset of the cluster with dose-responsive behavior and a ≥ 5-fold increase or decrease in signal intensity at the highest dose. These genes primarily encompassed CAR/RXR activation, P450 metabolism, liver hypertrophy- glutathione depletion, LPS/IL-1-mediated inhibition of RXR, and NRF2-mediated oxidative stress pathways. Median BMDL(T) estimates from the subset were concordant (within a factor of 2.4) with apical benchmark doses (BMDL(A)) for increased liver weight at 30 days for the 5 conazoles. The 30-day median BMDL(T) estimates were within one-half order of magnitude of the chronic BMDLA for hepatocellular tumors. Potency differences seen in the dose-responsive transcription of certain phase II metabolism, bile acid detoxification, and lipid oxidation genes mirrored each conazole's tumorigenic potency. The 30-day BMDL(T) corresponded to tumorigenic potency on a milligram per kilogram day basis with cyproconazole > epoxiconazole > propiconazole > triadimefon > myclobutanil (nontumorigenic). These results support the utility of measuring short-term gene expression changes to inform quantitative risk assessments from long-term exposures.

Transcriptional responses in rat brain associated with sub-chronic toluene inhalation are not predicted by effects of acute toluene inhalation.

A primary public health concern regarding environmental chemicals is the potential for persistent effects from long-term exposure, and approaches to estimate these effects from short-term exposures are needed. Toluene, a ubiquitous air pollutant, exerts well-documented acute and persistent CNS-mediated effects from a variety of exposure scenarios, and so provides a useful case for determining whether its persistent effects can be predicted from its acute effects on the CNS. We recently reported that acute inhalation of toluene produced transcriptional effects in rat brain 18 h following a single, acute 6-h exposure to toluene. The goal of the present study was to determine whether these acute effects are also evident after long-term (sub-chronic) exposure to toluene, and thereby provide a mechanistic basis for predicting its persistent effects from short-term exposures. Male Long-Evans rats were exposed to toluene via inhalation (0, 10, 100, 1000 ppm, n=5/dose), 6h/day for 64 days, excluding weekends. The day following the final exposure, total mRNA was extracted from the cerebral cortex and striatum, and gene expression evaluated using Affymetrix arrays. Principal component analysis using all samples showed a clear discrimination of tissues, with striatum having more within-group variance than cortex. Differentially-expressed genes (DEGs) whose expression was altered by toluene were identified in each tissue by ANOVA followed by mapping to pathways. Analysis of striatum revealed 22, 57, and 94 significant DEGs for the 10 ppm, 100 ppm, and 1000 ppm doses, respectively, far fewer than the 3352 DEGS previously observed after acute exposure. In addition, the direction of change in the 57 DEGs common to both exposures differed between acute and sub-chronic exposure scenarios. Thus, relative to acute toluene exposure, sub-chronic exposure yielded both quantitative and qualitative differences in transcriptional response. Based on the current data, long-term gene expression changes after toluene inhalation cannot be readily predicted by acute responses.

A multi-cell, multi-scale model of vertebrate segmentation and somite formation.

Somitogenesis, the formation of the body's primary segmental structure common to all vertebrate development, requires coordination between biological mechanisms at several scales. Explaining how these mechanisms interact across scales and how events are coordinated in space and time is necessary for a complete understanding of somitogenesis and its evolutionary flexibility. So far, mechanisms of somitogenesis have been studied independently. To test the consistency, integrability and combined explanatory power of current prevailing hypotheses, we built an integrated clock-and-wavefront model including submodels of the intracellular segmentation clock, intercellular segmentation-clock coupling via Delta/Notch signaling, an FGF8 determination front, delayed differentiation, clock-wavefront readout, and differential-cell-cell-adhesion-driven cell sorting. We identify inconsistencies between existing submodels and gaps in the current understanding of somitogenesis mechanisms, and propose novel submodels and extensions of existing submodels where necessary. For reasonable initial conditions, 2D simulations of our model robustly generate spatially and temporally regular somites, realistic dynamic morphologies and spontaneous emergence of anterior-traveling stripes of Lfng. We show that these traveling stripes are pseudo-waves rather than true propagating waves. Our model is flexible enough to generate interspecies-like variation in somite size in response to changes in the PSM growth rate and segmentation-clock period, and in the number and width of Lfng stripes in response to changes in the PSM growth rate, segmentation-clock period and PSM length.

Acute toluene exposure alters expression of genes in the central nervous system associated with synaptic structure and function.

Toluene is a volatile organic compound (VOC) and a ubiquitous air pollutant of interest to EPA regulatory programs. Whereas its acute functional effects are well described, several modes of action in the CNS have been proposed. Therefore, we sought to identify potential pathways mediating direct or indirect effects of VOCs by investigating the genomic response of the rat CNS to acutely-inhaled toluene. Adult male Long-Evans rats inhaled clean air or 1000 ppm toluene vapor for 6 h. Specific brain regions were collected from the rats either immediately after 6 h of treatment or 18 h after removal from the exposure chambers (n=6/group/time). Total mRNA was extracted from the striatum and hybridized to Rat 230A Affymetrix arrays. Statistical analyses showed 226 and 3352 transcripts altered in the toluene-exposed groups relative to controls at the 6 h time point and after the 18 h recovery period, respectively. Relative to controls, toluene exposure was associated with induction or repression of genes in pathways associated with synaptic plasticity, including long-term depression, GABA receptor signaling and mitochondrial function. In each of these pathways, responses were characterized by changes in a small number of transcripts following the 6 h toluene inhalation and with substantial increases in numbers of changed transcripts at 18 h recovery following termination of exposure. This report provides the first global genomic evidence that CNS pathways affected by toluene are strongly associated with neurological processes participating in synaptic transmission and plasticity.

Transcriptional profile of diuron-induced toxicity on the urinary bladder of male Wistar rats to inform mode of action.

Diuron (3-(3,4-dichlorophenyl)-1,1-dimethylurea) is a substituted urea herbicide that induces rat urinary bladder urothelial tumors at high dietary levels (2500 ppm). The specific mode of action and molecular alterations triggered by diuron, however, have not been clarified. The present study evaluated the dose-dependent effects of mucosal alterations and transcriptional changes in the urinary bladder of rats exposed to diuron. Six-week-old male Wistar rats were treated with 0, 60, 125, 1250, and 2500 ppm of diuron in the diet for 20 weeks. Histologic examination showed urothelial hyperplasia present in rats treated with either 1250 or 2500 ppm of diuron but not 60 or 125 ppm. Comprehensive gene expression analyses of urothelial cell RNA were conducted using Affymetrix microarrays. The numbers of differentially expressed transcripts between each treatment group and control increased with diuron dose. Based on similar histology and gene expression responses, the treatment groups were regrouped into a high-dose (1250 and 2500 ppm) and low-dose group (60 and 125 ppm). These data suggest that persistent exposure to high dietary concentrations of diuron induces oxidative stress, increases cellular metabolism, and enhances cell death that is associated with sustained urothelial hyperplasia.