Evolutionary dynamics at the tumor edge reveal metabolic imaging biomarkers

Spatial interactions modulate tumor growth and immune infiltration

Lenia, a cellular automata framework used in artificial life, provides a natural setting to implement mathematical models of cancer incorporating features such as morphogenesis, homeostasis, motility, reproduction, growth, stimuli response, evolvability, and adaptation. Historically, agent-based models of cancer progression have been constructed with rules that govern birth, death and migration, with attempts to map local rules to emergent global growth dynamics. In contrast, Lenia provides a flexible framework for considering a spectrum of local (cell-scale) to global (tumor-scale) dynamics by defining an interaction kernel governing density-dependent growth dynamics. Lenia can recapitulate a range of cancer model classifications including local or global, deterministic or stochastic, non-spatial or spatial, single or multi-population, and off or on-lattice. Lenia is subsequently used to develop data-informed models of 1) single-population growth dynamics, 2) multi-population cell-cell competition models, and 3) cell migration or chemotaxis. Mathematical modeling provides important mechanistic insights. First, short-range interaction kernels provide a mechanism for tumor cell survival under conditions with strong Allee effects. Next, we find that asymmetric interaction tumor-immune kernels lead to poor immune response. Finally, modeling recapitulates immune-ECM interactions where patterns of collagen formation provide immune protection, indicated by an emergent inverse relationship between disease stage and immune coverage.

Preoperative Three-Dimensional Morphological Tumor Features Predict Microvascular Invasion in Hepatocellular Carcinoma

RATIONALE AND OBJECTIVES

The study was designed to evaluate microvascular invasion (MVI) using three-dimensional (3D) morphological indicators prior to surgery.

MATERIALS AND METHODS

This retrospective study included 156 patients with hepatocellular carcinoma (HCC) at our hospital from 2017 to 2018. Through thin-layer CT scanning and 3D reconstruction, the tumor surface inclination angles can be quantitatively analyzed to determine the surface irregularity rate (SIR), which serves as a comprehensive assessment method for tumor irregularity based on preoperative 3D morphological evaluation. Univariate and multivariate logistic regression analyses were employed to investigate the correlation with MVI.

RESULTS

The SIR was related to MVI (OR: 10.667, P < 0.001). Multivariate logistic regression analysis showed that the SIR was an independent risk factor for MVI. The area under the receiver operating characteristic curve (ROC) of prediction model composed of the morphological indicator SIR was 0.831 (95% confidence interval: 0.759-0.895).

CONCLUSION

The preoperative 3D morphological indicator SIR of a tumor is an accurate predictor of MVI, providing a valuable tool in clinical decision-making.

Integrated modeling of labile and glycated hemoglobin with glucose for enhanced diabetes detection and short-term monitoring

Metabolic biomarkers, particularly glycated hemoglobin and fasting plasma glucose, are pivotal in the diagnosis and control of diabetes mellitus. Despite their importance, they exhibit limitations in assessing short-term glucose variations. In this study, we propose labile hemoglobin as an additional biomarker, providing insightful perspectives into these fluctuations. By utilizing datasets from 40,652 retrospective general participants and conducting glucose tolerance tests on 60 prospective pediatric subjects, we explored the relationship between plasma glucose and labile hemoglobin. A mathematical model was developed to encapsulate short-term glucose kinetics in the pediatric group. Applying dimensionality reduction techniques, we successfully identified participant subclusters, facilitating the differentiation between diabetic and non-diabetic individuals. Intriguingly, by integrating labile hemoglobin measurements with plasma glucose values, we were able to predict the likelihood of diabetes in pediatric subjects, underscoring the potential of labile hemoglobin as a significant glycemic biomarker for diabetes research.

Promising Candidate Prognostic Biomarkers in [18F]FDG PET Images: Evaluation in Independent Cohorts of Non-Small Cell Lung Cancer Patients

The normalized distances from the hot spot of radiotracer uptake (SUVmax) to the tumor centroid (NHOC) and to the tumor perimeter (NHOP) have recently been suggested as novel PET features reflecting tumor aggressiveness. These biomarkers characterizing the shift of SUVmax toward the lesion edge during tumor progression have been shown to be prognostic factors in breast and non-small cell lung cancer (NSCLC) patients. We assessed the impact of imaging parameters on NHOC and NHOP, their complementarity to conventional PET features, and their prognostic value for advanced-NSCLC patients. Methods: This retrospective study investigated baseline [18F]FDG PET scans: cohort 1 included 99 NSCLC patients with no treatment-related inclusion criteria (robustness study); cohort 2 included 244 NSCLC patients (survival analysis) treated with targeted therapy (93), immunotherapy (63), or immunochemotherapy (88). Although 98% of patients had metastases, radiomic features including SUVs were extracted from the primary tumor only. NHOCs and NHOPs were computed using 2 approaches: the normalized distance from the localization of SUVmax or SUVpeak to the tumor centroid or perimeter. Bland-Altman analyses were performed to investigate the impact of both spatial resolution (comparing PET images with and without gaussian postfiltering) and image sampling (comparing 2 voxel sizes) on feature values. The correlation of NHOCs and NHOPs with other features was studied using Spearman correlation coefficients (r). The ability of NHOCs and NHOPs to predict overall survival (OS) was estimated using the Kaplan-Meier method. Results: In cohort 1, NHOC and NHOP features were more robust to image filtering and to resampling than were SUVs. The correlations were weak between NHOCs and NHOPs (r ≤ 0.45) and between NHOCs or NHOPs and any other radiomic features (r ≤ 0.60). In cohort 2, the patients with short OS demonstrated higher NHOCs and lower NHOPs than those with long OS. NHOCs significantly distinguished 2 survival profiles in patients treated with immunotherapy (log-rank test, P < 0.01), whereas NHOPs stratified patients regarding OS in the targeted therapy (P = 0.02) and immunotherapy (P < 0.01) subcohorts. Conclusion: Our findings suggest that even in advanced NSCLC patients, NHOC and NHOP features pertaining to the primary tumor have prognostic potential. Moreover, these features appeared to be robust with respect to imaging protocol parameters and complementary to other radiomic features and are now available in LIFEx software to be independently tested by others.

Spatial interactions modulate tumor growth and immune infiltration

Direct observation of immune cell trafficking patterns and tumor-immune interactions is unlikely in human tumors with currently available technology, but computational simulations based on clinical data can provide insight to test hypotheses. It is hypothesized that patterns of collagen formation evolve as a mechanism of immune escape, but the exact nature of the interaction between immune cells and collagen is poorly understood. Spatial data quantifying the degree of collagen fiber alignment in squamous cell carcinomas indicates that late stage disease is associated with highly aligned fibers. Here, we introduce a computational modeling framework (called Lenia) to discriminate between two hypotheses: immune cell migration that moves 1) parallel or 2) perpendicular to collagen fiber orientation. The modeling recapitulates immune-ECM interactions where collagen patterns provide immune protection, leading to an emergent inverse relationship between disease stage and immune coverage. We also illustrate the capabilities of Lenia to model the evolution of tumor progression and immune predation. Lenia provides a flexible framework for considering a spectrum of local (cell-scale) to global (tumor-scale) dynamics by defining a kernel cell-cell interaction function that governs tumor growth dynamics under immune predation with immune cell migration. Mathematical modeling provides important mechanistic insights into cell interactions. Short-range interaction kernels provide a mechanism for tumor cell survival under conditions with strong Allee effects, while asymmetric tumor-immune interaction kernels lead to poor immune response. Thus, the length scale of tumor-immune interactions drives tumor growth and infiltration.

Radiomics and artificial intelligence for precision medicine in lung cancer treatment

Lung cancer is the leading cause of cancer-related deaths worldwide. It exhibits, at the mesoscopic scale, phenotypic characteristics that are generally indiscernible to the human eye but can be captured non-invasively on medical imaging as radiomic features, which can form a high dimensional data space amenable to machine learning. Radiomic features can be harnessed and used in an artificial intelligence paradigm to risk stratify patients, and predict for histological and molecular findings, and clinical outcome measures, thereby facilitating precision medicine for improving patient care. Compared to tissue sampling-driven approaches, radiomics-based methods are superior for being non-invasive, reproducible, cheaper, and less susceptible to intra-tumoral heterogeneity. This review focuses on the application of radiomics, combined with artificial intelligence, for delivering precision medicine in lung cancer treatment, with discussion centered on pioneering and groundbreaking works, and future research directions in the area.

PAK1 copy number in breast cancer-Associations with proliferation and molecular subtypes

INTRODUCTION

P21-activated kinase 1 (PAK1) is known to be overexpressed in several human tumour types, including breast cancer (BC). It is located on chromosome 11 (11q13.5-q14.1) and plays a significant role in proliferation in BC. In this study we aimed to assess PAK1 gene copy number (CN) in primary breast tumours and their corresponding lymph node metastases, and associations between PAK1 CN and proliferation status, molecular subtype, and prognosis. In addition, we aimed to study associations between CNs of PAK1 and CCND1. Both genes are located on the long arm of chromosome 11 (11q13).

METHODS

Fluorescence in situ hybridization for PAK1 and Chromosome enumeration probe (CEP)11 were used on tissue microarray sections from a series of 512 BC cases. Copy numbers were estimated by counting the number of fluorescent signals for PAK1 and CEP11 in 20 tumour cell nuclei. Pearson's x2 test was performed to assess associations between PAK1 CN and tumour features, and between PAK1 and CCND1 CNs. Cumulative risk of death from BC and hazard ratios were estimated in analysis of prognosis.

RESULTS

We found mean PAK1 CN ≥4<6 in 26 (5.1%) tumours, and CN ≥ 6 in 22 (4.3%) tumours. The proportion of cases with copy number increase (mean CN ≥4) was highest among HER2 type and Luminal B (HER2-) tumours. We found an association between PAK1 CN increase, and high proliferation, and high histological grade, but not prognosis. Of cases with PAK1 CN ≥ 6, 30% also had CCND1 CN ≥ 6.

CONCLUSIONS

PAK1 copy number increase is associated with high proliferation and high histological grade, but not with prognosis. PAK1 CN increase was most frequent in the HER2 type and Luminal B (HER2-) subtype. PAK1 CN increase is associated with CN increase of CCND1.

Deciphering tumor ecosystems at super resolution from spatial transcriptomics with TESLA

Cell populations in the tumor microenvironment (TME), including their abundance, composition, and spatial location, are critical determinants of patient response to therapy. Recent advances in spatial transcriptomics (ST) have enabled the comprehensive characterization of gene expression in the TME. However, popular ST platforms, such as Visium, only measure expression in low-resolution spots and have large tissue areas that are not covered by any spots, which limits their usefulness in studying the detailed structure of TME. Here, we present TESLA, a machine learning framework for tissue annotation with pixel-level resolution in ST. TESLA integrates histological information with gene expression to annotate heterogeneous immune and tumor cells directly on the histology image. TESLA further detects unique TME features such as tertiary lymphoid structures, which represents a promising avenue for understanding the spatial architecture of the TME. Although we mainly illustrated the applications in cancer, TESLA can also be applied to other diseases.

Agent-based methods facilitate integrative science in cancer

In this opinion, we highlight agent-based modeling as a key tool for exploration of cell-cell and cell-environment interactions that drive cancer progression, therapeutic resistance, and metastasis. These biological phenomena are particularly suited to be captured at the cell-scale resolution possible only within agent-based or individual-based mathematical models. These modeling approaches complement experimental work (in vitro and in vivo systems) through parameterization and data extrapolation but also feed forward to drive new experiments that test model-generated predictions.

Metabolic activity grows in human cancers pushed by phenotypic variability

Different evolutionary processes push cancers to increasingly aggressive behaviors, energetically sustained by metabolic reprogramming. The collective signature emerging from this transition is macroscopically displayed by positron emission tomography (PET). In fact, the most readily PET measure, the maximum standardized uptake value (SUV_max), has been found to have prognostic value in different cancers. However, few works have linked the properties of this metabolic hotspot to cancer evolutionary dynamics. Here, by analyzing diagnostic PET images from 512 patients with cancer, we found that SUV_max scales superlinearly with the mean metabolic activity (SUV_mean), reflecting a dynamic preferential accumulation of activity on the hotspot. Additionally, SUV_max increased with metabolic tumor volume (MTV) following a power law. The behavior from the patients data was accurately captured by a mechanistic evolutionary dynamics model of tumor growth accounting for phenotypic transitions. This suggests that non-genetic changes may suffice to fuel the observed sustained increases in tumor metabolic activity.

Stochastic Fluctuations Drive Non-genetic Evolution of Proliferation in Clonal Cancer Cell Populations

Evolutionary dynamics allows us to understand many changes happening in a broad variety of biological systems, ranging from individuals to complete ecosystems. It is also behind a number of remarkable organizational changes that happen during the natural history of cancers. These reflect tumour heterogeneity, which is present at all cellular levels, including the genome, proteome and phenome, shaping its development and interrelation with its environment. An intriguing observation in different cohorts of oncological patients is that tumours exhibit an increased proliferation as the disease progresses, while the timescales involved are apparently too short for the fixation of sufficient driver mutations to promote explosive growth. Here, we discuss how phenotypic plasticity, emerging from a single genotype, may play a key role and provide a ground for a continuous acceleration of the proliferation rate of clonal populations with time. We address this question by combining the analysis of real-time growth of non-small-cell lung carcinoma cells (N-H460) together with stochastic and deterministic mathematical models that capture proliferation trait heterogeneity in clonal populations to elucidate the contribution of phenotypic transitions on tumour growth dynamics.

A Head-to-Head Comparison of 18F-Fluorocholine PET/CT and Conventional MRI as Predictors of Outcome in IDH Wild-Type High-Grade Gliomas

(1) Aim: To study the associations between imaging parameters derived from contrast-enhanced MRI (CE-MRI) and 18F-fluorocholine PET/CT and their performance as prognostic predictors in isocitrate dehydrogenase wild-type (IDH-wt) high-grade gliomas. (2) Methods: A prospective, multicenter study (FuMeGA: Functional and Metabolic Glioma Analysis) including patients with baseline CE-MRI and 18F-fluorocholine PET/CT and IDH wild-type high-grade gliomas. Clinical variables such as performance status, extent of surgery and adjuvant treatments (Stupp protocol vs others) were obtained and used to discriminate overall survival (OS) and progression-free survival (PFS) as end points. Multilesionality was assessed on the visual analysis of PET/CT and CE-MRI images. After tumor segmentation, standardized uptake value (SUV)-based variables for PET/CT and volume-based and geometrical variables for PET/CT and CE-MRI were calculated. The relationships among imaging techniques variables and their association with prognosis were evaluated using Pearson’s chi-square test and the t-test. Receiver operator characteristic, Kaplan−Meier and Cox regression were used for the survival analysis. (3) Results: 54 patients were assessed. The median PFS and OS were 5 and 11 months, respectively. Significant strong relationships between volume-dependent variables obtained from PET/CT and CE-MRI were found (r > 0.750, p < 0.05). For OS, significant associations were found with SUVmax, SUVpeak, SUVmean and sphericity (HR: 1.17, p = 0.035; HR: 1.24, p = 0.042; HR: 1.62, p = 0.040 and HR: 0.8, p = 0.022, respectively). Among clinical variables, only Stupp protocol and age showed significant associations with OS and PFS. No CE-MRI derived variables showed significant association with prognosis. In multivariate analysis, age (HR: 1.04, p = 0.002), Stupp protocol (HR: 2.81, p = 0.001), multilesionality (HR: 2.20, p = 0.013) and sphericity (HR: 0.79, p = 0.027) derived from PET/CT showed independent associations with OS. For PFS, only age (HR: 1.03, p = 0.021) and treatment protocol (HR: 2.20, p = 0.008) were significant predictors. (4) Conclusions: 18F-fluorocholine PET/CT metabolic and radiomic variables were robust prognostic predictors in patients with IDH-wt high-grade gliomas, outperforming CE-MRI derived variables.

Cell type-specific inference of differential expression in spatial transcriptomics

A central problem in spatial transcriptomics is detecting differentially expressed (DE) genes within cell types across tissue context. Challenges to learning DE include changing cell type composition across space and measurement pixels detecting transcripts from multiple cell types. Here, we introduce a statistical method, cell type-specific inference of differential expression (C-SIDE), that identifies cell type-specific DE in spatial transcriptomics, accounting for localization of other cell types. We model gene expression as an additive mixture across cell types of log-linear cell type-specific expression functions. C-SIDE's framework applies to many contexts: DE due to pathology, anatomical regions, cell-to-cell interactions and cellular microenvironment. Furthermore, C-SIDE enables statistical inference across multiple/replicates. Simulations and validation experiments on Slide-seq, MERFISH and Visium datasets demonstrate that C-SIDE accurately identifies DE with valid uncertainty quantification. Last, we apply C-SIDE to identify plaque-dependent immune activity in Alzheimer's disease and cellular interactions between tumor and immune cells. We distribute C-SIDE within the R package https://github.com/dmcable/spacexr .

SUVmax to tumor perimeter distance: a robust radiomics prognostic biomarker in resectable non-small cell lung cancer patients

OBJECTIVE

The purpose of this study was to evaluate the prognostic value of novel geometric variables obtained from pre-treatment [¹⁸F]FDG PET/CT with respect to classical ones in patients with non-small cell lung cancer (NSCLC).

METHODS

Retrospective study including stage I-III NSCLC patients with baseline [¹⁸F]FDG PET/CT. Clinical, histopathologic, and metabolic parameters were obtained. After tumor segmentation, SUV and volume-based variables, global texture, sphericity, and two novel parameters, normalized SUVpeak to centroid distance (nSCD) and normalized SUVmax to perimeter distance (nSPD), were obtained. Early recurrence (ER) and short-term mortality (STM) were used as end points. Univariate logistic regression and multivariate logistic regression with respect to ER and STM were performed.

RESULTS

A cohort of 173 patients was selected. ER was detected in 49/104 of patients with recurrent disease. Additionally, 100 patients died and 53 had STM. Age, pathologic lymphovascular invasion, lymph nodal infiltration, TNM stage, nSCD, and nSPD were associated with ER, although only age (aOR = 1.06, p = 0.002), pathologic lymphovascular invasion (aOR = 3.40, p = 0.022), and nSPD (aOR = 0.02, p = 0.018) were significant independent predictors of ER in multivariate analysis. Age, lymph nodal infiltration, TNM stage, nSCD, and nSPD were predictors of STM. Age (aOR = 1.05, p = 0.006), lymph nodal infiltration (aOR = 2.72, p = 0.005), and nSPD (aOR = 0.03, p = 0.022) were significantly associated with STM in multivariate analysis. Coefficient of variation (COV) and SUVmean/SUVmax ratio did not show significant predictive value with respect to ER or STM.

CONCLUSION

The geometric variables, nSCD and nSPD, are robust biomarkers of the poorest outcome prediction of patients with NSCLC with respect to classical PET variables.

KEY POINTS

• In NSCLC patients, it is crucial to find prognostic parameters since TNM system alone cannot explain the variation in lung cancer survival. • Age, lymphovascular invasion, lymph nodal infiltration, and metabolic geometrical parameters were useful as prognostic parameters. • The displacement grade of the highest point of metabolic activity towards the periphery assessed by geometric variables obtained from [¹⁸F]FDG PET/CT was a robust biomarker of the poorest outcome prediction of patients with NSCLC.