Genome-wide association for methamphetamine sensitivity in an advanced intercross mouse line

The HNRNPF/H RNA binding proteins and disease

The members of the HNRNPF/H family of heterogeneous nuclear RNA proteins-HNRNPF, HNRNPH1, HNRNPH2, HNRNPH3, and GRSF1, are critical regulators of RNA maturation. Documented functions of these proteins include regulating splicing, particularly alternative splicing, 5' capping and 3' polyadenylation of RNAs, and RNA export. The assignment of these proteins to the HNRNPF/H protein family members relates to differences in the amino acid composition of their RNA recognition motifs, which differ from those of other RNA binding proteins (RBPs). HNRNPF/H proteins typically bind RNA sequences enriched with guanine (G) residues, including sequences that, in the presence of a cation, have the potential to form higher-order G-quadruplex structures. The need to further investigate members of the HNRNPF/H family of RBPs has intensified with the recent descriptions of their involvement in several disease states, including the pediatric tumor Ewing sarcoma and the hematological malignancy mantle cell lymphoma; newly described groups of developmental syndromes; and neuronal-related disorders, including addictive behavior. Here, to foster the study of the HNRNPF/H family of RBPs, we discuss features of the genes encoding these proteins, their structures and functions, and emerging contributions to disease. This article is categorized under: RNA in Disease and Development > RNA in Disease RNA Processing > Splicing Regulation/Alternative Splicing RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications.

Genetics of murine type 2 diabetes and comorbidities

ABSTRAC

Type 2 diabetes (T2D) is a polygenic and multifactorial complex disease, defined as chronic metabolic disorder. It's a major global health concern with an estimated 463 million adults aged 20-79 years with diabetes and projected to increase up to 700 million by 2045. T2D was reported to be one of the four leading causes of non-communicable disease (NCD) deaths in 2012. Environmental factors play a part in the development of polygenic forms of diabetes. Polygenic forms of diabetes often run-in families. Fortunately, T2D, which accounts for 90-95% of the entire four types of diabetes including, Type 1 diabetes (T1D), T2D, monogenic diabetes syndromes (MGDS), and Gestational diabetes mellitus, can be prevented or delayed through nutrition and lifestyle changes as well as through pharmacologic interventions. Typical symptom of the T2D is high blood glucose levels and comprehensive insulin resistance of the body, producing an impaired glucose tolerance. Impaired glucose tolerance of T2D is accompanied by extensive health complications, including cardiovascular diseases (CVD) that vary in morbidity and mortality among populations. The pathogenesis of T2D varies between populations and/or ethnic groupings and is known to be attributed extremely by genetic components and environmental factors. It is evident that genetic background plays a critical role in determining the host response toward certain environmental conditions, whether or not of developing T2D (susceptibility versus resistant). T2D is considered as a silent disease that can progress for years before its diagnosis. Once T2D is diagnosed, many metabolic malfunctions are observed whether as side effects or as independent comorbidity. Mouse models have been proven to be a powerful tool for mapping genetic factors that underline the susceptibility to T2D development as well its comorbidities. Here, we have conducted a comprehensive search throughout the published data covering the time span from early 1990s till the time of writing this review, for already reported quantitative trait locus (QTL) associated with murine T2D and comorbidities in different mouse models, which contain different genetic backgrounds. Our search has resulted in finding 54 QTLs associated with T2D in addition to 72 QTLs associated with comorbidities associated with the disease. We summarized the genomic locations of these mapped QTLs in graphical formats, so as to show the overlapping positions between of these mapped QTLs, which may suggest that some of these QTLs could be underlined by sharing gene/s. Finally, we reviewed and addressed published reports that show the success of translation of the identified mouse QTLs/genes associated with the disease in humans.

Systems genetic analysis of binge-like eating in a C57BL/6J x DBA/2J-F2 cross

Binge eating is a heritable trait associated with eating disorders and refers to the rapid consumption of a large quantity of energy-dense food that is, associated with loss of control and negative affect. Binge eating disorder is the most common eating disorder in the United States; however, the genetic basis is unknown. We previously identified robust mouse inbred strain differences between C57BL/6J and DBA/2J in binge-like eating of sweetened palatable food in an intermittent access, conditioned place preference paradigm. To map the genetic basis of changes in body weight and binge-like eating (BLE) and to identify candidate genes, we conducted quantitative trait locus (QTL) analysis in 128 C57BL/6J x DBA/2J-F2 mice combined with PheQTL and trait covariance analysis in GeneNetwork2 using legacy BXD-RI trait datasets. We identified a QTL on Chromosome 18 influencing changes in body weight across days in females (log of the odds [LOD] = 6.3; 1.5-LOD: 3-12 cM) that contains the candidate gene Zeb1. We also identified a sex-combined QTL influencing initial palatable food intake on Chromosome 5 (LOD = 5.8; 1.5-LOD: 21-28 cM) that contains the candidate gene Lcorl and a second QTL influencing escalated palatable food intake on Chromosome 6 in males (LOD = 5.4; 1.5-LOD: 50-59 cM) that contains the candidate genes Adipor2 and Plxnd1. Finally, we identified a suggestive QTL in females for slope of BLE on distal Chromosome 18 (LOD = 4.1; p = 0.055; 1.5-LOD: 23-35 cM). Future studies will use BXD-RI strains to fine map loci and support candidate gene nomination for gene editing.

Genetic Dissection of Growth Traits in a Unique Chicken Advanced Intercross Line

The advanced intercross line (AIL) that is created by successive generations of pseudo-random mating after the F₂ generation is a valuable resource, especially in agricultural livestock and poultry species, because it improves the precision of quantitative trait loci (QTL) mapping compared with traditional association populations by introducing more recombination events. The growth traits of broilers have significant economic value in the chicken industry, and many QTLs affecting growth traits have been identified, especially on chromosomes 1, 4, and 27, albeit with large confidence intervals that potentially contain dozens of genes. To promote a better understanding of the underlying genetic architecture of growth trait differences, specifically body weight and bone development, in this study, we report a nine-generation AIL derived from two divergent outbred lines: High Quality chicken Line A (HQLA) and Huiyang Bearded (HB) chicken. We evaluate the genetic architecture of the F₀, F₂, F₈, and F₉ generations of AIL and demonstrate that the population of the F₉ generation sufficiently randomized the founder genomes and has the characteristics of rapid linkage disequilibrium decay, limited allele frequency decline, and abundant nucleotide diversity. This AIL yielded a much narrower QTL than the F₂ generations, especially the QTL on chromosome 27, which was reduced to 120 Kb. An ancestral haplotype association analysis showed that most of the dominant haplotypes are inherited from HQLA but with fluctuation of the effects between them. We highlight the important role of four candidate genes (PHOSPHO1, IGF2BP1, ZNF652, and GIP) in bone growth. We also retrieved a missing QTL from AIL on chromosome 4 by identifying the founder selection signatures, which are explained by the loss of association power that results from rare alleles. Our study provides a reasonable resource for detecting quantitative trait genes and tracking ancestor history and will facilitate our understanding of the genetic mechanisms underlying chicken bone growth.

Differential genetic risk for methamphetamine intake confers differential sensitivity to the temperature-altering effects of other addictive drugs

Mice selectively bred for high methamphetamine (MA) drinking (MAHDR), compared with mice bred for low MA drinking (MALDR), exhibit greater sensitivity to MA reward and insensitivity to aversive and hypothermic effects of MA. Previous work identified the trace amine-associated receptor 1 gene (Taar1) as a quantitative trait gene for MA intake that also impacts thermal response to MA. All MAHDR mice are homozygous for the mutant Taar1 m1J allele, whereas all MALDR mice possess at least one copy of the reference Taar1 + allele. To determine if their differential sensitivity to MA-induced hypothermia extends to drugs of similar and different classes, we examined sensitivity to the hypothermic effect of the stimulant cocaine, the amphetamine-like substance 3,4-methylenedioxymethamphetamine (MDMA), and the opioid morphine in these lines. The lines did not differ in thermal response to cocaine, only MALDR mice exhibited a hypothermic response to MDMA, and MAHDR mice were more sensitive to the hypothermic effect of morphine than MALDR mice. We speculated that the μ-opioid receptor gene (Oprm1) impacts morphine response, and genotyped the mice tested for morphine-induced hypothermia. We report genetic linkage between Taar1 and Oprm1; MAHDR mice more often inherit the Oprm1 D2 allele and MALDR mice more often inherit the Oprm1 B6 allele. Data from a family of recombinant inbred mouse strains support the influence of Oprm1 genotype, but not Taar1 genotype, on thermal response to morphine. These results nominate Oprm1 as a genetic risk factor for morphine-induced hypothermia, and provide additional evidence for a connection between drug preference and drug thermal response.

5' UTR variants in the quantitative trait gene Hnrnph1 support reduced 5' UTR usage and hnRNP H protein as a molecular mechanism underlying reduced methamphetamine sensitivity

We previously identified a 210 kb region on chromosome 11 (50.37-50.58 Mb, mm10) containing two protein-coding genes (Hnrnph1, Rufy1) that was necessary for reduced methamphetamine-induced locomotor activity in C57BL/6J congenic mice harboring DBA/2J polymorphisms. Gene editing of a small deletion in the first coding exon supported Hnrnph1 as a quantitative trait gene. We have since shown that Hnrnph1 mutants also exhibit reduced methamphetamine-induced reward, reinforcement, and dopamine release. However, the quantitative trait variants (QTVs) that modulate Hnrnph1 function at the molecular level are not known. Nine single nucleotide polymorphisms and seven indels distinguish C57BL/6J from DBA/2J within Hnrnph1, including four variants within the 5' untranslated region (UTR). Here, we show that a 114 kb introgressed region containing Hnrnph1 and Rufy1 was sufficient to cause a decrease in MA-induced locomotor activity. Gene-level transcriptome analysis of striatal tissue from 114 kb congenics vs Hnrnph1 mutants identified a nearly perfect correlation of fold-change in expression for those differentially expressed genes that were common to both mouse lines, indicating functionally similar effects on the transcriptome and behavior. Exon-level analysis (including noncoding exons) revealed decreased 5' UTR usage of Hnrnph1 and immunoblot analysis identified a corresponding decrease in hnRNP H protein in 114 kb congenic mice. Molecular cloning of the Hnrnph1 5' UTR containing all four variants (but none of them individually) upstream of a reporter induced a decrease in reporter signal in both HEK293 and N2a cells, thus, identifying a set of QTVs underlying molecular regulation of Hnrnph1.

Genome-Wide Association Study in Two Cohorts from a Multi-generational Mouse Advanced Intercross Line Highlights the Difficulty of Replication Due to Study-Specific Heterogeneity

There has been extensive discussion of the "Replication Crisis" in many fields, including genome-wide association studies (GWAS). We explored replication in a mouse model using an advanced intercross line (AIL), which is a multigenerational intercross between two inbred strains. We re-genotyped a previously published cohort of LG/J x SM/J AIL mice (F34; n = 428) using a denser marker set and genotyped a new cohort of AIL mice (F39-43; n = 600) for the first time. We identified 36 novel genome-wide significant loci in the F34 and 25 novel loci in the F39-43 cohort. The subset of traits that were measured in both cohorts (locomotor activity, body weight, and coat color) showed high genetic correlations, although the SNP heritabilities were slightly lower in the F39-43 cohort. For this subset of traits, we attempted to replicate loci identified in either F34 or F39-43 in the other cohort. Coat color was robustly replicated; locomotor activity and body weight were only partially replicated, which was inconsistent with our power simulations. We used a random effects model to show that the partial replications could not be explained by Winner's Curse but could be explained by study-specific heterogeneity. Despite this heterogeneity, we performed a mega-analysis by combining F34 and F39-43 cohorts (n = 1,028), which identified four novel loci associated with locomotor activity and body weight. These results illustrate that even with the high degree of genetic and environmental control possible in our experimental system, replication was hindered by study-specific heterogeneity, which has broad implications for ongoing concerns about reproducibility.

Using Heterogeneous Stocks for Fine-Mapping Genetically Complex Traits

In this chapter we will review both the rationale and experimental design for using Heterogeneous Stock (HS) populations for fine-mapping of complex traits in mice and rats. We define an HS as an outbred population derived from an intercross between two or more inbred strains. HS have been used to perform genome-wide association studies (GWAS) for multiple behavioral, physiological, and gene expression traits. GWAS using HS require four key steps, which we review: selection of an appropriate HS population, phenotyping, genotyping, and statistical analysis. We provide advice on the selection of an HS, comment on key issues related to phenotyping, discuss genotyping methods relevant to these populations, and describe statistical genetic analyses that are applicable to genetic analyses that use HS.

Genetic differences in the behavioral organization of binge eating, conditioned food reward, and compulsive-like eating in C57BL/6J and DBA/2J strains

Binge eating (BE) is a heritable symptom of eating disorders associated with anxiety, depression, malnutrition, and obesity. Genetic analysis of BE could facilitate therapeutic discovery. We used an intermittent, limited access BE paradigm involving sweetened palatable food (PF) to examine genetic differences in BE, conditioned food reward, and compulsive-like eating between C57BL/6J (B6J) and DBA/2J (D2J) inbred mouse strains. D2J mice showed a robust escalation in intake and conditioned place preference for the PF-paired side. D2J mice also showed a unique style of compulsive-like eating in the light/dark conflict test where they rapidly hoarded and consumed PF in the preferred unlit environment. BE and compulsive-like eating exhibited narrow-sense heritability estimates between 56 and 73%. To gain insight into the genetic basis, we phenotyped and genotyped a small cohort of 133 B6J × D2J-F₂ mice at the peak location of three quantitative trait loci (QTL) previously identified in F₂ mice for sweet taste (chromosome 4: 156 Mb), bitter taste (chromosome 6: 133 Mb) and behavioral sensitivity to drugs of abuse (chromosome 11: 50 Mb). The D2J allele on chromosome 6 was associated with greater PF intake on training days and greater compulsive-like PF intake, but only in males, suggesting that decreased bitter taste may increase BE in males. The D2J allele on chromosome 11 was associated with an increase in final PF intake and slope of escalation across days. Future studies employing larger crosses and genetic reference panels comprising B6J and D2J alleles will identify causal genes and neurobiological mechanisms.

Heterogeneous Stock Populations for Analysis of Complex Traits

Heterogeneous Stock (HS) populations allow for fine-resolution genetic mapping of a variety of complex traits. HS mice and rats were created from breeding together eight inbred strains, followed by maintaining the colony in a manner that minimizes inbreeding. After 50 or more generations of breeding, the resulting animals' chromosomes represent a genetic mosaic of the founders' haplotypes, with the average distance between recombination events in the centiMorgan range. This allows for genetic mapping to only a few Mb, a much smaller region than what can be identified using traditional F2 intercross or backcross mapping strategies. HS animals have been used to fine-map a variety of complex traits including anxiety and fear behaviors, diabetes, asthma, and heart disease, among others. Once a quantitative trait locus (QTL) has been identified, founder sequence and expression analysis can be used to identify underlying causal genes. In the following review, we provide an overview of how HS rats and mice have been used to identify genetic loci, and in some cases the causal genes, underlying complex traits. We discuss the creation and breeding strategies for both HS rats and mice. We then discuss the statistical analyses used to identify genetic loci, as well as strategies to identify causal genes underlying these loci. We end the chapter by discussing limitations faced when using HS populations, including several statistical challenges that have not been fully resolved.

Identifying genes for neurobehavioural traits in rodents: progress and pitfalls

Identifying genes and pathways that contribute to differences in neurobehavioural traits is a key goal in psychiatric research. Despite considerable success in identifying quantitative trait loci (QTLs) associated with behaviour in laboratory rodents, pinpointing the causal variants and genes is more challenging. For a long time, the main obstacle was the size of QTLs, which could encompass tens if not hundreds of genes. However, recent studies have exploited mouse and rat resources that allow mapping of phenotypes to narrow intervals, encompassing only a few genes. Here, we review these studies, showcase the rodent resources they have used and highlight the insights into neurobehavioural traits provided to date. We discuss what we see as the biggest challenge in the field - translating QTLs into biological knowledge by experimentally validating and functionally characterizing candidate genes - and propose that the CRISPR/Cas genome-editing system holds the key to overcoming this obstacle. Finally, we challenge traditional views on inbred versus outbred resources in the light of recent resource and technology developments.

Dissection of Host Susceptibility to Bacterial Infections and Its Toxins

Infection is one of the leading causes of human mortality and morbidity. Exposure to microbial agents is obviously required. However, also non-microbial environmental and host factors play a key role in the onset, development and outcome of infectious disease, resulting in large of clinical variability between individuals in a population infected with the same microbe. Controlled and standardized investigations of the genetics of susceptibility to infectious disease are almost impossible to perform in humans whereas mouse models allow application of powerful genomic techniques to identify and validate causative genes underlying human diseases with complex etiologies. Most of current animal models used in complex traits diseases genetic mapping have limited genetic diversity. This limitation impedes the ability to create incorporated network using genetic interactions, epigenetics, environmental factors, microbiota, and other phenotypes. A novel mouse genetic reference population for high-resolution mapping and subsequently identifying genes underlying the QTL, namely the Collaborative Cross (CC) mouse genetic reference population (GRP) was recently developed. In this chapter, we discuss a variety of approaches using CC mice for mapping genes underlying quantitative trait loci (QTL) to dissect the host response to polygenic traits, including infectious disease caused by bacterial agents and its toxins.

Genome-wide association study of behavioral, physiological and gene expression traits in outbred CFW mice

Although mice are the most widely used mammalian model organism, genetic studies have suffered from limited mapping resolution due to extensive linkage disequilibrium (LD) that is characteristic of crosses among inbred strains. Carworth Farms White (CFW) mice are a commercially available outbred mouse population that exhibit rapid LD decay in comparison to other available mouse populations. We performed a genome-wide association study (GWAS) of behavioral, physiological and gene expression phenotypes using 1,200 male CFW mice. We used genotyping by sequencing (GBS) to obtain genotypes at 92,734 SNPs. We also measured gene expression using RNA sequencing in three brain regions. Our study identified numerous behavioral, physiological and expression quantitative trait loci (QTLs). We integrated the behavioral QTL and eQTL results to implicate specific genes, including Azi2 in sensitivity to methamphetamine and Zmynd11 in anxiety-like behavior. The combination of CFW mice, GBS and RNA sequencing constitutes a powerful approach to GWAS in mice.

The Genetic Basis of Baculum Size and Shape Variation in Mice

The rapid divergence of male genitalia is a preeminent evolutionary pattern. This rapid divergence is especially striking in the baculum, a bone that occurs in the penis of many mammalian species. Closely related species often display diverse baculum morphology where no other morphological differences can be discerned. While this fundamental pattern of evolution has been appreciated at the level of gross morphology, nearly nothing is known about the genetic basis of size and shape divergence. Quantifying the genetic basis of baculum size and shape variation has been difficult because these structures generally lack obvious landmarks, so comparing them in three dimensions is not straightforward. Here, we develop a novel morphometric approach to quantify size and shape variation from three-dimensional micro-CT scans taken from 369 bacula, representing 75 distinct strains of the BXD family of mice. We identify two quantitative trait loci (QTL) that explain ∼50% of the variance in baculum size, and a third QTL that explains more than 20% of the variance in shape. Together, our study demonstrates that baculum morphology may diverge relatively easily, with mutations at a few loci of large effect that independently modulate size and shape. Based on a combination of bioinformatic investigations and new data on RNA expression, we prioritized these QTL to 16 candidate genes, which have hypothesized roles in bone morphogenesis and may enable future genetic manipulation of baculum morphology.

Hnrnph1 Is A Quantitative Trait Gene for Methamphetamine Sensitivity

Psychostimulant addiction is a heritable substance use disorder; however its genetic basis is almost entirely unknown. Quantitative trait locus (QTL) mapping in mice offers a complementary approach to human genome-wide association studies and can facilitate environment control, statistical power, novel gene discovery, and neurobiological mechanisms. We used interval-specific congenic mouse lines carrying various segments of chromosome 11 from the DBA/2J strain on an isogenic C57BL/6J background to positionally clone a 206 kb QTL (50,185,512–50,391,845 bp) that was causally associated with a reduction in the locomotor stimulant response to methamphetamine (2 mg/kg, i.p.; DBA/2J < C57BL/6J)—a non-contingent, drug-induced behavior that is associated with stimulation of the dopaminergic reward circuitry. This chromosomal region contained only two protein coding genes—heterogeneous nuclear ribonucleoprotein, H1 (Hnrnph1) and RUN and FYVE domain-containing 1 (Rufy1). Transcriptome analysis via mRNA sequencing in the striatum implicated a neurobiological mechanism involving a reduction in mesolimbic innervation and striatal neurotransmission. For instance, Nr4a2 (nuclear receptor subfamily 4, group A, member 2), a transcription factor crucial for midbrain dopaminergic neuron development, exhibited a 2.1-fold decrease in expression (DBA/2J < C57BL/6J; p 4.2 x 10−¹⁵). Transcription activator-like effector nucleases (TALENs)-mediated introduction of frameshift deletions in the first coding exon of Hnrnph1, but not Rufy1, recapitulated the reduced methamphetamine behavioral response, thus identifying Hnrnph1 as a quantitative trait gene for methamphetamine sensitivity. These results define a novel contribution of Hnrnph1 to neurobehavioral dysfunction associated with dopaminergic neurotransmission. These findings could have implications for understanding the genetic basis of methamphetamine addiction in humans and the development of novel therapeutics for prevention and treatment of substance abuse and possibly other psychiatric disorders.

Both genetic and environmental factors can powerfully modulate susceptibility to substance use disorders. Quantitative trait locus (QTL) mapping is an unbiased discovery-based approach that is used to identify novel genetic factors and provide new mechanistic insight into phenotypic variation associated with disease. In this study, we focused on the genetic basis of variation in sensitivity to the acute locomotor stimulant response to methamphetamine which is a behavioral phenotype in rodents that is associated with stimulated dopamine release and activation of the brain reward circuitry involved in addiction. Using brute force monitoring of recombination events associated with changes in behavior, we fortuitously narrowed the genotype-phenotype association down to just two genes that we subsequently targeted using a contemporary genome editing approach. The gene that we validated–Hnrnph1 –is an RNA binding protein that did not have any previously known function in psychostimulant behavior or psychostimulant addiction. Our behavioral data combined with our gene expression results provide a compelling rationale for a new line of investigation regarding Hnrnph1 and its role in neural development and plasticity associated with the addictions and perhaps other dopamine-dependent psychiatric disorders.

High-resolution genetic mapping of complex traits from a combined analysis of F2 and advanced intercross mice

Genetic influences on anxiety disorders are well documented; however, the specific genes underlying these disorders remain largely unknown. To identify quantitative trait loci (QTL) for conditioned fear and open field behavior, we used an F2 intercross (n = 490) and a 34th-generation advanced intercross line (AIL) (n = 687) from the LG/J and SM/J inbred mouse strains. The F2 provided strong support for several QTL, but within wide chromosomal regions. The AIL yielded much narrower QTL, but the results were less statistically significant, despite the larger number of mice. Simultaneous analysis of the F2 and AIL provided strong support for QTL and within much narrower regions. We used a linear mixed-model approach, implemented in the program QTLRel, to correct for possible confounding due to familial relatedness. Because we recorded the full pedigree, we were able to empirically compare two ways of accounting for relatedness: using the pedigree to estimate kinship coefficients and using genetic marker estimates of "realized relatedness." QTL mapping using the marker-based estimates yielded more support for QTL, but only when we excluded the chromosome being scanned from the marker-based relatedness estimates. We used a forward model selection procedure to assess evidence for multiple QTL on the same chromosome. Overall, we identified 12 significant loci for behaviors in the open field and 12 significant loci for conditioned fear behaviors. Our approach implements multiple advances to integrated analysis of F2 and AILs that provide both power and precision, while maintaining the advantages of using only two inbred strains to map QTL.

Association of novelty-related behaviors and intravenous cocaine self-administration in Diversity Outbred mice

RATIONALE

The preference for and reaction to novelty are strongly associated with addiction to cocaine and other drugs. However, the genetic variants and molecular mechanisms underlying these phenomena remain largely unknown. Although the relationship between novelty- and addiction-related traits has been observed in rats, studies in mice have failed to demonstrate this association. New, genetically diverse, high-precision mouse populations including Diversity Outbred (DO) mice provide an opportunity to assess an expanded range of behavioral variation enabling detection of associations of novelty- and addiction-related traits in mice.

METHODS

To examine the relationship between novelty- and addiction-related traits, male (n = 51) and female (n = 47) DO mice were tested on open field exploration, hole board exploration, and novelty preference followed by intravenous cocaine self-administration (IVSA; ten 2-h sessions of fixed ratio 1 and one 6-h session of progressive ratio).

RESULTS

We observed high variation of cocaine IVSA in DO mice with 43 % reaching and 57 % not reaching conventional acquisition criteria. As a group, mice that did not reach these criteria still demonstrated significant lever discrimination. Mice experiencing catheter occlusion or other technical issues (n = 17) were excluded from the analysis. Novelty-related behaviors were positively associated with cocaine IVSA. Multivariate analysis of associations among novelty- and addiction-related traits revealed a large degree of shared variance (45 %).

CONCLUSIONS

Covariation among cocaine IVSA and novelty-related phenotypes in DO mice indicates that this relationship is amenable to genetic dissection. The high genetic precision and phenotypic diversity in the DO may facilitate discovery of previously undetectable mechanisms underlying predisposition to develop addiction disorders.

Discovery and refinement of muscle weight QTLs in B6 × D2 advanced intercross mice

The genes underlying variation in skeletal muscle mass are poorly understood. Although many quantitative trait loci (QTLs) have been mapped in crosses of mouse strains, the limited resolution inherent in these conventional studies has made it difficult to reliably pinpoint the causal genetic variants. The accumulated recombination events in an advanced intercross line (AIL), in which mice from two inbred strains are mated at random for several generations, can improve mapping resolution. We demonstrate these advancements in mapping QTLs for hindlimb muscle weights in an AIL (n = 832) of the C57BL/6J (B6) and DBA/2J (D2) strains, generations F8-F13. We mapped muscle weight QTLs using the high-density MegaMUGA SNP panel. The QTLs highlight the shared genetic architecture of four hindlimb muscles and suggest that the genetic contributions to muscle variation are substantially different in males and females, at least in the B6D2 lineage. Out of the 15 muscle weight QTLs identified in the AIL, nine overlapped the genomic regions discovered in an earlier B6D2 F2 intercross. Mapping resolution, however, was substantially improved in our study to a median QTL interval of 12.5 Mb. Subsequent sequence analysis of the QTL regions revealed 20 genes with nonsense or potentially damaging missense mutations. Further refinement of the muscle weight QTLs using additional functional information, such as gene expression differences between alleles, will be important for discerning the causal genes.

Fine-mapping QTLs in advanced intercross lines and other outbred populations

Quantitative genetic studies in model organisms, particularly in mice, have been extremely successful in identifying chromosomal regions that are associated with a wide variety of behavioral and other traits. However, it is now widely understood that identification of the underlying genes will be far more challenging. In the last few years, a variety of populations have been utilized in an effort to more finely map these chromosomal regions with the goal of identifying specific genes. The common property of these newer populations is that linkage disequilibrium spans relatively short distances, which permits fine-scale mapping resolution. This review focuses on advanced intercross lines (AILs) which are the simplest such population. As originally proposed in 1995 by Darvasi and Soller, an AIL is the product of intercrossing two inbred strains beyond the F2 generation. Unlike recombinant inbred strains, AILs are maintained as outbred populations; brother-sister matings are specifically avoided. Each generation of intercrossing beyond the F2 further degrades linkage disequilibrium between adjacent makers, which allows for fine-scale mapping of quantitative trait loci (QTLs). Advances in genotyping technology and techniques for the statistical analysis of AILs have permitted rapid advances in the application of AILs. We review some of the analytical issues and available software, including QTLRel, EMMA, EMMAX, GEMMA, TASSEL, GRAMMAR, WOMBAT, Mendel, and others.

QTL mapping in outbred populations: successes and challenges

Quantitative trait locus (QTL) mapping in animal populations has been a successful strategy for identifying genomic regions that play a role in complex diseases and traits. When conducted in an F2 intercross or backcross population, the resulting QTL is frequently large, often encompassing 30 Mb or more and containing hundreds of genes. To narrow the locus and identify candidate genes, additional strategies are needed. Congenic strains have proven useful but work less well when there are multiple tightly linked loci, frequently resulting in loss of phenotype. As an alternative, we discuss the use of highly recombinant outbred models for directly fine-mapping QTL to only a few megabases. We discuss the use of several currently available models such as the advanced intercross (AI), heterogeneous stocks (HS), the diversity outbred (DO), and commercially available outbred stocks (CO). Once a QTL has been fine-mapped, founder sequence and expression QTL mapping can be used to identify candidate genes. In this regard, the large number of alleles found in outbred stocks can be leveraged to identify causative genes and variants. We end this review by discussing some important statistical considerations when analyzing outbred populations. Fine-resolution mapping in outbred models, coupled with full genome sequence, has already led to the identification of several underlying causative genes for many complex traits and diseases. These resources will likely lead to additional successes in the coming years.

Out of the bottleneck: the Diversity Outcross and Collaborative Cross mouse populations in behavioral genetics research

The historical origins of classical laboratory mouse strains have led to a relatively limited range of genetic and phenotypic variation, particularly for the study of behavior. Many recent efforts have resulted in improved diversity and precision of mouse genetic resources for behavioral research, including the Collaborative Cross and Diversity Outcross population. These two populations, derived from an eight way cross of common and wild-derived strains, have high precision and allelic diversity. Behavioral variation in the population is expanded, both qualitatively and quantitatively. Variation that had once been canalized among the various inbred lines has been made amenable to genetic dissection. The genetic attributes of these complementary populations, along with advances in genetic and genomic technologies, makes a systems genetic analyses of behavior more readily tractable, enabling discovery of a greater range of neurobiological phenomena underlying behavioral variation.

Practical considerations regarding the use of genotype and pedigree data to model relatedness in the context of genome-wide association studies

Genome-wide association studies of complex traits often are complicated by relatedness among individuals. Ignoring or inappropriately accounting for relatedness often results in inflated type I error rates. Either genotype or pedigree data can be used to estimate relatedness for use in mixed-models when undertaking quantitative trait locus mapping. We performed simulations to investigate methods for controlling type I error and optimizing power considering both full and partial pedigrees and, similarly, both sparse and dense marker coverage; we also examined real data sets. (1) When marker density was low, estimating relatedness by genotype data alone failed to control the type I error rate; (2) this was resolved by combining both genotype and pedigree data. (3) When sufficiently dense marker data were used to estimate relatedness, type I error was well controlled and power increased; however, (4) this was only true when the relatedness was estimated using genotype data that excluded genotypes on the chromosome currently being scanned for a quantitative trait locus.

High-precision genetic mapping of behavioral traits in the diversity outbred mouse population

Historically our ability to identify genetic variants underlying complex behavioral traits in mice has been limited by low mapping resolution of conventional mouse crosses. The newly developed Diversity Outbred (DO) population promises to deliver improved resolution that will circumvent costly fine-mapping studies. The DO is derived from the same founder strains as the Collaborative Cross (CC), including three wild-derived strains. Thus the DO provides more allelic diversity and greater potential for discovery compared to crosses involving standard mouse strains. We have characterized 283 male and female DO mice using open-field, light–dark box, tail-suspension and visual-cliff avoidance tests to generate 38 behavioral measures. We identified several quantitative trait loci (QTL) for these traits with support intervals ranging from 1 to 3 Mb in size. These intervals contain relatively few genes (ranging from 5 to 96). For a majority of QTL, using the founder allelic effects together with whole genome sequence data, we could further narrow the positional candidates. Several QTL replicate previously published loci. Novel loci were also identified for anxiety- and activity-related traits. Half of the QTLs are associated with wild-derived alleles, confirming the value to behavioral genetics of added genetic diversity in the DO. In the presence of wild-alleles we sometimes observe behaviors that are qualitatively different from the expected response. Our results demonstrate that high-precision mapping of behavioral traits can be achieved with moderate numbers of DO animals, representing a significant advance in our ability to leverage the mouse as a tool for behavioral genetics

Methamphetamine-induced conditioned place preference in LG/J and SM/J mouse strains and an F45/F46 advanced intercross line

The conditioned place preference (CPP) test is frequently used to evaluate the rewarding properties of drugs of abuse in mice. Despite its widespread use in transgenic and knockout experiments, there are few forward genetic studies using CPP to identify novel genes contributing to drug reward. In this study, we tested LG/J and SM/J inbred strains and the parents/offspring of 10 families of an F₄₅/F₄₆ advanced intercross line (AIL) for methamphetamine-induced CPP (MA-CPP) once per week over 2 weeks. Both LG/J and SM/J mice exhibited significant MA-CPP that was not significantly different between the two strains. Furthermore, LG/J mice showed significantly less acute MA-induced locomotor activity as well as locomotor sensitization following subsequent MA injections. AIL mice (N = 105) segregating LG/J and SM/J alleles also demonstrated significant MA-CPP that was equal in magnitude between the first and second week of training. Importantly, MA-CPP in AIL mice did not correlate with drug-free or MA-induced locomotor activity, indicating that MA-CPP was not confounded by test session activity and implying that MA-CPP is genetically distinct from acute psychomotor sensitivity. We estimated the heritability of MA-CPP and locomotor phenotypes using midparent-offspring regression and maximum likelihood estimates derived from the kinship coefficients of the AIL pedigree. Heritability estimates of MA-CPP were low (0–0.21) and variable (SE = 0–0.33) which reflected our poor power to estimate heritability using only 10 midparent-offspring observations. In sum, we established a short-term protocol for MA-CPP in AIL mice that could reveal LG/J and SM/J alleles important for MA reward. The use of highly recombinant genetic populations like AIL should facilitate the identification of these genes and may have implications for understanding psychostimulant abuse in humans.

Csnk1e is a genetic regulator of sensitivity to psychostimulants and opioids

Csnk1e, the gene encoding casein kinase 1-epsilon, has been implicated in sensitivity to amphetamines. Additionally, a polymorphism in CSNK1E was associated with heroin addiction, suggesting that this gene may also affect opioid sensitivity. In this study, we first conducted genome-wide quantitative trait locus (QTL) mapping of methamphetamine (MA)-induced locomotor activity in C57BL/6J (B6) × DBA/2J (D2)-F(2) mice and a more highly recombinant F(8) advanced intercross line. We identified a QTL on chromosome 15 that contained Csnk1e (63-86 Mb; Csnk1e=79.25 Mb). We replicated this result and further narrowed the locus using B6.D2(Csnk1e) and D2.B6(Csnk1e) reciprocal congenic lines (78-86.8 and 78.7-81.6 Mb, respectively). This locus also affected sensitivity to the μ-opioid receptor agonist fentanyl. Next, we directly tested the hypothesis that Csnk1e is a genetic regulator of sensitivity to psychostimulants and opioids. Mice harboring a null allele of Csnk1e showed an increase in locomotor activity following MA administration. Consistent with this result, coadministration of a selective pharmacological inhibitor of Csnk1e (PF-4800567) increased the locomotor stimulant response to both MA and fentanyl. These results show that a narrow genetic locus that contains Csnk1e is associated with differences in sensitivity to MA and fentanyl. Furthermore, gene knockout and selective pharmacological inhibition of Csnk1e define its role as a negative regulator of sensitivity to psychostimulants and opioids.

Genome-wide association for fear conditioning in an advanced intercross mouse line

Fear conditioning (FC) may provide a useful model for some components of post-traumatic stress disorder (PTSD). We used a C57BL/6J × DBA/2J F(2) intercross (n = 620) and a C57BL/6J × DBA/2J F(8) advanced intercross line (n = 567) to fine-map quantitative trait loci (QTL) associated with FC. We conducted an integrated genome-wide association analysis in QTLRel and identified five highly significant QTL affecting freezing to context as well as four highly significant QTL associated with freezing to cue. The average percent decrease in QTL width between the F(2) and the integrated analysis was 59.2%. Next, we exploited bioinformatic sequence and expression data to identify candidate genes based on the existence of non-synonymous coding polymorphisms and/or expression QTLs. We identified numerous candidate genes that have been previously implicated in either fear learning in animal models (Bcl2, Btg2, Dbi, Gabr1b, Lypd1, Pam and Rgs14) or PTSD in humans (Gabra2, Oprm1 and Trkb); other identified genes may represent novel findings. The integration of F(2) and AIL data maintains the advantages of studying FC in model organisms while significantly improving resolution over previous approaches.