Textural measures have been widely explored as imaging biomarkers in cancer. However, their robustness under dynamic range and spatial resolution changes in brain 3D magnetic resonance images (MRI) ...has not been assessed. The aim of this work was to study potential variations of textural measures due to changes in MRI protocols.
Twenty patients harboring glioblastoma with pretreatment 3D T1-weighted MRIs were included in the study. Four different spatial resolution combinations and three dynamic ranges were studied for each patient. Sixteen three-dimensional textural heterogeneity measures were computed for each patient and configuration including co-occurrence matrices (CM) features and run-length matrices (RLM) features. The coefficient of variation was used to assess the robustness of the measures in two series of experiments corresponding to (i) changing the dynamic range and (ii) changing the matrix size.
No textural measures were robust under dynamic range changes. Entropy was the only textural feature robust under spatial resolution changes (coefficient of variation under 10% in all cases).
Textural measures of three-dimensional brain tumor images are not robust neither under dynamic range nor under matrix size changes. Standards should be harmonized to use textural features as imaging biomarkers in radiomic-based studies. The implications of this work go beyond the specific tumor type studied here and pose the need for standardization in textural feature calculation of oncological images.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Gliomas are the most frequent type of primary brain tumours. Low grade gliomas (LGGs, WHO grade II gliomas) may grow very slowly for the long periods of time, however they inevitably cause death due ...to the phenomenon known as the malignant transformation. This refers to the transition of LGGs to more aggressive forms of high grade gliomas (HGGs, WHO grade III and IV gliomas). In this paper we propose a mathematical model describing the spatio-temporal transition of LGGs into HGGs. Our modelling approach is based on two cellular populations with transitions between them being driven by the tumour microenvironment transformation occurring when the tumour cell density grows beyond a critical level. We show that the proposed model describes real patient data well. We discuss the relationship between patient prognosis and model parameters. We approximate tumour radius and velocity before malignant transformation as well as estimate the onset of this process.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Increasingly complex in silico modeling approaches offer a way to simultaneously access cancerous processes at different spatio-temporal scales. High-level models, such as those based on partial ...differential equations, are computationally affordable and allow large tumor sizes and long temporal windows to be studied, but miss the discrete nature of many key underlying cellular processes. Individual-based approaches provide a much more detailed description of tumors, but have difficulties when trying to handle full-sized real cancers. Thus, there exists a trade-off between the integration of macroscopic and microscopic information, now widely available, and the ability to attain clinical tumor sizes. In this paper we put forward a stochastic mesoscopic simulation framework that incorporates key cellular processes during tumor progression while keeping computational costs to a minimum. Our framework captures a physical scale that allows both the incorporation of microscopic information, tracking the spatio-temporal emergence of tumor heterogeneity and the underlying evolutionary dynamics, and the reconstruction of clinically sized tumors from high-resolution medical imaging data, with the additional benefit of low computational cost. We illustrate the functionality of our modeling approach for the case of glioblastoma, a paradigm of tumor heterogeneity that remains extremely challenging in the clinical setting.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Brain metastasis (BM) is one of the main complications of many cancers, and the most frequent malignancy of the central nervous system. Imaging studies of BMs are routinely used for diagnosis of ...disease, treatment planning and follow-up. Artificial Intelligence (AI) has great potential to provide automated tools to assist in the management of disease. However, AI methods require large datasets for training and validation, and to date there have been just one publicly available imaging dataset of 156 BMs. This paper publishes 637 high-resolution imaging studies of 75 patients harboring 260 BM lesions, and their respective clinical data. It also includes semi-automatic segmentations of 593 BMs, including pre- and post-treatment T1-weighted cases, and a set of morphological and radiomic features for the cases segmented. This data-sharing initiative is expected to enable research into and performance evaluation of automatic BM detection, lesion segmentation, disease status evaluation and treatment planning methods for BMs, as well as the development and validation of predictive and prognostic tools with clinical applicability.
Metastasis is the process through which cancer cells break away from a primary tumor, travel through the blood or lymph system, and form new tumors in distant tissues. One of the preferred sites for ...metastatic dissemination is the brain, affecting more than 20% of all cancer patients. This figure is increasing steadily due to improvements in treatments of primary tumors. Stereotactic radiosurgery (SRS) is one of the main treatment options for patients with a small or moderate number of brain metastases (BMs). A frequent adverse event of SRS is radiation necrosis (RN), an inflammatory condition caused by late normal tissue cell death. A major diagnostic problem is that RNs are difficult to distinguish from BM recurrences, due to their similarities on standard magnetic resonance images (MRIs). However, this distinction is key to choosing the best therapeutic approach since RNs resolve often without further interventions, while relapsing BMs may require open brain surgery. Recent research has shown that RNs have a faster growth dynamics than recurrent BMs, providing a way to differentiate the two entities, but no mechanistic explanation has been provided for those observations. In this study, computational frameworks were developed based on mathematical models of increasing complexity, providing mechanistic explanations for the differential growth dynamics of BMs relapse versus RN events and explaining the observed clinical phenomenology. Simulated tumor relapses were found to have growth exponents substantially smaller than the group in which there was inflammation due to damage induced by SRS to normal brain tissue adjacent to the BMs, thus leading to RN. ROC curves with the synthetic data had an optimal threshold that maximized the sensitivity and specificity values for a growth exponent β* = 1.05, very close to that observed in patient datasets.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Most physical and other natural systems are complex entities composed of a large number of interacting individual elements. It is a surprising fact that they often obey the so-called scaling laws ...relating an observable quantity with a measure of the size of the system. Here we describe the discovery of universal superlinear metabolic scaling laws in human cancers. This dependence underpins increasing tumour aggressiveness, due to evolutionary dynamics, which leads to an explosive growth as the disease progresses. We validated this dynamic using longitudinal volumetric data of different histologies from large cohorts of cancer patients. To explain our observations we put forward increasingly-complex biologically-inspired mathematical models that captured the key processes governing tumor growth. Our models predicted that the emergence of superlinear allometric scaling laws is an inherently three-dimensional phenomenon. Moreover, the scaling laws thereby identified allowed us to define a set of metabolic metrics with prognostic value, thus providing added clinical utility to the base findings.
Full text
Available for:
FZAB, GEOZS, IJS, IMTLJ, KILJ, KISLJ, MFDPS, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBJE, SBMB, SBNM, UKNU, UL, UM, UPUK, VKSCE, ZAGLJ
Low-grade gliomas are primary brain tumors that arise from glial cells and are usually treated with temozolomide (TMZ) as a chemotherapeutic option. They are often incurable, but patients have a ...prolonged survival. One of the shortcomings of the treatment is that patients eventually develop drug resistance. Recent findings show that persisters, cells that enter a dormancy state to resist treatment, play an important role in the development of resistance to TMZ. In this study we constructed a mathematical model of low-grade glioma response to TMZ incorporating a persister population. The model was able to describe the volumetric longitudinal dynamics, observed in routine FLAIR 3D sequences, of low-grade glioma patients acquiring TMZ resistance. We used the model to explore different TMZ administration protocols, first on virtual clones of real patients and afterwards on virtual patients preserving the relationships between parameters of real patients. In silico clinical trials showed that resistance development was deferred by protocols in which individual doses are administered after rest periods, rather than the 28-days cycle standard protocol. This led to median survival gains in virtual patients of more than 15 months when using resting periods between two and three weeks and agreed with recent experimental observations in animal models. Additionally, we tested adaptive variations of these new protocols, what showed a potential reduction in toxicity, but no survival gain. Our computational results highlight the need of further clinical trials that could obtain better results from treatment with TMZ in low grade gliomas.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Grade II gliomas are slowly growing primary brain tumors that affect mostly young patients. Cytotoxic therapies (radiotherapy and/or chemotherapy) are used initially only for patients having a bad ...prognosis. These therapies are planned following the "maximum dose in minimum time" principle, i. e. the same schedule used for high-grade brain tumors in spite of their very different behavior. These tumors transform after a variable time into high-grade gliomas, which significantly decreases the patient's life expectancy. In this paper we study mathematical models describing the growth of grade II gliomas in response to radiotherapy. We find that protracted metronomic fractionations, i.e. therapeutical schedules enlarging the time interval between low-dose radiotherapy fractions, may lead to a better tumor control without an increase in toxicity. Other non-standard fractionations such as protracted or hypoprotracted schemes may also be beneficial. The potential survival improvement depends on the tumor's proliferation rate and can be even of the order of years. A conservative metronomic scheme, still being a suboptimal treatment, delays the time to malignant progression by at least one year when compared to the standard scheme.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Purpose To evaluate the prognostic and predictive value of surface-derived imaging biomarkers obtained from contrast material-enhanced volumetric T1-weighted pretreatment magnetic resonance (MR) ...imaging sequences in patients with glioblastoma multiforme. Materials and Methods A discovery cohort from five local institutions (165 patients; mean age, 62 years ± 12 standard deviation; 43% women and 57% men) and an independent validation cohort (51 patients; mean age, 60 years ± 12; 39% women and 61% men) from The Cancer Imaging Archive with volumetric T1-weighted pretreatment contrast-enhanced MR imaging sequences were included in the study. Clinical variables such as age, treatment, and survival were collected. After tumor segmentation and image processing, tumor surface regularity, measuring how much the tumor surface deviates from a sphere of the same volume, was obtained. Kaplan-Meier, Cox proportional hazards, correlations, and concordance indexes were used to compare variables and patient subgroups. Results Surface regularity was a powerful predictor of survival in the discovery (P = .005, hazard ratio HR = 1.61) and validation groups (P = .05, HR = 1.84). Multivariate analysis selected age and surface regularity as significant variables in a combined prognostic model (P < .001, HR = 3.05). The model achieved concordance indexes of 0.76 and 0.74 for the discovery and validation cohorts, respectively. Tumor surface regularity was a predictor of survival for patients who underwent complete resection (P = .01, HR = 1.90). Tumors with irregular surfaces did not benefit from total over subtotal resections (P = .57, HR = 1.17), but those with regular surfaces did (P = .004, HR = 2.07). Conclusion The surface regularity obtained from high-resolution contrast-enhanced pretreatment volumetric T1-weighted MR images is a predictor of survival in patients with glioblastoma. It may help in classifying patients for surgery.
Glioblastoma (GBM) is one of the most lethal tumor types. Hypercellular regions, named pseudopalisades, are characteristic in these tumors and have been hypothesized to be waves of migrating ...glioblastoma cells. These "waves" of cells are thought to be induced by oxygen and nutrient depletion caused by tumor-induced blood vessel occlusion. Although the universal presence of these structures in GBM tumors suggests that they may play an instrumental role in GBM's spread and invasion, the recreation of these structures in vitro has remained challenging.
Here we present a new microfluidic model of GBM that mimics the dynamics of pseudopalisade formation. To do this, we embedded U-251 MG cells within a collagen hydrogel in a custom-designed microfluidic device. By controlling the medium flow through lateral microchannels, we can mimic and control blood-vessel obstruction events associated with this disease.
Through the use of this new system, we show that nutrient and oxygen starvation triggers a strong migratory process leading to pseudopalisade generation in vitro. These results validate the hypothesis of pseudopalisade formation and show an excellent agreement with a systems-biology model based on a hypoxia-driven phenomenon.
This paper shows the potential of microfluidic devices as advanced artificial systems capable of modeling in vivo nutrient and oxygen gradients during tumor evolution.