# Methylenetetrahydrofolate Reductase Gene rs1801133 and rs1801131 Polymorphisms and Essential Hypertension Risk: A Comprehensive Analysis

Authors: Yingchao Fan; Liting Wu; Wenfang Zhuang
Journal: Cardiovascular Therapeutics
 (2022)

---

## Abstract

Background. Essential hypertension (EH) is a common and multifactorial disorder that is likely to be influenced by multiple genes. The methylenetetrahydrofolate reductase (MTHFR) gene rs1801133 and rs1801131 polymorphisms influence MTHFR enzyme activity and plasma homocysteine concentration. In addition, variations in MTHFR functions likely play roles in the etiology of EH. Thus far, a large number of studies investigating the associations between the MTHFR polymorphisms and EH have provided controversial or inconclusive results. To better assess the purported relationship, we performed a comprehensive analysis of 52 published studies. Objective and Methods. Eligible studies were identified by searching the PubMed, Wanfang, and China National Knowledge Infrastructure (CNKI) databases. Odds ratios (ORs) with 95% confidence intervals (CIs) were estimated to assess the potential association between the MTHFR rs1801133 polymorphism and EH. Results. Overall, 10712 patients and 11916 controls were involved; we observed significantly increased association between the MTHFR rs1801133 polymorphism and EH risk (such as T vs. C: OR=1.38, 95% CI=1.25−1.54, P≤0.001), with similar results evident within race subgroups (such as Asian: T vs. C: OR=1.47, 95% CI=1.30−1.67, P≤0.001; compared to Chinese: T vs. C: OR=1.54, 95% CI=1.33−1.79, P≤0.001). Similar associations were also found in subgroups defined by the source of controls and genotype methods. To our regret, based on the limited studies, no association was detected for rs1801131 polymorphism. Conclusions. Our study provides evidence that the MTHFR rs1801133 null genotype may increase EH risk. Future studies with larger sample sizes are warranted to evaluate this association in more detail.


---

## Body

## 1. Introduction

Essential hypertension (EH) has a high prevalence rate worldwide and is considered to derive from complex interactions between diverse genes and environmental conditions [1, 2]. The present evidence-based treatment of EH is a critical intervention in reducing cardiovascular (CV) morbidity and mortality [3]. A contemporary meta-analysis of 123 studies with 613,815 hypertensive participants showed that, for every 10 mmHg reduction in systolic blood pressure, there is a significantly decreased risk of major CV disease events (relative risk 0.80 and 95% confidence interval (CI) 0.77–0.83), coronary heart disease (0.83 and 0.78–0.88), stroke (0.73 and 0.68–0.77), and heart failure (0.72 and 0.67–0.78) [4]. EH accounts for 90%-95% of all patients with hypertension, and approximately 20%-60% of its etiology is influenced by genetic factors [5].Most people with hypertension do not present with typical symptoms, making the condition easy to ignore and reducing the opportunities for early medical intervention. Therefore, the identification of high-risk patients is particularly meaningful, since identifying such patients with high blood pressure as early as possible allows for timely monitoring of blood pressure, regular follow-ups, and options for improving unhealthy lifestyles. In known high-risk patients, rising blood pressure can be immediately detected and treatment can be provided timeously, thus avoiding complications and, ultimately, reducing the incidence of hypertension and improving quality of life [3, 6, 7].Studies have shown that a high plasma concentration of homocysteine (Hcy) may injure the vascular endothelium, which results in hypertension and a predisposition toward atherosclerosis. Elevated Hcy levels have been identified as an independent risk factor for hypertension [8–11]. Methylenetetrahydrofolate reductase (MTHFR) plays an important role in the metabolism of Hcy. The rs1801133 polymorphism of the MTHFR gene is a C to T transition at nucleotide position 677 (C667T) in exon 4, which results in a change from alanine to valine at amino acid 222. This mutation may lead to a decrease in MTHFR activity and heat tolerance and thus metabolic damage of Hcy, which may moderately increase plasma Hcy levels [12, 13]; moreover, 1298 (A to C) also may cause a significant reduction in enzyme activity. This information suggests that MTHFR two polymorphisms may be related to EH development and susceptibility. It is therefore worthwhile to demonstrate whether there is an association between this polymorphism and EH risk, as it may provide guidance for the prevention and diagnosis of EH in the clinic. Thus far, numerous studies have reported the association between the MTHFR two common polymorphisms and EH risk. We therefore performed a comprehensive analysis of 52 different case-control studies to derive a convincing conclusion regarding this apparent association [14].


## 2. Materials and Methods

### 2.1. Identification of Eligible Studies

Searches were performed in the PubMed, Wanfang, and China National Knowledge Infrastructure (CNKI) databases (updated on Dec. 10, 2021) using the following related keywords: polymorphism/variant/mutation, hypertension/essential hypertension, and MTHFR/methylenetetrahydrofolate reductase. We included all studies that described a relationship between theMTHFR two polymorphisms and EH susceptibility. All of the included studies met the following criteria: (1) association between MTHFR two polymorphisms and EH risk; (2) case-control study; (3) each genotype frequency is shown in tables; and (4) genotype distributions of the control were consistent, with a Hardy-Weinberg equilibrium (HWE) more than 0.05. Otherwise, studies should be excluded with the following issues: (1) no control, (2) incomplete genotype frequency data, (3) duplication studies, and (4) not according with HWE in control groups.


### 2.2. Data Extraction

We collected the following information in our analysis: first author’s last name, year of publication, ethnic origin, sample size (cases/controls), study design, HWE of controls and genotype methods, and genotype frequencies in cases/controls.


### 2.3. Quality Assessment

In this meta-analysis, the quality was assessed using the Newcastle-Ottawa Scale (NOS) for cross-sectional study quality assessment. The methodological quality of each study (sampling strategy, response rate, and representativeness of the study), comparability, and outcome were checked using the NOS tool. Studies with a score of more than 7 out of 10 were considered as good quality. This cut-off point was determined after reviewing relevant meta-analyses from the literature [15–17]. Moreover, to assess the quality of each qualified articles by quality score to explore other potential sources of heterogeneity, another score of quality assessment was applied. The quality scores of the studies ranged from 0 (lowest) to 15 (highest). Studies with scores≤9 were categorized into low quality, while those with scores>9 were considered as high quality [18].


### 2.4. Statistical Analysis

The extracted data were imported into the Stata software program (version 10.0, Corporation, College Station, Texas) for analysis. Odds ratios (ORs) with 95% CIs were used to measure the strength of the association [19, 20]. The subgroup analysis stratified by race was performed first. Race was categorized as European, Asian, Mixed, Chinese, and non-Chinese (all people who are not Chinese) subtypes. The source of the control subgroups was defined based on two classifications: hospital-based (HB) and population-based (PB). For MTHFR rs1801133 and rs1801131, we investigated the relationship between genetic variants and EH risk in five different models (T-allele vs. C-allele or A-allele vs. C-allele, TT vs. CC or AA vs. CC, TC vs. CC or AC vs. CC, TT+TC vs. CC or AA+AC vs. CC, and TT vs. TC+CC or AA vs. AC+CC).The evaluation of heterogeneity within the included studies was assessed using Cochrane’sQ test (chi-square) and I2 (%) statistical analysis. A fixed effect model was applied when the effects were assumed to be homogenous (P>0.05, I2≤50%); otherwise, the random-effect model was adopted (P<0.05, I2≥50%) [21, 22]. When heterogeneity was observed, the source of the heterogeneity was explored via subgroup analysis, performed using the ethnicity, publication year, study design, and genotype methods.The presence of potential publication bias was determined using the Egger/Begg’s test and presented graphically in the form of a funnel plot [23]. In addition, the departure of MTHFR polymorphism frequencies from expectations under HWE were assessed in controls using the Pearson chi-square test [24]. Another sensitivity analysis was conducted to assess the stability of the results. Finally, the power and sample size analysis of our meta-analysis was calculated using a program called PS: Power and Sample Size Calculation (http://biostat.mc.vanderbilt.edu/wiki/Main/PowerSampleSize#Windows) [25].


### 2.5. Genotyping Methods

Genotyping to identify single nucleotide polymorphisms (SNPs) in theMTHFR gene was conducted using the following analyses: polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP), TaqMan, sequencing, PCR, amplification refractory mutation system-PCR (ARMS-PCR), and high-resolution melt (HRM) genotyping.


### 2.6. Metaregression

Random-effect metaregression analysis was performed to define the source of publication bias, with the subgroups of publication year, ethnicity, source of control, and genotype method set as independent variables and the log values considered as dependent variables [26].


### 2.7. Gene Interaction Network Analysis of theMTHFR Gene

In order to fully understand the role ofMTHFR and its potential functional partners in EH, we used the STRING online server (http://string-db.org/) to construct a MTHFR gene-gene interaction network [27, 28].


## 2.1. Identification of Eligible Studies

Searches were performed in the PubMed, Wanfang, and China National Knowledge Infrastructure (CNKI) databases (updated on Dec. 10, 2021) using the following related keywords: polymorphism/variant/mutation, hypertension/essential hypertension, and MTHFR/methylenetetrahydrofolate reductase. We included all studies that described a relationship between theMTHFR two polymorphisms and EH susceptibility. All of the included studies met the following criteria: (1) association between MTHFR two polymorphisms and EH risk; (2) case-control study; (3) each genotype frequency is shown in tables; and (4) genotype distributions of the control were consistent, with a Hardy-Weinberg equilibrium (HWE) more than 0.05. Otherwise, studies should be excluded with the following issues: (1) no control, (2) incomplete genotype frequency data, (3) duplication studies, and (4) not according with HWE in control groups.


## 2.2. Data Extraction

We collected the following information in our analysis: first author’s last name, year of publication, ethnic origin, sample size (cases/controls), study design, HWE of controls and genotype methods, and genotype frequencies in cases/controls.


## 2.3. Quality Assessment

In this meta-analysis, the quality was assessed using the Newcastle-Ottawa Scale (NOS) for cross-sectional study quality assessment. The methodological quality of each study (sampling strategy, response rate, and representativeness of the study), comparability, and outcome were checked using the NOS tool. Studies with a score of more than 7 out of 10 were considered as good quality. This cut-off point was determined after reviewing relevant meta-analyses from the literature [15–17]. Moreover, to assess the quality of each qualified articles by quality score to explore other potential sources of heterogeneity, another score of quality assessment was applied. The quality scores of the studies ranged from 0 (lowest) to 15 (highest). Studies with scores≤9 were categorized into low quality, while those with scores>9 were considered as high quality [18].


## 2.4. Statistical Analysis

The extracted data were imported into the Stata software program (version 10.0, Corporation, College Station, Texas) for analysis. Odds ratios (ORs) with 95% CIs were used to measure the strength of the association [19, 20]. The subgroup analysis stratified by race was performed first. Race was categorized as European, Asian, Mixed, Chinese, and non-Chinese (all people who are not Chinese) subtypes. The source of the control subgroups was defined based on two classifications: hospital-based (HB) and population-based (PB). For MTHFR rs1801133 and rs1801131, we investigated the relationship between genetic variants and EH risk in five different models (T-allele vs. C-allele or A-allele vs. C-allele, TT vs. CC or AA vs. CC, TC vs. CC or AC vs. CC, TT+TC vs. CC or AA+AC vs. CC, and TT vs. TC+CC or AA vs. AC+CC).The evaluation of heterogeneity within the included studies was assessed using Cochrane’sQ test (chi-square) and I2 (%) statistical analysis. A fixed effect model was applied when the effects were assumed to be homogenous (P>0.05, I2≤50%); otherwise, the random-effect model was adopted (P<0.05, I2≥50%) [21, 22]. When heterogeneity was observed, the source of the heterogeneity was explored via subgroup analysis, performed using the ethnicity, publication year, study design, and genotype methods.The presence of potential publication bias was determined using the Egger/Begg’s test and presented graphically in the form of a funnel plot [23]. In addition, the departure of MTHFR polymorphism frequencies from expectations under HWE were assessed in controls using the Pearson chi-square test [24]. Another sensitivity analysis was conducted to assess the stability of the results. Finally, the power and sample size analysis of our meta-analysis was calculated using a program called PS: Power and Sample Size Calculation (http://biostat.mc.vanderbilt.edu/wiki/Main/PowerSampleSize#Windows) [25].


## 2.5. Genotyping Methods

Genotyping to identify single nucleotide polymorphisms (SNPs) in theMTHFR gene was conducted using the following analyses: polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP), TaqMan, sequencing, PCR, amplification refractory mutation system-PCR (ARMS-PCR), and high-resolution melt (HRM) genotyping.


## 2.6. Metaregression

Random-effect metaregression analysis was performed to define the source of publication bias, with the subgroups of publication year, ethnicity, source of control, and genotype method set as independent variables and the log values considered as dependent variables [26].


## 2.7. Gene Interaction Network Analysis of theMTHFR Gene

In order to fully understand the role ofMTHFR and its potential functional partners in EH, we used the STRING online server (http://string-db.org/) to construct a MTHFR gene-gene interaction network [27, 28].


## 3. Results

### 3.1. Study Selection and Characteristics in our Meta-analysis

Using appropriate keywords (see Materials and Methods), we identified 364 articles from PubMed, 56 from CNKI, and 362 from the Wanfang database. In total, 552 articles were excluded after review of the abstract, leaving 230 articles for full article evaluation. Among them, 39 articles featured systematic analysis/meta-analysis/review; 24 covered only case groups; 17 articles were duplicates of other papers; 35 had no original numbers for case/control groups and presented only total numbers; 28 articles focused on H-type hypertension; 4 were related to aortic hypertension; and 27 covered hypertension in pregnancy (Figure1). After exclusion of the above studies by full article review, we were left with 55 articles covering 60 case-control studies, of which 8 case-control studies were not consistent with HWE and were excluded. In total, 49 case-control studies about rs1801133 and 3 case-control studies about rs1801131 were included in our current analysis. All essential information is listed in Table 1, including first author, publishing year, race, the numbers of cases and controls, HWE, genotype numbers in cases/controls, study design, and genotype methods. Our study comprised 9 European, 36 Asian, and 4 Mixed-race case-control studies. The T frequency was 43.59% in the Asian group, 35.2% among Europeans, and 47.7% in the Mixed-race group, which indicated that the European group had lower frequency of T-allele than the Asian and Mixed groups (P<0.05). The distribution of genotypes in all the controls was in agreement with HWE. In addition, we confirmed the minor allele frequency (MAF) reported for the seven main worldwide populations in the 1000 Genomes Browser [29] (https://www.ncbi.nlm.nih.gov/snp/rs1801133), as follows: Global (0.335), European (0.345), East Asian (0.328), South Asian (0.167), African (0.123), African American (0.125), and Asian (0.265) (Figure 2(a)). When we examined the frequency of T- and C-alleles both in the case and control groups, our analyses showed similar allele frequencies in each group (Figure 2(b)). Finally, we used The Cancer Genome Atlas (TCGA) database to search for trends in the frequency of rs1801133 polymorphisms. Our results indicated that the TT (AA) frequency was relatively low compared to other genotypes (Figure 2(c)). This polymorphism is associated with the coronary artery, rather than the aorta artery left ventricle and the tibial artery (https://www.gtexportal.org/home/) (Figure 2(d)). We also showed the corresponding information about rs1801131 polymorphism (Figures 3(a)–3(d)).Figure 1
Flowchart illustrating the search strategy used to identify association studies forMTHFR gene rs1801133 polymorphisms and EH risk.Table 1
Characteristics of studies ofMTHFR gene two common polymorphisms and essential hypertension risk included in our meta-analysis.
AuthorYearCountryEthnicityCaseControlCaseControlSOCHWEGenotypeScoreNOSTTTCCCTTTCCCGao1999ChinaAsian127170156844248462PB0.6PCR-RFLP77Wang2002ChinaAsian1054637511792314HB0.935PCR-RFLP76Sun2003ChinaAsian5517327226186986HB0.456ARMS-PCR66Liu2005ChinaAsian100100264529195031HB0.883PCR77Li2006ChinaAsian263026182721HB0.226PCR-RFLP76Hu2006ChinaAsian110115163955124261PB0.249PCR107Tang2007ChinaAsian2521952093139651138HB0.629PCR-RFLP86Lin2008ChinaAsian501234271964473HB0.847PCR76Zhang2012ChinaAsian189165853128741117PB0.175PCR-RFLP108Yin2012ChinaAsian6706826835824451309322PB0.096PCR-RFLP128Yao2013ChinaAsian150150496932226761HB0.607PCR-RFLP87Cai2014ChinaAsian200200629939508961PB0.129PCR107Liu2019ChinaAsian9341075295439200356505214HB0.151TaqMan107Ghogomu2016ChinaAsian4150142430545HB0.709PCR76Dai2015ChinaAsian114104325725264929PB0.562PCR-RFLP108Wen2015ChinaAsian17463476534585291258PB0.837PCR107Wu2016ChinaAsian123120113973104070PB0.223Gene Chip97Fan2016ChinaAsian2144947510237141234119HB0.493TaqMan96Zhao2017ChinaAsian200200479954299180HB0.705PCR-RFLP87Sui2020ChinaAsian102109315219224938HB0.397HRM76Yu2020ChinaAsian1371285478554281HB0.877TaqMan86Yu2020ChinaAsian163160317953277657HB0.845TaqMan87Nong2019ChinaAsian122110495815165935PB0.267PCR-RFLP98Zhao2016ChinaAsian100100235027154540HB0.689PCR-RFLP67Tian2018ChinaAsian743718203373167148370200HB0.333HRM97Sui2019ChinaAsian11373445019104122HB0.186PCR-RFLP77Zhang2020ChinaAsian74153816431326492268178PB0.603TaqMan128Nishio1996JapanAsian47825261694429HB0.201PCR57Nakata1998JapanAsian173184199163368365HB0.309PCR87Lwin2006JapanAsian1162191958393811764PB0.215PCR68Hui2007JapanAsian261271491298344123104HB0.454PCR-SSCP66Markan2007IndiaAsian153133840105028105PB0.174PCR-RFLP118Nassereddine2015MoroccoAsian10110214404734554PB0.074PCR-RFLP107Candrasatria2020IndonesiaAsian213202673134342157PB0.92TaqMan97Arina2019IndonesiaAsian53535163201043HB0.448PCR-RFLP67Kahleova2002CzechEuropean164173275582186986PB0.457ARMS-PCR88Rodriguez-Esparragon2003SpainEuropean23221534115832010095PB0.386PCR118Heux2004New ZealandEuropean247249351258725119105HB0.298PCR-RFLP77Tylicki2005Austria/PolandEuropean9090113940103842PB0.752PCR117Ilhan2008TurkeyEuropean7810010323622672HB0.844Real-time PCR47Deshmukh2009USAEuropean4211841622184852PB0.221Sequencing108Ng2009AustraliaEuropean388010141483240PB0.67PCR78Fowdar2012AustraliaEuropean3773933317417035183175PB0.186PCR-RFLP117Bayramoglu2013TurkeyEuropean1259922386553856HB0.654PCR87Fridman2013ArgentinaMixed7515064029156471PB0.917PCR-RFLP98Perez-Razo2015MexicoMixed3733918717411210120090HB0.637TaqMan96Perez-Razo2015MexicoMixed1991996798345610835HB0.168TaqMan87Vazquez-Alaniz2014MexicoMixed194194399362439754HB0.964TaqMan87Yadav2021IndiaAsian20725410581391165178PB0.116PCR-RFLP119HB: hospital-based; PB: population-based; SOC; source of control; PCR-RFLP: polymerase chain reaction followed by restriction fragment length polymorphism; HRM: high-resolution melt; ARMS-PCR: amplification refractory mutation system-PCR; HWE: Hardy-Weinberg equilibrium of control group.Figure 2
(a) The MAF of minor allele (mutant allele) forMTHFR gene rs1801133 polymorphism from the 1000 Genomes online database. (b) The frequency about T-allele or C-allele both in the case and control groups. (c) The distribution of each genotype from online GTEx Portal (https://www.gtexportal.org/home/). (d) The risk frequency of rs1801133 polymorphism to several diseases from TCGA database.
(a)(b)(c)(d)Figure 3
(a) The MAF of minor allele (mutant allele) forMTHFR gene rs1801131 polymorphism from the 1000 Genomes online database. (b) The frequency about T-allele or C-allele both in the case and control groups. (c) The distribution of each genotype from online GTEx Portal (https://www.gtexportal.org/home/). (d) The risk frequency of rs1801131 polymorphism to several diseases from TCGA database.
(a)(b)(c)(d)


### 3.2. Quantitative Data Synthesis

Table2 shows the summary of odds ratios of MTHFR, based on 10,712 EH cases and 11,916 matched controls. We observed an increased association between the MTHFR rs1801133 polymorphism and EH in total population groups (for example, T-allele vs. C-allele: OR=1.38, 95% CI=1.25−1.54, Ph<0.001, P<0.001, I2=81.3%, Figure 4). Similar trends were observed in the analyses of ethnic subgroups (T-allele vs. C-allele: OR=1.47, 95% CI=1.30−1.67, Ph<0.001, P<0.001, I2=83.2% for Asians, Figure 4; OR=1.28, 95% CI=1.05−1.57, Ph=0.004, P=0.014, I2=64.2% for Europeans, Figure 4; OR=1.54, 95% CI=1.33−1.79, Ph<0.001, P<0.001, I2=85.6% for Chinese, Figure 5; and OR=1.20, 95% CI=1.06−1.37, Ph<0.001, P=0.004, I2=66.3% for non-Chinese). In order to analyze the source of controls and determine the source of heterogeneity, the ratios were calculated for the HB and PB subgroups. The results showed significantly increased relationships in these groups (T-allele vs. C-allele: OR=1.49, 95% CI=1.28−1.75, Ph<0.001, P<0.001, I2=85.4% for HB and OR=1.29, 95% CI=1.13−1.47, Ph<0.001, P<0.001, I2=73.3% for PB) (Figure 6). As different methods for detecting this polymorphism were applied in all of the included studies, we considered whether positive results were associated with particular genotyping methods. Our analyses revealed some significant findings, such as PCR (T-allele vs. C-allele: OR=1.51, 95% CI=1.14−2.01, Ph<0.001, P=0.004, I2=86.1%), PCR-RFLP (T-allele vs. C-allele: OR=1.47, 95% CI=1.29−1.68, Ph<0.001, P<0.001, I2=64.3%, Figure 7), and HRM analyses (T-allele vs. C-allele: OR=1.32, 95% CI=1.15−1.51, Ph<0.001, P<0.001, I2=47.5%, Figure 8). To our regret, no positive associations were observed for rs1801131 polymorphism (Table 2).Table 2
Stratified analyses ofMTHFR gene rs1801133 polymorphism and essential hypertension risk.
VariablesNCase/controlT vs. CTC vs. CCTT vs. CCTT+TC vs. CCTT vs. TC+CCOR (95% CI) PhPI2OR (95% CI) PhPOR (95% CI) PhPrs1801133Total499613/107131.38 (1.25-1.54) 0.000 0.000 81.3%1.26 (1.13-1.41) 0.000 0.000 58.7%1.80 (1.47-2.20) 0.000 0.000 75.6%1.41 (1.25-1.59) 0.000 0.000 71.1%1.57 (1.33-1.85) 0.000 0.000 72.8%EthnicityAsian367379/82621.47 (1.30-1.67) 0.000 0.000 83.2%1.36 (1.19-1.55) 0.000 0.000 61.0%1.99 (1.56-2.55) 0.000 0.000 77.6%1.53 (1.32-1.78) 0.000 0.000 72.9%1.67 (1.36-2.05) 0.000 0.000 76.2%European91393/15171.28 (1.05-1.57) 0.004 0.014 64.2%1.11 (0.94-1.29) 0.2225 0.210 24.6%1.73 (1.15-2.59) 0.024 0.009 54.7%1.24 (0.99-1.54) 0.057 0.054 47.2%1.63 (1.13-2.35) 0.045 0.009 49.4%Mixed4841/9340.93 (0.82-1.06) 0.219 0.313 32.2%0.86 (0.69-1.08) 0.156 0.205 42.5%1.50 (0.88-2.56) 0.042 0.133 63.4%1.37 (0.90-2.09) 0.017 0.139 70.6%1.06 (0.92-1.24) 0.118 0.397 48.9%Chinese276055/67621.54 (1.33-1.79) 0.000 0.000 85.61.41 (1.19-1.66) 0.000 0.000 67.4%2.17 (1.65-2.85) 0.000 0.000 80.1%1.63 (1.36-1.95) 0.000 0.000 77.2%1.77 (1.41-2.25) 0.000 0.000 78.7%Non-Chinese223558/39511.20 (1.06-1.37) 0.000 0.004 66.3%1.12 (0.98-1.28) 0.048 0.091 36.0%1.35 (1.03-1.77) 0.000 0.030 58.9%1.19 (1.03-1.38) 0.002 0.018 52.6%1.28 (1.02-1.61) 0.002 0.035 53.3%SOCHB275235/57461.49 (1.28-1.75) 0.000 0.000 85.4%1.38 (1.16-1.64) 0.000 0.000 67.1%1.96 (1.47-2.67) 0.000 0.000 79.5%1.57 (1.28-1.91) 0.000 0.000 78.1%1.60 (1.29-1.99) 0.000 0.000 72.5%PB224378/49671.29 (1.13-1.47) 0.000 0.000 73.3%1.17 (1.03-1.33) 0.019 0.019 42.7%1.64 (1.24-2.17) 0.000 0.001 69.0%1.28 (1.12-1.48) 0.001 0.000 55.7%1.51 (1.16-1.97) 0.000 0.002 70.8%Genotype methodsPCR131496/21911.51 (1.14-2.01) 0.000 0.004 86.1%1.31 (0.97-1.75) 0.000 0.075 69.4%1.93 (1.18-3.17) 0.000 0.009 78.9%1.51 (1.10-2.07) 0.000 0.010 76.6%1.71 (1.07-2.73) 0.000 0.024 81.1%PCR-RFLP193381/33591.47 (1.29-1.68) 0.000 0.000 64.3%1.38 (1.24-1.54) 0.547 0.000 0.0%2.14 (1.58-2.91) 0.000 0.000 60.3%1.54 (1.33-1.78) 0.027 0.000 42.3%1.74 (1.34-2.25) 0.003 0.000 54.2%ARMS-PCR2219/3462.39 (0.55-10.45) 0.000 0.244 96.3%1.85 (0.87-2.47) 0.002 0.468 89.9%5.64 (0.43-73.32) 0.000 0.186 94.3%2.72 (0.34-21.96) 0.000 0.348 94.3%3.72 (0.78-17.71) 0.001 0.099 90.5%TaqMan93168/33811.06 (0.93-1.21) 0.006 0.358 63.0%1.01 (0.82-1.24) 0.008 0.941 61.1%1.05 (0.85-1.30) 0.094 0.639 41.0%1.04 (0.85-1.27) 0.005 0.691 63.9%1.07 (0.95-1.20) 0.230 0.247 24.0%HRM2845/8271.32 (1.15-1.51) 0.167 0.000 47.5%1.29 (1.02-1.64) 0.125 0.031 57.6%1.76 (1.34-2.32) 0.202 0.000 38.5%1.43 (1.15-1.78) 0.104 0.002 62.2%1.48 (1.18-1.86) 0.609 0.001 0.0%Other4504/6091.25 (0.79-1.97) 0.002 0.346 80.0%1.25 (0.82-1.91) 0.082 0.291 55.2%1.45 (0.63-3.32) 0.030 0.386 66.6%1.29 (0.78-2.15) 0.021 0.319 72.8%1.29 (0.65-2.56) 0.081 0.468 55.5%A vs. CAC vs. CCAA vs. CCAA+AC vs. CCAA vs. AC+CCOR (95% CI) PhP I-squaredOR (95% CI) PhPOR (95% CI) PhPrs180113131099/12030.93 (0.80-1.09) 0.710 0.384 0.0%0.97 (0.59-1.59) 0.990 0.901 0.0%0.90 (0.55-1.46) 0.912 0.667 0.0%0.92 (0.57-1.48) 0.979 0.730 0.0%0.92 (0.77-1.11) 0.632 0.382 0.0%Ph: value of Q test for heterogeneity test; P: Z test for the statistical significance of the OR.Figure 4
T-allele frequencies for theMTHFR gene rs1801133 polymorphism among cases/controls stratified by subgroups in T-allele vs. C-allele model in regular ethnic subgroup.Figure 5
T-allele frequencies for theMTHFR gene rs1801133 polymorphism among cases/controls stratified by subgroups in T-allele vs. C-allele model in Chinese/non-Chinese subgroup.Figure 6
T-allele frequencies for theMTHFR gene rs1801133 polymorphism among cases/controls stratified by subgroups in T-allele vs. C-allele model in source of control subgroup.Figure 7
T-allele frequencies for theMTHFR gene rs1801133 polymorphism among cases/controls stratified by subgroups in T-allele vs. C-allele model in genotype method subgroup by random-effect model.Figure 8
T-allele frequencies for theMTHFR gene rs1801133 polymorphism among cases/controls stratified by subgroups in T-allele vs. C-allele model in genotype method subgroup by fixed effect model.


### 3.3. Publication Bias and Sensitivity Analysis

Begg’s funnel plot and Egger’s test were performed to determine the publication bias of the included studies. Significant obvious evidence of publication bias was detected in five genetic model analyses (such as Figures9(a) and 9(b) regarding T-allele vs. C-allele) (Table 3).Figure 9
(a) Begg’s funnel plot for publication bias test (T-allele vs. C-allele). (b) Egger’s publication bias plot (T-allele vs. C-allele). (c) Sensitivity analysis (T-allele vs. C-allele).
(a)(b)(c)Table 3
Publication bias tests (Begg’s funnel plot and Egger’s test for publication bias test) forMTHFR gene rs1801133 polymorphism and essential hypertension risk.
Egger’s testBegg’s testGenetic typeCoefficientStandard errortP value95% CI of interceptzP valueT vs. C2.4680.7193.430.001(1.021-3.915)2.370.018TC vs. CC1.7260.4893.530.001(0.743-2.709)2.940.003TT vs. CC0.8310.2892.870.006(0.248-1.413)2.080.038TT+TC vs. CC1.8030.5193.470.001(0.757-2.849)2.960.003TT vs. TC+CC0.8580.3262.630.012(0.201-1.514)1.590.111To remove studies which may influence the power and stability of the current meta-analysis, sensitivity analysis was applied, but no sensitive case-control studies were found for this SNP among the above five models (such as Figure9(c) regarding T-allele vs. C-allele).


### 3.4. Metaregression

The metaregression analysis indicated that there is a significant relationship for the allele model (T-allele vs. C-allele) with respect to ethnicity, source of control, and genotype methods, with a regression coefficient of 0.001, 0.004, 0.010, and 0.002, respectively. There was no association with publication year, which suggests that the heterogeneity from the rs1801133 polymorphism in EH may be due to the ethnicity, source of control, and genotype method subgroups (Figures10(a)–10(e)).Figure 10
Random-effect metaregression of log odds ratio versus publication year (a), regular ethnicity (b), Chinese/non-Chinese (c), source of control (d), and genotype methods (e), respectively, in EH.
(a)(b)(c)(d)(e)


### 3.5. Gene-Gene Network Diagram and Interactions

Our analysis using the STRING online server indicated thatMTHFR interacts with numerous genes. The ten most significant genes from the network of gene-gene interactions have been listed in Figures 11(a) and 11(b). These ten genes are methionine (MTR), thymidylate synthase (TYMS), C-1-tetrahydrofolate synthase (MTHFD1), serine hydroxymethyltransferase 1 (SHMT1), serine hydroxymethyltransferase 2 (SHMT2), bifunctional methylenetetrahydrofolate dehydrogenase (MTHFD2), probable bifunctional methylenetetrahydrofolate dehydrogenase (MTHFD2L), aminomethyltransferase (AMT), and methionine synthase reductase (MTRR).Figure 11
Human MTHFR interaction network with other genes obtained from STRING server. At least 10 genes have been indicated to correlate with MTHFR gene. (a) The gene-gene interaction. (b) The detail of relative ten core genes.
(a)(b)


## 3.1. Study Selection and Characteristics in our Meta-analysis

Using appropriate keywords (see Materials and Methods), we identified 364 articles from PubMed, 56 from CNKI, and 362 from the Wanfang database. In total, 552 articles were excluded after review of the abstract, leaving 230 articles for full article evaluation. Among them, 39 articles featured systematic analysis/meta-analysis/review; 24 covered only case groups; 17 articles were duplicates of other papers; 35 had no original numbers for case/control groups and presented only total numbers; 28 articles focused on H-type hypertension; 4 were related to aortic hypertension; and 27 covered hypertension in pregnancy (Figure1). After exclusion of the above studies by full article review, we were left with 55 articles covering 60 case-control studies, of which 8 case-control studies were not consistent with HWE and were excluded. In total, 49 case-control studies about rs1801133 and 3 case-control studies about rs1801131 were included in our current analysis. All essential information is listed in Table 1, including first author, publishing year, race, the numbers of cases and controls, HWE, genotype numbers in cases/controls, study design, and genotype methods. Our study comprised 9 European, 36 Asian, and 4 Mixed-race case-control studies. The T frequency was 43.59% in the Asian group, 35.2% among Europeans, and 47.7% in the Mixed-race group, which indicated that the European group had lower frequency of T-allele than the Asian and Mixed groups (P<0.05). The distribution of genotypes in all the controls was in agreement with HWE. In addition, we confirmed the minor allele frequency (MAF) reported for the seven main worldwide populations in the 1000 Genomes Browser [29] (https://www.ncbi.nlm.nih.gov/snp/rs1801133), as follows: Global (0.335), European (0.345), East Asian (0.328), South Asian (0.167), African (0.123), African American (0.125), and Asian (0.265) (Figure 2(a)). When we examined the frequency of T- and C-alleles both in the case and control groups, our analyses showed similar allele frequencies in each group (Figure 2(b)). Finally, we used The Cancer Genome Atlas (TCGA) database to search for trends in the frequency of rs1801133 polymorphisms. Our results indicated that the TT (AA) frequency was relatively low compared to other genotypes (Figure 2(c)). This polymorphism is associated with the coronary artery, rather than the aorta artery left ventricle and the tibial artery (https://www.gtexportal.org/home/) (Figure 2(d)). We also showed the corresponding information about rs1801131 polymorphism (Figures 3(a)–3(d)).Figure 1
Flowchart illustrating the search strategy used to identify association studies forMTHFR gene rs1801133 polymorphisms and EH risk.Table 1
Characteristics of studies ofMTHFR gene two common polymorphisms and essential hypertension risk included in our meta-analysis.
AuthorYearCountryEthnicityCaseControlCaseControlSOCHWEGenotypeScoreNOSTTTCCCTTTCCCGao1999ChinaAsian127170156844248462PB0.6PCR-RFLP77Wang2002ChinaAsian1054637511792314HB0.935PCR-RFLP76Sun2003ChinaAsian5517327226186986HB0.456ARMS-PCR66Liu2005ChinaAsian100100264529195031HB0.883PCR77Li2006ChinaAsian263026182721HB0.226PCR-RFLP76Hu2006ChinaAsian110115163955124261PB0.249PCR107Tang2007ChinaAsian2521952093139651138HB0.629PCR-RFLP86Lin2008ChinaAsian501234271964473HB0.847PCR76Zhang2012ChinaAsian189165853128741117PB0.175PCR-RFLP108Yin2012ChinaAsian6706826835824451309322PB0.096PCR-RFLP128Yao2013ChinaAsian150150496932226761HB0.607PCR-RFLP87Cai2014ChinaAsian200200629939508961PB0.129PCR107Liu2019ChinaAsian9341075295439200356505214HB0.151TaqMan107Ghogomu2016ChinaAsian4150142430545HB0.709PCR76Dai2015ChinaAsian114104325725264929PB0.562PCR-RFLP108Wen2015ChinaAsian17463476534585291258PB0.837PCR107Wu2016ChinaAsian123120113973104070PB0.223Gene Chip97Fan2016ChinaAsian2144947510237141234119HB0.493TaqMan96Zhao2017ChinaAsian200200479954299180HB0.705PCR-RFLP87Sui2020ChinaAsian102109315219224938HB0.397HRM76Yu2020ChinaAsian1371285478554281HB0.877TaqMan86Yu2020ChinaAsian163160317953277657HB0.845TaqMan87Nong2019ChinaAsian122110495815165935PB0.267PCR-RFLP98Zhao2016ChinaAsian100100235027154540HB0.689PCR-RFLP67Tian2018ChinaAsian743718203373167148370200HB0.333HRM97Sui2019ChinaAsian11373445019104122HB0.186PCR-RFLP77Zhang2020ChinaAsian74153816431326492268178PB0.603TaqMan128Nishio1996JapanAsian47825261694429HB0.201PCR57Nakata1998JapanAsian173184199163368365HB0.309PCR87Lwin2006JapanAsian1162191958393811764PB0.215PCR68Hui2007JapanAsian261271491298344123104HB0.454PCR-SSCP66Markan2007IndiaAsian153133840105028105PB0.174PCR-RFLP118Nassereddine2015MoroccoAsian10110214404734554PB0.074PCR-RFLP107Candrasatria2020IndonesiaAsian213202673134342157PB0.92TaqMan97Arina2019IndonesiaAsian53535163201043HB0.448PCR-RFLP67Kahleova2002CzechEuropean164173275582186986PB0.457ARMS-PCR88Rodriguez-Esparragon2003SpainEuropean23221534115832010095PB0.386PCR118Heux2004New ZealandEuropean247249351258725119105HB0.298PCR-RFLP77Tylicki2005Austria/PolandEuropean9090113940103842PB0.752PCR117Ilhan2008TurkeyEuropean7810010323622672HB0.844Real-time PCR47Deshmukh2009USAEuropean4211841622184852PB0.221Sequencing108Ng2009AustraliaEuropean388010141483240PB0.67PCR78Fowdar2012AustraliaEuropean3773933317417035183175PB0.186PCR-RFLP117Bayramoglu2013TurkeyEuropean1259922386553856HB0.654PCR87Fridman2013ArgentinaMixed7515064029156471PB0.917PCR-RFLP98Perez-Razo2015MexicoMixed3733918717411210120090HB0.637TaqMan96Perez-Razo2015MexicoMixed1991996798345610835HB0.168TaqMan87Vazquez-Alaniz2014MexicoMixed194194399362439754HB0.964TaqMan87Yadav2021IndiaAsian20725410581391165178PB0.116PCR-RFLP119HB: hospital-based; PB: population-based; SOC; source of control; PCR-RFLP: polymerase chain reaction followed by restriction fragment length polymorphism; HRM: high-resolution melt; ARMS-PCR: amplification refractory mutation system-PCR; HWE: Hardy-Weinberg equilibrium of control group.Figure 2
(a) The MAF of minor allele (mutant allele) forMTHFR gene rs1801133 polymorphism from the 1000 Genomes online database. (b) The frequency about T-allele or C-allele both in the case and control groups. (c) The distribution of each genotype from online GTEx Portal (https://www.gtexportal.org/home/). (d) The risk frequency of rs1801133 polymorphism to several diseases from TCGA database.
(a)(b)(c)(d)Figure 3
(a) The MAF of minor allele (mutant allele) forMTHFR gene rs1801131 polymorphism from the 1000 Genomes online database. (b) The frequency about T-allele or C-allele both in the case and control groups. (c) The distribution of each genotype from online GTEx Portal (https://www.gtexportal.org/home/). (d) The risk frequency of rs1801131 polymorphism to several diseases from TCGA database.
(a)(b)(c)(d)


## 3.2. Quantitative Data Synthesis

Table2 shows the summary of odds ratios of MTHFR, based on 10,712 EH cases and 11,916 matched controls. We observed an increased association between the MTHFR rs1801133 polymorphism and EH in total population groups (for example, T-allele vs. C-allele: OR=1.38, 95% CI=1.25−1.54, Ph<0.001, P<0.001, I2=81.3%, Figure 4). Similar trends were observed in the analyses of ethnic subgroups (T-allele vs. C-allele: OR=1.47, 95% CI=1.30−1.67, Ph<0.001, P<0.001, I2=83.2% for Asians, Figure 4; OR=1.28, 95% CI=1.05−1.57, Ph=0.004, P=0.014, I2=64.2% for Europeans, Figure 4; OR=1.54, 95% CI=1.33−1.79, Ph<0.001, P<0.001, I2=85.6% for Chinese, Figure 5; and OR=1.20, 95% CI=1.06−1.37, Ph<0.001, P=0.004, I2=66.3% for non-Chinese). In order to analyze the source of controls and determine the source of heterogeneity, the ratios were calculated for the HB and PB subgroups. The results showed significantly increased relationships in these groups (T-allele vs. C-allele: OR=1.49, 95% CI=1.28−1.75, Ph<0.001, P<0.001, I2=85.4% for HB and OR=1.29, 95% CI=1.13−1.47, Ph<0.001, P<0.001, I2=73.3% for PB) (Figure 6). As different methods for detecting this polymorphism were applied in all of the included studies, we considered whether positive results were associated with particular genotyping methods. Our analyses revealed some significant findings, such as PCR (T-allele vs. C-allele: OR=1.51, 95% CI=1.14−2.01, Ph<0.001, P=0.004, I2=86.1%), PCR-RFLP (T-allele vs. C-allele: OR=1.47, 95% CI=1.29−1.68, Ph<0.001, P<0.001, I2=64.3%, Figure 7), and HRM analyses (T-allele vs. C-allele: OR=1.32, 95% CI=1.15−1.51, Ph<0.001, P<0.001, I2=47.5%, Figure 8). To our regret, no positive associations were observed for rs1801131 polymorphism (Table 2).Table 2
Stratified analyses ofMTHFR gene rs1801133 polymorphism and essential hypertension risk.
VariablesNCase/controlT vs. CTC vs. CCTT vs. CCTT+TC vs. CCTT vs. TC+CCOR (95% CI) PhPI2OR (95% CI) PhPOR (95% CI) PhPrs1801133Total499613/107131.38 (1.25-1.54) 0.000 0.000 81.3%1.26 (1.13-1.41) 0.000 0.000 58.7%1.80 (1.47-2.20) 0.000 0.000 75.6%1.41 (1.25-1.59) 0.000 0.000 71.1%1.57 (1.33-1.85) 0.000 0.000 72.8%EthnicityAsian367379/82621.47 (1.30-1.67) 0.000 0.000 83.2%1.36 (1.19-1.55) 0.000 0.000 61.0%1.99 (1.56-2.55) 0.000 0.000 77.6%1.53 (1.32-1.78) 0.000 0.000 72.9%1.67 (1.36-2.05) 0.000 0.000 76.2%European91393/15171.28 (1.05-1.57) 0.004 0.014 64.2%1.11 (0.94-1.29) 0.2225 0.210 24.6%1.73 (1.15-2.59) 0.024 0.009 54.7%1.24 (0.99-1.54) 0.057 0.054 47.2%1.63 (1.13-2.35) 0.045 0.009 49.4%Mixed4841/9340.93 (0.82-1.06) 0.219 0.313 32.2%0.86 (0.69-1.08) 0.156 0.205 42.5%1.50 (0.88-2.56) 0.042 0.133 63.4%1.37 (0.90-2.09) 0.017 0.139 70.6%1.06 (0.92-1.24) 0.118 0.397 48.9%Chinese276055/67621.54 (1.33-1.79) 0.000 0.000 85.61.41 (1.19-1.66) 0.000 0.000 67.4%2.17 (1.65-2.85) 0.000 0.000 80.1%1.63 (1.36-1.95) 0.000 0.000 77.2%1.77 (1.41-2.25) 0.000 0.000 78.7%Non-Chinese223558/39511.20 (1.06-1.37) 0.000 0.004 66.3%1.12 (0.98-1.28) 0.048 0.091 36.0%1.35 (1.03-1.77) 0.000 0.030 58.9%1.19 (1.03-1.38) 0.002 0.018 52.6%1.28 (1.02-1.61) 0.002 0.035 53.3%SOCHB275235/57461.49 (1.28-1.75) 0.000 0.000 85.4%1.38 (1.16-1.64) 0.000 0.000 67.1%1.96 (1.47-2.67) 0.000 0.000 79.5%1.57 (1.28-1.91) 0.000 0.000 78.1%1.60 (1.29-1.99) 0.000 0.000 72.5%PB224378/49671.29 (1.13-1.47) 0.000 0.000 73.3%1.17 (1.03-1.33) 0.019 0.019 42.7%1.64 (1.24-2.17) 0.000 0.001 69.0%1.28 (1.12-1.48) 0.001 0.000 55.7%1.51 (1.16-1.97) 0.000 0.002 70.8%Genotype methodsPCR131496/21911.51 (1.14-2.01) 0.000 0.004 86.1%1.31 (0.97-1.75) 0.000 0.075 69.4%1.93 (1.18-3.17) 0.000 0.009 78.9%1.51 (1.10-2.07) 0.000 0.010 76.6%1.71 (1.07-2.73) 0.000 0.024 81.1%PCR-RFLP193381/33591.47 (1.29-1.68) 0.000 0.000 64.3%1.38 (1.24-1.54) 0.547 0.000 0.0%2.14 (1.58-2.91) 0.000 0.000 60.3%1.54 (1.33-1.78) 0.027 0.000 42.3%1.74 (1.34-2.25) 0.003 0.000 54.2%ARMS-PCR2219/3462.39 (0.55-10.45) 0.000 0.244 96.3%1.85 (0.87-2.47) 0.002 0.468 89.9%5.64 (0.43-73.32) 0.000 0.186 94.3%2.72 (0.34-21.96) 0.000 0.348 94.3%3.72 (0.78-17.71) 0.001 0.099 90.5%TaqMan93168/33811.06 (0.93-1.21) 0.006 0.358 63.0%1.01 (0.82-1.24) 0.008 0.941 61.1%1.05 (0.85-1.30) 0.094 0.639 41.0%1.04 (0.85-1.27) 0.005 0.691 63.9%1.07 (0.95-1.20) 0.230 0.247 24.0%HRM2845/8271.32 (1.15-1.51) 0.167 0.000 47.5%1.29 (1.02-1.64) 0.125 0.031 57.6%1.76 (1.34-2.32) 0.202 0.000 38.5%1.43 (1.15-1.78) 0.104 0.002 62.2%1.48 (1.18-1.86) 0.609 0.001 0.0%Other4504/6091.25 (0.79-1.97) 0.002 0.346 80.0%1.25 (0.82-1.91) 0.082 0.291 55.2%1.45 (0.63-3.32) 0.030 0.386 66.6%1.29 (0.78-2.15) 0.021 0.319 72.8%1.29 (0.65-2.56) 0.081 0.468 55.5%A vs. CAC vs. CCAA vs. CCAA+AC vs. CCAA vs. AC+CCOR (95% CI) PhP I-squaredOR (95% CI) PhPOR (95% CI) PhPrs180113131099/12030.93 (0.80-1.09) 0.710 0.384 0.0%0.97 (0.59-1.59) 0.990 0.901 0.0%0.90 (0.55-1.46) 0.912 0.667 0.0%0.92 (0.57-1.48) 0.979 0.730 0.0%0.92 (0.77-1.11) 0.632 0.382 0.0%Ph: value of Q test for heterogeneity test; P: Z test for the statistical significance of the OR.Figure 4
T-allele frequencies for theMTHFR gene rs1801133 polymorphism among cases/controls stratified by subgroups in T-allele vs. C-allele model in regular ethnic subgroup.Figure 5
T-allele frequencies for theMTHFR gene rs1801133 polymorphism among cases/controls stratified by subgroups in T-allele vs. C-allele model in Chinese/non-Chinese subgroup.Figure 6
T-allele frequencies for theMTHFR gene rs1801133 polymorphism among cases/controls stratified by subgroups in T-allele vs. C-allele model in source of control subgroup.Figure 7
T-allele frequencies for theMTHFR gene rs1801133 polymorphism among cases/controls stratified by subgroups in T-allele vs. C-allele model in genotype method subgroup by random-effect model.Figure 8
T-allele frequencies for theMTHFR gene rs1801133 polymorphism among cases/controls stratified by subgroups in T-allele vs. C-allele model in genotype method subgroup by fixed effect model.


## 3.3. Publication Bias and Sensitivity Analysis

Begg’s funnel plot and Egger’s test were performed to determine the publication bias of the included studies. Significant obvious evidence of publication bias was detected in five genetic model analyses (such as Figures9(a) and 9(b) regarding T-allele vs. C-allele) (Table 3).Figure 9
(a) Begg’s funnel plot for publication bias test (T-allele vs. C-allele). (b) Egger’s publication bias plot (T-allele vs. C-allele). (c) Sensitivity analysis (T-allele vs. C-allele).
(a)(b)(c)Table 3
Publication bias tests (Begg’s funnel plot and Egger’s test for publication bias test) forMTHFR gene rs1801133 polymorphism and essential hypertension risk.
Egger’s testBegg’s testGenetic typeCoefficientStandard errortP value95% CI of interceptzP valueT vs. C2.4680.7193.430.001(1.021-3.915)2.370.018TC vs. CC1.7260.4893.530.001(0.743-2.709)2.940.003TT vs. CC0.8310.2892.870.006(0.248-1.413)2.080.038TT+TC vs. CC1.8030.5193.470.001(0.757-2.849)2.960.003TT vs. TC+CC0.8580.3262.630.012(0.201-1.514)1.590.111To remove studies which may influence the power and stability of the current meta-analysis, sensitivity analysis was applied, but no sensitive case-control studies were found for this SNP among the above five models (such as Figure9(c) regarding T-allele vs. C-allele).


## 3.4. Metaregression

The metaregression analysis indicated that there is a significant relationship for the allele model (T-allele vs. C-allele) with respect to ethnicity, source of control, and genotype methods, with a regression coefficient of 0.001, 0.004, 0.010, and 0.002, respectively. There was no association with publication year, which suggests that the heterogeneity from the rs1801133 polymorphism in EH may be due to the ethnicity, source of control, and genotype method subgroups (Figures10(a)–10(e)).Figure 10
Random-effect metaregression of log odds ratio versus publication year (a), regular ethnicity (b), Chinese/non-Chinese (c), source of control (d), and genotype methods (e), respectively, in EH.
(a)(b)(c)(d)(e)


## 3.5. Gene-Gene Network Diagram and Interactions

Our analysis using the STRING online server indicated thatMTHFR interacts with numerous genes. The ten most significant genes from the network of gene-gene interactions have been listed in Figures 11(a) and 11(b). These ten genes are methionine (MTR), thymidylate synthase (TYMS), C-1-tetrahydrofolate synthase (MTHFD1), serine hydroxymethyltransferase 1 (SHMT1), serine hydroxymethyltransferase 2 (SHMT2), bifunctional methylenetetrahydrofolate dehydrogenase (MTHFD2), probable bifunctional methylenetetrahydrofolate dehydrogenase (MTHFD2L), aminomethyltransferase (AMT), and methionine synthase reductase (MTRR).Figure 11
Human MTHFR interaction network with other genes obtained from STRING server. At least 10 genes have been indicated to correlate with MTHFR gene. (a) The gene-gene interaction. (b) The detail of relative ten core genes.
(a)(b)


## 4. Discussion

While the precise causes of hypertension are still unknown, the risk factors include genetic factors, age, and unhealthy lifestyle practices, with 70-80% of hypertensive cases resulting from unhealthy lifestyle practices. As the risk factors for high blood pressure accumulate, the risk of high blood pressure increases [30, 31].The detection of significant polymorphisms may be a suitable method to predict the risk of hypertension in susceptible individuals. Our current study focused on EH, a common type of hypertension, and included 10,712 patients with EH and 11,916 healthy individuals. In the overall analysis, we observed that individuals carrying TT or a T-allele may have an increased risk of developing EH compared to those with CC or a C-allele (between 37% and 89%). In other words, individuals carrying T-allele or TT genotype may have more possibility to suffer from hypertension in a current or in future time. This can give these people warnings, such as regular changes in blood pressure, changing bad habits, moderating physical exercise, or early medical intervention. Significant heterogeneity was indicated in all of the genetic models. To determine the source of the variation, we analyzed the associations in other subgroups, such as ethnicity, source of control, and genotype method. In parallel, significant relationships were also observed with the ethnicity, source of control, and genotype method subgroups, which provided further support for hypothesis that the rs1801133 polymorphism is a risk factor for EH. In addition, when we used metaregression analysis to evaluate the source of heterogeneity, the aforementioned three subgroups emerged as significant sources of variation. The power of our study was 1, suggesting that our conclusions were accurate.Several meta-analyses regarding the rs1801133 polymorphism and hypertension have been published to date. For example, Wu et al. included 30 case-control studies, and their findings supported a role for the rs1801133 polymorphism in the risk of developing EH [32]. Kosmas et al. identified 23 comparisons relating to hypertensive disease in pregnancy, and they concluded that the T-allele of the rs1801133 polymorphism may increase the risk of hypertension in pregnancy by 1.21-fold [33]. Qian et al. combined 26 published studies of both hypertension and hypertension in pregnancy, suggesting that the rs1801133 polymorphism may be one independent risk factor [34]. Finally, Yang et al. performed a meta-analysis of 27 studies, including 5,418 EH and 4,997 controls, and their findings supported the evidence that carriers of the rs180113 T-allele were susceptible to EH [35]. However, the above studies were subject to some disadvantages. First, in several studies, the HWE values were not consistently above 0.05, which may have increased the heterogeneity and reduced the power of the conclusions. Second, each study omitted other related case-control studies, whereas our current study is a relative comprehensive analysis. Third, some articles did not distinguish the types of hypertension, which may have introduced variability because different kinds of hypertension have different etiologies, pathogenesis, and genetic backgrounds. For these reasons, we focused on one form of hypertension for our analysis. Fourth, our analysis increased genotype subgroup and evaluation of power. Fifth, we analyzed gene-gene interactions between MTHFR and related genes to elucidate potential functional interactions. Despite these disadvantages, it should be noted that the conclusions from our current study are similar to those from previous meta-analyses.A key feature of our study was the identification of gene-gene interactions forMTHFR. The average scores of the ten most significant genes were more than 0.9. The top three genes that were suggested to interact with MTHFR were MTR (0.995), TYMS (0.992), and MTHFD1 (0.989). MTHFR and MTR both participate in homocysteine metabolism, regulating different stages of the process. MTHFR converts 5,10-methylene-THF into 5-methyl-THF, while MTR catalyzes the demethylation of 5-methyl-THF to THF and the remethylation of homocysteine to methionine [36, 37]. The MTR 2756 A/G polymorphism is also associated with the risk of hypertension [38].Some limitations in our meta-analysis should be considered. First, our analyses indicated that the heterogeneity identified in our study emerged from ethnicity, source of control, and genotype methods. Future studies should optimize the design of both retrospective and prospective research projects to overcome this deficiency. Second, EH is a complex disease including genetic and other factors (such as environment, diet, and concomitant disease) [39], and future studies should analyze the gene-gene or gene-environment interactions using larger sample sizes. Third, further meta-analyses should cover all kinds of hypertension, analyzing the associations and determining the genetic background of each type individually. Fourth, the specific mechanism underpinning the impact of the rs1801133 or rs1801131 polymorphism should be explored, for example, using animal models of EH.


## 5. Conclusion

Our meta-analysis provides evidence that theMTHFR 677T null genotype is associated with increased risk of EH, suggesting that further, well-designed large studies are necessary to confirm this relationship. Furthermore, it will be important to focus on the mechanism of action of the rs1801133T-allele to explain the complete chain of evidence for the prevention of EH in the future.



---

*Source: 2144443-2022-02-22.xml*