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.

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.

Characterizing cell type specific transcriptional differences between the living and postmortem human brain.

Single-nucleus RNA sequencing (snRNA-seq) is often used to define gene expression patterns characteristic of brain cell types as well as to identify cell type specific gene expression signatures of neurological and mental illnesses in postmortem human brains. As methods to obtain brain tissue from living individuals emerge, it is essential to characterize gene expression differences associated with tissue originating from either living or postmortem subjects using snRNA-seq, and to assess whether and how such differences may impact snRNA-seq studies of brain tissue. To address this, human prefrontal cortex single nuclei gene expression was generated and compared between 31 samples from living individuals and 21 postmortem samples. The same cell types were consistently identified in living and postmortem nuclei, though for each cell type, a large proportion of genes were differentially expressed between samples from postmortem and living individuals. Notably, estimation of cell type proportions by cell type deconvolution of pseudo-bulk data was found to be more accurate in samples from living individuals. To allow for future integration of living and postmortem brain gene expression, a model was developed that quantifies from gene expression data the probability a human brain tissue sample was obtained postmortem. These probabilities are established as a means to statistically account for the gene expression differences between samples from living and postmortem individuals. Together, the results presented here provide a deep characterization of both differences between snRNA-seq derived from samples from living and postmortem individuals, as well as qualify and account for their effect on common analyses performed on this type of data.

Divergent landscapes of A-to-I editing in postmortem and living human brain.

Adenosine-to-inosine (A-to-I) editing is a prevalent post-transcriptional RNA modification within the brain. Yet, most research has relied on postmortem samples, assuming it is an accurate representation of RNA biology in the living brain. We challenge this assumption by comparing A-to-I editing between postmortem and living prefrontal cortical tissues. Major differences were found, with over 70,000 A-to-I sites showing higher editing levels in postmortem tissues. Increased A-to-I editing in postmortem tissues is linked to higher ADAR1 and ADARB1 expression, is more pronounced in non-neuronal cells, and indicative of postmortem activation of inflammation and hypoxia. Higher A-to-I editing in living tissues marks sites that are evolutionarily preserved, synaptic, developmentally timed, and disrupted in neurological conditions. Common genetic variants were also found to differentially affect A-to-I editing levels in living versus postmortem tissues. Collectively, these discoveries illuminate the nuanced functions and intricate regulatory mechanisms of RNA editing within the human brain.

Divergent landscapes of A-to-I editing in postmortem and living human brain.

Adenosine-to-inosine (A-to-I) editing is a prevalent post-transcriptional RNA modification within the brain. Yet, most research has relied on postmortem samples, assuming it is an accurate representation of RNA biology in the living brain. We challenge this assumption by comparing A-to-I editing between postmortem and living prefrontal cortical tissues. Major differences were found, with over 70,000 A-to-I sites showing higher editing levels in postmortem tissues. Increased A-to-I editing in postmortem tissues is linked to higher ADAR and ADARB1 expression, is more pronounced in non-neuronal cells, and indicative of postmortem activation of inflammation and hypoxia. Higher A-to-I editing in living tissues marks sites that are evolutionarily preserved, synaptic, developmentally timed, and disrupted in neurological conditions. Common genetic variants were also found to differentially affect A-to-I editing levels in living versus postmortem tissues. Collectively, these discoveries offer more nuanced and accurate insights into the regulatory mechanisms of RNA editing in the human brain.

Multiomic foundations of human prefrontal cortex tissue function.

The prefrontal cortex (PFC) is a region of the brain that in humans is involved in the production of higher-order functions such as cognition, emotion, perception, and behavior. Neurotransmission in the PFC produces higher-order functions by integrating information from other areas of the brain. At the foundation of neurotransmission, and by extension at the foundation of higher-order brain functions, are an untold number of coordinated molecular processes involving the DNA sequence variants in the genome, RNA transcripts in the transcriptome, and proteins in the proteome. These "multiomic" foundations are poorly understood in humans, perhaps in part because most modern studies that characterize the molecular state of the human PFC use tissue obtained when neurotransmission and higher-order brain functions have ceased (i.e., the postmortem state). Here, analyses are presented on data generated for the Living Brain Project (LBP) to investigate whether PFC tissue from individuals with intact higher-order brain function has characteristic multiomic foundations. Two complementary strategies were employed towards this end. The first strategy was to identify in PFC samples obtained from living study participants a signature of RNA transcript expression associated with neurotransmission measured intracranially at the time of PFC sampling, in some cases while participants performed a task engaging higher-order brain functions. The second strategy was to perform multiomic comparisons between PFC samples obtained from individuals with intact higher-order brain function at the time of sampling (i.e., living study participants) and PFC samples obtained in the postmortem state. RNA transcript expression within multiple PFC cell types was associated with fluctuations of dopaminergic, serotonergic, and/or noradrenergic neurotransmission in the substantia nigra measured while participants played a computer game that engaged higher-order brain functions. A subset of these associations - termed the "transcriptional program associated with neurotransmission" (TPAWN) - were reproduced in analyses of brain RNA transcript expression and intracranial neurotransmission data obtained from a second LBP cohort and from a cohort in an independent study. RNA transcripts involved in TPAWN were found to be (1) enriched for RNA transcripts associated with measures of neurotransmission in rodent and cell models, (2) enriched for RNA transcripts encoded by evolutionarily constrained genes, (3) depleted of RNA transcripts regulated by common DNA sequence variants, and (4) enriched for RNA transcripts implicated in higher-order brain functions by human population genetic studies. In PFC excitatory neurons of living study participants, higher expression of the genes in TPAWN tracked with higher expression of RNA transcripts that in rodent PFC samples are markers of a class of excitatory neurons that connect the PFC to deep brain structures. TPAWN was further reproduced by RNA transcript expression patterns differentiating living PFC samples from postmortem PFC samples, and significant differences between living and postmortem PFC samples were additionally observed with respect to (1) the expression of most primary RNA transcripts, mature RNA transcripts, and proteins, (2) the splicing of most primary RNA transcripts into mature RNA transcripts, (3) the patterns of co-expression between RNA transcripts and proteins, and (4) the effects of some DNA sequence variants on RNA transcript and protein expression. Taken together, this report highlights that studies of brain tissue obtained in a safe and ethical manner from large cohorts of living individuals can help advance understanding of the multiomic foundations of brain function.

A study of gene expression in the living human brain.

A goal of medical research is to determine the molecular basis of human brain health and illness. One way to achieve this goal is through observational studies of gene expression in human brain tissue. Due to the unavailability of brain tissue from living people, most such studies are performed using tissue from postmortem brain donors. An assumption underlying this practice is that gene expression in the postmortem human brain is an accurate representation of gene expression in the living human brain. Here, this assumption - which, until now, had not been adequately tested - is tested by comparing human prefrontal cortex gene expression between 275 living samples and 243 postmortem samples. Expression levels differed significantly for nearly 80% of genes, and a systematic examination of alternative explanations for this observation determined that these differences are not a consequence of cell type composition, RNA quality, postmortem interval, age, medication, morbidity, symptom severity, tissue pathology, sample handling, batch effects, or computational methods utilized. Analyses integrating the data generated for this study with data from earlier landmark studies that used tissue from postmortem brain donors showed that postmortem brain gene expression signatures of neurological and mental illnesses, as well as of normal traits such as aging, may not be accurate representations of these gene expression signatures in the living brain. By using tissue from large cohorts living people, future observational studies of human brain biology have the potential to (1) determine the medical research questions that can be addressed using postmortem tissue as a proxy for living tissue and (2) expand the scope of medical research to include questions about the molecular basis of human brain health and illness that can only be addressed in living people (e.g., "What happens at the molecular level in the brain as a person experiences an emotion?").

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.

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.