Panel Of Inbred Strains Research Articles

Into the Wild: A novel wild-derived inbred strain resource expands the genomic and phenotypic diversity of laboratory mouse models.

The laboratory mouse has served as the premier animal model system for both basic and preclinical investigations for over a century. However, laboratory mice capture only a subset of the genetic variation found in wild mouse populations, ultimately limiting the potential of classical inbred strains to uncover phenotype-associated variants and pathways. Wild mouse populations are reservoirs of genetic diversity that could facilitate the discovery of new functional and disease-associated alleles, but the scarcity of commercially available, well-characterized wild mouse strains limits their broader adoption in biomedical research. To overcome this barrier, we have recently developed, sequenced, and phenotyped a set of 11 inbred strains derived from wild-caught Mus musculus domesticus. Each of these "Nachman strains" immortalizes a unique wild haplotype sampled from one of five environmentally distinct locations across North and South America. Whole genome sequence analysis reveals that each strain carries between 4.73-6.54 million single nucleotide differences relative to the GRCm39 mouse reference, with 42.5% of variants in the Nachman strain genomes absent from current classical inbred mouse strain panels. We phenotyped the Nachman strains on a customized pipeline to assess the scope of disease-relevant neurobehavioral, biochemical, physiological, metabolic, and morphological trait variation. The Nachman strains exhibit significant inter-strain variation in >90% of 1119 surveyed traits and expand the range of phenotypic diversity captured in classical inbred strain panels. These novel wild-derived inbred mouse strain resources are set to empower new discoveries in both basic and preclinical research.

Collaborative cross strain CC011/UncJ as a novel mouse model of T2-high, severe asthma

Among asthmatics, there is significant heterogeneity in the clinical presentation and underlying pathophysiological mechanisms, leading to the recognition of multiple disease endotypes (e.g., T2-high vs. T2-low). This heterogeneity extends to severe asthmatics, who may struggle to control symptoms even with high-dose corticosteroid treatment and other therapies. However, there are limited mouse models available to model the spectrum of severe asthma endotypes. We sought to identify a new mouse model of severe asthma by first examining responses to chronic allergen exposure among strains from the Collaborative Cross (CC) mouse genetics reference population, which contains greater genetic diversity than other inbred strain panels previously used for models of asthma. Mice from five CC strains and the often-used classical inbred strain BALB/cJ were chronically exposed to house dust mite (HDM) allergen for five weeks followed by measurements of airway inflammation. CC strain CC011/UncJ (CC011) exhibited extreme responses to HDM including high levels of airway eosinophilia, elevated lung resistance, and extensive airway wall remodeling, and even fatalities among ~ 50% of mice prior to study completion. Compared to BALB/cJ mice, CC011 mice had stronger Th2-mediated airway responses demonstrated by significantly elevated total and HDM-specific IgE and increased Th2 cytokines during tests of antigen recall, but not enhanced ILC2 activation. Airway eosinophilia in CC011 mice was completely dependent upon CD4+ T-cells. Notably, we also found that airway eosinophilia in CC011 mice was resistant to dexamethasone steroid treatment. Thus, the CC011 strain provides a new mouse model of T2-high, severe asthma driven by natural genetic variation likely acting through CD4+ T-cells. Future studies aimed at determining the genetic basis of this phenotype will provide new insights into mechanisms underlying severe asthma.

Natural genetic variation in a dopamine receptor is associated with variation in female fertility in Drosophila melanogaster.

Fertility is a major component of fitness but its genetic architecture remains poorly understood. Using a full diallel cross of 50 Drosophila Genetic Reference Panel inbred lines with whole genome sequences, we found substantial genetic variation in fertility largely attributable to females. We mapped genes associated with variation in female fertility by genome-wide association analysis of common variants in the fly genome. Validation of candidate genes by RNAi knockdown confirmed the role of the dopamine 2-like receptor (Dop2R) in promoting egg laying. We replicated the Dop2R effect in an independently collected productivity dataset and showed that the effect of the Dop2R variant was mediated in part by regulatory gene expression variation. This study demonstrates the strong potential of genome-wide association analysis in this diverse panel of inbred strains and subsequent functional analyses for understanding the genetic architecture of fitness traits.

A Mouse Systems Genetics Approach Reveals Common and Uncommon Genetic Modifiers of Hepatic Lysosomal Enzyme Activities and Glycosphingolipids.

Identification of genetic modulators of lysosomal enzyme activities and glycosphingolipids (GSLs) may facilitate the development of therapeutics for diseases in which they participate, including Lysosomal Storage Disorders (LSDs). To this end, we used a systems genetics approach: we measured 11 hepatic lysosomal enzymes and many of their natural substrates (GSLs), followed by modifier gene mapping by GWAS and transcriptomics associations in a panel of inbred strains. Unexpectedly, most GSLs showed no association between their levels and the enzyme activity that catabolizes them. Genomic mapping identified 30 shared predicted modifier genes between the enzymes and GSLs, which are clustered in three pathways and are associated with other diseases. Surprisingly, they are regulated by ten common transcription factors, and their majority by miRNA-340p. In conclusion, we have identified novel regulators of GSL metabolism, which may serve as therapeutic targets for LSDs and may suggest the involvement of GSL metabolism in other pathologies.

Mouse genome-wide association and systems genetics identifies Lhfp as a regulator of bone mass

Bone mineral density (BMD) is a strong predictor of osteoporotic fracture. It is also one of the most heritable disease-associated quantitative traits. As a result, there has been considerable effort focused on dissecting its genetic basis. Here, we performed a genome-wide association study (GWAS) in a panel of inbred strains to identify associations influencing BMD. This analysis identified a significant (P = 3.1 x 10−12) BMD locus on Chromosome 3@52.5 Mbp that replicated in two separate inbred strain panels and overlapped a BMD quantitative trait locus (QTL) previously identified in a F2 intercross. The association mapped to a 300 Kbp region containing four genes; Gm2447, Gm20750, Cog6, and Lhfp. Further analysis found that Lipoma HMGIC Fusion Partner (Lhfp) was highly expressed in bone and osteoblasts. Furthermore, its expression was regulated by a local expression QTL (eQTL), which overlapped the BMD association. A co-expression network analysis revealed that Lhfp was strongly connected to genes involved in osteoblast differentiation. To directly evaluate its role in bone, Lhfp deficient mice (Lhfp-/-) were created using CRISPR/Cas9. Consistent with genetic and network predictions, bone marrow stromal cells (BMSCs) from Lhfp-/- mice displayed increased osteogenic differentiation. Lhfp-/- mice also had elevated BMD due to increased cortical bone mass. Lastly, we identified SNPs in human LHFP that were associated (P = 1.2 x 10−5) with heel BMD. In conclusion, we used GWAS and systems genetics to identify Lhfp as a regulator of osteoblast activity and bone mass.

Networking in Biology: The Hybrid Rat Diversity Panel.

One of the most fruitful resources for systems genetic studies of nonhuman mammals is a panel of inbred strains that exhibits significant genetic diversity between strains but genetic stability (isogenicity) within strains. These characteristics allow for fine mapping of complex phenotypes (QTLs) and provide statistical power to identify loci which contribute nominally to the phenotype. This type of resource also allows the planning and performance of investigations using the same genetic backgrounds over several generations of the test animals. Often, rats are preferred over mice for physiologic and behavioral studies because of their larger size and more distinguishable anatomy (particularly for their central nervous system). The Hybrid Rat Diversity Panel (HRDP) is a panel of inbred rat strains, which combines two recombinant inbred panels (the HXB/BXH, 30 strains; the LEXF/FXLE, 34 strains and 35 more strains of inbred rats which were selected for genetic diversity, based on their fully sequenced genomes and/or thorough genotyping). The genetic diversity and statistical power of this panel for mapping studies rivals or surpasses currently available panels in mouse. The genetic stability of this panel makes it particularly suitable for collection of high-throughput omics data as relevant technology becomes available for engaging in truly integrative systems biology. The PhenoGen website ( http://phenogen.org ) is the repository for the initial transcriptome data, making the raw data, the processed data, and the analysis results, e.g., organ-specific protein coding and noncoding transcripts, isoform analysis, expression quantitative trait loci, and co-expression networks, available to the research public. The data sets and tools being developed will complement current efforts to analyze the human transcriptome and its genetic controls (the Genotype-Tissue Expression Project (GTEx)) and allow for dissection of genetic networks that predispose to particular phenotypes and gene-by-environment interactions that are difficult or even impossible to study in humans. The HRDP is an essential population for exploring truly integrative systems genetics.

High-Resolution Maps of Mouse Reference Populations.

Genetic reference panels are widely used to map complex, quantitative traits in model organisms. We have generated new high-resolution genetic maps of 259 mouse inbred strains from recombinant inbred strain panels (C57BL/6J × DBA/2J, ILS/IbgTejJ × ISS/IbgTejJ, and C57BL/6J × A/J) and chromosome substitution strain panels (C57BL/6J-Chr#<A/J>, C57BL/6J-Chr#<PWD/Ph>, and C57BL/6J-Chr#<MSM/Ms>). We genotyped all samples using the Affymetrix Mouse Diversity Array with an average intermarker spacing of 4.3 kb. The new genetic maps provide increased precision in the localization of recombination breakpoints compared to the previous maps. Although the strains were presumed to be fully inbred, we found residual heterozygosity in 40% of individual mice from five of the six panels. We also identified de novo deletions and duplications, in homozygous or heterozygous state, ranging in size from 21 kb to 8.4 Mb. Almost two-thirds (46 out of 76) of these deletions overlap exons of protein coding genes and may have phenotypic consequences. Twenty-nine putative gene conversions were identified in the chromosome substitution strains. We find that gene conversions are more likely to occur in regions where the homologous chromosomes are more similar. The raw genotyping data and genetic maps of these strain panels are available at http://churchill-lab.jax.org/website/MDA.

Network Analysis Implicates Alpha-Synuclein (Snca) in the Regulation of Ovariectomy-Induced Bone Loss.

The postmenopausal period in women is associated with decreased circulating estrogen levels, which accelerate bone loss and increase the risk of fracture. Here, we gained novel insight into the molecular mechanisms mediating bone loss in ovariectomized (OVX) mice, a model of human menopause, using co-expression network analysis. Specifically, we generated a co-expression network consisting of 53 gene modules using expression profiles from intact and OVX mice from a panel of inbred strains. The expression of four modules was altered by OVX, including module 23 whose expression was decreased by OVX across all strains. Module 23 was enriched for genes involved in the response to oxidative stress, a process known to be involved in OVX-induced bone loss. Additionally, module 23 homologs were co-expressed in human bone marrow. Alpha synuclein (Snca) was one of the most highly connected “hub” genes in module 23. We characterized mice deficient in Snca and observed a 40% reduction in OVX-induced bone loss. Furthermore, protection was associated with the altered expression of specific network modules, including module 23. In summary, the results of this study suggest that Snca regulates bone network homeostasis and ovariectomy-induced bone loss.

Genetic Links between Recombination and Speciation

Genetic Links between Recombination and Speciation

Acute adaptive immune response correlates with late radiation-induced pulmonary fibrosis in mice.

BackgroundThe lung response to radiation exposure can involve an immediate or early reaction to the radiation challenge, including cell death and an initial immune reaction, and can be followed by a tissue injury response, of pneumonitis or fibrosis, to this acute reaction. Herein, we aimed to determine whether markers of the initial immune response, measured within days of radiation exposure, are correlated with the lung tissue injury responses occurring weeks later.MethodsInbred strains of mice known to be susceptible (KK/HIJ, C57BL/6J, 129S1/SvImJ) or resistant (C3H/HeJ, A/J, AKR/J) to radiation-induced pulmonary fibrosis and to vary in time to onset of respiratory distress post thoracic irradiation (from 10–23 weeks) were studied. Mice were untreated (controls) or received 18 Gy whole thorax irradiation and were euthanized at 6 h, 1d or 7 d after radiation treatment. Pulmonary CD4+ lymphocytes, bronchoalveolar cell profile & cytokine level, and serum cytokine levels were assayed.ResultsThoracic irradiation and inbred strain background significantly affected the numbers of CD4+ cells in the lungs and the bronchoalveolar lavage cell differential of exposed mice. At the 7 day timepoint greater numbers of pulmonary Th1 and Th17 lymphocytes and reduced lavage interleukin17 and interferonγ levels were significant predictors of late stage fibrosis. Lavage levels of interleukin-10, measured at the 7 day timepoint, were inversely correlated with fibrosis score (R = −0.80, p = 0.05), while serum levels of interleukin-17 in control mice significantly correlated with post irradiation survival time (R = 0.81, p = 0.04). Lavage macrophage, lymphocyte or neutrophil counts were not significantly correlated with either of fibrosis score or time to respiratory distress in the six mouse strains.ConclusionSpecific cytokine and lymphocyte levels, but not strain dependent lavage cell profiles, were predictive of later radiation-induced lung injury in this panel of inbred strains.Electronic supplementary materialThe online version of this article (doi:10.1186/s13014-015-0359-y) contains supplementary material, which is available to authorized users.

Strain specific differences in the effects of TGF-beta on hepatocyte viability and CYP450 activity in mice

Primary hepatocytes represent a widely used experimental system in the field of toxicology research and drug studies. In this study, we compared hepatocytes from mouse inbred strains C57/BL6 and DBA/2J. These two lines are the parental strains of the BXD recombinant inbred strain panel. We performed a detailed comparison of isolated hepatocytes of the two strains for susceptibility to transforming growth factor (TGF)-beta exposure and induced effects. We detected strain specific differences upon TGF-beta exposure, namely a significant decrease in the viable cell number within five days of culture and an earlier onset of dedifferentiation in DBA/2J as compared to C57/BL6 hepatocytes. In addition, the analysis of the central carbon metabolism of both strains upon TGF-beta exposure and control conditions revealed no significant differences between the hepatocytes of both strains. Therefore, in contrast to the analyzed liver-specific functions, the central carbon metabolism appears to be robust against TGF-beta exposure. The presented results provide detailed insight into strain specific differences and thereby valuable information for the selection of strains for studies such as clinically relevant fibrosis mechanisms, drug metabolism and toxicity as well as studies dealing with the analysis of TGF-beta signaling.

Abstract 595: Genomic, Transcriptomic and MiRNomic Analysis of Triglyceride and Cholesterol Distribution into Lipoprotein Fractions in the Rat

Background: Dyslipidemia is central to the definition of metabolic syndrome (MS), one of most prevalent human diseases worldwide. We preformed genome-wide association study (GWAS) of triacylglycerols (TG) and cholesterol (C) distribution into lipoprotein fractions in the recombinant inbred strain panel PXO (segregating alleles of two MS models, SHR and PD strains, together with those of normolipidemic Brown Norway strain origin), followed by transcriptomic and miRNomic (microRNA profiling) analyses. Methods: We established morphometric and metabolic profile in adult male rats of 14 PXO strains and two progenitor strains (n=183) including glucose tolerance and TG and C concentrations in 20 lipoprotein fractions. GWAS utilizing &gt;20,000 SNPs was performed using MapManager, the significance validated by 2000 permutations per trait. The hepatic transcriptome and miRNome profiles of the most contrasting strains were generated using Affymetrix GeneAtlas system followed by network analysis (Ingenuity Pathway Analysis). Results: We have identified 14 haplotype blocks showing suggestive or significant linkage to studied traits. Except for LDL-TG loci on chromosomes 3 and 12, PXO strains carrying the SHR allele displayed significantly higher values of the lipid linked traits, e.g. LDL-C (21.2±0.4 vs. 12.5±0.4 mg/dl in PXO strains with SHR allele vs. BXH2 allele). C concentrations in large, medium and very small LDL particles were significantly associated to a single gene ( Lrp1b ). Subsequent transcriptomic comparison of phenotypically most contrasting identified series of dysregulated metabolic and signaling pathways including cholesterol biosynthesis and the key upstream regulators such as HNF1 , HNF4 and PPARA . Conclusion: We identified several novel variants associated to TG and C concentrations in lipoprotein fractions together with their transcriptomic and biological network correlates.

Mapping genetic determinants of kidney damage in rat models

During the last two decades, significant progress in our understanding of the development of kidney diseases has been achieved by unravelling the mechanisms underlying rare familial forms of human kidney diseases. Due to the genetic heterogeneity in human populations and the complex multifactorial pathogenesis of the disease phenotypes, the dissection of the genetic basis of common chronic kidney diseases (CKD) remains a difficult task. In this regard, several inbred rat models provide valuable complementary tools to uncover the genetic basis of complex renal disease phenotypes that are related to common forms of CKD. In this review, data obtained in nine experimental rat models, including the Buffalo (BUF), Dahl salt-sensitive (SS), Fawn-hooded hypertensive (FHH), Goto-Kakizaki (GK), Lyon hypertensive (LH), Munich Wistar Frömter (MWF), Sabra hypertension-prone (SBH), spontaneously hypertensive rat (SHR) and stroke-prone spontaneously hypertensive rat (SHRSP) inbred strains, that contributed to the genetic dissection of renal disease phenotypes are presented. In this panel of inbred strains, a large number of quantitative trait loci (QTL) linked to albuminuria/proteinuria and other functional or structural kidney abnormalities could be identified by QTL mapping analysis and follow-up studies including consomic and congenic rat lines. The comprehensive exploitation of the genotype-renal phenotype associations that are inherited in this panel of rat strains is suitable for making a significant contribution to the development of an integrated approach to the systems genetics of common CKD.

Joint genotyping on the fly: identifying variation among a sequenced panel of inbred lines.

High-throughput sequencing is enabling remarkably deep surveys of genomic variation. It is now possible to completely sequence multiple individuals from a single species, yet the identification of variation among them remains an evolving computational challenge. This challenge is compounded for experimental organisms when strains are studied instead of individuals. In response, we present the Joint Genotyper for Inbred Lines (JGIL) as a method for obtaining genotypes and identifying variation among a large panel of inbred strains or lines. JGIL inputs the sequence reads from each line after their alignment to a common reference. Its probabilistic model includes site-specific parameters common to all lines that describe the frequency of nucleotides segregating in the population from which the inbred panel was derived. The distribution of line genotypes is conditional on these parameters and reflects the experimental design. Site-specific error probabilities, also common to all lines, parameterize the distribution of reads conditional on line genotype and realized coverage. Both sets of parameters are estimated per site from the aggregate read data, and posterior probabilities are calculated to decode the genotype of each line. We present an application of JGIL to 162 inbred Drosophila melanogaster lines from the Drosophila Genetic Reference Panel. We explore by simulation the effect of varying coverage, sequencing error, mapping error, and the number of lines. In doing so, we illustrate how JGIL is robust to moderate levels of error. Supported by these analyses, we advocate the importance of modeling the data and the experimental design when possible.

Genetic analysis of complex traits in the emerging Collaborative Cross

The Collaborative Cross (CC) is a mouse recombinant inbred strain panel that is being developed as a resource for mammalian systems genetics. Here we describe an experiment that uses partially inbred CC lines to evaluate the genetic properties and utility of this emerging resource. Genome-wide analysis of the incipient strains reveals high genetic diversity, balanced allele frequencies, and dense, evenly distributed recombination sites-all ideal qualities for a systems genetics resource. We map discrete, complex, and biomolecular traits and contrast two quantitative trait locus (QTL) mapping approaches. Analysis based on inferred haplotypes improves power, reduces false discovery, and provides information to identify and prioritize candidate genes that is unique to multifounder crosses like the CC. The number of expression QTLs discovered here exceeds all previous efforts at eQTL mapping in mice, and we map local eQTL at 1-Mb resolution. We demonstrate that the genetic diversity of the CC, which derives from random mixing of eight founder strains, results in high phenotypic diversity and enhances our ability to map causative loci underlying complex disease-related traits.

The Role of Interleukin-1 in Wound Biology. Part I

Wound healing is a multistep, complex process that involves the coordinated action of multiple cell types. Conflicting results have been obtained when conventional methods have been used to study wound biology. Therefore, we analyzed the wound response in a mouse genetic model. We analyzed inflammatory mediators produced within incisional wounds induced in 16 inbred mouse strains. Computational haplotype-based genetic analysis of inter-strain differences in the level of production of 2 chemokines in wounds was performed. An in vitro experimental analysis system was developed to investigate whether interleukin (IL)-1 could affect chemokine production by 2 different types of cells that are present within wounds. The level of 2 chemokines, keratinocyte-derived chemokine (KC) and macrophage inflammatory protein 1α, exhibited very large (75- and 463-fold, respectively) interstrain differences within wound tissue across this inbred strain panel. Genetic variation within Nalp1, an inflammasome component that regulates IL-1 production, correlated with the interstrain differences in KC and macrophage inhibitory protein 1α production. Consistent with the genetic correlation, IL-1β was shown to stimulate KC production by murine keratinocyte and fibroblast cell lines in vitro. Genetic variation within Nalp1 could contribute to interstrain differences in wound chemokine production by altering the amount of IL-1 produced.

Strain Differences in Stress Responsivity Are Associated with Divergent Amygdala Gene Expression and Glutamate-Mediated Neuronal Excitability

Stress is a major risk factor for numerous neuropsychiatric diseases. However, susceptibility to stress and the qualitative nature of stress effects on behavior differ markedly among individuals. This is partly because of the moderating influence of genetic factors. Inbred mouse strains provide a relatively stable and restricted range of genetic and environmental variability that is valuable for disentangling gene-stress interactions. Here, we screened a panel of inbred strains for anxiety- and depression-related phenotypes at baseline (trait) and after exposure to repeated restraint. Two strains, DBA/2J and C57BL/6J, differed in trait and restraint-induced anxiety-related behavior (dark/light exploration, elevated plus maze). Gene expression analysis of amygdala, medial prefrontal cortex, and hippocampus revealed divergent expression in DBA/2J and C57BL/6J both at baseline and after repeated restraint. Restraint produced strain-dependent expression alterations in various genes including glutamate receptors (e.g., Grin1, Grik1). To elucidate neuronal correlates of these strain differences, we performed ex vivo analysis of glutamate excitatory neurotransmission in amygdala principal neurons. Repeated restraint augmented amygdala excitatory postsynaptic signaling and altered metaplasticity (temporal summation of NMDA receptor currents) in DBA/2J but not C57BL/6J. Furthermore, we found that the C57BL/6J-like changes in anxiety-related behavior after restraint were absent in null mutants lacking the modulatory NMDA receptor subunit Grin2a, but not the AMPA receptor subunit Gria1. Grin2a null mutants exhibited significant ( approximately 30%) loss of dendritic spines on amygdala principal neurons under nonrestraint conditions. Collectively, our data support a model in which genetic variation in glutamatergic neuroplasticity in corticolimbic circuitry underlies phenotypic variation in responsivity to stress.

Murine Warriors or Worriers: The Saga of Comt1, B2 SINE Elements, and the Future of Translational Genetics

Murine Warriors or Worriers: The Saga of Comt1, B2 SINE Elements, and the Future of Translational Genetics

A Common and Unstable Copy Number Variant Is Associated with Differences in Glo1 Expression and Anxiety-Like Behavior

Glyoxalase 1 (Glo1) has been implicated in anxiety-like behavior in mice and in multiple psychiatric diseases in humans. We used mouse Affymetrix exon arrays to detect copy number variants (CNV) among inbred mouse strains and thereby identified a ∼475 kb tandem duplication on chromosome 17 that includes Glo1 (30,174,390–30,651,226 Mb; mouse genome build 36). We developed a PCR-based strategy and used it to detect this duplication in 23 of 71 inbred strains tested, and in various outbred and wild-caught mice. Presence of the duplication is associated with a cis-acting expression QTL for Glo1 (LOD>30) in BXD recombinant inbred strains. However, evidence for an eQTL for Glo1 was not obtained when we analyzed single SNPs or 3-SNP haplotypes in a panel of 27 inbred strains. We conclude that association analysis in the inbred strain panel failed to detect an eQTL because the duplication was present on multiple highly divergent haplotypes. Furthermore, we suggest that non-allelic homologous recombination has led to multiple reversions to the non-duplicated state among inbred strains. We show associations between multiple duplication-containing haplotypes, Glo1 expression and anxiety-like behavior in both inbred strain panels and outbred CD-1 mice. Our findings provide a molecular basis for differential expression of Glo1 and further implicate Glo1 in anxiety-like behavior. More broadly, these results identify problems with commonly employed tests for association in inbred strains when CNVs are present. Finally, these data provide an example of biologically significant phenotypic variability in model organisms that can be attributed to CNVs.

Animal Models of Gene–Nutrient Interactions

Food intake of humans is governed by the food's nutritional value and pleasing taste, but also by other factors such as food cost and availability, cultural imperatives, and social status. The biological determinants of human food intake are not easily parsed from these other factors, making them hard to study against the whirligig aspects of human life in a modern age. The study of animals provides a useful alternative. Humans have a history of studying animal food intake, for agricultural reasons (e.g., pigs and cows), and for personal reasons (e.g., dogs and cats), and these practical concerns have been joined with the appreciation that other models can teach us the principles of behavior, genetics, and nutrition. Thus there is a steady use of the traditional animal models in this type of research, as well as growth in the use of other systems such as worms and flies. Rats and mice occupy a special niche as animal models for two reasons; first, they share with humans a love of the same types of food, and second, they are the target of a number of well-developed genetic tools. The available genetic tools that make mice a popular model include a well-annotated genome (Mouse Build 37), profiles of RNA expression from many tissues, a diverse panel of inbred strains, and the ability to manipulate genes in the whole animal, including removing a gene only in specific tissues (e.g., Cre-lox system). Mice have been harnessed to find genotypes that contribute to sweet-liking, and other studies are underway to understand how genetic variation might at least partially explain other puzzles of human appetites. Animal models provide a way to study the genetic determinants of food selection with experimental rigor and therefore complement human genetics studies.