Gut microbiome-brain-cirrhosis axis.

Cirrhosis is characterized by inflammation, degeneration, and fibrosis of liver tissue. Along with being the most common cause of liver failure and liver transplant, cirrhosis is a significant risk factor for several neuropsychiatric conditions. The most common of these is HE, which is characterized by cognitive and ataxic symptoms, resulting from the buildup of metabolic toxins with liver failure. However, cirrhosis patients also show a significantly increased risk for neurodegenerative diseases such as Alzheimer and Parkinson diseases, and for mood disorders such as anxiety and depression. In recent years, more attention has been played to communication between the ways the gut and liver communicate with each other and with the central nervous system, and the way these organs influence each other's function. This bidirectional communication has come to be known as the gut-liver-brain axis. The gut microbiome has emerged as a key mechanism affecting gut-liver, gut-brain, and brain-liver communication. Clinical studies and animal models have demonstrated the significant patterns of gut dysbiosis when cirrhosis is present, both with or without concomitant alcohol use disorder, and have provided compelling evidence that this dysbiosis also influences the cognitive and mood-related behaviors. In this review, we have summarized the pathophysiological and cognitive effects associated with cirrhosis, links to cirrhosis-associated disruption of the gut microbiome, and the current evidence from clinical and preclinical studies for the modulation of the gut microbiome as a treatment for cirrhosis and associated neuropsychiatric conditions.

Identification of candidate genes for nicotine withdrawal in C57BL/6J × DBA/2J recombinant inbred mice.

Nicotine is the reinforcing ingredient in tobacco. Following chronic exposure, sudden cessation of nicotine use produces negative symptoms of withdrawal that contribute to dependence. The molecular mechanisms underlying nicotine withdrawal behaviors, however, are poorly understood. Using recombinant inbred mice, chronic nicotine was delivered by minipump and withdrawal induced using mecamylamine. Somatic signs of withdrawal, and anxiety-like behavior using elevated plus maze, were then assessed. Interval mapping was used to identify associations between genetic variation and withdrawal behaviors, and with basal gene expression. Differential gene expression following nicotine exposure and withdrawal was also assessed in progenitor mice using microarrays. Quantitative trait loci mapping identified chromosome intervals with significant genetic associations to somatic signs of withdrawal or withdrawal-induced anxiety-like behavior. Using bioinformatics, and association with basal gene expression in nucleus accumbens, we implicated Rb1, Bnip3l, Pnma2, Itm2b, and Kif13b as candidate genes for somatic signs of withdrawal, and Galr1, which showed trans-regulation from a region of chromosome 14 that was associated with somatic signs of withdrawal. Candidate genes within the chromosome 9 region associated with anxiety-like withdrawal behavior included Dixdc1, Ncam1, and Sorl1. Bioinformatics identified six genes that were also significantly associated with nicotine or alcohol traits in recent human genome-wide association studies. Withdrawal-associated somatic signs and anxiety-like behavior had strong non-overlapping genetic associations, respectively, with regions of chromosome 14 and chromosome 9. Genetic, behavioral and gene expression correlations, and bioinformatics analysis identified several candidate genes that may represent novel molecular targets for modulating nicotine withdrawal symptoms.

Different Disease Endotypes in Phenotypically Similar Vasculitides Affecting Small-to-Medium Sized Blood Vessels.

Objectives: Chronic primary vasculitis describes a group of complex and rare diseases that are characterized by blood vessel inflammation. Classification of vasculitis subtypes is based predominantly on the size of the involved vessels and clinical phenotype. There is a recognized need to improve classification, especially for small-to-medium sized vessel vasculitides, that, ideally, is based on the underlying biology with a view to informing treatment. Methods: We performed RNA-Seq on blood samples from children (n = 41) and from adults (n = 11) with small-to-medium sized vessel vasculitis, and used unsupervised hierarchical clustering of gene expression patterns in combination with clinical metadata to define disease subtypes. Results: Differential gene expression at the time of diagnosis separated patients into two primary endotypes that differed in the expression of ~3,800 genes in children, and ~1,600 genes in adults. These endotypes were also present during disease flares, and both adult and pediatric endotypes could be discriminated based on the expression of just 20 differentially expressed genes. Endotypes were associated with distinct biological processes, namely neutrophil degranulation and T cell receptor signaling. Conclusions: Phenotypically similar subsets of small-to-medium sized vessel vasculitis may have different mechanistic drivers involving innate vs. adaptive immune processes. Discovery of these differentiating immune features provides a mechanistic-based alternative for subclassification of vasculitis.

Ivacaftor or lumacaftor/ivacaftor treatment does not alter the core CF airway epithelial gene response to rhinovirus.

BACKGROUND: Aberrant responses by the cystic fibrosis airway epithelium during viral infection may underly the clinical observations. Whether CFTR modulators affect antiviral responses by CF epithelia is presently unknown. We tested the hypothesis that treatment of CF epithelial cells with ivacaftor (Iva) or ivacaftor/lumacaftor (Iva/Lum) would improve control of rhinovirus infection. METHODS: Nineteen CF epithelial cultures (10 homozygous for p.Phe508del as CFTR Class 2, 9 p.Phe508del/p.Gly551Asp as Class 3) were infected with rhinovirus 1B at multiplicity of infection 12 for 24 h. Culture RNA and supernatants were harvested to assess gene and protein expression respectively. RESULTS: RNA-seq analysis comparing rhinovirus infected cultures to control identified 796 and 629 differentially expressed genes for Class 2 and Class 3, respectively. This gene response was highly conserved when cells were treated with CFTR modulators and were predicted to be driven by the same interferon-pathway transcriptional regulators (IFNA, IFNL1, IFNG, IRF7, STAT1). Direct comparisons between treated and untreated infected cultures did not yield any differentially expressed genes for Class 3 and only 68 genes for Class 2. Changes were predominantly related to regulators of lipid metabolism and inflammation, aspects of epithelial biology known to be dysregulated in CF. In addition, CFTR modulators did not affect viral copy number, or levels of pro-inflammatory cytokines produced post-infection. CONCLUSIONS: Though long-term clinical data is not yet available, results presented here suggest that first generation CFTR modulators do not interfere with core airway epithelial responses to rhinovirus infection. Future work should investigate the latest triple modulation therapies.

The Small RNAs PA2952.1 and PrrH as Regulators of Virulence, Motility, and Iron Metabolism in Pseudomonas aeruginosa.

Pseudomonas aeruginosa is a Gram-negative opportunistic pathogen that undergoes swarming motility in response to semisolid conditions with amino acids as a nitrogen source. With a genome encoding hundreds of potential intergenic small RNAs (sRNAs), P. aeruginosa can easily adapt to different conditions and stresses. We previously identified 20 sRNAs that were differentially expressed (DE) under swarming conditions. Here, these sRNAs were overexpressed in strain PAO1 and were subjected to an array of phenotypic screens. Overexpression of the PrrH sRNA resulted in decreased swimming motility, whereas a ΔprrH mutant had decreased cytotoxicity and increased pyoverdine production. Overexpression of the previously uncharacterized PA2952.1 sRNA resulted in decreased swarming and swimming motilities, increased gentamicin and tobramycin resistance under swarming conditions, and increased trimethoprim susceptibility. Transcriptome sequencing (RNA-Seq) and proteomic analysis were performed on the wild type (WT) overexpressing PA2952.1 compared to the empty vector control under swarming conditions, and these revealed the differential expression (absolute fold change [FC] ≥ 1.5) of 784 genes and the differential abundance (absolute FC ≥ 1.25) of 59 proteins. Among these were found 73 transcriptional regulators, two-component systems, and sigma and anti-sigma factors. Downstream effectors included downregulated pilus and flagellar genes, the upregulated efflux pump MexGHI-OpmD, and the upregulated arn operon. Genes involved in iron and zinc uptake were generally upregulated, and certain pyoverdine genes were upregulated. Overall, the sRNAs PA2952.1 and PrrH appeared to be involved in regulating virulence-related programs in P. aeruginosa, including iron acquisition and motility.IMPORTANCE Due to the rising incidence of multidrug-resistant (MDR) strains and the difficulty of eliminating P. aeruginosa infections, it is important to understand the regulatory mechanisms that allow this bacterium to adapt to and thrive under a variety of conditions. Small RNAs (sRNAs) are one regulatory mechanism that allows bacteria to change the amount of protein synthesized. In this study, we overexpressed 20 different sRNAs in order to investigate how this might affect different bacterial behaviors. We found that one of the sRNAs, PrrH, played a role in swimming motility and virulence phenotypes, indicating a potentially important role in clinical infections. Another sRNA, PA2952.1, affected other clinically relevant phenotypes, including motility and antibiotic resistance. RNA-Seq and proteomics of the strain overexpressing PA2952.1 revealed the differential expression of 784 genes and 59 proteins, with a total of 73 regulatory factors. This substantial dysregulation indicates an important role for the sRNA PA2952.1.

Unique inducible filamentous motility identified in pathogenic Bacillus cereus group species.

Active migration across semi-solid surfaces is important for bacterial success by facilitating colonization of unoccupied niches and is often associated with altered virulence and antibiotic resistance profiles. We isolated an atmospheric contaminant, subsequently identified as a new strain of Bacillus mobilis, which showed a unique, robust, rapid, and inducible filamentous surface motility. This flagella-independent migration was characterized by formation of elongated cells at the expanding edge and was induced when cells were inoculated onto lawns of metabolically inactive Campylobacter jejuni cells, autoclaved bacterial biomass, adsorbed milk, and adsorbed blood atop hard agar plates. Phosphatidylcholine (PC), bacterial membrane components, and sterile human fecal extracts were also sufficient to induce filamentous expansion. Screening of eight other Bacillus spp. showed that filamentous motility was conserved amongst B. cereus group species to varying degrees. RNA-Seq of elongated expanding cells collected from adsorbed milk and PC lawns versus control rod-shaped cells revealed dysregulation of genes involved in metabolism and membrane transport, sporulation, quorum sensing, antibiotic synthesis, and virulence (e.g., hblA/B/C/D and plcR). These findings characterize the robustness and ecological significance of filamentous surface motility in B. cereus group species and lay the foundation for understanding the biological role it may play during environment and host colonization.

Overexpression of the Small RNA PA0805.1 in Pseudomonas aeruginosa Modulates the Expression of a Large Set of Genes and Proteins, Resulting in Altered Motility, Cytotoxicity, and Tobramycin Resistance.

Pseudomonas aeruginosa is a motile species that initiates swarming motility in response to specific environmental cues, i.e., a semisolid surface with amino acids as a nitrogen source (relevant to the human lung). Swarming is an intricately regulated process, but to date posttranscriptional regulation has not been extensively investigated. Small noncoding RNAs (sRNAs) are hypothesized to play posttranscriptional regulatory roles, largely through suppression of translation, and we previously demonstrated 20 sRNA species that were dysregulated under swarming conditions. One of these, sRNA PA0805.1 (which was 5-fold upregulated under swarming conditions), when cloned, transformed into wild-type (WT) PAO1, and overexpressed, led to broad phenotypic changes, including reduced swarming, swimming, and twitching motilities, as well as increased adherence, cytotoxicity, and tobramycin resistance. A ΔPA0805.1 deletion mutant was more susceptible to tobramycin than the WT under swarming conditions. The strain overexpressing PA0805.1 was compared to the empty-vector strain by transcriptome sequencing (RNA-Seq) and proteomics under swarming conditions to determine sRNA targets. Broad transcriptional and proteomic profiles showed 1,121 differentially expressed genes and 258 proteins with significantly different abundance. Importantly, these included 106 transcriptional regulators, two-component regulatory systems, and sigma and anti-sigma factors. Downstream of these regulators were found downregulated type IV pilus genes, many upregulated adherence and virulence factors, and two multidrug efflux systems, mexXY and mexGHI-opmD Therefore, the sRNA PA0805.1 appears to be a global regulator that influences diverse bacterial lifestyles, most likely through a regulatory cascade.IMPORTANCE P. aeruginosa is an opportunistic pathogen of humans. With roughly 10% of its genes encoding transcriptional regulators, and hundreds of small noncoding RNAs (sRNAs) interspersed throughout the genome, P. aeruginosa is able to fine-tune its response to adapt and survive in the host and resist antimicrobial agents. Understanding mechanisms of genetic regulation is therefore crucial to combat pathogenesis. The previously uncharacterized sRNA PA0805.1 was overexpressed in P. aeruginosa strain PAO1, resulting in decreased motility, increased adherence, cytotoxicity, and tobramycin resistance. In contrast, a ΔPA0805.1 deletion mutant had increased susceptibility to tobramycin under swarming conditions. Omic approaches uncovered 1,121 transcriptomic and 258 proteomic changes in the overexpression strain compared with the empty-vector strain, which included 106 regulatory factors. Downstream of these regulators were upregulated adherence factors, multidrug efflux systems, and virulence factors in both transcriptomics and proteomics. This study provides insights into the role of the sRNA PA0805.1 in modulating bacterial adaptations.

Brain regional gene expression network analysis identifies unique interactions between chronic ethanol exposure and consumption.

Progressive increases in ethanol consumption is a hallmark of alcohol use disorder (AUD). Persistent changes in brain gene expression are hypothesized to underlie the altered neural signaling producing abusive consumption in AUD. To identify brain regional gene expression networks contributing to progressive ethanol consumption, we performed microarray and scale-free network analysis of expression responses in a C57BL/6J mouse model utilizing chronic intermittent ethanol by vapor chamber (CIE) in combination with limited access oral ethanol consumption. This model has previously been shown to produce long-lasting increased ethanol consumption, particularly when combining oral ethanol access with repeated cycles of intermittent vapor exposure. The interaction of CIE and oral consumption was studied by expression profiling and network analysis in medial prefrontal cortex, nucleus accumbens, hippocampus, bed nucleus of the stria terminalis, and central nucleus of the amygdala. Brain region expression networks were analyzed for ethanol-responsive gene expression, correlation with ethanol consumption and functional content using extensive bioinformatics studies. In all brain-regions studied the largest number of changes in gene expression were seen when comparing ethanol naïve mice to those exposed to CIE and drinking. In the prefrontal cortex, however, unique patterns of gene expression were seen compared to other brain-regions. Network analysis identified modules of co-expressed genes in all brain regions. The prefrontal cortex and nucleus accumbens showed the greatest number of modules with significant correlation to drinking behavior. Across brain-regions, however, many modules with strong correlations to drinking, both baseline intake and amount consumed after CIE, showed functional enrichment for synaptic transmission and synaptic plasticity.

Humoral immunity prevents clinical malaria during Plasmodium relapses without eliminating gametocytes.

Plasmodium relapses are attributed to the activation of dormant liver-stage parasites and are responsible for a significant number of recurring malaria blood-stage infections. While characteristic of human infections caused by P. vivax and P. ovale, their relative contribution to malaria disease burden and transmission remains poorly understood. This is largely because it is difficult to identify 'bona fide' relapse infections due to ongoing transmission in most endemic areas. Here, we use the P. cynomolgi-rhesus macaque model of relapsing malaria to demonstrate that clinical immunity can form after a single sporozoite-initiated blood-stage infection and prevent illness during relapses and homologous reinfections. By integrating data from whole blood RNA-sequencing, flow cytometry, P. cynomolgi-specific ELISAs, and opsonic phagocytosis assays, we demonstrate that this immunity is associated with a rapid recall response by memory B cells that expand and produce anti-parasite IgG1 that can mediate parasite clearance of relapsing parasites. The reduction in parasitemia during relapses was mirrored by a reduction in the total number of circulating gametocytes, but importantly, the cumulative proportion of gametocytes increased during relapses. Overall, this study reveals that P. cynomolgi relapse infections can be clinically silent in macaques due to rapid memory B cell responses that help to clear asexual-stage parasites but still carry gametocytes.

Cross-Species Co-analysis of Prefrontal Cortex Chronic Ethanol Transcriptome Responses in Mice and Monkeys.

Despite recent extensive genomic and genetic studies on behavioral responses to ethanol, relatively few new therapeutic targets for the treatment of alcohol use disorder have been validated. Here, we describe a cross-species genomic approach focused on identifying gene networks associated with chronic ethanol consumption. To identify brain mechanisms underlying a chronic ethanol consumption phenotype highly relevant to human alcohol use disorder, and to elucidate potential future therapeutic targets, we conducted a genomic study in a non-human primate model of chronic open-access ethanol consumption. Microarray analysis of RNA expression in anterior cingulate and subgenual cortices from rhesus macaques was performed across multiple cohorts of animals. Gene networks correlating with ethanol consumption or showing enrichment for ethanol-regulated genes were identified, as were major ethanol-related hub genes within these networks. A subsequent consensus module analysis was used to co-analyze monkey data with expression data from a chronic intermittent ethanol vapor-exposure and consumption model in C57BL/6J mice. Ethanol-related gene networks conserved between primates and rodents were enriched for genes involved in discrete biological functions, including; myelination, synaptic transmission, chromatin modification, Golgi apparatus function, translation, cellular respiration, and RNA processing. The myelin-related network, in particular, showed strong correlations with ethanol consumption behavior and displayed marked network reorganization between control and ethanol-drinking animals. Further bioinformatics analysis revealed that these networks also showed highly significant overlap with other ethanol-regulated gene sets. Altogether, these studies provide robust primate and rodent cross-species validation of gene networks associated with chronic ethanol consumption. Our results also suggest potential novel focal points for future therapeutic interventions in alcohol use disorder.

Correction: Time-Course Analysis of Brain Regional Expression Network Responses to Chronic Intermittent Ethanol and Withdrawal: Implications for Mechanisms Underlying Excessive Ethanol Consumption.

[This corrects the article DOI: 10.1371/journal.pone.0146257.].

Systems Biology-Based Investigation of Host-Plasmodium Interactions.

Malaria is a serious, complex disease caused by parasites of the genus Plasmodium. Plasmodium parasites affect multiple tissues as they evade immune responses, replicate, sexually reproduce, and transmit between vertebrate and invertebrate hosts. The explosion of omics technologies has enabled large-scale collection of Plasmodium infection data, revealing systems-scale patterns, mechanisms of pathogenesis, and the ways that host and pathogen affect each other. Here, we provide an overview of recent efforts using systems biology approaches to study host-Plasmodium interactions and the biological themes that have emerged from these efforts. We discuss some of the challenges in using systems biology for this goal, key research efforts needed to address those issues, and promising future malaria applications of systems biology.

Differential potassium channel gene regulation in BXD mice reveals novel targets for pharmacogenetic therapies to reduce heavy alcohol drinking.

Alcohol (ethanol) dependence is a chronic relapsing brain disorder partially influenced by genetics and characterized by an inability to regulate harmful levels of drinking. Emerging evidence has linked genes that encode KV7, KIR, and KCa2 K+ channels with variation in alcohol-related behaviors in rodents and humans. This led us to experimentally test relations between K+ channel genes and escalation of drinking in a chronic-intermittent ethanol (CIE) exposure model of dependence in BXD recombinant inbred strains of mice. Transcript levels for K+ channel genes in the prefrontal cortex (PFC) and nucleus accumbens (NAc) covary with voluntary ethanol drinking in a non-dependent cohort. Transcripts that encode KV7 channels covary negatively with drinking in non-dependent BXD strains. Using a pharmacological approach to validate the genetic findings, C57BL/6J mice were allowed intermittent access to ethanol to establish baseline consumption before they were treated with retigabine, an FDA-approved KV7 channel positive modulator. Systemic administration significantly reduced drinking, and consistent with previous evidence, retigabine was more effective at reducing voluntary consumption in high-drinking than low-drinking subjects. We evaluated the specific K+ channel genes that were most sensitive to CIE exposure and identified a gene subset in the NAc and PFC that were dysregulated in the alcohol-dependent BXD cohort. CIE-induced modulation of nine genes in the NAc and six genes in the PFC covaried well with the changes in drinking induced by ethanol dependence. Here we identified novel candidate genes in the NAc and PFC that are regulated by ethanol dependence and correlate with voluntary drinking in non-dependent and dependent BXD mice. The findings that Kcnq expression correlates with drinking and that retigabine reduces consumption suggest that KV7 channels could be pharmacogenetic targets to treat individuals with alcohol addiction.

The allostatic impact of chronic ethanol on gene expression: A genetic analysis of chronic intermittent ethanol treatment in the BXD cohort.

The transition from acute to chronic ethanol exposure leads to lasting behavioral and physiological changes such as increased consumption, dependence, and withdrawal. Changes in brain gene expression are hypothesized to underlie these adaptive responses to ethanol. Previous studies on acute ethanol identified genetic variation in brain gene expression networks and behavioral responses to ethanol across the BXD panel of recombinant inbred mice. In this work, we have performed the first joint genetic and genomic analysis of transcriptome shifts in response to chronic intermittent ethanol (CIE) by vapor chamber exposure in a BXD cohort. CIE treatment is known to produce significant and sustained changes in ethanol consumption with repeated cycles of ethanol vapor. Using Affymetrix microarray analysis of prefrontal cortex (PFC) and nucleus accumbens (NAC) RNA, we compared CIE expression responses to those seen following acute ethanol treatment, and to voluntary ethanol consumption. Gene expression changes in PFC and NAC after CIE overlapped significantly across brain regions and with previously published expression following acute ethanol. Genes highly modulated by CIE were enriched for specific biological processes including synaptic transmission, neuron ensheathment, intracellular signaling, and neuronal projection development. Expression quantitative trait locus (eQTL) analyses identified genomic loci associated with ethanol-induced transcriptional changes with largely distinct loci identified between brain regions. Correlating CIE-regulated genes to ethanol consumption data identified specific genes highly associated with variation in the increase in drinking seen with repeated cycles of CIE. In particular, multiple myelin-related genes were identified. Furthermore, genetic variance in or near dynamin3 (Dnm3) on Chr1 at ∼164 Mb may have a major regulatory role in CIE-responsive gene expression. Dnm3 expression correlates significantly with ethanol consumption, is contained in a highly ranked functional group of CIE-regulated genes in the NAC, and has a cis-eQTL within a genomic region linked with multiple CIE-responsive genes.

Genetic divergence in the transcriptional engram of chronic alcohol abuse: A laser-capture RNA-seq study of the mouse mesocorticolimbic system.

Genetic factors that influence the transition from initial drinking to dependence remain enigmatic. Recent studies have leveraged chronic intermittent ethanol (CIE) paradigms to measure changes in brain gene expression in a single strain at 0, 8, 72 h, and even 7 days following CIE. We extend these findings using LCM RNA-seq to profile expression in 11 brain regions in two inbred strains - C57BL/6J (B6) and DBA/2J (D2) - 72 h following multiple cycles of ethanol self-administration and CIE. Linear models identified differential expression based on treatment, region, strain, or interactions with treatment. Nearly 40% of genes showed a robust effect (FDR < 0.01) of region, and hippocampus CA1, cortex, bed nucleus stria terminalis, and nucleus accumbens core had the highest number of differentially expressed genes after treatment. Another 8% of differentially expressed genes demonstrated a robust effect of strain. As expected, based on similar studies in B6, treatment had a much smaller impact on expression; only 72 genes (p < 0.01) are modulated by treatment (independent of region or strain). Strikingly, many more genes (415) show a strain-specific and largely opposite response to treatment and are enriched in processes related to RNA metabolism, transcription factor activity, and mitochondrial function. Over 3 times as many changes in gene expression were detected in D2 compared to B6, and weighted gene co-expression network analysis (WGCNA) module comparison identified more modules enriched for treatment effects in D2. Substantial strain differences exist in the temporal pattern of transcriptional neuroadaptation to CIE, and these may drive individual differences in risk of addiction following excessive alcohol consumption.

Integrative Analysis of Genetic, Genomic, and Phenotypic Data for Ethanol Behaviors: A Network-Based Pipeline for Identifying Mechanisms and Potential Drug Targets.

Complex behavioral traits, such as alcohol abuse, are caused by an interplay of genetic and environmental factors, producing deleterious functional adaptations in the central nervous system. The long-term behavioral consequences of such changes are of substantial cost to both the individual and society. Substantial progress has been made in the last two decades in understanding elements of brain mechanisms underlying responses to ethanol in animal models and risk factors for alcohol use disorder (AUD) in humans. However, treatments for AUD remain largely ineffective and few medications for this disease state have been licensed. Genome-wide genetic polymorphism analysis (GWAS) in humans, behavioral genetic studies in animal models and brain gene expression studies produced by microarrays or RNA-seq have the potential to produce nonbiased and novel insight into the underlying neurobiology of AUD. However, the complexity of such information, both statistical and informational, has slowed progress toward identifying new targets for intervention in AUD. This chapter describes one approach for integrating behavioral, genetic, and genomic information across animal model and human studies. The goal of this approach is to identify networks of genes functioning in the brain that are most relevant to the underlying mechanisms of a complex disease such as AUD. We illustrate an example of how genomic studies in animal models can be used to produce robust gene networks that have functional implications, and to integrate such animal model genomic data with human genetic studies such as GWAS for AUD. We describe several useful analysis tools for such studies: ComBAT, WGCNA, and EW_dmGWAS. The end result of this analysis is a ranking of gene networks and identification of their cognate hub genes, which might provide eventual targets for future therapeutic development. Furthermore, this combined approach may also improve our understanding of basic mechanisms underlying gene x environmental interactions affecting brain functioning in health and disease.

Time-Course Analysis of Brain Regional Expression Network Responses to Chronic Intermittent Ethanol and Withdrawal: Implications for Mechanisms Underlying Excessive Ethanol Consumption.

Long lasting abusive consumption, dependence, and withdrawal are characteristic features of alcohol use disorders (AUD). Mechanistically, persistent changes in gene expression are hypothesized to contribute to brain adaptations leading to ethanol toxicity and AUD. We employed repeated chronic intermittent ethanol (CIE) exposure by vapor chamber as a mouse model to simulate the cycles of ethanol exposure and withdrawal commonly seen with AUD. This model has been shown to induce progressive ethanol consumption in rodents. Brain CIE-responsive expression networks were identified by microarray analysis across five regions of the mesolimbic dopamine system and extended amygdala with tissue harvested from 0-hours to 7-days following CIE. Weighted Gene Correlated Network Analysis (WGCNA) was used to identify gene networks over-represented for CIE-induced temporal expression changes across brain regions. Differential gene expression analysis showed that long-lasting gene regulation occurred 7-days after the final cycle of ethanol exposure only in prefrontal cortex (PFC) and hippocampus. Across all brain regions, however, ethanol-responsive expression changes occurred mainly within the first 8-hours after removal from ethanol. Bioinformatics analysis showed that neuroinflammatory responses were seen across multiple brain regions at early time-points, whereas co-expression modules related to neuroplasticity, chromatin remodeling, and neurodevelopment were seen at later time-points and in specific brain regions (PFC or HPC). In PFC a module containing Bdnf was identified as highly CIE responsive in a biphasic manner, with peak changes at 0 hours and 5 days following CIE, suggesting a possible role in mechanisms underlying long-term molecular and behavioral response to CIE. Bioinformatics analysis of this network and several other modules identified Let-7 family microRNAs as potential regulators of gene expression changes induced by CIE. Our results suggest a complex temporal and regional pattern of widespread gene network responses involving neuroinflammatory and neuroplasticity related genes as contributing to physiological and behavioral responses to chronic ethanol.