Search

A major barrier to remission from an alcohol use disorder (AUD) is the continued risk of relapse during abstinence. Assessing the neuroadaptations after chronic alcohol and repeated abstinence is important to identify mechanisms that may contribute to relapse. In this study, we used a rhesus macaque model of long-term alcohol use and repeated abstinence, providing a platform to extend mechanistic findings from rodents to primates. The central amygdala (CeA) displays elevated GABA release following chronic alcohol in rodents and in abstinent male macaques, highlighting this neuroadaptation as a conserved mechanism that may underlie excessive alcohol consumption. Here, we determined circulating interleukin-1β (IL-1β) levels, CeA transcriptomic changes, and the effects of IL-1β and corticotropin releasing factor (CRF) signaling on CeA GABA transmission in male controls and abstinent drinkers. While no significant differences in peripheral IL-1β or the CeA transcriptome were observed, pathway analysis identified several canonical immune-related pathways. We addressed this potential dysregulation of CeA immune signaling in abstient drinkers with an electrophysiological approach. We found that IL-1β decreased CeA GABA release in controls while abstinent drinkers were less sensitive to IL-1β’s effects, suggesting adaptations in the neuromodulatory role of IL-1β. In contrast, CRF enhanced CeA GABA release similarly in controls and abstinent drinkers, consistent with rodent studies. Notably, CeA CRF expression was inversely correlated with intoxication, suggesting that CRF levels during abstinence may predict future intoxication. Together, our findings highlight conserved and divergent actions of chronic alcohol on neuroimmune and stress signaling on CeA GABA transmission across rodents and macaques.

Circuit manipulation has been a staple technique in neuroscience to identify how the brain functions to control complex behaviors. Chemogenetics, including designer receptors exclusively activated by designer drugs (DREADDs), have proven to be a powerful tool for the reversible modulation of discrete brain circuitry without the need for implantable devices, thereby making them especially useful in awake and unrestrained animals. This study used a DREADD approach to query the role of the nucleus accumbens (NAc) in mediating the interoceptive effects of 1.0 g/kg ethanol (i.g.) in rhesus monkeys (n = 7) using a drug discrimination procedure. After training, stereotaxic surgery was performed to introduce an AAV carrying the human muscarinic 4 receptor DREADD (hM4Di) bilaterally into the NAc. The hypothesis was that decreasing the output of the NAc by activation of hM4Di with the DREADD actuator, clozapine-n-oxide (CNO), would potentiate the discriminative stimulus effect of ethanol (i.e., a leftward shift the ethanol dose discrimination curve). The results showed individual variability shifts of the ethanol dose-response determination under DREADD activation. Characterization of the expression and function of hM4Di with MRI, immunohistochemical, and electrophysiological techniques found the selectivity of NAc transduction was proportional to behavioral effect. Specifically, the proportion of hM4Di expression restricted to the NAc was associated with the potency of the discriminative stimulus effects of ethanol. Together, these experiments highlight the NAc in mediating the interoceptive effects of ethanol, provide a framework for validation of chemogenetic tools in primates, and underscore the importance of robust within-subjects examination of DREADD expression for interpretation of behavioral findings.

BACKGROUND Phosphatidylethanol (PEth) homologs are ethanol metabolites used to identify and monitor alcohol drinking in humans. In this study, we measured levels of the 2 most abundant homologs, PEth 16:0/18:1 and PEth 16:0/18:2, in whole blood samples from rhesus macaque monkeys that drank ethanol daily ad libitum to assess the relationship between PEth levels and recent ethanol exposure in this animal model. METHODS Blood samples were obtained from The Monkey Alcohol Tissue Research Resource. The monkeys were first induced to consume 4% (w/v) ethanol in water from a panel attached to their home cage. Then, monkeys were allowed to drink ethanol and water ad libitum 22 h daily for 12 months and the daily amount of ethanol each monkey consumed was measured. Whole, uncoagulated blood was collected from each animal at the end of the entire experimental procedure. PEth 16:0/18:1 and PEth 16:0/18:2 levels were analyzed by HPLC with tandem mass spectrometry, and the ethanol consumed during the preceding 14 days was measured. Combined PEth was the sum of the concentrations of both homologs. RESULTS Our results show that (1) PEth accumulates in the blood of rhesus monkeys after ethanol consumption; (2) PEth homolog levels were correlated with the daily average ethanol intake during the 14-day period immediately preceding blood collection; (3) the application of established human PEth 16:0/18:1 cutoff concentrations indicative of light social or no ethanol consumption (

Chronic heavy alcohol use is often associated with reduced bone mineral density and altered bone turnover. However, the dose response effects of ethanol on bone turnover have not been established. This study examined the effects of graded increases of ethanol consumption on biochemical markers of bone turnover in young adult male cynomolgus macaques (Macaca fascicularis). For this study, 6.6-year-old (95% CI: 6.5, 6.7) male macaques were subjected to three 30-day sessions of increased ethanol intake over a 90-day interval. During the first 30 days, the monkeys drank a predetermined volume of ethanol corresponding to 0.5 g/kg/day, followed by 1.0 g/kg/day and 1.5 g/kg/day. Osteocalcin, a marker of bone formation, and carboxyterminal cross-linking telopeptide of type 1 collagen (CTX), a marker of resorption, were measured during each 30-day session. In addition, the ratio of osteocalcin to CTX was determined as a surrogate measure of global turnover balance. Mean osteocalcin decreased by 2.6 ng/mL (1.8, 3.5) for each one-half unit (0.5 g/kg/day) increase in dose (p 0.001). Mean CTX decreased by 0.13 ng/mL (0.06, 0.20) for each one-half unit increase in dose (p 0.001). Furthermore, there was an inverse relationship between dose and the ratio of osteocalcin to CTX, such that the mean ratio decreased by 0.9 (0.3, 1.5) for each one-half unit increase in dose (p = 0.01). In summary, male cynomolgus macaques had decreased blood osteocalcin and CTX, and osteocalcin to CTX ratio during the 90-day interval of graded increases in ethanol consumption, indicative of reduced bone turnover and negative turnover balance, respectively. These findings suggest that over the range ingested, ethanol resulted in a linear decrease in bone turnover. Furthermore, the negative bone turnover balance observed is consistent with reported effects of chronic alcohol intake on the skeleton.

It is well established that chronic heavy alcohol drinking (CHD) results in significant organ damage, increased susceptibility to infections, and poor outcomes following injury. In contrast, chronic moderate drinking (CMD) has been associated with improved cardiovascular health and immunity. These differential outcomes have been linked to alterations in both innate and adaptive branches of the immune system; however, the mechanisms remain poorly understood. To address this question, we determined the impact of chronic drinking on the transcriptional and functional responses of peripheral blood mononuclear cells (PBMC) collected from male rhesus macaques classified as CMD or CHD after 12 months of voluntary ethanol self-administration. Our analysis suggests that chronic alcohol drinking, regardless of dose alters resting transcriptomes of PBMC, with the largest impact seen in innate immune cells. These transcriptional changes are partially explained by alterations in microRNA profiles. Additionally, chronic alcohol drinking is associated with a dose dependent heightened inflammatory profiled at resting and following LPS stimulation. Moreover, we observed a dose-dependent shift in the kinetics of transcriptional responses to LPS. These findings may explain the dichotomy in clinical and immunological outcomes observed with moderate versus heavy alcohol drinking.

The nucleus accumbens core (NAcc) has been repeatedly demonstrated to be a key component of the circuitry associated with excessive ethanol consumption. Previous studies have illustrated that in a nonhuman primate (NHP) model of chronic ethanol consumption, there is significant epigenetic remodeling of the NAcc. In the current study, RNA-Seq was used to examine genome-wide gene expression in eight each of control, low/binge (LD*), and high/very high (HD*) rhesus macaque drinkers. Using an FDR < 0.05, zero genes were significantly differentially expressed (DE) between LD* and controls, six genes between HD* and LD*, and 734 genes between HD* and controls. Focusing on HD* versus control DE genes, the upregulated genes (N = 366) were enriched in genes with annotations associated with signal recognition particle (SRP)-dependent co-translational protein targeting to membrane (FDR < 3 × 10-59 ), structural constituent of ribosome (FDR < 3 × 10-47 ), and ribosomal subunit (FDR < 5 × 10-48 ). Downregulated genes (N = 363) were enriched in annotations associated with behavior (FDR < 2 × 10-4 ), membrane organization (FDR < 1 × 10-4 ), inorganic cation transmembrane transporter activity (FDR < 2 × 10-3 ), synapse part (FDR < 4 × 10-10 ), glutamatergic synapse (FDR < 1 × 10-6 ), and GABAergic synapse (FDR < 6 × 10-4 ). Ingenuity Pathway Analysis (IPA) revealed that EIF2 signaling and mTOR pathways were significantly upregulated in HD* animals (FDR < 3 × 10-33 and

The efficacy of short-term treatment with mifepristone (MIFE), a high-affinity, nonselective glucocorticoid receptor antagonist, to reduce ethanol drinking was tested in a rhesus macaque model. Stable individual daily ethanol intakes were established, ranging from 1.6 to 4.0 g/kg per day (n = 9 monkeys). After establishment of chronic ethanol intake, a MIFE dosing regimen that modeled a study of rodent drinking and human alcohol craving was evaluated. Three doses of MIFE (17, 30, and 56 mg/kg per day) were each administered for four consecutive days. Both 30 and 56 mg/kg decreased ethanol intake compared with baseline drinking levels without a change in water intake. The dose of 56 mg/kg per day of MIFE produced the largest reduction in ethanol self-administration, with the average intake at 57% of baseline intakes. Cortisol was elevated during MIFE dosing, and a mediation analysis revealed that the effect on ethanol drinking was fully mediated through cortisol. During a forced abstinence phase, access to 1.5 g/kg ethanol resulted in relapse in all drinkers and was not altered by treatment with 56 mg/kg MIFE. Overall, these results show that during active drinking MIFE is efficacious in reducing heavy alcohol intake in a monkey model, an effect that was related to MIFE-induced increase in cortisol. However, MIFE treatment did not eliminate ethanol drinking. Further, cessation of MIFE treatment resulted in a rapid return to baseline intakes, and MIFE was not effective in preventing a relapse during early abstinence. SIGNIFICANCE STATEMENT: Mifepristone reliably decreases average daily ethanol self-administration in a nonhuman primate model. This effect was mediated by cortisol, was most effective during open-access conditions, and did not prevent or reduce relapse drinking.

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.

This is the first description of the relationship between chronic ethanol self-administration and the brain transcriptome in a non-human primate (rhesus macaque). Thirty-one male animals self-administered ethanol on a daily basis for over 12 months. Gene transcription was quantified with RNA-Seq in the central nucleus of the amygdala (CeA) and cortical Area 32. We constructed coexpression and cosplicing networks, and we identified areas of preservation and areas of differentiation between regions and network types. Correlations between intake and transcription included largely distinct gene sets and annotation categories across brain regions and between expression and splicing; positive and negative correlations were also associated with distinct annotation groups. Membrane, synaptic and splicing annotation categories were over-represented in the modules (gene clusters) enriched in positive correlations (CeA); our cosplicing analysis further identified the genes affected only at the exon inclusion level. In the CeA coexpression network, we identified Rab6b, Cdk18 and Igsf21 among the intake-correlated hubs, while in the Area 32, we identified a distinct hub set that included Ppp3r1 and Myeov2. Overall, the data illustrate that excessive ethanol self-administration is associated with broad expression and splicing mechanisms that involve membrane and synapse genes.

BACKGROUND: The Monkey Alcohol and Tissue Research Resource (MATRR) is a repository and analytics platform for detailed data derived from well-documented Non-Human Primate (NHP) alcohol self-administration studies. This macaque model has demonstrated categorical drinking norms reflective of human drinking populations, resulting in consumption pattern classifications of Very Heavy Drinking (VHD), Heavy Drinking (HD), Binge Drinking (BD), and Low Drinking (LD) individuals. Here we expand on previous findings that suggest ethanol drinking patterns during initial drinking to intoxication can reliably predict future drinking category assignment. METHODS: The classification strategy uses a machine learning approach to examine an extensive set of daily drinking attributes during 90 sessions of induction across seven cohorts of five to eight monkeys for a total of 50 animals. A Random Forest classifier is employed to accurately predict categorical drinking after 12 months of self-administration. RESULTS: Predictive outcome accuracy is approximately 78% when classes are aggregated into two groups, "LD and BD" and "HD and VHD". A subsequent two-step classification model distinguishes individual LD and BD categories with 90% accuracy and between HD and VHD categories with 95% accuracy. Average four-category classification accuracy is 74%, and provides putative distinguishing behavioral characteristics between groupings. CONCLUSION: We demonstrate that data derived from the induction phase of this ethanol self-administration protocol has significant predictive power for future ethanol consumption patterns. Importantly, numerous predictive factors are longitudinal, measuring the change of drinking patterns through three stages of induction. Factors during induction that predict future heavy drinkers include being younger at the time of first intoxication and developing a shorter latency to first ethanol drink. Overall, this analysis identifies predictive characteristics in future very heavy drinkers that optimize intoxication, such as having increasingly fewer bouts with more drinks. This analysis also identifies characteristic avoidance of intoxicating topographies in future low drinkers, such as increasing number of bouts and waiting longer before the first ethanol drink. This article is protected by copyright. All rights reserved.This article is protected by copyright. All rights reserved. Language: en

Background Chronic heavy alcohol consumption is an established risk factor for bone fracture, but comorbidities associated with alcohol intake may contribute to increased fracture rates in alcohol abusers. To address the specific effects of alcohol on bone, we used a nonhuman primate model and evaluated voluntary alcohol consumption on: (i) global markers of bone turnover in blood and (ii) cancellous bone mass, density, microarchitecture, turnover, and microdamage in lumbar vertebra. Methods Following a 4-month induction period, 6-year-old male rhesus macaques (Macaca mulatta, n = 13) voluntarily self-administered water or ethanol (EtOH; 4% w/v) for 22 h/d, 7 d/wk, for a total of 12 months. Control animals (n = 9) consumed an isocaloric maltose-dextrin solution. Tetracycline hydrochloride was administered orally 17 and 3 days prior to sacrifice to label mineralizing bone surfaces. Global skeletal response to EtOH was evaluated by measuring plasma osteocalcin and carboxyterminal collagen cross-links (CTX). Local response was evaluated in lumbar vertebra using dual-energy X-ray absorptiometry, microcomputed tomography, static and dynamic histomorphometry, and histological assessment of microdamage. Results Monkeys in the EtOH group consumed an average of 2.8 ± 0.2 (mean ± SE) g/kg/d of EtOH (30 ± 2% of total calories), resulting in an average blood EtOH concentration of 88.3 ± 8.8 mg/dl 7 hours after the session onset. Plasma CTX and osteocalcin tended to be lower in EtOH-consuming monkeys compared to controls. Significant differences in bone mineral density in lumbar vertebrae 1 to 4 were not detected with treatment. However, cancellous bone volume fraction (in cores biopsied from the central region of the third vertebral body) was lower in EtOH-consuming monkeys compared to controls. Furthermore, EtOH-consuming monkeys had lower osteoblast perimeter and mineralizing perimeter, no significant difference in osteoclast perimeter, and higher bone marrow adiposity than controls. No significant differences between groups were detected in microcrack density (2nd lumbar vertebra). Conclusions Voluntary chronic heavy EtOH consumption reduces cancellous bone formation in lumbar vertebra by decreasing osteoblast-lined bone perimeter, a response associated with an increase in bone marrow adiposity.

We identified a locus on mouse chromosome 10 that accounts for 60% of the genetic variance in methamphetamine intake in mice selectively bred for high versus low methamphetamine consumption. We nominated the trace amine-associated receptor 1 gene, Taar1, as the strongest candidate and identified regulation of the mu-opioid receptor 1 gene, Oprm1, as another contributor. This study exploited CRISPR-Cas9 to test the causal role of Taar1 in methamphetamine intake and a genetically-associated thermal response to methamphetamine. The methamphetamine-related traits were rescued, converting them to levels found in methamphetamine-avoiding animals. We used a family of recombinant inbred mouse strains for interval mapping and to examine independent and epistatic effects of Taar1 and Oprm1. Both methamphetamine intake and the thermal response mapped to Taar1 and the independent effect of Taar1 was dependent on genotype at Oprm1. Our findings encourage investigation of the contribution of Taar1 and Oprm1 variants to human methamphetamine addiction., eLife digest People who misuse drugs often do so partly in response to the environment they find themselves in, and partly because of their genetics. The genetic component of someone’s risk is influenced by many different genes, and most research has found that each gene has a small individual effect. A method called quantitative trait locus (QTL) analysis can help find parts of the genome that influence someone’s risk of misusing drugs. In 2013, researchers found one region on chromosome 10 in mice has a particularly large influence on how much methamphetamine an individual mouse will ingest if the drug was available in one of its two water bottles. A gene called Taar1 was particularly important in this region and another gene, called Oprm1, may also play a significant role. When the Taar1 gene is switched off, mice consume larger amounts of methamphetamine, have a heightened reward response from the drug, and are insensitive to the adverse effects – such as hypothermia. But whether Taar1 directly caused these effects, and whether Taar1 and Oprm1 interact, had not yet been determined. If these genes played a causal role, they could be useful targets for treatment of methamphetamine-use disorder. Stafford, Reed et al. – who include several of the researchers involved in the 2013 work – now report that when a particular variant of Taar1 was present in mice they consumed large amounts of methamphetamine. The variant codes for a faulty version of a receptor protein. When this variant was replaced with a working version using gene editing, the mice consumed less methamphetamine and also became sensitive to hypothermia induced by the drug. This confirms that this gene does play a causal role in methamphetamine consumption and hypothermia. Next, Stafford, Reed et al. tested mice with different combinations of variants of Oprm1 and Taar1 to see how the genes interacted. The results showed that the effects of Taar1 on both consumption of the drug and hypothermia depended on the Oprm1 variant present. The findings suggest that variants of these two genes in humans could influence an individual’s risk of addiction to methamphetamine. It is possible that in future the disorder could be treated by drugs that modify the brain activity impacted by these receptors. But first, it will be important to find out if these genes play a similar role in humans as they do in mice.

Background Genome‐wide profiling to examine brain transcriptional features associated with excessive ethanol (EtOH) consumption has been applied to a variety of species including rodents, nonhuman primates (NHPs), and humans. However, these data were obtained from cross‐sectional samples which are particularly vulnerable to individual variation when obtained from small outbred populations typical of human and NHP studies. In the current study, a novel within‐subject design was used to examine the effects of voluntary EtOH consumption on prefrontal cortex (PFC) gene expression in a NHP model. Methods Two cohorts of cynomolgus macaques (n = 23) underwent a schedule‐induced polydipsia procedure to establish EtOH self‐administration followed by 6 months of daily open access to EtOH (4% w/v) and water. Individual daily EtOH intakes ranged from an average of 0.7 to 3.7 g/kg/d. Dorsal lateral PFC area 46 (A46) brain biopsies were collected in EtOH‐naïve and control monkeys; contralateral A46 biopsies were collected from the same monkeys following the 6 months of fluid consumption. Gene expression changes were assessed using RNA‐Seq paired analysis, which allowed for correction of individual baseline differences in gene expression. Results A total of 675 genes were significantly down‐regulated following EtOH consumption; these were functionally enriched for immune response, cell adhesion, plasma membrane, and extracellular matrix. A total of 567 genes that were up‐regulated following EtOH consumption were enriched in microRNA target sites and included target sites associated with Toll‐like receptor pathways. The differentially expressed genes were also significantly enriched in transcription factor binding sites. Conclusions The data presented here are the first to use a longitudinal biopsy strategy to examine how chronic EtOH consumption affects gene expression in the primate PFC. Prominent effects were seen in both cell adhesion and neuroimmune pathways; the latter contained both pro‐ and antiinflammatory genes. The data also indicate that changes in miRNAs and transcription factors may be important epigenetic regulators of EtOH consumption., Walter et al. used a unique within subject design to detect changes in the brain prefrontal cortex transcriptome associated with chronic ethanol exposure in non‐human primates. The data obtained showed that broad categories of genes were affected, especially those associated with neuro‐immune function, cell adhesion and the extracellular matrix. The data suggest that chronic ethanol exposure has marked functional and structural effects on the prefrontal cortex.

INTRODUCTION: Gestational ethanol exposure is associated with multiple developmental abnormalities, collectively termed fetal alcohol spectrum disorder (FASD). While the majority of women abstain from ethanol following knowledge of pregnancy, one contributing factor to the high FASD prevalence is that pregnancy is not detected until 4–6 weeks. Thus, ethanol consumption continues during the initial stages of fetal development. METHODS: An experimental protocol is described in which rhesus macaques self-administer 1.5 g/kg/day ethanol (or isocaloric maltose-dextrin) prior to pregnancy and through the first 60 days of a 168-day gestation term. Menstrual cycles were monitored, including measurements of circulating estradiol and progesterone levels. The latency to consume 1.5g/kg ethanol and blood ethanol concentration (BEC) were measured. RESULTS: Twenty-eight fetuses (14 ethanol, 14 controls) were generated in this study. Ethanol did not affect menstrual cycles, or the probability of successful breeding. No ethanol-induced gross adverse effects on pregnancy were observed. Individual variability in latency to complete drinking translated into variability in BEC, measured 90 minutes following session start. Drinking latencies in controls and ethanol-drinkers were longer in the second gestational month than in the first. All pregnancies reached the planned experimental time point of G85, G110, or G135, when in utero MRIs were performed, fetuses were delivered by cesarean section, and brains were evaluated with ex vivo procedures, including slice electrophysiology. Fetal tissues have been deposited to the Monkey Alcohol Tissue Research Resource (MATRR). CONCLUSION: This FASD model takes advantage of the similarities between humans and rhesus macaques in gestational length relative to brain development, as well as similarities in ethanol self-administration and metabolism. The daily 1.5 g/kg dose of ethanol through the first trimester does not influence pregnancy success rates. However, pregnancy influences drinking behavior during the second month of pregnancy. Future publications using this model will describe the effect of early-gestational ethanol exposure on anatomical and functional brain development at subsequent gestational ages.

We previously identified a region on chromosome 1 that harbor quantitative trait loci (QTLs) with large effects on alcohol withdrawal risk using both chronic and acute models in mice. Here, using newly created and existing QTL interval-specific congenic (ISC) models, we report the first evidence that this region harbors two distinct alcohol withdrawal QTLs (Alcw11and Alcw12), which underlie 13% and 3–6%, respectively, of the genetic variance in alcohol withdrawal severity measured using the handling-induced convulsion. Our results also precisely localize Alcw11 and Alcw12 to discreet chromosome regions (syntenic with human 1q23.1–23.3) that encompass a limited number of genes with validated genotype-dependent transcript expression and/or non-synonymous sequence variation that may underlie QTL phenotypic effects. ISC analyses also implicate Alcw11and Alcw12 in withdrawal-induced anxiety-like behavior, representing the first evidence for their broader roles in alcohol withdrawal beyond convulsions; but detect no evidence for Alcw12 involvement in ethanol conditioned place preference (CPP) or consumption. Our data point to high-quality candidates for Alcw12, including genes involved in mitochondrial respiration, spatial buffering, and neural plasticity, and to Kcnj9 as a high-quality candidate for Alcw11. Our studies are the first to show, using two null mutant models on different genetic backgrounds, that Kcnj9 −/− mice demonstrate significantly less severe alcohol withdrawal than wildtype littermates using acute and repeated exposure paradigms. We also demonstrate that Kcnj9 −/− voluntarily consume significantly more alcohol (20%, two-bottle choice) than wildtype littermates. Taken together with evidence implicating Kcnj9 in ethanol CPP, our results support a broad role for this locus in ethanol reward and withdrawal phenotypes. In summary, our results demonstrate two distinct chromosome 1 QTLs that significantly affect risk for ethanol withdrawal, and point to their distinct unique roles in alcohol reward phenotypes.

Alcohol use disorder is a significant global burden. Stress has been identified as an etiological factor in the initiation and continuation of ethanol consumption. Understanding adaptations within stress circuitry is an important step toward novel treatment strategies. The effects of protracted abstinence following long-term ethanol self-administration on the central nucleus of the amygdala (CeA) and the hypothalamic paraventricular nucleus (PVN) were evaluated in male rhesus monkeys. Using whole-cell patch-clamp electrophysiology, inhibitory GABAergic transmission in the CeA and excitatory glutamatergic transmission in the PVN were measured. CeA neurons from abstinent drinkers displayed an elevated baseline spontaneous inhibitory postsynaptic current (sIPSC) frequency compared with controls, indicating increased presynaptic GABA release. Application of acute ethanol significantly increased the frequency of sIPSCs in controls, but not in abstinent drinkers, suggesting a tolerance to ethanol-enhanced GABA release in abstinent rhesus monkeys with a history of chronic ethanol self-administration and repeated abstinence. In the PVN, the frequency of spontaneous excitatory postsynaptic currents (sEPSC) was elevated in abstinent drinkers compared with controls, indicating increased presynaptic glutamate release. Notably, acute ethanol decreased presynaptic glutamate release onto parvocellular PVN neurons in both controls and abstinent drinkers, suggesting a lack of tolerance to acute ethanol among PVN neurons. These results are the first to demonstrate distinct synaptic adaptations and ethanol sensitivity in both the extrahypothalamic and hypothalamic stress circuits in abstinent rhesus males. Importantly, our findings describe adaptations in stress circuitry present in the brain at a state during abstinence, just prior to relapse to ethanol drinking.

Genetic factors significantly affect vulnerability to alcohol dependence (alcoholism). We previously identified quantitative trait loci on distal mouse chromosome 1 with large effects on predisposition to alcohol physiological dependence and associated withdrawal following both chronic and acute alcohol exposure in mice (Alcdp1 and Alcw1, respectively). We fine-mapped these loci to a 1.1-1.7 Mb interval syntenic with human 1q23.2-23.3. Alcw1/Alcdp1 interval genes show remarkable genetic variation among mice derived from the C57BL/6J and DBA/2J strains, the two most widely studied genetic animal models for alcohol-related traits. Here, we report the creation of a novel recombinant Alcw1/Alcdp1 congenic model (R2) in which the Alcw1/Alcdp1 interval from a donor C57BL/6J strain is introgressed onto a uniform, inbred DBA/2J genetic background. As expected, R2 mice demonstrate significantly less severe alcohol withdrawal compared to wild-type littermates. Additionally, comparing R2 and background strain animals, as well as reciprocal congenic (R8) and appropriate background strain animals, we assessed Alcw1/Alcdp1 dependent brain gene expression using microarray and quantitative PCR analyses. To our knowledge this includes the first Weighted Gene Co-expression Network Analysis using reciprocal congenic models. Importantly, this allows detection of co-expression patterns limited to one or common to both genetic backgrounds with high or low predisposition to alcohol withdrawal severity. The gene expression patterns (modules) in common contain genes related to oxidative phosphorylation, building upon human and animal model studies that implicate involvement of oxidative phosphorylation in alcohol use disorders. Finally, we demonstrate that administration of N-acetylcysteine, an FDA-approved antioxidant, significantly reduces symptoms of alcohol withdrawal (convulsions) in mice, thus validating a phenotypic role for this network. Taken together, these studies support the importance of mitochondrial oxidative homeostasis in alcohol withdrawal and identify this network as a valuable therapeutic target in human alcohol use disorders.

Association studies implicate the multiple PDZ domain protein (MUPP1/MPDZ) gene in risk for alcoholism in humans and alcohol withdrawal in mice. Although manipulation of the Mpdz gene by homologous recombination and bacterial artificial chromosome transgenesis has suggested that its expression affects alcohol withdrawal risk, the potential confounding effects of linked genes and developmental compensation currently limit interpretation. Here, using RNA interference, we directly test the impact of Mpdz expression on alcohol withdrawal severity and provide brain regional mechanistic information. Lentiviral-mediated delivery of Mpdz short hairpin RNA (shRNA) to the caudolateral substantia nigra pars reticulata significantly reduces Mpdz expression and exacerbates alcohol withdrawal convulsions compared to control mice delivered a scrambled shRNA. Neither baseline nor pentylenetetrazol enhanced convulsions differed between Mpdz shRNA and control animals, indicating that Mpdz expression in the caudolateral substantia nigra pars reticulata does not generally affect seizure susceptibility. To our knowledge, these represent the first in vivo Mpdz RNA interference analyses, and provide the first direct evidence that Mpdz expression impacts behavior. Our results confirm that Mpdz is a quantitative trait gene for alcohol withdrawal and demonstrate that its expression in the caudolateral substantia nigra pars reticulata is crucially involved in risk for alcohol withdrawal.

Attentional set-shifting ability is an executive function underling cognitive flexibility in humans and animals. In humans, this function is typically observed during a single experimental session where dimensions of playing cards are used to measure flexibility in the face of changing rules for reinforcement (i.e., the Wisconsin Card Sorting Test (WCST)). In laboratory animals, particularly non-human primates, variants of the WCST involve extensive training and testing on a series of dimensional discriminations, usually in social isolation. In the present study, a novel experimental approach was used to assess attentional set-shifting simultaneously in 12 rhesus monkeys. Specifically, monkeys living in individual cages but in the same room were trained at the same time each day in a set-shifting task in the same housing environment. As opposed to the previous studies, each daily session began with a simple single-dimension discrimination regardless of the animal's performance on the previous session. A total of eight increasingly difficult, discriminations (sets) were possible in each daily 45 min session. Correct responses were reinforced under a second-order schedule of flavored food pellet delivery, and criteria for completing a set was 12 correct trials out of a running total of 15 trials. Monkeys progressed through the sets at their own pace and abilities. The results demonstrate that all 12 monkeys acquired the simple discrimination (the first set), but individual differences in the ability to progress through all eight sets were apparent. A performance index (PI) that encompassed progression through the sets, errors and session duration was calculated and used to rank each monkey's performance in relation to each other. Overall, this version of a set-shifting task results in an efficient assessment of reliable differences in cognitive flexibility in a group of monkeys.

Mitochondrial dysfunction is implicated in the etiology and pathogenesis of numerous human disorders involving tissues with high energy demand. Murine models are widely used to elucidate genetic determinants of phenotypes relevant to human disease, with recent studies of C57BL/6J (B6), DBA/2J (D2) and B6xD2 populations implicating naturally occurring genetic variation in mitochondrial function/dysfunction. Using blue native polyacrylamide gel electrophoresis, immunoblots and in-gel activity analyses of complexes I, II, III, IV and V, our studies are the first to assess abundance, organization and catalytic activity of mitochondrial respiratory complexes and supercomplexes in mouse brain. Remarkable strain differences in supercomplex assembly and associated activity are evident, without differences in individual complexes I, II, III or IV. Supercomplexes I1 III2 IV2-3 exhibit robust complex III immunoreactivity and activities of complexes I and IV in D2, but with little detected in B6 for I1 III2 IV2 , and I1 III2 IV3 is not detected in B6. I1 III2 IV1 and I1 III2 are abundant and catalytically active in both strains, but significantly more so in B6. Furthermore, while supercomplex III2 IV1 is abundant in D2, none is detected in B6. In aggregate, these results indicate a shift toward more highly assembled supercomplexes in D2. Respiratory supercomplexes are thought to increase electron flow efficiency and individual complex stability, and to reduce electron leak and generation of reactive oxygen species. Our results provide a framework to begin assessing the role of respiratory complex suprastructure in genetic vulnerability and treatment for a wide variety of mitochondrial-related disorders.

Association studies implicate multiple PDZ domain protein (MPDZ/MUPP1) sequence and/or expression in risk for alcoholism in humans and ethanol withdrawal (EW) in mice, but confirmation has been hindered by the dearth of targeted genetic models. We report the creation of transgenic (MPDZ-TG) and knockout heterozygote (Mpdz(+/-) ) mice, with increased (2.9-fold) and decreased (53%) target expression, respectively. Both models differ in EW compared with wild-type littermates (P ≤ 0.03), providing compelling evidence for an inverse relationship between Mpdz expression and EW severity. Additionally, ethanol consumption is reduced up to 18% (P = 0.006) in Mpdz(+/-) , providing the first evidence implicating Mpdz in ethanol self-administration.

α-Synuclein (α-syn) protein and endocannabinoid CB1 receptors are primarily located in presynaptic terminals. An association between α-syn and CB1 receptors has recently been established in Parkinson's disease, but it is completely unknown whether there is an association between these two proteins in alcohol addiction. Therefore, we aimed to examine the α-syn mRNA transcript and protein expression levels in the prefrontal cortex, striatum, amygdala and hippocampus. These brain regions are the most frequently implicated in alcohol and other drug addiction. In these studies, we used C57BL/6 mice carrying a spontaneous deletion of the α-syn gene (C57BL/6Snca-/-) and their respective controls (C57BL/6Snca+/+). These animals were monitored for spontaneous alcohol consumption (3–10%) and their response to a hypnotic-sedative dose of alcohol (3 g kg−1) was also assessed. Compared with the C57BL/6Snca+/+ mice, we found that the C57BL/6Snca-/- mice exhibited a higher expression level of the CB1 mRNA transcript and CB1 receptor in the hippocampus and amygdala. Furthermore, C57BL/6Snca-/- mice showed an increase in alcohol consumption when offered a 10% alcohol solution. There was no significant difference in sleep time after the injection of 3 g/kg alcohol. These results are the first to reveal an association between α-syn and the CB1 receptor in the brain regions that are most frequently implicated in alcohol and other drug addictions. Synapse 00:000–000, 2013. © 2013 Wiley Periodicals, Inc.

Advances in next-generation sequencing suggest that RNA-Seq is poised to supplant microarray-based approaches for transcriptome analysis. This article briefly reviews the use of microarrays in the brain-behavior context and then illustrates why RNA-Seq is a superior strategy. Compared with microarrays, RNA-Seq has a greater dynamic range, detects both coding and noncoding RNAs, is superior for gene network construction, detects alternative spliced transcripts, detects allele specific expression and can be used to extract genotype information, e.g. nonsynonymous coding single nucleotide polymorphisms. Examples of where RNA-Seq has been used to assess brain gene expression are provided. Despite the advantages of RNA-Seq, some disadvantages remain. These include the high cost of RNA-Seq and the computational complexities associated with data analysis. RNA-Seq embraces the complexity of the transcriptome and provides a mechanism to understand the underlying regulatory code; the potential to inform the brain-behavior relationship is substantial.

The current study examined the changes in striatal gene network structure induced by short-term selective breeding from a heterogeneous stock for haloperidol response. Brain (striatum) gene expression data were obtained using the Illumina WG 8.2 array, and the datasets from responding and non-responding selected lines were independently interrogated using a weighted gene coexpression network analysis (WGCNA). We detected several gene modules (groups of coexpressed genes) in each dataset; the membership of the modules was found to be largely concordant, and a consensus network was constructed. Further validation of the network topology showed that using approximately 35 samples is sufficient to reliably infer the transcriptome network. An in-depth analysis showed significant changes in network structure and gene connectivity associated with the selected lines; these changes were validated using a bootstrapping procedure. The most dramatic changes were associated with a gene module richly annotated with neurobehavioral traits. The changes in network connectivity were concentrated in the links between this module and the rest of the network, in addition to changes within the module; this observation is consistent with recent results in protein and metabolic networks. These results suggest that a network-based strategy will help identify the genetic factors associated with haloperidol response.

Evidence for genetic linkage to alcohol and other substance dependence phenotypes in areas of the human and mouse genome have now been reported with some consistency across studies. However, the question remains as to whether the genes that underlie the alcohol-related behaviors seen in mice are the same as those that underlie the behaviors observed in human alcoholics. The aims of the current set of analyses were to identify a small set of alcohol-related phenotypes in human and in mouse by which to compare quantitative trait locus (QTL) data between the species using syntenic mapping. These analyses identified that QTLs for alcohol consumption and acute and chronic alcohol withdrawal on distal mouse chromosome 1 are syntenic to a region on human chromosome 1q where a number of studies have identified QTLs for alcohol-related phenotypes. Additionally, a QTL on human chromosome 15 for alcohol dependence severity/withdrawal identified in two human studies was found to be largely syntenic with a region on mouse chromosome 9, where two groups have found QTLs for alcohol preference. In both of these cases, while the QTLs were found to be syntenic, the exact phenotypes between humans and mice did not necessarily overlap. These studies demonstrate how this technique might be useful in the search for genes underlying alcohol-related phenotypes in multiple species. However, these findings also suggest that trying to match exact phenotypes in humans and mice may not be necessary or even optimal for determining whether similar genes influence a range of alcohol-related behaviors between the two species.

Here, we map a quantitative trait locus (QTL) with a large effect on predisposition to barbiturate (pentobarbital) withdrawal to a 0.44 Mb interval of mouse chromosome 1 syntenic with human 1q23.2. We report a detailed analysis of the genes within this interval and show that it contains 15 known and predicted genes, 12 of which demonstrate validated genotype-dependent transcript expression and/or nonsynonymous coding sequence variation that may underlie the influence of the QTL on withdrawal. These candidates are involved in diverse cellular functions including intracellular trafficking, potassium conductance and spatial buffering, and multimolecular complex dynamics, and indicate both established and novel aspects of neurobiological response to sedative-hypnotics. This work represents a substantial advancement toward identification of the gene(s) that underlie the phenotypic effects of the QTL. We identifyKcnj9as a particularly promising candidate and report the development of aKcnj9-null mutant model that exhibits significantly less severe withdrawal from pentobarbital as well as other sedative-hypnotics (zolpidem and ethanol) versus wild-type littermates. Reduced expression ofKcnj9, which encodes GIRK3 (Kir3.3), is associated with less severe sedative-hypnotic withdrawal. A multitude of QTLs for a variety of complex traits, including diverse responses to sedative-hypnotics, have been detected on distal chromosome 1 in mice, and as many as four QTLs on human chromosome 1q have been implicated in human studies of alcohol dependence. Thus, our results will be primary to additional efforts to identify genes involved in a wide variety of behavioral responses to sedative-hypnotics and may directly facilitate progress in human genetics.

We report here the confirmation of the quantitative trait locus for haloperidol-induced catalepsy on distal chromosome (Chr) 1. We determined that this quantitative trait locus was captured in the B6.D2-Mtv7a/Ty congenic mouse strain, whose introgressed genomic interval extends from approximately 169.1 to 191.3 Mb. We then constructed a group of overlapping interval-specific congenic strains to further break up the interval and remapped the locus between 177.5 and 183.4 Mb. We next queried single nucleotide polymorphism (SNP) data sets and identified three genes with nonsynonymous coding SNPs in the quantitative trait locus. We also queried two brain gene expression data sets and found five known genes in this 5.9-Mb interval that are differentially expressed in both whole brain and striatum. Three of the candidate quantitative trait genes were differentially expressed using quantitative real-time polymerase chain reaction analyses. Overall, the current study illustrates how multiple approaches, including congenic fine mapping, SNP analysis and microarray gene expression screens, can be integrated both to reduce the quantitative trait locus interval significantly and to detect promising candidate quantitative trait genes.

Cell adhesion molecules are involved in a diverse array of cellular processes. Recent data suggests that human immunodeficiency virus (HIV-1) co-opts their functions, in particular the properties of the intercellular cell adhesion molecules (ICAMs), to enhance viral infection and transmission. To investigate mechanisms that may underlie the non-progression that occurs in immunodeficiency virus-infected chimpanzees, we amplified the protein coding regions of multiple non-human primate ICAMs 1–5 and two ICAM ligands, leukocyte function-associated antigen-1 (LFA-1) and macrophage antigen 1 (Mac-1). We then employed a phylogenetic tree-based approach to comparative genomics, in order to screen for the presence of adaptive changes. Strong Darwinian positive selection in chimpanzee ICAMs 1, 2 and 3 was observed, most markedly in domains that are critical for the integrity and maintenance of ICAM-1 dimerization. As binding of ligands, including the attachment of virions, is influenced by the state of ICAM 1 dimerization, chimpanzee ICAMs may have evolved to modulate their own dimerization. In concert with previous evidence suggesting an ancient retroviral pandemic as a prominent selective force in chimpanzee evolution, adaptation of chimpanzee ICAMs may have effected a mechanism that explains the lack of immunosuppression observed following HIV-1 or simian immunodeficiency virus (SIVcpz) infection.

Various lines of evidence suggest that a polymorphism in the gene for proopiomelanocortin, Pomc1, might account for some portion of the genetic variance for open-field activity in the LS × SS RI strains. To test this hypothesis, approximately 1600 bp of Pomc1was sequenced from genomic DNA of seven of the LS × SS strains. Two distinct alleles containing a total of five single-base pair differences were detected. A substitution was found in the coding region causing a Pro-to-Ser conversion, two substitutions occurred in the 3′ untranslated region of the mRNA, and a substitution and a deletion were detected in the 3′ untranscribed flanking region. The fragment containing the coding region substitution was sequenced in an additional 15 of the LS × SS strains. A total of 12 strains contained one form of the allele, while 10 had the other. The strain distribution pattern of open-field activity scores between the two alleles suggests that these alleles do not contribute to the genetic variation of open-field activity in the LS × SS RI strains.

Background Data from C57BL/6J (B6) × DBA/2J (D2) F2 intercrosses (B6xD2 F2), standard and recombinant inbred strains, and heterogeneous stock mice indicate that a reciprocal (or inverse) genetic relationship exists between alcohol consumption and withdrawal severity. Furthermore, some genetic studies have detected reciprocal quantitative trait loci (QTLs) for these traits. We used a novel mouse model developed by simultaneous selection for both high alcohol consumption/low withdrawal and low alcohol consumption/high withdrawal and analyzed the gene expression and genome-wide genotypic differences. Methods Randomly chosen third selected generation (S3) mice (N = 24/sex/line), bred from a B6xD2 F2, were genotyped using the Mouse Universal Genotyping Array, which provided 2,760 informative markers. QTL analysis used a marker-by-marker strategy with the threshold for a significant log of the odds (LOD) set at 10. Gene expression in the ventral striatum was measured using the Illumina Mouse 8.2 array. Differential gene expression and the weighted gene co-expression network analysis (WGCNA) were implemented. Results Significant QTLs for consumption/withdrawal were detected on chromosomes (Chr) 2, 4, 9, and 12. A suggestive QTL mapped to Chr 6. Some of the QTLs overlapped with known QTLs mapped for 1 of the traits individually. One thousand seven hundred and forty-five transcripts were detected as being differentially expressed between the lines; there was some overlap with known withdrawal genes (e.g., Mpdz) located within QTL regions. WGCNA revealed several modules of co-expressed genes showing significant effects in both differential expression and intramodular connectivity; a module richly annotated with kinase-related annotations was most affected. Conclusions Marked effects of selection on expression and network structure were detected. QTLs overlapping with differentially expressed genes on Chr 2 (distal) and 4 suggest that these are cis-eQTLs (Chr 2: Kif3b, Kcnq2; Chr 4: Mpdz, Snapc3). Other QTLs identified were on Chr 2 (proximal), 9, and 12. Network results point to involvement of kinase-related mechanisms and outline the need for further efforts such as interrogation of noncoding RNAs.

The human genome is thought to harbor 50,000 to 100,000 genes, of which about half have been sampled to date in the form of expressed sequence tags. An international consortium was organized to develop and map gene-based sequence tagged site markers on a set of two radiation hybrid panels and a yeast artificial chromosome library. More than 16,000 human genes have been mapped relative to a framework map that contains about 1000 polymorphic genetic markers. The gene map unifies the existing genetic and physical maps with the nucleotide and protein sequence databases in a fashion that should speed the discovery of genes underlying inherited human disease. The integrated resource is available through a site on the World Wide Web at http://www.ncbi.nlm.nih.gov/SCIENCE96/.

Using our data set of 3,143 single pass sequences from human brain cDNA libraries, we have developed a strategy in which gene-based sequence-tagged-sites (STSs), derived from 3'untranslated regions of human cDNAs, are rapidly assigned to megabase-insert yeast artificial chromosomes and somatic cell hybrids to generate regional gene mapping data. Employing this approach, we have mapped 318 cDNAs, representing 308 human genes. Ninety-two of these mapped to regions implicated in human genetic diseases, identifying them as candidate genes. Extension of this strategy has the potential to result in virtually every human gene having, at its 3' end, its own associated STS, with each STS in turn specifying both a corresponding genomic clone and a specific regional location in the genome.

Methamphetamine (MA) and neurotransmitter precursors and metabolites such as tyramine, octopamine, and β-phenethylamine stimulate the G protein-coupled trace amine-associated receptor 1 (TAAR1). TAAR1 has been implicated in human conditions including obesity, schizophrenia, depression, fibromyalgia, migraine, and addiction. Additionally TAAR1 is expressed on lymphocytes and astrocytes involved in inflammation and response to infection. In brain, TAAR1 stimulation reduces synaptic dopamine availability and alters glutamatergic function. TAAR1 is also expressed at low levels in heart, and may regulate cardiovascular tone. Taar1 knockout mice orally self-administer more MA than wild type and are insensitive to its aversive effects. DBA/2J (D2) mice express a non-synonymous single nucleotide polymorphism (SNP) in Taar1 that does not respond to MA, and D2 mice are predisposed to high MA intake, compared to C57BL/6 (B6) mice. Here we demonstrate that endogenous agonists stimulate the recombinant B6 mouse TAAR1, but do not activate the D2 mouse receptor. Progeny of the B6XD2 (BxD) family of recombinant inbred (RI) strains have been used to characterize the genetic etiology of diseases, but contrary to expectations, BXDs derived 30–40 years ago express only the functional B6 Taar1 allele whereas some more recently derived BXD RI strains express the D2 allele. Data indicate that the D2 mutation arose subsequent to derivation of the original RIs. Finally, we demonstrate that SNPs in human TAAR1 alter its function, resulting in expressed, but functional, sub-functional and non-functional receptors. Our findings are important for identifying a predisposition to human diseases, as well as for developing personalized treatment options.

Background Allelic variation is the cornerstone of genetically determined differences in gene expression, gene product structure, physiology, and behavior. However, allelic variation, particularly cryptic (unknown or not annotated) variation, is problematic for follow up analyses. Polymorphisms result in a high incidence of false positive and false negative results in hybridization based analyses and hinder the identification of the true variation underlying genetically determined differences in physiology and behavior. Given the proliferation of mouse genetic models (e.g., knockout models, selectively bred lines, heterogeneous stocks derived from standard inbred strains and wild mice) and the wealth of gene expression microarray and phenotypic studies using genetic models, the impact of naturally-occurring polymorphisms on these data is critical. With the advent of next-generation, high-throughput sequencing, we are now in a position to determine to what extent polymorphisms are currently cryptic in such models and their impact on downstream analyses. Results We sequenced the two most commonly used inbred mouse strains, DBA/2J and C57BL/6J, across a region of chromosome 1 (171.6 – 174.6 megabases) using two next generation high-throughput sequencing platforms: Applied Biosystems (SOLiD) and Illumina (Genome Analyzer). Using the same templates on both platforms, we compared realignments and single nucleotide polymorphism (SNP) detection with an 80 fold average read depth across platforms and samples. While public datasets currently annotate 4,527 SNPs between the two strains in this interval, thorough high-throughput sequencing identified a total of 11,824 SNPs in the interval, including 7,663 new SNPs. Furthermore, we confirmed 40 missense SNPs and discovered 36 new missense SNPs. Conclusion Comparisons utilizing even two of the best characterized mouse genetic models, DBA/2J and C57BL/6J, indicate that more than half of naturally-occurring SNPs remain cryptic. The magnitude of this problem is compounded when using more divergent or poorly annotated genetic models. This warrants full genomic sequencing of the mouse strains used as genetic models.

As described in other articles in this Special Section, microarrays are widely used to evaluate gene expression at the genome scale. However, all too often the importance of data analysis at the level of the individual probe is overlooked. This is a particular problem when trying to detect differences in gene expression levels among genetically unique animals, across inbred animal strains, or among genetically modified animals. Of particular concern is the presence of small modifications in the DNA (i.e., single nucleotide polymorphisms [SNPs]) that occur naturally and differentiate one individual from the next. This article describes the potential impact of SNPs on analyses of gene expression differences and introduces an approach called SNP masking, which implements removal of SNP-affected probes. SNP masking is a valuable and feasible approach that can ameliorate these hybridization problems.

Physiological dependence and associated withdrawal episodes can constitute a powerful motivational force that perpetuates drug use and abuse. Using robust behavioral models of drug physiological dependence in mice, positional cloning, and sequence and expression analyses, we identified an addiction-relevant quantitative trait gene, Mpdz. Our findings provide a framework to define the protein interactions and neural circuit by which this gene's product (multiple PDZ domain protein) affects drug dependence, withdrawal and relapse.

We describe a method for the rapid identification and mapping of human genes, including those possibly contributing to disease and alcohol-related phenotypes. New human genes are identified from cDNA libraries through single-pass sequencing into the 3' untranslated (3'UT) regions of human brain cDNAs. Primers derived from the 3'UT region sequences [representing gene-based, sequence-tagged sites (STSs)] are used for polymerase chain reaction (PCR) analyses of the CEPH megabase insert yeast artificial chromosome (YAC) DNA pools. With this approach, approximately 18,000 megabase YACs can be screened and a single YAC identified using only 52 PCR reactions. The YAC localization in conjunction with other mapping techniques, such as PCR mapping to human chromosomes using somatic cell hybrids, allows identification of chromosomal band locations. In this manner, each gene can be associated with its own STS, which in turn specifies both a corresponding genomic clone and specific location in the genome. These locations can be compared with the purported locations of disease genes. The locations of the STSs can also be compared with those of Quantitative Trait Loci implicated for quantitative traits (e.g., alcohol-related phenotypes) on the basis of synteny between the mouse and human genes. Using this strategy, we found candidates for 78 human disease/syndrome genes among the first 220 genes mapped.

C57BL/6J (B6) and DBA/2J (D2) are two of the most commonly used inbred mouse strains in neuroscience research. However, the only currently available mouse genome is based entirely on the B6 strain sequence. Subsequently, oligonucleotide microarray probes are based solely on this B6 reference sequence, making their application for gene expression profiling comparisons across mouse strains dubious due to their allelic sequence differences, including single nucleotide polymorphisms (SNPs). The emergence of next-generation sequencing (NGS) and the RNA-Seq application provides a clear alternative to oligonucleotide arrays for detecting differential gene expression without the problems inherent to hybridization-based technologies. Using RNA-Seq, an average of 22 million short sequencing reads were generated per sample for 21 samples (10 B6 and 11 D2), and these reads were aligned to the mouse reference genome, allowing 16,183 Ensembl genes to be queried in striatum for both strains. To determine differential expression, 'digital mRNA counting' is applied based on reads that map to exons. The current study compares RNA-Seq (Illumina GA IIx) with two microarray platforms (Illumina MouseRef-8 v2.0 and Affymetrix MOE 430 2.0) to detect differential striatal gene expression between the B6 and D2 inbred mouse strains. We show that by using stringent data processing requirements differential expression as determined by RNA-Seq is concordant with both the Affymetrix and Illumina platforms in more instances than it is concordant with only a single platform, and that instances of discordance with respect to direction of fold change were rare. Finally, we show that additional information is gained from RNA-Seq compared to hybridization-based techniques as RNA-Seq detects more genes than either microarray platform. The majority of genes differentially expressed in RNA-Seq were only detected as present in RNA-Seq, which is important for studies with smaller effect sizes where the sensitivity of hybridization-based techniques could bias interpretation.

Background The current study focused on the extent genetic diversity within a species (Mus musculus) affects gene co-expression network structure. To examine this issue, we have created a new mouse resource, a heterogeneous stock (HS) formed from the same eight inbred strains that have been used to create the collaborative cross (CC). The eight inbred strains capture > 90% of the genetic diversity available within the species. For contrast with the HS-CC, a C57BL/6J (B6) × DBA/2J (D2) F2 intercross and the HS4, derived from crossing the B6, D2, BALB/cJ and LP/J strains, were used. Brain (striatum) gene expression data were obtained using the Illumina Mouse WG 6.1 array, and the data sets were interrogated using a weighted gene co-expression network analysis (WGCNA). Results Genes reliably detected as expressed were similar in all three data sets as was the variability of expression. As measured by the WGCNA, the modular structure of the transcriptome networks was also preserved both on the basis of module assignment and from the perspective of the topological overlap maps. Details of the HS-CC gene modules are provided; essentially identical results were obtained for the HS4 and F2 modules. Gene ontology annotation of the modules revealed a significant overrepresentation in some modules for neuronal processes, e.g., central nervous system development. Integration with known protein-protein interactions data indicated significant enrichment among co-expressed genes. We also noted significant overlap with markers of central nervous system cell types (neurons, oligodendrocytes and astrocytes). Using the Allen Brain Atlas, we found evidence of spatial co-localization within the striatum for several modules. Finally, for some modules it was possible to detect an enrichment of transcription binding sites. The binding site for Wt1, which is associated with neurodegeneration, was the most significantly overrepresented. Conclusions Despite the marked differences in genetic diversity, the transcriptome structure was remarkably similar for the F2, HS4 and HS-CC. These data suggest that it should be possible to integrate network data from simple and complex crosses. A careful examination of the HS-CC transcriptome revealed the expected structure for striatal gene expression. Importantly, we demonstrate the integration of anatomical and network expression data.

Background The Collaborative Cross (CC) is a large panel of genetically diverse recombinant inbred mouse strains specifically designed to provide a systems genetics resource for the study of complex traits. In part, the utility of the CC stems from the extensive genome-wide annotations of founder strain sequence and structural variation. Still missing, however, are transcriptome-specific annotations of the CC founder strains that could further enhance the utility of this resource. Results We provide a comprehensive survey of the splicing landscape of the 8 CC founder strains by leveraging the high level of alternative splicing within the brain. Using deep transcriptome sequencing, we found that a majority of the splicing landscape is conserved among the 8 strains, with ~65% of junctions being shared by at least 2 strains. We, however, found a large number of potential strain-specific splicing events as well, with an average of ~3000 and ~500 with ≥3 and ≥10 sequence read coverage, respectively, within each strain. To better understand strain-specific splicing within the CC founder strains, we defined criteria for and identified high-confidence strain-specific splicing events. These splicing events were defined as exon-exon junctions 1) found within only one strain, 2) with a read coverage ≥10, and 3) defined by a canonical splice site. With these criteria, a total of 1509 high-confidence strain-specific splicing events were identified, with the majority found within two of the wild-derived strains, CAST and PWK. Strikingly, the overwhelming majority, 94%, of these strain-specific splicing events are not yet annotated. Strain-specific splicing was also located within genomic regions recently reported to be over- and under-represented within CC populations. Conclusions Phenotypic characterization of CC populations is increasing; thus these results will not only aid in further elucidating the transcriptomic architecture of the individual CC founder strains, but they will also help in guiding the utilization of the CC populations in the study of complex traits. This report is also the first to establish guidelines in defining and identifying strain-specific splicing across different mouse strains. Electronic supplementary material The online version of this article (doi:10.1186/s12864-015-1267-0) contains supplementary material, which is available to authorized users.