- Research article
- Open Access
Positively correlated miRNA-mRNA regulatory networks in mouse frontal cortex during early stages of alcohol dependence
BMC Genomics volume 14, Article number: 725 (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
In order to evaluate the impact of alcohol consumption on miRNA and gene expression in mouse brain, we conducted miRNA and mRNA profiling studies on 32 samples (20 ethanol-treated mice and 12 matched controls) as described in Materials and Methods. Ethanol-treated mice consumed an average of 4.92 g/kg/3 h (range: 4.08-6.31 g/kg/3 h; stdev: 0.52). The miRNA expression analysis indicated that, as in prefrontal cortex (PFC) of human alcoholics, miRNAs appear predominantly upregulated in FCtx of alcohol-drinking mice, with 52 miRNA families upregulated in mouse brain (Table 1). Importantly, we found that a highly significant number of miRNA families that are upregulated in PFC of human alcoholics  are also upregulated in the FCtx of ethanol-treated mice (Figure 1). Fourteen out of 32 alcohol-induced human miRNA families changed expression in the ethanol-treated mice (P < 1 × 10-5 as determined after 10,000 Monte Carlo simulations). The fact that upregulation is predominant at the significance level chosen (FDR < 10%) among differentially expressed miRNAs in the mouse frontal cortex is consistent with previous results in human alcoholic PFC . The miRNA families that change expression in both mouse and human were: let-7, miR-7, miR-15, miR-101, miR-140, miR-152 (all validated by qPCR, P < 0.05), as well as miR-17, miR-34, miR-135, miR-144, miR-146, miR-301, miR-339, miR-368 (qPCR not performed). Other differentially expressed miRNAs specific to the mouse model were also validated by qPCR, including miR-195 (member of the miR-15 family) and miR-541 family members. Moreover, we found a significant match among differentially expressed mouse miRNA families and those reported in PFC of ethanol-treated rats . Seventeen out of 33 rat miRNA families (representing 41 differentially expressed rat miRNAs) were matched to our mouse data, P < 1 × 10-5. The miRNA families that change expression in both mice and rats were: mir-7, mir-9, mir-10, mir-15, mir-17, mir-26, mir-29, mir-30, mir-101, mir-130, mir-181, mir-204, mir-339, mir-340, mir-368, mir-434, mir-467. Overall, these results underscore the relevance of gene regulation by miRNAs in response to alcohol consumption and suggest conservation of alcohol-responsive miRNA regulatory pathways from rodent to human.
Ethanol-responsive genes in brains of mice significantly match respective homologs differentially expressed in brains of human alcoholics
Differential expression analysis of the Illumina microarray data identified 709 genes with expression changes in response to alcohol in the FCtx of ethanol-treated mice (Additional file 1: Table S1, FDR < 1%). A representative group of altered genes, including Syt11, Tom1, Atp2b1, and Fermt2, were validated by qPCR (P < 0.05) using TaqMan® assays as described in Materials and Methods. A significant number of the differentially expressed genes (Additional file 1: Figure S1A) also changed in PFC of human alcoholics, as reported by the Mayfield group [11, 14]. Such highly statistically significant matches (P < 1 × 10-5 as empirically estimated by Monte Carlo simulations) underscore the relevance of the commonly affected genes. However, it should be noted that 79% of the common gene changes are occurring in opposite directions when the mouse and human models are compared (Figure 2A). When our mouse gene expression dataset was compared to an additional profiled cohort of human alcoholics , we consistently found a highly significant match among differentially expressed genes between the mouse and human models (84 common differentially expressed genes, P < 1 × 10-7, (Additional file 1: Figure S1B). Moreover, we similarly observed that a relatively large percentage (40%) changed in opposite directions (Figure 2B). Furthermore, comparison to the human datasets indicated that a majority of the common differentially expressed genes (65% on average) are upregulated in the mouse brain, while only 34% on average are upregulated in the brain of human alcoholics. These results suggest that brain genes upregulated in early stages of development of alcohol dependence (using mouse model) might undergo downregulation in late stages of the disease (using human alcoholic model), possibly due to homeostatic mechanisms driven by miRNA regulation.
Integrated miRNA-mRNA networks reveal positively correlated action of brain miRNAs and persistent upregulation of overtargeted mRNAs in ethanol drinking mice
In order to better understand miRNA-mRNA regulatory relationships, we constructed respective interaction networks among negatively and positively correlated interaction groups, using the miRNA-mRNA interaction universe generated from miRecords predicted (consensus of at least 4 tools) and validated interactions (described in Materials and Methods). We found that the negatively correlated network of interactions among upregulated miRNAs and downregulated mRNAs (average number of neighbors: 2.95, Figure 3A) was no different from randomized networks generated as a control tool (average number of neighbors: 3.03, shown in Additional file 1: Figure S2A). On the other hand, the corresponding network containing the miRNA-mRNA interactions that were positively correlated (upregulated miRNA-upregulated mRNA network) provided unexpected evidence that this network was twice as interconnected (average number of neighbors: 4.97, Figure 3B) as expected by chance (average number of neighbors: 2.52, shown in Additional file 1: Figure S2B).
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
To determine the relevant alcohol-responsive pathways regulated by miRNAs, we investigated the organization of the transcriptome in the frontal cortex of the mouse brain by applying weighted gene coexpression network analysis to the 32 samples as described in Materials and Methods. The weighted network identified groups (modules) of genes with similar expression patterns and high topological overlap. Topological overlap is a pairwise measurement that describes the similarity of two genes’ coexpression relationships with all other genes in the network [16, 17] and has been shown to globally identify, for example, coexpression patterns that correlate with validated protein interactions in the human brain . Our WGCNA identified 10 mRNA modules (Figure 4) with highly significant coexpression patterns, which we graphically corroborated by plotting the expression levels for the top 10 genes with the strongest membership on alcohol-related modules (Additional file 1: Figure S3). All modules were further correlated using the module eigengenes, with the phenotypic traits of interest, namely average ethanol consumption per voluntary drinking session and expression pattern of differentially expressed miRNAs. Table 2 summarizes the WGCNA module detection results. We included correlation values to the ethanol consumption trait and the associated P value of these correlations, as well as results from a variety of enrichment analyses. Additional details for enrichment analyses are provided in Additional file 1: Tables S4-S7). We want to point out that the number of modules produced by our WGCNA strategy seems rather low and the number of unassigned genes (grey module) rather large. One of the reasons for this behavior is that we constrained the minimum number of genes per module to a relatively high number (100 genes) in an effort to generate larger modules that improve detection of functional enrichment after gene enrichment analysis. Moreover, we purposely chose a low cutting height of 0.99 (default in WGCNA R package is 0.995) while implementing the dynamic tree cut method of the blockwiseModules function, to increase the stringency of the module detection step. Cutting height sets a threshold for dissimilarity between expression patterns to be considered for module detection: the lower the threshold, the more similar the patterns need to be for assignment into modules. Yet another potentially influential factor is the relatively small difference in gene expression between control and ethanol-treated animals, which may make it more difficult to detect significantly different patterns at stringent levels of statistical significance. The relatively small changes in gene and miRNA expression observed in rodent and human models of alcohol consumption are a common phenomenon described by researchers in the field [11, 13, 19].
Five modules (yellow, red, turquoise, pink, and brown) presented significant average correlations with the ethanol-related traits (Table 2, P < 0.01; Figure 5A and 5B; P < 0.05), thus warranting further investigation. The alcohol-responsive modules also presented highly significant correlations between module membership (MM) and gene significance (GS) to the ethanol-consumption trait (Figure 5C), which strengthen the inference of biological functions of the alcohol-significant genes by analyzing ontology category enrichment in the respective modules. Moreover, the five modules significantly correlated with ethanol drinking were significantly enriched with differentially expressed genes (Table 2). We also detected significant correlations (FDR < 0.1) between individual gene expression patterns and individual differentially expressed miRNA expression patterns (Figure 4, Additional file 1: Table S8), which provide indirect experimental validation for a group of predicted miRNA-mRNA interactions. By correlating the module eigengenes with the differentially expressed miRNA expression patterns (Figure 6) and considering the individual significant correlations shown in Figure 4, three modules (red and brown –positively correlated with ethanol drinking, and turquoise –negatively correlated with ethanol drinking), were found to be preferentially targeted by the differentially expressed miRNAs. To assess the validity of this approach, we corroborated the positively correlated expression levels of interacting miRNAs and mRNAs (i.e., mmu-let-7 vs. Syt11 and mmu-let-7 vs. Tom1) by qPCR. Additional confirmation for WGCNA-suggested interactions was obtained from TarBase, the database for validated interactions, where we found validated interactions between miRNAs and mRNAs that were positively correlated in our study but negatively correlated in previous publications: mmu-miR-34a-5p --| ACTB, mmu-miR-200b-3p --| ZEB2, mmu-miR-30a-5p --| TNRC6A, mmu-miR-152-3p --| CAMK2A, mmu-miR-200c-3p --| FLT1, mmu-miR-20a-5p --| ZBTB7A. Also, we found evidence in TarBase for validated interactions between miRNAs and mRNAs that were negatively correlated in our study as well as in previous publications: mmu-miR-29b-5p --| COL4A2 and mmu-miR-30a-5p5p --| AK1. Apparent spurious correlation was observed among some differentially expressed miRNAs and the black module. However, this module did not show significant correlation with the alcohol consumption phenotype. Hence the correlations detected among the black module and differentially expressed miRNAs are understood as pleiotropic interactions, akin to off-target effects (in this case referring to the miRNA targeting of non-alcohol-relevant transcripts).
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)
To gain insight into the alcohol-responsive gene modules preferentially targeted by alcohol-induced miRNAs, we constructed gene interaction networks based on WGCNA correlation measurements for the red, brown, and turquoise modules. In the red module containing alcohol-upregulated genes, synaptotagmin-11 (Syt11), a member of the synaptotagmin gene family, which is involved in vesicle trafficking and synaptic transmission, is at the centre of the interaction network and the hub gene with higher connectivity (Figure 7A). 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 interleukin IL-1 receptor, among others [21–24]. Regarding miRNA targeting of these hub genes (Figure 7B), miR-34c-5p and let-7g-5p are main regulatory candidates based on significant expression pattern correlations (FDR < 0.10) with Syt11, Snx-7, Tom1 (both let-7g-5p and miR-34c-5p), and Bim1 (miR-34c-5p). This suggests that miR-34c-5p and let-7g-5p play an important role in the attempted regulation of this module. The central role played by miR-34c-5p and let-7g-5p while targeting additional genes in the red module is also evident in the expression correlation network shown in Additional file 1: Figure S4 and S5 (note the higher connectivity of these particular miRNAs).
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).
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-7
Farris SP, Miles MF: Ethanol Modulation of Gene Networks: Implications for Alcoholism. Neurobiol Dis. 2011, 1-7. doi:10.1016/j.nbd.2011.04.013
Ron D, Messing RO: Signaling Pathways Mediating Alcohol Effects. Current Topics in Behavioural Neurosciences. 2011, doi:10.1007/7854_2011_161
Nunez YO, Mayfield RD: Understanding Alcoholism Through microRNA Signatures in Brains of Human Alcoholics. Front Genet. 2012, 3: 43-doi:10.3389/fgene.2012.00043
Mayfield RD, Harris RA, Schuckit MA: Genetic Factors Influencing Alcohol Dependence. Br J Pharmacol. 2008, 154 (2): 275-287. doi:10.1038/bjp.2008.88
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.00271
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.00109
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.00077
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.00070
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.00106
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.x
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.026
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\
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.1300947
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.2012
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.1073374
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.1128
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.2207
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.019
Allen NJ, Barres BA: Neuroscience: Glia - More Than Just Brain Glue. Nature. 2009, doi:10.1038/457675a
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.
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/ncb1986
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.
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.x
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-4
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.
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.250
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.006
Jahn R, Fasshauer D: Molecular Machines Governing Exocytosis of Synaptic Vesicles. Nature. 2012, 490 (7419): 201-207. 10.1038/nature11320. doi:10.1038/nature11320
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.201111121
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/157015911795017164
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.200800202
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.061
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.
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.012
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.
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-y
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.x
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.
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.11
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.
Sanna PP, Simpson C, Lutjens R, Koob G: ERK Regulation in Chronic Ethanol Exposure and Withdrawal. Brain Res. 2002, 948 (1–2): 186-191.
Ö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.0000930
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.019
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.x
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.008
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/ddm368
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.
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.x
Hutt D, Balch WE: The Proteome in Balance. Science. 2010, 329 (5993): 766-767. 10.1126/science.1194160. doi:10.1126/science.1194160
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.1191542
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.427336
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.018
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.
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.
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.
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.
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.185124
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.2003
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.
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/38738
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.1094925
Lynch JW, Callister RJ: Glycine Receptors: a New Therapeutic Target in Pain Pathways. Curr Opin Investig Drugs. 2006, 7 (1): 48-53.
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.012
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.t003
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.004
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.1006301107
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/bhq027
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.022
Kalivas PW: The Glutamate Homeostasis Hypothesis of Addiction. Nat Rev Neurosci. 2009, 10 (8): 561-572. 10.1038/nrn2515. doi:10.1038/nrn2515
Koob GF: Drug Addiction: the Yin and Yang of Hedonic Homeostasis. Neuron. 1996, 16: 893-896. 10.1016/S0896-6273(00)80109-9.
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/S1461145708008584
Franklin K, Paxinos G: The Mouse Brain in Streotaxic Coordinates. 2007, San Diego, CA: Academic Press, 1-351.
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_23
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)
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-
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-0496
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/10543400903572753
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/gkr1161
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/gkn851
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.
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.
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.
The authors declare that they have no competing interests.
YON planned and carried out experiments, implemented bioinformatic algorithms for data analysis, and drafted the manuscript. JMT, GG, and ONP planned and carried out experiments and contributed to data analysis. YAB planned and supervised the mouse alcohol consumption studies. RAH supervised experiments and critically reviewed manuscript for important intellectual content. RDM planned and supervised all experiments, analysed data, critically reviewed and gave final approval of the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Figure S1: Random interaction networks generated as control networks. Differentially expressed genes in FCtx of ethanol-treated mice. Table S2. Overrepresentation analysis of miRNA targets among upregulated genes. Table S3. Overrepresentation analysis of miRNA targets among downregulated genes. Table S4. Cell type-specific modular enrichment. Table S5. Functional enrichment analysis for red module genes with high GS and high MM. Table S6. Functional enrichment analysis for brown module genes with high GS and high MM. Table S7. Functional enrichment analysis for turquoise module genes with high GS and high MM. Table S8. Prediction of miRNA-mRNA interactions based on expression correlation patterns between single mRNAs and single differentially expressed (alcohol-responsive) miRNAs. (DOCX 2 MB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Nunez, Y.O., Truitt, J.M., Gorini, G. et al. Positively correlated miRNA-mRNA regulatory networks in mouse frontal cortex during early stages of alcohol dependence. BMC Genomics 14, 725 (2013). https://doi.org/10.1186/1471-2164-14-725
- Expression profiling
- Network analysis
- Alcohol dependence
- Frontal cortex