Scott C. Wesselkamper

Based on existing data and previous work, a series of studies is proposed as a basis toward a pragmatic early step in transforming toxicity testing. These studies were assembled into a data-driven framework that invokes successive tiers of testing with margin of exposure (MOE) as the primary metric. The first tier of the framework integrates data from high-throughput in vitro assays, in vitro-to-in vivo extrapolation (IVIVE) pharmacokinetic modeling, and exposure modeling. The in vitro assays are used to separate chemicals based on their relative selectivity in interacting with biological targets and identify the concentration at which these interactions occur. The IVIVE modeling converts in vitro concentrations into external dose for calculation of the point of departure (POD) and comparisons to human exposure estimates to yield a MOE. The second tier involves short-term in vivo studies, expanded pharmacokinetic evaluations, and refined human exposure estimates. The results from the second tier studies provide more accurate estimates of the POD and the MOE. The third tier contains the traditional animal studies currently used to assess chemical safety. In each tier, the POD for selective chemicals is based primarily on endpoints associated with a proposed mode of action, whereas the POD for nonselective chemicals is based on potential biological perturbation. Based on the MOE, a significant percentage of chemicals evaluated in the first 2 tiers could be eliminated from further testing. The framework provides a risk-based and animal-sparing approach to evaluate chemical safety, drawing broadly from previous experience but incorporating technological advances to increase efficiency.

BACKGROUND: The small sample sizes often used for microarray experiments result in poor estimates of variance if each gene is considered independently. Yet accurately estimating variability of gene expression measurements in microarray experiments is essential for correctly identifying differentially expressed genes. Several recently developed methods for testing differential expression of genes utilize hierarchical Bayesian models to "pool" information from multiple genes. We have developed a statistical testing procedure that further improves upon current methods by incorporating the well-documented relationship between the absolute gene expression level and the variance of gene expression measurements into the general empirical Bayes framework. RESULTS: We present a novel Bayesian moderated-T, which we show to perform favorably in simulations, with two real, dual-channel microarray experiments and in two controlled single-channel experiments. In simulations, the new method achieved greater power while correctly estimating the true proportion of false positives, and in the analysis of two publicly-available "spike-in" experiments, the new method performed favorably compared to all tested alternatives. We also applied our method to two experimental datasets and discuss the additional biological insights as revealed by our method in contrast to the others. The R-source code for implementing our algorithm is freely available at http://eh3.uc.edu/ibmt. CONCLUSION: We use a Bayesian hierarchical normal model to define a novel Intensity-Based Moderated T-statistic (IBMT). The method is completely data-dependent using empirical Bayes philosophy to estimate hyperparameters, and thus does not require specification of any free parameters. IBMT has the strength of balancing two important factors in the analysis of microarray data: the degree of independence of variances relative to the degree of identity (i.e. t-tests vs. equal variance assumption), and the relationship between variance and signal intensity. When this variance-intensity relationship is weak or does not exist, IBMT reduces to a previously described moderated t-statistic. Furthermore, our method may be directly applied to any array platform and experimental design. Together, these properties show IBMT to be a valuable option in the analysis of virtually any microarray experiment.

The number of legacy chemicals without toxicity reference values combined with the rate of new chemical development is overwhelming the capacity of the traditional risk assessment paradigm. More efficient approaches are needed to quantitatively estimate chemical risks. In this study, rats were dosed orally with multiple doses of six chemicals for 5 days and 2, 4, and 13 weeks. Target organs were analyzed for traditional histological and organ weight changes and transcriptional changes using microarrays. Histological and organ weight changes in this study and the tumor incidences in the original cancer bioassays were analyzed using benchmark dose (BMD) methods to identify noncancer and cancer points of departure. The dose-response changes in gene expression were also analyzed using BMD methods and the responses grouped based on signaling pathways. A comparison of transcriptional BMD values for the most sensitive pathway with BMD values for the noncancer and cancer apical endpoints showed a high degree of correlation at all time points. When the analysis included data from an earlier study with eight additional chemicals, transcriptional BMD values for the most sensitive pathway were significantly correlated with noncancer (r = 0.827, p = 0.0031) and cancer-related (r = 0.940, p = 0.0002) BMD values at 13 weeks. The average ratio of apical-to-transcriptional BMD values was less than two, suggesting that for the current chemicals, transcriptional perturbation did not occur at significantly lower doses than apical responses. Based on our results, we propose a practical framework for application of transcriptomic data to chemical risk assessment.

The traditional approach for estimating noncancer and cancer reference values in quantitative chemical risk assessment is time and resource intensive. The extent and nature of the studies required under the traditional approach has limited the number of chemicals with published risk assessments. In this study, female mice were exposed for 13 weeks to multiple concentrations of five chemicals that were positive in a 2-year cancer bioassay. Traditional histological and organ weight changes were evaluated, and gene expression microarray analysis was performed on the target tissues. The histological, organ weight changes, and the original tumor incidences in the original cancer bioassay were analyzed using standard benchmark dose (BMD) methods to identify noncancer and cancer points of departure, respectively. The dose-related changes in gene expression were also analyzed using a BMD approach and the responses grouped based on cellular biological processes. A comparison of the transcriptional BMD values with those for the traditional noncancer and cancer apical endpoints showed a high degree of correlation for specific cellular biological processes. For chemicals with human exposure data, the transcriptional BMD values were also used to calculate a margin of exposure. The margins of exposure ranged from 1900 to 54,000. Both the correlation between the BMD values for the transcriptional and apical endpoints and the margin of exposure analysis suggest that transcriptional BMD values may be used as potential points of departure for noncancer and cancer risk assessment.