Polyendocrine metabolic ovarian syndrome (PMOS), formerly known as polycystic ovary syndrome (PCOS), is the most common endocrine disorder in women and is closely associated with complex diseases such as cardiovascular disease and type 2 diabetes. However, the mechanistic links between PMOS and its comorbidities remain poorly understood. Here, we present an integrative systems genetics platform that leverages genetic diversity in both mice and humans to dissect the drivers of PMOS and its associated complications. This framework uncovers conserved genetic and environmental factors underlying PMOS, identifies susceptible cell types and organs, and elucidates mechanisms linking PMOS to subsequent pathologies. For instance, we show that increased ovarian area contributes to both PMOS susceptibility and ovarian cancer progression, while specific ovary-heart signaling circuits modulate cardiac function with aging. We further identify ovarian SF3B1-mediated alternative splicing as a key mechanistic link between PMOS and metabolic traits. Pharmacologic inhibition of SF3B1 in mice reduced circulating testosterone, insulin and glucose levels, as well as fat mass expansion. Transcriptomics analysis of ovaries from mice and experiments using human cell lines localized these effects to exon skipping events in granulosa cells. Together, this study offers a mechanistic framework for modeling the diversity of PMOS pathologies and uncovers SF3B1-mediated splicing as a link between ovary function and systemic metabolism.

ABSTRACT Developmental coordination disorder (DCD) is a common neurodevelopmental condition characterized by impaired motor coordination and learning, yet its neurobiological and genetic bases remain poorly understood. Here, we leverage the BXD recombinant inbred mouse panel to model the polygenic architecture of DCD and link behaviour, brain structure, and genotype. High-resolution ex vivo MRI across 14 strains revealed that DCD-like mice have modestly reduced total brain volume, with a distinct neuroanatomical profile characterized by enlarged cortical regions alongside reduced cerebellar, thalamic, and other subcortical volumes. These structural differences closely mirror findings reported in human DCD. Across strains, variation in brain structure strongly correlated with motor behaviours, with coordinated patterns linking increased cortical and decreased subcortical volumes to poorer motor coordination, while more focal associations were observed for motor learning. Multivariate analysis identified a dominant brain-behaviour axis capturing this cortical-subcortical trade-off. Quantitative trait locus (QTL) mapping revealed multiple loci influencing regional brain volumes, including a prominent locus on chromosome 12 regulating cerebellar structures, but did not identify single loci driving the main multivariate brain-behaviour relationships, consistent with a distributed genetic architecture. Together, these findings demonstrate that DCD-like motor impairments arise from coordinated alterations across distributed brain systems under polygenic control. This work establishes a translational framework linking genetic variation to brain organization and motor function, and suggests that DCD reflects the extreme of a continuous spectrum of neurobiological variation rather than a discrete condition.

Exploring the molecular genetic cascades responsible for behavioral responses to opioids can improve our understanding of drug use initiation. We generated high-precision time-series data for 105 morphine- and naloxone-related traits across ∼700 young adult BXD mice (64 diverse strains and both sexes) for 3 hours after a single morphine injection. Variations in responses were mapped using high precision sequencing-based genotypes. The initial locomotor responses to morphine map precisely to the µ opioid receptor gene (MOR or Oprm1) on chromosome (Chr) 10 with a peak linkage of 12.4 (–log10P). The B allele inherited from C57BL/6J is associated with up to 60% higher activity. This effect climaxes at 75 min but is exhausted by 160 min. A second major modulator of opioid activation emerges after about 100 min and is located on Chr 16 with peak linkages of 10.6 (–log10P) in females, also associated with a high B allele. This locus includes only one compelling candidate—fibroblast growth factor 12 (Fgf12), a 600 Kb gene that controls sodium current kinetics at the axon hillock. A strong and transient epistatic interaction exists between the Oprm1 and Fgf12 loci during a short time window (45–75 min). The combination of a B haplotype at Oprm1 with a D haplotype from DBA/2J at Fgf12 is associated with unusually high activity. In a complementary study in heterogeneous stock rats we demonstrate that Oprm1 and Fgf12 are co-expressed in one specific subtype of Drd1+ medium spiny neuron. A Bayesian network analysis supports an Oprm1-to-Fgf12 network that involves a MAP kinase cascade—Mapk8ip2, Map3k11, and Map3k12—that we hypothesize modulates FGF12 phosphorylation, Nav1.2 sodium channel state, and locomotor activation. OPRM1 and FGF12 networks in human GWAS data highlight enrichment of signals associated with substance use disorder. This is the first demonstration of a time-dependent epistatic interaction modulating drug response in mammals and the first linkage of Fgf12 to opioid sensitivity and potentially to sodium channel activity.

DNA variants modulate mortality risks across an entire lifespan but their dynamic age-dependent effects have not been resolved in any species for either sex. Here we mapped variants that shape mortality using an actuarial approach, starting with a base population of 6,438 pubescent mice and ending with 559 survivors that lived beyond 1,100 days of age. Twenty-nine Vita loci influence lifespan with strong age- and sex-specific effects. Most act during distinct stages with polarities that often invert with age, but a minority have consistent age-dependent effects in one or both sexes. A separate set of 30 Soma loci influence correlations between body mass and life expectancy. Nineteen Soma loci mediate higher mortality in larger young mice, whereas 11 mediate lower mortality in larger old mice. All effects are stronger in male mice than in female mice. Vita and Soma loci form epistatic networks split strictly by sex. These findings provide a genetic bridge between evolutionary theories of ageing and molecular mechanisms that can guide interventions to extend healthy lifespan.

Alzheimer's disease has widespread effects on brain structure, function and behavior, but we lack a systematic dissection of its impact across hundreds of forebrain and brainstem regions. Here, using diffusion tensor MRI at 25 µm, we mapped the global consequences of mutations in APP and PSEN1 across 231 regions of interest (ROIs) in male and female 5×FAD BXD hybrid mice at 14 months. Over half of the ROIs change in volume along rostrocaudal and mediolateral axes of the CNS, with unexpected swelling in the neocortex, hippocampus and amygdala of up to 10%, counterbalanced by shrinkage in the thalamus, brainstem and most white matter tracts. Yet, total brain volume is unaltered. Variation in individual ROI volumes is highest in females. Differences in fear acquisition and contextual memory performance covary with volumes of several regions and can have opposite polarities between cases and controls. These structural benchmarks establish a foundation for testing therapeutic interventions in preclinical trials.

The fundamental skills for motor coordination and motor control emerge through development. Neurodevelopmental disorders such as developmental coordination disorder (DCD) lead to impaired acquisition of motor skills. This study investigated motor behaviors that reflect the core symptoms of human DCD through the use of BXD recombinant inbred strains of mice that are known to have divergent phenotypes in many behavioral traits, including motor activity. We sought to correlate behavior in basic motor control tasks with the known genotypes of these reference populations of mice using quantitative trait locus (QTL) mapping. We used 12 BXD strains with an average of 16 mice per group to assess the onset of reflexes during the early neonatal stage of life and differences in motor coordination using the tests for open field, rotarod, and gait behaviors during the adolescent/young adulthood period. Results indicated significant variability between strains in when neonatal reflexes appeared and significant strain differences for all measures of motor coordination. Five strains (BXD15, BXD27, BXD28, BXD75, BXD86) struggled with sensorimotor coordination as seen in gait analysis, rotarod, and open field, similar to human presentation of DCD. We identified three significant quantitative trait loci for gait on proximal Chr 3, Chr 4, and distal Chr 6. Based on expression, function, and polymorphism within the mapped QTL intervals, seven candidate genes (Gpr63, Spata5, Trpc3, Cntn6, Chl1, Grm7, Ogg1) emerged. This study offers new insights into mouse motor behavior, which promises to be a first murine model to explore the genetics and neural correlates of DCD.

Alzheimer's disease (AD) affects brain regions with remarkable heterogeneity, but the precise impact of this disease on hundreds of small cortical, subcortical, and brainstem regions remains poorly defined. Here, as a prelude to testing preclinical models to prevent AD, we systematically quantified effects of human AD mutations in APP and PSEN1 on 231 regions and comprehensively evaluated changes in volume with unprecedented resolution in genetically diverse mice as a function of sex and genetic background. We studied 34 5XFAD F1 hybrids and 23 sibling controls at 14 months, evaluating learning and memory behaviors, followed by ex vivo diffusion tensor images (DTI) at 25 micron resolution. We delineated 231 regions of interest (ROIs) bilaterally with high precision. Remarkably, we found bidirectional changes: marked volume increases (up to 10%) in neocortex, hippocampus, amygdala, and sensory nuclei, contrasted with decreases in striatum, pallidum, thalamus, hypothalamus, and most fiber tracts. These opposing effects are unrelated to amyloid load and are likely to reflect temporal gradients in susceptibility of ROIs. Effects are similar in both sexes but far more prominent in females. Genetic background strongly modulates penetrance of the human mutations with the AD-BXD77 F1 type having the greatest sensitivity. Light sheet microscopy and stereological analysis of NeuN+ neurons and amyloid beta; aggregates in 22 regions revealed up to 20% loss of cells in CA3 and in anterior and intralaminar parts of the thalamus. Volumetric changes correlated with impaired fear acquisition and memory, with cases and controls often showing opposite relations between performance and regional volumes. These findings reveal unprecedented regional heterogeneity in AD progression and suggest therapeutic efficacy may vary substantially across genetic backgrounds and between sexes. ### Competing Interest Statement The authors have declared no competing interest.

The dynamics of lifespan are shaped by DNA variants that exert effects at different ages. We have mapped genetic loci that modulate age-specific mortality using an actuarial approach. We started with an initial population of 6,438 pubescent siblings and ended with a survivorship of 559 mice that lived to at least 1100 days. Twenty-nine Vita loci dynamically modulate the mean lifespan of survivorships with strong age- and sex-specific effects. Fourteen have relatively steady effects on mortality while other loci act forcefully only early or late in life and with polarities of effects that invert. A distinct set of 19 Soma loci shape the negative correlation between weights of young adults with their life expectancies-much more strongly so in males than females. Another set of 11 Soma loci shape the positive correlation between weights at older ages with life expectancies. The Vita and Soma loci share 289 age-dependent epistatic interactions (LODs ≥3.8) but fewer than 4% are common to both sexes. We provide two examples of how to move from maps toward potential mechanisms. Our findings provide an empirical bridge between evolutionary theories on aging and genetic and molecular causes. These loci and their interactions are key to begin to understand the impact of interventions that may extend healthy lifespan in mice and even in humans.

Hundreds of inbred mouse strains and intercross populations have been used to characterize the function of genetic variants that contribute to disease. Thousands of disease-relevant traits have been characterized in mice and made publicly available. New strains and populations including consomics, the collaborative cross, expanded BXD, and inbred wild-derived strains add to existing complex disease mouse models, mapping populations, and sensitized backgrounds for engineered mutations. The genome sequences of inbred strains, along with dense genotypes from others, enable integrated analysis of trait-variant associations across populations, but these analyses are hampered by the sparsity of genotypes available. Moreover, the data are not readily interoperable with other resources. To address these limitations, we created a uniformly dense variant resource by harmonizing multiple data sets. Missing genotypes were imputed using the Viterbi algorithm with a data-driven technique that incorporates local phylogenetic information, an approach that is extendable to other model organisms. The result is a web- and programmatically accessible data service called GenomeMUSter, comprising single-nucleotide variants covering 657 strains at 106.8 million segregating sites. Interoperation with phenotype databases, analytic tools, and other resources enable a wealth of applications, including multitrait, multipopulation meta-analysis. We show this in cross-species comparisons of type 2 diabetes and substance use disorder meta-analyses, leveraging mouse data to characterize the likely role of human variant effects in disease. Other applications include refinement of mapped loci and prioritization of strain backgrounds for disease modeling to further unlock extant mouse diversity for genetic and genomic studies in health and disease.

We require experimental populations in which to test precision medicine interventions before they can be translated to humans. Mouse models can provide immortal genometypes that can be resampled at any stage, in either sex, and under any experimental conditions to model genome-by-environment interactions and test genome-phenome predictions. This allows for experimental precision medicine—systems in which causal models of complex interactions among DNA variants, phenotype variation at many levels, and innumerable environmental factors can be built. In this chapter, I review some of the challenges of precision medicine and how mouse populations can overcome these challenges. Lastly, I provide specific suggestions and tools for future development of experimental precision medicine.

Abstract Disclosure: L.M. Velez: None. C. Johnson: None. I. Tamburrini: None. M. Zhou: None. C. Viesi: None. N. Ujagar: None. D. Ashbrook: None. M. Nelson: None. A. Senior: None. D. James: None. R. Williams: None. D. Nicholas: None. M. Seldin: None. Polycystic ovary syndrome (PCOS) is the most common endocrinopathy in women, with a prevalence of ∼4-20% in women of reproductive age. The diagnosis of the syndrome generally occurs when the patient consults for fertility issues and is only based on a reproductive criterion, which includes (1) hyperandrogenism, (2) oligo-anovulation, and (3) polycystic ovary morphology. However, the overlap of PCOS with cardiometabolic diseases is significant. To put in numbers; up to 75% of PCOS women present some degree of insulin insensitivity, 38-88% present obesity or overweight, 20-50% develop type 2 diabetes by age 40, and PCOS women are at increased risk of cardiovascular disease. Despite these facts, shared reproductive/metabolic mechanisms are largely underexplored. Moreover, studies addressing the genetic architecture of PCOS are missing. Here, we induced a PCOS-like condition in 25 recombinant and classical inbred female strains and matched placebo controls over 6 weeks. Comprehensive in vivo and terminal reproductive/metabolic analyses were performed, as well as ovary and adipose RNA-Seq. These strains varied in PCOS response in a number of key metabolic and reproductive traits, including circulating hormone levels, glucose metabolism, and cardiac function. We applied a linear mixed-effects model to estimate heritability, and genetic (h2), PCOS, and gene-by-PCOS interactions. High h2 was observed for lean and fat mass, glucose, and AUC, whereas PCOS effects were high for the BW change, testosterone, and AUC. Substantial gene-by-PCOS interactions were found for reproductive hormones. Undirected network construction and centrality estimates showed that the reproductive hormones LH and LH/FSH ratio were the strongest central traits connecting metabolic phenotypes. We also showed select strains represent subtypes of human PCOS-metabolism interaction with varied susceptibilities to disease in a PCOS setting. Ovarian RNA-seq analysis of PCOS DEGs showed strong enrichments with human disease settings such as hyperandrogenism, inflammation, and pregnancy hypertension. Similar analyses in GWAT RNA-seq showed enrichments in weight gain, liposarcoma, inflammation, and reproductive diseases were at the top, with adipose genes connecting these diseases and potentially involved with PCOS. In conclusion, we established a PCOS model to study relevant mechanisms intersecting reproduction with metabolism in the context of genetic variation. Presentation: 6/2/2024

In aged humans and mice, aggregates of hypobranched glycogen molecules called polyglucosan bodies (PGBs) accumulate in hippocampal astrocytes. PGBs are known to drive cognitive decline in neurological diseases but remain largely unstudied in the context of typical brain aging. Here, we show that PGBs arise in autophagy-dysregulated astrocytes of the aged C57BL/6J mouse hippocampus. To map the genetic cause of age-related PGB accumulation, we quantified PGB burden in 32 fully sequenced BXD-recombinant inbred mouse strains, which display a 400-fold variation in hippocampal PGB burden at 16-18 months of age. A major modifier locus was mapped to chromosome 1 at 72-75 Mb, which we defined as the Pgb1 locus. To evaluate candidate genes and downstream mechanisms by which Pgb1 controls the aggregation of glycogen, extensive hippocampal transcriptomic and proteomic datasets were produced for aged mice of the BXD family. We utilized these datasets to identify Smarcal1 and Usp37 as potential regulators of PGB accumulation. To assess the effect of PGB burden on age-related cognitive decline, we performed phenome-wide association scans, transcriptomic analyses as well as conditioned fear memory and Y-maze testing. Importantly, we did not find any evidence suggesting a negative impact of PGBs on cognition. Taken together, our study demonstrates that the Pgb1 locus controls glycogen aggregation in astrocytes of the aged hippocampus without affecting age-related cognitive decline.

We created GNQA, a generative pre-trained transformer (GPT) knowledge base driven by a performant retrieval augmented generation (RAG) with a focus on aging, dementia, Alzheimer’s and diabetes. We uploaded a corpus of three thousand peer reviewed publications on these topics into the RAG. To address concerns about inaccurate responses and GPT ‘hallucinations’, we implemented a context provenance tracking mechanism that enables researchers to validate responses against the original material and to get references to the original papers. To assess the effectiveness of contextual information we collected evaluations and feedback from both domain expert users and ‘citizen scientists’ on the relevance of GPT responses.A key innovation of our study is automated evaluation by way of a RAG assessment system (RAGAS). RAGAS combines human expert assessment with AI-driven evaluation to measure the effectiveness of RAG systems. When evaluating the responses to their questions, human respondents give a “thumbs-up” 76% of the time. Meanwhile, RAGAS scores 90% on answer relevance on questions posed by experts. And when GPT-generates questions, RAGAS scores 74% on answer relevance. With RAGAS we created a benchmark that can be used to continuously assess the performance of our knowledge base.Full GNQA functionality is embedded in the freeGeneNetwork.orgweb service, an open-source system containing over 25 years of experimental data on model organisms and human. The code developed for this study is published under a free and open-source software license athttps://git.genenetwork.org/gn-ai/tree/README.md.

Pangenome graphs can represent all variation between multiple reference genomes, but current approaches to build them exclude complex sequences or are based upon a single reference. In response, we developed the PanGenome Graph Builder, a pipeline for constructing pangenome graphs without bias or exclusion. The PanGenome Graph Builder uses all-to-all alignments to build a variation graph in which we can identify variation, measure conservation, detect recombination events and infer phylogenetic relationships. PGGB is a modular framework for efficiently building unbiased pangenome graphs, supporting diverse downstream analyses.

Breast cancer (BC) is the most common cancer and the second cause of death in US women. Our lack of understanding of how genetic variants affect molecular mechanisms that mediate BC aggression poses a substantial obstacle to advancements in cancer diagnosis and therapy. To examine genetic variants on BC traits, a novel murine model was created with robust phenotypic and genomic variation. The FVB C3(1)-T-antigen (“C3Tag”) mouse develops spontaneous tumors in the mammary glands of female mice with a mean latency of 4-5 months of age. This genetically engineered mouse model (GEMM) is well established to resemble human basal-like TNBC. TNBC is an aggressive subtype with few clinical approaches and poor patient outcomes. Thus, to model human heterogeneity in BC outcomes, we systematically crossed the C3Tag GEMM into the BXD recombinant inbred family – the largest and best characterized genetic reference population. The new model is termed “BXD-BC” and F1 hybrids of the cross have isogenic genomes that are reproducible. BXD-BCs are a potent tool to determine the impact of genetic modifiers on BC tumor traits. We hypothesized that examination of BXD-BC GEMMs will enable the identification of susceptibility loci, candidate genes, and molecular networks that underlie variation of multiple BC phenotypes. Using N=29 BXD-BC strains, we demonstrated significant heritable variations in the severity of TNBC characteristics such as tumor latency, multiplicity, and survival. Interestingly, 2 BXD-BC strains never developed tumors out to 1 year of age. Thus, BXD-BC strains demonstrate variance in cancer susceptibility and progression compared to the parent C3Tag GEMM, indicating the presence of genetic modifiers. Through an unbiased systematic quantification of breast cancer severity across BXD-BC hybrids, we identified several significant quantitative trait loci (QTL) and candidate genes for specific tumor traits. In combination with public human GWAS datasets, we defined syntenic regions, candidate genes, and underlying networks through cross-species systems genetics analyses to demonstrate the translational validity of conserved, biologically relevant, and targetable candidates. Our findings suggest conserved candidates predicting TNBC patient survival. In sum, the BXD-BC resource is an innovative, reliable, and robust preclinical model that reflects robust genetic heterogeneity. Using cutting edge systems genetics, we have identified genetic modifiers of BC phenotypic variation that could be targeted to advance therapeutic limitations or as biomarkers of risk or response to therapy.### Competing Interest StatementThe authors have declared no competing interest.

AbstractGenetic differences among mammalian hosts andMycobacterium tuberculosis(Mtb) strains determine diverse tuberculosis (TB) patient outcomes. The advent of recombinant inbred mouse panels and next-generation transposon mutagenesis and sequencing approaches has enabled dissection of complex host- pathogen interactions. To identify host and pathogen genetic determinants ofMtbpathogenesis, we infected members of the BXD family of mouse strains with a comprehensive library ofMtbtransposon mutants (TnSeq). Members of the BXD family segregate forMtb-resistant C57BL/6J (B6 orB) andMtb-susceptible DBA/2J (D2 orD) haplotypes. The survival of each bacterial mutant was quantified within each BXD host, and we identified those bacterial genes that were differentially required forMtbfitness across BXD genotypes. Mutants that varied in survival among the host family of strains were leveraged as reporters for “endophenotypes”, each bacterial fitness profile directly probing specific components of the infection microenvironment. We conducted QTL mapping of these bacterial fitness endophenotypes and identified 140host-pathogenquantitative trait loci (hpQTL). We identified a QTL hotspot on chromosome 6 (75.97–88.58 Mb) associated with the genetic requirement of multipleMtbgenes;Rv0127(mak),Rv0359(rip2),Rv0955(perM), andRv3849(espR). Together, this screen reinforces the utility of bacterial mutant libraries as precise reporters of the host immunological microenvironment during infection and highlights specific host-pathogen genetic interactions for further investigation. To enable downstream follow-up for both bacterial and mammalian genetic research communities, all bacterial fitness profiles have been deposited into GeneNetwork.org and added into the comprehensive collection of TnSeq libraries in MtbTnDB.

ABSTRACT Hundreds of inbred laboratory mouse strains and intercross populations have been used to functionalize genetic variants that contribute to disease. Thousands of disease relevant traits have been characterized in mice and made publicly available. New strains and populations including the Collaborative Cross, expanded BXD and inbred wild-derived strains add to set of complex disease mouse models, genetic mapping resources and sensitized backgrounds against which to evaluate engineered mutations. The genome sequences of many inbred strains, along with dense genotypes from others could allow integrated analysis of trait – variant associations across populations, but these analyses are not feasible due to the sparsity of genotypes available. Moreover, the data are not readily interoperable with other resources. To address these limitations, we created a uniformly dense data resource by harmonizing multiple variant datasets. Missing genotypes were imputed using the Viterbi algorithm with a data-driven technique that incorporates local phylogenetic information, an approach that is extensible to other model organism species. The result is a web– and programmatically-accessible data service called GenomeMUSter ( https://muster.jax.org ), comprising allelic data covering 657 strains at 106.8M segregating sites. Interoperation with phenotype databases, analytic tools and other resources enable a wealth of applications including multi-trait, multi-population meta-analysis. We demonstrate this in a cross-species comparison of the meta-analysis of Type 2 Diabetes and of substance use disorders, resulting in the more specific characterization of the role of human variant effects in light of mouse phenotype data. Other applications include refinement of mapped loci and prioritization of strain backgrounds for disease modeling to further unlock extant mouse diversity for genetic and genomic studies in health and disease.

Of the nearly 1 million military personnel who participated in the 1990-1991 Gulf War, between 25% and 35% became ill with what now is referred to as Gulf War Illness (GWI) by the Department of Defense. Symptoms varied from gastrointestinal distress to lethargy, memory loss, inability to concentrate, depression, respiratory, and reproductive problems. The symptoms have persisted for 30 years in those afflicted but the basis of the illness remains largely unknown. Nerve agents and other chemical exposures in the war zone have been implicated but the long-term effects of these acute exposures have left few if any identifiable signatures. The major aim of this study is to elucidate the possible genomic basis for the persistence of symptoms, especially of the neurological and behavioral effects. To address this, we performed a whole genome epigenetic analysis of the proposed cause of GWI, viz., exposure to organophosphate neurotoxicants combined with high circulating glucocorticoids in two inbred mouse strains, C57BL/6J and DBA/2J. The animals received corticosterone in their drinking water for 7 days followed by injection of diisopropylfluorophosphate, a nerve agent surrogate. Six weeks after DFP injection, the animals were euthanized and medial prefrontal cortex harvested for genome-wide DNA methylation analysis using high-throughput sequencing. We observed 67 differentially methylated genes, notably among them, Ttll7, Akr1c14, Slc44a4, and Rusc2, all related to different symptoms of GWI. Our results support proof of principle of genetic differences in the chronic effects of GWI-related exposures and may reveal why the disease has persisted in many of the now aging Gulf War veterans.

Developmental Coordination Disorder (DCD) is a neurodevelopmental disorder of unknown etiology that affects one in 20 children. There is an indication that DCD has an underlying genetic component due to its high heritability. Therefore, we explored the use of a recombinant inbred family of mice known as the BXD panel to understand the genetic basis of complex traits (i.e., motor learning) through identification of quantitative trait loci (QTLs). The overall aim of this study was to utilize the QTL approach to evaluate the genome-to-phenome correlation in BXD strains of mice in order to better understand the human presentation of DCD. Results of this current study confirm differences in motor learning in selected BXD strains and strains with altered cerebellar volume. Five strains - BXD15, BXD27, BXD28, BXD75, and BXD86 - exhibited the most DCD-like phenotype when compared with other BXD strains of interest. Results indicate that BXD15 and BXD75 struggled primarily with gross motor skills, BXD28 primarily had difficulties with fine motor skills, and BXD27 and BXD86 strains struggled with both fine and gross motor skills. The functional roles of genes within significant QTLs were assessed in relation to DCD-like behavior. Only Rab3a (Ras-related protein Rab-3A) emerged as a high likelihood candidate gene for the horizontal ladder rung task. This gene is associated with brain and skeletal muscle development, but lacked nonsynonymous polymorphisms. This study along with Gill et al. (same issue) is the first studies to specifically examine the genetic linkage of DCD using BXD strains of mice.