The role of aldehyde dehydrogenase-1 (ALDH1A1) polymorphisms in harmful alcohol consumption in a Finnish population.

Liver cystolic aldehyde dehydrogenase 1 (ALDH1A1) has been previously associated with both alcohol dependence and alcohol consumption behaviour, and has been implicated in alcohol-induced flushing and alcohol sensitivity in Caucasians. The present study tested for association between ALDH1A1 and alcohol consumption behaviour and susceptibility to problem drinking or alcohol dependence in Finnish cohorts of unrelated male subjects recruited from alcoholism clinical treatment facilities ( n = 104) and from the general population ( n = 201). All participants completed the Alcohol Use Disorder Identification Test (AUDIT) and were genotyped for eight single nucleotide polymorphisms (SNPs) within or flanking ALDH1A1 . To test for association between alcohol consumption behaviour and these polymorphisms, we used generalised linear models and haplotypic analysis. Three SNPs were nominally associated (rs348449, p = 0.043; rs610529, p = 0.013; rs348479, p = 0.025) with the quantitative AUDIT score, which evaluates alcohol consumption behaviour. Two-locus (rs610529-rs2288087) haplotype analysis increased the strength of association with AUDIT score ( p = 0.0015). Additionally, rs348449 is highly associated with problem drinking (allelic odds ratio [OR] 7.87, 95 per cent confidence interval [CI] 1.67-37.01) but due to the low minor allele frequency (0.01 and 0.07 in controls and problem drinkers, respectively), more samples are required to validate this observation. Conversely, rs348479 ( p = 0.019) and rs610529 (allelic OR 0.65, 95 per cent CI 0.43-0.98; genotypic OR 0.32, 95 per cent CI 0.12-0.84) are implicated in alcohol dependence status. This study provides further evidence for a role for ALDH1A1 in alcohol consumption behaviour, including problem drinking and possibly alcohol dependence, in our Finnish population.

Introduction

Cystolic aldehyde dehydrogenase (ALDH1A1) has been associated with alcohol-induced flushing in Caucasians, and alcohol sensitivity and dependence [1,2], and functions in the detoxification of acetaldehyde, the first metabolite of ethanol oxidation [3]. Additionally, ALDH1A1 is involved in the degradation of dopamine (DA) in the ventral tegmental area (VTA) [4]. This role in the regulation of DA degradation may be important in the development and maintenance of excessive alcohol consumption because VTA DA function can be ethanol induced via increased DA neurone firing, and the VTA lies within the mesolimbic system, which functions in reward, motivation and various substance disorders, including alcohol dependence [5,6].

Association between genetic variation in ALDH1A1 (on chromosome 9q21.13) and inter-individual differences in the vulnerability to alcohol dependence and other alcohol-related phenotypes has been investigated in several studies. While there are three alternative splice patterns producing two ALDH1A1 isoforms [7], ALDH1A1 is highly conserved in humans. Two promoter polymorphisms have recently been correlated with the development of alcohol dependence in south-west Californian Indians [8] and other populations [1,9]. Spence and colleagues [1] sequenced the ALDH1A1 promoter region in Asian, Caucasian and African-American alcoholics and control subjects. They discovered two rare promoter polymorphisms, ALDH1A1*2 and ALDH1A1*3 (only detected in African Americans). While no association was observed between ALDH1A1*2 and alcohol dependence, there was a nominally significant difference in ALDH1A1*3 allelic frequency between alcoholics (6 per cent) and controls (0), and the ALDH1A1*3 polymorphism resulted in reduced ALDH1A1 expression. Ehlers et al. replicated this study in a south-west Californian Indian population of alcoholics and controls and observed a protective role for the ALDH1A1*2 polymorphism, whereby individuals with an ALDH1A1*2 allele reported lower rates of alcohol dependence, drank fewer drinks per session and reported a lower maximum number of drinks consumed in a 24-hour period. Conversely, the ALDH1A1*2 alleles were later reported to be associated with increased development of alcohol dependency in Indo-Trinidadians [9].

Two other studies have reported genetic linkage in the 9q21 region with phenotypic measures of alcohol consumption in predominantly Caucasian American samples from the Framingham Heart Study (FHS). Bergen et al. [10] observed nominally significant evidence of linkage in the 9q21.11 region between three loci and a square-root-transformed 'maximum alcohol consumed in a 24-hour period' variable, MAXAPD. Stronger evidence for linkage of MAXAPD on chromosome 9 was detected by Ma et al. [11] with a maximal logarithm of the odds ratio (LOD) score of 2.27 in the region of ALDH1A1.

Due to the recently established role of ALDH1A1 in alcohol-related phenotypes and the dual function of ALDH1A1 in both the alcohol and DA metabolism pathways, we screened eight polymorphisms across the 52.8-kilobase (kb) ALDH1A1 gene for association with alcohol consumption patterns and alcohol dependence in our Finnish population. Alcohol dependence and problem drinking were determined by alcoholism treatment status and the Alcohol Use Disorder Identification Test (AUDIT), a measure of alcohol consumption behaviour.

Materials and methods

Recruitment and assessment of subjects

The Finnish subjects (with Finnish ethnicity) consisted of 104 men enlisted during clinical treatment for alcoholism and a sample of 201 male volunteers from the general population (including alcohol-dependent and control subjects), either recruited through advertisements (n = 129) or drawn from the Finnish Population Register (n = 72) by the National Public Health Institute of Helsinki, Finland. All subjects were from the metropolitan area. Blood was obtained from all subjects and total genomic DNA was extracted using the PUREGENE® DNA Isolation Kit (Gentra, Inc., Minneapolis, MN, USA) according to the manufacturer's instructions.

Hazardous or harmful alcohol consumption behaviour and alcohol-related problems were assessed by the AUDIT. This instrument consists of ten items, which, when responded to in the affirmative, increase the likelihood that the subject is alcohol dependent. The questionnaire measures alcohol intake (quantity and frequency; items 1-3), dependence symptoms (items 4-6) and alcohol-related problems (items 7-10). All items are equally weighted and scored between 0 and 4 and the completed AUDIT is scored by summing the values for each item. The interpretation of the AUDIT, with respect to alcohol use disorders, has been well validated and reviewed by a number of groups [12-15]. Babor et al. [16] distinguished between hazardous and harmful drinking, whereby hazardous alcohol intake is indicative of an increased risk of harmful consequences for the user or others, and harmful alcohol use is defined as consumption behaviour that has resulted in physical or psychological damage. Additionally, Babor and colleagues [14] recommend further evaluation for alcohol dependence for those who score 20 or higher. Accordingly, our study population was subdivided into four groups based upon total AUDIT score and history of alcohol dependence treatment: (1) control drinkers scored ≤12 on the AUDIT but were not alcohol naïve (n = 105); (2) hazardous drinkers scored 13-19 on the AUDIT (excluded from further genetic analyses, n = 32); (3) problem drinkers had an AUDIT score ≥20 (n = 64); and (4) treated alcoholics were recruited from alcohol treatment facilities (n = 104).

Genotyping

Since genotype data were not available from the International HapMap Project public database [17] at the time this study was conceived (July 2003), eight single nucleotide polymorphisms (SNPs) in or near the ADLH1A1 locus were selected from the Applied Biosystems SNP database[18]. SNPs were initially chosen to be evenly distributed across ADLH1A1, including one SNP in both the 5' and 3' untranslated regions (UTRs), with an average inter-marker spacing of 10.7 kb.

All primers, probes and reagents were purchased from Applied Biosystems Inc. (ABI; Foster City, CA, USA). To ensure quality control, DNA samples from cases and controls were randomly distributed across a 384-well plate. SNPs were genotyped using TaqManTM fluorescence 5' exonuclease technology. Each 5 μl experiment contained 25 ng genomic DNA, 1.6X TaqMan assay primer/probe mix, 1X PCR buffer A, 2.5 mM MgCl2, 250 μM deoxyribonucleotide triphosphates (dNTPs) and 0.5 U AmpliTaq Gold polymerase. Thermocycling was performed as recommended by ABI using a Dual 384-Well GeneAmp® PCR System 9700 cycler and cycling conditions as follows: an initial denaturation stage at 95°C for 12 minutes, then 40 cycles of 95°C for 15 seconds and 60°C for one minute, followed by storage at 4°C in the dark. Genotypes were determined on an ABI 7900HT Fast Real-Time PCR System using the allelic discrimination mode.

Statistical analysis

All polymorphisms were first assessed to determine if the observed genotype frequencies were consistent with Hardy-Weinberg proportions using Pearson chi-square tests within Stata (version 8.1; Stata Corporation, College Station, TX, USA). Pair-wise marker-marker linkage disequilibrium (LD) was assessed in Haploview (version 4.0) [19]. To test for association between each SNP and the continuously distributed 'AUDIT score', we performed analysis of variance models (ANOVA) using PROC GLM within the SAS software system (version 8.0; SAS Institute, Cary, NC, USA). Genotype was first tested for general association (no mode of inheritance assumption) using a two-degrees of freedom F-test. Further testing of the SNP effects assuming genetic modes of inheritance for that SNP were then performed. For dominant or recessive genetic transmission models, a single indicator variable (eg Idom is an indicator variable taking the value 1 ifan individual has genotype 1/2 or 2/2, and 0 otherwise) was used. A variable taking the values -1 for genotype 1/1; 0 for genotype 1/2; and 1 for genotype 2/2 was used to test for additive genetic effects.

In addition, odds ratios (OR) and 95 per cent Wald confidence intervals (95 per cent CI) were estimated in Stata 8.1 for both the allelic and genotypic distribution by dividing the dataset into controls (n = 105, AUDIT score ≤12), problem drinkers (n = 64, AUDIT score ≥20) and treated alcoholics (n = 104). Genotype and allele distributions for the problem drinkers and treated alcoholic subjects were each compared with those in the controls as the reference group. Genotypic odds ratios were calculated by collapsing two genotypes (one heterozygous and one homozygous group) into one, and comparing the collapsed genotype count with that of the count of the remaining homozygous genotype.

To test for haplotypic association between the markers and AUDIT score, a generalised linear model for the quantitative trait 'AUDIT score' was run in Qtphase manufacturers. Qtphase is part of the UNPHASED [20] suite of programs for association analysis of multilocus haplotypes (an expectation-maximisation [E-M] algorithm is used to deal with uncertain haplotypes) from unphased genotype data [21]. All chi-square statistics and p-values are presented without correction for multiple testing.

Results

A total of 305 male participants provided blood samples and completed the AUDIT questionnaire. The mean age of the study population was 45.1 ± 11.1 years, ranging from 21 to 75 years. While the average age of the treated alcoholics (44.3 years), problem drinkers (47.5 years) and control groups (44.2 years) at the time of interview did not differ significantly (p ≥ 0.05), the average weight (kg, p = 0.026) but not body mass index (kg/m2, p = 0.19) of the treated alcoholics (77.6 kg, 24.7 kg/m2) and control (81.3 kg, 25.3 kg/m2) groups did. Among the general population (n = 201), 15.9 per cent displayed hazardous patterns of alcohol consumption and 31.8 per cent scored ≥ 20 on the AUDIT and were classified in our analyses as problem drinkers (since they generally experienced both dependence symptoms and alcohol-related problems) and were also included in the problem/dependent group.

The problem drinkers and treated alcoholics had marked differences (p ≤ 0.0001) in their mean responses for all ten AUDIT items compared with the controls (Table 1). Although all groups tended to consume alcoholic beverages on a weekly basis (item 1; 55 per cent of controls and an average of 77 per cent for the problem and dependent drinkers), the majority (82 per cent) of problem drinkers and treated subjects consumed an average of six or more drinks per session (item 2) compared with 13 per cent of controls (on average, the control group imbibed three or four drinks per session). This trend was similar for the frequency with which subjects consumed six or more drinks per session (item 3). While these significant differences in alcohol consumption behaviour (items 1-3) were observed between the control and problem/treated groups, an unambiguous distinction from the controls was seen with respect to those items related to dependence symptoms (items 4-6) and alcohol-related problems (items 7-10). For these seven items, the mean response was ≤0.5 and ≥2.2 for controls and problem/ treated groups, respectively, which on a 'occurs on a monthly or more' basis represents an affirmative response by ≤ 2 per cent and ≥ 63 per cent of the respective group subjects.

Table 1

Alcohol Use Disorder Identification Test (AUDIT) response, mean and standard deviation per item, in the Finnish control drinkers, problem drinkers and treated alcohol-dependent groups, based on AUDIT score and/or alcohol treatment status

Item	Assessment	AUDIT item descriptiona	Controlsb	Problem drinkersc	Treated alcoholicsd
--	Subjects(n)	--	105	64	104
--	Clinically treated for alcoholism	--	0	0	100%
1	Hazardous alcohol consumption	Frequency of drinking days	2.5 ± 1.0	3.4 ± 0.8	3.1 ± 1.0
2	Hazardous alcohol consumption	Average number of drinks per session	1.1 ± 1.1	3.4 ± 10.8	3.4 ± 1.1
3	Hazardous alcohol consumption	Frequency of six or more drinks	1.4 ± 0.8	3.4 ± 0.7	3.0 ± 1.0
4	Dependence symptoms	Impaired control over drinking	0.1 ± 0.3	2.6 ± 1.2	2.3 ± 1.4
5	Dependence symptoms	Increased salience of drinking	0.1 ± 0.3	2.0 ± 1.2	2.3 ± 1.3
6	Dependence symptoms	Morning drinking	0.1 ± 0.3	2.6 ± 1.3	2.7 ± 1.3
7	Alcohol-related problems	Guilt after drinking	0.5 ± 0.6	2.9 ± 1.2	2.7 ± 1.4
8	Alcohol-related problems	Alcohol-related blackouts	0.4 ± 0.5	2.4 ± 1.2	2.4 ± 1.3
9	Alcohol-related problems	Alcohol-related injuries	0.2 ± 0.7	2.2 ± 1.5	2.6 ± 1.5
10	Alcohol-related problems	Others concerned about drinking	0.0 ± 0.4	3.4 ± 0.9	3.6 ± 0.9
--	Hazardous/harmful alcohol usee	AUDIT score (sum of items)	6.3 ± 3.1	28.3 ± 5.4	27.7 ± 8.4

aAll AUDIT items, except questions 9 and 10, are confined to previous year.

bThe Finnish population consists of alcoholic subjects recruited from alcohol treatment facilities (n = 104) and volunteer subjects from the general population that includes control drinkers (n = 105), hazardous drinkers (n = 32) and problem drinkers, as defined by the AUDIT (n = 64). Control drinkers had an AUDIT score ≤12 and had consumed alcohol in the previous year.

cProblem drinkers have an AUDIT score ≥20 but were not recruited from alcohol treatment facilities.

dTreated alcoholics were recruited from alcohol treatment facilities.

eAn AUDIT score ≥12 indicates a pattern of hazardous consumption, while an AUDIT score ≥20 suggests that the subject is alcohol dependent13.

The physical locations of the eight SNPs across the ALDH1A1 gene are schematically presented in Figure 1. Frequencies of the ALDH1A1 genotypes were in general in Hardy-Weinberg equilibrium (HWE), except for rs595958 in intron 1 (p = 0.0008), and are listed in Table 2. Minor allele frequencies within the control sample were similar to those observed in the Centre d'Etude du Polymorphisme Humain (CEPH) sample (data not shown) of the International HapMap Project [22]. The extent of LD between the eight SNPs is presented in Figure 1. LD was strongest between markers rs610529 (intron 8) and rs13959 (exon 3), with one haplotype block structure detected by Haploview [19] between these SNPs. The effective number of independent SNPs was 6.0, as determined by single nucleotide polymorphism spectral decomposition (SNPSpD) [23,24].

Figure 1

A schematic representation of the human . The gene structure of ALDH1A1 is shown with exons numbered and relative exon size denoted by the width of the vertical bars. The eight SNPs analysed in this study are shown in relation to their location across ALDH1A1 and regions of low to high pair-wise LD, as measured by the r2 statistic, are represented by light grey to black shading, respectively.

Table 2

Association analysis of eight cystolic aldehyde dehydrogenase (ALDH1A1) polymorphisms with AUDIT score

Variant	Chromosomal location (bp)a	Gene location	Allelesb		MAFc	Allelic AUDIT score meansd		Genotypic AUDIT score meanse			Genotypic association (p value)		Mode of inheritance
			Allele 1	Allele 2		Allele 1	Allele 2	1/1	1/2	2/2	Generalf	Best modelg
rs348479	74702120	3' UTR	G	T	0.177	17.79	19.66	11.00	19.47	19.33	0.08	0.025	Dominant
rs348472	74710880	Intron 1 1	A	G	0.035	20.67	19.30	28.00	18.90	19.08	0.75	0.44	Dominant
rs610529	74723144	Intron 8	A	G	0.399	20.21	17.85	20.50	19.07	14.65	0.027	0.013	Recessive
rs2288087	74733858	Intron 4	A	T	0.475	18.94	19.71	17.53	19.61	19.34	0.48	0.23	Dominant
rs13959	74735702	Exon 3	A	G	0.472	19.65	19.00	19.34	19.48	17.83	0.62	0.33	Recessive
rs348449	74744802	Intron 2	A	G	0.033	19.18	24.20	18.78	24.44	22.00	0.13	0.043	Dominant
rs595958	74750694	Intron 1	C	T	0.354	19.52	19.01	19.19	19.70	17.44	0.53	0.28	Recessive
rs918836	74777539	5' UTR	C	G	0.278	19.03	19.48	18.87	18.79	19.38	0.91	0.66	Recessive

Abbreviations: AUDIT, Alcohol Use Disorder Identification Test; MAF, minor allele frequency.

abp = base-pair position on NCBI dbSNP Build 128.

bThe allele with the lowest frequency is italicised.

cMinor allele frequency of the italicised allele in the volunteers from the general population (n = 201).

dAllelic AUDIT score means were calculated in Qtphase17.

eGenotypic AUDIT score means were calculated in SAS (SAS Institute, version 8.0, Cary, NC) using a general test of association.

fProbability for general association between each SNP and the AUDIT score using analysis of variance models.

gLowest probability of association assuming a genetic mode of inheritance (given in Table) with respect to allele 2.

We observed nominally significant differences between the genotypic and/or allelic frequencies and mean AUDIT scores for three (rs348479 in the 3' UTR, rs610529 in intron 8 and rs348449 in intron 2) of the eight SNPs genotyped in ALDH1A1. The genotypic distribution and allelic frequencies of the three SNPs are presented in Table 3.

Table 3

Genotypic and allelic analysis of the ALDH1A1 rs348479, rs6l0529 and rs348449 polymorphisms in Finnish controls, problem drinkers and clinically treated alcoholic subjects

SNPa and groupb	Number	Group mean AUDIT score	MAF (%)	Genotype			Genotypic			Allelic
							p-valuec	OR (95% Cl)d	p valuec	OR (95% Cl)d
rs348479				GG	GT	TT		(TT and TG) vs GG
controls	105	6.26	18.1	6.7	22.9	70.4	Ref	Ref	Ref	Ref
problem	64	28.27	18.0	3.1	29.7	67.2	0.4	1.01 (0.57-1.79)	1.0	1.01 (0.57-1.79)
alcoholic	104	27.74	14.9	0	29.8	70.2	0.019	N/A	0.4	1.26 (0.75-2.12)
rs610529				AA	AG	GG		GG vs (GA and AA)
controls	105	6.26	39.5	37.1	46.7	16.2	Ref	Ref	Ref	Ref
problem	64	28.27	37.5	35.9	53.1	10.9	0.6	0.92 (0.58-1.44)	0.7	0.92 (0.58-1.44)
alcoholic	104	27.74	29.8	46.1	48.1	5.8	0.045	0.32 (0.12-0.84)	0.040	0.65 (0.43-0.98)
rs348449				AA	AG	GG		AA vs (AG and GG)
controls	105	6.26	1.0	98.1	1.9	0	Ref	Ref	Ref	Ref
problem	64	28.27	7.0	87.5	10.9	1.6	0.017	7.36 (1.51-35.84)	0.003	7.87 (1.67-37.01)
alcoholic	104	27.74	3.4	93.3	6.7	0	0.09	3.72 (0.75-18.33)	0.1	3.62 (0.74-17.64)

aThe population (n = 305) SNP allele frequencies are as follows: rs348479 G allele frequency = 0.17, rs610529 G allele frequency = 0.37 and rs348449 G allele frequency = 0.03.

bThe three groups are: (1) Controls (AUDIT score ≤ 12); (2) Problem drinkers (AUDIT score ≥20, not treated for alcoholism); and (3) Clinically treated alcoholics.

cGenotype and allele distributions of the problem drinkers and treated alcoholic subjects were each compared with those in the control group using two-sided Pearson's chi-square or Fisher exact tests, respectively.

dOdds ratios (OR) and Wald 95% confidence intervals (95% CI) were determined separately for the allele and genotype counts of problem drinkers and treated alcoholic subjects versus control drinkers. Genotype odds ratios for each SNP were calculated as follows: rs348479 (GG versus [GT and TT]), rs6l0529 ([AA and AG] versus GG) and rs348449 (AA versus [AG and GG]).

Abbreviations: MAF, minor allele frequency; N/A, not applicable; Ref, reference group.

The rs348479 T allele was more frequent in subjects with higher AUDIT scores, where the mean AUDIT scores observed for the rs348479 G and T alleles were 17.8 and 19.7, respectively (Table 2) and the frequency of the T allele was 82 per cent in controls compared with 85 per cent (Table 3) in the treated alcoholic drinkers. In the analysis of the entire sample (n = 305), the association of the T allele with the AUDIT score approached nominal significance (p = 0.08 for general association) and the estimated genotypic AUDIT score means were suggestive of a dominant genetic model (p = 0.025) with respect to the T allele, where the AUDIT score is highest in those individuals with one or two copies of the T allele compared with individuals with two copies of the G allele (mean AUDIT score of 11.0, 19.5 and 19.3 with zero, one and two rs348479 T alleles, respectively). The rs348479 genotypic frequencies were significantly different between the controls and treated alcoholics (p = 0.019); however, case-control odds ratio calculations did not indicate significant allele frequency differences between the controls and treated alcoholics (Table 3). Furthermore, the absence of significant allelic and genotypic differences between the controls and problem drinkers may indicate that rs348479 is involved in risk for alcohol dependence rather than increased or problem alcohol consumption. Nonetheless, the results of this study should also be interpreted in the context of multiple testing, since, correcting for the number of independent SNPs alone, a Bonferroni-corrected p-value of 0.0083 (0.05/6) would be required for an association to be declared significant.

The rs610529 polymorphism in intron 8 was associated (p = 0.027 and p = 0.013 for the general and best model of association, which is based upon a recessive mode of inheritance for the G allele) with the AUDIT score. The genotypic AUDIT score means were 20.5, 19.1 and 14.7 for zero, one and two rs610529 G alleles, respectively (Table 2). Although there were no significant differences in allelic frequencies between control and problem drinkers (OR = 0.92, 95 per cent CI 0.58-1.44), treated alcoholics were approximately 1.5 times less likely to have the G allele (OR = 0.65, 95 per cent CI 0.43-0.98) and three times more likely to be either a homozygous AA or a heterozygous AG (OR = 0.32, 95 per cent CI 0.12-0.84) (Table 3). As with rs348479, these analyses provide some evidence for an association between rs6105249 and risk for alcohol dependence rather than problematic alcohol consumption behaviour.

The rs348449 polymorphism was nominally associated (p = 0.043 for the dominant model of inheritance; p = 0.13 for general association) where the AUDIT score genotypic means were 18.9, 24.4 and 22.0 for zero, one and two G alleles, respectively (Table 2). Significant differences in both the rs348449 genotypic and allelic frequencies were detected between the controls and the problem drinkers (Table 3). The frequency of the minor allele G was 1 per cent in controls compared with 7 per cent in the problem drinkers and 3.4 per cent in treated alcoholics. No G/G homozygotes were observed in the control and treated alcoholics groups. Problem drinkers were 7.8 times more likely to have the rs348449 G allele (OR = 7.87, 95 per cent CI 1.67-37.01, p = 0.003) and 7.4 times more likely to have the GG or GA genotype (p = 0.017). By contrast, the allelic and genotypic distributions in the treated alcoholics and control groups were not significantly different. Therefore, rs348449 may play a role in the consumption of high levels of alcohol (which can lead to dependence symptoms and alcohol-related problems) rather than the development of alcohol dependence.

Since the age of the control subjects ranged between 24 and 69 years, it is possible that a subset of the younger controls could develop alcohol dependence later in life. Therefore, we divided the controls into an older cohort, aged over 40 years (n = 66), and repeated the allelic and genotypic analyses. With each SNP, we observed the same direction of, and similar levels of, association (data not shown).

Two- and three-locus haplotypic association between ALDH1A1 SNP genotypes and the quantitative AUDIT score is graphically presented in Figure 2. Using Qtphase and a sliding window approach, the significance of association increased from p = 0.014 for the single rs610529 SNP analysis to 0.0015 (likelihood ratio statistic [LRS] = 15.47, degrees of freedom [DF] = 3) for rs610529-rs2288087 haplotype analysis, where the two SNPs are in high LD (D' = 0.97, r2 = 0.60). Three common haplotypes (≥1 per cent frequency) were observed, totalling 99.4 per cent of the observed haplotypes. Two common haplotypes, '1-1' (10.4 per cent frequency) and '1-2' (35.9 per cent frequency) had mean AUDIT scores of 23.42 (χ2 = 8.37, p = 0.0038) and 17.64 (χ2 = 7.3, p = 0.0070), respectively. The third common haplotype ('2-1', 53.1 per cent) was not significantly associated (χ2 = 0.4, p = 0.55) with the AUDIT score (mean score = 19.59). The addition of a third neighbouring SNP, rs13959, did not increase the strength of association (p = 0.0024).

Figure 2

Plot of minus log of .

Discussion

Previous linkage and association studies of alcohol consumption or dependence measures have implicated the ALDH1A1 gene [1,2,8,10,11]. Although the majority of alcohol studies related to the aldehyde dehydrogenase family focus on ALDH2, particularly in Asian populations [25,26], recent studies[1,8,9] are shining light on the role that ALDH1A1 may play in several human alcohol-related traits. In the present study, we analysed eight SNPs located across ALDH1A1 (including one SNP in both the 50 and 30 UTRs) for an association with alcohol dependence status (using case-control methodology) and the quantitative AUDIT score in our male Finnish population.

A 1992 study of Finnish drinking habits that utilised the AUDIT questionnaire determined that 22 per cent of males (26 per cent of men aged 30-49) who completed the AUDIT displayed hazardous alcohol consumption (alcohol consumption indicative of an increased risk of harmful consequences for the user or others) and 42 per cent usually drank seven or more units of alcohol per session (item 2), and concluded that in the Finnish culture there was an intoxication-seeking nature to alcohol consumption [27]. In 2000, the 12-month prevalence of alcohol dependence and alcohol use disorders in a representative sample of male Finnish adults (>30 years) were 6.5 per cent and 7.3 per cent, respectively [28]. In this study, 15.9 per cent of subjects that were not recruited from alcoholism treatment facilities had experienced hazardous or harmful alcohol consumption in the previous year (AUDIT score of 12-19), which is similar to the findings of Holmila [27], while 31.8 per cent were categorised as problem drinkers (AUDIT score ≥ 20). It should be noted, however, that alcohol-dependent and control drinkers were both specifically recruited in the advertisement campaign and, therefore, the problem-drinking group may have contained a percentage of alcohol-dependent individuals.

Our findings suggest that both alcohol consumption behaviour and alcohol dependence status are influenced by genetic variation in ALDH1A1. The results from ANOVA analyses demonstrate that two intronic SNPs (rs610529 and rs348449) and one 3' UTR SNP (rs348479) in ALDH1A1 are nominally (p = 0.01-0.04) associated with alcohol consumption behaviour, as described by the AUDIT (Table 2). While rs348479 and rs610529 are associated with risk of alcohol dependence rather than problem drinking, rs348449 in intron 2 influences harmful alcohol consumption behaviour that leads to alcohol-related problems and alcohol dependence symptoms (as determined by odds ratio and χ2 contingency analyses of allele and genotype frequencies in control and problem drinkers).

Although the allelic and genotypic distribution of rs348449 was not significantly different between controls and treated alcohol-dependent subjects, it may play a role in alcohol dependence. If Babor et al. [14] are correct in their assertion that an AUDIT score ≥ 20 indicates alcohol dependence, rs348449 may be associated with the transition from harmful alcohol consumption to alcohol dependence: the mean AUDIT scores were 19.2 and 24.2 for the A and G alleles, respectively, and on average a subject's AUDIT score increased from 18.8 to 24.4 and 22.0 with zero, one and two G alleles, respectively. Additionally, Finns with a G allele were 5.2-fold more likely to have been treated for alcohol dependence or have scored ≥ 20 on the AUDIT. We have reservations with respect to the validity of the observed associations between this SNP and the two alcohol-related phenotypes, however; the low minor allele frequency (MAF) in the three drinking groups (1-7 per cent) and the absence of G/G homozygotes among the control and treated alcohol dependent drinkers means that there is a high risk of false-positive finding with this SNP.

The three SNPs that demonstrated association with alcohol consumption behaviour in this study have not previously been reported in the literature and their functionality is unclear; they are neither coding SNPs nor are they near alternative splice sites. Nonetheless, we believe that at the genotypic and/or allelic level, each SNP may influence transitions in drinking patterns and behaviour in our population, particularly the progression from problem to alcohol-dependent drinker. For example, the genotypic AUDIT score means of the 3' UTR SNP, rs348479, are 11.00 (these subjects are classified as control drinkers) and 19.33 (bordering between harmful and problem drinkers) for G/ G and T/T homozygotes, respectively. Furthermore, the three SNPs are in moderate LD (D' = 0.4 - 0.6), so it may be possible that the SNPs are in strong LD with an as yet ungenotyped polymorphism that is causal with respect to these alcohol-related phenotypes. An important question, with regard to the relevance of the present data, is the functionality of these polymorphisms; this needs to be settled in future studies.

In conclusion, while the associations detected in the present study are not significant after taking into account multiple testing, multiple studies reporting associations between ALDH1A1 and alcohol-related phenotypes have been published recently, and our results provide further evidence for a role of ALDH1A1 in the risk for problem alcohol consumption behaviour and possibly alcohol dependence. With this in mind, our findings need to be replicated in a larger Finnish or Caucasian population and a systematic search for, and genotyping of, additional ALDH1A1 variants that influence alcohol-related phenotypes is warranted.