Radiogenomics of diffuse intrinsic pontine gliomas (DIPGs): correlation of histological and biological characteristics with multimodal MRI features.

OBJECTIVES: The diffuse intrinsic pontine gliomas (DIPGs) are now defined by the type of histone H3 mutated at lysine 27. We aimed to correlate the multimodal MRI features of DIPGs, H3K27M mutant, with their histological and molecular characteristics. METHODS: Twenty-seven treatment-naïve children with histopathologically confirmed DIPG H3K27M mutant were prospectively included. MRI performed prior to biopsy included multi-b-value diffusion-weighted imaging, ASL, and dynamic susceptibility contrast (DSC) perfusion imaging. The ADC and cerebral blood flow (CBF) and blood volume (CBV) were measured at the biopsy site. We assessed quantitative histological data, including microvascular density, nuclear density, and H3K27M-positive nuclear density. Gene expression profiling was also assessed in the samples. We compared imaging and histopathological data according to histone subgroup. We correlated MRI quantitative data with histological data and gene expression. RESULTS: H3.1K27M mutated tumors showed higher ADC values (median 3151 µm2/s vs 1741 µm2/s, p = 0.003), and lower perfusion values (DSC-rCBF median 0.71 vs 1.43, p = 0.002, and DSC-rCBV median 1.00 vs 1.71, p = 0.02) than H3.3K27M ones. They had similar microvascular and nuclear density, but lower H3K27M-positive nuclear density (p = 0.007). The DSC-rCBV was positively correlated to the H3K27M-positive nuclear density (rho = 0.74, p = 0.02). ADC values were not correlated with nuclear density nor perfusion values with microvascular density. The expression of gated channel activity-related genes tended to be inversely correlated with ADC values and positively correlated with DSC perfusion. CONCLUSIONS: H3.1K27M mutated tumors have higher ADC and lower perfusion values than H3.3K27M ones, without direct correlation with microvascular or nuclear density. This may be due to tissular edema possibly related to gated channel activity-related gene expression. KEY POINTS: • H3.1K27M mutant DIPG had higher apparent diffusion coefficient (p = 0.003), lower α (p = 0.048), and lower relative cerebral blood volume (p = 0.02) than H3.3K27M mutant DIPG at their biopsy sites. • Biopsy samples obtained within the tumor's enhancing portion showed higher microvascular density (p = 0.03) than samples obtained outside the tumor's enhancing portion, but similar H3K27M-positive nuclear density (p = 0.84). • Relative cerebral blood volume measured at the biopsy site was significantly correlated with H3K27M-positive nuclear density (rho = 0.74, p = 0.02).

Experimental and modeling study of the formation of cell aggregates with differential substrate adhesion.

The study of cell aggregation in vitro has a tremendous importance these days. In cancer biology, aggregates and spheroids serve as model systems and are considered as pseudo-tumors that are more realistic than 2D cell cultures. Recently, in the context of brain tumors (gliomas), we developed a new poly(ethylene glycol) (PEG)-based hydrogel, with adhesive properties that can be controlled by the addition of poly(L-lysine) (PLL), and a stiffness close to the brain's. This substrate allows the motion of individual cells and the formation of cell aggregates (within one day), and we showed that on a non-adhesive substrate (PEG without PLL is inert for cells), the aggregates are bigger and less numerous than on an adhesive substrate (with PLL). In this article, we present new experimental results on the follow-up of the formation of aggregates on our hydrogels, from the early stages (individual cells) to the late stages (aggregate compaction), in order to compare, for two cell lines (F98 and U87), the aggregation process on the adhesive and non-adhesive substrates. We first show that a spaceless model of perikinetic aggregation can reproduce the experimental evolution of the number of aggregates, but not of the mean area of the aggregates. We thus develop a minimal off-lattice agent-based model, with a few simple rules reproducing the main processes that are at stack during aggregation. Our spatial model can reproduce very well the experimental temporal evolution of both the number of aggregates and their mean area, on adhesive and non-adhesive soft gels and for the two different cell lines. From the fit of the experimental data, we were able to infer the quantitative values of the speed of motion of each cell line, its rate of proliferation in aggregates and its ability to organize in 3D. We also found qualitative differences between the two cell lines regarding the ability of aggregates to compact. These parameters could be inferred for any cell line, and correlated with clinical properties such as aggressiveness and invasiveness.

WHO grade has no prognostic value in the pediatric high-grade glioma included in the HERBY trial.

BACKGROUND: The World Health Organization (WHO) adult glioma grading system is questionable in pediatric high-grade gliomas (pHGGs), which are biologically distinct from adult HGGs. We took advantage of the neuropathological review data obtained during one of the largest prospective randomized pHGG trials, namely HERBY (NCT01390948), to address this issue in children with newly diagnosed non-brainstem HGG. METHODS: HGG diagnosis was confirmed by pre-randomization, real-time central pathology review using WHO 2007 criteria, followed by a consensus review blinded to clinical factors and outcomes. We evaluated association between WHO 2007 grade and other clinical/radiological/biological characteristics and the prognostic value of WHO 2007 grade, midline location, and selected biomarkers (Ki-67 index/Olig2/CD34/EGFR/p53/H3F3A K27M mutation) on overall survival. RESULTS: Real-time central neuropathological review was feasible in a multicenter study, with a mean time of 2.4 days, and led to the rejection of HGG diagnosis in 20 of 163 cases (12.3%). The different grading criteria and resulting WHO grade were not significantly associated with overall survival in the entire population (n = 118) or in midline and non-midline subgroups. H3F3A K27M mutation was significantly associated with poor outcome. No significant prognostic value was observed for grade, even after regrading H3F3A K27M-mutated midline glioma as grade IV (WHO 2016). Midline location and a high Ki-67 index (≥20%) were associated with poor outcome (P = 0.004 and P = 0.04, respectively). A 10% increase in Ki-67 index was associated with a hazard ratio of 1.53 (95% CI: 1.27-1.83; P < 0.0001). CONCLUSION: Our findings suggest that WHO grade III versus IV has no prognostic value in pediatric HGG.

Modeling the dynamics of oligodendrocyte precursor cells and the genesis of gliomas.

Oligodendrocyte precursor cells (OPCs) have remarkable properties: they represent the most abundant cycling cell population in the adult normal brain and they manage to achieve a uniform and constant density throughout the adult brain. This equilibrium is obtained by the interplay of four processes: division, differentiation or death, migration and active self-repulsion. They are also strongly suspected to be at the origin of gliomas, when their equilibrium is disrupted. In this article, we present a model of the dynamics of OPCs, first in a normal tissue. This model is based on a cellular automaton and its rules are mimicking the ones that regulate the dynamics of real OPCs. The model is able to reproduce the homeostasis of the cell population, with the maintenance of a constant and uniform cell density and the healing of a lesion. We show that there exists a fair quantitative agreement between the simulated and experimental parameters, such as the cell velocity, the time taken to close a lesion, and the duration of the cell cycle. We present three possible scenarios of disruption of the equilibrium: the appearance of an over-proliferating cell, of a deadless/non-differentiating cell, or of a cell that lost any contact-inhibition. We show that the appearance of an over-proliferating cell is sufficient to trigger the growth of a tumor that has low-grade glioma features: an invasive behaviour, a linear radial growth of the tumor with a corresponding growth velocity of less than 2 mm per year, as well a cell density at the center which exceeds the one in normal tissue by a factor of less than two. The loss of contact inhibition leads to a more high-grade-like glioma. The results of our model contribute to the body of evidence that identify OPCs as possible cells of origin of gliomas.

Interactions between glioma and pregnancy: insight from a 52-case multicenter series.

OBJECTIVE The goal of this study was to provide insight into the influence of gliomas on gestational outcomes, the impact of pregnancy on gliomas, and the identification of patients at risk. METHODS In this multiinstitutional retrospective study, the authors identified 52 pregnancies in 50 women diagnosed with a glioma. RESULTS For gliomas known prior to pregnancy (n = 24), we found the following: 1) An increase in the quantified imaging growth rates occurred during pregnancy in 87% of cases. 2) Clinical deterioration occurred in 38% of cases, with seizures alone resolving after delivery in 57.2% of cases. 3) Oncological treatments were immediately performed after delivery in 25% of cases. For gliomas diagnosed during pregnancy (n = 28), we demonstrated the following: 1) The tumor was discovered during the second and third trimesters in 29% and 54% of cases, respectively, with seizures being the presenting symptom in 68% of cases. 2) The quantified imaging growth rates did not significantly decrease after delivery and before oncological treatment. 3) Clinical deterioration resolved after delivery in 21.4% of cases. 4) Oncological treatments were immediately performed after delivery in 70% of cases. Gliomas with a high grade of malignancy, negative immunoexpression of alpha-internexin, or positive immunoexpression for p53 were more likely to be associated with tumor progression during pregnancy. Deliveries were all uneventful (cesarean section in 54.5% of cases and vaginal delivery in 45.5%), and the infants were developmentally normal. CONCLUSIONS When a woman harboring a glioma envisions a pregnancy, or when a glioma is discovered in a pregnant patient, the authors suggest informing her and her partner that pregnancy may impact the evolution of the glioma clinically and radiologically. They strongly advise a multidisciplinary approach to management. ■ CLASSIFICATION OF EVIDENCE Type of question: association; study design: case series; evidence: Class IV.

Surgical Decision Making From Image-Based Biophysical Modeling of Glioblastoma: Not Ready for Primetime.

BACKGROUND: Biophysical modeling of glioma is gaining more interest for clinical practice. The most popular model describes aggressivity of tumor cells by two parameters: net proliferation rate (ρ) and propensity to migrate (D). The ratio ρ/D, which can be estimated from a single preoperative magnetic resonance imaging (MRI), characterizes tumor invasiveness profile (high ρ/D: nodular; low ρ/D: diffuse). A recent study reported, from a large series of glioblastoma multiforme (GBM) patients, that gross total resection (GTR) would improve survival only in patients with nodular tumors. OBJECTIVE: To replicate these results, that is to verify that benefit of GTR would be only observed for nodular tumors. METHODS: Between 2005 and 2012, we considered 234 GBM patients with pre- and postoperative MRI. Stereotactic biopsy (BST) was performed in 109 patients. Extent of resection was assessed on postoperative MRI and classified as GTR or partial resection (PR). Invasiveness ρ/D was estimated from the preoperative tumor volumes on T1-Gadolinium-enhanced and fluid-attenuated inversion recovery sequences. RESULTS: We demonstrate that patients with diffuse GBM (low ρ/D), as well as more nodular (mid and high ρ/D) GBM, presented significant survival benefit from GTR over PR/BST ( P < .001). CONCLUSION: Whatever the degree of tumor invasiveness, as estimated from MRI-driven biophysical modeling, GTR improves survival of GBM patients, compared to PR or BST. This conflicting result should motivate further studies.

Arterial Spin Labeling to Predict Brain Tumor Grading in Children: Correlations between Histopathologic Vascular Density and Perfusion MR Imaging.

Purpose To compare arterial spin labeling (ASL) data between low- and high-grade brain tumors in children to establish a cutoff to distinguish low- from high-grade neoplasms and to assess potential correlations between cerebral blood flow (CBF) and quantitative histologic microvascular data. Materials and Methods Approval was obtained from the regional review board. ASL data obtained in 129 children between 2011 and 2015 were retrospectively analyzed. CBF and relative CBF in the most perfused area of each neoplasm and contrast enhancement were quantified with a semiquantitative ratio. The correlation between CBF and microvascular density was analyzed in specimens stained with anti-CD34. Results were controlled in two validation cohorts with 1.5- and 3.0-T magnetic resonance (MR) imaging. Results Mean CBF was significantly higher for high-grade than for low-grade hemispheric (116 mL/min/100 g [interquartile range {IQR}, 73-131 mL/min/100 g] vs 29 mL/min/100 g [IQR, 23-35 29 mL/min/100 g], P < .001), thalamic (87 mL/min/100 g [IQR, 73-100 mL/min/100 g] vs 36 mL/min/100 g [IQR, 30-40 mL/min/100 g], P = .016), and posterior fossa (59 mL/min/100 g [IQR, 45-91 mL/min/100 g] vs 33 mL/min/100 g [IQR, 25-40 mL/min/100 g], P < .001) tumors. With a cutoff of 50 mL/min/100 g, sensitivity and specificity were 90% (95% confidence interval [CI]: 68, 100) and 93% (95% CI: 66, 100), respectively, for hemispheric tumors; 100% (95% CI: 48, 100) and 80% (95% CI: 28, 100), respectively, for thalamic tumors; and 65% (95% CI: 51, 78) and 94% (95% CI: 80, 99), respectively, for posterior fossa tumors. In posterior fossa tumors, additional use of the CBF-to-contrast enhancement ratio yielded sensitivity and specificity of 96% (95% CI: 87, 100) and 97% (95% CI: 84, 100), respectively. Use of a simple algorithm based on these values yielded an accuracy of 93% (95% CI: 87, 97). Validation sets yielded similar results, with grading accuracy of 88% (95% CI: 62, 98) with 1.5-T MR imaging and 77% (95% CI: 46, 95) with 3.0-T MR imaging. CBF was strongly correlated with microvascular density (R = 0.66, P < .001). Conclusion High-grade pediatric brain tumors display higher CBF than do low-grade tumors, and they may be accurately graded by using these values. CBF is correlated with tumor microvascular density. © RSNA, 2016 Online supplemental material is available for this article.

Towards better digital pathology workflows: programming libraries for high-speed sharpness assessment of Whole Slide Images.

BACKGROUND: Since microscopic slides can now be automatically digitized and integrated in the clinical workflow, quality assessment of Whole Slide Images (WSI) has become a crucial issue. We present a no-reference quality assessment method that has been thoroughly tested since 2010 and is under implementation in multiple sites, both public university-hospitals and private entities. It is part of the FlexMIm R&D project which aims to improve the global workflow of digital pathology. For these uses, we have developed two programming libraries, in Java and Python, which can be integrated in various types of WSI acquisition systems, viewers and image analysis tools. METHODS: Development and testing have been carried out on a MacBook Pro i7 and on a bi-Xeon 2.7GHz server. Libraries implementing the blur assessment method have been developed in Java, Python, PHP5 and MySQL5. For web applications, JavaScript, Ajax, JSON and Sockets were also used, as well as the Google Maps API. Aperio SVS files were converted into the Google Maps format using VIPS and Openslide libraries. RESULTS: We designed the Java library as a Service Provider Interface (SPI), extendable by third parties. Analysis is computed in real-time (3 billion pixels per minute). Tests were made on 5000 single images, 200 NDPI WSI, 100 Aperio SVS WSI converted to the Google Maps format. CONCLUSIONS: Applications based on our method and libraries can be used upstream, as calibration and quality control tool for the WSI acquisition systems, or as tools to reacquire tiles while the WSI is being scanned. They can also be used downstream to reacquire the complete slides that are below the quality threshold for surgical pathology analysis. WSI may also be displayed in a smarter way by sending and displaying the regions of highest quality before other regions. Such quality assessment scores could be integrated as WSI's metadata shared in clinical, research or teaching contexts, for a more efficient medical informatics workflow.

Analyzing huge pathology images with open source software.

BACKGROUND: Digital pathology images are increasingly used both for diagnosis and research, because slide scanners are nowadays broadly available and because the quantitative study of these images yields new insights in systems biology. However, such virtual slides build up a technical challenge since the images occupy often several gigabytes and cannot be fully opened in a computer's memory. Moreover, there is no standard format. Therefore, most common open source tools such as ImageJ fail at treating them, and the others require expensive hardware while still being prohibitively slow. RESULTS: We have developed several cross-platform open source software tools to overcome these limitations. The NDPITools provide a way to transform microscopy images initially in the loosely supported NDPI format into one or several standard TIFF files, and to create mosaics (division of huge images into small ones, with or without overlap) in various TIFF and JPEG formats. They can be driven through ImageJ plugins. The LargeTIFFTools achieve similar functionality for huge TIFF images which do not fit into RAM. We test the performance of these tools on several digital slides and compare them, when applicable, to standard software. A statistical study of the cells in a tissue sample from an oligodendroglioma was performed on an average laptop computer to demonstrate the efficiency of the tools. CONCLUSIONS: Our open source software enables dealing with huge images with standard software on average computers. They are cross-platform, independent of proprietary libraries and very modular, allowing them to be used in other open source projects. They have excellent performance in terms of execution speed and RAM requirements. They open promising perspectives both to the clinician who wants to study a single slide and to the research team or data centre who do image analysis of many slides on a computer cluster. VIRTUAL SLIDES: The virtual slide(s) for this article can be found here:http://www.diagnosticpathology.diagnomx.eu/vs/5955513929846272.

Quantitative characterization of the imaging limits of diffuse low-grade oligodendrogliomas.

BACKGROUND: Supratentorial diffuse low-grade gliomas in adults extend beyond maximal visible MRI-defined abnormalities, and a gap exists between the imaging signal changes and the actual tumor margins. Direct quantitative comparisons between imaging and histological analyses are lacking to date. However, they are of the utmost importance if one wishes to develop realistic models for diffuse glioma growth. METHODS: In this study, we quantitatively compared the cell concentration and the edema fraction from human histological biopsy samples (BSs) performed inside and outside imaging abnormalities during serial imaging-based stereotactic biopsy of diffuse low-grade gliomas. RESULTS: The cell concentration was significantly higher in BSs located inside (1189 ± 378 cell/mm(2)) than outside (740 ± 124 cell/mm(2)) MRI-defined abnormalities (P = .0003). The edema fraction was significantly higher in BSs located inside (mean, 45% ± 23%) than outside (mean, 5 %± 9%) MRI-defined abnormalities (P < .0001). At borders of the MRI-defined abnormalities, 20% of the tissue surface area was occupied by edema and only 3% by tumor cells. The cycling cell concentration was significantly higher in BSs located inside (10 ± 12 cell/mm(2)), compared with outside (0.5 ± 0.9 cell/mm(2)), MRI-defined abnormalities (P = .0001). CONCLUSIONS: We showed that the margins of T2-weighted signal changes are mainly correlated with the edema fraction. In 62.5% of patients, the cycling tumor cell fraction (defined as the ratio of the cycling tumor cell concentration to the total number of tumor cells) was higher at the limits of the MRI-defined abnormalities than closer to the center of the tumor. In the remaining patients, the cycling tumor cell fraction increased towards the center of the tumor.

Exclusion processes: short-range correlations induced by adhesion and contact interactions.

We analyze the out-of-equilibrium behavior of exclusion processes where agents interact with their nearest neighbors, and we study the short-range correlations which develop because of the exclusion and other contact interactions. The form of interactions we focus on, including adhesion and contact-preserving interactions, is especially relevant for migration processes of living cells. We show the local agent density and nearest-neighbor two-point correlations resulting from simulations on two-dimensional lattices in the transient regime where agents invade an initially empty space from a source and in the stationary regime between a source and a sink. We compare the results of simulations with the corresponding quantities derived from the master equation of the exclusion processes, and in both cases, we show that, during the invasion of space by agents, a wave of correlations travels with velocity v(t)~t(-1/2). The relative placement of this wave to the agent density front and the time dependence of its height may be used to discriminate between different forms of contact interactions or to quantitatively estimate the intensity of interactions. We discuss, in the stationary density profile between a full and an empty reservoir of agents, the presence of a discontinuity close to the empty reservoir. Then we develop a method for deriving approximate hydrodynamic limits of the processes. From the resulting systems of partial differential equations, we recover the self-similar behavior of the agent density and correlations during space invasion.

Pregnancy increases the growth rates of World Health Organization grade II gliomas.

Twelve pregnancies in 11 adult women harboring World Health Organization (WHO) grade II gliomas (GIIGs) prior to pregnancy were reviewed to address whether pregnancy affects tumor growth using a quantitative approach of the radiological velocity of diametric expansion (VDE) on successive magnetic resonance images. VDE was significantly increased during pregnancy as compared to prepregnancy (p < 0.001) and to postdelivery (p = 0.012) periods. Pregnancy increases the radiological growth rates of GIIGs. An increase in seizure frequency was observed concomitantly in 40% of cases and further oncological treatment was started after delivery in 25% of cases.

Modeling tumor cell migration: From microscopic to macroscopic models.

It has been shown experimentally that contact interactions may influence the migration of cancer cells. Previous works have modelized this thanks to stochastic, discrete models (cellular automata) at the cell level. However, for the study of the growth of real-size tumors with several million cells, it is best to use a macroscopic model having the form of a partial differential equation (PDE) for the density of cells. The difficulty is to predict the effect, at the macroscopic scale, of contact interactions that take place at the microscopic scale. To address this, we use a multiscale approach: starting from a very simple, yet experimentally validated, microscopic model of migration with contact interactions, we derive a macroscopic model. We show that a diffusion equation arises, as is often postulated in the field of glioma modeling, but it is nonlinear because of the interactions. We give the explicit dependence of diffusivity on the cell density and on a parameter governing cell-cell interactions. We discuss in detail the conditions of validity of the approximations used in the derivation, and we compare analytic results from our PDE to numerical simulations and to some in vitro experiments. We notice that the family of microscopic models we started from includes as special cases some kinetically constrained models that were introduced for the study of the physics of glasses, supercooled liquids, and jamming systems.

Field-theoretic approach to metastability in the contact process.

A "quantum" field-theoretic formulation of the dynamics of the contact process on a regular graph of degree z is introduced. A perturbative calculation in powers of 1/z of the effective potential for the density of particles phi(t) and an instantonic field psi(t) emerging from the formalism is performed. Corrections to the mean-field distribution of densities of particles in the out-of-equilibrium stationary state are derived in powers of 1/z. Results for typical (e.g., average density) and rare fluctuation (e.g. lifetime of the metastable state) properties are in very good agreement with numerical simulations carried out on D-dimensional hypercubic (z=2D) and Cayley lattices.