"Synthetic" DSC perfusion MRI with adjustable acquisition parameters in brain tumors using dynamic spin-and-gradient-echo echoplanar imaging.

BACKGROUND AND PURPOSE: Normalized relative cerebral blood volume (nrCBV) and percentage of signal recovery (PSR) computed from dynamic susceptibility contrast (DSC) perfusion imaging are useful biomarkers for differential diagnosis and treatment response assessment in brain tumors. However, their measurements are dependent on DSC acquisition factors, and CBV-optimized protocols technically differ from PSR-optimized protocols. This study aimed to generate "synthetic" DSC data with adjustable synthetic acquisition parameters using dual-echo gradient-echo (GE) DSC datasets extracted from dynamic spin-and-gradient-echo echoplanar imaging (dynamic SAGE-EPI). Synthetic DSC was aimed at: 1) simultaneously create nrCBV and PSR maps using optimal sequence parameters, 2) compare DSC datasets with heterogeneous external cohorts, and 3) assess the impact of acquisition factors on DSC metrics. MATERIALS AND METHODS: Thirty-eight patients with contrast-enhancing brain tumors were prospectively imaged with dynamic SAGE-EPI during a non-preloaded single-dose contrast injection and included in this cross-sectional study. Multiple synthetic DSC curves with desired pulse sequence parameters were generated using the Bloch equations applied to the dual-echo GE data extracted from dynamic SAGE-EPI datasets, with or without optional preload simulation. RESULTS: Dynamic SAGE-EPI allowed for simultaneous generation of CBV-optimized and PSR-optimized DSC datasets with a single contrast injection, while PSR computation from guideline-compliant CBV-optimized protocols resulted in rank variations within the cohort (Spearman's ρ=0.83-0.89, i.e. 31%-21% rank variation). Treatment-naïve glioblastoma exhibited lower parameter-matched PSR compared to the external cohorts of treatment-naïve primary CNS lymphomas (PCNSL) (p<0.0001), supporting a role of synthetic DSC for multicenter comparisons. Acquisition factors highly impacted PSR, and nrCBV without leakage correction also showed parameter-dependence, although less pronounced. However, this dependence was remarkably mitigated by post-hoc leakage correction. CONCLUSIONS: Dynamic SAGE-EPI allows for simultaneous generation of CBV-optimized and PSR-optimized DSC data with one acquisition and a single contrast injection, facilitating the use of a single perfusion protocol for all DSC applications. This approach may also be useful for comparisons of perfusion metrics across heterogeneous multicenter datasets, as it facilitates post-hoc harmonization. ABBREVIATIONS: DSC = dynamic susceptibility contrast; FA = flip angle; GBCA = gadolinium-based contrast agent; GBM = glioblastoma; GE = gradient echo; IDH = isocitrate dehydrogenase; IDHm = IDH-mutant; IDHwt = IDH-wild-type; 1p19qcod = 1p19q codeleted; 1p19qint = 1p19q intact; MRI = magnetic resonance imaging; PCNSL = primary CNS lymphoma; PSR = percentage of signal recovery; Rec = recurrent; SAGE-EPI = spin-and-gradient-echo echoplanar imaging; CBV = cerebral blood volume; nrCBV = normalized relative CBV; ROI = region of interest; TE = echo time; TN = treatment-naïve; TR = repetition time.

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.

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.

Seizure outcomes following single-fraction versus hypofractionated radiosurgery for brain metastases: a single-center experience.

OBJECTIVE: Although seizures are a relatively common phenomenon in the setting of brain metastases (BMs), there are no discrete recommendations regarding the use of antiepileptic drugs (AEDs) in this population, either in general or in the context of treatment. The authors' aim was to better understand the underlying pathological factors as well as the therapeutic techniques that may lead to seizures following the radiosurgical treatment of BMs with the goal of guiding appropriate AED prophylaxis. METHODS: Adult patients with BMs diagnosed from 2013 to 2020 at a single academic institution and treated with radiation therapy were included in this study. The authors evaluated factors associated with the incidence of seizures throughout the disease course, with a focus on seizures in the 90-day period following stereotactic radiosurgery (SRS). RESULTS: Four hundred forty-four patients with newly diagnosed BMs were identified, 10% of whom had seizures at the time of presentation and 28% of whom had a seizure at any point during the study period. Tumor histology was significantly associated with initial seizure risk. AED use was highly variable. In the 90-day post-SRS period, the summed total planning target volume (PTV) was independently predictive of post-SRS seizures, regardless of the fractionation scheme (single fraction vs hypofractionated) and other clinical factors. The number of supratentorial BMs was not predictive of post-SRS seizures. CONCLUSIONS: PTV is a superior predictor of post-SRS seizures relative to the number of supratentorial BMs, as it serves as a volumetric proxy for intracranial disease burden. A larger PTV, alongside tumor histology and prior seizure history, should be considered in the decision-making process for AED use following radiosurgery.

Change in volumetric tumor growth rate after cytotoxic therapy is predictive of overall survival in recurrent glioblastoma.

Background: Alterations in tumor growth rate (TGR) in recurrent glioblastoma (rGBM) after treatment may be useful for identifying therapeutic activity. The aim of this study was to assess the impact of volumetric TGR alterations on overall survival (OS) in rGBM treated with chemotherapy with or without radiation therapy (RT). Methods: Sixty-one rGBM patients treated with chemotherapy with or without concomitant radiation therapy (RT) at 1st or 2nd recurrence were retrospectively examined. Pre- and post-treatment contrast enhancing volumes were computed. Patients were considered "responders" if they reached progression-free survival at 6 months (PFS6) and showed a decrease in TGR after treatment and "non-responders" if they didn't reach PFS6 or if TGR increased. Results: Stratification by PFS6 and based on TGR resulted in significant differences in OS both for all patients and for patients without RT (P < 0.05). A decrease of TGR (P = 0.009), smaller baseline tumor volume (P = 0.02), O6-methylguanine-DNA methyltransferase promoter methylation (P = 0.048) and fewer number of recurrences (P = 0.048) were significantly associated with longer OS after controlling for age, sex and concomitant RT. Conclusion: A decrease in TGR in patients with PFS6, along with smaller baseline tumor volume, were associated with a significantly longer OS in rGBM treated with chemotherapy with or without radiation. Importantly, all patients that exhibited PFS6 also showed a measurable decrease in TGR.

Depth of Radiographic Response and Time to Tumor Regrowth Predicts Overall Survival Following Anti-VEGF Therapy in Recurrent Glioblastoma.

PURPOSE: Antiangiogenic therapies are known to cause high radiographic response rates due to reduction in vascular permeability resulting in a lower degree of contrast extravasation. In this study, we investigate the prognostic ability for model-derived parameters describing enhancing tumor volumetric dynamics to predict survival in recurrent glioblastoma treated with antiangiogenic therapy. EXPERIMENTAL DESIGN: N = 276 patients in two phase II trials were used as training data, including bevacizumab ± irinotecan (NCT00345163) and cabozantinib (NCT00704288), and N = 74 patients in the bevacizumab arm of a phase III trial (NCT02511405) were used for validation. Enhancing volumes were estimated using T1 subtraction maps, and a biexponential model was used to estimate regrowth (g) and regression (d) rates, time to tumor regrowth (TTG), and the depth of response (DpR). Response characteristics were compared to diffusion MR phenotypes previously shown to predict survival. RESULTS: Optimized thresholds occurred at g = 0.07 months-1 (phase II: HR = 0.2579, P = 5 × 10-20; phase III: HR = 0.2197, P = 5 × 10-5); d = 0.11 months-1 (HR = 0.3365, P < 0.0001; HR = 0.3675, P = 0.0113); TTG = 3.8 months (HR = 0.2702, P = 6 × 10-17; HR = 0.2061, P = 2 × 10-5); and DpR = 11.3% (HR = 0.6326, P = 0.0028; HR = 0.4785, P = 0.0206). Multivariable Cox regression controlling for age and baseline tumor volume confirmed these factors as significant predictors of survival. Patients with a favorable pretreatment diffusion MRI phenotype had a significantly longer TTG and slower regrowth. CONCLUSIONS: Recurrent glioblastoma patients with a large, durable radiographic response to antiangiogenic agents have significantly longer survival. This information is useful for interpreting activity of antiangiogenic agents in recurrent glioblastoma.

Application of a genetic signature of late GU toxicity in SCIMITAR, a Post-op SBRT trial.

Predictors of genitourinary toxicity after post-prostatectomy radiotherapy remain elusive. A previously defined germline DNA signature (PROSTOX) has shown predictive ability for late grade ≥ 2 GU toxicity after intact prostate stereotactic body radiotherapy. We explore whether PROSTOX would predict toxicity among patients receiving post-prostatectomy SBRT on a phase II clinical trial.

Germline biomarkers predict toxicity to anti-PD1/PDL1 checkpoint therapy.

BACKGROUND: There is great interest in finding ways to identify patients who will develop toxicity to cancer therapies. This has become especially pressing in the era of immune therapy, where toxicity can be long-lasting and life-altering, and primarily comes in the form of immune-related adverse effects (irAEs). Treatment with the first drugs in this class, anti-programmed death 1 (anti-PD1)/programmed death-ligand 1 (PDL1) checkpoint therapies, results in grade 2 or higher irAEs in up to 25%-30% of patients, which occur most commonly within the first 6 months of treatment and can include arthralgias, rash, pruritus, pneumonitis, diarrhea and/or colitis, hepatitis, and endocrinopathies. We tested the hypothesis that germline microRNA pathway functional variants, known to predict altered systemic stress responses to cancer therapies, would predict irAEs in patients across cancer types. METHODS: MicroRNA pathway variants were evaluated for an association with grade 2 or higher toxicity using four classifiers on 62 patients with melanoma, and then the panel's performance was validated on 99 patients with other cancer types. Trained classifiers included classification trees, LASSO-regularized logistic regression, boosted trees, and random forests. Final performance measures were reported on the training set using leave-one-out cross validation and validated on held-out samples. The predicted probability of toxicity was evaluated for its association, if any, with response categories to anti-PD1/PDL1 therapy in the melanoma cohort. RESULTS: A biomarker panel was identified that predicts toxicity with 80% accuracy (F1=0.76, area under the curve (AUC)=0.82) in the melanoma training cohort and 77.6% accuracy (F1=0.621, AUC=0.778) in the pan-cancer validation cohort. In the melanoma cohort, the predictive probability of toxicity was not associated with response categories to anti-PD1/PDL1 therapy (p=0.70). In the same cohort, the most significant biomarker of toxicity in RAC1, predicting a greater than ninefold increased risk of toxicity (p<0.001), was also not associated with response to anti-PD1/PDL1 therapy (p=0.151). CONCLUSIONS: A germline microRNA-based biomarker signature predicts grade 2 and higher irAEs to anti-PD1/PDL1 therapy, regardless of tumor type, in a pan-cancer manner. These findings represent an important step toward personalizing checkpoint therapy, the use of which is growing rapidly.

Germline variants disrupting microRNAs predict long-term genitourinary toxicity after prostate cancer radiation.

BACKGROUND AND PURPOSE: The purpose of this study was to determine whether single nucleotide polymorphisms disrupting microRNA targets (mirSNPs) can serve as predictive biomarkers for toxicity after radiotherapy for prostate cancer and whether these may be differentially predictive depending on radiation fractionation. MATERIALS AND METHODS: We identified 201 men treated with two forms of definitive radiotherapy for prostate cancer at two institutions: 108 men received conventionally-fractionated radiotherapy (CF-RT) and 93 received stereotactic body radiotherapy (SBRT). Germline DNA was evaluated for the presence of functional mirSNPs. Random forest, boosted trees and elastic net models were developed to predict late grade ≥2 GU toxicity by the RTOG scale. RESULTS: The crude incidence of late grade ≥2 GU toxicity was 16% after CF-RT and 15% after SBRT. An elastic net model based on 22 mirSNPs differentiated CF-RT patients at high risk (71.5%) versus low risk (7.5%) for toxicity, with an area under the curve (AUC) values of 0.76-0.81. An elastic net model based on 32 mirSNPs differentiated SBRT patients at high risk (64.7%) versus low risk (3.9%) for toxicity, with an area under the curve (AUC) values of 0.81-0.87. These models were specific to treatment type delivered. Prospective studies are warranted to further validate these results. CONCLUSION: Predictive models using germline mirSNPs have high accuracy for predicting late grade ≥2 GU toxicity after either CF-RT or SBRT, and are unique for each treatment, suggesting that germline predictors of late radiation sensitivity are fractionation-dependent. Prospective studies are warranted to further validate these results.

Liposomal Delivery of Mitoxantrone and a Cholesteryl Indoximod Prodrug Provides Effective Chemo-immunotherapy in Multiple Solid Tumors.

We developed a custom-designed liposome carrier for codelivery of a potent immunogenic cell death (ICD) stimulus plus an inhibitor of the indoleamine 2,3-dioxygenase (IDO-1) pathway to establish a chemo-immunotherapy approach for solid tumors in syngeneic mice. The carrier was constructed by remote import of the anthraquinone chemotherapeutic agent, mitoxantrone (MTO), into the liposomes, which were further endowed with a cholesterol-conjugated indoximod (IND) prodrug in the lipid bilayer. For proof-of-principle testing, we used IV injection of the MTO/IND liposome in a CT26 colon cancer model to demonstrate the generation of a robust immune response, characterized by the appearance of ICD markers (CRT and HMGB-1) as well as evidence of cytotoxic cancer cell death, mediated by perforin and granzyme B. Noteworthy, the cytotoxic effects involved natural killer (NK) cell, which suggests a different type of ICD response. The immunotherapy response was significantly augmented by codelivery of the IND prodrug, which induced additional CRT expression, reduced number of Foxp3+ Treg, and increased perforin release, in addition to extending animal survival beyond the effect of an MTO-only liposome. The outcome reflects the improved pharmacokinetics of MTO delivery to the cancer site by the carrier. In light of the success in the CT26 model, we also assessed the platform efficacy in further breast cancer (EMT6 and 4T1) and renal cancer (RENCA) models, which overexpress IDO-1. Encapsulated MTO delivery was highly effective for inducing chemo-immunotherapy responses, with NK participation, in all tumor models. Moreover, the growth inhibitory effect of MTO was enhanced by IND codelivery in EMT6 and 4T1 tumors. All considered, our data support the use of encapsulated MTO delivery for chemo-immunotherapy, with the possibility to boost the immune response by codelivery of an IDO-1 pathway inhibitor.

A Phase II Trial of 5-Day Neoadjuvant Radiotherapy for Patients with High-Risk Primary Soft Tissue Sarcoma.

PURPOSE: In a single-institution phase II study, we evaluated the safety of a 5-day dose-equivalent neoadjuvant radiotherapy (RT) regimen for high-risk primary soft tissue sarcoma. PATIENTS AND METHODS: Patients received neoadjuvant RT alone (30 Gy in five fractions) to the primary tumor with standard margins. The primary endpoint was grade ≥2 late-radiation toxicity. Major wound complications, local recurrences, and distant metastases were also examined. In exploratory analysis, we evaluated germline biomarkers for wound toxicity and the effects of the study on treatment utilization. RESULTS: Over 2 years, 52 patients were enrolled with median follow-up of 29 months. Seven of 44 evaluable patients (16%) developed grade ≥2 late toxicity. Major wound complications occurred in 16 of 50 patients (32%); a signature defined by 19 germline SNPs in miRNA-binding sites of immune and DNA damage response genes, in addition to lower extremity tumor location, demonstrated strong predictive performance for major wound complications. Compared with the preceding 2-year period, the number of patients treated with neoadjuvant RT alone at our institution increased 3-fold, with a concomitant increase in the catchment area. CONCLUSIONS: A shorter 5-day neoadjuvant RT regimen results in favorable rates of wound complications and grade ≥2 toxicity after 2-year follow-up. Five-day RT significantly increased utilization of neoadjuvant RT at our high-volume sarcoma center. With further validation, a putative germline biomarker for wound complications may guide safer RT utilization.

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.

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.

Ensemble of sparse classifiers for high-dimensional biological data.

Biological data are often high in dimension while the number of samples is small. In such cases, the performance of classification can be improved by reducing the dimension of data, which is referred to as feature selection. Recently, a novel feature selection method has been proposed utilising the sparsity of high-dimensional biological data where a small subset of features accounts for most variance of the dataset. In this study we propose a new classification method for high-dimensional biological data, which performs both feature selection and classification within a single framework. Our proposed method utilises a sparse linear solution technique and the bootstrap aggregating algorithm. We tested its performance on four public mass spectrometry cancer datasets along with two other conventional classification techniques such as Support Vector Machines and Adaptive Boosting. The results demonstrate that our proposed method performs more accurate classification across various cancer datasets than those conventional classification techniques.

Ribosomal Proteins RPS11 and RPS20, Two Stress-Response Markers of Glioblastoma Stem Cells, Are Novel Predictors of Poor Prognosis in Glioblastoma Patients.

Glioblastoma stem cells (GSC) co-exhibiting a tumor-initiating capacity and a radio-chemoresistant phenotype, are a compelling cell model for explaining tumor recurrence. We have previously characterized patient-derived, treatment-resistant GSC clones (TRGC) that survived radiochemotherapy. Compared to glucose-dependent, treatment-sensitive GSC clones (TSGC), TRGC exhibited reduced glucose dependence that favor the fatty acid oxidation pathway as their energy source. Using comparative genome-wide transcriptome analysis, a series of defense signatures associated with TRGC survival were identified and verified by siRNA-based gene knockdown experiments that led to loss of cell integrity. In this study, we investigate the prognostic value of defense signatures in glioblastoma (GBM) patients using gene expression analysis with Probeset Analyzer (131 GBM) and The Cancer Genome Atlas (TCGA) data, and protein expression with a tissue microarray (50 GBM), yielding the first TRGC-derived prognostic biomarkers for GBM patients. Ribosomal protein S11 (RPS11), RPS20, individually and together, consistently predicted poor survival of newly diagnosed primary GBM tumors when overexpressed at the RNA or protein level [RPS11: Hazard Ratio (HR) = 11.5, p<0.001; RPS20: HR = 4.5, p = 0.03; RPS11+RPS20: HR = 17.99, p = 0.001]. The prognostic significance of RPS11 and RPS20 was further supported by whole tissue section RPS11 immunostaining (27 GBM; HR = 4.05, p = 0.01) and TCGA gene expression data (578 primary GBM; RPS11: HR = 1.19, p = 0.06; RPS20: HR = 1.25, p = 0.02; RPS11+RPS20: HR = 1.43, p = 0.01). Moreover, tumors that exhibited unmethylated O-6-methylguanine-DNA methyltransferase (MGMT) or wild-type isocitrate dehydrogenase 1 (IDH1) were associated with higher RPS11 expression levels [corr (IDH1, RPS11) = 0.64, p = 0.03); [corr (MGMT, RPS11) = 0.52, p = 0.04]. These data indicate that increased expression of RPS11 and RPS20 predicts shorter patient survival. The study also suggests that TRGC are clinically relevant cells that represent resistant tumorigenic clones from patient tumors and that their properties, at least in part, are reflected in poor-prognosis GBM. The screening of TRGC signatures may represent a novel alternative strategy for identifying new prognostic biomarkers.

Estrogen withdrawal, increased breast cancer risk and the KRAS-variant.

The KRAS-variant is a biologically functional, microRNA binding site variant, which predicts increased cancer risk especially for women. Because external exposures, such as chemotherapy, differentially impact the effect of this mutation, we evaluated the association of estrogen exposures, breast cancer (BC) risk and tumor biology in women with the KRAS-variant. Women with BC (n = 1712), the subset with the KRAS-variant (n = 286) and KRAS-variant unaffected controls (n = 80) were evaluated, and hormonal exposures, KRAS-variant status, and pathology were compared. The impact of estrogen withdrawal on transformation of isogenic normal breast cell lines with or without the KRAS-variant was studied. Finally, the association and presentation characteristics of the KRAS-variant and multiple primary breast cancer (MPBC) were evaluated. KRAS-variant BC patients were more likely to have ovarian removal pre-BC diagnosis than non-variant BC patients (p = 0.033). In addition, KRAS-variant BC patients also appeared to have a lower estrogen state than KRAS-variant unaffected controls, with a lower BMI (P < 0.001). Finally, hormone replacement therapy (HRT) discontinuation in KRAS-variant patients was associated with a diagnosis of triple negative BC (P < 0.001). Biologically confirming our clinical findings, acute estrogen withdrawal led to oncogenic transformation in KRAS-variant positive isogenic cell lines. Finally, KRAS-variant BC patients had greater than an 11-fold increased risk of presenting with MPBC compared to non-variant patients (45.39% vs 6.78%, OR 11.44 [3.42-37.87], P < 0.001). Thus, estrogen withdrawal and a low estrogen state appear to increase BC risk and to predict aggressive tumor biology in women with the KRAS-variant, who are also significantly more likely to present with multiple primary breast cancer.

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.

The CREB-binding protein inhibitor ICG-001 suppresses pancreatic cancer growth.

Pancreatic ductal adenocarcinoma (PDAC) is a highly lethal cancer due in part to a lack of highly robust cytotoxic or molecular-based therapies. Recent studies investigating ligand-mediated Wnt/ß-catenin signaling have highlighted its importance in pancreatic cancer initiation and progression, as well as its potential as a therapeutic target in PDAC. The small-molecule ICG-001 binds cAMP-responsive element binding (CREB)-binding protein (CBP) to disrupt its interaction with ß-catenin and inhibit CBP function as a coactivator of Wnt/ß-catenin-mediated transcription. Given its ability to inhibit Wnt/ß-catenin-mediated transcription in vitro and in vivo, as well as its efficacy in preclinical models of colorectal cancer and other Wnt-driven diseases, we examined ICG-001 and its potential role as a therapeutic in PDAC. ICG-001 alone significantly inhibited anchorage-dependent and -independent growth of multiple PDAC lines, and augmented in vitro growth inhibition when used in combination with gemcitabine. ICG-001 had only variable modest effects on PDAC apoptosis and instead mediated PDAC growth inhibition primarily through robust induction of G1 cell-cycle arrest. These effects, however, seemed decoupled from its inhibition of Wnt/ß-catenin-mediated transcription. DNA microarrays performed on PDAC cells in the context of ICG-001 treatment revealed ICG-001 altered the expression of several genes with well-established roles in DNA replication and cell-cycle progression, including direct actions on SKP2 and CDKN1A. ICG-001 also significantly prolonged survival in an in vivo orthotopic xenograft model of PDAC, indicating ICG-001 or derived compounds that disrupt CBP activity are potentially useful small-molecule therapeutics for pancreatic cancer.

Lyophilized brain tumor specimens can be used for histologic, nucleic acid, and protein analyses after 1 year of room temperature storage.

Frozen tissue, a gold standard biospecimen, can yield well preserved nucleic acids and proteins after over a decade but is vulnerable to thawing and has substantial fiscal, spatial, and environmental costs. A long-term room temperature biospecimen storage alternative that preserves broad analytical utility can potentially empower tissue-based research. As there is scant data on the analytical utility of lyophilized brain tumor biospecimens, we evaluated lyophilized (freeze-dried) samples stored for 1 year at room temperature. Lyophilized tumor tissue processed into paraffin sections produced good histology. Yields of extracted DNA, RNA, and protein approximated those of frozen tissue. After 1 year, lyophilized samples yielded high molecular weight DNA that permitted copy number variation analysis, IDH 1 mutation detection, and MGMT promoter methylation PCR. A 27 % decrease in RIN scores over the 1 year suggests that RNA degradation was inhibited though incompletely. Nevertheless, RT-PCR studies on lyophilized tissue performed similarly to frozen tissue. In contrast to FFPE tissues where protein bands were absent or shifted to a lower molecular weight, lyophilized samples showed similar protein bands as frozen tissue on SDS-PAGE analysis. Lyophilized tissue performed similarly to frozen tissue for Western blots and enzyme activity assays. Immunohistochemistry of lyophilized tissue that were processed into FFPE blocks often required longer incubation times for staining than standard FFPE samples but generally provided robust antigen detection. This preliminary study suggests that lyophilization has promise for long-term room temperature storage while permitting varied tests; however, further work is required to better stabilize nucleic acids particularly RNA.