Patient-Specific, Mechanistic Models of Tumor Growth Incorporating Artificial Intelligence and Big Data.

Despite the remarkable advances in cancer diagnosis, treatment, and management over the past decade, malignant tumors remain a major public health problem. Further progress in combating cancer may be enabled by personalizing the delivery of therapies according to the predicted response for each individual patient. The design of personalized therapies requires the integration of patient-specific information with an appropriate mathematical model of tumor response. A fundamental barrier to realizing this paradigm is the current lack of a rigorous yet practical mathematical theory of tumor initiation, development, invasion, and response to therapy. We begin this review with an overview of different approaches to modeling tumor growth and treatment, including mechanistic as well as data-driven models based on big data and artificial intelligence. We then present illustrative examples of mathematical models manifesting their utility and discuss the limitations of stand-alone mechanistic and data-driven models. We then discuss the potential of mechanistic models for not only predicting but also optimizing response to therapy on a patient-specific basis. We describe current efforts and future possibilities to integrate mechanistic and data-driven models. We conclude by proposing five fundamental challenges that must be addressed to fully realize personalized care for cancer patients driven by computational models.

A Pilot Study on Patient-specific Computational Forecasting of Prostate Cancer Growth during Active Surveillance Using an Imaging-informed Biomechanistic Model.

Active surveillance (AS) is a suitable management option for newly diagnosed prostate cancer, which usually presents low to intermediate clinical risk. Patients enrolled in AS have their tumor monitored via longitudinal multiparametric MRI (mpMRI), PSA tests, and biopsies. Hence, treatment is prescribed when these tests identify progression to higher-risk prostate cancer. However, current AS protocols rely on detecting tumor progression through direct observation according to population-based monitoring strategies. This approach limits the design of patient-specific AS plans and may delay the detection of tumor progression. Here, we present a pilot study to address these issues by leveraging personalized computational predictions of prostate cancer growth. Our forecasts are obtained with a spatiotemporal biomechanistic model informed by patient-specific longitudinal mpMRI data (T2-weighted MRI and apparent diffusion coefficient maps from diffusion-weighted MRI). Our results show that our technology can represent and forecast the global tumor burden for individual patients, achieving concordance correlation coefficients from 0.93 to 0.99 across our cohort (n = 7). In addition, we identify a model-based biomarker of higher-risk prostate cancer: the mean proliferation activity of the tumor (P = 0.041). Using logistic regression, we construct a prostate cancer risk classifier based on this biomarker that achieves an area under the ROC curve of 0.83. We further show that coupling our tumor forecasts with this prostate cancer risk classifier enables the early identification of prostate cancer progression to higher-risk disease by more than 1 year. Thus, we posit that our predictive technology constitutes a promising clinical decision-making tool to design personalized AS plans for patients with prostate cancer. SIGNIFICANCE: Personalization of a biomechanistic model of prostate cancer with mpMRI data enables the prediction of tumor progression, thereby showing promise to guide clinical decision-making during AS for each individual patient.

Toward mechanistic medical digital twins: some use cases in immunology.

A fundamental challenge for personalized medicine is to capture enough of the complexity of an individual patient to determine an optimal way to keep them healthy or restore their health. This will require personalized computational models of sufficient resolution and with enough mechanistic information to provide actionable information to the clinician. Such personalized models are increasingly referred to as medical digital twins. Digital twin technology for health applications is still in its infancy, and extensive research and development is required. This article focuses on several projects in different stages of development that can lead to specific-and practical-medical digital twins or digital twin modeling platforms. It emerged from a two-day forum on problems related to medical digital twins, particularly those involving an immune system component. Open access video recordings of the forum discussions are available.

Systematic evaluation of MRI-based characterization of tumor-associated vascular morphology and hemodynamics via a dynamic digital phantom.

Purpose: Validation of quantitative imaging biomarkers is a challenging task, due to the difficulty in measuring the ground truth of the target biological process. A digital phantom-based framework is established to systematically validate the quantitative characterization of tumor-associated vascular morphology and hemodynamics based on dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI). Approach: A digital phantom is employed to provide a ground-truth vascular system within which 45 synthetic tumors are simulated. Morphological analysis is performed on high-spatial resolution DCE-MRI data (spatial/temporal resolution = 30 to 300 µm/60 s) to determine the accuracy of locating the arterial inputs of tumor-associated vessels (TAVs). Hemodynamic analysis is then performed on the combination of high-spatial resolution and high-temporal resolution (spatial/temporal resolution = 60 to 300 µm/1 to 10 s) DCE-MRI data, determining the accuracy of estimating tumor-associated blood pressure, vascular extraction rate, interstitial pressure, and interstitial flow velocity. Results: The observed effects of acquisition settings demonstrate that, when optimizing the DCE-MRI protocol for the morphological analysis, increasing the spatial resolution is helpful but not necessary, as the location and arterial input of TAVs can be recovered with high accuracy even with the lowest investigated spatial resolution. When optimizing the DCE-MRI protocol for hemodynamic analysis, increasing the spatial resolution of the images used for vessel segmentation is essential, and the spatial and temporal resolutions of the images used for the kinetic parameter fitting require simultaneous optimization. Conclusion: An in silico validation framework was generated to systematically quantify the effects of image acquisition settings on the ability to accurately estimate tumor-associated characteristics.

Forum on immune digital twins: a meeting report.

Medical digital twins are computational models of human biology relevant to a given medical condition, which are tailored to an individual patient, thereby predicting the course of disease and individualized treatments, an important goal of personalized medicine. The immune system, which has a central role in many diseases, is highly heterogeneous between individuals, and thus poses a major challenge for this technology. In February 2023, an international group of experts convened for two days to discuss these challenges related to immune digital twins. The group consisted of clinicians, immunologists, biologists, and mathematical modelers, representative of the interdisciplinary nature of medical digital twin development. A video recording of the entire event is available. This paper presents a synopsis of the discussions, brief descriptions of ongoing digital twin projects at different stages of progress. It also proposes a 5-year action plan for further developing this technology. The main recommendations are to identify and pursue a small number of promising use cases, to develop stimulation-specific assays of immune function in a clinical setting, and to develop a database of existing computational immune models, as well as advanced modeling technology and infrastructure.

A mathematical model for predicting the spatiotemporal response of breast cancer cells treated with doxorubicin.

Tumor heterogeneity contributes significantly to chemoresistance, a leading cause of treatment failure. To better personalize therapies, it is essential to develop tools capable of identifying and predicting intra- and inter-tumor heterogeneities. Biology-inspired mathematical models are capable of attacking this problem, but tumor heterogeneity is often overlooked in in-vivo modeling studies, while phenotypic considerations capturing spatial dynamics are not typically included in in-vitro modeling studies. We present a data assimilation-prediction pipeline with a two-phenotype model that includes a spatiotemporal component to characterize and predict the evolution of in-vitro breast cancer cells and their heterogeneous response to chemotherapy. Our model assumes that the cells can be divided into two subpopulations: surviving cells unaffected by the treatment, and irreversibly damaged cells undergoing treatment-induced death. MCF7 breast cancer cells were previously cultivated in wells for up to 1000 hours, treated with various concentrations of doxorubicin and imaged with time-resolved microscopy to record spatiotemporally-resolved cell count data. Images were used to generate cell density maps. Treatment response predictions were initialized by a training set and updated by weekly measurements. Our mathematical model successfully calibrated the spatiotemporal cell growth dynamics, achieving median [range] concordance correlation coefficients of > .99 [.88, >.99] and .73 [.58, .85] across the whole well and individual pixels, respectively. Our proposed data assimilation-prediction approach achieved values of .97 [.44, >.99] and .69 [.35, .79] for the whole well and individual pixels, respectively. Thus, our model can capture and predict the spatiotemporal dynamics of MCF7 cells treated with doxorubicin in an in-vitro setting.

Breast Multiparametric MRI for Prediction of Neoadjuvant Chemotherapy Response in Breast Cancer: The BMMR2 Challenge.

Purpose To describe the design, conduct, and results of the Breast Multiparametric MRI for prediction of neoadjuvant chemotherapy Response (BMMR2) challenge. Materials and Methods The BMMR2 computational challenge opened on May 28, 2021, and closed on December 21, 2021. The goal of the challenge was to identify image-based markers derived from multiparametric breast MRI, including diffusion-weighted imaging (DWI) and dynamic contrast-enhanced (DCE) MRI, along with clinical data for predicting pathologic complete response (pCR) following neoadjuvant treatment. Data included 573 breast MRI studies from 191 women (mean age [±SD], 48.9 years ± 10.56) in the I-SPY 2/American College of Radiology Imaging Network (ACRIN) 6698 trial (ClinicalTrials.gov: NCT01042379). The challenge cohort was split into training (60%) and test (40%) sets, with teams blinded to test set pCR outcomes. Prediction performance was evaluated by area under the receiver operating characteristic curve (AUC) and compared with the benchmark established from the ACRIN 6698 primary analysis. Results Eight teams submitted final predictions. Entries from three teams had point estimators of AUC that were higher than the benchmark performance (AUC, 0.782 [95% CI: 0.670, 0.893], with AUCs of 0.803 [95% CI: 0.702, 0.904], 0.838 [95% CI: 0.748, 0.928], and 0.840 [95% CI: 0.748, 0.932]). A variety of approaches were used, ranging from extraction of individual features to deep learning and artificial intelligence methods, incorporating DCE and DWI alone or in combination. Conclusion The BMMR2 challenge identified several models with high predictive performance, which may further expand the value of multiparametric breast MRI as an early marker of treatment response. Clinical trial registration no. NCT01042379 Keywords: MRI, Breast, Tumor Response Supplemental material is available for this article. © RSNA, 2024.

Designing clinical trials for patients who are not average.

The heterogeneity inherent in cancer means that even a successful clinical trial merely results in a therapeutic regimen that achieves, on average, a positive result only in a subset of patients. The only way to optimize an intervention for an individual patient is to reframe their treatment as their own, personalized trial. Toward this goal, we formulate a computational framework for performing personalized trials that rely on four mathematical techniques. First, mathematical models that can be calibrated with patient-specific data to make accurate predictions of response. Second, digital twins built on these models capable of simulating the effects of interventions. Third, optimal control theory applied to the digital twins to optimize outcomes. Fourth, data assimilation to continually update and refine predictions in response to therapeutic interventions. In this perspective, we describe each of these techniques, quantify their "state of readiness", and identify use cases for personalized clinical trials.

The ISMRM Open Science Initiative for Perfusion Imaging (OSIPI): Results from the OSIPI-Dynamic Contrast-Enhanced challenge.

PURPOSE: K trans $$ {K}^{\mathrm{trans}} $$ has often been proposed as a quantitative imaging biomarker for diagnosis, prognosis, and treatment response assessment for various tumors. None of the many software tools for K trans $$ {K}^{\mathrm{trans}} $$ quantification are standardized. The ISMRM Open Science Initiative for Perfusion Imaging-Dynamic Contrast-Enhanced (OSIPI-DCE) challenge was designed to benchmark methods to better help the efforts to standardize K trans $$ {K}^{\mathrm{trans}} $$ measurement. METHODS: A framework was created to evaluate K trans $$ {K}^{\mathrm{trans}} $$ values produced by DCE-MRI analysis pipelines to enable benchmarking. The perfusion MRI community was invited to apply their pipelines for K trans $$ {K}^{\mathrm{trans}} $$ quantification in glioblastoma from clinical and synthetic patients. Submissions were required to include the entrants' K trans $$ {K}^{\mathrm{trans}} $$ values, the applied software, and a standard operating procedure. These were evaluated using the proposed OSIP I gold $$ \mathrm{OSIP}{\mathrm{I}}_{\mathrm{gold}} $$ score defined with accuracy, repeatability, and reproducibility components. RESULTS: Across the 10 received submissions, the OSIP I gold $$ \mathrm{OSIP}{\mathrm{I}}_{\mathrm{gold}} $$ score ranged from 28% to 78% with a 59% median. The accuracy, repeatability, and reproducibility scores ranged from 0.54 to 0.92, 0.64 to 0.86, and 0.65 to 1.00, respectively (0-1 = lowest-highest). Manual arterial input function selection markedly affected the reproducibility and showed greater variability in K trans $$ {K}^{\mathrm{trans}} $$ analysis than automated methods. Furthermore, provision of a detailed standard operating procedure was critical for higher reproducibility. CONCLUSIONS: This study reports results from the OSIPI-DCE challenge and highlights the high inter-software variability within K trans $$ {K}^{\mathrm{trans}} $$ estimation, providing a framework for ongoing benchmarking against the scores presented. Through this challenge, the participating teams were ranked based on the performance of their software tools in the particular setting of this challenge. In a real-world clinical setting, many of these tools may perform differently with different benchmarking methodology.

Vascularized Hepatocellular Carcinoma on a Chip to Control Chemoresistance through Cirrhosis, Inflammation and Metabolic Activity.

Understanding the effects of inflammation and cirrhosis on the regulation of drug metabolism during the progression of hepatocellular carcinoma (HCC) is critical for developing patient-specific treatment strategies. In this work, we created novel three-dimensional vascularized HCC-on-a-chips (HCCoC), composed of HCC, endothelial, stellate, and Kupffer cells tuned to mimic normal or cirrhotic liver stiffness. HCC inflammation was controlled by tuning Kupffer macrophage numbers, and the impact of cytochrome P450-3A4 (CYP3A4) was investigated by culturing HepG2 HCC cells transfected with CYP3A4 to upregulate expression from baseline. This model allowed for the simulation of chemotherapeutic delivery methods such as intravenous injection and transcatheter arterial chemoembolization (TACE). We showed that upregulation of metabolic activity, incorporation of cirrhosis and inflammation, increase vascular permeability due to upregulated inflammatory cytokines leading to significant variability in chemotherapeutic treatment efficacy. Specifically, we show that further modulation of CYP3A4 activity of HCC cells by TACE delivery of doxorubicin provides an additional improvement to treatment response and reduces chemotherapy-associated endothelial porosity increase. The HCCoCs were shown to have utility in uncovering the impact of the tumor microenvironment (TME) during cancer progression on vascular properties, tumor response to therapeutics, and drug delivery strategies.

Predicting response to combination evofosfamide and immunotherapy under hypoxic conditions in murine models of colon cancer.

The goal of this study is to develop a mathematical model that captures the interaction between evofosfamide, immunotherapy, and the hypoxic landscape of the tumor in the treatment of tumors. Recently, we showed that evofosfamide, a hypoxia-activated prodrug, can synergistically improve treatment outcomes when combined with immunotherapy, while evofosfamide alone showed no effects in an in vivo syngeneic model of colorectal cancer. However, the mechanisms behind the interaction between the tumor microenvironment in the context of oxygenation (hypoxic, normoxic), immunotherapy, and tumor cells are not fully understood. To begin to understand this issue, we develop a system of ordinary differential equations to simulate the growth and decline of tumors and their vascularization (oxygenation) in response to treatment with evofosfamide and immunotherapy (6 combinations of scenarios). The model is calibrated to data from in vivo experiments on mice implanted with colon adenocarcinoma cells and longitudinally imaged with [18F]-fluoromisonidazole ([18F]FMISO) positron emission tomography (PET) to quantify hypoxia. The results show that evofosfamide is able to rescue the immune response and sensitize hypoxic tumors to immunotherapy. In the hypoxic scenario, evofosfamide reduces tumor burden by $ 45.07 \pm 2.55 $%, compared to immunotherapy alone, as measured by tumor volume. The model accurately predicts the temporal evolution of five different treatment scenarios, including control, hypoxic tumors that received immunotherapy, normoxic tumors that received immunotherapy, evofosfamide alone, and hypoxic tumors that received combination immunotherapy and evofosfamide. The average concordance correlation coefficient (CCC) between predicted and observed tumor volume is $ 0.86 \pm 0.05 $. Interestingly, the model values to fit those five treatment arms was unable to accurately predict the response of normoxic tumors to combination evofosfamide and immunotherapy (CCC = $ -0.064 \pm 0.003 $). However, guided by the sensitivity analysis to rank the most influential parameters on the tumor volume, we found that increasing the tumor death rate due to immunotherapy by a factor of $ 18.6 \pm 9.3 $ increases CCC of $ 0.981 \pm 0.001 $. To the best of our knowledge, this is the first study to mathematically predict and describe the increased efficacy of immunotherapy following evofosfamide.

Investigating tumor-host response dynamics in preclinical immunotherapy experiments using a stepwise mathematical modeling strategy.

Immunotherapies such as checkpoint blockade to PD1 and CTLA4 can have varied effects on individual tumors. To quantify the successes and failures of these therapeutics, we developed a stepwise mathematical modeling strategy and applied it to mouse models of colorectal and breast cancer that displayed a range of therapeutic responses. Using longitudinal tumor volume data, an exponential growth model was utilized to designate response groups for each tumor type. The exponential growth model was then extended to describe the dynamics of the quality of vasculature in the tumors via [18F] fluoromisonidazole (FMISO)-positron emission tomography (PET) data estimating tumor hypoxia over time. By calibrating the mathematical system to the PET data, several biological drivers of the observed deterioration of the vasculature were quantified. The mathematical model was then further expanded to explicitly include both the immune response and drug dosing, so that model simulations are able to systematically investigate biological hypotheses about immunotherapy failure and to generate experimentally testable predictions of immune response. The modeling results suggest elevated immune response fractions (> 30 %) in tumors unresponsive to immunotherapy is due to a functional immune response that wanes over time. This experimental-mathematical approach provides a means to evaluate dynamics of the system that could not have been explored using the data alone, including tumor aggressiveness, immune exhaustion, and immune cell functionality.

Predictive digital twin for optimizing patient-specific radiotherapy regimens under uncertainty in high-grade gliomas.

We develop a methodology to create data-driven predictive digital twins for optimal risk-aware clinical decision-making. We illustrate the methodology as an enabler for an anticipatory personalized treatment that accounts for uncertainties in the underlying tumor biology in high-grade gliomas, where heterogeneity in the response to standard-of-care (SOC) radiotherapy contributes to sub-optimal patient outcomes. The digital twin is initialized through prior distributions derived from population-level clinical data in the literature for a mechanistic model's parameters. Then the digital twin is personalized using Bayesian model calibration for assimilating patient-specific magnetic resonance imaging data. The calibrated digital twin is used to propose optimal radiotherapy treatment regimens by solving a multi-objective risk-based optimization under uncertainty problem. The solution leads to a suite of patient-specific optimal radiotherapy treatment regimens exhibiting varying levels of trade-off between the two competing clinical objectives: (i) maximizing tumor control (characterized by minimizing the risk of tumor volume growth) and (ii) minimizing the toxicity from radiotherapy. The proposed digital twin framework is illustrated by generating an in silico cohort of 100 patients with high-grade glioma growth and response properties typically observed in the literature. For the same total radiation dose as the SOC, the personalized treatment regimens lead to median increase in tumor time to progression of around six days. Alternatively, for the same level of tumor control as the SOC, the digital twin provides optimal treatment options that lead to a median reduction in radiation dose by 16.7% (10 Gy) compared to SOC total dose of 60 Gy. The range of optimal solutions also provide options with increased doses for patients with aggressive cancer, where SOC does not lead to sufficient tumor control.

Patient-specific, mechanistic models of tumor growth incorporating artificial intelligence and big data.

Despite the remarkable advances in cancer diagnosis, treatment, and management that have occurred over the past decade, malignant tumors remain a major public health problem. Further progress in combating cancer may be enabled by personalizing the delivery of therapies according to the predicted response for each individual patient. The design of personalized therapies requires patient-specific information integrated into an appropriate mathematical model of tumor response. A fundamental barrier to realizing this paradigm is the current lack of a rigorous, yet practical, mathematical theory of tumor initiation, development, invasion, and response to therapy. In this review, we begin by providing an overview of different approaches to modeling tumor growth and treatment, including mechanistic as well as data-driven models based on ``big data" and artificial intelligence. Next, we present illustrative examples of mathematical models manifesting their utility and discussing the limitations of stand-alone mechanistic and data-driven models. We further discuss the potential of mechanistic models for not only predicting, but also optimizing response to therapy on a patient-specific basis. We then discuss current efforts and future possibilities to integrate mechanistic and data-driven models. We conclude by proposing five fundamental challenges that must be addressed to fully realize personalized care for cancer patients driven by computational models.

Dynamic parameterization of a modified SEIRD model to analyze and forecast the dynamics of COVID-19 outbreaks in the United States.

The rapid spread of the numerous outbreaks of the coronavirus disease 2019 (COVID-19) pandemic has fueled interest in mathematical models designed to understand and predict infectious disease spread, with the ultimate goal of contributing to the decision making of public health authorities. Here, we propose a computational pipeline that dynamically parameterizes a modified SEIRD (susceptible-exposed-infected-recovered-deceased) model using standard daily series of COVID-19 cases and deaths, along with isolated estimates of population-level seroprevalence. We test our pipeline in five heavily impacted states of the US (New York, California, Florida, Illinois, and Texas) between March and August 2020, considering two scenarios with different calibration time horizons to assess the update in model performance as new epidemiologic data become available. Our results show a median normalized root mean squared error (NRMSE) of 2.38% and 4.28% in calibrating cumulative cases and deaths in the first scenario, and 2.41% and 2.30% when new data are assimilated in the second scenario, respectively. Then, 2-week (4-week) forecasts of the calibrated model resulted in median NRMSE of cumulative cases and deaths of 5.85% and 4.68% (8.60% and 17.94%) in the first scenario, and 1.86% and 1.93% (2.21% and 1.45%) in the second. Additionally, we show that our method provides significantly more accurate predictions of cases and deaths than a constant parameterization in the second scenario (p < 0.05). Thus, we posit that our methodology is a promising approach to analyze the dynamics of infectious disease outbreaks, and that our forecasts could contribute to designing effective pandemic-arresting public health policies. Supplementary Information: The online version contains supplementary material available at 10.1007/s00366-023-01816-9.

Comparing mechanism-based and machine learning models for predicting the effects of glucose accessibility on tumor cell proliferation.

Glucose plays a central role in tumor metabolism and development and is a target for novel therapeutics. To characterize the response of cancer cells to blockade of glucose uptake, we collected time-resolved microscopy data to track the growth of MDA-MB-231 breast cancer cells. We then developed a mechanism-based, mathematical model to predict how a glucose transporter (GLUT1) inhibitor (Cytochalasin B) influences the growth of the MDA-MB-231 cells by limiting access to glucose. The model includes a parameter describing dose dependent inhibition to quantify both the total glucose level in the system and the glucose level accessible to the tumor cells. Four common machine learning models were also used to predict tumor cell growth. Both the mechanism-based and machine learning models were trained and validated, and the prediction error was evaluated by the coefficient of determination (R2). The random forest model provided the highest accuracy predicting cell dynamics (R2 = 0.92), followed by the decision tree (R2 = 0.89), k-nearest-neighbor regression (R2 = 0.84), mechanism-based (R2 = 0.77), and linear regression model (R2 = 0.69). Thus, the mechanism-based model has a predictive capability comparable to machine learning models with the added benefit of elucidating biological mechanisms.

Co-Clinical Imaging Metadata Information (CIMI) for Cancer Research to Promote Open Science, Standardization, and Reproducibility in Preclinical Imaging.

Preclinical imaging is a critical component in translational research with significant complexities in workflow and site differences in deployment. Importantly, the National Cancer Institute's (NCI) precision medicine initiative emphasizes the use of translational co-clinical oncology models to address the biological and molecular bases of cancer prevention and treatment. The use of oncology models, such as patient-derived tumor xenografts (PDX) and genetically engineered mouse models (GEMMs), has ushered in an era of co-clinical trials by which preclinical studies can inform clinical trials and protocols, thus bridging the translational divide in cancer research. Similarly, preclinical imaging fills a translational gap as an enabling technology for translational imaging research. Unlike clinical imaging, where equipment manufacturers strive to meet standards in practice at clinical sites, standards are neither fully developed nor implemented in preclinical imaging. This fundamentally limits the collection and reporting of metadata to qualify preclinical imaging studies, thereby hindering open science and impacting the reproducibility of co-clinical imaging research. To begin to address these issues, the NCI co-clinical imaging research program (CIRP) conducted a survey to identify metadata requirements for reproducible quantitative co-clinical imaging. The enclosed consensus-based report summarizes co-clinical imaging metadata information (CIMI) to support quantitative co-clinical imaging research with broad implications for capturing co-clinical data, enabling interoperability and data sharing, as well as potentially leading to updates to the preclinical Digital Imaging and Communications in Medicine (DICOM) standard.

An Online Repository for Pre-Clinical Imaging Protocols (PIPs).

Providing method descriptions that are more detailed than currently available in typical peer reviewed journals has been identified as an actionable area for improvement. In the biochemical and cell biology space, this need has been met through the creation of new journals focused on detailed protocols and materials sourcing. However, this format is not well suited for capturing instrument validation, detailed imaging protocols, and extensive statistical analysis. Furthermore, the need for additional information must be counterbalanced by the additional time burden placed upon researchers who may be already overtasked. To address these competing issues, this white paper describes protocol templates for positron emission tomography (PET), X-ray computed tomography (CT), and magnetic resonance imaging (MRI) that can be leveraged by the broad community of quantitative imaging experts to write and self-publish protocols in protocols.io. Similar to the Structured Transparent Accessible Reproducible (STAR) or Journal of Visualized Experiments (JoVE) articles, authors are encouraged to publish peer reviewed papers and then to submit more detailed experimental protocols using this template to the online resource. Such protocols should be easy to use, readily accessible, readily searchable, considered open access, enable community feedback, editable, and citable by the author.

Toward Practical Integration of Omic and Imaging Data in Co-Clinical Trials.

Co-clinical trials are the concurrent or sequential evaluation of therapeutics in both patients clinically and patient-derived xenografts (PDX) pre-clinically, in a manner designed to match the pharmacokinetics and pharmacodynamics of the agent(s) used. The primary goal is to determine the degree to which PDX cohort responses recapitulate patient cohort responses at the phenotypic and molecular levels, such that pre-clinical and clinical trials can inform one another. A major issue is how to manage, integrate, and analyze the abundance of data generated across both spatial and temporal scales, as well as across species. To address this issue, we are developing MIRACCL (molecular and imaging response analysis of co-clinical trials), a web-based analytical tool. For prototyping, we simulated data for a co-clinical trial in "triple-negative" breast cancer (TNBC) by pairing pre- (T0) and on-treatment (T1) magnetic resonance imaging (MRI) from the I-SPY2 trial, as well as PDX-based T0 and T1 MRI. Baseline (T0) and on-treatment (T1) RNA expression data were also simulated for TNBC and PDX. Image features derived from both datasets were cross-referenced to omic data to evaluate MIRACCL functionality for correlating and displaying MRI-based changes in tumor size, vascularity, and cellularity with changes in mRNA expression as a function of treatment.