Daniel C. Edelman

MYC-driven accumulation of 2-hydroxyglutarate is associated with breast cancer prognosis

Metabolic profiling of cancer cells has recently been established as a promising tool for the development of therapies and identification of cancer biomarkers. Here we characterized the metabolomic profile of human breast tumors and uncovered intrinsic metabolite signatures in these tumors using an untargeted discovery approach and validation of key metabolites. The oncometabolite 2-hydroxyglutarate (2HG) accumulated at high levels in a subset of tumors and human breast cancer cell lines. We discovered an association between increased 2HG levels and MYC pathway activation in breast cancer, and further corroborated this relationship using MYC overexpression and knockdown in human mammary epithelial and breast cancer cells. Further analyses revealed globally increased DNA methylation in 2HG-high tumors and identified a tumor subtype with high tissue 2HG and a distinct DNA methylation pattern that was associated with poor prognosis and occurred with higher frequency in African-American patients. Tumors of this subtype had a stem cell-like transcriptional signature and tended to overexpress glutaminase, suggestive of a functional relationship between glutamine and 2HG metabolism in breast cancer. Accordingly, 13C-labeled glutamine was incorporated into 2HG in cells with aberrant 2HG accumulation, whereas pharmacologic and siRNA-mediated glutaminase inhibition reduced 2HG levels. Our findings implicate 2HG as a candidate breast cancer oncometabolite associated with MYC activation and poor prognosis.

DNA methylation profiling identifies global methylation differences and markers of adrenocortical tumors.

CONTEXT It is not known whether there are any DNA methylation alterations in adrenocortical tumors. OBJECTIVE The objective of the study was to determine the methylation profile of normal adrenal cortex and benign and malignant adrenocortical tumors. METHODS Genome-wide methylation status of CpG regions were determined in normal (n = 19), benign (n = 48), primary malignant (n = 8), and metastatic malignant (n = 12) adrenocortical tissue samples. An integrated analysis of genome-wide methylation and mRNA expression in benign vs. malignant adrenocortical tissue samples was also performed. RESULTS Methylation profiling revealed the following: 1) that methylation patterns were distinctly different and could distinguish normal, benign, primary malignant, and metastatic tissue samples; 2) that malignant samples have global hypomethylation; and 3) that the methylation of CpG regions are different in benign adrenocortical tumors by functional status. Normal compared with benign samples had the least amount of methylation differences, whereas normal compared with primary and metastatic adrenocortical carcinoma samples had the greatest variability in methylation (adjusted P ≤ 0.01). Of 215 down-regulated genes (≥2-fold, adjusted P ≤ 0.05) in malignant primary adrenocortical tumor samples, 52 of these genes were also hypermethylated. CONCLUSIONS Malignant adrenocortical tumors are globally hypomethylated as compared with normal and benign tumors. Methylation profile differences may accurately distinguish between primary benign and malignant adrenocortical tumors. Several differentially methylated sites are associated with genes known to be dysregulated in malignant adrenocortical tumors.

An Integrated Prognostic Classifier for Stage I Lung Adenocarcinoma Based on mRNA, microRNA, and DNA Methylation Biomarkers.

Introduction: Up to 30% stage I lung cancer patients suffer recurrence within 5 years of curative surgery. We sought to improve existing protein-coding gene and microRNA expression prognostic classifiers by incorporating epigenetic biomarkers. Methods: Genome-wide screening of DNA methylation and pyrosequencing analysis of HOXA9 promoter methylation were performed in two independently collected cohorts of stage I lung adenocarcinoma. The prognostic value of HOXA9 promoter methylation alone and in combination with mRNA and miRNA biomarkers was assessed by Cox regression and Kaplan–Meier survival analysis in both cohorts. Results: Promoters of genes marked by polycomb in embryonic stem cells were methylated de novo in tumors and identified patients with poor prognosis. The HOXA9 locus was methylated de novo in stage I tumors (p < 0.0005). High HOXA9 promoter methylation was associated with worse cancer-specific survival (hazard ratio [HR], 2.6; p = 0.02) and recurrence-free survival (HR, 3.0; p = 0.01), and identified high-risk patients in stratified analysis of stages IA and IB. Four protein-coding gene (XPO1, BRCA1, HIF1&agr;, and DLC1), miR-21 expression, and HOXA9 promoter methylation were each independently associated with outcome (HR, 2.8; p = 0.002; HR, 2.3; p = 0.01; and HR, 2.4; p = 0.005, respectively), and when combined, identified high-risk, therapy naive, stage I patients (HR, 10.2; p = 3 × 10−5). All associations were confirmed in two independently collected cohorts. Conclusion: A prognostic classifier comprising three types of genomic and epigenomic data may help guide the postoperative management of stage I lung cancer patients at high risk of recurrence.

Genome-Wide Methylation Patterns in Papillary Thyroid Cancer Are Distinct Based on Histological Subtype and Tumor Genotype

CONTEXT Aberrant DNA methylation is known to be a major factor in oncogenesis and cancer progression, but effects of methylation in papillary thyroid cancer (PTC) are not well defined. OBJECTIVE The objective of the study was to identify altered methylation patterns, which may be associated with PTC disease behavior. DESIGN This study was a genome-wide methylation analysis of PTC. SETTING The study was conducted at the National Institutes of Health Clinical Center. PATIENTS PTC tissue from 51 patients were analyzed and compared with normal thyroid tissue from seven patients. INTERVENTIONS CpG methylation status was assessed using advanced genome-wide methylation bead chips. OUTCOME MEASURES Altered methylation patterns in PTC were analyzed by stage, recurrence, histological subtype of tumor, and tumor genotype. RESULTS PTC is globally hypomethylated compared with normal thyroid with 2837 differentially methylated CpG sites. The follicular variant of PTC demonstrated less differential methylation with only 569 differentially methylated CpG sites. Tumors with mutations in BRAF, RET/PTC, and RAS demonstrated a 3.6-fold increase in the number of differentially methylated sites compared with wild-type tumors. The differentially methylated genes were associated with oncological pathways including cellular movement, growth, and proliferation. CONCLUSION PTC is epigenetically distinct from the follicular variant of PTC and by gene mutation status (BRAF, RET/PTC, and RAS).

Translational predictive biomarker analysis of the phase 1b sorafenib and bevacizumab study expansion cohort

Predictive biomarkers are needed to triage patients to the best therapy. We prospectively planned examination of sequential blood, biopsy, and functional imaging with which to confirm the mechanism and to identify potential predictive biomarkers in a phase Ib clinical trial expansion of patients with solid tumors receiving sorafenib/bevacizumab. The maximally tolerated doses of sorafenib at 200 mg twice daily with bevacizumab at 5 mg/kg every other week were given to biopsiable patients. Patients were randomized to receive either sorafenib or bevacizumab monotherapy for the first 28-day cycle with the second drug added with cycle 2. Biopsies, dynamic contrast-enhanced MRI, and fluorodeoxyglucose-proton emission tomography were done pre-therapy and at 2 and 6 weeks (2 weeks into combination therapy). Tumor and serum proteomics, Ras/Raf mutational analysis, and functional imaging results were examined individually and across the dataset to identify potential changes predictive of response to therapy and those that confirm the biochemical drug mechanism(s). Therapy with sorafenib/bevacizumab resulted in clinical benefit in 45% of this mixed solid tumor group. ERK activation and microvessel density were decreased with monotherapy treatment with sorafenib or bevacizumab, respectively; whereas a decreased signal over the group of total AKT, phospho(p)-VEGF receptor2, p-endothelial nitric-oxide synthase, b-RAF, and cleaved poly(ADP-ribose) polymerase was associated with earlier progression of disease. Tumor metabolic activity decreased in those patients with clinical benefits lasting longer than 4 months, and activity increased with progression of disease. Cleavage of caspase 3 and poly(ADP-ribose) polymerase was increased, and Ki67 expression decreased in patients with prolonged clinical benefits, consistent with decreased proliferation and increased apoptosis. The conglomerate analysis, incorporating pharmacodynamic and tumor biochemistry, demonstrated sorafenib/bevacizumab-targeted vascular activity in the tumor. Results suggest potential biomarkers for which changes, as a group, during early therapeutic exposure may predict clinical benefit.

Distinct methylation profiles characterize fusion-positive and fusion-negative rhabdomyosarcoma.

Rhabdomyosarcoma comprises two major subtypes, fusion positive (PAX3–FOXO1 or PAX7–FOXO1) and fusion negative. To investigate the significance of DNA methylation in these subtypes, we analyzed methylation profiles of 37 rhabdomyosarcoma tumors and 10 rhabdomyosarcoma cell lines, as well as 8 normal tissues. Unsupervised clustering of DNA methylation clearly distinguished the fusion-positive and fusion-negative subsets. The fusion-positive tumors showed substantially lower overall levels of methylation compared with fusion-negative tumors. Comparison with the methylation pattern of normal skeletal muscle and bone marrow indicates that fusion-negative rhabdomyosarcoma is more similar to these normal tissues compared with fusion-positive rhabdomyosarcoma, and suggests that many of the methylation differences between these subtypes arise from ‘aberrant’ hyper- and hypomethylation events in fusion-positive rhabdomyosarcoma. Integrative methylation and gene expression analysis revealed that methylation differences between fusion-positive and fusion-negative tumors could either be positively or negatively associated with mRNA expression. There was no significant difference in the distribution of PAX3–FOXO1-binding sites between genes with and without differential methylation. However, the finding that PAX3–FOXO1-binding sites were enriched among genes that were both differentially methylated and differentially expressed suggests that the fusion protein interacts with DNA methylation to regulate target gene expression. An 11-gene DNA methylation signature, classifying the rhabdomyosarcoma tumors into fusion-positive and fusion-negative subsets, was established and validated by pyrosequencing assays. Notably, EMILIN1 (part of the 11-gene signature) showed higher methylation and lower mRNA expression in fusion-positive compared with fusion-negative tumors, and demonstrated demethylation and re-expression in multiple fusion-positive cell lines after treatment with 5-aza-2′-deoxycytidine. In conclusion, our study demonstrates that fusion-positive and fusion-negative rhabdomyosarcoma tumors possess characteristic methylation profiles that contribute to the expression differences between these fusion subtypes. These findings indicate an important relationship between fusion status and epigenetic changes in rhabdomyosarcoma, present a novel approach for ascertaining fusion status, and may identify new therapeutic targets in rhabdomyosarcoma.

Telomere Length and Survival of Patients with Hepatocellular Carcinoma in the United States.

Background Telomere shortening is an important molecular event in hepatocellular carcinoma (HCC) initiation; however, its role in HCC progression and prognosis is less clear. Our study aimed to examine the association of telomere length with survival of patients with HCC. Methods We measured telomere length in tumor and adjacent non-tumor tissues from 126 persons with HCC in the United States (U.S.) who were followed for mortality outcomes. Relative telomere length (RTL) was measured by a monochrome multiplex quantitative polymerase chain reaction assay. Multivariable Cox proportional hazards modeling was used to calculate hazard ratios (HRs) and 95% CIs for the association between telomere length and all-cause mortality. We also examined associations between telomere length and patient characteristics using multiple linear regression. Results During a mean follow-up of 6.0 years, 79 deaths occurred among 114 individuals for whom survival data were available. The ratio of RTL in tumor relative to non-tumor tissue was greater for individuals with regional or distant stage tumors (0.97) than localized stage tumors (0.77), and for individuals with grade III or IV tumors (0.95) than grade II (0.88) or grade I (0.67) tumors. An RTL ratio ≥1 was not associated with survival (HR 0.92, 95% CI 0.55, 1.55) compared to a ratio <1, after adjusting for age at diagnosis, sex, tumor stage and tumor size. Similarly, RTL in the tumor and non-tumor tissue, respectively, were not associated with survival. Conclusions This U.S. based study found that telomeres may be longer in more aggressive HCCs. There was no evidence, however, that telomere length was associated with survival of patients with HCC. Future investigations are warranted to clarify the role of telomere length in HCC prognosis.

Aneuploidy, TP53 mutation, and amplification of MYC correlate with increased intratumor heterogeneity and poor prognosis of breast cancer patients

The clinical course of breast cancer varies from one patient to another. Currently, the choice of therapy relies on clinical parameters and histological and molecular tumor features. Alas, these markers are informative in only a subset of patients. Therefore, additional predictors of disease outcome would be valuable for treatment stratification. Extensive studies showed that the degree of variation of the nuclear DNA content, i.e., aneuploidy, determines prognosis. Our aim was to further elucidate the molecular basis of aneuploidy. We analyzed five diploid and six aneuploid tumors with more than 20 years of follow‐up. By performing FISH with a multiplexed panel of 10 probes to enumerate copy numbers in individual cells, and by sequencing 563 cancer‐related genes, we analyzed how aneuploidy is linked to intratumor heterogeneity. In our cohort, none of the patients with diploid tumors died of breast cancer during follow‐up in contrast to four of six patients with aneuploid tumors (mean survival 86.4 months). The FISH analysis showed markedly increased genomic instability and intratumor heterogeneity in aneuploid tumors. MYC gain was observed in only 20% of the diploid cancers, while all aneuploid cases showed a gain. The mutation burden was similar in diploid and aneuploid tumors, however, TP53 mutations were not observed in diploid tumors, but in all aneuploid tumors in our collective. We conclude that quantitative measurements of intratumor heterogeneity by multiplex FISH, detection of MYC amplification and TP53 mutation could augment prognostication in breast cancer patients.

Impact of telomere length on survival in classic and variant hairy cell leukemia

Telomeres, which protect the ends of chromosomes, are shortened in several hematologic malignancies, often with adverse prognostic implications, but their effect on prognosis of classic and variant hairy cell leukemia (HCL and HCLv) has not been reported. HCL/HCLv genomic DNA from 46 patients was studied by PCR to determine the ratio of telomere to single copy gene number (T/S). T/S was unrelated to diagnosis of HCL or HCLv (p=0.27), but shorter T/S was associated with unmutated immunoglobulin rearrangements (p=0.033) and age above the median at diagnosis (p=0.017). Low T/S was associated with shorter overall survival from diagnosis (OS), particularly T/S <0.655 (p=0.0064, adjusted p=0.019). Shorter OS was also associated with presence of unmutated (p<0.0001) or IGHV4-34+ (p<0.0001) rearrangements, or increasing age (p=0.0002). Multivariable analysis with Cox modeling showed that short T/S along with either unmutated or IGHV4-34+ rearrangements remained associated with reduced OS (p=0.0071, p=0.0024, respectively) after age adjustment. While T/S is relatively long in HCL and the disease usually indolent with excellent survival, shortened telomeres in HCL/HCLv are associated with decreased survival. Shortened T/S could represent a risk factor needing further investigation/intervention to determine if non-chemotherapy treatment options, in addition to or instead of chemotherapy, might be particularly useful.

Long-term Stability of Total RNA in RNAstable® as Evaluated by Expression Microarray

Storage of labile RNA in laboratories is accomplished through ultra-low freezing of the nucleic acids. This however requires expensive freezers, convenient storage, reliable electrical power, and increased shipping costs, thereby making it a less viable option. Biomatrica (San Diego, CA) has created RNAstable(®), a stabilization reagent that is used to store RNA in a dehydrated state at room temperature (RT) and protects the RNA from degradation. Our objective was to investigate the sequence integrity and suitability of RNA when stored in RNAstable at extended time periods and at varying temperatures through use of Illumina and Agilent RNA expression microarrays. We observed in Bioanalyzer electropherograms that total RNA extracted from 293 cells stored at RT in RNAstable for 4.5 and 11.5 months is similar in quality to RNA stored at -80°C. Illumina mRNA expression array QC metrics and gene expression patterns from RNAstable-protected RNA, in contrast to RNA stored without RNAstable, correlated well with those of freezer controls. Significantly, when RNA was stored in RNAstable at 45°C for 4.5 months, equivalent to 22 months RT storage, RNA quality, microarray probe signal intensities, probe detection rates, and expression profiles remained similar between RNAstable-protected RNA at RT and the -80°C controls. At 10.5 months, miRNA levels were compared among the storage conditions using miRNA expression arrays. Here too we found strong concordance between miRNA expression patterns when total RNA was stored in RNAstable or at -80°C. Further, Bioanalyzer electrophoresis of RNAstable-protected samples stored at RT for a relative total of 33 months or 50.5 months showed comparable integrity scores to those of -80°C controls. We conclude that use of RNAstable holds promise as an effective stabilization reagent for total RNA and should be useful in situations where shipping and storage options are limited resources.