Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.
Advertisement
Scientific Reports volume 12, Article number: 17630 (2022)
246
3
Metrics details
Non-alcoholic fatty liver disease is the most common chronic liver disease and is associated with chronic kidney disease. The fibrosis-4 index and non-alcoholic fatty liver disease score are widely used as non-invasive diagnostic methods for non-alcoholic fatty liver disease. However, the relationship between these markers and specific renal histopathologies in chronic kidney disease remain unclear. This study included 179 patients aged between 16 and 80 years who underwent renal biopsy. We examined the association between the fibrosis-4 index or non-alcoholic fatty liver disease score and change in estimated glomerular filtration rate 12 months after kidney biopsy for each renal histopathology. Renal histopathologies were determined by renal biopsy. Our results showed that there was a significant negative correlation between the fibrosis-4 index and estimated glomerular filtration rate. In nephrosclerosis, the non-alcoholic fatty liver disease score and estimated glomerular filtration rate tended to have a negative correlation, albeit without significance. In IgA nephropathy, both the fibrosis-4 index and non-alcoholic fatty liver disease score were significantly negatively correlated with estimated glomerular filtration rate. Furthermore, the fibrosis-4 index was not associated with urinary protein-to-creatinine ratio or renal function markers such as urinary b2 microglobulin and urinary N-acetyl-d-glucosamine. Our kidney biopsy-based study showed that the liver fibrosis markers fibrosis-4 index and non-alcoholic fatty liver disease score were negatively correlated with the estimated glomerular filtration rate in nephrosclerosis and IgA nephropathy.
Chronic kidney disease is associated with increased risks of mortality and morbidity. Furthermore, CKD is also associated with diabetes mellitus, including diabetic kidney disease (DKD)1. Obesity and insulin resistance are the most common causes of CKD2,3. It has been reported that non-alcoholic fatty liver disease (NAFLD) is a clinicopathological condition that encompasses a variety of chronic liver diseases ranging from asymptomatic hepatic fat accumulation to progressive liver disease, with clinical symptoms resembling alcoholic liver disease despite the absence of excessive alcohol consumption4,5. A recent study revealed that the prevalence of NAFLD diagnosed by ultrasound was estimated to be about 25%6. It is estimated that either NAFLD or non-alcoholic steatohepatitis (NASH) will be the most common indication for liver transplantation by 20307. NAFLD is a multisystem disease that is associated with various diseases such as diabetes mellitus, obesity, metabolic syndrome, cardiovascular diseases, and CKD8. The fibrosis-4 (FIB-4) index and NAFLD fibrosis score (NFS) are simple and validated diagnostic markers used to stratify the degree of liver fibrosis according to risk in patients with suspected NAFLD9. Moreover, these non-invasive fibrosis markers have been shown to be associated with the prevalence of CKD in recent studies10. However, there have been no reports on the association between these fibrosis markers and decreased estimated glomerular filtration rate (eGFR) by renal histopathology classification as determined by kidney biopsy. Thus, we performed this renal pathology-based study to elucidate whether the FIB-4 index and NFS can predict eGFR decline in CKD patients.
We determined the clinical characteristics of the 179 cases after excluding patients who underwent kidney biopsy. The median age of the study participants was 51 years, and 96 patients (53.6%) were female. The median body mass index was 22.7, and the median FIB-4 index was 0.96. The other median values are as follows: systolic blood pressure (SBP), 129 mmHg; diastolic blood pressure (DBP), 78 mmHg; serum creatinine (Cr), 76.02 mmol/L; blood urea nitrogen (BUN), 5 mmol/L; estimated glomerular filtration rate (eGFR), 64 mL/min/1.73 m2; urine protein-to-creatinine ratio, 0.538; N-acetyl-β-d-glucosaminidase (urinary NAG), 8.9 IU/L; and urinary β2 microglobulin (urinary β2 MG), 198 mg/L (Table 1). The underlying diseases are as follows: nephrosclerosis (10 cases, 5.6%), IgA nephropathy (74 cases, 41.3%), minimal change disease (24 cases, 13.4%), diabetic nephropathy (7 cases, 3.9%), lupus nephritis (12 cases, 6.7%), membranous nephropathy (18 cases, 10.1%), ANCA-associated glomerulonephritis (9 cases, 5.0%), interstitial nephritis (9 cases, 5.0%), and others (16 cases, 8.9%) (Table 2).
We calculated the FIB-4 index and eGFR before renal biopsy for all diseases. The results showed that there was a weak but negative correlation between the FIB-4 index and eGFR (R2 = 0.1458, P < 0.01, Fig. 1a). In contrast, the FIB-4 index was not correlated with the urinary protein-to-creatinine ratio (R2 = 0.0164, P = 0.08, Fig. 1b) and urinary β2 MG (R2 = 0.0819, P = 0.08, Fig. 1c). Furthermore, there was a negative correlation between the FIB-4 index and urinary NAG (R2 = 0.0483, P < 0.01, Fig. 1d).
Associations of the FIB-4 index and renal function markers. (a) Correlation between the FIB-4 index and estimated glomerular filtration rate. (b) Correlation between the FIB-4 index and urine protein-creatinine ratio. (c) Correlation between the FIB-4 index and urinary b2 microglobulin. (d) Correlation between the FIB-4 index and urinary N-acetyl-b-d-glucosaminidase.
Next, we examined the relationship between FIB-4 index and NFS, and the rate of change in eGFR before and 12 months after renal biopsy. In nephrosclerosis, we found that the FIB-4 index was significantly negatively correlated with the eGFR (R2 = 0.4362, P = 0.04, Fig. 2a). Additionally, the correlation between NFS and eGFR tended to be negative, albeit without statistical significance (R2 = 0.1503, P = 0.27, Fig. 2b). In IgA nephropathy, the FIB-4 index (R2 = 0.1282, P < 0.01, Fig. 2c) and NFS (R2 = 0.2245, P < 0.01, Fig. 2d) were significantly negatively correlated with the eGFR. The rate of change in eGFR was not negatively associated with either FIB-4 or NFS in lupus nephritis (R2 = 0.2939, P = 0.07; R2 = 0.4573, P = 0.02, respectively; Fig. 3a,b), diabetic nephropathy (R2 = 0.0854, P = 0.52; R2 = 0.1894, P = 0.33, respectively; Fig. 3c,d), membranous nephropathy (R2 = 0.0307, P = 0.49; R2 = 0.1213, P = 0.16, respectively; Fig. 3e,f), minimal change disease (R2 = 0.0541, P = 0.30; R2 = 0.0132, P = 0.61, respectively; Fig. 4a,b), ANCA-associated glomerulonephritis (R2 = 0.0022, P = 0.91; R2 = 0.0381, P = 0.61, respectively; Fig. 4c,d), or interstitial nephritis (R2 = 0.0933, P = 0.46; R2 = 0.721, P = 0.52, respectively; Fig. 4e,f).
Associations between the change in estimated glomerular filtration rate and fibrosis markers in nephrosclerosis and IgA nephropathy. (a,b) Correlation between the change in the estimated glomerular filtration rate and the FIB-4 index (a) and NFS (b) in nephrosclerosis. (c,d) Correlation between the change in the estimated glomerular filtration rate and the FIB-4 index (c) and NFS (d) in IgA nephropathy.
Associations between the change in estimated glomerular filtration rate and fibrosis markers in lupus nephritis, diabetic nephropathy, and membranous nephropathy. (a,b) Correlation between the change in the estimated glomerular filtration rate and the FIB-4 index (a) and NFS (b) in lupus nephritis. (c,d) Correlation between the change in the estimated glomerular filtration rate and the FIB-4 index (c) and NFS (d) in diabetic nephropathy. (e,f) Correlation between the change in the estimated glomerular filtration rate and the FIB-4 index (e) and NFS (f) in membranous nephropathy.
Associations between the change in estimated glomerular filtration rate and fibrosis markers in minimal change disease, ANCA-associated glomerulonephritis, and interstitial nephritis. (a,b) Correlation between the change in the estimated glomerular filtration rate and the FIB-4 index (a) and NFS (b) in minimal change disease. (c,d) Correlation between the change in the estimated glomerular filtration rate and the FIB-4 index (c) and NFS (d) in ANCA-associated glomerulonephritis. (e,f) Correlation between the change in the estimated glomerular filtration rate and the FIB-4 index (e) and NFS (f) in interstitial nephritis.
Next, we have checked in the high FIB-4 index group (≥ 1.3). In nephrosclerosis, we found that the FIB-4 index was more significantly negatively correlated with the rate of change in eGFR when examined in all cases including the low FIB-4 index (< 1.3) (R2 = 0.8115, P < 0.01, Fig. 5a). A similar trend was observed in interstitial nephritis, but the difference was not statistically significant (R2 = 0.4291, P = 0.29, Fig. 5b). However, other diseases did not show any changes when classified by high FIB-4 index levels.
Associations between the change in estimated glomerular filtration rate and FIB-4 index in nephrosclerosis and IgA nephropathy of high FIB-4 index group. (a) Nephrosclerosis in high FIB-4 index group. (b) IgA nephropathy in high FIB-4 index group.
The association between the rate of change in eGFR and the levels of Cr, BUN, and urinary β2 MG before kidney biopsy was revealed by multivariate analysis (Table 3). As shown in Table 4, multivariate analysis showed that the rate of change in eGFR and FIB-4 index were independently associated in both nephrosclerosis and IgA nephropathy, even adjusting for relevant factors.
In the present study, we examined the relationship between each renal histopathology as determined by biopsy and non-invasive fibrosis markers. In nephrosclerosis, we found a statistically significant negative correlation between the FIB-4 index and eGFR, whereas the NFS and eGFR tended to have a negative correlation, albeit insignificant. Furthermore, in IgA nephropathy, both the NFS and FIB-4 index had statistically significant negative correlations with eGFR. Meanwhile, the pre-renal biopsy eGFR and FIB-4 index were negatively correlated in all subjects. This is consistent with previous reports, since patients eligible for renal biopsy suffered from CKD.
Our results showed that there was an association, albeit very weak, between the FIB-4 index and changes in eGFR. This may be due to the wide range of eGFR (1–100 ml/min/1.73 m2 or even higher), resulting in greater variability.
Recent studies have indicated that NAFLD could affect CKD, including DKD11. The FIB-4 index and NFS are widely used hepatic fibrosis markers for the diagnosis of NAFLD or NASH. Furthermore, the FIB-4 index could be useful for the prediction of the development of cardiovascular diseases12. It has been reported that the FIB-4 index is an indicator of inflammation; FIB-4 index is calculated using the platelet count, which contributes to thrombotic inflammatory actions owing to their ability to functionally interact with activated endothelial cells, leukocytes, and coagulation-related proteins13. Several studies have shown that inflammation is involved in the development of CKD14,15. As previously reported, nuclear factor NF-kappa-B (NF-kB) acts via inflammation and protein kinase C to activate the expression of cytokines; NF-kB is transferred to the nucleus, promptly activating the subsequent transcription of tumor necrosis factor-α, vascular cell adhesion molecule 1, or interleukin-6 in the kidneys of CKD patients3,16,17.
Previous reports have indicated that the activation of the TGF-β/Smad pathway is associated with the FIB-4 index and NFS18. We reported that Smad1 signaling plays a significant role in increasing the extracellular matrix, such as type 4 collagen and α-smooth muscle actin, in mesangial cells19. Inflammatory processes inevitably cause damage; as part of the healing process, glomerulosclerosis increases the mesangial extracellular matrix, leading to a decreased glomerular filtration rate14,20. Our results showed that the FIB-4 index and NFS were negatively correlated with eGFR in nephrosclerosis and IgA nephropathy, which are both characterized by an increase in mesangial extracellular matrix.
Obesity is the most common risk factor for NAFLD and an independent risk for developing CKD21. Furthermore, we have reported the mechanisms of developing obesity-related CKD3,14,15,16,22. NAFLD is one of the most common causes of chronic liver disease in clinical practice in its entire spectrum of conditions, ranging from simple lipidosis to steatohepatitis and cirrhosis23,24,25. It has been reported that liver fibrosis and fatty liver are crucial in the development of DKD26. In contrast, another study using multivariate analysis showed that the presence of NAFLD was not a significant predictor of the development of DKD11. In this regard, the association between NAFLD and DKD remains inconsistent. Although patients in our study tended to be obese (median value of body mass index: 22.7), we could not find any correlation between the FIB-4 index and eGFR in patients with type 2 diabetes.
Unlike eGFR, there was no clear correlation between urinary protein or tubulointerstitial markers and FIB-4 index. It has been reported that the FIB-4 index may be associated with vascular endothelial function27, which is a mechanism of glomerulosclerosis caused by glomerular endothelial dysfunction3,17. Thus, the FIB-4 index may be more reflective of glomerulosclerotic lesions. In contrast, proteinuria, which is mainly caused by podocyte dysfunction, may not be associated with the FIB-4 index, which is associated with glomerular endothelial cells28,29.
This study has several limitations. First, liver diseases, such as NASH or NAFLD, were not evaluated by imaging and liver biopsy. Second, we conducted a retrospective analysis of a cohort from a single center. Finally, some diseases, especially nephrosclerosis, had a small number of cases.
The FIB-4 index and NFS, which are non-invasive fibrosis markers, could be predictive markers of reduced eGFR, especially in nephrosclerosis and IgA nephropathy.
All procedures performed in studies involving human participants were in accordance with the ethical standards of the institutional and/or national research committee and with the 1964 Helsinki Declaration and its later amendments or comparable ethical standards. The ethics committee of Osaka Medical and Pharmaceutical University approved this study (approval number: 2020-095). Written informed consent was obtained from all subjects. This study included a retrospective cohort of 179 cases who underwent kidney biopsy at Osaka Medical and Pharmaceutical University Hospital between April 2015 and July 2020. Data were collected and analyzed retrospectively by using electronic medical records maintained by the Department of Nephrology at Osaka Medical and Pharmaceutical University Hospital. Cr, BUN, eGFR, urine protein-to-creatinine ratio, urinary NAG, urinary β2 MG, and patient characteristics (e.g., underlying diseases, age, sex, blood pressure, and body mass index) were searched from the electronic medical records and kidney biopsy database. No diagnosis of fatty liver was diagnosed by liver biopsy or ultrasonography.
Kidney specimens were obtained with a 16-gauge biopsy needle (Bard, New Providence, NJ). Specimens were fixed in 10% formalin, and the prepared sections were stained with hematoxylin–eosin, Masson trichrome, periodic acid silver methenamine, or periodic acid-Schiff. At least three pathologists evaluated each specimen.
The degree of liver fibrosis was evaluated by the FIB-4 index and NFS. FIB-4 was calculated using the following formula: age (years) × aspartate aminotransferase (AST, IU/L) / platelet (109/L) × √ alanine aminotransferase (ALT, IU/L). High FIB-4 index group was defined as 1.3 or higher30. NFS was calculated using the following formula: − 1.675 + 0.037 × age (years) + 0.094 × body mass index (kg/m2) + 1.13 × impaired fasting glucose/diabetes (yes = 1, no = 0) + 0.99 × AST (IU/L) / ALT (IU/L)—0.013 × platelet (109/L)—0.66 × albumin (g/dL).
Levels of serum creatinine were measured in all samples using an enzymatic method in laboratories, and the values are represented using two decimal places. The eGFR of each patient was calculated using the following formula: 194 × serum creatinine − 1.094 × age − 0.287 × 0.739 (if female). The rate of change in eGFR was used to evaluate renal function in each period and was calculated using the following formula: ((eGFR 12 months after kidney biopsy—eGFR before kidney biopsy))/eGFR before kidney biopsy) × 100.
Continuous variables are presented as medians. Comparisons were made between groups using the Mann–Whitney U test; categorical variables were presented as numbers (percentage) and compared using Pearson’s Chi-squared test and Fisher’s exact test, as appropriate. All analyses were performed using StatView (SAS Institute, Cary, CA, USA) and Excel software. Statistical significance was defined as P < 0.05.
The data that support the findings of this study are available in the figshare at https://figshare.com/articles/dataset/raw_data_update_xlsx/19590235, but restrictions apply to the availability of these data, and so are not publicity available. Data are however available from the author upon reasonable request and with permission of figshare.
Diabetes, C. et al. The effect of intensive treatment of diabetes on the development and progression of long-term complications in insulin-dependent diabetes mellitus. N. Engl. J. Med. 329, 977–986 (1993).
Article Google Scholar
Kramer, H. J. et al. Increasing body mass index and obesity in the incident ESRD population. J. Am. Soc. Nephrol. 17, 1453–1459 (2006).
Article Google Scholar
Mima, A. et al. Glomerular-specific protein kinase C-beta-induced insulin receptor substrate-1 dysfunction and insulin resistance in rat models of diabetes and obesity. Kidney Int. 79, 883–896 (2011).
Article CAS Google Scholar
Fuster, D. & Samet, J. H. Alcohol use in patients with chronic liver disease. N. Engl. J. Med. 379, 2579 (2018).
Article Google Scholar
Carmiel-Haggai, M., Cederbaum, A. I. & Nieto, N. A high-fat diet leads to the progression of non-alcoholic fatty liver disease in obese rats. FASEB J. 19, 136–138 (2005).
Article CAS Google Scholar
Sorino, P. et al. Development and validation of a neural network for NAFLD diagnosis. Sci. Rep. 11, 20240 (2021).
Article ADS CAS Google Scholar
Cotter, T. G. & Rinella, M. Nonalcoholic fatty liver disease 2020: The state of the disease. Gastroenterology 158, 1851–1864 (2020).
Article CAS Google Scholar
Chon, Y. E. et al. Decrease in waist-to-hip ratio reduced the development of chronic kidney disease in non-obese non-alcoholic fatty liver disease. Sci. Rep. 10, 8996 (2020).
Article ADS CAS Google Scholar
Castera, L., Friedrich-Rust, M. & Loomba, R. Noninvasive assessment of liver disease in patients with nonalcoholic fatty liver disease. Gastroenterology 156, 1264–1281 (2019).
Article Google Scholar
Hara, M. et al. Non-invasive fibrosis assessments of non-alcoholic fatty liver disease associated with low estimated glomerular filtration rate among CKD patients: The Fukuoka Kidney disease Registry Study. Clin. Exp. Nephrol. 25, 822–834 (2021).
Article CAS Google Scholar
Saito, H. et al. High FIB4 index is an independent risk factor of diabetic kidney disease in type 2 diabetes. Sci. Rep. 11, 11753 (2021).
Article ADS CAS Google Scholar
Vieira Barbosa, J. et al. Fibrosis-4 index can independently predict major adverse cardiovascular events in nonalcoholic fatty liver disease. Am. J. Gastroenterol. 117, 453–461 (2022).
Article Google Scholar
Fraser, D. A. et al. A structurally engineered fatty acid, icosabutate, suppresses liver inflammation and fibrosis in NASH. J. Hepatol. 76, 800–811 (2021).
Article Google Scholar
Mima, A. Inflammation and oxidative stress in diabetic nephropathy: New insights on its inhibition as new therapeutic targets. J. Diabetes Res. 2013, 248563 (2013).
Article Google Scholar
Mima, A., Yasuzawa, T., King, G. L. & Ueshima, S. Obesity-associated glomerular inflammation increases albuminuria without renal histological changes. FEBS Open Bio 8, 664–670 (2018).
Article CAS Google Scholar
Mima, A. Diabetic nephropathy: Protective factors and a new therapeutic paradigm. J. Diabetes Compl. 27, 526–530 (2013).
Article Google Scholar
Mima, A. et al. Protective effects of GLP-1 on glomerular endothelium and its inhibition by PKCbeta activation in diabetes. Diabetes 61, 2967–2979 (2012).
Article CAS Google Scholar
Dawood, A. A., Saleh, A. A., Elbahr, O., Gohar, S. F. & Habieb, M. S. Inverse relationship between the level of miRNA 148a–3p and both TGF-beta1 and FIB-4 in hepatocellular carcinoma. Biochem. Biophys. Rep. 27, 101082 (2021).
CAS PubMed PubMed Central Google Scholar
Mima, A. et al. Activation of Src mediates PDGF-induced Smad1 phosphorylation and contributes to the progression of glomerulosclerosis in glomerulonephritis. PLoS ONE 6, e17929 (2011).
Article ADS CAS Google Scholar
Fogo, A. B. Mesangial matrix modulation and glomerulosclerosis. Exp. Nephrol. 7, 147–159 (1999).
Article CAS Google Scholar
Kiapidou, S., Liava, C., Kalogirou, M., Akriviadis, E. & Sinakos, E. Chronic kidney disease in patients with non-alcoholic fatty liver disease: What the Hepatologist should know?. Ann. Hepatol. 19, 134–144 (2020).
Article Google Scholar
Mima, A., Qi, W. & King, G. L. Implications of treatment that target protective mechanisms against diabetic nephropathy. Semin. Nephrol. 32, 471–478 (2012).
Article CAS Google Scholar
Adams, L. A., Angulo, P. & Lindor, K. D. Nonalcoholic fatty liver disease. CMAJ 172, 899–905 (2005).
Article Google Scholar
Day, C. P. Non-alcoholic fatty liver disease: Current concepts and management strategies. Clin. Med. (Lond) 6, 19–25 (2006).
Article Google Scholar
Tolman, K. G., Fonseca, V., Dalpiaz, A. & Tan, M. H. Spectrum of liver disease in type 2 diabetes and management of patients with diabetes and liver disease. Diabetes Care 30, 734–743 (2007).
Article CAS Google Scholar
Han, E. & Lee, Y. H. Non-alcoholic fatty liver disease: The emerging burden in cardiometabolic and renal diseases. Diabetes Metab. J. 41, 430–437 (2017).
Article CAS Google Scholar
Lambrecht, J. et al. A PDGFRbeta-based score predicts significant liver fibrosis in patients with chronic alcohol abuse, NAFLD and viral liver disease. EBioMedicine 43, 501–512 (2019).
Article Google Scholar
Mundel, P. & Reiser, J. Proteinuria: An enzymatic disease of the podocyte?. Kidney Int. 77, 571–580 (2010).
Article CAS Google Scholar
Mima, A. et al. Glomerular VEGF resistance induced by PKCdelta/SHP-1 activation and contribution to diabetic nephropathy. FASEB J. 26, 2963–2974 (2012).
Article CAS Google Scholar
Kim, M. C. et al. Serum zinc level and hepatic fibrosis in patients with nonalcoholic fatty liver disease. PLoS ONE 15, e0240195 (2020).
Article CAS Google Scholar
Download references
The authors thank Dr. Tatsuhiko Mori (Osaka Medical and Pharmaceutical University) for the discussions on the statistical analysis and renal pathology. The authors also thank to Ms. Yumiko Ohgami for preparation of the renal pathology data. This work is supported by grants to A.M. from JSPS KAKENHI Grant Number 22K08368, the Osaka Kidney Foundation Grant Number OKF21-0003 and Takeda Medical Foundation.
Department of Nephrology, Osaka Medical and Pharmaceutical University, Osaka, 569-8686, Japan
Akira Mima
You can also search for this author in PubMed Google Scholar
A.M. contributed to the conception and design of the study. Material preparation, data collection, and analyses were performed by A.M. The first draft of the manuscript was written by A.M.
Correspondence to Akira Mima.
A.M. received a speaker’s honorarium from Otsuka, Kyowa Kirin, Mitsubishi Tanabe, Torii, Kowa, Sumitomo Pharma, Bayer, Eli Lilly, Astellas, and Boehringer Ingelheim. A.M. received research grants from Kyowa Kirin, Sumitomo Pharma, Otsuka, Torii, Daiichi-Sankyo, Mitsubishi Tanabe, Boehringer Ingelheim, and Eli Lilly.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Reprints and Permissions
Mima, A. Prediction of decreased estimated glomerular filtration rate using liver fibrosis markers: a renal biopsy-based study. Sci Rep 12, 17630 (2022). https://doi.org/10.1038/s41598-022-22636-9
Download citation
Received: 12 April 2022
Accepted: 18 October 2022
Published: 21 October 2022
DOI: https://doi.org/10.1038/s41598-022-22636-9
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.
Provided by the Springer Nature SharedIt content-sharing initiative
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
Advertisement
© 2022 Springer Nature Limited
Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.