This study includes two steps. First, clinical and survival data of MPCC were obtained from the SEER database. Then, the R software (version 4.3.3) was used to analyze and construct a competing risk model. The workflow diagram of the study is presented in Figure 1. This study does not require ethical approval and consent to participate. The data used in this study was obtained from databases.
1. Data acquisition
<ol>
	<li>Download the SEER.Stat 8.4.3 software from the SEER database website (http://seer.cancer.gov/about/overview.html). Register and log in to SEER.Stat 8.4.3 to obtain relevant patient data.</li>
	<li>After logging in to SEER.Stat 8.4.3, click on Case Listing Session &#62; Data and select the Incidence SEER Research Data, 17 Registries, Nov 2022 Sub (2000-2020) database.</li>
	<li>Click on Selection &#62; Edit and choose {Race, Sex, Year Dx. Year of diagnosis} = &#39;2004&#39;, &#39;2005&#39;, &#39;2006&#39;, &#39;2007&#39;, &#39;2008&#39;, &#39;2009&#39;, &#39;2010&#39;, &#39;2011&#39;, &#39;2012&#39;, &#39;2013&#39;, &#39;2014&#39;, &#39;2015&#39; AND {Site and Morphology. Site recode ICD-O-3/WHO 2008} = &#39;Colon and Rectum&#39; AND {Site and Morphology. Diagnostic Confirmation} = &#39;Positive histology&#39; AND {Multiple Primary Fields. Sequence number}! = &#39;One primary only&#39;.</li>
	<li>Click on OK and save the selection. Click on Table, and in the Available variables interface, select Age recode with single ages and 85+, Sex, Site recode ICD-O-3/WHO 2008, CS tumor size, Grade Recode (thru 2017), Derived AJCC T, 6th ed (2004-2015), Derived AJCC N, 6th ed (2004-2015), Derived AJCC M, 6th ed (2004-2015), Radiation recode, Chemotherapy recode (yes, no/unk), SEER cause-specific death classification, SEER other cause of death classification, Survival months, Patient ID, RX Summ-Surg Prim Site (1998+), and click on Column.</li>
	<li>Click on Output, name the data, and click on Execute to output and save the data.</li>
	<li>The SEER database does not directly provide information on whether a patient has MPCC. After completing the data download as described above, use the patient ID to filter out patients with MPCC, that is, those diagnosed with two or more occurrences. After identifying MPCC patients, calculate the interval between tumor diagnoses based on their survival times. Using 6 months as the cutoff, classify patients into SMPCC and MMPCC.</li>
	<li>Based on their age at the time of initial diagnosis, categorize patients&#39; ages as = 65 years and&#62; 65 years. For tumor grade, grade I corresponds to good differentiation, grade II corresponds to moderate differentiation, grade III corresponds to low differentiation, and grade IV corresponds to undifferentiation. Classify tumor location based on the distribution of multiple tumors as right colon, left colorectum, and entire colorectum.
	<ol>
		<li>There are obvious differences in embryonic development and biological characteristics between the proximal and distal colorectum, divide the tumor location based on the splenic flexure23. The right colon is defined as proximal to the splenic flexure, including the cecum, ascending colon, hepatic flexure, and transverse colon, while the left colorectum includes the splenic flexure, descending colon, sigmoid colon, rectosigmoid junction, and rectum. If a patient&#39;s tumors are entirely located in the right colon, define the tumor location as right colon; if all tumors are located in the left colorectum, define the location as left colorectum; if there are tumors in both the right colon and left colorectum, define the location as entire colorectum.</li>
	</ol>
	</li>
	<li>Determine tumor size by selecting the largest tumor diameter among multiple tumors in the same patient. Based on the tumor diameter, categorize size as = 5 cm or &#62;5 cm. A total of 8,931 patients were included in the study. To create a training cohort and a validation cohort, randomly divide all cases in a 7:3 ratio.</li>
</ol>
2. Model construction
<ol>
	<li>Download RStudio (2023.12.1+402) and R software (4.3.3). Open RStudio to run R software. Click on New File and select R Script to create a new R programming interface. Enter the relevant code in the code editor and click on Run to execute the code. 
	NOTE: The R language and its functions offer a wide variety of parameters. Adding or modifying these parameters can enhance data analysis and visualization.</li>
	<li>Use the following code to perform univariate analysis and plot the CIF curve. 
	library(tidycmprsk) 
	library(gtsummary) 
	library(ggsurvfit) 
	library(ggprism) 
	data &#60;-read.csv(&#39;data.csv&#39;)&#34; 
	data$status&#60;-factor(data$status, levels=c(0,1,2),labels=c(&#34;0&#34;,&#34;CSS&#34;,&#34;OCS&#34;)) 
	CIF &#60;- tidycmprsk::cuminc(Surv(time, status) ~ Sex, data = data) 
	ggcuminc(CIF,outcome= c(&#34;CSS&#34;, &#34;OCS&#34;),size=1.5) 
	where data.csv is data obtained from the SEER database.</li>
	<li>After running the above code, click on Export, then click Save as Image, and finally click Save to save the image. The method for saving subsequent images will be the same as in this step. Replace Sex in the above code one by one with other factors to perform univariate analysis for all factors.</li>
	<li>Use the following code to perform BSR and multivariate analysis and visualization. 
	library(leaps) 
	library(riskRegression) 
	library(prodlim) 
	library(forestplot) 
	leaps&#60;-regsubsets(status==0~Sex+Size+Grade+T+N+M,data = data) 
	plot(leaps,scale=&#34;adjr2&#34;) 
	multi &#60;- FGR(Hist(time,status)~Sex + T + N + M +Grade+ Size,cause=1,data=data) 
	summary(multi) 
	Multi &#60;- read.csv(&#34;multi.csv&#34;,header = T) 
	forestplot(labeltext=as.matrix(Multi [,1:4]), mean= Multi $HR_mean, ower= Multi $HR_1, upper= Multi $HR_2) 
	where the data in multi.csv comes from the results of the previous code.</li>
	<li>Follow step 2.3. to save the image.</li>
	<li>Use the following code to plot the nomogram, ROC curve, calibration curve, and DCA curve. 
	library(QHScrnomo) 
	library(rms) 
	library(riskRegression) 
	library(prodlim) 
	data &#60;- read.csv(&#34;data.csv&#34;) 
	d &#60;- datadist(data) 
	options(datadist = &#34;d&#34;) 
	e &#60;-cph(Surv(time,status==1)~Sex + T + N + M +Grade,data = data, 
	x=T, y=T, surv=T,time.inc=60) 
	nomo &#60;- crr.fit(e,failcode=1,cencode = 0) 
	nomogram.crr(fit =nomo, lp = F, xfrac = 0.5, fun.at =seq(from=0, to=1, by= 0.1) , failtime =c(12,36,60),funlabel = c(&#34;1-year CSS Cumulative Incidence&#34;, &#34;3-year CSS Cumulative Incidence&#34;,&#34;5-year CSS Cumulative Incidence&#34;)) 
	set.seed(123) 
	data$pro &#60;- tenf.crr(m3,time = 60) 
	groupci(x=data$pro, ftime = data$time, fstatus = data$status, failcode = 1, cencode = 0, ci = TRUE) 
	f &#60;- CSC(Hist(time,status)~Sex + T + N + M +Grade,data = data) 
	x &#60;- Score(list(model1=f), Hist(time,status)~1, data=data, cause=1, times=c(12,36,60), se.fit=1L, plots=&#34;roc&#34;, metrics=&#34;auc&#34;) 
	g &#60;- as.data.frame(x$AUC$score) 
	h &#60;- g[g$times %in% c(12,36,60),] 
	col = c(&#34;darkcyan&#34;,&#34;tomato&#34;,&#34;purple&#34;) 
	plotROC(x, xlab=&#34;1-Specificity&#34;, ylab=&#34;Sensitivity&#34;,col=col[1], cex=1.5, legend=&#34;&#34;, auc.in.legend = F, times = 12) 
	plotROC(x,col=col[2],legend = &#39;&#39;, cex=1.5,times =36,auc.in.legend = F,add=T) 
	plotROC(x,col=col[3], times =60, add=T, cex=1.5, legend = &#39;&#39;,auc.in.legend = F) 
	leg &#60;- paste0(c(&#34;1-Year: &#34;,&#34;3-Year: &#34;,&#34;5-Year: &#34;),substr(a$AUC,1,5))</li>
	<li>Follow step 2.3. to save the image.</li>
	<li>Use the ggscidca package to plot the DCA curve. 
	library(ggscidca) 
	df_surv&#60;-read.csv(&#34;data.csv&#34;,header = T) 
	cox_model &#60;- coxph(Surv(time,status==1)~Sex + T + N+ M +Grade, data = df_surv) 
	cox_model1&#60;-newcrr(cox_model) 
	scidca(cox_model1, newdata= df_surv,threshold.line = T,threshold.text = T) 
	where data.csv is data obtained from the SEER database.</li>
	<li>Follow step 2.3. to save the image.</li>
</ol>

Nomogram Development for Survival Prediction in MPCC Patients

<ol>
	<li>Bray, F., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Global+cancer+statistics+2022:+GLOBOCAN+estimates+of+incidence+and+mortality+worldwide+for+36+cancers+in+185+countries.">Global cancer statistics 2022: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries.</a> CA Cancer J Clin. 74 (3), 229-263 (2024).</li><li>Derwinger, K., Gustavsson, B. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+study+of+aspects+on+gender+and+prognosis+in+synchronous+colorectal+cancer.">A study of aspects on gender and prognosis in synchronous colorectal cancer.</a> Clin Med Insights Oncol. 5, 259-264 (2011).</li><li>Chin, C. C., Kuo, Y. H., Chiang, J. M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Synchronous+colorectal+carcinoma:+predisposing+factors+and+characteristics.">Synchronous colorectal carcinoma: predisposing factors and characteristics.</a> Colorectal Dis. 21 (4), 432-440 (2019).</li><li>Nikoloudis, N., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Synchronous+colorectal+cancer.">Synchronous colorectal cancer.</a> Tech Coloproctol. 8 (Suppl 1), s177-s179 (2004).</li><li>Huang, C. S., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Synchronous+and+metachronous+colorectal+cancers:+Distinct+disease+entities+or+different+disease+courses.">Synchronous and metachronous colorectal cancers: Distinct disease entities or different disease courses.</a> Hepatogastroenterology. 62 (140), 838-842 (2015).</li><li>He, W., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Prognosis+of+synchronous+colorectal+carcinoma+compared+to+solitary+colorectal+carcinoma:+a+matched+pair+analysis.">Prognosis of synchronous colorectal carcinoma compared to solitary colorectal carcinoma: a matched pair analysis.</a> Eur J Gastroenterol Hepatol. 31 (12), 1489-1495 (2019).</li><li>Ranstam, J., Cook, J. A., Collins, G. S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Clinical+prediction+models.">Clinical prediction models.</a> Br J Surg. 103 (13), 1886 (2016).</li><li>Smith, T., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Comparison+of+prognostic+models+to+predict+the+occurrence+of+colorectal+cancer+in+asymptomatic+individuals:+a+systematic+literature+review+and+external+validation+in+the+EPIC+and+UK+Biobank+prospective+cohort+studies.">Comparison of prognostic models to predict the occurrence of colorectal cancer in asymptomatic individuals: a systematic literature review and external validation in the EPIC and UK Biobank prospective cohort studies.</a> Gut. 68 (4), 672-683 (2019).</li><li>Iasonos, A., Schrag, D., Raj, G. V., Panageas, K. S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=How+to+build+and+interpret+a+nomogram+for+cancer+prognosis.">How to build and interpret a nomogram for cancer prognosis.</a> J Clin Oncol. 26 (8), 1364-1370 (2008).</li><li>Conroy, T., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Five-year+outcomes+of+FOLFIRINOX+vs+Gemcitabine+as+adjuvant+therapy+for+pancreatic+cancer:+A+randomized+clinical+trial.">Five-year outcomes of FOLFIRINOX vs Gemcitabine as adjuvant therapy for pancreatic cancer: A randomized clinical trial.</a> JAMA Oncol. 8 (11), 1571-1578 (2022).</li><li>Zhong, M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Impact+of+lung+metastasis+versus+metastasis+of+bone,+brain,+or+liver+on+overall+survival+and+thyroid+cancer-specific+survival+of+thyroid+cancer+patients:+A+population-based+study.">Impact of lung metastasis versus metastasis of bone, brain, or liver on overall survival and thyroid cancer-specific survival of thyroid cancer patients: A population-based study.</a> Cancers. 14 (13), 3133 (2022).</li><li>Xu, Y., Wang, X., Huang, Y., Ye, D., Chi, P. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+LASSO-based+survival+prediction+model+for+patients+with+synchronous+colorectal+carcinomas+based+on+SEER.">A LASSO-based survival prediction model for patients with synchronous colorectal carcinomas based on SEER.</a> Transl Cancer Res. 11 (8), 2795-2809 (2022).</li><li>Zhang, X., Zhao, L., Hu, Y., Deng, K., Ren, W. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+novel+risk+prediction+nomogram+for+early+death+in+patients+with+resected+synchronous+multiple+primary+colorectal+cancer+based+on+the+SEER+database.">A novel risk prediction nomogram for early death in patients with resected synchronous multiple primary colorectal cancer based on the SEER database.</a> Int J Colorectal Dis. 38 (1), 130 (2023).</li><li>Zhang, X., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Developing+prognostic+nomograms+to+predict+overall+survival+and+cancer-specific+survival+in+synchronous+multiple+primary+colorectal+cancer+based+on+the+SEER+database.">Developing prognostic nomograms to predict overall survival and cancer-specific survival in synchronous multiple primary colorectal cancer based on the SEER database.</a> J Cancer Res Clin Oncol. 149 (15), 14057-14070 (2023).</li><li>Austin, P. C., Lee, D. S., Fine, J. P. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Introduction+to+the+analysis+of+survival+data+in+the+presence+of+competing+risks.">Introduction to the analysis of survival data in the presence of competing risks.</a> Circulation. 133 (6), 601-609 (2016).</li><li>Lau, B., Cole, S. R., Gange, S. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Competing+risk+regression+models+for+epidemiologic+data.">Competing risk regression models for epidemiologic data.</a> Am J Epidemiol. 170 (2), 244-256 (2009).</li><li>Nolan, E. K., Chen, H. Y. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+comparison+of+the+Cox+model+to+the+Fine-Gray+model+for+survival+analyses+of+re-fracture+rates.">A comparison of the Cox model to the Fine-Gray model for survival analyses of re-fracture rates.</a> Arch Osteoporos. 15 (1), 86 (2020).</li><li>Park, S. Y. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Nomogram:+An+analogue+tool+to+deliver+digital+knowledge.">Nomogram: An analogue tool to deliver digital knowledge.</a> J Thorac Cardiovasc Surg. 155 (4), 1793 (2018).</li><li>Balachandran, V. P., Gonen, M., Smith, J. J., DeMatteo, R. P. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Nomograms+in+oncology:+more+than+meets+the+eye.">Nomograms in oncology: more than meets the eye.</a> Lancet Oncol. 16 (4), e173-e180 (2015).</li><li>Liang, W., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Development+and+validation+of+a+nomogram+for+predicting+survival+in+patients+with+resected+non-small-cell+lung+cancer.">Development and validation of a nomogram for predicting survival in patients with resected non-small-cell lung cancer.</a> J Clin Oncol. 33 (8), 861-869 (2015).</li><li>Chen, D., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Predicting+postoperative+peritoneal+metastasis+in+gastric+cancer+with+serosal+invasion+using+a+collagen+nomogram.">Predicting postoperative peritoneal metastasis in gastric cancer with serosal invasion using a collagen nomogram.</a> Nat Commun. 12 (1), 179 (2021).</li><li>Niu, X., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+prognostic+nomogram+for+patients+with+newly+diagnosed+adult+thalamic+glioma+in+a+surgical+cohort.">A prognostic nomogram for patients with newly diagnosed adult thalamic glioma in a surgical cohort.</a> Neuro Oncol. 23 (2), 337-338 (2021).</li><li>Luo, S., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Comparison+of+left-+and+right-sided+colorectal+cancer+to+explore+prognostic+signatures+related+to+pyroptosis.">Comparison of left- and right-sided colorectal cancer to explore prognostic signatures related to pyroptosis.</a> Heliyon. 10 (7), e28091 (2024).</li><li>Keum, N., Giovannucci, E. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Global+burden+of+colorectal+cancer:+emerging+trends,+risk+factors+and+prevention+strategies.">Global burden of colorectal cancer: emerging trends, risk factors and prevention strategies.</a> Nat Rev Gastroenterol Hepatol. 16 (12), 713-732 (2019).</li><li>Foster, P. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Oestrogen+and+colorectal+cancer:+mechanisms+and+controversies.">Oestrogen and colorectal cancer: mechanisms and controversies.</a> Int J Colorectal Dis. 28 (6), 737-749 (2013).</li><li>Wang, L., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Male-biased+gut+microbiome+and+metabolites+aggravate+colorectal+cancer+development.">Male-biased gut microbiome and metabolites aggravate colorectal cancer development.</a> Adv Sci. 10 (25), e2206238 (2023).</li><li>Kr&#228;mer, H. U., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Type+2+diabetes+mellitus+and+gender-specific+risk+for+colorectal+neoplasia.">Type 2 diabetes mellitus and gender-specific risk for colorectal neoplasia.</a> Eur J Epidemiol. 27 (5), 341-347 (2012).</li><li>Aigner, E., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Gender-+and+site-specific+differences+of+colorectal+neoplasia+relate+to+vitamin+D.">Gender- and site-specific differences of colorectal neoplasia relate to vitamin D.</a> Aliment Pharmacol Ther. 40 (11-12), 1341-1348 (2014).</li><li>Dai, W., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Does+tumor+size+have+its+prognostic+role+in+colorectal+cancer?+Re-evaluating+its+value+in+colorectal+adenocarcinoma+with+different+macroscopic+growth+pattern.">Does tumor size have its prognostic role in colorectal cancer? Re-evaluating its value in colorectal adenocarcinoma with different macroscopic growth pattern.</a> Int J Surg. 45, 105-112 (2017).</li><li>R&#246;ssler, O., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Tumor+size,+tumor+location,+and+antitumor+inflammatory+response+are+associated+with+lymph+node+size+in+colorectal+cancer+patients.">Tumor size, tumor location, and antitumor inflammatory response are associated with lymph node size in colorectal cancer patients.</a> Mod Pathol. 30 (6), 897-904 (2017).</li><li>Zhang, Q., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Prognostic+impact+of+tumor+size+on+patients+with+metastatic+colorectal+cancer:+a+large+SEER-based+retrospective+cohort+study.">Prognostic impact of tumor size on patients with metastatic colorectal cancer: a large SEER-based retrospective cohort study.</a> Updates Surg. 75 (5), 1135-1147 (2023).</li><li>Greene, F. L. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=TNM:+our+language+of+cancer.">TNM: our language of cancer.</a> CA Cancer J Clin. 54 (3), 129-130 (2004).</li><li>Ueno, H., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Optimal+colorectal+cancer+staging+criteria+in+TNM+classification.">Optimal colorectal cancer staging criteria in TNM classification.</a> J Clin Oncol. 30 (13), 1519-1526 (2012).</li><li>Yahagi, M., Okabayashi, K., Hasegawa, H., Tsuruta, M., Kitagawa, Y. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+worse+prognosis+of+right-sided+compared+with+left-sided+colon+cancers:+a+systematic+review+and+meta-analysis.">The worse prognosis of right-sided compared with left-sided colon cancers: a systematic review and meta-analysis.</a> J Gastrointest Surg. 20 (3), 648-655 (2016).</li><li>Holch, J. W., Ricard, I., Stintzing, S., Modest, D. P., Heinemann, V. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+relevance+of+primary+tumour+location+in+patients+with+metastatic+colorectal+cancer:+A+meta-analysis+of+first-line+clinical+trials.">The relevance of primary tumour location in patients with metastatic colorectal cancer: A meta-analysis of first-line clinical trials.</a> Eur J Cancer. 70, 87-98 (2017).</li><li>Petrelli, F., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Prognostic+survival+associated+with+left-sided+vs+right-sided+colon+cancer:+A+systematic+review+and+meta-analysis.">Prognostic survival associated with left-sided vs right-sided colon cancer: A systematic review and meta-analysis.</a> JAMA Oncol. 3 (2), 211-219 (2017).</li><li>Zheng, C., Jiang, F., Lin, H., Li, S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Clinical+characteristics+and+prognosis+of+different+primary+tumor+location+in+colorectal+cancer:+a+population-based+cohort+study.">Clinical characteristics and prognosis of different primary tumor location in colorectal cancer: a population-based cohort study.</a> Clin Transl Oncol. 21 (11), 1524-1531 (2019).</li><li>Karim, S., Brennan, K., Nanji, S., Berry, S. R., Booth, C. M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Association+between+prognosis+and+tumor+laterality+in+early-stage+colon+cancer.">Association between prognosis and tumor laterality in early-stage colon cancer.</a> JAMA Oncol. 3 (10), 1386-1392 (2017).</li><li>Lindberg, L. J., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Risk+of+multiple+colorectal+cancer+development+depends+on+age+and+subgroup+in+individuals+with+hereditary+predisposition.">Risk of multiple colorectal cancer development depends on age and subgroup in individuals with hereditary predisposition.</a> Fam Cancer. 18 (2), 183-191 (2019).</li></ol>

The authors have nothing to disclose.

As a common digestive system tumor, factors influencing the prognosis of solitary colorectal cancer have been studied and confirmed in previous research24. However, there has been limited research on prognostic factors for MPCC patients. This study included 8,931 MPCC patients who underwent surgery from the SEER database between 2004 and 2015. We used a competing risk model to investigate the risk factors affecting CSS and constructed a predictive model. In this study, 39.1% of the deceased patients died due to causes other than cancer, emphasizing the necessity of using a competing risk model to analyze risk factors for CSS.
Previous study has found that advanced age is a risk factor affecting the OS of MPCC patients12. However, in this study, we found that age was not a risk factor for CSS. Elderly patients often exhibit poorer baseline health conditions and more comorbidities, which can lead to a shorter OS compared to younger patients. The elderly are also more likely to die from other causes, such as cardiovascular events and severe infections, making them more susceptible to these factors than younger patients. In this study, we focused on the competing relationship between deaths due to other causes and deaths due to cancer itself. By using a competing risk model and excluding the interference of deaths caused by other reasons, we found that under these conditions, advanced age was no longer a risk factor for CSS of MPCC. This indicates that MPCC patients across different age groups may face similar tumor burdens. Overall, young and elderly patients with this disease have similar CSS, which provides valuable guidance for developing personalized treatment strategies.
This study also found that the incidence and prognosis of MPCC exhibit gender-related tendencies, with a higher proportion and poorer CSS of male patients. This is consistent with the situation observed in solitary colorectal cancer. Previous studies have shown that solitary colorectal cancer is more common in males, and male patients have a worse prognosis than female patients. This may be due to the influence of estrogen on the occurrence and progression of colorectal cancer25. Additionally, other research has suggested that the gut microbiota and gut metabolites in male patients may be one of the reasons for the gender differences observed in colorectal cancer patients26. A study in Germany found that type 2 diabetes has a greater impact on colorectal cancer in female27. What&#39;s more, vitamin D has a protective effect against colorectal neoplasia in females, but no similar findings were found in males28.
Competing univariate analysis found that MPCC patients with larger tumor sizes had poorer CSS. In previous studies, tumor size has often been considered an indicator of tumor aggressiveness 29,30,31. However, the results from competing multivariate analyses did not support tumor size as an independent risk factor. This suggests that in hollow organs like the colorectum, tumor size may have limited ability to reflect tumor aggressiveness. This phenomenon may stem from various factors, including the complexity of tumor biology and the differing patterns of tumor growth and spread across different organs.
Poor tumor grade indicates stronger invasive and migratory abilities of tumor cells, which aligns with the conclusions of this study: a worse tumor grade suggests a poorer CSS for MPCC patients. TNM stage is the most commonly used clinical method for guiding patient treatment and predicting prognosis32,33. This study found that the deeper the tumor infiltration, the greater the number of metastatic lymph nodes, and the presence of organ metastasis, the worse the patient&#39;s cancer-specific survival. This is entirely consistent with clinical consensus.
In colorectal cancer, prognosis varies depending on the location of the tumor. Many studies have found that the right colon has a poorer prognosis compared to the left colorectum34,35,36,37. However, some research suggests that in resectable colorectal cancer, the side of the tumor location does not impact the long-term prognosis38. In this study, MPCC located in the right colon had worse CSS in the short term but a better long-term prognosis. Similarly, MPCC patients who received radiation and chemotherapy had better short-term CSS but worse long-term CSS. SCRC had worse short-term CSS compared to MCRC but better long-term CSS. The reasons for these dual effects on prognosis are still unclear.
This study has some limitations. Due to the inherent constraints of the SEER database, we could not obtain or analyze some known prognostic indicators, such as carcinoembryonic antigen (CEA)levels and microsatellite status. Moreover, we could not exclude patients with inflammatory bowel disease, hereditary nonpolyposis colorectal cancer, and familial adenomatous polyposis who have higher risks of MPCC39. Additionally, immunotherapy has become an increasingly important aspect of colorectal cancer treatment, but we were unable to obtain data on this. The retrospective nature of our analysis and reliance on a single data set may introduce inherent biases.
The study based on the SEER database found that male gender, poor tumor grade, and advanced TNM stage were associated with poorer CSS in MPCC patients after surgery. It is important to pay close attention to patients with these risk factors during diagnosis and treatment. Furthermore, this study developed a nomogram to predict the CSS of MPCC patients, which can accurately predict the prognosis and contribute to clinical treatment decision-making.

Colorectal cancer, as one of the most common digestive system tumors globally, has seen a continuous increase in its incidence over the past few decades. The persistently high incidence and mortality rates have drawn significant attention. According to the latest statistics, colorectal cancer ranks third among common malignant tumors worldwide, with its mortality rate ranking second globally1. Multiple primary colorectal cancer (MPCC) is a special subtype of colorectal cancer that has garnered increasing attention from researchers in recent years. It is defined as the diagnosis of two or more independent colorectal cancer lesions in the same patient, with a distance of more than 5 cm between the edges of the lesions. In MPCC, when multiple tumors are diagnosed simultaneously or within less than 6 months, it is defined as synchronous colorectal cancer (SCRC), while if the interval between diagnoses is greater than 6 months, it is defined as metachronous colorectal cancer (MCRC)2,3,4,5.
The proportion of MPCC within all colorectal cancers is relatively low, and reports on its incidence vary across different regions and studies. Recent studies have reported that MPCC accounts for between 2%-10% of colorectal cancers2,3,4,5. Compared to solitary colorectal cancer, MPCC has a worse prognosis3,6. Currently, clinical assessment of MPCC prognosis mainly relies on the tumor lymph node metastasis (TNM) staging system of the American Joint Committee of Cancer, which uses the most advanced stage of multiple lesions as the actual stage of MPCC. However, basing prognosis predictions solely on TNM staging is insufficient. There is still a lack of an effective tool to predict the prognosis of MPCC post-surgery. In the current era of precision medicine, clinical prediction models for quantifying risk are widely used in clinical decision-making and patient prognosis assessment7,8. In recent years, nomograms based on independent prognostic risk factors have been widely accepted for predicting tumor prognosis9. Nomograms can visualize complex statistical models, making them easier to apply in clinical settings. In the assessment of tumor prognosis, overall survival (OS) and cancer-specific survival (CSS) are commonly used outcome indicators9. OS refers to the time from the confirmation of a cancer in a patient to death from any cause. CSS refers to the time from the diagnosis of a tumor to death caused specifically by the cancer, offering a more precise reflection of the risk of death due to the cancer. When CSS is used as the outcome indicator, deaths caused by other factors can affect the probability of cancer-specific mortality, introducing a competing risk relationship between the two events10,11. Therefore, studies of tumor-specific survival should use competing risk models to eliminate the impact of competing events. Previous studies have constructed models to predict the prognosis of MPCC, but these have been limited to synchronous colorectal cancer and have not utilized competing risk models to account for the impact of competing events on CSS12,13,14.
In this study, we explored the competing independent risk factors affecting cancer-specific survival for MPCC after surgery. The rationale for employing a competing risks model stems from its ability to account for the possibility that patients may die from causes other than cancer, which is critical for obtaining unbiased survival estimates15,16. Traditional survival analysis techniques, such as the Cox model, may overestimate survival probabilities in the presence of competing events, making the use of competing risk models more appropriate in these scenarios17.
Based on the identified independent risk factors, we constructed a nomogram to predict survival probability and validated its performance. Nomograms have gained popularity in clinical settings because they provide a user-friendly graphical representation of complex statistical models, allowing clinicians to easily calculate individualized risk scores18,19. Unlike other predictive tools, nomograms incorporate multiple risk factors and provide more accurate, personalized prognosis estimates. This approach has been increasingly validated across various cancers, demonstrating superior performance over traditional staging systems20,21,22. Our tool aims to assist clinicians in making personalized and more accurate prognosis estimates during diagnosis and treatment, thereby enhancing decision-making in the management of MPCC.

<table border="1" fo:keep-together.within-page="1" fo:keep-with-next.within-page="always">
	<tbody>
		<tr>
			<td>forestplot package</td>
			<td>Comprehensive R Archive Network (CRAN)</td>
			<td>forest plot 3.1.3</td>
			<td>A forest plot that allows for multiple confidence intervals per row, custom fonts for each text element, custom confidence intervals, text mixed with expressions, and more.</td>
		</tr>
		<tr>
			<td>ggprism package</td>
			<td>Comprehensive R Archive Network (CRAN)</td>
			<td>ggprism 1.0.5</td>
			<td>The ggprism package provides various themes, palettes, and other useful functions to customise ggplots and give them the &#8216;GraphPad Prism&#8217; look.</td>
		</tr>
		<tr>
			<td>ggscidca package</td>
			<td>Comprehensive R Archive Network (CRAN)</td>
			<td>ggscidca package</td>
			<td>&#160;The &#39;ggscidca&#39; package adds coloured bars of discriminant relevance to the traditional decision curve. Improved practicality and aesthetics.&#160;</td>
		</tr>
		<tr>
			<td>ggsurvfit package</td>
			<td>Comprehensive R Archive Network (CRAN)</td>
			<td>ggsurvfit 1.1.0</td>
			<td>The ggsurvfit package eases the creation of time-to-event (aka survival) summary figures with ggplot2. The concise and modular code creates images that are ready for publication or sharing.&#160;</td>
		</tr>
		<tr>
			<td>gtsummary package</td>
			<td>Comprehensive R Archive Network (CRAN)</td>
			<td>gtsummary 2.0.0</td>
			<td>The gtsummary package provides an elegant and flexible way to create publication-ready analytical and summary tables using the R programming language. The {gtsummary} package summarizes data sets, regression models, and more, using sensible defaults with highly customizable capabilities.</td>
		</tr>
		<tr>
			<td>QHScrnomo package</td>
			<td>Comprehensive R Archive Network (CRAN)</td>
			<td>QHScrnomo 3.0.1</td>
			<td>The goal of QHScrnomo is to provide functionality to construct nomograms in the context of time-to-event (survival) analysis in the presence of competing risks. It also contains functions to build, validate, and summarize these models.</td>
		</tr>
		<tr>
			<td>R Software</td>
			<td>R Core Team</td>
			<td>R 4.3.3</td>
			<td>Free software environment for statistical computing and graphics</td>
		</tr>
		<tr>
			<td>riskRegression package</td>
			<td>Comprehensive R Archive Network (CRAN)</td>
			<td>riskRegression 1.3.7</td>
			<td>Risk Regression Models and Prediction Scores for Survival Analysis with Competing Risks</td>
		</tr>
		<tr>
			<td>rms package</td>
			<td>Comprehensive R Archive Network (CRAN)</td>
			<td>rms 6.8-1</td>
			<td>&#160;rms does regression modeling, testing, estimation, validation, graphics, prediction, and typesetting by storing enhanced model design attributes in the fit.</td>
		</tr>
		<tr>
			<td>RStudio</td>
			<td>RStudio, Public Benefit Corporation(PBC)</td>
			<td>Rstudio 2023.12.1+402</td>
			<td>Integrated Development Environment (IDE) used for running R scripts, data analysis, and model development. Provides a user-friendly interface for R programming with advanced features like script editing, debugging, and version control.</td>
		</tr>
		<tr>
			<td>SEERstat</td>
			<td>National Cancer Institute (NCI)</td>
			<td>SEERstat 8.4.3</td>
			<td>Software for statistical analysis of SEER and other cancer-related databases</td>
		</tr>
		<tr>
			<td>tidycmprsk package</td>
			<td>Comprehensive R Archive Network (CRAN)</td>
			<td>tidycmprsk 1.0.0</td>
			<td>The tidycmprsk package provides an intuitive interface for working with the competing risk endpoints</td>
		</tr>
	</tbody>
</table>

competing-risk nomogram for predicting cancer-specific survival in multiple primary colorectal cancer patients after surgery

Cancer-specific survival (CSS) in multiple primary colorectal cancer (MPCC) patients is competitively affected by death from other causes. This study aimed to investigate CSS and associated risk factors by competing risk analysis in MPCC patients. The data of this study is from the SEER database. Using univariable and multivariable analysis of the competitive risk model to weaken the impact of competitive events, it explores the risk factors of CSS and develops a nomogram model. Then, the performance of the model is verified by the ROC curve, calibration curve, and DCA. The study encompasses a total of 8931 patients, with 6255 in the training cohort and 2676 in the validation cohort. Univariable and multivariable analyses showed that sex, Tumor Node Metastasis (TNM) stage, and tumor grade are independent risk factors for cancer-specific survival in MPCC patients. Based on the risk factors, we developed a diagram model to predict CSS. ROC curve, calibration curve, and DCA also show good results. In conclusion, a nomogram is developed that serves as a valuable tool for predicting CSS in MPCC patients, providing clinicians with crucial insights for personalized treatment planning.

Author Spotlight: Analyzing Cancer-Specific Survival in Multiple Primary Colorectal Cancer Using Competing Risks

Competing-Risk Nomogram for Predicting Cancer-Specific Survival in Multiple Primary Colorectal Cancer Patients after Surgery

Patient characteristics 
A total of 8,931 patients were included in the study. In terms of sex distribution, there was a higher proportion of male patients (56%) compared to female patients (44%). Regarding tumor location distribution, the majority of tumors were distributed across the total colorectum, while the fewest tumors were located in the right colon (21%). In terms of tumor grade, the most common grade was grade II, with 6,251 patients, accounting for 70% of the total; this was followed by grade III, with 2,026 patients (23%). Grades I and IV had 355 and 299 patients, respectively, accounting for 4% and 3.3% of the total. In terms of tumor invasion, T3 and T4 stages were the most common. More than half of the patients did not have lymph node metastasis, while 17% of the patients experienced distant metastasis. Radiation therapy was received by 1,446 patients (16%), and chemotherapy was administered to 3,843 patients (43%). SCRC was more common, with 5,327 patients accounting for 60% of the total; there were 3,604 patients with metachronous colorectal cancer, accounting for 40%. The patients were randomly divided into training and validation sets in a 7:3 ratio. There were no statistically significant differences in baseline data between the training and validation cohorts, as shown in Table 1.
Univariable analysis 
After controlling the influence of competitive events, the results of the univariate analysis showed sex, tumor grade and size, TNM stage, radiation, chemotherapy, synchronous or metachronous status, and tumor location were the prognostic factors affecting CSS in MPCC patients. Only age is not a prognostic factor for CSS in MPCC patients. We note that there is a significant intersection of CIF in radiation, chemotherapy, location, and synchronous or metachronous status, indicating that the short-term and long-term prognostic effects of radiation on MPCC patients were different, as did chemotherapy, location, and synchronous or metachronous status. The cumulative risk curve of each subgroup is shown in Figure 2.
Multivariable analysis 
The prognostic factors obtained by univariable analysis were incorporated into the best subsets regression (BSR) and multivariable analysis of the competitive risk model. Among them, radiation, chemotherapy, location, and synchronous or metachronous status are excluded because of their dual effects on prognosis. The results of the BSR and multivariable analysis both showed that sex, TNM stage, and tumor grade were independent risk factors for CSS in MPCC patients. (Figure 2)
Construction and verification of Nomogram 
Based on the independent risk factors obtained by multivariable analysis, we construct a line nomogram to predict CSS and verify the performance of the prediction model (Figure 3). Then, we use the ROC curve, calibration curve, and DCA to evaluate the model. The ROC curve showed that the AUCs of 1-year, 3-year, and 5-year of the training cohort were 0.762, 0.742, and 0.734, and the AUCs of 1-year, 3-year, and 5-year of the verification cohort were 0.801, 0.740, and 0.743. In the training cohort and verification cohort, the calibration curve showed a high agreement between the projected probability and the actual data (Figure 4). In order to verify the performance of the model in clinical application, we used DCA to evaluate the clinical value of the model. The results show that the model shows a good net benefit (Figure 5).
<img alt="Figure 1" class="xfigimg" src="/files/ftp_upload/67433/67433fig01.jpg" /> 
Figure 1: Workflow diagram of the study. This study consists of two steps: first, data was obtained using SEER.Stat, and then data analysis and visualization were performed using R. <a href="https://www.jove.com/files/ftp_upload/67433/67433fig01large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Figure 2" class="xfigimg" src="/files/ftp_upload/67433/67433fig02.jpg" /> 
Figure 2: Analysis results. (A) CIF of subgroups. * indicates p &#60;0.05. (B) Best subsets regression. Under the best goodness of fit, sex, grade, and TNM stage were considered for inclusion (shown as black blocks at the top). (C) Multivariate analysis also showed that sex, grade, and TNM stage are independent risk factors. <a href="https://www.jove.com/files/ftp_upload/67433/67433fig02large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Figure 3" class="xfigimg" src="/files/ftp_upload/67433/67433fig03.jpg" /> 
Figure 3: Nomogram for CSS in MPCC patients. The patient&#39;s total score can be calculated by adding the scores corresponding to each factor. Based on the total score, the probability of cancer-specific death at 1, 3, and 5 years can be predicted <a href="https://www.jove.com/files/ftp_upload/67433/67433fig03large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Figure 4" class="xfigimg" src="/files/ftp_upload/67433/67433fig04.jpg" /> 
Figure 4:&#160;ROC curve and calibration curve. (A, B) The curves in training cohort (n=6255) and (C, D) validation cohort (n=2676). The closer the AUC value is to 1, the better the model&#39;s classification performance. The error bars show the 95% confidence interval for the probability of the actual event occurring. <a href="https://www.jove.com/files/ftp_upload/67433/67433fig04large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Figure 5" class="xfigimg" src="/files/ftp_upload/67433/67433fig05.jpg" /> 
Figure 5: DCA for 1-year, 3-year, and 5-year. (A, B, C) DCA in training cohort and (D, E, F) validation cohort. The green line represents the net benefit of all positives, the blue line represents the net benefit of none positives, and the red line indicates the net benefit of the model. The red area below represents the model&#39;s benefit exceeding that of all positives and none positives, indicating the actual benefit range of the model. <a href="https://www.jove.com/files/ftp_upload/67433/67433fig05large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
Table 1: Clinicopathologic and baseline characteristics of patients. <a href="https://www.jove.com/files/ftp_upload/67433/table 1_RE.xlsx" target="_blank">Please click here to download this Table.</a>

The study found that male gender, poor tumor grade, and advanced Tumor Node Metastasis stage were associated with poorer cancer-specific survival (CSS) in multiple primary colorectal cancer (MPCC) patients after surgery. We developed a nomogram to predict the CSS of MPCC patients and contribute to clinical treatment decision-making.

Cancer-specific survival (CSS) in multiple primary colorectal cancer (MPCC) patients is competitively affected by death from other causes. This study ...

Competing-Risk Nomogram for Predicting Cancer-Specific Survival in Multiple Primary Colorectal Cancer Patients after Surgery

Abstract