Novel computational biology modeling system can accurately forecast response to neoadjuvant therapy in early breast cancer
Breast Cancer Research volume 25, Article number: 54 (2023)
Generalizable population-based studies are unable to account for individual tumor heterogeneity that contributes to variability in a patient’s response to physician-chosen therapy. Although molecular characterization of tumors has advanced precision medicine, in early-stage and locally advanced breast cancer patients, predicting a patient’s response to neoadjuvant therapy (NAT) remains a gap in current clinical practice. Here, we perform a study in an independent cohort of early-stage and locally advanced breast cancer patients to forecast tumor response to NAT and assess the stability of a previously validated biophysical simulation platform.
A single-blinded study was performed using a retrospective database from a single institution (9/2014–12/2020). Patients included: ≥ 18 years with breast cancer who completed NAT, with pre-treatment dynamic contrast enhanced magnetic resonance imaging. Demographics, chemotherapy, baseline (pre-treatment) MRI and pathologic data were input into the TumorScope Predict (TS) biophysical simulation platform to generate predictions. Primary outcomes included predictions of pathological complete response (pCR) versus residual disease (RD) and final volume for each tumor. For validation, post-NAT predicted pCR and tumor volumes were compared to actual pathological assessment and MRI-assessed volumes. Predicted pCR was pre-defined as residual tumor volume ≤ 0.01 cm3 (≥ 99.9% reduction).
The cohort consisted of eighty patients; 36 Caucasian and 40 African American. Most tumors were high-grade (54.4% grade 3) invasive ductal carcinomas (90.0%). Receptor subtypes included hormone receptor positive (HR+)/human epidermal growth factor receptor 2 positive (HER2+, 30%), HR+/HER2− (35%), HR−/HER2+ (12.5%) and triple negative breast cancer (TNBC, 22.5%). Simulated tumor volume was significantly correlated with post-treatment radiographic MRI calculated volumes (r = 0.53, p = 1.3 × 10–7, mean absolute error of 6.57%). TS prediction of pCR compared favorably to pathological assessment (pCR: TS n = 28; Path n = 27; RD: TS n = 52; Path n = 53), for an overall accuracy of 91.2% (95% CI: 82.8% – 96.4%; Clopper–Pearson interval). Five-year risk of recurrence demonstrated similar prognostic performance between TS predictions (Hazard ratio (HR): − 1.99; 95% CI [− 3.96, − 0.02]; p = 0.043) and clinically assessed pCR (HR: − 1.76; 95% CI [− 3.75, 0.23]; p = 0.054).
We demonstrated TS ability to simulate and model tumor in vivo conditions in silico and forecast volume response to NAT across breast tumor subtypes.
Despite pharmacologic and genomic progress in oncology, tumor heterogeneity hinders the ability to optimize therapy . Molecular tumor characterization has advanced targeted therapeutics; despite this, expression of target biomarkers in any given patient does not necessarily ensure a predictable or durable response to therapy. Further compoudning challenges to optimizing therapy include the inability, a process of escalating or de-escalating therapy when it will result in better or comparable results with similar or fewer side effects, include the inability of large-scale trials with finite enrollment criteria to adequately address individual tumor characteristics (morphology [2, 3], vascularity [4, 5], location , tumor microenvironment ) and deliver personalized oncologic treatment recommendations that comprehensively address these variations in individual tumor conditions resulting from tumor heterogeneity.
Increasingly neoadjuvant therapy is being employed and International  and European  standard-of-care (SOC) guidelines for clinical practice continue to evolve. In the USA, these SOC regimens are designated by guidelines set forth by the National Comprehensive Cancer Network (NCCN) . A rising number of SOC regimens, including new drug approvals and combinations of existing drugs, have become available for a single indication. These expanded coverage options for breast cancer are largely driven by an accelerated awareness of the complexity of tumor biology and a quest to increase offerings to the medical oncology community. What remains unknown is which SOC treatment will engender the best response in an individual patient and provide the best opportunity to achieve pathological complete response or a reduction in residual tumor burden. While some ineffectiveness may be attributed to suboptimal-guideline adherence [11, 12], even with adherence using clinical decision support systems to guide the decision-making process, clinical outcomes can vary within any sub-cohort given the same regimen. 
An innovative approach is needed to forecast treatment response and identify the likelihood of success of any chosen treatment. Accurately predicting tumor response affords numerous opportunities for the clinician in terms of increasing confidence around treatment selection, by providing information at the time of diagnosis to assist in understanding the treatment plan. This understanding will be particularly meaningful for patients who may be exposed to additional toxicity but would respond to an escalated regimen or those who may benefit from treatment de-escalation with comparable disease outcomes .
Currently, most methodologies appropriate clinical outcomes by stratifying patients according to risk of recurrence (ROR) and overall survival (OS) based on risk factors identified on retrospective analyses or -omics data [15,16,17,18]. However, modern day computational biology now affords a unique opportunity to use cutting-edge mathematical modeling to approximate in vivo biological processes in silico with a high degree of accuracy. Much of the field of perfusion kinetics [19, 20] has been derived from advances in medical imaging which now makes it possible to acquire dynamic, high spatial resolution images that have advanced multi-scale tumor modeling potential. These multi-scale computational models are now able to simulate and closely mirror fundamental cancer biological processes and predict spatiotemporal changes . The next logical step in resolving the gap is in its predictive capacity to capture the dynamic environment on the time axis to forecast spatiotemporal changes of a patient’s tumor as it responds to chemotherapeutic agents.
Here, we perform an independent assessment of a previously validated biophysical tumor modeling platform , TumorScope Predict (TS). TS constructs a three-dimensional (3D), dynamic model of an individual’s tumor. By integrating modeling of tumor morphology and metabolism, the platform simulates biological processes and interactions that take place within the tumor microenvironment including vascularity, nutrient availability, drug delivery, sensitivity and resistance . The integration of these complex interactions enables multimodal forecasting of tumor response over time to a given NAT.
In the current study, we further validate the platform’s ability to accurately evaluate and predict pCR, a surrogate marker for long-term outcome , in an independent breast cancer cohort receiving NAT for early-stage or locally advanced breast cancer (hereafter, referred to collectively as “early stage”). Additionally, we have undertaken efforts to further our understanding of the predictive capacity of TS, spanning residual tumor morphology, and ROR for breast cancer subtypes.
Materials and methods
Patients included in the study were age 18 or older, diagnosed with breast cancer, from September 2014 through December 2020, treated at a single institution who completed NAT with a SOC chemotherapy regimen, and had pre-treatment T1-weighted dynamic contrast enhanced magnetic resonance imaging (DCE-MRI). Patients with all subtypes and histology of breast cancer were included as well as those with bilateral breast tumors. One bilateral breast cancer case was analyzed as two separate tumors in the model, each with a distinct volumetric and pCR prediction; characteristics for the two tumors were analyzed to determine underlying biology. Retrospective data encompassing imaging and baseline diagnostic information were separated temporally into two time points: pre-treatment and post-treatment. The pre-treatment data (timepoint 1) was used for all analysis and computation. The post-treatment data (timepoint 2) was used for validation. Outcomes were single-blinded and run prospectively using timepoint 1. As de-identified patient data was used for the study; an institutional review board (IRB) exemption was granted.
A total of 81 tumors from 80 patients were analyzed. Patients with metastatic disease, those receiving neoadjuvant endocrine therapy alone, or those on clinical trials with experimental therapies or regimens were excluded from the study. The American Joint Committee on Cancer (AJCC) 8th edition staging system was used for all clinical and pathological assessments reported herein .
Tumor segmentation and model design
Briefly, the previously described TS platform  combines artificial intelligence for tumor segmentation from the surrounding tissues with biophysical simulations to simulate the tumor response to therapy. It constructs a 3D virtual (in silico) tumor model of a patient (e.g., a virtual “twin”) that can be used to simulate how a patient will respond to a particular therapy using only baseline demographic, pathological, and medical imaging data. To create the 3D model, a pre-treatment DCE-MRI is processed using a semantic segmentation convolutional neural network (CNN)—a type of deep learning model with multiple successive layers that analyzes a neighborhood of voxels (e.g., pixels) and distils properties of them in increasingly abstract ways to assign a single class to each voxels —to create a 3D representation of the breast tissues. In this case, the CNN is a multi-class volumetric residual UNet (e.g., a ResVNet , that examines each voxel in the MRI and classifies it as comprising primarily of tumor, vasculature, fibroglandular tissue, adipose, skin, or chest, producing a grid of cubic voxels that are 0.5 mm3 (the “spatial model”; see Fig. 1C, D) which forms the basis of the virtual twin. Briefly, the ResVNet consists of three encoding tiers, a bottleneck, and three decoding tiers and takes as input a 48 × 48 × 96 cubic voxel region with three channels (pre-contrast, early and late post-contrast) of the DCE-MRI. It was trained using fivefold cross-validation on the training set with a cross-entropy loss function; then, each fold was combined via weighted geometric means to produce the final model weights. None of the patients in this study were used in training the CNN. A recent publication describes the model training and validation .
Once the 3D model is created, the virtual twin is further personalized by incorporating patient demographic (e.g., age, race/ethnicity), and pathological characteristics (e.g., T (tumor size) and N (lymph node spread) stage, estrogen receptor (ER) percent staining, progesterone receptor (PR) percent staining, HER2 status, grade, histology type). These are used to select a mathematical model representing the patient’s tumor biology. Tumor biology algorithms were developed to describe the tumor’s requirements for nutrients (e.g., how quickly they consume nutrients and produce byproduct chemicals), interaction with neighboring healthy tissues (e.g., competition for, and cross feeding of, nutrients with fatty and glandular tissues), intrinsic growth rate (e.g., how fast the tumor doubles in size when consuming differing amounts of available nutrients), and susceptibility to drugs (e.g., pharmacodynamics). The combination of the tumor biology algorithms and the 3D spatial model constitutes the whole patient-specific virtual twin. The details of this process are described previously .
A completed virtual twin can then be simulated using the previously described biophysical model  with a NAT drug regimen. The simulation progresses by dosing drugs at regular intervals as defined by the National Comprehensive Cancer Network (NCCN) guidelines . Drug dosing is captured by pharmacokinetic equations, which introduce the drug into the tissues according to a model based on the local vascular perfusion in a time dependent manner . Vascular perfusion at each voxel is computed via a Tofts-like model that is parameterized from the DCE-MRI . The simulation progresses by computing the local nutrient environment (i.e., concentration of metabolites like glucose, oxygen, alanine, etc.) and drugs at a specific time during treatment (see Fig. 1G). The growth/death rate of the tumor is computed as a function of the local nutrient and drug concentrations via the tumor biology model. The tumor size and shape are updated in response to the growth/death rate, and time is advanced by an increment of 1 h. This process is repeated until the entire course of treatment has been simulated (e.g., 60–180 days depending on drug regimen) [28,29,30,31]. The result is a 4D (3D + time) simulated picture of how the tumor responds to treatment.
In this study, the drug regimen prescribed for a patient by the physician was simulated. In total, 81 tumors in 80 patients were simulated. The final pre-surgery simulation tumor size and shape (see Fig. 1F) are used to determine whether the patient has pCR or RD as described in the next section.
Model volume prediction and validation
Simulation volume for the model was defined as the sum of all volume fractions across all voxels that have disease. This sum was then multiplied by the physical volume of the voxel (in mm3) to generate the predicted, post-treatment volume.
For the radiologist-assessed volume, the tumor region of interest (ROI) was depicted (segmented) in 3D. Then, the number of voxels included in the ROI was summed and multiplied by the physical volume of the voxel (in mm3). Statistical analysis comparing the two models was then carried out.
Virtual twins (of each individual tumor) were simulated with the physician-prescribed regimen. All analyses were performed on individual tumors; for patients with bilateral cancer, separate predictions were made for the tumor in each breast. Baseline demographics between the pCR and RD subgroups were compared using a chi-squared test (Table 1). Primary outcomes were the pCR/residual prediction along with predicted final volume for each tumor, which were both assessed at the timepoint 2 (e.g., the total simulated time duration was timepoint 2-timepoint 1). pCR predictions were compared to pathologist assessments reported in post-surgery pathology notes. Predicted volumes were compared to volumes segmented from post-treatment (timepoint 2, pre-surgery) MRIs. Segmented volumes were assessed by a board-certified radiologist (10 + years) specializing in breast cancer.
The primary outcome metrics of positive predictive value (PPV), negative predictive value (NPV), sensitivity, specificity, and accuracy for pCR in the overall population—as well as in individual molecular subtypes [HR+/HER2−, HR+/HER2+, HR−/HER2+, and triple negative breast cancer (TNBC)]—were computed along with 95% confidence intervals (CI) estimated using the Clopper-Pearson exact binomial interval . The accuracy of both residual tumor volume and percentage reduction in tumor volume as predictors of pCR was measured with the area under the receiver operating characteristic curve (AUROC). A pre-defined cutoff was used for binary pCR prediction; namely, a residual predicted tumor volume less than or equal to 0.01 cm3 or a greater than 99.9% reduction in predicted final tumor volume compared to diagnosis was considered pCR. This cutoff was pre-defined from feasibility studies of over 600 patients [22, 29, 33,34,35,36]. Predicted pCR was compared to standard pathologist assessments of pCR, defined as ypT0N0 per the AJCC staging system 8th edition . To explore the ability of TS-predicted pCR to discriminate event-free survival (EFS), a long-rank test between survival curves for predicted pCR and predicted residual disease (RD) patients was used.
TS predictions were also correlated with radiographic response, assessed with the Pearson’s correlation coefficient with significance tested using the Fisher transformation. Statistics were computed with the Python package statsmodels v0.12.2. All statistical testing was done with α = 0.05, and false discovery correction was performed separately for descriptive analysis of demographic data and analysis of outcome metrics (sensitivity, specificity, correlation of volume with predictions, and survival outcomes)—using the Benjamini–Hochberg method with a false discovery rate of 0.05.
A total of 89 breast cancer patients who received NAT and had corresponding pre-treatment DCE-MRIs were identified at a single institution. Nine patients were excluded due to missing pre-treatment DCE-MRIs (5 patients), missing cancer histology (2 patients), or missing pCR calls (2 patients) (Fig. 1).
Eighty patients were eligible for inclusion in the generated for analysis (Table 1, Fig. 1), of which 36 (45.0%) were self-identified as African American women and 40 (50.0%) were self-identified Caucasian women. The median patient age was 53 years. Most tumors were high grade (54.4% grade 3) invasive ductal carcinomas (90%). The most common receptor subtype was HR+/HER2− (35%); additional receptor subtypes in the cohort were HR+/HER2+ (30%), TNBC (22.5%), and HR−/HER2+ (12.5%). The majority of patients were T2 (71%) and had either no spread to nearby lymph nodes (N0, 30%) or spread to 1–3 axillary lymph nodes (N1, 60%). Within this cohort, TNBC patients experienced higher rates of pCR, while HR+/HER2− patients had higher rates of RD (p < 1 × 10–11). Most of the pCR responses were observed in high-grade tumors (Table 1), consistent with other studies [37, 38] showing that higher grade (e.g., highly proliferative) tumors respond better to chemotherapy. Achievement of pCR varied across the different drug regimens that the patients received (p < 1 × 10–8).
The demographic and cancer characteristics of the population in this study were statistically different (χ2 test) when compared to a previous  study. Patients differed in racial/ethnic compositions, with a lower proportion of self-reported Caucasian patients (p = 0.011), a lower proportion of ductal versus lobular cancers (p = 0.005), a higher proportion of HR+/HER2+ and lower proportion of TNBC patients (p < 0.001). Moreover, this population had a higher proportion of grade 1 and 2 tumors (p < 0.001), a lower proportion of T1 versus T2 tumors (p = 0.004), and a lower proportion of N0 versus N1 tumors (p < 0.001).
Previously acquired pre-treatment diagnostic data for these patients (demographic characteristics, drug regimen information, imaging (DCE-MRI), and pathology data) were input into TS (Fig. 2).
TS then modeled the weekly volumetric response throughout the specific treatment regimen for each patient up to the point of surgery. Predictions of tumor volume strongly correlated with radiographic assessment of tumor volume for follow-up (n = 53) MRIs obtained after treatment for each patient (r = 0.53, p < 1.3 × 10–7; Fig. 3A, Additional file 1: Figs. S1 and S3). It should be noted that follow-up MRIs were used solely for validation, not as input to TS.
The mean absolute error in predicted tumor volume change for all follow-up MRIs was 6.57% demonstrating high accuracy of TS predictions (Table 2, Fig. 3B top and C). While error in the volume is difficult to interpret without a reference error rate, agreement between radiologically-assessed tumor size and pathologically-assessed tumor size has been examined . Simulated volume errors were compared to the level of agreement (LOA) between breast MRIs and histopathology reported previously  (Fig. 3B bottom). For both inter-regimen and post-treatment MRI scans the mean volume error remained under 3.7% (Table 2) with most of the error distribution within the LOA.
Based on simulated tumor volumetric response, TS was then evaluated for its ability to predict pCR via the pre-defined tumor volume threshold 0.01 cm3 or percent reduction of 99.9% or greater. pCR was predicted in 28 patients, while RD was predicted in the remaining 52 patients (Fig. 4). The device predicted pCR/RD for 73 of 80 patients correctly, for an overall accuracy of 91.2 (95% CI 82.8–96.4%; Clopper-Pearson interval).
Device performance assessed via AUROC was found to be state-of-the-art (AUROC = 0.91; range 0.75–0.94). TS performance was robust throughout receptor subtypes, with the highest pCR prediction accuracy observed in TNBC patients at 93.8% (95% CI 55.5–99.8%) with a sensitivity of 90%, and the lowest seen in HR+/HER2− patients where the accuracy was 75% (95% CI 68–93.2%) corresponding to the lowest sensitivity of 50% (Table 3, Additional file 1: Fig. S2). The lower performance in HR+/HER2− patients was observed in a prior study , however the performance relative to other subtypes was not as extreme. The algorithm performance was also robust across chemotherapy regimens with a prediction accuracy of over 89% for both anthracycline (n = 33) and non-anthracycline (n = 47) regimens (Additional file 1: Table 1).
Next, probability of EFS was modeled to assess the prognostic ability of TS prediction of pCR versus RD compared to pathological assessment of pCR versus RD. TS-simulated pCR was associated with a 5-year ROR with a hazard ratio (HR) = − 1.99 [ − 3.96, − 0.02] (95% CI; p = 0.043) while the clinically assessed pCR data had a HR = − 1.76 [ − 3.75,0.23] (95% CI; p = 0.054) demonstrating similar prognostic performance between TS predictions and pathological assessment (Fig. 4).
Advances in breast cancer oncology within the last decade have been largely in genomic profiling [15, 40, 41]. There is burgeoning recognition of the intricacies of tumor heterogeneity [42,43,44] that confounds expected treatment response to targeted pharmaceutical management, and newer oncologic insights are manifesting at an unprecedented rate [45, 46]. Existing unimodal technologies, including single cell biology [47, 48], organoids [49, 50] and spatial biology [47, 51] platforms, attempt to recapitulate varying aspects of multi-omics, biomarker expression patterns, and the complexity that stems from cells organizing and interacting in the 3D tumor microenvironment. Here, we perform a study of an independent cohort using a previously validated biophysical platform, TumorScope Predict. Predicated on mathematical modeling algorithms that incorporate dynamic, high-resolution imaging [52,53,54,55,56], this technology confers capacity to identify and manage tumor heterogeneity, such as region-specific vascular density, perfusion, simulated nutrient availability and simulated drug delivery to forecast tumor response on a case-by-case basis (see Tumor Segmentation and Model Design).
While traditional radiomics, the quantitative analysis of medical images, is a field of active investigation and interest, it relies on static 2D or 3D images for analysis. For example, a multivariate study  using MRI to predict pCR found that for all patients receiving NAT, pCR prediction failed to reach an AUROC over 0.660 , suggesting only modest gains in prediction performance. Radiomic features were only considered for prediction which may also explain lower performance. The addition of pre-treatment clinical information with machine learning increases the accuracy of pCR predictions , however, this approach remains suboptimal as it still fails to account for tumor progression and response throughout therapy. Even in vitro methods, such as organoids and in vivo methods, such as xenografts, which explore the behavior of living cells in response to a variety of targeted therapies, still lack the fundamental capacity to analyze tumor heterogeneity and complexity in relation to the dynamic communication and interplay between genetically or transcriptionally distinct regions of a tumor. Instead, both organoids and xenografts approach tumors as essentially a series of clonal experiments, without incorporating the higher-order behavior of a structurally diverse tumor. Thus, these methods do not truly address the heterogenous nature of the tumor that is thought to be a critical driver of therapy response.
The TumorScope technology works by amalgamating the fields of perfusion kinetics and computational biology to simulate the biophysical interactions between the tumor, surrounding tissues, nutrients and drugs in a 3D model. The inputs for the technology include non-invasive high-resolution image acquisition (DCE-MRIs) and standard of care pre-treatment information including demographics, pathological report data (including available molecular markers). Once basal conditions are established, the platform simulates physician-chosen drug regimens with fundamental in vivo cancer biological processes (e.g., metabolic activity) in a multi-scale biophysical computational model over the axis of time to create an adaptive, 4D model (3D model over time).
Precedence for this technology exists in preclinical models; Zangooei et al. demonstrated that microscale computed tomography (microCT) enabled computational simulations (multi-scale modeling) to recapitulate the dynamics of tumor growth in a rodent model . In the current study, we found that TS-simulated tumors, as compared with radiographically-assessed tumors from MRIs, demonstrated highly accurate prediction of both tumor volume and percent change in tumor volume throughout therapy, with low median absolute error (Fig. 3, Table 2). Tumor volume has been used in several approaches to assess clinical outcomes either as an independent prognostic indicator  or in combination with other prognostic indicators such as pCR , residual cancer burden (RCB) , RFS [62, 63], progression-free survival (PFS)  and OS [64,65,66,67]. Previous studies have used initial volume parameters from 2D images obtained from medical images (MRIs and PET/CT scans) to predict long-term outcomes (pCR [68, 69] and RFS [62, 63]) in breast cancer [62, 63] lung [68, 69] and rectal cancer . However, the accuracy of these methodologies has only reached an AUROC ranging from 0.62  to 0.73 . The TS platform uses multi-scale multimodal approach to predict the tumor’s response over time; attained pCR predictions with predictive accuracy ranging from 0.893 to 0.944, associated with an overall AUROC of at least 0.91. Importantly, overall performance of TS was robust in all breast cancer subtypes, indicating the tool may offer value in the NAT decision-making process. This substantial improvement relative to existing technologies not only advances the field of precision oncology, but additionally increases the likelihood of meaningful impact on clinical care.
More importantly, these volume predictions are forecasted based on physician-chosen drug dose and timing of standard of care chemotherapeutic regimens; newer agents (immunotherapy, antibody–drug conjugates) are actively being added into the platform’s drug repertoire to accommodate the changing landscape of pharmaceutical management of breast cancer. TS demonstrated high fidelity and accuracy near or above 90% in prediction of tumor response to both anthracycline and non-anthracycline based regimens. As anthracyclines have been a major target for de-escalation of therapy, given their known cardiotoxicity and association with secondary leukemia, these results are encouraging to inform decision-making and assist with understanding response to physician-chosen regimen. The proposed biophysical modeling approach thus confers domain over an expanse of clinical support areas ranging from de-escalation to prognostic predictions related to tumor volume such as pCR, RD and associated RFS.
While pCR at the clinical trial cohort level may not be generally indicative of an overall superior therapy [23, 71], it has proved to have meaningful prognostic value for individual patients, particularly for the TNBC and HER2+ subtypes . TS-generated predictions of RFS derived from pCR mirrors RFS from clinically assessed pCR (Fig. 4). As TS predictions [23, 73, 74] rely on volume-based metrics (tumor volume rendered in 3D, changes in the tumor over time in 4D), these outputs provide insights on residual disease including morphology and volume. Thus, TS has the capacity to capture residual cancer burden (RCB)  in breast cancer, as well as potential for use in other cancers with the Response Evaluation Criteria in Solid Tumors (RECIST)  criteria. Given these outputs, volume-based tumor responses can be stratified into risk categories and assist in evaluating long-term prognosis.
There are certain limitations to our model. Foremost, prediction of pCR using TS is performed on analysis of the primary tumor only; no axillary lymph node disease is considered. However, studies have shown that lymph node histology is similar to that of the primary tumor; therefore, modeling the primary tumor could be representative of the actual tumor response . Indeed, Fayanju et al. found that in a study of over 20,000 patients, only 1.5% of patients achieved pCR in the primary tumor but not the lymph nodes , and thus, one might anticipate that the error in overpredicting pCR would be small. Also, it depends on the availability and quality of the medical images. Not all patients are indicated to receive a pre-treatment DCE-MRI scan. Given that roughly 5.6% (5/89) of patients were excluded due to a missing pre-treatment MRI, this is hurdle for use, and potentially a source of bias in the study. Additionally, MRI quality will depend on the machines and techniques used and can affect the precision of the technology. Our platform models NAT per standard-of-care recommended regimens, and in this assessment potential treatment interruptions, delays and changes in medication may not have been considered. Finally, the nature of a single site study with a limited number of patients has resulted in wide confidence intervals; however, our results remain consistent with TS performance observed in concurrent validation studies [22, 33] with a multi-center analysis that addresses a much larger cohort of patients, currently underway. One assumption of the study is that multiple tumors in a single patient are independent, which may bias the statistics. One African American patient with HR+/HER2− invasive ductal carcinoma had bilateral cancer. These tumors differed only in that the left tumor was grade 2, while the right was grade 3. Both tumors were predicted (correctly) to have residual disease, which may slightly skew the statistics for the HR+/HER2− cohort. Despite any limitations imposed by these constraints, TS performance was robust across all subtypes. While the focus of this validation study is breast cancer, the potential to use biophysical simulation modeling approaches for other solid tumors is an area of active investigation for our team with encouraging preliminary results.
This validation study demonstrates that TS accurately predicts patient-specific tumor volume response to NAT across all breast cancer subtypes. These capabilities allow the platform to predict pCR in a patient-specific manner, which could in turn be used to optimize chemotherapy regimens and escalation/de-escalation decisions, predict downstaging, and inform the physician–patient discussion. By predicting tumor volume across the timescale of cancer therapy, TS offers predictive capacity with regards to pCR, RCB and EFS and is a useful technology for adjuvant therapy planning. Changes in volume, morphology, and particularly size (diameter) derived from the simulations in the current study with associated prognostic metrics can be used in patients diagnosed with other types of solid tumors, advancing precision oncology across a spectrum of disease.
While this and previous studies demonstrate method performance in a limited set of independent centers, a broader multi-center validation consisting of a significantly larger patient population is necessary to account for patient and clinical care variability seen in current practice.
Availability of data and materials
The data generated in this study are available within the article and its supplementary data files.
Area under the receivers operating curve
Dynamic contrast enhanced magnetic resonance imaging
Human epidermal growth factor receptor 2
Hormone receptor positive
Hormone receptor negative
Institutional review board
Pathological complete response
Residual cancer burden
Response evaluation criteria in solid tumors
Risk of recurrence
Dagogo-Jack I, Shaw AT. Tumour heterogeneity and resistance to cancer therapies. Nat Rev Clin Oncol. 2017;15(2):81–94. https://doi.org/10.1038/nrclinonc.2017.166.
Beck AH, Sangoi AR, Leung S, et al. Imaging: systematic analysis of breast cancer morphology uncovers stromal features associated with survival. Sci Transl Med. 2011. https://doi.org/10.1126/SCITRANSLMED.3002564.
Laws A, Pastorello R, Dey T, et al. Impact of the histologic pattern of residual tumor after neoadjuvant chemotherapy on recurrence and survival in stage I–III breast cancer. Ann Surg Oncol. 2022;29(12):7726–36. https://doi.org/10.1245/S10434-022-12054-6.
Oshi M, Newman S, Tokumaru Y, et al. Intra-tumoral angiogenesis is associated with inflammation, immune reaction and metastatic recurrence in breast cancer. Int J Mol Sci. 2020;21(18):6708. https://doi.org/10.3390/IJMS21186708.
Jacquemier J, Charafe-Jauffret E, Monville F, et al. Association of GATA3, P53, Ki67 status and vascular peritumoral invasion are strongly prognostic in luminal breast cancer. Breast Cancer Res. 2009;11(2):1–11. https://doi.org/10.1186/BCR2249.
Ji F, Xiao WK, Yang CQ, et al. Tumor location of the central and nipple portion is associated with impaired survival for women with breast cancer. Cancer Manag Res. 2019;11:2915. https://doi.org/10.2147/CMAR.S186205.
Chitalia RD, Rowland J, McDonald ES, et al. Imaging phenotypes of breast cancer heterogeneity in preoperative breast dynamic contrast enhanced magnetic resonance imaging (DCE-MRI) scans predict 10-year recurrence. Clin Cancer Res. 2020;26(4):862–9. https://doi.org/10.1158/1078-0432.CCR-18-4067.
Burstein HJ, Curigliano G, Thürlimann B, et al. Customizing local and systemic therapies for women with early breast cancer: the St. Gallen International Consensus Guidelines for treatment of early breast cancer 2021. Ann Oncol Off J Eur Soc Med Oncol. 2021;32(10):1216–35. https://doi.org/10.1016/J.ANNONC.2021.06.023.
Guidelines | ESMO. Accessed March 2, 2022. https://www.esmo.org/guidelines.
Network NCC. Breast cancer (version 05.2020) (2020).
Rocque GB, Williams CP, Jackson BE, et al. Impact of nonconcordance with NCCN guidelines on resource utilization, cost, and mortality in de novo metastatic breast cancer. J Natl Compr Cancer Netw. 2018;16(9):1084–91. https://doi.org/10.6004/JNCCN.2018.7036.
Ricci-Cabello I, Vásquez-Mejía A, Canelo-Aybar C, et al. Adherence to breast cancer guidelines is associated with better survival outcomes: A systematic review and meta-analysis of observational studies in EU countries. BMC Health Serv Res. 2020;20(1):1–12. https://doi.org/10.1186/S12913-020-05753-X.
Klarenbeek SE, Weekenstroo HHA, Sedelaar JPM, Fütterer JJ, Prokop M, Tummers M. The effect of higher level computerized clinical decision support systems on oncology care: a systematic review. Cancers. 2020;12(4):1032. https://doi.org/10.3390/CANCERS12041032.
Piccart MJ, Hilbers FS, Bliss JM, et al. Road map to safe and well-designed de-escalation trials of systemic adjuvant therapy for solid tumors. J Clin Oncol. 2020;38(34):4120–9. https://doi.org/10.1200/JCO.20.01382.
Sparano JA, Paik S. Development of the 21-gene assay and its application in clinical practice and clinical trials. J Clin Oncol. 2008;26(5):721–8. https://doi.org/10.1200/JCO.2007.15.1068.
Mamounas EP, Tang G, Fisher B, et al. Association between the 21-gene recurrence score assay and risk of locoregional recurrence in node-negative, estrogen receptor-positive breast cancer: results from NSABP B-14 and NSABP B-20. J Clin Oncol. 2010;28(10):1677. https://doi.org/10.1200/JCO.2009.23.7610.
Rossi C, Cicalini I, Cufaro MC, et al. Breast cancer in the era of integrating “Omics” approaches. Oncogenesis. 2022;11(1):1–13. https://doi.org/10.1038/s41389-022-00393-8.
Osdoit M, Yau C, Symmans WF, et al. Association of residual ductal carcinoma in situ with breast cancer recurrence in the neoadjuvant I-SPY2 trial. JAMA Surg. 2022;157(11):1034–41. https://doi.org/10.1001/JAMASURG.2022.4118.
Johnson GB, Harms HJ, Johnson DR, Jacobson MS. PET imaging of tumor perfusion: a potential cancer biomarker? Semin Nucl Med. 2020;50(6):549–61. https://doi.org/10.1053/J.SEMNUCLMED.2020.07.001.
Baboli M, Zhang J, Kim SG. Advances in diffusion and perfusion MRI for quantitative cancer imaging. Curr Pathobiol Rep. 2019;7(4):129–41. https://doi.org/10.1007/S40139-019-00204-7.
Zangooei MH, Margolis R, Hoyt K. Multiscale computational modeling of cancer growth using features derived from microCT images. Sci Rep. 2021;11(1):1–17. https://doi.org/10.1038/s41598-021-97966-1.
Howard FM, He G, Peterson JR, et al. Highly accurate response prediction in high-risk early breast cancer patients using a biophysical simulation platform. Breast Cancer Res Treat. 2022;196(1):57–66. https://doi.org/10.1007/S10549-022-06722-0.
Conforti F, Pala L, Sala I, et al. Evaluation of pathological complete response as surrogate endpoint in neoadjuvant randomised clinical trials of early stage breast cancer: systematic review and meta-analysis. https://doi.org/10.1136/bmj-2021-066381.
Amin MB, Frederick H, Greene L, et al. The eighth edition AJCC cancer staging manual: continuing to build a bridge from a population-based to a more “personalized” approach to cancer staging. CA Cancer J Clin. 2017;67(2):93–9. https://doi.org/10.3322/CAAC.21388.
Yadav SS, Jadhav SM. Deep convolutional neural network based medical image classification for disease diagnosis. J Big Data. 2019;6(1):1–18. https://doi.org/10.1186/S40537-019-0276-2.
Milletari F, Navab N, Ahmadi SA. V-Net: fully convolutional neural networks for volumetric medical image segmentation. In: Proceedings of 2016 4th international conference 3D vision, 3DV 2016. Published online June 15, 2016, pp 565–571. https://doi.org/10.48550/arxiv.1606.04797.
Tofts PS, Kermode AG. Measurement of the blood-brain barrier permeability and leakage space using dynamic MR imaging. 1. Fundamental concepts. Magn Reson Med. 1991;17(2):357–367. https://doi.org/10.1002/MRM.1910170208
Cole JA, Peterson JR, Earnest TM, et al. Perfusion kinetics from clinical DCE mris increase the accuracy of predictions of tumor response to chemotherapy.. 2020;38(15_suppl):e12651–e12651. https://doi.org/10.1200/JCO.2020.38.15_SUPPL.E12651.
Cole JA, Peterson JR, Earnest TM, et al. Abstract P1–08–31: Simbiosys tumorscope: Biophysical modeling of patient-specific response to chemotherapy. Cancer Res. 2022;82(4_Supplement):P1-08-31. https://doi.org/10.1158/1538-7445.SABCS21-P1-08-31.
Jarrett AM, Kazerouni AS, Wu C, et al. Quantitative magnetic resonance imaging and tumor forecasting of breast cancer patients in the community setting. Nat Protoc. 2021;16(11):5309–38. https://doi.org/10.1038/s41596-021-00617-y.
Bowers HJ, Douglas E, Ansley K, Thomas A, Weis JA. Dynamic characterization of breast cancer response to neoadjuvant therapy using biophysical metrics of spatial proliferation. Sci Rep. 2022;12(1):11718. https://doi.org/10.1038/s41598-022-15801-7.
Clopper CJ, Pearson ES. The use of confidence or fiducial limits illustrated in the case of the binomial. Biometrika. 1934;26(4):404. https://doi.org/10.2307/2331986.
Howard FM, He G, Pfeiffer JR, et al. Abstract P4-05-03: Evaluation of the prognostic accuracy of SimBioSys TumorScope in early breast cancer. Cancer Res. 2022;82(4_Supplement):P4-05-03. https://doi.org/10.1158/1538-7445.SABCS21-P4-05-03.
Cole JA, Peterson JR, Earnest TM, Hallock MJ, Braun E. Abstract P106–04: SimBioSys TumorScope: Spatio-temporal modeling of the breast tumor microenvironment accurately predicts chemotherapeutic response. Cancer Res. 2020;80(4_Supplement):P1-06-04. https://doi.org/10.1158/1538-7445.SABCS19-P1-06-04.
Cole JA, Peterson JR, Earnest TM, et al. SimBioSys TumorScope: Spatio-temporal modeling of the tumor microenvironment to predict chemotherapeutic response. J Clin Oncol. 2020;38(15_suppl):e12650–e12650. https://doi.org/10.1200/jco.2020.38.15_suppl.e12650.
Mamillapalli CK, Markwell TK, Ellis JK, et al. Spatiotemporal modeling with SimBioSys TumorScope to predict chemotherapeutic response in breast tumor microenvironments. 2020;38(15_suppl):e12656–e12656. https://doi.org/10.1200/JCO.2020.38.15_SUPPL.E12656.
Masood S. Neoadjuvant chemotherapy in breast cancers. 2016;12(5):480–491. https://doi.org/10.1177/1745505716677139.
Alba E, Lluch A, Ribelles N, et al. High proliferation predicts pathological complete response to neoadjuvant chemotherapy in early breast cancer. Oncologist. 2016;21(2):150–5. https://doi.org/10.1634/THEONCOLOGIST.2015-0312.
Lobbes MBI, Lalji UC, Nelemans PJ, et al. The quality of tumor size assessment by contrast-enhanced spectral mammography and the benefit of additional breast MRI. J Cancer. 2015;6(2):144–50. https://doi.org/10.7150/JCA.10705.
Van’t Veer LJ, Dai H, Van de Vijver MJ, et al. Gene expression profiling predicts clinical outcome of breast cancer. Nature. 2002;415(6871):530–6. https://doi.org/10.1038/415530a.
Bernard PS, Parker JS, Mullins M, et al. Supervised risk predictor of breast cancer based on intrinsic subtypes. J Clin Oncol. 2009;27(8):1160–7. https://doi.org/10.1200/JCO.2008.18.1370.
Alizadeh AA, Aranda V, Bardelli A, Blanpain C, Bock C, Borowski C. Toward understanding and exploiting tumor heterogeneity. Nat Med. 2015;21(8):846–53. https://doi.org/10.1038/nm.3915.
Ellsworth RE, Blackburn HL, Shriver CD, Soon-Shiong P, Ellsworth DL. Molecular heterogeneity in breast cancer: state of the science and implications for patient care. Semin Cell Dev Biol. 2017;64:65–72. https://doi.org/10.1016/J.SEMCDB.2016.08.025.
Turashvili G, Brogi E. Tumor heterogeneity in breast cancer. Tumor Heterog Breast Cancer Front Med. 2017;4:227. https://doi.org/10.3389/fmed.2017.00227.
Zhao N, Rosen JM. Breast cancer heterogeneity through the lens of single-cell analysis and spatial pathologies. Semin Cancer Biol. 2022;82:3–10. https://doi.org/10.1016/j.semcancer.2021.07.010.
Keller L, Pantel K. Unravelling tumour heterogeneity by single-cell profiling of circulating tumour cells. Nat Rev Cancer. 2019;19(10):553–67. https://doi.org/10.1038/s41568-019-0180-2.
Wu SZ, Al-Eryani G, Roden DL, et al. A single-cell and spatially resolved atlas of human breast cancers. Nat Genet. 2021;53(9):1334–47. https://doi.org/10.1038/s41588-021-00911-1.
Chung W, Eum HH, Lee HO, et al. Single-cell RNA-seq enables comprehensive tumour and immune cell profiling in primary breast cancer. Nat Commun. 2017. https://doi.org/10.1038/ncomms15081.
Guillen KP, Fujita M, Butterfield AJ, et al. A human breast cancer-derived xenograft and organoid platform for drug discovery and precision oncology. Nat Cancer. 2022;3(2):232–50. https://doi.org/10.1038/s43018-022-00337-6.
Larsen BM, Kannan M, Langer LF, et al. A pan-cancer organoid platform for precision medicine. Cell Rep. 2021;36(4):109429. https://doi.org/10.1016/j.celrep.2021.109429.
Zheng B, Fang L. Spatially resolved transcriptomics provide a new method for cancer research. J Exp Clin Cancer Res. 2022. https://doi.org/10.1186/s13046-022-02385-3.
Wu C, Pineda F, Hormuth DA, Karczmar GS, Yankeelov TE. Quantitative analysis of vascular properties derived from ultrafast DCE-MRI to discriminate malignant and benign breast tumors. Magn Reson Med. 2019;81(3):2147–60. https://doi.org/10.1002/MRM.27529.
Hormuth DA, Al Feghali KA, Elliott AM, Yankeelov TE, Chung C. Image-based personalization of computational models for predicting response of high-grade glioma to chemoradiation. Sci Rep. 2021;11(1):1–14. https://doi.org/10.1038/s41598-021-87887-4.
Rockne RC, Hawkins-Daarud A, Swanson KR, et al. The 2019 mathematical oncology roadmap. Phys Biol. 2019;16(4):041005. https://doi.org/10.1088/1478-3975/AB1A09.
Yankeelov TE, Atuegwu N, Hormuth D, et al. Clinically relevant modeling of tumor growth and treatment response. Sci Transl Med. 2013. https://doi.org/10.1126/SCITRANSLMED.3005686.
Fertig EJ, Jaffee EM, Macklin P, Stearns V, Wang C. Forecasting cancer: from precision to predictive medicine. Med. 2021;2(9):1004–10. https://doi.org/10.1016/J.MEDJ.2021.08.007.
Cain EH, Saha A, Harowicz MR, Marks JR, Marcom PK, Mazurowski MA. Multivariate machine learning models for prediction of pathologic response to neoadjuvant therapy in breast cancer using MRI features: a study using an independent validation set. Breast Cancer Res Treat. 2018;173(2):455–63. https://doi.org/10.1007/S10549-018-4990-9.
Joo S, Ko ES, Kwon S, et al. Multimodal deep learning models for the prediction of pathologic response to neoadjuvant chemotherapy in breast cancer. Sci Rep. 2021;11(1):1–8. https://doi.org/10.1038/s41598-021-98408-8.
Alexander BM, Othus M, Caglar HB, Allen AM. Tumor volume is a prognostic factor in non-small-cell lung cancer treated with chemoradiotherapy. Int J Radiat Oncol. 2011;79(5):1381–7. https://doi.org/10.1016/J.IJROBP.2009.12.060.
Higuchi T, Fujimoto Y, Ozawa H, et al. Significance of metabolic tumor volume at baseline and reduction of mean standardized uptake value in 18 F-FDG-PET/CT imaging for predicting pathological complete response in breast cancers treated with preoperative chemotherapy. Ann Surg Oncol. 2019;26(7):2175–83. https://doi.org/10.1245/S10434-019-07325-8.
Henderson SA, Gowdh NM, Purdie CA, et al. Breast cancer: influence of tumour volume estimation method at MRI on prediction of pathological response to neoadjuvant chemotherapy. Br J Radiol. 2018. https://doi.org/10.1259/BJR.20180123.
Partridge SC, Gibbs JE, Lu Y, et al. MRI measurements of breast tumor volume predict response to neoadjuvant chemotherapy and recurrence-free survival. 2012;184(6):1774–1781. https://doi.org/10.2214/AJR.184.6.01841774.
Jafri NF, Newitt DC, Kornak J, Esserman LJ, Joe BN, Hylton NM. Optimized breast MRI functional tumor volume as a biomarker of recurrence-free survival following neoadjuvant chemotherapy. J Magn Reson Imaging. 2014;40(2):476–82. https://doi.org/10.1002/JMRI.24351.
Lee P, Bazan JG, Lavori PW, et al. Metabolic tumor volume is an independent prognostic factor in patients treated definitively for non–small-cell lung cancer. Clin Lung Cancer. 2012;13(1):52–8. https://doi.org/10.1016/J.CLLC.2011.05.001.
Kumasaka S, Nakajima T, Arisaka Y, et al. Prognostic value of metabolic tumor volume of pretreatment 18F-FAMT PET/CT in non-small cell lung Cancer. BMC Med Imaging. 2018;18(1):1–8. https://doi.org/10.1186/S12880-018-0292-2.
Seifert R, Kessel K, Schlack K, et al. PSMA PET total tumor volume predicts outcome of patients with advanced prostate cancer receiving [177Lu]Lu-PSMA-617 radioligand therapy in a bicentric analysis. Eur J Nucl Med Mol Imaging. 2021;48(4):1200–10. https://doi.org/10.1007/S00259-020-05040-1.
Seifert R, Herrmann K, Kleesiek J, et al. Semiautomatically quantified tumor volume using 68Ga-PSMA-11 PET as a biomarker for survival in patients with advanced prostate cancer. J Nucl Med. 2020;61(12):1786–92. https://doi.org/10.2967/JNUMED.120.242057.
Takahashi S, Go T, Anada M, et al. Correlation of pathological complete response with tumor volume reduction during neoadjuvant chemoradiotherapy in lung cancer. Anticancer Res. 2020;40:4327–30. https://doi.org/10.21873/anticanres.14435.
Pellegrino S, Fonti R, Mazziotti E, et al. 2019 Total metabolic tumor volume by 18F-FDG PET/CT for the prediction of outcome in patients with non-small cell lung cancer. Ann Nucl Med. 2019;33(12):937–44. https://doi.org/10.1007/S12149-019-01407-Z.
Lutsyk M, Awawda M, Gourevich K, Ben YR. Tumor volume as predictor of pathologic complete response following neoadjuvant chemoradiation in locally advanced rectal cancer. Am J Clin Oncol. 2021;44(9):482–6. https://doi.org/10.1097/COC.0000000000000846.
Cortazar P, Zhang L, Untch M, et al. Pathological complete response and long-term clinical benefit in breast cancer: the CTNeoBC pooled analysis. Lancet. 2014;384(9938):164–72. https://doi.org/10.1016/S0140-6736(13)62422-8.
Spring LM, Fell G, Arfe A, et al. Pathologic complete response after neoadjuvant chemotherapy and impact on breast cancer recurrence and survival: a comprehensive meta-analysis. Clin Cancer Res. 2020;26(12):2838–48. https://doi.org/10.1158/1078-0432.CCR-19-3492.
Yau C, Osdoit M, van der Noordaa M, et al. Residual cancer burden after neoadjuvant chemotherapy and long-term survival outcomes in breast cancer: a multicentre pooled analysis of 5161 patients. Lancet Oncol. 2022;23(1):149–60. https://doi.org/10.1016/S1470-2045(21)00589-1.
Valero V, Chavez-Macgregor M, Symmans WF, et al. Measurement of residual breast cancer burden to predict survival Aer neoadjuvant chemotherapy measurement of residual breast cancer burden to predict survival after neoadjuvant chemotherapy. J Clin Oncol. 2007;25:4414–22. https://doi.org/10.1200/JCO.2007.10.6823.
Nishino M. Tumor response assessment for precision cancer therapy: response evaluation criteria in solid tumors and beyond. Am Soc Clin Oncol Educ Book. 2018;38:1019–29. https://doi.org/10.1200/edbk_201441.
Samiei S, Van Nijnatten TJA, De Munck L, et al. Correlation between pathologic complete response in the breast and absence of axillary lymph node metastases after neoadjuvant systemic therapy. Ann Surg. 2020;271(3):574–80. https://doi.org/10.1097/SLA.0000000000003126.
Fayanju OM, Ren Y, Thomas SM, et al. The clinical significance of breast-only and node-only pathologic complete response (pCR) after neoadjuvant chemotherapy (NACT): a review of 20,000 breast cancer patients in the national cancer data base (NCDB). Ann Surg. 2018;268(4):591. https://doi.org/10.1097/SLA.0000000000002953.
Ethics approval and consent to participate
De-identified patient data were used for the study; thus, an institutional review board (IRB) exemption was granted.
Consent for publication
J.R. Peterson, J.A. Cole J.R.Pfeiffer, T. Pandey and A. K. Antony report ownership interest (stocks, stock options, patent or other intellectual property or other ownership interest excluding diversified mutual funds, receipt of intellectual property rights / patent holder, salary) from SimBioSys, Inc. G. Norris Y., Zhang D.A., Lopez-Ramos, and M. Biancalana report salary from SimBioSys, Inc. The other authors declare no potential conflicts of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Peterson, J.R., Cole, J.A., Pfeiffer, J.R. et al. Novel computational biology modeling system can accurately forecast response to neoadjuvant therapy in early breast cancer. Breast Cancer Res 25, 54 (2023). https://doi.org/10.1186/s13058-023-01654-z