CRC-SPIN Fred Hutch Cancer Center

Simulate the natural history of colorectal cancer, including diagnosis and survival after diagnosis, to evalute screening strategies.

Contact: Carolyn Rutter crutter@fredhutch.org

Summary

The Colorectal Cancer Simulated Population model for Incidence and Natural history (CRC-SPIN; Fred Hutchinson Cancer Center) describes discrete event state transitions in continuous time.1,2 State transitions can depend on the age and sex of the individual and lesion location within the large intestine. For each simulated individual, CRC-SPIN first generates a time of birth and a time of death. Next, CRC-SPIN generates adenomas within the individual using a non-homogeneous Poisson process. The risk of developing a new adenoma is based on a log-linear model that depends on age and sex. Each individual is assigned a baseline log-risk that is normally distributed, which represents their individual risk level. Once initiated, CRC-SPIN adenomas are randomly assigned a location in the large intestine.

CRC-SPIN simulates continuous adenoma size and models adenoma growth using a Richard’s growth model with a minimum detectable adenoma size equal to 1 mm and a maximum (but rarely reached) adenoma size of 50mm.3 For each adenoma, CRC-SPIN simulates a time to reach 10 mm using an inverse Weibull distribution. The time to reach 10 mm depends on adenoma location (colon or rectum). The adenoma growth rate is a closed function of the (randomly distributed) time to reach 10 mm. Adenomas of any size may progress to preclinical CRC, although most do not progress within an individual’s simulated lifetime and transition at a small size is rare. The size at transition is a function of sex, location, and age at adenoma initiation.

The model simulates the time spent in the preclinical cancer state (sojourn time) using a Weibull distribution with shape parameter fixed at 5 (to focus on distributions with a limited range of skewness) and scale parameter (which determines mean sojourn time) that depends on adenoma location (colon or rectum). Individuals can have multiple preclinical cancers, and, given location, the sojourn time of each preclinical cancer is independent.

Once a cancer becomes clinically detectable, the size and stage at clinical detection is simulated based on observed SEER data from 1975-1979 (a period before screening). Given cancer size at clinical detection, the size at any point during the preclinical period is determined using an exponential growth model, assuming a minimum cancer size of 0.5 mm and replacement of adenoma cells with cancer cells until the cancer overtakes the adenoma. Stage at screen detection is assigned based on tumor size, with the restriction that the screen-detected stage must be no later than the stage at clinical detection.

After detection, CRC-SPIN simulates the survival time to death from CRC using relative survival estimates from SEER, following an approach that is common across the three models.4 CRC survival time depends on age at detection, calendar year at detection, stage at detection, cancer location (colon or rectum), and sex. For individuals with synchronous CRCs at the time of diagnosis, CRC-SPIN uses the stage-specific survival of the cancer with the highest stage present. For individuals with CRC, their death date is set to the earliest simulated death date (either due to CRC death or other causes).

Risk Factor Modeling

The CRC-SPIN model includes a risk factor input that can be incorporated to increase or decrease cancer risk in the overall population or in subgroups of the population. Risk factor inputs can affect each state transition via multiplicative effects similar to accelerated failure time models. That is, risk factors can affect adenoma incidence, adenoma growth, transition to preclinical cancer, and sojourn time. This approach allows inclusion of basic (essential) risk factors, such as race, in the model. The model currently incorporates this risk information as an input (estimated outside the model) rather than a calibrated parameter (chosen in conjunction with other parameters to match target data). The CRC-SPIN model is currently being revised to incorporate differences in CRC risk by race/ethnicity.

References

  1. Rutter CM, Ozik J, DeYoreo M, Collier N. Microsimulation model calibration using incremental mixture approximate Bayesian computation. The Annals of Applied Statistics. 2019 Dec 1;13(4):2189–2212. PMCID: PMC8534811
  2. DeYoreo M, Rutter CM, Ozik J, Collier N. Sequentially calibrating a Bayesian microsimulation model to incorporate new information and assumptions. BMC Med Inform Decis Mak. 2022 Jan 12;22(1):12. PMCID: PMC8756687
  3. Tjørve E, Tjørve KMC. A unified approach to the Richards-model family for use in growth analyses: Why we need only two model forms. Journal of Theoretical Biology. Elsevier; 2010 Dec;267(3):417–425. PMID: 20831877
  4. Rutter CM, Johnson EA, Feuer EJ, Knudsen AB, Kuntz KM, Schrag D. Secular trends in colon and rectal cancer relative survival. J Natl Cancer Inst. 2013;105(23):1806-13.