Recombinant Inbred Strains Research Articles

Tail Tendon Break Time: A Biomarker of Aging?

Research has attempted to identify biomarkers of aging that are predictive of longevity and specific age-related changes during animal life span. Tail tendon break time (TTBT), one presumed biomarker, measures collagen cross-linking, known to increase with age. Significant differences in the rate of increase of TTBT with age have been reported between mouse strains and animal species. We measured both TTBT and longevity in C57BL/6J, DBA/2J, and 23 recombinant inbred (RI) strains (B×D RIs), with TTBT measured at 200, 500, and 800 days of age. Longevity demonstrated considerable variability among these strains (116-951 days). TTBT, also highly variable, increased significantly with age in both sexes and all genotypes. Neither TTBT nor its rate of change correlated significantly with life span. There were suggestive trends for rate of TTBT change to correlate with male longevity and strain longevity to correlate with female TTBT. We conclude that for the range of genetic variation found among these mouse genotypes, TTBT cannot be considered a robust biomarker of longevity.

Quantitative Trait Loci Analysis of Tail Tendon Break Time in Mice of C57BL/6J and DBA/2J Lineage

Tail tendon break time (TTBT), a measure of collagen cross-linking, shown to increase with age differs significantly among inbred strains of mice, indicating underlying genetic influences. This study was aimed to identify quantitative trait loci (QTLs) associated with tail tendon break time at three ages (200, 500, and 800 days of age) for 23 BxD recombinant inbred strains of mice and B6D2F(2) mice derived from C57BL/6J and DBA/2J strains. Heritability estimates were calculated, and QTL analyses were conducted using interval-mapping methods. Mean tail tendon break time values were higher in males and increased nonlinearly with age. Eight total QTLs were nominated in the B6D2F(2) mice at the three measured ages, with the QTL at 800 days confirmed in the recombinant inbred strains. Allelic effect modeling for the identified QTLs suggests differences in gene action between sexes. Candidate genes in the QTL regions include collagen genes and an advanced glycation end-product receptor. The QTLs identified demonstrate influence at some but not all ages.

Systems genetics analysis of mouse chondrocyte differentiation

One of the goals of systems genetics is the reconstruction of gene networks that underlie key processes in development and disease. To identify cartilage gene networks that play an important role in bone development, we used a systems genetics approach that integrated microarray gene expression profiles from cartilage and bone phenotypic data from two sets of recombinant inbred strains. Microarray profiles generated from isolated chondrocytes were used to generate weighted gene coexpression networks. This analysis resulted in the identification of subnetworks (modules) of coexpressed genes that then were examined for relationships with bone geometry and density. One module exhibited significant correlation with femur length (r = 0.416), anteroposterior diameter (r = 0.418), mediolateral diameter (r = 0.576), and bone mineral density (r = 0.475). Highly connected genes (n = 28) from this and other modules were tested in vitro using prechondrocyte ATDC5 cells and RNA interference. Five of the 28 genes were found to play a role in chondrocyte differentiation. Two of these, Hspd1 and Cdkn1a, were known previously to function in chondrocyte development, whereas the other three, Bhlhb9, Cugbp1, and Spcs3, are novel genes. Our integrative analysis provided a systems-level view of cartilage development and identified genes that may be involved in bone development. © 2011 American Society for Bone and Mineral Research.

Genetic Pathways Regulating Glutamate Levels in Retinal Müller Cells

Müller cells serve many functions including the regulation of extracellular glutamate levels. The product of two genes, Slc1a3 [aka solute carrier family 1 (glial high affinity glutamate transporter), member 3] and Glul (aka glutamine synthetase) are the primary role players that transport glutamate into the Müller cell and convert it into glutamine. In this study, we sought to identify the genetic regulation of both genes. Given their tightly coupled biological functions, we predicted that they would be similarly regulated. Using an array of 75 recombinant inbred strains of mice, we determined that Slc1a3 and Glul are differentially regulated by distinct chromosomal regions. Interestingly, despite their independent regulation, gene ontology analysis of tightly correlated genes reveals that the enriched and statistically significant molecular function categories of both directed acyclic graphs have substantial overlap, indicating that the shared functions of correlates of Slc1a3 and Glul include production and usage of ATP.

Identification of a Chr 11 quantitative trait locus that modulates proliferation in the rostral migratory stream of the adult mouse brain

Neuron production takes place continuously in the rostral migratory stream (RMS) of the adult mammalian brain. The molecular mechanisms that regulate progenitor cell division and differentiation in the RMS remain largely unknown. Here, we surveyed the mouse genome in an unbiased manner to identify candidate gene loci that regulate proliferation in the adult RMS. We quantified neurogenesis in adult C57BL/6J and A/J mice, and 27 recombinant inbred lines derived from those parental strains. We showed that the A/J RMS had greater numbers of bromodeoxyuridine-labeled cells than that of C57BL/6J mice with similar cell cycle parameters, indicating that the differences in the number of bromodeoxyuridine-positive cells reflected the number of proliferating cells between the strains. AXB and BXA recombinant inbred strains demonstrated even greater variation in the numbers of proliferating cells. Genome-wide mapping of this trait revealed that chromosome 11 harbors a significant quantitative trait locus at 116.75 +/- 0.75 Mb that affects cell proliferation in the adult RMS. The genomic regions that influence RMS proliferation did not overlap with genomic regions regulating proliferation in the adult subgranular zone of the hippocampal dentate gyrus. On the contrary, a different, suggestive locus that modulates cell proliferation in the subgranular zone was mapped to chromosome 3 at 102 +/- 7 Mb. A subset of genes in the chromosome 11 quantitative trait locus region is associated with neurogenesis and cell proliferation. Our findings provide new insights into the genetic control of neural proliferation and an excellent starting point to identify genes critical to this process.

A novel strategy for genetic dissection of complex traits: the population of specific chromosome substitution strains from laboratory and wild mice

The mouse is an irreplaceable model for understanding the precise genetic mechanisms of mammalian physiological pathways. Thousands of quantitative trait loci (QTLs) have been mapped onto the mouse genome during the last two decades. However, only a few genes' underlying complex traits have been successfully identified, and rapid fine mapping of QTL genes still remains a challenge for mouse geneticists. Currently, the Collaborative Cross (CC) has proceeded to the goal of establishing more than 1,000 recombinant inbred strains for the sub-centimorgan mapping resolution of QTL loci. In this article, a novel complementary strategy, designated as population of specific chromosome substitution strains or PSCSS, is proposed for rapid fine mapping of QTLs on the substituted chromosome. One specific chromosome (Chr 1) of recipient mouse strain C57BL/6 J has been substituted by homologous counterparts from five different inbred strains (C3H/He, FVB/N, AKR, NOD/LtJ, NZW/LacJ), an outbred line Kunmin mouse in China, and 23 wild mice captured in different localities. The primary genetic studies on the structure of these wild donor chromosomes (Chr 1) show that these donor chromosomes harbor extensive genetic polymorphisms, with a high density of SNP distribution, abundant variations of STR alleles, and a high level of historical recombination accumulation. These specific chromosome substitution strains eventually form a special population that has the identical genetic background of the recipient strain and differs only in the donor chromosomes. With simple association studies, known QTLs on the donor chromosome can be rapidly mapped in high resolution without requirement of further crosses. This approach, taking advantage of the extensive genetic polymorphisms of wild resources and chromosome substitution strategy, brings a new outlook for genetic dissection of complex traits.

Phenotype and Genetics of Progressive Sensorineural Hearing Loss (Snhl1) in the LXS Set of Recombinant Inbred Strains of Mice

Progressive sensorineural hearing loss is the most common form of acquired hearing impairment in the human population. It is also highly prevalent in inbred strains of mice, providing an experimental avenue to systematically map genetic risk factors and to dissect the molecular pathways that orchestrate hearing in peripheral sensory hair cells. Therefore, we ascertained hearing function in the inbred long sleep (ILS) and inbred short sleep (ISS) strains. Using auditory-evoked brain stem response (ABR) and distortion product otoacoustic emission (DPOAE) measurements, we found that ISS mice developed a high-frequency hearing loss at twelve weeks of age that progressed to lower frequencies by 26 weeks of age in the presence of normal endocochlear potentials and unremarkable inner ear histology. ILS mice exhibited milder hearing loss, showing elevated thresholds and reduced DPOAEs at the higher frequencies by 26 weeks of age. To map the genetic variants that underlie this hearing loss we computed ABR thresholds of 63 recombinant inbred stains derived from the ISS and ILS founder strains. A single locus was linked to markers associated with ISS alleles on chromosome 10 with a highly significant logarithm of odds (LOD) score of 15.8. The 2-LOD confidence interval spans ∼4 Megabases located at position 54–60 Mb. This locus, termed sensorineural hearing loss 1 (Snhl1), accounts for approximately 82% of the phenotypic variation. In summary, this study identifies a novel hearing loss locus on chromosome 10 and attests to the prevalence and genetic heterogeneity of progressive hearing loss in common mouse strains.

High-throughput sequencing of the DBA/2J mouse genome

BackgroundThe DBA/2J mouse is not only the oldest inbred strain,but also one of the most widely used strains. DBA/2Jexhibits many unique anatomical, physiological, andbehavior traits. In addition, DBA/2J is one parent of thelarge BXD family of recombinant inbred strains [1]. Thegenome of the other parent of this BXD family—C57BL/6J—has been sequenced and serves as themouse reference genome [2]. We sequenced the gen-ome of DBA/2J using SOLiD and Illumina highthroughput short read protocols to generate a compre-hensive set of ~5 million sequence variants segregatingin the BXD family that ultimately cause developmental,anatomical, functional and behavioral differences amongthese 80+ strains.ResultsWe generated approximately 13.2 and 38.9× whole-gen-ome short reads of DBA/2J females using Illumina GA2and ABI SOLiD massively parallel DNA sequencingplatforms. Comparing to the C57BL/6J reference gen-ome sequence, we identified over 4.5 million singlenucleotide polymorphisms (SNPs), including 84 non-sense and ~11,000 missense mutations, 78% of whichare novel. We also detected ~568,000 insertions anddeletions (indels) within single short reads and ~9,400between mate-paired reads. Approximately 300 inver-sions were detected by SOLiD mate-pair reads, 46 ofwhich span at least one exon. In addition, we identified~22,000 copy number variants (CNVs) in the range of1 Kb to 100 Kb (Figure 1).ConclusionOur study generates the first consensus sequence for theDBA/2J and creates a compendium of sequence andstructural variations that will be used by the communityof researchers who study complex traits in mouse mod-els.Thesequencedataprovideanovelresourcewithwhich to initiate reverse genetic analysis of complex

Correlation of Centromeric Heterochromatin C-band Polymorphism with Breeding Failure in Mice

The aim of the present study was to test the hypothesis about the relation between segregation of chromosomes 14 and 18 and the deterioration of mouse fertility and vitality. The analysis was possible because C-banding on chromosome 14 and chromosome 18 of the CBA/Kw and KE strains show size polymorphism. A small sized C-band on chromosome 14 is characteristic for the CBA/Kw mice, while the KE mice show small C-bands on chromosomes 18. Thus, if fertility parameters are affected in a centromere-dependent manner, we should observe non-random inheritance of both chromosome pairs in recombinant inbred (RI) strains. The results showed statistically significant preferential segregation of chromosomes 14 and 18 with small C-bands. Most of the RI strains inherited chromosome 14 from the CBA/Kw strain and chromosome 18 from the KE strain, and did not manifest a deterioration of fertility and vitality. On the contrary, RI strains that inherited chromosomes 14 and 18 from one of the parental strains, particularly the KE strain, stopped breeding or had difficulties in producing the next generation.

Genetically determined phenotype covariation networks control bone strength.

To identify genes affecting bone strength, we studied how genetic variants regulate components of a phenotypic covariation network that was previously shown to accurately characterize the compensatory trait interactions involved in functional adaptation during growth. Quantitative trait loci (QTLs) regulating femoral robustness, morphologic compensation, and mineralization (tissue quality) were mapped at three ages during growth using AXB/BXA Recombinant Inbred (RI) mouse strains and adult B6-iA Chromosome Substitution Strains (CSS). QTLs for robustness were identified on chromosomes 8, 12, 18, and 19 and confirmed at all three ages, indicating that genetic variants established robustness postnatally without further modification. A QTL for morphologic compensation, which was measured as the relationship between cortical area and body weight, was identified on chromosome 8. This QTL limited the amount of bone formed during growth and thus acted as a setpoint for diaphyseal bone mass. Additional QTLs were identified from the CSS analysis. QTLs for robustness and morphologic compensation regulated bone structure independently (ie, in a nonpleiotropic manner), indicating that each trait may be targeted separately to individualize treatments aiming to improve strength. Multiple regression analyses showed that variation in morphologic compensation and tissue quality, not bone size, determined femoral strength relative to body weight. Thus an individual inheriting slender bones will not necessarily inherit weak bones unless the individual also inherits a gene that impairs compensation. This systems genetic analysis showed that genetically determined phenotype covariation networks control bone strength, suggesting that incorporating functional adaptation into genetic analyses will advance our understanding of the genetic basis of bone strength. © 2010 American Society for Bone and Mineral Research.

Sex-specific gene expression in the BXD mouse liver.

Differences in clinical phenotypes between the sexes are well documented and have their roots in differential gene expression. While sex has a major effect on gene expression, transcription is also influenced by complex interactions between individual genetic variation and environmental stimuli. In this study, we sought to understand how genetic variation affects sex-related differences in liver gene expression by performing genetic mapping of genomewide liver mRNA expression data in a genetically defined population of naive male and female mice from C57BL/6J, DBA/2J, B6D2F1, and 37 C57BL/6J x DBA/2J (BXD) recombinant inbred strains. As expected, we found that many genes important to xenobiotic metabolism and other important pathways exhibit sexually dimorphic expression. We also performed gene expression quantitative trait locus mapping in this panel and report that the most significant loci that appear to regulate a larger number of genes than expected by chance are largely sex independent. Importantly, we found that the degree of correlation within gene expression networks differs substantially between the sexes. Finally, we compare our results to a recently released human liver gene expression data set and report on important similarities in sexually dimorphic liver gene expression between mouse and human. This study enhances our understanding of sex differences at the genome level and between species, as well as increasing our knowledge of the molecular underpinnings of sex differences in responses to xenobiotics.

Inferring combinatorial association logic networks in multimodal genome-wide screens

Motivation: We propose an efficient method to infer combinatorial association logic networks from multiple genome-wide measurements from the same sample. We demonstrate our method on a genetical genomics dataset, in which we search for Boolean combinations of multiple genetic loci that associate with transcript levels.Results: Our method provably finds the global solution and is very efficient with runtimes of up to four orders of magnitude faster than the exhaustive search. This enables permutation procedures for determining accurate false positive rates and allows selection of the most parsimonious model. When applied to transcript levels measured in myeloid cells from 24 genotyped recombinant inbred mouse strains, we discovered that nine gene clusters are putatively modulated by a logical combination of trait loci rather than a single locus. A literature survey supports and further elucidates one of these findings. Due to our approach, optimal solutions for multi-locus logic models and accurate estimates of the associated false discovery rates become feasible. Our algorithm, therefore, offers a valuable alternative to approaches employing complex, albeit suboptimal optimization strategies to identify complex models.Availability: The MATLAB code of the prototype implementation is available on: http://bioinformatics.tudelft.nl/ or http://bioinformatics.nki.nl/Contact: m.j.t.reinders@tudelft.nl; l.wessels@nki.nl

Susceptibility to Anthrax Lethal Toxin-Induced Rat Death Is Controlled by a Single Chromosome 10 Locus That Includes rNlrp1

Anthrax lethal toxin (LT) is a bipartite protease-containing toxin and a key virulence determinant of Bacillus anthracis. In mice, LT causes the rapid lysis of macrophages isolated from certain inbred strains, but the correlation between murine macrophage sensitivity and mouse strain susceptibility to toxin challenge is poor. In rats, LT induces a rapid death in as little as 37 minutes through unknown mechanisms. We used a recombinant inbred (RI) rat panel of 19 strains generated from LT-sensitive and LT-resistant progenitors to map LT sensitivity in rats to a locus on chromosome 10 that includes the inflammasome NOD-like receptor (NLR) sensor, Nlrp1. This gene is the closest rat homolog of mouse Nlrp1b, which was previously shown to control murine macrophage sensitivity to LT. An absolute correlation between in vitro macrophage sensitivity to LT-induced lysis and animal susceptibility to the toxin was found for the 19 RI strains and 12 additional rat strains. Sequencing Nlrp1 from these strains identified five polymorphic alleles. Polymorphisms within the N-terminal 100 amino acids of the Nlrp1 protein were perfectly correlated with LT sensitivity. These data suggest that toxin-mediated lethality in rats as well as macrophage sensitivity in this animal model are controlled by a single locus on chromosome 10 that is likely to be the inflammasome NLR sensor, Nlrp1.

Genetic dissection of dietary restriction in mice supports the metabolic efficiency model of life extension

Dietary restriction (DR) has been used for decades to retard aging in rodents, but its mechanism of action remains an enigma. A principal roadblock has been that DR affects many different processes, making it difficult to distinguish cause and effect. To address this problem, we applied a quantitative genetics approach utilizing the ILSXISS series of mouse recombinant inbred strains. Across 42 strains, mean female lifespan ranged from 380 to 1070 days on DR (fed 60% of ad libitum [AL]) and from 490 to 1020 days on an AL diet. Longevity under DR and AL is under genetic control, showing 34% and 36% heritability, respectively. There was no correlation between lifespans on DR and AL; thus different genes modulate longevity under the two regimens. DR lifespans are significantly correlated with female fertility after return to an AL diet after various periods of DR ( R = 0.44, P = 0.006). We assessed fuel efficiency (FE, ability to maintain growth and body weight independent of absolute food intake) using a multivariate approach and found it to be correlated with longevity and female fertility, suggesting possible causality. We found several quantitative trait loci responsible for these traits, mapping to chromosomes 7, 9, and 15. We present a metabolic model in which the anti-aging effects of DR are consistent with the ability to efficiently utilize dietary resources.

Calpain-10 is a component of the obesity-related quantitative trait locus Adip1

Calpain-10 is a component of the obesity-related quantitative trait locus Adip1

Genetic regulation of catecholamine synthesis, storage and secretion in the spontaneously hypertensive rat

Understanding catecholamine metabolism is crucial for elucidating the pathogenesis of hereditary hypertension. Here we integrated transcriptional and biochemical profiling with physiologic quantitative trait locus (eQTL and pQTL) mapping in adrenal glands of the HXB/BXH recombinant inbred (RI) strains, derived from the spontaneously hypertensive rat (SHR) and normotensive Brown Norway (BN.Lx). We found simultaneous down-regulation of five heritable transcripts in the catecholaminergic pathway in young (6 weeks) SHRs. We identified cis-acting eQTLs for Dbh, Pnmt (catecholamine biosynthesis) and Vamp1 (catecholamine secretion); enzymatic activities of Dbh and Pnmt paralleled transcripts, with pQTLs for activities mirroring eQTLs. We also detected trans-regulated expression of Vmat1 and Chga (both involved in catecholamine storage), with co-localization of these trans-eQTLs to the Pnmt locus. Pnmt re-sequencing revealed promoter polymorphisms that result in decreased response of the transfected SHR promoter to glucocorticoid, compared with BN.Lx. Of physiological pertinence, Dbh activity negatively correlated with systolic blood pressure in RI strains, whereas Pnmt activity was negatively correlated with heart rate. The finding of such cis- and trans-QTLs at an age before the onset of frank hypertension suggests that these heritable changes in biosynthetic enzyme expression represent primary genetic mechanisms for regulation of catecholamine action and blood pressure control in this widely studied model of hypertension.

Genomic Analysis of Variation in Hindlimb Musculature of Mice from the C57BL/6J and DBA/2J Lineage

The precise locations of attachment points of muscle to bone-the origin and insertion sites-are crucial anatomical and functional characteristics that influence locomotor performance. Mechanisms that control the development of these interactions between muscle, tendon, and bone are currently not well understood. In a subset of BXD recombinant inbred (RI) strains derived from the C57BL/6J and DBA/2J strains, we observed a soleus femoral attachment anomaly (SFAA) that was rare in both parental strains (Lionikas, Glover et al. 2006). The aim of the present study was to assess suitability of SFAA as a model to study the genetic mechanisms underlying variation in musculoskeletal anatomy. We scored the incidence of SFAA in 55 BXD strains (n = 9 to 136, median = 26, phenotyped animals per strain, for a total number of 2367). Seven strains (BXD1, 12, 38, 43, 48, 54, and 56) exhibited a high incidence of unilateral SFAA (47-89%), whereas 23 strains scored 0%. Exploration of the mechanisms underlying SFAA in 2 high incidence strains, BXD1 and BXD38, indicated that SFAA-relevant genes are to be found in both C57BL/6J and DBA/2J regions of the BXD1 genome. However, not all alleles relevant for the expression of the phenotype were shared between the 2 high-incidence BXD strains. In conclusion, the anatomical origin of the soleus muscle in mouse is controlled by a polygenic system. A panel of BXD RI strains is a useful tool in exploring the genetic mechanisms underlying SFAA and improving our understanding of musculoskeletal development.

High-throughput behavioral phenotyping in the expanded panel of BXD recombinant inbred strains.

Genetic reference populations, particularly the BXD recombinant inbred (BXD RI) strains derived from C57BL/6J and DBA/2J mice, are a valuable resource for the discovery of the bio-molecular substrates and genetic drivers responsible for trait variation and covariation. This approach can be profitably applied in the analysis of susceptibility and mechanisms of drug and alcohol use disorders for which many predisposing behaviors may predict the occurrence and manifestation of increased preference for these substances. Many of these traits are modeled by common mouse behavioral assays, facilitating the detection of patterns and sources of genetic coregulation of predisposing phenotypes and substance consumption. Members of the Tennessee Mouse Genome Consortium (TMGC) have obtained phenotype data from over 250 measures related to multiple behavioral assays across several batteries: response to, and withdrawal from cocaine, 3,4-methylenedioxymethamphetamine; “ecstasy” (MDMA), morphine and alcohol; novelty seeking; behavioral despair and related neurological phenomena; pain sensitivity; stress sensitivity; anxiety; hyperactivity and sleep/wake cycles. All traits have been measured in both sexes in approximately 70 strains of the recently expanded panel of BXD RI strains. Sex differences and heritability estimates were obtained for each trait, and a comparison of early (N = 32) and recent (N = 37) BXD RI lines was performed. Primary data are publicly available for heritability, sex difference and genetic analyses using the MouseTrack database, and are also available in GeneNetwork.org for quantitative trait locus (QTL) detection and genetic analysis of gene expression. Together with the results of related studies, these data form a public resource for integrative systems genetic analysis of neurobehavioral traits.

Systems genetics analyses predict a transcription role for P2P-R: Molecular confirmation that P2P-R is a transcriptional co-repressor

BackgroundThe 250 kDa P2P-R protein (also known as PACT and Rbbp6) was cloned over a decade ago and was found to bind both the p53 and Rb1 tumor suppressor proteins. In addition, P2P-R has been associated with multiple biological functions, such as mitosis, mRNA processing, translation and ubiquitination. In the current studies, the online GeneNetwork system was employed to further probe P2P-R biological functions. Molecular studies were then performed to confirm the GeneNetwork evaluations.ResultsGeneNetwork and associated gene ontology links were used to investigate the coexpression of P2P-R with distinct functional sets of genes in an adipocyte genetic reference panel of HXB/BXH recombinant strains of rats and an eye genetic reference panel of BXD recombinant inbred strains of mice. The results establish that biological networks of 75 and 135 transcription-associated gene products that include P2P-R are co-expressed in a genetically-defined manner in rat adipocytes and in the mouse eye, respectively. Of this large set of transcription-associated genes, >10% are associated with hormone-mediated transcription. Since it has been previously reported that P2P-R can bind the SRC-1 transcription co-regulatory factor (steroid receptor co-activator 1, [Ncoa1]), the possible effects of P2P-R on estrogen-induced transcription were evaluated. Estrogen-induced transcription was repressed 50-70% by the transient transfection of P2P-R plasmid constructs into four different cell types. In addition, knockdown of P2P-R expression using an antisense oligonucleotide increased estrogen-mediated transcription. Co-immunoprecipitation assays confirmed that P2P-R interacts with SRC-1 and also demonstrated that P2P-R interacts with estrogen receptor α.ConclusionsThe findings presented in this study provide strong support for the value of systems genetics, especially GeneNetwork, in discovering new functions of genes that can be confirmed by molecular analysis. More specifically, these data provide evidence that the expression of P2P-R co-varies in a genetically-defined manner with large transcription networks and that P2P-R can function as a co-repressor of estrogen-dependent transcription.

Identifying genetic loci and spleen gene coexpression networks underlying immunophenotypes in BXD recombinant inbred mice.

The immune system plays a pivotal role in the susceptibility to and progression of a variety of diseases. Due to a strong genetic basis, heritable differences in immune function may contribute to differential disease susceptibility between individuals. Genetic reference populations, such as the BXD (C57BL/6J × DBA/2J) panel of recombinant inbred (RI) mouse strains, provide unique models through which to integrate baseline phenotypes in healthy individuals with heritable risk for disease because of the ability to combine data collected from these populations across both multiple studies and time. We performed basic immunophenotyping (e.g., percentage of circulating B and T lymphocytes and CD4(+) and CD8(+) T cell subpopulations) in peripheral blood of healthy mice from 41 BXD RI strains to define the immunophenotypic variation in this strain panel and to characterize the genetic architecture that underlies these traits. Significant QTL models that explained the majority (50-77%) of phenotypic variance were derived for each trait and for the T:B cell and CD4(+):CD8(+) ratios. Combining QTL mapping with spleen gene expression data uncovered two quantitative trait transcripts, Ptprk and Acp1, as candidates for heritable differences in the relative abundance of helper and cytotoxic T cells. These data will be valuable in extracting genetic correlates of the immune system in the BXD panel. In addition, they will be a useful resource for prospective, phenotype-driven model selection to test hypotheses about differential disease or environmental susceptibility between individuals with baseline differences in the composition of the immune system.