Abstract
Alcohol dependence is a complex neuropsychiatric disorder. Numerous studies investigated association between MTHFR gene C677T polymorphism and alcohol dependence (AD), but the results of this association remain conflicting. Accordingly, authors conducted a meta-analysis to further investigate such an association. PubMed, Elsevier Science Direct and Springer Link databases were searched for studies on the association between the MTHFR C677T polymorphism and AD. Pooled odds ratio (OR) with 95% confidence interval (CI) was calculated using the fixed- or random-effects model. Statistical analysis was performed with the software program MetaAnayst and MIX.
A total of 11 articles were identified through a search of electronic databases, up to February 28, 2020. The results of the present meta-analysis did not show any association between MTHFR C677T polymorphisms and AD risk (for T vs. C: OR = 1.04, 95% CI = 0.88-1.24; CT vs. CC: OR=1.02, 95%CI= 0.62-1.68; for TT + CT vs. CC: OR = 1.10, 95% CI = 0.94-1.29; for TT vs. CC: OR = 1.01, 95% CI = 0.66-1.51; for TT vs. CT + CC: OR = 0.97, 95% CI = 0.66-1.40). Results of subgroup analysis showed no significant association between MTHFR C677T polymorphism with AD in Asian as well as in Caucasian population. In conclusion, C677T polymorphism is not a risk factor for alcohol dependence.
Introduction
Alcohol dependence or Alcoholism, also regarded as alcohol use disorder (AUD), is a complex and relapsing neuropsychiatric disorder (Koob 2003; Volkow et al., 2009). World Health Organization(WHO) reported that approximately 140 million individuals addicted to alcohol globally, resulting in to 2.5 million death each year (WHO,2011). AD is regarded as a “reward deficiency syndrome” that intemperately affects public health (Comings and Blum,2000; Parry et al., 2011). It has been found to be influenced by both genetic and environmental factors (Prescott and Kendler,1999; Kendler et al., 2007). Exact patho-physiological and molecular mechanism of AD is not known yet. However, molecular genetic studies support that multiple genes determine an individual"s predisposition to AD (Begleiter et al., 1999; Sokolov et al., 2003). Heritability of AD likely plays an important role in its development and is determined to be moderate to high (Goldman et al.,2005; Heath et al.,1997). It was reported frequently that alcohol consumption increased homocysteine (Hcy) concentration i.e hyperhomocysteinaemia (Bleich et al.,2005). However, inconsistent results of the combined effect of both positive and negative association have been reported between alcohol intake and Hcy (Ganji and Kafai 2003). Hyperhomocystenemia is already reported as risk factor for several diseases or disorders including neural tube defects, Alzheimer disease, schizophrenia, pregnancy complications, cardiovascular diseases, noninsulin dependent diabetes and end-stage renal disease as evidenced from several studies (Seshadri et al.,2002).
Homocysteine is a sulfur containing amino acid, several genetic and environmenta ris factor are reported for higher plasma concentration of homocysteine (Gudnason et al.,1998). Homocysteine (Hcy) is synthesized in methionine and folate cycle by demethylation of methionine. 5,10-methylenetetrahydrofolate reductase (MTHFR) enzyme of folate cycle plays an important role in homocysteine metabolism. MTHFR gene is present on chromosome 1p36.3. Numerous single nucleotide polymorphisms (SNP) are known in MTHFR gene like C677T and A1298C etc (Frosst et al. 1995; Goyette et al., 1995). The most clinically important and studies polymorphism is C677T (rs 1801133), in which cytosine (C) is substituted with thymine (T) at 677 nucleotide position and consequently alanine is replaced by Valine in MTHFR enzyme(Ala 222 Val) (Rozen et al., 1997; Chango et al., 2000). The variant MTHFR enzyme is thermolabile with reduced activity (‘∼70%) and it increased the plasma homocysteine concentrations (Frosst et al.,1995). Globally, frequency of mutant T allele varies greatly (Rady et al., 2002; Wilcken et al., 2003; Rai et al., 2010, 2012; Yadav et al., 2018). Yadav et al (2018) have conducted a comprehensive C667T polymorphism study and reported the highest frequency in European populations ranging from 24.1% to 64.3% and, lowest frequency from African population. Several studies revealed association of MTHFR gene C677T polymorphism with AD. However, findings showed inconsistent results (Lutz et al., 2006; Saffroy et al., 2008; Benyamina et al., 2009). To derive a more precise estimation of the relationship, authors performed a meta-analysis.
Methods
Meta-analysis of observational studies in epidemiology (MOOSE) guidelines (Stroup et al., 2006) are followed in present meta-analysis.
Retrieval strategy and selection criteria
Articles were retrieved through Pubmed, Google scholar, Springer Link, and Science Direct databases up to February 28, 2020, using following key words: ‘Methylenetetrahydrofoate reductase’ or ‘MTHFR’ or ‘C677T’ and ‘Alcohol dependence’ or ‘AD’ or ‘Addiction’.
Inclusion and exclusion criteria
Inclusion criteria were following: (1) MTHFR C677T polymorphism and alcohol dependence association was investigated in the study, (2) MTHFR C677T genotype/ allele numbers in alcohol dependence cases and controls were given in the study and (3) sufficient information for calculating the odds ratio (OR) with 95% confidence interval (CI). Major reasons for studies exclusion were as follows: (1) no alcohol dependence cases analyzed, (2) the C677T polymorphism details information missing, and (3) duplicate article.
Data extraction
Name of first author, country name, number of cases and controls, number of genotypes in cases and controls and journal name with full reference from each article were extracted.
Statistical Analysis
All analysis were done according to the method of Rai et al (2014). ORs with 95 % confidence intervals (CIs) were calculated using fixed effect and random effect models (Mantel and Haenszel, 1959; Dersimonian and Laird, 1986). A five genetic models viz. alelle contrast, co-dominant, homozygote, dominant and recessive models were calculated. Heterogeneity was investigated and quantified by I2 statistic (Higgins and Thompson,2002). Chi-squared analysis was used to determine whether the genotype distribution of control group was in Hardy–Weinberg equilibrium or not. Subgroup analysis was conducted by ethnicity. Publication bias was investigated by Egger"s regression intercept test (Egger et al., 1997). P value <0.05 was considered statistically significant. All calculations were done by softwares MIX version 1.7 (Bax et al., 2006) and MetaAnalyst (Wallace et al., 2013) program.
Results
Eligible studies
Following the exclusion criteria, 10 individual case-control studies with a total of 1676 cases and 1594 controls were included into this meta-analysis (Bonsch et al., 2006; Lutz et al.,2006; Lutz et al., 2007; Saffroy et al., 2008; Benyamina et al., 2009; Fabris et al., 2009; Shin et al., 2010; Supic et al.,2010; Singh et al.,2014; Singh et al.,2015). One author (Supic et al.,2010), reported their data in to two categories, we included both set of data as different studies. Hence, total number of included studies in present meta-analysis is eleven (Tabel 1).
Summary Statistics
Overall, eleven studies provided 1676/1594 cases/controls for MTHFR C677T polymorphism. The prevalence of C and T alleles in AD cases was 71.22% and 28.79% respectively. The percentage frequency of TT genotype among cases and controls was 9.43% and 12.98%, respectively whereas prevalence of CT heterozygote among AD cases was 38.72% and 39.21% in controls. The prevalence of CC homozygote among AD cases and controls was 51.85% and 47.80%, respectively. Genotypes were in Hardy-Weinberg equilibrium in all controls. In control group the percentage of C and T allele frequencies was 67.41% and 32.59% respectively (Figure 1).
Meta-analysis
No significant association was observed between the MTHFR C677T polymorphism and the susceptibility to AD in all the genetic models using random effect model (for T vs. C (allele contrast): OR = 1.04, 95% CI = 0.88-1.24; CT vs. CC (co-dominant): OR=1.02, 95%CI= 0.62-1.68; for TT + CT vs. CC (dominant): OR = 1.10, 95% CI = 0.94-1.29; for TT vs. CC (homozygote): OR = 1.01, 95% CI = 0.66-1.51; for TT vs. CT + CC (recessive): OR = 0.97, 95% CI = 0.66-1.40)(Table 2; Figures 2).
A true heterogeneity existed between studies for allele contrast (Pheterogeneity =0.02, Q= 20.64, I2= 51.56%, t2=0.04, z = 0.69), co-dominant genotype (Pheterogeneity <0.0001, Q= 86.64, I2=88.46%, t2=0.61, z = 4.29), homozygote genotype (Pheterogeneity =0.02, Q= 20.93, I2=52.24%, t2=0.24, z = 0.1), and recessive genotype (Pheterogeneity =0.02, Q= 21.00, I2= 52.40%, t2=0.20, z = 0.47) comparisons. The ‘I2’ value of more than 50% shows high level of true heterogeneity.
Subgroup analysis
Out of 11 studies included in the present meta-analysis, 3 studies were carried out in Asian countries, and 8 studies were carried out on Caucasian (Table 1). The subgroup analysis by ethnicity did not reveal any significant association between MTHFR C677T polymorphism and AD in Asian population (T vs. C: OR= 1.16; 95% CI= 0.93-1.44; p= 0.17; I2= 3.1%; Pheterogeneity= 0.65; TT vs. CC: OR= 1.16; 95% CI= 0.62-2.02; p= 0.69; I2= 3.1%; Pheterogeneity= 0.89; and TT+CT vs. CC: OR= 1.26; 95% CI= 0.96-1.67; p= 0.09; I2= 3.1%; Pheterogeneity=0.81); and Caucasian population (T vs. C: OR= 0.99; 95% CI= 0.86-1.14; p= 0.93; I2= 61.75%; Pheterogeneity= 0.01; TT vs. CC: OR= 0.95; 95% CI= 0.70-1.29; p= 0.75; I2= 65.63%; Pheterogeneity= 0.004; and TT+CT vs. CC: OR= 1.03; 95% CI= 0.85-1.25; p= 0.75; I2= 13.93%; Pheterogeneity=0.32)
Publication bias
Symmetrical shape of Funnel plots’ revealed absence of pubication bias. P values of Egger’s test were more than 0.05, also provided statistical evidence for the funnel plots’ symmetry (p= 0.62for T vs. C; p= 0.48 for TT vs CC; p= 0.26 for CT vs. CC; p= 0.48 for TT+AC vs. CC; p= 0.28 for TT vs. CT+CC) (Table 3; Figure 4).
Discussion
MTHFR enzyme is crucial for the methylation of homocysteine in to methionine and individulas with C677T polymorphism the plasma concentration of homocysteine is high. In vivo and in vitro studies has demonstrated that homocysteine has neurotoxic effects especially on dopamine neurons of reward pathway (Bleich et al. 2000). In addition, hyperhomocysteinemia is also reported in AD (Bleich et al 2000; Cravo and Camilo, 2000). According to deficit hypothesis of addiction, C677T polymorphism-dependent alteration of the reward system possibly leads to alcohol addiction.. Further, homovanillic acid (HVA) is a potential indicator of central dopaminergic neuronal activity (Amin et al., 1992) and experimentally, it was demonstrated that hyperhomocystein lowers the eve of HVA in rat stria region (Lee et al., 2005). On the basis of 11 studies providing data on MTHFR C677T genotype and AD risk in two ethnic populations, including over 3,205 subjects, our meta-analysis provides an evidence that TT and CT genotypes or T allele are not associated with AD risk. Hence the present meta-analysis indicated that C677T is a not a risk factor of AD.
Meta-analysis is a statistical tool to combine the information of independent case-control studies with similar target (Ioannidis et al., 2002). Several meta-analysis are published, which evaluated effects of gene polymorphisms in susceptibility of diseases/disorders-cleft lip and palate (Rai, 2017), Glucose-6-phosphate dehydrogenase deficiency (Kumar et al.,2016), down syndrome (Rai,2011; Rai et al., 2017; Rai and Kumar, 2018), male infertility (Rai and Kumar,2017), schizophrenia (Yadav et al., 2015; Rai et al., 2017), obsessive compulsive disorder (Kumar and Rai,2019), depression (Rai,2014,2017), epilepsy (Rai and Kumar,2018), Alzheimers disease (Rai, 2016), esophageal cancer (Kumar and Rai, 2018), prostate cancer (Yadav et al.,2016), breast cancer (Rai et al.,2017), digestive tract cancer (Yadav et al.,2018), ovary cancer (Rai, 2016), endometrial cancer (Kumar et al., 2018), uterine leiomyioma (Kumar and Rai, 2018), and MTRR gene frequency (Yadav et al.,2019).
Several limitations that should be acknowledged like (i)calculated crude Odds ratio, (ii) included the less number of available studies(10 studies) and the limited sample size of each included study, (iii) observed higher between study heterogeneity, (iv) considered only one gene polymorphism and (v) not considered other confounding factors like diet, gender etc. In addition to limitations, current meta-analysis has several strength also such as - higher study power and larger sample size in comparison to individual case control studies, and absence of publication bias etc.
In conclusion, pooled analysis of data from 11separate articles indicates that the MTHFR 677TT genotype is not a risk factor for AD. The results of present meta-analysis should be interpreted with certain cautions due to presence of higher heterogeneity and small number of included studies. Future large-scale, population-based association studies from different regions of the world are required to investigate potential gene–gene and gene–environment interactions involving the MTHFR C677T polymorphism in determining AD risk.
Data Availability
Yes data are availabe