Hnrnph1 Is A Quantitative Trait Gene for Methamphetamine Sensitivity

Atp1a2 and Kcnj9 Are Candidate Genes Underlying Sensitivity to Oxycodone-Induced Locomotor Activation and Withdrawal-Induced Anxiety-Like Behaviors in C57BL/6 Substrains

Opioid use disorder is heritable, yet its genetic etiology is largely unknown. C57BL/6J and C57BL/6NJ mouse substrains exhibit phenotypic diversity in the context of limited genetic diversity which together can facilitate genetic discovery. Here, we found C57BL/6NJ mice were less sensitive to oxycodone (OXY)-induced locomotor activation versus C57BL/6J mice in a conditioned place preference paradigm. Narrow-sense heritability of OXY-induced locomotor activity traits ranged from 0.22 to 0.31, implicating suitability for genetic analysis. Quantitative trait locus (QTL) mapping in an F2 cross identified a chromosome 1 QTL explaining 7%-12% of the variance in OXY locomotion and anxiety-like withdrawal in the elevated plus maze. A second QTL for EPM withdrawal behavior on chromosome 5 near Gabra2 (alpha-2 subunit of GABA-A receptor) explained 9% of the variance. To narrow the chromosome 1 locus, we generated recombinant lines spanning 163-181 Mb, captured the QTL for OXY locomotor traits and withdrawal, and fine-mapped a 2.45-Mb region (170.16-172.61 Mb). Transcriptome analysis identified five, localized striatal cis-eQTL transcripts and two were confirmed at the protein level (KCNJ9, ATP1A2). Kcnj9 codes for a potassium channel (GIRK3) that is a major effector of mu opioid receptor signaling. Atp1a2 codes for a subunit of a Na+/K+ ATPase enzyme that regulates neuronal excitability and shows functional adaptations following chronic opioid administration. To summarize, we identified two candidate genes underlying the physiological and behavioral properties of opioids, with direct preclinical relevance to investigators employing these widely used substrains and clinical relevance to human genetic studies of opioid use disorder.

Atp1a2 and Kcnj9 are candidate genes underlying sensitivity to oxycodone-induced locomotor activation and withdrawal-induced anxiety-like behaviors in C57BL/6 substrains

Opioid use disorder is heritable, yet its genetic etiology is largely unknown. C57BL/6J and C57BL/6NJ mouse substrains exhibit phenotypic diversity in the context of limited genetic diversity which together can facilitate genetic discovery. Here, we found C57BL/6NJ mice were less sensitive to oxycodone (OXY)-induced locomotor activation versus C57BL/6J mice in a conditioned place preference paradigm. Narrow-sense heritability was estimated at 0.22-0.31, implicating suitability for genetic analysis. Quantitative trait locus (QTL) mapping in an F2 cross identified a chromosome 1 QTL explaining 7-12% of the variance in OXY locomotion and anxiety-like withdrawal in the elevated plus maze. A second QTL for EPM withdrawal behavior on chromosome 5 near Gabra2 (alpha-2 subunit of GABA-A receptor) explained 9% of the variance. To narrow the chromosome 1 locus, we generated recombinant lines spanning 163-181 Mb, captured the QTL for OXY locomotor traits and withdrawal, and fine-mapped a 2.45-Mb region (170.16-172.61 Mb). Transcriptome analysis identified five, localized striatal cis-eQTL transcripts and two were confirmed at the protein level (KCNJ9, ATP1A2). Kcnj9 codes for a potassium channel (GIRK3) that is a major effector of mu opioid receptor signaling. Atp1a2 codes for a subunit of a Na+/K+ ATPase enzyme that regulates neuronal excitability and shows functional adaptations following chronic opioid administration. To summarize, we identified two candidate genes underlying the physiological and behavioral properties of opioids, with direct preclinical relevance to investigators employing these widely used substrains and clinical relevance to human genetic studies of opioid use disorder.

Decreased myelin-related gene expression in the nucleus accumbens during spontaneous neonatal opioid withdrawal in the absence of long-term behavioral effects in adult outbred CFW mice

Prenatal opioid exposure is a major health concern in the United States, with the incidence of neonatal opioid withdrawal syndrome (NOWS) escalating in recent years. NOWS occurs upon cessation of in utero opioid exposure and is characterized by increased irritability, disrupted sleep patterns, high-pitched crying, and dysregulated feeding. The main pharmacological strategy for alleviating symptoms is treatment with replacement opioids. The neural mechanisms mediating NOWS and the long-term neurobehavioral effects are poorly understood. We used a third trimester-approximate model in which neonatal outbred pups (Carworth Farms White; CFW) were administered once-daily morphine (15 mg/kg, s.c.) from postnatal day (P) day 1 through P14 and were then assessed for behavioral and transcriptomic adaptations within the nucleus accumbens (NAc) on P15. We also investigated the long-term effects of perinatal morphine exposure on adult learning and reward sensitivity. We observed significant weight deficits, spontaneous thermal hyperalgesia, and altered ultrasonic vocalization (USV) profiles following repeated morphine and during spontaneous withdrawal. Transcriptome analysis of NAc from opioid-withdrawn P15 neonates via bulk mRNA sequencing identified an enrichment profile consistent with downregulation of myelin-associated transcripts. Despite the neonatal behavioral and molecular effects, there were no significant long-term effects of perinatal morphine exposure on adult spatial memory function in the Barnes Maze, emotional learning in fear conditioning, or in baseline or methamphetamine-potentiated reward sensitivity as measured via intracranial self-stimulation. Thus, the once daily third trimester-approximate exposure regimen, while inducing NOWS model traits and significant transcriptomic effects in neonates, had no significant long-term effects on adult behaviors.

The HNRNPF/H RNA binding proteins and disease

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

Decreased myelin-related gene expression in the nucleus accumbens during spontaneous neonatal opioid withdrawal in the absence of long-term behavioral effects in adult outbred CFW mice

Prenatal opioid exposure is a major health concern in the United States, with the incidence of neonatal opioid withdrawal syndrome (NOWS) escalating in recent years. NOWS occurs upon cessation of in utero opioid exposure and is characterized by increased irritability, disrupted sleep patterns, high-pitched crying, and dysregulated feeding. The main pharmacological strategy for alleviating symptoms is treatment with replacement opioids. The neural mechanisms mediating NOWS and the long-term neurobehavioral effects are poorly understood. We used a third trimester-approximate model in which neonatal outbred pups (Carworth Farms White; CFW) were administered once-daily morphine (15 mg/kg, s.c.) from postnatal day (P) day 1 through P14 and were then assessed for behavioral and transcriptomic adaptations within the nucleus accumbens (NAc) on P15. We also investigated the long-term effects of perinatal morphine exposure on adult learning and reward sensitivity. We observed significant weight deficits, spontaneous thermal hyperalgesia, and altered ultrasonic vocalization (USV) profiles following repeated morphine and during spontaneous withdrawal. Transcriptome analysis of NAc from opioid-withdrawn P15 neonates via bulk mRNA sequencing identified an enrichment profile consistent with downregulation of myelin-associated transcripts. Despite the neonatal behavioral and molecular effects, there were no significant long-term effects of perinatal morphine exposure on adult spatial memory function in the Barnes Maze, emotional learning in fear conditioning, or in baseline or methamphetamine-potentiated reward sensitivity as measured via intracranial self-stimulation. Thus, the once daily third trimester-approximate exposure regimen, while inducing NOWS model traits and significant transcriptomic effects in neonates, had no significant long-term effects on adult behaviors.

HIGHLIGHTS

We replicated some NOWS model traits via 1x-daily morphine (P1-P14).We found a downregulation of myelination genes in nucleus accumbens on P15.There were no effects on learning/memory or reward sensitivity in adults.

Deletion of Cryab increases the vulnerability of mice to the addiction-like effects of the cannabinoid JWH-018 via upregulation of striatal NF-κB expression

Synthetic cannabinoids have exhibited unpredictable abuse liabilities, especially self-administration (SA) responses in normal rodent models, despite seemingly inducing addiction-like effects in humans. Thus, an efficient pre-clinical model must be developed to determine cannabinoid abuse potential in animals and describe the mechanism that may mediate cannabinoid sensitivity. The Cryab knockout (KO) mice were recently discovered to be potentially sensitive to the addictive effects of psychoactive drugs. Herein, we examined the responses of Cryab KO mice to JWH-018 using SA, conditioned place preference, and electroencephalography. Additionally, the effects of repeated JWH-018 exposure on endocannabinoid- and dopamine-related genes in various addiction-associated brain regions were examined, along with protein expressions involving neuroinflammation and synaptic plasticity. Cryab KO mice exhibited greater cannabinoid-induced SA responses and place preference, along with divergent gamma wave alterations, compared to wild-type (WT) mice, implying their higher sensitivity to cannabinoids. Endocannabinoid- or dopamine-related mRNA expressions and accumbal dopamine concentrations after repeated JWH-018 exposure were not significantly different between the WT and Cryab KO mice. Further analyses revealed that repeated JWH-018 administration led to possibly greater neuroinflammation in Cryab KO mice, which may arise from upregulated NF-κB, accompanied by higher expressions of synaptic plasticity markers, which might have contributed to the development of cannabinoid addiction-related behavior in Cryab KO mice. These findings signify that increased neuroinflammation via NF-κB may mediate the enhanced addiction-like responses of Cryab KO mice to cannabinoids. Altogether, Cryab KO mice may be a potential model for cannabinoid abuse susceptibility.

Hunting for Genes Underlying Emotionality in the Laboratory Rat: Maps, Tools and Traps

Scientists have systematically investigated the hereditary bases of behaviors since the 19th century, moved by either evolutionary questions or clinically-motivated purposes. The pioneer studies on the genetic selection of laboratory animals had already indicated, one hundred years ago, the immense complexity of analyzing behaviors that were influenced by a large number of small-effect genes and an incalculable amount of environmental factors. Merging Mendelian, quantitative and molecular approaches in the 1990s made it possible to map specific rodent behaviors to known chromosome regions. From that point on, Quantitative Trait Locus (QTL) analyses coupled with behavioral and molecular techniques, which involved in vivo isolation of relevant blocks of genes, opened new avenues for gene mapping and characterization. This review examines the QTL strategy applied to the behavioral study of emotionality, with a focus on the laboratory rat. We discuss the challenges, advances and limitations of the search for Quantitative Trait Genes (QTG) playing a role in regulating emotionality. For the past 25 years, we have marched the long journey from emotionality-related behaviors to genes. In this context, our experiences are used to illustrate why and how one should move forward in the molecular understanding of complex psychiatric illnesses. The promise of exploring genetic links between immunological and emotional responses are also discussed. New strategies based on humans, rodents and other animals (such as zebrafish) are also acknowledged, as they are likely to allow substantial progress to be made in the near future.

HNRNPH1 destabilizes the G-quadruplex structures formed by G-rich RNA sequences that regulate the alternative splicing of an oncogenic fusion transcript

In the presence of physiological monovalent cations, thousands of RNA G-rich sequences can form parallel G-quadruplexes (G4s) unless RNA-binding proteins inhibit, destabilize, or resolve the formation of such secondary RNA structures. Here, we have used a disease-relevant model system to investigate the biophysical properties of the RNA-binding protein HNRNPH1's interaction with G-rich sequences. We demonstrate the importance of two EWSR1-exon 8 G-rich regions in mediating the exclusion of this exon from the oncogenic EWS-FLI1 transcripts expressed in a subset of Ewing sarcomas, using complementary analysis of tumor data, long-read sequencing, and minigene studies. We determined that HNRNPH1 binds the EWSR1-exon 8 G-rich sequences with low nM affinities irrespective of whether in a non-G4 or G4 state but exhibits different kinetics depending on RNA structure. Specifically, HNRNPH1 associates and dissociates from G4-folded RNA faster than the identical sequences in a non-G4 state. Importantly, we demonstrate using gel shift and spectroscopic assays that HNRNPH1, particularly the qRRM1-qRRM2 domains, destabilizes the G4s formed by the EWSR1-exon 8 G-rich sequences in a non-catalytic fashion. Our results indicate that HNRNPH1's binding of G-rich sequences favors the accumulation of RNA in a non-G4 state and that this contributes to its regulation of RNA processing.

Genome-wide association mapping of ethanol sensitivity in the Diversity Outbred mouse population

BACKGROUND

A strong predictor for the development of alcohol use disorder (AUD) is altered sensitivity to the intoxicating effects of alcohol. Individual differences in the initial sensitivity to alcohol are controlled in part by genetic factors. Mice offer a powerful tool to elucidate the genetic basis of behavioral and physiological traits relevant to AUD, but conventional experimental crosses have only been able to identify large chromosomal regions rather than specific genes. Genetically diverse, highly recombinant mouse populations make it possible to observe a wider range of phenotypic variation, offer greater mapping precision, and thus increase the potential for efficient gene identification.

METHODS

We have taken advantage of the Diversity Outbred (DO) mouse population to identify and precisely map quantitative trait loci (QTL) associated with ethanol sensitivity. We phenotyped 798 male J:DO mice for three measures of ethanol sensitivity: ataxia, hypothermia, and loss of the righting response. We used high-density MegaMUGA and GigaMUGA to obtain genotypes ranging from 77,808 to 143,259 SNPs. We also performed RNA sequencing in striatum to map expression QTLs and identify gene expression-trait correlations. We then applied a systems genetic strategy to identify narrow QTLs and construct the network of correlations that exists between DNA sequence, gene expression values, and ethanol-related phenotypes to prioritize our list of positional candidate genes.

RESULTS

We observed large amounts of phenotypic variation with the DO population and identified suggestive and significant QTLs associated with ethanol sensitivity on chromosomes 1, 2, and 16. The implicated regions were narrow (4.5-6.9 Mb in size) and each QTL explained ~4-5% of the variance.

CONCLUSIONS

Our results can be used to identify alleles that contribute to AUD in humans, elucidate causative biological mechanisms, or assist in the development of novel therapeutic interventions.

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.

Differentially Expressed Genes in Period 2-Overexpressing Mice Striatum May Underlie Their Lower Sensitivity to Methamphetamine Addiction-Like Behavior

Previous reports have demonstrated that genetic mechanisms greatly mediate responses to drugs of abuse, including methamphetamine (METH). The circadian gene Period 2 (Per2) has been previously associated with differential responses towards METH in mice. While the behavioral consequences of eliminating Per2 have been illustrated previously, Per2 overexpression has not yet been comprehensively described; although, Per2-overexpressing (Per2 OE) mice previously showed reduced sensitivity towards METH-induced addiction-like behaviors. To further elucidate this distinct behavior of Per2 OE mice to METH, we identified possible candidate biomarkers by determining striatal differentially expressed genes (DEGs) in both drug-naïve and METH-treated Per2 OE mice relative to wild-type (WT), through RNA sequencing. Of the several DEGs in drug naïve Per2 OE mice, we identified six genes that were altered after repeated METH treatment in WT mice, but not in Per2 OE mice. These results, validated by quantitative real-time polymerase chain reaction, could suggest that the identified DEGs might underlie the previously reported weaker METH-induced responses of Per2 OE mice compared to WT. Gene network analysis also revealed that Asic3, Hba-a1, and Rnf17 are possibly associated with Per2 through physical interactions and predicted correlations, and might potentially participate in addiction. Inhibiting the functional protein of Asic3 prior to METH administration resulted in the partial reduction of METH-induced conditioned place preference in WT mice, supporting a possible involvement of Asic3 in METH-induced reward. Although encouraging further investigations, our findings suggest that these DEGs, including Asic3, may play significant roles in the lower sensitivity of Per2 OE mice to METH.

Inflammatory-miR-301a circuitry drives mTOR and Stat3-dependent PSC activation in chronic pancreatitis and PanIN

Activated pancreatic stellate cells (PSCs) are the main cells involved in chronic pancreatitis and pancreatic intraepithelial neoplasia lesion (PanIN). Fine-tuning the precise molecular targets in PSC activation might help the development of PSC-specific therapeutic strategies to tackle progression of pancreatic cancer-related fibrosis. miR-301a is a pro-inflammatory microRNA known to be activated by multiple inflammatory factors in the tumor stroma. Here, we show that miR-301a is highly expressed in activated PSCs in mice, sustained tissue fibrosis in caerulein-induced chronic pancreatitis, and accelerated PanIN formation. Genetic ablation of miR-301a reduced pancreatic fibrosis in mouse models with chronic pancreatitis and PanIN. Cell proliferation and activation of PSCs was inhibited by downregulation of miR-301a via two of its targets, Tsc1 and Gadd45g. Moreover, aberrant PSC expression of miR-301a and Gadd45g restricted the interplay between PSCs and pancreatic cancer cells in tumorigenesis. Our findings suggest that miR-301a activates two major cell proliferation pathways, Tsc1/mTOR and Gadd45g/Stat3, in vivo, to facilitate development of inflammatory-induced PanIN and maintenance of PSC activation and desmoplasia in pancreatic cancer.

Genetic basis of thermal nociceptive sensitivity and brain weight in a BALB/c reduced complexity cross

Thermal nociception involves the transmission of temperature-related noxious information from the periphery to the CNS and is a heritable trait that could predict transition to persistent pain. Rodent forward genetics complement human studies by controlling genetic complexity and environmental factors, analysis of end point tissue, and validation of variants on appropriate genetic backgrounds. Reduced complexity crosses between nearly identical inbred substrains with robust trait differences can greatly facilitate unbiased discovery of novel genes and variants. We found BALB/cByJ mice showed enhanced sensitivity on the 53.5°C hot plate and mechanical stimulation in the von Frey test compared to BALB/cJ mice and replicated decreased gross brain weight in BALB/cByJ versus BALB/cJ. We then identified a quantitative trait locus (QTL) on chromosome 13 for hot plate sensitivity (LOD = 10.7; p < 0.001; peak = 56 Mb) and a QTL for brain weight on chromosome 5 (LOD = 8.7; p < 0.001). Expression QTL mapping of brain tissues identified H2afy (56.07 Mb) as the top transcript with the strongest association at the hot plate locus (FDR = 0.0002) and spliceome analysis identified differential exon usage within H2afy associated with the same locus. Whole brain proteomics further supported decreased H2AFY expression could underlie enhanced hot plate sensitivity, and identified ACADS as a candidate for reduced brain weight. To summarize, a BALB/c reduced complexity cross combined with multiple-omics approaches facilitated identification of candidate genes underlying thermal nociception and brain weight. These substrains provide a powerful, reciprocal platform for future validation of candidate variants.

A quantitative trait variant in Gabra2 underlies increased methamphetamine stimulant sensitivity

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.

Sex Differences in Behavioral and Brainstem Transcriptomic Neuroadaptations following Neonatal Opioid Exposure in Outbred Mice

The opioid epidemic led to an increase in the number of neonatal opioid withdrawal syndrome (NOWS) cases in infants born to opioid-dependent mothers. Hallmark features of NOWS include weight loss, severe irritability, respiratory problems, and sleep fragmentation. Mouse models provide an opportunity to identify brain mechanisms that contribute to NOWS. Neonatal outbred Swiss Webster Cartworth Farms White (CFW) mice were administered morphine (15 mg/kg, s.c.) twice daily from postnatal day 1 (P1) to P14, an approximation of the third trimester of human gestation. Female and male mice underwent behavioral testing on P7 and P14 to determine the impact of opioid exposure on anxiety and pain sensitivity. Ultrasonic vocalizations (USVs) and daily body weights were also recorded. Brainstems containing pons and medulla were collected during morphine withdrawal on P14 for RNA sequencing. Morphine induced weight loss from P2 to P14, which persisted during adolescence (P21) and adulthood (P50). USVs markedly increased at P7 in females, emerging earlier than males. On P7 and P14, both morphine-exposed female and male mice displayed hyperalgesia on the hot plate and tail-flick assays, with females showing greater hyperalgesia than males. Morphine-exposed mice exhibited increased anxiety-like behavior in the open-field arena on P21. Transcriptome analysis of the brainstem, an area implicated in opioid withdrawal and NOWS, identified pathways enriched for noradrenergic signaling in females and males. We also found sex-specific pathways related to mitochondrial function and neurodevelopment in females and circadian entrainment in males. Sex-specific transcriptomic neuroadaptations implicate unique neurobiological mechanisms underlying NOWS-like behaviors.

Spontaneously Hypertensive Rat substrains show differences in model traits for addiction risk and cocaine self-administration: Implications for a novel rat reduced complexity cross

Forward genetic mapping of F2 crosses between closely related substrains of inbred rodents - referred to as a reduced complexity cross (RCC) - is a relatively new strategy for accelerating the pace of gene discovery for complex traits, such as drug addiction. RCCs to date were generated in mice, but rats are thought to be optimal for addiction genetic studies. Based on past literature, one inbred Spontaneously Hypertensive Rat substrain, SHR/NCrl, is predicted to exhibit a distinct behavioral profile as it relates to cocaine self-administration traits relative to another substrain, SHR/NHsd. Direct substrain comparisons are a necessary first step before implementing an RCC. We evaluated model traits for cocaine addiction risk and cocaine self-administration behaviors using a longitudinal within-subjects design. Impulsive-like and compulsive-like traits were greater in SHR/NCrl than SHR/NHsd, as were reactivity to sucrose reward, sensitivity to acute psychostimulant effects of cocaine, and cocaine use studied under fixed-ratio and tandem schedules of cocaine self-administration. Compulsive-like behavior correlated with the acute psychostimulant effects of cocaine, which in turn correlated with cocaine taking under the tandem schedule. Compulsive-like behavior also was the best predictor of cocaine seeking responses. Heritability estimates indicated that 22 %-40 % of the variances for the above phenotypes can be explained by additive genetic factors, providing sufficient genetic variance to conduct genetic mapping in F2 crosses of SHR/NCrl and SHR/NHsd. These results provide compelling support for using an RCC approach in SHR substrains to uncover candidate genes and variants that are of relevance to cocaine use disorders.

Intracranial self-stimulation and concomitant behaviors following systemic methamphetamine administration in Hnrnph1 mutant mice

RATIONALE

Methamphetamine (MA) addiction is a major public health issue in the USA, with a poorly understood genetic component. We previously identified heterogeneous nuclear ribonucleoprotein H1 (Hnrnph1; H1) as a quantitative trait gene underlying sensitivity to MA-induced behavioral sensitivity. Mice heterozygous for a frameshift deletion in the first coding exon of H1 (H1^+/-) showed reduced MA phenotypes including oral self-administration, locomotor activity, dopamine release, and dose-dependent differences in MA conditioned place preference. However, the effects of H1^+/- on innate and MA-modulated reward sensitivity are not known.

OBJECTIVES

We examined innate reward sensitivity and facilitation by MA in H1^+/- mice via intracranial self-stimulation (ICSS).

METHODS

We used intracranial self-stimulation (ICSS) of the medial forebrain bundle to assess shifts in reward sensitivity following acute, ascending doses of MA (0.5-4.0 mg/kg, i.p.) using a within-subjects design. We also assessed video-recorded behaviors during ICSS testing sessions.

RESULTS

H1^+/- mice displayed reduced normalized maximum response rates in response to MA. H1^+/- females had lower normalized M50 values compared to wild-type females, suggesting enhanced reward facilitation by MA. Finally, regardless of genotype, there was a dose-dependent reduction in distance to the response wheel following MA administration, providing an additional measure of MA-induced reward-driven behavior.

CONCLUSIONS

H1^+/- mice displayed a complex ICSS phenotype following MA, displaying indications of both blunted reward magnitude (lower normalized maximum response rates) and enhanced reward sensitivity specific to H1^+/- females (lower normalized M50 values).

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

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

Hnrnph1 is a novel regulator of alcohol reward

BACKGROUND

Hnrnph1 is a validated quantitative trait gene for methamphetamine behavioral sensitivity that encodes for heterogeneous nuclear ribonucleoprotein H1 (hnRNP H1). This RNA-binding protein is involved in all stages of RNA metabolism that impacts mesocorticolimbic dopamine neurotransmission to influence addiction-related behavior.

METHODS

We characterized the alcohol behavioral phenotypes of mice heterozygous for a deletion in the first coding exon of Hnrnph1 (Hnrnph1+/-). We examined alcohol intake under both continuous- and limited-access procedures, as well as alcohol-induced place-conditioning. Follow-up studies examined genotypic differences in the psychomotor-activating and sedative-hypnotic effects of acute and repeated alcohol, and a behavioral test battery was employed to determine the effects of Hnrnph1 deletion on the manifestation of negative affect during alcohol withdrawal.

RESULTS

Relative to wild-type (WT) controls, Hnrnph1+/- males exhibited blunted intake of high alcohol concentrations under both drinking procedures. Hnrnph1 deletion did not impact the conditioned rewarding properties of low-dose alcohol, but reversed the conditioned place-aversion elicited by higher alcohol doses (2 and 4 g/kg), with more robust effects in male versus female mice. No genotypic differences were observed for alcohol-induced locomotor activity. Hnrnph1+/- mice exhibited a modest increase in sensitivity to alcohol's sedative-hypnotic effects, but did not differ from WT mice with regard to tolerance to alcohol's sedative-hypnotic effects or alcohol metabolism, Inconsistent effects of Hnrnph1 deletion were observed in models for withdrawal-induced negative affect.

CONCLUSIONS

These data identify Hnrnph1 as a novel, male-selective, driver of alcohol consumption and high-dose alcohol aversion that is potentially relevant to the neurobiology of alcohol abuse and alcoholism.

Confirmation of a Causal Taar1 Allelic Variant in Addiction-Relevant Methamphetamine Behaviors

Sensitivity to rewarding and reinforcing drug effects has a critical role in initial use, but the role of initial aversive drug effects has received less attention. Methamphetamine effects on dopamine re-uptake and efflux are associated with its addiction potential. However, methamphetamine also serves as a substrate for the trace amine-associated receptor 1 (TAAR1). Growing evidence in animal models indicates that increasing TAAR1 function reduces drug self-administration and intake. We previously determined that a non-synonymous single nucleotide polymorphism (SNP) in Taar1 predicts a conformational change in the receptor that has functional consequences. A Taar1 ^m1J mutant allele existing in DBA/2J mice expresses a non-functional receptor. In comparison to mice that possess one or more copies of the reference Taar1 allele (Taar1 ^+/+ or Taar1 ^+/m1J ), mice with the Taar1 ^m1J/m1J genotype readily consume methamphetamine, express low sensitivity to aversive effects of methamphetamine, and lack sensitivity to acute methamphetamine-induced hypothermia. We used three sets of knock-in and control mice in which one Taar1 allele was exchanged with the alternative allele to determine if other methamphetamine-related traits and an opioid trait are impacted by the same Taar1 SNP proven to affect MA consumption and hypothermia. First, we measured sensitivity to conditioned rewarding and aversive effects of methamphetamine to determine if an impact of the Taar1 SNP on these traits could be proven. Next, we used multiple genetic backgrounds to study the consistency of Taar1 allelic effects on methamphetamine intake and hypothermia. Finally, we studied morphine-induced hypothermia to confirm prior data suggesting that a gene in linkage disequilibrium with Taar1, rather than Taar1, accounts for prior observed differences in sensitivity. We found that a single SNP exchange reduced sensitivity to methamphetamine conditioned reward and increased sensitivity to conditioned aversion. Profound differences in methamphetamine intake and hypothermia consistently corresponded with genotype at the SNP location, with only slight variation in magnitude across genetic backgrounds. Morphine-induced hypothermia was not dependent on Taar1 genotype. Thus, Taar1 genotype and TAAR1 function impact multiple methamphetamine-related effects that likely predict the potential for methamphetamine use. These data support further investigation of their potential roles in risk for methamphetamine addiction and therapeutic development.

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.

Facilitating Complex Trait Analysis via Reduced Complexity Crosses

Genetically diverse inbred strains are frequently used in quantitative trait mapping to identify sequence variants underlying trait variation. Poor locus resolution and high genetic complexity impede variant discovery. As a solution, we explore reduced complexity crosses (RCCs) between phenotypically divergent, yet genetically similar, rodent substrains. RCCs accelerate functional variant discovery via decreasing the number of segregating variants by orders of magnitude. The simplified genetic architecture of RCCs often permit immediate identification of causal variants or rapid fine-mapping of broad loci to smaller intervals. Whole-genome sequences of substrains make RCCs possible by supporting the development of array- and targeted sequencing-based genotyping platforms, coupled with rapid genome editing for variant validation. In summary, RCCs enhance discovery-based genetics of complex traits.

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

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

The effect of the demyelinating agent cuprizone on binge-like eating of sweetened palatable food in female and male C57BL/6 substrains

Binge eating is a heritable symptom of eating disorders with an unknown genetic etiology. Rodent models for binge-like eating (BLE) of palatable food permit the study of genetic and biological mechanisms. We previously genetically mapped a coding mutation in Cyfip2 associated with increased BLE of sweetened palatable food in the C57BL/6NJ versus C57BL/6J substrain. The increase in BLE in C57BL/6NJ mice was associated with a decrease in transcription of genes enriched for myelination in the striatum. Here, we tested the hypothesis that decreasing myelin levels with the demyelinating agent cuprizone would enhance BLE. Mice were treated with a 0.3% cuprizone home cage diet for two weeks. Cuprizone induced similar weight loss in both substrains and sexes that recovered within 48 h after removal of cuprizone. Following a three-week recovery period, mice were trained for BLE in an intermittent, limited access procedure. Surprisingly, cuprizone significantly reduced BLE in male but not female C57BL/6NJ mice while having no effect in C57BL/6J mice. Cuprizone also reduced myelin basic protein (MBP) at seven weeks post-cuprizone removal while having no effect on myelin-associated glycoprotein at this time point. C57BL/6NJ mice also showed less MBP than C57BL/6J mice. There were no statistical interactions of Treatment with Sex on MBP levels, indicating that differences in MBP reduction are unlikely to account for sex differences in BLE. To summarize, cuprizone induced an unexpected, significant reduction in BLE in C57BL/6NJ males, which could indicate genotype-dependent sex differences in the biological mechanisms of BLE.

A Mutation in Hnrnph1 That Decreases Methamphetamine-Induced Reinforcement, Reward, and Dopamine Release and Increases Synaptosomal hnRNP H and Mitochondrial Proteins

Individual variation in the addiction liability of amphetamines has a heritable genetic component. We previously identified Hnrnph1 (heterogeneous nuclear ribonucleoprotein H1) as a quantitative trait gene underlying decreased methamphetamine-induced locomotor activity in mice. Here, we showed that mice (both females and males) with a heterozygous mutation in the first coding exon of Hnrnph1 (H1^+/-) showed reduced methamphetamine reinforcement and intake and dose-dependent changes in methamphetamine reward as measured via conditioned place preference. Furthermore, H1^+/- mice showed a robust decrease in methamphetamine-induced dopamine release in the NAc with no change in baseline extracellular dopamine, striatal whole-tissue dopamine, dopamine transporter protein, dopamine uptake, or striatal methamphetamine and amphetamine metabolite levels. Immunohistochemical and immunoblot staining of midbrain dopaminergic neurons and their forebrain projections for TH did not reveal any major changes in staining intensity, cell number, or forebrain puncta counts. Surprisingly, there was a twofold increase in hnRNP H protein in the striatal synaptosome of H1^+/- mice with no change in whole-tissue levels. To gain insight into the mechanisms linking increased synaptic hnRNP H with decreased methamphetamine-induced dopamine release and behaviors, synaptosomal proteomic analysis identified an increased baseline abundance of several mitochondrial complex I and V proteins that rapidly decreased at 30 min after methamphetamine administration in H1^+/- mice. In contrast, the much lower level of basal synaptosomal mitochondrial proteins in WT mice showed a rapid increase. We conclude that H1^+/- decreases methamphetamine-induced dopamine release, reward, and reinforcement and induces dynamic changes in basal and methamphetamine-induced synaptic mitochondrial function.SIGNIFICANCE STATEMENT Methamphetamine dependence is a significant public health concern with no FDA-approved treatment. We discovered a role for the RNA binding protein hnRNP H in methamphetamine reward and reinforcement. Hnrnph1 mutation also blunted methamphetamine-induced dopamine release in the NAc, a key neurochemical event contributing to methamphetamine addiction liability. Finally, Hnrnph1 mutants showed a marked increase in basal level of synaptosomal hnRNP H and mitochondrial proteins that decreased in response to methamphetamine, whereas WT mice showed a methamphetamine-induced increase in synaptosomal mitochondrial proteins. Thus, we identified a potential role for hnRNP H in basal and dynamic mitochondrial function that informs methamphetamine-induced cellular adaptations associated with reduced addiction liability.

Cyfip1 Haploinsufficiency Increases Compulsive-Like Behavior and Modulates Palatable Food Intake in Mice: Dependence on Cyfip2 Genetic Background, Parent-of Origin, and Sex

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 (Cyfip2N/N) background and the BE-resistant C57BL/6J (Cyfip2J/J) background. Cyfip1+/- mice showed increased compulsive-like behavior on both backgrounds and increased PF intake on the Cyfip2N/N background. In contrast, maternal Cyfip1 haploinsufficiency on the BE-resistant Cyfip2J/J background induced a robust escalation in PF intake in wild-type Cyfip1J/J males while having no effect in Cyfip1J/- 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.

Regional Analysis of the Brain Transcriptome in Mice Bred for High and Low Methamphetamine Consumption

Transcriptome profiling can broadly characterize drug effects and risk for addiction in the absence of drug exposure. Modern large-scale molecular methods, including RNA-sequencing (RNA-Seq), have been extensively applied to alcohol-related disease traits, but rarely to risk for methamphetamine (MA) addiction. We used RNA-Seq data from selectively bred mice with high or low risk for voluntary MA intake to construct coexpression and cosplicing networks for differential risk. Three brain reward circuitry regions were explored, the nucleus accumbens (NAc), prefrontal cortex (PFC), and ventral midbrain (VMB). With respect to differential gene expression and wiring, the VMB was more strongly affected than either the PFC or NAc. Coexpression network connectivity was higher in the low MA drinking line than in the high MA drinking line in the VMB, oppositely affected in the NAc, and little impacted in the PFC. Gene modules protected from the effects of selection may help to eliminate certain mechanisms from significant involvement in risk for MA intake. One such module was enriched in genes with dopamine-associated annotations. Overall, the data suggest that mitochondrial function and glutamate-mediated synaptic plasticity have key roles in the outcomes of selective breeding for high versus low levels of MA intake.

miR-301a promotes lung tumorigenesis by suppressing Runx3

BACKGROUND

Our previous report demonstrated that genetic ablation of miR-301a reduces Kras-driven lung tumorigenesis in mice. However, the impact of miR-301a on host anti-tumor immunity remains unexplored. Here we assessed the underlying molecular mechanisms of miR-301a in the tumor microenvironment.

METHODS

The differentially expressed genes were identified by using deep sequencing. The immune cell counts, and cytokines expression were analyzed by realtime PCR, immunohistochemistry and flow cytometry. The role of miR-301a/Runx3 in lung tumor was evaluated on cell growth, migration and invasion. The function of miR-301a/Runx3 in regulating tumor microenvironment and tumor metastasis were evaluated in Kras transgenic mice and B16/LLC1 syngeneic xenografts tumor models.

RESULTS

In this work, we identified 1166 up-regulated and 475 down-regulated differentially expressed genes in lung tumor tissues between KrasLA2 and miR-301a-/-; KrasLA2 mice. Immune response and cell cycle were major pathways involved in the protective role of miR-301a deletion in lung tumorigenesis. Overexpression of the miR-301a target, Runx3, was an early event identified in miR-301a-/-; KrasLA2 mice compared to WT-KrasLA2 mice. We found that miR-301a deletion enhanced CD8+ T cell accumulation and IFN-γ production in the tumor microenvironment and mediated antitumor immunity. Further studies revealed that miR-301a deficiency in the tumor microenvironment effectively reduced tumor metastasis by elevating Runx3 and recruiting CD8+ T cells, whereas miR-301a knockdown in tumor cells themselves restrained cell migration by elevating Runx3 expression.

CONCLUSIONS

Our findings further underscore that miR-301a facilitates tumor microenvironment antitumor immunity by Runx3 suppression in lung tumorigenesis.

Single-Cell RNA Sequencing of Human Embryonic Stem Cell Differentiation Delineates Adverse Effects of Nicotine on Embryonic Development

Nicotine, the main chemical constituent of tobacco, is highly detrimental to the developing fetus by increasing the risk of gestational complications and organ disorders. The effects of nicotine on human embryonic development and related mechanisms, however, remain poorly understood. Here, we performed single-cell RNA sequencing (scRNA-seq) of human embryonic stem cell (hESC)-derived embryoid body (EB) in the presence or absence of nicotine. Nicotine-induced lineage-specific responses and dysregulated cell-to-cell communication in EBs, shedding light on the adverse effects of nicotine on human embryonic development. In addition, nicotine reduced cell viability, increased reactive oxygen species (ROS), and altered cell cycling in EBs. Abnormal Ca2+ signaling was found in muscle cells upon nicotine exposure, as verified in hESC-derived cardiomyocytes. Consequently, our scRNA-seq data suggest direct adverse effects of nicotine on hESC differentiation at the single-cell level and offer a new method for evaluating drug and environmental toxicity on human embryonic development in utero.

C57BL/6 substrain differences in inflammatory and neuropathic nociception and genetic mapping of a major quantitative trait locus underlying acute thermal nociception

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.

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

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

Changes in neuronal immunofluorescence in the C- versus N-terminal domains of hnRNP H following D1 dopamine receptor activation

RNA binding proteins are a diverse class of proteins that regulate all aspects of RNA metabolism. Accumulating studies indicate that heterogeneous nuclear ribonucleoproteins are associated with cellular adaptations in response to drugs of abuse. We recently mapped and validated heterogeneous nuclear ribonucleoprotein H1 (Hnrnph1) as a quantitative trait gene underlying differential behavioral sensitivity to methamphetamine. The molecular mechanisms by which hnRNP H1 alters methamphetamine behaviors are unknown but could involve pre- and/or post-synaptic changes in protein localization and function. Methamphetamine initiates post-synaptic D1 dopamine receptor signaling indirectly by binding to pre-synaptic dopamine transporters and vesicular monoamine transporters of midbrain dopaminergic neurons which triggers reverse transport and accumulation of dopamine at the synapse. Here, we examined changes in neuronal localization of hnRNP H in primary rat cortical neurons that express dopamine receptors that can be modulated by the D1 or D2 dopamine receptor agonists SKF38393 and (-)-Quinpirole HCl, respectively. Basal immunostaining of hnRNP H was localized primarily to the nucleus. D1 dopamine receptor activation induced an increase in hnRNP H nuclear immunostaining as detected by immunocytochemistry with a C-domain directed antibody containing epitope near the glycine-rich domain but not with an N-domain specific antibody. Although there was no change in hnRNP H protein in the nucleus or cytoplasm, there was a decrease in Hnrnph1 transcript following D1 receptor stimulation. Taken together, these results suggest that D1 receptor activation increases availability of the hnRNP H C-terminal epitope, which could potentially reflect changes in protein-protein interactions. Thus, D1 receptor signaling could represent a key molecular post-synaptic event linking Hnrnph1 polymorphisms to drug-induced behavior.

Intravenous cocaine self-administration in a panel of inbred mouse strains differing in acute locomotor sensitivity to cocaine

RATIONALE

Initial sensitivity to drugs of abuse often predicts subsequent use and abuse, but this relationship is not always observed in human studies. Moreover, studies examining the relationship between initial locomotor sensitivity and the rewarding and reinforcing effects of drugs in animal models have also been equivocal. Understanding the relationship between initial drug effects and propensity to continue use, potentially resulting in the development of a substance use disorder, may help to identify key targets for prevention and treatment.

OBJECTIVES

We examined intravenous cocaine self-administration in a set of mouse strains that were previously identified to be at the phenotypic extremes for cocaine-induced locomotor activation to determine if initial locomotor sensitivity predicted acquisition, extinction, dose response, or progressive ratio (PR) breakpoint.

METHODS

We selected eight inbred mouse strains based on locomotor sensitivity to 20 mg/kg cocaine. These strains, designated as low and high responders, were tested in an intravenous self-administration paradigm that included acquisition of 0.5 mg/(kg*inf) under a FR1 schedule, extinction, re-acquisition, dose response to 0.125, 0.25, 0.5, 1, and 2 mg/(kg*inf), and progressive ratio.

RESULTS

We observed overall differences in self-administration behavior between high and low responders. Low responders self-administered less cocaine and had lower breakpoints under the PR schedule. However, we also observed strain differences within each group. Self-administration in the low responder, LG/J, more closely resembled the behavior of the high-responding group, and the high responder, P/J, had self-administration behavior that more closely resembled the low-responding group.

CONCLUSIONS

We conclude that acute cocaine-induced locomotor activation does predict self-administration behavior, but in a strain-specific manner. These data support the idea that genetic background influences the relationship among addiction-related behaviors.

A Spontaneous Mutation in Taar1 Impacts Methamphetamine-Related Traits Exclusively in DBA/2 Mice from a Single Vendor

Major gene effects on traits associated with substance use disorders are rare. Previous findings in methamphetamine drinking (MADR) lines of mice, bred for high or low voluntary MA intake, and in null mutants demonstrate a major impact of the trace amine-associated receptor 1 (Taar1) gene on a triad of MA-related traits: MA consumption, MA-induced conditioned taste aversion and MA-induced hypothermia. While inbred strains are fundamentally genetically stable, rare spontaneous mutations can become fixed and result in new or aberrant phenotypes. A single nucleotide polymorphism in Taar1 that encodes a missense proline to threonine mutation in the second transmembrane domain (Taar1^m1J) has been identified in the DBA/2J strain. MA is an agonist at this receptor, but the receptor produced by Taar1^m1J does not respond to MA or endogenous ligands. In the present study, we used progeny of the C57BL/6J × DBA/2J F2 cross, the MADR lines, C57BL/6J × DBA/2J recombinant inbred strains, and DBA/2 mice sourced from four vendors to further examine Taar1-MA phenotype relations and to define the chronology of the fixation of the Taar1^m1J mutation. Mice homozygous for Taar1^m1J were found at high frequency early in selection for high MA intake in multiple replicates of the high MADR line, whereas Taar1^m1J homozygotes were absent in the low MADR line. The homozygous Taar1^m1J genotype is causally linked to increased MA intake, reduced MA-induced conditioned taste aversion, and reduced MA-induced hypothermia across models. Genotype-phenotype correlations range from 0.68 to 0.96. This Taar1 polymorphism exists in DBA/2J mice sourced directly from The Jackson Laboratory, but not DBA/2 mice sourced from Charles River (DBA/2NCrl), Envigo (formerly Harlan Sprague Dawley; DBA/2NHsd) or Taconic (DBA/2NTac). By genotyping archived samples from The Jackson Laboratory, we have determined that this mutation arose in 2001–2003. Our data strengthen the conclusion that the mutant Taar1^m1J allele, which codes for a non-functional receptor protein, increases risk for multiple MA-related traits, including MA intake, in homozygous Taar1^m1J individuals.

Casein kinase 1-epsilon deletion increases mu opioid receptor-dependent behaviors and binge eating1

Genetic and pharmacological studies indicate that casein kinase 1 epsilon (Csnk1e) contributes to psychostimulant, opioid, and ethanol motivated behaviors. We previously used pharmacological inhibition to demonstrate that Csnk1e negatively regulates the locomotor stimulant properties of opioids and psychostimulants. Here, we tested the hypothesis that Csnk1e negatively regulates opioid and psychostimulant reward using genetic inhibition and the conditioned place preference assay in Csnk1e knockout mice. Similar to pharmacological inhibition, Csnk1e knockout mice showed enhanced opioid-induced locomotor activity with the mu opioid receptor agonist fentanyl (0.2 mg/kg i.p.) as well as enhanced sensitivity to low-dose fentanyl reward (0.05 mg/kg). Interestingly, female knockout mice also showed a markedly greater escalation in consumption of sweetened palatable food - a behavioral pattern consistent with binge eating that also depends on mu opioid receptor activation. No difference was observed in fentanyl analgesia in the 52.5°C hot plate assay (0-0.4 mg/kg), naloxone conditioned place aversion (4 mg/kg), or methamphetamine conditioned place preference (0-4 mg/kg). To identify molecular adaptations associated with increased drug and food behaviors in knockout mice, we completed transcriptome analysis via mRNA sequencing of the striatum. Enrichment analysis identified terms associated with myelination and axon guidance and pathway analysis identified a differentially expressed gene set predicted to be regulated by the Wnt signaling transcription factor, Tcf7l2. To summarize, Csnk1e deletion increased mu opioid receptor-dependent behaviors, supporting previous studies indicating an endogenous negative regulatory role of Csnk1e in opioid behavior.

Methamphetamine Addiction Vulnerability: The Glutamate, the Bad, and the Ugly

BACKGROUND

The high prevalence and severity of methamphetamine (MA) abuse demands greater neurobiological understanding of its etiology.

METHODS

We conducted immunoblotting and in vivo microdialysis procedures in MA high/low drinking mice, as well as in isogenic C57BL/6J mice that varied in their MA preference/taking, to examine the glutamate underpinnings of MA abuse vulnerability. Neuropharmacological and Homer2 knockdown approaches were also used in C57BL/6J mice to confirm the role for nucleus accumbens (NAC) glutamate/Homer2 expression in MA preference/aversion.

RESULTS

We identified a hyperglutamatergic state within the NAC as a biochemical trait corresponding with both genetic and idiopathic vulnerability for high MA preference and taking. We also confirmed that subchronic subtoxic MA experience elicits a hyperglutamatergic state within the NAC during protracted withdrawal, characterized by elevated metabotropic glutamate 1/5 receptor function and Homer2 receptor-scaffolding protein expression. A high MA-preferring phenotype was recapitulated by elevating endogenous glutamate within the NAC shell of mice and we reversed MA preference/taking by lowering endogenous glutamate and/or Homer2 expression within this subregion.

CONCLUSIONS

Our data point to an idiopathic, genetic, or drug-induced hyperglutamatergic state within the NAC as a mediator of MA addiction vulnerability.

Cytoplasmic FMR1-Interacting Protein 2 Is a Major Genetic Factor Underlying Binge Eating

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 F₂ cross to identify quantitative trait loci associated with binge eating. We used gene targeting to validate candidate genetic factors. Finally, we used transcriptome analysis of the striatum via messenger RNA sequencing 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 quantitative trait locus on chromosome 11 (logarithm of the odds = 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.

Identification of a Putative Quantitative Trait Gene for Resistance to Obesity in Mice Using Transcriptome Analysis and Causal Inference Tests

It is still challenging to identify causal genes governing obesity. Pbwg1.5, a quantitative trait locus (QTL) for resistance to obesity, was previously discovered from wild Mus musculus castaneus mice and was fine-mapped to a 2.1-Mb genomic region of mouse chromosome 2, where no known gene with an effect on white adipose tissue (WAT) has been reported. The aim of this study was to identify a strong candidate gene for Pbwg1.5 by an integration approach of transcriptome analysis (RNA-sequencing followed by real-time PCR analysis) and the causal inference test (CIT), a statistical method to infer causal relationships between diplotypes, gene expression and trait values. Body weight, body composition and biochemical traits were measured in F2 mice obtained from an intercross between the C57BL/6JJcl strain and a congenic strain carrying Pbwg1.5 on the C57BL/6JJcl background. The F2 mice showed significant diplotype differences in 12 traits including body weight, WAT weight and serum cholesterol/triglyceride levels. The transcriptome analysis revealed that Ly75, Pla2r1, Fap and Gca genes were differentially expressed in the liver and that Fap, Ifih1 and Grb14 were differentially expressed in WAT. However, CITs indicated statistical evidence that only the liver Ly75 gene mediated between genotype and WAT. Ly75 expression was negatively associated with WAT weight. The results suggested that Ly75 is a putative quantitative trait gene for the obesity-resistant Pbwg1.5 QTL discovered from the wild M. m. castaneus mouse. The finding provides a novel insight into a better understanding of the genetic basis for prevention of obesity.

RNA-binding proteins, neural development and the addictions

Transcriptional and post-transcriptional regulation of gene expression defines the neurobiological mechanisms that bridge genetic and environmental risk factors with neurobehavioral dysfunction underlying the addictions. More than 1000 genes in the eukaryotic genome code for multifunctional RNA-binding proteins (RBPs) that can regulate all levels of RNA biogenesis. More than 50% of these RBPs are expressed in the brain where they regulate alternative splicing, transport, localization, stability and translation of RNAs during development and adulthood. Dysfunction of RBPs can exert global effects on their targetomes that underlie neurodegenerative disorders such as Alzheimer's and Parkinson's diseases as well as neurodevelopmental disorders, including autism and schizophrenia. Here, we consider the evidence that RBPs influence key molecular targets, neurodevelopment, synaptic plasticity and neurobehavioral dysfunction underlying the addictions. Increasingly well-powered genome-wide association studies in humans and mammalian model organisms combined with ever more precise transcriptomic and proteomic approaches will continue to uncover novel and possibly selective roles for RBPs in the addictions. Key challenges include identifying the biological functions of the dynamic RBP targetomes from specific cell types throughout subcellular space (e.g. the nuclear spliceome vs. the synaptic translatome) and time and manipulating RBP programs through post-transcriptional modifications to prevent or reverse aberrant neurodevelopment and plasticity underlying the addictions.