Multi-omic and multi-species meta-analyses of nicotine consumption

Identification of gene co-expression modules from zebrafish brain data: Applications in psychiatry illustrated through alcohol-related traits

Cumulative evidence suggests that zebrafish is a useful model in psychiatric research. Weighted Gene Co-expression Network Analysis (WGCNA) enables the reduction of genome-wide expression data to modules of highly co-expressed genes, which are hypothesized to interact within molecular networks. In this study, we first applied WGCNA to zebrafish brain expression data across different experimental conditions. Then, we characterized the different co-expression modules by gene-set enrichment analysis and hub gene-phenotype association. Finally, we analyzed association of polygenic risk scores (PRSs) based on genes of some interesting co-expression modules with alcohol dependence in 524 patients and 729 controls from Galicia, using competitive tests. Our approach revealed 34 co-expression modules in the zebrafish brain, with some showing enrichment in human synaptic genes, brain tissues, or brain developmental stages. Moreover, certain co-expression modules were enriched in psychiatry-related GWAS and comprised hub genes associated with psychiatry-related traits in both human GWAS and zebrafish models. Expression patterns of some co-expression modules were associated with the tested experimental conditions, mainly with substance withdrawal and cold stress. Notably, a PRS based on genes from co-expression modules exclusively associated with substance withdrawal in zebrafish showed a stronger association with human alcohol dependence than PRSs based on randomly selected brain-expressed genes. In conclusion, our analysis led to the identification of co-expressed gene modules that may model human brain gene networks involved in psychiatry-related traits. Specifically, we detected a cluster of co-expressed genes whose expression was exclusively associated with substance withdrawal in zebrafish, which significantly contributed to alcohol dependence susceptibility in humans.

Loci on chromosome 20 interact with rs16969968 to influence cigarettes per day in European ancestry individuals

BACKGROUND

The understanding of the molecular genetic contributions to smoking is largely limited to the additive effects of individual single nucleotide polymorphisms (SNPs), but the underlying genetic risk is likely to also include dominance, epistatic, and gene-environment interactions.

METHODS

To begin to address this complexity, we attempted to identify genetic interactions between rs16969968, the most replicated SNP associated with smoking quantity, and all SNPs and genes across the genome.

RESULTS

Using the UK Biobank European subsample, we found one SNP, rs1892967, and two genes, PCNA and TMEM230, that showed a significant genome-wide interaction with rs16969968 for log10 CPD and raw CPD, respectively, in a sample of 116 442 individuals who self-reported currently or previously smoking. We extended these analyses to individuals of South Asian descent and meta-analyzed the combined sample of 117 212 individuals of European and South Asian ancestry. We replicated the gene findings in a meta-analysis of five Finnish samples (N=40 140): FinHealth, FINRISK, Finnish Twin Cohort, GeneRISK, and Health-2000-2011.

CONCLUSIONS

To our knowledge, this represents the first reliable epistatic association between single nucleotide polymorphisms for smoking behaviors and provides a novel direction for possible future functional studies related to this interaction. Furthermore, this work demonstrates the feasibility of these analyses by pooling multiple datasets across various ancestries, which may be applied to other top SNPs for smoking and/or other phenotypes.

Testing associations between human anxiety and genes previously implicated by mouse anxiety models

Anxiety disorders are common and can be debilitating, with effective treatments remaining hampered by an incomplete understanding of the underlying genetic etiology. Improvements have been made in understanding the genetic influences on mouse behavioral models of anxiety, yet it is unclear the extent to which genes identified in these experimental systems contribute to genetic variation in human anxiety phenotypes. Leveraging new and existing large-scale human genome-wide association studies, we tested whether sets of genes previously identified in mouse anxiety-like behavior studies contribute to a range of human anxiety disorders. When tested as individual genes, 13 mouse-identified genes were associated with human anxiety phenotypes, suggesting an overlap of individual genes contributing to both mouse models of anxiety-like behaviors and human anxiety traits. When genes were tested as sets, we did identify 14 significant associations between mouse gene sets and human anxiety, but the majority of gene sets showed no significant association with human anxiety phenotypes. These few significant associations indicate a need to identify and develop more translatable mouse models by identifying sets of genes that "match" between model systems and specific human phenotypes of interest. We suggest that continuing to develop improved behavioral paradigms and finer-scale experimental data, for instance from individual neuronal subtypes or cell-type-specific expression data, is likely to improve our understanding of the genetic etiology and underlying functional changes in anxiety disorders.

Gene expression genetics of the striatum of Diversity Outbred mice

Brain transcriptional variation is a heritable trait that mediates complex behaviors, including addiction. Expression quantitative trait locus (eQTL) mapping reveals genomic regions harboring genetic variants that influence transcript abundance. In this study, we profiled transcript abundance in the striatum of 386 Diversity Outbred (J:DO) mice of both sexes using RNA-Seq. All mice were characterized using a behavioral battery of widely-used exploratory and risk-taking assays prior to transcriptional profiling. We performed eQTL mapping, incorporated the results into a browser-based eQTL viewer, and deposited co-expression network members in GeneWeaver. The eQTL viewer allows researchers to query specific genes to obtain allelic effect plots, analyze SNP associations, assess gene expression correlations, and apply mediation analysis to evaluate whether the regulatory variant is acting through the expression of another gene. GeneWeaver allows multi-species comparison of gene sets using statistical and combinatorial tools. This data resource allows users to find genetic variants that regulate differentially expressed transcripts and place them in the context of other studies of striatal gene expression and function in addiction-related behavior.

Polygenic influences on the behavioral effects of alcohol withdrawal in a mixed-ancestry population from the collaborative study on the genetics of alcoholism (COGA)

Alcohol withdrawal (AW) is a feature of alcohol use disorder that may occur in up to half of individuals with chronic, heavy alcohol consumption whenever alcohol use is abruptly stopped or significantly reduced. To date, few genes have been robustly associated with AW; this may be partly due to most studies defining AW as a binary construct despite the multiple symptoms and their range in severity from mild to severe. The current study examined the effects of genome-wide loci on a factor score for AW in high risk and community family samples in the Collaborative Study for the Genetics of Alcoholism (COGA). In addition, we tested whether differentially expressed genes associated with alcohol withdrawal in model organisms are enriched in human genome-wide association study (GWAS) effects. Analyses employed roughly equal numbers of males and females (mean age 35, standard deviation = 15; total N = 8009) and included individuals from multiple ancestral backgrounds. Genomic data were imputed to the HRC reference panel and underwent strict quality control procedures using Plink2. Analyses controlled for age, sex, and population stratification effects using ancestral principal components. We found support that AW is a polygenic disease (SNP-heritability = 0.08 [95 % CI = 0.01, 0.15; pedigree-based heritability = 0.12 [0.08,0.16]. We identified five single nucleotide variants that met genomewide significance, some of which have previously been associated with alcohol phenotypes. Gene-level analyses suggest a role for COL19A1 in AW; H-MAGMA analyses implicated 12 genes associated with AW. Cross-species enrichment analyses indicated that variation within genes identified in model organism studies explained <1 % of the phenotypic variability in human AW. Notably, the surrounding regulatory regions of model organism genes explained more variance than expected by chance, indicating that these regulatory regions and gene sets may be important for human AW. Lastly, when comparing the overlap in genes identified from the human GWAS and H-MAGMA analyses with the genes identified from the animal studies, there was modest overlap, indicating some convergence between the methods and organisms.

Gene expression genetics of the striatum of Diversity Outbred mice

Brain transcriptional variation is a heritable trait that mediates complex behaviors, including addiction. Expression quantitative trait locus (eQTL) mapping reveals genomic regions harboring genetic variants that influence transcript abundance. In this study, we profiled transcript abundance in the striatum of 386 Diversity Outbred (J:DO) mice of both sexes using RNA-Seq. All mice were characterized using a behavioral battery of widely-used exploratory and risk-taking assays prior to transcriptional profiling. We performed eQTL mapping, incorporated the results into a browser-based eQTL viewer, and deposited co-expression network members in GeneWeaver. The eQTL viewer allows researchers to query specific genes to obtain allelic effect plots, analyze SNP associations, assess gene expression correlations, and apply mediation analysis to evaluate whether the regulatory variant is acting through the expression of another gene. GeneWeaver allows multi-species comparison of gene sets using statistical and combinatorial tools. This data resource allows users to find genetic variants that regulate differentially expressed transcripts and place them in the context of other studies of striatal gene expression and function in addiction-related behavior.

Testing Association of Previously Implicated Gene Sets and Gene-Networks in Nicotine Exposed Mouse Models with Human Smoking Phenotypes

INTRODUCTION

Smoking behaviors are partly heritable, yet the genetic and environmental mechanisms underlying smoking phenotypes are not fully understood. Developmental nicotine exposure (DNE) is a significant risk factor for smoking and leads to gene expression changes in mouse models; however, it is unknown whether the same genes whose expression is impacted by DNE are also those underlying smoking genetic liability. We examined whether genes whose expression in D1-type striatal medium spiny neurons due to DNE in the mouse are also associated with human smoking behaviors.

METHODS

Specifically, we assessed whether human orthologs of mouse-identified genes, either individually or as a set, were genetically associated with five human smoking traits using MAGMA and S-LDSC while implementing a novel expression-based gene-SNP annotation methodology.

RESULTS

We found no strong evidence that these genes sets were more strongly associated with smoking behaviors than the rest of the genome, but ten of these individual genes were significantly associated with three of the five human smoking traits examined (p < 2.5e-6). Three of these genes have not been reported previously and were discovered only when implementing the expression-based annotation.

CONCLUSIONS

These results suggest the genes whose expression is impacted by DNE in mice are largely distinct from those contributing to smoking genetic liability in humans. However, examining a single mouse neuronal cell type may be too fine a resolution for comparison, suggesting that experimental manipulation of nicotine consumption, reward, or withdrawal in mice may better capture genes related to the complex genetics of human tobacco use.

IMPLICATIONS

Genes whose expression is impacted by DNE in mouse D1-type striatal medium spiny neurons were not found to be, as a whole, more strongly associated with human smoking behaviors than the rest of the genome, though ten individual mouse-identified genes were associated with human smoking traits. This suggests little overlap between the genetic mechanisms impacted by DNE and those influencing heritable liability to smoking phenotypes in humans. Further research is warranted to characterize how developmental nicotine exposure paradigms in mice can be translated to understand nicotine use in humans and their heritable effects on smoking.

The Genetically Informed Neurobiology of Addiction (GINA) model

Addictions are heritable and unfold dynamically across the lifespan. One prominent neurobiological theory proposes that substance-induced changes in neural circuitry promote the progression of addiction. Genome-wide association studies have begun to characterize the polygenic architecture undergirding addiction liability and revealed that genetic loci associated with risk can be divided into those associated with a general broad-spectrum liability to addiction and those associated with drug-specific addiction risk. In this Perspective, we integrate these genomic findings with our current understanding of the neurobiology of addiction to propose a new Genetically Informed Neurobiology of Addiction (GINA) model.

Meta-Analysis of Transcriptome-Wide Association Studies across 13 Brain Tissues Identified Novel Clusters of Genes Associated with Nicotine Addiction

Genome-wide association studies (GWAS) have identified and reproduced thousands of diseases associated loci, but many of them are not directly interpretable due to the strong linkage disequilibrium among variants. Transcriptome-wide association studies (TWAS) incorporated expression quantitative trait loci (eQTL) cohorts as a reference panel to detect associations with the phenotype at the gene level and have been gaining popularity in recent years. For nicotine addiction, several important susceptible genetic variants were identified by GWAS, but TWAS that detected genes associated with nicotine addiction and unveiled the underlying molecular mechanism were still lacking. In this study, we used eQTL data from the Genotype-Tissue Expression (GTEx) consortium as a reference panel to conduct tissue-specific TWAS on cigarettes per day (CPD) over thirteen brain tissues in two large cohorts: UK Biobank (UKBB; number of participants (N) = 142,202) and the GWAS & Sequencing Consortium of Alcohol and Nicotine use (GSCAN; N = 143,210), then meta-analyzing the results across tissues while considering the heterogeneity across tissues. We identified three major clusters of genes with different meta-patterns across tissues consistent in both cohorts, including homogenous genes associated with CPD in all brain tissues; partially homogeneous genes associated with CPD in cortex, cerebellum, and hippocampus tissues; and, lastly, the tissue-specific genes associated with CPD in only a few specific brain tissues. Downstream enrichment analyses on each gene cluster identified unique biological pathways associated with CPD and provided important biological insights into the regulatory mechanism of nicotine dependence in the brain.

Genes identified in rodent studies of alcohol intake are enriched for heritability of human substance use

BACKGROUND

Rodent paradigms and human genome-wide association studies (GWAS) on drug use have the potential to provide biological insight into the pathophysiology of addiction.

METHODS

Using GeneWeaver, we created rodent alcohol and nicotine gene-sets derived from 19 gene expression studies on alcohol and nicotine outcomes. We partitioned the SNP-heritability of these gene-sets using four large human GWAS: 1) alcoholic drinks per week, 2) problematic alcohol use, 3) cigarettes per day and 4) smoking cessation. We benchmarked our findings with curated human alcoholism and nicotine addiction gene-sets and performed specificity analyses using other rodent gene-sets (e.g., locomotor behavior) and other human GWAS (e.g., height).

RESULTS

The rodent alcohol gene-set was enriched for heritability of drinks per week, cigarettes per day, and smoking cessation, but not problematic alcohol use. However, the rodent nicotine gene-set was not significantly associated with any of these traits. Both rodent gene-sets showed enrichment for several non-substance use GWAS, and the extent of this relationship tended to increase as a function of trait heritability. In general, larger gene-sets demonstrated more significant enrichment. Finally, when evaluating human traits with similar heritabilities, both rodent gene-sets showed greater enrichment for substance use traits.

CONCLUSION

Our results suggest that rodent gene expression studies can help to identify genes that contribute to heritability of some substance use traits in humans, yet there was less specificity than expected. We outline various limitations, interpretations and considerations for future research.

Smoking: Pathogenetic Mechanisms, Possibilities of Medical Correction (Literature Review)

The smoking epidemic is one of the greatest threats to human health, as it is the leading cause of death worldwide, as well as many preventable diseases. Diseases associated with smoking include cardiovascular disease, atherosclerosis, type 2 diabetes, chronic pancreatitis, chronic obstructive pulmonary disease, cancer, osteoporosis and others. The aim is to analyze information on the pathogenetic mechanisms of tobacco exposure and methods of its drug correction. Among the 4,000 chemicals identified in tobacco smoke, nicotine is the main active ingredient in tobacco products. It is highly toxic and potentially lethal. Nicotine affects many systems, including the neurological, neuromuscular, cardiovascular, respiratory, immunological and gastrointestinal systems. Therefore, nicotine-containing drugs reduce the symptoms of nicotine withdrawal in regular smokers who abstain from smoking. The presence of different types of nicotinic acetylcholine receptors, their regulation and desensitization affect these complex physiological effects. Nicotine replacement drugs contain pure nicotine, in order to reduce the patient’s propensity to use tobacco, physiological and psychomotor withdrawal symptoms. They increase the level of nicotine in the blood, thus facilitating the transition from cigarette smoking to complete abstinence. Smokers can quit without using nicotine-containing drugs, but most those who try do not succeed on the first try. Nicotine replacement drugs only affect physical dependence. Nicotine is the main factor responsible for dependence on cigarette smoking, and its rate of absorption and entry into the brain are considered key factors responsible for the high potential for abuse of this drug. The advantage of nicotine-containing drugs is the patient’s independence from the method of administration, determination of abstinence, duration of treatment, level of supportive counseling, place of treatment (hospital, clinic or support group), as well as fixed, variable or reduced doses. There are a number of nicotine-containing drugs in various forms that effectively promote smoking cessation in adults.

Integration of evidence across human and model organism studies: A meeting report

The National Institute on Drug Abuse and Joint Institute for Biological Sciences at the Oak Ridge National Laboratory hosted a meeting attended by a diverse group of scientists with expertise in substance use disorders (SUDs), computational biology, and FAIR (Findability, Accessibility, Interoperability, and Reusability) data sharing. The meeting's objective was to discuss and evaluate better strategies to integrate genetic, epigenetic, and 'omics data across human and model organisms to achieve deeper mechanistic insight into SUDs. Specific topics were to (a) evaluate the current state of substance use genetics and genomics research and fundamental gaps, (b) identify opportunities and challenges of integration and sharing across species and data types, (c) identify current tools and resources for integration of genetic, epigenetic, and phenotypic data, (d) discuss steps and impediment related to data integration, and (e) outline future steps to support more effective collaboration-particularly between animal model research communities and human genetics and clinical research teams. This review summarizes key facets of this catalytic discussion with a focus on new opportunities and gaps in resources and knowledge on SUDs.

Mining the key genes for ventilator-induced lung injury using co-expression network analysis

Mechanical ventilation is extensively adopted in general anesthesia and respiratory failure management, but it can also induce ventilator-induced lung injury (VILI). Therefore, it is of great urgency to explore the mechanisms involved in the VILI pathogenesis, which might contribute to its future prevention and treatment. Four microarray datasets from the GEO database were selected in our investigation, and were subjected to the Weighted Gene Co-Expression Network Analysis (WGCNA) to identify the VILI-correlated gene modules. The limma package in R software was used to identify the differentially expressed genes (DEGs) between the VILI and control groups. WGCNA was constructed by merging the GSE9314, GSE9368, GSE11434 and GSE11662 datasets. A total of 49 co-expression network modules were determined as associated with VILI. The intersected genes between hub genes screened from DEGs for VILI and those identified using WGCNA were as follows: Tlr2, Hmox1, Serpine1, Mmp9, Il6, Il1b, Ptgs2, Fos and Atf3, which were determined to be key genes for VILI. Those key genes were validated by GSE86229 and quantitative PCR (qPCR) experiment to have significantly statistical difference in their expression between the VILI and control groups. In a nutshell, nine key genes with expression differences in VILI were screened by WGCNA by integrating multiple datasets.