Use of comparative genomics approaches to characterize interspecies differences in response to environmental chemicals: challenges, opportunities, and research needs.

A critical challenge for environmental chemical risk assessment is the characterization and reduction of uncertainties introduced when extrapolating inferences from one species to another. The purpose of this article is to explore the challenges, opportunities, and research needs surrounding the issue of how genomics data and computational and systems level approaches can be applied to inform differences in response to environmental chemical exposure across species. We propose that the data, tools, and evolutionary framework of comparative genomics be adapted to inform interspecies differences in chemical mechanisms of action. We compare and contrast existing approaches, from disciplines as varied as evolutionary biology, systems biology, mathematics, and computer science, that can be used, modified, and combined in new ways to discover and characterize interspecies differences in chemical mechanism of action which, in turn, can be explored for application to risk assessment. We consider how genetic, protein, pathway, and network information can be interrogated from an evolutionary biology perspective to effectively characterize variations in biological processes of toxicological relevance among organisms. We conclude that comparative genomics approaches show promise for characterizing interspecies differences in mechanisms of action, and further, for improving our understanding of the uncertainties inherent in extrapolating inferences across species in both ecological and human health risk assessment. To achieve long-term relevance and consistent use in environmental chemical risk assessment, improved bioinformatics tools, computational methods robust to data gaps, and quantitative approaches for conducting extrapolations across species are critically needed. Specific areas ripe for research to address these needs are recommended.

An experimental assessment of in silico haplotype association mapping in laboratory mice.

BACKGROUND: To assess the utility of haplotype association mapping (HAM) as a quantitative trait locus (QTL) discovery tool, we conducted HAM analyses for red blood cell count (RBC) and high density lipoprotein cholesterol (HDL) in mice. We then experimentally tested each HAM QTL using published crosses or new F2 intercrosses guided by the haplotype at the HAM peaks. RESULTS: The HAM for RBC, using 33 classic inbred lines, revealed 8 QTLs; 2 of these were true positives as shown by published crosses. A HAM-guided (C57BL/6J x CBA/J)F2 intercross we carried out verified 2 more as true positives and 4 as false positives. The HAM for HDL, using 81 strains including recombinant inbred lines and chromosome substitution strains, detected 46 QTLs. Of these, 36 were true positives as shown by published crosses. A HAM-guided (C57BL/6J x A/J)F2 intercross that we carried out verified 2 more as true positives and 8 as false positives. By testing each HAM QTL for RBC and HDL, we demonstrated that 78% of the 54 HAM peaks were true positives and 22% were false positives. Interestingly, all false positives were in significant allelic association with one or more real QTL. CONCLUSION: Because type I errors (false positives) can be detected experimentally, we conclude that HAM is useful for QTL detection and narrowing. We advocate the powerful and economical combined approach demonstrated here: the use of HAM for QTL discovery, followed by mitigation of the false positive problem by testing the HAM-predicted QTLs with small HAM-guided experimental crosses.

Effects of atherogenic diet on hepatic gene expression across mouse strains.

Diets high in fat and cholesterol are associated with increased obesity and metabolic disease in mice and humans. To study the molecular basis of the metabolic response to dietary fat, 10 inbred strains of mice were fed atherogenic high-fat and control low-fat diets. Liver gene expression and whole animal phenotypes were measured and analyzed in both sexes. The effects of diet, strain, and sex on gene expression were determined irrespective of complex processes, such as feedback mechanisms, that could have mediated the genomic responses. Global gene expression analyses demonstrated that animals of the same strain and sex have similar transcriptional profiles on a low-fat diet, but strains may show considerable variability in response to high-fat diet. Functional profiling indicated that high-fat feeding induced genes in the immune response, indicating liver damage, and repressed cholesterol biosynthesis. The physiological significance of the transcriptional changes was confirmed by a correlation analysis of transcript levels with whole animal phenotypes. The results found here were used to confirm a previously identified quantitative trait locus on chromosome 17 identified in males fed a high-fat diet in two crosses, PERA x DBA/2 and PERA x I/Ln. The gene expression data and phenotype data have been made publicly available as an online tool for exploring the effects of atherogenic diet in inbred mouse strains (http://cgd-array.jax.org/DietStrainSurvey).

Practical applications of the bioinformatics toolbox for narrowing quantitative trait loci.

Dissecting the genes involved in complex traits can be confounded by multiple factors, including extensive epistatic interactions among genes, the involvement of epigenetic regulators, and the variable expressivity of traits. Although quantitative trait locus (QTL) analysis has been a powerful tool for localizing the chromosomal regions underlying complex traits, systematically identifying the causal genes remains challenging. Here, through its application to plasma levels of high-density lipoprotein cholesterol (HDL) in mice, we demonstrate a strategy for narrowing QTL that utilizes comparative genomics and bioinformatics techniques. We show how QTL detected in multiple crosses are subjected to both combined cross analysis and haplotype block analysis; how QTL from one species are mapped to the concordant regions in another species; and how genomewide scans associating haplotype groups with their phenotypes can be used to prioritize the narrowed regions. Then we illustrate how these individual methods for narrowing QTL can be systematically integrated for mouse chromosomes 12 and 15, resulting in a significantly reduced number of candidate genes, often from hundreds to <10. Finally, we give an example of how additional bioinformatics resources can be combined with experiments to determine the most likely quantitative trait genes.