Genome‐wide Association Meta‐analyses Identify Novel Genetic Risk Loci and Polygenic Phenotype Associations for Heroin, Methamphetamine and Alcohol Dependences
Su-Hua Chang,Yan Sun,Fan Wang,Xiang-Wen Chang,Ying-Jian Zhang,Tian-Ye Jia,Hong-Qiang Sun,Wei-Hua Yue,Ping Wu,Lin Lu,Jie Shi
DOI: https://doi.org/10.1002/ctm2.659
IF: 8.554
2022-01-01
Clinical and Translational Medicine
Abstract:This work presented a significant correlation between heroin and methamphetamine dependence (HD and MD) which distinguished with alcohol dependence (AD) at the genome-wide level. Three novel risk loci were identified for HD and MD. The shared polygenic risk with cognition and attention-deficit hyperactivity disorder (ADHD) further profiled the similar genetic characteristics between HD and MD compared to AD. Substance dependencies (SD) are one of the leading public health concerns worldwide. Drug markets are in a constant state of flux. The combined use of different addictive substances, especially illegal drugs, is common among addiction patients. Genetic factors contribute to approximately 40%–70% of the variance in persistent SD.1 Capturing the shared and specific genetic mechanism of different substance dependences is crucial for coping with the changeable types of addiction. In this study, genome-wide association meta-analyses (GWMA) were performed independently for HD, MD and AD based on two data sets (DS) of substance-specific dependence patients (1028 HD, 1750 MD, 537 AD and 2862 shared controls for DS1, 980 HD, 701 MD, 224 AD and 1111 shared controls for DS2) (Supplementary Tables S1 and S2, Supplementary Figures S1 and S2, supplementary methods). Consist with our previous study in DS12, One significant locus at chr12 ANKS1B was identified for both HD (peaked in rs112706556, pmeta = 4.99e-8) and MD (peaked in rs58720542, pmeta = 4.401e-9) (Table 1, Figure 1A and B); two well-known loci in chr4 ADH cluster and chr12 ALDH cluster were verified for AD (Supplementary Figure S3). Genetic correlation analysis using LD score regression3 showed MD and HD were significantly correlated (rg = 0.8618, p = 7.361e-5), but not for HD and AD, MD and AD. The cross-dataset pairwise polygenic risk score (PRS) analysis further validated these relationship patterns (Supplementary Table S3). Based on the high genetic correlation between HD and MD, the combined GWMA for HD&MD versus controls were analysed, in which HD and MD were grouped as cases and compared with controls (Table 1 and Figure 1). In addition to the ANKS1B locus (peaked in rs140254085, pmeta = 6.355e-10), another two novel loci located in chr2 locus (peaked in rs74330628 downstream of NRXN1, pmeta = 3.84e-8) and chr7 locus (peaked in GTF2IRD1 intron rs76965632, pmeta = 1.41e-8) were identified (Supplementary Table S4). The functional annotations for the three significant loci are shown in Table 1. Several regulatory features were located in the loci (Supplementary Figure S4), and eQTL data showed the loci regulated the expression of ANKS1B and NRXN1 in brain tissues (Supplementary Table S5, Figure 2A and Supplementary Figure S5). By comparison, the three significant loci of HD&MD were not associated with AD, while the chr4 ADH and chr12 ALDH loci were not associated with HD and MD or had a different direction in HD and MD with AD (Supplementary Figure S6). These findings are consistent with previous twin study that showed a closer genetic correlation across illicit drug dependences compared to alcohol dependence.4 Next, the association with addiction characteristics of the significant loci for HD&MD was examined (Supplementary Table S6 and Figure 2B). The protective allele (A) of rs74330628 in the NRXN1 locus was associated with a lower frequency of heroin usage (padj = .0167). The protective allele (C) of rs76965632 in the GTF2IRD1 locus was negatively associated with craving of MA (padj = .0256). Using the Human Connectome Project data, the SNPs of ANKS1B was associated with risk of ever used illicit drug (EverDrugs) (Supplementary Table S7), which was associated with the volume of the left and right amygdala (Supplementary Table S8). Mediation analysis showed the ANKS1B locus had an indirect effect on EverDrugs by the mediation of the left amygdala and right amygdala (Figure 2C). ANKS1B protein may be involved in the neural plasticity of the amygdala during drug use by affecting glutamatergic neurotransmission.5 Phenome-wide association analysis based on the GWAS Atlas6 further validated the association of ANKS1B, NRXN1 and GTF2IRD1 with psychiatric traits (Supplementary Figure S7 and Supplementary Tables S9–S11). We next explored the associated genes and pathways for the shared risk of HD and MD. ANKS1B, GRM7, RBFOX1 and CDH13 were shared by HD, MD and HD&MD (Supplementary Table S12). Enrichment analysis showed the HD- and MD-related 22 unique genes were significantly enriched in brain-related tissues, GO term ‘modulation of chemical synaptic transmission’ and drug abuse and neurodevelopmental disorder-related diseases (Supplementary Figure S8). These associated genes and pathways provide new candidates and clues for understanding the genetic mechanism of drug dependences. Additionally, polygenic associations with phenotypes of addiction-related traits, risk behaviour, cognition and psychiatric disorders (Supplementary Table S13) were performed for HD, MD and AD. The three alcohol-related traits were significantly associated with AD, but not with HD and MD (Figure 3A and Supplementary Table S14). Sexual and smoke behaviour risks were positively associated with the MD and HD&MD, but not with AD (Figure 3B and Supplementary Table S15). Cognition performance, education attainment and IQ were negatively associated with HD, MD and HD&MD but not with AD (Figure 3C and Supplementary Table S16). Our findings support that pre-existing cognition disruption could increase the risk of illicit drug dependence.7 For the eight psychiatric disorders, only the PRS of ADHD was positively associated with HD, MD and HD&MD, but not with AD (Figure 3D and Supplementary Table S17). Mendelian randomisation analysis further highlighted the causal effect of ADHD on HD and HD&MD (Supplementary Table S18). Our previous study also demonstrated that ADHD-relevant childhood behaviour was a risk factor for MA-induced psychosis.8 First-degree relatives of ADHD probands have an increased risk for drug dependence.9 This suggests that ADHD and drug dependence have shared genetic factors. Cluster analysis based on the standardised coefficient matrix from the PRS association results showed HD, MD and HD&MD were in one sub-cluster and were different from AD, which further highlighted the closer polygenic correlations between HD and MD compared to AD. In conclusion, at the genome-wide loci, genes and polygenic levels, we identified significant genetic overlap between HD and MD, which distinguished with AD. Notably, the combined HD&MD GWAS identified three common risk loci, located on ANSK1B, NRXN1 and GTF2IRD1 genes. At the polygenic level, HD and MD were significantly associated with cognition deficiency and ADHD, which distinguished them from AD. Our results would help with the fine-mapping of the common and unique genetic mechanisms underlying drug dependences and alcohol dependence and may provide clues for their prevention. This work was supported by the National Natural Science Foundation of China (No. 31871259, 82171488, U1802283, 81821092), the National Key Research and Development Program of China (No. 2021YFF0306500), and the Beijing Municipal Science and Technology Commission (No. Z181100001518005). The authors declare no competing interest. Please note: The publisher is not responsible for the content or functionality of any supporting information supplied by the authors. Any queries (other than missing content) should be directed to the corresponding author for the article.