- State Key Laboratory of Ophthalmology, Zhongshan Ophthalmic Center, Sun Yat-sen University, Guangdong Provincial Key Laboratory of Ophthalmology and Visual Science, Guangdong Provincial Clinical Research Center for Ocular Diseases, Guangzhou, China
Purpose: To develop a computational method for oxygen-saturation-related functional parameter analysis of retinal vessels based on traditional color fundus photography, and to explore their characteristic alterations in type 2 diabetes mellitus (DM).
Methods: 50 type 2 DM patients with no-clinically detectable retinopathy (NDR) and 50 healthy subjects were enrolled in the study. An optical density ratio (ODR) extraction algorithm based on the separation of oxygen-sensitive and oxygen-insensitive channels in color fundus photography was proposed. With precise vascular network segmentation and arteriovenous labeling, ODRs were acquired from different vascular subgroups, and the global ODR variability (ODRv) was calculated. Student’s t-test was used to analyze the differences of the functional parameters between groups, and regression analysis and receiver operating characteristic (ROC) curves were used to explore the discrimination efficiency of DM patients from healthy subjects based on these functional parameters.
Results: There was no significant difference in the baseline characteristics between the NDR and healthy normal groups. The ODRs of all vascular subgroups except the micro venule were significantly higher (
Conclusion: A computational method to extract the retinal vascular oxygen-saturation-related optical density ratios (ODRs) with single color fundus photography was developed, and increased ODRs and decreased ODRv of retinal vessels could be new potential image biomarkers of DM.
1 Introduction
The eye is a structurally advanced optical organ with a complex vascular network. The maintenance of retinal function is critically dependent on the normal functioning of the blood-retinal barrier (Kaur et al., 2008; O’Leary and Campbell, 2023). Fundus diseases and related systematic diseases can cause functional changes in the fundus such as ischemia and hypoxia, which can damage the blood-retinal barrier, leading to retinal vascular permeability changes, hemorrhage, exudation, neovascularization, and other lesions (Pournaras et al., 2008; Ascunce et al., 2023; Monickaraj et al., 2023). Diabetes mellitus (DM) is a common endocrine metabolic disease with its microvascular complications in the eye that can result in diabetic retinopathy (DR), which is the leading cause of blindness in working-age people (Bourne et al., 2021; Steinmetz et al., 2021). Hypoxia is considered to be the core pathogenesis of DR besides of hyperglycemia in retinal blood vessels caused by DM (Hardarson, 2013). Hypoxia along with hyperglycemia stimulate vascular endothelium and eventually lead to vascular endothelial cell dysfunction, which in turn leads to the occurrence of DR. (Hardarson and Stefánsson, 2012; Khoobehi et al., 2013). Therefore, it is equally important to detect retinal oxygen saturation as well as blood glucose in DM patients, which can help to prompt the potential occurrence of DR and avoid further visual impairment.
Retinal vessel is the only vascular system in the body that can be observed in vivo (Hanssen et al., 2022), and non-invasive imaging technique is the first choice for retinal vascular oxygen saturation measurement. Nevertheless, there are only limited imaging instruments available to measure retinal oxygen saturation, in which dual-wavelength retinal oximetry is a promising technique. By measuring the light of different wavelengths reflected from the eye, it takes advantage of the light absorption variations between oxyhemoglobin and hemoglobin, and a calculation can be made that correlates directly to the retinal vascular oxygenation (Beach, 2014; Jeppesen and Bek, 2019; Garg et al., 2021). Previous studies using retinal oximetry have found that the oxygen saturation of retinal vessels in DR is significantly higher than that in normal subjects, and increases with the severity of DR (Hardarson and Stefánsson, 2012; Khoobehi et al., 2013; Jørgensen and Bek, 2014). Retinal oximetry has also been used to monitor retinal oxygen saturation in other eye diseases such as central retinal vein occlusion (CRVO) and glaucoma, showing varying degrees of increase in retinal oxygen saturation (Stefánsson et al., 2019). However, despite its potential, retinal oximetry is not widely used in ophthalmology clinics due to its recent commercialization.
Traditional color fundus cameras are the most commonly used equipment for fundus examination and disease screening (Karlsson et al., 2021). Their broad-spectrum cameras use complementary metal-oxide semiconductor (CMOS) or charge couple device (CCD) photoreceptors covered with red-green-blue (RGB) filters (Malvar et al., 2004). Thus, oxygen-sensitive and insensitive wavelengths reflected from retina vessels, which corresponds to different oxyhemoglobin and hemoglobin light absorptions, are generally received in the red and green channels, respectively, allowing for the estimation of blood oxygen function by extracting light absorption information from these channels.
In this study, we aim to develop a computational method for analyzing oxygen-saturation-related functional parameters with traditional color fundus photography based on precise vessel segmentation and light absorption analysis, providing an easy-to-apply image analysis algorithm for traditional color fundus cameras. With that, we also explore the characteristic alterations of retinal vascular oxygen function in type 2 DM patients to look for new potential image biomarkers for the monitoring of DM.
2 Methods
2.1 Ethics
This was a retrospective observational study. The project underwent a formal administrative review and was determined to be methodologically validated according to the institutional policy of Zhongshan Ophthalmic Center, Sun Yat-sen University. Therefore, this study was considered not to be a human subject’s study, but the study was still reviewed by Institutional Review Board to avoid information leakage during data processing (protocol number: 2017KYPJ104).
2.2 Population
We included 50 patients with type 2 diabetes mellitus (DM) and a normal group of 50 age- and sex-matched healthy normal subjects between 1 January 2019 and 15 May 2021. All the subjects were adults. We extracted the age and sex information of the subjects from the medical records for intergroup matching and determined whether the subjects are patients with type 2 DM or not. To assess the existence of diabetic retinopathy, standard seven-fields fundus photographs were obtained with a mydriasis-free digital fundus camera (RetiCam 3100, SYSEYE, China). Patients with clinically detectable retinopathy, or with a history of ocular disease, inflammation, trauma, or any intraocular surgery were excluded. Further exclusion criteria include related systemic diseases, such as Alzheimer’s disease, hypertension or lung disease. The healthy normal group should also meet the above exclusion criteria in addition to no diabetes.
2.3 Vessel segmentation and labeling
Only the 50° fundus photography centered on the macula are analyzed. Figure 1 shows the flowchart for the processing of the region of interest (ROI) of retinal vessels. To extract the retinal vascular network, we employed an intelligent automatic vascular segmentation method, details of which are available in our earlier publication (Wang et al., 2021; 2023). Briefly, the original fundus color photographs (Figure 1A) were processed with an enhanced U-net by employing a multipath attention network model (MA-net) for improved vessel segmentation, achieving an area under the curve (AUC) of 0.9838 for segment accuracy (Wang et al., 2023). Blood vessel junctions on the obtained retinal vascular binarization network (Figure 1B) were detected then through a convolution method on vessel skeleton, facilitating the vessel segments extraction (Figure 1C). The arteriovenous labeling of the segmented vessels was further processed based on an intra-image regularization classifier (Figure 1D) (Xu et al., 2017).
FIGURE 1. Region of interest extraction process. (A) Original color fundus photography. (B) Retinal vascular binarization network. (C) Retinal vessel segments map. (D) Arterial and venous labeling. (E,F) Green and red image channels extracted from the original fundus image. (G,H) Green and red channel of the vessel segments map. (I,J) Intravascular region of interest (IROI) and extravascular region of interest (EROI).
To ensure the calculation of the oxygen-saturated-related functional parameter of retinal vessels, the red and green channels (Figures 1E,F) composing the original fundus color photography were extracted through channel separation. By multiplying them with the binarized vessel segments map, oxygen-sensitive and insensitive wavelengths reflected from the retinal vascular areas were obtained (Figures 1G,H), of which each vessel segment was marked as the intravascular regions of interest (IROI, Figure 1I). To obtain the light reflection information of the retinal background around the vessel segment, a range of one vessel diameter is extended outward along the blood vessel segment outlines, which was marked as the extravascular region of interest (EROI, Figure 1J).
2.4 Optical density ratio extraction
According to the EROI and IROI of blood vessels, the optical density (OD) values were calculated for every three pixels along the vascular central line (Formula 1), the red channel is used as the oxygen sensitive channel, and the green channel represents the oxygen insensitive channel, and the optical density ratio (ODR) of the two is defined as the blood oxygen-saturation-related parameter (Formula 2) (Tiedeman et al., 1998; Beach et al., 1999).
In order to study the differences of different types of blood vessels, all branch vessels output ODR according to the labeling of arteries and veins, and distinguish between main and micro vessels referring to the setting of commercial retinal oximeter (Stefánsson et al., 2017), i.e., ≥6pixel diameter (1pixel≈12.69 μm in this study) was defined as main vessels and <6pixel diameter was defined as micro vessels. Finally, ODR outputs the results according to all vessels (All), all arteriole (A), main arteriole (main A), micro arteriole (mic A), all venule (V), main venule (main V), and micro venule (mic V).
In addition to the ODR of each vascular subgroup, we also average the absolute value of the difference of ODR between all adjacent blocks of pixels and define it as ODR variability (ODRv, Formula 3).
2.5 Statistical analysis
Student’s t-test or chi-square test used to compare the differences between Normal and DM groups. All data were recorded as mean ± standard deviation (SD). Data distribution was analyzed by the Kolmogorov–Smirnov test to determine normally distributed data (
3 Results
There is no difference in the age and sex between the 50 healthy volunteers in Normal group and the 50 patients in DM group (Table 1).
Table 1 alsodemonstrates the ODR for all vascular subgroups, ODRv, as well as the t-statistic and p-values. In both healthy normal group and DM group, A ODR is higher than V ODR (mean ± SD, A vs. V, Normal: 0.713 ± 0.057 vs. 0.708 ± 0.063; DM: 0.755 ± 0.083 vs. 0.754 ± 0.082), main A ODR is higher than main V ODR (Normal: 0.728 ± 0.056 vs. 0.717 ± 0.058; DM: 0.773 ± 0.083 vs. 0.761 ± 0.082) and micro A is higher than micro V ODR (Normal: 0.634 ± 0.068 vs. 0.599 ± 0.082; DM: 0.672 ± 0.096 vs. 0.627 ± 0.117), which means that ODR is positively correlated with blood oxygen saturation. Compared to Normal group, ODR is significantly increased in all vascular subgroups except micro venular in DM group (
As shown in the logistics regression results in Table 2, higher ODRs are significantly associated with the occurrence of DM, especially in main A (OR 1.953, 95% CI 1.252-3.046) and main V (OR 1.908, 95% CI 1.227-2.968), and also in micro A (OR 1.604, 1.047-2.455). While Lower ODRv (OR 0.381, 95% CI 0.210-0.693) is significantly associated with the occurrence of DM. The C-statistic for the discrimination of DM from the healthy normal group combining all the acquired oxygen-saturation-related functional parameters is 0.777 (95% CI 0.687-0.867,
4 Discussion
In this study, we present a computational approach for extracting oxygen-saturation-related functional parameters through precise vascular segmentation and light reflection channel analysis using single traditional color fundus photography. By applying this method to the comprehensive analysis of the fundus images acquired from normal individuals and patients with type 2 diabetes, we revealed that the retinal vascular ODRs were increased while the ODR variability was reduced significantly in DM patients, providing new potential image biomarkers for the management of DM.
The successful extraction of retinal vascular ODRs with color fundus photography offers a new protocol for retinal vascular function evaluation, providing additional quantitative functional information based on the most commonly used fundus imaging tools in clinical settings. Compared to the current existed retinal oximetry, which only reveals oxygen saturation level in main retinal vessels and requires subjects to withstand long light exposure time (Hardarson, 2013), our method achieves information in microvascular level without any additional hardware implementation, which ensures its convenience and accessibility with greatly reduced cost. These advantages suggest extensive potential application value in the diagnosis and screening of retinal diseases.
In our study, we have shown that the increase of proposed oxygen-saturation-related ODRs were statistically significant at the stage of DM, which could be related to the stronger affinity of glycosylated hemoglobin to oxygen in DM (Graham et al., 1980; van Kampen and Zijlstra, 1983). Indeed, former studies using retinal oximetry have found obvious increase in retinal oxygen saturation level during DR, especially in proliferative DR, but no significant change was shown in DM stage before (Khoobehi et al., 2013). This indicates that our newly proposed ODRs might be more sensitive parameters than that of commercial retinal oximetry in finding alterations in blood oxygen function. Moreover, our exploration of ODR in the microvascular groups revealing significant increase in micro arterioles of DM patients further proves the high sensitivity of our method.
We have proposed and analyzed the ODR variability for the first time on the basis of extracted ODR of the full retinal vascular network and found it to be significantly decreased in DM patients compared to normal subjects. While glycosylated hemoglobin has higher affinity to oxygen (Graham et al., 1980), it damps the oxygen-release capacity (Ditzel, 1972) of retinal blood, results in reduced oxygen saturation variations. Studies have also found that the retinal blood flow velocity increases in DM patients (Landa et al., 2011), meaning less oxygen release time of blood in retinal vessels, further reduced the ODR variability.
Our method suffers limited retinal vessel resolution and imaging depth due to the physical nature of fundus photography. Only the superficial retinal vessels can be analyzed without 3-dimensional (3D) illustration. The same limitations also exist in current commercial retinal oximetry. Recently, visible-light optical coherence tomography based techniques for retinal oxygen saturation measurement have been gradually proposed (Pi et al., 2020; Poddar and Basu, 2020), achieving high resolution 3D imaging of retinal vascular oxygen function information (Pi et al., 2020), although most of these studies are still in laboratorial stage.
Our attempt to distinguish DM from normal subjects with all the extracted functional ODR parameters showed a relative low C statistic result (AUC = 0.777), which is not sufficient for DM diagnosis or screening. Since our previous study found that the morphological parameters of retinal vessels in DM patients also changed significantly (Li et al., 2021), further combining functional and morphological parameters of retinal blood vessels and multi-center studies or follow-up studies may help us better understand the changes in retinal blood vessel characteristics during the occurrence and development of DM and even DR, so as to diagnose DM or DR More accurately. Furthermore, applying our method to other eye diseases related to retinal vascular oxygen function alterations such as glaucoma (Vandewalle et al., 2014), central retinal vein occlusion (Hardarson and Stefánsson, 2010), age-related macular degeneration (Geirsdottir et al., 2014), etc. in the future would further demonstrate the reliability and universality of our method.
Data availability statement
The raw data supporting the conclusions of this article are available from the corresponding authors upon reasonable request.
Ethics statement
The studies involving human participants were reviewed and approved by Institutional Review Board of Zhongshan Ophthalmic Center, Sun Yat-sen University. The patients/participants provided their written informed consent to participate in this study.
Author contributions
JZ and ZL: Conceptualization, methodology, software, and writing—original draft. GW and YH: methodology and writing—review and editing. KF, YL, and JL: Investigation and data collection and writing—review and editing. JY and PX: Conceptualization, writing—review and editing, and supervision. All authors read and approved the final manuscript. All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.
Funding
This work was supported by the National Natural Science Foundation of China (Nos. 82230033 and 82271133), the Department of Science and Technology of Guangdong Province (Nos. 2021TX06Y127 and 2021TQ06Y137), and the Basic and Applied Basic Research Foundation of Guangdong Province (No. 2022A1515011486).
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
Ascunce, K., Dhodapkar, R. M., Huang, D., and Hafler, B. P. (2023). Innate immune biology in age-related macular degeneration. Front. Cell Dev. Biol. 11, 1118524. doi:10.3389/fcell.2023.1118524
Beach, J. M., Schwenzer, K. J., Srinivas, S., Kim, D., and Tiedeman, J. S. (1999). Oximetry of retinal vessels by dual-wavelength imaging: Calibration and influence of pigmentation. J. Appl. Physiol. 86, 748–758. doi:10.1152/jappl.1999.86.2.748
Bourne, R., Steinmetz, J. D., Flaxman, S., Briant, P. S., Taylor, H. R., Resnikoff, S., et al. (2021). Trends in prevalence of blindness and distance and near vision impairment over 30 years: An analysis for the global burden of disease study. Lancet Glob. Health 9, e130–e143. doi:10.1016/S2214-109X(20)30425-3
Ditzel, J. (1972). Impaired oxygen release caused by alterations of the metabolism in the erythrocytes in diabetes. Lancet 1, 721–723. doi:10.1016/s0140-6736(72)90233-4
Garg, A. K., Knight, D., Lando, L., and Chao, D. L. (2021). Advances in retinal oximetry. Transl. Vis. Sci. Technol. 10, 5. doi:10.1167/tvst.10.2.5
Geirsdottir, A., Hardarson, S. H., Olafsdottir, O. B., and Stefánsson, E. (2014). Retinal oxygen metabolism in exudative age-related macular degeneration. Acta Ophthalmol. 92, 27–33. doi:10.1111/aos.12294
Graham, J. J., Ryall, R. G., and Wise, P. H. (1980). Glycosylated haemoglobin and relative polycythaemia in diabetes mellitus. Diabetologia 18, 205–207. doi:10.1007/BF00251917
Hanssen, H., Streese, L., and Vilser, W. (2022). Retinal vessel diameters and function in cardiovascular risk and disease. Prog. Retin Eye Res. 91, 101095. doi:10.1016/j.preteyeres.2022.101095
Hardarson, S. H., and Stefánsson, E. (2010). Oxygen saturation in central retinal vein occlusion. Am. J. Ophthalmol. 150, 871–875. doi:10.1016/j.ajo.2010.06.020
Hardarson, S. H., and Stefánsson, E. (2012). Retinal oxygen saturation is altered in diabetic retinopathy. Br. J. Ophthalmol. 96, 560–563. doi:10.1136/bjophthalmol-2011-300640
Jeppesen, S. K., and Bek, T. (2019). The retinal oxygen saturation measured by dual wavelength oximetry in larger retinal vessels is influenced by the linear velocity of the blood. Curr. Eye Res. 44, 46–52. doi:10.1080/02713683.2018.1524015
Jørgensen, C., and Bek, T. (2014). Increasing oxygen saturation in larger retinal vessels after photocoagulation for diabetic retinopathy. Invest. Ophthalmol. Vis. Sci. 55, 5365–5369. doi:10.1167/iovs.14-14811
Karlsson, R. A., Jonsson, B. A., Hardarson, S. H., Olafsdottir, O. B., Halldorsson, G. H., and Stefansson, E. (2021). Automatic fundus image quality assessment on a continuous scale. Comput. Biol. Med. 129, 104114. doi:10.1016/j.compbiomed.2020.104114
Kaur, C., Foulds, W. S., and Ling, E. A. (2008). Blood-retinal barrier in hypoxic ischaemic conditions: Basic concepts, clinical features and management. Prog. Retin Eye Res. 27, 622–647. doi:10.1016/j.preteyeres.2008.09.003
Khoobehi, B., Firn, K., Thompson, H., Reinoso, M., and Beach, J. (2013). Retinal arterial and venous oxygen saturation is altered in diabetic patients. Invest. Ophthalmol. Vis. Sci. 54, 7103–7106. doi:10.1167/iovs.13-12723
Landa, G., Amde, W., Haileselassie, Y., and Rosen, R. B. (2011). Cilioretinal arteries in diabetic eyes are associated with increased retinal blood flow velocity and occurrence of diabetic macular edema. Retina 31, 304–311. doi:10.1097/IAE.0b013e3181e91108
Li, M., Wang, G., Xia, H., Feng, Z., Xiao, P., and Yuan, J. (2021). Retinal vascular geometry detection as a biomarker in diabetes mellitus. Eur. J. Ophthalmol. 112067212110334, 1710–1719. doi:10.1177/11206721211033488
Malvar, H. S., He, L.-W., and Cutler, R. (2004). “High-quality linear interpolation for demosaicing of Bayer-patterned color images,” in 2004 IEEE International Conference on Acoustics, Speech, and Signal Processing, Montreal, Canada, May 2004 (IEEE), 485–488.
Monickaraj, F., Acosta, G., Cabrera, A. P., and Das, A. (2023). Transcriptomic profiling reveals chemokine CXCL1 as a mediator for neutrophil recruitment associated with blood-retinal barrier alteration in diabetic retinopathy. Diabetes, db220619. doi:10.2337/db22-0619
O’Leary, F., and Campbell, M. (2023). The blood-retina barrier in health and disease. FEBS J. 290, 878–891. doi:10.1111/febs.16330
Pi, S., Hormel, T. T., Wei, X., Cepurna, W., Wang, B., Morrison, J. C., et al. (2020). Retinal capillary oximetry with visible light optical coherence tomography. Proc. Natl. Acad. Sci. U.S.A. 117, 11658–11666. doi:10.1073/pnas.1918546117
Poddar, R., and Basu, M. (2020). Characterization and oxygen saturation study of human retinal blood vessels evaluated by spectroscopic Optical Coherence Tomography Angiography. Opt. Laser Technol. 122, 105886. doi:10.1016/j.optlastec.2019.105886
Pournaras, C. J., Rungger-Brändle, E., Riva, C. E., Hardarson, S. H., and Stefansson, E. (2008). Regulation of retinal blood flow in health and disease. Prog. Retin Eye Res. 27, 284–330. doi:10.1016/j.preteyeres.2008.02.002
Stefánsson, E., Olafsdottir, O. B., Einarsdottir, A. B., Eliasdottir, T. S., Eysteinsson, T., Vehmeijer, W., et al. (2017). Retinal oximetry discovers novel biomarkers in retinal and brain diseases. Invest. Ophthalmol. Vis. Sci. 58, BIO227–BIO233. doi:10.1167/iovs.17-21776
Stefánsson, E., Olafsdottir, O. B., Eliasdottir, T. S., Vehmeijer, W., Einarsdottir, A. B., Bek, T., et al. (2019). Retinal oximetry: Metabolic imaging for diseases of the retina and brain. Prog. Retin. Eye Res. 70, 1–22. doi:10.1016/j.preteyeres.2019.04.001
Steinmetz, J. D., Bourne, R. R. A., Briant, P. S., Flaxman, S. R., Taylor, H. R. B., Jonas, J. B., et al. (2021). Causes of blindness and vision impairment in 2020 and trends over 30 years, and prevalence of avoidable blindness in relation to VISION 2020: The right to sight: An analysis for the global burden of disease study. Lancet Glob. Health 9, e144–e160. doi:10.1016/S2214-109X(20)30489-7
Tiedeman, J. S., Kirk, S. E., Srinivas, S., and Beach, J. M. (1998). Retinal oxygen consumption during hyperglycemia in patients with diabetes without retinopathy. Ophthalmology 105, 31–36. doi:10.1016/s0161-6420(98)71029-1
Van Kampen, E. J., and Zijlstra, W. G. (1983). Spectrophotometry of hemoglobin and hemoglobin derivatives. Adv. Clin. Chem. 23, 199–257. doi:10.1016/s0065-2423(08)60401-1
Vandewalle, E., Abegão Pinto, L., Olafsdottir, O. B., De Clerck, E., Stalmans, P., Van Calster, J., et al. (2014). Oximetry in glaucoma: Correlation of metabolic change with structural and functional damage. Acta Ophthalmol. 92, 105–110. doi:10.1111/aos.12011
Wang, G., Huang, Y., Ma, K., Duan, Z., Luo, Z., Xiao, P., et al. (2023). Automatic vessel crossing and bifurcation detection based on multi-attention network vessel segmentation and directed graph search. Comput. Biol. Med. 155, 106647. doi:10.1016/j.compbiomed.2023.106647
Wang, G., Li, M., Yun, Z., Duan, Z., Ma, K., Luo, Z., et al. (2021). A novel multiple subdivision-based algorithm for quantitative assessment of retinal vascular tortuosity. Exp. Biol. Med. (Maywood) 246, 2222–2229. doi:10.1177/15353702211032898
Keywords: fundus photography, diabetes mellitus, retina vessel segmentation, oxygen saturation level, light reflection analysis
Citation: Zhang J, Luo Z, Wang G, Huang Y, Fei K, Liu Y, Li J, Yuan J and Xiao P (2023) Oxygen-saturation-related functional parameter as a biomarker for diabetes mellitus—extraction method and clinical validation. Front. Cell Dev. Biol. 11:1195873. doi: 10.3389/fcell.2023.1195873
Received: 29 March 2023; Accepted: 02 May 2023;
Published: 11 May 2023.
Edited by:
Weihua Yang, Jinan University, ChinaReviewed by:
Bilian Ke, Shanghai General Hospital, ChinaGangjin Kang, The Affiliated Hospital of Southwest Medical University, China
Fan Li, Zhongshan People’s Hospital (ZSPH), China
Copyright © 2023 Zhang, Luo, Wang, Huang, Fei, Liu, Li, Yuan and Xiao. 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: Jin Yuan, eXVhbmppbmNvcm5lYUAxMjYuY29t; Peng Xiao, eGlhb3BlbmdhZGRpc0Bob3RtYWlsLmNvbQ==
†These authors have contributed equally to this work