Skip to main content

SYSTEMATIC REVIEW article

Front. Neurol. , 20 March 2025

Sec. Dementia and Neurodegenerative Diseases

Volume 16 - 2025 | https://doi.org/10.3389/fneur.2025.1556535

This article is part of the Research Topic Cellular Senescence and its Impact on Neurological Disorders View all articles

Association between methylenetetrahydrofolate reductase C677T polymorphism and cerebral small vessel disease: a systematic review and meta-analysis

\r\nHao-tao Zheng,&#x;Hao-tao Zheng1,2Wen-wen Lai&#x;Wen-wen Lai3Jian-jun WangJian-jun Wang2Fan-xin KongFan-xin Kong2Hao-bin CaiHao-bin Cai2Song-jun LinSong-jun Lin2Xu WangXu Wang2Dong-bin CaiDong-bin Cai2Min Pi
Min Pi4*Xiu-de Qin
Xiu-de Qin2*
  • 1Department of Encephalopathy and Phycology, The Fourth Clinical Medical College of Guangzhou University of Chinese Medicine, Shenzhen, China
  • 2Department of Encephalopathy and Phycology, Shenzhen Traditional Chinese Medicine Hospital, Shenzhen, China
  • 3Department of Child Healthcare, Luohu District Maternal and Child Health Hospital, Shenzhen, China
  • 4Department of Preventive Healthcare and Hospital Infection Management, Shenzhen Hospital of Traditional Chinese Medicine, Shenzhen, China

Objective: This systematic review and meta-analysis aimed to evaluate the association between the methylenetetrahydrofolate reductase (5,10-methylenetetrahydrofolate reductase, MTHFR) cytosine (C)677thymine (T) polymorphism and cerebral small vessel disease (CSVD), addressing potential sources of heterogeneity and publication bias.

Methods: An extensive search of databases, including PubMed, the Excerpta Medical Database, and The Cochrane Database of Systematic Reviews, was conducted to identify studies assessing the prevalence of the MTHFR C677T variant associated with CSVD subtypes in humans. Random or fixed effects models were used to accommodate heterogeneity across the study results. Odds ratios (ORs) and weighted mean differences with 95% confidence intervals (CIs) were used for pooled analyses of the relationships between the MTHFR C677T variant associated and CSVD subtypes. Subgroup analyses and assessments of publication bias were performed using Stata software.

Results: Nineteen studies involving 12,441 participants were included. Significant associations were observed across all genetic models: recessive (OR = 1.33; 95%CI = 1.16, 1.52), dominant (OR = 1.25; 95%CI = 1.14, 1.37), allelic (OR = 1.24; 95%CI = 1.14, 1.35), TT vs. CC (OR = 1.42; 95%CI = 1.25, 1.61), and CT vs. CC (OR = 1.20; 95%CI = 1.09, 1.32). Subgroup analyses revealed stronger associations in CSVD-NOS. However, the trim-and-fill method indicated significant publication bias, with adjusted ORs becoming non-significant (recessive model: OR =1.10, 95% CI=0.81, 1.49). Heterogeneity was low to moderate across models (I2 = 14.2–32.4%).

Conclusion: This study highlights the significant association between MTHFR C677T genotyping and CSVD. Early assessment of MTHFR C677T genotyping during the clinical evaluation of elderly patients may improve patient management and reduce the adverse prognostic impact of the CSVD burden. However, further validation of these findings in large-scale, high-quality prospective studies is required.

Systematic review registration: https://www.crd.york.ac.uk/prospero/; identifier: CRD42023339320.

1 Introduction

The incidence of cerebral small vessel disease (CSVD) increases with increasing age, and CSVD is characterized by a range of pathological changes that affect the small arteries, arterioles, venules, and capillaries of the brain. It significantly contributes to age-associated cerebrovascular incidents, including stroke, and is implicated in the decline of vascular-related cognitive function, mood disorders, and mobility challenges (1). Epidemiological data indicate that CSVD occurs in ~50% of all patients with dementia and accounts for approximately 20% of stroke cases worldwide, including a quarter of ischemic stroke cases, approximately one-fifth of patients with stroke have persistent disabilities (2, 3). Traditional risk factors such as smoking, hypertension, diabetes, hyperlipidemia, and obstructive sleep apnea, though extensively investigated, account for only a few CSVD cases. Growing evidence from longitudinal cohort studies highlights the influence of genetic factors on CSVD incidence, suggesting a synergistic interplay between genetic and environmental elements in the etiology of the disease (4, 5). Homocysteine (Hcy), an endogenous amino acid, is associated with endothelial dysfunction and extracellular matrix proliferation, which may lead to vascular damage. Observational research has consistently associated elevated plasma total Hcy (tHcy) levels with an increased CSVD risk (68). Furthermore, studies have reinforced the hypothesis that heightened tHcy levels exacerbate vascular pathology (8). Hcy metabolism, which involves folic acid effects, crucial enzymatic activities, and gene polymorphisms, is a key area of study. In particular, the MTHFR gene variant, characterized by a cytosine (C) to thymine (T) substitution at position 677 (rs1801133), was identified as the most significant genetic factor correlated with elevated tHcy levels.

Recent studies have established a correlation between hyperhomocysteinemia (HHcy) and CSVD (7). However, concerns have been raised regarding the reliability of these findings, which is attributed to confounding factors that influence peripheral Hcy levels. Gene detection methodologies offer a strategy to circumvent these confounding effects. Emerging research has revealed an association between a mutation in the Hcy-related MTHFR C677T gene and CSVD pathogenesis (5, 9).

Emerging research has highlighted the association between the MTHFR C677T gene mutation, related to Hcy, and the pathogenesis of CSVD (6). This correlation has attracted research attention, with studies exploring the relationship between the MTHFR C677T polymorphism and both cerebral infarction and vascular cognitive impairment. These investigations have provided preliminary clues regarding the potential link between the MTHFR C677T genotype and CSVD; however, a comprehensive quantitative assessment is currently lacking.

Furthermore, existing research examining the relationship between the MTHFR C677T genotype and specific CSVD subtypes, as well as imaging burden, is limited, and the conclusions are inconsistent (1013). This inconsistency may be due to differences in sample size, study design, and the rigor of execution, which impact the consistency and reliability of the findings.

Therefore, this meta-analysis aims to clarify the association between the MTHFR C677T gene polymorphism and CSVD by analyzing existing studies. It will assess the impact of the MTHFR C677T gene polymorphism on the risk of CSVD, offering insights for further research into the pathophysiology of CSVD.

2 Methods

This study adhered rigorously to the established protocols delineated in the Meta-analysis of Observational Studies in Epidemiology guidelines (MOOSE) (14). Data extraction was conducted in strict accordance with the 2020 Preferred Reporting Items for Systematic Reviews and Meta-Analyses guidelines (15).

2.1 Data sources

An exhaustive online search was systematically performed to identify suitable studies using key electronic databases, including PubMed (MEDLINE), the Excerpta Medical Database, and The Cochrane Database of Systematic Reviews. The search spanned the period from the inception of each database until December, 2024. Comprehensive descriptions of the search methodologies and strategies applied to each database are provided in the Supplementary Information. We augmented our study collection by manually reviewing the bibliographies of pertinent research articles. In instances of duplicate publications, the study with the largest cohort was selected for inclusion in the analysis.

2.2 Inclusion criteria

The studies were required to meet the criteria as follow: (1) Population (P): Adults were diagnosed with CSVD or its subtypes by neuroimaging (CT or MRI). (2) Intervention/Exposure (I): Genotyping of MTHFR C677T polymorphism with validated methods (such as TaqMan, PCR-RFLP, etc.). (3) Comparison (C): Control groups without CSVD confirmed by neuroimaging (CT or MRI). (4) Outcomes (O): Report the frequencies of cases associated with the MTHFR C677T polymorphism and CSVD. (5) Study Design (S): Observational studies (case-control, cross-sectional, cohort) published in peer-reviewed journals. In cases in which the studies included overlapping participant cohorts, we included only the most extensive dataset available for analysis.

2.3 Exclusion criteria

In the present meta-analysis, exclusion criteria were (1) Repeatedly published data; (2) Literature with incomplete data or lacking target indicators; (3) Case reports, prevalence studies, review articles, letters, conference abstracts, and editorials.

2.4 Data abstraction

Two independent reviewers meticulously screened the articles and evaluated the titles, abstracts, and full texts of each article. They also assessed the eligibility of the studies and systematically extracted relevant information and data. In instances where a consensus was not reached, an independent expert was consulted for resolution. The inclusion criteria for the trials were studies that provided data on specific outcomes, including CSVD subtypes and the prevalence of the MTHFR C677T genotype in both the case and control groups. Three independently investigators extracted data from the selected articles. The extracted data included the publication year, study design, sample size, demographic details, CSVD subtype, and outcome indicators.

2.5 Evaluation of the quality of the studies

The Newcastle-Ottawa Quality Assessment Scale was used to appraise the methodological quality of observational case-control studies. This scale evaluates three primary domains: the selection process of cases and controls, comparability between groups, and accuracy of exposure or outcome ascertainment, with a scoring range of 0–8 points. A higher score indicates superior methodological quality.

The quality of the selected studies was independently assessed by two investigators, and the outcomes of these assessments were meticulously recorded in a dedicated spreadsheet. Discrepancies in ratings were reconciled either through consensus or through the intervention of a third reviewer. An average NOS score of 7.1 indicates the high methodological quality of the included studies. The level of agreement between the raters was remarkably high, as indicated by a Cohen's kappa coefficient of 0.694 (95% confidence interval [CI]: 0.42–0.97).

2.6 Outcomes measures

Studies were considered eligible for inclusion if they presented data addressing specific outcomes, such as CSVD subtypes and the prevalence of the MTHFR C677T genotype in both case and control cohorts. The analysis primarily aimed to discern variations in the aggregated prevalence of individuals exhibiting diverse MTHFR C677T polymorphisms, with a comparison between CSVD and control groups and among different neuroimaging CSVD subtypes.

CSVD encompasses a range of diseases identified through clinical and imaging observations, resulting from pathological processes that affect cerebral arterioles, capillaries, and venules. Neuroimaging features of SVD include recent small subcortical infarcts, lacunes, white matter hyperintensities, perivascular spaces, microbleeds, and brain atrophy (16).

Lacunar infarction (LI): LI referred to as lacunar strokes or lacunar syndromes, represent a morphological entity characterized by small areas of cerebral infarction. These infarctions typically range in size from 0.5 to 15.0 mm in diameter and are predominantly located in deep cerebral regions, including the subcortical white matter, basal ganglia, and pontine base. The central liquefaction in lacunar infarctions results in MRI features that mimicking cerebrospinal fluid, characterized by hypointensity on T1-weighted imaging (T1WI), hyperintensity on T2-weighted imaging (T2WI), and hypointensity on fluid-attenuated inversion recovery (FLAIR) sequences, occasionally surrounded by a hyperintense rim. On diffusion-weighted imaging (DWI), these infarctions exhibit isointense or hypointense signals, and hyperintensity on T2*/susceptibility-weighted imaging (SWI), with diameters typically ranging from 3 to 15 mm (17).

White Matter Hyperintensities (WMH): Vasogenic WMHs exhibit a range of pathological features, including demyelination, gliosis, axonal loss, and oligodendrocyte depletion, with confluent lesions spanning broader regions. The magnetic resonance imaging characteristics of periventricular and deep WMHs encompass thickening of vessel walls, enlargement of the perivascular space (PVS), reduced vascular density, and increased vascular tortuosity. These hyperintensities are observed on T2WI and FLAIR sequences as foci of variable size, with isointense or hypointense signals on T1WI and DWI, and hyperintensity on SWI, which is distinct from the cerebrospinal fluid signal associated with cavitation. They are frequently bilaterally symmetrical and situated adjacent to the lateral ventricles and within the deep white matter (17).

Cerebral Microbleeds (CMBs): Cerebral microbleeds are commonly identified as clusters of iron-laden pigment granules surrounded by macrophages, with associated perilesional tissue alterations such as rarefaction, destruction, and gliosis. These granules may also be present within the PVS or brain tissue without associated tissue damage. On MRI, microbleeds appear as small, round hypointensities on SWI, with isointense signals on T1WI, T2WI, FLAIR, and DWI sequences, typically measuring 2 to 5 mm in diameter, although they can occasionally extend up to 10 mm (17).

2.7 Statistical analysis

The chi-square test was used to determine whether the genotype distribution in the control group adhered to the Hardy–Weinberg equilibrium (HWE). Odds ratios (ORs) with 95% confidence intervals (CIs) were extracted or calculated to evaluate the relationship between the MTHFR C677T polymorphism and CSVD subtypes. Meta-analyses were conducted under five genetic models: recessive model (TT vs. CT + CC), dominant model (TT + CT vs. CC), allelic model (T vs. C), homozygote contrast (TT vs. CC), and heterozygote genotype vs. wild homozygote genotype contrast (CT vs. CC). Log-transformed ORs and CIs were pooled using a random-effects model to account for anticipated heterogeneity across studies. The extent of heterogeneity among study results was quantified using I2 statistics. An I2 value of 0%, < 25%, 25%−50%, and >50% indicated the absence of heterogeneity, low heterogeneity, moderate heterogeneity, and high heterogeneity, respectively. A fixed-effects model was applied only if both I2 ≤ 50% and Cochran's Q-test P >0.10 indicated homogeneity. Subgroup analyses and meta-regression: Predefined subgroups included disease subtype (CSVD-NOS, WMH, LI, CMBs), ethnicity (Asian, European), imaging criteria (MRI, CT/MRI), and HWE compliance (yes/no). Sensitivity analysis was performed to assess the impact of individual studies on the aggregate OR by sequentially excluding one study at a time and recalculating the combined ORs for the remaining studies. Egger's and Begg's tests were conducted to evaluate potential publication bias. These tests are recognized as standard procedures in meta-analyses for examining funnel plot asymmetry. If bias was detected, the trim-and-fill method adjusted for potential missing studies. Two studies (18, 19) were excluded from genotype-specific analyses due to incomplete reporting of CC/CT frequencies. Despite email correspondence with authors, raw data were unavailable. All analyses were performed in Stata 13.0, with statistical significance at two-tailed P < 0.05.

3 Results

3.1 Study selection

The database search yielded 7,557 citations. Following the screening process detailed in Figure 1, 19 full-length articles met inclusion criteria for detailed consideration. Consequently, the meta-analysis included 12,441 participants (3,676 cases and 8,765 controls) across four cerebrovascular phenotypes: 1,143 with CSVD-NOS, 1,125 with WMH, 1,247 with LI, 161 with CMBs. Table 1 presents a detailed overview of studies focusing on the association between the MTHFR C677T genotype and the clinical manifestations of CSVD.

Figure 1
www.frontiersin.org

Figure 1. The prism flowchart indicating the articles' screening and exclusion process up to final inclusion in the qualitative and quantitative analyses.

Table 1
www.frontiersin.org

Table 1. Demographics and clinical features of the case-control studies.

3.2 Study characteristics

The analysis incorporated 19 case-control studies (4, 5, 1833) (Table 1), including 8 in CSVD-NOS (4, 5, 20, 25, 28, 30, 31, 33), 6 in WMH (18, 23, 24, 27, 32, 34), 4 in LI (21, 22, 26, 29), and 1 in CMBs (19). Sample sizes ranged from 90 to 2,345 participants (total N = 12,441), with mean ages of 40–78 years. The geographical distribution of the studies was as follows: 14 studies were from Asia (China: 9; Japan: 3; South Korea: 2), and 5 studies were from Europe (Hungary: 3; Serbia: 1; United Kingdom: 1). The imaging modalities comprised MRI-only protocols (10 studies) (1820, 23, 24, 27, 28, 32, 34), and CT/MRI protocols (nine studies) (4, 5, 21, 22, 25, 26, 2931, 33). Thirteen studies conformed to HWE, while six showed deviations (Table 2). Eighteen studies reported the prevalence of hypertension, while only four studies reported levels of homocysteine (19, 24, 32).

Table 2
www.frontiersin.org

Table 2. Distribution of genotypes in the individual studies.

3.3 Association between the MTHFR C677T gene polymorphism and CSVD

The results demonstrated significant associations between the MTHFR C677T polymorphism and CSVD across all genetic models: recessive model (OR = 1.33; 95%CI = 1.16, 1.52; fixed effects), dominant model (OR = 1.25; 95%CI = 1.14, 1.37; fixed effects), allelic model (OR = 1.24; 95%CI= 1.14, 1.35; random effects), homozygote contrast (TT vs. CC: OR = 1.42; 95%CI = 1.25, 1.61; fixed effects), and heterozygote contrast (CT vs. CC: OR = 1.20; 95%CI = 1.09, 1.32; fixed effects) (Table 3, Figures 26).

Table 3
www.frontiersin.org

Table 3. Main results of the pooled ORs in meta-analysis of the C677T polymorphism.

Figure 2
www.frontiersin.org

Figure 2. Forest plot of MTHFR C677T and CSVD (recessive model: TT vs. CT + CC).

Figure 3
www.frontiersin.org

Figure 3. Forest plot of MTHFR C677T and CSVD (dominant model: TT + CT vs. CC).

Figure 4
www.frontiersin.org

Figure 4. Forest plot of MTHFR C677T and CSVD (allelic model: T vs. C).

Figure 5
www.frontiersin.org

Figure 5. Forest plot of MTHFR C677T and CSVD (homozygote contrast: TT vs. CC).

Figure 6
www.frontiersin.org

Figure 6. Forest plot of MTHFR C677T and CSVD (heterozygote genotype vs. wild homozygote genotype contrast: CT vs. CC).

3.4 Subgroup analysis and meta-regression

Meta-regression identified disease subtype as a significant contributor to heterogeneity in the recessive model (β= −0.20; 95% CI = −0.34, −0.06; P = 0.009), and homozygote contrast model (TT vs. CC: β= −0.20; 95% CI= −0.37, −0.03; P = 0.027). No significant heterogeneity sources were identified in other subgroups (Table 4).

Table 4
www.frontiersin.org

Table 4. Meta-regression of MTHFR C677T polymorphism and CSVD subgroup.

In disease subgroup analyses, notable associations between the MTHFR C677T polymorphism and CSVD-NOS emerged across all genetic models: recessive (OR= 1.73; 95% CI = 1.40, 2.14), dominant (OR = 1.38; 95% CI = 1.15, 1.67), allelic (OR = 1.38; 95% CI= 1.22, 1.57), homozygote contrast (TT vs. CC: OR = 2.04; 95% CI = 1.60, 2.60), and heterozygote contrast (CT vs. CC: OR = 1.23; 95% CI = 1.02, 1.47). For lacunar infarction, no significant associations were observed in most genetic models. White matter hyperintensity exhibited significant associations in the dominant (OR = 1.21; 95% CI = 1.02, 1.44), allelic (OR = 1.17; 95% CI = 1.04, 1.32), and homozygote contrast (TT vs. CC: OR = 1.37; 95% CI = 1.07, 1.74). Cerebral microbleed analyses were underpowered due to limited data, though a non-significant trend emerged in the recessive model (OR = 1.52; 95% CI = 0.99, 2.33) (Table 5, Supplementary Figures 15).

Table 5
www.frontiersin.org

Table 5. Pooled ORs of C677T polymorphism meta-analysis: subgroup-stratified main results.

Ethnic subgroup analyses revealed significant associations in European populations for recessive (OR = 1.41; 95% CI = 1.03, 1.91), dominant (OR = 1.28; 95% CI = 1.02, 1.62), allelic (OR = 1.30; 95% CI = 1.10, 1.54), homozygote contrast (TT vs. CC: OR = 1.78; 95% CI = 1.25, 2.53) models. Asian populations revealed stronger effect sizes across all models: recessive (OR = 1.32; 95% CI = 1.12, 1.54), dominant (OR = 1.31; 95% CI = 1.13, 1.51), allelic (OR = 1.24; 95% CI = 1.12, 1.36), homozygote contrast (TT vs. CC: OR = 1.49; 95% CI = 1.22–1.83), heterozygote contrast (CT vs. CC: OR = 1.24; 95% CI = 1.08–1.44) (Table 5, Supplementary Figures 15).

MRI-based studies showed significant associations in recessive (OR = 1.41, 95% CI = 1.17, 1.71), dominant (OR = 1.24; 95% CI = 1.07, 1.45), allelic (OR = 1.25; 95% CI = 1.12, 1.40), and homozygote contrast (TT vs. CC: OR = 1.57; 95% CI = 1.24, 1.97) models. CT/MRI studies maintained consistent significance across all models (OR range = 1.23, 1.49) (Table 5, Supplementary Figures 15).

HWE-compliant studies demonstrated robust associations: recessive (OR = 1.44; 95% CI = 1.22, 1.70), dominant (OR = 1.37; 95% CI = 1.16, 1.63), allelic (OR = 1.30; 95% CI = 1.16, 1.45), homozygote contrast (OR = 1.66; 95% CI = 1.34, 2.05), and heterozygote contrast (CT vs. CC: OR = 1.28; 95% CI = 1.07, 1.52). HWE-deviant studies exhibited attenuated effects: dominant (OR = 1.17; 95% CI = 1.02, 1.35), and allelic (OR = 1.12; 95% CI= 1.02, 1.23) (Table 5, Supplementary Figures 15).

3.5 Heterogeneity evaluation

The between study heterogeneity across genetic models was systematically evaluated through I2 statistics and Q-test (Table 3). The recessive model showed low-to-moderate heterogeneity (I2 = 29.3%, P = 0.113), while the dominant model exhibited minimal heterogeneity (I2 = 24.1%, P = 0.176). Comparable heterogeneity levels were observed in the allelic (I2 = 32.4%, P = 0.097) and homozygote contrast (TT vs. CC: I2 = 31.8%, P = 0.102) models, with the heterozygote model demonstrating the lowest heterogeneity (I2 = 14.2%, P = 0.287).

Disease subgroup analyses revealed minimal heterogeneity for CSVD-NOS (I2 ≤ 22.5%, P ≥ 0.250), significant heterogeneity in LI (dominant model: I2 = 68.0%, P = 0.025; heterozygote model: I2 = 70.5%, P = 0.017), and negligible heterogeneity in WMH (I2 = 0.0%, P ≥ 0.685).

Ethnic subgroups displayed differential heterogeneity patterns: European populations showed moderate heterogeneity in the recessive model (I2 = 27.1%, P = 0.241), other genetic models exhibited low heterogeneity. While Asian populations exhibited systematically higher heterogeneity across models (I2 = 33.4–41.2%, P = 0.060–0.115).

Imaging subgroups demonstrated low heterogeneity in MRI-based studies (I2 ≤ 23.4%, P ≥ 0.228) vs. moderate-to-high heterogeneity in CT/MRI studies (I2 = 33.2–56.3%, P = 0.019–0.152).

HWE stratification analyses showed consistently lower heterogeneity in HWE-compliant studies (I2 = 15.7–39.8%, P ≥ 0.075) compared to HWE-deviant studies (I2 ≤ 26.0%, P ≥ 0.248).

3.6 Publication bias and sensitivity analysis

To assess the potential effects of small studies, both Begg's rank correlation test and Egger's linear regression analysis were used, supplemented by a visual examination of comparison-adjusted funnel plots for each outcome (Supplementary Figures 1115). The heterozygote model (CT vs. CC) showed no statistically significant publication bias through either method (Begg's z = 1.94, P = 0.053; Egger's t = 1.64, P = 0.123). Conversely, substantial evidence of publication bias across the other four genetic models (Table 4).

The trim-and-fill method revealed significant publication bias in the meta-analysis of the MTHFR C677T polymorphism and CSVD. After adjustment, ORs decreased to non-significant levels (recessive model: OR = 1.10, 95%CI = 0.81–1.49; dominant model: OR = 1.12, 95%CI = 0.85–1.49; allelic model: OR = 1.10, 95%CI = 0.90–1.33; TT vs. CC: OR = 1.18, 95%CI = 0.79–1.76; CT vs. CC: OR = 1.09, 95%CI = 0.82–1.45), indicating publication bias likely inflated initial estimates.

For sensitivity analysis, we employed a sequential exclusion method for each study. Upon excluding any single study, the aggregated results of the remaining studies demonstrated no significant deviation from the original combined results indicating robust stability in our findings (Supplementary Figures 610).

4 Discussion

Earlier research endeavors have begun to delineate a potential link between the MTHFR C677T polymorphism and various CSVD phenotypes, including white matter hyperintensities, lacunar infarcts, and cerebral microbleeds, conclusive evidence has been elusive. This systematic review and meta-analysis of 20 case-control studies investigates the association between MTHFR C677T polymorphism and CSVD subtypes in humans. The design of case-control studies provided our analysis with a robust capacity for causal inference. An average bias risk score of 7.1 indicates the high methodological quality of the included studies, enhancing the credibility of our findings. Key findings demonstrate significant associations across genetic models, though publication bias attenuated effects after trim-and-fill adjustment. Sensitivity analysis using sequential exclusion confirmed the stability of pooled estimates, with no single study disproportionately influencing the results.

The etiology of CSVD encompasses a multifaceted array of mechanisms, such as endothelial dysfunction, compromised blood-brain barrier integrity, oxidative stress, mitochondrial impairment, and inflammatory processes. These factors contribute to the development of arteriolar sclerosis, the tortuous course of deep white matter microvessels, and a reduction in capillary density, culminating in a cascade of pathological brain tissue alterations (35, 36). These alterations manifest as white matter injury, lacunar infarction, cerebral microbleeds, and expanded perivascular spaces. Previous studies have demonstrated that small cerebral arteries exhibit heightened sensitivity to HHcy, suggesting that the impact of HHcy on CSVD may be particularly pronounced. A significant correlation has been established between HHcy and CSVD, with HHcy identified as an independent risk factor for CSVD (37). The MTHFR C677T polymorphism has been shown to diminish enzyme activity and thermal stability, with the enzyme activity in individuals with the CT and TT genotypes being 65% and 30% of that observed in the CC genotype, respectively. Furthermore, individuals harboring the TT genotype exhibit serum homocysteine levels that are 20% to 70% higher than those of the CC genotype (38). We propose that the MTHFR C677T polymorphism diminishes the enzymatic activity and thermal stability of the enzyme, leading to dysregulated homocysteine metabolism, potentially precipitating CSVD through multiple pathways, including endothelial dysfunction, oxidation of low-density lipoproteins, and heightened oxidative stress.

In our subgroup analysis, the MTHFR C677T genotype correlated with LI, WMH and CMBs, displaying a trend, yet failing to attain statistical significance. WMH, a prevalent subtype of CSVD previously considered a benign aspect of aging, have been reevaluated in the studies as a marker indicative of compromised cerebral and cardiovascular health. Studies have consistently demonstrated that WMH is associated with an elevated risk of stroke, cognitive impairment, depressive disorders, and mortality (39). The relationship between the MTHFR C677T polymorphism and WMH is controversial, with mixed findings from previous studies. Rutten-Jacobs' study suggests a correlation between the MTHFR C677T polymorphism and WMH volume. Another study, using Mendelian randomization, suggested a correlation between MTHFR C677T polymorphism and WMH (4042). However, several other studies have failed to replicate these findings (24, 32, 34, 43). Our analysis suggests a positive trend but without statistical significance. The possible reason may be that WMH is influenced by various factors such as age, lipids, fasting glucose, smoking, inflammatory response, vitamin B12 levels, etc., and most of the current studies did not incorporate the additional variables, which interfered with the effect of MTHFR gene polymorphisms on WMH.

Lacunar infarction, the predominant subtype of CSVD, implicates plasma homocysteine in its pathogenesis. The potential mechanisms underlying homocysteine-induced atherosclerotic thrombosis encompass endothelial dysfunction, proliferation of vascular smooth muscle cells, and oxidation of low-density lipoproteins (44, 45). Previous researches have proposed a possible association between the MTHFR C677T polymorphism and lacunar infarction, potentially through increased susceptibility to, or interaction with, high blood pressure (26, 42). However, a subset of studies has failed to establish significant correlations (21, 29). Our meta-analysis did not demonstrate a significant correlation between lacunar infarction and MTHFR genotype, which may be attributed to the omission of homocysteine and hypertension as covariates in the original studies.

CMBs research is in its initial stages. Elevated total homocysteine levels have been observed in patients exhibiting a high burden of CMBs, implying that serum tHcy levels may exert an independent influence on the occurrence of CMBs (46). The correlation between the MTHFR C677T TT genotype and CMBs may be attributed to homocysteine's role in coagulation, leading to arteriolar and venous infarctions that result in microbleeds and the initiation of microhemorrhages (47).

Although our meta-analysis identified a nominally significant association between the MTHFR C677T polymorphism and CSVD, we recognize the need to interpret these findings within clinical practice. The association of the T allele reduced MTHFR enzyme activity, and subsequent hyperhomocysteinemia provides a plausible biological basis, particularly in folate-deficient populations. This supports the potential application of MTHFR C677T genotyping in CSVD risk stratification, especially in areas with low dietary folate intake or in populations prone to micronutrient deficiencies, such as the elderly or those with malabsorption syndromes. Clinically, identification of T allele carriers may prompt targeted interventions such as homocysteine-lowering therapy and early neuroimaging surveillance. However, current guidelines do not universally recommend MTHFR genotyping for CSVD due to inconsistent evidence. Our findings suggest that such testing may be of value in specific subgroups, such as patients with early-onset CSVD who lack traditional risk factors.

This study has some limitations. First, it suffers from the potential for publication bias noted in the meta-analysis, which was the main concern. Despite our efforts to comprehensively identify and include all studies that met the inclusion criteria, the use of funnel plots and Egger's and Begg's regression analyses indicated a possible risk of publication bias. The trim-and-fill method revealed that publication bias likely influenced the initial results, and the corrected Ors became no longer significant in all genetic models. Second, the methods used to measure white matter hyperintensities varied significantly across studies, contributing to methodological heterogeneity. Notably, many studies reported high inter-rater reliability; however, the subjective nature of assessing white matter hyperintensities using visual rating scales cannot be overlooked. The inclusion of CT and MRI studies introduced further heterogeneity. The lower sensitivity of CT in detecting lacunar infarctions and white matter hyperintensities might have led to underrepresentation in such studies. However, sensitivity analyses revealed that the exclusion of these studies did not significantly alter the pooled log OR values. Consequently, we deemed their inclusion to be justified and appropriate.

5 Conclusion

This study highlights the significant association between MTHFR C677T genotyping and CSVD. Early assessment of MTHFR C677T genotyping during the clinical evaluation of elderly patients may improve patient management and reduce the adverse prognostic impact of the CSVD burden. However, further validation of these findings in large-scale, high-quality prospective studies is required.

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary material.

Author contributions

H-tZ: Data curation, Methodology, Software, Writing – original draft. W-wL: Data curation, Methodology, Software, Writing – original draft. J-jW: Formal analysis, Methodology, Supervision, Writing – original draft. F-xK: Supervision, Visualization, Writing – original draft. H-bC: Methodology, Supervision, Writing – review & editing. S-jL: Data curation, Visualization, Writing – original draft. XW: Data curation, Validation, Writing – original draft. D-bC: Formal analysis, Methodology, Supervision, Writing – original draft. MP: Funding acquisition, Supervision, Writing – review & editing. X-dQ: Supervision, Writing – review & editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This work was supported by Shenzhen's Sanming Project (No. SZZYSM 20211101), Shenzhen Science and Technology Program (No. JCYJ20220531092205012), and Provincial Natural Science Foundation of Guangdong (No. 2023A1515030265).

Acknowledgments

The authors would like to thank all study participants who were enrolled in this study.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Gen AI was used in the creation of this manuscript.

Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fneur.2025.1556535/full#supplementary-material

References

1. Duering M, Biessels GJ, Brodtmann A, Chen C, Cordonnier C, de Leeuw FE, et al. Neuroimaging standards for research into small vessel disease-advances since 2013. Lancet Neurol. (2023) 22:602–18. doi: 10.1016/S1474-4422(23)00131-X

PubMed Abstract | Crossref Full Text | Google Scholar

2. Huo Y, Li J, Qin X, Huang Y, Wang X, Gottesman RF, et al. Efficacy of folic acid therapy in primary prevention of stroke among adults with hypertension in China: the CSPPT randomized clinical trial. JAMA. (2015) 313:1325–35. doi: 10.1001/jama.2015.2274

PubMed Abstract | Crossref Full Text | Google Scholar

3. Valdes HMC, Morris Z, Dickie DA, Royle NA, Muñoz Maniega S, Aribisala BS, et al. Close correlation between quantitative and qualitative assessments of white matter lesions. Neuroepidemiology. (2013) 40:13–22. doi: 10.1159/000341859

PubMed Abstract | Crossref Full Text | Google Scholar

4. Hassan A. Homocysteine is a risk factor for cerebral small vessel disease, acting via endothelial dysfunction. Brain. (2004) 127:212–9. doi: 10.1093/brain/awh023

PubMed Abstract | Crossref Full Text | Google Scholar

5. Yuan H, Fu M, Yang X, Huang K, Ren X. Single nucleotide polymorphism ofMTHFR rs1801133 associated with elevated Hcy levels affects susceptibility to cerebral small vessel disease. PeerJ. (2020) 8:e8627. doi: 10.7717/peerj.8627

PubMed Abstract | Crossref Full Text | Google Scholar

6. Cao Y, Su N, Zhang D, Zhou L, Yao M, Zhang S, et al. Correlation between total homocysteine and cerebral small vessel disease: a Mendelian randomization study. Eur J Neurol. (2021) 28:1931–8. doi: 10.1111/ene.14708

PubMed Abstract | Crossref Full Text | Google Scholar

7. Feng C, Bai X, Xu Y, Hua T, Huang J, Liu X-Y. Hyperhomocysteinemia associates with small vessel disease more closely than large vessel disease. Int J Med Sci. (2013) 10:408–12. doi: 10.7150/ijms.5272

PubMed Abstract | Crossref Full Text | Google Scholar

8. Nam K-W, Kwon H-M, Jeong H-Y, Park J-H, Kwon H, Jeong S-M. Serum homocysteine level is related to cerebral small vessel disease in a healthy population. Neurology. (2019) 92:e317–25. doi: 10.1212/WNL.0000000000006816

PubMed Abstract | Crossref Full Text | Google Scholar

9. Annus Á, Juhász LÁ, Szabó E, Rárosi F, Szpisjak L, Vécsei L, et al. Connection between small vessel disease related stroke and the MTHFR C677T polymorphism in a Hungarian population. Heliyon. (2020) 6:e05305. doi: 10.1016/j.heliyon.2020.e05305

PubMed Abstract | Crossref Full Text | Google Scholar

10. Jacob MA, Cai M, van de Donk V, Bergkamp M, Marques J, Norris DG, et al. Cerebral small vessel disease progression and the risk of dementia: a 14-year follow-up study. Am J Psychiatry. (2023) 180:508–18. doi: 10.1176/appi.ajp.20220380

PubMed Abstract | Crossref Full Text | Google Scholar

11. Mishra A, Duplaa C, Vojinovic D, Suzuki H, Sargurupremraj M, Zilhão NR, et al. Gene-mapping study of extremes of cerebral small vessel disease reveals TRIM47 as a strong candidate. Brain. (2022) 145:1992–2007. doi: 10.1093/brain/awab432

PubMed Abstract | Crossref Full Text | Google Scholar

12. Rudilosso S, Mena L, Esteller D, Olivera M, Mengual JJ, Montull C, et al. Higher Cerebral Small Vessel Disease Burden in Patients with White Matter Recent Small Subcortical Infarcts. J Stroke Cerebrovasc Dis. (2021) 30:105824. doi: 10.1016/j.jstrokecerebrovasdis.2021.105824

PubMed Abstract | Crossref Full Text | Google Scholar

13. Sargurupremraj M, Suzuki H, Jian X, Sarnowski C, Evans TE, Bis JC, et al. Cerebral small vessel disease genomics and its implications across the lifespan. Nat Commun. (2020) 11:6285.

PubMed Abstract | Google Scholar

14. Stroup DF, Berlin JA, Morton SC, et al. Meta-analysis of observational studies in epidemiology: a proposal for reporting. Meta-analysis of observational studies in epidemiology (MOOSE) group. JAMA. (2000) 283:2008–12. doi: 10.1001/jama.283.15.2008

PubMed Abstract | Crossref Full Text | Google Scholar

15. Moher D, Liberati A, Tetzlaff J, Altman DG, PRISMA Group. Preferred reporting items for systematic reviews and meta-analyses: the PRISMA statement. Int J Surg. (2010) 8:336–41. doi: 10.1016/j.ijsu.2010.02.007

PubMed Abstract | Crossref Full Text | Google Scholar

16. Wardlaw JM, Smith EE, Biessels GJ, Cordonnier C, Fazekas F, Frayne R, et al. Neuroimaging standards for research into small vessel disease and its contribution to ageing and neurodegeneration. Lancet Neurol. (2013) 12:822–38. doi: 10.1016/S1474-4422(13)70124-8

PubMed Abstract | Crossref Full Text | Google Scholar

17. Shi Y, Wardlaw JM. Update on cerebral small vessel disease: a dynamic whole-brain disease. BMJ. (2016) 1:83–92. doi: 10.1136/svn-2016-000035

PubMed Abstract | Crossref Full Text | Google Scholar

18. Szolnoki Z, Somogyvári F, Kondacs A, Szabó M, Fodor L, Bene J, et al. Specific APO E genotypes in combination with the ACE D/D or MTHFR 677TT mutation yield an independent genetic risk of leukoaraiosis. Acta Neurol Scand. (2004) 109:222–7. doi: 10.1046/j.1600-0404.2003.00218.x

PubMed Abstract | Crossref Full Text | Google Scholar

19. Yoo JS, Ryu C-H, Kim YS, Kim H-J, Bushnell CD, Kim HY. Homocysteinemia is associated with the presence of microbleeds in cognitively impaired patients. J Stroke Cerebrovas Dis. (2020) 29:105302. doi: 10.1016/j.jstrokecerebrovasdis.2020.105302

PubMed Abstract | Crossref Full Text | Google Scholar

20. Choi BO, Kim NK, Kim SH, Kang MS, Lee S, Ahn JY, et al. Homozygous C677T mutation in the MTHFR gene as an independent risk factor for multiple small-artery occlusions. Thromb Res. (2003) 111:39–44. doi: 10.1016/j.thromres.2003.08.022

PubMed Abstract | Crossref Full Text | Google Scholar

21. Feng C, Yang Y, Yang S, Tu X, Wang Y, Song Y, et al. Effect of gene–gene and gene–environment interaction on the risk of first-ever stroke and poststroke death. Mol Genet Genomic Med. (2019) 7:846. doi: 10.1002/mgg3.846

PubMed Abstract | Crossref Full Text | Google Scholar

22. Kawamoto R, Kohara K, Oka Y, et al. An association of 5,10-Methylenetetrahydrofolate Reductase (MTHFR) gene polymorphism and ischemic stroke. J Stroke Cerebrovas Dis. (2005) 14:67–74. doi: 10.1016/j.jstrokecerebrovasdis.2004.12.003

PubMed Abstract | Crossref Full Text | Google Scholar

23. Kohara K, Fujisawa M, Ando F, Tabara Y, Niino N, Miki T, et al. MTHFR gene polymorphism as a risk factor for silent brain infarcts and white matter lesions in the Japanese general population. Stroke. (2003) 34:1130–5. doi: 10.1161/01.STR.0000069163.02611.B0

PubMed Abstract | Crossref Full Text | Google Scholar

24. Li M, Fu B, Dong W. Correlations between plasma homocysteine and MTHFR gene polymorphism and white matter lesions. Folia Neuropathol. (2018) 56:301–7. doi: 10.5114/fn.2018.80863

PubMed Abstract | Crossref Full Text | Google Scholar

25. Luo M, Li J, Sun X, Lai R, Wang Y, Xu X, et al. Interactions among candidate genes selected by meta-analyses resulting in higher risk of ischemic stroke in a Chinese population. PLoS ONE. (2015) 10:e0145399. doi: 10.1371/journal.pone.0145399

PubMed Abstract | Crossref Full Text | Google Scholar

26. Mao X, Han L. The relationship of methylenetetrahydrofolate reductase gene C677T polymorphism and ischemic stroke in Chinese Han population. Ann Clin Lab Sci. (2018) 48:242–7.

Google Scholar

27. Notsu Y, Nabika T, Park HY, Masuda J, Kobayashi S. Evaluation of genetic risk factors for silent brain infarction. Stroke. (1970) 30:1881–6. doi: 10.1161/01.STR.30.9.1881

PubMed Abstract | Crossref Full Text | Google Scholar

28. Pavlovic AM, Pekmezovic T, Obrenovic R, Novakovic I, Tomic G, Mijajlovic M, et al. Increased total homocysteine level is associated with clinical status and severity of white matter changes in symptomatic patients with subcortical small vessel disease. Clin Neurol Neurosurg. (2011) 113:711–5. doi: 10.1016/j.clineuro.2011.07.004

PubMed Abstract | Crossref Full Text | Google Scholar

29. Shen C, Zhang W, Sun K, Ling MX, Shen CC, Liu WG, et al. Interaction of genetic risk factors confers higher risk for thrombotic stroke in male Chinese: a multicenter case-control study. Ann Hum Genet. (2007) 71:620–9. doi: 10.1111/j.1469-1809.2007.00364.x

PubMed Abstract | Crossref Full Text | Google Scholar

30. Szolnoki Z, Maasz A, Magyari L, Horvatovich K, Farago B, Somogyvari F, et al. The combination of homozygous MTHFR 677T and angiotensin II type-1 receptor 1166C variants confers the risk of small-vessel-associated ischemic stroke. J Mol Neurosci. (2007) 31:201–7. doi: 10.1385/JMN:31:03:201

PubMed Abstract | Crossref Full Text | Google Scholar

31. Yeh P-S, Lin H-J, Li Y-H, Lin K-C, Cheng T-J, Chang C-Y, et al. Prognosis of young ischemic stroke in Taiwan: impact of prothrombotic genetic polymorphisms. Thromb Haemost. (2004) 92:583–9. doi: 10.1160/TH04-03-0099

PubMed Abstract | Crossref Full Text | Google Scholar

32. You M, Zhou X, Yin W, Wan K, Zhang W, Li C, et al. The influence of MTHFR polymorphism on gray matter volume in patients with amnestic mild cognitive impairment. Front Neurosci. (2021) 15:778123. doi: 10.3389/fnins.2021.778123

PubMed Abstract | Crossref Full Text | Google Scholar

33. Chen L, Wu C, Dong Z, Cao S, Ren N, Yan X. Methylenetetrahydrofolate reductase polymorphisms and elevated plasma homocysteine levels in small vessel disease. Brain Behav. (2023) 13:2960. doi: 10.1002/brb3.2960

PubMed Abstract | Crossref Full Text | Google Scholar

34. Szolnoki Z, Szaniszlo I, Szekeres M, Hitri K, Kondacs A, Mandi Y, et al. Evaluation of the MTHFR A1298C variant in leukoaraiosis. J Mol Neurosci. (2012) 46:492–6. doi: 10.1007/s12031-011-9621-4

PubMed Abstract | Crossref Full Text | Google Scholar

35. Li T, Huang Y, Cai W, Chen X, Men X, Lu T, et al. Age-related cerebral small vessel disease and inflammaging. Cell Death Dis. (2020) 11:932. doi: 10.1038/s41419-020-03137-x

PubMed Abstract | Crossref Full Text | Google Scholar

36. Martinez SS, Smith KJ. Understanding a role for hypoxia in lesion formation and location in the deep and periventricular white matter in small vessel disease and multiple sclerosis. Clin Sci (Lond). (2017) 131:2503–24. doi: 10.1042/CS20170981

PubMed Abstract | Crossref Full Text | Google Scholar

37. Li S, Li G, Luo X, Huang Y, Wen L, Li J. Endothelial dysfunction and hyperhomocysteinemia-linked cerebral small vessel disease: underlying mechanisms and treatment timing. Front Neurol. (2021) 12:736309. doi: 10.3389/fneur.2021.736309

PubMed Abstract | Crossref Full Text | Google Scholar

38. Rozen R. Genetic predisposition to hyperhomocysteinemia: deficiency of methylenetetrahydrofolate reductase (MTHFR). Thromb Haemost. (1997) 78:523–6. doi: 10.1055/s-0038-1657581

Crossref Full Text | Google Scholar

39. Ottavi TP, Pepper E, Bateman G, Fiorentino M, Brodtmann A. Consensus statement for the management of incidentally found brain white matter hyperintensities in general medical practice. Med J Aust. (2023) 219:278–84. doi: 10.5694/mja2.52079

PubMed Abstract | Crossref Full Text | Google Scholar

40. Hong ED, Taylor WD, McQuoid DR, Potter GG, Payne ME, Ashley-Koch A, et al. Influence of the MTHFR C677T polymorphism on magnetic resonance imaging hyperintensity volume and cognition in geriatric depression. Am J Geriatr Psychiatry. (2009) 17:847–55. doi: 10.1097/JGP.0b013e3181aad5b2

PubMed Abstract | Crossref Full Text | Google Scholar

41. Li Z, Wu X, Huang H, Xu F, Liang G, Lin C, et al. MTHFR C677T polymorphism and cerebrovascular lesions in elderly patients with CSVD: a correlation analysis. Front Genet. (2022) 13:987519. doi: 10.3389/fgene.2022.987519

PubMed Abstract | Crossref Full Text | Google Scholar

42. Rutten-Jacobs LCA, Traylor M, Adib-Samii P, Thijs V, Sudlow C, Rothwell PM, et al. Association of MTHFR C677T genotype with ischemic stroke is confined to cerebral small vessel disease subtype. Stroke. (2016) 47:3. doi: 10.1161/STROKEAHA.115.011545

PubMed Abstract | Crossref Full Text | Google Scholar

43. Rajagopalan P, Jahanshad N, Stein JL, Hua X, Madsen SK, Kohannim O, et al. Common folate gene variant, MTHFR C677T, is associated with brain structure in two independent cohorts of people with mild cognitive impairment. Neuroimage Clin. (2012) 1:179–87. doi: 10.1016/j.nicl.2012.09.012

PubMed Abstract | Crossref Full Text | Google Scholar

44. Fan H, Yang S, Li Y, Yin J, Qin W, Yang L, et al. Assessment of homocysteine as a diagnostic and early prognostic biomarker for patients with acute lacunar infarction. Eur Neurol. (2018) 79:54–62. doi: 10.1159/000484893

PubMed Abstract | Crossref Full Text | Google Scholar

45. Spence JD. Homocysteine lowering for stroke prevention: unravelling the complexity of the evidence. Int J Stroke. (2016) 11:744–7. doi: 10.1177/1747493016662038

PubMed Abstract | Crossref Full Text | Google Scholar

46. Wang X-B, Dong H, Qiu Y-G, Lou C-C, Huang D-Y, Zhang J, et al. Nomogram based on clinical and brain computed tomography characteristics for predicting more than 5 cerebral microbleeds in the hypertensive population. Front Neurol. (2022) 13:955378. doi: 10.3389/fneur.2022.955378

PubMed Abstract | Crossref Full Text | Google Scholar

47. Zhang DD, Cao Y, Mu JY, Liu YM, Gao F, Han F, et al. Inflammatory biomarkers and cerebral small vessel disease: a community-based cohort study. Stroke Vasc Neurol. (2022) 7:302–9. doi: 10.1136/svn-2021-001102

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: MTHFR C677T polymorphism, cerebral small vessel disease, meta-analysis, white matter hyperintensities, lacunar infarction, microbleeds

Citation: Zheng H-t, Lai W-w, Wang J-j, Kong F-x, Cai H-b, Lin S-j, Wang X, Cai D-b, Pi M and Qin X-d (2025) Association between methylenetetrahydrofolate reductase C677T polymorphism and cerebral small vessel disease: a systematic review and meta-analysis. Front. Neurol. 16:1556535. doi: 10.3389/fneur.2025.1556535

Received: 07 January 2025; Accepted: 28 February 2025;
Published: 20 March 2025.

Edited by:

José Alfonso Cruz Ramos, University of Guadalajara, Mexico

Reviewed by:

Jiawei Wang, Capital Medical University, China
Manabesh Nath, All India Institute of Medical Sciences, India
Osama Youssef, The University of Utah, United States

Copyright © 2025 Zheng, Lai, Wang, Kong, Cai, Lin, Wang, Cai, Pi and Qin. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Min Pi, cG0wMzA1QGd6dWNtLmVkdS5jbg==; Xiu-de Qin, cWlueGl1ZGVAZm94bWFpbC5jb20=

These authors have contributed equally to this work and share first authorship

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Research integrity at Frontiers

Man ultramarathon runner in the mountains he trains at sunset

95% of researchers rate our articles as excellent or good

Learn more about the work of our research integrity team to safeguard the quality of each article we publish.


Find out more