Genetic architecture of voluntary exercise in an advanced intercross line of mice.

Exercise is essential for health, yet the amount, duration, and intensity that individuals engage in are strikingly variable, even under prescription. Our focus was to identify the locations and effects of quantitative trait loci (QTL) controlling genetic predisposition for exercise-related traits, utilizing a large advanced intercross line (AIL) of mice. This AIL (G(4)) population originated from a reciprocal cross between mice with genetic propensity for increased voluntary exercise [high-runner (HR) line, selectively bred for increased wheel running] and the inbred strain C57BL/6J. After adjusting for family structure, we detected 32 significant and 13 suggestive QTL representing both daily running traits (distance, duration, average speed, and maximum speed) and the mean of these traits on days 5 and 6 (the selection criteria for HR) of a 6-day test conducted at 8 wk of age, with many co-localizing to similar genomic regions. Additionally, seven significant and five suggestive QTL were observed for the slope and intercept of a linear regression across all 6 days of running, some representing a combination of the daily traits. We also observed two significant and two suggestive QTL for body mass before exercise. These results, from a well-defined animal model, reinforce a genetic basis for the predisposition to engage in voluntary exercise, dissect this predisposition into daily segments across a continuous time period, and present unique QTL that may provide insight into the initiation, continuation, and temporal pattern of voluntary activity in mammals.

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.

Using gene expression databases for classical trait QTL candidate gene discovery in the BXD recombinant inbred genetic reference population: mouse forebrain weight.

BACKGROUND: Successful strategies for QTL gene identification benefit from combined experimental and bioinformatic approaches. Unique design aspects of the BXD recombinant inbred line mapping panel allow use of archived gene microarray expression data to filter likely from unlikely candidates. This prompted us to propose a simple five-filter protocol for candidate nomination. To filter more likely from less likely candidates, we required candidate genes near to the QTL to have mRNA abundance that correlated with the phenotype among the BXD lines as well as differed between the parental lines C57BL/6J and DBA/2J. We also required verification of mRNA abundance by an independent method, and finally we required either differences in protein levels or confirmed DNA sequence differences. RESULTS: QTL mapping of mouse forebrain weight in 34 BXD RI lines found significant association on chromosomes 1 and 11, with each C57BL/6J allele increasing weight by more than half a standard deviation. The intersection of gene lists that were within +/- 10 Mb of the strongest associated location, that had forebrain mRNA abundance correlated with forebrain weight among the BXD, and that had forebrain mRNA abundance differing between C57BL/6J and DBA/2J, produced two candidates, Tnni1 (troponin 1) and Asb3 (ankyrin repeat and SOCS box-containing protein 3). Quantitative RT-PCR confirmed the direction of an increased expression in C57BL/6J genotype over the DBA/2J genotype for both genes, a difference that translated to a 2-fold difference in Asb3 protein. Although Tnni1 protein differences could not be confirmed, a 273 bp indel polymorphism was discovered 1 Kb upstream of the transcription start site. CONCLUSION: Delivery of well supported candidate genes following a single quantitative trait locus mapping experiment is difficult. However, by combining available gene expression data with QTL mapping, we illustrated a five-filter protocol that nominated Asb3 and Tnni1 as candidates affecting increased mouse forebrain weight. We recommend our approach when (1) investigators are working with phenotypic differences between C57BL/6J and DBA/2J, and (2) gene expression data are available on http://www.genenetwork.org that relate to the phenotype of interest. Under these circumstances, measurement of the phenotype in the BXD lines will likely also deliver excellent candidate genes.

Genome Reshuffling for Advanced Intercross Permutation (GRAIP): simulation and permutation for advanced intercross population analysis.

BACKGROUND: Advanced intercross lines (AIL) are segregating populations created using a multi-generation breeding protocol for fine mapping complex trait loci (QTL) in mice and other organisms. Applying QTL mapping methods for intercross and backcross populations, often followed by naïve permutation of individuals and phenotypes, does not account for the effect of AIL family structure in which final generations have been expanded and leads to inappropriately low significance thresholds. The critical problem with naïve mapping approaches in AIL populations is that the individual is not an exchangeable unit. METHODOLOGY/PRINCIPAL FINDINGS: The effect of family structure has immediate implications for the optimal AIL creation (many crosses, few animals per cross, and population expansion before the final generation) and we discuss these and the utility of AIL populations for QTL fine mapping. We also describe Genome Reshuffling for Advanced Intercross Permutation, (GRAIP) a method for analyzing AIL data that accounts for family structure. GRAIP permutes a more interchangeable unit in the final generation crosses - the parental genome - and simulating regeneration of a permuted AIL population based on exchanged parental identities. GRAIP determines appropriate genome-wide significance thresholds and locus-specific P-values for AILs and other populations with similar family structures. We contrast GRAIP with naïve permutation using a large densely genotyped mouse AIL population (1333 individuals from 32 crosses). A naïve permutation using coat color as a model phenotype demonstrates high false-positive locus identification and uncertain significance levels, which are corrected using GRAIP. GRAIP also detects an established hippocampus weight locus and a new locus, Hipp9a. CONCLUSIONS AND SIGNIFICANCE: GRAIP determines appropriate genome-wide significance thresholds and locus-specific P-values for AILs and other populations with similar family structures. The effect of family structure has immediate implications for the optimal AIL creation and we discuss these and the utility of AIL populations.

A simple method for combining genetic mapping data from multiple crosses and experimental designs.

BACKGROUND: Over the past decade many linkage studies have defined chromosomal intervals containing polymorphisms that modulate a variety of traits. Many phenotypes are now associated with enough mapping data that meta-analysis could help refine locations of known QTLs and detect many novel QTLs. METHODOLOGY/PRINCIPAL FINDINGS: We describe a simple approach to combining QTL mapping results for multiple studies and demonstrate its utility using two hippocampus weight loci. Using data taken from two populations, a recombinant inbred strain set and an advanced intercross population we demonstrate considerable improvements in significance and resolution for both loci. 1-LOD support intervals were improved 51% for Hipp1a and 37% for Hipp9a. We first generate locus-wise permuted P-values for association with the phenotype from multiple maps, which can be done using a permutation method appropriate to each population. These results are then assigned to defined physical positions by interpolation between markers with known physical and genetic positions. We then use Fisher's combination test to combine position-by-position probabilities among experiments. Finally, we calculate genome-wide combined P-values by generating locus-specific P-values for each permuted map for each experiment. These permuted maps are then sampled with replacement and combined. The distribution of best locus-specific P-values for each combined map is the null distribution of genome-wide adjusted P-values. CONCLUSIONS/SIGNIFICANCE: Our approach is applicable to a wide variety of segregating and non-segregating mapping populations, facilitates rapid refinement of physical QTL position, is complementary to other QTL fine mapping methods, and provides an appropriate genome-wide criterion of significance for combined mapping results.

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.

Combining gene expression QTL mapping and phenotypic spectrum analysis to uncover gene regulatory relationships.

Gene expression QTL (eQTL) mapping can suggest candidate regulatory relationships between genes. Recent advances in mammalian phenotype annotation such as mammalian phenotype ontology (MPO) enable systematic analysis of the phenotypic spectrum subserved by many genes. In this study we combined eQTL mapping and phenotypic spectrum analysis to predict gene regulatory relationships. Five pairs of genes with similar phenotypic effects and potential regulatory relationships suggested by eQTL mapping were identified. Lines of evidence supporting some of the predicted regulatory relationships were obtained from biological literature. A particularly notable example is that promoter sequence analysis and real-time PCR assays support the predicted regulation of protein kinase C epsilon (Prkce) by cAMP responsive element binding protein 1 (Creb1). Our results show that the combination of gene eQTL mapping and phenotypic spectrum analysis may provide a valuable approach to uncovering gene regulatory relations underlying mammalian phenotypes.

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.

Genetic analysis of barrel field size in the first somatosensory area (SI) in inbred and recombinant inbred strains of mice.

We measured the combined area of posterior medial barrel subfield (PMBSF) and anterior lateral barrel subfield (ALBSF) areas in four common inbred strains (C3H/HeJ, A /J, C57BL /6J, DBA/2J), B6D2F1, and ten recombinant inbred (RI) strains generated from C57BL/6J and DBA/2J progenitors (BXD) as an initial attempt to examine the genetic influences underlying natural variation in barrel field size in adult mice. These two subfields are associated with the representation of the whisker pad and sinus hairs on the contralateral face. Using cytochrome oxidase labeling to visualize the barrel field, we measured the size of the combined subfields in each mouse strain. We also measured body weight and brain weight in each strain. We report that DBA/2J mice have a larger combined PMBSF/ALBSF area (6.15 +/- 0.10 mm(2), n = 7) than C57BL /6J (5.48 +/- 0.13 mm(2), n = 10), C3H/HeJ (5.37 +/- 0.16 mm(2), n = 10), and A/J mice (5.04 +/- 0.09 mm(2), n = 15), despite the fact that DBA/2J mice have smaller average brain and body sizes. This finding may reflect dissociation between systems that control brain size with those that regulate barrel field area. In addition, BXD strains (average n = 4) and parental strains showed considerable and continuous variation in PMBSF/ALBSF area, suggesting that this trait is polygenic. Furthermore, brain, body, and cortex weights have heritable differences between inbred strains and among BXD strains. PMBSF/ALBSF pattern appears similar among inbred and BXD strains, suggesting that somatosensory patterning reflects a common plan of organization. This data is an important first step in the quantitative genetic analysis of the parcellation of neocortex into diverse cytoarchitectonic zones that vary widely within and between species, and in identifying the genetic factors underlying barrel field size using quantitative trait locus (QTL) analyses.

A new set of BXD recombinant inbred lines from advanced intercross populations in mice.

BACKGROUND: Recombinant inbred (RI) strains are an important resource for mapping complex traits in many species. While large RI panels are available for Arabidopsis, maize, C. elegans, and Drosophila, mouse RI panels typically consist of fewer than 30 lines. This is a severe constraint on the power and precision of mapping efforts and greatly hampers analysis of epistatic interactions. RESULTS: In order to address these limitations and to provide the community with a more effective collaborative RI mapping panel we generated new BXD RI strains from two independent advanced intercrosses (AI) between C57BL/6J (B6) and DBA/2J (D2) progenitor strains. Progeny were intercrossed for 9 to 14 generations before initiating inbreeding, which is still ongoing for some strains. Since this AI base population is highly recombinant, the 46 advanced recombinant inbred (ARI) strains incorporate approximately twice as many recombinations as standard RI strains, a fraction of which are inevitably shared by descent. When combined with the existing BXD RI strains, the merged BXD strain set triples the number of previously available unique recombinations and quadruples the total number of recombinations in the BXD background. CONCLUSION: The combined BXD strain set is the largest mouse RI mapping panel. It is a powerful tool for collaborative analysis of quantitative traits and gene function that will be especially useful to study variation in transcriptome and proteome data sets under multiple environments. Additional strains also extend the value of the extensive phenotypic characterization of the previously available strains. A final advantage of expanding the BXD strain set is that both progenitors have been sequenced, and approximately 1.8 million SNPs have been characterized. This provides unprecedented power in screening candidate genes and can reduce the effective length of QTL intervals. It also makes it possible to reverse standard mapping strategies and to explore downstream effects of known sequence variants.