A systems genetic analysis of alcohol drinking by mice, rats and men: influence of brain GABAergic transmission

Genetic vitamin B6 deficiency and alcohol interaction in behavior and metabolism

Alcohol abuse is a leading cause of preventable deaths, affecting brain function and metabolism, including GABA transmission and vitamin B6 (VB6) levels. However, the interaction between genetic VB6 deficiency and alcohol consumption remains unexplored. Here, we utilized Drosophila models with mutations in pyridox(am)ine-5'-phosphate oxidase (PNPO), a key enzyme in VB6 metabolism, to examine this interaction at behavioral and biochemical levels. Our findings demonstrate that PNPO deficiency reduces alcohol aversion, increases consumption, and alters locomotor behavior. Biochemically, PNPO deficiency and alcohol exposure converge on amino acid metabolism, elevating inhibitory neurotransmitters GABA and glycine. Moreover, both PNPO deficiency and alcohol exposure lead to lethality with significant interaction, which can be rescued by VB6 supplementation. These results highlight a functional interaction between genetic VB6 deficiency and alcohol, suggesting potential therapeutic strategies for alcohol-related behaviors.

High ethanol preference and dissociated memory are co-occurring phenotypes associated with hippocampal GABAAR-δ receptor levels

Alcohol use disorder (AUD) frequently co-occurs with dissociative disorders and disorders with dissociative symptoms, suggesting a common neurobiological basis. It has been proposed that facilitated information processing under the influence of alcohol, resulting in the formation of dissociated memories, might be an important factor controlling alcohol use. Access to such memories is facilitated under the effect of alcohol, thus further reinforcing alcohol use. To interrogate possible mechanisms associated with these phenotypes, we used a mouse model of dissociative amnesia, combined with a high-alcohol preferring (HAP) model of AUD. Dissociated memory was induced by activation of hippocampal extrasynaptic GABA type A receptor delta subunits (GABA_AR-δ), which control tonic inhibition and to which ethanol binds with high affinity. Increased ethanol preference was associated with increased propensity to form dissociated memories dependent on GABA_AR-δ in the dorsal hippocampus (DH). Furthermore, the DH level of GABA_AR-δ protein, but not mRNA, was increased in HAP mice, and was inversely correlated to the level of miR-365-3p, suggesting an miRNA-mediated post-transcriptional mechanism contributing to elevated GABA_AR-δ. The observed changes of DH GABA_AR-δ were associated with a severe reduction of excitatory projections stemming from GABA_AR-δ-containing pyramidal neurons in the subiculum and terminating in the mammillary body. These results suggest that both molecular and circuit dysfunction involving hippocampal GABA_AR-δ receptors might contribute to the co-occurrence of ethanol preference and dissociated information processing.

Heritability of ethanol consumption and pharmacokinetics in a genetically diverse panel of collaborative cross mouse strains and their inbred founders

BACKGROUND

Interindividual variation in voluntary ethanol consumption and ethanol response is partially influenced by genetic variation. Discovery of the genes and allelic variants that affect these phenotypes may clarify the etiology and pathophysiology of problematic alcohol use, including alcohol use disorder. Genetically diverse mouse populations, which demonstrate heritable variation in ethanol consumption, can be utilized to discover the genes and gene networks that influence this trait. The Collaborative Cross (CC) recombinant inbred strains, Diversity Outbred (DO) population and their 8 founder strains are complementary mouse resources that capture substantial genetic diversity and can demonstrate expansive phenotypic variation in heritable traits. These populations may be utilized to discover candidate genes and gene networks that moderate ethanol consumption and other ethanol-related traits.

METHODS

We characterized ethanol consumption, preference, and pharmacokinetics in the 8 founder strains and 10 CC strains in 12-hour drinking sessions during the dark phase of the circadian cycle.

RESULTS

Ethanol consumption was substantially heritable, both early in ethanol access and over a chronic intermittent access schedule. Ethanol pharmacokinetics were also heritable; however, no association between strain-level ethanol consumption and pharmacokinetics was detected. The PWK/PhJ strain was the highest drinking strain, with consumption substantially exceeding that of the C57BL/6J strain, which is commonly used as a model of "high" or "binge" drinking. Notably, we found strong evidence that sex moderated genetic effects on voluntary ethanol drinking.

CONCLUSIONS

Collectively, this research serves as a foundation for expanded genetic study of ethanol consumption in the CC/DO and related populations. Moreover, we identified reference strains with extreme consumption phenotypes that effectively represent polygenic models of excessive ethanol use.

Systems genetics analysis of the LXS recombinant inbred mouse strains:Genetic and molecular insights into acute ethanol tolerance

We have been using the Inbred Long- and Short-Sleep mouse strains (ILS, ISS) and a recombinant inbred panel derived from them, the LXS, to investigate the genetic underpinnings of acute ethanol tolerance which is considered to be a risk factor for alcohol use disorders (AUDs). Here, we have used RNA-seq to examine the transcriptome of whole brain in 40 of the LXS strains 8 hours after a saline or ethanol "pretreatment" as in previous behavioral studies. Approximately 1/3 of the 14,184 expressed genes were significantly heritable and many were unique to the pretreatment. Several thousand cis- and trans-eQTLs were mapped; a portion of these also were unique to pretreatment. Ethanol pretreatment caused differential expression (DE) of 1,230 genes. Gene Ontology (GO) enrichment analysis suggested involvement in numerous biological processes including astrocyte differentiation, histone acetylation, mRNA splicing, and neuron projection development. Genetic correlation analysis identified hundreds of genes that were correlated to the behaviors. GO analysis indicated that these genes are involved in gene expression, chromosome organization, and protein transport, among others. The expression profiles of the DE genes and genes correlated to AFT in the ethanol pretreatment group (AFT-Et) were found to be similar to profiles of HDAC inhibitors. Hdac1, a cis-regulated gene that is located at the peak of a previously mapped QTL for AFT-Et, was correlated to 437 genes, most of which were also correlated to AFT-Et. GO analysis of these genes identified several enriched biological process terms including neuron-neuron synaptic transmission and potassium transport. In summary, the results suggest widespread genetic effects on gene expression, including effects that are pretreatment-specific. A number of candidate genes and biological functions were identified that could be mediating the behavioral responses. The most prominent of these was Hdac1 which may be regulating genes associated with glutamatergic signaling and potassium conductance.

Cross-species alcohol dependence-associated gene networks: Co-analysis of mouse brain gene expression and human genome-wide association data

Genome-wide association studies on alcohol dependence, by themselves, have yet to account for the estimated heritability of the disorder and provide incomplete mechanistic understanding of this complex trait. Integrating brain ethanol-responsive gene expression networks from model organisms with human genetic data on alcohol dependence could aid in identifying dependence-associated genes and functional networks in which they are involved. This study used a modification of the Edge-Weighted Dense Module Searching for genome-wide association studies (EW-dmGWAS) approach to co-analyze whole-genome gene expression data from ethanol-exposed mouse brain tissue, human protein-protein interaction databases and alcohol dependence-related genome-wide association studies. Results revealed novel ethanol-responsive and alcohol dependence-associated gene networks in prefrontal cortex, nucleus accumbens, and ventral tegmental area. Three of these networks were overrepresented with genome-wide association signals from an independent dataset. These networks were significantly overrepresented for gene ontology categories involving several mechanisms, including actin filament-based activity, transcript regulation, Wnt and Syndecan-mediated signaling, and ubiquitination. Together, these studies provide novel insight for brain mechanisms contributing to alcohol dependence.

Systems genetics of sensation seeking

Sensation seeking is a multifaceted, heritable trait which predicts the development of substance use and abuse in humans; similar phenomena have been observed in rodents. Genetic correlations among sensation seeking and substance use indicate shared biological mechanisms, but the genes and networks underlying these relationships remain elusive. Here, we used a systems genetics approach in the BXD recombinant inbred mouse panel to identify shared genetic mechanisms underlying substance use and preference for sensory stimuli, an intermediate phenotype of sensation seeking. Using the operant sensation seeking (OSS) paradigm, we quantified preference for sensory stimuli in 120 male and 127 female mice from 62 BXD strains and the C57BL/6J and DBA/2J founder strains. We used relative preference for the active and inactive levers to dissociate preference for sensory stimuli from locomotion and exploration phenotypes. We identified genomic regions on chromosome 4 (155.236‐155.742 Mb) and chromosome 13 (72.969‐89.423 Mb) associated with distinct behavioral components of OSS. Using publicly available behavioral data and mRNA expression data from brain regions involved in reward processing, we identified (a) genes within these behavioral QTL exhibiting genome‐wide significant cis‐eQTL and (b) genetic correlations among OSS phenotypes, ethanol phenotypes and mRNA expression. From these analyses, we nominated positional candidates for behavioral QTL associated with distinct OSS phenotypes including Gnb1 and Mef2c. Genetic covariation of Gnb1 expression, preference for sensory stimuli and multiple ethanol phenotypes suggest that heritable variation in Gnb1 expression in reward circuitry partially underlies the widely reported relationship between sensation seeking and substance use.

Predictive modeling of miRNA-mediated predisposition to alcohol-related phenotypes in mouse

BACKGROUND

MicroRNAs (miRNAs) are small non-coding RNAs that bind messenger RNAs and promote their degradation or repress their translation. There is increasing evidence of miRNAs playing an important role in alcohol related disorders. However, the role of miRNAs as mediators of the genetic effect on alcohol phenotypes is not fully understood. We conducted a high-throughput sequencing study to measure miRNA expression levels in alcohol naïve animals in the LXS panel of recombinant inbred (RI) mouse strains. We then combined the sequencing data with genotype data, microarry gene expression data, and data on alcohol-related behavioral phenotypes such as 'Drinking in the dark', 'Sleep time', and 'Low dose activation' from the same RI panel. SNP-miRNA-gene triplets with strong association within the triplet that were also associated with one of the 4 alcohol phenotypes were selected and a Bayesian network analysis was used to aggregate results into a directed network model.

RESULTS

We found several triplets with strong association within the triplet that were also associated with one of the alcohol phenotypes. The Bayesian network analysis found two networks where a miRNA mediates the genetic effect on the alcohol phenotype. The miRNAs were found to influence the expression of protein-coding genes, which in turn influences the quantitative phenotypes. The pathways in which these genes are enriched have been previously associated with alcohol-related traits.

CONCLUSION

This work enhances association studies by identifying miRNAs that may be mediating the association between genetic markers (SNPs) and the alcohol phenotypes. It suggests a mechanism of how genetic variants are affecting traits of interest through the modification of miRNA expression.

GABAA Receptor Subtype Mechanisms and the Abuse-Related Effects of Ethanol: Genetic and Pharmacological Evidence

Ethanol's reinforcing and subjective effects, as well as its ability to induce relapse, are powerful factors contributing to its widespread use and abuse. A significant mediator of these behavioral effects is the GABA_A receptor system. GABA_A receptors are the target for γ-aminobutyric acid (GABA), the major inhibitory neurotransmitter in the CNS. Structurally, they are pentameric, transmembrane chloride ion channels comprised of subunits from at least eight different families of distinct proteins. The contribution of different GABA_A subunits to ethanol's diverse abuse-related effects is not clear and remains an area of research focus. This chapter details the clinical and preclinical findings supporting roles for different α, β, γ, and δ subunit-containing GABA_A receptors in ethanol's reinforcing, subjective/discriminative stimulus, and relapse-inducing effects. The reinforcing properties of ethanol have been studied the most systematically, and convergent preclinical evidence suggests a key role for the α5 subunit in those effects. Regarding ethanol's subjective/discriminative stimulus effects, clinical and genetic findings support a primary role for the α2 subunit, whereas preclinical evidence implicates the α5 subunit. At present, too few studies investigating ethanol relapse exist to make any solid conclusions regarding the role of specific GABA_A subunits in this abuse-related effect.

From integrative genomics to systems genetics in the rat to link genotypes to phenotypes

Complementary to traditional gene mapping approaches used to identify the hereditary components of complex diseases, integrative genomics and systems genetics have emerged as powerful strategies to decipher the key genetic drivers of molecular pathways that underlie disease. Broadly speaking, integrative genomics aims to link cellular-level traits (such as mRNA expression) to the genome to identify their genetic determinants. With the characterization of several cellular-level traits within the same system, the integrative genomics approach evolved into a more comprehensive study design, called systems genetics, which aims to unravel the complex biological networks and pathways involved in disease, and in turn map their genetic control points. The first fully integrated systems genetics study was carried out in rats, and the results, which revealed conserved trans-acting genetic regulation of a pro-inflammatory network relevant to type 1 diabetes, were translated to humans. Many studies using different organisms subsequently stemmed from this example. The aim of this Review is to describe the most recent advances in the fields of integrative genomics and systems genetics applied in the rat, with a focus on studies of complex diseases ranging from inflammatory to cardiometabolic disorders. We aim to provide the genetics community with a comprehensive insight into how the systems genetics approach came to life, starting from the first integrative genomics strategies [such as expression quantitative trait loci (eQTLs) mapping] and concluding with the most sophisticated gene network-based analyses in multiple systems and disease states. Although not limited to studies that have been directly translated to humans, we will focus particularly on the successful investigations in the rat that have led to primary discoveries of genes and pathways relevant to human disease.

Genome characterization of the selected long- and short-sleep mouse lines

The Inbred Long- and Short-Sleep (ILS, ISS) mouse lines were selected for differences in acute ethanol sensitivity using the loss of righting response (LORR) as the selection trait. The lines show an over tenfold difference in LORR and, along with a recombinant inbred panel derived from them (the LXS), have been widely used to dissect the genetic underpinnings of acute ethanol sensitivity. Here we have sequenced the genomes of the ILS and ISS to investigate the DNA variants that contribute to their sensitivity difference. We identified ~2.7 million high-confidence SNPs and small indels and ~7000 structural variants between the lines; variants were found to occur in 6382 annotated genes. Using a hidden Markov model, we were able to reconstruct the genome-wide ancestry patterns of the eight inbred progenitor strains from which the ILS and ISS were derived, and found that quantitative trait loci that have been mapped for LORR were slightly enriched for DNA variants. Finally, by mapping and quantifying RNA-seq reads from the ILS and ISS to their strain-specific genomes rather than to the reference genome, we found a substantial improvement in a differential expression analysis between the lines. This work will help in identifying and characterizing the DNA sequence variants that contribute to the difference in ethanol sensitivity between the ILS and ISS and will also aid in accurate quantification of RNA-seq data generated from the LXS RIs.

Quantitative trait locus mapping of acute functional tolerance in the LXS recombinant inbred strains

BACKGROUND

We previously reported that acute functional tolerance (AFT) to the hypnotic effects of alcohol was significantly correlated with drinking in the dark (DID) in the LXS recombinant inbred panel, but only in mice that had been pretreated with alcohol. Here, we have conducted quantitative trait locus (QTL) mapping for AFT. DNA sequencing of the progenitor ILS and ISS strains and microarray analyses were also conducted to identify candidate genes and functional correlates.

METHODS

LXS mice were given either saline or alcohol (5 g/kg) on day 1 and then tested for loss of righting reflex AFT on day 2. QTLs were mapped using standard procedures. Two microarray analyses from brain were conducted: (i) naïve LXS mice and (ii) an alcohol treatment time course in the ILS and ISS. The full genomes of the ILS and ISS were sequenced to a depth of approximately 30×.

RESULTS

A significant QTL for AFT in the alcohol pretreatment group was mapped to distal chromosome 4; numerous suggestive QTLs were also mapped. Preference drinking and DID have previously been mapped to the chromosome 4 locus. The credible interval of the significant chromosome 4 QTL spanned 23 Mb and included 716 annotated genes of which 150 had at least 1 nonsynonymous single nucleotide polymorphism or small indel that differed between the ILS and ISS; expression of 48 of the genes was cis-regulated. Enrichment analysis indicated broad functional categories underlying AFT, including proteolysis, transcription regulation, chromatin modification, protein kinase activity, and apoptosis.

CONCLUSIONS

The chromosome 4 QTL is a key region containing possibly pleiotropic genes for AFT and drinking behavior. Given that the region contains many viable candidates and a large number of the genes in the interval fall into 1 or more of the enriched functional categories, we postulate that many genes of varying effect size contribute to the observed QTL effect.

Cytotoxicity effects of three-dimensional graphene in NIH-3T3 fibroblasts

We present an evaluation of the in vitro cytotoxicity of 3D graphene sheets fabricated by carbonization of polydopamine (PDA) films on a template of aligned nanopore arrays (NPAs) on a stainless steel surface.

Genes and Alcohol Consumption: Studies with Mutant Mice

In this chapter, we review the effects of global null mutant and overexpressing transgenic mouse lines on voluntary self-administration of alcohol. We examine approximately 200 publications pertaining to the effects of 155 mouse genes on alcohol consumption in different drinking models. The targeted genes vary in function and include neurotransmitter, ion channel, neuroimmune, and neuropeptide signaling systems. The alcohol self-administration models include operant conditioning, two- and four-bottle choice continuous and intermittent access, drinking in the dark limited access, chronic intermittent ethanol, and scheduled high alcohol consumption tests. Comparisons of different drinking models using the same mutant mice are potentially the most informative, and we will highlight those examples. More mutants have been tested for continuous two-bottle choice consumption than any other test; of the 137 mouse genes examined using this model, 97 (72%) altered drinking in at least one sex. Overall, the effects of genetic manipulations on alcohol drinking often depend on the sex of the mice, alcohol concentration and time of access, genetic background, as well as the drinking test.

The sequenced rat brain transcriptome--its use in identifying networks predisposing alcohol consumption

A quantitative genetic approach, which involves correlation of transcriptional networks with the phenotype in a recombinant inbred (RI) population and in selectively bred lines of rats, and determination of coinciding quantitative trait loci for gene expression and the trait of interest, has been applied in the present study. In this analysis, a novel approach was used that combined DNA-Seq data, data from brain exon array analysis of HXB/BXH RI rat strains and six pairs of rat lines selectively bred for high and low alcohol preference, and RNA-Seq data (including rat brain transcriptome reconstruction) to quantify transcript expression levels, generate co-expression modules and identify biological functions that contribute to the predisposition of consuming varying amounts of alcohol. A gene co-expression module was identified in the RI rat strains that contained both annotated and unannotated transcripts expressed in the brain, and was associated with alcohol consumption in the RI panel. This module was found to be enriched with differentially expressed genes from the selected lines of rats. The candidate genes within the module and differentially expressed genes between high and low drinking selected lines were associated with glia (microglia and astrocytes) and could be categorized as being related to immune function, energy metabolism and calcium homeostasis, as well as glial-neuronal communication. The results of the present study show that there are multiple combinations of genetic factors that can produce the same phenotypic outcome. Although no single gene accounts for predisposition to a particular level of alcohol consumption in every animal model, coordinated differential expression of subsets of genes in the identified pathways produce similar phenotypic outcomes.

DATABASE

The datasets supporting the results of the present study are available at http://phenogen.ucdenver.edu.

Influence of sex on genetic regulation of "drinking in the dark" alcohol consumption

The ILSXISS (LXS) recombinant inbred (RI) panel of mice is a valuable resource for genetic mapping studies of complex traits, due to its genetic diversity and large number of strains. Male and female mice from this panel were used to investigate genetic influences on alcohol consumption in the "drinking in the dark" (DID) model. Male mice (38 strains) and female mice (36 strains) were given access to 20% ethanol during the early phase of their circadian dark cycle for four consecutive days. The first principal component of alcohol consumption measures on days 2, 3, and 4 was used as a phenotype (DID phenotype) to calculate QTLs, using a SNP marker set for the LXS RI panel. Five QTLs were identified, three of which included a significant genotype by sex interaction, i.e., a significant genotype effect in males and not females. To investigate candidate genes associated with the DID phenotype, data from brain microarray analysis (Affymetrix Mouse Exon 1.0 ST Arrays) of male LXS RI strains were combined with RNA-Seq data (mouse brain transcriptome reconstruction) from the parental ILS and ISS strains in order to identify expressed mouse brain transcripts. Candidate genes were determined based on common eQTL and DID phenotype QTL regions and correlation of transcript expression levels with the DID phenotype. The resulting candidate genes (in particular, Arntl/Bmal1) focused attention on the influence of circadian regulation on the variation in the DID phenotype in this population of mice.

Genetics of gene expression characterizes response to selective breeding for alcohol preference

Numerous selective breeding experiments have been performed with rodents, in an attempt to understand the genetic basis for innate differences in preference for alcohol consumption. Quantitative trait locus (QTL) analysis has been used to determine regions of the genome that are associated with the behavioral difference in alcohol preference/consumption. Recent work suggests that differences in gene expression represent a major genetic basis for complex traits. Therefore, the QTLs are likely to harbor regulatory regions (eQTLs) for the differentially expressed genes that are associated with the trait. In this study, we examined brain gene expression differences over generations of selection of the third replicate lines of high and low alcohol-preferring (HAP3 and LAP3) mice, and determined regions of the genome that control the expression of these differentially expressed genes (de eQTLs). We also determined eQTL regions (rv eQTLs) for genes that showed a decrease in variance of expression levels over the course of selection. We postulated that de eQTLs that overlap with rv eQTLs, and also with phenotypic QTLs, represent genomic regions that are affected by the process of selection. These overlapping regions controlled the expression of candidate genes (that displayed differential expression and reduced variance of expression) for the predisposition to differences in alcohol consumption by the HAP3/LAP3 mice.

Genetic influences on brain gene expression in rats selected for tameness and aggression

Interindividual differences in many behaviors are partly due to genetic differences, but the identification of the genes and variants that influence behavior remains challenging. Here, we studied an F2 intercross of two outbred lines of rats selected for tame and aggressive behavior toward humans for >64 generations. By using a mapping approach that is able to identify genetic loci segregating within the lines, we identified four times more loci influencing tameness and aggression than by an approach that assumes fixation of causative alleles, suggesting that many causative loci were not driven to fixation by the selection. We used RNA sequencing in 150 F2 animals to identify hundreds of loci that influence brain gene expression. Several of these loci colocalize with tameness loci and may reflect the same genetic variants. Through analyses of correlations between allele effects on behavior and gene expression, differential expression between the tame and aggressive rat selection lines, and correlations between gene expression and tameness in F2 animals, we identify the genes Gltscr2, Lgi4, Zfp40, and Slc17a7 as candidate contributors to the strikingly different behavior of the tame and aggressive animals.

The multiMiR R package and database: integration of microRNA-target interactions along with their disease and drug associations

microRNAs (miRNAs) regulate expression by promoting degradation or repressing translation of target transcripts. miRNA target sites have been catalogued in databases based on experimental validation and computational prediction using various algorithms. Several online resources provide collections of multiple databases but need to be imported into other software, such as R, for processing, tabulation, graphing and computation. Currently available miRNA target site packages in R are limited in the number of databases, types of databases and flexibility. We present multiMiR, a new miRNA-target interaction R package and database, which includes several novel features not available in existing R packages: (i) compilation of nearly 50 million records in human and mouse from 14 different databases, more than any other collection; (ii) expansion of databases to those based on disease annotation and drug microRNAresponse, in addition to many experimental and computational databases; and (iii) user-defined cutoffs for predicted binding strength to provide the most confident selection. Case studies are reported on various biomedical applications including mouse models of alcohol consumption, studies of chronic obstructive pulmonary disease in human subjects, and human cell line models of bladder cancer metastasis. We also demonstrate how multiMiR was used to generate testable hypotheses that were pursued experimentally.

Identification of a QTL in Mus musculus for alcohol preference, withdrawal, and Ap3m2 expression using integrative functional genomics and precision genetics

Extensive genetic and genomic studies of the relationship between alcohol drinking preference and withdrawal severity have been performed using animal models. Data from multiple such publications and public data resources have been incorporated in the GeneWeaver database with >60,000 gene sets including 285 alcohol withdrawal and preference-related gene sets. Among these are evidence for positional candidates regulating these behaviors in overlapping quantitative trait loci (QTL) mapped in distinct mouse populations. Combinatorial integration of functional genomics experimental results revealed a single QTL positional candidate gene in one of the loci common to both preference and withdrawal. Functional validation studies in Ap3m2 knockout mice confirmed these relationships. Genetic validation involves confirming the existence of segregating polymorphisms that could account for the phenotypic effect. By exploiting recent advances in mouse genotyping, sequence, epigenetics, and phylogeny resources, we confirmed that Ap3m2 resides in an appropriately segregating genomic region. We have demonstrated genetic and alcohol-induced regulation of Ap3m2 expression. Although sequence analysis revealed no polymorphisms in the Ap3m2-coding region that could account for all phenotypic differences, there are several upstream SNPs that could. We have identified one of these to be an H3K4me3 site that exhibits strain differences in methylation. Thus, by making cross-species functional genomics readily computable we identified a common QTL candidate for two related bio-behavioral processes via functional evidence and demonstrate sufficiency of the genetic locus as a source of variation underlying two traits.

Genetics and genomics of alcohol sensitivity

Alcohol abuse and alcoholism incur a heavy socioeconomic cost in many countries. Both genetic and environmental factors contribute to variation in the inebriating effects of alcohol and alcohol addiction among individuals within and across populations. From a genetics perspective, alcohol sensitivity is a quantitative trait determined by the cumulative effects of multiple segregating genes and their interactions with the environment. This review summarizes insights from model organisms as well as human populations that represent our current understanding of the genetic and genomic underpinnings that govern alcohol metabolism and the sedative and addictive effects of alcohol on the nervous system.

Altered gamma-aminobutyric acid type B receptor subunit 1 splicing in alcoholics

BACKGROUND

Chronic alcohol exposure can change splice variant expression. The gamma-aminobutyric acid type B (GABAB) receptor undergoes splicing and is an alcoholism treatment target, but there is little information about splicing changes in this receptor in alcoholics. We studied GABAB receptor subunit 1 (GABAB1) splicing in alcoholic postmortem brains.

METHODS

To maximize GABAB1 splice junction identification, we combined gene specific libraries with RNA-seq. Splice junctions and mapped reads were also found from intronic and intergenic regions. We compared GABAB1 splice junctions in prefrontal cortices from 14 alcoholic and 15 control subjects and introduced new strategies, reads per kilobase of splice junction model per million mapped reads and reads per kilobase of gene model per million mapped reads, for quantitating splice junction and gene expression.

RESULTS

Novel splice junction detection indicated that the GABAB1 gene is at least two times longer than the previously reported gene length. GABAB1 exon and intron expression data showed low expression at the 5' end exons and exon grouping. This indicated that there are short splicing variants in addition to GABAB receptor subunit GABAB1a, the longest known major transcript. We found that chronic alcohol altered exon/intron expression and splice junction levels. Decreased expression of the gamma-aminobutyric acid binding site, a transmembrane domain and a microRNA binding site may decrease normal GABAB1 transcript population and thereby decrease normal signal transduction in alcoholics.

CONCLUSIONS

We discovered novel, complex splicing of GABAB1 in human brain and showed that chronic alcohol produces additional splicing complexity.

Progress in a replicated selection for elevated blood ethanol concentrations in HDID mice

Drinking in the dark (DID) is a limited access ethanol-drinking phenotype in mice. High Drinking in the Dark (HDID-1) mice have been bred for 27 selected generations (S27) for elevated blood ethanol concentrations (BECs) after a 4-h period of access to 20% ethanol. A second replicate line (HDID-2) was started later from the same founder population and is currently in S20. An initial report of response to selection in HDID-1 was published after S11. This article reports genetic and behavioral characteristics of both lines in comparison with the HS controls. Heritability is low in both replicates (h(2)  = 0.09) but the lines have shown 4-5 fold increases in BEC since S0; 80% of HDID-1 and 60% of HDID-2 mice reach BECs greater than 1.0 mg/ml. Several hours after a DID test, HDID mice show mild signs of withdrawal. Although not considered during selection, intake of ethanol (g/kg) during the DID test increased by approximately 80% in HDID-1 and 60% in HDID-2. Common genetic influences were more important than environmental influences in determining the similarity between BEC and intake for HDID mice. Analysis of the partitioning of intake showed that 60% of intake is concentrated in the last 2 h of the 4 h session. However, this has not changed during selection. Hourly BECs during the DID test reach peak levels after 3 or 4 h of drinking. HDID mice do not differ from HS mice in their rate of elimination of an acute dose of alcohol.

Introduction to sequencing the brain transcriptome

High-throughput next-generation sequencing is now entering its second decade. However, it was not until 2008 that the first report of sequencing the brain transcriptome appeared (Mortazavi, Williams, Mccue, Schaeffer, & Wold, 2008). These authors compared short-read RNA-Seq data for mouse whole brain with microarray results for the same sample and noted both the advantages and disadvantages of the RNA-Seq approach. While RNA-Seq provided exon level resolution, the majority of the reads were provided by a small proportion of highly expressed genes and the data analysis was exceedingly complex. Over the past 6 years, there have been substantial improvements in both RNA-Seq technology and data analysis. This volume contains 11 chapters that detail various aspects of sequencing the brain transcriptome. Some of the chapters are very methods driven, while others focus on the use of RNA-Seq to study such diverse areas as development, schizophrenia, and drug abuse. This chapter briefly reviews the transition from microarrays to RNA-Seq as the preferred method for analyzing the brain transcriptome. 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, and can be used to extract genotype information, e.g., nonsynonymous coding single nucleotide polymorphisms. RNA-Seq embraces the complexity of the brain transcriptome and provides a mechanism to understand the underlying regulatory code; the potential to inform the brain-behavior-disease relationships is substantial.

Hands-on workshops as an effective means of learning advanced technologies including genomics, proteomics and bioinformatics

Genomics and proteomics have emerged as key technologies in biomedical research, resulting in a surge of interest in training by investigators keen to incorporate these technologies into their research. At least two types of training can be envisioned in order to produce meaningful results, quality publications and successful grant applications: (1) immediate short-term training workshops and (2) long-term graduate education or visiting scientist programs. We aimed to fill the former need by providing a comprehensive hands-on training course in genomics, proteomics and informatics in a coherent, experimentally-based framework. This was accomplished through a National Heart, Lung, and Blood Institute (NHLBI)-sponsored 10-day Genomics and Proteomics Hands-on Workshop held at National Jewish Health (NJH) and the University of Colorado School of Medicine (UCD). The course content included comprehensive lectures and laboratories in mass spectrometry and genomics technologies, extensive hands-on experience with instrumentation and software, video demonstrations, optional workshops, online sessions, invited keynote speakers, and local and national guest faculty. Here we describe the detailed curriculum and present the results of short- and long-term evaluations from course attendees. Our educational program consistently received positive reviews from participants and had a substantial impact on grant writing and review, manuscript submissions and publications.

Gene expression within the extended amygdala of 5 pairs of rat lines selectively bred for high or low ethanol consumption

The objectives of this study were to determine innate differences in gene expression in 2 regions of the extended amygdala between 5 different pairs of lines of male rats selectively bred for high or low ethanol consumption: a) alcohol-preferring (P) vs. alcohol-non-preferring (NP) rats, b) high-alcohol-drinking (HAD) vs. low-alcohol-drinking (LAD) rats (replicate line-pairs 1 and 2), c) ALKO alcohol (AA) vs. nonalcohol (ANA) rats, and d) Sardinian alcohol-preferring (sP) vs. Sardinian alcohol-nonpreferring (sNP) rats, and then to determine if these differences are common across the line-pairs. Microarray analysis revealed up to 1772 unique named genes in the nucleus accumbens shell (AcbSh) and 494 unique named genes in the central nucleus of the amygdala (CeA) that significantly differed [False Discovery Rate (FDR) = 0.10; fold-change at least 1.2] in expression between the individual line-pairs. Analysis using Gene Ontology (GO) and Ingenuity Pathways information indicated significant categories and networks in common for up to 3 or 4 line-pairs, but not for all 5 line-pairs. However, there were almost no individual genes in common within these categories and networks. ANOVAs of the combined data for the 5 line-pairs indicated 1014 and 731 significant (p < 0.01) differences in expression of named genes in the AcbSh and CeA, respectively. There were 4-6 individual named genes that significantly differed across up to 3 line-pairs in both regions; only 1 gene (Gsta4 in the CeA) differed in as many as 4 line-pairs. Overall, the findings suggest that a) some biological categories or networks (e.g., cell-to-cell signaling, cellular stress response, cellular organization, etc.) may be in common for subsets of line-pairs within either the AcbSh or CeA, and b) regulation of different genes and/or combinations of multiple biological systems may be contributing to the disparate alcohol drinking behaviors of these line-pairs.

Positively correlated miRNA-mRNA regulatory networks in mouse frontal cortex during early stages of alcohol dependence

Background

Although the study of gene regulation via the action of specific microRNAs (miRNAs) has experienced a boom in recent years, the analysis of genome-wide interaction networks among miRNAs and respective targeted mRNAs has lagged behind. MicroRNAs simultaneously target many transcripts and fine-tune the expression of genes through cooperative/combinatorial targeting. Therefore, they have a large regulatory potential that could widely impact development and progression of diseases, as well as contribute unpredicted collateral effects due to their natural, pathophysiological, or treatment-induced modulation. We support the viewpoint that whole mirnome-transcriptome interaction analysis is required to better understand the mechanisms and potential consequences of miRNA regulation and/or deregulation in relevant biological models. In this study, we tested the hypotheses that ethanol consumption induces changes in miRNA-mRNA interaction networks in the mouse frontal cortex and that some of the changes observed in the mouse are equivalent to changes in similar brain regions from human alcoholics.

Results

miRNA-mRNA interaction networks responding to ethanol insult were identified by differential expression analysis and weighted gene coexpression network analysis (WGCNA). Important pathways (coexpressed modular networks detected by WGCNA) and hub genes central to the neuronal response to ethanol are highlighted, as well as key miRNAs that regulate these processes and therefore represent potential therapeutic targets for treating alcohol addiction. Importantly, we discovered a conserved signature of changing miRNAs between ethanol-treated mice and human alcoholics, which provides a valuable tool for future biomarker/diagnostic studies in humans. We report positively correlated miRNA-mRNA expression networks that suggest an adaptive, targeted miRNA response due to binge ethanol drinking.

Conclusions

This study provides new evidence for the role of miRNA regulation in brain homeostasis and sheds new light on current understanding of the development of alcohol dependence. To our knowledge this is the first report that activated expression of miRNAs correlates with activated expression of mRNAs rather than with mRNA downregulation in an in vivo model. We speculate that early activation of miRNAs designed to limit the effects of alcohol-induced genes may be an essential adaptive response during disease progression.

Transcriptome analysis of Inbred Long Sleep and Inbred Short Sleep mice

Many studies have utilized the Inbred Long Sleep and Inbred Short Sleep mouse strains to model the genetic influence on initial sensitivity to ethanol. The mechanisms underlying this divergent phenotype are still not completely understood. In this study, we attempt to identify genes that are differentially expressed between these two strains and to identify baseline networks of co-expressed genes, which may provide insight regarding their phenotypic differences. We examined the whole brain and striatal transcriptomes of both strains, using next generation RNA sequencing techniques. Many genes were differentially expressed between strains, including several in chromosomal regions previously shown to influence initial sensitivity to ethanol. These results are in concordance with a similar sample of striatal transcriptomes measured using microarrays. In addition to the higher dynamic range, RNA-Seq is not hindered by high background noise or polymorphisms in probesets as with microarray technology, and we are able to analyze exome sequence of abundant genes. Furthermore, utilizing Weighted Gene Co-expression Network Analysis, we identified several modules of co-expressed genes corresponding to strain differences. Several candidate genes were identified, including protein phosphatase 1 regulatory unit 1b (Ppp1r1b), prodynorphin (Pdyn), proenkephalin (Penk), ras association (RalGDS/AF-6) domain family member 2 (Rassf2), myosin 1d (Myo1d) and transthyretin (Ttr). In addition, we propose a role for potassium channel activity as well as map kinase signaling in the observed phenotypic differences between the two strains.

Genetic studies of acute tolerance, rapid tolerance, and drinking in the dark in the LXS recombinant inbred strains

BACKGROUND

We hypothesized that rapid tolerance (1-day tolerance) for the duration of the loss of righting reflex ("sleep time" [ST]) was mediated by an increase in acute functional tolerance (AFT). We also hypothesized that increased AFT would correspond to increased drinking. These questions were addressed using the LXS recombinant inbred mouse strain panel.

METHODS

Mice were given a pretreatment dose of either saline or 5 g/kg alcohol on day 1. On day 2, mice were tested for ST (4.1 g/kg) using a method with which it is possible to accurately assess AFT. Genetic correlation analysis was conducted among the ST-related variables and also with "drinking in the dark" (DID) which was previously measured by Saba and colleagues (2011).

RESULTS

Saline-pretreated mice showed a continuous distribution of ST ranging from ~40 minutes to over 3 hours. Of the 43 strains tested, 9 showed significantly decreased ST after alcohol pretreatment, while in 3 strains, ST was significantly increased. AFT scores ranged from 0 to over 200 mg% in the saline group, and in the alcohol group, 8 strains showed a significant increase in AFT and 2 strains showed significant decrease in AFT. In the saline group, AFT was significantly correlated with ST (r = -0.47), but not in the alcohol group (r = -0.22). DID was significantly correlated with only AFT in the alcohol pretreated group (r = 0.64).

CONCLUSIONS

The results suggest that AFT is an important component of the overall ST response, but that the alcohol pretreatment-induced change in AFT does not contribute to rapid ST tolerance. The significant correlation between DID and AFT in the alcohol group suggests that AFT may be a more relevant predictor of drinking behavior than the static measurement of ST. Moreover, preexposure to alcohol seems to change AFT in a way that makes it an even stronger predictor of drinking behavior.

Selection for drinking in the dark alters brain gene coexpression networks

BACKGROUND

Heterogeneous stock (HS/NPT) mice have been used to create lines selectively bred in replicate for elevated drinking in the dark (DID). Both selected lines routinely reach a blood ethanol (EtOH) concentration (BEC) of 1.00 mg/ml or greater at the end of the 4-hour period of access in Day 2. The mechanisms through which genetic differences influence DID are currently unclear. Therefore, the current study examines the transcriptome, the first stage at which genetic variability affects neurobiology. Rather than focusing solely on differential expression (DE), we also examine changes in the ways that gene transcripts collectively interact with each other, as revealed by changes in coexpression patterns.

METHODS

Naïve mice (N = 48/group) were genotyped using the Mouse Universal Genotyping Array, which provided 3,683 informative markers. Quantitative trait locus (QTL) analysis used a marker-by-marker strategy with the threshold for a significant logarithm of odds (LOD) set at 10.6. Gene expression in the ventral striatum was measured using the Illumina Mouse 8.2 array. Differential gene expression and the weighted gene coexpression network analysis (WGCNA) were implemented largely as described elsewhere.

RESULTS

Significant QTLs for elevated BECs after DID were detected on chromosomes 4, 14, and 16; the latter 2 were associated with gene-poor regions. None of the QTLs overlapped with known QTLs for EtOH preference drinking. Ninety-four transcripts were detected as being differentially expressed in both selected lines versus HS controls; there was no overlap with known preference genes. The WGCNA revealed 2 modules as showing significant effects of both selections on intramodular connectivity. A number of genes known to be associated with EtOH phenotypes (e.g., Gabrg1, Glra2, Grik1, Npy2r, and Nts) showed significant changes in connectivity.

CONCLUSIONS

We found marked and consistent effects of selection on coexpression patterns; DE changes were more modest and less concordant. The QTLs and differentially expressed genes detected here are distinct from the preference phenotype. This is consistent with behavioral data and suggests that the DID and preference phenotypes are markedly different genetically.

Genes, behavior and next-generation RNA sequencing

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.

Modeling the diagnostic criteria for alcohol dependence with genetic animal models

A diagnosis of alcohol dependence (AD) using the DSM-IV-R is categorical, based on an individual's manifestation of three or more symptoms from a list of seven. AD risk can be traced to both genetic and environmental sources. Most genetic studies of AD risk implicitly assume that an AD diagnosis represents a single underlying genetic factor. We recently found that the criteria for an AD diagnosis represent three somewhat distinct genetic paths to individual risk. Specifically, heavy use and tolerance versus withdrawal and continued use despite problems reflected separate genetic factors. However, some data suggest that genetic risk for AD is adequately described with a single underlying genetic risk factor. Rodent animal models for alcohol-related phenotypes typically target discrete aspects of the complex human AD diagnosis. Here, we review the literature derived from genetic animal models in an attempt to determine whether they support a single-factor or multiple-factor genetic structure. We conclude that there is modest support in the animal literature that alcohol tolerance and withdrawal reflect distinct genetic risk factors, in agreement with our human data. We suggest areas where more research could clarify this attempt to align the rodent and human data.

The hunt for gene effects pertinent to behavioral traits and psychiatric disorders: from mouse to human

The field of behavioral genetics was reviewed in the classic 1960 text by Fuller and Thompson. Since then, there has been remarkable progress in the genetic analysis of animal behavior. Many molecular genetic methods in common use today were not even anticipated in 1960. Animal models for many human psychiatric disorders have been discovered or created. In human behavior genetics, however, powerful new methods have failed to reveal even one bona fide, replicable gene effect pertinent to the normal range of variation in intelligence and personality. There is no explanatory or predictive value in that genetic information. For several psychiatric disorders, including autism and schizophrenia, many large genetic effects arise from de novo mutations. Genetically, the disorders are heterogeneous; different cases with the same diagnosis have different causes. The promises of the molecular genetic revolution have not been fulfilled in behavioral domains of most interest to human psychology.

Introduction to the special issue "Pharmacotherapies for the treatment of alcohol abuse and dependence" and a summary of patents targeting other neurotransmitter systems

This paper introduces the Special Section: Pharmacotherapies for the Treatment of Alcohol Abuse and Dependence and provides a summary of patents targeting neurotransmitter systems not covered in the other four chapters. The World Health Organization notes that alcoholic-type drinking results in 2.5 million deaths per year, and these deaths occur to a disproportionately greater extent among adolescents and young adults. Developing a pharmacological treatment targeting alcohol abuse and dependence is complicated by (a) the heterogeneous nature of the disease(s), (b) alcohol affecting multiple neurotransmitter and neuromodulator systems, and (c) alcohol affecting multiple organ systems which in turn influence the function of the central nervous system. Presently, the USA Federal Drug Administration has approved three pharmacotherapies for alcoholism: disulfiram, naltrexone, and acamprosate. This chapter provides a summary of the following systems, which are not covered in the accompanying chapters; alcohol and acetaldehyde metabolism, opioid, glycinergic, GABA-A, neurosteroid, dopaminergic, serotonergic, and endocannabinoid, as well as patents targeting these systems for the treatment of alcoholism. Finally, an overview is presented on the use of pharmacogenetics and pharmacogenomics in tailoring treatments for certain subpopulations of alcoholics, which is expected to continue in the future.

A gene co-expression network in whole blood of schizophrenia patients is independent of antipsychotic-use and enriched for brain-expressed genes

Despite large-scale genome-wide association studies (GWAS), the underlying genes for schizophrenia are largely unknown. Additional approaches are therefore required to identify the genetic background of this disorder. Here we report findings from a large gene expression study in peripheral blood of schizophrenia patients and controls. We applied a systems biology approach to genome-wide expression data from whole blood of 92 medicated and 29 antipsychotic-free schizophrenia patients and 118 healthy controls. We show that gene expression profiling in whole blood can identify twelve large gene co-expression modules associated with schizophrenia. Several of these disease related modules are likely to reflect expression changes due to antipsychotic medication. However, two of the disease modules could be replicated in an independent second data set involving antipsychotic-free patients and controls. One of these robustly defined disease modules is significantly enriched with brain-expressed genes and with genetic variants that were implicated in a GWAS study, which could imply a causal role in schizophrenia etiology. The most highly connected intramodular hub gene in this module (ABCF1), is located in, and regulated by the major histocompatibility (MHC) complex, which is intriguing in light of the fact that common allelic variants from the MHC region have been implicated in schizophrenia. This suggests that the MHC increases schizophrenia susceptibility via altered gene expression of regulatory genes in this network.

Understanding Alcoholism Through microRNA Signatures in Brains of Human Alcoholics

Advances in the fields of genomics and genetics in the last decade have identified a large number of genes that can potentially influence alcohol-drinking behavior in humans as well as animal models. Consequently, the task of identifying efficient molecular targets that could be used to develop effective therapeutics against the disease has become increasingly daunting. One of the reasons for this is the fact that each of the many alcohol-responsive genes only contributes a small effect to the overall mechanism and disease phenotype, as is characteristic of complex traits. Current research trends are hence shifting toward the analysis of gene networks rather than emphasizing individual genes. The discovery of microRNAs and their mechanisms of action on regulation of transcript level and protein translation have made evident the utility of these small non-coding RNA molecules that act as central coordinators of multiple cross-communicating cellular pathways. Cells exploit the fact that a single microRNA can target hundreds of mRNA transcripts and that a single mRNA transcript can be simultaneously targeted by distinct microRNAs, to ensure fine-tuned and/or redundant control over a large number of cellular functions. By the same token, we can use these properties of microRNAs to develop novel, targeted strategies to combat complex disorders. In this review, we will focus on recent discoveries of microRNA signatures in brain of human alcoholics supporting the hypothesis that changes in gene expression and regulation by microRNAs are responsible for long-term neuroadaptations occurring during development of alcoholism. We also discuss insights into the potential modulation of epigenetic regulators by a subset of microRNAs. Taken together, microRNA activity may be controlling many of the cellular mechanisms already known to be involved in the development of alcoholism, and suggests potential targets for the development of novel therapeutic interventions.

GABAergic gene expression in postmortem hippocampus from alcoholics and cocaine addicts; corresponding findings in alcohol-naïve P and NP rats

BACKGROUND

By performing identical studies in humans and rats, we attempted to distinguish vulnerability factors for addiction from neurobiological effects of chronic drug exposure. We focused on the GABAergic system within the hippocampus, a brain region that is a constituent of the memory/conditioning neuronal circuitry of addiction that is considered to be important in drug reinforcement behaviors in animals and craving and relapse in humans.

METHODOLOGY

Using RNA-Seq we quantified mRNA transcripts in postmortem total hippocampus from alcoholics, cocaine addicts and controls and also from alcohol-naïve, alcohol preferring (P) and non-preferring (NP) rats selectively bred for extremes of alcohol-seeking behavior that also show a general addictive tendency. A pathway-targeted analysis of 25 GABAergic genes encoding proteins implicated in GABA synthesis, metabolism, synaptic transmission and re-uptake was undertaken.

PRINCIPAL FINDINGS

Directionally consistent and biologically plausible overlapping and specific changes were detected: 14/25 of the human genes and 12/25 of the rat genes showed nominally significant differences in gene expression (global p values: 9×10⁻¹⁴, 7×10⁻¹¹ respectively). Principal FDR-corrected findings were that GABBR1 was down-regulated in alcoholics, cocaine addicts and P rats with congruent findings in NSF, implicated in GABAB signaling efficacy, potentially resulting in increased synaptic GABA. GABRG2, encoding the gamma2 subunit required for postsynaptic clustering of GABAA receptors together with GPHN, encoding the associated scaffolding protein gephryin, were both down-regulated in alcoholics and cocaine addicts but were both up-regulated in P rats. There were also expression changes specific to cocaine addicts (GAD1, GAD2), alcoholics (GABRA2) and P rats (ABAT, GABRG3).

CONCLUSIONS/SIGNIFICANCE

Our study confirms the involvement of the GABAergic system in alcoholism but also reveals a hippocampal GABA input in cocaine addiction. Congruent findings in human addicts and P rats provide clues to predisposing factors for alcohol and drug addiction. Finally, the results of this study have therapeutic implications.

MicroRNAs in addiction: adaptation's middlemen?

A central question in addiction is how drug-induced changes in synaptic signaling are converted into long-term neuroadaptations. Emerging evidence reveals that microRNAs (miRNAs) have a distinct role in this process through rapid response to cellular signals and dynamic regulation of local mRNA transcripts. Because each miRNA can target hundreds of mRNAs, relative changes in the expression of miRNAs can greatly impact cellular responsiveness, synaptic plasticity and transcriptional events. These diverse consequences of miRNA action occur through coordination with genes implicated in addictions, the most compelling of these being the neurotrophin BDNF, the transcription factor cAMP-responsive element-binding protein (CREB) and the DNA-binding methyl CpG binding protein 2 (MeCP2). In this study, we review the recent progress in the understanding of miRNAs in general mechanisms of plasticity and neuroadaptation and then focus on specific examples of miRNA regulation in the context of addiction. We conclude that miRNA-mediated gene regulation is a conserved means of converting environmental signals into neuronal response, which holds significant implications for addiction and other psychiatric illnesses.