Abstract

Background. Intrahepatic cholangiocarcinoma (ICC) is the second most common liver malignancy after hepatocellular carcinoma (HCC), with a dismal prognosis and high heterogeneity. The oncological advantages of anatomical resection (AR) and nonanatomical resection (NAR) in HCC have been studied, but surgical strategies for ICC remain controversial with insufficient investigations. Materials and Methods. From Jan 2013 to Dec 2016, 3880 consecutive patients were retrospectively reviewed from a single center. Patients with ICC undergoing AR or NAR have been enrolled according to inclusion and exclusion criteria. Propensity score matching (PSM) analysis was performed between two groups with a 1 : 1 ratio. The primary endpoint was overall survival (OS), and the secondary endpoints included disease-free survival (DFS), intraoperative patterns, postoperative morbidity, mortality, complications and recurrence. A prognostic nomogram was developed by a multivariate Cox proportion hazard model. Results. After PSM, 99 paired cases were selected from 276 patients enrolled in this study. Patients in the AR group achieved better 1-, 3-, and 5-year OS (70%, 46%, and 34%, respectively) and DFS (61%, 21%, and 10%, respectively) than patients in the NAR group with statistical significance after PSM analysis. The postoperative complications and recurrence patterns were comparable between the two groups. Multivariate analysis identified NAR, tumor size >5 cm, multiple tumors, and poor differentiation as independent risk factors for OS (). Selected patients can benefit most from AR, according to subgroup analysis. A prognostic nomogram based on six independent risk factors for OS and factors with clinical significance was constructed to predict OS in ICC patients. Conclusion. AR improved the long-term survival of ICC with comparable postoperative complications and similar recurrence patterns. AR is suggested in ICC patients with sufficient remnant liver volume. In addition to surgery strategy, malignant characteristics of tumors are risk factors for ICC prognosis.

1. Introduction

Cholangiocarcinoma (CCA) is a highly lethal hepatic malignancy with increasing incidence and mortality worldwide, with 0.3–6 cases per 100,000 inhabitants per year in Western and over 6 cases in some East Asian regions [13]. This heterogeneous cancer with aggressive invasiveness can be classified as intrahepatic CCA (iCCA or ICC), perihilar CCA (pCCA), or distal CCA (dCCA), according to the anatomic structure [46]. Though advances in neoadjuvant chemotherapy and targeting immunotherapy have brought scientific and clinical insights into treatments for ICC, there is still a lack of decisive evidence for their clinical application [7, 8].

Surgical resection still remains the only potentially curative treatment, with low resectability rates varying from 19% to 74% [9]. The surgery strategies and resection range are controversial for ICC in speaking of prognosis improvement. Resection margin and lymph node metastasis are significant prognostic factors, but it is debated whether wide resection margins or lymphadenectomy improved postoperative survival [10, 11]. Aggressive invasiveness characteristics of ICC present multifocality, lymph node metastasis, and vascular invasion with dismal outcomes [12]. The optimal therapeutic strategy for ICC has not been determined yet [13]. Appropriate preoperative surgical planning and subsequent treatments decided by multidisciplinary team (MDT) discussion are indispensable in achieving optimal outcomes for ICC patients [13].

Anatomical resection (AR) was firstly proposed by Japanese surgeon Makuuchi in the 1980s to eradicate potential micrometastasis [14]. Studies for oncological advantages of AR in primary liver cancer mainly focus on hepatocellular carcinoma (HCC), though the results have been long controversial [15, 16]. Compared with intrahepatic metastasis of HCC via vascular structures, tumor cells of ICC spread along with the biliary tree and lymph system, which supports the oncological and prognostic superiority of AR over nonanatomical resection (NAR) in theory [17]. The complete excision of tumor-bearing portal tributaries of AR in ICC may prevent microscopic intrahepatic metastasis by preventing tumor cells from spreading along the biliary tree or lymph system and reducing the rates of local recurrence and prolonging survival. The efficacy of AR for ICC may also vary according to clinicopathological factors, including tumor size, number, location, differentiation, preoperative liver function, and remnant liver parenchyma.

However, the impact of AR and NAR on short-term and long-term outcomes in ICC still lacks investigation [1820]. We conducted this retrospective study to evaluate the operative and postoperative characteristics of AR and NAR groups. Propensity score matching (PSM) analysis was applied to minimize the selection bias of the surgery strategy.

2. Materials and Methods

2.1. Patient Enrollment

From Jan 2013 to Dec 2016, 3880 consecutive patients were retrospectively reviewed in the Institute of Hepato-Pancreato-Biliary Surgery, Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology. The inclusion criteria were as follows: (1) patients with ICC confirmed pathologically; (2) patients with age ≥18 years; (3) patients initially treated with AR or NAR; (4) patients with R0 resection margin. The exclusion criteria were as follows: (1) patients with severe underlying diseases; (2) patients with other malignancies; (3) patients with metastasis prior to the resection; (4) patients who received previous antitumor treatments; (5) patients undergoing hepatectomy combined with transarterial chemoembolization (TACE) or ablation; (6) patients with microscopically positive (R1) resection or macroscopically positive (R2) resection margin; (7) patients who died within 30 days after surgery or due to other nontumor causes; (8) patients with insufficient clinical data.

This study was reviewed and approved by the Medical Ethics Committee of Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology. Written informed consent was obtained from the patients.

2.2. Liver Resection and Follow-Up

Surgery strategy of AR and NAR was determined based on a general evaluation of jaundice, serum level of total bilirubin, and remnant liver volume assessed by CT scanning and three-dimensional reconstruction techniques as well as indocyanine green retention at 15 min (ICGR15). AR was defined as the complete resection of Couinaud’s segments by prior ischemia or dye staining of indocyanine green (ICG) fluorescence, including segmentectomy, lobectomy, or hemihepatectomy. NAR was defined as incomplete resection of the portal tributaries of the lesion segment with a surgical margin of at least 1 cm or without exposing the tumor surface on the parenchymal transection, including partial resection or enucleation of the liver tumor [15, 17, 21]. Intraoperative ultrasound was routinely performed to evaluate tumor location, number, size and adjacent parenchyma, in addition to preoperative imaging. Pringle’s maneuver was routinely performed with interval clamping/unclamping of 10 min/5 min. Portal occlusion and/or inferior vena cava (IVC) occlusion were applied when necessary. Routine lymphadenectomy at the level of hepato-duodenal ligament has been performed according to guidelines [22]. All patients in this study have achieved R0 resection.

The first follow-up was carried out 1 month after the operation, and every 2-3 months within the first year, then every 6–12 months afterwards. Physical examination and measurements of liver and kidney function, serum level of tumor markers (carbohydrate antigen 19–9 (CA19-9), carcinoembryonic antigen (CEA), and alpha-fetoprotein (AFP)), and imaging examination (abdominal ultrasound, contrast-enhanced computed tomography (CT), magnetic resonance imaging (MRI), and chest X-ray) were performed during the follow-up. When recurrence occurred during the follow-up, reoperation, microwave ablation, transarterial chemoembolization (TACE), chemotherapy, radiotherapy, targeted therapy, immunotherapy or palliative care were given according to clinical guidelines, MDT discussion, and the patients' wishes [6, 2225]. The last follow-up date was Jan 2022.

2.3. Study Endpoints

The primary endpoint was overall survival (OS). The secondary endpoints included disease-free survival (DFS), intraoperative parameters (operation time, blood loss, blood transfusion, portal vein/IVC occlusion and laparoscopic/open surgery), postoperative morbidity, mortality, complications, and recurrence (intrahepatic, extrahepatic, or both). OS was defined as the time from the date of surgery to either the date of death or the last date of follow-up, while DFS was defined as the time from the date of surgery to either the date of disease recurrence or the last date of follow-up. Postoperative morbidity and mortality were defined as events that happened within the first 60 days after surgery. Complications were graded according to the Clavien–Dindo classification [26].

2.4. Propensity Score Matching Analysis and Nomogram Prediction

PSM analysis was introduced to reduce the bias of treatment selection. As previously described, the ICC patients in AR and NAR groups were matched by propensity score [27]. The propensity score for an individual was calculated given the covariates of tumor number, Child–Pugh classification, lymph node metastasis, and portal vein tumor thrombosis (PVTT) in pathology using a logistic regression model. Thereafter, 1 : 1 nearest neighbor matching with a calliper of 0.01 and without replacement was applied to ensure minimized conditional bias [28].

Independent risk factors selected by multivariate analysis for OS and the demographic characteristics with clinical significance were incorporated into the nomogram to predict 1-, 3-, and 5-year survival.

2.5. Statistical Analysis

Categorical variables are presented as numbers (%) and were compared by the Chi-squared test. Continuous variables are presented as mean ± standard deviation (SD). Normally distributed data were compared with Student’s t-test, and nonnormally distribution was analyzed by the Mann–Whitney U test. The Kaplan–Meier method was applied to depict the survival curves before and after PSM, compared using the Log-rank test. Candidate variables with statistical significance in univariate analysis were introduced into multivariable Cox proportional hazards regression analyses to determine the independent risk factors associated with OS and DFS in the PSM cohort. In order to evaluate the impact of surgical strategy (AR or NAR) on OS, the groups were included in the multivariable Cox regression, regardless of whether the value was statistically significant in univariate analysis. Statistical analysis was performed by IBM SPSS version 22.0 (SPSS Inc., Chicago, IL, USA) and R software version 4.1.2 with the “rms,” “survival,” “survminer,” “ggsci,” and “forestplot” packages. was considered to be statistically significant.

3. Results

3.1. Perioperative Characteristics

From Jan 2013 to Dec 2016, a total of 3880 consecutive patients with malignant liver diseases from our single HPB center were retrospectively reviewed and validated. In these cases, patients with hepatocellular carcinoma (n = 2903), metastatic liver cancer (n = 265) and other malignancies (n = 43) were excluded. Among the patients with cholangiocarcinoma (n = 669), patients with perihilar CCA (n = 189) and distal CCA (n = 102) and patients undergoing TACE (n = 48) and microwave ablation (n = 54) were excluded. 276 ICC patients with surgical treatment were initially included in the analytic cohort, with 137 patients in the NAR group and 139 patients in the AR group. PSM analysis determined 99 pairs of patients for further survival and risk factor analysis (Figure 1).

Perioperative characteristics are summarized in Table 1. The preoperative RBC count was lower in the AR group (). There were more patients in the NAR group with Child–Pugh class A (91.97% vs. 84.17%), compared with those in the AR group (). There is no significant difference among the other baseline characteristics in two groups. The operation time in the AR group was longer (). Patients in the AR group had a larger amount of intraoperative hemorrhage () and the percentage of transfusion was higher (). More patients in AR group had singular lesion (), lymph node metastasis () and PVTT (). The hospital stay was shorter in NAR group (). There was no statistical significance in the other operative and postoperative parameters between the two groups.

PSM analysis selected 99 matched pairs from each group. There were no significant differences in baseline characteristics in AR and NAR groups after PSM. However, more patients in the AR group had transfusions (), and the hospital was longer in the AR group () as well. There were no significant differences in the incidence of complication (), types (), or Clavien–Dindo grade () between the two groups in the PSM cohort (Table 2).

3.2. Survival Analysis

Before PSM analysis, the 1-, 3-, and 5-year OS in the AR group (68%, 44%, and 28%, respectively) were comparable to those in the NAR group (64%, 34%, and 22%, respectively), though the Kaplan–Meier curves were suggestive but not significant (Figure 2(a)). There was no significant difference of 1-, 3-, and 5-year DFS between AR group (58%, 19%, and 7%, respectively) and NAR group (50%, 13%, and 3%, respectively) (Figure 2(c)). Patients showed better 1-, 3-, and 5-year OS in AR group (70%, 46%, and 34%, respectively) than NAR group (60%, 28%, and 16%, respectively) () (Figure 2(b)) after PSM analysis. In the PSM cohort, the 1-, 3-, and 5-year DFS in the AR group (61%, 21%, and 10%, respectively) were significantly better than those in the NAR group (49%, 12%, and 4%, respectively) () (Figure 2(d)).

During the follow-up, 6 (6.06%) patients in the AR group and 1 (1.01%) patient in the NAR group did not report disease recurrence (Table 3). Recurrence pattern analysis showed there were no statistical differences in intrahepatic recurrence, extrahepatic recurrence, and both intra-/extrahepatic recurrence. Patients in the AR group tended to have an intrahepatic recurrence in distant segments (39.39% vs. 31.31%) and less possibility of recurrence in resection margin (2.02% vs. 5.05%) or adjacent segments (26.26% vs. 32.32%), though the difference was not statistically significant ().

Subgroup analysis of OS and DFS after PSM analysis showed that the surgical strategy of AR in selected patients had an associated improved prognosis (Figure 3). Specific subgroups of patients, age >60 years, male, without HBV infection, CA19-9 > 37 U/mL, without laparoscopic approach, without portal vein occlusion, with IVC occlusion, singular tumor, well/moderate differentiation, without microvascular invasion, regardless of the tumor size, may benefit more from AR in OS. Other subgroups of patients, age >60 years, CA19-9 > 37 U/mL, Child–Pugh class A, with cirrhosis, without laparoscopic approach, without portal occlusion, tumor size >5 cm, singular tumor, well/moderate differentiation, and without microvascular invasion, may benefit more from AR in DFS (Figure 4).

3.3. Risk Factor Analysis

In the PSM cohort, the univariate analysis identified NAR, tumor size >5 cm, multiple tumors, poor differentiation and lymph node metastasis as significant risk factors for OS () (Table 4). NAR, without laparoscopic approach, tumor size >5 cm, multiple tumors, poor differentiation and lymph node metastasis were considered significant risk factors for DFS () (Table 5). Multivariate analysis showed NAR, tumor size >5 cm, multiple tumors and poor differentiation were independent risk factors for OS (), NAR, tumor size >5 cm, and poor differentiation were independent risk factors for DFS () after PSM analysis.

3.4. Nomogram Prediction

A nomogram model predicting OS of patients with ICC undergoing hepatectomy is shown in Figure 5. The prognostic nomogram was developed based on the following six prognostic factors: age, sex (male or female), group (AR or NAR), tumor number (singular or multiple), tumor size (≤5 cm or >5 cm), and differentiation (well/moderate or poor). Each factor was ascribed a weighted point total that implied a survival prognosis.

4. Discussion

ICC is the second most common primary hepatic malignant tumor, with radical liver resection as the only curative option [29]. Extended liver resection and vascular reconstruction, together with systemic therapy and locoregional treatments, enabled increasing rates of resection and improved OS in selected ICC patients [30]. AR and NAR have been debated in HCC for decades with still controversial results, while clinical studies focusing on surgical strategy of AR or NAR for ICC are rare. In this study, we found that AR improved the 1-, 3-, and 5-year OS (70%, 46%, and 34%, respectively) and DFS (61%, 21%, and 10%, respectively) for ICC patients with statistical significance after PSM analysis. AR presented comparable complications and recurrence when compared with NAR. Multivariate analysis identified NAR, tumor size >5 cm, multiple tumors, and poor differentiation as independent risk factors for OS (). Selected patients can benefit most from AR, according to subgroup analysis. A nomogram based on independent risk factors for OS and factors with clinical significance was constructed to predict OS in ICC patients.

Primary liver cancer is the fifth most commonly diagnosed malignancy worldwide, with a high prevalence in Asia and Africa. HCC arising from hepatocytes and ICC from bile duct epithelium are major types of primary liver cancer [31, 32]. ICC is located in the second-order bile ducts in the hepatic parenchyma, extinguished with pCCA in left and right common hepatic ducts and dCCA in common bile duct [33]. ICC can be classified as mass-forming, periductal infiltrating, and intraductal growth types by morphology [34, 35]. Chronic biliary tract inflammation owing to choledocholithiasis, cholelithiasis, primary sclerosing cholangitis, or liver fluke infection is associated with CCA, while patients are usually asymptomatic and have no underlying liver diseases [36]. Patients with HCC usually have underlying diseases like HBV/HCV infection, steatohepatitis or cirrhosis, and inclined to metastasis with blood flow, while ICC characterized itself with jaundice caused by biliary obstruction and lymph node metastasis. HBV/HCV infection may also be involved in the carcinogenesis of ICC [37, 38]. Interestingly, HBV-associated ICC has been reported to have a favorable prognosis, probably due to early diagnosis [39, 40]. Capecitabine is now first-line adjuvant therapy after curative intent resection [29]. Consecutive therapy plans based on MDT discussion are needed for this aggressive cancer with distinct anatomic, molecular, and clinical characteristics [41, 42].

Radical surgical treatment is still the only therapy with curative potential for ICC. An aggressive surgical approach, including major liver resection, has been recommended in many centers to improve outcomes. Our study found that patients with ICC benefited from AR in OS and DFS after PSM, indicating complete removal of tumor-bearing segments plays a significant role in improving the survival outcomes. Shen reported better survival outcomes were associated with AR in ICC patients with stage IB or II tumors without vascular invasion [18]. The 1-, 3-, and 5-year OS were 72.9%, 45.7%, and 36.0% in the AR group and 62.0%, 30.8%, and 25.3% in the NAR group after PSM. However, Yang concluded that NAR was not inferior to AR in survival outcomes for primary solitary ICC without direct invasion of contiguous organs or extrahepatic metastasis and potential benefits exist in NAR [19]. In their study, the NAR group had a more positive surgical margin, but the surgical margin had no significant impact on OS or DFS before and after PSM analysis.

Resection margin status and length are supposed to be associated with the incidence of local recurrence in theory, which is a technical concern during ICC surgery [11]. Technically, AR obtains a larger distance between surgical margin and tumor lesion, which had a higher potential for negative resection margin than NAR. Previous studies found a residual tumor on the surgical margin is likely to grow and spread much more aggressively, leading to early recurrence and dismal survival [43, 44]. However, the prognostic value of a wide margin remains controversial [10]. Resection margin ≥1 cm was associated with improved survival. Intrahepatic recurrence is inclined to happen in resection margin or adjacent segments in the NAR group but distant segments in the AR group. The result was suggestive based on proportion, although there was no significant difference. Postoperative complications were comparable (Clavien–Dindo grade I/II) in two groups, indicating both AR and NAR are technically safe in hepatectomy for ICC. Though surgical treatments improve the survival of ICC patients, systematic surveillance among patients with high-risk factors is necessary to avoid a late diagnosis of ICC in intermediate/advanced stages [33].

Lymph node (LN) metastases are recognized as an extremely poor prognostic risk factor no matter whether curative resection is applied [34]. The essential of surgery in patients with LN metastases detected preoperatively and routine lymphadenectomy in ICC surgery remain controversial [45]. The 8thAmerican Joint Committee on Cancer (AJCC) staging system recommends six nodes need to be analyzed. Combined with resection margin and perineural invasion, lymph node ratio (LNR) 15 was reported to be an independent predictor of DFS, OS, early, local, and distal recurrence [46]. In our research, LN metastases were a relative risk factor for OS in univariate analysis, but it was not an independent risk factor in multivariate analysis. However, lymphadenectomy is still recommended in ICC patients undergoing hepatectomy by experienced surgeons with/without preoperation evidence of LN metastasis [22, 47]. Accurate staging determined by routine lymphadenectomy is crucial for predicting prognosis and providing options for following treatments [48]. In addition, complete removal of lymph nodes with metastasis potential helps to reduce rates of regional recurrence, as well as jaundice or pyloric obstruction induced by recurrence.

The diagnosis of cholangiocarcinoma is accomplished by the combination of clinical/biochemical features and imaging findings with nonspecific tumor markers for suggestive complementary [49]. There is a lack of prognostic and predictive tumor markers for ICC, particularly compared with the clinical significance indicated by elevated AFP level and prognostic efficiency of PIVKA-II in HCC [50, 51]. The combination of CEA, CA125, and CA19-9 had been reported to have diagnostic effects [52, 53]. In our research, about half of the patients presented elevated CA19-9 before or after PSM analysis, while only a very small proportion of patients presented abnormal serum CEA and AFP. In subgroup analysis, AR was superior to NAR in patients with elevated CA19-9, which shows low sensitivity in early stages but increased sensitivity in advanced diseases [2, 54]. Besides traditional serum biomarkers, biomarkers from extracellular vesicles, metabolites, and nucleic acids, as well as next-generation biomarkers detected by high-throughput omics-based approaches, have the clinical application potential [49, 55].

Although surgical resection is the only curative treatment for patients with ICC, most patients are ineligible for surgery treatment at the time of the first diagnosis due to metastasis or local advancement [36]. Early detection of ICC by screening is of vital importance in potential population with risk factors. Tumor size, R0 resection, lymph node metastasis, differentiation, adjuvant chemotherapy, CA19-9, T stage, PVTT, HBV infection/vaccination, and Eastern Cooperative Oncology Group performance status (ECOG-PS) have been reported as prognostic factors influencing survival outcomes of ICC in recent researches [2, 31, 37, 40, 5658]. We identified NAR, tumor size >5 cm, multiple tumors, and poor differentiation as independent risk factors for OS in ICC patients undergoing hepatectomy. The nomogram in our study was developed from independent risk factors and factors with clinical significance. The AJCC TNM system is the most commonly used staging system for ICC, while other systems have been proposed, including the staging system for mass-forming type by Okabayashi and the LCSGJ system by Wang. [59]. Many researchers have attempted to establish nomogram models to be used as alternative standards in staging ICC subgroups [5961].

The minimally invasive approach for curative surgery of ICC has not been well established [62]. Laparoscopy is associated with less intraoperative blood loss, faster recovery, and fewer complications for most surgeries in general, but the benefits of laparoscopic resection for ICC are controversial due to the challenges in major hepatectomy, vascular and biliary reconstruction, and extended lymphadenectomy. Complete resection (R0) with adequate remnant liver parenchyma is the aim of resection regardless of surgery types, which should be technically feasible for both open and minimally invasive approaches. We enrolled 19 patients and 17 patients in the PSM cohort undergoing NAR or AR with a laparoscopic approach, respectively. Without laparoscopic approach was considered a risk factor for DFS in univariate analysis, though multivariate analysis showed it cannot predict the prognosis independently. Subgroup analysis showed AR and NAR were comparable in the laparoscopic group, while AR achieved better OS in patients with open surgery, which might be the result of patient selection. Surgeons preferred to choose laparoscopy on patients with singular, relatively smaller tumors, and better general situation to ensure the operation safety, the outcomes of whom tended to be impacted more by tumor characteristics rather than surgical intervention. However, the advantages of AR can be observed in patients with open surgery. Robotic resection has not been further analyzed due to the limited cases in the entire cohort (2.92% and 2.28%) and in the PSM cohort (3.03% and 3.04%) for NAR and AR groups.

We acknowledge the potential limitations in this study of selection bias and sample size. Although PSM analysis was applied to reduce the selection bias in a new cohort with comparable baseline characteristics, the possibility of other unconsidered biases remains in a retrospective study. The entire cohort was selected from 3880 consecutive patients from a single center in the past 4 years according to the inclusion and exclusion criteria; however, further randomized clinical trials (RCTs) among multiple centers with large surgery volumes are still demanded.

In conclusion, AR improved the long-term survival of ICC with comparable postoperative complications and similar recurrence patterns. Multivariate analysis showed NAR, tumor size >5 cm, multiple tumors, and poor differentiation were independent risk factors for OS. AR is suggested in ICC patients with sufficient remnant liver volume.

Data Availability

The data for the current study are available from the corresponding authors on reasonable request.

Ethical Approval

All procedures followed were in accordance with the ethical standards of the Medical Ethics Committee of Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology, and with the Helsinki Declaration of 1975, as revised in 2000.

Informed consent was obtained from all patients included in the study.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this article.

Authors’ Contributions

Chao Wang was responsible for the investigation, methodology, funding acquisition, software, and writing the original draft. Pingcuo Ciren contributed to the data curation and investigation. Awang Danzeng was responsible for data curation and investigation. Cheng-Long Zeng contributed to data curation and investigation. Zhi-Wei Zhang, Zhi-Yong Huang, Wan-Guang Zhang, Yi-Fa Chen, and Bi-Xiang Zhang were responsible for the investigation, acquiring resources, and methodology. Bin-Hao Zhang and Xiao-Ping Chen contributed to the conceptualization, validation, formal analysis, funding acquisition, supervision, and reviewing and editing of the manuscript.

Acknowledgments

The authors thank Dr. Long Zhang from Ganzhou People’s Hospital of Jiangxi Province (Ganzhou Hospital Affiliated to Nanchang University), Dr. Jia-yu Shi from Tongji Tianyou Hospital, Dr. Yuan Zhu, and Dr. Chen-Liang Liu from the China University of Geosciences for their assistance with software and statistical analysis. This research was supported by the Chen Xiao-ping Foundation for the Development of Science and Technology of Hubei Province (Grant no. CXPJJH121001-2021004) to Dr. Bin-hao Zhang and Chen Xiao-ping Foundation for the Development of Science and Technology of Hubei Province (Grant no. CXPJJH12000002-2020058) to Dr. Chao Wang.