Establishing a living biobank of patient-derived organoids of intraductal papillary mucinous neoplasms of the pancreas.

Pancreatic cancer (PaCa) is the third leading cause of cancer-related deaths in the United States. There is an unmet need to develop strategies to detect PaCa at an early, operable stage and prevent its progression. Intraductal papillary mucinous neoplasms (IPMNs) are cystic PaCa precursors that comprise nearly 50% of pancreatic cysts detected incidentally via cross-sectional imaging. Since IPMNs can progress from low- and moderate-grade dysplasia to high-grade dysplasia and invasion, the study of these lesions offers a prime opportunity to develop early detection and prevention strategies. Organoids are an ideal preclinical platform to study IPMNs, and the objective of the current investigation was to establish a living biobank of patient-derived organoids (PDO) from IPMNs. IPMN tumors and adjacent normal pancreatic tissues were successfully harvested from 15 patients with IPMNs undergoing pancreatic surgical resection at Moffitt Cancer Center & Research Institute (Tampa, FL) between May of 2017 and March of 2019. Organoid cultures were also generated from cryopreserved tissues. Organoid count and size were determined over time by both Image-Pro Premier 3D Version 9.1 digital platform and Matlab application of a Circular Hough Transform algorithm, and histologic and genomic characterization of a subset of the organoids was performed using immunohistochemistry and targeted sequencing, respectively. The success rates for organoid generation from IPMN tumor and adjacent normal pancreatic tissues were 81% and 87%, respectively. IPMN organoids derived from different epithelial subtypes showed different morphologies in vitro, and organoids recapitulated histologic and genomic characteristics of the parental IPMN tumor. In summary, this preclinical model has the potential to provide new opportunities to unveil mechanisms of IPMN progression to invasion and to shed insight into novel biomarkers for early detection and targets for chemoprevention.

The Image Biomarker Standardization Initiative: Standardized Quantitative Radiomics for High-Throughput Image-based Phenotyping.

Background Radiomic features may quantify characteristics present in medical imaging. However, the lack of standardized definitions and validated reference values have hampered clinical use. Purpose To standardize a set of 174 radiomic features. Materials and Methods Radiomic features were assessed in three phases. In phase I, 487 features were derived from the basic set of 174 features. Twenty-five research teams with unique radiomics software implementations computed feature values directly from a digital phantom, without any additional image processing. In phase II, 15 teams computed values for 1347 derived features using a CT image of a patient with lung cancer and predefined image processing configurations. In both phases, consensus among the teams on the validity of tentative reference values was measured through the frequency of the modal value and classified as follows: less than three matches, weak; three to five matches, moderate; six to nine matches, strong; 10 or more matches, very strong. In the final phase (phase III), a public data set of multimodality images (CT, fluorine 18 fluorodeoxyglucose PET, and T1-weighted MRI) from 51 patients with soft-tissue sarcoma was used to prospectively assess reproducibility of standardized features. Results Consensus on reference values was initially weak for 232 of 302 features (76.8%) at phase I and 703 of 1075 features (65.4%) at phase II. At the final iteration, weak consensus remained for only two of 487 features (0.4%) at phase I and 19 of 1347 features (1.4%) at phase II. Strong or better consensus was achieved for 463 of 487 features (95.1%) at phase I and 1220 of 1347 features (90.6%) at phase II. Overall, 169 of 174 features were standardized in the first two phases. In the final validation phase (phase III), most of the 169 standardized features could be excellently reproduced (166 with CT; 164 with PET; and 164 with MRI). Conclusion A set of 169 radiomics features was standardized, which enabled verification and calibration of different radiomics software. © RSNA, 2020 Online supplemental material is available for this article. See also the editorial by Kuhl and Truhn in this issue.

Prediction of radiologic outcome-optimized dose plans and post-treatment magnetic resonance images: A proof-of-concept study in breast cancer brain metastases treated with stereotactic radiosurgery.

Background and purpose: Information in multiparametric Magnetic Resonance (mpMR) images is relatable to voxel-level tumor response to Radiation Treatment (RT). We have investigated a deep learning framework to predict (i) post-treatment mpMR images from pre-treatment mpMR images and the dose map ("forward models"), and, (ii) the RT dose map that will produce prescribed changes within the Gross Tumor Volume (GTV) on post-treatment mpMR images ("inverse model"), in Breast Cancer Metastases to the Brain (BCMB) treated with Stereotactic Radiosurgery (SRS). Materials and methods: Local outcomes, planning computed tomography (CT) images, dose maps, and pre-treatment and post-treatment Apparent Diffusion Coefficient of water (ADC) maps, T1-weighted unenhanced (T1w) and contrast-enhanced (T1wCE), T2-weighted (T2w) and Fluid-Attenuated Inversion Recovery (FLAIR) mpMR images were curated from 39 BCMB patients. mpMR images were co-registered to the planning CT and intensity-calibrated. A 2D pix2pix architecture was used to train 5 forward models (ADC, T2w, FLAIR, T1w, T1wCE) and 1 inverse model on 1940 slices from 18 BCMB patients, and tested on 437 slices from another 9 BCMB patients. Results: Root Mean Square Percent Error (RMSPE) within the GTV between predicted and ground-truth post-RT images for the 5 forward models, in 136 test slices containing GTV, were (mean ± SD) 0.12 ± 0.044 (ADC), 0.14 ± 0.066 (T2w), 0.08 ± 0.038 (T1w), 0.13 ± 0.058 (T1wCE), and 0.09 ± 0.056 (FLAIR). RMSPE within the GTV on the same 136 test slices, between the predicted and ground-truth dose maps, was 0.37 ± 0.20 for the inverse model. Conclusions: A deep learning-based approach for radiologic outcome-optimized dose planning in SRS of BCMB has been demonstrated.

Interactions between ploidy and resource availability shape clonal interference at initiation and recurrence of glioblastoma.

Glioblastoma (GBM) is the most aggressive form of primary brain tumor. Complete surgical resection of GBM is almost impossible due to the infiltrative nature of the cancer. While no evidence for recent selection events have been found after diagnosis, the selective forces that govern gliomagenesis are strong, shaping the tumor's cell composition during the initial progression to malignancy with late consequences for invasiveness and therapy response. We present a mathematical model that simulates the growth and invasion of a glioma, given its ploidy level and the nature of its brain tissue micro-environment (TME), and use it to make inferences about GBM initiation and response to standard-of-care treatment. We approximate the spatial distribution of resource access in the TME through integration of in-silico modelling, multi-omics data and image analysis of primary and recurrent GBM. In the pre-malignant setting, our in-silico results suggest that low ploidy cancer cells are more resistant to starvation-induced cell death. In the malignant setting, between first and second surgery, simulated tumors with different ploidy compositions progressed at different rates. Whether higher ploidy predicted fast recurrence, however, depended on the TME. Historical data supports this dependence on TME resources, as shown by a significant correlation between the median glucose uptake rates in human tissues and the median ploidy of cancer types that arise in the respective tissues (Spearman r = -0.70; P = 0.026). Taken together our findings suggest that availability of metabolic substrates in the TME drives different cell fate decisions for cancer cells with different ploidy and shapes GBM disease initiation and relapse characteristics.

In Vivo Imaging of Rat Vascularity with FDG-Labeled Erythrocytes.

Microvascular disease is frequently found in major pathologies affecting vital organs, such as the brain, heart, and kidneys. While imaging modalities, such as ultrasound, computed tomography, single photon emission computed tomography, and magnetic resonance imaging, are widely used to visualize vascular abnormalities, the ability to non-invasively assess an organ's total vasculature, including microvasculature, is often limited or cumbersome. Previously, we have demonstrated proof of concept that non-invasive imaging of the total mouse vasculature can be achieved with 18F-fluorodeoxyglucose (18F-FDG)-labeled human erythrocytes and positron emission tomography/computerized tomography (PET/CT). In this work, we demonstrate that changes in the total vascular volume of the brain and left ventricular myocardium of normal rats can be seen after pharmacological vasodilation using 18F-FDG-labeled rat red blood cells (FDG RBCs) and microPET/CT imaging. FDG RBC PET imaging was also used to approximate the location of myocardial injury in a surgical myocardial infarction rat model. Finally, we show that FDG RBC PET imaging can detect relative differences in the degree of drug-induced intra-myocardial vasodilation between diabetic rats and normal controls. This FDG-labeled RBC PET imaging technique may thus be useful for assessing microvascular disease pathologies and characterizing pharmacological responses in the vascular bed of interest.

Deep-learning and MR images to target hypoxic habitats with evofosfamide in preclinical models of sarcoma.

Rationale: Hypoxic regions (habitats) within tumors are heterogeneously distributed and can be widely variant. Hypoxic habitats are generally pan-therapy resistant. For this reason, hypoxia-activated prodrugs (HAPs) have been developed to target these resistant volumes. The HAP evofosfamide (TH-302) has shown promise in preclinical and early clinical trials of sarcoma. However, in a phase III clinical trial of non-resectable soft tissue sarcomas, TH-302 did not improve survival in combination with doxorubicin (Dox), possibly due to a lack of patient stratification based on hypoxic status. Therefore, we used magnetic resonance imaging (MRI) to identify hypoxic habitats and non-invasively follow therapies response in sarcoma mouse models. Methods: We developed deep-learning (DL) models to identify hypoxia, using multiparametric MRI and co-registered histology, and monitored response to TH-302 in a patient-derived xenograft (PDX) of rhabdomyosarcoma and a syngeneic model of fibrosarcoma (radiation-induced fibrosarcoma, RIF-1). Results: A DL convolutional neural network showed strong correlations (>0.76) between the true hypoxia fraction in histology and the predicted hypoxia fraction in multiparametric MRI. TH-302 monotherapy or in combination with Dox delayed tumor growth and increased survival in the hypoxic PDX model (p<0.05), but not in the RIF-1 model, which had a lower volume of hypoxic habitats. Control studies showed that RIF-1 resistance was due to hypoxia and not other causes. Notably, PDX tumors developed resistance to TH-302 under prolonged treatment that was not due to a reduction in hypoxic volumes. Conclusion: Artificial intelligence analysis of pre-therapy MR images can predict hypoxia and subsequent response to HAPs. This approach can be used to monitor therapy response and adapt schedules to forestall the emergence of resistance.

Multi-window CT based Radiomic signatures in differentiating indolent versus aggressive lung cancers in the National Lung Screening Trial: a retrospective study.

BACKGROUND: We retrospectively evaluated the capability of radiomic features to predict tumor growth in lung cancer screening and compared the performance of multi-window radiomic features and single window radiomic features. METHODS: One hundred fifty lung nodules among 114 screen-detected, incident lung cancer patients from the National Lung Screening Trial (NLST) were investigated. Volume double time (VDT) was calculated as the difference between continuous two scans and used to define indolent and aggressive lung cancers. Lung nodules were semi-automatically segmented using lung and mediastinal windows separately, and subtracting the mediastinal window region from the lung window region generated the difference region. 364 radiomic features were separately exacted from nodules using the lung window, the mediastinal window and the difference region. Multivariable models were conducted to identify the most predictive features in predicting tumor growth. Clinical information was also obtained from the database. RESULTS: Based on our definition, 26% of the cases were indolent lung cancer. The tumor growth pattern could be predicted by radiomic models constructed using features obtained in the lung window, the difference region, and by combining features obtained in both the lung window and difference regions with areas under the receiver operator characteristic (AUROCs) of 0.799, 0.819, and 0.846, respectively. The multi-window feature model showed better performance compared to single window features (P < 0.001). Incorporating clinical factors into the multi-window feature models showed improvement, yielding an accuracy of 84.67% and AUROC of 0.855 for distinguishing indolent from aggressive disease. CONCLUSIONS: Multi-window CT based radiomics features are valuable predictors of indolent lung cancers and out performed single CT window setting. Combining clinical information improved predicting performance.

Multiparametric MRI and Coregistered Histology Identify Tumor Habitats in Breast Cancer Mouse Models.

It is well-recognized that solid tumors are genomically, anatomically, and physiologically heterogeneous. In general, more heterogeneous tumors have poorer outcomes, likely due to the increased probability of harboring therapy-resistant cells and regions. It is hypothesized that the genomic and physiologic heterogeneity are related, because physiologically distinct regions will exert variable selection pressures leading to the outgrowth of clones with variable genomic/proteomic profiles. To investigate this, methods must be in place to interrogate and define, at the microscopic scale, the cytotypes that exist within physiologically distinct subregions ("habitats") that are present at mesoscopic scales. MRI provides a noninvasive approach to interrogate physiologically distinct local environments, due to the biophysical principles that govern MRI signal generation. Here, we interrogate different physiologic parameters, such as perfusion, cell density, and edema, using multiparametric MRI (mpMRI). Signals from six different acquisition schema were combined voxel-by-voxel into four clusters identified using a Gaussian mixture model. These were compared with histologic and IHC characterizations of sections that were coregistered using MRI-guided 3D printed tumor molds. Specifically, we identified a specific set of MRI parameters to classify viable-normoxic, viable-hypoxic, nonviable-hypoxic, and nonviable-normoxic tissue types within orthotopic 4T1 and MDA-MB-231 breast tumors. This is the first coregistered study to show that mpMRI can be used to define physiologically distinct tumor habitats within breast tumor models. SIGNIFICANCE: This study demonstrates that noninvasive imaging metrics can be used to distinguish subregions within heterogeneous tumors with histopathologic correlation.