- 1Wuxi School of Medicine, Jiangnan University, Wuxi, China
- 2Department of Urology, Jiangnan University Medical Center, Wuxi No. 2 People’s Hospital, Wuxi, China
- 3School of Medicine, Nanjing Medical University, Nanjing, China
- 4School of Medicine, Nantong University, Nantong, China
Background: Kidney stones are a common urologic disease with an increasing incidence year by year, and there are similar influences between gout status and kidney stone incidence. Therefore the contribution of gout status to the incidence of kidney stones is unclear. The aim of this study was to investigate the relationship between gout status and kidney stones and to further explore the causal relationship by Mendelian randomization (MR) analysis.
Method: An epidemiologic study of 49,693 participants in the 2009–2018 National Health and Nutrition Examination Survey (NHANES) was conducted to examine the association between the two. The causal relationship between gout status and kidney stones was assessed by Mendelian randomization analysis of data from the GWAS database.
Result: A total of 28,742 participants were included in the NHANES analysis. We found that gout status was associated with an increased risk of kidney stones [odds ratio (OR) = 1.45 (95%CI, 1.243–1.692); p < 0.001]. In the MR analysis, we found a causal relationship between gout status and the risk of developing kidney stones (OR = 1.047, 95%CI, 1.011–1.085, p = 0.009).
Conclusion: There may be an association between gout status and kidney stone risk. This finding requires further large-sample studies and adequate follow-up.
Introduction
Kidney stones are a kidney disease with a global and increasing incidence (Peerapen and Thongboonkerd, 2023). It has been reported that the incidence of kidney stones in the United States has increased from 3.2% to 9.6% from 1980 to the present (Di X et al., 2023). Kidney stones are a huge economic burden, costing approximately billions of dollars annually (Litwin et al., 2005). Symptomatic kidney stones are more prevalent in men compared with women (Kittanamongkolchai et al., 2018). According to epidemiologic surveys, the highest incidence of kidney stones was found in men after the age of 60, and in women at the age of 50 (Lieske et al., 2006). This suggests a potential correlation between the causes of kidney stones and gender and age. The occurrence of kidney stones not only reduces the quality of life of the patient but can also be life-threatening (Bryant et al., 2012; Whitehurst et al., 2019). Although there are many treatment strategies available for kidney stones, it is more important to prevent their occurrence.
It is well known that there is a significant correlation between gout and kidney stones (Kim et al., 2017). Gout is a chronic inflammatory disease caused by excessive deposition of uric acid crystals in joints and non-joint structures (Dalbeth et al., 2021). Meanwhile, gout is the most common inflammatory arthritis (Keller and Mandell, 2021). And the prevalence and disability of gout worldwide is increasing yearly (McCormick et al., 2022). There is a correlation between gout and the development of several kidney diseases, and a correlation has been reported between elevated urate levels and the development of kidney stones (Ramos and Goldfarb, 2022). This suggests that gout may play a crucial role in the prevention of kidney stones.
Mendelian randomization (MR) is an emerging and valuable approach to causality research. It uses single nucleotide polymorphisms (SNPs) as instrumental variables to the causal relationship between inflammatory exposures (e.g., gout) and outcomes (e.g., kidney stones). Because gene alleles are randomly assigned during meiosis independent of environmental factors, genetic associations observed from MR analyses are unlikely to be affected by confounding bias and risk of reverse causation (Li et al., 2023). Because the method is effective in avoiding confounding factors from interfering with the analyzed results, MR is widely used for causal hypothesis studies in epidemiology.
In this study, we investigate the causal relationship between gout and kidney stones by combining the National Health and Nutrition Examination Survey (NHANES) and MR analysis.
Materials and methods
Study population in NHANES
The National Health and Nutrition Examination Survey (NHANES) is a 2-year cross-sectional study designed to investigate the health status of the U.S. population. Data on sociodemographic characteristics, diseases, and related aspects are obtained through questionnaires and laboratory data. Gout was labeled in the database as “Doctor ever told you that you had gout?”. Kidney stones were labeled as “Ever had kidney stones? “To make our findings more accurate, we selected data from 2009–2018 (N = 49,693) for analysis. In this study, we excluded patients with “don't know” and “missing data” from the gout data, and then excluded patients with “don't know” from the kidney stones, and finally included the remaining 28,742 patients in our study (Figure 1).
Figure 1. Flow chart for inclusion of participants based on inclusion and exclusion criteria. Schematic diagram showing study participants included for the present analysis from the 2009 to 2018 NHANES. A total of 28,742 participants were included. NHANES, National Health and Nutrition Examination Survey.
Covariates used in NHANES
To control for potential confounders, we chose the following demographic characteristics: sex, age, race, education level, BMI (body mass index), smoking, alcohol consumption, sleep disorders, uric acid levels, and water consumption. These variables were chosen because of their possible association with outcomes. Participants were categorized as smokers and nonsmokers based on whether they had smoked fewer than 100 cigarettes in their lifetime. Participants were categorized as drinkers, nondrinkers, refused to answer, and don't know based on whether they had more than 12 drinks in their lifetime. Participants were categorized as sleep-disordered and non-sleep-disordered based on whether they had a sleep disorder. Uric acid levels were categorized as < 4 mg/dL, 4–7 mg/dL, and > 7 mg/dL based on circulating uric acid levels.
Data sources for exposure and outcomes
Genetic variants strongly associated with gout were obtained from the IEU OPEN GWAS, and we chose the gout number “ukb-a-107”. A total of 337,159 investigators (332,352 controls and 4,807 cases) were included in this data. Genetic variants associated with kidney stones were obtained from UK biobank, number “ukb-a-72”, which contained 337,159 investigators (334,465 controls and 2,694 cases).
Selection of genetic instruments
Instrumental variables (single nucleotide polymorphisms, SNPs) were screened by the following criteria: 1) significant gene-wide correlation (p < 5 × 10−8); 2) pruned by linkage disequilibrium (r2 < 0.001 and within 10,000 kb from the index variant).
Statistical analysis
In our observational study, we selected appropriate data in the NHANES database and used multivariate logistic regression to assess the association between gout status and kidney stone risk. In addition, we used the sample weights provided by NHANES in our observational study.
The validity of the results of the MR analysis depends on the key assumption of no pleiotropy (Lawlor et al., 2008). First, we applied the MR-Egger method, which is based on the Instrument Strength Independent of Direct Effect assumption. In this method, the effect of SNPs on gout is charted against the effect of SNPs on kidney stone, and an intercept dissimilar from that of the source is considered as proof of pleiotropic effects (Bowden et al., 2015; Bowden et al., 2015). This method yields bias-free estimations even if all the selected SNPs are unfounded. Furthermore, we used three more MR models that are more robust under horizontal pleiotropy conditions: simple mode, weighted mode, weighted median (Burgess et al., 2016; Li et al., 2023; Lin et al., 2023). We used Cochrane’s Q test to test for possible SNP heterogeneity. We also performed omission analysis and MR-Egger intercept tests to assess significant estimates of horizontal pleiotropy. MR analyses were performed using the R package “TwoSampleMR”, recognizing outliers with the “MR-PRESSO” package (Burgess et al., 2016; Li et al., 2023; Lin et al., 2023). All results are reported as odd ratios, with 95% confdence intervals (CI). All statistical analysis and data visualization were performed in R software 4.2.0 (https://www.r-project.org/).
Result
Observations on gout and kidney stones in NHANES
By selection, 28,742 participants were included in our study (Figure 1). Table 1 shows the Baseline characteristics of participants by gout in the NHANES 2009–2018. The age range of the participants was 20–80 years old, with 51.6% females and 48.4% males. Gout was diagnosed by doctors in 3.2% of the men and 1.5% of the women. The occurrence of gout was associated with age (64 vs. 49 years; p < 0.001) and was more likely to be associated with males (67.6% vs. 32.4%; p < 0.001) and non-Hispanic whites (46.4% vs. 38.8%; p < 0.001), and with a BMI > 30 kg/m2 (53.8% vs. 38.0%; p < 0.001). In addition, the prevalence of kidney stones (16.6% vs. 9.0%; p = 0.004) and whether sleep difficulties (39% vs. 25.2%; p < 0.001) were significantly higher in gout patients than in non-gout patients among all participants. Gout patients had higher rates of uric acid > 7 mg/dL (p < 0.001), smoking (p < 0.001), alcohol consumption (p < 0.001) and sleep difficulties (p < 0.001) compared to non-gout patients. Gout patients had higher rates of uric acid > 7 mg/dL (p < 0.001), smoking (p < 0.001), alcohol consumption (p < 0.001), and sleep difficulties (p < 0.001) compared to healthy individuals. Among patients with gout, 16.6% had a history of kidney stones. In contrast, the prevalence of gout among patients reporting kidney stones was 8.4%.
In the age-adjusted model, kidney stones were significantly associated with gout in the study population [odds ratio (OR) = 1.45; 95%CI, 1.243–1.692] (Table 2). In the multivariate model, we found that gout was significantly and positively associated with kidney stones (OR = 1.284; 95%CI, 1.084–1.522) (Table 3). However, in the multivariate model, the result for educational level was not statistically significant (p = 0.738 > 0.05), so the difference in educational level was not significant among gout patients.
Causal association between gout and kidney stones in MR
After data merging preprocessing, we selected a total of 22 SNPs for MR analysis (Figure 2A). We performed 5 MR analyses, in which the results of IVW, Weighted media, and Weighted mode were significant (p < 0.05) (Table 4). IVW showed that there was a causal relationship between gout and kidney stones (OR = 1.047; 95%CI,1.011–1.085, p = 0.009). Weighted media (OR = 1.065; 95%CI,1.016–1.116, p = 0.007) and Weighted mode (OR = 1.067; 95%CI, 1.010–1.128, p = 0.030) as sensitivity analyses similarly confirmed this. In contrast, MR Egger (OR = 1.044; 95%CI, 0.983–1.110, p = 0.174) and Simple mode (OR = 1.072, 95%CI, 0.968–1.187, p = 0.190) turned out to be non-significant (p > 0.05). However, the IVW, which provided a precise causal estimate, had a p = 0.009 < 0.05, and the beta values of all five tests were >0 in the same direction (Figure 2B), i.e., individuals who developed gout were at increased risk of kidney stones, which suggests that there is indeed a causal relationship between gout and kidney stones. In order to exclude the heterogeneity of SNPs affecting the results causing bias, we subsequently performed a heterogeneity analysis (Table 5), using Cochran’s Q test to assess the degree of heterogeneity of individual effects on the results, in which the Q test P for IVW (Q = 27, p = 0.145) as well as for MR-Egger (Q = 27, p = 0.113) values were > 0.05 suggesting that there is no heterogeneity,the effect of gout on kidney stones does not vary from individual to individual. Funnel plot results show that SNPs are symmetric, suggesting stable results (Figure 2C). In the test of multiplicity (Table 6), we verified the MR-Egger (except = 1.26e-05, p = 0.910), p = 0.910 > 0.05 proved the weak possibility of horizontal multiplicity, i.e., SNP directly affects the risk of renal stone but not gout, which can be ignored. After passing the tests for heterogeneity and pleiotropy, we performed leave-one-out analysis (Figure 2D) and observed that the gradual elimination of each SNP had little effect on the final results, with little change in the overall error line, establishing that it was not a single SNP that was responsible for the association between gout and kidney stone risk.
Figure 2. Causal relationship between gout and kidney stones. (A) Forest map of gout and kidney stones. (B) Scatter plot of gout and kidney stones. (C) Funnel plot of gout and kidney stones. (D) Results of the “leave-one-out” method in gout and kidney stones.
Discussion
In this study, we conducted the first large-scale investigation of the association between gout status and kidney stone risk based on the NHANES database and MR analysis. We investigated the relationship between gout status and kidney stone risk and the potential causal relationship between them. Our results suggest that gout status is positively associated with kidney stone risk and that there is a causal relationship between them.
Previously, some researchers have explored the causal relationship between serum uric acid and the risk of urolithiasis. Ravi K. Narang et al. conducted a large-scale MR analysis comprising 359,827 participants, and the results showed that the existence of a causal relationship between serum uric acid and the risk of urolithiasis was not supported (Narang et al., 2021). However, in our cross-sectional study (2009–2018) based on the NHANES database, we found a significant positive association between gout status and kidney stone risk (p = 0.004). A significant correlation between patients with urinary tract stones and gout was found based on a previous case-finding study, which is consistent with our findings (Marchini et al., 2013).
The diagnosis of gout relies heavily on circulating levels of uric acid, and the most common components of urinary stones are: calcium oxalate, carbapatite, urate, struvite, and brushite (Ludwig and Matlaga, 2018; Ye et al., 2020). In the present study, our large investigative study in NHANES and MR analyses went on to explore whether gout status is a risk factor for kidney stones. Although the correlation between gouty status and kidney stones was demonstrated in our study, correlation studies between the two are scarce and further research is needed. In our MR analysis, more than half of the analytical methods supported a causal relationship between gout status and high risk of kidney stones. In our study, the result of MR-Egger (OR = 1.044, 95%CI, 0.983–1.110, p = 0.174) was not significant (p > 0.05), which suggests that horizontal pleiotropy does not exist in this study (Luo et al., 2021). Although the result of Simple mode (OR 1.072, 95%CI, 0.968–1.187, p = 0.190) was also not significant (p > 0.05), the results of our study were based on IVW(Yang et al., 2023).
Although the definitive mechanism between gout and kidney stones is not known, we discuss several possible hypotheses. First, gout patients have circulating levels of uric acid that are significantly higher than those of healthy individuals, and acidification of the urine due to high uric acid levels favors the formation of kidney stones (Ferraro et al., 2020). Second, the disease state of gout patients is accompanied by an unpleasant state, and the persistence or frequency of the unpleasant state leads to the development of kidney stones (Wang et al., 2022). A previous study conducted a meta-analysis of the correlation between gout and depression and anxiety, which showed a positive correlation between gout and mental health disorders (Howren et al., 2021). Meanwhile, a study reported a positive correlation between the degree of depression and the risk of kidney stones (Wang et al., 2022). However, the clear mechanism of the related studies still needs further research.
The main strength of this study was the use of observational studies from 2009–2018 in the NHANES database combined with MR analysis. The combined collation of many factors and the large sample size enabled us to adjust the analysis of confounders in the regression model in a timely manner and with sufficient computational power to analyze the causal relationship between gout and kidney stones. Observational studies alone can only illustrate the correlation between observed factors, whereas MR analysis further analyzes the causal relationship between observed factors on the basis of correlation analysis, and MR analysis can avoid the bias of results caused by confounding factors (Zuber et al., 2022). However, there are some limitations of our study. First, gout-related data in the NHANWS database were incomplete in recent years due to the impact of neocoronary pneumonia, which affected the further refinement of the results analyzed. Second, the study population was limited to the Americas and Europe, which may have had an impact on the generalizability of the results.
Conclusion
In conclusion, our study demonstrated a strong positive correlation between gout status and the risk of developing kidney stones. And our MR analysis results also proved the causal relationship between gout status and kidney stones. However, our findings still need to be further validated and their underlying mechanisms need to be further investigated.
Data availability statement
The data that support the findings of this study are openly available at https://wwwn.cdc.gov/nchs/nhanes/Default.aspx and http://gwas.mrcieu.ac.uk
Author contributions
SJ: Data curation, Investigation, Validation, Visualization, Writing–original draft. HG: Formal Analysis, Supervision, Validation, Writing–review and editing. YL: Data curation, Validation, Writing–review and editing. YhZ: Investigation, Validation, Writing–review and editing. JL: Data curation, Project administration, Software, Writing–review and editing. CF: Validation, Writing–review and editing. YwZ: Supervision, Writing–review and editing. ML: Methodology, Supervision, Writing–review and editing. FN: Conceptualization, Funding acquisition, Writing–review and editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This study was supported by the National Natural Science Foundation of China (Item No: 82370777, 82302654), Major Project of Jiangsu Commission of Health: Basic research and clinical application of biomarkers for early diagnosis of prostate cancer (Item No: ZD2021002), National Natural Science Foundation of Jiangsu Province (No. SBK2023044153), “San Ming” Strategic Talent Project of Wuxi Municipal Health Commission (No. HB2023032), Science and Technology Research Project of Wuxi Science and Technology Bureau (No K20221029).
Acknowledgments
We thank all the contributors and participants in the National Health and Nutrition Examination Survey and the original genome-wide association studies.
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.
References
Bowden, J., Davey, S. G., and Burgess, S. (2015). Mendelian randomization with invalid instruments: effect estimation and bias detection through egger regression. Int. J. Epidemiol. 44 (2), 512–525. doi:10.1093/ije/dyv080
Bryant, M., Angell, J., Tu, H., Goodman, M., Pattaras, J., and Ogan, K. (2012). Health related quality of life for stone formers. J. Urol. 188 (2), 436–440. doi:10.1016/j.juro.2012.04.015
Burgess, S., Dudbridge, F., and Thompson, S. G. (2016). Combining information on multiple instrumental variables in mendelian randomization: comparison of allele score and summarized data methods. Stat. Med. 35 (11), 1880–1906. doi:10.1002/sim.6835
Dalbeth, N., Gosling, A. L., Gaffo, A., and Abhishek, A. (2021). Gout. Lancet. 397 (10287), 1843–1855. doi:10.1016/S0140-6736(21)00569-9
Di, X., Liu, S., Xiang, L., and Jin, X. (2023). Association between the systemic immune-inflammation index and kidney stone: a cross-sectional study of nhanes 2007-2018. Front. Immunol. 14, 1116224. doi:10.3389/fimmu.2023.1116224
Ferraro, P. M., Bargagli, M., Trinchieri, A., and Gambaro, G. (2020). Risk of kidney stones: influence of dietary factors, dietary patterns, and vegetarian-vegan diets. Nutrients 12 (3), 779. doi:10.3390/nu12030779
Howren, A., Bowie, D., Choi, H. K., Rai, S. K., and De Vera, M. A. (2021). Epidemiology of depression and anxiety in gout: a systematic review and metaanalysis. J. Rheumatol. 48 (1), 129–137. doi:10.3899/jrheum.190974
Keller, S. F., and Mandell, B. F. (2021). Management and cure of gouty arthritis. Med. Clin. N. Am. 105 (2), 297–310. doi:10.1016/j.mcna.2020.09.013
Kim, S., Chang, Y., Yun, K. E., Jung, H. S., Lee, S. J., Shin, H., et al. (2017). Development of nephrolithiasis in asymptomatic hyperuricemia: a cohort study. Am. J. Kidney Dis. 70 (2), 173–181. doi:10.1053/j.ajkd.2017.01.053
Kittanamongkolchai, W., Vaughan, L. E., Enders, F. T., Dhondup, T., Mehta, R. A., Krambeck, A. E., et al. (2018). The changing incidence and presentation of urinary stones over 3 decades. Mayo Clin. Proc. 93 (3), 291–299. doi:10.1016/j.mayocp.2017.11.018
Lawlor, D. A., Harbord, R. M., Sterne, J. A., Timpson, N., and Davey, S. G. (2008). Mendelian randomization: using genes as instruments for making causal inferences in epidemiology. Stat. Med. 27 (8), 1133–1163. doi:10.1002/sim.3034
Li, W., Zheng, Q., Xu, M., Zeng, C., and Deng, X. (2023). Association between circulating 25-hydroxyvitamin d metabolites and periodontitis: results from the nhanes 2009-2012 and mendelian randomization study. J. Clin. Periodontol. 50 (2), 252–264. doi:10.1111/jcpe.13736
Lieske, J. C., Pena, D. L. V. L., Slezak, J. M., Bergstralh, E. J., Leibson, C. L., Ho, K. L., et al. (2006). Renal stone epidemiology in rochester, Minnesota: an update. Kidney Int. 69 (4), 760–764. doi:10.1038/sj.ki.5000150
Lin, Z., Huang, J., Xie, S., Zheng, Z., Tang, K., Li, S., et al. (2023). The association between insulin use and asthma: an epidemiological observational analysis and mendelian randomization study. Lung 201 (2), 189–199. doi:10.1007/s00408-023-00611-z
Litwin, M. S., Saigal, C. S., Yano, E. M., Avila, C., Geschwind, S. A., Hanley, J. M., et al. (2005). Urologic diseases in America project: analytical methods and principal findings. J. Urol. 173 (3), 933–937. doi:10.1097/01.ju.0000152365.43125.3b
Ludwig, W. W., and Matlaga, B. R. (2018). Urinary stone disease: diagnosis, medical therapy, and surgical management. Med. Clin. N. Am. 102 (2), 265–277. doi:10.1016/j.mcna.2017.10.004
Luo, J., le Cessie, S., van Heemst, D., and Noordam, R. (2021). Diet-derived circulating antioxidants and risk of coronary heart disease: a mendelian randomization study. J. Am. Coll. Cardiol. 77 (1), 45–54. doi:10.1016/j.jacc.2020.10.048
Marchini, G. S., Sarkissian, C., Tian, D., Gebreselassie, S., and Monga, M. (2013). Gout, stone composition and urinary stone risk: a matched case comparative study. J. Urol. 189 (4), 1334–1339. doi:10.1016/j.juro.2012.09.102
McCormick, N., Lu, N., Yokose, C., Joshi, A. D., Sheehy, S., Rosenberg, L., et al. (2022). Racial and sex disparities in gout prevalence among us adults. JAMA Netw. Open 5 (8), e2226804. doi:10.1001/jamanetworkopen.2022.26804
Narang, R. K., Gamble, G. G., Topless, R., Cadzow, M., Stamp, L. K., Merriman, T. R., et al. (2021). Assessing the relationship between serum urate and urolithiasis using mendelian randomization: an analysis of the UK biobank. Am. J. Kidney Dis. 78 (2), 210–218. doi:10.1053/j.ajkd.2020.11.018
Peerapen, P., and Thongboonkerd, V. (2023). Kidney stone prevention. Adv. Nutr. 14 (3), 555–569. doi:10.1016/j.advnut.2023.03.002
Ramos, G. K., and Goldfarb, D. S. (2022). Update on uric acid and the kidney. Curr. Rheumatol. Rep. 24 (5), 132–138. doi:10.1007/s11926-022-01069-3
Wang, M., Jian, Z., Ma, Y., Jin, X., Li, H., and Wang, K. (2022). Depression increases the risk of kidney stone: results from the national health and nutrition examination survey 2007-2018 and mendelian randomization analysis. J. Affect. Disord. 312, 17–21. doi:10.1016/j.jad.2022.06.008
Whitehurst, L., Jones, P., and Somani, B. K. (2019). Mortality from kidney stone disease (ksd) as reported in the literature over the last two decades: a systematic review. World J. Urol. 37 (5), 759–776. doi:10.1007/s00345-018-2424-2
Yang, H., Song, J., Li, A., Lv, L., Sun, X., Mao, Y., et al. (2023). Genetically predicted levels of folate, vitamin b(12), and risk of autoimmune diseases: a mendelian randomization study. Front. Immunol. 14, 1139799. doi:10.3389/fimmu.2023.1139799
Ye, Z., Zeng, G., Yang, H., Li, J., Tang, K., Wang, G., et al. (2020). The status and characteristics of urinary stone composition in china. BJU Int. 125 (6), 801–809. doi:10.1111/bju.14765
Keywords: gout, kidney stone, national health and nutrition examination survey, cross-sectional research, mendelian randomization
Citation: Jin S, Geng H, Lu Y, Zhou Y, Lv J, Fu C, Zhang Y, Li M and Feng N (2024) Association between gout and kidney stone: results from mendelian randomization and the NHANES study. Front. Genet. 15:1417663. doi: 10.3389/fgene.2024.1417663
Received: 19 April 2024; Accepted: 15 October 2024;
Published: 24 October 2024.
Edited by:
Hong Weng, Wuhan University, ChinaReviewed by:
Yu Wang, Vanderbilt University Medical Center, United StatesRicardo Adrian Nugraha, Airlangga University, Indonesia
Copyright © 2024 Jin, Geng, Lu, Zhou, Lv, Fu, Zhang, Li and Feng. 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: Ninghan Feng, bi5mZW5nQG5qbXUuZWR1LmNu; Menglu Li, bG1sYncxMjMyQDE2My5jb20=; Yuwei Zhang, emhhbmd5dzAxMjNAbmptdS5lZHUuY24=
†These authors have contributed equally to this work and share first authorship