Deep neural network to locate and segment brain tumors outperformed the expert technicians who created the training data

Purpose: Deep learning (DL) algorithms have shown promising results for brain tumor segmentation in MRI. However, validation is required prior to routine clinical use. We report the first randomized and blinded comparison of DL and trained technician segmentations.

Approach: We compiled a multi-institutional database of 741 pretreatment MRI exams. Each contained a postcontrast T1-weighted exam, a T2-weighted fluid-attenuated inversion recovery exam, and at least one technician-derived tumor segmentation. The database included 729 unique patients (470 males and 259 females). Of these exams, 641 were used for training the DL system, and 100 were reserved for testing. We developed a platform to enable qualitative, blinded, controlled assessment of lesion segmentations made by technicians and the DL method. On this platform, 20 neuroradiologists performed 400 side-by-side comparisons of segmentations on 100 test cases. They scored each segmentation between 0 (poor) and 10 (perfect). Agreement between segmentations from technicians and the DL method was also evaluated quantitatively using the Dice coefficient, which produces values between 0 (no overlap) and 1 (perfect overlap).

Results: The neuroradiologists gave technician and DL segmentations mean scores of 6.97 and 7.31, respectively (p<0.00007). The DL method achieved a mean Dice coefficient of 0.87 on the test cases.

Conclusions: This was the first objective comparison of automated and human segmentation using a blinded controlled assessment study. Our DL system learned to outperform its “human teachers” and produced output that was better, on average, than its training data.

Consensus recommendations for a dynamic susceptibility contrast MRI protocol for use in high-grade gliomas

Despite the widespread clinical use of dynamic susceptibility contrast (DSC) MRI, DSC-MRI methodology has not been standardized, hindering its utilization for response assessment in multicenter trials. Recently, the DSC-MRI Standardization Subcommittee of the Jumpstarting Brain Tumor Drug Development Coalition issued an updated consensus DSC-MRI protocol compatible with the standardized brain tumor imaging protocol (BTIP) for high-grade gliomas that is increasingly used in the clinical setting and is the default MRI protocol for the National Clinical Trials Network. After reviewing the basis for controversy over DSC-MRI protocols, this paper provides evidence-based best practices for clinical DSC-MRI as determined by the Committee, including pulse sequence (gradient echo vs spin echo), BTIP-compliant contrast agent dosing (preload and bolus), flip angle (FA), echo time (TE), and post-processing leakage correction. In summary, full-dose preload, full-dose bolus dosing using intermediate (60°) FA and field strength-dependent TE (40-50 ms at 1.5 T, 20-35 ms at 3 T) provides overall best accuracy and precision for cerebral blood volume estimates. When single-dose contrast agent usage is desired, no-preload, full-dose bolus dosing using low FA (30°) and field strength-dependent TE provides excellent performance, with reduced contrast agent usage and elimination of potential systematic errors introduced by variations in preload dose and incubation time.

Identifying the spatial and temporal dynamics of molecularly-distinct glioblastoma sub-populations

Glioblastomas (GBMs) are the most aggressive primary brain tumours and have no known cure. Each individual tumour comprises multiple sub-populations of genetically-distinct cells that may respond differently to targeted therapies and may contribute to disappointing clinical trial results. Image-localized biopsy techniques allow multiple biopsies to be taken during surgery and provide information that identifies regions where particular sub-populations occur within an individual GBM, thus providing insight into their regional genetic variability. These sub-populations may also interact with one another in a competitive or cooperative manner; it is important to ascertain the nature of these interactions, as they may have implications for responses to targeted therapies. We combine genetic information from biopsies with a mechanistic model of interacting GBM sub-populations to characterise the nature of interactions between two commonly occurring GBM sub-populations, those with EGFR and PDGFRA genes amplified. We study population levels found across image-localized biopsy data from a cohort of 25 patients and compare this to model outputs under competitive, cooperative and neutral interaction assumptions. We explore other factors affecting the observed simulated sub-populations, such as selection advantages and phylogenetic ordering of mutations, which may also contribute to the levels of EGFR and PDGFRA amplified populations observed in biopsy data.

Evaluating the Use of rCBV as a Tumor Grade and Treatment Response Classifier Across NCI Quantitative Imaging Network Sites: Part II of the DSC-MRI Digital Reference Object (DRO) Challenge

We have previously characterized the reproducibility of brain tumor relative cerebral blood volume (rCBV) using a dynamic susceptibility contrast magnetic resonance imaging digital reference object across 12 sites using a range of imaging protocols and software platforms. As expected, reproducibility was highest when imaging protocols and software were consistent, but decreased when they were variable. Our goal in this study was to determine the impact of rCBV reproducibility for tumor grade and treatment response classification. We found that varying imaging protocols and software platforms produced a range of optimal thresholds for both tumor grading and treatment response, but the performance of these thresholds was similar. These findings further underscore the importance of standardizing acquisition and analysis protocols across sites and software benchmarking.

Sex-specific impact of patterns of imageable tumor growth on survival of primary glioblastoma patients

Background

Sex is recognized as a significant determinant of outcome among glioblastoma patients, but the relative prognostic importance of glioblastoma features has not been thoroughly explored for sex differences.

Methods

Combining multi-modal MR images, biomathematical models, and patient clinical information, this investigation assesses which pretreatment variables have a sex-specific impact on the survival of glioblastoma patients (299 males and 195 females).

Results

Among males, tumor (T1Gd) radius was a predictor of overall survival (HR = 1.027, p = 0.044). Among females, higher tumor cell net invasion rate was a significant detriment to overall survival (HR = 1.011, p < 0.001). Female extreme survivors had significantly smaller tumors (T1Gd) (p = 0.010 t-test), but tumor size was not correlated with female overall survival (p = 0.955 CPH). Both male and female extreme survivors had significantly lower tumor cell net proliferation rates than other patients (M p = 0.004, F p = 0.001, t-test).

Conclusion

Despite similar distributions of the MR imaging parameters between males and females, there was a sex-specific difference in how these parameters related to outcomes.

Image-based metric of invasiveness predicts response to adjuvant temozolomide for primary glioblastoma

Background

Temozolomide (TMZ) has been the standard-of-care chemotherapy for glioblastoma (GBM) patients for more than a decade. Despite this long time in use, significant questions remain regarding how best to optimize TMZ therapy for individual patients. Understanding the relationship between TMZ response and factors such as number of adjuvant TMZ cycles, patient age, patient sex, and image–based tumor features, might help predict which GBM patients would benefit most from TMZ, particularly for those whose tumors lack O⁶–methylguanine–DNA methyltransferase (MGMT) promoter methylation.

Methods and findings

Using a cohort of 90 newly–diagnosed GBM patients treated according to the standard of care, we examined the relationships between several patient and tumor characteristics and volumetric and survival outcomes during adjuvant chemotherapy. Volumetric changes in MR imaging abnormalities during adjuvant therapy were used to assess TMZ response. T1Gd volumetric response is associated with younger patient age, increased number of TMZ cycles, longer time to nadir volume, and decreased tumor invasiveness. Moreover, increased adjuvant TMZ cycles corresponded with improved volumetric response only among more nodular tumors, and this volumetric response was associated with improved survival outcomes. Finally, in a subcohort of patients with known MGMT methylation status, methylated tumors were more diffusely invasive than unmethylated tumors, suggesting the improved response in nodular tumors is not driven by a preponderance of MGMT methylated tumors.

Conclusions

Our finding that less diffusely invasive tumors are associated with greater volumetric response to TMZ suggests patients with these tumors may benefit from additional adjuvant TMZ cycles, even for those without MGMT methylation.

Performance of Standardized Relative CBV for Quantifying Regional Histologic Tumor Burden in Recurrent High-Grade Glioma: Comparison against Normalized Relative CBV Using Image-Localized Stereotactic Biopsies

BACKGROUND AND PURPOSE

Perfusion MR imaging measures of relative CBV can distinguish recurrent tumor from posttreatment radiation effects in high-grade gliomas. Currently, relative CBV measurement requires normalization based on user-defined reference tissues. A recently proposed method of relative CBV standardization eliminates the need for user input. This study compares the predictive performance of relative CBV standardization against relative CBV normalization for quantifying recurrent tumor burden in high-grade gliomas relative to posttreatment radiation effects.

MATERIALS AND METHODS

We recruited 38 previously treated patients with high-grade gliomas (World Health Organization grades III or IV) undergoing surgical re-resection for new contrast-enhancing lesions concerning for recurrent tumor versus posttreatment radiation effects. We recovered 112 image-localized biopsies and quantified the percentage of histologic tumor content versus posttreatment radiation effects for each sample. We measured spatially matched normalized and standardized relative CBV metrics (mean, median) and fractional tumor burden for each biopsy. We compared relative CBV performance to predict tumor content, including the Pearson correlation (r), against histologic tumor content (0%-100%) and the receiver operating characteristic area under the curve for predicting high-versus-low tumor content using binary histologic cutoffs (≥50%; ≥80% tumor).

RESULTS

Across relative CBV metrics, fractional tumor burden showed the highest correlations with tumor content (0%-100%) for normalized (r = 0.63, P < .001) and standardized (r = 0.66, P < .001) values. With binary cutoffs (ie, ≥50%; ≥80% tumor), predictive accuracies were similar for both standardized and normalized metrics and across relative CBV metrics. Median relative CBV achieved the highest area under the curve (normalized = 0.87, standardized = 0.86) for predicting ≥50% tumor, while fractional tumor burden achieved the highest area under the curve (normalized = 0.77, standardized = 0.80) for predicting ≥80% tumor.

CONCLUSIONS

Standardization of relative CBV achieves similar performance compared with normalized relative CBV and offers an important step toward workflow optimization and consensus methodology.

Imaging of intratumoral heterogeneity in high-grade glioma

High-grade glioma (HGG), and particularly Glioblastoma (GBM), can exhibit pronounced intratumoral heterogeneity that confounds clinical diagnosis and management. While conventional contrast-enhanced MRI lacks the capability to resolve this heterogeneity, advanced MRI techniques and PET imaging offer a spectrum of physiologic and biophysical image features to improve the specificity of imaging diagnoses. Published studies have shown how integrating these advanced techniques can help better define histologically distinct targets for surgical and radiation treatment planning, and help evaluate the regional heterogeneity of tumor recurrence and response assessment following standard adjuvant therapy. Application of texture analysis and machine learning (ML) algorithms has also enabled the emerging field of radiogenomics, which can spatially resolve the regional and genetically distinct subpopulations that coexist within a single GBM tumor. This review focuses on the latest advances in neuro-oncologic imaging and their clinical applications for the assessment of intratumoral heterogeneity.

Days gained response discriminates treatment response in patients with recurrent glioblastoma receiving bevacizumab-based therapies

Background

Accurate assessments of patient response to therapy are a critical component of personalized medicine. In glioblastoma (GBM), the most aggressive form of brain cancer, tumor growth dynamics are heterogenous across patients, complicating assessment of treatment response. This study aimed to analyze days gained (DG), a burgeoning model-based dynamic metric, for response assessment in patients with recurrent GBM who received bevacizumab-based therapies.

Methods

DG response scores were calculated using volumetric tumor segmentations for patients receiving bevacizumab with and without concurrent cytotoxic therapy (N = 62). Kaplan–Meier and Cox proportional hazards analyses were implemented to examine DG prognostic relationship to overall (OS) and progression-free survival (PFS) from the onset of treatment for recurrent GBM.

Results

In patients receiving concurrent bevacizumab and cytotoxic therapy, Kaplan–Meier analysis showed significant differences in OS and PFS at DG cutoffs consistent with previously identified values from newly diagnosed GBM using T1-weighted gadolinium-enhanced magnetic resonance imaging (T1Gd). DG scores for bevacizumab monotherapy patients only approached significance for PFS. Cox regression showed that increases of 25 DG on T1Gd imaging were significantly associated with a 12.5% reduction in OS hazard for concurrent therapy patients and a 4.4% reduction in PFS hazard for bevacizumab monotherapy patients.

Conclusion

DG has significant meaning in recurrent therapy as a metric of treatment response, even in the context of anti-angiogenic therapies. This provides further evidence supporting the use of DG as an adjunct response metric that quantitatively connects treatment response and clinical outcomes.

ENvironmental Dynamics Underlying Responsive Extreme Survivors (ENDURES) of Glioblastoma: A Multidisciplinary Team-based, Multifactorial Analytical Approach

Although glioblastoma (GBM) is a fatal primary brain cancer with short median survival of 15 months, a small number of patients survive >5 years after diagnosis; they are known as extreme survivors (ES). Because of their rarity, very little is known about what differentiates these outliers from other patients with GBM. For the purpose of identifying unknown drivers of extreme survivorship in GBM, the ENDURES consortium (ENvironmental Dynamics Underlying Responsive Extreme Survivors of GBM) was developed. This consortium is a multicenter collaborative network of investigators focused on the integration of multiple types of clinical data and the creation of patient-specific models of tumor growth informed by radiographic and histologic parameters. Leveraging our combined resources, the goals of the ENDURES consortium are 2-fold: (1) to build a curated, searchable, multilayered repository housing clinical and outcome data on a large cohort of ES patients with GBM; and (2) to leverage the ENDURES repository for new insights into tumor behavior and novel targets for prolonging survival for all patients with GBM. In this article, the authors review the available literature and discuss what is already known about ES. The authors then describe the creation of their consortium and some preliminary results.

Integration of machine learning and mechanistic models accurately predicts variation in cell density of glioblastoma using multiparametric MRI

Glioblastoma (GBM) is a heterogeneous and lethal brain cancer. These tumors are followed using magnetic resonance imaging (MRI), which is unable to precisely identify tumor cell invasion, impairing effective surgery and radiation planning. We present a novel hybrid model, based on multiparametric intensities, which combines machine learning (ML) with a mechanistic model of tumor growth to provide spatially resolved tumor cell density predictions. The ML component is an imaging data-driven graph-based semi-supervised learning model and we use the Proliferation-Invasion (PI) mechanistic tumor growth model. We thus refer to the hybrid model as the ML-PI model. The hybrid model was trained using 82 image-localized biopsies from 18 primary GBM patients with pre-operative MRI using a leave-one-patient-out cross validation framework. A Relief algorithm was developed to quantify relative contributions from the data sources. The ML-PI model statistically significantly outperformed (p < 0.001) both individual models, ML and PI, achieving a mean absolute predicted error (MAPE) of 0.106 ± 0.125 versus 0.199 ± 0.186 (ML) and 0.227 ± 0.215 (PI), respectively. Associated Pearson correlation coefficients for ML-PI, ML, and PI were 0.838, 0.518, and 0.437, respectively. The Relief algorithm showed the PI model had the greatest contribution to the result, emphasizing the importance of the hybrid model in achieving the high accuracy.

Moving Toward a Consensus DSC-MRI Protocol: Validation of a Low-Flip Angle Single-Dose Option as a Reference Standard for Brain Tumors

BACKGROUND AND PURPOSE

DSC-MR imaging using preload, intermediate (60°) flip angle and postprocessing leakage correction has gained traction as a standard methodology. Simulations suggest that DSC-MR imaging with flip angle = 30° and no preload yields relative CBV practically equivalent to the reference standard. This study tested this hypothesis in vivo.

MATERIALS AND METHODS

Eighty-four patients with brain lesions were enrolled in this 3-institution study. Forty-three patients satisfied the inclusion criteria. DSC-MR imaging (3T, single-dose gadobutrol, gradient recalled-echo-EPI, TE = 20-35 ms, TR = 1.2-1.63 seconds) was performed twice for each patient, with flip angle = 30°-35° and no preload (P-), which provided preload (P+) for the subsequent intermediate flip angle = 60°. Normalized relative CBV and standardized relative CBV maps were generated, including postprocessing with contrast agent leakage correction (C+) and without (C-) contrast agent leakage correction. Contrast-enhancing lesion volume, mean relative CBV, and contrast-to-noise ratio obtained with 30°/P-/C-, 30°/P-/C+, and 60°/P+/C- were compared with 60°/P+/C+ using the Lin concordance correlation coefficient and Bland-Altman analysis. Equivalence between the 30°/P-/C+ and 60°/P+/C+ protocols and the temporal SNR for the 30°/P- and 60°/P+ DSC-MR imaging data was also determined.

RESULTS

Compared with 60°/P+/C+, 30°/P-/C+ had closest mean standardized relative CBV (P = .61), highest Lin concordance correlation coefficient (0.96), and lowest Bland-Altman bias (μ = 1.89), compared with 30°/P-/C- (P = .02, Lin concordance correlation coefficient = 0.59, μ = 14.6) and 60°/P+/C- (P = .03, Lin concordance correlation coefficient = 0.88, μ = -10.1) with no statistical difference in contrast-to-noise ratios across protocols. The normalized relative CBV and standardized relative CBV were statistically equivalent at the 10% level using either the 30°/P-/C+ or 60°/P+/C+ protocols. Temporal SNR was not significantly different for 30°/P- and 60°/P+ (P = .06).

CONCLUSIONS

Tumor relative CBV derived from low-flip angle, no-preload DSC-MR imaging with leakage correction is an attractive single-dose alternative to the higher dose reference standard.

Evaluating Multisite rCBV Consistency from DSC-MRI Imaging Protocols and Postprocessing Software Across the NCI Quantitative Imaging Network Sites Using a Digital Reference Object (DRO)

Relative cerebral blood volume (rCBV) cannot be used as a response metric in clinical trials, in part, because of variations in biomarker consistency and associated interpretation across sites, stemming from differences in image acquisition and postprocessing methods (PMs). This study leveraged a dynamic susceptibility contrast magnetic resonance imaging digital reference object to characterize rCBV consistency across 12 sites participating in the Quantitative Imaging Network (QIN), specifically focusing on differences in site-specific imaging protocols (IPs; n = 17), and PMs (n = 19) and differences due to site-specific IPs and PMs (n = 25). Thus, high agreement across sites occurs when 1 managing center processes rCBV despite slight variations in the IP. This result is most likely supported by current initiatives to standardize IPs. However, marked intersite disagreement was observed when site-specific software was applied for rCBV measurements. This study's results have important implications for comparing rCBV values across sites and trials, where variability in PMs could confound the comparison of therapeutic effectiveness and/or any attempts to establish thresholds for categorical response to therapy. To overcome these challenges and ensure the successful use of rCBV as a clinical trial biomarker, we recommend the establishment of qualifying and validating site- and trial-specific criteria for scanners and acquisition methods (eg, using a validated phantom) and the software tools used for dynamic susceptibility contrast magnetic resonance imaging analysis (eg, using a digital reference object where the ground truth is known).

Accurate Patient-Specific Machine Learning Models of Glioblastoma Invasion Using Transfer Learning

BACKGROUND AND PURPOSE

MR imaging-based modeling of tumor cell density can substantially improve targeted treatment of glioblastoma. Unfortunately, interpatient variability limits the predictive ability of many modeling approaches. We present a transfer learning method that generates individualized patient models, grounded in the wealth of population data, while also detecting and adjusting for interpatient variabilities based on each patient's own histologic data.

MATERIALS AND METHODS

We recruited patients with primary glioblastoma undergoing image-guided biopsies and preoperative imaging, including contrast-enhanced MR imaging, dynamic susceptibility contrast MR imaging, and diffusion tensor imaging. We calculated relative cerebral blood volume from DSC-MR imaging and mean diffusivity and fractional anisotropy from DTI. Following image coregistration, we assessed tumor cell density for each biopsy and identified corresponding localized MR imaging measurements. We then explored a range of univariate and multivariate predictive models of tumor cell density based on MR imaging measurements in a generalized one-model-fits-all approach. We then implemented both univariate and multivariate individualized transfer learning predictive models, which harness the available population-level data but allow individual variability in their predictions. Finally, we compared Pearson correlation coefficients and mean absolute error between the individualized transfer learning and generalized one-model-fits-all models.

RESULTS

Tumor cell density significantly correlated with relative CBV (r = 0.33, P < .001), and T1-weighted postcontrast (r = 0.36, P < .001) on univariate analysis after correcting for multiple comparisons. With single-variable modeling (using relative CBV), transfer learning increased predictive performance (r = 0.53, mean absolute error = 15.19%) compared with one-model-fits-all (r = 0.27, mean absolute error = 17.79%). With multivariate modeling, transfer learning further improved performance (r = 0.88, mean absolute error = 5.66%) compared with one-model-fits-all (r = 0.39, mean absolute error = 16.55%).

CONCLUSIONS

Transfer learning significantly improves predictive modeling performance for quantifying tumor cell density in glioblastoma.

Is the blood-brain barrier really disrupted in all glioblastomas? A critical assessment of existing clinical data

The blood-brain barrier (BBB) excludes the vast majority of cancer therapeutics from normal brain. However, the importance of the BBB in limiting drug delivery and efficacy is controversial in high-grade brain tumors, such as glioblastoma (GBM). The accumulation of normally brain impenetrant radiographic contrast material in essentially all GBM has popularized a belief that the BBB is uniformly disrupted in all GBM patients so that consideration of drug distribution across the BBB is not relevant in designing therapies for GBM. However, contrary to this view, overwhelming clinical evidence demonstrates that there is also a clinically significant tumor burden with an intact BBB in all GBM, and there is little doubt that drugs with poor BBB permeability do not provide therapeutically effective drug exposures to this fraction of tumor cells. This review provides an overview of the clinical literature to support a central hypothesis: that all GBM patients have tumor regions with an intact BBB, and cure for GBM will only be possible if these regions of tumor are adequately treated.

Optimization of Acquisition and Analysis Methods for Clinical Dynamic Susceptibility Contrast MRI Using a Population-Based Digital Reference Object

BACKGROUND AND PURPOSE

The accuracy of DSC-MR imaging CBV maps in glioblastoma depends on acquisition and analysis protocols. Multisite protocol heterogeneity has challenged standardization initiatives due to the difficulties of in vivo validation. This study sought to compare the accuracy of routinely used protocols using a digital reference object.

MATERIALS AND METHODS

The digital reference object consisted of approximately 10,000 simulated voxels recapitulating typical signal heterogeneity encountered in vivo. The influence of acquisition and postprocessing methods on CBV reliability was evaluated across 6912 parameter combinations, including contrast agent dosing schemes, pulse sequence parameters, field strengths, and postprocessing methods. Accuracy and precision were assessed using the concordance correlation coefficient and coefficient of variation.

RESULTS

Across all parameter space, the optimal protocol included full-dose contrast agent preload and bolus, intermediate (60°) flip angle, 30-ms TE, and postprocessing with a leakage-correction algorithm (concordance correlation coefficient = 0.97, coefficient of variation = 6.6%). Protocols with no preload or fractional dose preload and bolus using these acquisition parameters were generally less robust. However, a protocol with no preload, full-dose bolus, and low (30°) flip angle performed very well (concordance correlation coefficient = 0.93, coefficient of variation = 8.7% at 1.5T and concordance correlation coefficient = 0.92, coefficient of variation = 8.2% at 3T).

CONCLUSIONS

Schemes with full-dose preload and bolus maximize CBV accuracy and reduce variability, which could enable smaller sample sizes and more reliable detection of CBV changes in clinical trials. When a lower total contrast agent dose is desired, use of a low flip angle, no preload, and full-dose bolus protocol may provide an attractive alternative.

Neuroimaging abnormalities in patients with Cowden syndrome: Retrospective single-center study

Background

We retrospectively reviewed the neuroimaging findings of patients with Cowden syndrome and determined their frequency in a single cohort.

Methods

Electronic medical records were queried from January 1999 to January 2017 to identify patients who fit the clinical criteria for diagnosis of Cowden syndrome with or without a documented PTEN mutation. Patients with brain MRI examinations were then identified.

Results

We retrospectively identified 44 patients with Cowden syndrome, 22 of whom had neuroimaging for review. Eleven (50%) had Lhermitte-Duclos disease, 4 (18.1%) had meningiomas, 13 (59.1%) had at least one developmental venous anomaly, 3 had cavernous malformations, 2 had evidence of dural arteriovenous fistula, 7 had increased white matter signal abnormalities relative to age (31.8%), 4 had prominent perivascular spaces, cerebellar tonsillar ectopia was present in 7 of 21 (33.3%), and 1 had cortical malformation.

Conclusions

It is important to recognize that in addition to Lhermitte-Duclos disease, other intracranial findings such as multiple venous anomalies, meningiomas, greater than expected white matter signal abnormality, prominent perivascular spaces, and cortical malformations may warrant a thorough evaluation for Cowden syndrome in the appropriate clinical setting. We further recommend that this broader spectrum of intracranial abnormalities be considered for addition to the Cowden syndrome diagnostic criteria at the time of next revision.

Characterizing the Influence of Preload Dosing on Percent Signal Recovery (PSR) and Cerebral Blood Volume (CBV) Measurements in a Patient Population With High-Grade Glioma Using Dynamic Susceptibility Contrast MRI

With DSC-MRI, contrast agent leakage effects in brain tumors can either be leveraged for percent signal recovery (PSR) measurements or be adequately resolved for accurate relative cerebral blood volume (rCBV) measurements. Leakage effects can be dimished by administration of a preload dose before imaging and/or specific postprocessing steps. This study compares the consistency of both PSR and rCBV measurements as a function of varying preload doses in a retrospective analysis of 14 subjects with high-grade gliomas. The scans consisted of 6 DSC-MRI scans during 6 sequential bolus injections (0.05 mmol/kg). Mean PSR was calculated for tumor and normal-appearing white matter regions of interest. DSC-MRI data were corrected for leakage effects before computing mean tumor rCBV. Statistical differences were seen across varying preloads for tumor PSR (P value = 4.57E-24). Tumor rCBV values did not exhibit statistically significant differences across preloads (P value = .14) and were found to be highly consistent for clinically relevant preloads (intraclass correlation coefficient = 0.93). For a 0.05 mmol/kg injection bolus and pulse sequence parameters used, the highest PSR contrast between normal-appearing white matter and tumor occurs when no preload is used. This suggests that studies using PSR as a biomarker should acquire DSC-MRI data without preload. The finding that leakage-corrected rCBV values do not depend on the presence or dose of preload contradicts that of previous studies with dissimilar acquisition protocols. This further confirms the sensitivity of rCBV to preload dosing schemes and pulse sequence parameters and highlights the importance of standardization efforts for achieving multisite rCBV consistency.

Dynamic Susceptibility Contrast-MRI Quantification Software Tool: Development and Evaluation

Relative cerebral blood volume (rCBV) is a magnetic resonance imaging biomarker that is used to differentiate progression from pseudoprogression in patients with glioblastoma multiforme, the most common primary brain tumor. However, calculated rCBV depends considerably on the software used. Automating all steps required for rCBV calculation is important, as user interaction can lead to increased variability and possible inaccuracies in clinical decision-making. Here, we present an automated tool for computing rCBV from dynamic susceptibility contrast-magnetic resonance imaging that includes leakage correction. The entrance and exit bolus time points are automatically calculated using wavelet-based detection. The proposed tool is compared with 3 Food and Drug Administration-approved software packages, 1 automatic and 2 requiring user interaction, on a data set of 43 patients. We also evaluate manual and automated white matter (WM) selection for normalization of the cerebral blood volume maps. Our system showed good agreement with 2 of the 3 software packages. The intraclass correlation coefficient for all comparisons between the same software operated by different people was >0.880, except for FuncTool when operated by user 1 versus user 2. Little variability in agreement between software tools was observed when using different WM selection techniques. Our algorithm for automatic rCBV calculation with leakage correction and automated WM selection agrees well with 2 out of the 3 FDA-approved software packages.

Radiogenomics to characterize regional genetic heterogeneity in glioblastoma

BACKGROUND

Glioblastoma (GBM) exhibits profound intratumoral genetic heterogeneity. Each tumor comprises multiple genetically distinct clonal populations with different therapeutic sensitivities. This has implications for targeted therapy and genetically informed paradigms. Contrast-enhanced (CE)-MRI and conventional sampling techniques have failed to resolve this heterogeneity, particularly for nonenhancing tumor populations. This study explores the feasibility of using multiparametric MRI and texture analysis to characterize regional genetic heterogeneity throughout MRI-enhancing and nonenhancing tumor segments.

METHODS

We collected multiple image-guided biopsies from primary GBM patients throughout regions of enhancement (ENH) and nonenhancing parenchyma (so called brain-around-tumor, [BAT]). For each biopsy, we analyzed DNA copy number variants for core GBM driver genes reported by The Cancer Genome Atlas. We co-registered biopsy locations with MRI and texture maps to correlate regional genetic status with spatially matched imaging measurements. We also built multivariate predictive decision-tree models for each GBM driver gene and validated accuracies using leave-one-out-cross-validation (LOOCV).

RESULTS

We collected 48 biopsies (13 tumors) and identified significant imaging correlations (univariate analysis) for 6 driver genes: EGFR, PDGFRA, PTEN, CDKN2A, RB1, and TP53. Predictive model accuracies (on LOOCV) varied by driver gene of interest. Highest accuracies were observed for PDGFRA (77.1%), EGFR (75%), CDKN2A (87.5%), and RB1 (87.5%), while lowest accuracy was observed in TP53 (37.5%). Models for 4 driver genes (EGFR, RB1, CDKN2A, and PTEN) showed higher accuracy in BAT samples (n = 16) compared with those from ENH segments (n = 32).

CONCLUSION

MRI and texture analysis can help characterize regional genetic heterogeneity, which offers potential diagnostic value under the paradigm of individualized oncology.

Genetic diversity in the germplasm of black pepper determined by EST-SSR markers

This study aimed to assess genetic diversity in the germplasm of black pepper from around the world using SSR markers from EST. In total, 13 markers were selected and successfully amplified the target loci across the black pepper germplasm. All the EST-SSR markers showed high levels of polymorphisms with an average polymorphism information content of 0.93. The genetic similarity coefficients among all accessions ranged from 0.724 to 1.000, with an average of 0.867. These results indicated that black pepper germplasms possess a complex genetic background and high genetic diversity. Based on a cluster analysis, 148 black pepper germplasms were grouped in two major clades: the Neotropics and the Asian tropics. Peperomia pellucida was grouped separately and distantly from all other accessions. These results generally agreed with the genetic and geographic distances. However, the Asian tropics clade did not cluster according to their geographic origins. In addition, compared with the American accessions, the Asian wild accessions and cultivated accessions grouped together, indicating a close genetic relationship. This verified the origin of black pepper. The newly developed EST-SSRs are highly valuable resources for the conservation of black pepper germplasm diversity and for black pepper breeding.

Multi-Parametric MRI and Texture Analysis to Visualize Spatial Histologic Heterogeneity and Tumor Extent in Glioblastoma

Background

Genetic profiling represents the future of neuro-oncology but suffers from inadequate biopsies in heterogeneous tumors like Glioblastoma (GBM). Contrast-enhanced MRI (CE-MRI) targets enhancing core (ENH) but yields adequate tumor in only ~60% of cases. Further, CE-MRI poorly localizes infiltrative tumor within surrounding non-enhancing parenchyma, or brain-around-tumor (BAT), despite the importance of characterizing this tumor segment, which universally recurs. In this study, we use multiple texture analysis and machine learning (ML) algorithms to analyze multi-parametric MRI, and produce new images indicating tumor-rich targets in GBM.

Methods

We recruited primary GBM patients undergoing image-guided biopsies and acquired pre-operative MRI: CE-MRI, Dynamic-Susceptibility-weighted-Contrast-enhanced-MRI, and Diffusion Tensor Imaging. Following image coregistration and region of interest placement at biopsy locations, we compared MRI metrics and regional texture with histologic diagnoses of high- vs low-tumor content (≥80% vs <80% tumor nuclei) for corresponding samples. In a training set, we used three texture analysis algorithms and three ML methods to identify MRI-texture features that optimized model accuracy to distinguish tumor content. We confirmed model accuracy in a separate validation set.

Results

We collected 82 biopsies from 18 GBMs throughout ENH and BAT. The MRI-based model achieved 85% cross-validated accuracy to diagnose high- vs low-tumor in the training set (60 biopsies, 11 patients). The model achieved 81.8% accuracy in the validation set (22 biopsies, 7 patients).

Conclusion

Multi-parametric MRI and texture analysis can help characterize and visualize GBM’s spatial histologic heterogeneity to identify regional tumor-rich biopsy targets.

Impact of Software Modeling on the Accuracy of Perfusion MRI in Glioma

BACKGROUND AND PURPOSE

Relative cerebral blood volume, as measured by T2*-weighted dynamic susceptibility-weighted contrast-enhanced MRI, represents the most robust and widely used perfusion MR imaging metric in neuro-oncology. Our aim was to determine whether differences in modeling implementation will impact the correction of leakage effects (from blood-brain barrier disruption) and the accuracy of relative CBV calculations as measured on T2*-weighted dynamic susceptibility-weighted contrast-enhanced MR imaging at 3T field strength.

MATERIALS AND METHODS

This study included 52 patients with glioma undergoing DSC MR imaging. Thirty-six patients underwent both non-preload dose- and preload dose-corrected DSC acquisitions, with 16 patients undergoing preload dose-corrected acquisitions only. For each acquisition, we generated 2 sets of relative CBV metrics by using 2 separate, widely published, FDA-approved commercial software packages: IB Neuro and nordicICE. We calculated 4 relative CBV metrics within tumor volumes: mean relative CBV, mode relative CBV, percentage of voxels with relative CBV > 1.75, and percentage of voxels with relative CBV > 1.0 (fractional tumor burden). We determined Pearson (r) and Spearman (ρ) correlations between non-preload dose- and preload dose-corrected metrics. In a subset of patients with recurrent glioblastoma (n = 25), we determined receiver operating characteristic area under the curve for fractional tumor burden accuracy to predict the tissue diagnosis of tumor recurrence versus posttreatment effect. We also determined correlations between rCBV and microvessel area from stereotactic biopsies (n = 29) in 12 patients.

RESULTS

With IB Neuro, relative CBV metrics correlated highly between non-preload dose- and preload dose-corrected conditions for fractional tumor burden (r = 0.96, ρ = 0.94), percentage > 1.75 (r = 0.93, ρ = 0.91), mean (r = 0.87, ρ = 0.86), and mode (r = 0.78, ρ = 0.76). These correlations dropped substantially with nordicICE. With fractional tumor burden, IB Neuro was more accurate than nordicICE in diagnosing tumor versus posttreatment effect (area under the curve = 0.85 versus 0.67) (P < .01). The highest relative CBV-microvessel area correlations required preload dose and IB Neuro (r = 0.64, ρ = 0.58, P = .001).

CONCLUSIONS

Different implementations of perfusion MR imaging software modeling can impact the accuracy of leakage correction, relative CBV calculation, and correlations with histologic benchmarks.

Biopsy validation of 18F-DOPA PET and biodistribution in gliomas for neurosurgical planning and radiotherapy target delineation: results of a prospective pilot study

BACKGROUND

Delineation of glioma extent for surgical or radiotherapy planning is routinely based on MRI. There is increasing awareness that contrast enhancement on T1-weighted images (T1-CE) may not reflect the entire extent of disease. The amino acid tracer (18)F-DOPA (3,4-dihydroxy-6-[18F] fluoro-l-phenylalanine) has a high tumor-to-background signal and high sensitivity for glioma imaging. This study compares (18)F-DOPA PET against conventional MRI for neurosurgical biopsy targeting, resection planning, and radiotherapy target volume delineation.

METHODS

Conventional MR and (18)F-DOPA PET/CT images were acquired in 10 patients with suspected malignant brain tumors. One to 3 biopsy locations per patient were chosen in regions of concordant and discordant (18)F-DOPA uptake and MR contrast enhancement. Histopathology was reviewed on 23 biopsies. (18)F-DOPA PET was quantified using standardized uptake values (SUV) and tumor-to-normal hemispheric tissue (T/N) ratios.

RESULTS

Pathologic review confirmed glioma in 22 of 23 biopsy specimens. Thirteen of 16 high-grade biopsy specimens were obtained from regions of elevated (18)F-DOPA uptake, while T1-CE was present in only 6 of those 16 samples. Optimal (18)F-DOPA PET thresholds corresponding to high-grade disease based on histopathology were calculated as T/N > 2.0. In every patient, (18)F-DOPA uptake regions with T/N > 2.0 extended beyond T1-CE up to a maximum of 3.5 cm. SUV was found to correlate with grade and cellularity.

CONCLUSIONS

(18)F-DOPA PET SUV(max) may more accurately identify regions of higher-grade/higher-density disease in patients with astrocytomas and will have utility in guiding stereotactic biopsy selection. Using SUV-based thresholds to define high-grade portions of disease may be valuable in delineating radiotherapy boost volumes.