Tensor GSVD of Patient- and Platform-Matched Tumor and Normal DNA Copy-Number Profiles Uncovers Chromosome Arm-Wide Patterns of Tumor-Exclusive Platform-Consistent Alterations Encoding for Cell Transformation and Predicting Ovarian Cancer Survival

The number of large-scale high-dimensional datasets recording different aspects of a single disease is growing, accompanied by a need for frameworks that can create one coherent model from multiple tensors of matched columns, e.g., patients and platforms, but independent rows, e.g., probes. We define and prove the mathematical properties of a novel tensor generalized singular value decomposition (GSVD), which can simultaneously find the similarities and dissimilarities, i.e., patterns of varying relative significance, between any two such tensors. We demonstrate the tensor GSVD in comparative modeling of patient- and platform-matched but probe-independent ovarian serous cystadenocarcinoma (OV) tumor, mostly high-grade, and normal DNA copy-number profiles, across each chromosome arm, and combination of two arms, separately. The modeling uncovers previously unrecognized patterns of tumor-exclusive platform-consistent co-occurring copy-number alterations (CNAs). We find, first, and validate that each of the patterns across only 7p and Xq, and the combination of 6p+12p, is correlated with a patient’s prognosis, is independent of the tumor’s stage, the best predictor of OV survival to date, and together with stage makes a better predictor than stage alone. Second, these patterns include most known OV-associated CNAs that map to these chromosome arms, as well as several previously unreported, yet frequent focal CNAs. Third, differential mRNA, microRNA, and protein expression consistently map to the DNA CNAs. A coherent picture emerges for each pattern, suggesting roles for the CNAs in OV pathogenesis and personalized therapy. In 6p+12p, deletion of the p21-encoding CDKN1A and p38-encoding MAPK14 and amplification of RAD51AP1 and KRAS encode for human cell transformation, and are correlated with a cell’s immortality, and a patient’s shorter survival time. In 7p, RPA3 deletion and POLD2 amplification are correlated with DNA stability, and a longer survival. In Xq, PABPC5 deletion and BCAP31 amplification are correlated with a cellular immune response, and a longer survival.     

Focal chromosomal copy number aberrations in cancer—Needles in a genome haystack

The extent of focal chromosomal copy number aberrations (CNAs) in cancer has been uncovered through technical innovations, and this discovery has been critical for the identification of new cancer driver genes in genomics projects such as TCGA and ICGC. Unlike constitutive copy number variations (CNVs), focal CNAs are the result of many selection events during the evolution of cancer genomes. Therefore, it is possible that a single gene in a focal CNA gives the tumor a selective growth advantage. This concept has been instrumental in the discovery of new cancer driver genes. However, focal CNAs lack a consensus definition; therefore, we propose one based on pragmatic considerations. We also describe different strategies to identify focal CNAs and procedures to distinguish them from large CNAs and CNVs.     

自噬基因CTSL表达与胶质母细胞瘤患者预后相关

本研究旨在探讨自噬基因CTSL对胶质母细胞瘤(GBM)患者的预后影响。利用癌症基因组图谱(TCGA)、人类自噬数据库(HADB)、中国脑胶质瘤基因组图谱(CGGA)数据库、基因表达谱分析(GEPIA)获取数据信息,通过筛选差异表达基因及单因素和多因素COX分析确定GBM的独立预后危险因素,同时通过基因本体论(GO)、基因组百科全书途径(KEGG)、临床病理相关性、基因集富集分析(GSEA)、自噬基因网络分析CTSL的相关作用机制。结果显示:(1)富集分析显示胶质母细胞瘤中差异自噬基因(ARG)与自噬体的形成、细胞凋亡、血管生成、细胞化疗等相关;(2)GBM中CTSL的mRNA水平明显高于正常组织样本;(3)多因素COX回归分析显示自噬基因CTSL的高表达为GBM预后的独立危险因素,STUPP治疗(术后替莫唑胺[Tmz]同步放化疗+Tmz辅助化疗)为独立保护因素;(4)自噬基因CTSL在非GCIMP(CpG岛甲基化)型、间质型、IDH野生型、1p/19q无缺失型胶质母细胞瘤及化疗后表达量更高。综上所述,本研究分析了自噬基因在GBM中的作用,并表明自噬基因CTSL的过表达预示胶质母细胞瘤患者不良预后,显示自噬基因CTSL有作为有效靶标的潜质。     

SAMSN1 Is Highly Expressed and Associated with a Poor Survival in Glioblastoma Multiforme

Objectives

To study the expression pattern and prognostic significance of SAMSN1 in glioma.

Methods

Affymetrix and Arrystar gene microarray data in the setting of glioma was analyzed to preliminarily study the expression pattern of SAMSN1 in glioma tissues, and Hieratical clustering of gene microarray data was performed to filter out genes that have prognostic value in malignant glioma. Survival analysis by Kaplan-Meier estimates stratified by SAMSN1 expression was then made based on the data of more than 500 GBM cases provided by The Cancer Genome Atlas (TCGA) project. At last, we detected the expression of SAMSN1 in large numbers of glioma and normal brain tissue samples using Tissue Microarray (TMA). Survival analysis by Kaplan-Meier estimates in each grade of glioma was stratified by SAMSN1 expression. Multivariate survival analysis was made by Cox proportional hazards regression models in corresponding groups of glioma.

Results

With the expression data of SAMSN1 and 68 other genes, high-grade glioma could be classified into two groups with clearly different prognoses. Gene and large sample tissue microarrays showed high expression of SAMSN1 in glioma particularly in GBM. Survival analysis based on the TCGA GBM data matrix and TMA multi-grade glioma dataset found that SAMSN1 expression was closely related to the prognosis of GBM, either PFS or OS (P<0.05). Multivariate survival analysis with Cox proportional hazards regression models confirmed that high expression of SAMSN1 was a strong risk factor for PFS and OS of GBM patients.

Conclusion

SAMSN1 is over-expressed in glioma as compared with that found in normal brains, especially in GBM. High expression of SAMSN1 is a significant risk factor for the progression free and overall survival of 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.     

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.     

Genome-wide copy number profiling on high-density bacterial artificial chromosomes, single-nucleotide polymorphisms, and oligonucleotide microarrays: a platform comparison based on statistical power analysis.

Recently, comparative genomic hybridization onto bacterial artificial chromosome (BAC) arrays (array-based comparative genomic hybridization) has proved to be successful for the detection of submicroscopic DNA copy-number variations in health and disease. Technological improvements to achieve a higher resolution have resulted in the generation of additional microarray platforms encompassing larger numbers of shorter DNA targets (oligonucleotides). Here, we present a novel method to estimate the ability of a microarray to detect genomic copy-number variations of different sizes and types (i.e. deletions or duplications). We applied our method, which is based on statistical power analysis, to four widely used high-density genomic microarray platforms. By doing so, we found that the high-density oligonucleotide platforms are superior to the BAC platform for the genome-wide detection of copy-number variations smaller than 1 Mb. The capacity to reliably detect single copy-number variations below 100 kb, however, appeared to be limited for all platforms tested. In addition, our analysis revealed an unexpected platform-dependent difference in sensitivity to detect a single copy-number loss and a single copy-number gain. These analyses provide a first objective insight into the true capacities and limitations of different genomic microarrays to detect and define DNA copy-number variations.     

Development and validation of a CTNNB1-associated metabolic prognostic model for hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is a heterogeneous malignancy closely related to metabolic reprogramming. We investigated how CTNNB1 mutation regulates the HCC metabolic phenotype and thus affects the prognosis of HCC. We obtained the mRNA expression profiles and clinicopathological data from The Cancer Genome Atlas (TCGA), the International Cancer Genomics Consortium (ICGC) and the Gene Expression Omnibus database ( GSE14520 and GSE116174 ). We conducted gene set enrichment analysis on HCC patients with and without mutant CTNNB1 through TCGA dataset. The Kaplan-Meier analysis and univariate Cox regression analysis assisted in screening metabolic genes related to prognosis, and the prognosis model was constructed using the Lasso and multivariate Cox regression analysis. The prognostic model showed good prediction performance in both the training cohort (TCGA) and the validation cohorts (ICGC, GSE14520 , GSE116174 ), and the high-risk group presented obviously poorer overall survival compared with low-risk group. Cox regression analysis indicated that the risk score can be used as an independent predictor for the overall survival of HCC. The immune infiltration in different risk groups was also evaluated in this study to explore underlying mechanisms. This study is also the first to describe an metabolic prognostic model associated with CTNNB1 mutations and could be implemented for determining the prognoses of individual patients in clinical practice.     

A ten-microRNA expression signature predicts survival in glioblastoma

Glioblastoma (GBM) is the most common and aggressive primary brain tumor with very poor patient median survival. To identify a microRNA (miRNA) expression signature that can predict GBM patient survival, we analyzed the miRNA expression data of GBM patients (n=222) derived from The Cancer Genome Atlas (TCGA) dataset. We divided the patients randomly into training and testing sets with equal number in each group. We identified 10 significant miRNAs using Cox regression analysis on the training set and formulated a risk score based on the expression signature of these miRNAs that segregated the patients into high and low risk groups with significantly different survival times (hazard ratio [HR]=2.4; 95% CI=1.4-3.8; p<0.0001). Of these 10 miRNAs, 7 were found to be risky miRNAs and 3 were found to be protective. This signature was independently validated in the testing set (HR=1.7; 95% CI=1.1-2.8; p=0.002). GBM patients with high risk scores had overall poor survival compared to the patients with low risk scores. Overall survival among the entire patient set was 35.0% at 2 years, 21.5% at 3 years, 18.5% at 4 years and 11.8% at 5 years in the low risk group, versus 11.0%, 5.5%, 0.0 and 0.0% respectively in the high risk group (HR=2.0; 95% CI=1.4-2.8; p<0.0001). Cox multivariate analysis with patient age as a covariate on the entire patient set identified risk score based on the 10 miRNA expression signature to be an independent predictor of patient survival (HR=1.120; 95% CI=1.04-1.20; p=0.003). Thus we have identified a miRNA expression signature that can predict GBM patient survival. These findings may have implications in the understanding of gliomagenesis, development of targeted therapy and selection of high risk cancer patients for adjuvant therapy.     

A robust two-gene signature for glioblastoma survival prediction

Glioblastoma multiforme (GBM) is a highly malignant brain tumor. We explored the prognostic gene signature in 443 GBM samples by systematic bioinformatics analysis, using GSE16011 with microarray expression and corresponding clinical data from Gene Expression Omnibus as the training set. Meanwhile, patients from The Chinese Glioma Genome Atlas database (CGGA) were used as the test set and The Cancer Genome Atlas database (TCGA) as the validation set. Through Cox regression analysis, Kaplan-Meier analysis, t-distributed Stochastic Neighbor Embedding algorithm, clustering, and receiver operating characteristic analysis, a two-gene signature (GRIA2 and RYR3) associated with survival was selected in the GSE16011 dataset. The GRIA2-RYR3 signature divided patients into two risk groups with significantly different survival in the GSE16011 dataset (median: 0.72, 95% confidence interval [CI]: 0.64-0.98, vs median: 0.98, 95% CI: 0.65-1.61 years, logrank test P < .001), the CGGA dataset (median: 0.84, 95% CI: 0.70-1.18, vs median: 1.21, 95% CI: 0.95-2.94 years, logrank test P = .0017), and the TCGA dataset (median: 1.03, 95% CI: 0.86-1.24, vs median: 1.23, 95% CI: 1.04-1.85 years, logrank test P = .0064), validating the predictive value of the signature. And the survival predictive potency of the signature was independent from clinicopathological prognostic features in multivariable Cox analysis. We found that after transfection of U87 cells with small interfering RNA, GRIA2 and RYR3 influenced the biological behaviors of proliferation, migration, and invasion of glioblastoma cells. In conclusion, the two-gene signature was a robust prognostic model to predict GBM survival.     

Detection of DNA copy number abnormality by microarray expression analysis

Gene copy-number abnormalities (CNAs) are characteristic of solid tumors and are found in association with developmental abnormalities and/or mental retardation. The ultimate impact of CNAs is exerted by the altered expression of encoded genes. We have utilized high-density oligonucleotide arrays from Affymetrix to identify DNA CNAs via their impact on mRNA expression levels. In these studies, we have used three different trisomic cell lines (trisomy 9, trisomy 18, trisomy 21) as models of CNAs and have compared mRNA expression in those trisomic cells with that observed in diploid cell lines of matched tissue origin. Our data clearly show that genes from CNA chromosome regions are substantially over-represented (P<0.000001 by chi-square analysis) in the differentially expressed subset from comparisons of all three trisomic cell lines with normal matching cells. In addition, we have been able to detect the origin of the duplication by a statistical scan for over-expressed genes. These data show that microarray detection of differential mRNA expression can be used to identify significant DNA CNAs.     

The global DNA methylation surrogate LINE-1 methylation is correlated with MGMT promoter methylation and is a better prognostic factor for glioma

Gliomas are the most frequently occurring primary brain tumor in the central nervous system of adults. Glioblastoma multiformes (GBMs, WHO grade 4) have a dismal prognosis despite the use of the alkylating agent, temozolomide (TMZ), and even low grade gliomas (LGGs, WHO grade 2) eventually transform to malignant secondary GBMs. Although GBM patients benefit from promoter hypermethylation of the O(6)-methylguanine-DNA methyltransferase (MGMT) that is the main determinant of resistance to TMZ, recent studies suggested that MGMT promoter methylation is of prognostic as well as predictive significance for the efficacy of TMZ. Glioma-CpG island methylator phenotype (G-CIMP) in the global genome was shown to be a significant predictor of improved survival in patients with GBM. Collectively, we hypothesized that MGMT promoter methylation might reflect global DNA methylation. Additionally in LGGs, the significance of MGMT promoter methylation is still undetermined. In the current study, we aimed to determine the correlation between clinical, genetic, and epigenetic profiles including LINE-1 and different cancer-related genes and the clinical outcome in newly diagnosed 57 LGG and 54 GBM patients. Here, we demonstrated that (1) IDH1/2 mutation is closely correlated with MGMT promoter methylation and 1p/19q codeletion in LGGs, (2) LINE-1 methylation levels in primary and secondary GBMs are lower than those in LGGs and normal brain tissues, (3) LINE-1 methylation is proportional to MGMT promoter methylation in gliomas, and (4) higher LINE-1 methylation is a favorable prognostic factor in primary GBMs, even compared to MGMT promoter methylation. As a global DNA methylation marker, LINE-1 may be a promising marker in gliomas.     

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.     

Prognostic power of a lipid metabolism gene panel for diffuse gliomas

Lipid metabolism reprogramming plays important role in cell growth, proliferation, angiogenesis and invasion in cancers. However, the diverse lipid metabolism programmes and prognostic value during glioma progression remain unclear. Here, the lipid metabolism‐related genes were profiled using RNA sequencing data from The Cancer Genome Atlas (TCGA) and Chinese Glioma Genome Atlas (CGGA) database. Gene ontology (GO) and gene set enrichment analysis (GSEA) found that glioblastoma (GBM) mainly exhibited enrichment of glycosphingolipid metabolic progress, whereas lower grade gliomas (LGGs) showed enrichment of phosphatidylinositol metabolic progress. According to the differential genes of lipid metabolism between LGG and GBM, we developed a nine‐gene set using Cox proportional hazards model with elastic net penalty, and the CGGA cohort was used for validation data set. Survival analysis revealed that the obtained gene set could differentiate the outcome of low‐ and high‐risk patients in both cohorts. Meanwhile, multivariate Cox regression analysis indicated that this signature was a significantly independent prognostic factor in diffuse gliomas. Gene ontology and GSEA showed that high‐risk cases were associated with phenotypes of cell division and immune response. Collectively, our findings provided a new sight on lipid metabolism in diffuse gliomas.     

GISTIC2.0 facilitates sensitive and confident localization of the targets of focal somatic copy-number alteration in human cancers

We describe methods with enhanced power and specificity to identify genes targeted by somatic copy-number alterations (SCNAs) that drive cancer growth. By separating SCNA profiles into underlying arm-level and focal alterations, we improve the estimation of background rates for each category. We additionally describe a probabilistic method for defining the boundaries of selected-for SCNA regions with user-defined confidence. Here we detail this revised computational approach, GISTIC2.0, and validate its performance in real and simulated datasets.     

Overexpression of oncostatin M receptor regulates local immune response in glioblastoma

Glioblastoma (GBM) is the most lethal cancer in central nervous system. It is urgently needed to look for novel therapeutics for GBM. Oncostatin M receptor (OSMR) is a cytokine receptor gene of IL-6 family and has been reported to be involved in regulating GBM tumorigenesis. However, the role of OSMR regulating the disrupted immune response in GBM need to be further investigated. Three gene expression profiles, Chinese Glioma Genome Atlas (CGGA), The Cancer Genome Atlas (TCGA), and Gene Expression Omnibus (GEO) data set (GSE16011), were enrolled in our study and used for OSMR expression and survival analysis. The expression of OSMR was further verified with immunohistochemistry and western blot analysis in glioma tissues. Microenvironment cell populations-counter (MCP-counter) was applied for analyzing the relationship between OSMR expression and nontumor cells. The functions of OSMR in GBM was investigated by Gene Ontology, Gene set enrichment analysis (GSEA), gene set variation analysis and so on. The analysis of cytokine receptor activity-related genes in glioma identifies OSMR as a gene with an independent predictive factor for progressive malignancy in GBM. Furthermore, OSMR expression is a prognostic marker in the response prediction to radiotherapy and chemotherapy. OSMR contributes to the regulation of local immune response and extracellular matrix process in GBM. Our findings define an important role of OSMR in the regulation of local immune response in GBM, which may suggest OSMR as a possible biomarker in developing new therapeutic immune strategies in GBM.