- Department of Radiology, Huadong Hospital, Fudan University, Shanghai, China
Purpose: This study aims to explore the value of the delta-radiomics (DelRADx) model in predicting the invasiveness of lung adenocarcinoma manifesting as radiological part-solid nodules (PSNs).
Methods: A total of 299 PSNs histopathologically confirmed as lung adenocarcinoma (training set, n = 209; validation set, n = 90) in our hospital were retrospectively analyzed from January 2017 to December 2021. All patients underwent diagnostic noncontrast-enhanced CT (NCECT) and contrast-enhanced CT (CECT) before surgery. After image preprocessing and ROI segmentation, 740 radiomic features were extracted from NCECT and CECT, respectively, resulting in 740 DelRADx. A DelRADx model was constructed using the least absolute shrinkage and selection operator logistic (LASSO-logistic) algorithm based on the training cohort. The conventional radiomics model based on NCECT was also constructed following the same process for comparison purposes. The prediction performance was assessed using area under the ROC curve (AUC). To provide an easy-to-use tool, a radiomics-based integrated nomogram was constructed and evaluated by integrated discrimination increment (IDI), calibration curves, decision curve analysis (DCA), and clinical impact plot.
Results: The DelRADx signature, which consisted of nine robust selected features, showed significant differences between the AIS/MIA group and IAC group (p < 0.05) in both training and validation sets. The DelRADx signature showed a significantly higher AUC (0.902) compared to the conventional radiomics model based on NCECT (AUC = 0.856) in the validation set. The IDI was significant at 0.0769 for the integrated nomogram compared with the DelRADx signature. The calibration curve of the integrated nomogram demonstrated favorable agreement both in the training set and validation set with a mean absolute error of 0.001 and 0.019, respectively. Decision curve analysis and clinical impact plot indicated that if the threshold probability was within 90%, the integrated nomogram showed a high clinical application value.
Conclusion: The DelRADx method has the potential to assist doctors in predicting the invasiveness for patients with PSNs. The integrated nomogram incorporating the DelRADx signature with the radiographic features could facilitate the performance and serve as an alternative way for determining management.
Introduction
With the popularity of low-dose CT for lung cancer screening in recent years, a large number of early-stage lung adenocarcinomas radiographically manifesting as partial solid nodules (PSNs) have been screened out (1). Given the distinct surgical management strategy and prognosis for disease-free survival (DFS) and overall survival (OS), differentiation of the invasive adenocarcinoma (IAC) from the adenocarcinoma in situ (AIS) or minimally invasive adenocarcinoma (MIA) for this distinct subtype has been confirmed to be of great clinical significance (2). The malignant potential of PSNs in CT images has not been clarified; however, the evaluation of invasive qualities was primarily imaging- or clinic-based at the present stage, which is nonquantifiable and subjective (3, 4). Both radiologists and surgeons desire a precise and practical way to help guide clinical decision-making for the invasiveness of PSNs when faced with a such diagnostic quandary.
Histological evidence has revealed a tight relationship between the solid component of PSNs and the invasive component of adenocarcinomas (5, 6). The extra blood supply information of solid components may benefit from contrast-enhanced CT (CECT), which also brings out tumor heterogeneity (7, 8). Furthermore, a new high-throughput radiomics analysis method called delta-radiomics (DelRADx), which deals with a string of eigenvalue changes within different modes, has reportedly been linked to the efficacy or prognosis for malignancies such as colorectal, liver, or lung cancer (9–11). However, the use of DelRADx for the prediction of invasive adenocarcinomas manifested as PSNs is rarely reported. The purpose of this study is to develop a DelRADx model based on CECT and noncontrast-enhanced CT (NCECT) data in order to provide patients with PSNs for better decision support.
Materials and methods
Patients
Patients who underwent surgical resections for lung cancer at our hospital between January 2017 and December 2021 were reviewed retrospectively. The inclusion criteria were as follows: (1) pathologically confirmed as lung adenocarcinoma; (2) radiographically represented as PSNs (diameter, ≤3 cm) in axial CT at lung window setting; (3) NCECT and CECT images were obtained at one examination; and (4) thin-slice CT images (1–1.25 mm) could be obtained. The exclusion criteria were as follows: (1) nodule diameter less than 6 mm, due to the inherent calibration error of solid components within the threshold; (2) poor CT quality, such as severe motion artifacts; and (3) failure to extract the radiomic features for unknown reasons. For the multiple PSNs in one patient, only the lesion with the pathologically conclusive result was included.
Consequently, 299 cases with 299 PSNs who met the principles were enrolled in the cohort for the current analysis (detailed in Figure S1). According to the histopathological examination, 159 cases were diagnosed with AIS/MIA, while, 140 cases were diagnosed with IAC. The cases were randomly divided into training and validation sets in a ratio of 7:3. This study design was approved by the institutional research ethics board of our institution, and the informed consent requirement was waived for the retrospective research with anonymous data.
Image acquisition and radiographic feature assessment
NCECT examinations were conducted in the supine position with the arms up after deep inspiration. The CT data were acquired from one of the two scanners: Somatom Definition flash (Siemens Medical Solutions, Germany) and GE Discovery CT750 HD scanner (GE Healthcare, USA). The scanning parameters are detailed in Supplementary Table S1. Subsequently, the CECT was performed at 35 to 60s after injecting a dose of 80–100 ml nonionic IV contrast material (350 mg/ml, Omnipaque, GE Healthcare) mixed with isotonic saline into the ulnar vein using a high-pressure syringe at a rate of 3.0–4.0 ml/s.
The clinical and radiographic features were reviewed by two radiologists (H.W. and D.B.M.) with more than 10 years of experience in chest CT interpretation in a blinded fashion. We used the electronic caliper in our picture archiving and communication system to measure the maximum diameter of PSNs (nodule_max) and the maximum diameter of the corresponding solid component (solid_max). The consolidation-to-tumor ratio (CTR) was subsequently determined by dividing the solid_max by the nodule_max. Any discrepancies in describing the radiographic features were settled by consensus reading.
Nodule segmentation
The regions of interest (ROI) were manually contoured slicewise by one radiologist (Z.X.M.) with 5 years of experience to achieve three-dimensional segmentation using the open-source medical image processing and navigation software 3D slicer (version 4.8.0, Brigham and Women’s Hospital). Another radiologist (Y.L.Y.) with 6 years of experience segmented a random set of 20 nodules independently to assess the interobserver robustness of radiomic features. All ROI were exported in Nrrd (desensitization format) for the following analysis.
Image preprocessing and radiomics feature extraction
All images were isotropically resampled with 1.0 mm at X/Y/Z-spacing using linear interpolation to standardize. The Gaussian filter was used to reduce the noise influence of the voxel on radiomic features. Feature extraction was performed with pyRadiomics (https://doi.org/10.1158/0008-5472.CAN-17-0339). The DelRADx was defined as the change of radiomic features from NCECT to CECT: DelRADx = FeatureCECT − FeatureNCECT.
Feature selection and modeling
The training set was used for feature selection. Robust features with an interclass correlation coefficient (ICC) of >0.8 were chosen for further analysis. The discriminative ability was first evaluated using the Mann–Whitney U test. Features with statistical significance were taken for the following analysis. The correlation matrix with the pair-wise Spearman correlation analysis was built to eliminate the redundant features (correlation coefficient >0.90). The least absolute shrinkage and selection operator (LASSO) logistic regression algorithm with fivefold cross-validation was used to develop a predictive DelRADx signature by linear summing the core DelRADx multiplied by their coefficient. Another conventional radiomics model based on the NCECT features was constructed by using the same flow process to evaluate the optimal predictive model.
Statistical analysis
The feature selection, modeling, and statistical analysis were performed with R software (version 3.6.2; http://www.Rproject.org) or SPSS 21.0 (IBM, Chicago, IL, USA). The “irr” package was used for ICC. LASSO analysis was performed with the “glmnet” package. The nomogram and calibration curve were plotted based on the “rms” package. Decision curve analysis (DCA) and clinical impact plot were done with the “rmda” package. Multivariate binary logistic regression analysis was performed with an input parameter strategy. The model performance between the two radiomics models was evaluated by the ROC analysis. The significant difference was evaluated using the DeLong method. Integrated discrimination increment (IDI) was used to analyze the improvement of the integrated nomogram. A p-value of < 0.05 was considered statistically significant. All p-values were two sided in this study.
Results
Patients’ general characteristics
The general characteristics of 299 patients in the training and validation sets are summarized in Table 1. In the univariate analysis, there were no significant differences in gender, age, location, margin, air bronchogram, and vacuole sign (p >0.05). Statistically significant differences could be observed in the nodule_max, solid_max, and CTR both in the training and validation sets (p < 0.05). Consequently, these three quantitative parameters were chosen to establish a radiographic model. According to the univariate logistic regression analysis, only solid_max could independently predict the invasiveness of PSNs (Table S2).
The AUC of the radiographic model based on the binary logistic regression analysis was 0.890 (95% CI, 0.843–0.938) in the training set and 0.835 (95% CI, 0.745–0.926) in the validation set. The associated criterion of solid_max was 4.53 mm.
Feature selection and DelRADx signature building
A total of 740 radiomic features were extracted from the ROI of NCECT and CECT, respectively, resulting in 740 DelRADx. The schematic depiction of radiomics modeling was illustrated in Figure 1. After the process of ICC, ANOVA/MW, Spearman rank correlation test, and LASSO logistic regression analysis (detailed in Figure S2), nine robust DelRADx were ultimately selected. Based on the coefficients, the DelRADx signature was calculated for each patient. The signature formulas are provided in Figure S3.
Both the DelRADx and conventional radiomics signatures were significantly different between IAC and the AIS/MIA in the training set (0.78 ± 0.23 vs. 0.22 ± 0.25, 0.74 ± 0.24 vs. 0.26 ± 0.26) and the validation set (0.64 ± 1.25 vs. −1.99 ± 1.96, 0.44 ± 0.82 vs. −0.84 ± 0.81) (p-values < 0.05).
Performance comparison
The DelRADx model exhibited good performance with an AUC of 0.925 (95% CI, 0.888–0.962) in the training set and 0.902 (95% CI, 0.838–0.966) in the validation set, respectively. For the conventional radiomics models, the AUC was 0.894 (95% CI, 0.853~0.936) in the training set and 0.856 (95% CI, 0.767–0.945) in the validation set, respectively. Compared with the conventional radiomics model, the DelRADx model showed significantly higher AUC in the validation set (p < 0.05 of the DeLong test) (Figure 2).
Figure 2 The predictive performance of the DelRADx signature and conventional radiomics signature for each patient in the training set (A) and validation set (B).
Radiomics nomogram building and evaluation
To provide an easy-to-use tool, a delta-radiomics-based integrated nomogram incorporating the DelRADx signature and the solid_max was constructed using the multivariable logistical regression analysis (as shown in Figure 3). According to the univariate logistic regression analysis, both the DelRADx signature and the solid_max could independently predict the invasiveness of PSNs (Table 2). Compared with the DelRADx signature, the total IDI was significant at 0.0769 for the integrated nomogram (95% CI, 0.0394–0.1144, p < 0.001).
Table 2 Multivariate analysis of radiomics-based integrated nomogram for differentiating invasive adenocarcinoma.
The calibration curve of the integrated nomogram demonstrated favorable agreement with actual observation in the training cohort and was confirmed in the validation set with a mean absolute error of 0.001 and 0.019, respectively (Figure 4). Decision curve analysis and clinical impact plot indicated that if the threshold probability of a patient was within 90%, using the integrated nomogram to predict IAC added more benefit than either the treat-all-patient scheme or the treat-none scheme (Figures 5A, B).
Figure 4 The calibration curve of the integrated nomogram in the training dataset (A) and validation dataset (B).
Figure 5 Decision curve analysis (A) and clinical impact plot (B) for the integrated nomogram in the validation set. The x-axis indicates threshold probability, and the y-axis indicates the net benefit. The red line represents the integrated nomogram.
Discussion
In this study, we developed and validated a diagnostic DelRADx signature in the individualized evaluation of invasive adenocarcinomas for patients with PSNs. Our results demonstrated that the proposed DelRADx model outperformed the conventional radiomics model based on NCECT with an AUC of 0.925 in the training set and 0.902 in the validation set, respectively. Furthermore, our easy-to-use nomogram incorporating the DelRADx signature and radiographic independent predictor facilitated the performance compared with the DelRADx model alone. Our findings indicate that the DelRADx signature could provide additional information for invasive adenocarcinomas and, consequently, help to provide better support for decision-making when treating patients with PSNs.
The association between the solid component and pathological invasion has been extensively researched (12, 13). Several studies suggest good prediction of IACs if the threshold of the solid component is >5 mm (14). Our data are in keeping with previous observations that the diameter of solid components is a significant independent differentiator for IACs. While in our investigation, the prediction threshold was preshifted to 4.5 mm. The observed discrepancy could be explained in part by the fact that the current study included not only the subcohort of MIA but also a sample of AIS, which contained fewer solid components. Similarly, Luo et al. (15) developed a model with pleural indentation, solid component size, and solid component proportion for differentiating IAC from non-IAC in patients with PSNs and achieved an AUC of 0.85. Weng et al. (16) used lesion shape and solid component size to create a prediction model for PSNs with an AUC of 0.76.
Despite the predictive power of the radiographic features, the accurate identification of IACs with subjective assessment remains challenging in clinical practice. Radiomics is a high-throughput and automatic analysis technique for medical images. Numerous studies have proved that the development of radiomics represents a significant breakthrough in overcoming the limitation of subjective evaluation (17). Delta radiomics have shown potential advantages in earlier diagnosis and prognosis estimation in a variety of tumors.
In this study, we constructed an individualized delta-radiomics model to predict the invasiveness of PSNs. The results indicated that our delta-radiomics signature derived from CECT and NCECT outperformed the conventional radiomics signature. One of the potential explanations is that the images obtained from CECT could reveal the density of microvessels, which is closely linked with tumor invasion (18). Additionally, delta-radiomics could provide supplementary information on intratumoral heterogeneity, which is also supporting evidence for tumor invasion.
Previous studies support our interpretation. Son et al. (19) found that CECT imaging metrics could add value in distinguishing invasive adenocarcinoma from AIS/MIA that manifested as PSNs. For the DelRADx analysis, Wang et al. investigated a delta-radiomics model derived from the positioning CT and resetting CT after radiation therapy, and the result demonstrated that the delta-radiomics signature could be a promising image biomarker for the prediction of severe acute radiation pneumonitis. In line with our study, Saeed et al. (11) indicated that the Rider features in the machine learning-based delta-radiomics model could improve the performance of lung cancer screening. Our study found that the delta-radiomic features of the gray-level run-length matrix (GLRLM) had the added value in differentiating the invasiveness of lung adenocarcinoma. As is well-known, GLRLM is used to describe the distribution of texture changes between neighboring pixels. In our research, the change of GLRLM for IACs was significantly higher than for AIS/MIA, which verified the higher heterogeneity or asymmetry for IACs.
We further developed a novel delta radiomics-based nomogram by integrating the DelRADx signature with a radiographic independent predictor for convenient access in clinical applications. Our proposed nomogram showed facilitated performance with an IDI of 0.077 compared with the DelRADx model alone. Meanwhile, the developed tool revealed favorable calibration performance in the calibration curve analysis. Additionally, the decision curve analysis and clinical impact plot illustrated the potential clinical application value of our nomogram. Our findings concurred with the observations that delta radiomics incorporated with clinical data could improve the performance of the prediction model. Khorrami et al. (20) incorporated DelRADx and perinodular morphological characteristics into an individualized decision-making model for patients with advanced NSCLC, and the results showed that the noninvasive approach could effectively identify the candidates who benefit from immunotherapy. Notably, the data used in our model was all clinically accessible, requiring no more investigations and expenditures.
There are several limitations to this study. Firstly, inherent biases are inevitable in this retrospective, single-center analysis. For instance, the PSNs that were examined without CECT were excluded from this study. Secondly, although manual segmentation is regarded as the golden criterion for ROI, it is a labor- and time-intensive method that limits its application to a larger cohort. Thirdly, the reproducibility of radiomic features may face additional challenges due to the variance from different scanner parameters. Last but not the least, the stability of radiomic features may be affected by different acquisition phases of CECT scanning. A multicenter study with a prospective design is planned for our future study.
In conclusion, the delta-radiomics signatures can significantly improve the performance in the differential diagnosis of IACs from AIS/MIAs in patients with PSNs. The delta-radiomics-based nomogram coupled with radiographic features may serve as a convenient way in providing highly informative data for clinical decision support.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.
Ethics statement
This study design was approved by the institutional research ethics board of Huadong hospital, and the informed consent requirement was waived for the retrospective research with anonymous data.
Author contributions
Data curation: DBM, HW. Formal analysis: WFC, ZXM. Funding acquisition: ML. Methodology: RZW. Software: WFC, ZXM. Supervision: YQH. Validation: CL. Visualization: YLY. Writing - original draft: WFC. Writing - review and editing: ML. All authors contributed to the article and approved the submitted version.
Funding
This study was funded by the National Natural Science Foundation of China (grant no.: 61976238) and Research Project Plan of Shanghai Municipal Health Commission (grant no.: 20214Y0309).
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/fonc.2022.927974/full#supplementary-material
References
1. Zhang Y, Fu F, Chen H. Management of ground-glass opacities in the lung cancer spectrum. Ann Thorac Surg (2020) 110:1796–804. doi: 10.1016/j.athoracsur.2020.04.094
2. Yip R, Li K, Liu L, Xu D, Tam K, Yankelevitz DF, et al. Controversies on lung cancers manifesting as part-solid nodules. Eur Radiol (2018) 28:747–59. doi: 10.1007/s00330-017-4975-9
3. Nair A, Bartlett EC, Walsh S, Wells AU, Navani N, Hardavella G, et al. Variable radiological lung nodule evaluation leads to divergent management recommendations. Eur Respir J (2018) 52(6):1801359. doi: 10.1183/13993003.01359-2018
4. Revel MP, Mannes I, Benzakoun J, Guinet C, Léger T, Grenier P, et al. Subsolid lung nodule classification: A CT criterion for improving interobserver agreement. Radiology (2018) 286:316–25. doi: 10.1148/radiol.2017170044
5. Lee KH, Goo JM, Park SJ, Wi JY, Chung DH, Go H, et al. Correlation between the size of the solid component on thin-section CT and the invasive component on pathology in small lung adenocarcinomas manifesting as ground-glass nodules. J Thorac Oncol (2014) 9:74–82. doi: 10.1097/JTO.0000000000000019
6. Ko JP, Suh J, Ibidapo O, Escalon JG, Li J, Pass H, et al. Lung adenocarcinoma: Correlation of quantitative CT findings with pathologic findings. Radiology (2016) 280:931–9. doi: 10.1148/radiol.2016142975
7. Lucia F, Visvikis D, Desseroit MC, Miranda O, Malhaire JP, Robin P, et al. Prediction of outcome using pretreatment (18)F-FDG PET/CT and MRI radiomics in locally advanced cervical cancer treated with chemoradiotherapy. Eur J Nucl Med Mol Imaging (2018) 45:768–86. doi: 10.1007/s00259-017-3898-7
8. Lovinfosse P, Polus M, Van Daele D, Martinive P, Daenen F, Hatt M, et al. FDG PET/CT radiomics for predicting the outcome of locally advanced rectal cancer. Eur J Nucl Med Mol Imaging (2018) 45:365–75. doi: 10.1007/s00259-017-3855-5
9. Boldrini L, Cusumano D, Chiloiro G, Casà C, Masciocchi C, Lenkowicz J, et al. Delta radiomics for rectal cancer response prediction with hybrid 0.35 T magnetic resonance-guided radiotherapy (MRgRT): A hypothesis-generating study for an innovative personalized medicine approach. Radiol Med (2019) 124:145–53. doi: 10.1007/s11547-018-0951-y
10. Rao SX, Lambregts DM, Schnerr RS, Beckers RC, Maas M, Albarello F, et al. CT texture analysis in colorectal liver metastases: A better way than size and volume measurements to assess response to chemotherapy? United Eur Gastroenterol J (2016) 4:257–63. doi: 10.1177/2050640615601603
11. Alahmari SS, Cherezov D, Goldgof D, Hall L, Gillies RJ, Schabath MB. Delta radiomics improves pulmonary nodule malignancy prediction in lung cancer screening. IEEE Access (2018) 6:77796–806. doi: 10.1109/ACCESS.2018.2884126
12. Kamiya S, Iwano S, Umakoshi H, Ito R, Shimamoto H, Nakamura S, et al. Computer-aided volumetry of part-solid lung cancers by using CT: Solid component size predicts prognosis. Radiology (2018) 287:1030–40. doi: 10.1148/radiol.2018172319
13. Travis WD, Asamura H, Bankier AA, Beasley MB, Detterbeck F, Flieder DB, et al. The IASLC lung cancer staging project: Proposals for coding t categories for subsolid nodules and assessment of tumor size in part-solid tumors in the forthcoming eighth edition of the TNM classification of lung cancer. J Thorac Oncol (2016) 11:1204–23. doi: 10.1016/j.jtho.2016.03.025
14. Lee SM, Goo JM, Lee KH, Chung DH, Koh J, Park CM. CT findings of minimally invasive adenocarcinoma (MIA) of the lung and comparison of solid portion measurement methods at CT in 52 patients. Eur Radiol (2015) 25:2318–25. doi: 10.1007/s00330-015-3616-4
15. Luo T, Xu K, Zhang Z, Zhang L, Wu S. Radiomic features from computed tomography to differentiate invasive pulmonary adenocarcinomas from non-invasive pulmonary adenocarcinomas appearing as part-solid ground-glass nodules. Chin J Cancer Res (2019) 31:329–38. doi: 10.21147/j.issn.1000-9604.2019.02.07
16. Weng Q, Zhou L, Wang H, Hui J, Chen M, Pang P, et al. A radiomics model for determining the invasiveness of solitary pulmonary nodules that manifest as part-solid nodules. Clin Radiol (2019) 74:933–43. doi: 10.1016/j.crad.2019.07.026
17. Basler L, Gabryś HS, Hogan SA, Pavic M, Bogowicz M, Vuong D, et al. Radiomics, tumor volume, and blood biomarkers for early prediction of pseudoprogression in patients with metastatic melanoma treated with immune checkpoint inhibition. Clin Cancer Res (2020) 26:4414–25. doi: 10.1158/1078-0432.CCR-20-0020
18. Vinod SK, Hau E. Radiotherapy treatment for lung cancer: Current status and future directions. Respirology (2020) 25 Suppl 2:61–71. doi: 10.1111/resp.13870
19. Son JY, Lee HY, Kim JH, Han J, Jeong JY, Lee KS, et al. Quantitative CT analysis of pulmonary ground-glass opacity nodules for distinguishing invasive adenocarcinoma from non-invasive or minimally invasive adenocarcinoma: The added value of using iodine mapping. Eur Radiol (2016) 26:43–54. doi: 10.1007/s00330-015-3816-y
20. Khorrami M, Prasanna P, Gupta A, Patil P, Velu PD, Thawani R, et al. Changes in CT radiomic features associated with lymphocyte distribution predict overall survival and response to immunotherapy in non-small cell lung cancer. Cancer Immunol Res (2020) 8:108–19. doi: 10.1158/2326-6066.CIR-19-0476
Keywords: PSNs, lung adenocarcinoma, CECT, delta-radiomics, LASSO-logistic
Citation: Chen W, Wang R, Ma Z, Hua Y, Mao D, Wu H, Yang Y, Li C and Li M (2022) A delta-radiomics model for preoperative prediction of invasive lung adenocarcinomas manifesting as radiological part-solid nodules. Front. Oncol. 12:927974. doi: 10.3389/fonc.2022.927974
Received: 25 April 2022; Accepted: 24 October 2022;
Published: 16 November 2022.
Edited by:
Jianfei Shen, Wenzhou Medical University, ChinaReviewed by:
Ying Ji, Beijing Chaoyang Hospital, Capital Medical University, ChinaZhiyong Yuan, Tianjin Medical University Cancer Institute and Hospital, China
Copyright © 2022 Chen, Wang, Ma, Hua, Mao, Wu, Yang, Li 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: Wufei Chen, Y2hlbnd1ZmVpXzIwMDhAMTYzLmNvbQ==; Ming Li, bWlubGk3N0AxNjMuY29t