Genetic Underpinnings of the Transition From Alcohol Consumption to Alcohol Use Disorder: Shared and Unique Genetic Architectures in a Cross-Ancestry Sample
Publication: American Journal of Psychiatry
Abstract
Objective:
Recent genome-wide association studies (GWASs) of alcohol-related phenotypes have uncovered key differences in the underlying genetic architectures of alcohol consumption and alcohol use disorder (AUD), with the two traits having opposite genetic correlations with psychiatric disorders. Understanding the genetic factors that underlie the transition from heavy drinking to AUD has important theoretical and clinical implications.
Methods:
The authors used longitudinal data from the cross-ancestry Million Veteran Program sample to identify 1) novel loci associated with AUD and alcohol consumption (measured by the score on the consumption subscale of the Alcohol Use Disorders Identification Test [AUDIT-C]), 2) the impact of phenotypic variation on genetic discovery, and 3) genetic variants with direct effects on AUD that are not mediated through alcohol consumption.
Results:
The authors identified 26 loci associated with AUD and 22 loci associated with AUDIT-C score, including ancestry-specific and novel loci. In secondary GWASs that excluded individuals who report abstinence, the authors identified seven additional loci for AUD and eight additional loci for AUDIT-C score. Although the heterogeneity of the abstinent group biases the GWAS findings, unique variance between alcohol consumption and disorder remained after the abstinent group was excluded. Finally, using mediation analysis, the authors identified a set of variants with effects on AUD that are not mediated through alcohol consumption.
Conclusions:
Differences in genetic architecture between alcohol consumption and AUD are consistent with their having different biological contributions. Genetic variants with direct effects on AUD are potentially relevant to understanding the transition from heavy alcohol consumption to AUD and may be targets for translational prevention and treatment efforts.
Heavy alcohol use is common in the United States and is associated with multiple adverse consequences (1). Regular heavy drinking is the major epidemiological risk factor for alcohol use disorder (AUD), a problematic pattern of alcohol use accompanied by clinically significant impairment or distress (2). Recent genome-wide association studies (GWASs) have identified more than 30 genome-wide-significant variants contributing to the risk of AUD, alcohol dependence, or problematic alcohol use (3–5), and more than 100 genome-wide-significant variants have been identified as contributing to alcohol consumption measures (4, 6, 7) (e.g., score on the consumption subscale of the Alcohol Use Disorders Identification Test [AUDIT-C] [8], based on the first three questions of the AUDIT [9]). Nearly all GWAS loci have been identified in individuals of European ancestry, because comparable data from other population groups are not available.
Although the genetic variation contributing to alcohol consumption and AUD partially overlap, multiple variants are associated with one, but not the other, trait (10). Estimates of the genetic correlation between alcohol consumption and either AUD or problematic alcohol use are moderate (rg values, 0.52–0.67) (4, 5), suggesting that despite alcohol use being necessary for AUD, each trait has unique genetic liability. Furthermore, genetic correlations and phenome-wide association studies (PheWASs) of genetic liability for AUD and alcohol consumption show differential patterns of association with physical and mental health (4–6).
The genetic divergence between alcohol consumption and AUD may result from confounding, or it may indicate a true difference in the traits’ biology (11). Potential biases include confounding of alcohol consumption frequency by socioeconomic status (12) and the heterogeneity of participants reporting alcohol abstinence (13). However, differences in genetic liability between these traits are observed even after accounting for these confounders. Thus, a latent factor may underlie psychopathology and substance use disorders, but not substance use (14). As AUD is defined by several diagnostic criteria other than the frequency or quantity of alcohol consumption, individuals can exhibit the same level of alcohol consumption while differing with respect to an AUD diagnosis. Consequently, understanding the factors underlying AUD, separate from those that influence normative alcohol consumption, has important theoretical and clinical implications.
We sought to identify genetic variants that contribute uniquely to the development of AUD by conducting multiple cross-ancestry GWASs in the Million Veteran Program (MVP) (15). To elucidate the role of phenotypic variation in genetic discovery, we conducted additional GWASs varying the heterogeneity of the phenotype definition. To examine the transition from alcohol consumption to AUD, we conducted analyses to identify genetic variants that are not mediated through alcohol consumption, but rather directly affect the risk of AUD.
Methods
Overview of Analyses
In this study, we had three aims. First, we sought to identify novel loci associated with AUD and alcohol consumption by performing separate ancestry-specific primary GWASs of AUD and AUDIT-C score, followed by cross-ancestry meta-analyses. Second, we elucidated the impact of phenotypic variation on gene discovery by performing three additional secondary GWASs using 1) a less stringent definition of AUD, 2) AUD with abstainers excluded, and 3) AUDIT-C score with abstainers excluded. Third, we evaluated the mediating effect of alcohol consumption on genetic risk for AUD by conducting a causal mediation analysis that restricted AUD cases to those with an AUDIT-C score that predated their receipt of an AUD diagnosis and removed abstainers from the case and control groups. Figure S1 in the online supplement provides an overview of the phenotype definitions and sample counts for each analysis.
MVP Cohort and Phenotype Construction
The MVP is a large observational cohort study and biobank developed by the Department of Veterans Affairs (VA) (15). AUD diagnostic codes (based on ICD-9 and ICD-10) and AUDIT-C scores were obtained from the VA electronic health record (EHR) system. The AUDIT-C (8) comprises the first three questions of the 10-item AUDIT, a valid, reliable screening instrument for hazardous or harmful drinking (9). We used the maximum AUDIT-C score as a measure of alcohol consumption. The data set comprised all MVP participants with at least one AUDIT-C measure between 2007 and 2018.
Genome-Wide Association Analysis
For ancestry-specific analyses, we tested imputed single-nucleotide polymorphisms (SNPs) that passed quality control (see the Supplemental Methods section in the online supplement ) for association with five alcohol-related phenotypes: 1) AUD–stringent (primary GWAS), 2) AUDIT-C score (primary GWAS), 3) AUD only among individuals with an AUDIT-C score >0 (secondary GWAS), 4) AUDIT-C score only among individuals with an AUDIT-C score >0 (secondary GWAS), and 5) AUD–less stringent (secondary GWAS). The analyses used logistic or linear regression models, as appropriate, in PLINK, version 1.9 (16). Covariates included age at maximum AUDIT-C score, sex, and the first 10 genetic principal components within each ancestry. Cross-ancestry meta-analyses were performed in METAL (17) using the sample size–weighted method. Within each region, independent variants were identified by conditional analyses using GCTA-COJO (18) (for additional details, see the Supplemental Methods section in the online supplement ).
Gene-Based Analyses, Heritability Analyses, and Genetic Correlation
Gene-based association analyses for the primary AUD and AUDIT-C GWAS were performed using MAGMA (19) implemented in FUMA (20). Default settings were used to map input SNPs to 19,082 protein-coding genes. Genome-wide significance was defined as a p value ≥2.62×10−6 (0.05/19,082). SNP-based heritability (h2) for all phenotypes was estimated using linkage disequilibrium score regression (LDSC) (21). Genetic correlation analyses were performed using LDSC (21) and POPCORN (22) (see the Supplemental Methods section in the online supplement ).
Polygenic Risk Scores and Phenome-Wide Association Study
We performed validation analyses in the Vanderbilt Biobank (BioVU; African ancestry, N=12,384; European ancestry, N=66,903) and UK Biobank (African ancestry, N=1,182; European ancestry, N=123,561) by constructing polygenic risk scores (PRSs) for AUD and AUDIT-C score using PRS–continuous shrinkage (PRS-CS) (23) and testing the PRSs against the primary phenotypes available in each data set. The BioVU data set was also used to conduct PheWAS analyses. In each data set, summary statistics from the matched ancestry were used to calculate PRSs; details are provided in the Supplemental Methods section in the online supplement .
Mediation
Causal mediation analyses require that the exposure be associated with the mediator and/or the outcome. Therefore, we selected all variants associated with either AUD–stringent or AUDIT-C score at a suggestive p threshold of 1×10−5 in the cross-ancestry meta-analyses. Among variants present in all three ancestries, 3,092 were selected for AUD and 3,019 for AUDIT-C score. We performed linkage disequilibrium clumping using a range of 3,000 kb and an r2 >0.1 to create a set of 215 variants for AUD and 208 variants for AUDIT-C score, yielding a total set of 372 unique variants associated with either or both traits.
Phenotypic definitions for mediation analysis are provided in full in the Supplemental Methods section in the online supplement . To minimize the possibility of reverse causality, all AUD cases had an AUDIT-C score that predated the AUD diagnosis. For each variant in the set, we estimated the association with AUD and AUDIT-C score using logistic and linear regression models, respectively. To enhance interpretability, the variant associated with greater risk for AUD was used. The direct and indirect effects of each variant on AUD—with AUDIT-C score as the mediator—were estimated using the counterfactual method with natural effect models (24) implemented in the medflex package in R (25, 26). Data expansion was performed using the imputation-based approach. For each variant, two mediation models were implemented—one assuming no interaction between the exposure (variant) and the mediator (AUDIT-C score) and a variant–by–AUDIT-C score interaction model. Analyses were conducted separately for each ancestral group before calculating cross-ancestry meta-analyses in METAL. Variants that passed a Bonferroni-corrected p threshold of 1.34×10−4 (0.05/372) in the AUD–AUDIT-C score >0 analysis (i.e., were associated with the outcome) were explored to partition the total effect into a direct effect and an indirect effect. The proportion mediated (PM) on the log-odds scale was calculated as suggested for a binary outcome: PM=log(ORindirect)/log(ORtotal), where OR=odds ratio and ORtotal=ORdirect×ORindirect (27).
Results
We included 409,630 MVP participants from three ancestral groups: European American (N=296,989), African American (N=80,764), and Hispanic American (N=31,877) (see Table S1 in the online supplement ). Most of the population was male (91.5%). Nearly one-quarter (24.7%) of participants had at least one ICD-9 or ICD-10 code for AUD (i.e., AUD–less stringent; the rate was 21.0% among European Americans, 35.2% among African Americans, and 29.2% among Hispanic Americans), and 20.8% had at least one inpatient or two outpatient codes for AUD (i.e., AUD–stringent; the rate was 18.0% among European Americans, 31.0% among African Americans, and 24.6% among Hispanic Americans). The mean number of annual AUDIT-C screenings was 8.1 (SD=3.8; the mean was 8.1 [SD=3.7] among European Americans, 8.5 [SD=3.9] among African Americans, and 7.6 [SD=3.7] among Hispanic Americans). The mean number of AUDIT-C screenings that predated an AUD diagnosis was 3.3 (SD=2.8; the mean was 3.4 [SD=2.9] among European Americans, 3.2 [SD=2.7] among African Americans, and 3.1 [SD=2.6] among Hispanic Americans). The average maximum AUDIT-C score for all individuals was 3.0 (AUD cases, 5.6; control subjects, 2.2) and occurred at an average age of 59.8 years (see Figure S2 in the online supplement ).
Identification of Cross-Ancestry and Ancestry-Specific Novel Loci for AUD and AUDIT-C
In the cross-ancestry meta-analysis for AUD, 84 lead variants (mapping to 21 loci) were genome-wide significant (p<5×10−8) (Figure 1; see also Table S2 and Figure S3 in the online supplement ). After conditioning on the lead variants within each locus, 26 variants were independently associated with AUD, of which four are novel (no previous reported associations with alcohol-related phenotypes within 1 Mb; see Table S4 in the online supplement ), including two intronic (in ZNF804A and MLN) and two intergenic (near NICN1 and MIR5694). Figure S4 and Table S3 in the online supplement show the 27 lead variants (14 loci) in European Americans, five lead variants (three loci) in African Americans, and four lead variants (one locus) in Hispanic Americans. Of the lead variants identified in European Americans, four variants in three loci were genome-wide significant in the ancestry-specific analysis only (see Figure S5 in the online supplement ), of which one is novel (near MIR129-2; see Table S4 in the online supplement ). All European American ancestry–specific variants were polymorphic in all populations (minor allele frequency [MAF] >10%), except rs34517574, which was absent in the African American sample due to low imputation quality. In African Americans, two lead variants in two loci were ancestry specific (see Figure S6 in the online supplement ), with one novel (in an intron of IL7R; see Table S4 in the online supplement ). Although one of the variants (rs73792114) was rare in European Americans and Hispanic Americans (MAF: African American, 5.9%; European American, 0.04%; Hispanic American, 0.08%), the novel variant (rs7701111) was polymorphic in all populations (MAF: African American, 20.4%; European American, 32.7%; Hispanic American, 42.1%).
FIGURE 1. Manhattan plot for cross-ancestry meta-analysis of AUD and AUDIT-C scorea
aSample size–weighted meta-analyses were performed using METAL (17). In the upper plot, for AUD, Ns are 85,391 for case subjects and 308,488 for control subjects. In the lower plot, for AUDIT-C score, N=409,630. AUD=alcohol use disorder; AUDIT-C=consumption subscale of the Alcohol Use Disorders Identification Test.
In the cross-ancestry meta-analysis for AUDIT-C, 87 lead variants (mapping to 19 loci; see Figure 1; see also Table S5 and Figure S3 in the online supplement ) were identified. Conditional analysis yielded 24 independent variants, including six novel ones (see Table S7 in the online supplement ). Two of the novel variants are intronic (in NEGR1 and CACNG7) and four intergenic (near GBE1, LOC101927021, POT1, HMG20A). In ancestry-specific analyses (see Figure S7 and Table S6 in the online supplement ), 37 variants (14 loci) were genome-wide significant in European Americans, three variants (two loci) were genome-wide significant in African Americans, and four variants (one locus) were genome-wide significant in Hispanic Americans. In European Americans, two variants in two loci were ancestry specific, with one novel and polymorphic in all populations (see Figure S8 and Table S7 in the online supplement ; near COBLL1). In African Americans, there was one ancestry-specific variant, which was novel (see Figure S9 and Table S7 in the online supplement ; near OLFM1); although rare in African Americans, it was almost monomorphic in European Americans and absent in Hispanic Americans (MAF: African American, 0.002%; European American, <0.001%). None of the lead variants in Hispanic Americans were ancestry specific.
We conducted three additional GWASs: 1) AUD in individuals with AUDIT-C score >0; 2) AUDIT-C score in individuals with AUDIT-C score >0; and 3) AUD–less stringent (see Figure S3, Figures S10–S12, and Tables S8–S13 in the online supplement ). In all analyses, correlations between effect sizes for genome-wide significant loci identified in the primary GWAS and the secondary GWAS were high (r2=0.99). Effect sizes between AUD and AUDIT-C traits were highly correlated (Figure 2A).
FIGURE 2. Heritability and genetic correlation analyses for AUD–less stringent, AUD, AUD–AUDIT-C >0, AUDIT-C, and AUDIT-C >0a
aPanel A shows effect sizes from the five cross-ancestry meta-analyses. All are correlated with AUD–stringent. In panel B, heritability is calculated for European Americans and African Americans using linkage disequilibrium score regression (LDSC). Error bars indicate standard error. Panel C shows genetic correlations with other phenotypes in European Americans calculated using LDSC. Nongray bars are significant at a Bonferroni-corrected p value of 1.95×10−4. Error bars indicate 95% confidence interval. AUD=alcohol use disorder; AUDIT-C=consumption subscale of the Alcohol Use Disorders Identification Test.
Gene-Based Analyses
Gene-based association analyses for AUD identified 58 genome-wide significant (p<2.69×10–6) genes in European Americans, 10 in African Americans, and seven in Hispanic Americans (see Figure S13 in the online supplement ). Genes not identified in the SNP-based analyses of AUD (European American, N=47; African American, N=8; Hispanic American, N=4) include ADH4, which was genome-wide significant in both European Americans and African Americans, ADH5, which was genome-wide significant in African Americans, and ADH7, which was genome-wide significant in Hispanic Americans. Gene-based association analyses for AUDIT-C score identified 45 genome-wide significant genes in European Americans, seven in African Americans, and six in Hispanic Americans (see Figure S14 in the online supplement ). In European Americans, 38 of these were not identified in SNP-based analyses, including CACNA1C; six in African Americans, including METAP1; and three in Hispanic Americans, including TRMT10A.
Heritability and Genetic Correlation of AUD and AUDIT-C Phenotypes
SNP-based heritability (h2SNP) for the five related GWAS phenotypes was estimated using LDSC for both European Americans and African Americans (Figure 2B; see also Table S14 in the online supplement ). Heritability values, which were higher for African Americans than for European Americans and similar to previous analyses in MVP participants using smaller samples (4), increased with the specificity of the phenotype. Pairwise genetic correlations (rg) between the five GWAS phenotypes were high within both European Americans and African Americans (see Figure S15 in the online supplement ). Correlations were strongest within the AUD (rg range, 0.92–1) and AUDIT-C (rg range, 0.92–1) traits and weakest between AUD and AUDIT-C (European American rg range, 0.71–0.78; African American rg range, 0.93–0.97), but increased when individuals with AUDIT-C scores of 0 were removed (European American rg range, 0.86–0.92, African American rg range, 0.98–1). Genetic correlations between European Americans and African Americans within traits were estimated using POPCORN (rg range, 0.57–0.74; see Table S15 in the online supplement ).
Genetic Correlations With Other Phenotypes
The genetic correlations between the five GWAS phenotypes and 257 other phenotypes were estimated in the European American population using LDSC (Figure 2C; see also Table S16 in the online supplement ). As in our previous work (4), the magnitude and direction of genetic correlations differed significantly between AUD and AUDIT-C score (see Table S17 in the online supplement ). AUD was correlated with 22 phenotypes, including positive correlations with several psychiatric disorders, insomnia, and neuroticism, and negative correlations with years of schooling and age at first birth. In contrast, AUDIT-C score was correlated with 20 phenotypes, including positive correlations with HDL cholesterol and lung function and negative correlations with anthropometric phenotypes and type 2 diabetes. Phenotypes significantly genetically correlated with AUD in individuals with AUDIT-C score >0 and AUD–less stringent were generally the same as those correlated with AUD, with no significant differences in magnitude of correlation. Fifteen of the phenotypes that were significantly genetically correlated with AUDIT-C score were also correlated with AUDIT-C score in individuals with AUDIT-C score >0, with no significant differences in magnitude. However, AUDIT-C score in individuals with AUDIT-C score >0 was also significantly positively correlated with cross-disorder and schizophrenia phenotypes, and these correlations were not significantly different from those found for AUD.
Polygenic Risk Scores
In the UK Biobank, PRSs for AUD and AUDIT-C score were significantly associated with scores on the AUDIT-C, the AUDIT-P (the problem items in AUDIT), and the AUDIT in individuals of European ancestry, but not African ancestry (see Table S18 in the online supplement ). In BioVU, both PRSs were associated with alcohol-related disorders among individuals of African and European ancestry (see Tables S19–S22 in the online supplement ). Among individuals of European ancestry, PheWAS analyses across 1,335 phenotypes identified an additional 35 phenotypes associated with AUD PRS, including positive associations with “tobacco use disorder” and “mood disorders,” and 42 phenotypes associated with AUDIT-C PRS, including negative associations with “type 2 diabetes” and “obesity” (see Tables S21 and S22 and Figure S16 in the online supplement ).
Mediation
Mediation analysis identified total effects of variants on AUD and AUDIT-C score and direct and indirect effects of variants on AUD (Figure 3; see also Tables S23–S27 in the online supplement ). In the cross-ancestry meta-analysis, 55 variants were associated with AUD at a Bonferroni-corrected p threshold of 1.34×10−4, all with a concordant direction of effect on AUDIT-C score. Twelve of these showed no evidence of mediation (i.e., p values >0.05 for an indirect effect and small values for the proportion of mediation), indicating that AUD susceptibility at these loci was not driven by a genetic predisposition for alcohol consumption. These included variants near LILRB1, LOC157273, CDC5L, LOC101927282, HS3ST3B1, THEG5, BRD3, and LINC01249 and variants in CSMD1, ZBTB16, and SPECC1L-ADORA2A.
FIGURE 3. Mediation analysis of variants associated with AUDa
aOdds ratios for direct and indirect effects are shown for variants significantly associated with AUD. The color of the points denotes the proportion mediated on the log scale (yellow=high mediation, green=moderate mediation, purple=low mediation). Size of the points denotes total effect. Variants are annotated with the nearest gene. AUD=alcohol use disorder.
Forty-three variants showed evidence for mediation by alcohol consumption on their effect on AUD. The strongest evidence for mediation was found for two variants in FTO (in both cases with odds ratiodirect=1.006, pdirect=0.037; odds ratioindirect=1.022, pindirect=1.59×10−13 and 2.70×10−13), followed by a variant near GCKR (odds ratiodirect=1.009, pdirect=0.002; odds ratioindirect=1.020, pindirect=1.05×10−11). The well-known variant rs1229984 in ADH1B showed both a strong direct effect (odds ratiodirect=1.08, pdirect=2.30×10−26) and a strong indirect effect (odds ratioindirect=1.19, pindirect=4.45×10−117) on AUD, with an indirect effect higher than that for the association with AUD (odds ratio=1.14). This variant also showed evidence for an exposure-by-mediator interaction in European Americans (b=−0.03, p=0.025; see Table S27 in the online supplement ), reflecting a decreasing effect of the SNP on AUD as the level of alcohol consumption increases. Other variants with both strong direct and indirect effects on AUD included those near RAP1GDS1 and in SLC39A8.
Discussion
In this study, we performed cross-ancestry GWASs of AUD and AUDIT-C score in 409,630 individuals, including the largest sample of non-European-ancestry individuals for these traits to date. These populations were 40%–50% larger than in our previous GWAS for AUDIT-C score (4) and AUD (4), and we identified novel loci in both cross-ancestry meta-analysis and the ancestry-specific analyses. We conducted mediation analysis to establish a set of variants directly associated with AUD, independent of alcohol consumption. We propose that these variants are high-priority candidates for in vitro study of the biology of the transition from alcohol use to AUD, potentially helping to account for differences in risk of AUD among individuals with the same level of alcohol consumption.
We identified 26 loci associated with AUD, including four novel loci in the cross-ancestry meta-analysis and two novel loci in ancestry-specific analyses. Many of the loci, although novel for alcohol-related phenotypes, have been associated with psychiatric disorders or drug use or use disorders. One novel intronic variant in ZNF804A (Zinc Finger Protein 804A) is located in intron 2, close to a variant previously associated with increased susceptibility to a broad psychosis phenotype (including schizophrenia and bipolar disorder) and heroin addiction (28). The second novel intronic variant is located in MLN (Motilin), previously associated with smoking (7) and depression (29). A novel intergenic variant located upstream of NICN1 (Nicolin 1), which has been associated with depression (30) and intelligence (31), is a known expression quantitative trait locus for multiple other genes (32). The novel African American–specific variant is located in the intron of IL7R (Interleukin 7 Receptor), which plays a critical role in the development of immune cells (33). Chronic alcohol consumption is known to impair immune responses, and concentrations of IL-7 are altered during detoxification (34).
Twenty-two loci were associated with AUDIT-C score, including six novel loci in the cross-ancestry meta-analysis and two novel variants in the ancestry-specific analyses. No prior associations with psychiatric disorders have been identified for these novel loci, consistent with the observed divergence in the genetic architectures of AUDIT-C and AUD (10). Both novel intronic variants are located in genes with prior associations with intelligence and educational attainment (in NEGR1 [Neuronal Growth Regulator 1] and CACNG7 [Calcium Voltage-Gated Channel Auxiliary Subunit Gamma 7]) (31, 35). The novel European American–specific variant near COBLL1 has been associated with multiple phenotypes, including type 2 diabetes (36), whereas variants located nearby have shown associations with HDL cholesterol and triglyceride measures in current drinkers (37). The novel African American–specific variant is located near OLFM1, a gene associated with smoking behavior (7).
We refined both the AUD and AUDIT-C phenotypes by removing individuals who never reported alcohol consumption in the VA EHR. In the UK Biobank and the Veterans Aging Cohort Study (13, 38), where participants’ specific reasons for discontinuing drinking was available, the group not reporting alcohol consumption is heterogeneous, consisting of both lifelong abstainers and former drinkers, some who may have stopped drinking because of health problems, thereby confounding this variable. Confirming this, our analyses demonstrate that 15% of abstainers have an AUD diagnosis in their EHR. Removing them from the analyses did not increase the number of genome-wide significant loci, likely in part because this group makes up almost one-quarter of our original sample. However, by refining the phenotypes to include only individuals who report alcohol consumption, we observed an increase in the SNP heritability for both AUD and AUDIT-C score, and a stronger genetic correlation between AUD and AUDIT-C score. Furthermore, in genetic correlation analyses with other phenotypes, the magnitude of differences between AUD and AUDIT-C score was reduced, revealing positive correlations of AUDIT-C score with psychiatric phenotypes. These findings confirm previous results (12, 13) in which confounding in measures of alcohol consumption contributed in part to the difference in genetic architecture seen between alcohol consumption and AUD.
Evaluating the mediating effect of AUDIT-C score on the genetic risk of AUD identified variants in three genes involved in alcohol metabolism (alcohol dehydrogenase 1B [ADH1B], 1C [ADH1C], and 7 [ADH7]) with both direct and indirect effects on AUD risk. Variants in these genes alter the rate of alcohol metabolism and are consistently associated with protective effects against excessive alcohol consumption, thereby reducing AUD risk (39). The direct (nonmediated) effect of these variants on AUD could result from their effects on alcohol consumption earlier in life, prior to the AUDIT-C assessment, or convey risk of AUD via an independent pathway. Another variant with both strong direct and indirect effects on AUD was in SLC39A8 (solute carrier family 39 member 8; rs13107325), which has been associated with over 60 phenotypes, including anthropometric phenotypes, such as body mass index (BMI) (40), and psychiatric phenotypes, such as schizophrenia (41). Thus, this variant may affect risk for AUD both via alcohol consumption and via effects on the psychopathological dimension of addiction.
The effect of variants in FTO, GCKR, and BANK1 on AUD were shown to be largely mediated by their effect on alcohol consumption. FTO has multiple prior associations with obesity and BMI, with studies suggesting that the mechanism of effect could be through food intake and preference (42). Likewise, both GCKR and BANK1 have prior associations with anthropometric and metabolic phenotypes (43, 44). Thus, their association with AUD may be due to metabolic effects via increased alcohol consumption.
We identified 12 variants with no evidence of mediation through alcohol consumption on AUD risk. We hypothesize that these variants contribute to an addictive dimension of AUD that is not directly influenced by levels of consumption. Some of these were in or near genes previously associated with schizophrenia (CSMD1 [CUB And Sushi Multiple Domains 1]) (41), chronotype (LINC01249 [Long Intergenic Non-Protein Coding RNA 1249]) (45), lung function (CDC5L [Cell Division Cycle 5 Like]) (46), and caffeine consumption (SPECC1L-ADORA2A [Readthrough (NMD Candidate)]) (47). These variants warrant further examination as potential contributors to the transition from heavy drinking to AUD.
Our study has several limitations. First, our measure of alcohol consumption, the AUDIT-C, was self-reported and is subject to recall bias. Second, we systematically captured consumption reported in the VA health system only since 2007; if an individual’s maximum alcohol consumption occurred before then, this measure may not be accurate. Third, different responses to the three AUDIT-C items could result in the same overall score, yet each item may be differentially associated with AUD. In particular, the first item (frequency) is subject to bias and confounding (12). Future studies could examine each AUDIT-C item separately. Finally, the MVP sample is mostly male (>90%), precluding sex-specific analyses. As AUDIT-C scores vary by sex, this will be an important direction for future studies.
In summary, our study identifies novel associations for AUD and AUDIT-C score, including novel ancestry-specific associations. We addressed some of the differences between the genetic architectures of the traits by reducing heterogeneity in the data set. Although genetic variants that contribute to alcohol consumption are important in our understanding of AUD, there is unique variance that remains after refinement of the two traits. Further, we identified a set of variants with direct effects on AUD that are not mediated through alcohol consumption. We propose these as targets for research aimed at understanding the transition from heavy alcohol consumption to AUD, elucidation of the biology of which could inform prevention and treatment efforts.
Supplementary Material
File (appi.ajp.21090892.ds001.pdf)
- View/Download
- 190.33 KB
File (appi.ajp.21090892.ds002.xlsx)
- Download
- 950.04 KB
File (appi.ajp.21090892.ds003.pdf)
- View/Download
- 7.46 MB
References
1.
GBD 2016 Alcohol Collaborators: Alcohol use and burden for 195 countries and territories, 1990–2016: a systematic analysis for the Global Burden of Disease Study 2016. Lancet 2018; 392:1015–1035
2.
American Psychiatric Association: Diagnostic and Statistical Manual of Mental Disorders, 5th ed. Washington, DC, American Psychiatric Association, 2013
3.
Clarke TK, Adams MJ, Davies G, et al: Genome-wide association study of alcohol consumption and genetic overlap with other health-related traits in UK Biobank (N=112 117). Mol Psychiatry 2017; 22:1376–1384
4.
Kranzler HR, Zhou H, Kember RL, et al: Genome-wide association study of alcohol consumption and use disorder in 274,424 individuals from multiple populations. Nat Commun 2019; 10:1499
5.
Zhou H, Sealock JM, Sanchez-Roige S, et al: Genome-wide meta-analysis of problematic alcohol use in 435,563 individuals yields insights into biology and relationships with other traits. Nat Neurosci 2020; 23:809–818
6.
Sanchez-Roige S, Palmer AA, Fontanillas P, et al: Genome-wide association study meta-analysis of the Alcohol Use Disorders Identification Test (AUDIT) in two population-based cohorts. Am J Psychiatry 2019; 176:107–118
7.
Liu M, Jiang Y, Wedow R, et al: Association studies of up to 1.2 million individuals yield new insights into the genetic etiology of tobacco and alcohol use. Nat Genet 2019; 51:237–244
8.
Bush K, Kivlahan DR, McDonell MB, et al: The AUDIT alcohol consumption questions (AUDIT-C): an effective brief screening test for problem drinking. Arch Intern Med 1998; 158:1789–1795
9.
Saunders JB, Aasland OG, Babor TF, et al: Development of the Alcohol Use Disorders Identification Test (AUDIT): WHO collaborative project on early detection of persons with harmful alcohol consumption–II. Addiction 1993; 88:791–804
10.
Sanchez-Roige S, Palmer AA, Clarke TK: Recent efforts to dissect the genetic basis of alcohol use and abuse. Biol Psychiatry 2020; 87:609–618
11.
Mallard TT, Sanchez-Roige S: Dimensional phenotypes in psychiatric genetics: lessons from genome-wide association studies of alcohol use phenotypes. Complex Psychiatry 2021; 7:45–48
12.
Mallard TT, Savage JE, Johnson EC, et al: Item-level genome-wide association study of the Alcohol Use Disorders Identification Test in three population-based cohorts. Am J Psychiatry 2022; 179:58–70
13.
Dao C, Zhou H, Small A, et al: The impact of removing former drinkers from genome-wide association studies of AUDIT-C. Addiction 2021 Nov; 116:3044–3054
14.
Hatoum AS, Johnson EC, Colbert SMC, et al: The addiction risk factor: a unitary genetic vulnerability characterizes substance use disorders and their associations with common correlates. Neuropsychopharmacology 2022; 47:1739–1745
15.
Gaziano JM, Concato J, Brophy M, et al: Million Veteran Program: a mega-biobank to study genetic influences on health and disease. J Clin Epidemiol 2016; 70:214–223
16.
Chang CC, Chow CC, Tellier LC, et al: Second-generation PLINK: rising to the challenge of larger and richer datasets. GigaScience 2015; 4. https://doi.org/10.1186/s13742-015-0047-8
17.
Willer CJ, Li Y, Abecasis GR: METAL: fast and efficient meta-analysis of genomewide association scans. Bioinformatics 2010; 26:2190–2191
18.
Yang J, Ferreira T, Morris AP, et al: Conditional and joint multiple-SNP analysis of GWAS summary statistics identifies additional variants influencing complex traits. Nat Genet 2012; 44:369–375
19.
de Leeuw CA, Mooij JM, Heskes T, et al: MAGMA: generalized gene-set analysis of GWAS data. PLOS Comput Biol 2015; 11:e1004219
20.
Watanabe K, Taskesen E, van Bochoven A, et al: Functional mapping and annotation of genetic associations with FUMA. Nat Commun 2017; 8:1826
21.
Bulik-Sullivan BK, Loh PR, Finucane HK, et al: LD score regression distinguishes confounding from polygenicity in genome-wide association studies. Nat Genet 2015; 47:291–295
22.
Brown BC, Asian Genetic Epidemiology Network Type 2 Diabetes Consortium, Ye CJ, et al: Transethnic genetic-correlation estimates from summary statistics. Am J Hum Genet 2016; 99:76–88
23.
Ge T, Chen CY, Ni Y, et al: Polygenic prediction via Bayesian regression and continuous shrinkage priors. Nat Commun 2019; 10:1776
24.
VanderWeele TJ: Mediation analysis: a practitioner’s guide. Annu Rev Public Health 2016; 37:17–32
25.
Steen J, Loeys T, Moerkerke B, et al. medflex: an R package for flexible mediation analysis using natural effect models. J Stat Softw 2017; 76:1–46
26.
Lange T, Vansteelandt S, Bekaert M: A simple unified approach for estimating natural direct and indirect effects. Am J Epidemiol 2012; 176:190–195
27.
Vanderweele TJ, Vansteelandt S: Odds ratios for mediation analysis for a dichotomous outcome. Am J Epidemiol 2010; 172:1339–1348
28.
Chang H, Xiao X, Li M: The schizophrenia risk gene ZNF804A: clinical associations, biological mechanisms, and neuronal functions. Mol Psychiatry 2017; 22:944–953
29.
Coleman JRI, Peyrot WJ, Purves KL, et al: Genome-wide gene-environment analyses of major depressive disorder and reported lifetime traumatic experiences in UK Biobank. Mol Psychiatry 2020; 25:1430–1446
30.
Adewuyi EO, Mehta D, Sapkota Y, et al: Genetic analysis of endometriosis and depression identifies shared loci and implicates causal links with gastric mucosa abnormality. Hum Genet 2021; 140:529–552
31.
Hill WD, Marioni RE, Maghzian O, et al: A combined analysis of genetically correlated traits identifies 187 loci and a role for neurogenesis and myelination in intelligence. Mol Psychiatry 2019; 24:169–181
32.
GTEx Consortium: The GTEx Consortium atlas of genetic regulatory effects across human tissues. Science 2020; 369:1318–1330
33.
Kang J, Coles M: IL-7: the global builder of the innate lymphoid network and beyond, one niche at a time. Semin Immunol 2012; 24:190–197
34.
Nikou T, Ioannidis A, Zoga M, et al: Alteration in the concentrations of interleukin-7 (IL-7), interleukin-10 (IL-10), and granulocyte colony stimulating factor (G-CSF) in alcohol-dependent individuals without liver disease, during detoxification therapy. Drug Alcohol Depend 2016; 163:77–83
35.
Lee JJ, Wedow R, Okbay A, et al: Gene discovery and polygenic prediction from a genome-wide association study of educational attainment in 1.1 million individuals. Nat Genet 2018; 50:1112–1121
36.
Morris AP, Voight BF, Teslovich TM, et al: Large-scale association analysis provides insights into the genetic architecture and pathophysiology of type 2 diabetes. Nat Genet 2012; 44:981–990
37.
de Vries PS, Brown MR, Bentley AR, et al: Multiancestry genome-wide association study of lipid levels incorporating gene-alcohol interactions. Am J Epidemiol 2019; 188:1033–1054
38.
Gordon KS, McGinnis K, Dao C, et al: Differentiating types of self-reported alcohol abstinence. AIDS Behav 2020; 24:655–665
39.
Edenberg HJ, McClintick JN: Alcohol dehydrogenases, aldehyde dehydrogenases, and alcohol use disorders: a critical review. Alcohol Clin Exp Res 2018; 42:2281–2297
40.
Speliotes EK, Willer CJ, Berndt SI, et al: Association analyses of 249,796 individuals reveal 18 new loci associated with body mass index. Nat Genet 2010; 42:937–948
41.
Pardiñas AF, Holmans P, Pocklington AJ, et al: Common schizophrenia alleles are enriched in mutation-intolerant genes and in regions under strong background selection. Nat Genet 2018; 50:381–389
42.
Loos RJF, Yeo GSH: The bigger picture of FTO: the first GWAS-identified obesity gene. Nat Rev Endocrinol 2014; 10:51–61
43.
Sakaue S, Kanai M, Tanigawa Y, et al: A cross-population atlas of genetic associations for 220 human phenotypes. Nat Genet 2021; 53:1415–1424
44.
Onuma H, Tabara Y, Kawamoto R, et al: The GCKR rs780094 polymorphism is associated with susceptibility of type 2 diabetes, reduced fasting plasma glucose levels, increased triglycerides levels, and lower HOMA-IR in Japanese population. J Hum Genet 2010; 55:600–604
45.
Jones SE, Lane JM, Wood AR, et al: Genome-wide association analyses of chronotype in 697,828 individuals provides insights into circadian rhythms. Nat Commun 2019; 10:343
46.
Shrine N, Guyatt AL, Erzurumluoglu AM, et al: New genetic signals for lung function highlight pathways and chronic obstructive pulmonary disease associations across multiple ancestries. Nat Genet 2019; 51:481–493
47.
Said MA, van de Vegte YJ, Verweij N, et al: Associations of observational and genetically determined caffeine intake with coronary artery disease and diabetes mellitus. J Am Heart Assoc 2020; 9:e016808
Information & Authors
Information
Published In
History
Received: 8 September 2021
Revision received: 10 November 2022
Revision received: 2 March 2023
Accepted: 17 March 2023
Published online: 7 June 2023
Published in print: August 01, 2023
Keywords
Authors
Funding Information
Supported by Merit Review Award I01 BX003341 from the U.S. Department of Veterans Affairs Biomedical Laboratory Research and Development Service; the VISN 4 Mental Illness Research, Education, and Clinical Center; and National Institute on Alcohol Abuse and Alcoholism grant K01-AA028292 (to Dr. Kember).The views expressed in this article are those of the authors and do not necessarily represent the position or policy of the U.S. Department of Veterans Affairs or the U.S. government.Dr. Kranzler is a member of advisory boards for Clearmind Medicine, Dicerna Pharmaceuticals, Enthion Pharmaceuticals, and Sophrosyne Pharmaceuticals; he has served as a consultant for Sobrera Pharmaceuticals; he has received research funding and medication supplies for an investigator-initiated study from Alkermes; and he is a member of the American Society of Clinical Psychopharmacology’s Alcohol Clinical Trials Initiative, which was supported in the past 3 years by Alkermes, Dicerna, Ethypharm, Lundbeck, Mitsubishi, Otsuka, and Pear Therapeutics. Drs. Gelernter and Kranzler are holders of U.S. patent 10,900,082 (“Genotype-guided dosing of opioid agonists”). The other authors report no financial relationships with commercial interests.
Metrics & Citations
Metrics
Citations
Export Citations
If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click Download.
For more information or tips please see 'Downloading to a citation manager' in the Help menu.
There are no citations for this item
View Options
View options
PDF/ePub
View PDF/ePubGet Access
Login options
Already a subscriber? Access your subscription through your login credentials or your institution for full access to this article.
Personal login Institutional Login Open Athens loginNot a subscriber?
PsychiatryOnline subscription options offer access to the DSM-5-TR® library, books, journals, CME, and patient resources. This all-in-one virtual library provides psychiatrists and mental health professionals with key resources for diagnosis, treatment, research, and professional development.
Need more help? PsychiatryOnline Customer Service may be reached by emailing [email protected] or by calling 800-368-5777 (in the U.S.) or 703-907-7322 (outside the U.S.).