Multivariate spatiotemporal functional principal component analysis for modeling hospitalization and mortality rates in the dialysis population.

Dialysis patients experience frequent hospitalizations and a higher mortality rate compared to other Medicare populations, in whom hospitalizations are a major contributor to morbidity, mortality, and healthcare costs. Patients also typically remain on dialysis for the duration of their lives or until kidney transplantation. Hence, there is growing interest in studying the spatiotemporal trends in the correlated outcomes of hospitalization and mortality among dialysis patients as a function of time starting from transition to dialysis across the United States Utilizing national data from the United States Renal Data System (USRDS), we propose a novel multivariate spatiotemporal functional principal component analysis model to study the joint spatiotemporal patterns of hospitalization and mortality rates among dialysis patients. The proposal is based on a multivariate Karhunen-Loéve expansion that describes leading directions of variation across time and induces spatial correlations among region-specific scores. An efficient estimation procedure is proposed using only univariate principal components decompositions and a Markov Chain Monte Carlo framework for targeting the spatial correlations. The finite sample performance of the proposed method is studied through simulations. Novel applications to the USRDS data highlight hot spots across the United States with higher hospitalization and/or mortality rates and time periods of elevated risk.

Modeling intra-individual inter-trial EEG response variability in autism.

Autism spectrum disorder (autism) is a prevalent neurodevelopmental condition characterized by early emerging impairments in social behavior and communication. EEG represents a powerful and non-invasive tool for examining functional brain differences in autism. Recent EEG evidence suggests that greater intra-individual trial-to-trial variability across EEG responses in stimulus-related tasks may characterize brain differences in autism. Traditional analysis of EEG data largely focuses on mean trends of the trial-averaged data, where trial-level analysis is rarely performed due to low neural signal to noise ratio. We propose to use nonlinear (shape-invariant) mixed effects (NLME) models to study intra-individual inter-trial EEG response variability using trial-level EEG data. By providing more precise metrics of response variability, this approach could enrich our understanding of neural disparities in autism and potentially aid the identification of objective markers. The proposed multilevel NLME models quantify variability in the signal's interpretable and widely recognized features (e.g., latency and amplitude) while also regularizing estimation based on noisy trial-level data. Even though NLME models have been studied for more than three decades, existing methods cannot scale up to large data sets. We propose computationally feasible estimation and inference methods via the use of a novel minorization-maximization (MM) algorithm. Extensive simulations are conducted to show the efficacy of the proposed procedures. Applications to data from a large national consortium find that children with autism have larger intra-individual inter-trial variability in P1 latency in a visual evoked potential (VEP) task, compared to their neurotypical peers.

Bayesian analysis of longitudinal and multidimensional functional data.

Multi-dimensional functional data arises in numerous modern scientific experimental and observational studies. In this article, we focus on longitudinal functional data, a structured form of multidimensional functional data. Operating within a longitudinal functional framework we aim to capture low dimensional interpretable features. We propose a computationally efficient nonparametric Bayesian method to simultaneously smooth observed data, estimate conditional functional means and functional covariance surfaces. Statistical inference is based on Monte Carlo samples from the posterior measure through adaptive blocked Gibbs sampling. Several operative characteristics associated with the proposed modeling framework are assessed comparatively in a simulated environment. We illustrate the application of our work in two case studies. The first case study involves age-specific fertility collected over time for various countries. The second case study is an implicit learning experiment in children with autism spectrum disorder.

Multilevel hybrid principal components analysis for region-referenced functional electroencephalography data.

Electroencephalography experiments produce region-referenced functional data representing brain signals in the time or the frequency domain collected across the scalp. The data typically also have a multilevel structure with high-dimensional observations collected across multiple experimental conditions or visits. Common analysis approaches reduce the data complexity by collapsing the functional and regional dimensions, where event-related potential (ERP) features or band power are targeted in a pre-specified scalp region. This practice can fail to portray more comprehensive differences in the entire ERP signal or the power spectral density (PSD) across the scalp. Building on the weak separability of the high-dimensional covariance process, the proposed multilevel hybrid principal components analysis (M-HPCA) utilizes dimension reduction tools from both vector and functional principal components analysis to decompose the total variation into between- and within-subject variance. The resulting model components are estimated in a mixed effects modeling framework via a computationally efficient minorization-maximization algorithm coupled with bootstrap. The diverse array of applications of M-HPCA is showcased with two studies of individuals with autism. While ERP responses to match vs mismatch conditions are compared in an audio odd-ball paradigm in the first study, short-term reliability of the PSD across visits is compared in the second. Finite sample properties of the proposed methodology are studied in extensive simulations.

A study of longitudinal trends in time-frequency transformations of EEG data during a learning experiment.

EEG experiments yield high-dimensional event-related potential (ERP) data in response to repeatedly presented stimuli throughout the experiment. Changes in the high-dimensional ERP signal throughout the duration of an experiment (longitudinally) is the main quantity of interest in learning paradigms, where they represent the learning dynamics. Typical analysis, which can be performed in the time or the frequency domain, average the ERP waveform across all trials, leading to the loss of the potentially valuable longitudinal information in the data. Longitudinal time-frequency transformation of ERP (LTFT-ERP) is proposed to retain information from both the time and frequency domains, offering distinct but complementary information on the underlying cognitive processes evoked, while still retaining the longitudinal dynamics in the ERP waveforms. LTFT-ERP begins by time-frequency transformations of the ERP data, collected across subjects, electrodes, conditions and trials throughout the duration of the experiment, followed by a data driven multidimensional principal components analysis (PCA) approach for dimension reduction. Following projection of the data onto leading directions of variation in the time and frequency domains, longitudinal learning dynamics are modeled within a mixed effects modeling framework. Applications to a learning paradigm in autism depict distinct learning patterns throughout the experiment among children diagnosed with Autism Spectrum Disorder and their typically developing peers. LTFT-ERP time-frequency joint transformations are shown to bring an additional level of specificity to interpretations of the longitudinal learning patterns related to underlying cognitive processes, which is lacking in single domain analysis (in the time or the frequency domain only). Simulation studies show the efficacy of the proposed methodology.

Hybrid principal components analysis for region-referenced longitudinal functional EEG data.

Electroencephalography (EEG) data possess a complex structure that includes regional, functional, and longitudinal dimensions. Our motivating example is a word segmentation paradigm in which typically developing (TD) children, and children with autism spectrum disorder (ASD) were exposed to a continuous speech stream. For each subject, continuous EEG signals recorded at each electrode were divided into one-second segments and projected into the frequency domain via fast Fourier transform. Following a spectral principal components analysis, the resulting data consist of region-referenced principal power indexed regionally by scalp location, functionally across frequencies, and longitudinally by one-second segments. Standard EEG power analyses often collapse information across the longitudinal and functional dimensions by averaging power across segments and concentrating on specific frequency bands. We propose a hybrid principal components analysis for region-referenced longitudinal functional EEG data, which utilizes both vector and functional principal components analyses and does not collapse information along any of the three dimensions of the data. The proposed decomposition only assumes weak separability of the higher-dimensional covariance process and utilizes a product of one dimensional eigenvectors and eigenfunctions, obtained from the regional, functional, and longitudinal marginal covariances, to represent the observed data, providing a computationally feasible non-parametric approach. A mixed effects framework is proposed to estimate the model components coupled with a bootstrap test for group level inference, both geared towards sparse data applications. Analysis of the data from the word segmentation paradigm leads to valuable insights about group-region differences among the TD and verbal and minimally verbal children with ASD. Finite sample properties of the proposed estimation framework and bootstrap inference procedure are further studied via extensive simulations.

Covariate-adjusted region-referenced generalized functional linear model for EEG data.

Electroencephalography (EEG) studies produce region-referenced functional data in the form of EEG signals recorded across electrodes on the scalp. It is of clinical interest to relate the highly structured EEG data to scalar outcomes such as diagnostic status. In our motivating study, resting-state EEG is collected on both typically developing (TD) children and children with autism spectrum disorder (ASD) aged 2 to 12 years old. The peak alpha frequency (PAF), defined as the location of a prominent peak in the alpha frequency band of the spectral density, is an important biomarker linked to neurodevelopment and is known to shift from lower to higher frequencies as children age. To retain the most amount of information from the data, we consider the oscillations in the spectral density within the alpha band, rather than just the peak location, as a functional predictor of diagnostic status (TD vs ASD), adjusted for chronological age. A covariate-adjusted region-referenced generalized functional linear model is proposed for modeling scalar outcomes from region-referenced functional predictors, which utilizes a tensor basis formed from one-dimensional discrete and continuous bases to estimate functional effects across a discrete regional domain while simultaneously adjusting for additional nonfunctional covariates, such as age. The proposed methodology provides novel insights into differences in neural development of TD and ASD children. The efficacy of the proposed methodology is investigated through extensive simulation studies.

Exposure to the BPA-Substitute Bisphenol S Causes Unique Alterations of Germline Function.

Concerns about the safety of Bisphenol A, a chemical found in plastics, receipts, food packaging and more, have led to its replacement with substitutes now found in a multitude of consumer products. However, several popular BPA-free alternatives, such as Bisphenol S, share a high degree of structural similarity with BPA, suggesting that these substitutes may disrupt similar developmental and reproductive pathways. We compared the effects of BPA and BPS on germline and reproductive functions using the genetic model system Caenorhabditis elegans. We found that, similarly to BPA, BPS caused severe reproductive defects including germline apoptosis and embryonic lethality. However, meiotic recombination, targeted gene expression, whole transcriptome and ontology analyses as well as ToxCast data mining all indicate that these effects are partly achieved via mechanisms distinct from BPAs. These findings therefore raise new concerns about the safety of BPA alternatives and the risk associated with human exposure to mixtures.

Segmentation of High-Cost Adults in an Integrated Healthcare System Based on Empirical Clustering of Acute and Chronic Conditions.

BACKGROUND: High-cost patients are a frequent focus of improvement projects based on primary care and other settings. Efforts to characterize high-cost, high-need patients are needed to inform care planning, but such efforts often rely on a priori assumptions, masking underlying complexities of a heterogenous population. OBJECTIVE: To define recognizable subgroups of patients among high-cost adults based on clinical conditions, and describe their survival and future spending. DESIGN: Retrospective observational cohort study. PARTICIPANTS: Within a large integrated delivery system with 2.7 million adult members, we selected the top 1% of continuously enrolled adults with respect to total healthcare expenditures during 2010. MAIN MEASURES: We used latent class analysis to identify clusters of alike patients based on 53 hierarchical condition categories. Prognosis as measured by healthcare spending and survival was assessed through 2014 for the resulting classes of patients. RESULTS: Among 21,183 high-cost adults, seven clinically distinctive subgroups of patients emerged. Classes included end-stage renal disease (12% of high-cost population), cardiopulmonary conditions (17%), diabetes with multiple comorbidities (8%), acute illness superimposed on chronic conditions (11%), conditions requiring highly specialized care (14%), neurologic and catastrophic conditions (5%), and patients with few comorbidities (the largest class, 33%). Over 4 years of follow-up, 6566 (31%) patients died, and survival in the classes ranged from 43 to 88%. Spending regressed to the mean in all classes except the ESRD and diabetes with multiple comorbidities groups. CONCLUSIONS: Data-driven characterization of high-cost adults yielded clinically intuitive classes that were associated with survival and reflected markedly different healthcare needs. Relatively few high-cost patients remain persistently high cost over 4 years. Our results suggest that high-cost patients, while not a monolithic group, can be segmented into few subgroups. These subgroups may be the focus of future work to understand appropriateness of care and design interventions accordingly.

Chronic exposure to inhaled, traffic-related nitrogen dioxide and a blunted cortisol response in adolescents.

BACKGROUND: Chronic health effects of traffic-related air pollution, like nitrogen dioxide (NO2), are well-documented. Animal models suggested that NO2 exposures dysregulate cortisol function. OBJECTIVES: We evaluated the association between traffic-related NO2 exposure and adolescent human cortisol concentrations, utilizing measures of the cortisol diurnal slope. METHODS: 140 adolescents provided repeated salivary cortisol samples throughout one day. We built a land use regression model to estimate chronic NO2 exposures based on home and school addresses. We then generated model-based estimates of the association between cortisol and NO2 exposure one year prior to cortisol sampling, examining changes in cortisol diurnal slope. The final model was adjusted other criteria pollutants, measures of psychosocial stress, anthropometry, and other demographic and covariates. RESULTS: We observed a decrease in diurnal slope in cortisol for adolescents exposed to the estimated 75th percentile of ambient NO2 (high exposure) relative to those exposed at the 25th percentile (low exposure). For a highly exposed adolescent, the log cortisol was lower by 0.06 µg/dl at waking (95% CI: -0.15, 0.02), 0.07 µg/dl at 30 min post waking (95% CI: -0.15, 0.02), and higher by 0.05 µg/dl at bedtime (95% CI: 0.05, 0.15), compared to a low exposed adolescent. For an additional interquartile range of exposure, the model-based predicted diurnal slope significantly decreased by 0.12 (95% CI: -0.23, -0.01). CONCLUSIONS: In adolescents, we found that increased, chronic exposure to NO2 and the mixture of pollutants from traffic sources was associated with a flattened diurnal slope of cortisol, a marker of an abnormal cortisol response which we hypothesize may be a mechanism through which air pollution may affect respiratory function and asthma in adolescents.

Robust functional clustering of ERP data with application to a study of implicit learning in autism.

Motivated by a study on visual implicit learning in young children with Autism Spectrum Disorder (ASD), we propose a robust functional clustering (RFC) algorithm to identify subgroups within electroencephalography (EEG) data. The proposed RFC is an iterative algorithm based on functional principal component analysis, where cluster membership is updated via predictions of the functional trajectories obtained through a non-parametric random effects model. We consider functional data resulting from event-related potential (ERP) waveforms representing EEG time-locked to stimuli over the course of an implicit learning experiment, after applying a previously proposed meta-preprocessing step. This meta-preprocessing is designed to increase the low signal-to-noise ratio in the raw data and to mitigate the longitudinal changes in the ERP waveforms which characterize the nature and speed of learning. The resulting functional ERP components (peak amplitudes and latencies) inherently exhibit covariance heterogeneity due to low data quality over some stimuli inducing the averaging of different numbers of waveforms in sliding windows of the meta-preprocessing step. The proposed RFC algorithm incorporates this known covariance heterogeneity into the clustering algorithm, improving cluster quality, as illustrated in the data application and extensive simulation studies. ASD is a heterogeneous syndrome and identifying subgroups within ASD children is of interest for understanding the diverse nature of this complex disorder. Applications to the implicit learning paradigm identify subgroups within ASD and typically developing children with diverse learning patterns over the course of the experiment, which may inform clinical stratification of ASD.

The Genetic Heterogeneity among Different Mouse Strains Impacts the Lung Injury Potential of Multiwalled Carbon Nanotubes.

Genetic variation constitutes an important variable impacting the susceptibility to inhalable toxic substances and air pollutants, as reflected by epidemiological studies in humans and differences among animal strains. While multiwalled carbon nanotubes (MWCNTs) are capable of causing lung fibrosis in rodents, it is unclear to what extent the genetic variation in different mouse strains influence the outcome. Four inbred mouse strains, including C57Bl/6, Balb/c, NOD/ShiLtJ, and A/J, to test the pro-fibrogenic effects of a library of MWCNTs in vitro and in vivo are chosen. Ex vivo analysis of IL-1ß production in bone marrow-derived macrophages (BMDMs) as molecular initiating event (MIE) is performed. The order of cytokine production (Balb/c > A/J > C57Bl/6 > NOD/ShiLtJ) in BMDMs is also duplicated during assessment of IL-1ß production in the bronchoalveolar lavage fluid of the same mouse strains 40 h after oropharyngeal instillation of a representative MWCNT. Animal test after 21 d also confirms a similar hierarchy in TGF-ß1 production and collagen deposition in the lung. Statistical analysis confirms a correlation between IL-1ß production in BMDM and the lung fibrosis. All considered, these data demonstrate that genetic background indeed plays a major role in determining the pro-fibrogenic response to MWCNTs in the lung.

Incidence, survival, pathology, and genetics of adult Latino Americans with glioblastoma.

Latino Americans are a rapidly growing ethnic group in the United States but studies of glioblastoma in this population are limited. We have evaluated characteristics of 21,184 glioblastoma patients from the Surveillance, Epidemiology, and End Results (SEER) Program of the National Cancer Institute. This SEER data from 2001 to 2011 draws from 28% of the U.S. POPULATION: Latinos have a lower incidence of GBM and present slightly younger than non-Latino Whites. Cubans present at an older age than other Latino sub-populations. Latinos have a higher incidence of giant cell glioblastoma than non-Latino Whites while the incidence of gliosarcoma is similar. Despite lower rates of radiation therapy and greater rates of sub-total resection than non-Latino Whites, Latinos have better 1 and 5 year survival rates. SEER does not record chemotherapy data. Survivals of Latino sub-populations are similar with each other. Age, extent of resection, and the use of radiation therapy are associated with improved survival but none of these variables are sufficient in a multivariate analysis to explain the improved survival of Latinos relative to non-Latino Whites. As molecular data is not available in SEER records, we studied the MGMT and IDH status of 571 patients from a UCLA database. MGMT methylation and IDH1 mutation rates are not statistically significantly different between non-Latino Whites and Latinos. For UCLA patients with available information, chemotherapy and radiation rates are similar for non-Latino White and Latino patients, but the latter have lower rates of gross total resection and present at a younger age.

A multi-dimensional functional principal components analysis of EEG data.

The electroencephalography (EEG) data created in event-related potential (ERP) experiments have a complex high-dimensional structure. Each stimulus presentation, or trial, generates an ERP waveform which is an instance of functional data. The experiments are made up of sequences of multiple trials, resulting in longitudinal functional data and moreover, responses are recorded at multiple electrodes on the scalp, adding an electrode dimension. Traditional EEG analyses involve multiple simplifications of this structure to increase the signal-to-noise ratio, effectively collapsing the functional and longitudinal components by identifying key features of the ERPs and averaging them across trials. Motivated by an implicit learning paradigm used in autism research in which the functional, longitudinal, and electrode components all have critical interpretations, we propose a multidimensional functional principal components analysis (MD-FPCA) technique which does not collapse any of the dimensions of the ERP data. The proposed decomposition is based on separation of the total variation into subject and subunit level variation which are further decomposed in a two-stage functional principal components analysis. The proposed methodology is shown to be useful for modeling longitudinal trends in the ERP functions, leading to novel insights into the learning patterns of children with Autism Spectrum Disorder (ASD) and their typically developing peers as well as comparisons between the two groups. Finite sample properties of MD-FPCA are further studied via extensive simulations.

Detecting differential patterns of interaction in molecular pathways.

We consider statistical inference for potentially heterogeneous patterns of association characterizing the expression of bio-molecular pathways across different biologic conditions. We discuss a modeling approach based on Gaussian-directed acyclic graphs and provide computational and methodological details needed for posterior inference. Our application finds motivation in reverse phase protein array data from a study on acute myeloid leukemia, where interest centers on contrasting refractory versus relapsed patients. We illustrate the proposed method through both synthetic and case study data.

Bayesian inference for latent biologic structure with determinantal point processes (DPP).

We discuss the use of the determinantal point process (DPP) as a prior for latent structure in biomedical applications, where inference often centers on the interpretation of latent features as biologically or clinically meaningful structure. Typical examples include mixture models, when the terms of the mixture are meant to represent clinically meaningful subpopulations (of patients, genes, etc.). Another class of examples are feature allocation models. We propose the DPP prior as a repulsive prior on latent mixture components in the first example, and as prior on feature-specific parameters in the second case. We argue that the DPP is in general an attractive prior model for latent structure when biologically relevant interpretation of such structure is desired. We illustrate the advantages of DPP prior in three case studies, including inference in mixture models for magnetic resonance images (MRI) and for protein expression, and a feature allocation model for gene expression using data from The Cancer Genome Atlas. An important part of our argument are efficient and straightforward posterior simulation methods. We implement a variation of reversible jump Markov chain Monte Carlo simulation for inference under the DPP prior, using a density with respect to the unit rate Poisson process.

Identifying longitudinal trends within EEG experiments.

Differential brain response to sensory stimuli is very small (a few microvolts) compared to the overall magnitude of spontaneous electroencephalogram (EEG), yielding a low signal-to-noise ratio (SNR) in studies of event-related potentials (ERP). To cope with this phenomenon, stimuli are applied repeatedly and the ERP signals arising from the individual trials are averaged at the subject level. This results in loss of information about potentially important changes in the magnitude and form of ERP signals over the course of the experiment. In this article, we develop a meta-preprocessing step utilizing a moving average of ERP across sliding trial windows, to capture such longitudinal trends. We embed this procedure in a weighted linear mixed effects model to describe longitudinal trends in features such as ERP peak amplitude and latency across trials while adjusting for the inherent heteroskedasticity created at the meta-preprocessing step. The proposed unified framework, including the meta-processing and the weighted linear mixed effects modeling steps, is referred to as MAP-ERP (moving-averaged-processed ERP). We perform simulation studies to assess the performance of MAP-ERP in reconstructing existing longitudinal trends and apply MAP-ERP to data from young children with autism spectrum disorder (ASD) and their typically developing counterparts to examine differences in patterns of implicit learning, providing novel insights about the mechanisms underlying social and/or cognitive deficits in this disorder.

Toxicity of metal oxide nanoparticles in Escherichia coli correlates with conduction band and hydration energies.

Metal oxide nanoparticles (MOx NPs) are used for a host of applications, such as electronics, cosmetics, construction, and medicine, and as a result, the safety of these materials to humans and the environment is of considerable interest. A prior study of 24 MOx NPs in mammalian cells revealed that some of these materials show hazard potential. Here, we report the growth inhibitory effects of the same series of MOx NPs in the bacterium Escherichia coli and show that toxicity trends observed in E. coli parallel those seen previously in mammalian cells. Of the 24 materials studied, only ZnO, CuO, CoO, Mn2O3, Co3O4, Ni2O3, and Cr2O3 were found to exert significant growth inhibitory effects; these effects were found to relate to membrane damage and oxidative stress responses in minimal trophic media. A correlation of the toxicological data with physicochemical parameters of MOx NPs revealed that the probability of a MOx NP being toxic increases as the hydration enthalpy becomes less negative and as the conduction band energy approaches those of biological molecules. These observations are consistent with prior results observed in mammalian cells, revealing that mechanisms of toxicity of MOx NPs are consistent across two very different taxa. These results suggest that studying nanotoxicity in E. coli may help to predict toxicity patterns in higher organisms.

Small pretreatment lesion size and high sphericity as favorable prognostic factors after laser interstitial thermal therapy in brain metastases.

OBJECTIVE: The objective of this study was to identify baseline clinical and radiological characteristics of brain metastases (BMs) associated with a higher probability of lesion-specific progression-free survival (PFS-L) after laser interstitial thermal therapy (LITT). METHODS: A total of 47 lesions in 42 patients with BMs treated with LITT were retrospectively examined, including newly diagnosed BM, suspected recurrent BM, and suspected radiation necrosis. The association of baseline clinical and radiological features with PFS-L was assessed using survival analyses. Radiological features included lesion size measurements, diffusion and perfusion metrics, and sphericity, which is a radiomic feature ranging from 1 (perfect sphere) to 0. RESULTS: The probability of PFS-L for the entire cohort was 88.0% at 3 months, 70.6% at 6 months, 67.4% at 1 and 2 years, and 62.2% at 3 years. For lesions progressing after LITT (n = 13), the median time to progression was 3.9 months, and most lesions (n = 11) progressed within 6 months after LITT. In lesions showing response to LITT (n = 17), the median time to response was 12.1 months. All 3 newly diagnosed BMs showed a long-term response. The mean (± SD) follow-up duration for all censored lesions (n = 34) was 20.7 ± 19.4 months (range 12 days to 6.1 years). The mean pretreatment enhancing volume was 2.68 cm3 and the mean sphericity was 0.70. Pretreatment small enhancing volume (p = 0.003) and high sphericity (p = 0.024) computed from lesion segmentation predicted a longer PFS-L after LITT. Lesions meeting optimal cutoffs of either enhancing volume < 2.5 cm3 (adjusted p = 0.004) or sphericity ≥ 0.705 (adjusted p = 0.019) had longer PFS-L, and their probability of PFS-L was 86.8% at 3 years. Lesions meeting both cutoffs showed a cumulative benefit (p < 0.0001), with a 100% probability of PFS-L at 3 years, which was unchanged at the end of follow-up (4.1 years). Manually computed estimates of lesion size (maximal axial diameter, p = 0.011) and sphericity (p = 0.043) were also predictors of PFS-L. Optimal cutoffs of diameter < 2 cm (adjusted p = 0.035) or manual sphericity ≥ 0.91 (adjusted p = 0.092) identified lesions with longer PFS-L, and lesions meeting both cutoffs showed a cumulative benefit (p = 0.0023). Baseline diffusion imaging did not predict PFS-L. A subset of lesions (n = 7) with highly perfused hotspots had worse PFS-L (adjusted p = 0.010), but perfusion signal contamination from vessels and cortex and underlying size differences were possible confounders. CONCLUSIONS: Small size and high sphericity are ideal baseline features for lesions considered for LITT treatment, with a cumulative PFS-L benefit when both features are present, that could aid patient selection.

Quantification of T2-FLAIR Mismatch in Nonenhancing Diffuse Gliomas Using Digital Subtraction.

BACKGROUND AND PURPOSE: The T2-FLAIR mismatch sign on MR imaging is a highly specific imaging biomarker of isocitrate dehydrogenase (IDH)-mutant astrocytomas, which lack 1p/19q codeletion. However, most studies using the T2-FLAIR mismatch sign have used visual assessment. This study quantified the degree of T2-FLAIR mismatch using digital subtraction of fluid-nulled T2-weighted FLAIR images from non-fluid-nulled T2-weighted images in human nonenhancing diffuse gliomas and then used this information to assess improvements in diagnostic performance and investigate subregion characteristics within these lesions. MATERIALS AND METHODS: Two cohorts of treatment-naïve, nonenhancing gliomas with known IDH and 1p/19q status were studied (n = 71 from The Cancer Imaging Archive (TCIA) and n = 34 in the institutional cohort). 3D volumes of interest corresponding to the tumor were segmented, and digital subtraction maps of T2-weighted MR imaging minus T2-weighted FLAIR MR imaging were used to partition each volume of interest into a T2-FLAIR mismatched subregion (T2-FLAIR mismatch, corresponding to voxels with positive values on the subtraction maps) and nonmismatched subregion (T2-FLAIR nonmismatch corresponding to voxels with negative values on the subtraction maps). Tumor subregion volumes, percentage of T2-FLAIR mismatch volume, and T2-FLAIR nonmismatch subregion thickness were calculated, and 2 radiologists assessed the T2-FLAIR mismatch sign with and without the aid of T2-FLAIR subtraction maps. RESULTS: Thresholds of ≥42% T2-FLAIR mismatch volume classified IDH-mutant astrocytoma with a specificity/sensitivity of 100%/19.6% (TCIA) and 100%/31.6% (institutional); ≥25% T2-FLAIR mismatch volume showed 92.0%/32.6% and 100%/63.2% specificity/sensitivity, and ≥15% T2-FLAIR mismatch volume showed 88.0%/39.1% and 93.3%/79.0% specificity/sensitivity. In IDH-mutant astrocytomas with ≥15% T2-FLAIR mismatch volume, T2-FLAIR nonmismatch subregion thickness was negatively correlated with the percentage T2-FLAIR mismatch volume (P < .0001) across both cohorts. The percentage T2-FLAIR mismatch volume was higher in grades 3-4 compared with grade 2 IDH-mutant astrocytomas (P < .05), and ≥15% T2-FLAIR mismatch volume IDH-mutant astrocytomas were significantly larger than <15% T2-FLAIR mismatch volume IDH-mutant astrocytoma (P < .05) across both cohorts. When evaluated by 2 radiologists, the additional use of T2-FLAIR subtraction maps did not show a significant difference in interreader agreement, sensitivity, or specificity compared with a separate evaluation of T2-FLAIR and T2-weighted MR imaging alone. CONCLUSIONS: T2-FLAIR digital subtraction maps may be a useful, automated tool to obtain objective segmentations of tumor subregions based on quantitative thresholds for classifying IDH-mutant astrocytomas using the percentage T2 FLAIR mismatch volume with 100% specificity and exploring T2-FLAIR mismatch/T2-FLAIR nonmismatch subregion characteristics. Conversely, the addition of T2-FLAIR subtraction maps did not enhance the sensitivity or specificity of the visual T2-FLAIR mismatch sign assessment by experienced radiologists.