BXD Recombinant Inbred Strains Research Articles

A Highly Polymorphic Meiotic Recombination Mouse Hot Spot Exhibits Incomplete Repair

The recent mapping of recombination hot spots in the human genome has demonstrated that crossover is a nonrandom process that occurs at well-defined positions along chromosomes. However, the mechanisms that direct hot-spot turnover in complex mammalian genomes are poorly understood. Analyses of the human genome are impaired by the inability to genetically dissect and molecularly manipulate recombinogenic regions to test their roles in regulating hot spots. Here, using the BXD recombinant inbred strains as a crossover library, three new recombination hot spots have been identified on mouse chromosome 19. Analyses of a highly polymorphic recombination hot spot (HS22) revealed that approximately 4% of recombinant molecules display complex and incomplete repair with discontinuous conversion tracts, as well as persistent heteroduplex DNA at crossover sites in mature spermatozoa. Also, sequence analysis of the wild house mouse revealed instability at the center of this hot spot. This suggests that complete repair is not required for completion of mammalian meiosis, a scenario that leaves duplex DNA containing mismatches at crossover sites.

IN THIS ISSUE

[⇓][1] Whereas all C57BL/6J (B6) mice survive infection by the obligate intracellular bacterium Chlamydia psittaci , DBA/2J (D2) mice are susceptible and die. Miyairi et al. (p. [1814][2] ) found that three of 24 BXD recombinant inbred strains, derived from intercrossing the two parent strains,

The p47 GTPases Iigp2 and Irgb10 Regulate Innate Immunity and Inflammation to Murine Chlamydia psittaci Infection

C57BL/6J mice were 10(5)-fold more resistant to Chlamydia psittaci infection than DBA/2J mice by LD(100) determinations. Linkage analysis using BXD recombinant inbred strains revealed a single effector locus at a 1.5-Mbp region on chromosome 11 encoding a cluster of three p47 GTPases (Irgb10, Igtp, and Iigp2). Western blots of infected tissue showed that Irgb10 was elevated in resistant mice and one of the two possible Iigp2 protein isoforms was preferentially expressed in susceptible mice. The BXD39 strain, susceptible at Irgb10 and resistant at Iigp2, had an intermediate phenotype implicating the nonredundant role of these p47 GTPases. C57BL/6J and DBA/2J exhibited a difference in IFN-gamma-dependent chlamydial control, which was reversible by Iigp2 small interfering RNA knockdown. Microarrays of infected peritoneal lavage revealed >10-fold up-regulation of neutrophil-recruiting chemokines in susceptible mice and >100-fold increase in macrophage differentiation genes in resistant mice, indicating that the susceptibility pattern involves the stimulation of different inflammatory cell-recruiting pathways. Massive neutrophil recruitment was seen in susceptible mice by histology and flow cytometry, and neutrophil chemokine receptor (CXCR2) knockout mice on a susceptible background survived a lethal challenge, confirming that neutrophil recruitment was required for susceptibility. Congenic Igtp knockout mice also susceptible at Irgb10 and Iigp2 on a resistant background recruited neutrophils and succumbed to infection. We conclude that Irgb10 and Iigp2 act together to confer differential susceptibility against murine chlamydial infection. Data indicate that these p47 GTPases have cell-autonomous effects that result in vastly different inflammatory stimulations, leading to either recovery or death.

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

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. 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.

An integrative genomics strategy for systematic characterization of genetic loci modulating phenotypes

Naturally occurring genetic variations may affect certain phenotypes through influencing transcript levels of the genes that are causally related to those phenotypes. Genomic regions harboring common sequence variants that modulate gene expression can be mapped as quantitative trait loci (QTLs) using a newly developed genetical genomics approach. This enables a new strategy for systematically mapping novel genetic loci underlying various phenotypes. In this work, we started from a seed set of genes with variants that are known to affect behavioral and neurological phenotypes (as recorded in Mammalian Phenotype Ontology Database) and used microarrays to analyze their expression levels in brain samples of a panel of BXD recombinant inbred mouse strains. We then systematically mapped the QTLs controlling the expression of these genes. Candidate causal genes in the QTL intervals were evaluated for evidence of functional genetic polymorphisms. Using this method, we were able to predict novel genetic loci and causal genes for a number of behavioral and neurological phenotypes. Lines of independent evidence supporting some of our results were provided by transcription factor binding site analysis and by biomedical literature. This strategy integrates gene-phenotype relations from decades of experimental mutagenesis studies and new genomic resources to provide an approach to rapidly expand knowledge on genetic loci modulating phenotypes.

Quantitative trait loci linked to thalamus and cortex gray matter volumes in BXD recombinant inbred mice

To investigate whether there are separate or shared genetic influences on the development of the thalamus and cerebral cortex, we identified quantitative trait loci (QTLs) for relevant structural volumes in BXD recombinant inbred (RI) strains of mice. In 34 BXD RI strains and two parental strains (C57BL/6J and DBA/2J), we measured the volumes of the entire thalamus and cortex gray matter using point counting and Cavalieri's rule. Heritability was calculated using analysis of variance (ANOVA), and QTL analysis was carried out using WebQTL (http://www.genenetwork.org). The heritability of thalamus volume was 36%, and three suggestive QTLs for thalamus volume were identified on chromosomes 10, 11 and 16. The heritability of cortical gray matter was 43%, and four suggestive QTLs for cortex gray matter volume were identified on chromosomes 2, 8, 16 and 19. The genetic correlation between thalamus and cortex gray matter volumes was 0.64. Also, a single QTL on chromosome 16 (D16Mit100) was identified for thalamus volume, cortex gray matter volume and Morris water maze search-time preference (r=0.71). These results suggest that there are separate and shared genetic influences on the development of the thalamus and cerebral cortex.

QTL Analysis of Trabecular Bone in BXD F2 and RI Mice

A sample of 693 mice was used to identify regions of the mouse genome associated with trabecular bone architecture as measured using microCT. QTLs for bone in the proximal tibial metaphysis were identified on several chromosomes indicating regions containing genes that regulate properties of trabecular bone. Age-related osteoporosis is a condition of major concern because of the morbidity and mortality associated with osteoporotic fractures in humans. Osteoporosis is characterized by reduced bone density, strength, and altered trabecular architecture, all of which are quantitative traits resulting from the actions of many genes working in concert with each other and the environment over the lifespan. microCT gives accurate measures of trabecular bone architecture providing phenotypic data related to bone volume and trabecular morphology. The primary objective of this research was to identify chromosomal regions called quantitative trait loci (QTLs) that contain genes influencing trabecular architecture as measured by microCT. The study used crosses between C57BL/6J (B6) and DBA/2J (D2) as progenitor strains of a second filial (F2) generation (n = 141 males and 148 females) and 23 BXD recombinant inbred (RI) strains (n approximately 9 of each sex per strain). The proximal tibial metaphyses of the 200-day-old mice were analyzed by microCT to assess phenotypic traits characterizing trabecular bone, including bone volume fraction, trabecular connectivity, and quantitative measures of trabecular orientation and anisotropy. Heritabilities were calculated and QTLs were identified using composite interval mapping. A number of phenotypes were found to be highly heritable. Heritability values for measured phenotypes using RI strains ranged from 0.15 for degree of anisotropy in females to 0.51 for connectivity density in females and total volume in males. Significant and confirmed QTLs, with LOD scores 4.3 in the F2 cohort and 1.5 in the corresponding RI cohort were found on chromosomes 1 (43 cM), 5 (44 cM), 6 (20 cM), and 8 (49 cM). Other QTLs with LOD scores ranging from 2.8 to 6.9 in the F2 analyses were found on chromosomes 1, 5, 6, 8, 9, and 12. QTLs were identified using data sets comprised of both male and female quantitative traits, suggesting similar genetic action in both sexes, whereas others seemed to be associated exclusively with one sex or the other, suggesting the possibility of sex-dependent effects. Identification of the genes underlying these QTLs may lead to improvements in recognizing individuals most at risk for developing osteoporosis and in the design of new therapeutic interventions.

QTL influencing baseline hematocrit in the C57BL/6J and DBA/2J lineage: age-related effects

Baseline serum hematocrit varies substantially in the population. While additive genetic factors account for a large part of this variability, little is known about the genetic architecture underlying the trait. Because hematocrit levels vary with age, it is plausible that quantitative trait loci (QTL) that influence the phenotype also show an age-specific profile. To investigate this possibility, hematocrit was measured in three different age cohorts of mice (150, 450, and 750 days) of the C57BL/6J (B6) and the DBA2/J (D2) lineage. QTL were searched in the B6D2F(2) intercross and the BXD recombinant inbred (RI) strains. The effects of these QTL were explored across the different age groups. On the phenotypic level, baseline serum hematocrit declines with age in a sex-specific manner. In the B6D2F(2) intercross, suggestive QTL that influence the phenotype were located on Chromosomes (Chr) 1, 2, 7, 11, 13, and 16. With the exception of the QTL on Chr 2, all of these QTL exerted their largest effect at 750 days. The QTL on Chr 1, 2, 7, 11 and 16 were confirmed in the BXD RIs in a sex- and age-specific manner. Linkage analysis in the BXD RIs revealed an additional significant QTL on Chr 19. Baseline serum hematocrit is influenced by several QTL that appear to vary with the age and sex of the animal. These QTL primarily overlap with QTL that have been shown to regulate hematopoietic stem cell phenotypes.

Genetic determinants of weight of fast- and slow-twitch skeletal muscles in old mice

The main goal of the study was to explore the genetic architecture underlying muscle weight in old mice. Weight of soleus, tibialis anterior (TA), extensor digitorum longus (EDL), and gastrocnemius muscles was measured in the C57BL/6J (B6) and DBA/2J (D2) strains and derivative generations: a panel of the BXD recombinant inbred (RI) strains and a B6D2 F(2) intercross at the age of 800 days. The between-strain difference in muscle weight (B6 > D2) ranged between 16% and 38%. Linkage analysis identified suggestive quantitative trait loci (QTL) on Chromosomes (Chr) 2, 6, 7, 8, 19, and X that influenced muscle weight in the 800-day-old group. Comparison of weights at 200, 500, and 800 days revealed a variable effect of age among the four muscles. Linkage analysis in the B6D2 F(2) population combined across the three different age groups identified muscle-, sex-, and age-specific QTL on Chr 1, 2, 3, 5, 6, 8, 9, 11, 13, 17, X, and Y. Genetic factors that influence the rate of weight change (within-strain weight difference at two ages) over the lifespan of BXD RIs were mapped to the markers D2Mit369 and D3Mit130 at the genome-wide p < 0.05 for TA muscle in males (between 200 and 800 days) and females (between 500 and 800 days), respectively. Analysis of all age groups supported previous findings that the genetic effects may be muscle-, age-, and sex-specific.

Dopamine D2 receptor binding, Drd2 expression and the number of dopamine neurons in the BXD recombinant inbred series: genetic relationships to alcohol and other drug associated phenotypes.

It has not been established to what extent the natural variation in dopamine systems contribute to the variation in ethanol response. The current study addresses this issue by measuring D2 dopamine (DA) receptor binding, the expression of Drd2, the number of midbrain DA neurons in the BXD recombinant inbred (RI) series and then compares these strain means with those previously reported for a variety of ethanol and other drug-related phenotypes. Data were collected for 21 to 23 of the BXD RI strains and the parental strains. D2 DA receptor autoradiography was performed using 125I-epidepride as the ligand [ Kanes S, Dains K, Cipp L, Gatley J, Hitzemann B, Rasmussen E, Sanderson S, Silverman S, Hitzemann R (1996) Mapping the genes for haloperidol-induced catalepsy. J Pharmacol Exp Ther 277:1016-1025]. Drd2 expression was measured using the Affymetrix oligoarray system. Immunocytochemical techniques were used to determine the number of midbrain DA neurons [Hitzemann B, Dains K, Hitzemann R (1994) Further studies on the relationship between dopamine cell density and haloperidol response. J Pharmacol Exp Ther 271:969-976]. The range of difference in receptor binding for the RI strains was approximately 2-fold in all regions examined, the core, the shell of the nucleus accumbens (NAc) and the dorsomedial caudate-putamen (CPu); heritability in all regions was moderate--(h2 approximately 0.35). Drd2 expression in forebrain samples from the RI and parental strains ranged 1.5- to h2-fold and was moderate-0.47. Variation in the number of tyrosine hydroxylase (TH) positive neurons was moderate, 41% and 26% and h2 was low--0.19 and 0.15 for the ventral tegmental area (VTA) and substantia nigra compacta (SNc), respectively. Significant correlations were found between D2 DA receptor binding and the low dose (1.33 g/kg) ethanol stimulant response. (p < 0.002) and between expression and conditioned place preference (CPP) (p < 0.0005). No significant correlations were detected between ethanol preference and either receptor binding or Drd2 expression; however, a significant correlation was found between preference and Ncam expression. Ncam is approximately 0.2 Mb from Drd2. Overall, the data suggest ethanol preference and CPP are associated with the expression of Drd2 or closely linked genetic loci.

Sensitivity to the locomotor‐stimulant effects of ethanol and allopregnanolone: a quantitative trait locus study of common genetic influence

Previous studies have suggested that common genetic mechanisms influence sensitivity to the locomotor-stimulant effects of ethanol and allopregnanolone. We conducted two quantitative trait locus (QTL) studies to identify chromosomal regions that harbor genes that influence locomotor response to ethanol (2 g/kg) and allopregnanolone (17 mg/kg) using F2 crosses between C57BL/6J and DBA/2J mice. Because our previous data from the BXD recombinant inbred strains had indicated that chromosome 2 contained QTL for sensitivity to the locomotor-stimulant effects of both ethanol and allopregnanolone, we also tested reciprocal chromosome 2 congenic strains for sensitivity to the locomotor-stimulant effects of both drugs. The F2 analysis for ethanol sensitivity identified significant QTL on chromosomes 1 and 2 and suggestive QTL on chromosomes 5 and 9. The analysis of the allopregnanolone F2 study identified suggestive QTL on chromosomes 3, 5 and 12. Suggestive evidence for a female-specific QTL on chromosome 2 was also found. The studies of congenic mouse strains indicated that both the congenic strains captured one or more QTL for sensitivity to the locomotor-stimulant effects of both ethanol (2 g/kg) and allopregnanolone (17 mg/kg). When Fisher's method was used to combine the P values for the RI, F2 and congenic studies of the chromosome 2 QTL, cumulative probability scores of 9.6 x 10(-15) for ethanol and 7.7 x 10(-7) for allopregnanolone were obtained. These results confirm the presence of QTL for ethanol and allopregnanolone sensitivity in a common region of chromosome 2 and suggest possible pleiotropic genetic influence on sensitivity to these drugs.

Anomaly of anatomical origin of soleus muscle: A mouse model

In the laboratory mouse, the soleus muscle arises at the head of the fibula and inserts via the Achilles tendon on the tuber calcanei together with the gastrocnemius muscle. During routine dissection of mice from the BXD recombinant inbred (RI) strains, we found that the soleus often originated from the lateral epicondyle of the femur instead of the head of the fibula. This soleus femoral attachment anomaly (SFAA) changes the soleus from being a single-joint to a two-joint muscle. The incidence of SFAA was 45% in the BXD38 RI strain. Bilateral inspection indicated that SFAA may be present unilaterally or bilaterally within an individual mouse. We explored the effect of SFAA on muscle weight in mice with unilateral expression. The weight of SFAA soleii was significantly less (P < 0.01) than that of the soleii with normal attachment by 6% (females) and 14% (males). Similar anatomical anomalies of the soleus muscle have been noted in humans. The mouse model will provide the means to explore the physiological consequences and genetic basis for such anomalies.

Serum alkaline phosphatase activity is regulated by a chromosomal region containing the alkaline phosphatase 2 gene (Akp2) in C57BL/6J and DBA/2J mice

Quantitative trait locus (QTL) analyses were conducted to identify chromosomal regions that contribute to variability in serum alkaline phosphatase (AP) enzyme activity in mice derived from the C57BL/6J (B6) and DBA/2J (D2) inbred strains. Serum AP was measured in 400 B6D2 F2 mice at 5 mo and 400 B6D2 F2 mice at 15 mo of age that were genotyped at 96 microsatellite markers, and in 19 BXD recombinant inbred (RI) strains at 5 mo of age. A QTL on the distal end of chromosome 4 was present in all sex- and age-specific analyses with a peak logarithm of odds (LOD) score of 20.36 at 58.51 cM. The Akp2 gene, which encodes the major serum AP isozyme, falls within this QTL region at 70.2 cM where the LOD score reached 13.2 (LOD significance level set at 4.3). Serum AP activity was directly related to the number of D2 alleles of a single nucleotide polymorphism in the 5'-flanking region of the Akp2 gene, although no strain-related differences in hepatic expression of Akp2 RNA were found. A variety of sequence polymorphisms in this chromosomal region could be responsible for the differences in serum AP activity; the Akp2 gene, however, with several known amino acid substitutions between protein sequences of the B6 and D2 strains, is a leading candidate.

The Syntaxin Binding Protein 1 Gene (Stxbp1 ) Is a Candidate for an Ethanol Preference Drinking Locus on Mouse Chromosome 2

We previously mapped a quantitative trait locus (QTL) for ethanol preference drinking to mouse chromosome 2 (mapped with high confidence, LOD = 15.5, p = 3 x 10(-16)). The specific gene(s) in the QTL interval responsible for phenotypic variation in ethanol preference drinking has not been identified. In the current study, we investigated the association of the syntaxin binding protein 1 gene (Stxbp1) with ethanol preference drinking and other ethanol traits using a panel of B6 x D2 (BXD) recombinant inbred (RI) strains derived from the C57BL/6J (B6) and DBA/2J (D2) inbred mouse strains. Confirmation analyses for ethanol consumption and withdrawal were performed using a large B6D2 F2 cross, short-term selected lines derived from the B6 and D2 progenitor strains, and standard inbred strains. BXD RI strain analysis detected provisional associations between Stxbp1 molecular variants and ethanol consumption, as well as severity of acute ethanol withdrawal, ethanol-conditioned taste aversion, and ethanol-induced hypothermia. Confirmation analyses using three independent genetic models supported the involvement of Stxbp1 in ethanol preference drinking but not in ethanol withdrawal. Stxbp1 encodes a Sec1/Munc18-type protein essential for vesicular neurotransmitter release. The present study provides supporting evidence for the involvement of Stxbp1 in ethanol preference drinking.

WebQTL: An Internet resource for the integrative genetic analysis of gene expression and pain related phenotypes

Understanding individual differences in pain related phenotypes requires integration of data across levels of biological scale—a challenging task made readily possible by deployment of genetic reference populations such as the BXD recombinant inbred strains. These mice, derived from the C57BL/6J and DBA/2J strains, form a reference population that has been studied and expanded throughout the past 20 years. Recently, we have quantified mRNA expression levels in these strains and integrated molecular, neuropharmacological, and behavioral traits into a free, public, web-accessible interface called WebQTL (www.webQTL.org).

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.

Genetic determinants of weight of fast- and slow-twitch skeletal muscle in 500-day-old mice of the C57BL/6J and DBA/2J lineage

C57BL/6J (B6) and DBA/2J (D2) strains and two derivative populations, BXD recombinant inbred strains (BXD RIs) and B6D2F2, were used to explore genetic basis for variation in muscle weight at 500 days of age. In parallel with findings in 200-day-old mice (Lionikas A, Blizard DA, Vandenbergh DJ, Glover MG, Stout JT, Vogler GP, McClearn GE, and Larsson L. Physiol Genomics 16: 141-152, 2003), weight of slow-twitch soleus, mixed gastrocnemius, and fast-twitch tibialis anterior (TA) and extensor digitorum longus (EDL) muscles was 13-22% greater (P < 0.001) in B6 than in D2. Distribution of BXD RI strain means indicated that genetic influence on muscle weight (strain effect P < 0.001, all muscles) was of polygenic origin, and effect of genetic factors differed between males and females (strain-by-sex interaction: P < 0.01 for soleus, EDL; P < 0.05 for TA, gastrocnemius). Linkage analyses in B6D2F2 population identified QTL affecting muscle weight on Chr 1, 2, 6, and 9. Pleiotropic influences were observed for QTL on Chr 1 (soleus, TA), 2 (TA, EDL, gastrocnemius), and 9 (soleus, TA, EDL) and were not related to muscle type (fast/slow-twitch) or function (flexor/extensor). Effect of QTL on Chr 9 on soleus muscle was male specific. QTL on Chr 2 and 6 were previously observed at 200 days of age, whereas QTL on Chr 1 and 9 are novel muscle weight QTL. In summary, muscle weight in B6/D2 lineage is affected by a polygenic system that has variable influences at different ages, between males and females, and across muscles in a manner independent of muscle type.

Genetic Segregation of Spontaneous Erosive Arthritis and Generalized Autoimmune Disease in the BXD2 Recombinant Inbred Strain of Mice

The BXD2 strain of mice is one of approximately 80 BXD recombinant inbred (RI) mouse strains derived from an intercross between C57BL/6J (B6) and DBA/2J (D2) strains. We have discovered that adult BXD2 mice spontaneously develop generalized autoimmune disease, including glomerulonephritis (GN), increased serum titres of rheumatoid factor (RF) and anti-DNA antibody, and a spontaneous erosive arthritis characterized by mononuclear cell infiltration, synovial hyperplasia, and bone and cartilage erosion. The features of lupus and arthritis developed by the BXD2 mice segregate in F2 mice generated by crossing BXD2 mice with the parental B6 and D2 strains. Genetic linkage analysis of the serum levels of anti-DNA and RF by using the BXD RI strains shows that the serum titers of anti-DNA and RF were influenced by a genetic locus on mouse chromosome (Chr) 2 near the marker D2Mit412 (78 cm, 163 Mb) and on Chr 4 near D4Mit146 (53.6 cm, 109 Mb), respectively. Both loci are close to the B-cell hyperactivity, lupus or GN susceptibility loci that have been identified previously. The results of our study suggest that the BXD2 strain of mice is a novel model for complex autoimmune disease that will be useful in identifying the mechanisms critical for the immunopathogenesis and genetic segregation of lupus and erosive arthritis.

P2P‐R expression is genetically coregulated with components of the translation machinery and with PUM2, a translational repressor that associates with the P2P‐R mRNA

P2P-R is a nuclear protein with potential functional roles in the control of gene expression and mitosis. The P2P-R protein also interacts with the p53 and Rb1 tumor suppressor proteins. To search for additional functional associations of P2P-R, we employed the WebQTL database that contains the results of cDNA microarray analysis on forebrain, cerebellum, and hematopoietic stem cell (HSC) specimens of multiple BXD recombinant inbred strains of mice. Using WebQTL, gene products were identified that show genetically based coexpression with P2P-R. Initial studies identified general groups of mRNAs that share common functional roles and high covariation in expression with P2P-R. These functional groups involved the regulation of transcription, nucleotide binding, translation control, and ion transport. The findings related to translational mechanisms were further evaluated. In HSCs, expression of P2P-R mRNA demonstrates an impressive expression correlation with a group of gene products associated with translation; high expression of P2P-R specifically was associated with decreased expression of 29 ribosomal protein mRNAs. In all three tissues that were screened using the WebQTL database, a strong positive expression covariance between P2P-R and the Pum2 gene product also was observed. PUM2 is a member of the highly conserved Puf family of RNA binding proteins that often function as gene-specific translation regulators. The ability of Puf proteins to repress translation is mediated by their binding to specific elements located in the 3' untranslated region (UTR) of their target mRNAs. To assess the functional significance of the strong genetic correlation in expression of P2P-R and PUM2, the 3' UTR of the P2P-R mRNA was analyzed and found to contain one perfect consensus and two near-perfect consensus PUM2 binding sequences. PUM2 pull-down methods combined with reverse transcription and RT-PCR confirmed that PUM2 does indeed bind P2P-R mRNA. These results suggest that P2P-R expression may be translationally regulated by PUM2 and that P2P-R may modulate translation by influencing ribosomal protein gene expression. This study represents the first description of a RNA target for mammalian Puf proteins and the first molecular confirmation of information obtained using the WebQTL database.

Convergent analysis of cDNA and short oligomer microarrays, mouse null mutants and bioinformatics resources to study complex traits

Gene expression data sets have recently been exploited to study genetic factors that modulate complex traits. However, it has been challenging to establish a direct link between variation in patterns of gene expression and variation in higher order traits such as neuropharmacological responses and patterns of behavior. Here we illustrate an approach that combines gene expression data with new bioinformatics resources to discover genes that potentially modulate behavior. We have exploited three complementary genetic models to obtain convergent evidence that differential expression of a subset of genes and molecular pathways influences ethanol-induced conditioned taste aversion (CTA). As a first step, cDNA microarrays were used to compare gene expression profiles of two null mutant mouse lines with difference in ethanol-induced aversion. Mice lacking a functional copy of G protein-gated potassium channel subunit 2 (Girk2) show a decrease in the aversive effects of ethanol, whereas preproenkephalin (Penk) null mutant mice show the opposite response. We hypothesize that these behavioral differences are generated in part by alterations in expression downstream of the null alleles. We then exploited the WebQTL databases to examine the genetic covariance between mRNA expression levels and measurements of ethanol-induced CTA in BXD recombinant inbred (RI) strains. Finally, we identified a subset of genes and functional groups associated with ethanol-induced CTA in both null mutant lines and BXD RI strains. Collectively, these approaches highlight the phosphatidylinositol signaling pathway and identify several genes including protein kinase C beta isoform and preproenkephalin in regulation of ethanol- induced conditioned taste aversion. Our results point to the increasing potential of the convergent approach and biological databases to investigate genetic mechanisms of complex traits.