BackgroundThe current study focused on the extent genetic diversity within a species (Mus musculus) affects gene co-expression network structure. To examine this issue, we have created a new mouse resource, a heterogeneous stock (HS) formed from the same eight inbred strains that have been used to create the collaborative cross (CC). The eight inbred strains capture > 90% of the genetic diversity available within the species. For contrast with the HS-CC, a C57BL/6J (B6) × DBA/2J (D2) F2 intercross and the HS4, derived from crossing the B6, D2, BALB/cJ and LP/J strains, were used. Brain (striatum) gene expression data were obtained using the Illumina Mouse WG 6.1 array, and the data sets were interrogated using a weighted gene co-expression network analysis (WGCNA).ResultsGenes reliably detected as expressed were similar in all three data sets as was the variability of expression. As measured by the WGCNA, the modular structure of the transcriptome networks was also preserved both on the basis of module assignment and from the perspective of the topological overlap maps. Details of the HS-CC gene modules are provided; essentially identical results were obtained for the HS4 and F2 modules. Gene ontology annotation of the modules revealed a significant overrepresentation in some modules for neuronal processes, e.g., central nervous system development. Integration with known protein-protein interactions data indicated significant enrichment among co-expressed genes. We also noted significant overlap with markers of central nervous system cell types (neurons, oligodendrocytes and astrocytes). Using the Allen Brain Atlas, we found evidence of spatial co-localization within the striatum for several modules. Finally, for some modules it was possible to detect an enrichment of transcription binding sites. The binding site for Wt1, which is associated with neurodegeneration, was the most significantly overrepresented.ConclusionsDespite the marked differences in genetic diversity, the transcriptome structure was remarkably similar for the F2, HS4 and HS-CC. These data suggest that it should be possible to integrate network data from simple and complex crosses. A careful examination of the HS-CC transcriptome revealed the expected structure for striatal gene expression. Importantly, we demonstrate the integration of anatomical and network expression data.

Although hundreds if not thousands of quantitative trait loci (QTL) have been described for a wide variety of complex traits, only a very small number of these QTLs have been reduced to quantitative trait genes (QTGs) and quantitative trait nucleotides (QTNs). A strategy, Multiple Cross Mapping (MCM), is described for detecting QTGs and QTNs that is based on leveraging the information contained within the haplotype structure of the mouse genome. As described in the current report, the strategy utilizes the six F(2) intercrosses that can be formed from the C57BL/6J (B6), DBA/2J (D2), BALB/cJ (C), and LP/J (LP) inbred mouse strains. Focusing on the phenotype of basal locomotor activity, it was found that in all three B6 intercrosses, a QTL was detected on distal Chromosome (Chr) 1; no QTL was detected in the other three intercrosses, and thus, it was assumed that at the QTL, the C, D2, and LP strains had functionally identical alleles. These intercross data were used to form a simple algorithm for interrogating microsatellite, single nucleotide polymorphism (SNP), brain gene expression, and sequence databases. The results obtained point to Kcnj9 (which has a markedly lower expression in the B6 strain) as being the likely QTG. Further, it is suggested that the lower expression in the B6 strain results from a polymorphism in the 5'-UTR that disrupts the binding of at least three transcription factors. Overall, the method described should be widely applicable to the analysis of QTLs.

Although improvements are needed in the expression databases, the integration of QTL and gene expression analyses seems to have potential as a high-throughput strategy for moving from QTL to QTG.

This study examines the use of multiple cross mapping (MCM) to reduce the interval for an ethanol response QTL on mouse chromosome 1. The phenotype is the acute locomotor response to a 1.5-g/kg i.p. dose of ethanol. The MCM panel consisted of the six unique intercrosses that can be obtained from the C57BL/6J (B6), DBA/2J (D2), BALB/cJ (C) and LP/J (LP) inbred mouse strains (N Ն 600/cross). Ethanol response QTL were detected only with the B6xD2 and B6xC intercrosses. For both crosses, the D2 and C alleles were dominant and decreased ethanol response. The QTL information was used to develop an algorithm for sorting and editing the chromosome 1 Mit microsatellite marker set (http://www.jax.org). This process yielded a cluster of markers between 82 and 85 cM (MGI). Evidence that the QTL was localized in or near this interval was obtained by the analysis of a sample (n Ω 550) of advanced cross heterogenous stock animals. In addition, it was observed that one of the BXD recombinant inbred strains (BXD-32) had a recombination in the interval of interest which produced the expected change in behavior. Overall, the data obtained suggest that the information available within existing genetic maps coupled with MCM data can be used to reduce the QTL interval. In addition, the MCM data set can be used to interrogate gene expression data to estimate which polymorphisms within the interval of interest are relevant to the QTL.

Rationale Previous studies have suggested that there is an inverse genetic relationship between ethanol consumption (two-bottle choice, continuous access) and ethanol withdrawal (e.g., Metten et al., Behav Brain Res 95:113–122, 1998a). Objectives The current study used short-term selective breeding from heterogeneous stock (HS) animals to examine this relationship. The primary goal of the current study was to determine if reciprocal quantitative trait loci (QTLs) could be found in the selectively bred lines. The advantage of detecting QTLs in HS animals is that it is possible to extract a haplotype signature for the QTL, which in turn can be used to narrow the number of candidate genes generated from gene expression and sequence databases (see, e.g., Hitzemann et al., Mamm Genome 14:733–747, 2003). Results Seven reciprocal QTLs were detected on chromosomes (Chr) 1 (two), 3, 6, 11, 16, and 17 that exceeded the nominal LOD threshold of 10; genetic drift, which occurs during selection, dramatically increases the LOD threshold. The proximal Chr 1 QTL was examined in some detail. The haplotype structure of the QTL was such that the LP/J allele was associated with low withdrawal and high consumption. The QTL appears to be located in a gene-poor region between 170 and 173 Mbp. Based on available sequence data, two plausible candidate genes emerge—Nos1ap and Atf6α. Conclusions The data presented here confirm some aspects of the negative genetic relationship between acute ethanol withdrawal and ethanol consumption. The QTL data point to the potential involvement of NO signaling and/or the unfolded protein response.