Positively correlated miRNA-mRNA regulatory networks in mouse frontal cortex during early stages of alcohol dependence
© Nunez et al.; licensee BioMed Central Ltd. 2013
Received: 10 May 2013
Accepted: 4 October 2013
Published: 22 October 2013
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.
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.
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.
Alcohol dependence is a brain disorder of complex etiology characterized by brain-wide pathophysiological alterations under continuous crosstalk, with approximately 40-60% contribution from genetic factors . As suggested by Farris and Miles, interacting neuronal and glial networks across distinct brain regions likely control the variety of alcohol endophenotypes, which are ultimately controlled by the regulation of multiple gene networks expressed within individual neurons or glial cells . Support for these statements is found in the observation that chronic consumption of alcohol induces long-term changes in brain gene and protein expression, which allow brain cells to adapt through homeostatic alterations in distinct signaling pathways [3, 4]. Although the complex regulatory mechanisms governing these changes are not fully understood, there is mounting evidence that transcriptional reprogramming is brain area-specific and may reflect both pre-existing differences in gene expression and alterations in response to alcohol consumption .
MicroRNAs (miRNAs), a type of short non-coding RNA with well characterized post-transcriptional regulatory functions, have been recently implicated in cellular responses to multiple drugs of abuse including alcohol [4, 6–10]. Moreover, miRNAs have extensive regulatory capacity given that a single miRNA can simultaneously target multiple genes, and multiple miRNAs can cooperatively function while targeting a single gene, therefore allowing for fine-tuned regulation of targeted gene expression. We have previously shown that miRNAs, which localize to and display neurotransmitter signaling-related activities at neuronal synapses, are also capable of eliciting distinct and specific activities in other cell types and/or compartments of the cell, such as innate immunity- and epigenetic-related functions in neuronal and glial cells . This underscores the impact that molecular efficiency, signaling crosstalk, and cellular economy play in the adaptation and evolution of cellular systems and highlights the significance of miRNAs as efficient molecular targets for complex diseases such as addiction.
We have previously reported that in brains of human alcoholics, miRNAs appear mostly upregulated, with an overrepresented number of targets among downregulated mRNAs. Approximately 80% of those targets appeared to be combinatorially regulated by multiple miRNAs . Many of these upregulated miRNAs appear to be central regulators of epigenetic-, synaptic signaling-, and neuroimmune-related processes (Nunez and Mayfield ). Here we describe new findings in the frontal cerebral cortex (FCtx) of a mouse model for high voluntary ethanol consumption. Mouse models like ours, implementing drinking in the dark (DID) protocols that produce pharmacologically relevant levels of ethanol in blood, have been developed to examine binge-like ethanol consumption . Binge ethanol consumption in humans is “a pattern of behaviour that may emerge prior to, and contribute to the development of, ethanol dependence” . Our goal is to identify sets of miRNAs and relevant regulatory networks that contribute to our understanding of the mechanisms related to alcohol dependence in humans. We tested the hypotheses that alcohol induces changes in miRNA-mRNA interaction networks in mouse frontal cortex and that some of the changes observed in the mouse brain are equivalent to changes occurring in similar brain regions from human alcoholics. We report on the integrative analysis of genome-wide miRNA-mRNA expression profiles from mouse FCtx samples and on the implementation of weighted gene coexpression network analysis (WGCNA). Using WGCNA, we generated correlations among coexpressed gene modules and differentially expressed miRNAs and provided evidence for relevant regulatory networks responsive to alcohol actions. We unexpectedly uncovered striking positive correlation patterns between upregulated miRNAs and upregulated predicted, validated, and/or correlated mRNA targets.
Differential expression analysis underscores a major role for upregulation of miRNAs in FCtx of ethanol-treated mice
Differentially expressed miRNAs in mouse frontal cortex
Ethanol-responsive genes in brains of mice significantly match respective homologs differentially expressed in brains of human alcoholics
Integrated miRNA-mRNA networks reveal positively correlated action of brain miRNAs and persistent upregulation of overtargeted mRNAs in ethanol drinking mice
In addition when overrepresentation analyses were conducted on the interaction networks, we found that the majority of upregulated miRNAs (38/48 = 79%) have an overrepresented number of targets among upregulated genes (Additional file 1: Table S2, FDR < 0.10), while only about one quarter (10/41 = 24%) had an overrepresented number of targets among downregulated genes (Additional file 1: Table S3, FDR < 0.10). The expected number of miRNAs with overrepresented number of targets is 22/45 = 49%, as determined by averaging the proportions obtained for two random networks. Therefore, the proportion detected in the positively correlated (upregulated miRNA-upregulated mRNA) network is highly significantly enriched (79% vs. 49%, Chi-squared P < 1 × 10-4) in miRNAs that have a higher than expected number of targets (“overtargeting” miRNAs). On the other hand, the proportion detected in the negatively correlated (upregulated miRNA-downregulated mRNA) network is significantly smaller than expected by chance (24% vs. 49%, Chi-squared P = 0.0016), indicating that the negatively correlated network is depleted of overtargeting miRNAs. This is consistent with the observation that out of 413 upregulated genes, 265 (64%) represented putative upregulated miRNA targets, while only 119 out of 296 downregulated genes (40%) represented putative downregulated miRNA targets. These results suggest that the function of the alcohol-induced miRNAs is to target alcohol-induced genes.
WGCNA identifies brain gene modules perturbed by ethanol
Modular enrichment analysis
Module correlation to EtOH consumption
P value of correlation to consumption
Differentially expressed (DE) genes in module
P value of enrichment in DE genes
Functional category enrichment
Cell type-specific enrichment
Vesicle-mediated transport, Endocytosis
Transport, Synaptic transmission, Toll-like receptor & chemokine signaling
Cellular component organization or biogenesis, endosomal transport
Synaptic transmission, Neurotransmitter secretion
Neurons, Astrocytes, Oligodendrocytes
Regulation of developmental process, Neurogenesis
Astrocytes, Neurons, Oligodendrocytes
Blood vessel development, Angiogenesis, BMP signaling
Oxidative phosphorylation, Ribosome biogenesis
Synaptic transmission, Chemokine signaling, Toll-like receptor signaling
Electron transport chain, Translation, Ribosome biogenesis
A major theme emerging from these analyses is that the alcohol-responsive modules are enriched in genes involved in synaptic transmission and cytokine signaling pathways (Table 2, Additional file 1: Table S5-S7). This underscores the importance of these pathways in the transduction of alcohol actions. These modules showed enrichment for cell-type specific genes present in neurons, astrocytes, and oligodendrocytes, but not in microglia-specific genes, which highlight a major effect of alcohol on cells of the neuronal and macroglial lineages. These lineages share a common origin in precursor cells derived from the embryonic germ layer known as the neuroectoderm, which differ from the hematopoietic origin of microglia .
To assess potential network topology changes impinged by ethanol consumption, we conducted an alternative network analysis where we implemented WGCNA on each animal group independently. Calculation of the overlap between coexpressed gene modules from both networks identified significant differences in the modular organization of the independent networks (Additional file 1: Figure S4). Panel C in Additional file 1: Figure S4 reflects the reorganization of specific modules as the network “evolves” from a non-treated stage (network from control group) into an ethanol-treated stage (network from ethanol-treated group). For example, the control C_turquoise module diverted into two distinct modules in the alcohol network (A_brown and A_green). In addition, the alcohol (A_green) module included genes that had originally segregated into three distinct modules in the control network (C_turquoise, C greenyellow, and C_green). Our miRNA-mRNA integrative analysis suggests that such gene network reorganization events could be mediated, at least in part, by the activation of alcohol-responsive miRNAs.
Core of interactions in ethanol-relevant, miRNA-targeted network modules highlights important regulators of alcohol action (hub genes)
In the brown module, which also contains alcohol-upregulated genes, protein phosphatase Mg(2+)/Mn(2+)-dependent 1A (Ppm1a) and heat shock 70kDa protein 8 (Hspa8) are the main hub genes. In addition, another heat shock protein (Hsp90ab1), GDP dissociation inhibitor 1 (Gdi1), and the glycine receptor beta (Glrb) (Additional file 1: Figure S5A) were among the highly interconnected genes in the module. The MAPK signaling pathway was the most significantly enriched pathway in this module (including genes Hspa8, Ppm1a, Ppp3ca, Ppp5c, Fgf13, Rasgrf1, Akt3, Mapk1, Mapk9; FDR = 1.2 × 10-3 Additional file 1: Table S6). MicroRNA miR-34c-5p appeared again as a main modulator of hub genes in this module, namely Ppm1a and Hsp90ab1, as detected by significant correlation among respective miRNA and mRNA expression profiles. The central role of miR-34c-5p while targeting additional genes in the brown module is again evident in the expression correlation network shown in Additional file 1: Figure S5B (note the higher connectivity of this particular miRNA).
In the turquoise module (Additional file 1: Figure S6A) containing alcohol-downregulated genes, we detected a larger core of hub genes known to be enriched in synaptic transmission, endocytosis, and cytokine signaling (Additional file 1: Table S7). The CXCR4-mediated signaling pathway appeared as one of the most enriched. Concordantly, the CXCR4 gene was previously reported as one of the genes more significantly overtargeted by miRNAs in PFC of human alcoholics . Regarding miRNA regulation of this module, the central role of miR-34c-5p and let-7g-5p is yet again evident in the expression correlation network shown in Additional file 1: Figure S6B (note the higher connectivity of these particular miRNAs).
The goal of this study was to gain insights into the molecular mechanisms that are induced by alcohol consumption and governed by miRNA regulation in an early stage of development of dependence. This would allow for better understanding of how molecular progression of the disease occurs. Using a well-characterized mouse model developed in our lab [25, 26], we evaluated whole transcriptome profiles of miRNA and mRNA expression from the same samples (“paired profiles”).
We found that a significant number of miRNAs that are upregulated in the FCtx of ethanol-treated mice are also upregulated in the PFC of human alcoholics [4, 11]. Similarly, we found that many genes are differentially expressed in the FCtx of the treated animals, and a significant number of these are also differentially expressed in PFC of human alcoholics [14, 15]. Interestingly, about 40-80% of the common differentially expressed genes change in opposite directions when the human and mouse models are compared, with a majority of differentially expressed genes being upregulated in the ethanol-treated mouse brains but downregulated in the brains of alcoholics. Our results demonstrate first, a high degree of conservation of alcohol-relevant genes and pathways from mouse to human and second, unique differences in the direction of gene expression changes between the two models. Such distinctions may be a consequence of the duration of ethanol exposure, which lasted for most of the adult life in the case of the human alcoholics (long term adaptations), but only 20 days during early adulthood in the case of our treated mice (early stages of adaptation). These differences suggest a dynamic gene regulatory network that changes as subjects advance towards alcohol dependence. We acknowledge that there may be undetermined effects that influence our comparison of results across the two species: 1) our expression analysis in mice represents a snapshot at 24 hours after the drinking protocol was completed; therefore, the animals could have been experiencing acute withdrawal; while 2) the post-mortem human brain samples originated from human alcoholics that were considered active heavy drinkers until the time of death, thus negating effects of acute withdrawal in this group. It would be beneficial to conduct additional expression studies after longer drinking protocols to determine a time sequence of expression changes in the alcohol-treated mouse as it progresses towards long-term ethanol consumption. It is also important to conduct expression profiling after longer time points of ethanol withdrawal, to determine whether withdrawal may influence the direction of mRNA and miRNA expression changes in the mouse.
The mouse miRNA-mRNA interaction networks presented here provide evidence for some striking features in the alcohol-specific network topologies. Contrary to expectations and to what happens in the brain of chronic human alcoholics, we did not find differences between the negatively correlated interaction network (upregulated miRNAs targeting downregulated mRNAs) and a random network, which indicates no selective pressure under the alcohol treatment for this type of interaction. Interestingly, we discovered that the positively correlated network (upregulated miRNAs targeting upregulated mRNAs) displayed an unexpected behavior, being twice as interconnected than expected by chance, and with a highly dense core of interactions. When overrepresentation of miRNA targets was evaluated, we found that a majority of upregulated miRNAs (over 80%) have a significantly higher number of upregulated targets than is expected by chance (enrichment for over-targeting), while only 24% of the miRNAs displayed the behavior towards downregulated genes (depletion for over-targeting). We speculate that miRNA targeting of upregulated genes might be in an “uncompensated" state in our mouse model, where miRNAs have not yet been able to downregulate the expression of their relevant targets. The extended and highly statistically significant over-targeting of upregulated genes in general, might indicate a signal for transcript downregulation. We speculate that a selective pressure to target ethanol-induced (upregulated) genes is impinged by chronic, binge alcohol consumption in an apparent adaptive response from the cellular environment to recover homeostasis. To our knowledge, this is the first report identifying positively correlated miRNA-mRNA networks in an animal model, which we propose as an adaptive mechanism to reinstate cellular homeostasis. Our current working model suggests that alcohol drinking activates a variety of brain genes and pathways that are consequently targeted for down-regulation by mature miRNAs that are activated in response to the aforementioned gene activation. Over-targeting of red and brown module genes by several miRNAs described in this work are examples of the behavior. Our assumption is that activated miRNAs in the mouse frontal cortex had not yet been able to overcome the induced upregulation of alcohol-responsive genes, at the experimental time point chosen to collect the mouse brains and conduct the expression profiling studies. We also suggest, based on expression patterns in similar brain regions from human alcoholics, that miRNAs eventually dominate the regulatory interactions with targeted genes and consequently induce a downregulated state in those alcohol-responsive gene.
These results underscore the biological importance of targeting alcohol-induced genes by alcohol-induced miRNAs and suggest that induction of miRNA targeting in the frontal cortex is a conserved (from rodent to human), adaptive response to counter an initial activation of transcription induced by high ethanol consumption. One possible explanation for the unexpected correlation (between upregulated miRNAs and upregulated target mRNAs) in the early stages of adaptation (binge mouse model) is the temporal inability of miRNA expression levels to cope with the levels of expression of target mRNAs. It is unlikely that indirect activities (e.g., activation of certain ethanol-responsive transcription factors that concomitantly activate certain miRNAs and respective target genes) are responsible for such over-targeting of upregulated mRNAs by a majority of the upregulated miRNAs.
The coexpression network analysis underscores the main biological processes affected by alcohol and miRNA regulation in the mouse model, such as synaptic vesicle-mediated transport and endocytic recycling, as well as neuroimmune signaling mediated by chemokine and Toll-like receptor signaling. The fact that signaling pathways that are upregulated by alcohol and over-targeted by upregulated miRNAs also show up enriched in downregulated modules, confirms the importance of these pathways in transducing alcohol’s actions during early stages of chronic alcohol consumption and their role as triggers of a miRNA adaptive response aimed at counteracting such activation. MicroRNAs miR-34c-5p and let-7g-5p, in particular, appear to be central regulators of hub genes in the most significant ethanol-responsive gene modules. Similarly, these two miRNAs were also upregulated in PFC of human alcoholics [4, 11], and a let-7 family member was demonstrated to alter cocaine-induced conditioned place preference behavior when manipulated in the nucleus accumbens of rats . We are enthusiastic about future prospects of using let-7 modulation in mouse brain to induce behavioural changes related to alcohol consumption and development of dependence.
One of the miRNA-targeted hub genes highlighted by the coexpression network analysis was synaptotagmin 11 (red module). The synaptotagmin family members are known to function by cooperating with SNARE proteins and accelerating membrane fusion [28, 29]. Several other genes in the red module, such as sorting nexin-17 (Snx17), target of myb1 (Tom1), and bridging integrator 1 (Bin1), are also involved in membrane trafficking and sorting processes relevant to synaptic transmission as well as the recycling of adaptor components involved in signaling of Toll receptors and the interleukin IL-1 receptor [21–24, 30]. These results underscore the relevance of synaptic vesicle transport and endocytic recycling in the transduction of alcohol’s actions. Our group has previously gathered proteomics/genomics evidence that support the role of these pathways [31–33].
Among the hub genes in the brown module, Ppm1a is strongly expressed in the brain and involved in the regulation of the stress response by acting on mitogen-activated protein (MAP) kinase signaling and by regulating the immune response [34, 35]. Ethanol is well known to differentially modulate MAP kinase signaling cascades depending on the cell type, brain region, and ethanol treatment paradigm [3, 36–42]. PPM1A phosphatase has been reported to play an important role in the termination of TNFα-mediated NFκB activation through dephosphorylation and inactivation of IKKβ . Accordingly, the NF-κB pathway has been widely implicated in the transduction of alcohol effects, and activation of the pathway is associated with increased ethanol consumption [43–49]. Stress-induced chaperones Hspa8 (an Hsp70 family member) and Hsp90ab were also among brown module’s highly connected genes. These chaperons exert a variety of cellular functions including proteostasis maintenance [50, 51] and endotoxin-like effects through Toll-like receptor 4 (TLR4) [52, 53]. Upregulation of genes encoding several classes of chaperons has been previously reported after ethanol treatment [54–56]. Hsp70 and Hsp90 have also been reported as co-clustering factors of the LPS-sensing complex, which also contains CD14, CXCR4, GDF5, and TLR4 . Importantly, LPS-triggered signaling through this complex has been recently reported to increase ethanol consumption and to alter certain aspects of alcohol reward and aversion in mice . Glrb, another hub gene in the brown module, codes for the beta subunit of the glycine receptor (a neurotransmitter-gated pentameric ion channel) that is essential for targeting the receptor to the synapse. Glycine receptors contain an integral chloride channel and are known to have modulatory sites for anesthetics, neurosteroids, inhaled drugs of abuse, and ethanol [58–61]. Interestingly, α3 subunit-containing glycine receptors expressed in spinal cord dorsal horn synapses have been found to be specifically inhibited by inflammatory mediators [62, 63], which raises the possibility that ethanol and ethanol-induced inflammatory signals may initiate converging activities through glycine receptor signaling in the forebrain. Overall, this module emphasizes the primary role that activation of stress and inflammatory responses, as well as intracellular vesicle transport and neurotransmitter-mediated signaling, play in the establishment of ethanol-induced behaviors.
In the turquoise module, the CXCR4-mediated signaling pathway appeared as one of the most enriched. Concordantly, the CXCR4 gene was previously reported as one of the genes more significantly overtargeted by miRNAs in PFC of human alcoholics . This result underscores the relevance of this cytokine receptor pathway in the transduction of alcohol actions. CXCR4 has been reported to interact with the LPS-sensing complex that also includes CD14, Hsp70, Hsp90, GDF5, and TLR4 . As mentioned earlier, LPS-triggered signaling through this complex has been recently reported to increase ethanol consumption and to alter certain aspects of alcohol reward and aversion in mice . CXCR4 has also been implicated in opiate-induced hypernociception  and polymorphisms in SDF1 (aka CXCL12, the CXCR4 ligand) associated with several phenotypes, including alcohol consumption . In addition, CXCR4 has been implicated in glutamate exocytosis from astrocytes , differentiation of oligodendrocytes progenitors and remyelination , and appears to be required for proper distribution of GABAergic interneurons and the establishment of functional cortical circuitry in certain cortical regions and layers . Moreover, the CXCR4 ligand (CXCL12) induced presynaptic enhancement of GABA and glutamate release at serotonin dorsal raphe nucleus neurons .
The evidence shown here suggests that one of alcohol’s primary effects, in the early stages of development of dependence (represented here by the binge-drinking mouse model), is the transcriptional activation of genes involved in synaptic signaling, endocytic transport, and inflammatory response. A limitation of our study is the temporal resolution of the gene expression signal, which is a snapshot at the end of the drinking protocol. We speculate that our results describe the “intermediate state” of a dynamic process of mRNA regulation by miRNAs. Future studies need to explore gene expression profiles at additional time points during the drinking protocol and for an extended period of time and validate the regulatory relationship between miRNA/mRNA pairs. Based on our current model, we speculate that with time, the brain under ethanol insult mounts an adaptive miRNA response that counteracts the effect of the initial transcriptional activation. The well-adapted human alcoholic brain might represent this latter state. This hypothesis is consistent with processes of adaptive homeostasis generally triggered by addiction disorders [44, 70, 71]. The level of mechanistic conservation we observe between the mouse and human models is impressive, thereby warranting further investigation of the regulatory networks identified in our mouse model.
By simultaneously studying miRNA and mRNA expression profiles in the frontal cortex of ethanol-treated mice and comparing these with similar studies conducted in brains of human alcoholics, we provide evidence that alcohol drinking induces extensive activation of miRNA expression. We suggest that this is an adaptive response to modulate an initial broad activation of genes primarily involved in synaptic signaling, endocytic transport, and inflammatory response. Importantly, conserved subsets of miRNAs and mRNAs demonstrate changing expression levels in response to alcohol in both mice and humans. Our results provide the first evidence for positively correlated miRNA-mRNA interaction networks based on expression correlation analyses that suggest an adaptive miRNA response secondary to the activation of gene expression in brains under ethanol insult. This study underscores the value of integrative miRNA-mRNA analyses that could be critical for the future of addiction research. Further analysis should be conducted with independent expression profiling datasets to validate the use of these conserved miRNA and mRNA expression signatures as biomarkers for ethanol consumption. In vivo studies (i.e., manipulation of specific miRNA levels in mouse brain using direct stereotaxic injection in mutant mice) are currently in progress in our lab to elaborate on our findings and to discover specific miRNAs or relevant miRNA combinations that could alter drinking behaviors. Our results identify miRNAs as novel and promising therapeutic candidates for the treatment of alcoholism and other complex psychiatric disorders.
Mice and drinking protocol
Female hybrid F1 mice from reciprocal intercrosses of C57BL/6J x FVB/NJ F1 and FVB/NJ x C57BL/6J F1 (maternal strain x paternal strain) were used. These mice have been shown to drink the most ethanol at high concentrations and to achieve behaviorally significant blood ethanol concentrations . The B6 and FVB breeders were purchased from The Jackson Laboratory (Bar Harbor, ME) and mated at age 8 weeks in the Texas Genetic Animal Core of the INIA (Integrated Neuroscience Initiative on Alcohol) at the University of Texas at Austin. General methods implemented were as previously described . Alcohol drinking protocol followed a Two-Bottle Choice (2BC)-Drinking in the Dark (DID) paradigm. In this paradigm, mice were allowed to voluntarily drink 20% ethanol or water from two separate bottles (bottle positions were alternated daily) during a 3-hour time span starting three hours into the dark cycle. This modification of the DID test was published by Blednov and Harris . The drinking protocol was maintained for a total of 20 days and mice sacrificed by cervical dislocation 20 hours after the last drinking session. Brains were quickly removed and snap frozen in liquid nitrogen. Thirty-two mice were used in this study: 20 ethanol-treated and 12 matched controls. This study is in compliance with animal research guidelines and was approved by the Institutional Review Board of the University of Texas at Austin.
Total RNA extraction from brain tissue and microarray expression profiling
Frontal cortex (FCtx) tissue was dissected using a brain slicer (Zivic Instruments, Pittsburg, PA) to produce a 2-mm coronal section from the most rostral portion of the mouse brain devoid of olfactory bulbs (coordinates Bregma +1.56 to +3.56). The dorsal part of this coronal section, cut immediately above the forceps minor of the corpus callosum as the anatomical landmark, was used for RNA extraction. This section of the cortex is mostly composed of frontal associated cortex (FrA), cingulate cortex area 1 (Cg1), prelimbic cortex (PrL), and primary (M1) and secondary (M2) motor cortices, as depicted in the mouse brain atlas . Samples from both alcohol-treated and control groups were always included in each batch of extracted RNA. Total RNA was extracted using the mirVana® miRNA Isolation kit (Ambion, Austin, TX) according to the manufacturer’s instructions. Yield and quality of the total RNA preparation was determined using the Agilent 2100 Bioanalyzer (Agilent, Palo Alto, CA). For mRNA expression profiling, biotin-labeled cRNA was prepared using Illumina TotalPrep RNA Amplification kit (Ambion, Austin, TX) and then hybridized to Illumina MouseRef-8 v2.0 Expression BeadChips (Illumina, San Diego, CA). The quality of the Illumina bead summary data was assessed using the Bioconductor packages Lumi and arrayQualityMetrics. Data preprocessing included variance stabilization and quantile normalization using the Lumi package. Statistical analysis comparing ethanol-treated and control groups was performed using the Bioconductor package limma, which implements an empirical Bayes approach in R . For miRNA expression profiling, Exiqon miRCURY LNA microRNA Arrays 5th generation (Exiqon, Vedbaek, Denmark) were used for hybridization and scanning at the Moffitt Cancer Center Microarray Facility (Tampa, FL). Data analysis was performed using the limma package. Data preprocessing included minimum background correction and scale normalization between arrays. As each probe was spotted as four replicates on the arrays, within-array replication was assessed using the limma duplicate correlation function. False discovery rate (FDR) was assessed using the Benjamini-Hochberg method . Our list of differentially expressed genes was compared to the list of differentially expressed genes from human alcoholics reported previously , and also compared to the list of differentially expressed genes from a mouse model of alcohol dependence reported previously . The statistical significance of the matches was empirically evaluated by implementing Monte Carlo simulations in the R environment (R Core Team 2012. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. ISBN 3-900051-07-0, URL http://www.R-project.org/).
Validation by real-time PCR analysis
For miRNA validation, single-stranded cDNA was synthesized from total RNA using the TaqMan® MiRNA Reverse Transcription kit (Applied Biosystems, Foster City, CA). Following reverse transcription, quantitative real time PCR (qPCR) was performed in triplicate, using TaqMan® MiRNA Assays together with the TaqMan® Universal PCR Master Mix (Applied Biosystems), as per manufacturer’s instructions. TaqMan® miRNA assays used were: mmu-miR-7a-5p (ID: 000268), mmu-miR-15b-5p (ID: 000390), mmu-miR-101a-3p (ID: 002253), mmu-miR-101b-3p (ID: 002531), mmu-miR-140-5p (ID: 001187), mmu-miR-152-3p (ID: 000475), mmu-miR-195a-5p (ID: 000494), mmu-miR-541-5p (ID: 002562), mmu-let-7c-5p (ID: 000379), mmu-let-7f-5p (ID: 000382), and mmu-let-7g-5p (ID: 002282). Assays used for endogenous control were: snoRNA142 (ID: 001231) and snoRNA234 (ID: 001234), which were chosen out of five endogenous control assays tested in our samples. Reactions were carried out in a 7900HT Fast Real-Time PCR System and data collected using SDS software (Applied Biosystems). For mRNA validation, single-stranded cDNA was synthesized from total RNA using the TaqMan® High Capacity cDNA Reverse Transcription Kit (Applied Biosystems). Following reverse transcription, qPCR was performed in triplicate, using TaqMan® Gene Expression Assays together with the TaqMan® Universal PCR Master Mix (Applied Biosystems), as per manufacturer’s instructions. TaqMan® Gene expression assays used were: Syt11 (ID: Mm00444517_m1), Atp2b1(ID: Mm01245805_m1), and Fermt2 (ID: Mm00600590_m1). Assays Hprt (ID: Mm00446968_m1) and Gusb (ID: Mm00446953_m1) were used as endogenous controls after being selected as the most invariable control assays tested in our samples. GenEx software (MultiD Analyses AB, Gothenburg, Sweden) was used for analysis of real-time PCR data, including selection of endogenous controls using the GeNorm  and NormFinder  algorithms.
Weighted Gene Coexpression Network Analysis (WGCNA)
WGCNA was conducted as described previously  with modifications of a few parameters mentioned below. The general framework of WGCNA has been described in detail elsewhere . In short, we constructed a signed network by calculating the Pearson’s correlations for all pairs of genes and the signed similarity (Sij) matrix derived from Sij = (1 + cor(xi,xj))/2, where gene expression profiles xi and xj consist of the expression of genes i and j across multiple microarray samples. In the signed network, the direction of the changes in expression profiles can be inferred. The signed similarity (Sij) was then raised to the power β to represent the connection strength (aij): aij = Sijβ. This step aims to emphasize strong correlations and reduce the emphasis of weak correlations on an exponential scale . We chose a power of β = 12 so that the resulting networks exhibited approximate scale-free topology (Soft.R.sq = 0.88). All genes were hierarchically clustered based on a dissimilarity measure of topological overlap which measures interconnectedness for a pair of genes . The resulting gene dendrogram was used for module detection using the dynamic tree cut method and blockwiseModules function parameters: block size = 15000, minimum module size = 100, cutting height = 0.99, and deepSplit = 4. Gene modules corresponding to the branches cut off of the gene tree were labeled in unique colors. Unassigned genes were assigned to the grey module. Interaction networks were constructed for select modules. Only genes with high module membership (MM > 0.5) and high gene significance (GS > 0.5) were included in these networks. The trait-based gene significance measure is defined in WGCNA as the absolute value of the correlation between a specific node (gene) profile and the sample trait (ethanol consumption). The higher the value of GS, the more biologically significant the specific gene. On the other hand, module membership, also known as eigengene-based connectivity, correlates the gene expression profile of the specific gene with the module eigengene of a given module. Highly connected intramodular hub genes tend to have high module membership values to the respective module . The WGCNA analysis was first conducted on all mice (controls and alcohol-treated) and included all genes on the chips that were present (expressed over noise levels) in 75% or more of the samples. This first network is not unique to alcohol but represent an admixture of alcohol and control expression levels, which are later correlated with ethanol consumption to make inferences about the changes in expression being driven by ethanol consumption. A second analysis was conducted on each group independently in order to assess whether changes in network topology occurred between the coexpression network in control animals and the coexpression network in alcohol-treated animals.
Construction of miRNA-mRNA interaction database
A universe of miRNA-mRNA interactions for all reported miRNAs was constructed as a composite of all validated interactions downloaded from TarBase (http://diana.cslab.ece.ntua.gr/tarbase)  and predicted targets reported by miRecords (http://miRecords.biolead.org) . The predicted component integrates the predictions produced by 11 established miRNA target prediction programs (DIANA-microT, MicroInspector, miRanda, MirTarget2, miTarget, NBmiRTar, PicTar, PITA, RNA22, RNAhybrid, and TargetScan⁄TargertScanS). Predictions were filtered to only consider those targets predicted by at least 4 of 11 prediction algorithms. The TarBase validated interactions were then added to the filtered list and duplicated interactions eliminated (validated interactions were given preference) to avoid overcounting during overrepresentation analyses. miRecords’ predicted target information for every mouse miRNA accessible through the website was obtained using the RCurl package. MicroRNA family information was downloaded from miRBase (http://miRBase.org)  and included into the miRNA-mRNA interaction universe.
Integrative analysis of miRNA and mRNA microarray data
Statistical tests for comparison of multiple experimental proportions extracted from the integrated data set containing information about differentially expressed miRNAs, differentially expressed mRNAs, and putative/validated miRNA targets were implemented in the R environment. Specifically, hypergeometric tests for each differentially expressed miRNA were conducted to determine whether corresponding targets were overrepresented among inversely and directly correlated differentially expressed mRNAs from the same samples. Hypergeometric tests for each differentially expressed mRNA that is predicted to be targeted by more than one miRNA were conducted to determine whether any of such mRNAs could be subjected to “over-targeting” by miRNAs (an indication of redundant regulation or potential combinatorial action of miRNAs). When comparisons were made between differentially expressed sets of miRNAs from distinct species (i.e., mouse vs. human and mouse vs. rat) the species-specific miRNA IDs were converted to family IDs, which are consistent across species. MicroRNA families were downloaded from miRBase and ID conversion and matching code automated in R. When comparisons were implemented between differentially expressed sets of genes from distinct species (i.e., mouse vs. human), the human gene symbols were first converted into mouse symbols based on homology matching. The homology conversion table was downloaded from the Homologene database (http://ftp.ncbi.nih.gov/pub/HomoloGene/current/homologene.data) and ID conversion and matching code automated in R.
Functional enrichment analyses
Enrichment analysis of gene ontology annotations were conducted using lists of differentially expressed mRNAs and module-specific mRNAs using ToppGene functional analysis tool and the R packages GO.db and GOstats. To determine which cell types could have contributed to the detected gene expression changes and coexpression patterns in each module, cell-type enrichment analysis was conducted as described by Ponomarev and collaborators . Gene sets known to be preferentially expressed in mouse neurons, oligodendrocytes, astrocytes, and microglia were used. Statistical analysis was performed using the hypergeometric test and corrections for multiple comparisons using the Benjamini-Hochberg method.
Availability of supporting data
The data discussed in this publication have been deposited in NCBI's Gene Expression Omnibus  and are accessible through GEO Series accession number GSE50427 (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE50427).
We thank Dr. Jody Mayfield for comments on the manuscript and Drs. Sean Farris and Igor Ponomarev for useful discussions while conducting data analyses. This work was supported by NIAAA grants 5R01AA012404, 5P20AA017838, 5U01AA013520, and P01AA020683. YON was funded by postdoctoral fellowship under NIAAA Institutional Research Training Grant 5T32AA007471-24/25.
- Gorini G, Bell RL, Mayfield RD: Molecular Targets of Alcohol: Action Translational Research for Pharmacotherapy Development and Screening. Prog Mol Biol Transl Sci. 2011, 98: 293-347. doi:10.1016/B978-0-12-385506-0.00007-7View ArticlePubMedGoogle Scholar
- Farris SP, Miles MF: Ethanol Modulation of Gene Networks: Implications for Alcoholism. Neurobiol Dis. 2011, 1-7. doi:10.1016/j.nbd.2011.04.013Google Scholar
- Ron D, Messing RO: Signaling Pathways Mediating Alcohol Effects. Current Topics in Behavioural Neurosciences. 2011, doi:10.1007/7854_2011_161Google Scholar
- Nunez YO, Mayfield RD: Understanding Alcoholism Through microRNA Signatures in Brains of Human Alcoholics. Front Genet. 2012, 3: 43-doi:10.3389/fgene.2012.00043PubMed CentralView ArticlePubMedGoogle Scholar
- Mayfield RD, Harris RA, Schuckit MA: Genetic Factors Influencing Alcohol Dependence. Br J Pharmacol. 2008, 154 (2): 275-287. doi:10.1038/bjp.2008.88PubMed CentralView ArticlePubMedGoogle Scholar
- Pietrzykowski AZ: Coinciding Revolutions: How Discovery of Non-Coding DNA and RNA Can Change Our Understanding of Addiction. Front Genet. 2012, 3: 271-doi:10.3389/fgene.2012.00271PubMed CentralView ArticlePubMedGoogle Scholar
- Eipper-Mains JE, Eipper BA, Mains RE: Global Approaches to the Role of miRNAs in Drug-Induced Changes in Gene Expression. Front Genet. 2012, 3 (109): doi:10.3389/fgene.2012.00109Google Scholar
- Miranda RC: MicroRNAs and Fetal Brain Development: Implications for Ethanol Teratology During the Second Trimester Period of Neurogenesis. Front Genet. 2012, 3: 77-doi:10.3389/fgene.2012.00077PubMed CentralView ArticlePubMedGoogle Scholar
- Reilly M: Role of Non-Coding RNAs in the Neuroadaptation to Alcoholism and Fetal Alcohol Exposure. Front Genet. 2012, 3: 70-doi:10.3389/fgene.2012.00070PubMed CentralView ArticlePubMedGoogle Scholar
- Sartor GC, Laurent GS, Wahlestedt C: The Emerging Role of Non-Coding RNAs in Drug Addiction. Front Genet. 2012, 3: 106-doi:10.3389/fgene.2012.00106PubMed CentralView ArticlePubMedGoogle Scholar
- Lewohl JM, Nunez YO, Dodd PR, Tiwari GR, Harris RA, Mayfield RD: Up-Regulation of MicroRNAs in Brain of Human Alcoholics. Alcohol Clin Exp Res. 2011, 35 (11): 1928-1937. 10.1111/j.1530-0277.2011.01544.x. doi:10.1111/j.1530-0277.2011.01544.xPubMed CentralView ArticlePubMedGoogle Scholar
- Sprow GM, Thiele TE: The Neurobiology of Binge-Like Ethanol Drinking: Evidence From Rodent Models. Physiology & Behaviour. 2012, 106 (3): 325-331. 10.1016/j.physbeh.2011.12.026. doi:10.1016/j.physbeh.2011.12.026View ArticleGoogle Scholar
- Tapocik JD, Solomon M, Flanigan M, Meinhardt M, Barbier E, Schank JR, Schwandt M, Sommer WH, Heilig M: Coordinated Dysregulation of mRNAs and microRNAs in the Rat Medial Prefrontal Cortex Following a History of Alcohol Dependence. 2012, 1-11. doi:10.1038/tpj.2012.17\Google Scholar
- Liu J, Lewohl JM, Harris RA, Iyer VR, Dodd PR, Randall PK, Mayfield RD: Patterns of Gene Expression in the Frontal Cortex Discriminate Alcoholic From Nonalcoholic Individuals. Neuropsychopharmacology. 2006, 31 (7): 1574-1582. 10.1038/sj.npp.1300947. doi:10.1038/sj.npp.1300947View ArticlePubMedGoogle Scholar
- Ponomarev I, Wang S, Zhang L, Harris RA, Mayfield RD: Gene Coexpression Networks in Human Brain Identify Epigenetic Modifications in Alcohol Dependence. J Neurosci. 2012, 32 (5): 1884-1897. 10.1523/JNEUROSCI.3136-11.2012. doi:10.1523/JNEUROSCI.3136-11.2012PubMed CentralView ArticlePubMedGoogle Scholar
- Ravasz E, Somera AL, Mongru DA, Oltvai ZN, Barabási AL: Hierarchical Organization of Modularity in Metabolic Networks. Science. 2002, 297 (5586): 1551-1555. 10.1126/science.1073374. doi:10.1126/science.1073374View ArticlePubMedGoogle Scholar
- Zhang B, Horvath S: A General Framework for Weighted Gene Co-Expression Network Analysis. Stat Appl Genet Mol Biol. 2005, 8 (17): 8-doi:10.2202/1544-6115.1128Google Scholar
- Oldham MC, Konopka G, Iwamoto K, Langfelder P, Kato T, Horvath T, Geschwind DH: Functional Organization of the Transcriptome in Human Brain. Nat Neurosci. 2008, 11 (11): 1271-1282. 10.1038/nn.2207. doi:10.1038/nn.2207PubMed CentralView ArticlePubMedGoogle Scholar
- Saba LM, Bennett B, Hoffman PL, Barcomb K, Ishii T, Kechris K, Tabakoff B: A systems genetic analysis of alcohol drinking by mice, rats and men: Influence of brain GABAergic transmission. Neuropharmacology. 2011, 60 (7–8): 1269-1280. doi:10.1016/j.neuropharm.2010.12.019PubMed CentralView ArticlePubMedGoogle Scholar
- Allen NJ, Barres BA: Neuroscience: Glia - More Than Just Brain Glue. Nature. 2009, doi:10.1038/457675aGoogle Scholar
- Cullen PJ: Endosomal Sorting and Signalling: an Emerging Role for Sorting Nexins. Nat Rev Mol Cell Biol. 2008, 9 (7): 574-582. 10.1038/nrm2427.View ArticlePubMedGoogle Scholar
- Pant S, Sharma M, Patel K, Caplan S, Carr CM, Grant BD: AMPH-1/Amphiphysin/Bin1 Functions with RME-1/Ehd1 in Endocytic Recycling. Nat Cell Biol. 2009, 11 (12): 1399-1410. 10.1038/ncb1986. doi:10.1038/ncb1986PubMed CentralView ArticlePubMedGoogle Scholar
- Stockinger W, Sailler B, Strasser V, Recheis B, Fasching D, Kahr L, Schneider WJ, Nimpf J: The PX-Domain Protein SNX17 Interacts with Members of the LDL Receptor Family and Modulates Endocytosis of the LDL Receptor. EMBO J. 2002, 21 (16): 4259-4267. 10.1093/emboj/cdf435.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang T, Liu NS, Seet L, Hong W: The Emerging Role of VHS Domain-Containing Tom1, Tom1L1 and Tom1L2 in Membrane Trafficking. Traffic (Copenhagen, Denmark). 2010, 11 (9): 1119-1128. 10.1111/j.1600-0854.2010.01098.x. doi:10.1111/j.1600-0854.2010.01098.xView ArticleGoogle Scholar
- Blednov YA, Ozburn AR, Walker D, Ahmed S, Belknap JK, Harris RA: Hybrid Mice as Genetic Models of High Alcohol Consumption. Behaviour Genetics. 2010, 40 (1): 93-110. 10.1007/s10519-009-9298-4. doi:10.1007/s10519-009-9298-4View ArticleGoogle Scholar
- Blednov YA, Metten P, Finn DA, Rhodes JS, Bergeson SE, Harris RA, Crabbe JC: Hybrid C57BL/6J X FVB/NJ Mice Drink More Alcohol Than Do C57BL/6J Mice. Alcohol Clin Exp Res. 2005, 29 (11): 1949-1958. 10.1097/01.alc.0000187605.91468.17.PubMed CentralView ArticlePubMedGoogle Scholar
- Chandrasekar V, Dreyer J: Regulation of MiR-124, Let-7d, and MiR-181a in the Accumbens Affects the Expression, Extinction, and Reinstatement of Cocaine-Induced Conditioned Place Preference. Neuropsychopharmacology. 2011, 36 (6): 1149-1164. 10.1038/npp.2010.250. doi:10.1038/npp.2010.250PubMed CentralView ArticlePubMedGoogle Scholar
- Rizo J, Chen X, Arac D: Unraveling the Mechanisms of Synaptotagmin and SNARE Function in Neurotransmitter Release. Trends Cell Biol. 2006, 16 (7): 339-350. 10.1016/j.tcb.2006.04.006. doi:10.1016/j.tcb.2006.04.006View ArticlePubMedGoogle Scholar
- Jahn R, Fasshauer D: Molecular Machines Governing Exocytosis of Synaptic Vesicles. Nature. 2012, 490 (7419): 201-207. 10.1038/nature11320. doi:10.1038/nature11320PubMed CentralView ArticlePubMedGoogle Scholar
- Steinberg F, Heesom KJ, Bass MD, Cullen PJ: SNX17 Protects Integrins From Degradation by Sorting Between Lysosomal and Recycling Pathways. J Cell Biol. 2012, 197 (2): 219-230. 10.1083/jcb.201111121. doi:10.1083/jcb.201111121PubMed CentralView ArticlePubMedGoogle Scholar
- Etheridge N, Mayfield RD, Harris RA, Dodd PR: Identifying Changes in the Synaptic Proteome of Cirrhotic Alcoholic Superior Frontal Gyrus. Curr Neuropharmacol. 2011, 9 (1): 122-128. 10.2174/157015911795017164. doi:10.2174/157015911795017164PubMed CentralView ArticlePubMedGoogle Scholar
- Etheridge N, Lewohl JM, Mayfield RD, Harris RA, Dodd PR: Synaptic Proteome Changes in the Superior Frontal Gyrus and Occipital Cortex of the Alcoholic Brain. Proteomics - Clinical Applications. 2009, 3 (6): 730-742. 10.1002/prca.200800202. doi:10.1002/prca.200800202PubMed CentralView ArticlePubMedGoogle Scholar
- Gorini G, Ponomareva O, Shores KS, Person MD, Harris RA, Mayfield RD: Dynamin-1 Co-Associates with Native Mouse Brain BKCa Channels: Proteomics Analysis of Synaptic Protein Complexes. FEBS Lett. 2010, 584 (5): 845-851. 10.1016/j.febslet.2009.12.061. doi:10.1016/j.febslet.2009.12.061PubMed CentralView ArticlePubMedGoogle Scholar
- Malki K, Uher R, Paya-Cano J, Binder E, Rietschel M, Zobel A, Mors O: Convergent Animal and Human Evidence Suggests a Role of PPM1A Gene in Response to Antidepressants. Biol Psychiatry. 2011, 69 (4): 360-365. 10.1016/j.biopsych.2010.08.011.View ArticlePubMedGoogle Scholar
- Sun W, Yu Y, Dotti G, Shen T, Tan X, Savoldo B, Pass AK: PPM1A and PPM1B Act as IKKbeta Phosphatases to Terminate TNFalpha-Induced IKKbeta-NF-kappaB Activation. Cell Signal. 2009, 21 (1): 95-102. 10.1016/j.cellsig.2008.09.012. doi:10.1016/j.cellsig.2008.09.012PubMed CentralView ArticlePubMedGoogle Scholar
- Aroor AR, Shukla SD: MAP Kinase Signaling in Diverse Effects of Ethanol. Life Sci. 2004, 74 (19): 2339-2364. 10.1016/j.lfs.2003.11.001.View ArticlePubMedGoogle Scholar
- Moonat S, Starkman BG, Sakharkar A, Pandey SC: Neuroscience of Alcoholism: Molecular and Cellular Mechanisms. Cell Mol Life Sci. 2010, 67 (1): 73-88. 10.1007/s00018-009-0135-y. doi:10.1007/s00018-009-0135-yPubMed CentralView ArticlePubMedGoogle Scholar
- Aroor AR, James TT, Jackson DE, Shukla SD: Differential Changes in MAP Kinases, Histone Modifications, and Liver Injury in Rats Acutely Treated with Ethanol. Alcohol Clin Exp Res. 2010, 34: 1543-1551. 10.1111/j.1530-0277.2010.01239.x. doi:10.1111/j.1530-0277.2010.01239.xPubMed CentralView ArticlePubMedGoogle Scholar
- Chandler LJ, Sutton G: Acute Ethanol Inhibits Extracellular Signal-Regulated Kinase, Protein Kinase B, and Adenosine 3':5’-Cyclic Monophosphate Response Element Binding Protein Activity in an Age- and Brain Region-Specific Manner. Alcohol Clin Exp Res. 2005, 29 (4): 672-682. 10.1097/01.ALC.0000158935.53360.5F.View ArticlePubMedGoogle Scholar
- Hoffman PL, Miles M, Edenberg HJ, Sommer W, Tabakoff B, Wehner JM, Lewohl J: Gene Expression in Brain: a Window on Ethanol Dependence, Neuroadaptation, and Preference. 2003, 27: 155-168. doi:10.1097/01.ALC.0000060101.89334.11Google Scholar
- Mulligan MK, Ponomarev I, Hitzemann RJ, Belknap JK, Tabakoff B, Harris RA, Crabbe JC: Toward Understanding the Genetics of Alcohol Drinking Through Transcriptome Meta-Analysis. Proc Natl Acad Sci USA. 2006, 103 (16): 6368-6373. 10.1073/pnas.0510188103.PubMed CentralView ArticlePubMedGoogle Scholar
- Sanna PP, Simpson C, Lutjens R, Koob G: ERK Regulation in Chronic Ethanol Exposure and Withdrawal. Brain Res. 2002, 948 (1–2): 186-191.View ArticlePubMedGoogle Scholar
- Ökvist A, Johansson S, Kuzmin A, Bazov I, Merino-Martinez R, Igor Ponomarev R, Mayfield D: Neuroadaptations in Human Chronic Alcoholics: Dysregulation of the NF-kappaB System. PLoS ONE. 2007, 2 (9): e930-10.1371/journal.pone.0000930. doi:10.1371/journal.pone.0000930PubMed CentralView ArticlePubMedGoogle Scholar
- Yakovleva T, Bazov I, Watanabe H, Hauser KF, Bakalkin G: Transcriptional Control of Maladaptive and Protective Responses in Alcoholics: a Role of the NF-κB System. Brain Behaviour and Immunity. 2011, 25 (1): S29-S38. doi:10.1016/j.bbi.2010.12.019View ArticleGoogle Scholar
- Zou J, Crew F: Induction of Innate Immune Gene Expression Cascades in Brain Slice Cultures by Ethanol: Key Role of NF-κB and Proinflammatory Cytokines. Alcohol Clin Exp Res. 2010, 34 (5): 777-789. 10.1111/j.1530-0277.2010.01150.x. doi:10.1111/j.1530-0277.2010.01150.xView ArticlePubMedGoogle Scholar
- Blednov YA, Benavidez JM, Geil C, Perra S, Morikawa H, Harris RA: Activation of Inflammatory Signaling by Lipopolysaccharide Produces a Prolonged Increase of Voluntary Alcohol Intake in Mice. Brain Behaviour and Immunity. 2011, 25 (Suppl): S92-S105. doi:10.1016/j.bbi.2011.01.008View ArticleGoogle Scholar
- Edenberg HJ, Xuei X, Wetherill LF, Bierut L, Bucholz K, Dick DM, Hesselbrock V: Association of NFKB1, Which Encodes a Subunit of the Transcription Factor NF-kappaB, with Alcohol Dependence. Hum Mol Genet. 2008, 17 (7): 963-970. doi:10.1093/hmg/ddm368View ArticlePubMedGoogle Scholar
- Liu J, Yang AR, Kelly T, Puche A, Esoga C, June HL, Elnabawi A: Binge Alcohol Drinking Is Associated with GABAA Alpha2-Regulated Toll-Like Receptor 4 (TLR4) Expression in the Central Amygdala. Proc Natl Acad Sci USA. 2011, 108 (11): 4465-4470. 10.1073/pnas.1019020108.PubMed CentralView ArticlePubMedGoogle Scholar
- Marcos M, Pastor I, Gonzalez-Sarmiento R, Laso FJ: A Functional Polymorphism of the NFKB1 Gene Increases the Risk for Alcoholic Liver Cirrhosis in Patients with Alcohol Dependence. Alcohol Clin Exp Res. 2009, 33 (11): 1857-1862. 10.1111/j.1530-0277.2009.01023.x. doi:10.1111/j.1530-0277.2009.01023.xView ArticlePubMedGoogle Scholar
- Hutt D, Balch WE: The Proteome in Balance. Science. 2010, 329 (5993): 766-767. 10.1126/science.1194160. doi:10.1126/science.1194160PubMed CentralView ArticlePubMedGoogle Scholar
- Okiyoneda T, Barriere H, Bagdany M, Rabeh WM, Du K, Hohfeld J, Young JC, Lukacs GL: Peripheral Protein Quality Control Removes Unfolded CFTR From the Plasma Membrane. Science. 2010, 329 (5993): 805-810. 10.1126/science.1191542. doi:10.1126/science.1191542View ArticlePubMedGoogle Scholar
- Gupta A, Cooper ZA, Tulapurkar ME, Potla R, Maity T, Hasday JD, Singh IS: Toll-Like Receptor Agonists and Febrile Range Hyperthermia Synergize to Induce Heat Shock Protein 70 Expression and Extracellular Release. J Biol Chem. 2013, 288 (4): 2756-2766. 10.1074/jbc.M112.427336. doi:10.1074/jbc.M112.427336PubMed CentralView ArticlePubMedGoogle Scholar
- Lee K-H, Jeong J, Yoo C-G: Positive Feedback Regulation of Heat Shock Protein 70 (Hsp70) Is Mediated Through Toll-Like Receptor 4-PI3K/Akt-Glycogen Synthase Kinase-3β Pathway. Exp Cell Res. 2013, 319 (1): 88-95. 10.1016/j.yexcr.2012.09.018. doi:10.1016/j.yexcr.2012.09.018View ArticlePubMedGoogle Scholar
- Miles MF, Diaz JE, DeGuzman VS: Mechanisms of Neuronal Adaptation to Ethanol. Ethanol Induces Hsc70 Gene Transcription in NG108-15 Neuroblastoma X Glioma Cells. J Biol Chem. 1991, 266 (4): 2409-2414.PubMedGoogle Scholar
- Pignataro L, Varodayan FP, Tannenholz LE, Harrison NL: The Regulation of Neuronal Gene Expression by Alcohol. Pharmacol Ther. 2009, 124 (3): 324-335. 10.1016/j.pharmthera.2009.09.002.PubMed CentralView ArticlePubMedGoogle Scholar
- Wilke N, Sganga MW, Gayer GG, Hsieh KP, Miles MF: Characterization of Promoter Elements Mediating Ethanol Regulation of Hsc70 Gene Transcription. J Pharmacol Exp Ther. 2000, 292 (1): 173-180.PubMedGoogle Scholar
- Triantafilou M, Triantafilou K: Lipopolysaccharide Recognition: CD14, TLRs and the LPS-Activation Cluster. Trends Immunol. 2002, 23 (6): 301-304. 10.1016/S1471-4906(02)02233-0.View ArticlePubMedGoogle Scholar
- Blednov YA, Benavidez JM, Homanics GE, Harris RA: Behavioural Characterization of Knockin Mice with Mutations M287L and Q266I in the Glycine Receptor Α1 Subunit. J Pharmacol Exp Ther. 2012, 340 (2): 317-329. 10.1124/jpet.111.185124. doi:10.1124/jpet.111.185124PubMed CentralView ArticlePubMedGoogle Scholar
- Lynch JW: Molecular Structure and Function of the Glycine Receptor Chloride Channel. Physiol Rev. 2004, 84 (4): 1051-1095. 10.1152/physrev.00042.2003. doi:10.1152/physrev.00042.2003View ArticlePubMedGoogle Scholar
- McCracken LM, Blednov YA, Trudell JR, Benavidez JM, Betz H, Harris RA: Mutation of a Zinc-Binding Residue in the Glycine Receptor 1 Subunit Changes Ethanol Sensitivity in Vitro and Alcohol Consumption in Vivo. J Pharmacol Exp Ther. 2013, 344 (2): 489-500. 10.1124/jpet.112.197707.PubMed CentralView ArticlePubMedGoogle Scholar
- Mihic SJ, Ye Q, Wick MJ, Koltchine VV, Krasowski MD, Finn SE, Mascia MP: Sites of Alcohol and Volatile Anaesthetic Action on GABA(a) and Glycine Receptors. Nature. 1997, 389 (6649): 385-389. 10.1038/38738. doi:10.1038/38738View ArticlePubMedGoogle Scholar
- Harvey RJ, Depner UB, Wässle H, Ahmadi S, Heindl C, Reinold H, Smart TG: GlyR Alpha3: an Essential Target for Spinal PGE2-Mediated Inflammatory Pain Sensitization. Science. 2004, 304 (5672): 884-887. 10.1126/science.1094925. doi:10.1126/science.1094925View ArticlePubMedGoogle Scholar
- Lynch JW, Callister RJ: Glycine Receptors: a New Therapeutic Target in Pain Pathways. Curr Opin Investig Drugs. 2006, 7 (1): 48-53.PubMedGoogle Scholar
- White F, Wilson N: Opiate-Induced Hypernociception and Chemokine Receptors. Neuropharmacology. 2010, 58 (1): 35-37. 10.1016/j.neuropharm.2009.07.012. doi:10.1016/j.neuropharm.2009.07.012PubMed CentralView ArticlePubMedGoogle Scholar
- Xiao Q, Ye S, Oberhollenzer F, Mayr A, Jahangiri M, Willeit J, Kiechl S, Xu Q: SDF1 Gene Variation Is Associated with Circulating SDF1alpha Level and Endothelial Progenitor Cell Number: the Bruneck Study. PLoS ONE. Edited by: Qingzhong X, Shu Y, Friedrich O, Agnes M, Marjan J, Johann W, Stefan K, Qingbo X. 2008, (12) (December 30): e4061. doi:10.1371/journal.pone.0004061.t003Google Scholar
- Calì C, Bezzi P: CXCR4-Mediated Glutamate Exocytosis From Astrocytes. J Neuroimmunol. 2010, 224 (1–2): 13-21. doi:10.1016/j.jneuroim.2010.05.004View ArticlePubMedGoogle Scholar
- Patel JR, McCandless EE, Dorsey D, Klein RS: CXCR4 Promotes Differentiation of Oligodendrocyte Progenitors and Remyelination. Proc Natl Acad Sci USA. 2010, 107 (24): 11062-11067. 10.1073/pnas.1006301107. doi:10.1073/pnas.1006301107PubMed CentralView ArticlePubMedGoogle Scholar
- Tanaka DH, Mikami S, Nagasawa T, Miyazaki JI, Nakajima K, Murakami F: CXCR4 Is Required for Proper Regional and Laminar Distribution of Cortical Somatostatin-, Calretinin-, and Neuropeptide Y-Expressing GABAergic Interneurons. Cerebral Cortex (New York, N.Y.: 1991). 2010, 20 (12): 2810-2817. 10.1093/cercor/bhq027. doi:10.1093/cercor/bhq027View ArticleGoogle Scholar
- Heinisch S, Kirby LG: SDF-1alpha/CXCL12 Enhances GABA and Glutamate Synaptic Activity at Serotonin Neurons in the Rat Dorsal Raphe Nucleus. Neuropharmacology. 2010, 58 (2): 501-514. 10.1016/j.neuropharm.2009.08.022. doi:10.1016/j.neuropharm.2009.08.022PubMed CentralView ArticlePubMedGoogle Scholar
- Kalivas PW: The Glutamate Homeostasis Hypothesis of Addiction. Nat Rev Neurosci. 2009, 10 (8): 561-572. 10.1038/nrn2515. doi:10.1038/nrn2515View ArticlePubMedGoogle Scholar
- Koob GF: Drug Addiction: the Yin and Yang of Hedonic Homeostasis. Neuron. 1996, 16: 893-896. 10.1016/S0896-6273(00)80109-9.View ArticlePubMedGoogle Scholar
- Blednov YA, Harris RA: Metabotropic Glutamate Receptor 5 (mGluR5) Regulation of Ethanol Sedation, Dependence and Consumption: Relationship to Acamprosate Actions. Int J Neuropsychopharmacol. 2008, 11 (6): 775-793. doi:10.1017/S1461145708008584PubMed CentralView ArticlePubMedGoogle Scholar
- Franklin K, Paxinos G: The Mouse Brain in Streotaxic Coordinates. 2007, San Diego, CA: Academic Press, 1-351.Google Scholar
- Smyth GK: “Limma: Linear Models for Microarray Data”. In … And Computational Biology Solutions Using R and …, 397–420. Statistics for Biology and Health. 2005, New York: Springer-Verlag, doi:10.1007/0-387-29362-0_23Google Scholar
- Benjamini Y, Hochberg Y: Controlling the False Discovery Rate: a Practical and Powerful Approach to Multiple Testing. J R Stat Soc. 1995, 57( (1): 289-300. Series B (Methodological)Google Scholar
- Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, Speleman F: Accurate Normalization of Real-Time Quantitative RT-PCR Data by Geometric Averaging of Multiple Internal Control Genes. Genome Biol. 2002, 3 (7): RESEARCH0034-PubMed CentralView ArticlePubMedGoogle Scholar
- Andersen CL, Jensen JL, Ørntoft TF: Normalization of Real-Time Quantitative Reverse Transcription-PCR Data: a Model-Based Variance Estimation Approach to Identify Genes Suited for Normalization, Applied to Bladder and Colon Cancer Data Sets. Cancer Res. 2004, 64 (15): 5245-5250. 10.1158/0008-5472.CAN-04-0496. doi:10.1158/0008-5472.CAN-04-0496View ArticlePubMedGoogle Scholar
- Zhao W, Langfelder P, Fuller T, Dong J, Li A, Hovarth S: Weighted Gene Coexpression Network Analysis: State of the Art. J Biopharm Stat. 2010, 20 (2): 281-300. 10.1080/10543400903572753. doi:10.1080/10543400903572753View ArticlePubMedGoogle Scholar
- Vergoulis T, Vlachos IS, Alexiou P, Georgakilas G, Maragkakis M, Reczko M, Gerangelos S, Koziris N, Dalamagas T, Hatzigeorgiou AG: TarBase 6.0: Capturing the Exponential Growth of miRNA Targets with Experimental Support. Nucleic Acids Res. 2012, 40 (Database issue): 9-doi:10.1093/nar/gkr1161Google Scholar
- Xiao F, Zuo Z, Cai G, Kang S, Gao X, Li T: miRecords: an Integrated Resource for microRNA-Target Interactions. Nucleic Acids Res. 2009, 37 (Database): D105-D110. 10.1093/nar/gkn851. doi:10.1093/nar/gkn851PubMed CentralView ArticlePubMedGoogle Scholar
- Griffiths-Jones S, Grocock RJ, van Dongen S, Bateman A, Enrigh AJ: miRBase: microRNA Sequences, Targets and Gene Nomenclature. Nucleic Acids Res. 2006, 34 (Database issue): D140-D144.PubMed CentralView ArticlePubMedGoogle Scholar
- Edgar R, Domrachev M, Lash AE: Gene Expression Omnibus: NCBI Gene Expression and Hybridization Array Data Repository. Nucleic Acids Res. 2002, 30: 207-210. 10.1093/nar/30.1.207.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.