KAT2B polymorphism identified for drug abuse in African Americans with regulatory links to drug abuse pathways in human prefrontal cortex

Common and distinguishing genetic factors for substance use behavior and disorder: an integrated analysis of genomic and transcriptomic studies from both human and animal studies

BACKGROUND AND AIMS

Genomic and transcriptomic findings greatly broaden the biological knowledge regarding substance use. However, systematic convergence and comparison evidence of genome-wide findings is lacking for substance use. Here, we combined all the genome-wide findings from both substance use behavior and disorder (SUBD) and identified common and distinguishing genetic factors for different SUBDs.

METHODS

Systemic literature search for genome-wide association (GWAS) and RNA-seq studies of alcohol/nicotine/drug use behavior (partially meets or not reported diagnostic criteria) and alcohol use behavior and disorder (AUBD), nicotine use behavior and disorder (NUBD) and drug use behavior and disorder (DUBD) was performed using PubMed and the GWAS catalog. Drug use was focused upon cannabis, opioid, cocaine and methamphetamine use. GWAS studies required case-control or case/cohort samples. RNA-seq studies were based on brain tissues. The genes which contained significant single nucleotide polymorphism (P ≤ 1 × 10^-6 ) in GWAS and reported as significant in RNA-seq studies were extracted. Pathway enrichment was performed by using Metascape. Gene interaction networks were identified by using the Protein Interaction Network Analysis database.

RESULTS

Total SUBD-related 2910 genes were extracted from 75 GWAS studies (2 773 889 participants) and 17 RNA-seq studies. By overlapping the genes and pathways of AUBD, NUBD and DUBD, four shared genes (CACNB2, GRIN2B, PLXDC2 and PKNOX2), four shared pathways [two Gene Ontology (GO) terms of 'modulation of chemical synaptic transmission', 'regulation of trans-synaptic signaling', two Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways of 'dopaminergic synapse', 'cocaine addiction'] were identified (significantly higher than random, P < 1 × 10^-5 ). The top shared KEGG pathways (Benjamini-Hochberg-corrected P-value < 0.05) in the pairwise comparison of AUBD versus DUBD, NUBD versus DUBD, AUBD versus NUBD were 'Epstein-Barr virus infection', 'protein processing in endoplasmic reticulum' and 'neuroactive ligand-receptor interaction', respectively. We also identified substance-specific genetic factors: i.e. ADH1B and ALDH2 were unique for AUBD, while CHRNA3 and CHRNA4 were unique for NUBD.

CONCLUSIONS

This systematic review identifies the shared and unique genes and pathways for alcohol, nicotine and drug use behaviors and disorders at the genome-wide level and highlights critical biological processes for the common and distinguishing vulnerability of substance use behaviors and disorders.

GAWMerge expands GWAS sample size and diversity by combining array-based genotyping and whole-genome sequencing

Genome-wide association studies (GWAS) have made impactful discoveries for complex diseases, often by amassing very large sample sizes. Yet, GWAS of many diseases remain underpowered, especially for non-European ancestries. One cost-effective approach to increase sample size is to combine existing cohorts, which may have limited sample size or be case-only, with public controls, but this approach is limited by the need for a large overlap in variants across genotyping arrays and the scarcity of non-European controls. We developed and validated a protocol, Genotyping Array-WGS Merge (GAWMerge), for combining genotypes from arrays and whole-genome sequencing, ensuring complete variant overlap, and allowing for diverse samples like Trans-Omics for Precision Medicine to be used. Our protocol involves phasing, imputation, and filtering. We illustrated its ability to control technology driven artifacts and type-I error, as well as recover known disease-associated signals across technologies, independent datasets, and ancestries in smoking-related cohorts. GAWMerge enables genetic studies to leverage existing cohorts to validly increase sample size and enhance discovery for understudied traits and ancestries.

Molecular genetics of cocaine use disorders in humans

Drug addiction, one of the major health problems worldwide, is characterized by the loss of control in drug intake, craving, and withdrawal. At the individual level, drugs of abuse produce serious consequences on health and have a negative impact on the family environment and on interpersonal and work relationships. At a wider scale, they have significant socio-economic and public health consequences and they cause delinquency and citizen insecurity. Cocaine, a psychostimulant substance, is one of the most used illicit drugs, especially in America, Western Europe, and Australia. Cocaine use disorders (CUD) are complex multifactorial conditions driven by both genetic and environmental influences. Importantly, not all people who use cocaine develop CUD, and this is due, at least in part, to biological factors that are encoded in the genome of individuals. Acute and repeated use of cocaine induces epigenetic and gene expression changes responsible for the neuronal adaptations and the remodeling of brain circuits that lead to the transition from use to abuse or dependence. The purpose of this review is to delineate such factors, which should eventually help to understand the inter-individual variability in the susceptibility to cocaine addiction. Heritability estimates for CUD are high and genetic risk factors for cocaine addiction have been investigated by candidate gene association studies (CGAS) and genome-wide association studies (GWAS), reviewed here. Also, the high comorbidity that exists between CUD and several other psychiatric disorders is well known and includes phenotypes like schizophrenia, aggression, antisocial or risk-taking behaviors. Such comorbidities are associated with a worse lifetime trajectory, and here we report shared genetic factors that may contribute to them. Gene expression changes and epigenetic modifications induced by cocaine use and chronic abuse in humans are addressed by reviewing transcriptomic studies performed on neuronal cells and on postmortem brains. We report some genes which expression is altered by cocaine that also bear genetic risk variants for the disorder. Finally, we have a glance to the pharmacogenetics of CUD treatments, still in early stages. A better understanding of the genetic underpinnings of CUD will foster the search of effective treatments and help to move forward to personalized medicine.

Exploring allele specific methylation in drug dependence susceptibility

Drug dependence is a neuropsychiatric condition that involves genetic, epigenetic and environmental factors. Allele-specific methylation (ASM) is a common and stable epigenetic mechanism that involves genetic variants correlating with differential levels of methylation at CpG sites. We selected 182 single-nucleotide polymorphisms (SNPs) described to influence cis ASM in human brain regions to evaluate their possible contribution to drug dependence susceptibility. We performed a case-control association study in a discovery sample of 578 drug-dependent patients (including 428 cocaine-dependent subjects) and 656 controls from Spain, and then, we followed-up the significant associations in an independent sample of 1119 cases (including 589 cocaine-dependent subjects) and 1092 controls. In the discovery sample, we identified five nominal associations, one of them replicated in the follow-up sample (rs6020251). The pooled analysis revealed an association between drug dependence and rs6020251 but also rs11585570, both overcoming the Bonferroni correction for multiple testing. We performed the same analysis considering only cocaine-dependent patients and obtained similar results. The rs6020251 variant correlates with differential methylation levels of cg17974185 and lies in the first intron of the CTNNBL1 gene, in a genomic region with multiple histone marks related to enhancer and promoter regions in brain. Rs11585570 is an eQTL in brain and blood for the SCP2 and ECHDC2 genes and correlates with differential methylation of cg27535305 and cg13461509, located in the promoter regions of both genes. To conclude, using an approach that combines genetic and epigenetic data, we highlighted the CTNNBL1, SCP2 and ECHDC2 genes as potential contributors to drug dependence susceptibility.

Integrating genomics and transcriptomics: Towards deciphering ADHD

Attention-deficit/hyperactivity disorder (ADHD) is a highly heritable condition that represents the most common neurodevelopmental disorder in childhood, persisting into adulthood in around 40-65% of the cases. ADHD is characterised by age-inappropriate symptoms of inattention, impulsivity, and hyperactivity. Mounting evidence points towards ADHD having a strong genetic component and the first genome-wide significant findings have recently been reported. However, the functional characterization of variants unravelled by genome-wide association studies (GWAS) is challenging. Likewise, gene expression profiling studies have also been undertaken and novel integrative approaches combining genomic and transcriptomic data are starting to be conducted, which offers an exciting way that might provide a more informative insight towards the genetic architecture of ADHD. In this review, we summarised current knowledge on genomics, transcriptomics and integrative approaches in ADHD, focusing on GWAS and GWAS meta-analyses (GWAS-MA)- as genomics analyses- microarray and RNA-seq- as transcriptomics analyses-, and studies integrating genomics and transcriptomics data. In addition, current strengths and limitations of such approaches are discussed and further research avenues are proposed in order to face unsolved issues. Although important progress has been made, there is still a long way ahead to elucidate the biological mechanisms of ADHD, which eventually may lead to more personalized approaches in the future. Large- scale research efforts and new technological and statistical approaches are envisaged as important means towards deciphering ADHD in the upcoming years.

Different biases in meta-analyses of case-control and cohort studies: an example from genomics and precision medicine

PURPOSE

Meta-analyses of observational studies reduce the role of chance but also introduce bias because the individual component studies are not randomized. Further, it is plausible that the bias may be different in case-control and cohort studies. We explored these issues in meta-analyses of observational studies of Opioid Use Disorder (OUD).

METHODS

From a systematic literature review of 152 published meta-analyses, 11 fulfilled the initial inclusion criteria of observational studies of OUD. Of these, 9 were meta-analyses of case-control studies and 2 were meta-analyses of cohort studies but only 4 (3 case-control and 1 cohort) targeted more than one specific chromosomal location.

RESULTS

The meta-analyses of the 3 case-control studies, which included 13 individual studies, identified 12 different single nucleotide polymorphisms on 6 different genes on 5 different chromosomes. None was the same as the gene on Chromosome 15 identified from the meta-analysis of the cohort studies.

CONCLUSIONS

These data, from genetic studies, suggest biases are different in meta-analyses of case-control and cohort studies, perhaps due to greater selection bias in case-control studies. These observations have potential importance in the application of meta-analyses to many common and serious diseases, as well as genomics and precision medicine, including OUD.

Association of Stress, Glucocorticoid Receptor, and FK506 Binding Protein Gene Polymorphisms With Internalizing Disorders Among HIV-Infected Children and Adolescents From Kampala and Masaka Districts-Uganda

Children and adolescents living with human immunodeficiency virus (CA-HIV) suffer a considerable burden of internalizing disorders (IDs; depressive and anxiety disorders). Environmental and genetic factors have been reported to influence the vulnerability to IDs in western settings; however, their role among African populations remains inadequately explored. We investigated the individual and interactive effects of stress and single-nucleotide polymorphisms within the FK506 binding protein 5 (rs1360780) and glucocorticoid receptor (rs10482605) genes on ID status in a cohort of CA-HIV in Uganda. We genotyped rs10482605 (309 cases and 315 controls) and rs1360780 (350 cases and 335 controls) among CA-HIV with and without IDs using Kompetitive Allele-Specific PCR. Socio-demographic variables, as well as allele and genotype distributions, were compared between cases and controls using chi-square tests. Genotypes were assessed for Hardy-Weinberg equilibrium. Composite indices of recent and chronic stress classes were also generated. A hierarchical cluster analysis was used to generate cutoff points within each of the indices of recent and chronic stress. Logistic regression was used to assess the association between IDs and each of recent stress, chronic stress, and the investigated genotypes. The interaction effect of chronic/recent stress on the association between each of the polymorphisms and IDs was determined using a likelihood ratio test. We observed no significant association between IDs and rs1360780 and rs10482605 polymorphisms within the FKBP5 and glucocorticoid receptor genes, respectively (P > 0.050). Severe recent stress increased the vulnerability to IDs among CA-HIV (P = 0.001). We did not observe any gene-environment effect on vulnerability to IDs in this population. These findings support the currently held opinion that polymorphisms at single genetic loci only contribute a very small effect to the genetic vulnerability to IDs.

Genome-wide association meta-analysis of cocaine dependence: Shared genetics with comorbid conditions

Cocaine dependence is a complex psychiatric disorder that is highly comorbid with other psychiatric traits. Twin and adoption studies suggest that genetic variants contribute substantially to cocaine dependence susceptibility, which has an estimated heritability of 65-79%. Here we performed a meta-analysis of genome-wide association studies of cocaine dependence using four datasets from the dbGaP repository (2085 cases and 4293 controls, all of them selected by their European ancestry). Although no genome-wide significant hits were found in the SNP-based analysis, the gene-based analysis identified HIST1H2BD as associated with cocaine-dependence (10% FDR). This gene is located in a region on chromosome 6 enriched in histone-related genes, previously associated with schizophrenia (SCZ). Furthermore, we performed LD Score regression analysis with comorbid conditions and found significant genetic correlations between cocaine dependence and SCZ, ADHD, major depressive disorder (MDD) and risk taking. We also found, through polygenic risk score analysis, that all tested phenotypes are significantly associated with cocaine dependence status: SCZ (R² = 2.28%; P = 1.21e-26), ADHD (R² = 1.39%; P = 4.5e-17), risk taking (R² = 0.60%; P = 2.7e-08), MDD (R² = 1.21%; P = 4.35e-15), children's aggressive behavior (R² = 0.3%; P = 8.8e-05) and antisocial behavior (R² = 1.33%; P = 2.2e-16). To our knowledge, this is the largest reported cocaine dependence GWAS meta-analysis in European-ancestry individuals. We identified suggestive associations in regions that may be related to cocaine dependence and found evidence for shared genetic risk factors between cocaine dependence and several comorbid psychiatric traits. However, the sample size is limited and further studies are needed to confirm these results.

Prediction of causal genes and gene expression analysis of attention-deficit hyperactivity disorder in the different brain region, a comprehensive integrative analysis of ADHD

Recent genome-wide association study (GWAS) identified 12 independent loci for Attention-deficit hyperactivity disorder (ADHD). However, the causal genes expression and pathways of ADHD is still vague. We integrated GWAS, eQTL and genes expression data to find the causal genes, genes expression, and genes prioritization in the different brain tissues and whole blood cells. Overall 47 genes were prioritized, the most promising genes were LSG1, HYAL3, PIDD, PNPLA2, BLOC1S2, PLK1S1, CALN1, KAT2B, CTNNB1 and WDR11. Whereas, the CALN1, KAT2B, and WDR11 were previously associated with schizophrenia (SZ), bipolar (BP) and drug abuse. Gene ontology analysis shows that the glutamate receptor signaling pathway (P = 8.009E-07, with false discovery rate (FDR) < 5%), GRIK5 sub network (P = 2.887E-06, FDR < 5%), abnormal gait (P = 3.657E-06, FDR < 5%), REACTOME_SIGNALING_BY_ERBB2 (P = 5.161E-06, FDR < 5%), and abnormal nervous system physiology (P = 5.239E-06, FDR < 5%) were associated with ADHD. These causal genes were highly expressed in Fetal Astrocytes, Neurons, and Microglia/Macrophage. This study illustrates the comprehensive GWAS integrative approach of ADHD. However, further genetic and functional studies are required to validate the role of these genes in the etiology of ADHD, which should provide novel insights into the understanding of this disease.

A case-control genome wide association study of substance use disorder (SUD) identifies novel variants on chromosome 7p14.1 in patients from the United Arab Emirates (UAE)

Genome wide association studies (GWASs) have provided insights into the molecular basis of the disorder in different population. This study presents the first GWAS of substance use disorder (SUD) in patients from the United Arab Emirates (UAE). The aim was to identify genetic association(s) that may provide insights into the molecular basis of the disorder. The GWAS discovery cohort consisted of 512 (250 cases and 262 controls) male participants from the UAE. Controls with no prior history of SUD were available from the Emirates family registry. The replication cohort consisted of 520 (415 cases and 105 controls) Australian male Caucasian participants. The GWAS discovery samples were genotyped for 4.6 million single nucleotide polymorphism (SNP). The replication cohort was genotyped using TaqMan assay. The GWAS association analysis identified three potential SNPs rs118129027 (p-value = 6.24 × 10^-8 ), rs74477937 (p-value = 8.56 × 10^-8 ) and rs78707086 (p-value = 8.55 × 10^-8 ) on ch7p14.1, that did not meet the GWAS significance threshold but were highly suggestive. In the replication cohort, the association of the three top SNPs did not reach statistical significance. In a meta-analysis of the discovery and the replication cohorts, there were no strengthen evidence for association of the three SNPs. The top identified rs118129027 overlaps with a regulatory factor (enhancer) region that targets three neighboring genes LOC105375237, LOC105375240, and YAE1D1. The YAE1D1, which represents a potential locus that is involved in regulating translation initiation pathway. Novel associations that require further confirmation were identified, suggesting a new insight to the genetic basis of SUD.

Evaluation of previous substance dependence genome-wide significant findings in a Spanish sample

BACKGROUND

Substance dependence is a chronic and relapsing disorder explained by genetic and environmental risk factors. The aim of our study is to replicate previous genome-wide significant (GWS) hits identified in substance dependence in general or in cocaine dependence in particular using an independent sample from Spain.

METHODS

We evaluated, in a Spanish sample of 1711 subjects with substance dependence (1011 of them cocaine dependent) and 1719 control individuals, three SNPs identified as GWS in previous studies: rs1868152 and rs2952621 (located near LINC02052 and LINC01854, respectively), associated with substance dependence, and rs2629540 (in the first intron of FAM53B), associated with cocaine dependence.

RESULTS

We replicated the association between rs2952621 and substance dependence under the dominant model (P = 0.020), with the risk allele (T) being the same in our sample and in those two reported previously. We then performed a meta-analysis of the two samples used in the original study that reported the association of rs2952621 with substance dependence (Collaborative Studies on Genetics of Alcoholism (COGA) and Study of Addiction: Genetics and Environment (SAGE)) together with our Spanish sample. The meta-analysis of 3747 cases and 4043 controls confirmed the association (OR = 1.26, 95% CI = 1.15-1.39).

CONCLUSIONS

The rs2952621 variant, located downstream from the yet uncharacterized gene LINC01854, is associated with substance dependence in our Spanish sample. Further research is needed to understand its contribution to the susceptibility to substance dependence.

Human Genetics of Addiction: New Insights and Future Directions

PURPOSE OF REVIEW

With the advent of the genome-wide association study (GWAS), our understanding of the genetics of addiction has made significant strides forward. Here, we summarize genetic loci containing variants identified at genome-wide statistical significance (P < 5 × 10-8) and independently replicated, review evidence of functional or regulatory effects for GWAS-identified variants, and outline multi-omics approaches to enhance discovery and characterize addiction loci.

RECENT FINDINGS

Replicable GWAS findings span 11 genetic loci for smoking, eight loci for alcohol, and two loci for illicit drugs combined and include missense functional variants and noncoding variants with regulatory effects in human brain tissues traditionally viewed as addiction-relevant (e.g., prefrontal cortex [PFC]) and, more recently, tissues often overlooked (e.g., cerebellum). GWAS analyses have discovered several novel, replicable variants contributing to addiction. Using larger sample sizes from harmonized datasets and new approaches to integrate GWAS with multiple 'omics data across human brain tissues holds great promise to significantly advance our understanding of the biology underlying addiction.

Relationships between cell cycle pathway gene polymorphisms and risk of hepatocellular carcinoma

AIM: To investigate the associiations between the polymorphisms of cell cycle pathway genes and the risk of hepatocellular carcinoma (HCC).

METHODS: We enrolled 1127 cases newly diagnosed with HCC from the Tumor Hospital of Guangxi Medical University and 1200 non-tumor patients from the First Affiliated Hospital of Guangxi Medical University. General demographic characteristics, behavioral information, and hematological indices were collected by unified questionnaires. Genomic DNA was isolated from peripheral venous blood using Phenol-Chloroform. The genotyping was performed using the Sequenom MassARRAY iPLEX genotyping method. The association between genetic polymorphisms and risk of HCC was shown by P-value and the odd ratio (OR) with 95% confidence interval (CI) using the unconditional logistic regression after adjusting for age, sex, nationality, smoking, drinking, family history of HCC, and hepatitis B virus (HBV) infection. Moreover, stratified analysis was conducted on the basis of the status of HBV infection, smoking, and alcohol drinking.

RESULTS: The HCC risk was lower in patients with the MCM4 rs2305952 CC (OR = 0.22, 95%CI: 0.08-0.63, P = 0.01) and with the CHEK1 rs515255 TC, TT, TC/TT (OR = 0.73, 95%CI: 0.56-0.96, P = 0.02; OR = 0.67, 95%CI: 0.46-0.97, P = 0.04; OR = 0.72, 95%CI: 0.56-0.92, P = 0.01, respectively). Conversely, the HCC risk was higher in patients with the KAT2B rs17006625 GG (OR = 1.64, 95%CI: 1.01-2.64, P = 0.04). In addition, the risk was markedly lower for those who were carriers of MCM4 rs2305952 CC and were also HBsAg-positive and non-drinking and non-smoking (P < 0.05, respectively) and for those who were carriers of CHEK1 rs515255 TC, TT, TC/TT and were also HBsAg-negative and non-drinking (P < 0.05, respectively). Moreover, the risk was higher for those who were carriers of KAT2B rs17006625 GG and were also HBsAg-negative (P < 0.05).

CONCLUSION: Of 12 cell cycle pathway genes, MCM4, CHEK1 and KAT2B polymorphisms may be associated with the risk of HCC.

Genetic background and epigenetic modifications in the core of the nucleus accumbens predict addiction-like behavior in a rat model

This study provides a demonstration in the rat of a clear genetic difference in the propensity for addiction-related behaviors following prolonged cocaine self-administration. It relies on the use of selectively bred high-responder (bHR) and low-responder (bLR) rat lines that differ in several characteristics associated with "temperament," including novelty-induced locomotion and impulsivity. We show that bHR rats exhibit behaviors reminiscent of human addiction, including persistent cocaine-seeking and increased reinstatement of cocaine seeking. To uncover potential underlying mechanisms of this differential vulnerability, we focused on the core of the nucleus accumbens and examined expression and epigenetic regulation of two transcripts previously implicated in bHR/bLR differences: fibroblast growth factor (FGF2) and the dopamine D2 receptor (D2). Relative to bHRs, bLRs had lower FGF2 mRNA levels and increased association of a repressive mark on histones (H3K9me3) at the FGF2 promoter. These differences were apparent under basal conditions and persisted even following prolonged cocaine self-administration. In contrast, bHRs had lower D2 mRNA under basal conditions, with greater association of H3K9me3 at the D2 promoter and these differences were no longer apparent following prolonged cocaine self-administration. Correlational analyses indicate that the association of H3K9me3 at D2 may be a critical substrate underlying the propensity to relapse. These findings suggest that low D2 mRNA levels in the nucleus accumbens core, likely mediated via epigenetic modifications, may render individuals more susceptible to cocaine addiction. In contrast, low FGF2 levels, which appear immutable even following prolonged cocaine exposure, may serve as a protective factor.

A Review of Genome-Wide Association Studies of Stimulant and Opioid Use Disorders

Substance use disorders (SUD) are a major contributor to disability and disease burden worldwide. Risk for developing SUDs is influenced by variation in the genome. Identifying the genetic variants that influence SUD risk may help us to understand the biological mechanisms for the disorders and improve treatments. Genome-wide association studies (GWAS) have been successful in identifying many regions of the genome associated with common human disorders. Here, findings from recent GWAS of SUDs that involve illicit substances will be reviewed. Several GWAS have been reported, including studies on opioid and stimulant use disorder (cocaine and methamphetamine). Several of these GWAS report associations that are biologically interesting and statistically robust. Replication of the associations in independent samples and functional studies to understand the basis for the statistical associations will be important next steps.