Radiomic advances in the transarterial chemoembolization related therapy for hepatocellular carcinoma

Tian-You Chen, Zong-Guo Yang, Ying Li, Mao-Quan Li

Tian-You Chen, Department of Interventional Medicine, Shanghai Public Health Clinical Center,Fudan University, Shanghai 201508, China

Zong-Guo Yang, Department of Integrative Medicine, Shanghai Public Health Clinical Center,Fudan University, Shanghai 201508, China

Ying Li, Department of Radiology, Jinshan Hospital, Fudan University, Shanghai 201508, China

Mao-Quan Li, Department of Interventional & Vascular Surgery, Tenth People's Hospital of Tongji University, Tongji University, Shanghai 200433, China

Abstract Radiomics is a hot topic in the research on customized oncology treatment, efficacy evaluation, and tumor prognosis prediction. To achieve the goal of mining the heterogeneity information within the tumor tissue, the image features concealed within the tumoral images are turned into quantifiable data features. This article primarily describes the research progress of radiomics and clinicalradiomics combined model in the prediction of efficacy, the choice of treatment modality, and survival in transarterial chemoembolization (TACE) and TACE combination therapy for hepatocellular carcinoma.

Key Words: Transarterial chemoembolization; Hepatocellular carcinoma; Radiomics;Magnetic resonance imaging; Computed tomography

INTRODUCTION

Approximately 90% of all liver tumors are primary hepatocellular carcinomas (HCC), which are the most frequent primary liver tumors[1]. Among the most prevalent malignant tumors worldwide, HCC is the sixth most prevalent. HCC ranks third in the mortality rate of malignant tumors, which has a significant impact on people's quality of life and physical well-being[2]. According to statistics from the World Health Organization 2018, there are 840000 new instances of HCC each year, of which 780000 cases die[3].

The occurrence and progression of HCC is a complex process with multi-factor participation and multi-step interconnection. The majority of patients would lose their opportunity for radical treatment, such as surgical resection and liver transplantation, due to its insidious onset, atypical symptoms and signs, low early diagnosis rate, rapid disease progression, tumor multicentric origin, abundant blood supply, easy to form cancer thrombus, cause metastasis, highly heterogeneous tumor, and other malignant biological behaviors. The selection of HCC treatments is also constrained by characteristics other than the tumor itself, including the degree of cirrhosis, liver function reserve, and portal hypertension. The effectiveness of a single surgical procedure has demonstrated a "ceiling impact" in light of the HCC's biological nature, which is very aggressive.

With the continuous development of interventional oncology therapeutics, interventional therapy for HCC has developed from straightforward transhepatic arterial therapy to a diversified, combined treatment mode, covering early-, intermediate- and advanced-stage HCC. In the field of minimally invasive treatment of HCC, transarterial chemoembolization (TACE) is the most commonly used technique, which uses lipiodol and/or drug-eluting beads. In process of TACE, cytotoxic drugs like doxorubicin or cisplatin are injected while being emulsified in the radiopaque oil compound lipiodol. After that, the embolic agents, like a gelatin sponge, are given. Lipiodol promotes embolization of the tumor microcirculation and delivers cytotoxic drugs directly to the tumor[4]. Nevertheless, there is no acceptable approach to accurately assess the efficiency of interventional therapy due to the significant individual discrepancies presently.

Radiomics is a technique that combines high-throughput extraction of advanced quantitative features to objectively and quantitatively define tumor behaviors. It was first proposed in 2012 by Dutch researcher Lambin[5]. Radiomics extract a huge number of features from medical pictures in high throughput and turn imaging data into a high-resolution mining data space by automatic or semiautomatic analysis methods[5,6]. It is possible to comprehensively, non-invasively, and quantitatively monitor the spatial and temporal heterogeneity of tumors. Radiomics can be divided into below parts: (1) Image acquisition; (2) Region of interest (ROI) segmentation; (3) Feature Extraction; (4) Feature Selection; (5) Model establishment and evaluation; and (6) Clinical application[5,6]. Radiomics has no restrictions on the sorts of medical pictures used for sample collection, and the objects of analysis can include ultrasound, X-ray, CT, magnetic resonance imaging (MRI), positron emission tomography, and other medical images. The ROI can be altered depending on the clinical situation, whether it is normal tissue, the full tumor area, or a metastatic lesion. Generally, both manual and automatic segmentation of the ROI is used in radiomics. However, the segmentation effect of these automatic segmentation techniques is unstable. Currently, manual ROI segmentation by skilled imaging specialists is the approach most frequently used in medical imaging toolkits (such as ITK-SNAP segmentation software). Nevertheless, manual segmentation takes a lot of time, costs a lot of money, and is susceptible to subjective factors despite having a high accuracy rate. Extraction of high-throughput features for quantitative analysis of the important ROI attributes is the first stage in radiomics. The following three categories generally correspond to the extracted features: (1) Shape feature: Reflecting the information of tumor morphology, size, and regularity; (2) Histogram features: The most basic statistical descriptor is based on a global histogram of grayscale values, which includes the grayscale mean, maximum, minimum, variance, and percentile. These features are also known as first-order features because they are based on single-pixel or monomer element analysis; and (3) Texture feature: characteristics can measure details that are hard for the eye to easily see, such as tissue distribution inside the tumor or textural patterns. Such high-dimensional abstract features like the Baud sign may play varied roles to capture clinical information that is not readily perceptible by the vision for clinical problems that are difficult to convey by simple visual features of tumor images. Many hundred to tens of thousands of features may be obtained by feature extraction for feature screening, but not all of them are connected to the clinical issues that need to be resolved. On the other hand, because there are so few samples and so many characteristics in practice, it is simple to cause the overfitting phenomena in later models. Hence, choosing the extracted features is a crucial step in the study of radiomics. The commonly used feature selection methods include filtering, packaging, and embedding. The last step is model construction and model evaluation. To establish the key features of the above feature screening based on clinical research questions, or to further combine features other than imaging omics (such as clinical signs, pathology, and genetic test data) to establish predictive models (Figure 1). According to previous studies[5-8], radiomics properties are directly associated with tumor aggressiveness and heterogeneity indices at the cellular level. Radiomics is believed to be able to forecast clinical endpoints including survival and therapeutic response. Second, radiomics data are mineable, which means that, in sufficiently big data groups, they may be utilized to find previously unidentified indicators and patterns of illness evolution, progression, and treatment response.

Figure 1 Radiomics flowchart. Based on the imaging data obtained through scanning, the region of interest is delineated on the image, corresponding imaging features are extracted, features are filtered and dimensioned, and imageomics labels are constructed to obtain the rad-score of the image. Then, combined with clinical features, a joint prediction model is established, and a series of clinical applications are conducted by constructing nomogram, calibration curve, and receiver operating characteristic curve. ROI: Region of interest; ROC: Receiver operating characteristic; FPT: Ferroptosis; MRI: Magnetic resonance imaging.

The aim of this review was to summarize the predictive value of radiomics in TACE alone and combination therapy for HCC and its guiding significance for the formulation of treatment regimens.

PREDICTION OF EFFICACY

Tumors less than 5 cm

Kuanget al[9] retrospectively collected MRI images and clinical data of 153 patients from 3 hospitals, and divided the subjects into 1 internal training group (n= 113) and 1 external validation group (n= 40) to establish MRI radiomics nomograms that can forecast a patient's short-term response to TACE for an HCC with a diameter under 5 cm. By examining pictures [enhanced computed tomography (CT)/magnetic resonance (MR) scans] of all patients 3/4 mo after their initial TACE treatment, the response of HCC patients to TACE was assessed using the modified Response Evaluation Criteria in Solid Tumors (mRECIST). All cases were split into two groups based on the evaluation results: those who responded well (complete response, CR, and partial response, PR), and those who did not (stable disease, SD, and progressive disease, PD). Images from dynamic enhanced MRI (DCE-MRI) arterial phase and T2-weighted imaging (T2WI) were used to establish radiomics model, clinical model, and radiomics nomogram (by combining the radiomics and clinical model). Radiomics selected the image features using the minimum absolute contraction and selection operator (LASSO) regression and the maximum correlation-minimum redundancy (mRMR) algorithm and then selected 11 of the best subsets of features based on the arterial phase to calculate the radscore. The clinical model based on T2WI was ultimately developed by combining platelet count (PLT), pseudocapsule, boundary, and peritumoral augmentation, and well-response was suggested by a normal PLT value (＞ 100 × 109/L), pseudocapsule, and obvious border on the T2WI sequence. Another clinical model based on DCE-MRI arterial phase (AP) was developed by combining Child-Pugh class, border, and peritumoral augmentation. Low ChildPugh class (class A), clear border, and no peritumoral enhancement on AP all point to wellresponse. In contrast with radiomics and clinical models alone, the results showed that the nomogram had superior predictive potential for postoperative response to TACE. Clinics models and Radiomics were obtained which could predict the postoperative response of TACE, and the maximum area under the curve (AUC) value was 0.76, and 0.78 respectively. When combing Clinics models with Radiomics models, predictive power peaked at 0.84. The T2WI-Nomogram and AP-Nomogram were the models with the highest clinical decision effectiveness, according to a comparison of each model's clinical decision effectiveness. Furthermore, T2WI-based and DCE-MRI-based nomograms had approximately comparable predictive abilities. It is demonstrated that the T2WI sequence may effectively represent the heterogeneity of lesions and that adopting a T2WI sequence can result in the formulation of a stable prediction model, eliminate the need for enhanced scanning, lower the side effects of unnecessary contrast injections and medical costs. Patients included in this study were under MRI scanning from different manufacturers (GE and Siemens) with various field strengths (1.5T and 3.0T) and with various sequence parameters. It was discovered through analyzing and processing data that various machines and magnetic field strengths had a mild influence on the radiomics findings. This result indicated that the radiomics has a high degree of repeatability and stability.

Barcelona clinical liver cancer B

TACE offers a considerable survival advantage over supportive therapy for patients with stage B HCC who are initially diagnosed with Barcelona clinical liver cancer (BCLC). These patients had significant survival heterogeneity following TACE despite getting comparable therapies. Some pre-existing prognostic models have been put forth to forecast clinical outcomes following TACE, including Six-andtwelve score[10], Four-and-seven score[11], Hepatoma arterial-embolization prognostic (HAP score)[12], modified hepatoma arterial embolization prognostic (mHAP score)[13], modified hepatoma arterial embolization II (mHAP-II score)[14], modified hepatoma arterial embolization III (mHAP-III score)[15] and albumin-bilirubin(ALBI) grade[16] for arterial embolization of liver cancer. Using minimum absolute contraction and selection operator approaches, the author of this study creates radiomic signatures for survival (rad signatures). The most predictive 6 radiomic features were ultimately chosen. The remaining 4 features were generated from tumor variants of interest (VOI) imaging in the portal vein phase, leaving 2 of the 6 features based on tumor VOI and peritumoral VOI imaging in the arterial phase. Univariate Cox regression analyses were used to ascertain prognostic clinical factors. With multivariate Cox regression analyses, a combined radiomics-clinic (CRC) model was developed using the Rad-signature and clinical factors with a potential association with overall survival (OS). And then contrasted with a prognostic model that had previously been published. Researchers[17] conducted a retrospective multicenter analysis, and the outcomes revealed that the CRC model which was finally established with tumor number (＜ 4vs≥ 4) and the Rad-signature performed better than the other previously mentioned seven established prognostic models, with Harrell’s concordance-index (C-index) of 0.73 (95%CI: 0.68–0.79) and 0.70 (95%CI: 0.62–0.82) in the training and testing cohort, respectively. Among the seven models tested, the six-and-12 score and four-and-seven criteria performed better than the other models, with C-indices of 0.64 (95%CI: 0.58–0.70) and 0.65 (95%CI: 0.55–0.75) in the testing cohort, respectively. The authors created customized risk scores by linearly combining the radscore and the number of tumors (＜ 4vs≥ 4), weighted by their respective coefficients in the multivariate Cox regression model, in order to aid clinical practice. Patients were separated into two groups based on their median risk scores for the training cohort (0.0214): stratum 1, with a risk score of ＜ -0.0214, and stratum 2, with a risk score of ＞ -0.0214. In the training cohort, stratum 1 patients had a considerably longer median survival (31.3 mo) than stratum 2 patients (12.5 mo), with a hazard ratio of 3.63 (95%CI: 2.36-5.60, log-rank testP＜ 0.0001). Applying the same cutoff to the testing cohort the hazard ratio was 2.43 (95%CI: 1.91-4.98,P= 0.0014), and the median survival for the two groups was 30.9 and 17.0 mo, respectively. The CRC model showed improved survival predictive performance, and researchers believe that the CT radiomics signature represents an independent biomarker of survival in patients with HCC undergoing TACE.

Sunet al[18] used five distinct MR machines with various field strengths for traditional T1WI, T2WI, and T2 grease pressing as well as diffusion-weighted imaging (DWI) (b = 0 and b = 500) scanning. A total of 84 BCLC B HCC patients were enrolled, of whom 67 were enrolled in the training group and 17 in the test group. The mRECIST was used to assign patients into the disease progression (PD) group and nonprogressive disease (N-PD) group if they progressed 6 mo following TACE treatment. In the test group, the calculated AUC of the DWI feature-based model for TACE efficacy prediction is 0.79 (b = 0) and 0.73 (b = 500), followed by T2WI (0.729) and ADC (0.71). In comparison to any single MRI-based radiomics feature, the multiparametric MRI-based radiomics features' AUC rises to 0.80. While clinical data including age, α-fetoprotein (AFP), and tumor size were not statistically significant, a multivariate logistic analysis could identify the independent prognostic indicators of the radiomics signature of PD.

According to mRECIST 1.1, Baiet al[19] aim to evaluate the effectiveness of radiomics features based on preoperative contrast-enhanced computed tomography (CECT) in predicting response to transarterial chemotherapy (TACE), categorized 111 patients with intermediate HCC who received CE-CT in the arterial phase and venous phase before and following TACE into objective response group (n= 38) and non-response group (n= 73) groups. Radiomic feature extraction from CECT pictures. The best monophasic radiomic features of AP and VP in the training set were discovered using two feature ranking methods and three classifiers. In parallel, decision level fusion and feature level fusion were used to combine the pictures of the two CECT phases and establish the multiphase radiomic properties. For the AP signature, the combination of minimum redundancy maximum relevance (MRMR), and support vector machine (SVM) showed the best performance (AUC = 0.814). For VP signatures, the best performance was obtained by MRMR and LASSO (AUC = 0.861). For the performance of multiphase radiomic signatures, DLF signatures had the highest AUC value of 0.883 among all radiomic signatures by using the features selected by MRMR and SVM classifiers. Eventually, a nomogram was constructed by combining two common features (tumor size and tumor number) and radiomics signature, using multivariate logistic regression, and its prediction ability was assessed by AUC on the test data. To determine the accuracy, sensitivity, and specificity values, the cut-off point closest to the upper left corner of the training receiver operating characteristic (ROC) was used. Scores below or above the cutoff were considered to be either objective-response or non-response. It was discovered that multiphase radiomics features (AUC = 0.883) outperformed the best single-phase radiomics signature (AUC = 0.861) in predicting response to TACE treatment. The nomogram in the test dataset and training dataset showed better performance than any radiomics signatures. The authors conclude that the radiomic model will be helpful in assessing the therapeutic advantages of TACE treatment.

Barcelona clinical liver cancer B+C

Konget al[20] enrolled 99 patients with intermediate and advanced HCC (BCLC B and C) treated with TACE. Among them, 69 cases were enrolled in the training group and 30 cases in the test group. MRI was performed before TACE, and efficacy was evaluated according to mRECIST criteria 3 mo after TACE. The patients were divided into response groups (CR and PR) and non-response groups (PD and SD). There were significant differences in AFP values, Child-Pugh scores, and BCLC staging between the two groups. Before TACE, 396 radiomics characteristics were extracted from T2WI, and 6 features were chosen using minimal absolute contraction and selection operator (LASSO) regression. These features are variance, inverse difference moment angle90 offset7, long run emphasis angle90 offset7, short run emphasis all direction offset 4 SD, short run emphasis angle 135 offset1, and high intensity large area emphasis. They are combined to provide radscore to create response prediction models. In the training and validation cohorts, the AUCs of the ROC curve based on Radscore were 0.812 and 0.866, respectively. While the prediction performance of the radiomics model was greatly enhanced when the Child-Pugh class, BCLC stage, and AFP level were added. The AUC of the training group was 0.861 (95%CI: 0.774–0.949) with specificity and sensitivity as high as 0.811 and 0.844, respectively. The AUC of the validation group was 0.884 (95%CI: 0.764–1.000), and the specificity and sensitivity were 0.75 and 1.00, respectively. Based on the above conclusion, the authors believe that the quantitative nomogram model based on radscore and clinical predictors can accurately predict the response to TACE in intermediate and advanced HCC and can also be used as an auxiliary diagnosis for clinical prognosis.

Barcelona clinical liver cancer A+B

Retrospectively, 122 HCC (BCLC A and B) patients who had MRI scanning prior to initial TACE (response = 63; non-response = 59) were recruited by Ying Zhaoet al[21] and randomly divided into training group (n= 85) and validation group (n= 37). The manual segmentation of arteries, veins, and delayed periods on CE-MRI yielded a total of 2367 radiomics characteristics. Of all radiomics features, the three-phase radiomics model performed better in the training group, with an AUC of 0.838. The AUC of radiomics nomogram combining three-phase radiomics score and clinical-radiological risk factors (total bilirubin, tumor shape, and tumor encapsulation) were 0.878 and 0.833, respectively, which showed good calibration and prediction with better predictive power. However, other studies[22,23] have demonstrated that while total bilirubin can predict survival in HCC patients treated with TACE, it cannot be utilized as a clinical risk factor to evaluate the effectiveness of treatment.

CHOICE OF TREATMENT MODALITY

Almost all the guidelines advocate transitioning to systemic medication or combination therapy when a patient's response to TACE is poor. To tackle this problem, the concept of "TACE refractoriness" is put forth[24,25]. According to the Japanese Liver Association[26], vascular invasion, distant metastases, persistently increased tumor markers or an ongoing partial response to at least two TACEs are indicative of intractable resolution of TACE. New vascular invasion or distant metastases can operate as stand-in markers for refractory TACE during TACE therapy, according to a Korean study[27]. It's critical to anticipate TACE refractoriness while considering a potentially better treatment. The ability to predict TACE refractoriness before surgery has been demonstrated by some researchers utilizing the radiomics nomogram approach on arterial CT scans of 80 HCC patients who did not have extrahepatic metastases or macrovascular involvement[28]. The authors believe that the radscore, which consists of the Gray-Level Zone Length Matrix (GLZLM)-Long-Zone Low Gray-Level Emphasis and GLZLM-Gray-Level Non-Uniformity, T-stage, log AFP, and bilobar distribution is significantly related to TACE refractoriness. Multivariate logistic regression model reveals that these factors, such as radscore+T-stage (AUC = 0.95); radscore+bilobar distribution (AUC = 0.91); and radscore+logAFP (AUC = 0.91), perform well in predicting TACE refractoriness. Additionally, according to Niuet al[29], the radiomics nomogram of CT can be used to forecast the pretreatment prediction of TACE refractoriness and offer better guidance for choosing the next course of TACE treatment.

SURVIVAL ANALYSIS

A multicenter retrospective analysis included 140 patients with TACE HCC (BCLC A+B+C)[30]. In order to predict tumor response, response-related radscores were created using T2WI and DCE-MRI, respectively. The T2WI outperformed 4 MRI sequences with an AUC of 0.75. The T2WI radscore, BCLC stage, and albumin-bilirubin grade were combined to generate a nomogram using the logistic regression model. The nomogram yielded AUCs of 0.81 and 0.78 for the train and test groups for predicting tumor response. Progression-free survival (PFS) and OS in survival analyses showed substantial differences between responders and non-responders of TACE treatment. Albumin bilirubin grading, satellite nodes, and BCLC staging are independent predictors of OS. PFS and OS were independently predicted by albumin bilirubin grading and BCLC staging. According to the Cox model, PFS and OS's consistency indices were 0.70 and 0.73, respectively.

Another study extracted radiomics features from intratumoral and perineoplastic dilatation (1, 3, and 5 mm, respectively) from 184 patients with HCC treated with TACE. The radiomics model was constructed using LASSO[31], which was used to investigate the relationship between preoperatively CE-MRI and recurrence-free survival (RFS) in HCC patients regarding the RFS of HCC patients after TACE. Seven models of radiomics (APETV, PVPETV, PVPB1, PVPB3, and PVPB5, in which the TV stands for the total tumor volume and B represents millimeters beyond the tumor boundary), clinicalradiology parameters, and a combination model (nomogram) were generated. The results showed that sex, AFP, BCLC stage, IM (defined as a tumor with a budding portion at its periphery in both the AP and PVP), and APE (which was interpreted as the contrast enhanced shadow around the tumor in the AP, and isointensity compared with the background liver parenchyma in the PVP and delayed phase) were independent risk factors related to RFS after TACE. The performance of the PVPB3 model was 0.714 (C-index), which was equivalent to that of the portal stage intratumor model (C-index, 0.727,P= 0.409). The nomogram, which had a C-index of 0.802 with the best performance, incorporated five independent risk factors from PVPETV characteristics and clinical-radiological parameters. It was more useful for assessing RFS in HCC patients following TACE treatment. The combined model's cutoff relative to the median (1.74) allowed for a precise division of these individuals into high- and low-risk groupings, according to the Kaplan-Meier analysis. A study[29] based on radiomics features by CE-CT in the arterial phase had also revealed that radscore, BCLC stage, irregular tumor margin, largest tumor size, and tumor number were independent influencing factors for OS. CT-based radiomics nomogram could classify patients into high-risk (＞ 3.5) and low-risk (≤ 3.5) groups with significantly different prognoses between the two groups (OS: 12.3 movs23.6 mo,P＜ 0.001). Another study came to the same conclusion that the combined model (nomogram) included five radiomics features (surface area-tovolume ratio, kurtosis, median, gray-level co-occurrence matrix contrast, and size zone variability) and three clinical factors (Child-Pugh score, α-fetoprotein level, and tumor size) was better in predicting survival than either the clinical scoring model or the radscore model (hazard ratio 19.88;P＜ 0.001).

COMBINATION THERAPY

TACE + radio frequency ablation (RFA)

TACE and RFA combination therapy have synergistic cytotoxic effects for intermediate and advanced HCC, compared to TACE alone, which can better control local tumor recurrence and lengthen survival[4]. TACE lessens the "heat sink" that blood flow creates during radiofrequency ablation. It aggravates the edema of tissues, and clearly displays the outlines of lesions induced by iodine oil deposition on CT, making it easier to determine the RFA range. These elements may increase the severity of coagulation necrosis and, as a result, improve OS prominently[32]. However, there are also significant differences in the efficacy of TACE combined with RFA. Some scholars[33] combined radiomics and clinical factors to construct a predictive nomogram model, which can be used as a new strategy for predicting advanced HCC PFS in the treatment of TACE+RFA.

TACE+lenvatinib

Lenvatinib is just as efficient as sorafenib as a multikinase inhibitor with antiangiogenic effects in the first-line treatment of advanced HCC[34]. Combination therapy with TACE and lenvatinib has been shown to be effective in recent studies[35,36]. According to the findings, combination therapy outperformed lenvatinib monotherapy and TACE by a wide margin. Luoet al[37] studied 61 HCC patients and extracted radiomics characteristics from T1WI, T1WI arterial phase, T1WI portal phase, T1WI delay phase, T2WI, DWI (b = 800), and ADC maps. The results showed that the number of tumors [risk ratio, HR = 4.64, 95%CI: 1.03-20.88,P= 0.045] and arterial phase intensity HR = 0.24, (95%CI: 0.09-0.64) with a sensitivity of 0.99, specificity of 0.95, and accuracy of 0.71. The combined model significantly increased the clinical model's performance in predicting disease progression.

CURRENT LIMITATIONS

We discussed the use of radiomics in predicting TACE treatment in HCC, which quantifies the heterogeneity of HCC on medical images. These radiomics features can be further used as predictors of clinical outcomes to determine the efficacy of TACE, as well as provide individualized guidance for patient treatment options. However, there are also some difficulties need to be further resolved. First, it might be challenging by using a computer to pinpoint the tumor boundary automatically, and hand outlining may vulnerable to significant subjective variation. Second, considering radiomics is still in its infancy and different researchers employ different imaging tools and image algorithms, it is impossible for other researchers to evaluate and put the results to clinical use[38]. Last but not the least, it is challenging to create a worldwide uniform paradigm for radiomics in HCC research because the underlying etiology of the disease varies between different regions[39].

CONCLUSION

Radiomics can forecast the efficacy of TACE and TACE combined therapy regimens before the treatment conducting. Tailored counsel for patient treatment based on radiomics and increased comprehensive prediction abilities would be achieved when combined with clinical data.

FOOTNOTES

Author contributions:Chen TY, Yang ZG, Li Y, Li MQ equally contributed to this paper with conception and design of the study, literature review and analysis, drafting and critical revision and editing, and final approval of the final version.

Conflict-of-interest statement:All the authors declare that they have no conflicting interests.

Open-Access:This article is an open-access article that was selected by an in-house editor and fully peer-reviewed by external reviewers. It is distributed in accordance with the Creative Commons Attribution NonCommercial (CC BYNC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited and the use is noncommercial. See: https://creativecommons.org/Licenses/by-nc/4.0/

Country/Territory of origin:China

ORCID number:Tian-You Chen 0000-0002-8668-6650; Zong-Guo Yang 0000-0002-6623-4841; Ying Li 0000-0001-5062-0230; Mao-Quan Li 0000-0001-7849-7350.

S-Editor:Liu JH

L-Editor:A

P-Editor:Zhao S