{'fraction_non_alphanumeric': 0.035262723171738994, 'fraction_numerical': 0.017853739119384968, 'mean_word_length': 4.996190476190476, 'pattern_counts': {'":': 0, '<': 7, '<?xml version=': 0, '>': 2, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 1, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 1, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Nomogram for predicting occurrence of synchronous liver metastasis in colorectal cancer: a single-center retrospective study based on pathological factors

Purpose The purpose of this study was to explore the risk factors for synchronous liver metastasis (LM) of colorectal cancer (CRC) and to construct a nomogram for predicting the occurrence of synchronous LM based on baseline and pathological information. Methods The baseline and pathological information of 3190 CRC patients were enrolled in the study from the Department of Colorectal Surgery, the Second Affiliated Hospital of Harbin Medical University between 2012 and 2020. All patients were divided into development and validation cohorts with the 1:1 ratio. The characters of LM and none-LM patients in newly diagnosed colorectal cancer were utilized to explore the risk factors for synchronous LM with the univariate and multivariate logistic regression analyses. A predictive nomogram was constructed by using an R tool. In addition, receiver operating characteristic (ROC) curves was calculated to describe the discriminability of the nomogram. A calibration curve was plotted to compare the predicted and observed results of the nomogram. Decision-making curve analysis (DCA) was used to evaluate the clinical effect of nomogram. Results The nomogram consisted of six features including tumor site, vascular invasion (VI), T stage, N stage, preoperative CEA, and CA-199 level. ROC curves for the LM nomogram indicated good discrimination in the development (AUC = 0.885, 95% CI 0.854–0.916) and validation cohort (AUC = 0.857, 95% CI 0.821–0.893). The calibration curve showed that the prediction results of the nomogram were in good agreement with the actual observation results. Moreover, the DCA curves determined the clinical application value of predictive nomogram. Conclusions The pathologic-based nomogram could help clinicians to predict the occurrence of synchronous LM in postoperative CRC patients and provide a reference to perform appropriate metastatic screening plans and rational therapeutic options for the special population.


Introduction
Colorectal cancer (CRC) is the third most common malignancy in the world, with high incidence and mortality globally. Metastases is one of the most common contributors for death in CRC, of which liver metastasis (LM) is the most fatal [1]. At present, research on LM in Open Access *Correspondence: guiywang@163.com † Yunxiao Liu and Yuliuming Wang contributed equally to this work and should be considered as co-first authors. Department of Colorectal Surgery, the Second Affiliated Hospital of Harbin Medical University, Harbin, China CRC patients is continuing [2][3][4]. It is important to note that LM are found in more than 25% of patients with CRC at the time of their first diagnosis and occur in up to 25% of patients following removal of the primary tumor [5,6]. The 5-year survival rate for patients with CRC is approximately 56%, but can be significantly reduced when distant metastases are detected [7,8].
LM of CRC can be divided into synchronous and metachronous LM [9]. Synchronous LM was defined as LM detected at the time of diagnosis of the primary tumor or within 6 months after diagnosis. Jennie et al found that up to 18.3% of patients developed LM after radical CRC surgery. Therefore, early detection of highrisk postoperative populations for synchronous LM can help physicians to improve survival by targeting screening and individualizing treatment. With the development of scientific research capacity, more and more risk factors affecting LM of CRC have been discovered such as T stage, N stage, tumor site, preoperative CEA level, and so on [10,11]. However, most studies did not include vascular invasion (VI) as a pathological factor. Both the Association of Directors of Anatomic and Surgical Pathology [12] and the College of American Pathologists [13] emphasis is placed on documenting VI during routine pathological examination of cancer specimens. These institutions emphasized that VI was an independent predictor of poor prognosis and increased risk of LM because VI necessarily increased the risk of tumor cell entry into the bloodstream. Many studies have reached similar conclusions [14,15]. In addition, VI is an indication for adjuvant chemotherapy in stage II patients. Therefore, there are sufficient and necessary reasons to consider VI as a risk predictor for synchronous LM. Nomogram is a kind of graphical prediction model with friendly interface, which has strong clinical application value. By assigning points, we can not only observe the influence of certain parameters but also predict the probability of a particular event from the total score. To our knowledge, this is the first nomogram including VI to predict synchronous LM of CRC.

Patients
Three thousand one hundred ninety CRC cases were collected from the department of colorectal surgery, the Second Affiliated Hospital of Harbin Medical University between 2012 and 2020.
Inclusion criteria included (1) patients diagnosed with CRC and underwent surgery; (2) aged ≥ 18 years old; (3) Patients with complete pathological information; (4) CRC was the only primary malignancy. Exclusion criteria included (1) patients received preoperative neoadjuvant chemoradiotherapy; (2) the baseline and pathological information of the patient was incomplete; (3) distant metastases other than LM, and (4) postoperative resection margin of tumor was positive.
The colorectal cancer patients included in this study were treated with surgery only before the diagnosis of liver metastasis.

Variables
According to our study, age was regrouped into < 60, 60-74, and ≥ 75 years old; sex was classified as male and female; BMI was recorded as < 25 and ≥ 25; tumor size was divided into two groups: ≤ 5 cm and > 5 cm. The tumor site was grouped into right-sided colon (cecum, ascending colon, hepatic flexure, and transverse colon) and left-sided colon (splenic flexure, descending colon and sigmoid colon). The histology variable was classified as "adenocarcinoma, " "mucinous adenocarcinoma, " or "others"; the grade variable was classified as I/II and III/ IV ("well-differentiated/moderately differentiated, " and "poorly differentiated/undifferentiated") stage and the tumor type variable was classified as "ulcer type, " "uplift type, " and "infiltrating type. " Similarly, T stages, N stages, and lymph nodes examined (LNE) are grouped. The most prominent variable that VI, nerve invasion, and lymphatic invasion were classified as "yes" or "no. " The preoperative CEA level variable was classified as "positive" (≥ 5 ng/ml) and "negative" (< 5 ng/ml). In the same way, CA199 level variable was classified as "positive" (≥ 37 U/ ml) and "negative" (< 37 U/ml).

Statistical analysis
All statistical analyses were performed by R tool and SPSS 22.0. In this study, all patients were randomly (1:1 ratio) divided into development and validation cohorts and summarized by number and percentage. The characters of LM and none-LM patients in newly diagnosed colorectal cancer were utilized to explore the risk factors for synchronous LM with the univariate and multivariate logistic regression analyses. Nomogram was constructed based on the results of multivariate regression model and its performance was further evaluated by calibration and ROC. The DCA curve was used to evaluate the clinical decision ability of the model. In addition, all variables were assigned and the best cutoff value was calculated based on the total score by Youden's index. The difference was considered statistically significant for a two-sided P < 0.05.

Patients' characteristics
A total of 3190 patients were divided evenly into the development and validation cohorts. There were 104(6.5%) LM patients in development cohort and 107(6.7%) LM patients in validation cohort. In all patients, the majority of patients were men, aged 60-74 and had a BMI of less than 25. Overall, the main proportions of the patients were associated with rectum, tumor size ≤ 5 cm, ulcer type, grade I/II, adenocarcinoma, preoperative CEA and CA-199 level negative, T4 stage, N0 stage, and LNE ≥ 12. In the encroachment around, VI accounted for 28.8% and 30.5% in the development and validation cohorts. The detailed data was summarized in Table 1.
Using the nomogram derived scores, all LM patients were classified into two subgroup low-risk (risk score ≤ 193) and high-risk groups (risk score > 193) (Fig. 5). And we found there were significant differences in the occurrence of LM between the high-and low-risk groups in development cohort (P < 0.001) and validation cohort (P < 0.001) (Fig. 6).

Discussion
CRC is one of the major causes of cancer morbidity and mortality both in the worldwide. Distant metastasis disease is the main cause of poor prognosis in patients with CRC and the liver is the most common organ for metastasis. Over the course of the disease, a total of 50% of patients with CRC are likely to develop LM [16]. And the death rate remains high [17]. This also explains the necessity of this study. Nowadays, many studies have identified many independent factors for LM in patients with CRC but few studies developed predictive nomogram and the factors included are incomplete, this also has led to questions about the accuracy of those nomograms [18,19]. None of these studies had complete pathological information, such as VI, nerve invasion and lymphatic invasion, although nerve invasion and lymphatic invasion were not an independent risk factor in this study. Our study included not only the baseline information of the patients, but also the postoperative pathological information of the patients such as the VI. This greatly increases the sensitivity and specificity of this model. This nomogram can help identify high-risk groups for synchronous LM after surgery and help clinicians conduct targeted screening and individualized treatment. According to previous studies, the prevalence of CRC LM is more common than other metastases, such as the brain or bone [20,21] and the prevalence of LM was more than 20% [22,23]. However, the prevalence in this study was less than 20% which was consistent with one study conducted by Manfredi and his colleagues in France [24]. In this study, the incidence of elderly patients was significantly lower than that of young patients, with an incidence of only 13.2%. Because in the past few years, as colonoscopy screening has become more widespread, more and more early signs of cancer have been detected, leading to higher rates of disease in younger people and reducing the prevalence of LM. Tumor size has been shown to be an important factor in the development of LM from CRC [11], while little effect of tumor size was found in this study which might be partly attributed to sample size. Therefore, factors contributing to the development of LM after CRC should be identified and a screening method developed to determine the risk of LM in patients with postoperative CRC.
This study showed that tumor site, T stage, N stage, VI, preoperative CEA level, and CA-199 level were significantly associated with LM development, which were seldom reported before. In our nomogram, colon cancer (CC) is found to be more likely to metastasize than rectal cancer. The most common mechanism is that the metastatic pattern is different. In CC, most mesenteric drainage enters the hepatic portal vein system and the rectal venous-collected blood flows into the systemic circulation [25]. And right-sited CC is found to be more likely to metastasize than left-sited CC, which is consistent with previous articles [11,23,24], it was reported that the discrepancy was caused by molecular biological differences [26,27]. In addition, it is well known that high T stage, N stage, and tumor marker levels represent the malignancy of the tumor, indicating a high likelihood of LM. To some extent, it may explain the difference in prognosis, but this complex issue needs further research.
Since VI is closely related to progression and recurrence of the disease, both the Association of Directors of Anatomic and Surgical Pathology and the College of American Pathologists emphasis is placed on documenting VI during routine pathological examination of cancer specimens. It is reported that the pathologic detection rate of VI was 23% [28]. The proliferation of tumor cells requires energy and nutrition, tumor and other cells secrete VEGF, angiopoietin-like protein, and corresponding inflammatory cells to promote the formation of new blood vessels, thus providing an important channel for these functions [29]. Invasion of blood vessels by tumor cells is a key step in LM. Fujii T and Xie W et al. found that VI is closely related to the depth of tumor invasion and degree of differentiation [30,31]. A previous study has shown that VI is significantly associated with metastasis and high recurrence rates [32]. Furthermore, VI has been shown to be associated with reduced overall survival (OS) and disease-free survival (DFS), serving as a strong prognostic indicator [33]. These studies indicate that VI is an essential factor in LM of CRC. The inclusion of this factor improves the accuracy and persuasiveness of our model.    In addition, our research has certain research significance and advantages. Our model mainly identified patients at high risk of developing LM after surgery. For example, a right-sited (40 points) CC patient with VI (32 points), T4 stage (100 points), N2 (62 points), CEA level is positive (58 points), and CA199 level is negative (0 points) has a total of 292 points (high-risk group), resulting the diagnostic possibility is 0.38. In high-risk groups, where there may be micrometastases that were not detected before surgery, the enhanced computed tomography and monitoring should be justified after surgery. In addition, for high-risk patients, positron emission tomography computed tomography (PET-CT) and invasive procedures such as needle biopsy are also recommended for definite diagnosis if there is an undefined low-density liver lesion, because they are expensive and are not justified for general screening. If a positive liver result is found in a targeted examination, we can change the treatment plan in time so as not to delay the disease such as radiofrequency ablation and surgery [34][35][36]. Of course, this study only provides a reference rather than a guide, and specific decisions should be made based on clinical practice.
The study has several limitations. First of all, since this study is a retrospective study, involving 3190 patients from 2012 to 2020, there may be a possibility of inaccuracy due to the small amount of data, and inevitably there is observer and confusion bias. The current results require further validation in prospective clinical studies. Second, although our model has internal validation, there is a lack of external validation to further determine the model's accuracy. Third, some underlying factors such as gene status are unknown [37]. The inclusion of these important factors can further improve the effectiveness of the nomogram.
In conclusion, this study performed prediction of synchronous LM in patients undergoing CRC surgery. Therefore, if patients with preoperative diagnosis of colorectal cancer, no LM and underwent surgery are identified as a high-risk patient by nomogram, we should enhance post-operative imaging of the liver, such as enhanced computed tomography, MRI and positron emission tomography computed tomography (PET-CT). It can help clinicians timely detect the  disease progression of patients and take effective interventions to improve the quality of life of patients. This is crucial.