INTRODUCTION

Alcohol, cocaine, heroin, and methamphetamine addiction constitute major public health issues. The cost of dealing with drug abuse has grown to approximately one trillion dollars per year in the United States (Califano, 2007). An estimated 35.3 million Americans aged 12 years and older reported having used cocaine, and 2.4 million Americans were current users (according to the National Survey on Drug Use and Health). Despite a large amount of research on the etiopathology of substance use disorders, there remains a great deal of work to be done to identify specific genetic and environmental risk factors. Alcohol and drug addictions are familial and share some common genetic factors (Fu et al, 2002; True et al, 1999; Xian et al, 2008). Genetic association studies (Bierut et al, 2007; Edenberg et al, 2010; Zlojutro et al, 2011), eg, Collaborative Study on the Genetics of Alcoholism, have reported that a wide number of candidate genes contribute to the risk for alcoholism and other substance use disorders.

One such gene is the serotonin (5-hydroxytryptamine) transporter gene (SLC6A4 or 5-HTT). SLC6A4 is one of the most studied candidate genes and is the focus of this report. The serotonin transporter protein (SERT) is the presynaptic neuronal reuptake site for serotonin and a site of action for several drugs with central nervous system effects. SERT is a 630 amino-acid protein with 12 transmembrane domains and is a member of the Na+/Cl−-dependent transporter family (Ramamoorthy et al, 1993). The SLC6A4 gene, spanning 37 809 base pairs (bp) and consisting of 14 exons, is located on 17q11.1-q12(Gelernter et al, 1995). SLC6A4 has been linked to alcohol, heroin, cocaine, and methamphetamine dependence and abuse, such as linkage studies (Gelernter et al, 2006; Glatt et al, 2006) have identified the long arm of chromosome 17 as a susceptibility region to heroin dependence. Two important SLC6A4 variable number of tandem repeat (VNTR) polymorphisms, 5-HTTLPR and STin2, have been widely studied. The 5-HTTLPR polymorphism is located in the 5′ regulatory region and the two most common alleles are the ‘long’ (L) 16-repeat and ‘short’ (S) 14-repeat alleles (other alleles are observed in various populations (Gelernter et al, 1997)). This variation correlates with differential expression of the SERT protein in cell lines (Lesch et al, 1996). Another VNTR polymorphism, STin2, consists of the presence of a variable number of repeats (usually 9, 10, or 12 repeats) of a 17 bp segment that maps to intron 2. In vitro studies suggested that the short variants of STin2 decreased promoter activity and further the mRNA and protein concentration, therefore, cell lines carrying short variants showed a decreased serotonin uptake efficiency (Lesch et al, 1994). Study showed 5-HTTLPR may modulate the gene expression possibly via a combined effect with the STin2 polymorphism (Hranilovic et al, 2004). STin2 and 5-HTTLPR have been also highlighted to be associated with several other disorders, eg, affective disorders (Collier et al, 1996), obsessive-compulsive disorder (McDougle et al, 1998), suicidal behavior (Li and He, 2007), attention-deficit hyperactivity disorder (Gizer et al, 2009), amygdala activation (Munafo et al, 2008; Murphy et al, 2013), stress and depression (Karg et al, 2011), and antidepressant response (Kato and Serretti, 2010).

To date, there have been a number of case–control and family-based association studies that examined the relative risks of one or both of 5-HTTLPR and STin2 in alcohol or drug dependence. Some of these studies reported positive findings (Konishi et al, 2004; Patkar et al, 2001; Wu et al, 2008), but many others found no evidence of association (Supplementary Table 1). There are some reasons for lack of replication, such as insufficient sample size, inadequate statistical power, and failure to control for population variations. The aim of this study was to combine all of the available genotype data from prior case–control and family-based association studies in a multi-cultural meta-analysis of the SLC6A4 5-HTTLPR and STin2 polymorphisms with alcohol and drug (heroin, cocaine, and methamphetamine) abuse in the European, Asian, African, and Mexican populations.

MATERIALS AND METHODS

Literature Search

Published reports were selected from Scopus, PubMed, and Chinese Academic Journals database with keywords ‘SLC6A4', ‘5-HTT’, ‘serotonin transporter’, ‘association’, ‘associated’, ‘drug’, ‘substance’, ‘alcoholism’, ‘alcohol’, ‘alcoholics’, ‘heroin’, ‘cocaine’, ‘opiate’, ‘opioid’, ‘methamphetamine’, ‘morphine’, ‘opium’ and the specific names or abbreviations of the gene (ie, SERT'). Both English and Chinese keywords were used in searching the Chinese academic journals. All references cited in these studies and in published reviews were examined in order to identify additional works not indexed by the databases. The analyzed data cover all identified English and Chinese publications up to July 2012.

Inclusion Criteria

Eligible studies had to meet all of the following criteria: they (i) were published in peer-reviewed journals; (ii) contained original and independent data; (ii) presented sufficient samples to calculate the OR with confidence interval (CI) and P-value; (iii) were association studies investigating one or two of the polymorphisms using either case–control or family-based approach; (iv) described or referenced appropriate genotyping methods, primers, machines, or protocols; (v) investigated one or more of the following: alcohol, heroin, cocaine, methamphetamine, or more generally drug dependence (two studies (Hallikainen et al, 1999; Shin et al, 2009) included both alcohol dependence subjects and alcohol abusers and one study (Li et al, 2002) included heroin abusers); and (vi) used unrelated individuals with no explicit description of any of these disorders (some studies recruited healthy normal subjects, whereas other studies used random or general population) as controls for case–control study.

Phenotype Inclusion Criteria

We included studies that diagnosed the patients according to the World Health Organization’s International Statistical Classification of Diseases and Related Health Problems (ICD; World Health Organization) or American Psychiatric Association Diagnostic and Statistical Manual of Mental Disorder (DSM) system. We also included one study using Michigan Alcoholism Screening Test (Selzer, 1971), one study using Feighner Diagnostic Criteria (Feighner et al, 1972), and one study using Fagerstrom Test for Nicotine Dependence (Heatherton et al, 1991). The procedure of ‘extended-quality score’ (Li et al, 2006), which scores each paper categorizing it as having ‘high’, ‘median’, or ‘poor’ quality, was applied to assist the assessment of quality of the association studies. Authors were contacted in cases where we determined it would be useful to have additional information regarding their studies.

Statistical Analyses

Studies were classified according to design as case–control or family-based study, and the latter further subdivided according to statistical methodology into haplotype relative risk (HRR), pedigree disequilibrium test (PDT), and transmission disequilibrium test (TDT). Studies were also subdivided by ethnicity, ie, European ancestries, Asian ancestries, African ancestries, and Mexican ancestries. For studies that contained data from multiple populations, each was considered effectively as an independent study. Data from the case–control, HRR, and PDT studies were summarized by two-by-two tables and TDT studies were summarized by two-by-one tables. The two types of studies were statistically combined by the method described in the previous studies (Li and He, 2007; Lohmueller et al, 2003) to join population-based and family-based studies into a single meta-analysis.

From each table a log-OR and its sampling variance were calculated. The Cochran’s χ2-based Q statistic test was computed in order to assess heterogeneity to ensure that each group of studies was suitable for meta-analysis. Where heterogeneity was found, the random effects model, which yields a wider CI, was adopted; otherwise, the fixed effects model was adopted. Heterogeneity Q tests were also performed for differences in OR between subject ethnicities or phenotypes (eg, Europeans vs Asians or alcohol abuse vs heroin abuse). The Egger’s funnel plot asymmetry (Egger et al, 1997) was used to assess evidence for publication bias. The test uses a linear regression approach to measure funnel plot asymmetry on the natural logarithm of the OR. The larger the deviation of each study from the funnel curve, the more pronounced the asymmetry. Results from small studies will scatter widely at the bottom of the graph, with the spread narrowing among larger studies. The significance of the intercept was evaluated using the T-test. Begg and Mazumdar rank correlation (Begg and Mazumdar, 1994) was also employed to evaluate potential bias where the P-value for Kendall’s tau was computed. The ‘Duval and Tweedie’s Trim and Fill’ procedure (Duval and Tweedie, 2000) was adopted to impute the number of potentially missing studies if significant publication bias was found. The Trim and Fill procedure imputes the missing studies, adds them to the analysis, and then re-computes the adjusted overall effect size.

ORs were pooled using the method of DerSimonian and Laird (DerSimonian and Laird, 1986), and 95% CIs were constructed using Woolf’s method (Woolf, 1955). The significance of the overall OR was determined using the Z-test. To measure sensitivity of our analysis results, each study was removed in turn from the total, and the remainder then reanalyzed. This procedure was used to ensure that no individual study was entirely responsible for the combined results. Retrospective analysis was performed to better understand the potential effect of the year of publication upon the results. The type I error rate was set at 0.05. The tests were two-tailed. Haplotype construction, counting, and linkage disequilibrium (LD) block defining over a broader genomic region that include SLC6A4 was performed separately using HapMap samples. The multiallelic D′ and maximum likelihood haplotype blocks were calculated using the methods described in our previous study (Li et al, 2011). Each ethnic population and each phenotype (alcohol, cocaine, heroin, and methamphetamine dependence and abuse) as well as differently combined sub-groups were analyzed.

RESULTS

The combined search yielded 2394 references. After discarding overlapping references and those which clearly did not meet the inclusion criteria, 66 studies remained. These studies were then filtered to ensure conformity with the inclusion criteria. Four studies (Foley et al, 2004; Galeeva et al, 2002; Pastorelli et al, 2001) were excluded because no diagnostic criteria was described explicitly; one study (Munafo et al, 2005) because it was investigating social ‘drinkers’ rather than alcohol-dependent or abusing subjects; one study (Rasmussen et al, 2009) because it investigated alcohol and cigarette consumption rather than dependence or abuse; three studies (Budde et al, 2010; Mingione et al, 2012; Thompson et al, 2010) because no matched control data were described; and two studies (Wang et al, 2011a; Yang et al, 2012) because the genotype data were same as those in two other studies (Deng et al, 2008; Wang et al, 2011b), respectively. In the end, 55 studies (from 51 references; each of the three references (Gelernter et al, 1998; Kranzler et al, 2002; Nellissery et al, 2003) included two independent ethnic populations, and one reference (Saiz et al, 2009) included data for both alcohol and heroin dependence) composed of 50 case–control studies (Choi et al, 2006; Chu et al, 2010; Deng et al, 2008; Drago et al, 2009; Ezaki et al, 2008; Gelernter et al, 1998; Gelernter et al, 1997; Gerra et al, 2004; Gokturk et al, 2008; Gorwood et al, 2000; Grochans et al, 2011; Hallikainen et al, 1999; Hammoumi et al, 1999; Hong et al, 2003; Ishiguro et al, 1999; Johann et al, 2003; Kohnke et al, 2006; Konishi et al, 2004; Kotler et al, 1999; Kranzler et al, 2002; Lee et al, 2009; Li et al, 2002; Marques et al, 2006; Matsushita et al, 2001; Mokrovic et al, 2008; Namkoong et al, 2008; Nellissery et al, 2003; Parsian and Cloninger, 2001; Patkar et al, 2002; Patkar et al, 2001; Patkar et al, 2004; Philibert et al, 2008; Preuss et al, 2001; Reese et al, 2010; Saiz et al, 2008; Saiz et al, 2009; Sander et al, 1998; Sander et al, 1997; Shin et al, 2009; Stoltenberg et al, 2002; Tan et al, 1999; Thompson et al, 2000; Wang et al, 2012; Wang et al, 2011b; Wu et al, 2008; Yang et al, 2012), three TDT studies (Edenberg et al, 1998; Lichtermann et al, 2000; Samochowiec et al, 2006), one PDT study (Dick et al, 2007), and one HRR study (Hill et al, 2002), met our criteria for inclusion. These studies included 32 studies for European populations (Dick et al, 2007; Drago et al, 2009; Edenberg et al, 1998; Gelernter et al, 1998; Gelernter et al, 1997; Gerra et al, 2004; Gokturk et al, 2008; Gorwood et al, 2000; Grochans et al, 2011; Hallikainen et al, 1999; Hammoumi et al, 1999; Hill et al, 2002; Johann et al, 2003; Kohnke et al, 2006; Kotler et al, 1999; Kranzler et al, 2002; Lichtermann et al, 2000; Marques et al, 2006; Mokrovic et al, 2008; Nellissery et al, 2003; Parsian and Cloninger, 2001; Philibert et al, 2008; Preuss et al, 2001; Reese et al, 2010; Saiz et al, 2008; Saiz et al, 2009; Samochowiec et al, 2006; Sander et al, 1998; Sander et al, 1997; Stoltenberg et al, 2002; Thompson et al, 2000); 16 studies for Asian populations (Choi et al, 2006; Chu et al, 2010; Deng et al, 2008; Ezaki et al, 2008; Hong et al, 2003; Ishiguro et al, 1999; Lee et al, 2009; Li et al, 2002; Matsushita et al, 2001; Namkoong et al, 2008; Shin et al, 2009; Tan et al, 1999; Wang et al, 2012; Wang et al, 2011b; Wu et al, 2008; Yang et al, 2012), 6 studies for African Americans (Gelernter et al, 1998; Kranzler et al, 2002; Nellissery et al, 2003; Patkar et al, 2002; Patkar et al, 2001; Patkar et al, 2004), and 1 for Mexican Americans (Konishi et al, 2004). Among the 55 studies, eight studies (Deng et al, 2008; Gerra et al, 2004; Kotler et al, 1999; Li et al, 2002; Saiz et al, 2008; Saiz et al, 2009; Tan et al, 1999; Yang et al, 2012) investigated heroin dependence or abuse; three studies (Patkar et al, 2002; Patkar et al, 2001; Patkar et al, 2004) investigated cocaine dependence; two studies (Ezaki et al, 2008; Hong et al, 2003) investigated methamphetamine dependence; one study (Chu et al, 2010) investigated nicotine dependence; three studies (Gelernter et al, 1998; Gokturk et al, 2008) investigated alcohol dependence, drug dependence, or both; and the other 38 studies (Choi et al, 2006; Dick et al, 2007; Drago et al, 2009; Edenberg et al, 1998; Gelernter et al, 1997; Gorwood et al, 2000; Grochans et al, 2011; Hallikainen et al, 1999; Hammoumi et al, 1999; Hill et al, 2002; Ishiguro et al, 1999; Johann et al, 2003; Kohnke et al, 2006; Konishi et al, 2004; Kranzler et al, 2002; Lee et al, 2009; Lichtermann et al, 2000; Marques et al, 2006; Matsushita et al, 2001; Mokrovic et al, 2008; Namkoong et al, 2008; Nellissery et al, 2003; Parsian and Cloninger, 2001; Philibert et al, 2008; Preuss et al, 2001; Reese et al, 2010; Saiz et al, 2009; Samochowiec et al, 2006; Sander et al, 1998; Sander et al, 1997; Shin et al, 2009; Stoltenberg et al, 2002; Thompson et al, 2000; Wang et al, 2012; Wang et al, 2011b; Wu et al, 2008) investigated alcohol dependence or abuse. Among them, eight studies investigated antisocial alcoholism. For the study by Wu et al (2008), some new data, unavailable in the published paper, were provided by the authors. These studies included 7999 cases, 8264 controls, and 676 families or parent-offspring trios (Supplementary Table 1). The flow chart of literature search strategy is shown in Figure 1. The results for each polymorphism are detailed below.

Figure 1
figure 1

Flow chart of literature search strategy.

PowerPoint slide

5-HTTLPR

The frequency of the long variant allele (L) varied widely across the populations, high in European normal populations 57% (48–71%) and patients 55% (46–67%), but low in Asian normal populations 27% (15–71%) and patients 25% (17–63%). Figure 2 shows the average allele frequencies for the four populations. The frequencies were consistent with those observed in our previous study (Li and He, 2007), in which the L allele were 57% and 27% in European and Asian normal populations, respectively. Of the 51 studies included for this polymorphism, 30 studies showed lower frequency in cases than in controls (or less transmissions of the L allele in families), regardless of ethnicity and sample size (Supplementary Table 2).

Figure 2
figure 2

Average allele frequencies of the 5-HTTLPR ‘L’ allele.

PowerPoint slide

Alcohol dependence/abuse

The combined studies of alcohol dependence and abuse showed an overall allelic P-value of 0.02 (OR=0.91 (0.84, 0.99)) under the random effects model. Evidence of significant association was also found in European (P=0.048) and combined European, Asian, and Mexican populations (P=0.019). The overall P-value was still significant (P=0.037 and OR=0.91 (0.84, 0.99)) after the studies investigating alcohol abusers were excluded from the meta-analysis (only the subjects described as alcohol dependence were analyzed). The dominant model showed evidence of more significant associations, eg, P=0.009 and OR=0.83 (0.72, 0.95) in the combined European, African, and Mexican populations (Table 1).

Table 1 Results of the Overall and Subgroup Studies for the 5-HTTLPR Polymorphism

We also tested the association between 5-HTTLPR and type II alcoholism with antisocial behavior and the association between 5-HTTLPR and severe alcoholics (eg, in some of the included studies, the patients were described specifically with severe withdrawal symptoms, delirium tremens, and (or) seizure). The meta-analysis showed significant association with severe alcoholics in the combined European and Asian populations (P=0.007) under the dominant model. Significant association was also found in European populations (P=0.0058). However, the meta-analysis of the studies investigating type II alcoholism (or alcoholics) revealed no evidence of significant association.

Heroin dependence/abuse

Significant association was found in the combined studies of heroin dependence and abuse with an allelic P-value of 0.02, which was more significant in European populations with P=0.0089 and OR=0.82 (0.7, 0.95). The combined studies of heroin dependence revealed more significant results, eg, in the European populations the allelic P-value was 0.0009 (OR=0.73 (0.61, 0.88)), the P-values being 0.005 and 0.013 under the recessive and dominant models, respectively (Table 1).

Cocaine dependence and methamphetamine dependence

The studies of cocaine dependence and those of methamphetamine dependence only investigated African Americans and Asians, respectively, and the meta-analytic results showed evidence of significant associations with allelic P-values of 0.018 (OR=1.38 (1.06, 1.81)) and 0.04 (OR=0.75 (0.57, 0.99)), respectively, which were also significant under the recessive model (ORs=1.57 (1.05, 2.35) and 0.46 (0.23, 0.92), respectively).

Combined studies

When all the studies of different substance were combined, the allelic analysis showed evidence of significant association (overall P-value=0.02 and OR=0.92 (0.86, 0.99)) under the random effects model due to evidence of heterogeneity between studies (Table 1). The meta-analyses also showed associations in each major population, eg, in the combined European (P=0.03 and OR=0.91 (0.84, 0.99)), European and Asian (P=0.0047 and OR=0.91 (0.85, 0.97)), and non-African (European, Asian, and Mexican; P=0.002 and OR=0.9 (0.84, 0.96)) populations. The African populations also showed significant with a P-value of 0.0028 but in the opposite direction (OR=1.41 (1.13, 1.78)). The recessive model (LL vs LS plus SS) also revealed evidence of associations, eg, in the combined non-African (P=0.018 and OR=0.89 (0.81, 0.98)) and African (P=0.029 and OR=1.57 (1.05, 2.35)) populations. The results under the dominant model (LL plus LS vs SS) showed evidence of stronger association, for example, the overall P-value was 0.015 with OR of 0.87 (0.78, 0.97), and it was more significant in the combined European (P=0.0069 and OR=0.84 (0.74, 0.95)), European and Asian (P=0.002 and OR=0.87 (0.79, 0.95)), and non-African (P=0.0006 and OR=0.86 (0.78, 0.94)) populations.

Diagnosis criteria and control selection

When the studies that employed either the DSM or the ICD system for diagnosis were meta-analyzed, the results showed consistent evidence of significant association. For instance, the combined European and Asian populations showed evidence of significant association with substance use disorders (eg, P=0.0063 and OR=0.88 (0.8, 0.96) under the dominant model). The control subjects had no explicit description of alcohol or drug dependence or abuse. However, some studies explicitly described their controls as ‘healthy’ or ‘normal’ subjects. We also analyzed these ‘super controls’ separately. Evidence of significant association was also found, eg, the allelic P-value was 0.004 (OR=0.82 (0.72, 0.94)) in the non-European populations.

Between-group heterogeneity

There was no evidence of significant heterogeneity between Asian studies and European studies, between Asian studies and the others or between Chinese studies and the others, and between each pair of alcohol, heroin, and methamphetamine dependence/abuse for either the allelic or genotypic analyses (P(Q)>0.1). However, heterogeneity was observed for cocaine dependence (Table 2), which might be partially due to small sample size (all the subjects were Africans). The forest plots of the 5-HTTLPR polymorphism are shown in Figure 3 and supplementary Figure 1 for the allelic analysis and dominant model, respectively.

Table 2 Results of Heterogeneity Estimation Based on Ethnicities and Phenotypes
Figure 3
figure 3

Forest plots of ln(OR) with 95% CI for the 5-HTTLPR allelic analysis. Black squares indicate the ln(OR), with the size of the square inversely proportional to its variance, and horizontal lines represent the 95% CIs. The pooled results are indicated by the unshaded black diamond. For the results of meta-analysis, only the subgroups with P<0.05 are shown.

PowerPoint slide

STin2 VNTR

The 10-allele was high in European normal populations 34.6% (25–54%) and patients 35.8% (27–46%), but had an exceedingly low frequency in Asian normal populations 8% (5–10%) and patients 10% (9–15%), which were consistent with the frequencies that we reported previously (35 and 9% on average, respectively; Li and He, 2007). No evidence of significant association was found for alcohol dependence/abuse. However, the European studies of heroin dependence and abuse showed week association (P=0.02 for the 10/12 genotype). When all the studies of different phenotypes were combined, evidence of significant association was found in the combined Asian populations (P=0.009 and OR=1.47 (1.1, 1.95) for the 10-allele and P=0.02 and OR=0.71 (0.54, 0.95) for the 12-allele). The results are shown in Supplementary Table 3.

Publication Bias and Fail-safe Analyses

In the present meta-analysis, no evidence of significant publication bias was found in the meta-analyses of alcohol dependence/abuse, heroin dependence/abuse, combined drugs, or all the combined studies. The P-values were >0.05 for all these tests based on both Egger’s regression intercept and Begg’s rank correlation. For the 5-HTTLPR polymorphism, the classic fail-safe analysis showed that at least 38 and 3 assumed nonsignificant association studies would be required to bring the P-values to >0.05 for alcohol and heroin dependence/abuse, respectively. When the phenotypes were combined, the associations showed stronger: at least 42 and 49 assumed nonsignificant association studies would be required to bring the P-values to >0.05 for the allelic analysis and dominant model, respectively; for the meta-analysis of non-African (European, Asian, and Mexican) populations at least 114 and 71 assumed nonsignificant studies would be required to bring the P-value to >0.05 for the allelic analysis and dominant model, respectively. The results further supported the significant associations detected in the meta-analyses.

The funnel plots of the 5-HTTLPR studies of alcohol dependence/abuse, heroin dependence/abuse, and combined drugs are shown in Supplementary Figures 2-4, respectively; the plots of the combined studies (combined phenotypes) of the non-African populations are shown for the allelic analysis and dominant model in Supplementary Figures 5 and 6, respectively; the plots of all the combined studies are shown for the allelic analysis and dominant model in Supplementary Figures 7 and 8, respectively.

Sensitivity Analyses

The results of sensitivity analysis of 5-HTTLPR showed that no individual study included in the meta-analyses biased the significant association of heroin dependence. For instance, the studies of heroin dependence (European populations) showed that the P-values were never >0.0082 in the allelic analysis, regardless of the data set removed (not shown). Other major findings also showed consistency with the P-values <0.05. For example, when the phenotypes were combined, the studies of non-African populations showed consistency, regardless of the data set removed, with the P-values never >0.0064 and never >0.0062 for the allelic analysis and dominant model (Supplementary Tables 4 and 5), respectively.

Retrospective Analyses

The analyses in retrospect of 5-HTTLPR based on publication years showed that the cumulative results, as represented by the asymptote lines on the plots, have tended to be stable since 2006 for the meta-analysis. The plots of the alcohol dependence/abuse are shown for the allelic analysis and dominant model in Supplementary Figures 9 and 10, respectively. For an overall trend, the plots of the combined phenotypes of the European, Asian, and Mexican populations are shown for the allelic analysis and dominant model in Figure 4 and Supplementary Figure 11, respectively.

Figure 4
figure 4

Retrospective analysis for the 5-HTTLPR allelic analysis in non-African populations (combined European, Asian, and Mexican populations). Analysis in retrospect was based on publication year since 1997.

PowerPoint slide

LD and Haplotype Structure Analyses

The STin2 polymorphism was in a large haplotype block structure, whereas 5-HTTLPR was located outside of this block, in the gap between this haplotype block and an up-stream block. The plots are shown for the European, Asian, and African populations in Supplementary Figures 12–14, respectively. These structures appears to be consistent with the current results, and implies that the association of 5-HTTLPR may not be due to strong LD with a very close polymorphism. The same findings were observed in two other meta-analysis studies (Fan and Sklar, 2005; Li and He, 2007) dealing with between the two polymorphisms and schizophrenia and suicidal behavior, respectively, although the STin2 polymorphism was shown significant in the study by Fan and Sklar (2005).

DISCUSSION

Alcohol and drug dependence and abuse are multifactorial disorders, and the genetic contribution to vulnerability to develop the disorders is 40–70%, suggesting a complex inheritance mode in which multiple genes and polymorphisms exert a small effect (Gelernter and Kranzler, 2009; Kendler et al, 2007; Uhl et al, 2008). The meta-analyses found the associations between 5-HTTLPR and alcohol, heroin, cocaine, and methamphetamine dependence and abuse. For example, the smallest P-values were 0.0058 with OR=0.54 (0.35, 0.84); 0.0024 with OR=0.77 (0.66, 0.91); 0.018 with OR=1.38 (1.06, 1.81); and 0.028 with OR=0.46 (0.23, 0.92) for alcohol, heroin, cocaine, and methamphetamine dependence/abuse, respectively. When all the phenotypes are combined, the P-value was 0.0006 and OR was 0.86 (0.78, 0.94) in the combined European, Asian, and Mexican populations, whereas P was 0.0028 and OR was 1.41 (1.13, 1.78) in the African populations regarding the ‘L’ allele. Evidence of significant association was also observed in additional subgroup analyses regarding differently combined substance and populations. The effect sizes were comparable among the European, Asian, and Mexican populations, however, the risk ‘S’ allele was significantly more frequent in Asians (73%) than in Europeans (43%) and Mexicans (53%). Based on the between-phenotype heterogeneity analysis, the opposite directions of risk allele of African population vs non-African populations (also explained by ethnic heterogeneity) might be driven by the opposite directions of cocaine dependence vs other substance. The latter could be due to small sample size.

The individual association studies performed by different research groups have contradictory results, the similar phenomenon existing in a previous meta-analysis between SLC6A4 and suicidal behavior (Li and He, 2007). As shown in Figure 2, the 5-HTTLPR allele frequencies vary significantly in different ethnic populations, thus, a difference in sampling methods could differentiate the results. The discrepancy may also be due to insufficient sample size and low statistical power of individual study. In this meta-analysis, the random effects model, which yields larger P-values and wider CIs than the fixed effect model, was applied when heterogeneity was found. Evidence of significant association of 5-HTTLPR was identified with alcohol and drug dependence in the overall and subgroup analyses.

Alcohol and drug dependence are often comorbid with psychiatric disorders or behavior problems. For instance, type I alcoholism (Cloninger et al, 1981) was found to have both environmental and genetic risk factors; and type II alcoholism, the severe form of alcoholism, was found to have a more emphasized genetic etiological factor (Sigvardsson et al, 1996). According to Cloninger’s neurogenetic tripartite theory of personality, serotonin was hypothesized to be the major neuromodulator of harm avoidance (Cloninger, 1987). However, this meta-analysis found no evidence of significant association of SLC6A4 with antisocial behavior in alcoholism partially because of insufficient data published.

Compared with previous meta-analyses of alcohol dependence (Feinn et al, 2005; McHugh et al, 2010), which included 17 and 22 studies, respectively, mainly from European populations, and reported weak or marginal association (eg, P=0.03), the present study is a comprehensive meta-analysis combining (and also separately analyzing) alcohol, heroin, cocaine, and methamphetamine dependence and abuse as well as European, Asian, African, and Mexican populations from 55 case–control and family-based studies by using systematic approaches. However, there are some caveats in this study, for example, stress is a major risk factor in addiction and SLC6A4 modulates stress response and emotionality, but there are no stress exposure data available for meta-analysis. For future studies, it may be interesting to investigate other polymorphisms on SLC6A4 or nearby genes (eg, BLMH), including the noncoding regions as those polymorphisms may influence the gene functions according to the Encyclopedia of DNA Elements (ENCODE) project (Bernstein et al, 2012).

As the first association report between this polymorphism and affective disorders (Collier et al, 1996), 5-HTTLPR has been widely studied with a great number of neuropsychiatric disorders. The SLC6A4 gene protein transports the neurotransmitter serotonin from synaptic spaces into presynaptic neurons. 5-HTTLPR is a noncoding polymorphism that impacts its gene transcription. Various findings support that 5-HTTLPR has important role in the pathogenesis and/or etiology of brain diseases and psychiatric disorders. For example, 5-HTTLPR affects the rate of serotonin uptake and might influence gray matter in anterior cingulate brain region (Pezawas et al, 2005); and the short ‘S’ allele might drive amygdala hyper-reactivity (Hariri et al, 2005). Because of the fundamental roles, 5-HTTLPR is expected to have ‘pleiotropy effect’, ie, a same mutation allele affects multiple-related diseases and traits.

To conclude, our meta-analysis using existing genotype data supports that the association of SLC6A4 5-HTTLPR varies depending on substances (alcohol, heroin, cocaine, and methamphetamine). The ‘S’ allele was the risk allele in the European, Asian, and Mexican populations, whereas the ‘L’ allele was the risk allele in the African populations. The effect sizes were comparable in the European, Asian, and Mexican populations, but the risk ‘S’ allele was more frequent in Asians than Europeans or Mexicans. The opposite directions of the African populations might be driven by the opposite directions of cocaine dependence. Further studies using larger sample size are warranted.

Electronic-database information

Accession Numbers and URLs for data in this article are as follows:

GenBank, http://www.ncbi.nlm.nih.gov/Genbank/ for genomic structure of SLC6A4;

Online Mendelian Inheritance in Man (OMIM), http://www.ncbi.nlm.nih.gov/Omim for SLC6A4; Genotype data, http://www.hapmap.org/ for SLC6A4; Genome data, http://genome.ucsc.edu/ for SLC6A4.