A method to integrate benchmark dose estimates with genomic data to assess the functional effects of chemical exposure

The use of genomic technology for assessing health risks associated with chemical exposure has significant potential, but its direct application has proven to be challenging for the toxicology and risk assessment communities. In this study, a method was established for analyzing dose-response microarray data using benchmark dose (BMD) calculations and gene ontology (GO) classification. Gene expression changes in the rat nasal epithelium following acute formaldehyde exposure were used as a case study. The gene expression data were first analyzed using a one-way ANOVA to identify genes that showed significant dose-response behavior. These genes were then fit to a series of four statistical models (linear, second-degree polynomial, third-degree polynomial, and power models) and the least complex model that best described the data was selected. The genes were matched to their associated GO categories, and the average BMD and benchmark dose lower confidence limit (BMDL) were calculated for each GO category. The results were used to identify doses at which individual cellular processes were altered. For the formaldehyde exposures, the BMD estimates for the GO categories related to cell proliferation and DNA damage were similar to those measured in previous studies using cell labeling indices and DNA-protein cross-links and consistent with the BMD estimated for rat nasal tumors. The method represents a significant advance in applying genomic information to risk assessment by allowing a comprehensive survey of molecular changes associated with chemical exposure and providing the capability to identify reference doses at which particular cellular processes are altered.