High-throughput behavioral phenotyping in the expanded panel of BXD recombinant inbred strains

Identification of ethanol analgesia quantitative trait loci and candidate genes in BXD recombinant inbred mouse lines

Alcohol consumption produces acute analgesic effects, and people experiencing pain conditions may drink alcohol to alleviate discomfort. However, tolerance to the analgesic properties of alcohol could prompt escalating consumption and dependence. Both nociception and alcohol-induced analgesia are under significant genetic control. Understanding the genetic architecture of these processes could inform better treatment options for people with pain conditions. This study aims to identify quantitative trait loci (QTL) driving variation in ethanol-induced analgesia across BXD recombinant inbred mouse lines. Male and female mice from 62 BXD strains received ethanol or saline oral gavage for five days and were tested for hot plate (HP) latency at baseline, Day 1, and Day 5. QTL mapping of HP phenotypes identified a significant provisional QTL on chromosome 17 for Day 1 HP latency in mice receiving ethanol. An additional highly suggestive QTL was present on chromosome 9 for the difference in pre- and post-ethanol thermal nociception. Candidate genes within QTL support intervals were provisionally identified using HP phenotypic correlations to transcriptomic database, expression QTL analysis, and other bioinformatics inquiries. The combined behavioral and bioinformatic analyses yielded strong ethanol analgesia candidate genes, specifically Myo6 . Thus, the results of this genetic study of ethanol-induced analgesia in BXD mouse strains may contribute significantly to our understanding of the molecular basis for individual variation in the analgesic response to acute ethanol.

Complementation testing identifies genes mediating effects at quantitative trait loci underlying fear-related behavior

Knowing the genes involved in quantitative traits provides an entry point to understanding the biological bases of behavior, but there are very few examples where the pathway from genetic locus to behavioral change is known. To explore the role of specific genes in fear behavior, we mapped three fear-related traits, tested fourteen genes at six quantitative trait loci (QTLs) by quantitative complementation, and identified six genes. Four genes, Lamp, Ptprd, Nptx2, and Sh3gl, have known roles in synapse function; the fifth, Psip1, was not previously implicated in behavior; and the sixth is a long non-coding RNA, 4933413L06Rik, of unknown function. Variation in transcriptome and epigenetic modalities occurred preferentially in excitatory neurons, suggesting that genetic variation is more permissible in excitatory than inhibitory neuronal circuits. Our results relieve a bottleneck in using genetic mapping of QTLs to uncover biology underlying behavior and prompt a reconsideration of expected relationships between genetic and functional variation.

Lipophilic analogues of D-cysteine prevent and reverse physical dependence to fentanyl in male rats

We examined whether co-injections of the cell-permeant D-cysteine analogues, D-cysteine ethyl ester (D-CYSee) and D-cysteine ethyl amide (D-CYSea), prevent acquisition of physical dependence induced by twice-daily injections of fentanyl, and reverse acquired dependence to these injections in freely-moving male Sprague Dawley rats. Injection of the opioid receptor antagonist, naloxone HCl (NLX, 1.5 mg/kg, IV), elicited a series of withdrawal phenomena that included cardiorespiratory and behavioral responses, and falls in body weight and body temperature, in rats that received 5 or 10 injections of fentanyl (125 μg/kg, IV), and the same number of vehicle co-injections. Regarding the development of physical dependence, the NLX-precipitated withdrawal phenomena were markedly reduced in fentanyl-injected rats that had received co-injections of D-CYSee (250 μmol/kg, IV) or D-CYSea (100 μmol/kg, IV), but not D-cysteine (250 μmol/kg, IV). Regarding reversal of established dependence to fentanyl, the NLX-precipitated withdrawal phenomena in rats that had received 10 injections of fentanyl (125 μg/kg, IV) was markedly reduced in rats that received co-injections of D-CYSee (250 μmol/kg, IV) or D-CYSea (100 μmol/kg, IV), but not D-cysteine (250 μmol/kg, IV), starting with injection 6 of fentanyl. This study provides evidence that co-injections of D-CYSee and D-CYSea prevent the acquisition of physical dependence, and reverse acquired dependence to fentanyl in male rats. The lack of effect of D-cysteine suggests that the enhanced cell-penetrability of D-CYSee and D-CYSea into cells, particularly within the brain, is key to their ability to interact with intracellular signaling events involved in acquisition to physical dependence to fentanyl.

Discrete prefrontal neuronal circuits determine repeated stress-induced behavioral phenotypes in male mice

Chronic stress is a major risk factor for psychiatric disorders, including depression. Although depression is a highly heterogeneous syndrome, it remains unclear how chronic stress drives individual differences in behavioral responses. In this study, we developed a subtyping-based approach wherein stressed male mice were divided into four subtypes based on their behavioral patterns of social interaction deficits and anhedonia, the core symptoms of psychiatric disorders. We identified three prefrontal cortical neuronal projections that regulate repeated stress-induced behavioral phenotypes. Among them, the medial prefrontal cortex (mPFC)→anterior paraventricular thalamus (aPVT) pathway determines the specific behavioral subtype that exhibits both social deficits and anhedonia. Additionally, we identified the circuit-level molecular mechanism underlying this subtype: KDM5C-mediated epigenetic repression of Shisa2 transcription in aPVT projectors in the mPFC led to social deficits and anhedonia. Thus, we provide a set of biological aspects at the cellular, molecular, and epigenetic levels that determine distinctive stress-induced behavioral phenotypes.

Differential genetic expression within reward-specific ensembles in mice

Maladaptive reward seeking is a hallmark of cocaine use disorder. To develop therapeutic targets, it is critical to understand the neurobiological changes specific to cocaine-seeking without altering the seeking of natural rewards, e.g., sucrose. The prefrontal cortex (PFC) and the nucleus accumbens core (NAcore) are known regions associated with cocaine- and sucrose-seeking ensembles, i.e., a sparse population of co-activated neurons. Within ensembles, transcriptomic alterations in the PFC and NAcore underlie the learning and persistence of cocaine- and sucrose-seeking behavior. However, transcriptomes exclusively driving cocaine seeking independent from sucrose seeking have not yet been defined using a within-subject approach. Using Ai14:cFos-TRAP2 transgenic mice in a dual cocaine and sucrose self-administration model, we fluorescently sorted (FACS) and characterized (RNAseq) the transcriptomes defining cocaine- and sucrose-seeking ensembles. We found reward- and region-specific transcriptomic changes that will help develop clinically relevant genetic approaches to decrease cocaine-seeking behavior without altering non-drug reward-based positive reinforcement.

Systems genetics analysis reveals the common genetic basis for pain sensitivity and cognitive function

BACKGROUND

There is growing evidence of a strong correlation between pain sensitivity and cognitive function under both physiological and pathological conditions. However, the detailed mechanisms remain largely unknown. In the current study, we sought to explore candidate genes and common molecular mechanisms underlying pain sensitivity and cognitive function with a transcriptome-wide association study using recombinant inbred mice from the BXD family.

METHODS

The pain sensitivity determined by Hargreaves' paw withdrawal test and cognition-related phenotypes were systematically analyzed in 60 strains of BXD mice and correlated with hippocampus transcriptomes, followed by quantitative trait locus (QTL) mapping and systems genetics analysis.

RESULTS

The pain sensitivity showed significant variability across the BXD strains and co-varies with cognitive traits. Pain sensitivity correlated hippocampual genes showed a significant involvement in cognition-related pathways, including glutamatergic synapse, and PI3K-Akt signaling pathway. Moreover, QTL mapping identified a genomic region on chromosome 4, potentially regulating the variation of pain sensitivity. Integrative analysis of expression QTL mapping, correlation analysis, and Bayesian network modeling identified Ring finger protein 20 (Rnf20) as the best candidate. Further pathway analysis indicated that Rnf20 may regulate the expression of pain sensitivity and cognitive function through the PI3K-Akt signaling pathway, particularly through interactions with genes Ppp2r2b, Ppp2r5c, Col9a3, Met, Rps6, Tnc, and Kras.

CONCLUSIONS

Our study demonstrated that pain sensitivity is associated with genetic background and Rnf20-mediated PI3K-Akt signaling may involve in the regulation of pain sensitivity and cognitive functions.

Complementation testing identifies causal genes at quantitative trait loci underlying fear related behavior

Knowing the genes involved in quantitative traits provides a critical entry point to understanding the biological bases of behavior, but there are very few examples where the pathway from genetic locus to behavioral change is known. Here we address a key step towards that goal by deploying a test that directly queries whether a gene mediates the effect of a quantitative trait locus (QTL). To explore the role of specific genes in fear behavior, we mapped three fear-related traits, tested fourteen genes at six QTLs, and identified six genes. Four genes, Lsamp, Ptprd, Nptx2 and Sh3gl, have known roles in synapse function; the fifth gene, Psip1, is a transcriptional co-activator not previously implicated in behavior; the sixth is a long non-coding RNA 4933413L06Rik with no known function. Single nucleus transcriptomic and epigenetic analyses implicated excitatory neurons as likely mediating the genetic effects. Surprisingly, variation in transcriptome and epigenetic modalities between inbred strains occurred preferentially in excitatory neurons, suggesting that genetic variation is more permissible in excitatory than inhibitory neuronal circuits. Our results open a bottleneck in using genetic mapping of QTLs to find novel biology underlying behavior and prompt a reconsideration of expected relationships between genetic and functional variation.

Discovery and validation of genes driving drug-intake and related behavioral traits in mice

Substance use disorders are heritable disorders characterized by compulsive drug use, the biological mechanisms for which remain largely unknown. Genetic correlations reveal that predisposing drug-naïve phenotypes, including anxiety, depression, novelty preference and sensation seeking, are predictive of drug-use phenotypes, thereby implicating shared genetic mechanisms. High-throughput behavioral screening in knockout (KO) mice allows efficient discovery of the function of genes. We used this strategy in two rounds of candidate prioritization in which we identified 33 drug-use candidate genes based upon predisposing drug-naïve phenotypes and ultimately validated the perturbation of 22 genes as causal drivers of substance intake. We selected 19/221 KO strains (8.5%) that had a difference from control on at least one drug-naïve predictive behavioral phenotype and determined that 15/19 (~80%) affected the consumption or preference for alcohol, methamphetamine or both. No mutant exhibited a difference in nicotine consumption or preference which was possibly confounded with saccharin. In the second round of prioritization, we employed a multivariate approach to identify outliers and performed validation using methamphetamine two-bottle choice and ethanol drinking-in-the-dark protocols. We identified 15/401 KO strains (3.7%, which included one gene from the first cohort) that differed most from controls for the predisposing phenotypes. 8 of 15 gene deletions (53%) affected intake or preference for alcohol, methamphetamine or both. Using multivariate and bioinformatic analyses, we observed multiple relations between predisposing behaviors and drug intake, revealing many distinct biobehavioral processes underlying these relationships. The set of mouse models identified in this study can be used to characterize these addiction-related processes further.

L-cysteine ethyl ester prevents and reverses acquired physical dependence on morphine in male Sprague Dawley rats

The molecular mechanisms underlying the acquisition of addiction/dependence on morphine may result from the ability of the opioid to diminish the transport of L-cysteine into neurons via inhibition of excitatory amino acid transporter 3 (EAA3). The objective of this study was to determine whether the co-administration of the cell-penetrant L-thiol ester, L-cysteine ethyl ester (L-CYSee), would reduce physical dependence on morphine in male Sprague Dawley rats. Injection of the opioid-receptor antagonist, naloxone HCl (NLX; 1.5 mg/kg, IP), elicited pronounced withdrawal phenomena in rats which received a subcutaneous depot of morphine (150 mg/kg) for 36 h and were receiving a continuous infusion of saline (20 μL/h, IV) via osmotic minipumps for the same 36 h period. The withdrawal phenomena included wet-dog shakes, jumping, rearing, fore-paw licking, 360° circling, writhing, apneas, cardiovascular (pressor and tachycardia) responses, hypothermia, and body weight loss. NLX elicited substantially reduced withdrawal syndrome in rats that received an infusion of L-CYSee (20.8 μmol/kg/h, IV) for 36 h. NLX precipitated a marked withdrawal syndrome in rats that had received subcutaneous depots of morphine (150 mg/kg) for 48 h) and a co-infusion of vehicle. However, the NLX-precipitated withdrawal signs were markedly reduced in morphine (150 mg/kg for 48 h)-treated rats that began receiving an infusion of L-CYSee (20.8 μmol/kg/h, IV) at 36 h. In similar studies to those described previously, neither L-cysteine nor L-serine ethyl ester (both at 20.8 μmol/kg/h, IV) mimicked the effects of L-CYSee. This study demonstrates that 1) L-CYSee attenuates the development of physical dependence on morphine in male rats and 2) prior administration of L-CYSee reverses morphine dependence, most likely by intracellular actions within the brain. The lack of the effect of L-serine ethyl ester (oxygen atom instead of sulfur atom) strongly implicates thiol biochemistry in the efficacy of L-CYSee. Accordingly, L-CYSee and analogs may be a novel class of therapeutics that ameliorate the development of physical dependence on opioids in humans.

Developmental coordination disorder: What can we learn from RI mice using motor learning tasks and QTL analysis

Developmental Coordination Disorder (DCD) is a neurodevelopmental disorder of unknown etiology that affects one in 20 children. There is an indication that DCD has an underlying genetic component due to its high heritability. Therefore, we explored the use of a recombinant inbred family of mice known as the BXD panel to understand the genetic basis of complex traits (i.e., motor learning) through identification of quantitative trait loci (QTLs). The overall aim of this study was to utilize the QTL approach to evaluate the genome-to-phenome correlation in BXD strains of mice in order to better understand the human presentation of DCD. Results of this current study confirm differences in motor learning in selected BXD strains and strains with altered cerebellar volume. Five strains - BXD15, BXD27, BXD28, BXD75, and BXD86 - exhibited the most DCD-like phenotype when compared with other BXD strains of interest. Results indicate that BXD15 and BXD75 struggled primarily with gross motor skills, BXD28 primarily had difficulties with fine motor skills, and BXD27 and BXD86 strains struggled with both fine and gross motor skills. The functional roles of genes within significant QTLs were assessed in relation to DCD-like behavior. Only Rab3a (Ras-related protein Rab-3A) emerged as a high likelihood candidate gene for the horizontal ladder rung task. This gene is associated with brain and skeletal muscle development, but lacked nonsynonymous polymorphisms. This study along with Gill et al. (same issue) is the first studies to specifically examine the genetic linkage of DCD using BXD strains of mice.

Modeling Brain Gene Expression in Alcohol Use Disorder with Genetic Animal Models

Animal genetic models have and will continue to provide important new information about the behavioral and physiological adaptations associated with alcohol use disorder (AUD). This chapter focuses on two models, ethanol preference and drinking in the dark (DID), their usefulness in interrogating brain gene expression data and the relevance of the data obtained to interpret AUD-related GWAS and TWAS studies. Both the animal and human data point to the importance for AUD of changes in synaptic transmission (particularly glutamate and GABA transmission), of changes in the extracellular matrix (specifically including collagens, cadherins and protocadherins) and of changes in neuroimmune processes. The implementation of new technologies (e.g., cell type-specific gene expression) is expected to further enhance the value of genetic animal models in understanding AUD.

Adaptor protein complex 2 in the orbitofrontal cortex predicts alcohol use disorder

Alcohol use disorder (AUD) is a life-threatening disease characterized by compulsive drinking, cognitive deficits, and social impairment that continue despite negative consequences. The inability of individuals with AUD to regulate drinking may involve functional deficits in cortical areas that normally balance actions that have aspects of both reward and risk. Among these, the orbitofrontal cortex (OFC) is critically involved in goal-directed behavior and is thought to maintain a representation of reward value that guides decision making. In the present study, we analyzed post-mortem OFC brain samples collected from age- and sex-matched control subjects and those with AUD using proteomics, bioinformatics, machine learning, and reverse genetics approaches. Of the 4,500+ total unique proteins identified in the proteomics screen, there were 47 proteins that differed significantly by sex that were enriched in processes regulating extracellular matrix and axonal structure. Gene ontology enrichment analysis revealed that proteins differentially expressed in AUD cases were involved in synaptic and mitochondrial function, as well as transmembrane transporter activity. Alcohol-sensitive OFC proteins also mapped to abnormal social behaviors and social interactions. Machine learning analysis of the post-mortem OFC proteome revealed dysregulation of presynaptic (e.g., AP2A1) and mitochondrial proteins that predicted the occurrence and severity of AUD. Using a reverse genetics approach to validate a target protein, we found that prefrontal Ap2a1 expression significantly correlated with voluntary alcohol drinking in male and female genetically diverse mouse strains. Moreover, recombinant inbred strains that inherited the C57BL/6J allele at the Ap2a1 interval consumed higher amounts of alcohol than those that inherited the DBA/2J allele. Together, these findings highlight the impact of excessive alcohol consumption on the human OFC proteome and identify important cross-species cortical mechanisms and proteins that control drinking in individuals with AUD.

The ameliorative effects of choline on ethanol-induced cell death in the neural tube of susceptible BXD strains of mice

Introduction

Fetal alcohol spectrum disorders (FASD) are the leading preventable cause of intellectual disability, providing the impetus for evaluating various potential treatments to ameliorate ethanol's teratogenic effects, particularly in the nervous system. One treatment is the dietary supplement choline which has been shown to mitigate at least some of ethanol's teratogenic effects. The present study was designed to investigate the effects of genetics on choline's efficacy in ameliorating cell death in the developing neural tube. Previously, we examined BXD recombinant inbred mice, and their parental C57BL/6 J (B6) and DBA/2 J strains, and identified strains that were sensitive to ethanol's teratogenic actions. Thus, we used these strains to identify response to choline treatment.

Materials and methods

Timed pregnant mice from 4 strains (B6, BXD51, BXD73, BXD2) were given either ethanol or isocaloric maltose-dextrin (5.8 g/kg in two administrations separated by 2 h) with choline at one of 3 doses: 0, 100 or 250 mg/kg. Subjects were exposed via intragastric gavage on embryonic day 9 and embryos were collected 7 h after the initial ethanol administrations. Cell death was analyzed using TUNEL staining in the developing forebrain and brainstem.

Results

Choline ameliorated the ethanol-induced cell death across all 4 strains without causing enhanced cell death in control mice. Choline was effective in both the developing telencephalon and in the brainstem. Both doses diminished cell death, with some differences across strains and brain regions, although the 100 mg/kg dose was most consistent in mitigating ethanol-related cell death. Comparisons across strains showed that there was an effect of strain, particularly in the forebrain at the higher dose.

Discussion

These results show that choline is effective in ameliorating ethanol-induced cell death at this early stage of nervous system development. However, there were some strain differences in its efficacy, especially at the high dose, providing further evidence of the importance of genetics in influencing the ability of choline to protect against prenatal alcohol exposure.

Adaptor protein complex 2 in the orbitofrontal cortex predicts alcohol use disorder

Alcohol use disorder (AUD) is a life-threatening disease characterized by compulsive drinking, cognitive deficits, and social impairment that continue despite negative consequences. The inability of individuals with AUD to regulate drinking may involve functional deficits in cortical areas that normally balance actions that have aspects of both reward and risk. Among these, the orbitofrontal cortex (OFC) is critically involved in goal-directed behavior and is thought to maintain a representation of reward value that guides decision making. In the present study, we analyzed post-mortem OFC brain samples collected from age- and sex-matched control subjects and those with AUD using proteomics, bioinformatics, machine learning, and reverse genetics approaches. Of the 4,500+ total unique proteins identified in the proteomics screen, there were 47 proteins that differed significantly by sex that were enriched in processes regulating extracellular matrix and axonal structure. Gene ontology enrichment analysis revealed that proteins differentially expressed in AUD cases were involved in synaptic and mitochondrial function, as well as transmembrane transporter activity. Alcohol-sensitive OFC proteins also mapped to abnormal social behaviors and social interactions. Machine learning analysis of the post-mortem OFC proteome revealed dysregulation of presynaptic (e.g., AP2A1) and mitochondrial proteins that predicted the occurrence and severity of AUD. Using a reverse genetics approach to validate a target protein, we found that prefrontal Ap2a1 expression significantly correlated with voluntary alcohol drinking in male and female genetically diverse mouse strains. Moreover, recombinant inbred strains that inherited the C57BL/6J allele at the Ap2a1 interval consumed higher amounts of alcohol than those that inherited the DBA/2J allele. Together, these findings highlight the impact of excessive alcohol consumption on the human OFC proteome and identify important cross-species cortical mechanisms and proteins that control drinking in individuals with AUD.

Effects of genetics and sex on adolescent behaviors following neonatal ethanol exposure in BXD recombinant inbred strains

Introduction

Fetal alcohol spectrum disorders (FASD) are the leading preventable neurodevelopmental disorders and two hallmark symptoms of FASD are abnormal behavior, and cognitive and learning deficits. The severity of alcohol's teratogenic effects on the developing brain is influenced by genetics and sex. We previously identified recombinant inbred BXD mouse strains that show differential vulnerability to ethanol-induced cell death in the developing hippocampus, a brain region important in learning and memory. The present study aimed to test the hypothesis that strains with increased vulnerability to ethanol-induced cell death in the hippocampus have concomitant deficits in multiple hippocampal-related behaviors during adolescence.

Methods

The current study evaluated the effects of developmental ethanol exposure on adolescent behavior in two BXD strains that show high cell death (BXD48a, BXD100), two that show low cell death (BXD60, BXD71), and the two parental strains (C57BL/6 J (B6), DBA/2 J (D2)). On postnatal day 7, male and female neonatal pups were treated with ethanol (5.0 g/kg) or saline given in two equal doses 2 h apart. Adolescent behavior was assessed across multiple behavioral paradigms including the elevated plus maze, open field, Y-maze, and T-maze.

Results

Our results demonstrate that the effects of developmental ethanol exposure on adolescent behavioral responses are highly dependent on strain. The low cell death strains, BXD60 and BXD71, showed minimal effect of ethanol exposure on all behavioral measures but did present sex differences. The parental -B6 and D2-strains and high cell death strains, BXD48a and BXD100, showed ethanol-induced effects on activity-related or anxiety-like behaviors. Interestingly, the high cell death strains were the only strains that showed a significant effect of postnatal ethanol exposure on hippocampal-dependent spatial learning and memory behaviors.

Discussion

Overall, we identified effects of ethanol exposure, strain, and/or sex on multiple behavioral measures. Interestingly, the strains that showed the most effects of postnatal ethanol exposure on adolescent behavior were the BXD strains that show high ethanol-induced cell death in the neonatal hippocampus, consistent with our hypothesis. Additionally, we found evidence for interactions among strain and sex, demonstrating that these factors have a complex effect on alcohol responses and that both are important considerations.

DISCOVERY AND VALIDATION OF GENES DRIVING DRUG-INTAKE AND RELATED BEHAVIORAL TRAITS IN MICE

Substance use disorders (SUDs) are heritable disorders characterized by compulsive drug use, but the biological mechanisms driving addiction remain largely unknown. Genetic correlations reveal that predisposing drug-naïve phenotypes, including anxiety, depression, novelty preference, and sensation seeking, are predictive of drug-use phenotypes, implicating shared genetic mechanisms. Because of this relationship, high-throughput behavioral screening of predictive phenotypes in knockout (KO) mice allows efficient discovery of genes likely to be involved in drug use. We used this strategy in two rounds of screening in which we identified 33 drug-use candidate genes and ultimately validated the perturbation of 22 of these genes as causal drivers of substance intake. In our initial round of screening, we employed the two-bottle-choice paradigms to assess alcohol, methamphetamine, and nicotine intake. We identified 19 KO strains that were extreme responders on at least one predictive phenotype. Thirteen of the 19 gene deletions (68%) significantly affected alcohol use three methamphetamine use, and two both. In the second round of screening, we employed a multivariate approach to identify outliers and performed validation using methamphetamine two-bottle choice and ethanol drinking-in-the-dark protocols. We identified 15 KO strains that were extreme responders across the predisposing drug-naïve phenotypes. Eight of the 15 gene deletions (53%) significantly affected intake or preference for three alcohol, eight methamphetamine or three both (3). We observed multiple relations between predisposing behaviors and drug intake, revealing many distinct biobehavioral processes underlying these relationships. The set of mouse models identified in this study can be used to characterize these addiction-related processes further.

Behavioral and biochemical effects of alcohol withdrawal in female C3H/HeNRj and C57BL/6JRj mice

Background

Alcohol use disorder (AUD) is a major problem of our society and is often characterized and worsened by relapse. Prolonged alcohol exposure leads to numerous biochemical alterations that, upon cessation of alcohol intake, cause an array of immediate and lasting withdrawal symptoms. Acute withdrawal and neuroinflammation can be harmful in themselves, and lasting withdrawal symptoms contribute to relapse. Here, we conducted an initial feasibility study assessing several behavioral and neurochemical factors in female C3H/HeNRj (C3H) and C57BL/6JRj (B6) mice to determine which strain showed the clearest alcohol withdrawal symptoms during long-term abstinence and neurochemical alterations following re-exposure.

Methods

Female C3H and B6 mice (n = 12 per group/strain) were intermittently exposed to alcohol-containing or control liquid diets for 3 weeks. Acute and prolonged withdrawal symptoms were assessed over a period of 3 weeks using a battery of behavioral test, comprised of alcohol self-administration, anhedonia, hyperalgesia, anxiety-like and depressive-like disturbances. Brain inflammation was measured by multiplex cytokine assay. Monoamine levels in the hippocampus and striatum, as well as exploratory analyses of cations levels in the cerebellum, were assessed by High-Performance Liquid Chromatography (HPLC).

Results

Both C3H and B6 alcohol-exposed mice displayed decreased saccharin intake or preference and higher stress levels assessed by ultrasonic vocalizations (USVs) recordings. B6 but not C3H alcohol-exposed mice also exhibited a slower decline of alcohol oral self-administration (OSA), hyperalgesia, elevated brain TNF-α and elevated serotonin turnover.

Conclusion

Our findings highlight the suitability of the B6 strain to study the behavioral and neurochemical alterations caused by alcohol withdrawal and the potential efficacy of experimental treatments, not only in early detoxification, but also in prolonged abstinence. The feasibility of these assays is important because long-lasting withdrawal symptoms are often the main cause of relapse in alcohol-dependent patients.

Genetic Mapping of Behavioral Traits Using the Collaborative Cross Resource

The complicated interactions between genetic background, environment and lifestyle factors make it difficult to study the genetic basis of complex phenotypes, such as cognition and anxiety levels, in humans. However, environmental and other factors can be tightly controlled in mouse studies. The Collaborative Cross (CC) is a mouse genetic reference population whose common genetic and phenotypic diversity is on par with that of humans. Therefore, we leveraged the power of the CC to assess 52 behavioral measures associated with locomotor activity, anxiety level, learning and memory. This is the first application of the CC in novel object recognition tests, Morris water maze tasks, and fear conditioning tests. We found substantial continuous behavioral variations across the CC strains tested, and mapped six quantitative trait loci (QTLs) which influenced these traits, defining candidate genetic variants underlying these QTLs. Overall, our findings highlight the potential of the CC population in behavioral genetic research, while the identified genomic loci and genes driving the variation of relevant behavioral traits provide a foundation for further studies.

Zhx2 Is a Candidate Gene Underlying Oxymorphone Metabolite Brain Concentration Associated with State-Dependent Oxycodone Reward

Understanding the pharmacogenomics of opioid metabolism and behavior is vital to therapeutic success, as mutations can dramatically alter therapeutic efficacy and addiction liability. We found robust, sex-dependent BALB/c substrain differences in oxycodone behaviors and whole brain concentration of oxycodone metabolites. BALB/cJ females showed robust state-dependent oxycodone reward learning as measured via conditioned place preference when compared with the closely related BALB/cByJ substrain. Accordingly, BALB/cJ females also showed a robust increase in brain concentration of the inactive metabolite noroxycodone and the active metabolite oxymorphone compared with BALB/cByJ mice. Oxymorphone is a highly potent, full agonist at the mu opioid receptor that could enhance drug-induced interoception and state-dependent oxycodone reward learning. Quantitative trait locus (QTL) mapping in a BALB/c F2 reduced complexity cross revealed one major QTL on chromosome 15 underlying brain oxymorphone concentration that explained 32% of the female variance. BALB/cJ and BALB/cByJ differ by fewer than 10,000 variants, which can greatly facilitate candidate gene/variant identification. Hippocampal and striatal cis-expression QTL (eQTL) and exon-level eQTL analysis identified Zhx2, a candidate gene coding for a transcriptional repressor with a private BALB/cJ retroviral insertion that reduces Zhx2 expression and sex-dependent dysregulation of cytochrome P450 enzymes. Whole brain proteomics corroborated the Zhx2 eQTL and identified upregulated CYP2D11 that could increase brain oxymorphone in BALB/cJ females. To summarize, Zhx2 is a highly promising candidate gene underlying brain oxycodone metabolite levels. Future studies will validate Zhx2 and its site of action using reciprocal gene editing and tissue-specific viral manipulations in BALB/c substrains. SIGNIFICANCE STATEMENT: Our findings show that genetic variation can result in sex-specific alterations in whole brain concentration of a bioactive opioid metabolite after oxycodone administration, reinforcing the need for sex as a biological factor in pharmacogenomic studies. The cooccurrence of female-specific increased oxymorphone and state-dependent reward learning suggests that this minor yet potent and efficacious metabolite of oxycodone could increase opioid interoception and drug-cue associative learning of opioid reward, which has implications for cue-induced relapse of drug-seeking behavior and for precision pharmacogenetics.

Cocaine-Induced Locomotor Activation Differs Across Inbred Mouse Substrains

Cocaine use disorders (CUD) are devastating for affected individuals and impose a significant societal burden, but there are currently no FDA-approved therapies. The development of novel and effective treatments has been hindered by substantial gaps in our knowledge about the etiology of these disorders. The risk for developing a CUD is influenced by genetics, the environment and complex interactions between the two. Identifying specific genes and environmental risk factors that increase CUD risk would provide an avenue for the development of novel treatments. Rodent models of addiction-relevant behaviors have been a valuable tool for studying the genetics of behavioral responses to drugs of abuse. Traditional genetic mapping using genetically and phenotypically divergent inbred mice has been successful in identifying numerous chromosomal regions that influence addiction-relevant behaviors, but these strategies rarely result in identification of the causal gene or genetic variant. To overcome this challenge, reduced complexity crosses (RCC) between closely related inbred mouse strains have been proposed as a method for rapidly identifying and validating functional variants. The RCC approach is dependent on identifying phenotypic differences between substrains. To date, however, the study of addiction-relevant behaviors has been limited to very few sets of substrains, mostly comprising the C57BL/6 lineage. The present study expands upon the current literature to assess cocaine-induced locomotor activation in 20 inbred mouse substrains representing six inbred strain lineages (A/J, BALB/c, FVB/N, C3H/He, DBA/2 and NOD) that were either bred in-house or supplied directly by a commercial vendor. To our knowledge, we are the first to identify significant differences in cocaine-induced locomotor response in several of these inbred substrains. The identification of substrain differences allows for the initiation of RCC populations to more rapidly identify specific genetic variants associated with acute cocaine response. The observation of behavioral profiles that differ between mice generated in-house and those that are vendor-supplied also presents an opportunity to investigate the influence of environmental factors on cocaine-induced locomotor activity.

New Insights on Gene by Environmental Effects of Drugs of Abuse in Animal Models Using GeneNetwork

Gene-by-environment interactions are important for all facets of biology, especially behaviour. Families of isogenic strains of mice, such as the BXD strains, are excellently placed to study these interactions, as the same genome can be tested in multiple environments. BXD strains are recombinant inbred mouse strains derived from crossing two inbred strains—C57BL/6J and DBA/2J mice. Many reproducible genometypes can be leveraged, and old data can be reanalysed with new tools to produce novel insights. We obtained drug and behavioural phenotypes from Philip et al. Genes, Brain and Behaviour 2010, and reanalysed their data with new genotypes from sequencing, as well as new models (Genome-wide Efficient Mixed Model Association (GEMMA) and R/qtl2). We discovered QTLs on chromosomes 3, 5, 9, 11, and 14, not found in the original study. We reduced the candidate genes based on their ability to alter gene expression or protein function. Candidate genes included Slitrk6 and Cdk14. Slitrk6, in a Chromosome14 QTL for locomotion, was found to be part of a co-expression network involved in voluntary movement and associated with neuropsychiatric phenotypes. Cdk14, one of only three genes in a Chromosome5 QTL, is associated with handling induced convulsions after ethanol treatment, that is regulated by the anticonvulsant drug valproic acid. By using families of isogenic strains, we can reanalyse data to discover novel candidate genes involved in response to drugs of abuse.

BXD Recombinant Inbred Mice as a Model to Study Neurotoxicity

BXD recombinant inbred (RI) lines represent a genetic reference population derived from a cross between C57BL/6J mice (B6) and DBA/2J mice (D2), which through meiotic recombination events possesses recombinant chromosomes containing B6 or D2 haplotype segments. The quantitative trait loci (QTLs) are the locations of segregating genetic polymorphisms and are fundamental to understanding genetic diversity in human disease susceptibility and severity. QTL mapping represents the typical approach for identifying naturally occurring polymorphisms that influence complex phenotypes. In this process, genotypic values at markers of known genomic locations are associated with phenotypic values measured in a segregating population. Indeed, BXD RI strains provide a powerful tool to study neurotoxicity induced by different substances. In this review, we describe the use of BXD RI lines to understand the underlying mechanisms of neurotoxicity in response to ethanol and cocaine, as well as metals and pesticide exposures.

Genes identified in rodent studies of alcohol intake are enriched for heritability of human substance use

BACKGROUND

Rodent paradigms and human genome-wide association studies (GWAS) on drug use have the potential to provide biological insight into the pathophysiology of addiction.

METHODS

Using GeneWeaver, we created rodent alcohol and nicotine gene-sets derived from 19 gene expression studies on alcohol and nicotine outcomes. We partitioned the SNP-heritability of these gene-sets using four large human GWAS: 1) alcoholic drinks per week, 2) problematic alcohol use, 3) cigarettes per day and 4) smoking cessation. We benchmarked our findings with curated human alcoholism and nicotine addiction gene-sets and performed specificity analyses using other rodent gene-sets (e.g., locomotor behavior) and other human GWAS (e.g., height).

RESULTS

The rodent alcohol gene-set was enriched for heritability of drinks per week, cigarettes per day, and smoking cessation, but not problematic alcohol use. However, the rodent nicotine gene-set was not significantly associated with any of these traits. Both rodent gene-sets showed enrichment for several non-substance use GWAS, and the extent of this relationship tended to increase as a function of trait heritability. In general, larger gene-sets demonstrated more significant enrichment. Finally, when evaluating human traits with similar heritabilities, both rodent gene-sets showed greater enrichment for substance use traits.

CONCLUSION

Our results suggest that rodent gene expression studies can help to identify genes that contribute to heritability of some substance use traits in humans, yet there was less specificity than expected. We outline various limitations, interpretations and considerations for future research.

Heritable variation in locomotion, reward sensitivity and impulsive behaviors in a genetically diverse inbred mouse panel

Drugs of abuse, including alcohol and stimulants like cocaine, produce effects that are subject to individual variability, and genetic variation accounts for at least a portion of those differences. Notably, research in both animal models and human subjects point toward reward sensitivity and impulsivity as being trait characteristics that predict relatively greater positive subjective responses to stimulant drugs. Here we describe use of the eight collaborative cross (CC) founder strains and 38 (reversal learning) or 10 (all other tests) CC strains to examine the heritability of reward sensitivity and impulsivity traits, as well as genetic correlations between these measures and existing addiction-related phenotypes. Strains were all tested for activity in an open field and reward sensitivity (intake of chocolate BOOST®). Mice were then divided into two counterbalanced groups and underwent reversal learning (impulsive action and waiting impulsivity) or delay discounting (impulsive choice). CC and founder mice show significant heritability for impulsive action, impulsive choice, waiting impulsivity, locomotor activity, and reward sensitivity, with each impulsive phenotype determined to be non-correlating, independent traits. This research was conducted within the broader, inter-laboratory effort of the Center for Systems Neurogenetics of Addiction (CSNA) to characterize CC and DO mice for multiple, cocaine abuse related traits. These data will facilitate the discovery of genetic correlations between predictive traits, which will then guide discovery of genes and genetic variants that contribute to addictive behaviors.

Genetic Modulation of Initial Sensitivity to Δ9-Tetrahydrocannabinol (THC) Among the BXD Family of Mice

Cannabinoid receptor 1 activation by the major psychoactive component in cannabis, Δ9-tetrahydrocannabinol (THC), produces motor impairments, hypothermia, and analgesia upon acute exposure. In previous work, we demonstrated significant sex and strain differences in acute responses to THC following administration of a single dose (10 mg/kg, i.p.) in C57BL/6J (B6) and DBA/2J (D2) inbred mice. To determine the extent to which these differences are heritable, we quantified acute responses to a single dose of THC (10 mg/kg, i.p.) in males and females from 20 members of the BXD family of inbred strains derived by crossing and inbreeding B6 and D2 mice. Acute THC responses (initial sensitivity) were quantified as changes from baseline for: 1. spontaneous activity in the open field (mobility), 2. body temperature (hypothermia), and 3. tail withdrawal latency to a thermal stimulus (antinociception). Initial sensitivity to the immobilizing, hypothermic, and antinociceptive effects of THC varied substantially across the BXD family. Heritability was highest for mobility and hypothermia traits, indicating that segregating genetic variants modulate initial sensitivity to THC. We identified genomic loci and candidate genes, including Ndufs2, Scp2, Rps6kb1 or P70S6K, Pde4d, and Pten, that may control variation in THC initial sensitivity. We also detected strong correlations between initial responses to THC and legacy phenotypes related to intake or response to other drugs of abuse (cocaine, ethanol, and morphine). Our study demonstrates the feasibility of mapping genes and variants modulating THC responses in the BXDs to systematically define biological processes and liabilities associated with drug use and abuse.

Integration of evidence across human and model organism studies: A meeting report

The National Institute on Drug Abuse and Joint Institute for Biological Sciences at the Oak Ridge National Laboratory hosted a meeting attended by a diverse group of scientists with expertise in substance use disorders (SUDs), computational biology, and FAIR (Findability, Accessibility, Interoperability, and Reusability) data sharing. The meeting's objective was to discuss and evaluate better strategies to integrate genetic, epigenetic, and 'omics data across human and model organisms to achieve deeper mechanistic insight into SUDs. Specific topics were to (a) evaluate the current state of substance use genetics and genomics research and fundamental gaps, (b) identify opportunities and challenges of integration and sharing across species and data types, (c) identify current tools and resources for integration of genetic, epigenetic, and phenotypic data, (d) discuss steps and impediment related to data integration, and (e) outline future steps to support more effective collaboration-particularly between animal model research communities and human genetics and clinical research teams. This review summarizes key facets of this catalytic discussion with a focus on new opportunities and gaps in resources and knowledge on SUDs.

Differential effects by sex with Kmt5b loss

Lysine methyl transferase 5B (KMT5B) has been recently highlighted as a risk gene in genetic studies of neurodevelopmental disorders (NDDs), specifically, autism spectrum disorder (ASD) and intellectual disability (ID); yet, its role in the brain is not known. The goal of this work was to neurodevelopmentally characterize the effect(s) of KMT5B haploinsufficiency using a mouse model. A Kmt5b gene-trap mouse line was obtained from the Knockout Mouse Project. Wild type (WT) and heterozygous (HET) mice were subjected to a comprehensive neurodevelopmental test battery to assess reflexes, motor behavior, learning/memory, social behavior, repetitive movement, and common ASD comorbidities (obsessive compulsion, depression, and anxiety). Given the strong sex bias observed in the ASD patient population, we tested both a male and female cohort of animals and compared differences between genotypes and sexes. HET mice were significantly smaller than WT littermates starting at postnatal day 10 through young adulthood which was correlated with smaller brain size (i.e., microcephaly). This was more severe in males than females. HET male neonates also had delayed eye opening and significantly weaker reflexes than WT littermates. In young adults, significant differences between genotypes relative to anxiety, depression, fear, and extinction learning were observed. Interestingly, several sexually dimorphic differences were noted including increased repetitive grooming behavior in HET females and an increased latency to hot plate response in HET females versus a decreased latency in HET males. LAY SUMMARY: Lysine methyl transferase 5B (KMT5B) has been recently highlighted as a risk gene in neurodevelopmental disorders (NDDs), specifically, autism spectrum disorder (ASD) and intellectual disability (ID); yet its role in the brain is not known. Our study indicates that mice lacking one genomic copy of Kmt5b show deficits in neonatal reflexes, sociability, repetitive stress-induced grooming, changes in thermal pain sensing, decreased depression and anxiety, increased fear, slower extinction learning, and lower body weight, length, and brain size. Furthermore, several outcomes differed by sex, perhaps mirroring the sex bias in ASD.

Systems genetic analysis of nicotine withdrawal deficits in hippocampus-dependent learning

Cognitive deficits, such as disrupted learning, are a major symptom of nicotine withdrawal. These deficits are heritable, yet their genetic basis is largely unknown. Our lab has developed a mouse model of nicotine withdrawal deficits in learning, using chronic nicotine exposure via osmotic minipumps and fear conditioning. Here, we utilized the BXD genetic reference panel to identify genetic variants underlying nicotine withdrawal deficits in learning. Male and female mice (n = 6-11 per sex per strain, 31 strains) received either chronic saline or nicotine (6.3 mg/kg per day for 12 days), and were then tested for hippocampus-dependent learning deficits using contextual fear conditioning. Quantitative trait locus (QTL) mapping analyses using GeneNetwork identified a significant QTL on Chromosome 4 (82.13 Mb, LRS = 20.03, p < 0.05). Publicly available hippocampal gene expression data were used to identify eight positional candidates (Snacpc3, Mysm1, Rps6, Plaa, Lurap1l, Slc24a2, Hacd4, Ptprd) that overlapped with our behavioral QTL and correlated with our behavioral data. Overall, this study demonstrates that genetic factors impact cognitive deficits during nicotine withdrawal in the BXD recombinant inbred panel and identifies candidate genes for future research.

A platform for experimental precision medicine: The extended BXD mouse family

The challenge of precision medicine is to model complex interactions among DNA variants, phenotypes, development, environments, and treatments. We address this challenge by expanding the BXD family of mice to 140 fully isogenic strains, creating a uniquely powerful model for precision medicine. This family segregates for 6 million common DNA variants-a level that exceeds many human populations. Because each member can be replicated, heritable traits can be mapped with high power and precision. Current BXD phenomes are unsurpassed in coverage and include much omics data and thousands of quantitative traits. BXDs can be extended by a single-generation cross to as many as 19,460 isogenic F1 progeny, and this extended BXD family is an effective platform for testing causal modeling and for predictive validation. BXDs are a unique core resource for the field of experimental precision medicine.

Multi-omic and multi-species meta-analyses of nicotine consumption

Cross-species translational approaches to human genomic analyses are lacking. The present study uses an integrative framework to investigate how genes associated with nicotine use in model organisms contribute to the genetic architecture of human tobacco consumption. First, we created a model organism geneset by collecting results from five animal models of nicotine exposure (RNA expression changes in brain) and then tested the relevance of these genes and flanking genetic variation using genetic data from human cigarettes per day (UK BioBank N = 123,844; all European Ancestry). We tested three hypotheses: (1) DNA variation in, or around, the 'model organism geneset' will contribute to the heritability to human tobacco consumption, (2) that the model organism genes will be enriched for genes associated with human tobacco consumption, and (3) that a polygenic score based off our model organism geneset will predict tobacco consumption in the AddHealth sample (N = 1667; all European Ancestry). Our results suggested that: (1) model organism genes accounted for ~5-36% of the observed SNP-heritability in human tobacco consumption (enrichment: 1.60-31.45), (2) model organism genes, but not negative control genes, were enriched for the gene-based associations (MAGMA, H-MAGMA, SMultiXcan) for human cigarettes per day, and (3) polygenic scores based on our model organism geneset predicted cigarettes per day in an independent sample. Altogether, these findings highlight the advantages of using multiple species evidence to isolate genetic factors to better understand the etiological complexity of tobacco and other nicotine consumption.

Highlights from the Era of Open Source Web-Based Tools

High digital connectivity and a focus on reproducibility are contributing to an open science revolution in neuroscience. Repositories and platforms have emerged across the whole spectrum of subdisciplines, paving the way for a paradigm shift in the way we share, analyze, and reuse vast amounts of data collected across many laboratories. Here, we describe how open access web-based tools are changing the landscape and culture of neuroscience, highlighting six free resources that span subdisciplines from behavior to whole-brain mapping, circuits, neurons, and gene variants.

The influence of sex and strain on trace element dysregulation in the brain due to diet-induced obesity

BACKGROUND

The objective of this study was to identify interaction effects between diet, sex, and strain on trace element dysregulation and gene expression alterations due to diet-induced obesity (DIO) in the hippocampus, striatum, and midbrain.

METHODS

Male and female C57BL/6 J (B6 J) and DBA/2 J (D2 J) mice were fed either a low fat (10 % kcal) diet (LFD) or high fat (60 % kcal) diet (HFD) for 16 weeks, then assessed for trace element concentrations and gene expression patterns in the brain.

RESULTS

In the hippocampus, zinc was significantly increased by 48 % in D2 J males but decreased by 44 % in D2 J females, and divalent metal transporter 1 was substantially upregulated in B6 J males due to DIO. In the striatum, iron was significantly elevated in B6 J female mice, and ceruloplasmin was significantly upregulated in D2 J female mice due to DIO. In the midbrain, D2 J males fed a HFD had a 48 % reduction in Cu compared to the LFD group, and D2 J females had a 37 % reduction in Cu compared to the control group.

CONCLUSIONS

The alteration of trace element homeostasis and gene expression due to DIO was brain-region dependent and was highly influenced by sex and strain. A significant three-way interaction between diet, sex, and strain was discovered for zinc in the hippocampus (for mice fed a HFD, zinc increased in male D2 Js, decreased in female D2 Js, and had no effect in B6 J mice). A significant diet by sex interaction was observed for iron in the striatum (iron increased only in female mice fed a HFD). A main effect of decreased copper in the midbrain was found for the D2 J strain fed a HFD. These results emphasize the importance of considering sex and genetics as biological factors when investigating potential associations between DIO and neurodegenerative disease.

Turning strains into strengths for understanding psychiatric disorders

There is a paucity in the development of new mechanistic insights and therapeutic approaches for treating psychiatric disease. One of the major challenges is reflected in the growing consensus that risk for these diseases is not determined by a single gene, but rather is polygenic, arising from the action and interaction of multiple genes. Canonically, experimental models in mice have been designed to ascertain the relative contribution of a single gene to a disease by systematic manipulation (e.g., mutation or deletion) of a known candidate gene. Because these studies have been largely carried out using inbred isogenic mouse strains, in which there is no (or very little) genetic diversity among subjects, it is difficult to identify unique allelic variants, gene modifiers, and epigenetic factors that strongly affect the nature and severity of these diseases. Here, we review various methods that take advantage of existing genetic diversity or that increase genetic variance in mouse models to (1) strengthen conclusions of single-gene function; (2) model diversity among human populations; and (3) dissect complex phenotypes that arise from the actions of multiple genes.

Bioinformatics identification and pharmacological validation of Kcnn3/KCa2 channels as a mediator of negative affective behaviors and excessive alcohol drinking in mice

Mood disorders are often comorbid with alcohol use disorder (AUD) and play a considerable role in the development and maintenance of alcohol dependence and relapse. Because of this high comorbidity, it is necessary to determine shared and unique genetic factors driving heavy drinking and negative affective behaviors. In order to identify novel pharmacogenetic targets, a bioinformatics analysis was used to quantify the expression of amygdala K+ channel genes that covary with anxiety-related phenotypes in the well-phenotyped and fully sequenced family of BXD strains. We used a model of stress-induced escalation of drinking in alcohol-dependent mice to measure negative affective behaviors during abstinence. A pharmacological approach was used to validate the key bioinformatics findings in alcohol-dependent, stressed mice. Amygdalar expression of Kcnn3 correlated significantly with 40 anxiety-associated phenotypes. Further examination of Kcnn3 expression revealed a strong eigentrait for anxiety-like behaviors and negative correlations with binge-like and voluntary alcohol drinking. Mice treated with chronic intermittent alcohol exposure and repeated swim stress consumed more alcohol in their home cages and showed hypophagia on the novelty-suppressed feeding test during abstinence. Pharmacologically targeting Kcnn gene products with the KCa2 (SK) channel-positive modulator 1-EBIO decreased drinking and reduced feeding latency in alcohol-dependent, stressed mice. Collectively, these validation studies provide central nervous system links into the covariance of stress, negative affective behaviors, and AUD in the BXD strains. Further, the bioinformatics discovery tool is effective in identifying promising targets (i.e., KCa2 channels) for treating alcohol dependence exacerbated by comorbid mood disorders.

Machine learning-based automated phenotyping of inflammatory nocifensive behavior in mice

The discovery and development of new and potentially nonaddictive pain therapeutics requires rapid, yet clinically relevant assays of nociception in preclinical models. A reliable and scalable automated scoring system for nocifensive behavior of mice in the formalin assay would dramatically lower the time and labor costs associated with experiments and reduce experimental variability. Here, we present a method that exploits machine learning techniques for video recordings that consists of three components: key point detection, per frame feature extraction using these key points, and classification of behavior using the GentleBoost algorithm. This approach to automation is flexible as different model classifiers or key points can be used with only small losses in accuracy. The adopted system identified the behavior of licking/biting of the hind paw with an accuracy that was comparable to a human observer (98% agreement) over 111 different short videos (total 284 min) at a resolution of 1 s. To test the system over longer experimental conditions, the responses of two inbred strains, C57BL/6NJ and C57BL/6J, were recorded over 90 min post formalin challenge. The automated system easily scored over 80 h of video and revealed strain differences in both response timing and amplitude. This machine learning scoring system provides the required accuracy, consistency, and ease of use that could make the formalin assay a feasible choice for large-scale genetic studies.

Whole-Genome Sequencing of Inbred Mouse Strains Selected for High and Low Open-Field Activity

We conducted whole-genome sequencing of four inbred mouse strains initially selected for high (H1, H2) or low (L1, L2) open-field activity (OFA), and then examined strain distribution patterns for all DNA variants that differed between their BALB/cJ and C57BL/6J parental strains. Next, we assessed genome-wide sharing (3,678,826 variants) both between and within the High and Low Activity strains. Results suggested that about 10% of these DNA variants may be associated with OFA, and clearly demonstrated its polygenic nature. Finally, we conducted bioinformatic analyses of functional genomics data from mouse, rat, and human to refine previously identified quantitative trait loci (QTL) for anxiety-related measures. This combination of sequence analysis and genomic-data integration facilitated refinement of previously intractable QTL findings, and identified possible genes for functional follow-up studies.

Genetic variation regulates opioid-induced respiratory depression in mice

In the U.S., opioid prescription for treatment of pain nearly quadrupled from 1999 to 2014. The diversion and misuse of prescription opioids along with increased use of drugs like heroin and fentanyl, has led to an epidemic in addiction and overdose deaths. The most common cause of opioid overdose and death is opioid-induced respiratory depression (OIRD), a life-threatening depression in respiratory rate thought to be caused by stimulation of opioid receptors in the inspiratory-generating regions of the brain. Studies in mice have revealed that variation in opiate lethality is associated with strain differences, suggesting that sensitivity to OIRD is genetically determined. We first tested the hypothesis that genetic variation in inbred strains of mice influences the innate variability in opioid-induced responses in respiratory depression, recovery time and survival time. Using the founders of the advanced, high-diversity mouse population, the Diversity Outbred (DO), we found substantial sex and genetic effects on respiratory sensitivity and opiate lethality. We used DO mice treated with morphine to map quantitative trait loci for respiratory depression, recovery time and survival time. Trait mapping and integrative functional genomic analysis in GeneWeaver has allowed us to implicate Galnt11, an N-acetylgalactosaminyltransferase, as a gene that regulates OIRD.

Sex and Strain Variation in Initial Sensitivity and Rapid Tolerance to Δ9-Tetrahydrocannabinol

Background and Objectives: For cannabis and other drugs of abuse, initial response and/or tolerance to drug effects can predict later dependence and problematic use. Our objective is to identify sex and genetic (strain) differences in initial response and rapid tolerance to Δ9–tetrahydrocannabinol (THC), the main psychoactive ingredient in cannabis, between highly genetically divergent inbred mouse strains—C57BL/6J (B6) and DBA/2J (D2).

Experimental Approach: Sex and strain responses relative to baseline were quantified following daily exposure (i.p.) to 10 mg/kg THC or vehicle (VEH) over the course of 5 days. Dependent measures included hypothermia (decreased body temperature) and ataxia (decreased spontaneous activity in the open field), and antinociception (increase in tail withdrawal latency to a thermal stimulus). Initial sensitivity to THC was defined as the difference in response between baseline and day 1. Rapid tolerance to THC was defined as the difference in response between days 1 and 2.

Results: B6 exhibited greater THC-induced motor activity suppression and initial sensitivity to ataxia relative to the D2 strain. Females demonstrated greater levels of THC-induced hypothermia and initial sensitivity relative to males. Higher levels of THC-induced antinociception and initial sensitivity were observed for D2 relative to B6. Rapid tolerance to THC was observed for hypothermia and antinociception. Much less tolerance was observed for THC-induced ataxia. D2 exhibited rapid tolerance to THC-induced hypothermia and antinociception at time points associated with peak THC initial response. Likewise, at the peak initial THC response time point, females demonstrated greater levels of rapid tolerance to hypothermic effects relative to males.

Conclusions: Both sex and genetic factors drive variation in initial response and rapid tolerance to the ataxic, antinociceptive, and hypothermic effects of THC. As these traits directly result from THC activation of the cannabinoid receptor 1, gene variants between B6 and D2 in cannabinoid signaling pathways are likely to mediate strain differences in response to THC.

Neuroinflammation, Pain and Depression: An Overview of the Main Findings

Chronic pain is a serious public health problem with a strong affective-motivational component that makes it difficult to treat. Most patients with chronic pain suffer from severe depression; hence, both conditions coexist and exacerbate one another. Brain inflammatory mediators are critical for maintaining depression-pain syndrome and could be substrates for it. The goal of our paper was to review clinical and preclinical findings to identify the neuroinflammatory profile associated with the cooccurrence of pain and depression. In addition, we aimed to explore the regulatory effect of neuronal reorganization on the inflammatory response in pain and depression. We conducted a quantitative review supplemented by manual screening. Our results revealed inflammatory signatures in different preclinical models and clinical articles regarding depression-pain syndrome. We also identified that improvements in depressive symptoms and amelioration of pain can be modulated through direct targeting of inflammatory mediators, such as cytokines and molecular inhibitors of the inflammatory cascade. Additionally, therapeutic targets that improve and regulate the synaptic environment and its neurotransmitters may act as anti-inflammatory compounds, reducing local damage-associated molecular patterns and inhibiting the activation of immune and glial cells. Taken together, our data will help to better elucidate the neuroinflammatory profile in pain and depression and may help to identify pharmacological targets for effective management of depression-pain syndrome.

Do gene expression findings from mouse models of cocaine use recapitulate human cocaine use disorder in reward circuitry?

Animal models of drug use have investigated possible mechanisms governing human substance use traits for over 100 years. Most cross-species research on drug use/addiction examines behavioral overlap, but studies assessing neuromolecular (e.g. RNA) correspondence are lacking. Our study utilized transcriptome-wide data from the hippocampus and ventral tegmental area (VTA)/midbrain from a total of 35 human males with cocaine use disorder/controls and 49 male C57BL/6J cocaine/saline administering/exposed mice. We hypothesized differential expressed genes and systems of co-expressed genes (gene networks) would show appreciable overlap across mouse cocaine self-administration and human cocaine use disorder. We found modest, but significant relationships between differentially expressed genes associated with cocaine self-administration (short access) and cocaine use disorder within reward circuitry. Differentially expressed genes underlying models of acute cocaine exposure (cocaine), context re-exposure and cocaine + context re-exposure were not consistently associated with human CUD across brain regions. Investigating systems of co-expressed genes, we found several validated gene networks with weak to moderate conservation between cocaine/saline self-administering mice and disordered cocaine users/controls. The most conserved hippocampal and VTA gene networks demonstrated substantial overlap (2029 common genes) and included both novel and previously implicated targets for cocaine use/addiction. Lastly, we conducted (expression-based) phenome-wide association studies of the nine common hub genes across conserved gene networks. Common hub genes were associated with dopamine/serotonin function, cocaine self-administration and other relevant mouse traits. Overall, our study pinpointed and characterized conserved brain-related RNA patterns across mouse cocaine self-administration and human cocaine use disorder. We offer recommendations for future research and add to the dialogue surrounding pre-clinical animal research for human disease.

Discovery of a Role for Rab3b in Habituation and Cocaine Induced Locomotor Activation in Mice Using Heterogeneous Functional Genomic Analysis

Substance use disorders are prevalent and present a tremendous societal cost but the mechanisms underlying addiction behavior are poorly understood and few biological treatments exist. One strategy to identify novel molecular mechanisms of addiction is through functional genomic experimentation. However, results from individual experiments are often noisy. To address this problem, the convergent analysis of multiple genomic experiments can discern signal from these studies. In the present study, we examine genetic loci that modulate the locomotor response to cocaine identified in the recombinant inbred (BXD RI) genetic reference population. We then applied the GeneWeaver software system for heterogeneous functional genomic analysis to integrate and aggregate multiple studies of addiction genomics, resulting in the identification of Rab3b as a functional correlate of the locomotor response to cocaine in rodents. This gene encodes a member of the RAB family of Ras-like GTPases known to be involved in trafficking of secretory and endocytic vesicles in eukaryotic cells. The convergent evidence for a role of Rab3b includes co-occurrence in previously published genetic mapping studies of cocaine related behaviors; methamphetamine response and cocaine- and amphetamine-regulated transcript prepropeptide (Cartpt) transcript abundance; evidence related to other addictive substances; density of polymorphisms; and its expression pattern in reward pathways. To evaluate this finding, we examined the effect of RAB3 complex perturbation in cocaine response. B6;129-Rab3btm1Sud Rab3ctm1sud Rab3dtm1sud triple null mice (Rab3bcd -/-) exhibited significant deficits in habituation, and increased acute and repeated cocaine responses. This previously unidentified mechanism of the behavioral predisposition and response to cocaine is an example of many that can be identified and validated using aggregate genomic studies.

Using a mouse model to gain insights into developmental coordination disorder

Motor impairments are a common feature of many neurodevelopmental disorders; in fact, over 50% of children with Attentional Deficit Hyperactivity Disorder or Autism Spectrum Disorder may have a co-occurring diagnosis of developmental coordination disorder (DCD). DCD is a neurodevelopmental disorder of unknown etiology that affects motor coordination and learning, significantly impacting a child's ability to carry out everyday activities. Animal models play an important role in scientific investigation of behaviour and the mechanisms and processes that are involved in control of motor actions. The purpose of this paper is to present an approach in the mouse directed to gain behavioral and genetic insights into DCD that is designed with high face validity, construct validity and predictive validity. Pre-clinical and clinical expertise is used to establish a set of scientific criteria that the model will meet in order to investigate the potential underlying causes of DCD.

Prospects for finding the mechanisms of sex differences in addiction with human and model organism genetic analysis

Despite substantial evidence for sex differences in addiction epidemiology, addiction-relevant behaviors and associated neurobiological phenomena, the mechanisms and implications of these differences remain unknown. Genetic analysis in model organism is a potentially powerful and effective means of discovering the mechanisms that underlie sex differences in addiction. Human genetic studies are beginning to show precise risk variants that influence the mechanisms of addiction but typically lack sufficient power or neurobiological mechanistic access, particularly for the discovery of the mechanisms that underlie sex differences. Our thesis in this review is that genetic variation in model organisms are a promising approach that can complement these investigations to show the biological mechanisms that underlie sex differences in addiction.

Mouse Systems Genetics as a Prelude to Precision Medicine

Mouse models have been instrumental in understanding human disease biology and proposing possible new treatments. The precise control of the environment and genetic composition of mice allows more rigorous observations, but limits the generalizability and translatability of the results into human applications. In the era of precision medicine, strategies using mouse models have to be revisited to effectively emulate human populations. Systems genetics is one promising paradigm that may promote the transition to novel precision medicine strategies. Here, we review the state-of-the-art resources and discuss how mouse systems genetics helps to understand human diseases and to advance the development of precision medicine, with an emphasis on the existing resources and strategies.

A Cross-Species Systems Genetics Analysis Links APBB1IP as a Candidate for Schizophrenia and Prepulse Inhibition

Background: Prepulse inhibition (PPI) of the startle response is a highly conserved form of sensorimotor gating, disruption of which is found in schizophrenia patients and their unaffected first-degree relatives. PPI can be measured in many species, and shows considerable phenotypic variation between and within rodent models. This makes PPI a useful endophenotype. Genome-wide association studies (GWAS) have been carried out to identify genetic variants underlying schizophrenia, and these suggest that schizophrenia is highly polygenic. GWAS have been unable to account for the high heritability of schizophrenia seen in family studies, partly because of the low power of GWAS due to multiple comparisons. By contrast, complementary mouse model linkage studies often have high statistical power to detect variants for behavioral traits but lower resolution, producing loci that include tens or hundreds of genes. To capitalize on the advantages of both GWAS and genetic mouse models, our study uses a cross-species approach to identify novel genes associated with PPI regulation, which thus may contribute to the PPI deficits seen in schizophrenia.

Results: Using experimental data from the recombinant inbred (RI) mouse panel BXD, we identified two significant loci affecting PPI. These genomic regions contain genetic variants which influence PPI in mice and are therefore candidates that may be influencing aspects of schizophrenia in humans. We next investigated these regions in whole-genome data from the Psychiatric Genomics Consortium (PGC) schizophrenia GWAS and identify one novel candidate gene (ABPP1IP) that was significantly associated with PPI in mice and risk of schizophrenia in humans. A systems genetics approach demonstrates that APBB1IP coexpresses with several other genes related to schizophrenia in several brain regions. Gene coexpression and enrichment analysis shows clear links between APBB1IP and the immune system.

Conclusion: The combination of human GWAS and mouse quantitative trait loci (QTL) from some of the largest study systems available has enabled us to identify a novel gene, APBB1IP, which influences schizophrenia in humans and PPI in mice.

Discovery of early life stress interacting and sex-specific quantitative trait loci impacting cocaine responsiveness

BACKGROUND AND PURPOSE

Addiction vulnerability involves complex gene X environment interactions leading to a pathological response to drugs. Identification of the genes involved in these interactions is an important step in understanding the underlying neurobiology and rarely have such analyses examined sex-specific influences. To dissect this interaction, we examined the impact of prenatal stress (PNS) on cocaine responsiveness in male and female mice of the BXD recombinant inbred panel.

EXPERIMENTAL APPROACH

BXD strains were subjected to timed mating and assigned to PNS or control groups. PNS dams were subjected to restraint stress (1-hr restraint, three times daily) starting between embryonic day (E) 11 and 14 and continued until parturition. Adult male and female, control and PNS offspring were tested for locomotor response to initial and repeated cocaine injections (sensitization) as well as cocaine-induced conditioned place preference (CPP).

KEY RESULTS

Strain, PNS, and sex interacted to modulate initial and sensitized cocaine-induced locomotion, as well as CPP. Moreover, a quantitative trait locus (QTL) interacting with PNS regulating initial locomotor response to cocaine (chromosome X, 37.91 to 50.95 Mb) was identified. Also PNS-independent, female-specific QTLs regulating CPP (chromosome 11, 65.50 to 81.31 Mb) and sensitized cocaine-induced locomotion (chromosome 16, 95.79 to 98.32 Mb) were identified. Publicly available mRNA expression data were utilized to identify cis-eQTL and transcript covariation with the behavioural phenotype to prioritize candidate genes; including Aifm1.

CONCLUSIONS AND IMPLICATIONS

These QTL encompass genes that may moderate genetic susceptibility to PNS and interact with sex to determine adult responsiveness to cocaine and addiction vulnerability.

LINKED ARTICLES

This article is part of a themed section on The Importance of Sex Differences in Pharmacology Research. To view the other articles in this section visit http://onlinelibrary.wiley.com/doi/10.1111/bph.v176.21/issuetoc.

Exploring the involvement of Tac2 in the mouse hippocampal stress response through gene networking

Tachykinin 2 (Tac2) is expressed in a number of areas throughout the brain, including the hippocampus. However, knowledge about its function has been only well explored in the hypothalamus in the context of reproductive health. In this study, we identified and validated increased hippocampal Tac2 mRNA expression in response to chronic mild stress in mice. Expression quantitative trait locus (eQTL) analysis showed Tac2 is cis-regulated in the hippocampus. Using a systems genetics approach, we constructed a Tac2 co-expression network to better understand the relationship between Tac2 and the hippocampal stress response. Our network identified 69 total genes associated with Tac2, several of which encode major neuropeptides involved in hippocampal stress signaling as well as critical genes for producing neural plasticity, indicating that Tac2 is involved in these processes. Pathway analysis for the member of Tac2 gene network revealed a strong connection between Tac2 and neuroactive ligand-receptor interaction, calcium signaling pathway, as well as cardiac muscle contraction. In addition, we also identified 46 stress-related phenotypes, specifically fear conditioning response, that were significantly correlated with Tac2 expression. Our results provide evidence for Tac2 as a strong candidate gene who likely plays a role in hippocampal stress processing and neural plasticity.

Knockdown of heterochromatin protein 1 binding protein 3 recapitulates phenotypic, cellular, and molecular features of aging

Identifying genetic factors that modify an individual's susceptibility to cognitive decline in aging is critical to understanding biological processes involved and mitigating risk associated with a number of age‐related disorders. Recently, heterochromatin protein 1 binding protein 3 (Hp1bp3) was identified as a mediator of cognitive aging. Here, we provide a mechanistic explanation for these findings and show that targeted knockdown of Hp1bp3 in the hippocampus by 50%–75% is sufficient to induce cognitive deficits and transcriptional changes reminiscent of those observed in aging and Alzheimer's disease brains. Specifically, neuroinflammatory‐related pathways become activated following Hp1bp3 knockdown in combination with a robust decrease in genes involved in synaptic activity and neuronal function. To test the hypothesis that Hp1bp3 mediates susceptibility to cognitive deficits via a role in neuronal excitability, we performed slice electrophysiology demonstrate transcriptional changes after Hp1bp3 knockdown manifest functionally as a reduction in hippocampal neuronal intrinsic excitability and synaptic plasticity. In addition, as Hp1bp3 is a known mediator of miRNA biogenesis, here we profile the miRNA transcriptome and identify mir‐223 as a putative regulator of a portion of observed mRNA changes, particularly those that are inflammatory‐related. In summary, work here identifies Hp1bp3 as a critical mediator of aging‐related changes at the phenotypic, cellular, and molecular level and will help inform the development of therapeutics designed to target either Hp1bp3 or its downstream effectors in order to promote cognitive longevity.

Genome wide association analysis in a mouse advanced intercross line

The LG/J x SM/J advanced intercross line of mice (LG x SM AIL) is a multigenerational outbred population. High minor allele frequencies, a simple genetic background, and the fully sequenced LG and SM genomes make it a powerful population for genome-wide association studies. Here we use 1,063 AIL mice to identify 126 significant associations for 50 traits relevant to human health and disease. We also identify thousands of cis- and trans-eQTLs in the hippocampus, striatum, and prefrontal cortex of ~200 mice. We replicate an association between locomotor activity and Csmd1, which we identified in an earlier generation of this AIL, and show that Csmd1 mutant mice recapitulate the locomotor phenotype. Our results demonstrate the utility of the LG x SM AIL as a mapping population, identify numerous novel associations, and shed light on the genetic architecture of mammalian behavior.

Cross-species molecular dissection across alcohol behavioral domains

This review summarizes the proceedings of a symposium presented at the "Alcoholism and Stress: A Framework for Future Treatment Strategies" conference held in Volterra, Italy on May 9-12, 2017. Psychiatric diseases, including alcohol-use disorders (AUDs), are influenced through complex interactions of genes, neurobiological pathways, and environmental influences. A better understanding of the common neurobiological mechanisms underlying an AUD necessitates an integrative approach, involving a systematic assessment of diverse species and phenotype measures. As part of the World Congress on Stress and Alcoholism, this symposium provided a detailed account of current strategies to identify mechanisms underlying the development and progression of AUDs. Dr. Sean Farris discussed the integration and organization of transcriptome and postmortem human brain data to identify brain regional- and cell type-specific differences related to excessive alcohol consumption that are conserved across species. Dr. Brien Riley presented the results of a genome-wide association study of DSM-IV alcohol dependence; although replication of genetic associations with alcohol phenotypes in humans remains challenging, model organism studies show that COL6A3, KLF12, and RYR3 affect behavioral responses to ethanol, and provide substantial evidence for their role in human alcohol-related traits. Dr. Rob Williams expanded upon the systematic characterization of extensive genetic-genomic resources for quantifying and clarifying phenotypes across species that are relevant to precision medicine in human disease. The symposium concluded with Dr. Robert Hitzemann's description of transcriptome studies in a mouse model selectively bred for high alcohol ("binge-like") consumption and a non-human primate model of long-term alcohol consumption. Together, the different components of this session provided an overview of systems-based approaches that are pioneering the experimental prioritization and validation of novel genes and gene networks linked with a range of behavioral phenotypes associated with stress and AUDs.