Identification of Global DNA Methylation Signatures in Glioblastoma-Derived Cancer Stem Cells

Glioblastoma(GBM)is the most common and most aggressive primary brain tumor in adults.The existence of a small population of stem-like tumor cells that efficiently propagate tumors and resist cytotoxic therapy is one proposed mechanism leading to the resilient behavior of tumor cells and poor prognosis.In this study,we performed an in-depth analysis of the DNA methylation landscape in GBMderived cancer stem cells(GSCs).Parallel comparisons of primary tumors and GSC lines derived from these tumors with normal controls(a neural stem cell(NSC)line and normal brain tissue)identified groups of hyper- and hypomethylated genes that display a trend of either increasing or decreasing methylation levels in the order of controls,primary GBMs,and their counterpart GSC lines,respectively.Interestingly,concurrent promoter hypermethylation and gene body hypomethylation were observed in a subset of genes including MGMT,AJAP1 and PTPRN2.These unique DNA methylation signatures were also found in primary GBM-derived xenograft tumors indicating that they are not tissue culture-related epigenetic changes.Integration of GSC-specific epigenetic signatures with gene expression analysis further identified candidate tumor suppressor genes that are frequently down-regulated in GBMs such as SPINT2,NEFM and PENK.Forced re-expression of SPINT2 reduced glioma cell proliferative capacity,anchorage independent growth,cell motility,and tumor sphere formation in vitro.The results from this study demonstrate that GSCs possess unique epigenetic signatures that may play important roles in the pathogenesis of GBM.     

Integrative analysis identifies potential DNA methylation biomarkers for pan-cancer diagnosis and prognosis

DNA methylation status is closely associated with diverse diseases, and is generally more stable than gene expression, thus abnormal DNA methylation could be important biomarkers for tumor diagnosis, treatment and prognosis. However, the signatures regarding DNA methylation changes for pan-cancer diagnosis and prognosis are less explored. Here we systematically analyzed the genome-wide DNA methylation patterns in diverse TCGA cancers with machine learning. We identified seven CpG sites that could effectively discriminate tumor samples from adjacent normal tissue samples for 12 main cancers of TCGA (1216 samples, AUC > 0.99). Those seven potential diagnostic biomarkers were further validated in the other 9 different TCGA cancers and 4 independent datasets (AUC > 0.92). Three out of the seven CpG sites were correlated with cell division, DNA replication and cell cycle. We also identified 12 CpG sites that can effectively distinguish 26 different cancers (7605 samples), and the result was repeatable in independent datasets as well as two disparate tumors with metastases (micro-average AUC > 0.89). Furthermore, a series of potential signatures that could significantly predict the prognosis of tumor patients for 7 different cancer were identified via survival analysis (p-value < 1e-4). Collectively, DNA methylation patterns vary greatly between tumor and adjacent normal tissues, as well as among different types of cancers. Our identified signatures may aid the decision of clinical diagnosis and prognosis for pan-cancer and the potential cancer-specific biomarkers could be used to predict the primary site of metastatic breast and prostate cancers.     

Integrative epigenetic and genetic pan-cancer somatic alteration portraits

Genetic and epigenetic alterations are required for carcinogenesis and the mutation burden across tumor types has been investigated. Here, we investigate epigenetic alterations with a novel measure of global DNA methylation dysregulation, the methylation dysregulation index (MDI), across 14 cancer types in The Cancer Genome Atlas (TCGA) database. DNA methylation data—obtained using Illumina HumanMethylation450 BeadChip—was accessed from TCGA. We calculated the MDI in 14 tumor types (n = 5,592 tumors), using adjacent normal tissues (n = 701) from each tumor site. Copy number alteration, and mutation burden were retrieved from cBioportal (n = 5,152). We tested the relation of subject MDI across tumors and with age, gender, tumor stage, estimated tumor purity, and copy number alterations for both overall MDI and genomic-context-specific MDI. We also investigated the top most dysregulated loci shared across tumor types. There was a broad range of extent in methylation dysregulation across tumor types (P < 2.2E-16). However, a consistent pattern of methylation dysregulation stratified by genomic context was observed across tumor types where the highest dysregulation occurred at non-CpG island regions. Considering other summary measures of somatic alteration, MDI was correlated with copy number alterations but not with mutation burden. Using the top dysregulated CpG sites in common across tumors, 4 classes of cancer types were observed, and the functional consequences of these alterations to gene expression were confirmed. This work identified the global DNA methylation dysregulation patterns across 14 cancer types showing a higher impact for the non-CpG island areas. The most dysregulated loci across cancer types identified common clusters across cancer types that may have implications for future treatment and prevention measures.     

Network Signatures of Survival in Glioblastoma Multiforme

To determine a molecular basis for prognostic differences in glioblastoma multiforme (GBM), we employed a combinatorial network analysis framework to exhaustively search for molecular patterns in protein-protein interaction (PPI) networks. We identified a dysregulated molecular signature distinguishing short-term (survival<225 days) from long-term (survival>635 days) survivors of GBM using whole genome expression data from The Cancer Genome Atlas (TCGA). A 50-gene subnetwork signature achieved 80% prediction accuracy when tested against an independent gene expression dataset. Functional annotations for the subnetwork signature included “protein kinase cascade,” “IκB kinase/NFκB cascade,” and “regulation of programmed cell death” – all of which were not significant in signatures of existing subtypes. Finally, we used label-free proteomics to examine how our subnetwork signature predicted protein level expression differences in an independent GBM cohort of 16 patients. We found that the genes discovered using network biology had a higher probability of dysregulated protein expression than either genes exhibiting individual differential expression or genes derived from known GBM subtypes. In particular, the long-term survivor subtype was characterized by increased protein expression of DNM1 and MAPK1 and decreased expression of HSPA9, PSMD3, and CANX. Overall, we demonstrate that the combinatorial analysis of gene expression data constrained by PPIs outlines an approach for the discovery of robust and translatable molecular signatures in GBM.     

Comparative analyses identify molecular signature of MRI-classified SVZ-associated glioblastoma

Glioblastoma (GBM) is a highly aggressive brain cancer with limited therapeutic options. While efforts to identify genes responsible for GBM have revealed mutations and aberrant gene expression associated with distinct types of GBM, patients with GBM are often diagnosed and classified based on MRI features. Therefore, we seek to identify molecular representatives in parallel with MRI classification for group I and group II primary GBM associated with the subventricular zone (SVZ). As group I and II GBM contain stem-like signature, we compared gene expression profiles between these 2 groups of primary GBM and endogenous neural stem progenitor cells to reveal dysregulation of cell cycle, chromatin status, cellular morphogenesis, and signaling pathways in these 2 types of MRI-classified GBM. In the absence of IDH mutation, several genes associated with metabolism are differentially expressed in these subtypes of primary GBM, implicating metabolic reprogramming occurs in tumor microenvironment. Furthermore, histone lysine methyltransferase EZH2 was upregulated while histone lysine demethylases KDM2 and KDM4 were downregulated in both group I and II primary GBM. Lastly, we identified 9 common genes across large data sets of gene expression profiles among MRI-classified group I/II GBM, a large cohort of GBM subtypes from TCGA, and glioma stem cells by unsupervised clustering comparison. These commonly upregulated genes have known functions in cell cycle, centromere assembly, chromosome segregation, and mitotic progression. Our findings highlight altered expression of genes important in chromosome integrity across all GBM, suggesting a common mechanism of disrupted fidelity of chromosome structure in GBM.     

Molecular Predictors of Long-Term Survival in Glioblastoma Multiforme Patients

Glioblastoma multiforme (GBM) is the most common and aggressive adult primary brain cancer, with <10% of patients surviving for more than 3 years. Demographic and clinical factors (e.g. age) and individual molecular biomarkers have been associated with prolonged survival in GBM patients. However, comprehensive systems-level analyses of molecular profiles associated with long-term survival (LTS) in GBM patients are still lacking. We present an integrative study of molecular data and clinical variables in these long-term survivors (LTSs, patients surviving >3 years) to identify biomarkers associated with prolonged survival, and to assess the possible similarity of molecular characteristics between LGG and LTS GBM. We analyzed the relationship between multivariable molecular data and LTS in GBM patients from the Cancer Genome Atlas (TCGA), including germline and somatic point mutation, gene expression, DNA methylation, copy number variation (CNV) and microRNA (miRNA) expression using logistic regression models. The molecular relationship between GBM LTS and LGG tumors was examined through cluster analysis. We identified 13, 94, 43, 29, and 1 significant predictors of LTS using Lasso logistic regression from the somatic point mutation, gene expression, DNA methylation, CNV, and miRNA expression data sets, respectively. Individually, DNA methylation provided the best prediction performance (AUC = 0.84). Combining multiple classes of molecular data into joint regression models did not improve prediction accuracy, but did identify additional genes that were not significantly predictive in individual models. PCA and clustering analyses showed that GBM LTS typically had gene expression profiles similar to non-LTS GBM. Furthermore, cluster analysis did not identify a close affinity between LTS GBM and LGG, nor did we find a significant association between LTS and secondary GBM. The absence of unique LTS profiles and the lack of similarity between LTS GBM and LGG, indicates that there are multiple genetic and epigenetic pathways to LTS in GBM patients.     

G-Cimp Status Prediction Of Glioblastoma Samples Using mRNA Expression Data

Glioblastoma Multiforme (GBM) is a tumor with high mortality and no known cure. The dramatic molecular and clinical heterogeneity seen in this tumor has led to attempts to define genetically similar subgroups of GBM with the hope of developing tumor specific therapies targeted to the unique biology within each of these subgroups. Recently, a subset of relatively favorable prognosis GBMs has been identified. These glioma CpG island methylator phenotype, or G-CIMP tumors, have distinct genomic copy number aberrations, DNA methylation patterns, and (mRNA) expression profiles compared to other GBMs. While the standard method for identifying G-CIMP tumors is based on genome-wide DNA methylation data, such data is often not available compared to the more widely available gene expression data. In this study, we have developed and evaluated a method to predict the G-CIMP status of GBM samples based solely on gene expression data.     

Bone Marrow Stromal Antigen 2 (BST-2) DNA Is Demethylated in Breast Tumors and Breast Cancer Cells

Background

Bone marrow stromal antigen 2 (BST-2) is a known anti-viral gene that has been recently identified to be overexpressed in many cancers, including breast cancer. BST-2 is critical for the invasiveness of breast cancer cells and the formation of metastasis in vivo. Although the regulation of BST-2 in immune cells is unraveling, it is unknown how BST-2 expression is regulated in breast cancer. We hypothesized that meta-analyses of BST-2 gene expression and BST-2 DNA methylation profiles would illuminate mechanisms regulating elevated BST-2 expression in breast tumor tissues and cells.

Materials and Methods

We performed comprehensive meta-analyses of BST-2 gene expression and BST-2 DNA methylation in The Cancer Genome Atlas (TCGA) and various Gene Expression Omnibus (GEO) datasets. BST-2 expression levels and BST-2 DNA methylation status at specific CpG sites on the BST-2 gene were compared for various breast tumor molecular subtypes and breast cancer cell lines.

Results

We show that BST-2 gene expression is inversely associated with the methylation status at specific CpG sites in primary breast cancer specimens and breast cancer cell lines. BST-2 demethylation is significantly more prevalent in primary tumors and cancer cells than in normal breast tissues or normal mammary epithelial cells. Demethylation of the BST-2 gene significantly correlates with its mRNA expression. These studies provide the initial evidence that significant differences exist in BST-2 DNA methylation patterns between breast tumors and normal breast tissues, and that BST-2 expression patterns in tumors and cancer cells correlate with hypomethylated BST-2 DNA.

Conclusion

Our study suggests that the DNA methylation pattern and expression of BST-2 may play a role in disease pathogenesis and could serve as a biomarker for the diagnosis of breast cancer.     

The Proneural Molecular Signature Is Enriched in Oligodendrogliomas and Predicts Improved Survival among Diffuse Gliomas

The Cancer Genome Atlas Project (TCGA) has produced an extensive collection of ‘-omic’ data on glioblastoma (GBM), resulting in several key insights on expression signatures. Despite the richness of TCGA GBM data, the absence of lower grade gliomas in this data set prevents analysis genes related to progression and the uncovering of predictive signatures. A complementary dataset exists in the form of the NCI Repository for Molecular Brain Neoplasia Data (Rembrandt), which contains molecular and clinical data for diffuse gliomas across the full spectrum of histologic class and grade. Here we present an investigation of the significance of the TCGA consortium''s expression classification when applied to Rembrandt gliomas. We demonstrate that the proneural signature predicts improved clinical outcome among 176 Rembrandt gliomas that includes all histologies and grades, including GBMs (log rank test p = 1.16e-6), but also among 75 grade II and grade III samples (p = 2.65e-4). This gene expression signature was enriched in tumors with oligodendroglioma histology and also predicted improved survival in this tumor type (n = 43, p = 1.25e-4). Thus, expression signatures identified in the TCGA analysis of GBMs also have intrinsic prognostic value for lower grade oligodendrogliomas, and likely represent important differences in tumor biology with implications for treatment and therapy. Integrated DNA and RNA analysis of low-grade and high-grade proneural gliomas identified increased expression and gene amplification of several genes including GLIS3, TGFB2, TNC, AURKA, and VEGFA in proneural GBMs, with corresponding loss of DLL3 and HEY2. Pathway analysis highlights the importance of the Notch and Hedgehog pathways in the proneural subtype. This demonstrates that the expression signatures identified in the TCGA analysis of GBMs also have intrinsic prognostic value for low-grade oligodendrogliomas, and likely represent important differences in tumor biology with implications for treatment and therapy.     

DNA Methylation Patterns Can Estimate Nonequivalent Outcomes of Breast Cancer with the Same Receptor Subtypes

Breast cancer has various molecular subtypes and displays high heterogeneity. Aberrant DNA methylation is involved in tumor origin, development and progression. Moreover, distinct DNA methylation patterns are associated with specific breast cancer subtypes. We explored DNA methylation patterns in association with gene expression to assess their impact on the prognosis of breast cancer based on Infinium 450K arrays (training set) from The Cancer Genome Atlas (TCGA). The DNA methylation patterns of 12 featured genes that had a high correlation with gene expression were identified through univariate and multivariable Cox proportional hazards models and used to define the methylation risk score (MRS). An improved ability to distinguish the power of the DNA methylation pattern from the 12 featured genes (p = 0.00103) was observed compared with the average methylation levels (p = 0.956) or gene expression (p = 0.909). Furthermore, MRS provided a good prognostic value for breast cancers even when the patients had the same receptor status. We found that ER-, PR- or Her2- samples with high-MRS had the worst 5-year survival rate and overall survival time. An independent test set including 28 patients with death as an outcome was used to test the validity of the MRS of the 12 featured genes; this analysis obtained a prognostic value equivalent to the training set. The predict power was validated through two independent datasets from the GEO database. The DNA methylation pattern is a powerful predictor of breast cancer survival, and can predict outcomes of the same breast cancer molecular subtypes.     

Comparison of glioma stem cells to neural stem cells from the adult human brain identifies dysregulated Wnt- signaling and a fingerprint associated with clinical outcome

Glioblastoma is the most common brain tumor. Median survival in unselected patients is <10 months. The tumor harbors stem-like cells that self-renew and propagate upon serial transplantation in mice, although the clinical relevance of these cells has not been well documented. We have performed the first genome-wide analysis that directly relates the gene expression profile of nine enriched populations of glioblastoma stem cells (GSCs) to five identically isolated and cultivated populations of stem cells from the normal adult human brain. Although the two cell types share common stem- and lineage-related markers, GSCs show a more heterogeneous gene expression. We identified a number of pathways that are dysregulated in GSCs. A subset of these pathways has previously been identified in leukemic stem cells, suggesting that cancer stem cells of different origin may have common features. Genes upregulated in GSCs were also highly expressed in embryonic and induced pluripotent stem cells. We found that canonical Wnt-signaling plays an important role in GSCs, but not in adult human neural stem cells. As well we identified a 30-gene signature highly overexpressed in GSCs. The expression of these signature genes correlates with clinical outcome and demonstrates the clinical relevance of GSCs.     

Hypoxia-induced phenotypic transition from highly invasive to less invasive tumors in glioma stem-like cells: Significance of CD44 and osteopontin as therapeutic targets in glioblastoma

The poor prognosis of glioblastoma multiforme (GBM) is primarily due to highly invasive glioma stem-like cells (GSCs) in tumors. Upon GBM recurrence, GSCs with highly invasive and highly migratory activities must assume a less-motile state and proliferate to regenerate tumor mass. Elucidating the molecular mechanism underlying this transition from a highly invasive phenotype to a less-invasive, proliferative tumor could facilitate the identification of effective molecular targets for treating GBM. Here, we demonstrate that severe hypoxia (1% O₂) upregulates CD44 expression via activation of hypoxia-inducible factor (HIF-1α), inducing GSCs to assume a highly invasive tumor. In contrast, moderate hypoxia (5% O₂) upregulates osteopontin expression via activation of HIF-2α. The upregulated osteopontin inhibits CD44-promoted GSC migration and invasion and stimulates GSC proliferation, inducing GSCs to assume a less-invasive, highly proliferative tumor. These data indicate that the GSC phenotype is determined by interaction between CD44 and osteopontin. The expression of both CD44 and osteopontin is regulated by differential hypoxia levels. We found that CD44 knockdown significantly inhibited GSC migration and invasion both in vitro and in vivo. Mouse brain tumors generated from CD44-knockdown GSCs exhibited diminished invasiveness, and the mice survived significantly longer than control mice. In contrast, siRNA-mediated silencing of the osteopontin gene decreased GSC proliferation. These results suggest that interaction between CD44 and osteopontin plays a key role in tumor progression in GBM; inhibition of both CD44 and osteopontin may represent an effective therapeutic approach for suppressing tumor progression, thus resulting in a better prognosis for patients with GBM.