Next Article in Journal
Adsorption of Trace Estrogens in Ultrapure and Wastewater Treatment Plant Effluent by Magnetic Graphene Oxide
Previous Article in Journal
Land Use Regression Modelling of Outdoor NO2 and PM2.5 Concentrations in Three Low Income Areas in the Western Cape Province, South Africa
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJERPH
Volume 15
Issue 7
10.3390/ijerph15071453
ijerph-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Meta-Prediction of Methylenetetrahydrofolate-Reductase Polymorphisms and Air Pollution Increased the Risk of Ischemic Heart Diseases Worldwide

by
Zhao-Feng Chen
1,
Lufei Young
2,
Chong Ho Yu
3 and
S. Pamela K. Shiao
2,*
1
Department of Nursing, Yuanpei University of Medical Technology, Hsinchu 30015, Taiwan
2
College of Nursing, Augusta University, Augusta, GA 30912, USA
3
Department of Psychology, Azusa Pacific University, Azusa, CA 91702, USA
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2018, 15(7), 1453; https://doi.org/10.3390/ijerph15071453
Submission received: 29 May 2018 / Revised: 1 July 2018 / Accepted: 7 July 2018 / Published: 10 July 2018
Browse Figures
Versions Notes

Abstract

:
Ischemic heart disease (IHD) is among the leading causes of death worldwide. Methylenetetrahydrofolate reductase (MTHFR) polymorphisms have been associated with IHD risk, but the findings presented with heterogeneity. The purpose of the present meta-analysis was to provide an updated evaluation by integrating machine-learning based analytics to examine the potential source of heterogeneity on the associations between MTHFR polymorphisms and the risk of various subtypes of IHD, as well as the possible impact of air pollution on MTHFR polymorphisms and IHD risks. A comprehensive search of various databases was conducted to locate 123 studies (29,697 cases and 31,028 controls) for MTHFR C677T, and 18 studies (7158 cases and 5482 controls) for MTHFR A1298C. Overall, MTHFR 677 polymorphisms were risks for IHD (TT: Risk ratio (RR) = 1.23, p < 0.0001; CT: RR = 1.04, p = 0.0028, and TT plus CT: RR = 1.09, p < 0.0001). In contrast, MTHFR 677 CC wildtype was protective against IHD (RR = 0.91, p < 0.00001) for overall populations. Three countries with elevated IHD risks from MTHFR C677T polymorphism with RR >2 included India, Turkey, and Tunisia. Meta-predictive analysis revealed that increased air pollution was associated with increased MTHFR 677 TT and CT polymorphisms in both the case and control group (p < 0.05), with the trend of increased IHD risk resulting from increased air pollution. These results associate the potential inflammatory pathway with air pollution and the folate pathway with MTHFR polymorphism. Future intervention studies can be designed to mitigate MTHFR enzyme deficiencies resulting from gene polymorphisms to prevent IHDs for at-risk populations.

1. Introduction

Ischemic heart disease (IHD) remains the leading cause of death and disability around the world and in the United States (U.S.), accounting for 45.1% of deaths in the U.S. [1,2,3,4,5]. The estimated cost of IHD in 2012 to 2013 was $199.6 billion [4]. By the year, 2030, the medical costs of IHD are projected to increase by about 100% [4]. Coronary artery disease (CAD), a common cardiovascular disease (CVD), is also known as IHD, which includes a group of diseases across the spectrum of chronic to acute stages, including angina, acute myocardial infarction (MI), and sudden cardiac death [6]. Thus, the underlying mechanism of developing IHDs involves atherosclerosis of the blood vessels [3]. Among various risk factors, elevated plasma homocysteine level (hyperhomocysteinemia) is an independent predictor for IHD and subsequent death [7,8]. Homocysteine is positively associated with accelerating atherosclerosis, deterioration of endothelial progenitor cells, deregulating lipid metabolism, altering platelet activity, stimulating vascular smooth-muscle-cell proliferation, and inducing thrombosis [9,10,11]. The most investigated genetic variations affecting the homocysteine levels are the methylenetetrahydrofolate reductase (MTHFR) polymorphisms, which have two common variations: C677T (rs 1801133) and A1298C (rs 1801131) [12]. These two polymorphism-mutations have been associated with impaired or inactivated MTHFR enzyme and aberrant folate metabolism, leading to increased homocysteine and DNA hypomethylation [12]. Multiple meta-analysis studies also showed that MTHFR polymorphisms had a large effect on hyperhomocysteinemia in geographic regions of the world with low folate intakes [13,14,15].
Additional studies have supported the impact of air pollutants on the development of IHDs [16,17,18,19]. For instance, pollutants could cause the rupture of existing atherosclerotic plaque in the coronary vessels, triggering acute MI and stroke [20]. Further studies presented increased IHD-caused mortality and hospitalizations in areas with high concentrations of pollutants [17,20]. Fine particulate matters smaller than 2.5 microns (PM2.5) in air pollution can be inhaled into the lungs and pass through the lungs to reach cardio-pulmonary blood vessels to trigger inflammations [21]. These pollutants are mostly generated from traffic-related fuel combustion [18,20,22]. Over time, pollutants advance atherosclerosis [16,19] and consequent IHDs, which may lead to acute CAD and MI events.
Several meta-analysis studies have presented a significant association between MTHFR polymorphisms and IHD risk [13,14,15], while others failed to identify the association [23,24,25,26]. The discordance might be a result of the limited number of studies [23,24] or limited study populations from certain geographical areas [26,27]. Therefore, we performed an updated meta-analysis, with the addition of air pollution as a contributing factor to gene polymorphism and disease risk, by using meta-prediction analytic techniques to approach the heterogeneity of the findings and to fill the gap in the existing knowledge. In summary, we conducted a meta-prediction study using machine-learning analytics to examine (1) the associations for various IHD disease subtypes in adult populations, and (2) the impact of air pollution on the MTHFR polymorphisms and associated IHD risks.

2. Materials and Methods

2.1. Study Selection and Data Extraction

We conducted a comprehensive literature search from various databases, including MEDLINE through PubMed, Cochrane Library, Embase, EBSCO, Airiti Library, and WangFang databases for studies published between 1996 and 2017. This study was conducted based on the guidelines for Preferred Reporting Items of Systematic Reviews and Meta-Analysis (PRISMA) [28] and Systematic Reviews of Genetic Association Studies [29]. The search keywords and subject terms included “methylenetetrahydrofolate reductase” or “MTHFR” or “methylenetetrahydrofolate reductase gene polymorphisms” or “MTHFR polymorphisms” or “MTHFR variant” and “ischemic heart disease” or “IHD” or “heart disease” or “HD” or “coronary artery disease” or “CAD” or “myocardial infarction” or “MI” and “case control” or “case-control” or “meta-analysis”. We also performed cross-reference checks from previous meta-analyses (see meta-analysis reference list following Supplementary Materials Table S1). All databases were searched repeatedly at 3 different times at least 3 months apart until no additional eligible papers were identified. Two raters conducted the data collection and data extraction. The discrepancy between raters was reduced to zero prior to the data analyses.
An article was eligible if it: (1) presented the association of the MTHFR C677T and A1298C polymorphisms and IHD risk, with three genotype-allele counts in both the case and control groups; (2) was written in English, or (3) was written in non-English, but provided tables with genotype allele counts for both the case and control groups. IHD types, as presented in the literature for this meta-analysis, included CAD for chronic IHD and MI for acute stage IHD. The article was excluded if it: (1) did not have complete genotype allele counts in both the case and control groups, or (2) was written in non-English without tables listing the genotype allele counts.
Based on the inclusion/exclusion criteria, we first identified 469 articles. A total of 242 articles were excluded because they lacked the data for the genotypes per case and control groups. Another 113 articles were eliminated due to missing or incomplete genotype allele counts. We also removed two articles because each of them reused datasets that were in another included article [30,31]. As a result, we included 112 articles, with 123 studies in the final analysis (Figure 1). Eleven additional studies are described in the following section. Four articles (Supplementary Materials reference list, 37, 51, 69, and 79, following Supplementary Materials Table S1) included both gender groups of male and female, yielding 4 additional study groups; an article (Supplementary Materials reference list, 78) included 2 IHD types (CAD and MI), yielding 1 additional study group; an article (Supplementary Materials reference list, 89) included racial-ethnic groups of African and Caucasian, and gender groups of male and female, yielding 3 additional study groups; an article (Supplementary Materials reference list, 91) with MI included two diabetes mellitus (DM) groups (with DM and non-DM), yielding 1 additional study group; an article (Supplementary Materials reference list, 93) with MI included gender groups of male and female, yielding 1 additional study group; and an article (Supplementary Materials reference list, 130) with CAD included gender groups of male and female, yielding 1 additional study group.

2.2. Quality Assessment

The quality of each selected study was evaluated using a quality-assessment tool based on criteria appropriate for meta-analyses [28,32]. The examined domains included external validity (10 items, score range of 0–11), internal validity (12 items, score range of 0–12), and quality of reporting (6 items, score range from 0–6). The total scores ranged from 0 to 29, with a higher score indicating higher quality [32]. The quality score of all included studies ranged from 15 to 28. Therefore, all studies reached a minimal score of 15, above the midpoint of the total possible score for the quality assessment (Supplementary Materials Table S1).

2.3. Data Synthesis and Analysis

The pooled risk ratios (RR), odds ratio (OR), and 95% CI for the associations of MTHFR polymorphisms with IHD risk were calculated. The comparison between pooled RRs and ORs indicated that RRs were more conservative, with a lower type I error [32]. Additionally, RR, as a standardized ratio, is needed for the multi-dimensional gene-environment interaction analysis in this study [32]. Therefore, RRs were used to demonstrate the associations between MTHFR polymorphisms and IHD risk in this report. The StatsDirect, version 3.0.193 (StatsDirect Ltd., Cambridge, UK), was used to perform pooled analyses. All p values were two-tailed, with a significance level at 0.05. The sensitivity analyses were performed for studies that presented significant results on the Hardy-Weinberg Equilibrium (HWE) in the genotype distribution (Supplementary Materials Table S1). For those studies that did not meet the HWE (Supplementary Materials reference list, 47 (Czech), 58 (Netherlands), 77 (U.S.), 86 (Costa Rica), 112 (China), and 128 (India)), we performed sensitivity analyses for the HWE subgroups, and found similar trends and the same directions for the risk effects. We further examined the sensitivity of including data per source of control for the healthy control versus other hospital patients, and for the quality score. Studies that met the HWE compared to those that did not, studies with a low quality score compared to those with a high score, and studies with healthy controls compared to those with other hospital patients had presented the same direction with similar risk effects. Therefore, all studies were included in the final meta-analysis, as the intent of a meta-analysis is to provide summative evidence with pooled analyses, which is the recommended approach [32,33]. An Egger’s test and funnel plots were used to detect publication bias [34,35]. A random-effects model was used in the tests of association when the heterogeneity tests were significant, with p < 0.05. To identify the sources of heterogeneity in the associations between MTHFR polymorphisms and IHD risk, we conducted subgroup analyses by geographic regions, ethnic groups, gender, IHD sub-types, air pollution as an environmental factor, sources of control, and quality score. We used the JMP 13 Pro program (SAS Institute, Cary, NC, USA) to examine the gene-air pollution interaction on IHD risk [36,37], and generated the geographic information system (GIS) maps to visualize the geographic distribution and pattern of MTHFR polymorphisms and their associations with IHD risks [38]. These GIS maps were drawn by country distribution (through the graph builder in the JMP Pro program per outcomes of interest, such as percent polymorphism or IHD risks) and were helpful for visually identifying geographic patterns, and to manage the geospatial data.
We applied recursive partition trees in the JMP Pro 13 program to examine how an independent variable (e.g., air pollution) can make a decisive split of the data by partitioning the groups (such as air pollution levels) into the pairs of subgroups with reference to the dependent variable (the percentage of polymorphisms and IHD risks). The recursive partition tree does this by exhaustively searching all possible groupings [36]. These partitions of the data are done recursively, forming a tree of decision rules until the desired fit is reached. The partition process is driven by the Gini impurity criterion, which is based upon information theory [39,40]. Gini is a measure of group impurity, which is the inverse of group homogeneity, and, thus, a smaller Gini is better. The goodness of the partition result can be judged by using the Akaike’s information criterion correction (AICc). A smaller AICc suggests less complexity and a better model for optimization (a balance between fitness and parsimony) [41,42,43]. We entered the air-quality data for various countries using the guidelines from the World Health Organization on air quality measures, the death rates from air pollution (AP death) (Level 1: <50 deaths per million, Level 2: 51–100 deaths per million, Level 3: 101–250 deaths per million, Level 4: 251–400 deaths per million, and Level 5: >401 deaths per million of population) [44,45]. We further verified these levels with current scales on air pollution data [46,47,48,49], and used the most complete and current scaled air pollution data for the analyses. Only one study from Ireland (Supplementary Materials reference list, 72) presented a Level 1 air pollution level, thus, we combined Level 1 and Level 2 together for the grouping analysis.
Both GIS maps and recursive partition trees are common machine-learning analytical techniques for handling multidimensional and/or large-scale datasets. Different from conventional hypothesis testing, machine-learning or big-data analytics does not start with a pre-determined hypothesis. Rather, data-driven pattern recognition plays the central role. To triangulate and cross-validate the findings, we used machine-learning based validation analytics, including partition trees, and nonlinear-association fit modeling to explore the sources of heterogeneity in addition to the conventional approaches [32,36]. We used a bi-variate nonlinear fit plot to visualize the distribution patterns based on air pollution levels across countries on the pooled percentages of MTHFR genotypes per case and control groups, and IHD risks to cross-validate the results [32]. To compare AICc results with the partition trees, we used the Tukey’s test [36] to examine whether partition trees and Tukey tests concurred with each other. The aim of the meta-predictive analysis was to generate more precise predictions by integrating data from multidimensional sources. The main purpose for using both conventional statistical and machine learning (e.g., recursive partition trees) methods is to verify the results by cross-validation, including AICc, yielding a more accurate meta-prediction.

3. Results

3.1. Sample Characteristics and Genotype Frequency

The most investigated racial or ethnic populations in these studies were Caucasian (60 studies), followed by East Asian (32 studies), South Asian (13 studies), Middle-Eastern (7 studies), African (6 studies), mixed groups (3 studies), and Hispanic (2 studies). Among the 123 case-control studies, 93 studies were conducted in CAD patients and 30 in MI patients (Table 1 and Table 2). No significant association was found between MTHFR 1298 polymorphism and IHD risk. The frequencies of the MTHFR 677 homozygous TT genotype were highest in Hispanic populations (28.90%), followed by East Asian (20.23%), African (12.01%), Caucasian (11.10%), Middle Eastern (10.28%), mixed group (9.78%), and South Asian (3.76%). For pooled analyses in MTHFR 1298, the frequencies of the homozygous CC genotype were highest in the mixed group (27.36%), South Asian (16.55%), Caucasian (10.40%), African (6.86%), Middle Eastern (6.4%), and East Asian (3.70%) populations (Supplementary Materials Figure S1a). The frequencies of MTHFR 1298 polymorphisms are presented in the Supplementary Materials Figure S1b.

3.2. Pooled Meta-Analysis

Table 1 presents the schema of the significant findings across IHD types with CAD and MI subtypes for the MTHFR 677 genotypes and risk of IHD; Table 2 presents the results of the detailed analysis. For analyses of the IHD subtypes of CAD and MI, detailed results are presented in the Supplementary Materials Tables S2 and S3. For all included study groups, the pooled analysis presented MTHFR 677 polymorphisms as risks for IHD (TT: Risk ratio (RR) = 1.23, p < 0.0001; CT: RR = 1.04, p = 0.0028, and TT plus CT: RR = 1.09, p < 0.0001), with the T allele having an 11% greater risk of developing IHD (RR = 1.11, p < 0.0001). In contrast, MTHFR 677 CC wildtype (RR = 0.91, p < 00001) was protective against IHD, and holders of the C allele would have a 4% lesser chance of developing IHD (RR = 0.96, p <0.0001) (Table 2). The results of the pooled odds ratio (OR) for the MTHFR 677 polymorphisms presented greater ratios, as compared to the RRs, for IHD risk (TT: OR = 1.28, p < 0.0001; CT: OR = 1.08, p = 0.0012, and TT plus CT: OR = 1.22, p < 0.0001), with the T allele having a 17% greater odds of developing IHD (OR = 1.17, p < 0.0001) versus the CC wildtype (OR = 0.82, p < 00001), which presented a lesser odds as compared to RR against IHD, with holders of the C allele having a 15% lesser odds of developing IHD (OR = 0.85, p < 0.0001). Hence, the comparison between the pooled RRs and ORs indicated that RRs were more conservative, with a lower type I error [32]. As the RR is used for the standardized ratio in the multi-dimensional gene-environment interaction analysis [32], we used RRs to demonstrate the associations between MTHFR polymorphisms and IHD risk in this report.

3.3. Subgroup Analysis by Ethnicity and Countries

To identify sources of heterogeneity, subgroup analysis per ethnic groups showed that MTHFR 677 TT and CT genotypes were associated with IHD risk. The ranking of IHD risk with the MTHFR 677 TT genotype from high to low according to ethnic groups was Middle Eastern (RR = 2.62), African (RR = 2.14), South Asian (RR = 1.51), and East Asian (RR = 1.31); and with the CT genotype, the ranking of risks was South Asian (RR = 1.32) and Middle Eastern (RR = 1.18) (all p < 0.01). It is worthy to point out that the Middle Eastern and African samples had RR > 2, and, as in biological studies, RR > 2 could infer causality [33,50,51,52]. On the other hand, the MTHFR 677 CC wildtype played a protective role against IHD and ranked as follows: Middle Eastern (RR = 0.76), East Asian (RR = 0.84), and South Asian (RR = 0.91) (all p < 0.01) (Table 2). For the subgroup analysis by countries, we further examined whether MTHFR 677 TT genotypes posed a risk (RR > 1), were protective (RR < 1), or had a mixed effect (RR ~ 1) (Table 2). The countries that had MTHFR 677 TT as a risk genotype included Australia, European countries (Poland, Russia, Slovakia, Sweden, Croatia, Czech Republic, Netherlands, Switzerland, France, Portugal, United Kingdom, and Ireland), Canada, South America (Brazil and Mexico), Asian countries (Japan, China, and India), Middle-Eastern countries (Israel, Turkey, Iran, Egypt, Saudi Arabia, and Lebanon), and African countries (Tunisia and Morocco) (Supplementary Materials Figure S2a). It is worthy to point out the countries with RR > 2 because, as in biological studies, RR > 2 could infer causality [33,50,51,52]. These countries with elevated risks included India, Turkey, and Tunisia (Figure 2, forest plots). In Figure 2, six studies from India presented a very wide range of the 95% confidence interval (CI) (Supplementary Materials reference list, 125–127, 129, and 130, with both gender groups of male and female), with four studies presenting zero counts on the MTHFR 677 TT genotype in their control groups, and two studies (Supplementary Materials reference, 130, included two gender groups of male and female) presenting zero counts on the MTHFR TT genotype in the case groups. The countries that presented MTHFR 677 TT as a protective genotype (RR < 1) included three European countries (Hungary, Italy, and Spain), one Central American country (Costa Rica), two Asian countries (South Korea and Pakistan), and South Africa (study population consisted of descendants of Indian immigrants) (Supplementary Materials Figure S2b). For Germany, United States, and Taiwan, the MTHFR 677 TT genotype presented a mixed IHD risk (Supplementary Materials Figure S2c).

3.4. Subgroup Analysis by IHD Subtypes

We also analyzed subgroups by the IHD subtypes of CAD and MI to present the source of heterogeneity. For CAD subtype, risk genotypes were MTHFR 677 TT (RR = 1.27, p < 0.0001), CT (RR = 1.04, p = 0.0074), and TT plus CT (RR = 1.09, p < 0.0001). In contrast, the MTHFR 677 CC genotype played a protective role against CAD (RR = 0.91, p < 0.0001). Likewise, the MTHFR 677 T allele increased the risk of CAD by 12% (RR = 1.12, p < 0.0001), while its C allele could reduce CAD risk by 5% (RR = 0.95, p < 0.0001). Further subgroup analyses per ethnic subgroups revealed that the MTHFR 677 TT and CT genotypes were associated with CAD risk in Caucasian, East Asian, South Asian, Middle-Eastern, and African groups, with RRs ranging from 1.11 to 2.56 for TT, and 1.18 to 1.36 for CT. On the other hand, the MTHFR 677 CC wildtype played a protective role against CAD in East Asian, South Asian, and Middle-Eastern groups, with RRs ranging from 0.76 to 0.89 (Supplementary Materials Table S2). For the MI subtype, MTHFR 677 TT plus CT were risk genotypes (RR = 1.06, p = 0.0187), while the CC genotype played a protective role against MI (RR = 0.93, p = 0.0148). Individuals with the MTHFR 677 T allele had a 6% higher risk of having MI (RR = 1.06, p < 0.0073), while those with the C allele had a 3% lower chance of developing MI (RR = 0.97, p = 0.0223). Further subgroup analyses by ethnic groups showed that East Asian groups with the TT genotype had an increased risk of developing MI (RR = 1.30, p < 0.0001) (Supplementary Materials Table S3).

3.5. Heterogeneous Findings by GIS Map

To validate the heterogeneous findings, we further utilized a geographic information system to visualize the variability in regional distributions [53]. The global maps presented variations in the distribution of MTHFR 677 polymorphism rates and their roles in IHD risk across regions (Supplementary Materials Figure S3a,b). In the first two GIS maps, the continuous color spectrum from yellow to red represented the increasing rates of TT and TT plus CT genotypes, and, in the third map, the red-to-green color spectrum demonstrated a risk gradient, with red indicating IHD risk and green indicating protective effects. Consistent with the aforementioned subgroup analyses, GIS maps presented that MTHFR 677 TT and TT plus CT increased the risk of IHD in most countries except Italy, Spain, and South Korea.

3.6. Meta-Prediction: MTHFR Polymorphisms and Air Pollution Associated with Risk of IHD

For meta-prediction, we performed both a partition tree analysis and a Tukey’s test to examine the potential interaction between independent variables, such as levels of air pollution as measured by the death rates associated with air pollution per country (AP), source of controls, and the quality score, with outcome variables of MTHFR polymorphisms and IHD risk. We did not find significance with other factors except for the AP on the interaction between polymorphisms and IHD risks. The annual AP by country was reported and classified by the World Health Organization (WHO) using the following criteria: (1) Level 2 = 51–100 deaths/million population (DMP); (2) Level 3 = 101–250 DMP; and (3) Level 4 = 251–400 or greater DMP [49]. We present the partition tree (split groups) and Tukey’s test results side by side for the MTHFR 677 genotypes and AP (Table 3). The meta-prediction using the partition tree and the Tukey’s tests were not performed for the MTHFR 1298 genotypes because that genotype was rarely studied to date. There were significant increases in the percentages of MTHFR 677 TT and CT polymorphism between AP levels 2 and 3 (p < 0.05), and levels 2 and 4 (p < 0.01) for both the case and control groups. The partition tree results, depicted in Table 3, indicate a clear and consistent pattern: No matter what the variable is (TT%ct, TT%ca, CT%ct, CT%ca, etc.), the mean percent of polymorphisms of TT or CT were lower with lower levels of AP death, and increased with higher levels of AP death. In addition, strikingly higher polymorphisms were noted as AP rates increased for the IHD case group as compared to the corresponding control group. In contrast, the percentages of the MTHFR 677 CC wildtype decreased with increased levels of AP death for both the case and control groups. The risks of IHD from MTHFR 677 TT plus CT polymorphisms presented a greater risk with lower levels of AP death, and was significantly higher between AP levels 2 and 3 (p < 0.001), and levels 2 and 4 (p < 0.05). It is important to point out that, unlike ANOVA, in which levels are grouped by the analyst, the partition done by the decision tree is data-driven i.e., the algorithm determines what partition could result in better purity (homogeneity), meaning that the associational pattern is less likely to be a statistical artifact capitalized on chance. The nonlinear curves were further revealing of the differences on the percentage of MTHFR polymorphisms at different AP levels (Figure 3). With a change in AP levels from low (Level 2) to high (Level 3 and 4).
There was a substantial increase in the percentages of MTHFR 677 polymorphisms for both TT and TT plus CT genotypes in both case and control groups.

4. Discussion

In summary, we have updated and performed the pooled analyses with 123 studies of 29,697 IHD cases and 31,028 controls for MTHFR 677, and 18 studies of 7158 IHD cases and 5482 controls for MTHFR 1298 polymorphisms, adding the environmental air pollution factor. We extended the findings of previous meta-analyses [26,54,55,56], with a significant association between MTHFR 677 polymorphism and the risk of developing IHD identified. A great deal of heterogeneity was noted across geographic areas and populations based on the IHD subtypes of CAD and MI. To expand the findings from previous meta-analyses, we further conducted a meta-prediction to examine a potential interaction between air pollution and MTHFR polymorphisms on IHD risk. The nonlinear curve plots demonstrated that the percentages of MTHFR 677 polymorphism in both the case and control groups increased substantially in regions with higher levels of air pollution. Compared to control groups, the case groups showed higher percentages of MTHFR 677 TT and CT polymorphisms where the level of air pollution was higher. These findings are consistent with prior studies that found significantly increased MTHFR gene polymorphisms with increased air pollution levels in various diseases, including various cancers, Alzeimer’s disease, and hypertensive disorders during pregnancy (HDP) [32,57,58,59,60]. Also, an increased air pollution level was associated with MTHFR polymorphism associated with HDP risk [57].
To examine possible sources of heterogeneity, we examined data using other factors that became available with this study, including air pollution, HWE status (met or unmet), quality score (high or low), and the sources of controls (healthy control or inpatients as controls). We observed similar trends of these factors on the polymorphism and IHD risks, with the same directions for the risk effects. For the meta-prediction analysis, we examined the effects of air pollution, sources of controls, and quality score status on the polymorphism and IHD risks. Continuing the findings from previous studies [32,57,58,59,60], air pollution, as measured by the AP, was the only significant factor for the polymorphism and disease risks. It is worth noting the limitations of the differences found in Table 3 that may be caused by additional factors, such as cultural, economic, or industrial factors, but not from air pollution. In addition, despite the potential risk of ecological fallacy, in some cases, global data must be used because collecting individual data is impossible. Air pollution is a typical example. Although techniques of monitoring how much pollutant is absorbed by each individual has been under development by taking the immediate surroundings, individuals’ biophysical characteristics, and individuals’ space-time activities, such as culture, economic, or industrial factors, into account [61], these methods are not prevalent across countries.
Previous studies have shown that air pollution directly harms the cardiovascular system [18,62]. The association of environmental pollution with increased polymorphism rates is alarming for human health. The potential detriments of this association could be exacerbated by the greenhouse effect from air pollution because heat can further diminish the functioning of MTHFR enzymes. Decreased MTHFR enzyme functioning can compromise methylation pathways and result in elevation of plasma homocysteine, an independent risk factor for IHD [63,64]. Furthermore, both MTHFR 677 polymorphisms and air pollution may compound the detriments of hyperhomocysteinemia to compromise the epigenetic health status for persons with IHD.
For the meta-predictive analysis, meta-regression using linear modeling, such as Pearson’s correlation and ordinary least squares regression, has been criticized as overly simplistic. When the underlying data structure is nonlinear, such as the association between death from air pollution with polymorphism and IHD risk, then nonlinear modeling is more appropriate [64,65,66]. Therefore, we further demonstrated the association of MTHFR 677 polymorphism with IHD risks using machine learning meta-predictive analytics, including GIS maps, to present the heterogeneity of polymorphisms worldwide. The great heterogeneity in the relationship between MTHFR 677 polymorphism and IHD risk across ethnic groups and regions may be attributed to (1) the human migration paths and variations in the percentages of MTHFR 677 polymorphism among different ethnic groups across countries; (2) various IHD disease subtypes; (3) various levels of air pollution across different geographic areas; and (4) differences in lifestyles (e.g., eating behaviors and food choices). For instance, in countries where people are acculturated to eat plenty of vegetables, diet may have modified the effects of MTHFR 677 polymorphisms. Korea has a well-known taste for fermented cabbage, and Mediterranean cuisines, such as Italian or Spanish, include plants, such as peppers, tomatoes, eggplant, and zucchini; in such cultures, MTHFR 677 polymorphism was surprisingly protective against IHD [67,68]. MTHFR 677 polymorphism had a significant association with CAD, but not acute MI, in this analysis. This lack of significance for MI conditions could be due to a smaller number of studies (30 studies); perhaps because it might also be more feasible to conduct studies during the chronic disease state, such as CAD (93 studies), than during the acute MI stage. An additional potential explanation has been posited that homocysteine is considered an endothelial toxin that promotes thrombotic tendency and atherosclerosis for the pathogenesis of CAD, but not for that of acute MI [13]. To solidify the effects of MTHFR 677 polymorphism risk across populations, additional studies with purposeful designs and interventions can be helpful.
Despite the clear strengths of our study, some additional limitations merit some consideration. First, the data for potential confounding factors (e.g., cultural, economic, or industrial factors) were not provided in original studies, therefore, further investigations are needed to determine whether the differences in the magnitude of association between air pollution with MTHFR polymorphism and IHD risk are attributed to air pollution or other confounding factors. Second, to provide comprehensive pooled analyses, we included studies with small sample sizes and a wide range of 95% CI, with some studies that did not meet the HWE. However, based on our sensitivity analysis by the HWE status (met vs. unmet), including underpowered studies had no effect on the direction of effects based on MTHFR polymorphisms and IHD risks. Therefore, we included all studies in the meta-analyses, as the intent of meta-analysis is to provide summative filtered evidence with pooled analyses. Given the limitation and quality of original studies, future studies are needed to continue examining the effects of gene-environment interactions from air pollution and gene polymorphisms on health outcomes across the disease spectrum. The results of heightened disease risks from air pollution and gene polymorphisms for IHD could bring attention to the development of health policy for clean air environment and interventions to mitigate enzyme deficiency in folate metabolism pathways to prevent IHD.

5. Conclusions

MTHFR polymorphisms were risks for IHD and for gene-environment interactions, with the 677 wild type plying a greater protective role in countries with higher air pollution levels than in those countries with lower air pollution levels. Epigenetic factors, including environmental toxicants from air pollution and healthy lifestyles, such as food choices to detox, may affect the development of IHDs by modifying gene expressions in methylation pathways. Proactive strategies could be implemented in geographic regions with significant air pollution levels to mitigate the effects of MTHFR enzyme deficiencies that are a result of polymorphisms to prevent IHDs and to promote the health of susceptible, at-risk populations.

Supplementary Materials

The following are available online at https://www.mdpi.com/1660-4601/15/7/1453/s1, Table S1: Summary of participant characteristics and MTHFR 677 and 1298 loci distributions for studies on ischemic heart disease (IHD) included in this review, by geographic location (112 articles, 123 study groups based on counts available from both control and IHD groups). Table S2: Pooled analysis: MTHFR C677T genotypes and risks of coronary artery disease (CAD) by ethnicity (93 Studies). Table S3: Pooled analysis: MTHFR C677T genotypes and risks of myocardial infarction (MI) by ethnicity (30 Studies). Figure S1a. MTHFR C677T percentage of polymorphism per control and IHD case groups; Figure S1b. MTHFR A1298C percentage of polymorphism per control and IHD case groups. Figure S2a. Forest plot for meta-analysis of MTHFR 677 polymorphism by TT genotype, countries of Caucasian with risks >1. Figure S2b. Forest plot for meta-analysis of MTHFR 677 polymorphism by TT genotype, countries with risks <1. Figure S2c. Forest plot for meta-analysis of MTHFR 677 polymorphism by TT genotype, countries with risks varied around 1. Figure S3a. Geographic information maps for percentages of MTHFR 677 TT genotype per control and IHD case groups, and its association with IHD risks. Figure S3b. Geographic information maps for percentages of MTHFR 677 TT plus CT genotypes per control and ischemic heart disease (IHD) case groups, and their associations with IHD risks.

Author Contributions

All authors contributed to the conception of the research idea, coding of the data, data analyses, writing the manuscript, as well as critically analyzing the content and improving the quality of the manuscript.

Funding

This research was supported in part by the Augusta University Research Start up fund awarded to the corresponding author on the software and the APC.

Acknowledgments

The authors acknowledge the assistance from Suh-Mian Wu, who retrieved literature and rechecked data coding in earlier stage of study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Go, A.S.; Mozaffarian, D.; Roger, V.L.; Benjamin, E.J.; Berry, J.D.; Blaha, M.J.; Dai, S.; Ford, E.S.; Fox, C.S.; Franco, S.; et al. Executive summary: Heart disease and stroke statistics—2014 update: A report from the American Heart Association. Circulation 2014, 129, 399–410. [Google Scholar] [CrossRef] [PubMed]
  2. Nichols, M.; Townsend, N.; Scarborough, P.; Rayner, M. Cardiovascular disease in Europe 2014: Epidemiological update. Eur. Heart J. 2014, 35, 2950–2959. [Google Scholar] [CrossRef] [PubMed]
  3. Nabel, E.G.; Braunwald, E. A tale of coronary artery disease and myocardial infarction. N. Engl. J. Med. 2012, 366, 54–63. [Google Scholar] [CrossRef] [PubMed]
  4. Benjamin, E.J.; Blaha, M.J.; Chiuve, S.E.; Cushman, M.; Das, S.R.; Deo, R.; de Ferranti, S.D.; Floyd, J.; Fornage, M.; Gillespie, C.; et al. Heart disease and stroke statistics—2017 update: A report from the American Heart Association. Circulation 2017, 135, e146–e603. [Google Scholar] [CrossRef] [PubMed]
  5. Finegold, J.A.; Asaria, P.; Francis, D.P. Mortality from ischaemic heart disease by country, region, and age: Statistics from World Health Organisation and United Nations. Int. J. Cardiol. 2013, 168, 934–945. [Google Scholar] [CrossRef] [PubMed]
  6. Abubakar, I.I.; Tillmann, T.; Banerjee, A. Global, regional, and national age-sex specific all-cause and cause-specific mortality for 240 causes of death, 1990–2013: A systematic analysis for the global burden of disease study 2013. Lancet 2015, 385, 117–171. [Google Scholar] [CrossRef]
  7. Peng, H.Y.; Man, C.F.; Xu, J.; Fan, Y. Elevated homocysteine levels and risk of cardiovascular and all-cause mortality: A meta-analysis of prospective studies. J. Zhejiang Univ. Sci. B 2015, 16, 78–86. [Google Scholar] [CrossRef] [PubMed]
  8. Kim, Y.; Je, Y. Dietary fiber intake and total mortality: A meta-analysis of prospective cohort studies. Am. J. Epidemiol. 2014, 180, 565–573. [Google Scholar] [CrossRef] [PubMed]
  9. Zhang, D.; Sun, X.; Liu, J.; Xie, X.; Cui, W.; Zhu, Y. Homocysteine accelerates senescence of endothelial cells via DNA hypomethylation of human telomerase reverse transcriptase. Arterioscler. Thromb. Vasc. Biol. 2015, 35, 71–78. [Google Scholar] [CrossRef] [PubMed]
  10. Kones, R. Molecular sources of residual cardiovascular risk, clinical signals, and innovative solutions: Relationship with subclinical disease, undertreatment, and poor adherence: Implications of new evidence upon optimizing cardiovascular patient outcomes. Vasc. Health Risk Manag. 2013, 9, 617–670. [Google Scholar] [CrossRef] [PubMed]
  11. Krishna, S.M.; Dear, A.; Craig, J.M.; Norman, P.E.; Golledge, J. The potential role of homocysteine mediated DNA methylation and associated epigenetic changes in abdominal aortic aneurysm formation. Atheroscler. 2013, 228, 295–305. [Google Scholar] [CrossRef] [PubMed]
  12. Gil Nunez, A.C. Homocysteine: Vascular risk factor? Neurologia 2004, 19, 1–4. [Google Scholar] [PubMed]
  13. Luo, H.; Liu, B.; Hu, J.; Wang, X.; Zhan, S.; Kong, W. Hyperhomocysteinemia and methylenetetrahydrofolate reductase polymorphism in cervical artery dissection: A meta-analysis. Cerebrovasc. Dis. 2014, 37, 313–322. [Google Scholar] [CrossRef] [PubMed]
  14. Yadav, S.; Hasan, N.; Marjot, T.; Khan, M.S.; Prasad, K.; Bentley, P.; Sharma, P. Detailed analysis of gene polymorphisms associated with ischemic stroke in south Asians. PLoS ONE 2013, 8, e57305. [Google Scholar] [CrossRef] [PubMed]
  15. Klerk, M.; Verhoef, P.; Clarke, R.; Blom, H.J.; Kok, F.J.; Schouten, E.G. Mthfr 677c→t polymorphism and risk of coronary heart disease: A meta-analysis. JAMA 2002, 288, 2023–2031. [Google Scholar] [CrossRef] [PubMed]
  16. Langrish, J.P.; Li, X.; Wang, S.; Lee, M.M.; Barnes, G.D.; Miller, M.R.; Cassee, F.R.; Boon, N.A.; Donaldson, K.; Li, J.; et al. Reducing personal exposure to particulate air pollution improves cardiovascular health in patients with coronary heart disease. Environ. Health Perspect. 2012, 120, 367–372. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Gan, W.Q.; Koehoorn, M.; Davies, H.W.; Demers, P.A.; Tamburic, L.; Brauer, M. Long-term exposure to traffic-related air pollution and the risk of coronary heart disease hospitalization and mortality. Environ. Health Perspect. 2011, 119, 501–507. [Google Scholar] [CrossRef] [PubMed]
  18. Brook, R.D.; Rajagopalan, S.; Pope, C.A., 3rd; Brook, J.R.; Bhatnagar, A.; Diez-Roux, A.V.; Holguin, F.; Hong, Y.; Luepker, R.V.; Mittleman, M.A.; et al. Particulate matter air pollution and cardiovascular disease: An update to the scientific statement from the American Heart Association. Circulation 2010, 121, 2331–2378. [Google Scholar] [CrossRef] [PubMed]
  19. Brook, R.D.; Franklin, B.; Cascio, W.; Hong, Y.; Howard, G.; Lipsett, M.; Luepker, R.; Mittleman, M.; Samet, J.; Smith, S.C.; et al. Air pollution and cardiovascular disease: A statement for healthcare professionals from the expert panel on population and prevention science of the American Heart Association. Circulation 2004, 109, 2655–2671. [Google Scholar] [CrossRef] [PubMed]
  20. Dominici, F.; Peng, R.D.; Bell, M.L.; Pham, L.; McDermott, A.; Zeger, S.L.; Samet, J.M. Fine particulate air pollution and hospital admission for cardiovascular and respiratory diseases. JAMA 2006, 295, 1127–1134. [Google Scholar] [CrossRef] [PubMed]
  21. Ruckerl, R.; Ibald-Mulli, A.; Koenig, W.; Schneider, A.; Woelke, G.; Cyrys, J.; Heinrich, J.; Marder, V.; Frampton, M.; Wichmann, H.E.; et al. Air pollution and markers of inflammation and coagulation in patients with coronary heart disease. Am. J. Respir. Crit. Care Med. 2006, 173, 432–441. [Google Scholar] [CrossRef] [PubMed]
  22. Peters, A.; Dockery, D.W.; Muller, J.E.; Mittleman, M.A. Increased particulate air pollution and the triggering of myocardial infarction. Circulation 2001, 103, 2810–2815. [Google Scholar] [CrossRef] [PubMed]
  23. Van Meurs, J.B.; Pare, G.; Schwartz, S.M.; Hazra, A.; Tanaka, T.; Vermeulen, S.H.; Cotlarciuc, I.; Yuan, X.; Mälarstig, A.; Bandinelli, S.; et al. Common genetic loci influencing plasma homocysteine concentrations and their effect on risk of coronary artery disease. Am. J. Clin. Nutr. 2013, 98, 668–676. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Wald, D.S.; Morris, J.K.; Wald, N.J. Reconciling the evidence on serum homocysteine and ischaemic heart disease: A meta-analysis. PLoS ONE 2011, 6, e16473. [Google Scholar] [CrossRef] [PubMed]
  25. Clarke, R.; Bennett, D.A.; Parish, S.; Verhoef, P.; Dötsch-Klerk, M.; Lathrop, M.; Xu, P.; Nordestgaard, B.G.; Holm, H.; Hopewell, J.C.; et al. Homocysteine and coronary heart disease: Meta-analysis of MTHFR case-control studies, avoiding publication bias. PLoS Med. 2012, 9, e1001177. [Google Scholar] [CrossRef] [PubMed]
  26. Chen, L.; Liu, L.; Hong, K.; Hu, J.; Cheng, X. Three genetic polymorphisms of homocysteine-metabolizing enzymes and risk of coronary heart disease: A meta-analysis based on 23 case-control studies. DNA Cell Biol. 2012, 31, 238–249. [Google Scholar] [CrossRef] [PubMed]
  27. Jee, S.H.; Beaty, T.H.; Suh, I.; Yoon, Y.; Appel, L.J. The methylenetetrahydrofolate reductase gene is associated with increased cardiovascular risk in Japan, but not in other populations. Atherosclerosis 2000, 153, 161–168. [Google Scholar] [CrossRef]
  28. Moher, D.; Liberati, A.; Tetzlaff, J.; Altman, D.G. Preferred reporting items for systematic reviews and meta-analyses: The prisma statement. J. Clin. Epidemiol. 2009, 62, 1006–1012. [Google Scholar] [CrossRef] [PubMed]
  29. Sagoo, G.S.; Little, J.; Higgins, J.P. Systematic reviews of genetic association studies. Human genome epidemiology network. PLoS Med. 2009, 6, e28. [Google Scholar] [CrossRef] [PubMed]
  30. Girelli, D.; Friso, S.; Trabetti, E.; Olivieri, O.; Russo, C.; Pessotto, R.; Faccini, G.; Pignatti, P.F.; Mazzucco, A.; Corrocher, R. Methylenetetrahydrofolate reductase C677T mutation, plasma homocysteine, and folate in subjects from northern Italy with or without angiographically documented severe coronary atherosclerotic disease: Evidence for an important genetic-environmental interaction. Blood 1998, 91, 4158–4163. [Google Scholar] [PubMed]
  31. Morita, H.; Taguchi, J.I.; Kurihara, H.; Kitaoka, M.; Kaneda, H.; Kurihara, Y.; Maemura, K.; Minamino, T.; Ohno, M.; Yamaoki, K.; et al. Genetic polymorphism of 5,10-methylenetetrahydrofolate reductase (MTHFR) as a risk factor for coronary artery disease. Circulation 1997, 95, 2032–2036. [Google Scholar] [CrossRef] [PubMed]
  32. Shiao, S.P.K.; Yu, C.H. Meta-prediction of MTHFR gene polymorphism mutations and associated risk for colorectal cancer. Biol. Res. Nurs. 2016, 18, 357–369. [Google Scholar] [CrossRef] [PubMed]
  33. Deeks, J.J.; Higgins, J.; Altman, D.G. Analysing data and undertaking meta-analyses. In Cochrane Handbook for Systematic Reviews of Interventions: Cochrane Book Series; Green, S., Higgins, J., Eds.; John Wiley & Sons, Ltd.: Chichester, UK, 2008; pp. 243–296. [Google Scholar]
  34. Duval, S.; Tweedie, R. Trim and fill: A simple funnel-plot-based method of testing and adjusting for publication bias in meta-analysis. Biometrics 2000, 56, 455–463. [Google Scholar] [CrossRef] [PubMed]
  35. Egger, M.; Davey Smith, G.; Schneider, M.; Minder, C. Bias in meta-analysis detected by a simple, graphical test. BMJ 1997, 315, 629–634. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Speybroeck, N. Classification and regression trees. Int. J. Pub. Health 2012, 57, 243–246. [Google Scholar] [CrossRef] [PubMed]
  37. Jaccard, J.; Becker, M.A.; Wood, G. Pairwise multiple comparison procedures: A review. Psychol. Bull. 1984, 96, 589. [Google Scholar] [CrossRef]
  38. Albrecht, J. Key Concepts and Techniques in GIS; Sage: Thousand Oaks, CA, USA, 2007. [Google Scholar]
  39. Breiman, L.; Friedman, J.H.; Olshen, R.A.; Stone, C.J. Classification and Regression Trees; Chapman and Hall/CRC: Boca Raton, FL, USA, 1984. [Google Scholar]
  40. Bala, M.; Agrawal, R.K. Statistical measures to determine optimal structure of decision tree: One versus one support vector machine. Def. Sci. J. 2010, 60, 399–404. [Google Scholar] [CrossRef]
  41. Akaike, H. Akaike's Information Criterion; Springer: Berlin/Heidelberg, Germany, 2011. [Google Scholar]
  42. Akaike, H. Information theory and an extension of the maximum likelihood principle. In International Symposium on Information Theory; Petrov, B., Csaki, F., Eds.; Akademia Kiado: Budapest, Hungary, 1973; pp. 267–281. [Google Scholar]
  43. Symonds, M.R.; Moussalli, A. A brief guide to model selection, multimodel inference and model averaging in behavioural ecology using Akaike’s information criterion. Behav. Ecol. Sociobiol. 2011, 65, 13–21. [Google Scholar] [CrossRef]
  44. Kenworthy, J.; Laube, F. Urban transport patterns in a global sample of cities & their linkages to transport infrastructure, land use, economics & environment. World Transp. Policy Prac. 2002, 8, 5–20. [Google Scholar]
  45. World Health Organization. The Urban Environment. 2015. Available online: http://www.who.int/heli/risks/urban/urbanenv/en/ (accessed on 27 February 2018).
  46. US Environmental Protection Agency. Air Quality Index Basics. Available online: https://www.airnow.gov/index.cfm?action=aqibasics.aqi (accessed on 2 January 2018).
  47. World Health Organization. Deaths Attributable to Urban Air Pollution. 2004. Available online: http://www.who.int/heli/risks/urban/en/uapmap.1.pdf?ua¼1 (accessed on 2 January 2018).
  48. World Health Organization. Global Health Risks. 2009. Available online: http://www.who.int/healthinfo/global_burden_disease/GlobalHealthRisks_report_full.pdf (accessed on 2 January 2018).
  49. World Health Organization. Global Health Risks. 2012. Available online: https://commons.wikimedia.org/wiki/File:Deaths_from_air_pollution.png (accessed on 2 January 2018).
  50. American Institute for Cancer Research. How AICR Recommendations Cuts Colorectal Cancer Risk for Both Men and Women. Available online: http://www.aicr.org/cancer-research-update/2016/11_02/cru-how-AICR-recommendations-cuts-colorectal-cancer-risk-for-men-and-women.html (accessed on 16 March 2018).
  51. Viera, A.J. Odds ratios and risk ratios: What’s the difference and why does it matter? South Med. J. 2008, 101, 730–734. [Google Scholar] [CrossRef] [PubMed]
  52. Scott, I. Interpreting risks and ratios in therapy trials. Aust. Prescr. 2008, 31, 12–16. [Google Scholar] [CrossRef] [Green Version]
  53. Team, Q.D.; Quantum GIS Geographic Information System. Open Source Geospatial Foundation Project. 2012. Available online: http://qgis.osgeo.org (accessed on 2 January 2018).
  54. Bentley, P.; Peck, G.; Smeeth, L.; Whittaker, J.; Sharma, P. Causal relationship of susceptibility genes to ischemic stroke: Comparison to ischemic heart disease and biochemical determinants. PLoS ONE 2010, 5, e9136. [Google Scholar] [CrossRef] [PubMed]
  55. Drenos, F.; Whittaker, J.C.; Humphries, S.E. The use of meta-analysis risk estimates for candidate genes in combination to predict coronary heart disease risk. Ann. Hum. Genet. 2007, 71, 611–619. [Google Scholar] [CrossRef] [PubMed]
  56. Hou, X.; Chen, X.; Shi, J. Genetic polymorphism of MTHFR C677T and premature coronary artery disease susceptibility: A meta-analysis. Gene 2015, 565, 39–44. [Google Scholar] [CrossRef] [PubMed]
  57. Yang, Y.L.; Yang, H.L.; Shiao, S.P. Meta-prediction of MTHFR gene polymorphisms and air pollution on the risk of hypertensive disorders in pregnancy worldwide. Int. J. Environ. Res. Public Health 2018, 15, E326. [Google Scholar] [CrossRef] [PubMed]
  58. Wu, S.M.; Chen, Z.F.; Young, L.; Shiao, S.P.K. Meta-prediction of the effects of methylenetetrahydrofolate reductase gene polymorphisms and air pollution on risk of Alzheimer’s disease. Int. J. Environ. Res. Public Health 2017, 14, E63. [Google Scholar] [CrossRef] [PubMed]
  59. Gonzales, M.C.; Yu, P.J.; Shiao, S.P.K. Meta-prediction of MTHFR gene polymorphism-mutations and air pollution as risk factors for breast cancer. Nurs. Res. 2017, 66, 152–163. [Google Scholar] [CrossRef] [PubMed]
  60. Lien, S.Y.A.; Young, L.; Gau, B.S.; Shiao, S.P.K. Meta-prediction of MTHFR gene polymorphism-mutations, air pollution, and risks of leukemia among world populations. Oncotarget 2017, 8, 4387–4398. [Google Scholar] [CrossRef] [PubMed]
  61. Lu, Y.; Fang, T.B. Examining personal air pollution exposure, intake, and health danger zone using time geography and 3D geovisualization. ISPRS Int. J. Geo-Inf. 2015, 4, 32–46. [Google Scholar] [CrossRef]
  62. Leclerc, D.; Sibani, S.; Rozen, R. Molecular Biology of Methylenetetrahydrofolate Reductase (MTHFR) and Overview of Mutations/Polymorphisms. In MTHFR Polymorphisms and Disease; CRC Press: Boca Raton, FL, USA, 2005; pp. 15–34. [Google Scholar]
  63. Lee, B.J.; Kim, B.; Lee, K. Air pollution exposure and cardiovascular disease. Toxicol. Res. 2014, 30, 71–75. [Google Scholar] [CrossRef] [PubMed]
  64. Packard, G. Review article: On the use of log-transformation vs. nonlinear regression for analyzing biological power laws. Biol. J. Linn. Soc. 2014, 113, 1167–1178. [Google Scholar] [CrossRef]
  65. Xiao, X.; White, E.P.; Hooten, M.B.; Durham, S.L. On the use of log-transformation vs. nonlinear regression for analyzing biological power laws. Ecology 2011, 92, 1887–1894. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Yu, C.H.; Lee, H.S.; Gan, S.; Brown, E. Nonlinear modeling with big data in SAS and JMP. In Proceedings of the Western Users of SAS Software Conference, Long Beach, CA, USA, 20–22 September 2017. [Google Scholar]
  67. Choi, I.H.; Noh, J.S.; Han, J.S.; Kim, H.J.; Han, E.S.; Song, Y.O. Kimchi, a fermented vegetable, improves serum lipid profiles in healthy young adults: Randomized clinical trial. J. Med. Food. 2013, 16, 223–229. [Google Scholar] [CrossRef] [PubMed]
  68. Delgado-Lista, J.; Perez-Martinez, P.; Garcia-Rios, A.; Alcala-Diaz, J.F.; Perez-Caballero, A.I.; Gomez-Delgado, F.; Fuentes, F.; Quintana-Navarro, G.; Lopez-Segura, F.; Ortiz-Morales, A.M.; et al. CORonary Diet Intervention with Olive oil and cardiovascular PREVention study (the CORDIOPREV study): Rationale, methods, and baseline characteristics: A clinical trial comparing the efficacy of a Mediterranean diet rich in olive oil versus a low-fat diet on cardiovascular disease in coronary patients. Am. Heart J. 2016, 177, 42–50. [Google Scholar] [CrossRef] [PubMed]
Ijerph 15 01453 g001 550
Figure 1. Progression of the selection of studies for the meta-analysis.
Figure 1. Progression of the selection of studies for the meta-analysis.
Ijerph 15 01453 g001
Ijerph 15 01453 g002 550
Figure 2. Forest plot for the meta-analysis of MTHFR 677 polymorphism by TT genotype, countries of non-Caucasian ethnicity with risks >1.
Figure 2. Forest plot for the meta-analysis of MTHFR 677 polymorphism by TT genotype, countries of non-Caucasian ethnicity with risks >1.
Ijerph 15 01453 g002
Ijerph 15 01453 g003 550
Figure 3. Nonlinear curve fits of MTHFR 677 polymorphism TT (left) and TT plus CT (right) genotypes of control and ischemic heart disease (IHD) groups with annual death rates from air pollution (AP death) per million population (levels: 2 = 51–100 deaths/million, 3 = 101–250 deaths, 4 = 251–400 + deaths); case = red, control = blue. Note. TT%7ct: Percentage of MTHFR 677 TT genotype in control group; TT%7ca: Percentage of TT genotype in IHD group; TTCT%7ct: Percentage of TT plus CT genotypes in control group; TTCT%7ca: Percentage of TT plus CT genotypes in IHD group.
Figure 3. Nonlinear curve fits of MTHFR 677 polymorphism TT (left) and TT plus CT (right) genotypes of control and ischemic heart disease (IHD) groups with annual death rates from air pollution (AP death) per million population (levels: 2 = 51–100 deaths/million, 3 = 101–250 deaths, 4 = 251–400 + deaths); case = red, control = blue. Note. TT%7ct: Percentage of MTHFR 677 TT genotype in control group; TT%7ca: Percentage of TT genotype in IHD group; TTCT%7ct: Percentage of TT plus CT genotypes in control group; TTCT%7ca: Percentage of TT plus CT genotypes in IHD group.
Ijerph 15 01453 g003
Table
Table 1. Schema of significant findings across studies on MTHFR 677 genotypes and risk of ischemic heart disease (IHD).
Table 1. Schema of significant findings across studies on MTHFR 677 genotypes and risk of ischemic heart disease (IHD).
MTHFR 677IHDCADMI
(n Case/n Control)(n Case/n Control)(n Case/n Control)
123 Studies (29,697/31,028)93 Studies (22,994/20,221)30 Studies (6709/10,807)
OverallRisk Type: TT, CTRisk Type: TT, CT
Protective: CCProtective: CCProtective: CC
Subgroups
Caucasian60 Studies47 Studies 13 Studies
(22,042/14,925)(13,817/10,702)(2850/4223)
Risk trend: TTRisk Type: TTNS
Protective trend: CCProtective: CC
East Asian32 Studies25 Studies7 Studies
(7112/9471)(4855/4931)(2257/4540)
Risk Type: TTRisk Type: TTRisk Type: TT
Protective: CCProtective: CC
South Asian13 Studies8 Studies5 Studies
(2475/2763)(1447/1621)(1028/1142)
Risk Type: TT, CTRisk Type: TT, CTNS
Protective: CCProtective: CC
Mixed group3 Studies2 Studies1 Study
(951/852)(882/514)(69/338)
NSNS--
Middle Eastern7 Studies6 Studies1 Study
(1265/1723)(1169/1623)(99/100)
Risk Type: TT, CTRisk Type: TT, CT--
Protective: CCProtective: CC
Hispanic2 Studies--2 Studies
(353/364) (353/364)
NS NS
African6 Studies5 Studies1 Study
(874/930)(824/648)(50/282)
Risk Type: TTRisk Type: TT, CT--
Note. NS: No statistical significance; --: Not available or not applicable; CAD: Coronary artery disease, MI: Myocardial infarction.
Table
Table 2. Pooled analyses of MTHFR 677 genotypes and risks of ischemic heart disease (IHD) by ethnicity (123 Studies).
Table 2. Pooled analyses of MTHFR 677 genotypes and risks of ischemic heart disease (IHD) by ethnicity (123 Studies).
Genotype per Ethnicity (n Studies)IHD Cases
(n = 29,697)
n (%)		Controls
(n = 31,028)
n (%)		Test of HeterogeneityTest of Association
Qp *I2Risk Ratio
(95% CI)		p
TT (123)3812(12.84)3330(10.73)238.56<0.000148.9%1.23 [1.14, 1.33]<0.0001
 Caucasian (60)1850(11.10)1479(9.91)59.380.46150.6%1.07 [1.00, 1.14]0.0528
 East Asian (32)1439(20.23)1443(15.24)60.750.001149%1.31 [1.17, 1.47]<0.0001
 South Asian (13)93(3.76)71(2.57)20.190.063540.6%1.51 [1.13, 2.01]0.0048
 Mixed Group (3)93(9.78)88(10.33)0.050.97710%0.86 [0.63, 1.17]0.33
 Middle Eastern (7)130(10.28)89(5.17)27.380.000178.1%2.62 [1.37, 5.01]0.0034
 Hispanic (2)102(28.90)106(29.12)1.500.220433.4%0.99 [0.79, 1.25]0.94
 African (6)105(12.01)54(5.81)2.030.84540%2.14 [1.56, 2.94]<0.0001
CT (123)12,777(43.02)12,815(41.30)174.490.001330.1%1.04 [1.01, 1.07]0.0028
 Caucasian (60)7402(44.41)6448(43.20)47.630.85550%1.01 [0.98, 1.04]0.48
 East Asian (32)3338(46.93)4439(46.87)53.810.006742.4%1.03 [0.98, 1.09]0.27
 South Asian (13)634(25.62)531(19.22)27.600.006356.5%1.32 [1.12, 1.57]0.0011
 Mixed Group (3)415(43.64)327(38.38)1.440.48740%1.02 [0.90, 1.15]0.78
 Middle Eastern (7)489(38.66)593(34.42)5.670.46070%1.18 [1.07, 1.30]0.0011
 Hispanic (2)165(46.74)157(43.13)1.880.170946.7%1.08 [0.92, 1.27]0.34
 African (6)334(38.22)320(34.41)9.310.097546.3%1.13 [1.00, 1.28]0.0547
CC (123)13,108(44.14)14,883(47.97)269.85<0.000154.8%0.91 [0.89, 0.94]<0.0001
 Caucasian (60)7415(44.49)6998(46.89)63.450.32257%0.98 [0.95, 1.00]0.058
 East Asian (32)2335(32.83)3589(37.89)89.99<0.000165.6%0.84 [0.77, 0.91]<0.0001
 South Asian (13)1748(70.63)2161(78.21)44.14<0.000172.8%0.91 [0.86, 0.96]0.0011
 Mixed Group (3)443(46.58)437(51.29)1.090.58040%1.02 [0.92, 1.12]0.76
 Middle Eastern (7)646(51.07)1041(60.42)16.220.012663%0.76 [0.66, 0.87]<0.0001
 Hispanic (2)86(24.36)101(27.75)0.040.84870%0.88 [0.69, 1.13]0.31
 African (6)435(49.77)556(59.78)11.880.036557.9%0.87 [0.76, 1.00]0.0501
TT + CT (123)16,589(55.86)16,145(52.03)283.27<0.000156.9%1.09 [1.06, 1.11]<0.0001
CC + CT (123)25,885(87.16)27,698(89.27)304.35<0.000159.9%0.98 [0.97, 0.99]<0.0001
T (123)10,200(34.35)9737(31.38)180.110.000532.3%1.11 [1.07, 1.14]<0.0001
C (123)19,496(65.65)21,290(68.62)172.230.001929.2%0.96 [0.94, 0.97]<0.0001
Subgroups
TT Risk > 1 (87) 118,032(60.72)22,413(72.23)
 TT2515(13.95)2351(10.49)163.42<0.000147.4%1.40 [1.28, 1.53]<0.0001
 CT7720(42.81)9243(41.24)148.57<0.000142.1%1.05 [1.02, 1.09]0.0029
 CC7797(43.24)10,819(48.27)208.20<0.000158.7%0.88 [0.85, 0.91]<0.0001
 TT + CT10,235(56.76)11,594(51.73)225.28<0.000161.8%1.13 [1.09, 1.17]<0.0001
 CC + CT15,517(86.05)20,062(89.51)267.73<0.000167.9%0.96 [0.95, 0.98]<0.0001
 T6375(35.35)6972(31.11)134.060.000735.8%1.17 [1.12, 1.22]<0.0001
 C11,657(64.65)15,440(68.89)128.350.002133%0.94 [0.92, 0.95]<0.0001
TT Risk < 1 (12) 22731(9.20)2156(6.95)
 TT322(11.79)299(13.87)7.890.72290%0.86 [0.75, 1.00]0.051
 CT1151(42.15)848(39.33)8.710.64890%1.04 [0.97, 1.12]0.23
 CC1258(46.06)1009(46.80)11.750.38286.4%1.00 [0.95, 1.06]0.90
 TT + CT1473(53.94)1147(53.20)11.240.42352.1%1.00 [0.95, 1.05]0.90
 CC + CT2409(88.21)1857(86.13)13.300.274517.3%1.02 [1.00, 1.04]0.052
 T897.5(32.86)723(33.53)6.660.82620%0.97 [0.90, 1.05]0.44
 C1833.5(67.14)1433(66.47)7.300.77420%1.02 [0.98, 1.06]0.44
TT Risk~1 (24) 38934(30.08)6459(20.82)
 TT975(10.91)680(10.53)18.020.75620%0.99 [0.90, 1.09]0.86
 CT3906(43.72)2724(42.17)16.190.84680%1.01 [0.97, 1.05]0.73
 CC4053(45.37)3055(47.30)17.430.78750%1.00 [0.96, 1.03]0.81
 TT + CT4881(54.63)3404(52.70)16.510.83250%1.00 [0.97, 1.04]0.81
 CC + CT7959(89.09)5779(89.47)18.510.72950%1.00 [0.99, 1.01]0.86
 T2928(32.77)2042(31.61)8.890.99630%1.00 [0.95, 1.05]0.95
 C6006(67.23)4417(68.39)9.190.99530%1.00 [0.98, 1.02]0.95
Note. Q = Cochran’s Q; CI = Confidence interval; p * for random effects model; 1 countries with TT risks >1: Australia, Poland, Russia, Slovakia, Sweden, Croatia, Czech Republic, Netherlands, Switzerland, France, Portugal, United Kingdom, Canada, Ireland, Mexico, Brazil, Japan, China, India, Iran, Egypt, Israel, Lebanon, Morocco, Saudi Arabia, Tunisia, Turkey; 2 countries with TT with risks <1: Hungary, Italy, Spain, Costa Rica, South Korea, Pakistan, South Africa (Indian descendants); 3 countries with TT risks varied around 1: Germany, United States, Taiwan.
Table
Table 3. Meta-prediction: Death from air pollution (AP death) on percentages of MTHFR 677 genotypes for control (Ct) and ischemic heart disease (IHD) cases (Ca) and IHD risks.
Table 3. Meta-prediction: Death from air pollution (AP death) on percentages of MTHFR 677 genotypes for control (Ct) and ischemic heart disease (IHD) cases (Ca) and IHD risks.
VariablePartition TreeTukey’s Test
AICcAP Death LevelsCountMeanSDLevels ComparedDifferenceSE DifferenceLower CIUpper CIp
TT%ct8022346.67.04/25.11.41.88.40.0011
3 and 489116.03/24.51.51.18.00.0064
4/30.61.3−2.63.70.91
TT%ca8572 and 374116.84/26.11.71.9100.0019
449169.04/34.41.70.59.30.0251
3/21.71.8−2.66.00.63
CT%ct89423432133/2112.16.316<0.0001
3 and 489437.04/29.92.05.115<0.0001
3/41.41.9−3.26.00.75
CT%ca88723439124/28.02.03.3130.0002
3 and 489446.93/26.92.12.0120.0031
4/31.11.9−3.35.60.82
CC%ct98323461182/3163.08.623<0.0001
3 and 48946112/4152.98.122<0.0001
4/30.82.8−5.77.40.95
CC%ca99723453182/4143.06.921<0.0001
3 and 48942122/38.63.21.1160.0211
3/45.52.9−1.3120.14
TT + CT%ct98323439183/2163.08.623<0.0001
3 and 48954114/2152.98.122<0.0001
3/40.82.8−5.77.40.9499
TT + CT%ca99723447184/2143.06.921<0.0001
8958123/28.63.21.1160.0211
4/35.52.9−1.3120.1410
RRTT3732291.81.72/30.60.3−0.11.20.11
3 and 4891.40.92/40.30.3−0.31.00.42
4/30.20.3−0.40.80.60
RRCT1372341.30.82/30.30.10.10.60.0031
3 and 4891.10.22/40.30.10.030.50.0244
4/30.10.1−0.10.30.65
RRTT + CT1562341.40.82/30.40.10.20.60.0007
3 and 4891.10.22/40.30.10.040.50.0200
4/30.10.1−0.10.40.40
RRCC−842 and 3740.90.23/40.20.040.10.20.0001
4490.90.23/20.10.040.030.20.0041
2/40.020.04−0.10.10.81
Note. CI = confidence interval; AICc = Akaike’s information criterion correction; AP death levels = annual death rates from air pollution levels per million (Levels 2 = 51–100, 3 = 101–250, 4 = 251–400 and greater); RR = risk ratio; ct = controls; ca = IHD case; TT%ct = percentages of MTHFR 677 TT genotype in control group; TT%ca = percentages of TT genotype in IHD cases; CT%ct = percentages of CT genotype in control group; CT%ca = percentages of CT genotype in IHD cases; CC%ct = percentages of CC wildtype in control group; CC%ca = percentages of CC wildtype in IHD cases; TTCT%ct = percentages of TT plus CT genotypes in control group; TTCT%ca = percentages of TT plus CT genotypes in IHD cases; RRTT = risk ratio of TT; RRCT = risk ratio of CT; RRTT + CT = risk ratio of TT plus CT; RRCC = risk ratio of CC.

Share and Cite

MDPI and ACS Style

Chen, Z.-F.; Young, L.; Yu, C.H.; Shiao, S.P.K. A Meta-Prediction of Methylenetetrahydrofolate-Reductase Polymorphisms and Air Pollution Increased the Risk of Ischemic Heart Diseases Worldwide. Int. J. Environ. Res. Public Health 2018, 15, 1453. https://doi.org/10.3390/ijerph15071453

AMA Style

Chen Z-F, Young L, Yu CH, Shiao SPK. A Meta-Prediction of Methylenetetrahydrofolate-Reductase Polymorphisms and Air Pollution Increased the Risk of Ischemic Heart Diseases Worldwide. International Journal of Environmental Research and Public Health. 2018; 15(7):1453. https://doi.org/10.3390/ijerph15071453

Chicago/Turabian Style

Chen, Zhao-Feng, Lufei Young, Chong Ho Yu, and S. Pamela K. Shiao. 2018. "A Meta-Prediction of Methylenetetrahydrofolate-Reductase Polymorphisms and Air Pollution Increased the Risk of Ischemic Heart Diseases Worldwide" International Journal of Environmental Research and Public Health 15, no. 7: 1453. https://doi.org/10.3390/ijerph15071453

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop