Wu, Teresa

Background The current subclassification of migraine is according to headache frequency and aura status. The variability in migraine symptoms, disease course, and response to treatment suggest the presence of additional heterogeneity or subclasses within migraine. Objective The study objective was to subclassify migraine via a data‐driven approach, identifying latent factors by jointly exploiting multiple sets of brain structural features obtained via magnetic resonance imaging (MRI). Methods Migraineurs ( n = 66) and healthy controls ( n = 54) had brain MRI measurements of cortical thickness, cortical surface area, and volumes for 68 regions. A multimodality factor mixture model was used to subclassify MRIs and to determine the brain structural factors that most contributed to the subclassification. Clinical characteristics of subjects in each subgroup were compared. Results Automated MRI classification divided the subjects into two subgroups. Migraineurs in subgroup #1 had more severe allodynia symptoms during migraines (6.1 ± 5.3 vs. 3.6 ± 3.2, P = .03), more years with migraine (19.2 ± 11.3 years vs 13 ± 8.3 years, P = .01), and higher Migraine Disability Assessment (MIDAS) scores (25 ± 22.9 vs 15.7 ± 12.2, P = .04). There were not differences in headache frequency or migraine aura status between the two subgroups. Conclusions Data‐driven subclassification of brain MRIs based upon structural measurements identified two subgroups. Amongst migraineurs, the subgroups differed in allodynia symptom severity, years with migraine, and migraine‐related disability. Since allodynia is associated with this imaging‐based subclassification of migraine and prior publications suggest that allodynia impacts migraine treatment response and disease prognosis, future migraine diagnostic criteria could consider allodynia when defining migraine subgroups.

Abstract Purpose: Investigate the utility of patient-specific spatial predictions of tumor cell density from a bio-mathematical model. Introduction: Glioblastomas (GBMs) are the most malignant of all primary brain tumors. While it is known there is always a non-detectable portion of the tumor, current techniques of monitoring GBM progression, imaging and initial histological assessment, are not able to reliably estimate the tumor invasion past the enhancing region on T2-Weighted (T2W) imaging. Over the last two decades, a large effort has been made to create a simple patient-specific mathematical model of gliomas. The resulting model, referred to as the Proliferation-Invasion (PI) model, is based on two key parameters, the net growth rate, ρ, and the dispersal coefficient, D. In this model, the ratio of D/ρ is related to degree of invasion and the product D*ρ, is related to the speed of growth. The intuitive understanding provided by this model has been able to provide patient-specific understanding of disease kinetics enabling prediction of outcomes following surgical resection, radiation and the development of a prognostic response metric. Previous literature utilizing this model has been based on the assumption that what is seen on the pretreatment T1-Weighted contrast-enhanced (T1Gd) and T2W, images correspond to an 80% and 16% tumor cell density threshold respectively. This assumption allows for an estimate of D/ρ from a single time point of imaging. While these values were based on extensive experience, for ethical and technical reasons, they have never been rigorously investigated histologically. Recent technological advances have made it possible for surgeons to use an MRI to guide the acquisition of tissue making it possible to know with a good degree of accuracy where on the MR image the histological specimen comes from. Methods: Model Calibration : To estimate D/ρ for each patient, we assume abnormalities on the T1Gd and T2W images correspond to an 80% and 16% tumor cell density threshold respectively. We then utilize a Bayesian calibration approach based on adaptive grid refinement while holding the velocity constant to find the most likely value of D/ρ to match the observed radial measurements. Three-Dimensional Density Maps : Given a gray/white segmentation and an estimate for D/ρ, we can build a tumor cell density prediction in the patient's anatomy using the Eikonal equations and the modified Fast Marching Method (FMM) algorithm presented by Konukoglu et al. Patient Cohort : Eighteen patients were recruited with clinically suspected GBM undergoing preoperative stereotactic MRI for surgical resection with IRB approval Barrow Neurological Institute and Mayo Clinic in Arizona. Surgical Biopsy : Pre-operative conventional MRI, including T1Gd and T2W, was utilized to guide stereotactic biopsies. An average of 5–6 tissue specimens were acquired from each tumor by using stereotactic surgical localization, following the smallest possible diameter craniotomies to minimize brain shift. Histological Analysis : 4 μm tissue sections were stained with hematoxylin and eosin (H&E) for neuropathology review. H&E slides were reviewed by a neuropathologist. The percent tumor nuclei in the field was estimated. Results and Conclusion: Our analysis showed that the intuitive ordering of diffuse to nodular tumor profiles from the PI model is consistent with the histologic data. Further, utilizing the histologic data, limitations to a universal threshold cutoff of the T2/FLAIR region are specified and a better threshold cutoff value is suggested. We hope that this paper will lead to more studies of this kind resulting in more accurate ways to predict the spatial distribution of GBM tumor cells from medical imaging. Such advances could revolutionize surgical procedures and radiation therapy and may enable additional insights into tumor kinetic differences meaningful for treatment. Citation Format: Andrea Hawkins-Daarud, Lauren DeGirolamo, Joshua Jacobs, Kamala Clark-Swanson, Jennifer M. Eschbacher, Kris A. Smith, Peter Nakaji, Leslie C. Baxter, John P. Karis, Teresa Wu, J. Ross Mitchell, Jing Li, Leland Hu, Kristin R. Swanson. Histologic evidence for a bio-mathematical model of glioblastoma invasion. [abstract]. In: Proceedings of the AACR Special Conference on Engineering and Physical Sciences in Oncology; 2016 Jun 25-28; Boston, MA. Philadelphia (PA): AACR; Cancer Res 2017;77(2 Suppl):Abstract nr A08.

AbstractBackgroundGlioblastoma (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.MethodsWe 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).ResultsWe 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).ConclusionMRI and texture analysis can help characterize regional genetic heterogeneity, which offers potential diagnostic value under the paradigm of individualized oncology.

Background This study used machine-learning techniques to develop discriminative brain-connectivity biomarkers from resting-state functional magnetic resonance neuroimaging ( rs-fMRI) data that distinguish between individual migraine patients and healthy controls. Methods This study included 58 migraine patients (mean age = 36.3 years; SD = 11.5) and 50 healthy controls (mean age = 35.9 years; SD = 11.0). The functional connections of 33 seeded pain-related regions were used as input for a brain classification algorithm that tested the accuracy of determining whether an individual brain MRI belongs to someone with migraine or to a healthy control. Results The best classification accuracy using a 10-fold cross-validation method was 86.1%. Resting functional connectivity of the right middle temporal, posterior insula, middle cingulate, left ventromedial prefrontal and bilateral amygdala regions best discriminated the migraine brain from that of a healthy control. Migraineurs with longer disease durations were classified more accurately (>14 years; 96.7% accuracy) compared to migraineurs with shorter disease durations (≤14 years; 82.1% accuracy). Conclusions Classification of migraine using rs-fMRI provides insights into pain circuits that are altered in migraine and could potentially contribute to the development of a new, noninvasive migraine biomarker. Migraineurs with longer disease burden were classified more accurately than migraineurs with shorter disease burden, potentially indicating that disease duration leads to reorganization of brain circuitry.

Background The International Classification of Headache Disorders provides criteria for the diagnosis and subclassification of migraine. Since there is no objective gold standard by which to test these diagnostic criteria, the criteria are based on the consensus opinion of content experts. Accurate migraine classifiers consisting of brain structural measures could serve as an objective gold standard by which to test and revise diagnostic criteria. The objectives of this study were to utilize magnetic resonance imaging measures of brain structure for constructing classifiers: (1) that accurately identify individuals as having chronic vs episodic migraine vs being a healthy control; and (2) that test the currently used threshold of 15 headache days/month for differentiating chronic migraine from episodic migraine. Methods Study participants underwent magnetic resonance imaging for determination of regional cortical thickness, cortical surface area, and volume. Principal components analysis combined structural measurements into principal components accounting for 85% of variability in brain structure. Models consisting of these principal components were developed to achieve the classification objectives. Tenfold cross validation assessed classification accuracy within each of the 10 runs, with data from 90% of participants randomly selected for classifier development and data from the remaining 10% of participants used to test classification performance. Headache frequency thresholds ranging from 5‐15 headache days/month were evaluated to determine the threshold allowing for the most accurate subclassification of individuals into lower and higher frequency subgroups. Results Participants were 66 migraineurs and 54 healthy controls, 75.8% female, with an average age of 36 +/− 11 years. Average classifier accuracies were: (1) 68% for migraine (episodic + chronic) vs healthy controls; (2) 67.2% for episodic migraine vs healthy controls; (3) 86.3% for chronic migraine vs healthy controls; and (4) 84.2% for chronic migraine vs episodic migraine. The classifiers contained principal components consisting of several structural measures, commonly including the temporal pole, anterior cingulate cortex, superior temporal lobe, entorhinal cortex, medial orbital frontal gyrus, and pars triangularis. A threshold of 15 headache days/month allowed for the most accurate subclassification of migraineurs into lower frequency and higher frequency subgroups. Conclusions Classifiers consisting of cortical surface area, cortical thickness, and regional volumes were highly accurate for determining if individuals have chronic migraine. Furthermore, results provide objective support for the current use of 15 headache days/month as a threshold for dividing migraineurs into lower frequency (ie, episodic migraine) and higher frequency (ie, chronic migraine) subgroups.