- 1Emergency Department, China-Japan Union Hospital of Jilin University, Changchun, China
- 2Pathology Department, Jilin Cancer Hospital, Changchun, China
- 3Hematology Department, The First Hospital of Jilin University, Changchun, China
- 4Stroke Center, Department of Neurology, The Affiliated Hospital of Beihua University, Changchun, China
Aims: The metabolic syndrome (MetS), a collection of conditions that heighten the risk of disease development and impose economic burdens on patients. However, the causal relationship between education and MetS was uncertain. In this study, the Mendelian randomization (MR) method was employed to elucidate the potential causal link between education and the MetS and its components.
Method: Single nucleotide polymorphisms (SNPs) associated with education, MetS, and its components were sourced from a public database, with the inverse variance-weighted (IVW) method utilized for analysis.
Results: Education demonstrated a significant negative correlation with the risk of MetS (OR = 0.55, 95% CI = 0.48–0.63, p = 2.18E−51), waist circumference(OR = 0.80, 95% CI = 0.76–0.83, p = 4.98E-33), hypertension (OR = 0.96, 95% CI = 0.95–0.97; p = 4.54E-10), Fasting blood glucose (OR = 0.94, 95% CI = 0.91–0.97, p = 7.58E-6) and triglycerides (OR = 0.83, 95% CI = 0.79–0.87, p = 7.87E-18) while showing a positive association with high-density lipoprotein cholesterol (OR = 1.22, 95% CI = 1.18–1.25, p = 1.45E-31).
Conclusion: The findings of this study suggest that education can decrease the incidence of MetS.
1 Introduction
In the global context, the prevalence of metabolic syndrome (MetS) had ranged from 20 to 25% and exhibited gradually increased over the past several years (1, 2). MetS had been associated with an elevated risk of developing various illnesses, including a 29% increase in the risk of major coronary events, a 30% increase in the risk of recurrent myocardial infarction, and a 12% increase in the risk of dementia (3, 4). Meanwhile, MetS had been linked to a rising mortality rate, with a 24% increase in all-cause mortality and a 44% increase in mortality due to heart disease (5).
The level of education, as an indicator of socioeconomic status, has the potential to impact both the prevalence of cardiovascular disease and the prevalence of MetS (6, 7). Individuals with a college education demonstrated a 57% lower risk of developing MetS than those with only a primary school education (8). Recent research findings had indicated that while several studies had explored the association between education and MetS, a definitive causal relationship between the two had not been found.
Mendelian randomization (MR) was initially introduced by Katan in 1986, utilizing genome-wide association study(GWAS) data to employ single nucleotide polymorphisms(SNPs) as instrumental variables(IVs) for investigating causal relationships between exposure and outcome (9). This method distinguishes itself from others by its independence from confounding variables and ability to avoid reverse causation (10).
This study utilized a 2-sample MR method to estimate the association between education and MetS.
2 Method
2.1 Overview of the study design
As the assessment criteria outlined by the International Diabetes Federation (IDF) and the American Heart Association/National Heart, Lung, and Blood Institute (AHA/NHLBI) (11), five components (waist circumference (WC), hypertension, fasting blood glucose(FBG), triglycerides(TG), and high-density lipoprotein cholesterol(HDL-C)) were identified as elements of MetS.
In this research, a 2-sample MR analysis was conducted. The initial phase involved investigating the impact of education on MetS and WC, hypertension, FBG, TG, and HDL-C. The selection of IVs shown in Figure 1 was guided by three key assumptions: (1) the genetic instruments are strongly correlated with the exposure variable; (2) the genetic instruments are not affected by any confounders; and (3) the genetic instruments solely influence the outcome through the exposure variable. The study design is depicted in Figure 1.
Figure 1. Overview of the study design in this MR study. MR, Mendelian randomization; SNP, single-nucleotide polymorphism; MetS, metabolic syndrome, FBG, fasting blood glucose, TG, triglycerides, WC, waist circumference, HDL-C, high-density lipoprotein cholesterol.
2.2 Selection of genetic instruments for MR analyses
The criteria for selecting IVs in this study included the following: (1): only SNPs with p values less than 5 × 10–8 were included in the study; (2); SNPs with linkage disequilibrium (LD) based on a specified R2 threshold of 0.001 and a clumping distance of 10,000 kb were removed; and (3) the F statistic was calculated to assess the predictive potential of the IVs, with a threshold of F statistic >10 indicating strong predictive ability (Figure 2).
2.3 Data sources for education
The largest published GWAS for education was obtained from the Social Science Genetic Association Consortium (SSGAC), encompassing a total of 765,283 individuals of European ancestry (12). A total of 418 SNPs associated with WC with p values less than 5 × 10–8 were included in the study.
2.4 Data sources for MetS and its components
The genetic variation data for MetS were acquired from the GWAS, which includes a total of 291,107 individuals of European descent (comprising 59,677 cases and 231,430 controls) (13).
A summary of GWASs for WC was obtained from the Medical Research Council (MRC) Integrative Epidemiology Unit (IEU) GWAS pipeline, which uses a sample size of 407,661 individuals (14).
For hypertension, genetic data were extracted from the MRC-IEU UK Biobank GWAS pipeline, and a total of 484,598 subjects of European ancestry (129,909 cases and 354,689 controls) were included in this study (15).
The summary statistics for FBG, also taken from the MRC-IEU UK Biobank pipeline, included 58,074 individuals of European ancestry (16).
For TG, the IVs were obtained from the MRC-IEU UK Biobank pipeline, which included data from a cohort of over 115,082 individuals (17).
Similar to HDL-C, summary-level statistics were extracted from the MRC-IEU UK Biobank pipeline, which included data from a cohort of over 357,810 subjects of white British descent (14).
2.5 MR analysis method
Inverse variance-weighted (IVW), MR–Egger, weighted median, simple mode, and weighted mode were employed in 2-sample MR analysis to assess the association between education and MetS. The IVW method was applied in the primary analysis, involving the combination of Wald ratios of individual SNPs to derive an overall causal estimate. MR–Egger, weighted median, simple mode, and weighted mode methods were utilized for alternative analyses. MR egger was used for directional horizontal pleiotropy. MR–Egger, IVW methods, Cochran’s Q statistics and funnel plots were employed for the assessment of heterogeneity (10). Sensitivity analysis was performed through leave-one-out (LOO) analysis. TwoSample MR packages in R (version 4.3.3) were used to perform all the statistical analyses. The web tool1 was utilized to calculate the power for MR. p values <0.004 (0.05/6 = 0.008), as determined by Bonferroni correction, were used in this analysis. Significance was defined as p < 0.05 in the MR–Egger test and heterogeneity test.
3 Results
Among the 418 education-associated variants, 3 SNPs were not available for MetS; one SNP was not available for WC, hypertension or HDL-C; and 76 SNPs were not available for FBG, which is a no missing variant in triglycerides. Palindromic SNPs were excluded from the analysis. A total of 404 SNPs were utilized as genetic instruments for MetS, 356 for WC, 406 for hypertension, 333 for FBG, 407 for TG and 356 for HDL-C in this study. The R2 and F statistics were calculated to assess the efficacy of the genetic instruments (Table 1).
The results were shown in Table 2 and Figure 1. The scatter plots were presented in Supplementary Figure S3. The forest plots were presented in Supplementary Figure S4. The results of the IVW analysis [odds ratio (OR) = 0.55, 95% confidence interval (CI) = 0.48–0.63, p = 2.18E−51], MR–Egger regression (OR = 0.70, 95% CI = 0.41–0.99, p = 1.66E−2) and weighted median analysis (OR = 0.58, 95% CI = 0.49–0.67, p = 6.28E−34) indicated each one-s.d.-higher education level was negatively association with MetS (Table 2). The IVW analysis yielded Cochran’s Q test p values (p < 0.001), suggesting substantial heterogeneity in the results. However, the funnel plot did not display asymmetry (Supplementary Figure S1), indicating no evidence of publication bias. Furthermore, MR–Egger intercept values (intercept = −0.003; p = 0.098) indicate the absence of significant horizontal pleiotropy in this study. Furthermore, the LOO analysis demonstrated no significant changes upon the removal of any IVs (Supplementary Figure S2).
Table 2. Geneticly predicted effect of education on the risk of MetS and its components in the MR analysis.
In terms of its constituent elements, each one-s.d.-higher education level was negatively correlated with WC, hypertension, FBG, and triglycerides but positively correlated with HDL-C. In the IVW model, the ORs with 95% CIs per log-odds increase in genetically predicted education level were 0.80 (95% CI: 0.76–0.83; p = 4.98E–33) for WC, 0.96 (95% CI: 0.95–0.97; p = 4.54E–10) for hypertension, 0.94 (95% CI: 0.91–0.97; p = 7.58E–6) for FBG, 0.83 (95% CI: 0.79–0.87; p = 7.87E–18) for triglycerides, and 1.22 (95% CI, 1.18–1.25; p = 1. 45E–31) for HDL-C (see Table 2). The findings were consistent across the various models. Potential horizontal pleiotropy was not found in the IVW model (all p values >0.05), as shown in Table 2. The results of the Cochran Q test revealed significant heterogeneity in WC, hypertension, TG, and HDL-C, but not FBG. Additionally, analysis of the funnel plot (Supplementary Figure S1) did not detect any visible asymmetry, suggesting a lack of evidence for heterogeneity. Furthermore, the stability of the results remained consistent even after the exclusion of any single SNP through LOO analysis. (Supplementary Figure S2).
4 Discussion
This is the first article to elucidate the causal relationship between education and metabolic syndrome, using Mendelian randomization method. In this 2-sample MR study, each addition 4.2 years of education decreased 45% risk of MetS, had a negative association with WC, hypertension, FBG, and TG, a positive association with HDL-C.
Higher education levels are linked to improved health habits and lower rates of MetS, as well as a decreased risk of cardiovascular and cerebrovascular diseases such as heart disease and stroke (18–20). A cross-sectional analysis revealed a negative association between education and MetS (21), which is consistent with the findings of our own study. A study conducted in Korea demonstrated that education level has both direct and indirect impacts on MetS, with higher levels of education indirectly influencing food choices to decrease the prevalence of MetS. Specifically, individuals with higher education levels were more likely to adopt healthier food habits, such as consuming fruits, vegetables, and milk, with this effect being more pronounced in females (21). Previous research has revealed sex differences in the onset and progression of MetS. Interestingly, these differences also extend to the association between education level and MetS, with a more pronounced effect observed in females (8). A study focusing on a cohort of middle-aged women revealed that higher education levels were associated with a decreased prevalence of MetS. Specifically, individuals who had completed college or higher were 2.7 times more likely to have MetS (22). Prior research has generally posited that education is among the factors influencing metabolic syndrome. However, our study indicates that higher levels of education are associated with a reduced incidence of metabolic syndrome.
As the most important component of MetS, WC can also be influenced by education level. In the European cohort, higher levels of education were found to have a significant effect on reducing WC (23). A cross-sectional study conducted on elderly adults in Taipei further supports this finding, showing a positive correlation between obesity and lower education levels, particularly among females, with no significant association with income level (24). These results are consistent with our own study, which revealed a negative relationship between years of education and WC. In a prior study conducted by Jose P. Lopez-Lopez et al., the prevalence of hypertension was significantly higher among populations of African origin than among other populations. Additionally, the study identified an association between the risk of hypertension and lower levels of education (25). Another MR study also identified a lower level of education as a risk factor for hypertension (26). This finding is consistent with the results of our study. A study conducted by Braverman-Bronstein et al. demonstrated that the association between education level and diabetes is modulated by sex and region (27). Specifically, a negative correlation was detected in females, whereas males presented varying correlations across different regions (27). The study revealed a negative correlation in countries such as Argentina, Brazil, Colombia, Chile, and Mexico, whereas a positive correlation was observed in Peru, Panama, and El Salvador (27). A propensity score-based analysis conducted in Bangladesh revealed no significant influence of education level on the prevalence of diabetes and hypertension (28). This finding conflicts with our study. In our study, we observed that a longer education duration was associated with a lower risk of diabetes and hypertension. The discrepancies in findings across studies may be attributed to variations in study populations. A study conducted by Christopher R. Stephens demonstrated that an increase in educational duration was associated with a reduction in WC, blood pressure, FBG, and TG, alongside an increase in HDL-C. However, subgroup analysis revealed that WC increased with education duration exclusively in females (20). Additionally, a cross-sectional analysis of an Australian cohort by L. A. Simons indicated that TG and HDL-C levels were negatively correlated with educational attainment in females but not in males (29). In our study, we observed that extended educational duration was associated with increased WC, FBG, TG, and hypertension risk, as well as decreased HDL-C levels.
5 Strengths
First, MR was used in this study, reducing the influence of confounders. Second, GWAS data include comprehensive and summary-level data for exposure (education) and outcome (MetS, WC, hypertension, FBG, TG, HDL-C) and therefore have high investigative power. Second, multiple statistical methods were used in this study, increasing the credibility of the results.
6 Limitations
First, the influence of IVs on outcomes remains uncertain. Second, owing to the unavailability of data (For example, GWAS data for education is available only for overall population, not separated by sex), subgroup analyses, such as examinations of sex differences, could not be conducted. Finally, significant heterogeneity was identified in this study on the basis of the Cochran’s Q value, but LOO analysis revealed that the result was credible.
7 Conclusion
This study employed MR methodology to investigate the association between levels of education and MetS. A low-level education population should be focused on reducing the incidence of MetS, obesity, hypertension, FBG, TG and HDL-C.
Data availability statement
Publicly available datasets were analyzed in this study. This data can be found at: https://thessgac.com/papers/3; https://gwas.mrcieu.ac.uk/datasets/?trait__icontains=ebi-a-GCST90014020; https://gwas.mrcieu.ac.uk/datasets/ebi-a-GCST90038604/; https://gwas.mrcieu.ac.uk/datasets/ebi-a-GCST005186/; https://gwas.mrcieu.ac.uk/datasets/ebi-a-GCST90092992/; https://gwas.mrcieu.ac.uk/datasets/ebi-a-GCST90014007/; https://www.liebertpub.com/doi/10.1089/met.2019.0070.
Ethics statement
Ethical approval was not required for the study involving humans in accordance with the local legislation and institutional requirements. Written informed consent to participate in this study was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and the institutional requirements.
Author contributions
DL: Data curation, Writing – original draft, Writing – review & editing. Z-XX: Data curation, Formal analysis, Writing – original draft. X-lL: Data curation, Writing – review & editing. H-lY: Writing – review & editing. LW: Formal analysis, Writing – review & editing. YL: Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by the Jilin Provincial Department of Education (JJKH20210062KJ belonging to LW) the Health Commission of Jilin Province (2024A026 belonging to X-lL).
Acknowledgments
All the authors would like to thank all the collaborators involved 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.
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/fnut.2024.1477537/full#supplementary-material
Abbreviations
MetS, Metabolic syndrome; WC, Waist circumference; FBG, Fasting blood glucose; TG, Triglycerides; HDL-C, High density lipoprotein cholesterol levels; MR, Mendelian randomization; SNP, Single nucleotide polymorphism; IVW, inverse variance weighted; IVs, instrumental variables; GWAS, genome-wide association study; IDF, International Diabetes Federation; AHA/NHLBI, American Heart Association/National Heart, Lung, and Blood Institute; LD, linkage disequilibrium; SSGAC, Social Science Genetic Association Consortium; MRC, Medical Research Council; IEU, Integrative Epidemiology Unit; LOO, leave-one-out.
Footnotes
References
1. Zimmet, P, Alberti, KGM, and Ríos, MS. A new international diabetes federation (IDF) worldwide definition of the metabolic syndrome: the rationale and the results. Elsevier Doyma. (2005) 58:1371–5. doi: 10.1016/S1885-5857(06)60742-1
2. Alkhulaifi, F, and Darkoh, C. Meal timing, meal frequency and metabolic syndrome. Nutrients. (2022) 14:1719. doi: 10.3390/nu14091719
3. Cavallari, I, Cannon, CP, Braunwald, E, Goodrich, EL, Im, K, Lukas, MA, et al. Metabolic syndrome and the risk of adverse cardiovascular events after an acute coronary syndrome. Eur J Prev Cardiol. (2018) 25:830–8. doi: 10.1177/2047487318763897
4. Qureshi, D, Collister, J, Allen, NE, Kuźma, E, and Littlejohns, T. Association between metabolic syndrome and risk of incident dementia in UK biobank. Alzheimers Dement. (2024) 20:447–58. doi: 10.1002/alz.13439
5. Li, W, Chen, D, Peng, Y, Lu, Z, Kwan, MP, and Tse, LA. Association between metabolic syndrome and mortality: prospective cohort study. JMIR Public Health Surveill. (2023) 9:e44073. doi: 10.2196/44073
6. Seo, SR, Kim, SY, Lee, S-Y, Yoon, T-H, Park, H-G, Lee, SE, et al. The incidence of stroke by socioeconomic status, age, sex, and stroke subtype: a nationwide study in Korea. J Prev Med Public Health. (2014) 47:104–12. doi: 10.3961/jpmph.2014.47.2.104
7. Winkleby, MA, Jatulis, DE, Frank, E, and Fortmann, SP. Socioeconomic status and health: how education, income, and occupation contribute to risk factors for cardiovascular disease. Am J Public Health. (1992) 82:816–20. doi: 10.2105/AJPH.82.6.816
8. Li, FE, Zhang, FL, Zhang, P, Liu, D, Liu, HY, Guo, ZN, et al. Sex-based differences in and risk factors for metabolic syndrome in adults aged 40 years and above in Northeast China: results from the cross-sectional China national stroke screening survey. BMJ Open. (2021) 11:e038671. doi: 10.1136/bmjopen-2020-038671
9. Katan, M. Apoupoprotein e isoforms, serum cholesterol, and cancer. Lancet. (1986) 327:507–8. doi: 10.1016/S0140-6736(86)92972-7
10. Zhang, M, Chen, J, Yin, Z, Wang, L, and Peng, L. The association between depression and metabolic syndrome and its components: a bidirectional two-sample Mendelian randomization study. Transl Psychiatry. (2021) 11:633. doi: 10.1038/s41398-021-01759-z
11. Alberti, KG, Eckel, RH, Grundy, SM, Zimmet, PZ, Cleeman, JI, Donato, KA, et al. Harmonizing the metabolic syndrome: a joint interim statement of the international diabetes federation task force on epidemiology and prevention; National Heart, Lung, and Blood Institute; American Heart Association; world heart federation; international atherosclerosis society; and International Association for the Study of obesity. Circulation. (2009) 120:1640–5. doi: 10.1161/CIRCULATIONAHA.109.192644
12. Okbay, A, Wu, Y, Wang, N, Jayashankar, H, Bennett, M, Nehzati, SM, et al. Polygenic prediction of educational attainment within and between families from genome-wide association analyses in 3 million individuals. Nat Genet. (2022) 54:437–49. doi: 10.1038/s41588-022-01016-z
13. Lind, L. Genome-wide association study of the metabolic syndrome in UK biobank. Metab Syndr Relat Disord. (2019) 17:505–11. doi: 10.1089/met.2019.0070
14. Mbatchou, J, Barnard, L, Backman, J, Marcketta, A, Kosmicki, JA, Ziyatdinov, A, et al. Computationally efficient whole-genome regression for quantitative and binary traits. Nat Genet. (2021) 53:1097–103. doi: 10.1038/s41588-021-00870-7
15. Dönertaş, HM, Fabian, DK, Valenzuela, MF, Partridge, L, and Thornton, JM. Common genetic associations between age-related diseases. Nat Aging. (2021) 1:400–12. doi: 10.1038/s43587-021-00051-5
16. Manning, AK, Hivert, MF, Scott, RA, Grimsby, JL, Bouatia-Naji, N, Chen, H, et al. A genome-wide approach accounting for body mass index identifies genetic variants influencing fasting glycemic traits and insulin resistance. Nat Genet. (2012) 44:659–69. doi: 10.1038/ng.2274
17. Richardson, TG, Leyden, GM, Wang, Q, Bell, JA, Elsworth, B, Davey Smith, G, et al. Characterising metabolomic signatures of lipid-modifying therapies through drug target mendelian randomisation. PLoS Biol. (2022) 20:e3001547. doi: 10.1371/journal.pbio.3001547
18. Hu, M, Yang, T, and Yang, Y. Causal associations of education level with cardiovascular diseases, cardiovascular biomarkers, and socioeconomic factors. Am J Cardiol. (2024) 213:76–85. doi: 10.1016/j.amjcard.2023.06.044
19. Moustakim, R, Mziwira, M, El Ayachi, M, and Belahsen, R. Association of Metabolic Syndrome and Chronic Kidney Disease in Moroccan adult population. Metab Syndr Relat Disord. (2021) 19:460–8. doi: 10.1089/met.2020.0117
20. Stephens, CR, Easton, JF, Robles-Cabrera, A, Fossion, R, de la Cruz, L, Martínez-Tapia, R, et al. The impact of education and age on metabolic disorders. Front Public Health. (2020) 8:180. doi: 10.3389/fpubh.2020.00180
21. Kim, OY, Kwak, SY, Kim, B, Kim, YS, Kim, HY, and Shin, MJ. Selected food consumption mediates the association between education level and metabolic syndrome in Korean adults. Ann Nutr Metab. (2017) 70:122–31. doi: 10.1159/000470853
22. Wamala, SP, Lynch, J, Horsten, M, Mittleman, MA, Schenck-Gustafsson, K, and Orth-Gomér, K. Education and the metabolic syndrome in women. Diabetes Care. (1999) 22:1999–2003. doi: 10.2337/diacare.22.12.1999
23. Hermann, S, Rohrmann, S, Linseisen, J, May, AM, Kunst, A, Besson, H, et al. The association of education with body mass index and waist circumference in the EPIC-PANACEA study. BMC Public Health. (2011) 11:169. doi: 10.1186/1471-2458-11-169
24. Hsieh, TH, Lee, JJ, Yu, EW, Hu, HY, Lin, SY, and Ho, CY. Association between obesity and education level among the elderly in Taipei, Taiwan between 2013 and 2015: a cross-sectional study. Sci Rep. (2020) 10:20285. doi: 10.1038/s41598-020-77306-5
25. Lopez-Lopez, JP, Cohen, DD, Alarcon-Ariza, N, Mogollon-Zehr, M, Ney-Salazar, D, Chacon-Manosalva, MA, et al. Ethnic differences in the prevalence of hypertension in Colombia: association with education level. Am J Hypertens. (2022) 35:610–8. doi: 10.1093/ajh/hpac051
26. Wang, Y, Ye, C, Kong, L, Zheng, J, Xu, M, Xu, Y, et al. Independent associations of education, intelligence, and cognition with hypertension and the mediating effects of Cardiometabolic risk factors: a Mendelian randomization study. Hypertension. (2023) 80:192–203. doi: 10.1161/HYPERTENSIONAHA.122.20286
27. Braverman-Bronstein, A, Hessel, P, González-Uribe, C, Kroker, MF, Diez-Canseco, F, Langellier, B, et al. Association of education level with diabetes prevalence in Latin American cities and its modification by city social environment. J Epidemiol Community Health. (2021) 75:874–80. doi: 10.1136/jech-2020-216116
28. Gupta, RD, Chakraborty, PA, and Hossain, MB. Association of household wealth and education level with hypertension and diabetes among adults in Bangladesh: a propensity score-based analysis. Trop Med Int Health. (2021) 26:1047–56. doi: 10.1111/tmi.13625
Keywords: education, metabolic syndrome, waist circumference, hypertension, fasting blood glucose, triglycerides, high-density lipoprotein cholesterol
Citation: Liu D, Xu Z-X, Liu X-l, Yang H-L, Wang L-l and Li Y (2024) Education and metabolic syndrome: a Mendelian randomization study. Front. Nutr. 11:1477537. doi: 10.3389/fnut.2024.1477537
Edited by:
Yifan Zhang, Shanghai Jiao Tong University, ChinaReviewed by:
Marilena Stoian, Carol Davila University of Medicine and Pharmacy, RomaniaLannie O’Keefe, Victoria University, Australia
Copyright © 2024 Liu, Xu, Liu, Yang, Wang and Li. 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: Yan Li, eWxpMDFAamx1LmVkdS5jbg==; Ling-ling Wang, NzkyMjU4MzRAMTYzLmNvbQ==
†These authors have contributed equally to this work and share last authorship