ORIGINAL RESEARCH article

Front. Med., 31 January 2024

Sec. Gastroenterology

Volume 11 - 2024 | https://doi.org/10.3389/fmed.2024.1301487

New findings in prognostic factor assessment for adenocarcinoma of transverse colon: a comparison study between competing-risk and COX regression analysis

  • 1. Department of General Surgery, The First Affiliated Hospital of Jinan University, Guangzhou, Guangdong, China

  • 2. Department of Health Statistics, School of Public Health, Shanxi Medical University, Taiyuan, China

  • 3. Section of Occupational Medicine, Department of Special Medicine, Shanxi Medical University, Taiyuan, Shanxi, China

  • 4. Department of Clinical Research, The First Affiliated Hospital of Jinan University, Guangzhou, Guangdong, China

  • 5. Guangdong Provincial Key Laboratory of Traditional Chinese Medicine Informatization, Guangzhou, Guangdong, China

Abstract

Purpose:

Competing-risk analysis was used to accurately assess prognostic factors for cancer-specific death in patients with adenocarcinoma of transverse colon (ATC), and the results were compared with those from a conventional Cox regression analysis.

Materials and Methods:

Patients diagnosed with ATC between 2000 and 2019 were selected from the Surveillance, Epidemiology, and End Results database. The crude mortality rates of patients with ATC were calculated and their differences were tested using the Gray’s test, respectively. In performing multivariate analysis, the Cox regression model and the subdistribution hazard function (SD) in competing risk analysis were utilized, respectively.

Results:

This study included 21,477 eligible patients. The SD model indicated that age, etc. are actual independent prognostic factors. In contrast to previous recognition, the results of the Cox regression showed false-positives for sex and Carcinoembryonic antigen, and underestimated point-estimates in the stage and American Joint Committee on Cancer stage due to competing events. A detailed comparison of treatment revealed that the larger surgical scopes were prognostic risk factors compared with the smaller scope of local tumor excision, partial colectomy, or segmental resection. Patients treated with external proton beam radiotherapy had an increased risk compared with those with no radiotherapy and internal radiotherapy.

Conclusions:

After comparing the results of the two methods and mitigating the significant bias introduced by Cox regression, we found independent factors that really affect the prognosis of ATC. On the other hand, in terms of ATC, a larger surgical scope and external proton beam radiotherapy may not improve the long-term survival of patients. Therefore, when faced with ATC patients, these differences should be noted and treated differently from common colorectal cancer patients. Thus, clinicians are able to give more targeted treatment plans and prognostic assessments.

Introduction

Cancer has become the impediment to human longevity and high quality of life. According to the latest global cancer statistics from GLOBOCAN 2020, colorectal cancer ranked third in terms of new cases, with 1.93 million accounting for 10% of all new cancers, and second in terms of deaths, with 940,000 accounting for 9.4% of all deaths due to cancer in that year ().

Transverse colon cancer has been reported to account for about 10% of colorectal cancers (). The transverse colon is located in a special high position, in the middle and anterior part of the entire colon, between the ascending and descending colon, excluding the hepatic and splenic flexures. It has a maximum length of about 50 cm. Since the transverse colon differs from the rest of the colon in terms of embryonic development, anatomical structure, blood supply, and pathogenetic characteristics, it is necessary to clearly delineate the different segments of the colon and to provide precise and individualized treatments according to the specific characteristics of the transverse colon, which is also in line with contemporary medical concepts. However, most studies on colorectal cancer have focused on the ascending and descending colon, which have obvious differences. The transverse colon, which is the link between the two, has received little attention in research. Adenocarcinoma arises from the glandular epithelium, ducts, or secretory epithelium, and is characterized by adenoid structure formation. It is the most common clinical type of colon cancer, accounting for 90–95% of cases, and has a better prognosis than other pathological types. Research on adenocarcinoma of transverse colon (ATC) would therefore be helpful for improving the clinical outcomes.

Kaplan-Meier (KM) analyses and Cox proportional-hazards models are the classical statistical analysis methods used in investigation of the prognostic factors for colorectal cancer. However, advances in medicine and statistical methods are increasing the demand for more-accurate results. It needs to be remembered that dying from cancer is only one of the causes of death for cancer patients, since deaths due to noncancerous diseases and accidents also account for a significant proportion of the causes of death (). It is therefore necessary to consider cancer and non-cancer factors separately when estimating patient mortality. Non-cancer deaths are often considered competing events, and their presence makes the Cox proportional-hazards model inaccurate. Therefore, when analyzing the factors affecting the prognosis of patients with ATC, using a competing-risk analysis will reduce bias and increase the accuracy of the results, thereby more accurately reflecting the true situation.

This study extracted data from the Surveillance, Epidemiology, and End Results (SEER) database on patients diagnosed with ATC (), performed a competing-risk analysis, and compared the results with those of a Cox regression analysis. This protocol allowed for a more-accurate determination of the factors affecting the prognosis of ATC.

Materials and methods

Data collection and patient selection

Surveillance, Epidemiology, and End Results is one of the most-authoritative large oncology registry databases (). We used SEER*Stat software (version 8.4.0) and selected the “Incidence–SEER Research Plus Data, 17 Registries, Nov 2021 Sub (2000–2019)” database, which is derived from 17 registration stations and covers 26.5% of the entire US population. Basic and medical information on patients diagnosed with ATC was extracted, specifically demographic information such as age and race, clinical information such as type of pathology and surgical modality, and survival status. The following inclusion criteria were applied: age ≥ 18 years, diagnosed between 2000 and 2019, tumor located in the transverse colon (C18.4, excluding the hepatic flexure and splenic flexure), microscopy confirmation, and adenocarcinoma (including signet-ring-cell carcinoma). The exclusion criteria were no surgery, survival time 0 months or unknown, multiple primary malignant tumors, not the first tumor, not a primary malignant tumor, or too many incomplete variables (Figure 1). Since the causes of death and survival status of patients are documented in detail in the SEER database, we classify all patients into colon cancer-specific deaths, competing events (other causes of death), and survival (). Applying the inclusion and exclusion criteria resulted in 21,477 patients being included in this study.

FIGURE 1

Statistical analyses

The variables in the baseline data were described using number (N) and percentage values. In the competing-risk analysis, death from ATC and death from other causes was regarded as a competitive relationship. The cumulative-risk rate was estimated in the single-factor analysis using the cumulative incidence function (CIF) described by CIFk(t) = Pr(Tt, D = k), where function CIFk(t) represents the probability of the k-th event occurring before time t and other class events, and D represents the type of events that occur (). Gray’s test was used to perform between-group comparisons (). When no competing events exist, we used Cox regression for the multivariate analyses with the formula log[λ (t)] = log[λ 0 (t)]+ χβ 2, where λ 0(t) is the net risk and λ(t) is the baseline risk function; that is, the risk function when the covariate vector is 0, which can be written as λ (t) = λ 0(t)exp(χβ) ().

When competing events are present and the deletion-independence condition is not satisfied (), we provide the results of both the subdistribution hazard function (SD) and cause-specific (CS) hazard function belonging to competing-risk models (). The formula used for the SD model was , where SD represents the instantaneous probability of the occurrence of the k-th event being observed in the individual at time t (). The formula used for the CS model was , where CS represents the instantaneous probability of a class-k event being observed in the individual who did not experience any event at time t (). The interpretation and usage differ between the above two models, and so the results of both models need to be provided at the same time (). Lau et al. and Koller et al. proposed that the SD model only focus on the absolute incidence of events of interest (, ), and tends to be used in individual risk-prediction studies to estimate the risk and prognosis of a disease, and is suitable for establishing clinical prediction models and risk scores (); The CS model favors etiological studies, with regression coefficients reflecting the relative effects of covariates on the increased incidence of events of interest in subjects in the event-free risk set (). Therefore, the present study mainly adopted the conclusions from the SD model.λk(t) in the formulas for the SD and CS models is the crude risk rate, which is not equal to the net risk rate when competing events are present, and so the hazards ratio (HR) value and 95% confidence interval (CI) obtained using traditional Cox analysis are biased (). In view of this, we also compared the results of the competing-risk analysis with those of classical Cox analysis ().

R Studio software (version 2022.02.3) was used for all statistical analyses. All statistical tests were two-sided, with a probability value of P < 0.05 considered to indicate statistical significance (). The SEER database can be accessed free of charge, and this study was exempted from the need to obtain informed consent from the included patients by the institutional research committee of the First Affiliated Hospital of Jinan University.

Results

Clinical characteristics of patients and survival outcomes

Among the 21,477 patients with ATC who were screened, 10,867 (50.6%) died: (5,923 from ATC and 4,944 from other causes). The survival time was 64.67 ± 58.68 months.

For continuous variables, the average age of the sample was 68.08 ± 13.93 years, number of lymph nodes examined 17.78 ± 11.44, number of positive lymph nodes 1.63 ± 3.44. Table 1 lists the classification variables.

TABLE 1

VariableAll patients (%)Censored (%)Concerned (%)Competition (%)
Total number21,47710,6105,9234,944
Race
White people16,934 (78.85)8,160 (76.91)4,630 (78.17)4,144 (83.82)
Black people2,613 (12.17)1,274 (12.01)840 (14.18)499 (10.09)
Other1,930 (8.99)1,176 (11.08)301 (5.08)453 (9.16)
Sex
Male10,105 (47.05)5,126 (48.31)2,771 (46.78)2,208 (44.66)
Female11,372 (52.95)5,484 (51.69)3,152 (53.22)2,736 (55.34)
Marital status
Married11,316 (52.69)6,177 (58.22)2,945 (49.72)2,194 (44.38)
Unmarried10,161 (47.31)4,433 (41.78)2,978 (50.28)2,750 (55.62)
Grade
I1,774 (8.26)1,045 (9.85)283 (4.78)446 (9.02)
II14,460 (67.33)7,528 (70.95)3,521 (59.45)3,411 (68.99)
III4,243 (19.76)1,611 (15.18)1,761 (29.73)871 (17.62)
IV496 (2.31)191 (1.8)213 (3.6)92 (1.86)
Unknown504 (2.35)235 (2.21)145 (2.45)124 (2.51)
Stage
Localized7,801 (36.32)4,874 (45.94)681 (11.5)2,246 (45.43)
Regional10,624 (49.47)5,253 (49.51)2,871 (48.47)2,500 (50.57)
Distant3,048 (14.19)483 (4.55)2,370 (40.01)195 (3.94)
Unknown4 (0.02)0 (0)1 (0.02)3 (0.06)
AJCC
0186 (0.87)127 (1.2)11 (0.19)48 (0.97)
I3,787 (17.63)2,489 (23.46)208 (3.51)1,090 (22.05)
II8,068 (37.57)4,515 (42.55)1,220 (20.6)2,333 (47.19)
III6,433 (29.95)3,010 (28.37)2,140 (36.13)1,283 (25.95)
IV2,995 (13.95)468 (4.41)2,343 (39.56)184 (3.72)
Unknown8 (0.04)1 (0.01)1 (0.02)6 (0.12)
LymphExcision
Yes21,149 (98.47)10,456 (98.55)5,826 (98.36)4,867 (98.44)
No313 (1.46)151 (1.42)89 (1.5)73 (1.48)
Unknown15 (0.07)3 (0.03)8 (0.14)4 (0.08)
OthersiteSurgery
No19,932 (92.81)10,062 (94.84)5,141 (86.8)4,729 (95.65)
Single resection1,151 (5.36)400 (3.77)612 (10.33)139 (2.81)
Combination81 (0.38)27 (0.25)43 (0.73)11 (0.22)
Surgery, NOS231 (1.08)106 (1.0)82 (1.38)43 (0.87)
Unknown82 (0.38)15 (0.14)45 (0.76)22 (0.44)
PrimSiteSurg
LPS8,617 (40.12)4,111 (38.75)2,310 (39)2,196 (44.42)
SCH11,977 (55.77)6,074 (57.25)3,313 (55.93)2,590 (52.39)
TP504 (2.35)287 (2.7)142 (2.4)75 (1.52)
Surgery, NOS379 (1.76)138 (1.3)158 (2.67)83 (1.68)
Radiotherapy
No21,300 (99.18)10,560 (99.53)5,814 (98.16)4,926 (99.64)
Beam radiation144 (0.67)37 (0.35)94 (1.59)13 (0.26)
Internal radiotherapy9 (0.04)4 (0.04)4 (0.07)1 (0.02)
Radiotherapy, NOS7 (0.03)2 (0.02)3 (0.05)2 (0.04)
Unknown17 (0.08)7 (0.07)8 (0.14)2 (0.04)
Chemotherapy
Yes7,212 (33.58)3,645 (34.35)2,773 (46.82)794 (16.06)
No/unknown14,265 (66.42)6,965 (65.65)3,150 (53.18)4,150 (83.94)
MetsBone
Yes35 (0.16)3 (0.03)31 (0.52)1 (0.02)
No10,992 (51.18)7,395 (69.7)2,217 (37.43)1,380 (27.91)
Unknown10,450 (48.66)3,212 (30.27)3,675 (62.05)3,563 (72.07)
MetsBrain
Yes7 (0.03)1 (0.01)5 (0.08)1 (0.02)
No11,018 (51.30)7,397 (69.72)2,242 (37.85)1,379 (27.89)
Unknown10,452 (48.67)3,212 (30.27)3,676 (62.06)3,564 (72.09)
MetsLiver
Yes938 (4.37)229 (2.16)662 (11.18)47 (0.95)
No10,102 (47.04)7,168 (67.56)1,599 (27.0)1,335 (27.0)
Unknown10,437 (48.60)3,213 (30.28)3,662 (61.83)3,562 (72.05)
MetsLung
Yes179 (0.83)28 (0.26)144 (2.43)7 (0.14)
No10,849 (50.51)7,370 (69.46)2,107 (35.57)1,372 (27.75)
Unknown10,449 (48.65)3,212 (30.27)3,672 (62.0)3,565 (72.11)
MetsDistLN
Yes55 (0.26)25 (0.24)29 (0.49)1 (0.02)
No4,421 (20.58)3,738 (35.23)476 (8.04)207 (4.19)
Unknown17,001 (79.16)6,847 (64.53)5,418 (91.47)4,736 (95.79)
MetsOther
Yes4,313 (20.08)3,680 (34.68)429 (7.24)204 (4.13)
No161 (0.75)82 (0.77)75 (1.27)4 (0.08)
Unknown17,003 (79.17)6,848 (64.54)5,419 (91.49)4,736 (95.79)
TumorSize
≤4 cm10,378 (48.32)5,511 (51.94)2,255 (38.07)2,612 (52.83)
>4 cm11,099 (51.68)5,099 (48.06)3,668 (61.93)2,332 (47.17)
CEA
Negative7,092 (33.02)4,130 (38.93)1,612 (27.22)1,350 (27.31)
Border60 (0.28)31 (0.29)16 (0.27)13 (0.26)
Positive5,701 (26.54)2,945 (27.76)1,648 (27.82)1,108 (22.41)
Unknown8,624 (40.15)3,504 (33.03)2,647 (44.69)2,473 (50.02)
NeuroInvasion
Yes1,539 (7.17)729 (6.87)583 (9.84)227 (4.59)
No12,347 (57.49)7,411 (69.85)2,711 (45.77)2,225 (45.0)
Unknown7,591 (35.34)2,470 (23.28)2,629 (44.39)2,492 (50.4)

Basic characteristics of patients in this study.

Censored: Patients who are alive; Concerned: Patients who died of ATC; Competition: Patients who died of competing events; Unmarried, including single, widowed, divorced and separated. Stage: the most basic way of categorizing how far a cancer has spread from its point of origin. AJCC: American Joint Committee on Cancer; LymphExcision: if regional lymph nodes was removed; Single resection: removal of only regional metastases, or distant lymph nodes, or distant metastases; combination: Any combination of Single resection; Surgery, NOS: Surgery was performed but method is unknown; PrimSiteSurg: Surgery of Primary Site. LPS: Local tumor excision or Partial colectomy or Segmental resection; SCH: Subtotal colectomy or hemicolectomy; TP: Total colectomy or proctocolectomy; Internal radiotherapy: Radioactive implants or Radioisotopes; Radiotherapy, NOS: Radiotherapy was performed but method is unknown; Mets: metastases; MetsBone: bone metastases; CEA: carcinoembryonic antigen.

Univariate analyses

In the presence of competing-risk, we performed univariate analyses using the CIF and Gray’s test, with the results showing that race, marital status, grade of differentiation, summary stage, AJCC stage, surgery of primary site and other site, radiotherapy and chemotherapy status, tumor size, CEA, status of bone metastasis, brain metastasis, liver metastasis, lung metastasis, distant lymph node metastasis, other site metastasis, and neuro invasion exerted significant effects on the prognosis of ATC (P < 0.05). Analyzing the CIF of each variab le at 1 year, 3 years, and 5 years revealed that the CIF of almost all variables increased over time. Among the ordinal classification variables such as grade, the CIF also increased gradually with increasing degree. Detailed data are provided in Table 2. Because age at diagnosis, the number of lymph nodes examined or positive are continuous variables, they were directly included in the multivariate analyses rather than in univariate analyses.

TABLE 2

VariableGray’s testp-valueCIF
1 year3 years5 years
Race37.77<0.0001
White people0.1070.2180.267
Black people0.1020.2490.321
Other0.0780.1820.233
Sex0.020.8825
Male0.0980.2120.273
Female0.1090.2240.269
Marital status41.62<0.0001
Married0.0930.2010.253
Unmarried0.1170.2380.290
Grade743.54<0.0001
I0.0450.1070.152
II0.0770.1810.235
III0.2070.3690.417
IV0.2310.3730.4
Unknown0.1020.2350.283
Stage7027.36<0.0001
Localized0.0220.0480.071
Regional0.0840.2000.261
Distant0.3830.7130.805
Unknown0.1970.250.311
AJCC7418.21<0.0001
00.0160.0360.044
I0.0130.0270.042
II0.0430.0970.133
III0.1050.2530.329
IV0.3850.7180.811
Unknown0.1130.1460.213
LymphExcision4.050.132
Yes0.1040.2180.271
No0.1050.2090.262
Unknown0.2670.40.478
OthersiteSurgery517.84<0.0001
No0.0970.2020.252
Single resection0.1970.4540.537
Combination0.1980.4510.534
Surgery, NOS0.1200.2880.329
Unknown0.2980.4860.499
PrimSiteSurg42.48<0.0001
LPS0.1010.2090.258
SCH0.1040.2210.274
TP0.1070.2350.302
Surgery, NOS0.1730.3290.401
Radiotherapy128.16<0.0001
No0.1030.2160.268
Beam radiation0.3240.5690.636
Internal radiotherapy0.4670.8020.893
Radiotherapy, NOS0.1430.3560.429
Unknown0.1180.3970.548
Chemotherapy619.72<0.0001
Yes0.1040.3010.386
No/unknown0.1040.1770.213
MetsBone290.15<0.0001
Yes0.7710.90.95
No0.0820.1850.237
Unknown0.1240.2460.298
MetsBrain133.84<0.0001
Yes0.8570.8620.875
No0.0840.1870.23866
Unknown0.1240.2470.29829
MetsLiver1,829<0.0001
Yes0.3430.7070.8069
No0.0600.1390.18581
Unknown0.1230.2460.29749
MetsLung654.06<0.0001
Yes0.4730.8300.90761
No0.0780.1770.22818
Unknown0.1240.2460.29798
MetsDistLN142.06<0.0001
Yes0.2960.8180.874
No0.0690.1580.235
Unknown0.1120.2280.279
MetsOther245.65<0.0001
Yes0.2830.7340.875
No0.0640.1450.237
Unknown0.1120.2280.279
TumorSize402.69<0.0001
≤4 cm0.0670.1580.209
>4 cm0.1390.2750.329
CEA71.57<0.0001
Negative0.0850.1840.235
Border0.1240.2300.276
Positive0.1140.2460.301
Unknown0.1130.2280.279
NeuroInvasion252.55<0.0001
Yes0.1490.3420.425
No0.0850.1860.234
Unknown0.1260.2450.296

Univariate analysis of prognostic factors in patients with ATC.

CIF, cumulative incidence function.

Multivariate analyses

All variables that were statistically significant in the univariate analyses (P < 0.05) were entered into a Cox regression analysis and a competing-risk analysis for the multivariate analyses (). Since sex and regional lymphadenectomy were important demographic and surgical information, respectively, they also needed to be added to the multivariate analyses even though they were not statistically significant in univariate analyses ().

Both the Cox regression analysis and the SD model indicated that age, race, grade, stage, AJCC, lymph node excision status, surgery of primary site, radiotherapy and chemotherapy status, marital status, tumor size, the number of lymph nodes examined or positive, status of bone metastasis, lung metastasis, liver metastasis, and neuro invasion were independent factors affecting the prognosis of ATC, and the results for each subgroup were basically consistent in different models. When analyzing the surgery of primary site, the HR values of subtotal colectomy or hemicolectomy, and total colectomy or proctocolectomy gradually increased with the expansion of surgical scope. The same phenomenon was observed in patients who received external proton-beam radiation. Beam radiation became a risk factor, but there was no statistically significant difference between those who did and did not receive internal radiotherapy (). It was particularly interesting that being female and CEA-positive were risk factors in the Cox regression analysis. But in the SD model, there was no statistical difference between male and female and CEA-positive was also not a risk factor. There was no significant difference between patients with CEA-borderline and those who were CEA-negative. Most of the results from the CS model were consistent with those from the SD model. Since the CS model is often used to explore etiological issues, it only played an auxiliary role in the competing-risk analysis of this study, and so is not considered in detail (see Table 3 for more details).

TABLE 3

COX modelSD modelCS model
Prognostic factorsp-valueHR95% CIp-valueHR95% CIp-valueHR95% CI
Age<0.00011.041.0381.042<0.00011.0121.011.015<0.00011.0181.0161.021
Sex
MaleReferenceReferenceReference
Female<0.00010.8280.7960.8620.08970.9540.9031.0070.01120.9330.8840.984
Race
White peopleReferenceReferenceReference
Black people0.00091.1051.0421.1720.01571.0991.0181.1870.01081.1031.0231.19
Other<0.00010.7970.740.8590.02550.8930.8080.9860.00380.8660.7860.955
Marital status
MarriedReferenceReferenceReference
Unmarried<0.00011.2661.2161.3180.00321.0881.0291.15<0.00011.1431.0831.207
Grade
IReferenceReferenceReference
II0.151.0580.981.1440.14671.0940.9691.2350.16981.0890.9641.231
III<0.00011.291.1851.404<0.00011.4871.3071.693<0.00011.5081.3261.714
IV<0.00011.5181.3241.739<0.00011.6521.3681.996<0.00011.7351.4482.079
Unknown0.00741.2121.0531.3950.00021.4831.2031.828<0.00011.5131.2381.851
Stage
LocalizedReferenceReferenceReference
Regional<0.00011.1861.1091.267<0.00011.4891.3261.672<0.00011.5251.3581.712
Distant<0.00011.8481.3772.481<0.00013.1492.1044.714<0.00013.252.2794.634
Unknown0.83690.8540.1913.828<0.00013800.637442.0832674.750.90624105.52705.01E+63
AJCC
0ReferenceReferenceReference
I0.17611.1990.9221.5590.98261.0070.5551.8260.91671.0330.5621.898
II0.00651.4411.1071.8750.01052.1571.1973.8870.00772.2631.2414.129
III<0.00012.0541.5682.691<0.00014.2972.3687.796<0.00014.772.5988.756
IV<0.00014.4352.9896.58<0.00017.543.70515.347<0.00019.1164.55518.246
Unknown0.06522.9830.9349.529<0.00010.00100.0030.91990.00109.97E+56
LymphExcision
YesReferenceReferenceReference
No0.02761.1911.0191.3920.02551.2751.031.5790.00671.3371.0841.65
Unknown0.00242.4621.3754.408<0.00014.1572.4277.121<0.00014.512.2319.117
PrimSiteSurg
LPSReferenceReferenceReference
SCH0.0171.051.0091.0930.00041.1081.0471.172<0.00011.1261.0651.189
TP<0.00011.3581.1831.5590.00011.4441.1981.741<0.00011.4961.261.776
Surgery, NOS0.04171.1461.0051.3050.09891.1610.9721.3860.01291.2291.0451.447
Radiotherapy
noReferenceReferenceReference
Beam radiation0.00191.3671.1221.667<0.00011.5151.2391.8540.00021.4891.2061.839
Internal radiotherapy0.83191.10.4552.6580.62010.8340.4081.7070.80440.8830.3292.369
Radiotherapy, NOS0.3880.6790.2821.6360.85440.8680.1913.9470.69290.7960.2562.475
Unknown0.90540.9630.5171.7930.52051.3150.573.0330.58331.2150.6062.436
Chemotherapy
YesReferenceReferenceReference
No/unknown<0.00011.4711.41.545<0.00011.2441.1661.327<0.00011.4611.3751.552
LymphExamed<0.00010.9830.9810.985<0.00010.9830.980.986<0.00010.9790.9760.982
LymphPositive<0.00011.0661.0611.072<0.00011.0591.0521.066<0.00011.0691.0631.075
MetsBone
YesReferenceReferenceReference
No0.04710.6790.4630.9950.0340.6680.460.970.04680.6730.4550.994
Unknown0.17020.5980.2871.2470.87960.9420.4322.0540.39580.7240.3441.525
MetsLiver
YesReferenceReferenceReference
No<0.00010.7930.720.8730.00040.820.7350.915<0.00010.8130.7320.903
Unknown0.47240.8440.5311.340.87650.9560.5391.6950.49630.8460.5221.371
MetsLung
YesReferenceReferenceReference
No<0.00010.7030.5910.8370.00910.7510.6060.931<0.00010.6940.580.831
Unknown0.11520.6660.4021.1040.08740.5640.2921.0880.02650.5530.3270.933
TumorSize
≤4 cmReferenceReferenceReference
>4 cm<0.00011.1381.091.187<0.00011.1541.0881.223<0.00011.1941.1261.265
CEA
NegativeReferenceReferenceReference
Border0.30521.2110.841.7480.351.2590.7772.0390.43151.2190.7441.999
Positive0.00021.1041.0481.1630.0941.0630.991.1430.01711.0881.0151.166
Unknown0.00241.081.0281.1340.36341.0330.9641.1060.08361.0620.9921.136
NeuroInvasion
YesReferenceReferenceReference
No0.00020.8630.7990.9310.03140.8960.810.990.00380.8720.7950.957
Unknown0.00060.8630.7940.9390.10450.910.8111.020.03070.8890.7990.989

Multivariate analysis of 3 models of prognostic factors in patients with ATC.

HR, hazard radio; CI, confidence interval; CS, cause-specific hazard; SD, subdistribution hazard function.

Discussion

Colorectal cancer is becoming more and more common. In terms of anatomical sites, there have been many studies about the ascending colon, descending colon, and sigmoid colon, while research into the transverse colon is rare, and so many aspects of this part are still uncertain. Adenocarcinoma is one of the most common types. Therefore, it is necessary to perform further related investigations.

The findings of the present study were consistent with previous studies finding that age, etc., are independent factors influencing the prognosis (). The HR values and p-values for tumor size were consistent in the three models, which also confirms the results of previous studies, indicating that patients with tumors larger than 4 cm have a relatively poorer prognosis (). There is general agreement that racial disparities in health insurance and medical care result in black people having a higher risk. The conclusion that other races have a better prognosis than white people in this study may differ from the results of some other studies, which may be caused by genetic differences under similar medical conditions. The grade of differentiation also significantly affected the prognosis. In addition to moderate differentiation (grade II), grades III and IV had significant effects on the survival rate compared with well differentiation (grade I). Some previous studies on a similar topic have applied competing-risk analyses, and found a significant difference between grades I and II in SD and CS models, but not in the Cox model, which can be interpreted as a false-negative result (). Whereas, the results obtained in the present study for the three models indicated no significant difference between grades II and I. We believe that this discrepancy is attributable to differences in the characteristics or sizes of the selected samples, and it could also be due to inherent characteristics of ATC itself.

The statistical results of above variables are consistent in the range of P-values in the three models of this study. The correlation direction of risk factors and results was consistent. HR and 95% CI in Cox regression were slightly lower but still basically similar. It is worth noting that the point estimation and interval estimation of some factors differed markedly among the three models. For example, Cox regression analysis underestimated the risk of each level of AJCC by almost half, compared with the HR values for the SD and CS models, and the 95% CI were also correspondingly lower. This echoes what other researchers have found (). Compared with the results of Cox regression and competing-risk analysis in the “stage-distant” group, the degree of risk underestimation was surprisingly similar. Obviously, due to the existence of competing events (), Cox regression deviates markedly in both point estimation and interval estimation, whereas the results of the competing-risk analysis are more accurate ().

All of the independent related factors mentioned above have been widely recognized in previous studies. Our study not only further supports these conclusions, but also has produced some new findings by using competing-risk analysis. Many previous studies of colorectal cancer considered being female as a protective factor, and the present Cox regression analysis produced the same result. Based on the large population, this was not the case for the SD model, which indicated no difference in prognosis between males and females with ATC. Our analysis suggests that the incidence is slightly lower in females than males, but no difference was found in prognosis. This is consistent with the findings of Cheung et al. (). We therefore consider this to be a false-positive result caused by competing events.

In terms of treatment, surgery has previously been considered a protective factor, but differences among specific surgical methods have not been discussed. Based on the records in the SEER database, we divided the surgical procedures into three categories. For all three models, SCH and TP with a larger surgical scope did not appear to benefit patients with ATC, but the risk increased as the surgical scope expanded relative to LPS covering a smaller scope. A larger surgical scope may tend to be applied more when the disease advanced, and is also associated with a greater tissue damage, resulting in a poor prognosis (). Compared with the ascending colon and descending colon, the operation method of transverse colon is not so fixed, since a small deviation of the scope of hemicolectomy forward or backward can involve the ascending or descending colon (). Sometimes the prognosis varies markedly depending on whether or not the middle colon artery is preserved (). The results of the present competing-risk analysis further confirmed the results of the Cox regression analysis, hence making the results more robust and reliable.

The same was true for the variable of radiotherapy. Previous studies have only cared about whether patients received radiotherapy or not, and their conclusions are controversial. Some researchers have suggested that radiotherapy improves survival (). However, studies that found radiotherapy to not be beneficial for long-term survival did not explore the underlying reasons. In contrast, the present study analyzed detailed information on radiotherapy, and found that the increased risk is mainly due to the use of external proton beam, which is the most common type of external radiotherapy, while internal radiotherapy such as I125 did not affect the prognosis. This may be due to external radiotherapy killing some tumor cells, reducing the tumor load, relieving symptoms, and improving short-term survival. However, due to the side effects of external radiotherapy, such as perforation, bleeding, and pancytopenia, and the generally late stage of patients who need radiotherapy, the long-term prognosis might not be improved (). The present competing-risk analysis confirmed the results of the Cox regression analysis again.

Carcinoembryonic antigen is considered to be an important tumor marker for the diagnosis and monitoring of recurrence and metastasis, which has become a broad clinical consensus, but it often presents false-positive or false-negative results. Regarding the effect of CEA on prognosis, there is still a lot of controversy. The Cox regression analysis and the competing-risk analysis performed in this study produced different results. Elevated CEA was considered a risk factor for a poor prognosis in Cox regression, but no such difference was observed in the SD model. CEA assessment plays an important role in postoperative follow-up, in terms of detecting recurrence and metastasis early, thus improving the excision rate. However, this does not mean a higher or lower survival rate. It often needs to be combined with other indicators to evaluate the prognosis. Relevant studies by Ohlsson et al. () in Sweden and Kjeldsen et al. () in Denmark have confirmed this conclusion. Based on our results, particularly the findings from the SD model, it is evident that CEA cannot be used for prognostic prediction. This aligns with the prevailing views among most clinical experts and scholars. The clinical utility of CEA is therefore limited to the two points mentioned above. Our results once again validate the conclusions drawn by the majority of scholars. Relying solely on the results of Cox regression, however, can lead to significant errors. The Cox regression incorrectly estimates P-values due to the existence of competing events, and the results of competing-risk analysis will be more consistent with reality. It is well known that the location of colorectal cancer is closely related to the biological characteristics, genetic and epigenetic characteristics, pathological characteristics, and the prognosis of tumor cells (). Tumors vary by site over time, and various features of tumors will also change gradually along the colon segment.

The transverse colon is quite special because it is located between the ascending colon and descending colons, traverses the upper abdomen, and entangled by the transverse mesocolon, which is an internal peritoneal organ. It abuts many important organs which differ from that of the ascending colon and descending colon. In addition, it is mainly supplied by the middle colon artery, and the type of surgical method is directly related to whether the artery is preserved or not. Regional lymph nodes are mainly distributed along the middle colonic artery (), which makes lymph node dissection in this region a considerable surgical challenge. Therefore, the above-mentioned new findings in this study cannot only be attributed to the advantages of competing-risk analysis, but also reflect potential differences between ATC and cancers in other parts of the colon.

The classic KM and Cox methods rely on the assumption that censoring time and failure time are independent (); that is, there is a single endpoint without competing events. However, in reality, clinical research data often contains a substantial amount of right-censored data due to loss of follow-up and other reasons, leading to multiple outcomes with competitive relationships. Using classical KM for univariate analyses without considering these biases will overestimate cumulative mortality, while using Cox regression for multivariate analyses can lead to further biases. (). Our study has confirmed this theory. Some Cox regression results may have false-positives and biased effect estimates due to serious bias at a proportion of competing events exceeding 10%, with lower proportions also potentially resulting in false-positive or false-negative results (). Competing-risk analysis can overcome these shortcomings by establishing the dependence between correlation degree and covariates, which can better and more accurately explain the effect of covariates and standardize the distribution function for different types of competing risk (). It is noteworthy that the HR and 95% CI values obtained from the SD model were similar to those derived from the CS model, with consistent direction of associations and effect sizes that are basically consistent with the theory of SD ≤ CS; however, some results were inconsistent, which is also consistent with some studies. For instance, P-values for some factors in the SD and CS models were consistent, while others were not. This highlights the importance of obtaining results using both models in a competing-risk analysis, which helps to further differentiate the role of risk factors (). Broadly speaking, CS addresses upstream epidemiological questions related to disease etiology, while SD focuses on downstream clinical event rates. The SD model is primarily used for prognostic analysis, risk scoring, and clinical prediction modeling. Therefore, it is crucial to use competing-risk analysis when analyzing the prognostic risk factors of patients where competing events are present ().

As far as we know, this study is the first to report on the competing-risk analysis of prognostic factors for ATC, specifically using the largest number of samples and variables. It is based on the high-quality and large SEER database (), allowing for the identification of accurate prognostic factors for specific diseases like ATC. Moreover, it provides valuable insights for clinicians in assessing prognosis and avoiding harmful treatment strategies. Additionally, it serves as a reminder to researchers about the significant inaccuracies associated with the use of Cox regression. The study also has some limitations. First, we only selected patients from 2000 to 2019, which may have introduced bias due to the short time span. Second, although the SEER database contains a significant amount of variable information, it does not cover all information that may affect patient survival, such as gene expression (). Thus, further research is still needed to address these limitations.

Conclusion

Upon comparative analysis of the two methodologies, it provided conclusive evidence that age, etc. are the actual prognostic factors for ATC. However, sex and CEA do not qualify as independent prognostic factors. When analyzing prognostic factors with multiple endpoints, competing-risk analysis is more accurate and reliable than COX regression, which is prone to significant bias in the presence of competing events. Additionally, larger surgical scopes and external proton-beam radiotherapy may not improve long-term survival outcomes for patients with ATC. Therefore, clinicians should take note of these differences when treating ATC patients and may need to approach them differently from common cases of colorectal cancer. This study specifically examined ATC patients in detail, in contrast to previous crude analyses of prognostic factors for colorectal cancer. These results will contribute to a deeper understanding of ATC, as well as to diagnostic, therapeutic, and clinical decision-making processes.

Statements

Data availability statement

Publicly available datasets were analyzed in this study. This data can be found here: seer.cancer.gov/.

Ethics statement

Ethical approval was not required for the study involving humans in accordance with the local legislation and institutional requirements. Written informed consent to participate in this study was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and the institutional requirements.

Author contributions

HS: Data curation, Formal analysis, Software, Visualization, Writing—original draft, Writing—review and editing. SX: Software, Writing—review and editing. SW: Writing—review and editing. LH: Investigation, Visualization, Writing—review and editing. JL: Conceptualization, Writing—review and editing. YP: Conceptualization, Writing—review and editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by the Guangzhou Science and Technology Plan City-School Joint Funding Project (202201020084 and 202201020065), the Clinical Frontier Technology Program of the First Affiliated Hospital of Jinan University, China (No. JNU1AF-CFTP-2022a01223), and the Guangdong Provincial Key Laboratory of Traditional Chinese Medicine Informatization (2021B1212040007).

Acknowledgments

We thank the Surveillance, Epidemiology, and End Results (SEER) program staff for providing open access to the database.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

  • 1.

    SungHFerlayJSiegelRLaversanneMSoerjomataramIJemalAet alGlobal Cancer Statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 Countries.CA. (2021) 71:20949.

  • 2.

    WestNKobayashiHTakahashiKPerrakisAWeberKHohenbergerWet alUnderstanding optimal colonic cancer surgery: comparison of Japanese D3 resection and European complete mesocolic excision with central vascular ligation.J Clin Oncol. (2012) 30(15):17639.

  • 3.

    BaudinoT. Targeted cancer therapy: the next generation of cancer treatment.Curr Drug Discov Technol. (2015) 12:320.

  • 4.

    YangJPanZHeYZhaoFFengXLiuQet alCompeting-risks model for predicting the prognosis of penile cancer based on the SEER database.Cancer Med. (2019) 8:78819.

  • 5.

    LeiSGeYTianSCaiBGaoXWangNet alColorectal cancer metastases to brain or bone and the relationship to primary tumor location: a population-based study.J. Gastrointest Surg. (2020) 24: 183342.

  • 6.

    WuWLiYFengALiLHuangTXuAet alData mining in clinical big data: the frequently used databases, steps, and methodological models.Military Med Res. (2021) 8:44.

  • 7.

    WuWYangJLiDHuangQZhaoFFengXet alCompetitive risk analysis of prognosis in patients with cecum cancer: a population-based study.Cancer Control. (2021) 28:1073274821989316.

  • 8.

    HallerBSchmidtGUlmK. Applying competing risks regression models: an overview.Lifetime Data Analysis. (2013) 19:3358.

  • 9.

    MerrillRCapocacciaRFeuerEMariottoA. Cancer prevalence estimates based on tumour registry data in the Surveillance, Epidemiology, and End Results (SEER) Program.Int J Epidemiol. (2000) 29:197207.

  • 10.

    GoodwinJSametJHuntW. Determinants of survival in older cancer patients.J Natl Cancer Inst. (1996) 88:10318.

  • 11.

    WongCLiewTHanzoL. Turbo coded burst by burst adaptive wideband modulation with blind modem mode detection.Proceedings of the IEEE Transactions on Vehicular Technology.Sorrento: IEEE (1999). p. 3038.

  • 12.

    YangJLiYLiuQLiLFengAWangTet alBrief introduction of medical database and data mining technology in big data era.J Evid Based Med. (2020) 13:5769.

  • 13.

    EsnaolaNGebregziabherMFinneyCFordM. Underuse of surgical resection in black patients with nonmetastatic colorectal cancer: location, location, location.Ann Surg. (2009) 250:54957.

  • 14.

    LauBColeSGangeS. Competing risk regression models for epidemiologic data.Am J Epidemiol. (2009) 170:24456.

  • 15.

    KollerMRaatzHSteyerbergEWolbersM. Competing risks and the clinical community: irrelevance or ignorance?Stat Med. (2012) 31:108997.

  • 16.

    LiangCHeagertyPJ. A risk-based measure of time-varying prognostic discrimination for survival models.Biometrics. (2017) 73:72534.

  • 17.

    IslamiFKamangarFBoffettaP. Grand challenges in cancer epidemiology and prevention.Front Oncol. (2011) 1:3. 10.3389/fonc.2011.00003

  • 18.

    AvisP. Pressure homogenisation of mammalian cells.2nd ed. In: BirnieGeditor.Subcellular Components.Oxford: Butterworth-Heinemann (1972). p. 113.

  • 19.

    HashimDWeiderpassE. Cancer survival and survivorship.3rd ed. In: BoffettaPHainautPeditors. Encyclopedia of Cancer.Oxford: Academic Press (2019). p. 2509.

  • 20.

    FlickKSubletteCYip-SchneiderMMaatmanTColgateCSoufiMet alInsurance type and marital status impact hospital length of stay after pancreatoduodenectomy.J Surg Res. (2021) 257:58792.

  • 21.

    DistasoAAbatangeloLMagliettaRCreanzaTPiepoliACarellaMet alBiological and functional analysis of statistically significant pathways deregulated in colon cancer by using gene expression profiles.Int J Biol Sci. (2008) 4:36878.

  • 22.

    LongoCDe RobertisRValdoMMontemezziS. A case of acute bilateral iodine-induced submandibular sialadenitis: ultrasound findings.J Clin Ultrasound. (2022) 50:703.

  • 23.

    ChottanapundSChamroonsawasdiKTunyasitthisundhornPAekplakornWSilpasuwanPAnantachotiPet alModifiable factors and colon cancer risk in Thai population.Asian Pac J Cancer Prevent. (2021) 22: 3743.

  • 24.

    LiangYLiQHeDChenYLiJ. Tumor size improves the accuracy of the prognostic prediction of T4a stage colon cancer.Sci Rep. (2021) 11:16264.

  • 25.

    HuCCaoJZengLLuoYFanH. Prognostic factors for squamous cervical carcinoma identified by competing-risks analysis: a study based on the SEER database.Medicine. (2022) 101:e30901.

  • 26.

    LiangPChenTGiovannucciE. Cigarette smoking and colorectal cancer incidence and mortality: systematic review and meta-analysis.Int J Cancer. (2009) 124:240615.

  • 27.

    ZhaoXKallakuryBChahineJHartmannDZhangYChenYet alSurgical resection of SCLC: prognostic factors and the tumor microenvironment.J Thoracic Oncol. (2019) 14:91423.

  • 28.

    CheungWShiQO’ConnellMCassidyJBlankeCKerrDet alThe predictive and prognostic value of sex in early-stage colon cancer: a pooled analysis of 33,345 patients from the ACCENT database.Clin Colorectal Cancer. (2013) 12:17987.

  • 29.

    SuHWangSXieSHuangLPanYLyuJ. Prediction of death probability in adenocarcinoma of the transverse colon: competing-risk nomograms based on 21,469 patients.J Cancer Res Clin Oncol. (2023) 149:1043552.

  • 30.

    MurphyCHarlanLLundJLynchCGeigerA. Patterns of colorectal cancer care in the United States: 1990-2010.J Natl Cancer Inst. (2015) 107:djv198.

  • 31.

    AndersenBStimecBEdwinBKazaryanAMaziarzPIgnjatovicD. Re-interpreting mesenteric vascular anatomy on 3D virtual and/or physical models: positioning the middle colic artery bifurcation and its relevance to surgeons operating colon cancer.Surg Endoscopy. (2022) 36:1008.

  • 32.

    GuanXJiangZMaTLiuZHuHZhaoZet alRadiotherapy dose led to a substantial prolongation of survival in patients with locally advanced rectosigmoid junction cancer: a large population based study.Oncotarget. (2016) 7:2840819.

  • 33.

    RijkmansEvan TriestBNoutRKerkhofEBuijsenJRozemaTet alEvaluation of clinical and endoscopic toxicity after external beam radiotherapy and endorectal brachytherapy in elderly patients with rectal cancer treated in the HERBERT study.Radiother Oncol. (2018) 126:41723.

  • 34.

    OhlssonBBrelandUEkbergHGraffnerHTranbergK. Follow-up after curative surgery for colorectal carcinoma. Randomized comparison with no follow-up.Dis Colon Rectum. (1995) 38:61926.

  • 35.

    KjeldsenBKronborgOFengerCJørgensenOD. A prospective randomized study of follow-up after radical surgery for colorectal cancer.Br J Surg. (1997) 84:6669.

  • 36.

    YoungSGolzarianJ. Primary tumor location in colorectal cancer: comparison of right- and left-sided colorectal cancer characteristics for the interventional radiologist.Cardiovasc Intervent Radiol. (2018) 41:181925.

  • 37.

    LangmanGPatelABowleyD. Size and distribution of lymph nodes in rectal cancer resection specimens.Dis Colon Rectum. (2015) 58:40614.

  • 38.

    GiacconeGSoriaJ-C.Targeted Therapies in Oncology.2nd ed. Boca Raton, FL: CRC Press (2013).

  • 39.

    NieZ. A new perspective of clinical survival data: competitive risk model.Chin J Epidemiol. (2017) 38:112731.

  • 40.

    JungKWonYHongSKongHLeeE. Prediction of cancer incidence and mortality in Korea, 2020.Cancer Res Treat. (2020) 52:3518.

  • 41.

    ZhangJYangJWangHPanZYanXHuCet alDevelopment and validation of a nomogram for osteosarcoma-specific survival: a population-based study.Medicine. (2019) 98:e15988.

  • 42.

    SubramaniamDLiuSGiacconeG. Novel approaches in cancer immunotherapy.Discov Med. (2016) 21:26774.

Summary

Keywords

adenocarcinoma of transverse colon, competing-risk analysis, SEER, Cox regression, subdistribution hazard function

Citation

Su H, Xie S, Wang S, Huang L, Lyu J and Pan Y (2024) New findings in prognostic factor assessment for adenocarcinoma of transverse colon: a comparison study between competing-risk and COX regression analysis. Front. Med. 11:1301487. doi: 10.3389/fmed.2024.1301487

Received

25 September 2023

Accepted

08 January 2024

Published

31 January 2024

Volume

11 - 2024

Edited by

Muhammad Suleman, Qatar University, Qatar

Reviewed by

Pralay Mukhopadhyay, GlaxoSmithKline, United States

Zafar Ali, University of Swat, Pakistan

Updates

Copyright

*Correspondence: Yunlong Pan, Jun Lyu,

†These authors have contributed equally to this work

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics