Unexpected complexity in the haplotypes of commonly used inbred strains of laboratory mice

Large differences in collateral blood vessel abundance among individuals arise from multiple genetic variants

Collateral blood flow varies greatly among humans for reasons that remain unclear, resulting in significant differences in ischemic tissue damage. A similarly large variation has also been found in mice that is caused by genetic background-dependent differences in the extent of collateral formation, termed collaterogenesis-a unique angiogenic process that occurs during development and determines collateral number and diameter in the adult. Previous studies have identified several quantitative trait loci (QTL) linked to this variation. However, understanding has been hampered by the use of closely related inbred strains that do not model the wide genetic variation present in the "outbred" human population. The Collaborative Cross (CC) multiparent mouse genetic reference panel was developed to address this limitation. Herein we measured the number and average diameter of cerebral collaterals in 60 CC strains, their 8 founder strains, 8 F1 crosses of CC strains selected for abundant versus sparse collaterals, and 2 intercross populations created from the latter. Collateral number evidenced 47-fold variation among the 60 CC strains, with 14% having poor, 25% poor-to-intermediate, 47% intermediate-to-good, and 13% good collateral abundance, that was associated with large differences in post-stroke infarct volume. Collateral number in skeletal muscle and intestine of selected high- and low-collateral strains evidenced the same relative abundance as in brain. Genome-wide mapping demonstrated that collateral abundance is a highly polymorphic trait. Subsequent analysis identified: 6 novel QTL circumscribing 28 high-priority candidate genes harboring putative loss-of-function polymorphisms (SNPs) associated with low collateral number; 335 predicted-deleterious SNPs present in their human orthologs; and 32 genes associated with vascular development but lacking protein coding variants. Six additional suggestive QTL (LOD > 4.5) were also identified in CC-wide QTL mapping. This study provides a comprehensive set of candidate genes for future investigations aimed at identifying signaling proteins within the collaterogenesis pathway whose variants potentially underlie genetic-dependent collateral insufficiency in brain and other tissues.

Large differences in collateral blood vessel abundance among individuals arise from multiple genetic variants

Collateral blood flow varies greatly among humans for reasons that remain unclear, resulting in significant differences in ischemic tissue damage. A similarly large variation has also been found in mice that is caused by genetic background-dependent differences in the extent of collateral formation, termed collaterogenesis-a unique angiogenic process that occurs during development and determines collateral number and diameter in the adult. Previous studies have identified several quantitative trait loci (QTL) linked to this variation. However, understanding has been hampered by the use of closely related inbred strains that do not model the wide genetic variation present in the "outbred" human population. The Collaborative Cross (CC) multiparent mouse genetic reference panel was developed to address this limitation. Herein we measured the number and average diameter of cerebral collaterals in 60 CC strains, their 8 founder strains, 8 F1 crosses of CC strains selected for abundant versus sparse collaterals, and 2 intercross populations created from the latter. Collateral number evidenced 47-fold variation among the 60 CC strains, with 14% having poor, 25% poor-to-intermediate, 47% intermediate-to-good, and 13% good collateral abundance, that was associated with large differences in post-stroke infarct volume. Genome-wide mapping demonstrated that collateral abundance is a highly polymorphic trait. Subsequent analysis identified: 6 novel QTL circumscribing 28 high-priority candidate genes harboring putative loss-of-function polymorphisms (SNPs) associated with low collateral number; 335 predicted-deleterious SNPs present in their human orthologs; and 32 genes associated with vascular development but lacking protein coding variants. This study provides a comprehensive set of candidate genes for future investigations aimed at identifying signaling proteins within the collaterogenesis pathway whose variants potentially underlie genetic-dependent collateral insufficiency in brain and other tissues.

Hearing loss and tinnitus: association studies for complex-hearing disorders in mouse and man

Genome-wide association studies (GWAS) provide an unbiased first look at genetic loci involved in aging and noise-induced sensorineural hearing loss and tinnitus. The hearing phenotype, whether audiogram-based or self-report, is regressed against genotyped information at representative single nucleotide polymorphisms (SNPs) across the genome. Findings include the fact that both hearing loss and tinnitus are polygenic disorders, with up to thousands of genes, each of effect size of < 0.02. Smaller human GWAS' were able to use objective measures and identified a few loci; however, hundreds of thousands of participants have been required for the statistical power to identify significant variants, and GWAS is unable to assess rare variants with mean allele frequency < 1%. Animal studies are required as well because of inability to access the human cochlea. Mouse GWAS builds on linkage techniques and the known phenotypic differences in auditory function between inbred strains. With the advantage that the laboratory environment can be controlled for noise and aging, the Hybrid Mouse Diversity Panel (HDMP) combines 100 strains sequenced at high resolution. Lift-over regions between mice and humans have identified over 17,000 homologous genes. Since most significant SNPs are either intergenic or in introns, and binding sites between species are poorly preserved between species, expression quantitative trait locus information is required to bring humans and mice into agreement. Transcriptome-wide analysis studies (TWAS) can prioritize putative causal genes and tissues. Diverse species, each making a distinct contribution, carry a synergistic advantage in the quest for treatment and ultimate cure of sensorineural hearing difficulties.

What Has a Century of Quantitative Genetics Taught Us About Nature's Genetic Tool Kit?

The complexity of heredity has been appreciated for decades: Many traits are controlled not by a single genetic locus but instead by polymorphisms throughout the genome. The importance of complex traits in biology and medicine has motivated diverse approaches to understanding their detailed genetic bases. Here, we focus on recent systematic studies, many in budding yeast, which have revealed that large numbers of all kinds of molecular variation, from noncoding to synonymous variants, can make significant contributions to phenotype. Variants can affect different traits in opposing directions, and their contributions can be modified by both the environment and the epigenetic state of the cell. The integration of prospective (synthesizing and analyzing variants) and retrospective (examining standing variation) approaches promises to reveal how natural selection shapes quantitative traits. Only by comprehensively understanding nature's genetic tool kit can we predict how phenotypes arise from the complex ensembles of genetic variants in living organisms.

Collaborative Cross Mice Yield Genetic Modifiers for Pseudomonas aeruginosa Infection in Human Lung Disease

Human genetics influence a range of pathological and clinical phenotypes in respiratory infections; however, the contributions of disease modifiers remain underappreciated. We exploited the Collaborative Cross (CC) mouse genetic-reference population to map genetic modifiers that affect the severity of Pseudomonas aeruginosa lung infection. Screening for P. aeruginosa respiratory infection in a cohort of 39 CC lines exhibits distinct disease phenotypes ranging from complete resistance to lethal disease. Based on major changes in the survival times, a quantitative-trait locus (QTL) was mapped on murine chromosome 3 to the genomic interval of Mb 110.4 to 120.5. Within this locus, composed of 31 protein-coding genes, two candidate genes, namely, dihydropyrimidine dehydrogenase (Dpyd) and sphingosine-1-phosphate receptor 1 (S1pr1), were identified according to the level of genome-wide significance and disease gene prioritization. Functional validation of the S1pr1 gene by pharmacological targeting in C57BL/6NCrl mice confirmed its relevance in P. aeruginosa pathophysiology. However, in a cohort of Canadian patients with cystic fibrosis (CF) disease, regional genetic-association analysis of the syntenic human locus on chromosome 1 (Mb 97.0 to 105.0) identified two single-nucleotide polymorphisms (rs10875080 and rs11582736) annotated to the Dpyd gene that were significantly associated with age at first P. aeruginosa infection. Thus, there is evidence that both genes might be implicated in this disease. Our results demonstrate that the discovery of murine modifier loci may generate information that is relevant to human disease progression.IMPORTANCE Respiratory infection caused by P. aeruginosa is one of the most critical health burdens worldwide. People affected by P. aeruginosa infection include patients with a weakened immune system, such as those with cystic fibrosis (CF) genetic disease or non-CF bronchiectasis. Disease outcomes range from fatal pneumonia to chronic life-threatening infection and inflammation leading to the progressive deterioration of pulmonary function. The development of these respiratory infections is mediated by multiple causes. However, the genetic factors underlying infection susceptibility are poorly known and difficult to predict. Our study employed novel approaches and improved mouse disease models to identify genetic modifiers that affect the severity of P. aeruginosa lung infection. We identified candidate genes to enhance our understanding of P. aeruginosa infection in humans and provide a proof of concept that could be exploited for other human pathologies mediated by bacterial infection.

Strong conservation of inbred mouse strain microRNA loci but broad variation in brain microRNAs due to RNA editing and isomiR expression

Diversity in the structure and expression of microRNAs, important regulators of gene expression, arises from SNPs, duplications followed by divergence, production of isomiRs, and RNA editing. Inbred mouse strains and crosses using them are important reference populations for genetic mapping, and as models of human disease. We determined the nature and extent of interstrain miRNA variation by (i) identifying miRNA SNPs in whole-genome sequence data from 36 strains, and (ii) examining miRNA editing and expression in hippocampus (Hpc) and frontal cortex (FCx) of six strains, to facilitate the study of miRNAs in neurobehavioral phenotypes. miRNA loci were strongly conserved among the 36 strains, but even the highly conserved seed region contained 16 SNPs. In contrast, we identified RNA editing in 58.9% of miRNAs, including 11 consistent editing events in the seed region. We confirmed the functional significance of three conserved edits in the miR-379/410 cluster, demonstrating that edited miRNAs gained novel target mRNAs not recognized by the unedited miRNAs. We found significant interstrain differences in miRNA and isomiR expression: Of 779 miRNAs expressed in Hpc and 719 in FCx, 262 were differentially expressed (190 in Hpc, 126 in FCx, 54 in both). We also identified 32 novel miRNA candidates using miRNA prediction tools. Our studies provide the first comprehensive analysis of SNP, isomiR, and RNA editing variation in miRNA loci across inbred mouse strains, and a detailed catalog of expressed miRNAs in Hpc and FCx in six commonly used strains. These findings will facilitate the molecular analysis of neurological and behavioral phenotypes in this model organism.

Mapping Causal Variants with Single-Nucleotide Resolution Reveals Biochemical Drivers of Phenotypic Change

Understanding the sequence determinants that give rise to diversity among individuals and species is the central challenge of genetics. However, despite ever greater numbers of sequenced genomes, most genome-wide association studies cannot distinguish causal variants from linked passenger mutations spanning many genes. We report that this inherent challenge can be overcome in model organisms. By pushing the advantages of inbred crossing to its practical limit in Saccharomyces cerevisiae, we improved the statistical resolution of linkage analysis to single nucleotides. This "super-resolution" approach allowed us to map 370 causal variants across 26 quantitative traits. Missense, synonymous, and cis-regulatory mutations collectively gave rise to phenotypic diversity, providing mechanistic insight into the basis of evolutionary divergence. Our data also systematically unmasked complex genetic architectures, revealing that multiple closely linked driver mutations frequently act on the same quantitative trait. Single-nucleotide mapping thus complements traditional deletion and overexpression screening paradigms and opens new frontiers in quantitative genetics.

Genetics of metabolic syndrome: potential clues from wild-derived inbred mouse strains

The metabolic syndrome (MetS) is a complex constellation of metabolic abnormalities including obesity, abnormal glucose metabolism, dyslipidemia, and elevated blood pressure that together substantially increase risk for cardiovascular disease and Type 2 diabetes. Both genetic and environmental factors contribute to the development of MetS, but this process is still far from understood. Human studies have revealed only part of the underlying basis. Studies in mice offer many strengths that can complement human studies to help elucidate the etiology and pathophysiology of MetS. Here we review the ways mice can contribute to MetS research. In particular, we focus on the information that can be obtained from studies of the inbred strains, with specific focus on the phenotypes of the wild-derived inbred strains. These are newly derived inbred strains that were created from wild-caught mice. They contain substantial genetic variation that is not present in the classical inbred strains, have phenotypes of relevance for MetS, and various mouse strain resources have been created to facilitate the mining of this new genetic variation. Thus studies using wild-derived inbred strains hold great promise for increasing our understanding of MetS.

Whole exome sequencing of wild-derived inbred strains of mice improves power to link phenotype and genotype

The house mouse is a powerful model to dissect the genetic basis of phenotypic variation, and serves as a model to study human diseases. Despite a wealth of discoveries, most classical laboratory strains have captured only a small fraction of genetic variation known to segregate in their wild progenitors, and existing strains are often related to each other in complex ways. Inbred strains of mice independently derived from natural populations have the potential to increase power in genetic studies with the addition of novel genetic variation. Here, we perform exome-enrichment and high-throughput sequencing (~8× coverage) of 26 wild-derived strains known in the mouse research community as the "Montpellier strains." We identified 1.46 million SNPs in our dataset, approximately 19% of which have not been detected from other inbred strains. This novel genetic variation is expected to contribute to phenotypic variation, as they include 18,496 nonsynonymous variants and 262 early stop codons. Simulations demonstrate that the higher density of genetic variation in the Montpellier strains provides increased power for quantitative genetic studies. Inasmuch as the power to connect genotype to phenotype depends on genetic variation, it is important to incorporate these additional genetic strains into future research programs.

The Genetic Architecture of Noise-Induced Hearing Loss: Evidence for a Gene-by-Environment Interaction

The discovery of environmentally specific genetic effects is crucial to the understanding of complex traits, such as susceptibility to noise-induced hearing loss (NIHL). We describe the first genome-wide association study (GWAS) for NIHL in a large and well-characterized population of inbred mouse strains, known as the Hybrid Mouse Diversity Panel (HMDP). We recorded auditory brainstem response (ABR) thresholds both pre and post 2-hr exposure to 10-kHz octave band noise at 108 dB sound pressure level in 5–6-wk-old female mice from the HMDP (4–5 mice/strain). From the observation that NIHL susceptibility varied among the strains, we performed a GWAS with correction for population structure and mapped a locus on chromosome 6 that was statistically significantly associated with two adjacent frequencies. We then used a “genetical genomics” approach that included the analysis of cochlear eQTLs to identify candidate genes within the GWAS QTL. In order to validate the gene-by-environment interaction, we compared the effects of the postnoise exposure locus with that from the same unexposed strains. The most significant SNP at chromosome 6 (rs37517079) was associated with noise susceptibility, but was not significant at the same frequencies in our unexposed study. These findings demonstrate that the genetic architecture of NIHL is distinct from that of unexposed hearing levels and provide strong evidence for gene-by-environment interactions in NIHL.

Nutrition and the science of disease prevention: a systems approach to support metabolic health

Progress in nutritional science, genetics, computer science, and behavioral economics can be leveraged to address the challenge of noncommunicable disease. This report highlights the connection between nutrition and the complex science of preventing disease and discusses the promotion of optimal metabolic health, building on input from several complementary disciplines. The discussion focuses on (1) the basic science of optimal metabolic health, including data from gene-diet interactions, microbiome, and epidemiological research in nutrition, with the goal of defining better targets and interventions, and (2) how nutrition, from pharma to lifestyle, can build on systems science to address complex issues.

Collaborative cross mice in a genetic association study reveal new candidate genes for bone microarchitecture

Background

The microstructure of trabecular bone is a composite trait governed by a complex interaction of multiple genetic determinants. Identifying these genetic factors should significantly improve our ability to predict of osteoporosis and its associated risks. Genetic mapping using collaborative cross mice (CC), a genetically diverse recombinant inbred mouse reference panel, offers a powerful tool to identify causal loci at a resolution under one mega base-pairs, with a relatively small cohort size.

Here, we utilized 31 CC lines (160 mice of both sexes in total) to perform genome-wide haplotype mapping across 77,808 single-nucleotide polymorphisms (SNPs). Haplotype scans were refined by imputation with the catalogue of sequence variation segregating in the CC to suggest potential candidate genes. Trabecular traits were obtained following microtomographic analysis, performed on 10-μm resolution scans of the femoral distal metaphysis. We measured the trabecular bone volume fraction (BV/TV), number (Tb.N), thickness (Tb.Th), and connectivity density (Conn.D).

Results

Heritability of these traits ranged from 0.6 to 0.7. In addition there was a significant (P < 0.01) sex effect in all traits except Tb.Th. Our haplotype scans yielded six quantitative trait loci (QTL) at 1 % false discovery rate; BV/TV and Tb.Th produced two proximal loci each, on chromosome 2 and 7, respectively, and Tb.N and Conn.D yielded one locus on chromosomes 8 and 14, respectively. We identified candidate genes with previously-reported functions in bone biology, and implicated unexpected genes whose function in bone biology has yet to be assigned. Based on the literature, among the genes that ranked particularly high in our analyses (P < 10^-6) and which have a validated causal role in skeletal biology, are Avp, Oxt, B2m (associated with BV/TV), Cnot7 (with Tb.N), Pcsk6, Rgma (with Tb.Th), Rb1, and Cpb2 (with Conn.D). Other candidate genes strongly suggested by our analyses are Sgcz, Fgf20 (associated with Tb.N), and Chd2 (with Tb.Th).

Conclusion

We have demonstrated for the first time genome-wide significant association between several genetic loci and trabecular microstructural parameters for genes with previously reported experimental observations, as well as proposing a role for new candidate genes with no previously characterized skeletal function.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2213-x) contains supplementary material, which is available to authorized users.

The Genetic Architecture of Hearing Impairment in Mice: Evidence for Frequency-Specific Genetic Determinants

Genome-wide association studies (GWAS) have been successfully applied in humans for the study of many complex phenotypes. However, identification of the genetic determinants of hearing in adults has been hampered, in part, by the relative inability to control for environmental factors that might affect hearing throughout the lifetime, as well as a large degree of phenotypic heterogeneity. These and other factors have limited the number of large-scale studies performed in humans that have identified candidate genes that contribute to the etiology of this complex trait. To address these limitations, we performed a GWAS analysis using a set of inbred mouse strains from the Hybrid Mouse Diversity Panel. Among 99 strains characterized, we observed approximately two-fold to five-fold variation in hearing at six different frequencies, which are differentiated biologically from each other by the location in the cochlea where each frequency is registered. Among all frequencies tested, we identified a total of nine significant loci, several of which contained promising candidate genes for follow-up study. Taken together, our results indicate the existence of both genes that affect global cochlear function, as well as anatomical- and frequency-specific genes, and further demonstrate the complex nature of mammalian hearing variation.

Identification of a novel gene for diabetic traits in rats, mice, and humans

The genetic basis of type 2 diabetes remains incompletely defined despite the use of multiple genetic strategies. Multiparental populations such as heterogeneous stocks (HS) facilitate gene discovery by allowing fine mapping to only a few megabases, significantly decreasing the number of potential candidate genes compared to traditional mapping strategies. In the present work, we employed expression and sequence analysis in HS rats (Rattus norvegicus) to identify Tpcn2 as a likely causal gene underlying a 3.1-Mb locus for glucose and insulin levels. Global gene expression analysis on liver identified Tpcn2 as the only gene in the region that is differentially expressed between HS rats with glucose intolerance and those with normal glucose regulation. Tpcn2 also maps as a cis-regulating expression QTL and is negatively correlated with fasting glucose levels. We used founder sequence to identify variants within this region and assessed association between 18 variants and diabetic traits by conducting a mixed-model analysis, accounting for the complex family structure of the HS. We found that two variants were significantly associated with fasting glucose levels, including a nonsynonymous coding variant within Tpcn2. Studies in Tpcn2 knockout mice demonstrated a significant decrease in fasting glucose levels and insulin response to a glucose challenge relative to those in wild-type mice. Finally, we identified variants within Tpcn2 that are associated with fasting insulin in humans. These studies indicate that Tpcn2 is a likely causal gene that may play a role in human diabetes and demonstrate the utility of multiparental populations for positionally cloning genes within complex loci.

Susceptibility to Klebsiella pneumonaie infection in collaborative cross mice is a complex trait controlled by at least three loci acting at different time points

BACKGROUND

Klebsiella pneumoniae (Kp) is a bacterium causing severe pneumonia in immunocompromised hosts and is often associated with sepsis. With the rise of antibiotic resistant bacteria, there is a need for new effective and affordable control methods; understanding the genetic architecture of susceptibility to Kp will help in their development. We performed the first quantitative trait locus (QTL) mapping study of host susceptibility to Kp infection in immunocompetent Collaborative Cross mice (CC). We challenged 328 mice from 73 CC lines intraperitoneally with 104 colony forming units of Kp strain K2. Survival and body weight were monitored for 15 days post challenge. 48 of the CC lines were genotyped with 170,000 SNPs, with which we mapped QTLs.

RESULTS

CC lines differed significantly (P < 0.05) in mean survival time, between 1 to 15 days post infection, and broad sense heritability was 0.45. Distinct QTL were mapped at specific time points during the challenge. A QTL on chromosome 4 was found only on day 2 post infection, and QTL on chromosomes 8 and 18, only on day 8. By using the sequence variations of the eight inbred strain founders of the CC to refine QTL localization we identify several candidate genes.

CONCLUSION

Host susceptibility to Kp is a complex trait, controlled by multiple genetic factors that act sequentially during the course of infection.

The genetics of gene expression in complex mouse crosses as a tool to study the molecular underpinnings of behavior traits

Complex Mus musculus crosses provide increased resolution to examine the relationships between gene expression and behavior. While the advantages are clear, there are numerous analytical and technological concerns that arise from the increased genetic complexity that must be considered. Each of these issues is discussed, providing an initial framework for complex cross study design and planning.

Deconstructing Mus gemischus: advances in understanding ancestry, structure, and variation in the genome of the laboratory mouse

The laboratory mouse is an artificial construct with a complex relationship to its natural ancestors. In 2002, the mouse became the first mammalian model organism with a reference genome. Importantly, the mouse genome sequence was assembled from data on a single inbred laboratory strain, C57BL/6. Several large-scale genetic variant discovery efforts have been conducted, resulting in a catalog of tens of millions of SNPs and structural variants. High-density genotyping arrays covering a subset of those variants have been used to produce hundreds of millions of genotypes in laboratory stocks and a small number of wild mice. These landmark resources now enable us to determine relationships among laboratory mice, assign local ancestry at fine scale, resolve important controversies, and identify a new set of challenges-most importantly, the troubling scarcity of genetic data on the very natural populations from which the laboratory mouse was derived. Our aim with this review is to provide the reader with an historical context for the mouse as a model organism and to explain how practical decisions made in the past have influenced both the architecture of the laboratory mouse genome and the design and execution of current large-scale resources. We also provide examples on how the accomplishments of the past decade can be used by researchers to streamline the use of mice in their experiments and correctly interpret results. Finally, we propose future steps that will enable the mouse community to extend its successes in the decade to come.

Cofilin-1: a modulator of anxiety in mice

The genes involved in conferring susceptibility to anxiety remain obscure. We developed a new method to identify genes at quantitative trait loci (QTLs) in a population of heterogeneous stock mice descended from known progenitor strains. QTLs were partitioned into intervals that can be summarized by a single phylogenetic tree among progenitors and intervals tested for consistency with alleles influencing anxiety at each QTL. By searching for common Gene Ontology functions in candidate genes positioned within those intervals, we identified actin depolymerizing factors (ADFs), including cofilin-1 (Cfl1), as genes involved in regulating anxiety in mice. There was no enrichment for function in the totality of genes under each QTL, indicating the importance of phylogenetic filtering. We confirmed experimentally that forebrain-specific inactivation of Cfl1 decreased anxiety in knockout mice. Our results indicate that similarity of function of mammalian genes can be used to recognize key genetic regulators of anxiety and potentially of other emotional behaviours.

A better prognosis for genetic association studies in mice

Although inbred mouse strains have been the premier model organism used in biomedical research, multiple studies and analyses have indicated that genome-wide association studies (GWAS) cannot be productively performed using inbred mouse strains. However, there is one type of GWAS in mice that has successfully identified the genetic basis for many biomedical traits of interest: haplotype-based computational genetic mapping (HBCGM). Here, we describe how the methodological basis for a HBCGM study significantly differs from that of a conventional murine GWAS, and how an integrative analysis of its output within the context of other 'omic' information can enable genetic discovery. Consideration of these factors will substantially improve the prognosis for the utility of murine genetic association studies for biomedical discovery.

A modifier gene alleviates hypothyroidism-induced hearing impairment in Pou1f1dw dwarf mice

Thyroid hormone has pleiotropic effects on cochlear development, and genomic variation influences the severity of associated hearing deficits. DW/J-Pou1f1dw/dw mutant mice lack pituitary thyrotropin, which causes severe thyroid hormone deficiency and profound hearing impairment. To assess the genetic complexity of protective effects on hypothyroidism-induced hearing impairment, an F1 intercross was generated between DW/J-Pou1f1dw/+ carriers and an inbred strain with excellent hearing derived from Mus castaneus, CAST/EiJ. Approximately 24% of the (DW/J×CAST/EiJ) Pou1f1dw/dw F2 progeny had normal hearing. A genome scan revealed a locus on chromosome 2, named modifier of dw hearing, or Mdwh, that rescues hearing despite persistent hypothyroidism. This chromosomal region contains the modifier of tubby hearing 1 (Moth1) locus that encodes a protective allele of the microtubule-associated protein MTAP1A. DW/J-Pou1f1dw/+ carriers were crossed with the AKR strain, which also carries a protective allele of Mtap1a, and we found that AKR is not protective for hearing in the (DW/J×AKR) Pou1f1dw/dw F2 progeny. Thus, protective alleles of Mtap1a are not sufficient to rescue DW/J-Pou1f1dw/dw hearing. We expect that identification of protective modifiers will enhance our understanding of the mechanisms of hypothyroidism-induced hearing impairment.

Distribution of parental genome blocks in recombinant inbred lines

We consider recombinant inbred lines obtained by crossing two given homozygous parents and then applying multiple generations of self-crossings or full-sib matings. The chromosomal content of any such line forms a mosaic of blocks, each alternatively inherited identically by descent from one of the parents. Quantifying the statistical properties of such mosaic genomes has remained an open challenge for many years. Here, we solve this problem by taking a continuous chromosome picture and assuming crossovers to be noninterfering. Using a continuous-time random walk framework and Markov chain theory, we determine the statistical properties of these identical-by-descent blocks. We find that successive block lengths are only very slightly correlated. Furthermore, the blocks on the ends of chromosomes are larger on average than the others, a feature understandable from the nonexponential distribution of block lengths.

Collaborative Cross mice and their power to map host susceptibility to Aspergillus fumigatus infection

The Collaborative Cross (CC) is a genetic reference panel of recombinant inbred lines of mice, designed for the dissection of complex traits and gene networks. Each line is independently descended from eight genetically diverse founder strains such that the genomes of the CC lines, once fully inbred, are fine-grained homozygous mosaics of the founder haplotypes. We present an analysis of 120 CC lines, from a cohort of the CC bred at Tel Aviv University in collaboration with the University of Oxford, which at the time of this study were between the sixth and 12th generations of inbreeding and substantially homozygous at 170,000 SNPs. We show how CC genomes decompose into mosaics, and we identify loci that carry a deficiency or excess of a founder, many being deficient for the wild-derived strains WSB/EiJ and PWK/PhJ. We phenotyped 371 mice from 66 CC lines for a susceptibility to Aspergillus fumigatus infection. The survival time after infection varied significantly between CC lines. Quantitative trait locus (QTL) mapping identified genome-wide significant QTLs on chromosomes 2, 3, 8, 10 (two QTLs), 15, and 18. Simulations show that QTL mapping resolution (the median distance between the QTL peak and true location) varied between 0.47 and 1.18 Mb. Most of the QTLs involved contrasts between wild-derived founder strains and therefore would not segregate between classical inbred strains. Use of variation data from the genomes of the CC founder strains refined these QTLs further and suggested several candidate genes. These results support the use of the CC for dissecting complex traits.

Commercially available outbred mice for genome-wide association studies

Genome-wide association studies using commercially available outbred mice can detect genes involved in phenotypes of biomedical interest. Useful populations need high-frequency alleles to ensure high power to detect quantitative trait loci (QTLs), low linkage disequilibrium between markers to obtain accurate mapping resolution, and an absence of population structure to prevent false positive associations. We surveyed 66 colonies for inbreeding, genetic diversity, and linkage disequilibrium, and we demonstrate that some have haplotype blocks of less than 100 Kb, enabling gene-level mapping resolution. The same alleles contribute to variation in different colonies, so that when mapping progress stalls in one, another can be used in its stead. Colonies are genetically diverse: 45% of the total genetic variation is attributable to differences between colonies. However, quantitative differences in allele frequencies, rather than the existence of private alleles, are responsible for these population differences. The colonies derive from a limited pool of ancestral haplotypes resembling those found in inbred strains: over 95% of sequence variants segregating in outbred populations are found in inbred strains. Consequently it is possible to impute the sequence of any mouse from a dense SNP map combined with inbred strain sequence data, which opens up the possibility of cataloguing and testing all variants for association, a situation that has so far eluded studies in completely outbred populations. We demonstrate the colonies' potential by identifying a deletion in the promoter of H2-Ea as the molecular change that strongly contributes to setting the ratio of CD4+ and CD8+ lymphocytes.

Mapping quantitative traits and strategies to find quantitative trait genes

In 1999 a meeting took place at the Jackson Laboratory, a large mouse research centre in Bar Harbor, Maine, to consider the value of systematically collecting phenotypes on inbred strains of mice (Paigen and Eppig (2000) [1]). The group concluded that cataloguing the extensive phenotypic diversity present among laboratory mice, and in particular providing the research community with data from cohorts of animals, phenotyped according to standardized protocols, was essential if we were to take advantage of the possibilities of mouse genetics. Beginning with the collection of basic physiological, biochemical and behavioral data on nine commonly used inbred strains, the project has expanded so that by the beginning of 2010 data for 178 strains had been collected, with 105 phenotype projects yielding over 2000 different measurements (Bogue et al. (2007) [2].

Strain screen and haplotype association mapping of wheel running in inbred mouse strains

Previous genetic association studies of physical activity, in both animal and human models, have been limited in number of subjects and genetically homozygous strains used as well as number of genomic markers available for analysis. Expansion of the available mouse physical activity strain screens and the recently published dense single-nucleotide polymorphism (SNP) map of the mouse genome (approximately 8.3 million SNPs) and associated statistical methods allowed us to construct a more generalizable map of the quantitative trait loci (QTL) associated with physical activity. Specifically, we measured wheel running activity in male and female mice (average age 9 wk) in 41 inbred strains and used activity data from 38 of these strains in a haplotype association mapping analysis to determine QTL associated with activity. As seen previously, there was a large range of activity patterns among the strains, with the highest and lowest strains differing significantly in daily distance run (27.4-fold), duration of activity (23.6-fold), and speed (2.9-fold). On a daily basis, female mice ran further (24%), longer (13%), and faster (11%). Twelve QTL were identified, with three (on Chr. 12, 18, and 19) in both male and female mice, five specific to males, and four specific to females. Eight of the 12 QTL, including the 3 general QTL found for both sexes, fell into intergenic areas. The results of this study further support the findings of a moderate to high heritability of physical activity and add general genomic areas applicable to a large number of mouse strains that can be further mined for candidate genes associated with regulation of physical activity. Additionally, results suggest that potential genetic mechanisms arising from traditional noncoding regions of the genome may be involved in regulation of physical activity.

Identifying quantitative trait loci (QTLs) and genes (QTGs) for alcohol-related phenotypes in mice

Alcoholism is a complex clinical disorder with genetic and environmental contributions. Although no animal model duplicates alcoholism, models for specific factors, such as the withdrawal syndrome, are useful to identify potential genetic determinants of liability in humans. Murine models have been invaluable to identify quantitative trait loci (QTLs) that influence a variety of alcohol responses. However, the QTL regions are typically large, at least initially, and contain numerous genes, making identification of the causal quantitative trait gene(s) (QTGs) challenging. Here, we present QTG identification strategies currently used in the field of alcohol genetics and discuss relevance to alcoholic human populations.

Haplotype association mapping in mice

Haplotype Association Mapping (HAM) is a novel phenotype-driven approach to identify genetic loci and was originally developed for mice. This method, which is similar to Genome-Wide Association (GWA) studies in humans, looks for associations between the phenotype and the haplotypes of mouse inbred strains, treating inbred strains as individuals. Although this approach is still in development, we review the current literature, present the different methods and applications that are in use, and provide a glimpse of what is to come in the near future.

An integrative genetics approach to identify candidate genes regulating BMD: combining linkage, gene expression, and association

Numerous quantitative trait loci (QTLs) affecting bone traits have been identified in the mouse; however, few of the underlying genes have been discovered. To improve the process of transitioning from QTL to gene, we describe an integrative genetics approach, which combines linkage analysis, expression QTL (eQTL) mapping, causality modeling, and genetic association in outbred mice. In C57BL/6J x C3H/HeJ (BXH) F(2) mice, nine QTLs regulating femoral BMD were identified. To select candidate genes from within each QTL region, microarray gene expression profiles from individual F(2) mice were used to identify 148 genes whose expression was correlated with BMD and regulated by local eQTLs. Many of the genes that were the most highly correlated with BMD have been previously shown to modulate bone mass or skeletal development. Candidates were further prioritized by determining whether their expression was predicted to underlie variation in BMD. Using network edge orienting (NEO), a causality modeling algorithm, 18 of the 148 candidates were predicted to be causally related to differences in BMD. To fine-map QTLs, markers in outbred MF1 mice were tested for association with BMD. Three chromosome 11 SNPs were identified that were associated with BMD within the Bmd11 QTL. Finally, our approach provides strong support for Wnt9a, Rasd1, or both underlying Bmd11. Integration of multiple genetic and genomic data sets can substantially improve the efficiency of QTL fine-mapping and candidate gene identification.

Human-mouse quantitative trait locus concordance and the dissection of a human neuroticism locus

BACKGROUND

Exploiting synteny between mouse and human disease loci has been proposed as a cost-effective method for the identification of human susceptibility genes. Here we explore its utility in an analysis of a human personality trait, neuroticism, which can be modeled in mice by tests of emotionality. We investigated a mouse emotionality locus on chromosome 1 that contains no annotated genes but abuts four regulators of G protein signaling, one of which (rgs2) has been previously identified as a quantitative trait gene for emotionality. This locus is syntenic with a human region that has been consistently implicated in the genetic aetiology of neuroticism.

METHODS

The functional candidacy of 29 murine sequence variants was tested by a combination of gel shift and transient transfection assays. Murine sequences that contained functional variants and exhibited significant cross-species conservation were prioritized for investigation in humans. Genetic association with neuroticism was tested in 1869 high and 2032 low unrelated individuals scored for neuroticism, selected from the extremes of 88,141 people from southwest England.

RESULTS

Fifteen sequence variants contributed to variation in the expression of rgs18, the gene lying at the edge of the quantitative trait loci (QTL) interval. There was no evidence of association between neuroticism and single nucleotide polymorphisms (SNPs) lying in the human regions homologous to those of mouse functional variants. One SNP, rs6428058, in a region of sequence conservation 644 kb upstream of RGS18, showed significant association (p = .000631).

CONCLUSIONS

It is unlikely that a single variant is responsible for the mouse emotionality locus on chromosome 1. This level of underlying genetic complexity means that although cross-species QTL concordance may be invaluable for the identification of human disease loci, it is unlikely to be as informative in the identification of human disease-causing variants.

Nucleotide variation in wild and inbred mice

The house mouse is a well-established model organism, particularly for studying the genetics of complex traits. However, most studies of mice use classical inbred strains, whose genomes derive from multiple species. Relatively little is known about the distribution of genetic variation among these species or how variation among strains relates to variation in the wild. We sequenced intronic regions of five X-linked loci in large samples of wild Mus domesticus and M. musculus, and we found low levels of nucleotide diversity in both species. We compared these data to published data from short portions of six X-linked and 18 autosomal loci in wild mice. We estimate that M. domesticus and M. musculus diverged <500,000 years ago. Consistent with this recent divergence, some gene genealogies were reciprocally monophyletic between these species, while others were paraphyletic or polyphyletic. In general, the X chromosome was more differentiated than the autosomes. We resequenced classical inbred strains for all 29 loci and found that inbred strains contain only a small amount of the genetic variation seen in wild mice. Notably, the X chromosome contains proportionately less variation among inbred strains than do the autosomes. Moreover, variation among inbred strains derives from differences between species as well as from differences within species, and these proportions differ in different genomic regions. Wild mice thus provide a reservoir of additional genetic variation that may be useful for mapping studies. Together these results suggest that wild mice will be a valuable complement to laboratory strains for studying the genetics of complex traits.

An imputed genotype resource for the laboratory mouse

We have created a high-density SNP resource encompassing 7.87 million polymorphic loci across 49 inbred mouse strains of the laboratory mouse by combining data available from public databases and training a hidden Markov model to impute missing genotypes in the combined data. The strong linkage disequilibrium found in dense sets of SNP markers in the laboratory mouse provides the basis for accurate imputation. Using genotypes from eight independent SNP resources, we empirically validated the quality of the imputed genotypes and demonstrated that they are highly reliable for most inbred strains. The imputed SNP resource will be useful for studies of natural variation and complex traits. It will facilitate association study designs by providing high-density SNP genotypes for large numbers of mouse strains. We anticipate that this resource will continue to evolve as new genotype data become available for laboratory mouse strains. The data are available for bulk download or query at http://cgd.jax.org /.

Fine haplotype structure of a chromosome 17 region in the laboratory and wild mouse

Extensive linkage disequilibrium among classical laboratory strains represents an obstacle in the high-resolution haplotype mapping of mouse quantitative trait loci (QTL). To determine the potential of wild-derived mouse strains for fine QTL mapping, we constructed a haplotype map of a 250-kb region of the t-complex on chromosome 17 containing the Hybrid sterility 1 (Hst1) gene. We resequenced 33 loci from up to 80 chromosomes of five mouse (sub)species. Trans-species single-nucleotide polymorphisms (SNPs) were rare between Mus m. musculus (Mmmu) and Mus m. domesticus (Mmd). The haplotypes in Mmmu and Mmd differed and therefore strains from these subspecies should not be combined for haplotype-associated mapping. The haplotypes of t-chromosomes differed from all non-t Mmmu and Mmd haplotypes. Half of the SNPs and SN indels but only one of seven longer rearrangements found in classical laboratory strains were useful for haplotype mapping in the wild-derived M. m. domesticus. The largest Mmd haplotype block contained three genes of a highly conserved synteny. The lengths of the haplotype blocks deduced from 36 domesticus chromosomes were in tens of kilobases, suggesting that the wild-derived Mmd strains are suitable for fine interval-specific mapping.

Quantitative trait locus analysis: multiple cross and heterogeneous stock mapping

Until well into the 1990s, both preclinical and clinical research focused on finding "the" gene for human diseases, including alcoholism. This focus was reinforced by the emergence of technologies to either inactivate (i.e., knock out) a gene or add extra copies of an existing gene in a living organism, which clearly demonstrated that over- or underexpressing a single gene could have a profound effect on behavior. However, a small but vocal group of scientists, including many alcohol researchers, argued that behaviors, including alcohol-related behaviors, were complex traits and therefore no one gene likely would have a large effect. This view was consistent with a large body of genetic research conducted in plants and fruit flies (e.g., Paterson et al. 1988) indicating that, for example, even a presumably simple characteristic, such as the size of a tomato, was determined by several genes. However, it was difficult to convince the scientific community that, in terms of its genetic determination, behavior was similar to the size of a tomato. Only with the advent of new genetic tools did it become possible to prove that many different genes contribute to complex behavioral characteristics.

From genetics to mechanism of disease liability

With each advance in genomic technology, new statistical methods have regularly emerged to test genetic hypotheses in complex inheritance, as evidenced throughout this book. Notwithstanding the approach used, the greatest challenge in the genetics of complex traits remains the identification of the gene(s) and the molecular variant(s) accounting for a genetic inference based on statistical testing. We take the example of quantitative trait locus (QTL) mapping for blood pressure (BP) and related phenotypes in rodents to review the current landscape. Traditional approaches to refined mapping are typically hampered by the small effect and the small proportion of the variance attached to individual QTLs. The alternative of functional screens in intact animals, whether by chemical mutagenesis or gene targeting, remains a daunting undertaking. Such limitations account for the slow progress to date of inferences from QTL to gene(s). We select a QTL for differential sodium sensitivity between two mouse inbred lines to propose an approach that can be used in relatively large genomic regions (1) by optimizing the selection of candidate genes and (2) by subjecting such genes to high-throughput functional screens. While this is still work in progress, we think it abundantly illustrates what is ahead of us in delineating genetic variation that underlie complex disease.

Glutamate receptor metabotropic 7 is cis-regulated in the mouse brain and modulates alcohol drinking

Alcoholism is a heritable disease that afflicts about 8% of the adult population. Its development and symptoms, such as craving, loss of control, physical dependence, and tolerance, have been linked to changes in mesolimbic, mesocortical neurotransmitter systems utilizing biogenic amines, GABA, and glutamate. Identification of genes predisposing to alcoholism, or to alcohol-related behaviors in animal models, has been elusive because of variable interactions of multiple genes with relatively small individual effect size and sensitivity of the predisposing genotype to lifestyle and environmental factors. Here, using near-isogenic advanced animal models with reduced genetic background interactions, we integrate gene mapping and gene mRNA expression data in segregating and congenic mice and identify glutamate receptor metabotropic 7 (Grm7) as a cis-regulated gene for alcohol consumption. Traditionally, the mesoaccumbal dopamine reward hypothesis of addiction and the role of the ionotropic glutamate receptors have been emphasized. Our results lend support to an emerging direction of research on the role of metabotropic glutamate receptors in alcoholism and drug addiction. These data suggest for the first time that Grm7 is a risk factor for alcohol drinking and a new target in addiction therapy.

Meiotic crossover hotspots contained in haplotype block boundaries of the mouse genome

Fertility requires successful chromosome segregation in meiosis, which in most sexual organisms depends on the formation of appropriately placed crossovers. The nonrandom genome-wide distributions of meiotic recombination events have been examined at the molecular level experimentally in yeast and by inference from linkage disequilibrium patterns in humans. Thus far, no method has existed for pinpointing sites of crossing-over on a genome-wide scale in an experimentally tractable animal whose genome size and complexity models that of humans. Here, we present a genomic approach to identify mouse crossover hotspots, based on targeting haplotype block boundaries. This represents a previously undescribed method potentially applicable to large-scale mouse hotspot identification. Using this method, we have successfully predicted the location of two previously uncharacterized crossover hotspots in male mice. As increasing amounts of single-nucleotide polymorphism data emerge, this approach will be useful for investigating the recombination landscape of the mouse genome.

A sequence-based variation map of 8.27 million SNPs in inbred mouse strains

A dense map of genetic variation in the laboratory mouse genome will provide insights into the evolutionary history of the species and lead to an improved understanding of the relationship between inter-strain genotypic and phenotypic differences. Here we resequence the genomes of four wild-derived and eleven classical strains. We identify 8.27 million high-quality single nucleotide polymorphisms (SNPs) densely distributed across the genome, and determine the locations of the high (divergent subspecies ancestry) and low (common subspecies ancestry) SNP-rate intervals for every pairwise combination of classical strains. Using these data, we generate a genome-wide haplotype map containing 40,898 segments, each with an average of three distinct ancestral haplotypes. For the haplotypes in the classical strains that are unequivocally assigned ancestry, the genetic contributions of the Mus musculus subspecies--M. m. domesticus, M. m. musculus, M. m. castaneus and the hybrid M. m. molossinus--are 68%, 6%, 3% and 10%, respectively; the remaining 13% of haplotypes are of unknown ancestral origin. The considerable regional redundancy of the SNP data will facilitate imputation of the majority of these genotypes in less-densely typed classical inbred strains to provide a complete view of variation in additional strains.

On the subspecific origin of the laboratory mouse

The genome of the laboratory mouse is thought to be a mosaic of regions with distinct subspecific origins. We have developed a high-resolution map of the origin of the laboratory mouse by generating 25,400 phylogenetic trees at 100-kb intervals spanning the genome. On average, 92% of the genome is of Mus musculus domesticus origin, and the distribution of diversity is markedly nonrandom among the chromosomes. There are large regions of extremely low diversity, which represent blind spots for studies of natural variation and complex traits, and hot spots of diversity. In contrast with the mosaic model, we found that most of the genome has intermediate levels of variation of intrasubspecific origin. Finally, mouse strains derived from the wild that are supposed to represent different mouse subspecies show substantial intersubspecific introgression, which has strong implications for evolutionary studies that assume these are pure representatives of a given subspecies.

Mesencephalic dopamine neuron number and tyrosine hydroxylase content: Genetic control and candidate genes

The mesotelencephalic dopamine system shows substantial genetic variation which fundamentally affects normal and pathological behaviors related to motor function, motivation, and learning. Our earlier radioenzyme assay studies demonstrated significantly higher activity of tyrosine hydroxylase (TH), the first and rate limiting enzyme in the biosynthesis of catecholamine neurotransmitters, in the substantia nigra-ventral tegmental area of BALB/cJ mice in comparison with that of C57BL/6ByJ mice. Here, using quantitative immunoblotting and immunocytochemistry, we tested the hypothesis that mesencephalic TH protein content and number of nigral TH-positive neurons show strain-dependent differences in C57BL/6ByJ and BALB/cJ parallel to those observed in the TH activity studies. Immunoblotting experiments detected significantly higher mesencephalic TH protein content in BALB/cJ in comparison to C57BL/6ByJ (P<0.05). Immunocytochemical studies demonstrated that the number of TH-positive cells in substantia nigra was 31.3% higher in BALB/cJ than that in C57BL/6ByJ (P<0.01), while the average dopamine neuron volume was not significantly different. In a search for candidate genes that modulate TH content and the size of mesencephalic dopamine neuron populations we also studied near-isogenic mouse sublines derived from the C57BL/6ByJ and BALB/cJ progenitor strains. A whole-genome scan with 768 single nucleotide polymorphism markers indicated that two sublines, C4A6/N and C4A6/B, were genetically very similar (98.3%). We found significantly higher mesencephalic TH protein content in C4A6/B in comparison to C4A6/N (P=0.01), and a tendency for higher number of dopamine neurons in the substantia nigra in C4A6/B in comparison to C4A6/N, which, however, did not reach statistical significance. To identify the genetic source of the TH content difference we analyzed the single nucleotide polymorphism (SNP) genotype data of the whole-genome scan, and detected two small differential chromosome segments on chr. 13 and chr. 14. Microarray gene expression studies and bioinformatic analysis of the two differential regions implicated two cis-regulated genes (Spock1 and Cxcl14, chr. 13), and two growth factor genes [bone morphogenetic protein 6 (Bmp6) (chr. 13), and fibroblast growth factor 14 (Fgf14) (chr. 14)]. Taken together, the results suggest that (1) nigral dopamine neuron number and TH protein content may be genetically associated but further studies are needed to establish unequivocally this linkage, and (2) Spock1, Cxcl14, Bmp6, and Fgf14 are novel candidates for modulating the expression and maintenance of TH content in mesencephalic dopamine neurons in vivo.

Multiple trait measurements in 43 inbred mouse strains capture the phenotypic diversity characteristic of human populations

The breadth of genetic and phenotypic variation among inbred strains is often underappreciated because assessments include only a limited number of strains. Evaluation of a larger collection of inbred strains provides not only a greater understanding of this variation but collectively mimics much of the variation observed in human populations. We used a high-throughput phenotyping protocol to measure females and males of 43 inbred strains for body composition (weight, fat, lean tissue mass, and bone mineral density), plasma triglycerides, high-density lipoprotein and total cholesterol, glucose, insulin, and leptin levels while mice consumed a high-fat, high-cholesterol diet. Mice were fed a chow diet until they were 6-8 wk old and then fed the high-fat diet for an additional 18 wk. As expected, broad phenotypic diversity was observed among these strains. Significant variation between the sexes was also observed for most traits measured. Additionally, the response to the high-fat diet differed considerably among many strains. By the testing of such a large set of inbred strains for many traits, multiple phenotypes can be considered simultaneously and thereby aid in the selection of certain inbred strains as models for complex human diseases. These data are publicly available in the web-accessible Mouse Phenome Database (http://www.jax.org/phenome), an effort established to promote systematic characterization of biochemical and behavioral phenotypes of commonly used and genetically diverse inbred mouse strains. Data generated by this effort builds on the value of inbred mouse strains as a powerful tool for biomedical research.

Prospects for association mapping in classical inbred mouse strains

The collection of classical inbred mouse strains displays heritable variation in a large number of complex traits. Many generations of historical recombination have contributed to the panel of classical strain genomes, raising the possibility that quantitative trait loci could be located with high resolution by correlating strain genotypes and phenotypes. Although this association mapping framework has been successful in several empirical applications, its expected performance remains unclear. We used computer simulations based on a publicly available, dense single-nucleotide polymorphism (SNP) map to measure the power and false-positive rate of association mapping on a genomic scale across 30 commonly used classical inbred strains. Expected power is (i) often low for phenotypic effect sizes that are realistic for complex traits, (ii) highly variable across the genome, and (iii) correlated with linkage disequilibrium, aspects of the allele frequency distribution, and haplotype characteristics, as predicted by theory. Simulations also demonstrate clear potential for spurious associations to be generated by unequal relatedness among the strains. These findings suggest that association mapping in the classical strains is best applied in combination with other procedures, such as QTL mapping.

The genetic landscape of type 2 diabetes in mice

Inbred mouse strains provide genetic diversity comparable to that of the human population. Like humans, mice have a wide range of diabetes-related phenotypes. The inbred mouse strains differ in the response of their critical physiological functions, such as insulin sensitivity, insulin secretion, beta-cell proliferation and survival, and fuel partitioning, to diet and obesity. Most of the critical genes underlying these differences have not been identified, although many loci have been mapped. The dramatic improvements in genomic and bioinformatics resources are accelerating the pace of gene discovery. This review describes how mouse genetics can be used to discover diabetes-related genes, summarizes how the mouse strains differ in their diabetes-related phenotypes, and describes several examples of how loci identified in the mouse may directly relate to human diabetes.

Quantitative Trait Locus Mapping for Acute Functional Tolerance to Ethanol in the L × S Recombinant Inbred Panel

BACKGROUND

Acute functional tolerance (AFT) develops shortly after ethanol administration, and is determined as the change in brain or blood ethanol concentration (BEC) measured at 2 behavioral or physiological endpoints. Acute functional tolerance studies in some rodent strains support a long-held hypothesis that more sensitive strains develop more within-session tolerance. We used the new, 74-strain L x S recombinant inbred (RI) panel, developed from inbred long-sleep (ILS) and inbred short-sleep (ISS) strains, to revisit this hypothesis and to map quantitative trait loci (QTLs) for AFT. We report replication of QTL regions reported by earlier studies of AFT and preliminary application of a coarse single nucleotide polymorphism map analysis to limit QTL intervals for subsequent candidate gene hypotheses.

METHODS

Acute functional tolerance was assayed using a test of ataxia: loss and regain of balance on a stationary wooden dowel. Following an initial dose of 1.75 g/kg, BEC was measured at initial loss (BEC(0)) and regain of balance (BEC(1)). A second injection (2.0 g/kg) was administered and blood taken at the second regain of balance (BEC(2)). Acute functional tolerance was calculated as a difference score in 2 ways: (1) between BEC at the 2 successive regains of balance (AFT(1)), or (2) as the difference in BEC at final regain and at initial loss of balance (AFT(2)). We mapped QTLs for BEC(0), a measure of initial sensitivity, and both AFT scores.

RESULTS

All 4 parental strains (LS, SS, ILS, and ISS) developed tolerance, replicating previous published reports. There were significant sex effects for 3 of these strains. The L x S panel showed a 128-fold range in tolerance, with a few strains showing negative tolerance (sensitization). The ISS surpassed the next highest RI strain by 55% and was more than 4 times greater than SS. Heritability estimates for both AFT measures were close to 0.25 for both sexes. One significant QTL accounting for approximately 18% of phenotypic variance (V(P)), on chromosome 12 (AFT(1)), and 1 suggestive QTL (16% V(P)), on chromosome 16 (AFT(2)), were identified. These QTLs replicated regions reported in other studies. A multiple QTL model incorporating the effects of all significant interacting QTLs was developed, explaining almost 60% of V(P). The chromosome 12 region was further investigated by haplotype analysis, which identified many nonpolymorphic regions within the confidence interval, and possible candidate genes in the polymorphic regions.

CONCLUSIONS

Both SS and ISS developed greater AFT, assessed by both methods, than LS and ILS; this difference was significant in virtually all sex by strain comparisons. In the L x S RI, there was no correlation between initial sensitivity, measured by BEC at initial loss of balance, and either measure of AFT, on a stationary dowel. These results indicate that in this model system, initial sensitivity does not predict tolerance. Several QTLs for tolerance were identified; candidates in the narrowed chromosome 12 region, which has been reported in 2 other mapping studies, merit additional study.

Fine mapping of collagen-induced arthritis quantitative trait loci in an advanced intercross line

The generation of advanced intercross lines (AIL) is a powerful approach for high-resolution fine mapping of quantitative trait loci (QTLs), because they accumulate much more recombination events compared with conventional F2 intercross and N2 backcross. However, the application of this approach is severely hampered by the requirements of excessive resources to maintain such crosses, i.e., in terms of animal care, space, and time. Therefore, in this study, we produced an AIL to fine map collagen-induced arthritis (CIA) QTLs using comparatively limited resources. We used only 308 (DBA/1 x FVB/N)F11/12 AIL mice to refine QTLs controlling the severity and onset of arthritis as well as the Ab response and T cell subset in CIA, namely Cia2, Cia27, and Trmq3. These QTLs were originally identified in (DBA/1 x FVB/N)F2 progeny. The confidence intervals of the three QTLs were refined from 40, 43, and 48 Mb to 12, 4.1, and 12 Mb, respectively. The data were complemented by the use of another QTL fine-mapping approach, haplotype analysis, to further refine Cia2 into a 2-Mb genomic region. To aid in the search for candidate genes for the QTLs, genome-wide expression profiling was performed to identify strain-specific differentially expressed genes within the confidence intervals. Of the 1396 strain-specific differentially expressed genes, 3, 3, and 12 genes were within the support intervals of the Cia2, Cia27, and Trmq3, respectively. In addition, this study revealed that Cia27 and Trmq3 controlling anti-CII IgG2a Ab and CD4:CD8 T cell ratio, respectively, also regulated CIA clinical phenotypes.

Searching the genomes of inbred mouse strains for incompatibilities that reproductively isolate their wild relatives

Identification of the genes that underlie reproductive isolation provides important insights into the process of speciation. According to the Dobzhansky-Muller model, these genes suffer disrupted interactions in hybrids due to independent divergence in separate populations. In hybrid populations, natural selection acts to remove the deleterious heterospecific combinations that cause these functional disruptions. When selection is strong, this process can maintain multilocus associations, primarily between conspecific alleles, providing a signature that can be used to locate incompatibilities. We applied this logic to populations of house mice that were formed by hybridization involving two species that show partial reproductive isolation, Mus domesticus and Mus musculus. Using molecular markers likely to be informative about species ancestry, we scanned the genomes of 1) classical inbred strains and 2) recombinant inbred lines for pairs of loci that showed extreme linkage disequilibria. By using the same set of markers, we identified a list of locus pairs that displayed similar patterns in both scans. These genomic regions may contain genes that contribute to reproductive isolation between M. domesticus and M. musculus. This hypothesis can now be tested using laboratory crosses and surveys of introgression in the wild.

Further studies on using multiple-cross mapping (MCM) to map quantitative trait loci

We have completed whole-genome scans for quantitative trait loci (QTLs) associated with acute ethanol-induced activation in the six F(2) intercrosses that can be formed from the C57BL/6J (B6), DBA/2J (D2) , BALB/cJ (C), and LP/J (LP) inbred strains. The goal was to test the hypothesis that given the relatively simple structure of the laboratory mouse genome, the same QTLs will be detected in multiple crosses which in turn will provide support for the strategy of multiple-cross mapping (MCM). QTLs with LOD scores greater than 4 were detected on Chrs 1, 2, 3, 8, 9, 13, 14, and 16. Only for the QTL on distal Chr 1 was there convincing evidence that the same or at least a very similar QTL was detected in multiple crosses. We also mapped the Chr 2 QTL directly in heterogeneous stock (HS) animals derived from the four inbred strains. At G(19) the QTL was mapped to an approximately 3-Mbp interval and this interval was associated with a haplotype block with a largely biallelic structure: B6-L:C-D2. We conclude that mapping in HS animals not only provides significantly greater QTL resolution, at least in some cases it provides significantly more information about the QTL haplotype structure.

Identification of diabetes susceptibility loci in db mice by combined quantitative trait loci analysis and haplotype mapping

To identify the disease-susceptibility genes of type 2 diabetes, we performed quantitative trait loci (QTL) analysis in F(2) populations generated from a BKS.Cg-m+/+Lepr(db) and C3H/HeJ intercross, taking advantage of genetically determined obesity and diabetes traits associated with the db gene. A genome-wide scan in the F(2) populations divided by sex and db genotypes identified 14 QTLs in total and 3 major QTLs on chromosome (Chr) 3 (LOD 5.78) for fat pad weight, Chr 15 (LOD 6.64) for body weight, and Chr 16 (LOD 8.15) for blood glucose concentrations. A linear-model-based genome scan using interactive covariates allowed us to consider sex- or sex-by db-specific effects of each locus. For the most significant QTL on Chr 16, the high-resolution haplotype comparison between BKS and C3H strains reduced the critical QTL interval from 20 to 4.6 Mb by excluding shared haplotype regions and identified 11 nonsynonymous single-nucleotide polymorphisms in six candidate genes.