Background Eating disorders are lethal and heritable; however, the underlying genetic factors are unknown. Binge eating is a highly heritable trait associated with eating disorders that is comorbid with mood and substance use disorders. Therefore, understanding its genetic basis will inform therapeutic development that could improve several comorbid neuropsychiatric conditions. Methods We assessed binge eating in closely related C57BL/6 mouse substrains and in an F2 cross to identify quantitative trait loci (QTL) associated with binge eating. We used gene targeting to validated candidate genetic factors. Finally, we used transcriptome analysis of the striatum via mRNA sequencing (RNA-seq) to identify the premorbid transcriptome and the binge-induced transcriptome to inform molecular mechanisms mediating binge eating susceptibility and establishment. Results C57BL/6NJ but not C57BL/6J mice showed rapid and robust escalation in palatable food consumption. We mapped a single genome-wide significant QTL on chromosome 11 (LOD=7.4) to a missense mutation in cytoplasmic FMR1-interacting protein 2 (Cyfip2). We validated Cyfip2 as a major genetic factor underlying binge eating in heterozygous knockout mice on a C57BL/6N background that showed reduced binge eating toward a wild-type C57BL/6J-like level. Transcriptome analysis of premorbid genetic risk identified the enrichment terms “morphine addiction” and “retrograde endocannabinoid signaling” whereas binge eating resulted in the downregulation of a gene set enriched for decreased myelination, oligodendrocyte differentiation, and expression. Conclusions We identified Cyfip2 as a major significant genetic factor underlying binge eating and provide a behavioral paradigm for future genome-wide association studies in populations with increased genetic complexity.

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−15). 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.

Drug liking vs. drug disliking is a subjective motivational measure in humans that assesses the addiction liability of drugs. Variation in this trait is hypothesized to influence vulnerability vs. resilience toward substance abuse disorders and likely contains a genetic component. In rodents and humans, conditioned place preference (CPP)/aversion (CPA) is a Pavlovian conditioning paradigm whereby a learned preference for the drug-paired environment is used to infer drug liking whereas a learned avoidance or aversion is used to infer drug disliking. C57BL/6 inbred mouse substrains are nearly genetically identical, yet demonstrate robust differences in addiction-relevant behaviors, including locomotor sensitization to cocaine and consumption of ethanol. Here, we tested the hypothesis that B6 substrains would demonstrate differences in the rewarding properties of the mu opioid receptor agonist oxycodone (5 mg/kg, i.p.) and the aversive properties of the opioid receptor antagonist naloxone (4 mg/kg, i.p.). Both substrains showed similar degrees of oxycodone-induced CPP; however, there was a three-fold enhancement of naloxone-induced CPA in agonist-naïve C57BL/6J relative to C57Bl/6NJ mice. Exploratory factor analysis of CPP and CPA identified unique factors that explain variance in behavioral expression of reward vs. aversion. “Conditioned Opioid-Like Behavior” was a reward-based factor whereby drug-free locomotor variables resembling opioid treatment co-varied with the degree of CPP. “Avoidance and Freezing” was an aversion-based factor, whereby the increase in the number of freezing bouts co-varied with the degree of aversion. These results provide new insight into the behavioral architecture of the motivational properties of opioids. Future studies will use quantitative trait locus mapping in B6 substrains to identify novel genetic factors that contribute to the marked strain difference in NAL-CPA.

Sensitivity to different pain modalities has a genetic basis that remains largely unknown. Employing closely related inbred mouse substrains can facilitate gene mapping of nociceptive behaviors in preclinical pain models. We previously reported enhanced sensitivity to acute thermal nociception in C57BL/6J ( B6J ) versus C57BL/6N ( B6N ) substrains. Here, we expanded on nociceptive phenotypes and observed an increase in formalin-induced inflammatory nociceptive behaviors and paw diameter in B6J versus B6N mice (Charles River Laboratories). No strain differences were observed in mechanical or thermal hypersensitivity or in edema following the Complete Freund’s Adjuvant model of inflammatory pain, indicating specificity in the inflammatory nociceptive stimulus. In the chronic constrictive nerve injury, a model of neuropathic pain, no strain differences were observed in baseline mechanical threshold or in mechanical hypersensitivity up to one month post-chronic constrictive nerve injury. We replicated the enhanced thermal nociception in the 52.5°C hot plate test in B6J versus B6N mice from The Jackson Laboratory. Using a B6J × B6N-F2 cross (N = 164), we mapped a major quantitative trait locus underlying hot plate sensitivity to chromosome 7 that peaked at 26 Mb (log of the odds [LOD] = 3.81, p < 0.01; 8.74 Mb-36.50 Mb) that was more pronounced in males. Genes containing expression quantitative trait loci associated with the peak nociceptive marker that are implicated in pain and inflammation include Ryr1 , Cyp2a5 , Pou2f2 , Clip3 , Sirt2 , Actn4 , and Ltbp4 (false discovery rate < 0.05). Future studies involving positional cloning and gene editing will determine the quantitative trait gene(s) and potential pleiotropy of this locus across pain modalities.

Psychostimulant (methamphetamine, cocaine) use disorders have a genetic component that remains mostly unknown. We conducted genome‐wide quantitative trait locus (QTL) analysis of methamphetamine stimulant sensitivity. To facilitate gene identification, we employed a Reduced Complexity Cross between closely related C57BL/6 mouse substrains and examined maximum speed and distance traveled over 30 min following methamphetamine (2 mg/kg, i.p.). For maximum methamphetamine‐induced speed following the second and third administration, we identified a single genome‐wide significant QTL on chromosome 11 that peaked near the Cyfip2 locus (LOD = 3.5, 4.2; peak = 21 cM [36 Mb]). For methamphetamine‐induced distance traveled following the first and second administration, we identified a genome‐wide significant QTL on chromosome 5 that peaked near a functional intronic indel in Gabra2 coding for the alpha‐2 subunit of the GABA‐A receptor (LOD = 3.6–5.2; peak = 34–35 cM [66–67 Mb]). Striatal cis‐expression QTL mapping corroborated Gabra2 as a functional candidate gene underlying methamphetamine‐induced distance traveled. CRISPR/Cas9‐mediated correction of the mutant intronic deletion on the C57BL/6J background to the wild‐type C57BL/6NJ allele was sufficient to reduce methamphetamine‐induced locomotor activity toward the wild‐type C57BL/6NJ‐like level, thus validating the quantitative trait variant (QTV). These studies show the power and efficiency of Reduced Complexity Crosses in identifying causal variants underlying complex traits. Functionally restoring Gabra2 expression decreased methamphetamine stimulant sensitivity and supports preclinical and human genetic studies implicating the GABA‐A receptor in psychostimulant addiction‐relevant traits. Importantly, our findings have major implications for studying psychostimulants in the C57BL/6J strain—the gold standard strain in biomedical research.

Binge eating (BE) is a heritable trait associated with eating disorders and involves episodes of rapid, large amounts of food consumption. We previously identified cytoplasmic FMR1-interacting protein 2 ( Cyfip2 ) as a genetic factor underlying compulsive-like BE in mice. CYFIP2 is a homolog of CYFIP1 which is one of four paternally-deleted genes in patients with Type I Prader-Willi Syndrome (PWS), a neurodevelopmental disorder whereby 70% of cases involve paternal 15q11-q13 deletion. PWS symptoms include hyperphagia, obesity (if untreated), cognitive deficits, and obsessive-compulsive behaviors. We tested whether Cyfip1 haploinsufficiency (+/−) would enhance compulsive-like behavior and palatable food (PF) intake in a parental origin- and sex-dependent manner on two Cyfip2 genetic backgrounds, including the BE-prone C57BL/6N ( Cyfip2 N/N ) background and the BE-resistant C57BL/6J ( Cyfip2 J/J ) background. Cyfip1 +/− mice showed increased compulsive-like behavior on both backgrounds and increased PF intake on the Cyfip2 N/N background. In contrast, maternal Cyfip1 haploinsufficiency on the BE-resistant Cyfip2 J/J background induced a robust escalation in PF intake in wild-type Cyfip1 J/J males while having no effect in Cyfip1 J/- males. Notably, induction of behavioral phenotypes in wild-type males following maternal Fmr1 +/− has previously been reported. In the hypothalamus, there was a paternally-enhanced reduction in CYFIP1 protein whereas in the nucleus accumbens, there was a maternally-enhanced reduction in CYFIP1 protein. Nochange in FMR1 protein (FMRP) was observed in Cyfip1 +/− mice, regardless of parental origin. To summarize, Cyfip1 haploinsufficiency increased compulsive-like behavior and induced genetic background-dependent, sex-dependent, and parent-of-origin-dependent effects on PF consumption and CYFIP1 expression that could have relevance for neurodevelopmental and neuropsychiatric disorders.

Using chromosome substitution strains (CSS), we previously identified a large quantitative trait locus (QTL) for conditioned fear (CF) on mouse chromosome 10. Here, we used an F2 cross between CSS-10 and C57BL/6J (B6) to localize that QTL to distal chromosome 10. That QTL accounted for all of the difference between CSS-10 and B6. We then produced congenic strains to fine-map that interval. We identified two congenic strains that captured some or all of the QTL. The larger congenic strain (Line 1; 122.387121 – 129.068 Mb; build 37) appeared to account for all of the difference between CSS-10 and B6. The smaller congenic strain (Line 2; 127.277–129.068 Mb) was intermediate between CSS-10 and B6. We used haplotype mapping followed by qPCR to identify one gene that was differentially expressed in both lines relative to B6 (Rnf41) and one that was differentially expressed between only Line 1 and B6 (Shmt2). These cis-eQTLs may cause the behavioral QTLs; however, further studies are required to validate these candidate genes. More generally, our observation that a large QTL mapped using CSS and F2 crosses can be dissected into multiple smaller QTLs demonstrates a weaknesses of two-stage approaches that seek to use coarse mapping to identify large regions followed by fine-mapping. Indeed, additional dissection of these congenic strains might result in further subdivision of these QTL regions. Despite these limitations we have successfully fine mapped two QTLs to small regions and identified putative candidate genes, demonstrating that the congenic approach can be effective for fine mapping QTLs.

Psychostimulant (methamphetamine, cocaine) use disorders have a genetic component that remains mostly unknown. Here, we conducted genome-wide quantitative trait locus (QTL) analysis of methamphetamine stimulant sensitivity. To facilitate gene identification, we employed a Reduced Complexity Cross between closely related C57BL/6 mouse substrains and examined maximum speed and distance traveled over 30 min following methamphetamine (2 mg/kg, i.p.). For maximum methamphetamine-induced speed following the second and third administration, we identified a single genome-wide significant QTL on chromosome 11 that peaked near the Cyfip2 locus [LOD = 3.5, 4.2; peak = 21 cM (36 Mb)]. For methamphetamine-induced distance traveled, we identified a single genome-wide significant QTL on chromosome 5 that peaked near a functional intronic indel in Gabra2 that codes for the alpha-2 subunit of the GABA-A receptor [LOD = 5.2; peak = 35 cM (67 Mb)]. Striatal cis-expression QTL mapping corroborated Gabra2 as a functional candidate gene underlying methamphetamine-induced distance traveled. CRISPR/Cas9-mediated correction of the mutant intronic deletion on the C57BL/6J background to the wild-type C57BL/6NJ allele was sufficient to reduce methamphetamine-induced locomotor activity toward the wild-type C57BL/6NJ-like level, thus validating the quantitative trait variant (QTV). These studies demonstrate the power and efficiency of Reduced Complexity Crosses in identifying causal genes and variants underlying complex traits. Functionally restoring Gabra2 expression decreased methamphetamine stimulant sensitivity and supports preclinical and human genetic studies implicating the GABA-A receptor in psychostimulant addiction-relevant traits. Importantly, our findings have major implications for investigators studying psychostimulants in the C57BL/6J strain - the gold standard strain in biomedical research.