Quantitative trait loci mapping and gene network analysis implicate protocadherin‐15 as a determinant of brain serotonin transporter expression

Presynaptic serotonin (5‐hydroxytryptamine, 5‐HT) transporters (SERT) regulate 5‐HT signaling via antidepressant‐sensitive clearance of released neurotransmitter. Polymorphisms in the human SERT gene (SLC6A4) have been linked to risk for multiple neuropsychiatric disorders, including depression, obsessive‐compulsive disorder and autism. Using BXD recombinant inbred mice, a genetic reference population that can support the discovery of novel determinants of complex traits, merging collective trait assessments with bioinformatics approaches, we examine phenotypic and molecular networks associated with SERT gene and protein expression. Correlational analyses revealed a network of genes that significantly associated with SERT mRNA levels. We quantified SERT protein expression levels and identified region‐ and gender‐specific quantitative trait loci (QTLs), one of which associated with male midbrain SERT protein expression, centered on the protocadherin‐15 gene (Pcdh15), overlapped with a QTL for midbrain 5‐HT levels. Pcdh15 was also the only QTL‐associated gene whose midbrain mRNA expression significantly associated with both SERT protein and 5‐HT traits, suggesting an unrecognized role of the cell adhesion protein in the development or function of 5‐HT neurons. To test this hypothesis, we assessed SERT protein and 5‐HT traits in the Pcdh15 functional null line (Pcdh15^av‐3J), studies that revealed a strong, negative influence of Pcdh15 on these phenotypes. Together, our findings illustrate the power of multidimensional profiling of recombinant inbred lines in the analysis of molecular networks that support synaptic signaling, and that, as in the case of Pcdh15, can reveal novel relationships that may underlie risk for mental illness .     

Complex control of GABA(A) receptor subunit mRNA expression: variation, covariation, and genetic regulation

GABA type-A receptors are essential for fast inhibitory neurotransmission and are critical in brain function. Surprisingly, expression of receptor subunits is highly variable among individuals, but the cause and impact of this fluctuation remains unknown. We have studied sources of variation for all 19 receptor subunits using massive expression data sets collected across multiple brain regions and platforms in mice and humans. Expression of Gabra1, Gabra2, Gabrb2, Gabrb3, and Gabrg2 is highly variable and heritable among the large cohort of BXD strains derived from crosses of fully sequenced parents--C57BL/6J and DBA/2J. Genetic control of these subunits is complex and highly dependent on tissue and mRNA region. Remarkably, this high variation is generally not linked to phenotypic differences. The single exception is Gabrb3, a locus that is linked to anxiety. We identified upstream genetic loci that influence subunit expression, including three unlinked regions of chromosome 5 that modulate the expression of nine subunits in hippocampus, and that are also associated with multiple phenotypes. Candidate genes within these loci include, Naaa, Nos1, and Zkscan1. We confirmed a high level of coexpression for subunits comprising the major channel--Gabra1, Gabrb2, and Gabrg2--and identified conserved members of this expression network in mice and humans. Gucy1a3, Gucy1b3, and Lis1 are novel and conserved associates of multiple subunits that are involved in inhibitory signaling. Finally, proximal and distal regions of the 3' UTRs of single subunits have remarkably independent expression patterns in both species. However, corresponding regions of different subunits often show congruent genetic control and coexpression (proximal-to-proximal or distal-to-distal), even in the absence of sequence homology. Our findings identify novel sources of variation that modulate subunit expression and highlight the extraordinary capacity of biological networks to buffer 4-100 fold differences in mRNA levels.     

Identification of quantitative trait loci and candidate genes for an anxiolytic‐like response to ethanol in BXD recombinant inbred strains

Genetic differences in acute behavioral responses to ethanol contribute to the susceptibility to alcohol use disorder and the reduction of anxiety is a commonly reported motive underlying ethanol consumption among alcoholics. Therefore, we studied the genetic variance in anxiolytic‐like responses to ethanol across the BXD recombinant inbred (RI) mouse panel using the light–dark transition model of anxiety. Strain‐mean genetic mapping and a mixed‐model quantitative trait loci (QTL) analysis replicated several previously published QTL for locomotor activity and identified several novel anxiety‐related loci. Significant loci included a chromosome 11 saline anxiety‐like QTL (Salanq1) and a chromosome 12 locus (Etanq1) influencing the anxiolytic‐like response to ethanol. Etanq1 was successfully validated by studies with BXD advanced intercross strains and fine‐mapped to a region comprising less than 3.5 Mb. Through integration of genome‐wide mRNA expression profiles of the mesocorticolimbic reward circuit (prefrontal cortex, nucleus accumbens and ventral midbrain) across the BXD RI panel, we identified high priority candidate genes within Etanq1, the strongest of which was Ninein (Nin), a Gsk3β‐interacting protein that is highly expressed in the brain.     

Systems genetic analysis of hippocampal neuroanatomy and spatial learning in mice

Variation in hippocampal neuroanatomy correlates well with spatial learning ability in mice. Here, we have studied both hippocampal neuroanatomy and behavior in 53 isogenic BXD recombinant strains derived from C57BL/6J and DBA/2J parents. A combination of experimental, neuroinformatic and systems genetics methods was used to test the genetic bases of variation and covariation among traits. Data were collected on seven hippocampal subregions in CA3 and CA4 after testing spatial memory in an eight‐arm radial maze task. Quantitative trait loci were identified for hippocampal structure, including the areas of the intra‐ and infrapyramidal mossy fibers (IIPMFs), stratum radiatum and stratum pyramidale, and for a spatial learning parameter, error rate. We identified multiple loci and gene variants linked to either structural differences or behavior. Gpc4 and Tenm2 are strong candidate genes that may modulate IIPMF areas. Analysis of gene expression networks and trait correlations highlight several processes influencing morphometrical variation and spatial learning.     

Genetic loci for diet-induced atherosclerotic lesions and plasma lipids in mice

Genetic factors independent of those affecting plasma lipid levels are a major contributor to risk for atherosclerosis in humans, yet the basis for these is poorly understood. This study examined plasma lipids and diet-induced atherosclerosis in 16-month-old female mice of strains C56BL/6J and DBA/2J. Mice of the parental strains, from recombinant inbred strains derived from these (BXD RI), and F₂ progeny were fed an atherogenic diet for 16 weeks, beginning at 1 year of age. This induced atherosclerotic lesion formation in both parental strains, accompanied by increased plasma LDL levels. However, individual BXD RI strains and the BXD F₂ mice demonstrated a range of atherosclerotic lesion formation that was not or at best weakly correlated with plasma lipid levels. Quantitative trait locus (QTL) analysis of the BXD F₂ mice identified a locus with significant linkage (lod 4.5) for aortic lesion size on Chromosome (Chr) 10 that was independent of plasma lipids. Other loci with suggestive or significant linkage for various plasma lipid measures were identified on Chr 2, 3, 4, 5, 6, 7, 11, and 17. In this intercross, the genes primarily influencing atherosclerosis are distinct from those controlling plasma lipid levels.     

Systems Genetics of Metabolism: The Use of the BXD Murine Reference Panel for Multiscalar Integration of Traits

Metabolic homeostasis is achieved by complex molecular and cellular networks that differ significantly among individuals and are difficult to model with genetically engineered lines of mice optimized to study single gene function. Here, we systematically acquired metabolic phenotypes by using the EUMODIC EMPReSS protocols across a large panel of isogenic but diverse strains of mice (BXD type) to study the genetic control of metabolism. We generated and analyzed 140 classical phenotypes and deposited these in an open-access web service for systems genetics (www.genenetwork.org). Heritability, influence of sex, and genetic modifiers of traits were examined singly and jointly by using quantitative-trait locus (QTL) and expression QTL-mapping methods. Traits and networks were linked to loci encompassing both known variants and novel candidate genes, including alkaline phosphatase (ALPL), here linked to hypophosphatasia. The assembled and curated phenotypes provide key resources and exemplars that can be used to dissect complex metabolic traits and disorders.     

Systems genetics of sensation seeking

Sensation seeking is a multifaceted, heritable trait which predicts the development of substance use and abuse in humans; similar phenomena have been observed in rodents. Genetic correlations among sensation seeking and substance use indicate shared biological mechanisms, but the genes and networks underlying these relationships remain elusive. Here, we used a systems genetics approach in the BXD recombinant inbred mouse panel to identify shared genetic mechanisms underlying substance use and preference for sensory stimuli, an intermediate phenotype of sensation seeking. Using the operant sensation seeking (OSS) paradigm, we quantified preference for sensory stimuli in 120 male and 127 female mice from 62 BXD strains and the C57BL/6J and DBA/2J founder strains. We used relative preference for the active and inactive levers to dissociate preference for sensory stimuli from locomotion and exploration phenotypes. We identified genomic regions on chromosome 4 (155.236‐155.742 Mb) and chromosome 13 (72.969‐89.423 Mb) associated with distinct behavioral components of OSS. Using publicly available behavioral data and mRNA expression data from brain regions involved in reward processing, we identified (a) genes within these behavioral QTL exhibiting genome‐wide significant cis‐eQTL and (b) genetic correlations among OSS phenotypes, ethanol phenotypes and mRNA expression. From these analyses, we nominated positional candidates for behavioral QTL associated with distinct OSS phenotypes including Gnb1 and Mef2c. Genetic covariation of Gnb1 expression, preference for sensory stimuli and multiple ethanol phenotypes suggest that heritable variation in Gnb1 expression in reward circuitry partially underlies the widely reported relationship between sensation seeking and substance use.     

Genotyping new BXD recombinant inbred mouse strains and comparison of BXD and consensus maps

Nine additional BXD recombinant inbred (RI) strains have been developed from the F₂ cross of C57BL/6J and DBA/2J mouse strains. A tenth line stopped breeding in the F₁₂ generation. F₂₀ generation breeding pairs from the nine surviving strains and an F₁₂ pair from the extinct line were genotyped at 319 genetic markers (primarily microsatellites) spanning most of the genome. Where typing data were lacking, the established set of 26 BXD strains also were genotyped at these same loci. The availability of these additional nine strains enhances the value of the BXD RI set for analysis of complex phenotypic traits. The proportion of loci still segregating at the F₂₀ generation was found to closely approximate expectation, suggesting that selection favoring the retention of heterozygosity is not a strong factor. However, the number of crossovers between adjacent markers was frequently less than predicted from consensus map distances. A significant deficiency of recombinants was observed on Chrs 3, 4, 14, and X. On Chr 14, the estimated cumulative BXD map distance between the most proximal and distal markers was only 30.2 cM, compared with a distance of 60.0 cM in the consensus map. On the X Chr, the estimated and predicted cumulative distances were 38.8 and 69.5 cM, respectively. Over all chromosomes, the BXD RI map is 14.5% shorter than predicted from the consensus map. It is suggested that distances in some of the consensus maps are inflated. Alternatively, recombinant genotypes could be selected against during inbreeding owing to allelic interactions affecting fitness. The latter interpretation implies that relatively strong intrachromosomal epistasis is common. Received: 2 October 1998 / Accepted: 15 December 1998     

Genetic analysis of the neurosteroid deoxycorticosterone and its relation to alcohol phenotypes: identification of QTLs and downstream gene regulation

Background

Deoxycorticosterone (DOC) is an endogenous neurosteroid found in brain and serum, precursor of the GABAergic neuroactive steroid (3α,5α)-3,21-dihydroxypregnan-20-one (tetrahydrodeoxycorticosterone, THDOC) and the glucocorticoid corticosterone. These steroids are elevated following stress or ethanol administration, contribute to ethanol sensitivity, and their elevation is blunted in ethanol dependence.

Methodology/Principal Findings

To systematically define the genetic basis, regulation, and behavioral significance of DOC levels in plasma and cerebral cortex we examined such levels across 47 young adult males from C57BL/6J (B6)×DBA/2J (D2) (BXD) mouse strains for quantitative trait loci (QTL) and bioinformatics analyses of behavior and gene regulation. Mice were injected with saline or 0.075 mg/kg dexamethasone sodium salt at 8:00 am and were sacrificed 6 hours later. DOC levels were measured by radioimmunoassay. Basal cerebral cortical DOC levels ranged between 1.4 and 12.2 ng/g (8.7-fold variation, p<0.0001) with a heritability of ∼0.37. Basal plasma DOC levels ranged between 2.8 and 12.1 ng/ml (4.3-fold variation, p<0.0001) with heritability of ∼0.32. QTLs for basal DOC levels were identified on chromosomes 4 (cerebral cortex) and 14 (plasma). Dexamethasone-induced changes in DOC levels showed a 4.4-fold variation in cerebral cortex and a 4.1-fold variation in plasma, but no QTLs were identified. DOC levels across BXD strains were further shown to be co-regulated with networks of genes linked to neuronal, immune, and endocrine function. DOC levels and its responses to dexamethasone were associated with several behavioral measures of ethanol sensitivity previously determined across the BXD strains by multiple laboratories.

Conclusions/Significance

Both basal and dexamethasone-suppressed DOC levels are positively correlated with ethanol sensitivity suggesting that the neurosteroid DOC may be a putative biomarker of alcohol phenotypes. DOC levels were also strongly correlated with networks of genes associated with neuronal function, innate immune pathways, and steroid metabolism, likely linked to behavioral phenotypes.     

Melatonin synthesis in the optic lobes and midbrain of the grasshopper Oedipoda caerulescens

The pathways of insect melatonin (MEL) biosynthesis apparently follow the same routes as those identified in vertebrates but information on MEL synthesis variations related with serotonin (5‐HT), 5‐hydroxy‐indole acetic acid (5HIAA), and N‐acetylserotonin (NAS) levels, as well as 5‐HT N‐acetyltransferase (NAT) activity throughout the day, is very limited in the insect nervous system. In the present study, the levels of MEL, metabolites (5‐HT, NAS, and 5‐HIAA) and enzyme NAT were determined in the optic lobes and the midbrain of the grasshopper Oedipoda caerulescens, in conditions of light and darkness. In both tissues, a different pattern of MEL synthesis was observed over the light/dark cycle. Variations in the levels of 5‐HT, NAS and NAT activity related to the synthesis of cerebral MEL follow a pattern very similar to that observed in the pineal of mammals, with a peak of synthesis in the first half of the scotophase. Also, we observed differences in the metabolism of 5‐HT between the optic lobes and the midbrain light/dark‐dependent.     

Persistent sexual behavior in castrated, recombinant inbred mice.

In contrast to other male mice, many hybrid B6D2F1 males continue to mate after being castrated in adulthood. This persistent sexual behavior in the castrated hybrids presumably is based upon heterozygosity that might occur within individual loci or between different loci. Those possibilities were discriminated by studying BXD recombinant inbred mice. In BXD strains, individual loci are homozygous, but heterozygosity has arisen between separate loci. Males in two of six BXD strains persisted in copulating for 5 mo after being castrated. Thus, it is concluded that retention of masculine sexual behavior in castrated mice is a function of heterozygosity between loci.     

Identifying human disease genes through cross-species gene mapping of evolutionary conserved processes

Background

Understanding complex networks that modulate development in humans is hampered by genetic and phenotypic heterogeneity within and between populations. Here we present a method that exploits natural variation in highly diverse mouse genetic reference panels in which genetic and environmental factors can be tightly controlled. The aim of our study is to test a cross-species genetic mapping strategy, which compares data of gene mapping in human patients with functional data obtained by QTL mapping in recombinant inbred mouse strains in order to prioritize human disease candidate genes.

Methodology

We exploit evolutionary conservation of developmental phenotypes to discover gene variants that influence brain development in humans. We studied corpus callosum volume in a recombinant inbred mouse panel (C57BL/6J×DBA/2J, BXD strains) using high-field strength MRI technology. We aligned mouse mapping results for this neuro-anatomical phenotype with genetic data from patients with abnormal corpus callosum (ACC) development.

Principal Findings

From the 61 syndromes which involve an ACC, 51 human candidate genes have been identified. Through interval mapping, we identified a single significant QTL on mouse chromosome 7 for corpus callosum volume with a QTL peak located between 25.5 and 26.7 Mb. Comparing the genes in this mouse QTL region with those associated with human syndromes (involving ACC) and those covered by copy number variations (CNV) yielded a single overlap, namely HNRPU in humans and Hnrpul1 in mice. Further analysis of corpus callosum volume in BXD strains revealed that the corpus callosum was significantly larger in BXD mice with a B genotype at the Hnrpul1 locus than in BXD mice with a D genotype at Hnrpul1 (F = 22.48, p<9.87*10⁻⁵).

Conclusion

This approach that exploits highly diverse mouse strains provides an efficient and effective translational bridge to study the etiology of human developmental disorders, such as autism and schizophrenia.     

New murine polymorphisms detected by random amplified polymorphic DNA (RAPD) PCR and mapped by use of recombinant inbred strains

Oligonucleotide primers of random sequence that were 12 bases in length, 58% in GC content, and lacking internal palindromes were designed. By random amplified polymorphic DNA (RAPD) PCR, these primers were used to survey for DNA variations between the progenitors of the mouse AXB and BXA recombinant inbred sets (A/J and C57BL/6J). We identified 17 DNA variants detected by 10 primers. Map positions for these variants were determined by comparing their strain distribution patterns in the AXB, BXA recombinant inbred sets with strain distribution patterns of previously published loci. When necessary, BXD and NXSM recombinant inbred sets were also used. These 17 new loci mapped to 12 chromosomes. The 10 primers were also used to survey 20 inbred mouse strains including the progenitors of other recombinant inbred sets and four mouse strains recently inbred from the wild (CAST/Ei, MOLF/Ei, PERA/Ei, and SPRET/Ei).     

Segregation patterns of endogenous mouse mammary tumor viruses in five recombinant inbred strain sets.

We identified mouse mammary tumor proviral loci in the AKR/J, C3H/HeJ, C57BL/6J, C57L/J, DBA/2J, and SWR/J inbred mouse strains and determined their segregation patterns in the AKXD, AKXL, BXD, BXH, and SWXL recombinant inbred strain sets. Two new Mtv loci, Mtv-29 and Mtv-30, were identified. Mtv-30 was genetically mapped to chromosome 12. Additionally, two previously identified Mtv loci, Mtv-14 and Mtv-23, were genetically mapped to chromosome 4 and chromosome 6, respectively.     

Genetic regulation of endotoxin-induced airway disease

To identify novel genes regulating the biologic response to lipopolysaccharide (LPS), we used a combination of quantitative trait locus (QTL) analysis and microarray-based gene expression studies of C57BL/6J x DBA/2J(BXD) F2 and recombinant inbred (RI) mice. A QTL affecting pulmonary TNF-alpha production was identified on chromosome 2, and a region affecting both polymorphonuclear leukocyte recruitment and TNF-alpha levels was identified on chromosome 11. Microarray analyses of unchallenged and LPS-challenged BXD RI strains identified approximately 500 genes whose expression was significantly changed by inhalation of LPS. Of these genes, 28 reside within the chromosomal regions identified by the QTL analyses, implicating these genes as high priority candidates for functional studies. Additional high priority candidate genes were identified based on their differential expression in mice having high and low responses to LPS. Functional studies of these genes are expected to reveal important molecular mechanisms regulating the magnitude of biologic responses to LPS.     

The physical separation of Lps and Ifa loci in BXH recombinant inbred mice

Several reports in the literature suggest that many of the phenotypic defects of LPS-hyporesponsive C3H/HeJ mice may be attributed to decreased IFN production by their macrophages. The physical proximity on chromosome 4 of the gene which encodes sensitivity to LPS (Lps) and the structural gene cluster which encodes IFN-alpha (Ifa), suggests the possibility that the Lps gene, whose product is unknown, may actually be a part of the Ifa gene cluster. The C57BL/6J and C3H/HeJ mouse strains carry distinct alleles at both the LPs and the Ifa loci. In this study, we have analyzed these parental strains, as well as 12 recombinant inbred strains derived from these parental strains (e.g., BXH strains), for inheritance of these distinct alleles. The results show the segregation of these two loci: in 5 of 12 BXH RI strains, the IFN-alpha restriction fragment length polymorphism characteristic of one parental strain was discordant with the predicted LPS response phenotype. Therefore, we conclude that the Lps and the Ifa genes are physically distinct despite the apparent cause and effect relationship which is observed phenotypically.     

The effect of the Cmv-1 resistance gene, which is linked to the natural killer cell gene complex, is mediated by natural killer cells.

The resistance of mice to lethal infection by murine CMV (MCMV) is under complex host genetic control with contributions from both H-2 and non-H-2 genes. We have previously shown that an autosomal, non-MHC encoded gene, Cmv-1, controls MCMV replication in the spleen. We have investigated the mechanism by which the Cmv-1 resistance gene confers protection against MCMV infection. Using H-2 compatible irradiation bone marrow chimeras, the enhanced resistance to MCMV infection that is associated with the Cmv-1l allele in the C57BL background was shown to be mediated by an irradiation-sensitive bone marrow-derived cell population, or a factor produced by these cells. The lack of correlation between serum IFN titers and the strain distribution pattern of Cmv-1 in CXB recombinant inbred mouse strains suggests that IFN does not mediate resistance conferred by this gene. Similarly, the lack of effect of in vivo depletion of mature CD4+ and CD8+ T cells on virus replication in C57BL/6J mice indicates that T cells are unlikely to be involved. In contrast, in vivo depletion of NK cells by injection of the anti-NK1.1 mAb PK136 abrogated restricted splenic virus replication in C57BL/6J----BALB.B chimeric mice and in the Cmv-1l CXB strains. These data indicate that the effect of the Cmv-1 gene is mediated by NK cells. The significant augmentation in NK cell activity after MCMV infection of the susceptible Cmv-1h strains (BALB/cBy), CXBG/By, CXBH/By, CXBI/By, and CXBK/By) indicates the existence in these mice of NK cells that are functionally and phenotypically distinct from those in Cmv-1l strains. NK cells present in the Cmv-1h strains are unable to restrict efficiently splenic MCMV replication in vivo, possibly due to a lack of specificity for virus-infected target cells. Finally, flow cytometric analysis of NK1-1 expression in CXB and BXD RI mice together with MCMV replication studies in the BXD RI strains indicate that Cmv-1 is closely linked to NK1.1 and other loci that reside on a distal segment of murine chromosome 6 in a region that has recently been defined as the natural killer complex.     

Spontaneous mutations in recombinant inbred mice: mutant toll-like receptor 4 (Tlr4) in BXD29 mice

Recombinant inbred (RI) mice are frequently used to identify QTL that underlie differences in measurable phenotypes between two inbred strains of mice. Here we show that one RI strain, C57BL/6J x DBA/2J (BXD29), does not develop an inflammatory response following inhalation of LPS. Approximately 25% of F2 mice [F1(BXD29 x DBA/2J) x F1] are also unresponsive to inhaled LPS, suggesting the presence of a recessive mutation in the BXD29 strain. A genomic scan of these F2 mice revealed that unresponsive animals, but not responsive animals, are homozygous for C57BL/6J DNA at a single locus on chromosome 4 close to the genomic location of Tlr4. All progeny between BXD29 and gene-targeted Tlr4-deficient mice are unresponsive to inhaled LPS, suggesting that the mutation in the BXD29 strain is allelic with Tlr4. Moreover, the intact Tlr4 receptor is not displayed on the cell surface of BXD29 macrophages. Finally, a molecular analysis of the Tlr4 gene in BXD29 mice revealed that it is interrupted by a large insertion of repetitive DNA. These findings explain the unresponsiveness of BXD29 mice to LPS and suggest that data from BXD29 mice should not be included when using BXD mice to study phenotypes affected by Tlr4 function. Our results also suggest that the frequency of such unidentified, spontaneously occurring mutations is an issue that should be considered when RI strains are used to identify QTL.     

Genetic maps of mouse chromosome 17 including 12 new anonymous DNA loci and 25 anchor loci.

An interspecific backcross between lab mice and Mus spretus was used to construct a multilocus map of Chromosome 17 consisting of 12 new anonymous loci and 9 anchor loci. In addition, 7 anonymous DNA loci were added to the Chr 17 map for the BXD strains. Although we were able to identify readily the most likely gene order in the interspecific backcross, we found no evidence for an unambiguous gene order using the BXD recombinant inbred strains. Comparison of the interspecific backcross map and the BXD RI strain map revealed evidence in the interspecific backcross for a longer total genetic length, enhanced recombination distal to H-2, a segment showing suppressed recombination, and strong interference.     

Systems genetics analysis of iron regulation in the brain

Iron imbalances in the brain, including excess accumulation and deficiency, are associated with neurological disease and dysfunction; yet, their origins are poorly understood. Using systems genetics analysis, we have learned that large individual differences exist in brain iron concentrations, even in the absence of neurological disease. Much of the individual differences can be tied to the genetic makeup of the individual. This genetic-based differential regulation can be modeled in genetic reference populations of rodents. The work in our laboratory centers on iron regulation in the brain and our animal model consists of 25 BXD/Ty recombinant inbred mouse strains. By studying naturally occurring variation in iron phenotypes, such as tissue iron concentration, we can tie that variability to one or more genes by way of quantitative trait loci (QTL) analysis. Moreover, we can conduct genetic correlation analyses between our phenotypes and others previously measured in the BXD/Ty strains. We have observed several suggestive QTL related to ventral midbrain iron content, including one on chromosome 17 that contains btbd9, a gene that in humans has been associated with restless legs syndrome and serum ferritin. We have also observed gene expression correlations with ventral midbrain iron, including btbd9 expression and dopamine receptor expression. In addition, we have observed significant correlations between ventral midbrain iron content and dopamine-related phenotypes. The following is a discussion of iron regulation in the brain and the contributions a systems genetics approach can make toward understanding the genetic underpinnings and relation to neurological disease.