Genetic analysis in the Collaborative Cross breeding population.

Genetic reference populations in model organisms are critical resources for systems genetic analysis of disease related phenotypes. The breeding history of these inbred panels may influence detectable allelic and phenotypic diversity. The existing panel of common inbred strains reflects historical selection biases, and existing recombinant inbred panels have low allelic diversity. All such populations may be subject to consequences of inbreeding depression. The Collaborative Cross (CC) is a mouse reference population with high allelic diversity that is being constructed using a randomized breeding design that systematically outcrosses eight founder strains, followed by inbreeding to obtain new recombinant inbred strains. Five of the eight founders are common laboratory strains, and three are wild-derived. Since its inception, the partially inbred CC has been characterized for physiological, morphological, and behavioral traits. The construction of this population provided a unique opportunity to observe phenotypic variation as new allelic combinations arose through intercrossing and inbreeding to create new stable genetic combinations. Processes including inbreeding depression and its impact on allelic and phenotypic diversity were assessed. Phenotypic variation in the CC breeding population exceeds that of existing mouse genetic reference populations due to both high founder genetic diversity and novel epistatic combinations. However, some focal evidence of allele purging was detected including a suggestive QTL for litter size in a location of changing allele frequency. Despite these inescapable pressures, high diversity and precision for genetic mapping remain. These results demonstrate the potential of the CC population once completed and highlight implications for development of related populations.

Complex trait analysis of gene expression uncovers polygenic and pleiotropic networks that modulate nervous system function.

Patterns of gene expression in the central nervous system are highly variable and heritable. This genetic variation among normal individuals leads to considerable structural, functional and behavioral differences. We devised a general approach to dissect genetic networks systematically across biological scale, from base pairs to behavior, using a reference population of recombinant inbred strains. We profiled gene expression using Affymetrix oligonucleotide arrays in the BXD recombinant inbred strains, for which we have extensive SNP and haplotype data. We integrated a complementary database comprising 25 years of legacy phenotypic data on these strains. Covariance among gene expression and pharmacological and behavioral traits is often highly significant, corroborates known functional relations and is often generated by common quantitative trait loci. We found that a small number of major-effect quantitative trait loci jointly modulated large sets of transcripts and classical neural phenotypes in patterns specific to each tissue. We developed new analytic and graph theoretical approaches to study shared genetic modulation of networks of traits using gene sets involved in neural synapse function as an example. We built these tools into an open web resource called WebQTL that can be used to test a broad array of hypotheses.

Uncovering regulatory pathways that affect hematopoietic stem cell function using 'genetical genomics'.

We combined large-scale mRNA expression analysis and gene mapping to identify genes and loci that control hematopoietic stem cell (HSC) function. We measured mRNA expression levels in purified HSCs isolated from a panel of densely genotyped recombinant inbred mouse strains. We mapped quantitative trait loci (QTLs) associated with variation in expression of thousands of transcripts. By comparing the physical transcript position with the location of the controlling QTL, we identified polymorphic cis-acting stem cell genes. We also identified multiple trans-acting control loci that modify expression of large numbers of genes. These groups of coregulated transcripts identify pathways that specify variation in stem cells. We illustrate this concept with the identification of candidate genes involved with HSC turnover. We compared expression QTLs in HSCs and brain from the same mice and identified both shared and tissue-specific QTLs. Our data are accessible through WebQTL, a web-based interface that allows custom genetic linkage analysis and identification of coregulated transcripts.

Status and access to the Collaborative Cross population.

The Collaborative Cross (CC) is a panel of recombinant inbred lines derived from eight genetically diverse laboratory inbred strains. Recently, the genetic architecture of the CC population was reported based on the genotype of a single male per line, and other publications reported incompletely inbred CC mice that have been used to map a variety of traits. The three breeding sites, in the US, Israel, and Australia, are actively collaborating to accelerate the inbreeding process through marker-assisted inbreeding and to expedite community access of CC lines deemed to have reached defined thresholds of inbreeding. Plans are now being developed to provide access to this novel genetic reference population through distribution centers. Here we provide a description of the distribution efforts by the University of North Carolina Systems Genetics Core, Tel Aviv University, Israel and the University of Western Australia.

Combined expression trait correlations and expression quantitative trait locus mapping.

Coordinated regulation of gene expression levels across a series of experimental conditions provides valuable information about the functions of correlated transcripts. The consideration of gene expression correlation over a time or tissue dimension has proved valuable in predicting gene function. Here, we consider correlations over a genetic dimension. In addition to identifying coregulated genes, the genetic dimension also supplies us with information about the genomic locations of putative regulatory loci. We calculated correlations among approximately 45,000 expression traits derived from 60 individuals in an F2 sample segregating for obesity and diabetes. By combining the correlation results with linkage mapping information, we were able to identify regulatory networks, make functional predictions for uncharacterized genes, and characterize novel members of known pathways. We found evidence of coordinate regulation of 174 G protein-coupled receptor protein signaling pathway expression traits. Of the 174 traits, 50 had their major LOD peak within 10 cM of a locus on Chromosome 2, and 81 others had a secondary peak in this region. We also characterized a Riken cDNA clone that showed strong correlation with stearoyl-CoA desaturase 1 expression. Experimental validation confirmed that this clone is involved in the regulation of lipid metabolism. We conclude that trait correlation combined with linkage mapping can reveal regulatory networks that would otherwise be missed if we studied only mRNA traits with statistically significant linkages in this small cross. The combined analysis is more sensitive compared with linkage mapping alone.

An informatics approach to systems neurogenetics.

We outline the theory behind complex trait analysis and systems genetics and describe web-accessible resources including GeneNetwork (GN) that can be used for rapid exploratory analysis and hypothesis testing. GN, in particular, is a tightly integrated suite of bioinformatics tools and data sets, which supports the investigation of complex networks of gene variants, molecules, and cellular processes that modulate complex traits, including behavior and disease susceptibility. Using various statistical tools, users are able to analyze gene expression in various brain regions and tissues, map loci that modulate these traits, and explore genetic covariance among traits. Taken together, these tools enable the user to begin to assess complex interactions of gene networks, and facilitate analysis of traits using a systems approach.

Quantitative trait locus analysis using recombinant inbred intercrosses: theoretical and empirical considerations.

We describe a new approach, called recombinant inbred intercross (RIX) mapping, that extends the power of recombinant inbred (RI) lines to provide sensitive detection of quantitative trait loci (QTL) responsible for complex genetic and nongenetic interactions. RIXs are generated by producing F1 hybrids between all or a subset of parental RI lines. By dramatically extending the number of unique, reproducible genomes, RIXs share some of the best properties of both the parental RI and F2 mapping panels. These attributes make the RIX method ideally suited for experiments requiring analysis of multiple parameters, under different environmental conditions and/or temporal sampling. However, since any pair of RIX genomes shares either one or no parental RIs, this cross introduces an unusual population structure requiring special computational approaches for analysis. Herein, we propose an efficient statistical procedure for QTL mapping with RIXs and describe a novel empirical permutation procedure to assess genome-wide significance. This procedure will also be applicable to diallel crosses. Extensive simulations using strain distribution patterns from CXB, AXB/BXA, and BXD mouse RI lines show the theoretical power of the RIX approach and the analysis of CXB RIXs demonstrates the limitations of this procedure when using small RI panels.

Genetic networks controlling retinal injury.

PURPOSE: The present study defines genomic loci underlying coordinate changes in gene expression following retinal injury. METHODS: A group of acute phase genes expressed in diverse nervous system tissues was defined by combining microarray results from injury studies from rat retina, brain, and spinal cord. Genomic loci regulating the brain expression of acute phase genes were identified using a panel of BXD recombinant inbred (RI) mouse strains. Candidate upstream regulators within a locus were defined using single nucleotide polymorphism databases and promoter motif databases. RESULTS: The acute phase response of rat retina, brain, and spinal cord was dominated by transcription factors. Three genomic loci control transcript expression of acute phase genes in brains of BXD RI mouse strains. One locus was identified on chromosome 12 and was highly correlated with the expression of classic acute phase genes. Within the locus we identified the inhibitor of DNA binding 2 (Id2) as a candidate upstream regulator. Id2 was upregulated as an acute phase transcript in injury models of rat retina, brain, and spinal cord. CONCLUSIONS: We defined a group of transcriptional changes associated with the retinal acute injury response. Using genetic linkage analysis of natural transcript variation, we identified regulatory loci and candidate regulators that control transcript levels of acute phase genes.

An ENU-induced mutation in Rs1h causes disruption of retinal structure and function.

PURPOSE: The 44TNJ mutant mouse was generated by the Tennessee Mouse Genome Consortium (TMGC) using an ENU-based mutagenesis screen to produce recessive mutations that affect the eye and brain. Herein we present its retinal phenotype and genetic basis. METHODS: Fourth generation offspring (G4) and confirmed mutants were examined using slit lamp biomicroscopy, funduscopy, histology, immunohistochemistry, and electroretinography (ERG). 44TNJ mutant mice were crossed to C3BLiA or DBA/2 mice for chromosomal mapping purposes. Linkage analysis by PCR-based microsatellite marker genotyping was used to identify the disease locus. The Rs1h cDNA and its genomic DNA were sequenced directly. RESULTS: The 44TNJ pedigree was the first mutant pedigree identified by the ocular phenotyping domain of the TMGC. Examination of the fundus revealed numerous small and homogeneous intraretinal microflecks in the peripapillary region, which became courser and more irregular in the periphery. Males were typically more affected than females. Histology and immunohistochemistry revealed a disruption of the lamination of the retina, particularly at both margins of the outer nuclear layer, along with reduced calbindin immunostaining. ERG analyses revealed reduced amplitudes of both a-waves and b-waves. Linkage analysis mapped the 44TNJ mutation to the X chromosome close to the marker DXMit117. Sequence analysis of the positional candidate gene Rs1h revealed a T->C exchange at the second base of intron 2 of the Rs1h gene. CONCLUSIONS: We have generated and characterized a mutant mouse line that was produced using ENU-based mutagenesis. The 44TNJ pedigree manifests with photoreceptor dysfunction and concurrent structural and functional aberrations at the post-receptoral level. Genetic analysis revealed a mutation in Rs1h, making this the first murine model of X-linked retinoschisis in which the gene is expressed.

WebQTL: web-based complex trait analysis.

WebQTL is a website that combines databases of complex traits with fast software for mapping quantitative trait loci (QTLs) and for searching for correlations among traits. WebQTL also includes well-curated genotype data for five sets of mouse recombinant inbred (RI) lines. Thus, to identify QTLs, users need provide only quantitative trait data from one of the supported populations. The WebQTL databases include both biological traits--neuroanatomical, pharmacological, and behavioral traits--and microarray-based gene expression data from BXD RI lines. A search function finds correlations between RNA expression and biological traits, and mapping functions find QTLs for either type of trait. The WebQTL service is available at http://www.webqtl.org/.

Genetic correlates of gene expression in recombinant inbred strains: a relational model system to explore neurobehavioral phenotypes.

Full genome sequencing, high-density genotyping, expanding sets of microarray assays, and systematic phenotyping of neuroanatomical and behavioral traits are producing a wealth of data on the mouse central nervous system (CNS). These disparate resources are still poorly integrated. One solution is to acquire these data using a common reference population of isogenic lines of mice, providing a point of integration between the data types. Recombinant inbred (RI) mice, derived through inbreeding of progeny from an inbred cross, are a powerful tool for complex trait mapping and analysis of the challenging phenotypes of neuroscientific interest. These isogenic RI lines are a retrievable genetic resource that can be repeatedly studied using a wide variety of assays. Diverse data sets can be related through fixed and known genomes, using tools such as the interactive web-based system for complex trait analysis, www.WebQTL.org. In this report, we demonstrate the use of WebQTL to explore complex interactions among a wide variety of traits--from mRNA transcripts to the impressive behavioral and pharmacological variation among RI strains. The relational approach exploiting a common set of strains facilitates study of multiple effects of single genes (pleiotropy) without a priori hypotheses required. Here we demonstrate the power of this technique through genetic correlation of gene expression with a database of neurobehavioral phenotypes collected in these strains of mice through more than 20 years of experimentation. By repeatedly studying the same panel of mice, early data can be re-examined in light of technological advances unforeseen at the time of their initial collection.

Informatics center for mouse genomics: the dissection of complex traits of the nervous system.

In recent years, there has been an explosion in the number of tools and techniques available to researchers interested in exploring the genetic basis of all aspects of central nervous system (CNS) development and function. Here, we exploit a powerful new reductionist approach to explore the genetic basis of the very significant structural and molecular differences between the brains of different strains of mice, called either complex trait or quantitative trait loci (QTL) analysis. Our specific focus has been to provide universal access over the web to tools for the genetic dissection of complex traits of the CNS--tools that allow researchers to map genes that modulate phenotypes at a variety of levels ranging from the molecular all the way to the anatomy of the entire brain. Our website, The Mouse Brain Library (MBL; http://mbl.org) is comprised of four interrelated components that are designed to support this goal: The Brain Library, iScope, Neurocartographer, and WebQTL. The centerpiece of the MBL is an image database of histologically prepared museum-quality slides representing nearly 2000 mice from over 120 strains--a library suitable for stereologic analysis of regional volume. The iScope provides fast access to the entire slide collection using streaming video technology, enabling neuroscientists to acquire high-magnification images of any CNS region for any of the mice in the MBL. Neurocartographer provides automatic segmentation of images from the MBL by warping precisely delineated boundaries from a 3D atlas of the mouse brain. Finally, WebQTL provides statistical and graphical analysis of linkage between phenotypes and genotypes.

The Collaborative Cross at Oak Ridge National Laboratory: developing a powerful resource for systems genetics.

Complex traits and disease comorbidity in humans and in model organisms are the result of naturally occurring polymorphisms that interact with each other and with the environment. To ensure the availability of resources needed to investigate biomolecular networks and systems-level phenotypes underlying complex traits, we have initiated breeding of a new genetic reference population of mice, the Collaborative Cross. This population has been designed to optimally support systems genetics analysis. Its novel and important features include a high level of genetic diversity, a large population size to ensure sufficient power in high-dimensional studies, and high mapping precision through accumulation of independent recombination events. Implementation of the Collaborative Cross has been ongoing at the Oak Ridge National Laboratory (ORNL) since May 2005. Production has been systematically managed using a software-assisted breeding program with fully traceable lineages, performed in a controlled environment. Currently, there are 650 lines in production, and close to 200 lines are now beyond their seventh generation of inbreeding. Retired breeders enter a high-throughput phenotyping protocol and DNA samples are banked for analyses of recombination history, allele drift and loss, and population structure. Herein we present a progress report of the Collaborative Cross breeding program at ORNL and a description of the kinds of investigations that this resource will support.

Genome-level analysis of genetic regulation of liver gene expression networks.

UNLABELLED: The liver is the primary site for the metabolism of nutrients, drugs, and chemical agents. Although metabolic pathways are complex and tightly regulated, genetic variation among individuals, reflected in variations in gene expression levels, introduces complexity into research on liver disease. This study dissected genetic networks that control liver gene expression through the combination of large-scale quantitative mRNA expression analysis with genetic mapping in a reference population of BXD recombinant inbred mouse strains for which extensive single-nucleotide polymorphism, haplotype, and phenotypic data are publicly available. We profiled gene expression in livers of naive mice of both sexes from C57BL/6J, DBA/2J, B6D2F1, and 37 BXD strains using Agilent oligonucleotide microarrays. These data were used to map quantitative trait loci (QTLs) responsible for variations in the expression of about 19,000 transcripts. We identified polymorphic local and distant QTLs, including several loci that control the expression of large numbers of genes in liver, by comparing the physical transcript position with the location of the controlling QTL. CONCLUSION: The data are available through a public web-based resource (www.genenetwork.org) that allows custom data mining, identification of coregulated transcripts and correlated phenotypes, cross-tissue, and cross-species comparisons, as well as testing of a broad array of hypotheses.

Integrative genetic analysis of transcription modules: towards filling the gap between genetic loci and inherited traits.

Genetic loci that regulate inherited traits are routinely identified using quantitative trait locus (QTL) mapping methods. However, the genotype-phenotype associations do not provide information on the gene expression program through which the genetic loci regulate the traits. Transcription modules are 'self-consistent regulatory units' and are closely related to the modular components of gene regulatory network [Ihmels, J., Friedlander, G., Bergmann, S., Sarig, O., Ziv, Y. and Barkai, N. (2002) Revealing modular organization in the yeast transcriptional network. Nat. Genet., 31, 370-377; Segal, E., Shapira, M., Regev, A., Pe'er, D., Botstein, D., Koller, D. and Friedman, N. (2003) Module networks: identifying regulatory modules and their condition-specific regulators from gene expression data. Nat. Genet., 34, 166-176]. We used genome-wide genotype and gene expression data of a genetic reference population that consists of mice of 32 recombinant inbred strains to identify the transcription modules and the genetic loci regulating them. Twenty-nine transcription modules defined by genetic variations were identified. Statistically significant associations between the transcription modules and 18 classical physiological and behavioral traits were found. Genome-wide interval mapping showed that major QTLs regulating the transcription modules are often co-localized with the QTLs regulating the associated classical traits. The association and the possible co-regulation of the classical trait and transcription module indicate that the transcription module may be involved in the gene pathways connecting the QTL and the classical trait. Our results show that a transcription module may associate with multiple seemingly unrelated classical traits and a classical trait may associate with different modules. Literature mining results provided strong independent evidences for the relations among genes of the transcription modules, genes in the regions of the QTLs regulating the transcription modules and the keywords representing the classical traits.

How replicable are mRNA expression QTL?

Applying quantitative trait analysis methods to genome-wide microarray-derived mRNA expression phenotypes in segregating populations is a valuable tool in the attempt to link high-level traits to their molecular causes. The massive multiple-testing issues involved in analyzing these data make the correct level of confidence to place in mRNA abundance quantitative trait loci (QTL) a difficult problem. We use a unique resource to directly test mRNA abundance QTL replicability in mice: paired recombinant inbred (RI) and F(2) data sets derived from C57BL/6J (B6) and DBA/2J (D2) inbred strains and phenotyped using the same Affymetrix arrays. We have one forebrain and one striatum data set pair. We describe QTL replication at varying stringencies in these data. For instance, 78% of mRNA expression QTL (eQTL) with genome-wide adjusted p < or = 0.0001 in RI data replicate at a genome-wide adjusted p < 0.05 or better. Replicated QTL are disproportionately putatively cis-acting, and approximately 75% have higher apparent expression levels associated with B6 genotypes, which may be partly due to probe set generation using B6 sequence. Finally, we note that while trans-acting QTL do not replicate well between data sets in general, at least one cluster of trans-acting QTL on distal Chr 1 is notably preserved between data sets.

Reliability of statistical associations between genes and disease.

Many statistical associations between a disease and alleles of specific genes have proven to be irreproducible. In part, this irreproducibility can be attributed to a lack of replication before publication and the fact that, until recently, the relationship between statistical significance and various measures of reproducibility was not widely understood. This review proposes a classification system, the Better Associations for Disease and GEnes (BADGE) system, for describing genetic associations. The BADGE classes, first class through fifth class, are based on the P value of the association. A first-class association, with P < 2 x 10(-7), is expected to be reproducible even in the absence of other evidence supporting the association. A fifth-class association corresponds to conventional statistical significance (P < 5 x 10(-2)), which provides almost no assurance of reproducibility. Three intervening classes, described as second-, third-, and fourth-class associations, are defined by P values separated by factors of 20 or 25 from these extremes.

Weighting by heritability for detection of quantitative trait loci with microarray estimates of gene expression.

Heritable differences in transcribed RNA levels can be mapped as quantitative trait loci (QTLs). Transcribed RNA levels are often measured by hybridization to microarrays of oligonucleotide probes, in which each transcript is represented by multiple probes. The use of recombinant inbred lines allows an estimate of the heritability of expression measured by individual probes. This heritability varies greatly. We have tested heritability-weighted averages to define expression of a transcript and found that these allow detection of more QTLs than previously described methods.