Genetic influences on alcohol flushing in East Asian populations

Background Although it is known that variation in the aldehyde dehydrogenase 2 (ALDH2) gene family influences the East Asian alcohol flushing response, knowledge about other genetic variants that affect flushing symptoms is limited. Methods We performed a genome-wide association study meta-analysis and heritability analysis of alcohol flushing in 15,105 males of East Asian ancestry (Koreans and Chinese) to identify genetic associations with alcohol flushing. We also evaluated whether self-reported flushing can be used as an instrumental variable for alcohol intake. Results We identified variants in the region of ALDH2 strongly associated with alcohol flushing, replicating previous studies conducted in East Asian populations. Additionally, we identified variants in the alcohol dehydrogenase 1B (ADH1B) gene region associated with alcohol flushing. Several novel variants were identified after adjustment for the lead variants (ALDH2-rs671 and ADH1B-rs1229984), which need to be confirmed in larger studies. The estimated SNP-heritability on the liability scale was 13% (S.E. = 4%) for flushing, but the heritability estimate decreased to 6% (S.E. = 4%) when the effects of the lead variants were controlled for. Genetic instrumentation of higher alcohol intake using these variants recapitulated known associations of alcohol intake with hypertension. Using self-reported alcohol flushing as an instrument gave a similar association pattern of higher alcohol intake and cardiovascular disease-related traits (e.g. stroke). Conclusion This study confirms that ALDH2-rs671 and ADH1B-rs1229984 are associated with alcohol flushing in East Asian populations. Our findings also suggest that self-reported alcohol flushing can be used as an instrumental variable in future studies of alcohol consumption. Supplementary Information The online version contains supplementary material available at 10.1186/s12864-023-09721-7.


Background
Alcohol flushing is a heritable condition in which a person develops flushes on the face or skin after drinking alcohol.Whilst pronounced alcohol flushing is rarely observed in Europeans, approximately 36% of East Asians experience alcohol flushing as well as other unpleasant symptoms (e.g.nausea and tachycardia) [1].Previous genome-wide association studies (GWAS) identified two key genes associated with alcohol flushing, alcohol dehydrogenase 2 (ALDH2) and aldehyde dehydrogenase 1B (ADH1B) [2][3][4].These genes encode enzymes that metabolize alcohol into acetaldehyde (ADH1B) and acetaldehyde into acetate (ALDH2).Genetic variants in ALDH2 and ADH1B alter alcohol metabolism leading to prolonged, elevated levels of acetaldehyde.The excess acetaldehyde leads to physiological responses to alcohol consumption, including erythema on the face, nausea, and rapid heart rate [5,6].
Most previous GWAS have focused on genetic associations with alcohol drinking status, rather than alcohol-induced responses, such as alcohol flushing [7,8].Candidate gene association studies have provided evidence for the association of ALDH2 or ADH1B with alcohol flushing [9], but it is unclear whether there are loci other than ALDH2 or ADH1B at which genetic variation appreciably influences flushing symptoms.Furthermore, investigations of putative causal genes for alcohol-related physiological responses have been conducted almost exclusively in individuals of European ancestry to date [7,10], which risks missing variants with very low frequencies in European populations.Genetic biobanks from East Asian populations are growing in number, and with alcohol flushing highly prevalent amongst those participants there is an opportunity to improve our understanding of the relevant risk variants for the condition.
Recently, alcohol flushing has been proposed as a phenotypic instrumental variable (IV) for examining the health impacts of alcohol consumption in East Asian populations [11,12].Alcohol flushing is associated with lower levels of alcohol consumption and is assumed to be independent of confounders [13].Considering the ease of including alcohol flushing questions in surveys compared with collecting genetic information, using flushing as an IV may be beneficial, enabling IV analysis in a simple, cost-effective, and non-invasive manner.Therefore, it would be helpful to fully understand the effects of genetic variants on alcohol flushing and to further characterise its utility as an IV.
In this study, we perform the largest GWAS of alcohol flushing to date, using 15,016 male individuals of East Asian ancestry from the China Kadoorie Biobank (CKB; N = 13,456) and the Korean Genome and Epidemiology Study (KoGES; N = 1,560).We also estimated the SNPbased heritability of alcohol flushing.Furthermore, we examined whether self-reported alcohol flushing can be used as a phenotypic IV for alcohol intake, comparing estimates with results from the genotypic IV (rs671 in ALDH2).

Study population
This study was performed on two datasets, CKB (discovery set) and KoGES (replication set).CKB is a prospective study that recruited participants between 2004 and 2008.At baseline, 512,726 adults aged 30-79 years were recruited from 10 geographically defined regions of China (5 urban and 5 rural areas).All participants provide a 10mL blood sample which was processed into aliquots of buffy coat and plasma and stored at -70 °C.Participants were prospectively followed up for causespecific morbidity and mortality through linkage to death and disease registries and to the national health insurance system.Detailed information on the CKB is provided elsewhere [14,15].For the current analyses, we excluded individuals who were not genotyped or nondrinkers for whom information on alcohol flushing was not collected (Fig. 1).Individuals with non-local ancestry were excluded from region-stratified GWAS analyses.Analyses were limited to male participants only since female participants' alcohol intake is very low in China [16] and South Korea [17].In total, 13,456 male CKB participants were included in regional GWAS analyses.For the meta-analysis, data for a total of 1,560 Korean men were obtained from KoGES [18].For the IV analysis, we included 23,020 males from CKB who have information on alcohol flushing, alcohol intake amount and the known genetic instrument for alcohol (rs671 in ALDH2; Fig. 1).All participants provided written informed consent approved by relevant local, national, and international ethics committees.Detailed information on the samples is provided in Supplementary Data.

Assessment of alcohol flushing and drinking patterns
In CKB, alcohol drinking patterns were investigated using interviewer-administered questionnaires.Participants were asked how often they had drunk alcohol during the previous 12 months (never or almost never; occasionally; only at certain seasons; every month but less than weekly; usually at least once a week).Based on the questionnaire, individuals who reported alcohol consumption in most weeks in the past year were identified as current drinkers.Current drinkers were asked further questions including types of beverage consumed, amount of alcohol drunk, and experience of flushing after drinking.Total alcohol intake (g/day) was calculated using the average alcohol content of each type of alcoholic beverage.Detailed information on the assessment of alcohol intake is available elsewhere [16,19].To investigate the presence of alcohol flushing symptoms among current drinkers, the following question was used: "Do you usually experience hot flushes or dizziness after drinking?"Participants were offered four options: "Yes, immediately"; "Yes, after a small amount of alcohol"; "Yes, but only after drinking a large amount of alcohol", and "No".Participants who experienced flushing immediately after drinking alcohol and those who flushed after a small amount of alcohol were classified as alcohol flushers.For sensitivity analyses, we defined alcohol flushing using different criteria (main, relaxed, strict, and continuous; see the Methods section in Supplementary Data for more details).All questionnaires were provided in Mandarin.The definition of flushing for KoGES is described in Supplementary Data.

DNA sampling and genotyping
DNA was extracted from the buffy coat and was genotyped using the custom Affymetrix Axiom arrays and Illumina Golden Gate platform at BGI (Shenzhen, China), as previously described [15].Data for a total of 100,706 individuals passed quality control criteria (call rate ≥ 95%, no sex mismatch, heterozygosity F statistic SD score < + 3, no XY aneuploidy, no non-East Asian ancestry).Following variant QC (call rate > 0.98, no batch or plate effect, Hardy-Weinberg equilibrium P > 10 − 6 ), imputation was performed using SHAPEITv3/ IMPUTEv4 and the 1000 Genomes Project Phase 3 reference panel.After imputation, SNPs were removed if the MAF was low (< 0.01) or INFO was < 0.3.After QC, 8,001,732 autosomal SNPs were used for association testing.Detailed information on the genotyping method and QC for KoGES is provided in Supplementary Data.

Genome-wide association analyses
In CKB, genetic loci associated with flushing were investigated using BOLT-LMM v2.3.2 [20].Three models were constructed.The first model was adjusted for age, age squared, the first ten genetic principal components (PCs), and genotyping array version (Model 1).We performed second and third GWAS analyses adjusting for the dosages of the SNPs that are known to be strongly associated with alcohol metabolism -rs671 in ALDH2 (Model 2) and additionally rs1229984 in ADH1B (Model 3) [12].We performed further GWAS analyses using different definitions of alcohol flushing (Supplementary Data).Each of the GWA analyses described above was performed separately for each geographical region (10 study areas).Within each region, SNPs with a low minor allele count (MAC < 6) or with Hardy-Weinberg equilibrium test values of P < 1 × 10 − 6 were excluded.Betas and standard errors (S.E.) obtained from BOLT-LMM were converted to log-odds ratios (OR) using log(OR) = β/(µ(1 − µ)), where µ is the case-control ratio, following which region-level association statistics were combined using a fixed-effect inverse-variance-weighted meta-analysis using METAL [21].One region (region 46, Liuzhou; n = 682) was excluded from the meta-analysis since the heritability estimate in this region was close to 0. We did not apply genomic control correction to the meta-analysis data because there was little evidence for inflation (all λ < 1.02, Fig. 2).
In KoGES, association tests were performed using PLINK 1.90 (available at https://www.cog-genomics.org/plink2).The GWA analysis of alcohol flushing was conducted using logistic regression assuming an additive genetic model using the three constructed models described above (Supplementary Data).SNPs with a low minor allele count (MAC < 20) were excluded.
For the GWAS meta-analysis of CKB and KoGES, we performed a fixed-effect inverse variance-weighted metaanalysis of the GWAS summary statistics from the CKB and KoGES using METAL [21].
For all GWAS analyses, a genome-wide significance threshold of 5.0 × 10 − 8 was applied.We presented variants that were identified to be independent after linkage disequilibrium (LD) clumping (Supplementary Data).The distributions of the observed P-values of given SNPs were plotted against the theoretical distribution of expected P-values to yield a quantile-quantile (QQ) plot for flushing (Fig. 2).

Single nucleotide polymorphism-heritability analysis
The SNP heritability of alcohol flushing in the CKB sample was calculated using BOLT-REML, which provides a fast algorithm for multi-component modelling to partition SNP-heritability [22].Heritability (h 2 g ) was estimated using the restricted maximum likelihood estimation method implemented in BOLT-REML.Since we defined alcohol flushing as a binary trait, we transformed the heritability on the observed scale to that on the liability scale (h 2 l ) [23].Analyses were adjusted for the covariates used in the GWAS analyses.SNP heritability in KOGES was estimated using the bivariate restricted maximum likelihood analysis implemented in GCTA [24,25].Detailed methods are described in the Supplementary Data.

Mendelian randomisation analysis of alcohol flushing and disease outcomes
The causal effect of alcohol intake on blood pressure and cardiovascular diseases and related traits was evaluated using IV analyses with a two-stage least squares estimation method.A total of 23,020 individuals were included in the IV analyses (Fig. 1).Self-reported alcohol flushing and the rs671 genotype were used as the phenotypic and genotypic instruments, respectively.We used the strict definition of flushing (i.e., immediately after consuming alcohol) as our IV.The magnitude of the association of alcohol intake (g/week) was scaled into a 280 g/week effect, as in a previous study [16].For binary outcomes (i.e.stroke, myocardial infarction, coronary heart disease, hypertension, and diabetes), a two-stage logistic model was used.In the first stage, alcohol intake was instrumented by alcohol flushing or the rs671 genotype with adjustment for age, region, PCs (1-10), and genotyping array, using a linear regression model.In the second stage, the effect of alcohol on the risk of disease was estimated by fitting the alcohol intake value from the first stage, under a logistic regression model with adjustment for the same confounders as in the first stage.For continuous traits (i.e., aspartate aminotransferase [AST], gamma-glutamyl transferase [GTT], cholesterol, triglycerides, blood glucose, and blood pressure), a two-stage linear model was applied, similarly adjusting for confounders.Region-stratified analyses followed by metaanalysis gave similar results.
The values were reported as ORs per 280 g/week alcohol intake with 95% CIs for the binary outcomes and β-coefficients with 95% CIs for the continuous outcomes.We examined the strength and validity of each instrument using the F-statistic of the association of each instrument with alcohol intake (with an F-statistic > 10 indicating adequate strength).Statistical significance (at the 5% level) was evaluated using a P-value threshold of 0.05.The difference of estimates between instruments (alcohol flushing and rs671) was assessed using a difference of two means test [26] (P value threshold for significance = 0.05).

General characteristics of the study population
The baseline characteristics of the study subjects according to flushing status are presented in Supplementary Tables 1 and 2. In the CKB cohort, among 13,456 men with both alcohol flushing and genotype information, 17.9% reported flushing (i.e., flushing immediately after drinking alcohol or after drinking a small amount of alcohol).The mean weekly alcohol intake of non-flushers was 304.5 ± 259.0 g/week (mean ± standard deviation [SD]).Flushers had a lower mean weekly alcohol intake (228.1 ± 259.0 g/week) compared to non-flushers.Flushers had a higher proportion of rs671 A allele carriers (45.5% of flushers vs. 8.7% of non-flushers) as well as rs1229984 A allele carriers (90.3% of flushers vs. 87.3% of non-flushers) than non-flushers.The characteristics of 1,560 KoGES samples are described in Supplementary Table 2. Similar to the CKB, flushers in KoGES had a lower proportion of current drinkers who consumed relatively small amounts of alcohol compared to nonflushers.Also, flushers in KoGES had a higher proportion of rs671 A allele carriers (68.4% of flushers vs. 9.1% of non-flushers) and rs1229984 A allele carriers (95.5% of flushers vs. 93.2% of non-flushers) than non-flushers.
GWA analyses using different criteria for defining flushing showed no difference in the top signals for Models 1 and 2 across the different definitions of flushing (see Supplementary Methods) although the P-values for the lead SNPs varied (Table 1; Supplementary Figs.1-3; Supplementary Tables 10-16); The P values for the strongest signals became less significant for the relaxed flushing definition (i.e., flushing after drinking any amount of alcohol) (Table 1; Supplementary Tables 10-11).For the relaxed flushing definition, Model 2 identified additional signals on chromosome 2 (rs532522882 HPCAL1; P = 1.29 × 10 − 8 ) along with the signal at ADH1B on chromosome 4 (Table 1; Supplementary Table 11).For the strict flushing definition (i.e., flushing immediately after drinking alcohol), Model 3 identified a few rare variants (MAF < = 0.01; Table 1 and Supplementary Table 14) that reached genome-wide significance including rs150099059 in KCNH1 (P = 9.4 × 10 − 9 ), rs1011755 on chromosome 11 (P = 1.6 × 10 − 8 ), and rs142761523 in CNTN (P = 2.6 × 10 − 8 ).For each flushing definition, Model 3 also identified further suggestive associations marginally below the genome-wide significance threshold.These include rs148407052 in LOC105375361 (P = 5.1 × 10 − 7 ) for the relaxed flushing definition; and rs2903308 in SHISA9 (P = 1.4 × 10 − 7 ) for the continuous flushing definition.However, we were not able to replicate these findings in KoGES: either the association of these variants was strongly attenuated towards the null, or they were not available in KoGES (Supplementary Table 6).
The GWAS results from an independent Korean cohort (KoGES) are presented in Supplementary Tables 3 and   4. The GWAS identified strong association signals on chromosome 12 including rs671.In KoGES, ADH1B rs1229984 did not reach genome-wide significance across models 1-2.An apparent independent association at the chromosome 12 locus harbouring the ALDH2 gene was identified after adjusting for rs671 (rs2074356, beta = 2.85, S.E = 0.26, 2.7 × 10 − 28 ; Model 2; Supplementary Fig. 4 and Supplementary Table 4), or adjusting for rs12231737, which was the top signal obtained from Model 1 (rs2074356, beta = 2.26, S.E = 0.28, 2.9 × 10 − 16 ; Model 4; Supplementary Table 4).To explore the obtained signals further, we conducted fine mapping using SuSiE which returned a single credible set.The credible set suggested that the conditionally independent signals are likely due to measurement error induced by relatively low imputation quality around the rs671 locus (data available on request).
A summary of the strongest association signals from the meta-analysis is presented in Supplementary Tables 3 and 17-19.

SNP heritability for alcohol flushing in the CKB and KoGES
SNP heritability of alcohol flushing among drinkers was estimated to be 12.6% (SE = 4.0%) on the liability scale (h 2 l ) It decreased to 8.4% (S.E.= 4.2%) when we controlled Flushing and SNPs were regarded as dependent and independent variables, respectively. 3The beta estimates, standard errors, and P-values were obtained from the linear regression for rs671 in ALDH2 (Supplementary Table 5), and decreased further when we also controlled for rs1229984 in ADH1B (h 2 l = 6.3%;S.E.= 4.2%), suggesting that rs671 and rs1229984 together explain half of the common variant genetic variance in alcohol flushing in Chinese males.SNP heritability estimates of alcohol flushing amongst drinkers and non-drinkers in the Korean population were imprecise due to the relatively small sample size but showed a pattern consistent with that seen in CKB.

Using self-reported flushing as an instrumental variable
IV analyses among 23,020 men in CKB with flushing data showed that higher alcohol intake (as instrumented by absence of self-reported alcohol flushing) was nominally associated with a higher risk of intracerebral haemorrhage (OR per 280 g/week increase in alcohol intake = 3.28; 95% CI = 1.58-6.81),and total stroke (OR per 280 g/day increase in alcohol intake = 1.89; 95% CI = 1.28-6.81)as well as higher levels of AST, GGT, HDL cholesterol, log-transformed random blood glucose, and diastolic blood pressure (DBP; beta per 280 g/ day increase in alcohol intake = 2.3 mm Hg; 95% CI = 0.9-3.7;Table 2).These associations were generally consistent in direction and magnitude, although the estimates were more precise when using the rs671 genotype as an IV, which also provided evidence that higher alcohol intake caused a higher risk of hypertension and higher levels of systolic blood pressure (SBP), as well as increased risk of stroke types, coronary heart disease, and diabetes.

Discussion
In this study, we investigated genetic variation associated with alcohol flushing and estimated the heritability of flushing in Chinese and Korean male populations.Strong signals were detected in ALDH2 (Supplementary Table 3) in both populations, supporting the previous evidence [27].The SNP-based heritability estimate on the liability scale was 13% for flushing and decreased by 6% when the key variants (rs671 and rs1229984) were accounted for.The decrease in heritability supports the role of ALDH2 and ADH1B as major contributors to the self-reported alcohol flushing response in the Chinese and Korean populations.
In both cohorts (CKB and KoGES), a small proportion of non-flushers were carriers of ALDH2-rs671 A, whilst some flushers were not A allele carriers, suggesting that other genetic variants may play a role in alcohol flushing metabolism.Therefore, we adjusted for the ALDH2 1 ORs and beta coefficients by instrumental variable (IV) estimation were obtained from IV regressions with a two-stage least squares estimation method (in logistic and linear regression models, respectively), using alcohol flushing as an instrument for alcohol intake.To predict the amount of alcohol intake (280 g/week), nonflushers were regarded as a reference group. 2 P values were derived from IV regression analysis with adjustment for age, region, and PC scores 1-10.The difference of estimates between instruments (alcohol flushing and rs671) was assessed using a difference of two means test rs671 genotype to identify other variants that may influence alcohol flushing: this revealed a strong association of ADH1B rs1229984 with alcohol flushing: this revealed a strong association of ADH1B rs1229984 with alcohol flushing.rs1229984 is a missense variant that has been extensively reported to be associated with alcohol consumption phenotypes such as alcohol intake status, and alcohol use disorders, including in European populations where the variant is present at low-frequency [28][29][30].
There has been some disagreement relating to the association of ADH1B with alcohol flushing.A low-dose alcohol challenge followed by a metabolite screen in Han Chinese men suggested that ADH1B did not associate with elevated blood acetaldehyde [31].However, in a candidate gene study involving ALDH2 and ADH1B in a sample of Japanese individuals with alcohol dependence, ADH1B did associate with flushing [32].In CKB, the power to detect the ADH1B association is improved by reducing the residual variance after conditioning on rs671.However, the ADH1B association did not reach statistical significance in the Korean population.One theoretical explanation for that result is collider bias [33], in which flushing and ADH1B each influence alcohol dependence independently [32], and amongst cases become associated.Here, the ADH1B association is unlikely to arise due to this form of technical issue, because the association replicates in KoGES (albeit not at genome-wide significance) which has no alcohol consumption-related sample selection.Further GWAS in larger samples are required given the sample size of KoGES.
Several low-frequency variants were associated with different definitions of alcohol flushing in CKB (Table 1; Supplementary Tables 9-16), after controlling for the known variants (ALDH2 rs671 and ADH1B rs1229984).These include PTPRG rs1508403 (MAF = 0.013) for the main flushing definition (Supplementary Table 9), HPCAL1 rs532522882 (MAF = 0.004) and rs181957632 (MAF = 0.004) for the relaxed flushing definition (Supplementary Table 11), and KCNH1 rs150099059 (MAF = 0.01), and rs142761523 (MAF = 0.01) and rs144350123 in CNTN (MAF = 0.01) for the strict flushing definition (Supplementary Table 13).A GWAS study in 3,838 individuals of European-and African-American ancestry reported that the activities of PTPRG were associated with alcohol dependence [34].A study in mice reported that the expression of HPCAL1 was associated with alcohol consumption [35].Furthermore, a study in rats reported that the KCNH1 gene, which encodes potassium voltage-gated channels, is differentially expressed in binge drinking groups [36].The CNTN family has been suggested to be associated with alcohol independence by GWAS studies in European populations [37,38].Further studies with larger samples will be needed to replicate these findings.SNP-based heritability analyses estimated that around 13% of the phenotypic variation in flushing is explained by common genetic variants.The heritability estimates decreased substantially when ALDH2 rs671 was controlled for illustrating the strong effect of ALDH2 on flushing in the Chinese population.These heritability estimates for flushing were much lower than all previous estimates for alcohol consumption [39].One reason could be that our study only included regular drinkers.In this study, the subjects were asked about their experience of flushing based on their alcohol drinking status.This can be a source of selection bias where a sample can contain only those who report drinking.For example, individuals from CKB who do not regularly drink due to their knowledge of flushing are likely excluded from the current analysis.Also, individuals who drink regardless of their flushing symptom may have developed compensatory feedback mechanisms [40], which can possibly contribute to weaker flushing symptoms.Consequently, this may lead to lower variance in flushing severity in the study subjects that could lead to lower heritability estimates in Chinese population.
The IV results demonstrated that self-reported alcohol flushing can be used as an IV for alcohol consumption levels among drinkers.The pattern of associations of alcohol and disease traits was similar to a previous study in the Korean population that suggested the possibility of using self-reported alcohol flushing as an IV [11,41].However, we observed that the power to detect causal effects was generally attenuated in CKB when using selfreported flushing compared with the genetic IV, whereas the previous study by Cho et al. [41] demonstrated using self-reported alcohol flushing as an IV gave similar results to the use of the ALDH2 rs671 variant as an IV.One major difference between the two studies is that CKB only had data available on alcohol flushing amongst individuals who self-reported regular drinking.Such structured sample selection can induce collider bias [33].Indeed, in the CKB, the participants who regularly consumed alcohol had a lower prevalence of hypertension and lower BP levels than non-drinkers or ex-drinkers (Supplementary Table 7).This suggests that the IV analysis in CKB may have been affected by collider bias.For example, if higher levels of BP and flushing are both causally related to drinking, the association between alcohol intake and higher BP may be distorted (Supplementary Fig. 7), given non-drinkers who flush were excluded from the current study.In this case, the genetic instrument (e.g.rs671) for the overall population is likely to be more reliable than a questionnaire as the genotypes are distributed completely randomly within the whole sample, regardless of their drinking status.By contrast, the selfreported IV based on the questionnaire is more likely to be subject to individuals' drinking status.This study has several other limitations.First, despite this being the largest genome-wide study of alcohol flushing to date, it is possible that there was limited statistical power to detect influential loci other than ALDH2 and ADH1B.Second, our analyses included flushers who regularly drink, due to the design of the questionnaire used in CKB.Therefore, there is a possibility that those who do not drink alcohol due to their response to alcohol were not included in the current study.Nonetheless, results for our top loci are confirmed in two independent samples (Chinese and Koreans) showing that the identified genetic variants are likely to be strongly involved in flushing.Further GWAS and SNP heritability analyses are required in other East Asian populations.Third, some variants identified in CKB were relatively rare, and we could not test their association in KoGES, leaving the possibility that these variants were detected by chance.Fourth, although the variants used for GWAS were filtered to have high imputation scores (INFO > = 0.8), imputation accuracy using the 1000 genomes reference panel in Korean samples as was done for KoGES may still lead to measurement error.This is because, although the panel includes East Asian samples (Han Chinese and Japanese), it does not include Korean samples.It has been reported that the Korean population is genetically homogeneous due to geopolitical isolation, thus, Koreans genetically clustered distinctly from other East Asian populations [42].Therefore, it could be speculated that while rs671 associated very strongly with flushing, it was not detected as the top signal at the ALDH2 locus due to inaccuracy in imputation.Fifth, the use of alcohol flushing as an instrument may only reflect an effect of alcohol intake from a specific period of the life course (e.g. in adulthood) since alcohol flushing only occurs after an individual has started drinking (e.g. during adulthood).

Conclusions
Despite these limitations, the results have epidemiologic and public health implications.Our findings underline the importance of additive genetic effects in modifying alcohol consumption behaviour and support the use of flushing or genetic variants (e.g.rs671 in ALDH2) as proxies for alcohol consumption in East Asian populations.To the best of our knowledge, this is the first GWAS to investigate putative causal variants for alcohol flushing and estimate the heritability of the condition in East-Asian populations.

Fig. 1
Fig. 1 Flowchart of study population selection

Table 2
Associations of alcohol intake and disease traits using alcohol flushing or rs671 as instrumental variables