Integrating mRNA and miRNA Weighted Gene Co-Expression Networks with eQTLs in the Nucleus Accumbens of Subjects with Alcohol Dependence

Space radiation damage rescued by inhibition of key spaceflight associated miRNAs

Our previous research revealed a key microRNA signature that is associated with spaceflight that can be used as a biomarker and to develop countermeasure treatments to mitigate the damage caused by space radiation. Here, we expand on this work to determine the biological factors rescued by the countermeasure treatment. We performed RNA-sequencing and transcriptomic analysis on 3D microvessel cell cultures exposed to simulated deep space radiation (0.5 Gy of Galactic Cosmic Radiation) with and without the antagonists to three microRNAs: miR-16-5p, miR-125b-5p, and let-7a-5p (i.e., antagomirs). Significant reduction of inflammation and DNA double strand breaks (DSBs) activity and rescue of mitochondria functions are observed after antagomir treatment. Using data from astronaut participants in the NASA Twin Study, Inspiration4, and JAXA missions, we reveal the genes and pathways implicated in the action of these antagomirs are altered in humans. Our findings indicate a countermeasure strategy that can potentially be utilized by astronauts in spaceflight missions to mitigate space radiation damage.

Mapping microRNA expression quantitative trait loci in the prenatal human brain implicates miR-1908-5p expression in bipolar disorder and other brain-related traits

MicroRNA (miRNA) are small non-coding RNA involved in post-transcriptional gene regulation. Given their known involvement in early neurodevelopment processes, we here sought to identify common genetic variants associated with altered miRNA expression in the prenatal human brain. We performed small RNA sequencing on brain tissue from 112 genome-wide genotyped fetuses from the second trimester of gestation, identifying high-confidence (false discovery rate < 0.05) expression quantitative trait loci for 30 mature miRNA. Integrating our findings with genome-wide association study data for brain-related disorders, we implicate increased prenatal expression of miR-1908-5p as a risk mechanism for bipolar disorder and find that predicted mRNA targets of miR-1908-5p that are expressed in the fetal brain are enriched for common variant genetic association with the condition. Extending these analyses to other brain-related traits, we find that common genetic variation associated with increased miR-1908-5p expression in fetal brain is additionally associated with depressive symptoms, irritability, increased right cerebellum exterior volume and increased sleep duration in the general population. Our findings provide support to the view that altered miRNA expression can influence susceptibility to neuropsychiatric illness and suggest an early neurodevelopmental risk mechanism for bipolar disorder.

Addictive drugs modify neurogenesis, synaptogenesis and synaptic plasticity to impair memory formation through neurotransmitter imbalances and signaling dysfunction

Drug abuse changes neurophysiological functions at multiple cellular and molecular levels in the addicted brain. Well-supported scientific evidence suggests that drugs negatively affect memory formation, decision-making and inhibition, and emotional and cognitive behaviors. The mesocorticolimbic brain regions are involved in reward-related learning and habitual drug-seeking/taking behaviors to develop physiological and psychological dependence on the drugs. This review highlights the importance of specific drug-induced chemical imbalances resulting in memory impairment through various neurotransmitter receptor-mediated signaling pathways. The mesocorticolimbic modifications in the expression levels of brain-derived neurotrophic factor (BDNF) and the cAMP-response element binding protein (CREB) impair reward-related memory formation following drug abuse. The contributions of protein kinases and microRNAs (miRNAs), along with the transcriptional and epigenetic regulation have also been considered in memory impairment underlying drug addiction. Overall, we integrate the research on various types of drug-induced memory impairment in distinguished brain regions and provide a comprehensive review with clinical implications addressing the upcoming studies.

Analysis of the brain transcriptome for substance-associated genes: An update on large-scale genome-wide association studies

Substance use disorder (SUD) arises from the initiation to subsequent regular, irregular and harmful use of substances such as alcohol, tobacco/nicotine and cannabis. While thousands of genetic variants have been identified from recent large-scale genome-wide association studies (GWAS), understanding their functions in substance involvement and investigating the mechanisms by which they act in the addiction circuits remains challenging. In this study, we re-analysed the brain regional transcriptome data from the most comprehensive postmortem transcriptomic study, with a focus on up- or down-regulation of substance-associated protein-coding genes in the addiction circuit-related brain regions (AddictRegions), including six cortical and 11 subcortical regions. We found that substance-associated genes were overrepresented in AddictRegions. Interestingly, we observed a greater degree of genetic overlap between initiation and use and between use and SUD than between initiation and SUD. Moreover, substance initiation, use and SUD-associated genes tend to shift their enriched AddictRegions from the cortical for initiation and, to a lesser extent, substance use to subcortical regions for SUD (e.g., thalamus, substantia nigra and ventral tegmental area). We also uncovered a pattern of coordinated cortical up-regulation and subcortical down-regulation for the genes associated with tobacco initiation and alcohol use. Moreover, the Gene Ontology terms of glutamate receptor activity and neurotransmitter binding were most significantly overrepresented in AddictRegion-upregulated, substance-associated genes, with the strongest enrichment for those involved in common substance use behaviours. Overall, our analysis provides a resource of AddictRegion-related transcriptomes for studying substance-associated genes and generates intriguing insights into the genetic control of substance initiation, use and SUD.

The Effects of Transcranial Focused Ultrasound Stimulation of Nucleus Accumbens on Neuronal Gene Expression and Brain Tissue in High Alcohol-Preferring Rats

We investigated the effect of low-intensity focused ultrasound (LIFU) on gene expression related to alcohol dependence and histological effects on brain tissue. We also aimed at determining the miRNA-mRNA relationship and their pathways in alcohol dependence-induced expression changes after focused ultrasound therapy. We designed a case-control study for 100 days of observation to investigate differences in gene expression in the short-term stimulation group (STS) and long-term stimulation group (LTS) compared with the control sham group (SG). The study was performed in our Experimental Research Laboratory. 24 male high alcohol-preferring rats 63 to 79 days old, weighing 270 to 300 g, were included in the experiment. LTS received 50-day LIFU and STS received 10-day LIFU and 40-day sham stimulation, while the SG received 50-day sham stimulation. In miRNA expression analysis, it was found that LIFU caused gene expression differences in NAc. Significant differences were found between the groups for gene expression. Compared to the SG, the expression of 454 genes in the NAc region was changed in the STS while the expression of 382 genes was changed in the LTS. In the LTS, the expression of 32 genes was changed in total compared to STS. Our data suggest that LIFU targeted on NAc may assist in the treatment of alcohol dependence, especially in the long term possibly through altering gene expression. Our immunohistochemical studies verified that LIFU does not cause any tissue damage. These findings may lead to new studies in investigating the efficacy of LIFU for the treatment of alcohol dependence and also for other psychiatric disorders.

Proteomics and weighted gene correlated network analysis reveal glutamatergic synapse signaling in diazepam treatment of alcohol withdrawal

Background: Alcohol use disorder (AUD) is characterized by chronic excessive alcohol consumption, often alternating with periods of abstinence known as alcohol withdrawal syndrome (AWS). Diazepam is the preferred benzodiazepine for treatment of alcohol withdrawal syndrome under most circumstances, but the specific mechanism underlying the treatment needs further research. Methods: We constructed an animal model of two-bottle choices and chronic intermittent ethanol exposure. LC-MS/MS proteomic analysis based on the label-free and intensity-based quantification approach was used to detect the protein profile of the whole brain. Weighted gene correlated network analysis was applied for scale-free network topology analysis. We established a protein-protein interaction network based on the Search Tool for the Retrieval of Interacting Genes (STRING) database and Cytoscape software and identified hub proteins by CytoHubba and MCODE plugins of Cytoscape. The online tool Targetscan identified miRNA-mRNA pair interactions. Results: Seven hub proteins (Dlg3, Dlg4, Shank3, Grin2b, Camk2b, Camk2a and Syngap1) were implicated in alcohol withdrawal syndrome or diazepam treatment. In enrichment analysis, glutamatergic synapses were considered the most important pathway related to alcohol use disorder. Decreased glutamatergic synapses were observed in the late stage of withdrawal, as a protective mechanism that attenuated withdrawal-induced excitotoxicity. Diazepam treatment during withdrawal increased glutamatergic synapses, alleviating withdrawal-induced synapse inhibition. Conclusion: Glutamatergic synapses are considered the most important pathway related to alcohol use disorder that may be a potential molecular target for new interventional strategies.

Integrative analysis of microRNAs and mRNAs reveals the regulatory networks of triterpenoid saponin metabolism in Soapberry (Sapindus mukorossi Gaertn.)

Triterpenoid saponin are important secondary metabolites and bioactive constituents of soapberry (Sapindus mukorossi Gaertn.) and are widely used in medicine and toiletry products. However, little is known about the roles of miRNAs in the regulation of triterpenoid saponin biosynthesis in soapberry. In this study, a total of 3036 miRNAs were identified, of which 1372 miRNAs were differentially expressed at different stages of pericarp development. Important KEGG pathways, such as terpenoid backbone biosynthesis, sesquiterpenoid and triterpenoid biosynthesis, and basal transcription factors were highlighted, as well the roles of some key miRNAs, such as ath-miR5021, han-miR3630-3p, and ppe-miR858, which may play important roles in regulating triterpenoid saponin biosynthesis. In addition, 58 miRNAs might participate in saponin biosynthesis pathways by predicting the targets of those miRNAs to 53 saponin biosynthesis structural genes. And 75 miRNAs were identified to potentially play vital role in saponin accumulation by targeting transcript factor genes, bHLH, bZIP, ERF, MYB, and WRKY, respectively, which are candidate regulatory genes in the pathway of saponin biosynthesis. The results of weighted gene coexpression network analysis (WGCNA) suggested that two saponin-specific miRNA modules and 10 hub miRNAs may participate in saponin biosynthesis. Furthermore, multiple miRNA-mRNA regulatory networks potentially involved in saponin biosynthesis were generated, e.g., ath-miR5021-SmIDI2/SmGPS5/SmbAS1/SmCYP71D-3/SmUGT74G-2, han-miR3630-3p-SmCYP71A-14/SmbHLH54/SmMYB135/SmWRKY32, and ppe-miR858-SmMYB5/SmMYB32. qRT-PCR analysis validated the expression patterns of nine miRNAs and 12 corresponding target genes. This study represents the first comprehensive analysis of miRNAs in soapberry and lays the foundation for further understanding of miRNA-based regulation in triterpenoid saponin biosynthesis.

Functional roles, regulatory mechanisms and theranostics applications of ncRNAs in alcohol use disorder

Alcohol use disorder (AUD) is one of the most prevalent neuropsychological disorders worldwide, and its pathogenesis is convoluted and poorly understood. There is considerable evidence demonstrating significant associations between multiple heritable factors and the onset and progression of AUD. In recent years, a substantial body of research conducted by emerging biotechnologies has increasingly highlighted the crucial roles of noncoding RNAs (ncRNAs) in the pathophysiology of mental diseases. As in-depth understanding of ncRNAs and their mechanisms of action, they have emerged as prospective diagnostic indicators and preclinical therapeutic targets for a variety of psychiatric illness, including AUD. Of note, dysregulated expression of ncRNAs such as circRNAs, lncRNAs and miRNAs was routinely found in AUD individuals, and besides, exogenous regulation of partial ncRNAs has also been shown to be effective in ameliorating alcohol preference and excessive alcohol consumption. However, the exact molecular mechanism still remains elusive. Herein, we systematically summarized current knowledge regarding alterations in the expression of certain ncRNAs as well as their-mediated regulatory mechanisms in individuals with AUD. And finally, we detailedly reviewed the potential theranostics applications of gene therapy agents targeting ncRNAs in AUD mice. Overall, a deeper comprehension of functional roles and biological mechanisms of ncRNAs may make significant contributions to the accurate diagnosis and effective treatment of AUD.

Noncoding RNA therapeutics for substance use disorder

Although noncoding RNAs (ncRNAs) have been shown to regulate maladaptive neuroadaptations that drive compulsive drug use, ncRNA-targeting therapeutics for substance use disorder (SUD) have yet to be clinically tested. Recent advances in RNA-based drugs have improved many therapeutic issues related to immune response, specificity, and delivery, leading to multiple successful clinical trials for other diseases. As the need for safe and effective treatments for SUD continues to grow, novel nucleic acid-based therapeutics represent an appealing approach to target ncRNA mechanisms in SUD. Here, we review ncRNA processes implicated in SUD, discuss recent therapeutic approaches for targeting ncRNAs, and highlight potential opportunities and challenges of ncRNA-targeting therapeutics for SUD.

MicroRNAs and Synaptic Plasticity: From Their Molecular Roles to Response to Therapy

Synaptic plasticity is the ability of synapses to weaken or strengthen over time, in response to changes in the activity of the neurons. It is orchestrated by a variety of genes, proteins, and external and internal factors, especially epigenetic factors. MicroRNAs (miRNAs) are well-acknowledged epigenetic modulators that regulate the translation and degradation of target genes in the nervous system. Increasing evidence has suggested that a number of miRNAs play important roles in modulating various aspects of synaptic plasticity. The deregulation of miRNAs could be associated with pathological alterations in synaptic plasticity, which could lead to different CNS-related diseases. Herein, we provide an update on the role of miRNAs in governing synaptic plasticity. In addition, we also summarize recent researches on the role of miRNAs in drug addiction, and their targets and mechanism of action. Understanding of the way in which miRNAs contribute to synaptic plasticity provides rational clues in establishing the novel biomarkers and new therapeutic strategies for the diagnosis and treatment of plasticity-related diseases and drug addiction.

The value of the MicroRNAs on alcoholic addicts: A meta-analytic review

INTRODUCTION

A growing body of evidence indicates that repeated alcohol exposure or withdrawal from alcohol can result in persistent molecular and cellular adaptations. One molecular adaptation that occurs is the regulation of gene expression, which is believed to lead to functional alterations that characterize addiction. MicroRNAs (miRs) have been recently identified as master regulators of gene expression through post-transcriptional regulation. The aim of this meta-analytic review was to evaluate the regulatory forms of miRs during alcoholism.

METHODS

We used several databases such as PubMed, Scopus, and Web of Science without limitations on publication time. All studies were analyzed by Comprehensive Meta-Analysis software.

RESULTS AND DISCUSSION

Six clinical papers with 243 alcoholic patients and 162 controls were included. In this study, 1680 articles were initially reviewed and eventually, six clinical studies were included in the meta-analysis. The results of the meta-analysis showed that according to the random model, the difference between the upregulation and downregulation of central addiction targets was statistically significant, indicating that most dopamine- or gamma-aminobutyric acid receptor subunit (GABA)-related miRs are upregulated in alcoholics (P: 0.00, CI: 0.149-0.439).

CONCLUSION

This study strongly suggests that dopamine- or GABA-related miRs were mostly upregulated in alcoholism. Our findings revealed that about 9% of miRs were downregulated in alcoholism, including miR-567, miR-126, miR-1, miR-432, and miR-153. To identify other or specific miRs as potential biomarkers in alcoholics, large-scale studies and more clinical work are required.

An Augmented High-Dimensional Graphical Lasso Method to Incorporate Prior Biological Knowledge for Global Network Learning

Biological networks are often inferred through Gaussian graphical models (GGMs) using gene or protein expression data only. GGMs identify conditional dependence by estimating a precision matrix between genes or proteins. However, conventional GGM approaches often ignore prior knowledge about protein-protein interactions (PPI). Recently, several groups have extended GGM to weighted graphical Lasso (wGlasso) and network-based gene set analysis (Netgsa) and have demonstrated the advantages of incorporating PPI information. However, these methods are either computationally intractable for large-scale data, or disregard weights in the PPI networks. To address these shortcomings, we extended the Netgsa approach and developed an augmented high-dimensional graphical Lasso (AhGlasso) method to incorporate edge weights in known PPI with omics data for global network learning. This new method outperforms weighted graphical Lasso-based algorithms with respect to computational time in simulated large-scale data settings while achieving better or comparable prediction accuracy of node connections. The total runtime of AhGlasso is approximately five times faster than weighted Glasso methods when the graph size ranges from 1,000 to 3,000 with a fixed sample size (n = 300). The runtime difference between AhGlasso and weighted Glasso increases when the graph size increases. Using proteomic data from a study on chronic obstructive pulmonary disease, we demonstrate that AhGlasso improves protein network inference compared to the Netgsa approach by incorporating PPI information.

Epigenetic modification in alcohol use disorder and alcoholic cardiomyopathy: From pathophysiology to therapeutic opportunities

Alcohol consumption prompts detrimental psychological, pathophysiological and health issues, representing one of the major causes of death worldwide. Alcohol use disorder (AUD), which is characterized by compulsive alcohol intake and loss of control over alcohol usage, arises from a complex interplay between genetic and environmental factors. More importantly, long-term abuse of alcohol is often tied with unfavorable cardiac remodeling and contractile alterations, a cadre of cardiac responses collectively known as alcoholic cardiomyopathy (ACM). Recent evidence has denoted a pivotal role for ethanol-triggered epigenetic modifications, the interface between genome and environmental cues, in the organismal and cellular responses to ethanol exposure. To-date, three major epigenetic mechanisms (DNA methylation, histone modifications, and RNA-based mechanisms) have been identified for the onset and development of AUD and ACM. Importantly, these epigenetic changes induced by alcohol may be detectable in the blood, thus offering diagnostic, therapeutic, and prognostic promises of epigenetic markers for AUD and alcoholic complications. In addition, several epigenetic drugs have shown efficacies in the management of alcohol abuse, loss of control for alcohol usage, relapse, drinking-related anxiety and behavior in withdrawal. In this context, medications targeting epigenetic modifications may hold promises for pharmaceutical management of AUD and ACM.

Identifying a novel biological mechanism for alcohol addiction associated with circRNA networks acting as potential miRNA sponges

Our lab and others have shown that chronic alcohol use leads to gene and miRNA expression changes across the mesocorticolimbic (MCL) system. Circular RNAs (circRNAs) are noncoding RNAs that form closed-loop structures and are reported to alter gene expression through miRNA sequestration, thus providing a potentially novel neurobiological mechanism for the development of alcohol dependence (AD). Genome-wide expression of circRNA was assessed in the nucleus accumbens (NAc) from 32 AD-matched cases/controls. Significant circRNAs (unadj. p ≤ 0.05) were identified via regression and clustered in circRNA networks via weighted gene co-expression network analysis (WGCNA). CircRNA interactions with previously generated mRNA and miRNA were detected via correlation and bioinformatic analyses. Significant circRNAs (N = 542) clustered in nine significant AD modules (FWER p ≤ 0.05), within which we identified 137 circRNA hubs. We detected 23 significant circRNA-miRNA-mRNA interactions (FDR ≤ 0.10). Among these, circRNA-406742 and miR-1200 significantly interact with the highest number of mRNA, including genes associated with neuronal functioning and alcohol addiction (HRAS, PRKCB, HOMER1, and PCLO). Finally, we integrate genotypic information that revealed 96 significant circRNA expression quantitative trait loci (eQTLs) (unadj. p ≤ 0.002) that showed significant enrichment within recent alcohol use disorder (AUD) and smoking genome-wide association study (GWAS). To our knowledge, this is the first study to examine the role of circRNA in the neuropathology of AD. We show that circRNAs impact mRNA expression by interacting with miRNA in the NAc of AD subjects. More importantly, we provide indirect evidence for the clinical importance of circRNA in the development of AUD by detecting a significant enrichment of our circRNA eQTLs among GWAS of substance abuse.

Genetics of substance use disorders in the era of big data

Substance use disorders (SUDs) are conditions in which the use of legal or illegal substances, such as nicotine, alcohol or opioids, results in clinical and functional impairment. SUDs and, more generally, substance use are genetically complex traits that are enormously costly on an individual and societal basis. The past few years have seen remarkable progress in our understanding of the genetics, and therefore the biology, of substance use and abuse. Various studies - including of well-defined phenotypes in deeply phenotyped samples, as well as broadly defined phenotypes in meta-analysis and biobank samples - have revealed multiple risk loci for these common traits. A key emerging insight from this work establishes a biological and genetic distinction between quantity and/or frequency measures of substance use (which may involve low levels of use without dependence), versus symptoms related to physical dependence.

Targeting Epigenetic Mechanisms to Treat Alcohol Use Disorders (AUD)

BACKGROUND

The impact of abusive alcohol consumption on human health is remarkable. According to the World Health Organization (WHO), approximately 3.3 million people die annually because of harmful alcohol consumption (the figure represents around 5.9% of global deaths). Alcohol Use Disorder (AUD) is a chronic disease where individuals exhibit compulsive alcohol drinking and present negative emotional states when they do not drink. In the most severe manifestations of AUD, the individuals lose control over intake despite a decided will to stop drinking. Given the multiple faces and the specific forms of this disease, the term AUD often appears in the plural (AUDs). Since only a few approved pharmacological treatments are available to treat AUD and they do not apply to all individuals or AUD forms, the search for compounds that may help to eliminate the burden of the disease and complement other therapeutical approaches is necessary.

METHODS

This work reviews recent research focused on the involvement of epigenetic mechanisms in the pathophysiology of AUD. Excessive drinking leads to chronic and compulsive consumption that eventually damages the organism. The central nervous system is a key target and is the focus of this study. The search for the genetic and epigenetic mechanisms behind the intricated dysregulation induced by ethanol will aid researchers in establishing new therapy approaches.

CONCLUSION

Recent findings in the field of epigenetics are essential and offer new windows for observation and research. The study of small molecules that inhibit key epienzymes involved in nucleosome architecture dynamics is necessary in order to prove their action and specificity in the laboratory and to test their effectivity and safety in clinical trials with selected patients bearing defined alterations caused by ethanol.

Exploration of alcohol use disorder-associated brain miRNA-mRNA regulatory networks

Transcriptomic changes in specific brain regions can influence the risk of alcohol use disorder (AUD), but the underlying mechanism is not fully understood. We investigated AUD-associated miRNA-mRNA regulatory networks in multiple brain regions by analyzing transcriptomic changes in two sets of postmortem brain tissue samples and ethanol-exposed human embryonic stem cell (hESC)-derived cortical interneurons. miRNA and mRNA transcriptomes were profiled in 192 tissue samples (Set 1) from eight brain regions (amygdala, caudate nucleus, cerebellum, hippocampus, nucleus accumbens, prefrontal cortex, putamen, and ventral tegmental area) of 12 AUD and 12 control European Australians. Nineteen differentially expressed miRNAs (fold-change>2.0 & P < 0.05) and 97 differentially expressed mRNAs (fold-change>2.0 & P < 0.001) were identified in one or multiple brain regions of AUD subjects. AUD-associated miRNA-mRNA regulatory networks in each brain region were constructed using differentially expressed and negatively correlated miRNA-mRNA pairs. AUD-relevant pathways (including CREB Signaling, IL-8 Signaling, and Axonal Guidance Signaling) were potentially regulated by AUD-associated brain miRNA-mRNA pairs. Moreover, miRNA and mRNA transcriptomes were mapped in additional 96 tissue samples (Set 2) from six of the above eight brain regions of eight AUD and eight control European Australians. Some of the AUD-associated miRNA-mRNA regulatory networks were confirmed. In addition, miRNA and mRNA transcriptomes were analyzed in hESC-derived cortical interneurons with or without ethanol exposure, and ethanol-influenced miRNA-mRNA regulatory networks were constructed. This study provided evidence that alcohol could induce concerted miRNA and mRNA expression changes in reward-related or alcohol-responsive brain regions. We concluded that altered brain miRNA-mRNA regulatory networks might contribute to AUD development.

PaIRKAT: A pathway integrated regression-based kernel association test with applications to metabolomics and COPD phenotypes

High-throughput data such as metabolomics, genomics, transcriptomics, and proteomics have become familiar data types within the "-omics" family. For this work, we focus on subsets that interact with one another and represent these "pathways" as graphs. Observed pathways often have disjoint components, i.e., nodes or sets of nodes (metabolites, etc.) not connected to any other within the pathway, which notably lessens testing power. In this paper we propose the Pathway Integrated Regression-based Kernel Association Test (PaIRKAT), a new kernel machine regression method for incorporating known pathway information into the semi-parametric kernel regression framework. This work extends previous kernel machine approaches. This paper also contributes an application of a graph kernel regularization method for overcoming disconnected pathways. By incorporating a regularized or "smoothed" graph into a score test, PaIRKAT can provide more powerful tests for associations between biological pathways and phenotypes of interest and will be helpful in identifying novel pathways for targeted clinical research. We evaluate this method through several simulation studies and an application to real metabolomics data from the COPDGene study. Our simulation studies illustrate the robustness of this method to incorrect and incomplete pathway knowledge, and the real data analysis shows meaningful improvements of testing power in pathways. PaIRKAT was developed for application to metabolomic pathway data, but the techniques are easily generalizable to other data sources with a graph-like structure.

Brain-derived neurotrophic factor predominantly regulates the expression of synapse-related genes in the striatum: Insights from in vitro transcriptomics

Aim

The striatum, a main component of the basal ganglia, is a critical part of the motor and reward systems of the brain. It consists of GABAergic and cholinergic neurons and receives projections of dopaminergic, glutamatergic, and serotonergic neurons from other brain regions. Brain‐derived neurotrophic factor (BDNF) plays multiple roles in the central nervous system, and striatal BDNF has been suggested to be involved in psychiatric and neurodegenerative disorders. However, the transcriptomic impact of BDNF on the striatum remains largely unknown. In the present study, we performed transcriptomic profiling of striatal cells stimulated with BDNF to identify enriched gene sets (GSs) and their novel target genes in vitro.

Methods

We carried out RNA sequencing (RNA‐Seq) of messenger RNA extracted from primary dissociated cultures of rat striatum stimulated with BDNF and conducted Generally Applicable Gene‐set Enrichment (GAGE) analysis on 10599 genes. Significant differentially expressed genes (DEGs) were determined by differential expression analysis for sequence count data 2 (DESeq2).

Results

GAGE analysis identified significantly enriched GSs that included GSs related to regulation and dysregulation of synaptic functions, such as synaptic vesicle cycle and addiction to nicotine and morphine, respectively. It also detected GSs related to various types of synapses, including not only GABAergic and cholinergic synapses but also dopaminergic and glutamatergic synapses. DESeq2 revealed 72 significant DEGs, among which the highest significance was observed in the apolipoprotein L domain containing 1 (Apold1).

Conclusions

The present study indicates that BDNF predominantly regulates the expression of synaptic‐function‐related genes and that BDNF promotes synaptogenesis in various subtypes of neurons in the developing striatum. Apold1 may represent a unique target gene of BDNF in the striatum.

In the present study, we performed transcriptomic profiling of striatal cells stimulated with BDNF to identify enriched gene sets (GSs) in vitro. Generally Applicable Gene‐set Enrichment (GAGE) analysis followed by differential expression analysis for sequence count data 2 (DESeq2) suggested that BDNF predominantly regulates the expression of synaptic‐function‐related genes and that BDNF promotes synaptogenesis in various subtypes of neurons in the developing striatum.

An integrated analysis of mRNAs, lncRNAs, and miRNAs based on weighted gene co-expression network analysis involved in bovine endometritis

In dairy cattle, endometritis is a severe infectious disease that occurs following parturition. It is clear that genetic factors are involved in the etiology of endometritis, however, the molecular pathogenesis of endometritis is not entirely understood. In this study, a system biology approach was used to better understand the molecular mechanisms underlying the development of endometritis. Forty transcriptomic datasets comprising of 20 RNA-Seq (GSE66825) and 20 miRNA-Seq (GSE66826) were obtained from the GEO database. Next, the co-expressed modules were constructed based on RNA-Seq (Rb-modules) and miRNA-Seq (mb-modules) data, separately, using a weighted gene co-expression network analysis (WGCNA) approach. Preservation analysis was used to find the non-preserved Rb-modules in endometritis samples. Afterward, the non-preserved Rb-modules were assigned to the mb-modules to construct the integrated regulatory networks. Just highly connected genes (hubs) in the networks were considered and functional enrichment analysis was used to identify the biological pathways associated with the development of the disease. Furthermore, additional bioinformatic analysis including protein-protein interactions network and miRNA target prediction were applied to enhance the reliability of the results. Thirty-five Rb-modules and 10 mb-modules were identified and 19 and 10 modules were non-preserved, respectively, which were enriched in biological pathways related to endometritis like inflammation and ciliogenesis. Two non-preserved Rb-modules were significantly assigned to three mb-modules and three and two important sub-networks in the Rb-modules were identified, respectively, including important mRNAs, lncRNAs and miRNAs genes like IRAK1, CASP3, CCDC40, CCDC39, ZMYND10, FOXJ1, TLR4, IL10, STAT3, FN1, AKT1, CD68, ENSBTAG00000049936, ENSBTAG00000050527, ENSBTAG00000051242, ENSBTAG00000049287, bta-miR-449, bta-miR-484, bta-miR-149, bta-miR-30b and bta-miR-423. The potential roles of these genes have been previously demonstrated in endometritis or related pathways, which reinforced putative functions of the suggested integrated regulatory networks in the endometritis pathogenesis. These findings may help further elucidate the underlying mechanisms of bovine endometritis.

Effect of chronic ethanol consumption in rhesus macaques on the nucleus accumbens core transcriptome

The nucleus accumbens core (NAcc) has been repeatedly demonstrated to be a key component of the circuitry associated with excessive ethanol consumption. Previous studies have illustrated that in a nonhuman primate (NHP) model of chronic ethanol consumption, there is significant epigenetic remodeling of the NAcc. In the current study, RNA-Seq was used to examine genome-wide gene expression in eight each of control, low/binge (LD*), and high/very high (HD*) rhesus macaque drinkers. Using an FDR < 0.05, zero genes were significantly differentially expressed (DE) between LD* and controls, six genes between HD* and LD*, and 734 genes between HD* and controls. Focusing on HD* versus control DE genes, the upregulated genes (N = 366) were enriched in genes with annotations associated with signal recognition particle (SRP)-dependent co-translational protein targeting to membrane (FDR < 3 × 10-59 ), structural constituent of ribosome (FDR < 3 × 10-47 ), and ribosomal subunit (FDR < 5 × 10-48 ). Downregulated genes (N = 363) were enriched in annotations associated with behavior (FDR < 2 × 10-4 ), membrane organization (FDR < 1 × 10-4 ), inorganic cation transmembrane transporter activity (FDR < 2 × 10-3 ), synapse part (FDR < 4 × 10-10 ), glutamatergic synapse (FDR < 1 × 10-6 ), and GABAergic synapse (FDR < 6 × 10-4 ). Ingenuity Pathway Analysis (IPA) revealed that EIF2 signaling and mTOR pathways were significantly upregulated in HD* animals (FDR < 3 × 10-33 and <2 × 10-16 , respectively). Overall, the data supported our working hypothesis; excessive consumption would be associated with transcriptional differences in GABA/glutamate-related genes.

Unveiling the Pathogenesis of Psychiatric Disorders Using Network Models

Psychiatric disorders are complex brain disorders with a high degree of genetic heterogeneity, affecting millions of people worldwide. Despite advances in psychiatric genetics, the underlying pathogenic mechanisms of psychiatric disorders are still largely elusive, which impedes the development of novel rational therapies. There has been accumulating evidence suggesting that the genetics of complex disorders can be viewed through an omnigenic lens, which involves contextualizing genes in highly interconnected networks. Thus, applying network-based multi-omics integration methods could cast new light on the pathophysiology of psychiatric disorders. In this review, we first provide an overview of the recent advances in psychiatric genetics and highlight gaps in translating molecular associations into mechanistic insights. We then present an overview of network methodologies and review previous applications of network methods in the study of psychiatric disorders. Lastly, we describe the potential of such methodologies within a multi-tissue, multi-omics approach, and summarize the future directions in adopting diverse network approaches.

Role of microRNAs in the pathophysiology of addiction

Addiction is a chronic and relapsing brain disorder characterized by compulsive seeking despite adverse consequences. There are both heritable and epigenetic mechanisms underlying drug addiction. Emerging evidence suggests that non-coding RNAs (ncRNAs) such as microRNAs (miRNAs), long non-coding RNAs, and circular RNAs regulate synaptic plasticity and related behaviors caused by substances of abuse. These ncRNAs modify gene expression and may contribute to the behavioral phenotypes of addiction. Among the ncRNAs, the most widely researched and impactful are miRNAs. The goal in this systematic review is to provide a detailed account of recent research involving the role of miRNAs in addiction. This article is categorized under: RNA Interactions with Proteins and Other Molecules > Small Molecule-RNA Interactions RNA in Disease and Development > RNA in Disease.

Gene Expression Differences Between Young Adults Based on Trauma History and Post-traumatic Stress Disorder

Background: The purpose of this study was to identify gene expression differences associated with post-traumatic stress disorder (PTSD) and trauma exposure (TE) in a three-group study design comprised of those with and without trauma exposure and PTSD. Methods: We conducted gene expression and gene network analyses in a sample (n = 45) composed of female subjects of European Ancestry (EA) with PTSD, TE without PTSD, and controls. Results: We identified 283 genes differentially expressed between PTSD-TE groups. In an independent sample of Veterans (n = 78) a small minority of these genes were also differentially expressed. We identified 7 gene network modules significantly associated with PTSD and TE (Bonferroni corrected p ≤ 0.05), which at a false discovery rate (FDR) of q ≤ 0.2, were significantly enriched for biological pathways involved in focal adhesion, neuroactive ligand receptor interaction, and immune related processes among others. Conclusions: This study uses gene network analyses to identify significant gene modules associated with PTSD, TE, and controls. On an individual gene level, we identified a large number of differentially expressed genes between PTSD-TE groups, a minority of which were also differentially expressed in the independent sample. We also demonstrate a lack of network module preservation between PTSD and TE, suggesting that the molecular signature of PTSD and trauma are likely independent of each other. Our results provide a basis for the identification of likely disease pathways and biomarkers involved in the etiology of PTSD.

Chronic Alcohol Use Induces Molecular Genetic Changes in the Dorsomedial Thalamus of People with Alcohol-Related Disorders

The Mediodorsal (MD) thalamus that represents a fundamental subcortical relay has been underrepresented in the studies focusing on the molecular changes in the brains of subjects with alcohol use disorder (AUD). In the current study, MD thalamic regions from AUD subjects and controls were analyzed with Affymetrix Clariom S human microarray. Long-term alcohol use induced a significant (FDR ≤ 0.05) upregulation of 2802 transcripts and downregulation of 1893 genes in the MD thalamus of AUD subjects. A significant upregulation of GRIN1 (glutamate receptor NMDA type 1) and FTO (alpha-ketoglutarate dependent dioxygenase) was confirmed in western blot analysis. Immunohistochemical staining revealed similar heterogenous distribution of GRIN1 in the thalamic nuclei of both AUD and control subjects. The most prevalent functional categories of upregulated genes were related to glutamatergic and GABAergic neurotransmission, cellular metabolism, and neurodevelopment. The prevalent gene cluster among down-regulated genes was immune system mediators. Forty-two differentially expressed genes, including FTO, ADH1B, DRD2, CADM2, TCF4, GCKR, DPP6, MAPT and CHRH1, have been shown to have strong associations (FDR p < 10-8) with AUD or/and alcohol use phenotypes in recent GWA studies. Despite a small number of subjects, we were able to detect robust molecular changes in the mediodorsal thalamus caused by alcohol emphasizing the importance of deeper brain structures such as diencephalon, in the development of AUD-related dysregulation of neurocircuitry.

Exploration of the Role of Serine Proteinase Inhibitor A3 in Alcohol Dependence Using Gene Expression Omnibus Database

Background: Alcohol dependence is an overall health-related challenge; however, the specific mechanisms underlying alcohol dependence remain unclear. Serine proteinase inhibitor A3 (SERPINA3) plays crucial roles in multiple human diseases; however, its role in alcohol dependence clinical practice has not been confirmed. Methods: We screened Gene Expression Omnibus (GEO) expression profiles, and identified differentially expressed genes (DEGs). Protein-protein interaction (PPI) networks were generated using STRING and Cytoscape, and the key clustering module was identified using the MCODE plugin. SERPINA3-based target microRNA prediction was performed using online databases. Functional enrichment analysis was performed. Fifty-eight patients with alcohol dependence and 20 healthy controls were recruited. Clinical variables were collected and follow-up was conducted for 8 months for relapse. Results: SERPINA3 was identified as a DEG. ELANE and miR-137 were identified after PPI analysis. The enriched functions and pathways included acute inflammatory response, response to stress, immune response, and terpenoid backbone biosynthesis. SERPINA3 concentrations were significantly elevated in the alcohol dependence group than in healthy controls (P < 0.001). According to the median value of SERPINA3 expression level in alcohol dependence group, patients were divided into high SERPINA3 (≥2677.33 pg/ml, n = 29) and low SERPINA3 groups (<2677.33 pg/ml, n = 29). Binary logistic analysis indicated that IL-6 was statistically significant (P = 0.015) Kaplan-Meier survival analysis did not indicate any difference in event-free survival between patients with low and high SERPINA3 levels (P = 0.489) after 8 months of follow-up. Receiver characteristic curve analysis revealed that SERPINA3 had an area under the curve of 0.921 (P < 0.0001), with a sensitivity and specificity of 93.1 and 80.0%, respectively. Cox regression analysis revealed that aspartate transaminase level was a negative predictor of relapse (β = 0.003; hazard ratio = 1.003; P = 0.03). Conclusions: SERPINA3 level was remarkably elevated in patients with alcohol dependence than healthy controls, indicating that SERPINA3 is correlated with alcohol dependence. However, SERPINA3 may not be a potential predictive marker of relapse with patients in alcohol dependence.

Network preservation reveals shared and unique biological processes associated with chronic alcohol abuse in NAc and PFC

Chronic alcohol abuse has been linked to the disruption of executive function and allostatic conditioning of reward response dysregulation in the mesocorticolimbic pathway (MCL). Here, we analyzed genome-wide mRNA and miRNA expression from matched cases with alcohol dependence (AD) and controls (n = 35) via gene network analysis to identify unique and shared biological processes dysregulated in the prefrontal cortex (PFC) and nucleus accumbens (NAc). We further investigated potential mRNA/miRNA interactions at the network and individual gene expression levels to identify the neurobiological mechanisms underlying AD in the brain. By using genotyped and imputed SNP data, we identified expression quantitative trait loci (eQTL) uncovering potential genetic regulatory elements for gene networks associated with AD. At a Bonferroni corrected p≤0.05, we identified significant mRNA (NAc = 6; PFC = 3) and miRNA (NAc = 3; PFC = 2) AD modules. The gene-set enrichment analyses revealed modules preserved between PFC and NAc to be enriched for immune response processes, whereas genes involved in cellular morphogenesis/localization and cilia-based cell projection were enriched in NAc modules only. At a Bonferroni corrected p≤0.05, we identified significant mRNA/miRNA network module correlations (NAc = 6; PFC = 4), which at an individual transcript level implicated miR-449a/b as potential regulators for cellular morphogenesis/localization in NAc. Finally, we identified eQTLs (NAc: mRNA = 37, miRNA = 9; PFC: mRNA = 17, miRNA = 16) which potentially mediate alcohol's effect in a brain region-specific manner. Our study highlights the neurotoxic effects of chronic alcohol abuse as well as brain region specific molecular changes that may impact the development of alcohol addiction.

Assessing the Role of Long Noncoding RNA in Nucleus Accumbens in Subjects With Alcohol Dependence

BACKGROUND

Long noncoding RNA (lncRNA) have been implicated in the etiology of alcohol use. Since lncRNA provide another layer of complexity to the transcriptome, assessing their expression in the brain is the first critical step toward understanding lncRNA functions in alcohol use and addiction. Thus, we sought to profile lncRNA expression in the nucleus accumbens (NAc) in a large postmortem alcohol brain sample.

METHODS

LncRNA and protein-coding gene (PCG) expressions in the NAc from 41 subjects with alcohol dependence (AD) and 41 controls were assessed via a regression model. Weighted gene coexpression network analysis was used to identify lncRNA and PCG networks (i.e., modules) significantly correlated with AD. Within the significant modules, key network genes (i.e., hubs) were also identified. The lncRNA and PCG hubs were correlated via Pearson correlations to elucidate the potential biological functions of lncRNA. The lncRNA and PCG hubs were further integrated with GWAS data to identify expression quantitative trait loci (eQTL).

RESULTS

At Bonferroni adj. p-value ≤ 0.05, we identified 19 lncRNA and 5 PCG significant modules, which were enriched for neuronal and immune-related processes. In these modules, we further identified 86 and 315 PCG and lncRNA hubs, respectively. At false discovery rate (FDR) of 10%, the correlation analyses between the lncRNA and PCG hubs revealed 3,125 positive and 1,860 negative correlations. Integration of hubs with genotype data identified 243 eQTLs affecting the expression of 39 and 204 PCG and lncRNA hubs, respectively.

CONCLUSIONS

Our study identified lncRNA and gene networks significantly associated with AD in the NAc, coordinated lncRNA and mRNA coexpression changes, highlighting potentially regulatory functions for the lncRNA, and our genetic (cis-eQTL) analysis provides novel insights into the etiological mechanisms of AD.

Unsupervised discovery of phenotype-specific multi-omics networks

MOTIVATION

Complex diseases often involve a wide spectrum of phenotypic traits. Better understanding of the biological mechanisms relevant to each trait promotes understanding of the etiology of the disease and the potential for targeted and effective treatment plans. There have been many efforts towards omics data integration and network reconstruction, but limited work has examined the incorporation of relevant (quantitative) phenotypic traits.

RESULTS

We propose a novel technique, sparse multiple canonical correlation network analysis (SmCCNet), for integrating multiple omics data types along with a quantitative phenotype of interest, and for constructing multi-omics networks that are specific to the phenotype. As a case study, we focus on miRNA-mRNA networks. Through simulations, we demonstrate that SmCCNet has better overall prediction performance compared to popular gene expression network construction and integration approaches under realistic settings. Applying SmCCNet to studies on chronic obstructive pulmonary disease (COPD) and breast cancer, we found enrichment of known relevant pathways (e.g. the Cadherin pathway for COPD and the interferon-gamma signaling pathway for breast cancer) as well as less known omics features that may be important to the diseases. Although those applications focus on miRNA-mRNA co-expression networks, SmCCNet is applicable to a variety of omics and other data types. It can also be easily generalized to incorporate multiple quantitative phenotype simultaneously. The versatility of SmCCNet suggests great potential of the approach in many areas.

AVAILABILITY AND IMPLEMENTATION

The SmCCNet algorithm is written in R, and is freely available on the web at https://cran.r-project.org/web/packages/SmCCNet/index.html.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Interplay among miR-29 family, mineral metabolism, and gene regulation in Bos indicus muscle

An interplay between gene expression, mineral concentration, and beef quality traits in Bos indicus muscle has been reported previously under a network approach. However, growing evidence suggested that miRNAs not only modulate gene expression but are also involved with mineral homeostasis. To our knowledge, understanding of the miRNA-gene expression-mineral concentration relationship in mammals is still minimal. Therefore, we carried out a miRNA co-expression and multi-level miRNA-mRNA integration analyses to predict the putative drivers (miRNAs and genes) associated with muscle mineral concentration in Nelore steers. In this study, we identified calcium and iron to be the pivotal minerals associated with miRNAs and gene targets. Furthermore, we identified the miR-29 family (miR-29a, -29b, -29c, -29d-3p, and -29e) as the putative key regulators modulating mineral homeostasis. The miR-29 family targets genes involved with AMPK, insulin, mTOR, and thyroid hormone signaling pathways. Finally, we reported an interplay between miRNAs and minerals acting cooperatively to modulate co-expressed genes and signaling pathways both involved with mineral and energy homeostasis in Nelore muscle. Although we provided some evidence to understand this complex relationship, future work should determine the functional implications of minerals for miRNA levels and their feedback regulation system.\\An interplay between gene expression, mineral concentration, and beef quality traits in Bos indicus muscle has been reported previously under a network approach. However, growing evidence suggested that miRNAs not only modulate gene expression but are also involved with mineral homeostasis. To our knowledge, understanding of the miRNA-gene expression-mineral concentration relationship in mammals is still minimal. Therefore, we carried out a miRNA co-expression and multi-level miRNA-mRNA integration analyses to predict the putative drivers (miRNAs and genes) associated with muscle mineral concentration in Nelore steers. In this study, we identified calcium and iron to be the pivotal minerals associated with miRNAs and gene targets. Furthermore, we identified the miR-29 family (miR-29a, -29b, -29c, -29d-3p, and -29e) as the putative key regulators modulating mineral homeostasis. The miR-29 family targets genes involved with AMPK, insulin, mTOR, and thyroid hormone signaling pathways. Finally, we reported an interplay between miRNAs and minerals acting cooperatively to modulate co-expressed genes and signaling pathways both involved with mineral and energy homeostasis in Nelore muscle. Although we provided some evidence to understand this complex relationship, future work should determine the functional implications of minerals for miRNA levels and their feedback regulation system.

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

BACKGROUND

Genetic factors significantly affect alcohol consumption and vulnerability to withdrawal. Furthermore, some genetic models showing predisposition to severe withdrawal are also predisposed to low ethanol (EtOH) consumption and vice versa, even when tested independently in naïve animals.

METHODS

Beginning with a C57BL/6J × DBA/2J F2 intercross founder population, animals were simultaneously selectively bred for both high alcohol consumption and low acute withdrawal (SOT line), or vice versa (NOT line). Using randomly chosen fourth selected generation (S4) mice (N = 18-22/sex/line), RNA-Seq was employed to assess genome-wide gene expression in ventral striatum. The MegaMUGA array was used to detect genome-wide genotypic differences. Differential gene expression and the weighted gene co-expression network analysis were implemented as described elsewhere (Genes Brain Behav 16, 2017, 462).

RESULTS

The new selection of the SOT and NOT lines was similar to that reported previously (Alcohol Clin Exp Res 38, 2014, 2915). One thousand eight hundred and sixteen transcripts were detected as differentially expressed between the lines. For genes more highly expressed in the SOT line, there was enrichment in genes associated with cell adhesion, synapse organization, and postsynaptic membrane. The genes with a cell adhesion annotation included 23 protocadherins, Mpdz and Dlg2. Genes with a postsynaptic membrane annotation included Gabrb3, Gphn, Grid1, Grin2b, Grin2c, and Grm3. The genes more highly expressed in the NOT line were enriched in a network module (red) with annotations associated with mitochondrial function. Several of these genes were module hub nodes, and these included Nedd8, Guk1, Elof1, Ndufa8, and Atp6v1f.

CONCLUSIONS

Marked effects of selection on gene expression were detected. The NOT line was characterized by higher expression of hub nodes associated with mitochondrial function. Genes more highly expressed in the SOT aligned with previous findings, for example, Colville and colleagues (Genes Brain Behav 16, 2017, 462) that both high EtOH preference and consumption are associated with effects on cell adhesion and glutamate synaptic plasticity.

Chronic Voluntary Ethanol Drinking in Cynomolgus Macaques Elicits Gene Expression Changes in Prefrontal Cortical Area 46

BACKGROUND

Genome-wide profiling to examine brain transcriptional features associated with excessive ethanol (EtOH) consumption has been applied to a variety of species including rodents, nonhuman primates (NHPs), and humans. However, these data were obtained from cross-sectional samples which are particularly vulnerable to individual variation when obtained from small outbred populations typical of human and NHP studies. In the current study, a novel within-subject design was used to examine the effects of voluntary EtOH consumption on prefrontal cortex (PFC) gene expression in a NHP model.

METHODS

Two cohorts of cynomolgus macaques (n = 23) underwent a schedule-induced polydipsia procedure to establish EtOH self-administration followed by 6 months of daily open access to EtOH (4% w/v) and water. Individual daily EtOH intakes ranged from an average of 0.7 to 3.7 g/kg/d. Dorsal lateral PFC area 46 (A46) brain biopsies were collected in EtOH-naïve and control monkeys; contralateral A46 biopsies were collected from the same monkeys following the 6 months of fluid consumption. Gene expression changes were assessed using RNA-Seq paired analysis, which allowed for correction of individual baseline differences in gene expression.

RESULTS

A total of 675 genes were significantly down-regulated following EtOH consumption; these were functionally enriched for immune response, cell adhesion, plasma membrane, and extracellular matrix. A total of 567 genes that were up-regulated following EtOH consumption were enriched in microRNA target sites and included target sites associated with Toll-like receptor pathways. The differentially expressed genes were also significantly enriched in transcription factor binding sites.

CONCLUSIONS

The data presented here are the first to use a longitudinal biopsy strategy to examine how chronic EtOH consumption affects gene expression in the primate PFC. Prominent effects were seen in both cell adhesion and neuroimmune pathways; the latter contained both pro- and antiinflammatory genes. The data also indicate that changes in miRNAs and transcription factors may be important epigenetic regulators of EtOH consumption.

A Pathway-Based Genomic Approach to Identify Medications: Application to Alcohol Use Disorder

Chronic, excessive alcohol use alters brain gene expression patterns, which could be important for initiating, maintaining, or progressing the addicted state. It has been proposed that pharmaceuticals with opposing effects on gene expression could treat alcohol use disorder (AUD). Computational strategies comparing gene expression signatures of disease to those of pharmaceuticals show promise for nominating novel treatments. We reasoned that it may be sufficient for a treatment to target the biological pathway rather than lists of individual genes perturbed by AUD. We analyzed published and unpublished transcriptomic data using gene set enrichment of Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways to identify biological pathways disrupted in AUD brain and by compounds in the Library of Network-based Cellular Signatures (LINCS L1000) and Connectivity Map (CMap) databases. Several pathways were consistently disrupted in AUD brain, including an up-regulation of genes within the Complement and Coagulation Cascade, Focal Adhesion, Systemic Lupus Erythematosus, and MAPK signaling, and a down-regulation of genes within the Oxidative Phosphorylation pathway, strengthening evidence for their importance in AUD. Over 200 compounds targeted genes within those pathways in an opposing manner, more than twenty of which have already been shown to affect alcohol consumption, providing confidence in our approach. We created a user-friendly web-interface that researchers can use to identify drugs that target pathways of interest or nominate mechanism of action for drugs. This study demonstrates a unique systems pharmacology approach that can nominate pharmaceuticals that target pathways disrupted in disease states such as AUD and identify compounds that could be repurposed for AUD if sufficient evidence is attained in preclinical studies.

Silencing synaptic MicroRNA-411 reduces voluntary alcohol consumption in mice

Chronic alcohol consumption alters the levels of microRNAs and mRNAs in the brain, but the specific microRNAs and processes that target mRNAs to affect cellular function and behavior are not known. We examined the in vivo manipulation of previously identified alcohol-responsive microRNAs as potential targets to reduce alcohol consumption. Silencing of miR-411 by infusing antagomiR-411 into the prefrontal cortex of female C57BL/6J mice reduced alcohol consumption and preference, without altering total fluid consumption, saccharin consumption, or anxiety-related behaviors. AntagomiR-411 reduced alcohol consumption when given to mice exposed to a chronic alcohol drinking paradigm but did not affect the acquisition of consumption in mice without a history of alcohol exposure, suggesting that antagomiR-411 has a neuroadaptive, alcohol-dependent effect. AntagomiR-411 decreased the levels of miR-411, as well as the association of immunoprecipitated miR-411 with Argonaute2; and, it increased levels of Faah and Ppard mRNAs. Moreover, antagomiR-411 increased the neuronal expression of glutamate receptor AMPA-2 protein, a known alcohol target and a predicted target of miR-411. These results suggest that alcohol and miR-411 function in a homeostatic manner to regulate synaptic mRNA and protein, thus reversing alcohol-related neuroadaptations and reducing chronic alcohol consumption.

Molecular windows into the human brain for psychiatric disorders

Delineating the pathophysiology of psychiatric disorders has been extremely challenging but technological advances in recent decades have facilitated a deeper interrogation of molecular processes in the human brain. Initial candidate gene expression studies of the postmortem brain have evolved into genome wide profiling of the transcriptome and the epigenome, a critical regulator of gene expression. Here, we review the potential and challenges of direct molecular characterization of the postmortem human brain, and provide a brief overview of recent transcriptional and epigenetic studies with respect to neuropsychiatric disorders. Such information can now be leveraged and integrated with the growing number of genome-wide association databases to provide a functional context of trait-associated genetic variants linked to psychiatric illnesses and related phenotypes. While it is clear that the field is still developing and challenges remain to be surmounted, these recent advances nevertheless hold tremendous promise for delineating the neurobiological underpinnings of mental diseases and accelerating the development of novel medication strategies.

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

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

Postmortem brain tissue as an underutilized resource to study the molecular pathology of neuropsychiatric disorders across different ethnic populations

In recent years, large scale meta-analysis of genome-wide association studies (GWAS) have reliably identified genetic polymorphisms associated with neuropsychiatric disorders such as schizophrenia (SCZ), bipolar disorder (BPD) and major depressive disorder (MDD). However, the majority of disease-associated single nucleotide polymorphisms (SNPs) appear within functionally ambiguous non-coding genomic regions. Recently, increased emphasis has been placed on identifying the functional relevance of disease-associated variants via correlating risk polymorphisms with gene expression levels in etiologically relevant tissues. For neuropsychiatric disorders, the etiologically relevant tissue is brain, which requires robust postmortem sample sizes from varying genetic backgrounds. While small sample sizes are of decreasing concern, postmortem brain databases are composed almost exclusively of Caucasian samples, which significantly limits study design and result interpretation. In this review, we highlight the importance of gene expression and expression quantitative loci (eQTL) studies in clinically relevant postmortem tissue while addressing the current limitations of existing postmortem brain databases. Finally, we introduce future collaborations to develop postmortem brain databases for neuropsychiatric disorders from Chinese and Asian subpopulations.

Beyond genome-wide significance: integrative approaches to the interpretation and extension of GWAS findings for alcohol use disorder

Alcohol use disorder (AUD) is a heritable complex behavior. Due to the highly polygenic nature of AUD, identifying genetic variants that comprise this heritable variation has proved to be challenging. With the exception of functional variants in alcohol metabolizing genes (e.g. ADH1B and ALDH2), few other candidate loci have been confidently linked to AUD. Genome-wide association studies (GWAS) of AUD and other alcohol-related phenotypes have either produced few hits with genome-wide significance or have failed to replicate on further study. These issues reinforce the complex nature of the genetic underpinnings for AUD and suggest that both GWAS studies with larger samples and additional analysis approaches that better harness the nominally significant loci in existing GWAS are needed. Here, we review approaches of interest in the post-GWAS era, including in silico functional analyses; functional partitioning of single nucleotide polymorphism heritability; aggregation of signal into genes and gene networks; and validation of identified loci, genes and gene networks in postmortem brain tissue and across species. These integrative approaches hold promise to illuminate our understanding of the biological basis of AUD; however, we recognize that the main challenge continues to be the extremely polygenic nature of AUD, which necessitates large samples to identify multiple loci associated with AUD liability.

Analysis of whole genome-transcriptomic organization in brain to identify genes associated with alcoholism

Alcohol exposure triggers changes in gene expression and biological pathways in human brain. We explored alterations in gene expression in the Pre-Frontal Cortex (PFC) of 65 alcoholics and 73 controls of European descent, and identified 129 genes that showed altered expression (FDR < 0.05) in subjects with alcohol dependence. Differentially expressed genes were enriched for pathways related to interferon signaling and Growth Arrest and DNA Damage-inducible 45 (GADD45) signaling. A coexpression module (thistle2) identified by weighted gene co-expression network analysis (WGCNA) was significantly correlated with alcohol dependence, alcohol consumption, and AUDIT scores. Genes in the thistle2 module were enriched with genes related to calcium signaling pathways and showed significant downregulation of these pathways, as well as enrichment for biological processes related to nicotine response and opioid signaling. A second module (brown4) showed significant upregulation of pathways related to immune signaling. Expression quantitative trait loci (eQTLs) for genes in the brown4 module were also enriched for genetic associations with alcohol dependence and alcohol consumption in large genome-wide studies included in the Psychiatric Genetic Consortium and the UK Biobank's alcohol consumption dataset. By leveraging multi-omics data, this transcriptome analysis has identified genes and biological pathways that could provide insight for identifying therapeutic targets for alcohol dependence.

miRNA-373 promotes urinary bladder cancer cell proliferation, migration and invasion through upregulating epidermal growth factor receptor

MicroRNA (miRNA)-373 has been demonstrated to be involved in several types of cancer, whereas its involvement in urinary bladder cancer and the mechanism of its function remains poorly understood. The present study aimed to investigate the functionality of miRNA-373 in urinary bladder cancer. Tumor tissues and adjacent healthy tissues were collected from patients with urinary bladder cancer (n=55), and blood samples were collected from patients with urinary bladder cancer and healthy controls (n=45). The expression of miRNA-373 in these tissues was detected by reverse transcription quantitative polymerase chain reaction. The diagnostic value of serum miRNA-373 for urinary bladder cancer was investigated by receiver operating characteristic curve analysis and survival curve analysis, respectively. miRNA-373 mimics were transfected into urinary bladder cancer cells, and the effects on cancer cell proliferation, migration and invasion, and on epidermal growth factor receptor (EGFR) expression was assessed by Cell Counting kit-8 assay, Transwell migration and invasion assays, and western blot analysis. It was identified that the miRNA-373 expression level was increased in tumor tissues compared with adjacent healthy tissues. The serum level of miRNA-373 was increased in patients with cancer compared with the healthy controls. Serum miRNA-373 may be used to accurately predict urinary bladder cancer. miRNA-373 overexpression promoted tumor cell proliferation, migration and invasion, and resulted in upregulated EGFR expression in urinary bladder cancer cells. It was concluded that miRNA-373 overexpression may promote urinary bladder cancer cell proliferation, migration and invasion by upregulating EGFR.

An integrative transcriptome analysis indicates regulatory mRNA-miRNA networks for residual feed intake in Nelore cattle

Residual Feed Intake (RFI) is an economically relevant trait in beef cattle. Among the molecular regulatory mechanisms, microRNAs (miRNAs) are an important dimension in post-transcriptional regulation and have been associated with different biological pathways. Here, we performed differential miRNAs expression and weighted gene co-expression network analyses (WGCNA) to better understand the complex interactions between miRNAs and mRNAs expressed in bovine skeletal muscle and liver. MiRNA and mRNA expression data were obtained from Nelore steers that were genetically divergent for RFI (N = 10 [low RFI or feed efficient]; N = 10 [high RFI or feed inefficient]). Differentially expressed and hub miRNAs such as bta-miR-486, bta-miR-7, bta-miR15a, bta-miR-21, bta-miR 29, bta- miR-30b, bta-miR-106b, bta-miR-199a-3p, bta-miR-204, and bta-miR 296 may have a potential role in variation of RFI. Functional enrichment analysis of differentially expressed (DE) miRNA's target genes and miRNA-mRNA correlated modules revealed that insulin, lipid, immune system, oxidative stress and muscle development signaling pathways might potentially be involved in RFI in this population. Our study identified DE miRNAs, miRNA - mRNA regulatory networks and hub miRNAs related to RFI. These findings suggest a possible role of miRNAs in regulation of RFI, providing new insights into the potential molecular mechanisms that control feed efficiency in Nelore cattle.

Deep sequencing and miRNA profiles in alcohol-induced neuroinflammation and the TLR4 response in mice cerebral cortex

Alcohol abuse can induce brain injury and neurodegeneration, and recent evidence shows the participation of immune receptors toll-like in the neuroinflammation and brain damage. We evaluated the role of miRNAs as potential modulators of the neuroinflammation associated with alcohol abuse and the influence of the TLR4 response. Using mice cerebral cortex and next-generation sequencing (NGS), we identified miRNAs that were differentially expressed in the chronic alcohol-treated versus untreated WT or TLR4-KO mice. We observed a differentially expression of miR-183 Cluster (C) (miR-96/-182/-183), miR-200a and miR-200b, which were down-regulated, while mirR-125b was up-regulated in alcohol-treated WT versus (vs.) untreated mice. These miRNAs modulate targets genes related to the voltage-gated sodium channel, neuron hyperexcitability (Nav1.3, Trpv1, Smad3 and PP1-γ), as well as genes associated with innate immune TLR4 signaling response (Il1r1, Mapk14, Sirt1, Lrp6 and Bdnf). Functional enrichment of the miR-183C and miR-200a/b family target genes, revealed neuroinflammatory pathways networks involved in TLR4 signaling and alcohol abuse. The changes in the neuroinflammatory targets genes associated with alcohol abuse were mostly abolished in the TLR4-KO mice. Our results show the relationship between alcohol intake and miRNAs expression and open up new therapeutically targets to prevent deleterious effects of alcohol on the brain.

Genetics of Alcohol Use Disorder: A Role for Induced Pluripotent Stem Cells?

Alcohol use disorder (AUD) affects millions of people and costs nearly 250 billion dollars annually. Few effective FDA-approved treatments exist, and more are needed. AUDs have a strong heritability, but only a few genes have been identified with a large effect size on disease phenotype. Genomewide association studies (GWASs) have identified common variants with low effect sizes, most of which are in noncoding regions of the genome. Animal models frequently fail to recapitulate key molecular features of neuropsychiatric disease due to the polygenic nature of the disease, partial conservation of coding regions, and significant disparity in noncoding regions. By contrast, human induced pluripotent stem cells (hiPSCs) derived from patients provide a powerful platform for evaluating genes identified by GWAS and modeling complex interactions in the human genome. hiPSCs can be differentiated into a wide variety of human cells, including neurons, glia, and hepatic cells, which are compatible with numerous functional assays and genome editing techniques. In this review, we focus on current applications and future directions of patient hiPSC-derived central nervous system cells for modeling AUDs in addition to highlighting successful applications of hiPSCs in polygenic neuropsychiatric diseases.

An impaired hepatic clock system effects lipid metabolism in rats with nephropathy

Hyperlipidemia is a key clinical feature in patients with nephrotic syndrome (NS) that is associated with the incidence of cardiovascular events. Recent studies have suggested that the disorders of triglycerides, gluconeogenesis and liver glucose metabolism are associated with the abnormal transcription of clock genes. However, changes to the circadian rhythm of blood lipids in NS require further exploration, and the effects of NS on the hepatic clock system remain to be elucidated. In the present study, the impaired diurnal rhythm of the hepatic core clock genes (BMAL1, CLOCK, CRY1, CRY2, PER1 and PER2) significantly induced circadian rhythm abnormalities in liver-specific clock-controlled genes (LXR, CYP7A1, SREBP-1, ABCA1, DEC1 and DEC2; all P<0.05), which were significantly associated with the abnormal diurnal rhythms of triglyceride, total cholesterol, aspartate aminotransferase and alanine aminotransferase (all P<0.05) in rats with Adriamycin-induced nephropathy. Furthermore, a protein-protein interaction network was identified. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses based on the human database was conducted to obtain signaling pathway and correlation prediction analyses of overall human clock and clock-controlled gene correlations. Strong correlations of the aforementioned clock genes were detected (avg. local clustering coefficient, 0.849) which suggested significant enrichment in circadian rhythm signaling. The present results indicated that damage to hepatic clock systems may impact blood lipid circadian rhythm disorders in NS, and offer a starting point for understanding the crosstalk between peripheral organs and peripheral clock systems.

Cerebral white matter sex dimorphism in alcoholism: a diffusion tensor imaging study

Excessive alcohol consumption is associated with brain aberrations, including abnormalities in frontal and limbic brain regions. In a prior diffusion tensor magnetic resonance imaging (dMRI) study of neuronal circuitry connecting the frontal lobes and limbic system structures, we demonstrated decreases in white matter fractional anisotropy in abstinent alcoholic men. In the present study, we examined sex differences in alcoholism-related abnormalities of white matter connectivity and their association with alcoholism history. The dMRI scans were acquired from 49 abstinent alcoholic individuals (26 women) and 41 nonalcoholic controls (22 women). Tract-based spatial statistical tools were used to estimate regional FA of white matter tracts and to determine sex differences and their relation to measures of alcoholism history. Sex-related differences in white matter connectivity were observed in association with alcoholism: Compared to nonalcoholic men, alcoholic men had diminished FA in portions of the corpus callosum, the superior longitudinal fasciculi II and III, and the arcuate fasciculus and extreme capsule. In contrast, alcoholic women had higher FA in these regions. Sex differences also were observed for correlations between corpus callosum FA and length of sobriety. Our results suggest that sexual dimorphism in white matter microstructure in abstinent alcoholics may implicate underlying differences in the neurobehavioral liabilities for developing alcohol abuse disorders, or for sequelae following abuse.

miRNA-regulated transcription associated with mouse strains predisposed to hypnotic effects of ethanol

INTRODUCTION

Studying innate sensitivity to ethanol can be an important first step toward understanding alcohol use disorders. In brain, we investigated transcripts, with evidence of miRNA modulation related to a predisposition to the hypnotic effect of ethanol, as measured by loss of righting reflex (LORR).

METHODS

Expression of miRNAs (12 samples) and expression of mRNAs (353 samples) in brain were independently analyzed for an association with LORR in mice from the LXS recombinant inbred panel gathered across several small studies. These results were then integrated via a meta-analysis of miRNA-mRNA target pairs identified in miRNA-target interaction databases.

RESULTS

We found 112 significant miRNA-mRNA pairs where a large majority of miRNAs and mRNAs were highly interconnected. Most pairs indicated a pattern of increased levels of miRNAs and reduced levels of mRNAs being associated with more alcohol-sensitive strains. For example, CaMKIIn1 was targeted by multiple miRNAs associated with LORR. CAMK2N1 is an inhibitor of CAMK2, which among other functions, phosphorylates, or binds to GABAA and NMDA receptors.

CONCLUSIONS

Our results suggest a novel role of miRNA-mediated regulation of an inhibitor of CAMK2 and its downstream targets including the GABAA and NMDA receptors, which have been previously implicated to have a role in ethanol-induced sedation and sensitivity.

Transcriptome analysis of alcohol-treated microglia reveals downregulation of beta amyloid phagocytosis

BACKGROUND

Microglial activation contributes to the neuropathology associated with chronic alcohol exposure and withdrawal, including the expression of inflammatory and anti-inflammatory genes. In the current study, we examined the transcriptome of primary rat microglial cells following incubation with alcohol alone, or alcohol together with a robust inflammatory stimulus.

METHODS

Primary microglia were prepared from mixed rat glial cultures. Cells were incubated with 75 mM ethanol alone or with proinflammatory cytokines ("TII": IL1β, IFNγ, and TNFα). Isolated mRNA was used for RNAseq analysis and qPCR. Effects of alcohol on phagocytosis were determined by uptake of oligomeric amyloid beta.

RESULTS

Alcohol induced nitrite production in control cells and increased nitrite production in cells co-treated with TII. RNAseq analysis of microglia exposed for 24 h to alcohol identified 312 differentially expressed mRNAs ("Alc-DEs"), with changes confirmed by qPCR analysis. Gene ontology analysis identified phagosome as one of the highest-ranking KEGG pathways including transcripts regulating phagocytosis. Alcohol also increased several complement-related mRNAs that have roles in phagocytosis, including C1qa, b, and c; C3; and C3aR1. RNAseq analysis identified over 3000 differentially expressed mRNAs in microglia following overnight incubation with TII; and comparison to the group of Alc-DEs revealed 87 mRNAs modulated by alcohol but not by TII, including C1qa, b, and c. Consistent with observed changes in phagocytosis-related mRNAs, the uptake of amyloid beta1-42, by primary microglia, was reduced by alcohol.

CONCLUSIONS

Our results define alterations that occur to microglial gene expression following alcohol exposure and suggest that alcohol effects on phagocytosis could contribute to the development of Alzheimer's disease.

Integrating Genetic and Gene Co-expression Analysis Identifies Gene Networks Involved in Alcohol and Stress Responses

Although the link between stress and alcohol is well recognized, the underlying mechanisms of how they interplay at the molecular level remain unclear. The purpose of this study is to identify molecular networks underlying the effects of alcohol and stress responses, as well as their interaction on anxiety behaviors in the hippocampus of mice using a systems genetics approach. Here, we applied a gene co-expression network approach to transcriptomes of 41 BXD mouse strains under four conditions: stress, alcohol, stress-induced alcohol and control. The co-expression analysis identified 14 modules and characterized four expression patterns across the four conditions. The four expression patterns include up-regulation in no restraint stress and given an ethanol injection (NOE) but restoration in restraint stress followed by an ethanol injection (RSE; pattern 1), down-regulation in NOE but rescue in RSE (pattern 2), up-regulation in both restraint stress followed by a saline injection (RSS) and NOE, and further amplification in RSE (pattern 3), and up-regulation in RSS but reduction in both NOE and RSE (pattern 4). We further identified four functional subnetworks by superimposing protein-protein interactions (PPIs) to the 14 co-expression modules, including γ-aminobutyric acid receptor (GABA) signaling, glutamate signaling, neuropeptide signaling, cAMP-dependent signaling. We further performed module specificity analysis to identify modules that are specific to stress, alcohol, or stress-induced alcohol responses. Finally, we conducted causality analysis to link genetic variation to these identified modules, and anxiety behaviors after stress and alcohol treatments. This study underscores the importance of integrative analysis and offers new insights into the molecular networks underlying stress and alcohol responses.