New Insights Into a Classification-Based Microvascular Invasion Prediction Model in Hepatocellular Carcinoma: A Multicenter Study

Background and Aims Most microvascular invasion (MVI)-predicting models have not considered MVI classification, and thus do not reflect true MVI effects on prognosis of patients with hepatocellular carcinoma (HCC). We aimed to develop a novel MVI-predicting model focused on MVI classification, hoping to provide useful information for clinical treatment strategy decision-making. Methods A retrospective study was conducted with data from two Chinese medical centers for 800 consecutive patients with HCC (derivation cohort) and 250 matched patients (external validation cohort). MVI-associated variables were identified by ordinal logistic regression. Predictive models were constructed based on multivariate analysis results and validated internally and externally. The models’ discriminative ability and calibration ability were examined. Results Four factors associated independently with MVI: tumor diameter, tumor number, serum lactate dehydrogenase (LDH) ≥ 176.58 U/L, and γ-glutamyl transpeptidase (γ-GGT). Area under the curve (AUC)s for our M2, M1, and M0 nomograms were 0.864, 0.648, and 0.782. Internal validation of all three models was confirmed with AUC analyses in D-sets (development datasets) and V-sets (validation datasets) and C-indices for each cohort. GiViTI calibration belt plots and Hosmer-Lemeshow (HL) chi-squared calibration values demonstrated good consistency between observed frequencies and predicted probabilities for the M2 and M0 nomograms. Although the M1 nomogram was well calibrated, its discrimination was poor. Conclusion We developed and validated MVI prediction models in patients with HCC that differentiate MVI classification and may provide useful guidance for treatment planning.


INTRODUCTION
Primary liver cancer is the sixth most commonly diagnosed cancer and the third leading cause of cancer death worldwide, with hepatocellular carcinoma (HCC) comprising 75~85% of primary liver cancer cases (1). Chronic hepatitis B virus (HBV) infection has been identified as the predominant determinant of HCC risk in China (1,2). Hepatectomy and liver transplantation are first-line potentially curative therapies for HCC (3). Although there have been advances in HCC management over the past decade, overall survival has not improved significantly, and longterm prognoses remain poor for most patients (4).
The concept of microvascular invasion (MVI) describes the histologic presence of malignant cells in tumor-adjacent microvessels. MVI is thought to capture local and distant metastasis processes and has been shown to be a significant risk factor for both early recurrence and poor prognosis after tumor resection or liver transplantation (5). Thus, MVI diagnosis could, theoretically, affect therapeutic planning, including decisions regarding the width of margins to be taken during hepatectomy, microwave/radiofrequency ablation, and postoperative adjunct treatment modalities. However, currently, MVI can only be confirmed definitively by postoperative pathology, limiting the role of preoperative evaluation in surgical planning and patient prognosis (6,7).
Although the notion that MVI affects prognosis adversely has gained acceptance over time, there is not a consensus regarding the precise relationship between MVI and patient outcomes (8,9). Our previous findings indicated that a simple distinction of MVI presence or absence, as has been most commonly done in MVI-related prediction modeling studies (6,7), is inadequate for predicting early recurrence after curative hepatectomy for HCC (10). Some recent prediction models in the literature have classified MVI in some detail (11)(12)(13)(14)(15)(16)(17), but there is not yet a universally recognized consensus MVI classification scheme.
The purpose of this study was to identify clinical variables that are significantly associated with particular classifications of MVI. We employed data from independent clinical cohorts of patients receiving hepatectomies in one of two Chinese Medical Centers to develop prediction models for predicting MVI classes. The models were validated externally.

Study Design and Participants
A retrospective study was conducted with data from two Chinese Medical Centers: Hunan Provincial People's Hospital (HPPH) and Affiliated Hospital of North Sichuan Medical College (NSMC). The institutional review boards of both medical centers approved the study (HPPH approval no: 2020-018; NSMC approval no: 2020-055). Given the purely observational nature of the study and no patients being contacted for the purpose of this study, the need for written informed consent was waived by the Institutional Review Board. Patient privacy was ensured, and the data were anonymized or maintained with confidentiality. All procedures were performed in accordance with the 1964 Helsinki Declaration and its later amendments or comparable ethical standards. The study was conducted in accordance with the Transparent Reporting of a Multivariable Prediction Model for Individual Prognosis or Diagnosis statement for reporting multivariable prediction model derivation and validation (18). Neither patients nor the public were involved in the design, conduct, reporting, or dissemination plans of this research.
A derivation cohort was formed from 800 eligible consecutive patients treated at Hunan Provincial People's Hospital with the inclusion criterion of a histological diagnosis of HCC following an initial hepatectomy performed between January 2015 and September 2019. An external validation cohort was formed from 250 consecutive matched patients operated on at the Affiliated Hospital of North Sichuan Medical College between January 2018 and September 2019 with otherwise the same inclusion criterion. The exclusionary criteria were: (a) a pathology finding other than HCC; (b) hepatectomy of recurrent HCC; (c) a prior hepatectomy for HCC at another hospital;(d) previous history of other type of cancer; and (e) other treatment (such as transcatheter arterial chemoembolization, radiotherapy) before hepatectomy. The clinical characteristics of the patients are summarized in Table 1.

Preoperative Assessment
Patients' demographic and clinicopathological data were collected according to routine post-admission practices. The latest preoperative laboratory and imaging examination results were used in the present analyses. Portal hypertension was diagnosed according to published criteria (19) and classified as mild, moderate, or severe, as described in detail previously (10). Portal vein tumor thrombosis (PVTT) and hepatic vein tumor thrombosis (HVTT) were graded based on published clinical criteria (20). Bile duct tumor thrombus (BDTT) was classified as previously reported (21).

Tumor Staging and Pathology
Resected tumors were staged postoperatively in accordance with the American Joint Committee on Cancer's TNM staging system (8 th edition) (22) and the Barcelona Clinic Liver Cancer system (23). Tumor size was defined as the maximum diameter of the pathology specimen. Tumor number was classified as 1, 2, 3, or ≥4. Satellite nodules were defined as microscopic HCC nodules separated from the tumor by at least 2 cm of uninvolved liver parenchyma and included in tumor counts. Tumor cell differentiation was classified according to the Edmondson-Steiner system (24). HCC pathology samples were taken in accordance with a 7-point baseline sample collection protocol (25). MVI was defined as the presence of a cancer cell nest in vessels lined by endothelial cells on microscopy and classified based on the 2015 update of the Evidence-Based Practice Guidelines for Standardized Pathological Diagnosis of Primary Liver Cancer in China (25) as follows: M0 (MVI absent); M1 (≤5 MVI sites, all located within adjacent peritumoral liver tissues ≤ 1 cm from the tumor margin); or M2 (>MVI sites, or the presence of any MVI site(s) within adjacent peritumoral liver tissues > 1 cm away from the tumor margin). MVIs in all tissue specimens were counted by examination of serial sections under a light microscopic. Liver tissue inflammation and fibrosis in non-tumor areas were graded according to the METAVIR scoring system (26). Two senior pathologists with more than 10 years of hepatic pathology experience performed pathological examinations of the surgical specimens. A consensus was reached by discussion for discordant cases.

Surgical Therapies
Surgery was performed by senior hepatic specialists with ≥ 15 years of clinical experience. HCCs were removed by open or laparoscopic hepatectomy, according to surgical methods as we reported previously (10). PVTT and/or HVTT were removed by tumor thrombectomy or by en-bloc-resection of the thrombotic vein(s), with or without post-resection portal/hepatic vein reconstruction. BDTT was addressed with tumor thrombectomy or by en-bloc-resection inclusive of the bile duct. An inferior vena cava tumor thrombosis was resected concomitantly with the primary tumor resection after establishing well-controlled total hepatic vascular exclusion.

Statistical Analysis
Categorical variables, which are reported as numbers of cases with percentages, were analyzed with chi-square or Fisher's exact tests. Continuous variables, which are reported as medians with minimum-maximum value ranges, were analyzed with Mann-Whitney tests. To assess the predictive values of MVI data, univariate ordinal logistic regression analyses were conducted with derivation cohort data. Continuous variables that were determined to be significant in the univariate analyses were included in follow-up statistical analyses. Otherwise, cut-off values were determined based on the maximum Youden index from receiver operating characteristic (ROC) analysis and then dichotomized into univariate analysis. Variables with p values < 0.05 in the univariate analysis were submitted to multivariate analysis in accordance with a stepwise backward-elimination procedure (threshold p < 0.05). An MVI prediction nomogram was developed based on the multivariate ordinal logistic regression analysis results.
The predictive performance of the presently constructed MVI prediction models was validated internally with derivation cohort data by 2000 bootstrap re-samplings of 75% of the original data and then validated externally. Discrimination ability was assessed by constructing an area under the ROC curve (AUROC) and estimating an associated C-statistic. AUROC values are reported with 95% confidence intervals (CIs). C-statistic values were calculated as described previously (27). Agreement between the MVI classification model predictions and empirically observed frequencies were assessed with a GiViTI calibration belt plot and the Hosmer-Lemeshow (HL) chi square test of goodness of fit. All statistical analyses were performed in SPSS software (version 23.0; IBM Corporation, Armonk, NY) and R software (version 3.3.0) with the rms package (version 5.1-1; http://www.R-project.org). In all cases, p < 0.05 was considered significant.

Patient and Clinicopathologic Characteristics
A study overview flow diagram is provided in Figure 1. The baseline characteristics of all included patients, including the derivation (N = 800) and external validation (N = 250) cohorts, are presented in Table 1. All other clinical characteristics were statistically similar between the two groups, with the exceptions of hepatitis virus infection type, tumor maximum diameter and METAVIR grade of background liver fibrosis.

Predictor Variables for MVI
The cut-off values of continuous variables determined by ROC curves are shown in Table 2. The results of univariate ordinal logistic regression analysis of clinical features in the derivation cohort are shown in Table 3. In the multivariate analysis of putative predictive factors found to be significant in the univariate analyses, risk factors that independently associated with the presence of MVI were tumor size, number of tumors, serum LDH ≥ 176.58 U/L, and serum g-GGT (U/L) ( Table 4).

Establishment and Validation of an M2 Predictive Nomogram
A predictive nomogram for M2 developed based on our multivariate analysis results is presented in Figure 2A. M2 prediction scores calculated from this nomogram consisted of the total points for each patient obtained by adding the points for each of the five analyzed factors. The AUC of the nomogram to predict M2 were 0.864 (95% CI, 0.8081-0.9196) ( Figure 2B). In bootstrap analysis for internal validation of the derivation cohort data, we obtained AUCs for M2 prediction of 0.865 (95% CI, 0.8038-0.9234) for the development dataset (D-set) and 0.854 (95% CI, 0.7541-0.9534) for the validation dataset (V-set) ( Figure 2C), with a c-index of 0.817 (95% CI, 0.774-0.860) indicating that the nomogram had good performance for distinguishing between M2 and non-M2 patients. A GiViTI calibration plot showed good consistency between the observed frequency and predicted probability for M2 among patients in the derivation cohort ( Figure 2D). The HL chi-square calibration value of the M2 model was 7.601 for the derivation cohort (p = 0.473).
In the bootstrap analysis for internal validation performed with the external validation cohort data, AUCs for M2 prediction

Establishment and Validation of an M1 Predictive Nomogram
A predictive nomogram was constructed for predicting M1 status based on the results of our multivariate analysis ( Figure 3A). The AUC of the predictive model for M1 diagnosis was 0.648 (95% CI, 0.5773-0.7197) ( Figure 3B). In bootstrap analysis for internal validation, the AUCs for M1 prediction in the D-set and V-set were 0.631 (95% CI, 0.5473-0.7149) and 0.602 (95% CI, 0.5412-0.7533), respectively ( Figure 3C). A GiViTI calibration plot showed good consistency between the observed frequency and predicted probability for M1 among patients in the derivation cohort ( Figure 3D), the c-index was 0.623 (95% CI, 0.551-0.696). The HL chi-squared calibration values for M1 prediction were 5.284 (derivation cohort, p = 0.727).
In the bootstrap analysis for internal validation performed with the external validation cohort data, AUCs for M1 prediction were 0.626 (95% CI, 0.5026-0.7109) for the D-set and 0.587 (95% CI, 0.4827-0.6854) for the V-set ( Figure 3E), with a c-index of 0.623 (95% CI, 0.551-0.696). A GiViTI calibration plot affirmed that the model calibration yielded good consistency between the observed frequency and predicted probability for among patients with M1-class MVI in the external validation cohort ( Figure 3F), with an HL chi-squared calibration M1 prediction value of 4.523 for the external validation cohort (p =0.807).
These results indicate that there was fair discrimination, and the risk of disease was relatively close across groups and not readily distinguishable among M1-class and non-M1 MVI class, although the model has satisfactory calibration. These data do not support the effectiveness of the M1 predictive nomogram in clinical applications.

Establishment and Validation of an M0 Predictive Nomogram
A predictive nomogram was constructed for predicting M0 status based on the results of our multivariate analysis ( Figure 4A). The AUC of the M0 predictive nomogram was 0.782 (95% CI, 0.7044-0.8597) ( Figure 4B). In bootstrap analysis for internal validation, the AUCs for M0 prediction in the D-set and V-set were 0.758 (95% CI, 0.6921-0.8269) and 0.761 (95% CI, 0.6348-0.8877) ( Figure 4C), respectively, the c-index was 0.751(95% CI, 0.706-0.797). A GiViTI calibration belt plot was generated and demonstrated model calibration with good consistency between the observed frequency and predicted probability of patients with an M0 status in the derivation cohort ( Figure 4D). The HL chi-squared calibration value for M0 prediction was 9.213 for the derivation cohort (p = 0.325). In bootstrap analysis for internal validation of the external validation cohort, the AUCs for M0 prediction in the D-set and Vset were 0.780 (95% CI, 0.6978-0.8428) and 0.755 (95%CI, 0.6688-0.8574) ( Figure 4E), respectively, with a c-index was 0.732 (95% CI, 0.685-0.779). A GiViTI calibration plot was generated and demonstrated that the model calibration showed good      consistency between the observed frequency and predicted probability of M0 status among patients in the external validation cohort ( Figure 4F). The HL chi-square calibration value of the M0 model was 5.690 for the external validation cohort (p = 0.682). The M0 predictive nomogram showed good performance in distinguishing patients with versus not with M0 status.

DISCUSSION
The present multicenter study provides a comprehensive review of clinical data from hepatectomized HCC patients treated at two Chinese hepatobiliary centers. This work addresses the necessity to differentiate MVI classification before making treatment strategy in this patient population. Extensive studies have been conducted to explore useful clinical/tumor/radiological factors for prediction of MVI status. Most have been single-center cohort studies with relatively small sample sizes, characteristics which may impede their reliability and limit the generalizability of their conclusions (5). It is noteworthy that Lei et al. published a large cohort study (N = 1004) of patients with HCC (diagnosed in accordance with the Milan criteria) (28); their cohort consisted mainly of patients with HBV-related HCC and had characteristics similar to ours.  In another recent relatively large two-center retrospective cohort study conducted in China, Zhang et al. (29) adopted MVI classification criteria that were the same as ours. In their study, they predicted MVI status based on a combination of computed tomography radiomics data and two clinical factors: age and alpha-fetoprotein test result (positive/negative). They obtained AUCs of 0.806, 0.803, and 0.796 in their training, test, and independent validation cohorts, respectively. Although they used MVI classification criteria delineating M0, M1, and M2, they did not provide data on prediction of M1 or M0. Based on AUCs for prediction, our results for M2 are likely comparable with Zhang et al.'s MVI risk prediction for MVI positivity/ negativity; C-index values were not reported in Zhang's study. Importantly, simple MVI positivity cannot be regarded as being indicative of M2 status, and there were no robust data differentiating M0 or M1 that could be inferred from Zhang's model. MVI classification has not received sufficient attention in prognosis evaluation for patients with HCC. In a recent multicenter study, researchers from the Liver Cancer Pathology Group of China (30) used MVI classification data to analyze postoperative prognosis in 2,573 HCC patients who underwent curative hepatectomy. They found that the 3-year postoperative recurrence rates in their M0, M1, and M2 groups were 62.5%, 71.6%, and 86.1%, respectively (p < 0.001), with corresponding 3-year overall survival rates of 94.1%, 87.5% and 67.0%, respectively (p < 0.001). M1 grade was associated with early recurrence, whereas M2 grade was associated with both early and late recurrence. Our current study expands upon the prior data with a focus on MVI classification and further addresses the importance of MVI classification in evaluating the potential importance of MVI in HCC patient prognosis.
The main strength of the present study is that we developed models for predicting M0, M1, and M2 disease classifications. We found that only the M0 and M2 models could be validated as clinically applicable for use in patients undergoing hepatectomy for HCC. Our nomogram models demonstrated good prediction abilities for M0 and M2 and yielded good consistency between predicted probabilities of MVI and empirical frequency of MVIs in both internal and external validation cohorts. According to the MVI classification criteria used in our study, the clinicopathological factors selected by multivariate analysis showed unsatisfactory discrimination for the M1 prediction model, despite the model showing good calibration. Also, very few clinicopathological variables harvested meaningful cutoff point values determined by ROC analysis for M1 diagnosis. Given the unsatisfactory performance of the M1 nomogram in the present study, and our previous findings showing no difference in early recurrence after curative hepatectomy between M0 and M1 class HCC patients (10), it may be that MVI only has measurable effects on prognosis (postoperative recurrence, overall survival) when it reaches M2 class severity. The present study addresses the question of whether the M1 class distinction provides any advantage. When MVI is classified simply as absent (M0) or present (M1/M2), patients with M1 disease might receive overtreatment, such as unnecessary postoperative adjunct therapies, and inaccurate prognostic information given that clinical demonstrations have shown that adjuvant transarterial chemoembolization can improve overall and disease-free survival in HCC patients with MVI (31). The significance of MVI class is also relevant for HCC management with respect to deciding whether a particular patient should have local ablation, non-anatomical surgical resection, or anatomical surgical resection, and appropriate margin widths. A notable strength of our study is that we included complex and diverse HCC cases, augmenting the reliability and applicability of the results. It seems reasonable to presume that the incidence of MVI, especially M2, should be positively related to an advanced state of HCC with more invasive features. However, our current findings do not support such a view. We speculate that HCC biological characteristics related to MVI appear to differ from those related to macrovascular invasion and thus may have a distinct influence on prognosis. The molecular biological mechanisms of such a distinct influence remain to be clarified. Interestingly, given the previous studies which would suggest otherwise (16,32), we were surprised that preoperative HBV-DNA level was not found to be significantly associated with MVI in our HCC patients.
To date, various MVI classification schema have been proposed (11)(12)(13)(14)(15)(16)(17) and no consensus on an MVI classification scheme has been reached. We favor the strategy of establishing an MVI classification schema that takes into account the severity of MVI burden, which may include the number of malignant cells in an embolus, the presence or absence of vessel wall invasion, and distance from the primary lesion. Tumor size and number have been identified as risk factors in previous reports, albeit with different cut-off values (6). The presently obtained ROC curve analysis-based tumor size cut-off values for M1 (5.55cm) and M2 (6.95 cm) were quite close. A previous meta-analysis found that MVI was detected in about 25% of small liver cancer cases (<2 cm) (33). In the present study, among cases with a solitary HCC with a diameter ≤ 2 cm, M1 and M2 classification rates for the derivation cohort were 15.8% (9/57) and 3.5% (2/57), while M1 and M2 classification rates for the external validation group were 9.1% (1/ 11) and 0 (0/11), respectively. In our multivariate analysis, tumor size was inputted as a continuous variable to minimize information loss, thereby reflecting the value of tumor size for MVI prediction. With respect to tumor number, we were only able to harvest a clinically useful cut-off value from our ROC analysis for M2 prediction (≥1.5, effectively ≥2).
In this study, we also considered the number of satellite neoplastic lesions, which may be more reliable for preoperative imaging evaluation. HCC satellite lesions are detached micrometastatic nodules of HCC cells that have become embedded within the hepatic parenchyma, differing from MVI lesions in that they are thought to represent direct tissue invasion by malignant cells. Hence, MVI prediction models based on radiology examination alone may be intrinsically limited.
Several findings that have not been reported in previous studies are presented in this research report. Historically, serum LDH has been an indicator of cell injury and necrosis, and thus may reflect tumor burden. In a meta-analysis, preoperative LDH was shown to be significantly associated with poor prognosis in patients with HCC (34). However, further analysis showed that serum LDH levels correlated inversely with checkpoint inhibitor treatment responsivity, suggesting that the prognostic attributes of elevated LDH levels may go beyond tumor burden alone (35). High LDH levels could be consequent to elevated glycolytic activity of tumor cells and hypoxia-associated tumor necrosis; both glycolysis and hypoxia contribute to the immuno-suppressiveness of microenvironments (36). Elevation of the glycolysis metabolite lactate due to increased serum LDH can influence the complex metabolic crosstalk between tumor cells and tumor-infiltrating immune cells in the tumor microenvironment in multiple ways (37,38). Additionally, lactate increases acidify the tumor microenvironment, thereby promoting multiple critical oncogenic processes, including angiogenesis, tumor cell invasion, and metastasis, while suppressing anticancer immunity (33,39). Regarding the present finding showing that an LDH serum level ≥ 176.58 U/L was independently associated with MVI, it is interesting to note that this cut-off value is within the normal range for LDH (100.0-240.0 U/L). Moreover, we did not assess different LDH isoforms separately because serum LDH isoform analysis is not regularly performed in this patient population. In addition, the relationship between serum LDH levels and tumor LDH expression remains unclear, and the mechanism by which either may influence MVI are not known. Notwithstanding, our findings indicate that LDH levels may play an important role in MVI prediction in patients with HCC. g-GGT alone or combination with other markers has been extensively explored as a prognostic factor in HCC patients (40)(41)(42)(43)(44). Here, we found that serum g-GGT level was a significant predictor of MVI, a finding that has not been reported previously. High serum levels of g-GGT prior to treatment have been reported to correlate with poor survival as well as with unfavorable clinicopathological characteristics in patients diagnosed with HCC (45). For a detailed discussion of the potential mechanisms by which serum g-GGT level may have prognostic value in patients with HCC, the reader is referred to Sun and colleagues' metaanalysis of 18 studies from China and 2 studies from western countries (45). Since Sun's review was published, subsequent g-GGT studies of HCC have been conducted primarily in HBVendemic areas (40,41,43,44). It would be of interest to determine whether g-GGT may be a more valuable marker in HBV-related HCC than in non-HBV-related HCC. We speculate that g-GGT level might reflect the formation of a microenvironment within which oxidative stress promotes tumor cells to produce g-GGT, contributing, to some extent, to MVI formation. However, the exact direct mechanisms of elevated g-GGT in cancer initiation and progression have not been resolved. The question of whether such effects of g-GGT level on MVI occur with HCC arising from other etiologies deserves further investigation. Prior studies that have examined preoperative radiological evaluation of MVI status have been criticized fairly for having relatively small sample sizes, using varying imaging modalities and radiological parameters, and using advanced imaging technologies with limited utility in clinical practice. These drawbacks have precluded consensus and the establishment of valuable information to guide preoperative accurate MVI prediction in clinical practice. Despite being far from consensus, imaging evaluation approaches remain promising for MVI prediction and worthy of further research (46). In a recent meta-analysis of pooled data (46), diagnostic AUC values for radiomics and non-radiomics factors were similar (0.8550 and 0.8601, respectively). Our proposed model for M2 prediction outperformed these radiomics or non-radiomics methods slightly in terms of the simple index of AUC (0.864).
A major strength of our study is the inclusion of prediction models for M0 and M1. This valuable information cannot be acquired from the studies included in Huang's meta-analysis (46) because included studies dichotomized MVI status as present or absent. Overlooking these factors could affect the judgment of MVI status based on preoperative imaging tools. We have noticed several shortcomings when radiological and/or radiomics findings alone are used for MVI prediction model construction. Firstly, tumor characteristics inferred from preoperative imaging, including tumor size and number, may deviate from pathological findings. Secondly, the accuracy and practicality of radiomics models are questionable due to the lack of standardization in radiomics. Depending on the subjective judgment of particular diagnostic radiologists, some parameters may be overly specialized.
With respect to the clinical relevance of the present findings, it is important to note that our models involve variables that are easily detectable. Moreover, the fact that a large dataset was gathered from more than one center is an added strength in that it makes our conclusions more generalizable. The medical data that we analyzed were collected under real-world clinical conditions, providing high external validity and high generalizability. Our nomograms allowed reliable prediction of M0 and M2 for patients undergoing hepatectomy, and thus can be used to support preoperative treatment planning for patients with HCC.

STUDY LIMITATIONS
Firstly, due to the retrospective nature of the study, potential selection bias might exist, causing an unbalanced number of patients in different MVI classifications. To overcome this limitation to some extent, cases were selected consecutively to reduce bias. Secondly, this study was conducted in a region where HBV is endemic. Therefore, the applicability of the results may be limited to a certain extent. Despite these limitations, this study provides notable insights that can help surgeons and clinicians determine appropriate treatment strategies for patients with HCC.

CONCLUSION
We have developed and validated an MVI prediction model for patients with HCC. The model demonstrates excellent discriminative ability for M0 and M2 classes of MVI diagnosis. The present findings provide guidance information that clinicians can use in practice to improve preoperative treatment strategy decisions.

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
The studies involving human participants were reviewed and approved by Hunan Provincial People's Hospital (HPPH) and Affiliated Hospital of North Sichuan Medical College (NSMC). (HPPH approval no: 2020-018; NSMC approval no: 2020-055). The ethics committee waived the requirement of written informed consent for participation.

AUTHOR CONTRIBUTIONS
Study concept and design: WX and JL. Acquisition of data: YW, ZY, RL, and FL. Analysis and interpretation of data: WX, YW, ZY, and JL. Drafting of the manuscript: WX, YW, ZY, and JL. Statistical analysis: ZY, RL, and FL. Study supervision: WX and JL. All authors contributed to the article and approved the submitted version.

FUNDING
This study was supported by funds from the Health Commission (project no. 20200074) and Education Department (project no. 20A313) of Hunan Province.