RNA-Seq Analysis of Genetic and Transcriptome Network Effects of Dual-Trait Selection for Ethanol Preference and Withdrawal Using SOT and NOT Genetic Models

Identification of novel genetic loci and candidate genes for progressive ethanol consumption in diversity outbred mice

Mouse behavioral genetic mapping studies can identify genomic intervals modulating complex traits under well-controlled environmental conditions and have been used to study ethanol behaviors to aid in understanding genetic risk and the neurobiology of alcohol use disorder (AUD). However, historically such studies have produced large confidence intervals, thus complicating identification of potential causal candidate genes. Diversity Outbred (DO) mice offer the ability to perform high-resolution quantitative trait loci (QTL) mapping on a very genetically diverse background, thus facilitating identification of candidate genes. Here, we studied a population of 636 male DO mice with four weeks of intermittent ethanol access via a three-bottle choice procedure, producing a progressive ethanol consumption phenotype. QTL analysis identified 3 significant (Chrs 3, 4, and 12) and 13 suggestive loci for ethanol-drinking behaviors with narrow confidence intervals (1-4 Mbp for significant QTLs). Results suggested that genetic influences on initial versus progressive ethanol consumption were localized to different genomic intervals. A defined set of positional candidate genes were prioritized using haplotype analysis, identified coding polymorphisms, prefrontal cortex transcriptomics data, human GWAS data and prior rodent gene set data for ethanol or other misused substances. These candidates included Car8, the lone gene with a significant cis-eQTL within a Chr 4 QTL for week four ethanol consumption. These results represent the highest-resolution genetic mapping of ethanol consumption behaviors in mice to date, providing identification of novel loci and candidate genes for study in relation to the neurobiology of AUD.

Effects of repeated alcohol abstinence on within-subject prefrontal cortical gene expression in rhesus macaques

Male rhesus monkeys (n = 24) had a biopsy of prefrontal cortical area 46 prior to chronic ethanol self-administration (n = 17) or caloric control (n = 7). Fourteen months of daily self-administration (water vs. 4% alcohol, 22 h access/day termed "open-access") was followed by two cycles of prolonged abstinence (5 weeks) each followed by 3 months of open-access alcohol and a final abstinence followed by necropsy. At necropsy, a biopsy of Area 46, contralateral to the original biopsy, was obtained. Gene expression data (RNA-Seq) were collected comparing biopsy/necropsy samples. Monkeys were categorized by drinking status during the final post-abstinent drinking phase as light (LD), binge (BD), heavy (HD) and very heavy (VHD drinkers). Comparing pre-ethanol to post-abstinent biopsies, four animals that converted from HD to VHD status had significant ontology enrichments in downregulated genes (necropsy minus biopsy n = 286) that included immune response (FDR < 9 × 10-7) and plasma membrane changes (FDR < 1 × 10-7). Genes in the immune response category included IL16 and 18, CCR1, B2M, TLR3, 6 and 7, SP2 and CX3CR1. Upregulated genes (N = 388) were particularly enriched in genes associated with the negative regulation of MAP kinase activity (FDR < 3 × 10-5), including DUSP 1, 4, 5, 6 and 18, SPRY 2, 3, and 4, SPRED2, BMP4 and RGS2. Overall, these data illustrate the power of the NHP model and the within-subject design of genomic changes due to alcohol and suggest new targets for treating severe escalated drinking following repeated alcohol abstinence attempts.

Inflammatory Factor IL1α Induces Aberrant Astrocyte Proliferation in Spinal Cord Injury Through the Grin2c/Ca2+/CaMK2b Pathway

Spinal cord injury (SCI) is one of the most devastating traumas, and the aberrant proliferation of astrocytes usually causes neurological deficits. However, the mechanism underlying astrocyte over-proliferation after SCI is unclear. Grin2c (glutamate ionotropic receptor type 2c) plays an essential role in cell proliferation. Our bioinformatic analysis indicated that Grin2c and Ca2+ transport functions were inhibited in astrocytes after SCI. Suppression of Grin2c stimulated astrocyte proliferation by inhibiting the Ca2+/calmodulin-dependent protein kinase 2b (CaMK2b) pathway in vitro. By screening different inflammatory factors, interleukin 1α (IL1α) was further found to inhibit Grin2c/Ca2+/CaMK2b and enhance astrocyte proliferation in an oxidative damage model. Blockade of IL1α using neutralizing antibody resulted in increased Grin2c expression and the inhibition of astrocyte proliferation post-SCI. Overall, this study suggests that IL1α promotes astrocyte proliferation by suppressing the Grin2c/Ca2+/CaMK2b pathway after SCI, revealing a novel pathological mechanism of astrocyte proliferation, and may provide potential targets for SCI repair.

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

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

Brain gene expression differences related to ethanol preference in the collaborative cross founder strains

The collaborative cross (CC) founder strains include five classical inbred laboratory strains [129S1/SvlmJ (S129), A/J (AJ), C57BL/6J (B6), NOD/ShiLtJ (NOD), and NZO/HILtJ (NZO)] and three wild-derived strains [CAST/EiJ (CAST), PWK/PhJ (PWK), and WSB/EiJ (WSB)]. These strains encompass 89% of the genetic diversity available in Mus musculus and ∼10-20 times more genetic diversity than found in Homo sapiens. For more than 60 years the B6 strain has been widely used as a genetic model for high ethanol preference and consumption. However, another of the CC founder strains, PWK, has been identified as a high ethanol preference/high consumption strain. The current study determined how the transcriptomes of the B6 and PWK strains differed from the 6 low preference CC strains across 3 nodes of the brain addiction circuit. RNA-Seq data were collected from the central nucleus of the amygdala (CeA), the nucleus accumbens core (NAcc) and the prelimbic cortex (PrL). Differential expression (DE) analysis was performed in each of these brain regions for all 28 possible pairwise comparisons of the CC founder strains. Unique genes for each strain were identified by selecting for genes that differed significantly [false discovery rate (FDR) < 0.05] from all other strains in the same direction. B6 was identified as the most distinct classical inbred laboratory strain, having the highest number of total differently expressed genes (DEGs) and DEGs with high log fold change, and unique genes compared to other CC strains. Less than 50 unique DEGs were identified in common between B6 and PWK within all three brain regions, indicating the strains potentially represent two distinct genetic signatures for risk for high ethanol-preference. 338 DEGs were found to be commonly different between B6, PWK and the average expression of the remaining CC strains within all three regions. The commonly different up-expressed genes were significantly enriched (FDR < 0.001) among genes associated with neuroimmune function. These data compliment findings showing that neuroimmune signaling is key to understanding alcohol use disorder (AUD) and support use of these 8 strains and the highly heterogeneous mouse populations derived from them to identify alcohol-related brain mechanisms and treatment targets.

Neurobeachin, a promising target for use in the treatment of alcohol use disorder

Hazardous, heavy drinking increases risk for developing alcohol use disorder (AUD), which affects ~7% of adult Americans. Thus, understanding the molecular mechanisms promoting risk for heavy drinking is essential to developing more effective AUD pharmacotherapies than those currently approved by the FDA. Using genome-wide bisulfate sequencing, we identified DNA methylation (DNAm) signals within the nucleus accumbens core (NAcC) that differentiate nonheavy and heavy ethanol-drinking rhesus macaques. One differentially DNAm region (D-DMR) located within the gene neurobeachin (NBEA), which promotes synaptic membrane protein trafficking, was hypermethylated in heavy drinking macaques. A parallel study identified a similar NBEA D-DMR in human NAcC that distinguished alcoholic and nonalcoholic individuals. To investigate the role of NBEA in heavy ethanol drinking, we engineered a viral vector carrying a short hairpin RNA (shRNA) to reduce the expression of NBEA. Using two murine models of ethanol consumption: 4 days of drinking-in-the-dark and 4 weeks of chronic intermittent access, the knockdown of NBEA expression did not alter average ethanol consumption in either model. However, it did lead to a significant increase in the ethanol preference ratio. Following withdrawal, whole-cell patch clamp electrophysiological experiments revealed that Nbea knockdown led to an increase in spontaneous excitatory postsynaptic current amplitude with no alteration in spontaneous inhibitory postsynaptic currents, suggesting a specific role of NBEA in trafficking of glutamatergic receptors. Together, our findings suggest that NBEA could be targeted to modulate the preference for alcohol use.

On the Use of Heterogeneous Stock Mice to Map Transcriptomes Associated With Excessive Ethanol Consumption

We and many others have noted the advantages of using heterogeneous (HS) animals to map genes and gene networks associated with both behavioral and non-behavioral phenotypes. Importantly, genetically complex Mus musculus crosses provide substantially increased resolution to examine old and new relationships between gene expression and behavior. Here we report on data obtained from two HS populations: the HS/NPT derived from eight inbred laboratory mouse strains and the HS-CC derived from the eight collaborative cross inbred mouse strains that includes three wild-derived strains. Our work has focused on the genes and gene networks associated with risk for excessive ethanol consumption, individual variation in ethanol consumption and the consequences, including escalation, of long-term ethanol consumption. Background data on the development of HS mice is provided, including advantages for the detection of expression quantitative trait loci. Examples are also provided of using HS animals to probe the genes associated with ethanol preference and binge ethanol consumption.

Alcohol Dependence in Rats Is Associated with Global Changes in Gene Expression in the Central Amygdala

Alcohol dependence is associated with adverse consequences of alcohol (ethanol) use and is evident in most severe cases of alcohol use disorder (AUD). The central nucleus of the amygdala (CeA) plays a critical role in the development of alcohol dependence and escalation of alcohol consumption in dependent subjects. Molecular mechanisms underlying the CeA-driven behavioral changes are not well understood. Here, we examined the effects of alcohol on global gene expression in the CeA using a chronic intermittent ethanol (CIE) vapor model in rats and RNA sequencing (RNA-Seq). The CIE procedure resulted in robust changes in CeA gene expression during intoxication, as the number of differentially expressed genes (DEGs) was significantly greater than those expected by chance. Over-representation analysis of cell types, functional groups and molecular pathways revealed biological categories potentially important for the development of alcohol dependence in our model. Genes specific for astrocytes, myelinating oligodendrocytes, and endothelial cells were over-represented in the DEG category, suggesting that these cell types were particularly affected by the CIE procedure. The majority of the over-represented functional groups and molecular pathways were directly related to the functions of glial and endothelial cells, including extracellular matrix (ECM) organization, myelination, and the regulation of innate immune response. A coordinated regulation of several ECM metalloproteinases (e.g., Mmp2; Mmp14), their substrates (e.g., multiple collagen genes and myelin basic protein; Mbp), and a metalloproteinase inhibitor, Reck, suggests a specific mechanism for ECM re-organization in response to chronic alcohol, which may modulate neuronal activity and result in behavioral changes, such as an escalation of alcohol drinking. Our results highlight the importance of glial and endothelial cells in the effects of chronic alcohol exposure on the CeA, and demonstrate further insight into the molecular mechanisms of alcohol dependence in rats. These molecular targets may be used in future studies to develop therapeutics to treat AUD.