The joint role of methylation and immune-related lncRNAs in ovarian cancer: Defining molecular subtypes and developing prognostic signature

IntroductionComplex outcome of ovarian cancer (OC) stems from the tumor immune microenvironment (TIME) influenced by genetic and epigenetic factors. This study aimed to comprehensively explored the subclasses of OC through lncRNAs related to both N6-methyladenosine (m6A)/N1-methyladenosine (m1A)/N7-methylguanosine (m7G)/5-methylcytosine (m5C) in terms of epigenetic variability and immune molecules and develop a new set of risk predictive systems.Material and methodsThe lncRNA data of OC were collected from TCGA. Spearman correlation analysis on lncRNA data of OC with immune-related gene expression and with m6A/m5C/m1A/m7G were respectively conducted. The m6A/m5C/m1A/m7G-related m6A/m5C/m1A/m7G related immune lncRNA subtypes were identified on the basis of the prognostic lncRNAs. Heterogeneity among subtypes was evaluated by tumor mutation analysis, tumor microenvironment (TME) component analysis, response to immune checkpoint blocked (ICB) and chemotherapeutic drugs. A risk predictive system was developed based on the results of Cox regression analysis and random survival forest analysis of the differences between each specific cluster and other clusters.ResultsThree m6A/m5C/m1A/m7G-related immune lncRNA subtypes of OC showing distinct differences in prognosis, mutation pattern, TIME components, immunotherapy and chemotherapy response were identified. A set of risk predictive system consisting of 10 lncRNA for OC was developed, according to which the risk score of samples in each OC dataset was calculated and risk type was defined.ConclusionsThis study classified three m6A/m5C/m1A/m7G-related immune lncRNA subtypes with distinct heterogeneous mutation patterns, TME components, ICB therapy and immune response, and provided a set of risk predictive system consisted of 10 lncRNA for OC.     

Identification of immune-enhanced molecular subtype associated with BRCA1 mutations,immune checkpoints and clinical outcome in ovarian carcinoma

Ovarian carcinoma has the highest mortality among the malignant tumours in gynaecology, and new treatment strategies are urgently needed to improve the clinical status of ovarian carcinoma patients. The Cancer Genome Atlas (TCGA) cohort were performed to explore the immune function of the internal environment of tumours and its clinical correlation with ovarian carcinoma. Finally, four molecular subtypes were obtained based on the global immune-related genes. The correlation analysis and clinical characteristics showed that four subtypes were all significantly related to clinical stage; the immune scoring results indicated that most immune signatures were upregulated in C3 subtype, and the majority of tumour-infiltrating immune cells were upregulated in both C3 and C4 subtypes. Compared with other subtypes, C3 subtype had a higher BRCA1 mutation, higher expression of immune checkpoints, and optimal survival prognosis. These findings of the immunological microenvironment in tumours may provide new ideas for developing immunotherapeutic strategies for ovarian carcinoma.     

Identification of immune microenvironment subtypes that predicted the prognosis of patients with ovarian cancer

Ovarian cancer (OC) is associated with high mortality rate. However, the correlation between immune microenvironment and prognosis of OC remains unclear. This study aimed to explore prognostic significance of OC tumour microenvironment. The OC data set was selected from the cancer genome atlas (TCGA), and 307 samples were collected. Hierarchical clustering was performed according to the expression of 756 genes. The immune and matrix scores of all immune subtypes were determined, and Kruskal-Wallis test was used to analyse the differences in the immune and matrix scores between OC samples with different immune subtypes. The model for predicting prognosis was constructed based on the expression of immune-related genes. TIDE platform was applied to predict the effect of immunotherapy on patients with OC of different immune subtypes. The 307 OC samples were classified into three immune subtypes A-C. Patients in subtype B had poorer prognosis and lower survival rate. The infiltration of helper T cells and macrophages in microenvironment indicated significant differences between immune subtypes. Enrichment analyses of immune cell molecular pathways showed that JAK–STAT3 pathway changed significantly in subtype B. Furthermore, predictive response to immunotherapy in subtype B was significantly higher than that in subtype A and C. Immune subtyping can be used as an independent predictor of the prognosis of OC patients, which may be related to the infiltration patterns of immune cells in tumour microenvironment. In addition, patients in immune subtype B have superior response to immunotherapy, suggesting that patients in subtype B are suitable for immunotherapy.     

Identification of an m6A-related lncRNA prognostic signature in colorectal cancer

Data sets of colorectal cancer (CRC) were obtained from The Cancer Genome Atlas (TCGA), three N6-methyladenosine (m6A) subtypes were identified using 21 m6A-related long noncoding RNAs (lncRNAs) and differential m6A subtypes of different CRC tumors were determined in this study to evaluate the m6A expression and the prognosis of patients with CRC. Subsequently, eight key lncRNAs were identified based on co-expression with 21 m6A-related genes in CRC tumors using the single-factor Cox and least absolute shrinkage and selection operator. Finally, an m6A-related lncRNA risk score model of CRC tumor was established using multifactor Cox regression based on the eight important lncRNAs and found to have a better performance in evaluating the prognosis of patients in the TCGA-CRC data set. TCGA-CRC tumor samples were divided based on the risk scores: high and low. Then, the clinical characteristics, tumor mutation load, and tumor immune cell infiltration difference between the high- and low-risk-score groups were explored, and the predictive ability of the risk score was assessed for immunotherapeutic benefits. We found that the risk score model can determine the overall survival, be a relatively independent prognostic indicator, and better evaluate the immunotherapeutic benefits for patients with CRC. This study provides data support for accurate immunotherapy in CRC.     

Multi-Omics Analysis Elucidates The Immune And Intratumor Microbes Characteristics Of Ubiquitination Subtypes In Lung Adenocarcinoma

Ubiquitination modification is closely related to cancer and participates in the regulation of tumor microenvironment. However, the role of ubiquitination modification in the immune response and prognosis of lung adenocarcinoma has not been elucidated. This study aims to establish a disease classification associated with ubiquitination and reveal the landscape of intratumor microbes in patients with lung adenocarcinoma for the first time. A total of 1314 patients with lung adenocarcinoma in the GEO and TCGA databases were included in our study. We constructed a ubiquitination scoring model using WGCNA and constructed ubiquitination subtypes using unsupervised clustering, analyzed the clinical characteristics, immune characteristics, and intratumor microbes characteristics, and screened out the relevant gene signatures, which were verified by RT-qPCR in human cancer cells. The results showed that the high ubiquitination subtype had poor prognosis, low degree of immune infiltration, high index of tumor stemness, and poor effect of immunotherapy. The subtypes with lower ubiquitination scores have better prognosis, higher tumor microenvironment score and better immunotherapy effect. The C2 subtype has high level of immune infiltration, lower intratumor microbes diversity and abundance, and good prognosis. The C3 subtype has low level of immune infiltration, higher intratumor microbes diversity and abundance, and poor prognosis. The C1 subtype has characteristics between C2 and C3. In summary, this paper constructs a scoring system and several subtypes based on ubiquitination genes, and analyzed the characteristics, which can help provide new methods for clinical treatment.     

Five key lncRNAs considered as prognostic targets for predicting pancreatic ductal adenocarcinoma

Pancreatic ductal adenocarcinoma (PDAC) has a poor prognosis, and the 5‐year survival rate was only 7.7%. To improve prognosis, a screening biomarker for early diagnosis of pancreatic cancer is in urgent need. Long non‐coding RNA (lncRNA) expression profiles as potential cancer prognostic biomarkers play critical roles in development of tumorigenesis and metastasis of cancer. However, lncRNA signatures in predicting the survival of a patient with PDAC remain unknown. In the current study, we try to identify potential lncRNA biomarkers and their prognostic values in PDAC. LncRNAs expression profiles and corresponding clinical information for 182 cases with PDAC were acquired from The Cancer Genome Atlas (TCGA). A total of 14 470 lncRNA were identified in the cohort, and 175 PDAC patients had clinical variables. We obtained 108 differential expressed lncRNA via R packages. Univariate and multivariate Cox proportional hazards regression, lasso regression was performed to screen the potential prognostic lncRNA. Five lncRNAs have been recognized to significantly correlate with OS. We established a linear prognostic model of five lncRNA (C9orf139, MIR600HG, RP5‐965G21.4, RP11‐436K8.1, and CTC‐327F10.4) and divided patients into high‐ and low‐risk group according to the prognostic index. The five lncRNAs played independent prognostic biomarkers of OS of PDAC patients and the AUC of the ROC curve for the five lncRNAs signatures prediction 5‐year survival was 0.742. In addition, targeted genes of MIR600HG, C9orf139, and CTC‐327F10.4 were explored and functional enrichment was also conducted. These results suggested that this five‐lncRNAs signature could act as potential prognostic biomarkers in the prediction of PDAC patient's survival.     

Dysregulated ceRNA network modulated by copy number variation-driven lncRNAs in breast cancer: A comprehensive analysis

Breast cancer is a malignancy harmful to physical and mental health in women, with quite high mortality. Copy number variations (CNVs) are vital factors affecting the progression of breast cancer. Detecting CNVs in breast cancer to predict the prognosis of patients has become a promising approach to accurate treatment in recent years. The differential analysis was performed on CNVs of long noncoding RNAs (lncRNAs) as well as the expression of lncRNAs, microRNAs (miRNAs) and mRNAs in normal tissue and breast tumor tissue based on The Cancer Genome Atlas (TCGA) database. The CNV-driven lncRNAs were identified by the Kruskal–Wallis test. Meanwhile, a competitive endogenous RNA (ceRNA) network regulated by CNV-driven lncRNA was constructed. As the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses revealed, the mRNAs in the dysregulated ceRNA network were mainly enriched in the biological functions and signaling pathways, including the Focal Adhesion-PI3K-Akt–mTOR-signaling pathway, the neuronal system, metapathway biotransformation Phase I and II and blood circulation, etc. The relationship between the CNVs of five lncRNAs and their gene expression in the ceRNA network was analyzed via a chi-square test, which confirmed that except for LINC00243, the expression of four lncRNAs was notably correlated with the CNVs. The survival analysis revealed that only the copy number gain of LINC00536 was evidently related to the poor prognosis of patients. The CIBERSORT algorithm showed that five lncRNAs were correlated with the abundance of immune cell infiltration and immune checkpoints. In a word, by analyzing CNV-driven lncRNAs and the ceRNA network regulated by these lncRNAs, this study explored the mechanism of breast cancer and provided novel insights into new biomarkers.     

Ovarian Cancer Cell Line Panel (OCCP): Clinical Importance of In Vitro Morphological Subtypes

Epithelial ovarian cancer is a highly heterogeneous disease and remains the most lethal gynaecological malignancy in the Western world. Therapeutic approaches need to account for inter-patient and intra-tumoural heterogeneity and detailed characterization of in vitro models representing the different histological and molecular ovarian cancer subtypes is critical to enable reliable preclinical testing. There are approximately 100 publicly available ovarian cancer cell lines but their cellular and molecular characteristics are largely undescribed. We have characterized 39 ovarian cancer cell lines under uniform conditions for growth characteristics, mRNA/microRNA expression, exon sequencing, drug response for clinically-relevant therapeutics and collated all available information on the original clinical features and site of origin. We tested for statistical associations between the cellular and molecular features of the lines and clinical features. Of the 39 ovarian cancer cell lines, 14 were assigned as high-grade serous, four serous-type, one low-grade serous and 20 non-serous type. Three morphological subtypes: Epithelial (n = 21), Round (n = 7) and Spindle (n = 12) were identified that showed distinct biological and molecular characteristics, including overexpression of cell movement and migration-associated genes in the Spindle subtype. Comparison with the original clinical data showed association of the spindle-like tumours with metastasis, advanced stage, suboptimal debulking and poor prognosis. In addition, the expression profiles of Spindle, Round and Epithelial morphologies clustered with the previously described C1-stromal, C5-mesenchymal and C4 ovarian subtype expression profiles respectively. Comprehensive profiling of 39 ovarian cancer cell lines under controlled, uniform conditions demonstrates clinically relevant cellular and genomic characteristics. This data provides a rational basis for selecting models to develop specific treatment approaches for histological and molecular subtypes of ovarian cancer.     

A prognostic 11 long noncoding RNA expression signature for breast invasive carcinoma

Breast cancer, the most common cancer in women worldwide, is associated with high mortality. The long non-coding RNAs (lncRNAs) with a little capacity of coding proteins is playing an increasingly important role in the cancer paradigm. Accumulating evidences demonstrate that lncRNAs have crucial connections with breast cancer prognosis while the studies of lncRNAs in breast cancer are still in its primary stage. In this study, we collected 1052 clinical patient samples, a comparatively large sample size, including 13 159 lncRNA expression profiles of breast invasive carcinoma (BRCA) from The Cancer Genome Atlas database to identify prognosis-related lncRNAs. We randomly separated all of these clinical patient samples into training and testing sets. In the training set, we performed univariable Cox regression analysis for primary screening and played the model for Robust likelihood-based survival for 1000 times. Then 11 lncRNAs with a frequency more than 600 were selected for prediction of the prognosis of BRCA. Using the analysis of multivariate Cox regression, we established a signature risk-score formula for 11 lncRNA to identify the relationship between lncRNA signatures and overall survival. The 11 lncRNA signature was validated both in the testing and the complete set and could effectively classify the high-/low-risk group with different OS. We also verified our results in different stages. Moreover, we analyzed the connection between the 11 lncRNAs and the genes of ESR1, PGR, and Her2, of which protein products (ESR, PGR, and HER2) were used to classify the breast cancer subtypes widely. The results indicated correlations between 11 lncRNAs and the gene of PGR and ESR1. Thus, a prognostic model for 11 lncRNA expression was developed to classify the BRAC clinical patient samples, providing new avenues in understanding the potential therapeutic methods of breast cancer.     

Clinical and molecular characteristics of bladder urothelial carcinoma subtypes

Bladder urothelial carcinoma (BLCA) is a common malignancy with high heterogeneity. A reasonable molecular subtyping can facilitate biological study and personalized therapy of BLCA. In this study, unsupervised consensus clustering was used to acquire the molecular subtypes of BLCA based on messenger RNA (mRNA) and microRNA (miRNA) data. Gene signature markers and canonical signaling pathways were compared between different subtypes. The Database for Annotation, Visualization and Integrated Discovery (DAVID) was used for the functional annotation of overexpressed genes in different subtypes for Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Three molecular subtypes were identified including C1 (luminal-P53 like), C2 (luminal-other), and C3 (basal-immune-squamous). C2 was different from C1 and C3 in clinical characteristics, including younger, better prognosis, and a higher proportion of papillary, Asian, low-grade, early-stage, lymph node negative, and complete remission patients (P < 0.05). Three molecular subtypes also showed distinct mRNA and miRNA expression patterns. luminal and P53-like markers were highly expressed in subtype C1, luminal markers were highly expressed in subtype C2, and basal, EMT/claudin-low, immune and squamous-differentiation markers were highly expressed in subtype C3. In addition, highly expressed genes in C1 were involved in extracellular signal transduction and cell-cell interaction, highly expressed genes in C2 were associated with oxygen transport, energy and steroid metabolism, and highly expressed genes in C3 were related with inflammatory, immune, cytokine, and signal transduction. BLCA in different molecular subtypes showed different clinical and molecular characteristics and personalized therapy needed to be adopted according to the molecular subtypes.     

Prognostic value of the tumor-specific ceRNA network in epithelial ovarian cancer

Epithelial ovarian cancer is one of the leading causes of cancer-related death worldwide. Growing evidence indicates that multiple complex altered pathways play important regulatory roles in the development and progression of a variety of cancers, including epithelial ovarian cancer. However, the underlying mechanisms remain unclear. First, we identified differentially expressed messenger RNAs (mRNAs), long noncoding RNAs (lncRNAs), and microRNAs (miRNAs) in epithelial ovarian cancer by comparing the expression profiles between epithelial ovarian cancer samples and normal tissue samples in different GEO datasets. Then, GO- and KEGG-pathway-enrichment analyses were applied to investigate the primary functions of the overlapped differentially expressed mRNAs. Moreover, the primary enriched genes were used to construct the signal-network with Cytoscape software. In addition, we integrated the relationship among lncRNAs-miRNAs-mRNAs to create a competing endogenous RNA network. Finally, mRNAs that were associated with patient prognosis in epithelial ovarian cancer were selected using univariate Cox regression analysis. A total of 2,225 mRNAs, 336 lncRNAs, and 14 miRNAs were shown to be differentially expressed in epithelial ovarian cancer compared with normal tissues. The dysregulated mRNAs were primarily enriched in cell division and signal transduction, according to Gene Ontology, whereas, according to KEGG, they were primarily enriched in metabolic pathways and pathways in cancer. A total of 10 mRNAs were associated with patient prognosis in ovarian cancer. This study identifies a novel lncRNA–miRNA–mRNA network, which may suggest potential molecular mechanisms underlying the development of epithelial ovarian cancer, providing new insights for survival prediction and interventional strategies for epithelial ovarian cancer.     

Genome-wide screening and cohorts validation identifying novel lncRNAs as prognostic biomarkers for clear cell renal cell carcinoma

Clear cell renal cell carcinoma (ccRCC) is the main subtype of renal cell carcinoma with varied prognosis. We aimed to identify and assess the possible prognostic long noncoding RNA (lncRNA) biomarkers. LncRNAs expression data and corresponding clinical information of 619 ccRCC patients were downloaded from The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC) databases. Differentially expressed genes analysis, univariate Cox regression, the least absolute shrinkage and selection operator Cox regression model were utilized to identify hub lncRNAs. Multivariate Cox regression was used to establish the risk model. Statistical analysis was performed using R 3.5.3. The expression value of five lncRNAs and the risk-score levels were significantly associated with a survival prognosis of ccRCC patients (all P < .001). In the TCGA validation cohort, the area under the curve (AUC) for the integrated nomogram was 0.905 and 0.91 for 3-, 5-year prediction separately. The AUC reached up to 0.757 in an independent ICGC cohort. Besides, the calibration plots also illustrated well curve-fitting between observation values and predictive values. Weighted gene co-expression network analysis and subsequent pathway analysis revealed that the PI3K-Akt-mTOR and hypoxia-inducible factor signaling crosstalk might function as the most essential mechanisms related to the five-lncRNAs signature. Our study suggested that lncRNA AC009654.1, AC092490.2, LINC00524, LINC01234, and LINC01885 were significantly associated with ccRCC prognosis. The prognostic model based on this five lncRNA may predict the overall survival of ccRCC.     

基于TCGA和GEO数据库的胃癌lncRNA筛选与ceRNA网络构建

胃癌（gastric cancer,GC）是我国最常见的恶性肿瘤之一,严重危害人类健康。胃癌发病机制复杂,缺乏特异性预后生物标志物。长链非编码RNA（long non?coding RNA,lncRNA）可作为竞争性内源RNA（competing endogenous RNA,ceRNA）,影响microRNA（miRNA）与mRNA的结合,从而影响胃癌的发生、发展。基于TCGA和GEO数据库的转录组数据,筛选GC中差异表达的lncRNAs,并构建基于6条lncRNAs（HAGLROS、TMEM92?AS1、LINC01745、HOXC?AS3、SEMA3B?AS1、FEZF1?AS1）的lncRNA?miRNA?mRNA网络。网络核心基因的KEGG/GO富集和蛋白质互作分析结果显示,lncRNA可能通过miRNA海绵吸附作用,调控胃癌的发生、发展与转移。AC011352.1、AC087636.1、AC093627.1、GAS1RR与胃癌患者的预后相关性具有统计学意义（P<0.05）,并可能成为胃癌患者潜在的预后生物标志物。     

Identification of lncRNA–miRNA–mRNA regulatory network associated with epithelial ovarian cancer cisplatin-resistant

To construct a long noncoding RNA (lncRNA)–microRNA (miRNA)–messenger RNA (mRNA) regulatory network related to epithelial ovarian cancer (EOC) cisplatin-resistant, differentially expressed genes (DEGs), differentially expressed lncRNAs (DELs), and differentially expressed miRNAs (DEMs) between MDAH and TOV-112D cells lines were identified. Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were conducted to analyze the biological functions of DEGs. Downstream mRNAs or upstream lncRNAs for miRNAs were analyzed at miRTarBase 7.0 or DIANA-LncBase V2, respectively. A total of 485 significant DEGs, 85 DELs, and 5 DEMs were identified. Protein–protein interaction (PPI) network of DEGs contrains 81 nodes and 141 edges was constructed, and 25 hub genes related to EOC cisplatin-resistant were identified. Subsequently, a lncRNA–miRNA–mRNA regulatory network contains 4 lncRNAs, 4 miRNAs, and 35 mRNAs was established. Taken together, our study provided evidence concerning the alteration genes involved in EOC cisplatin-resistant, which will help to unravel the mechanisms underlying drug resistant.     

Cancer Specific Long Noncoding RNAs Show Differential Expression Patterns and Competing Endogenous RNA Potential in Hepatocellular Carcinoma

Long noncoding RNAs (lncRNAs) regulate gene expression by acting with microRNAs (miRNAs). However, the roles of cancer specific lncRNA and its related competitive endogenous RNAs (ceRNA) network in hepatocellular cell carcinoma (HCC) are not fully understood. The lncRNA profiles in 372 HCC patients, including 372 tumor and 48 adjacent non-tumor liver tissues, from The Cancer Genome Atlas (TCGA) and NCBI GEO omnibus (GSE65485) were analyzed. Cancer specific lncRNAs (or HCC related lncRNAs) were identified and correlated with clinical features. Based on bioinformatics generated from miRcode, starBase, and miRTarBase, we constructed an lncRNA-miRNA-mRNA network (ceRNA network) in HCC. We found 177 cancer specific lncRNAs in HCC (fold change ≥ 1.5, P < 0.01), 41 of them were also discriminatively expressed with gender, race, tumor grade, AJCC tumor stage, and AJCC TNM staging system. Six lncRNAs (CECR7, LINC00346, MAPKAPK5-AS1, LOC338651, FLJ90757, and LOC283663) were found to be significantly associated with overall survival (OS, log-rank P < 0.05). Collectively, our results showed the lncRNA expression patterns and a complex ceRNA network in HCC, and identified a complex cancer specific ceRNA network, which includes 14 lncRNAs and 17 miRNAs in HCC.     

Prognostic significance of immune landscape in tumour microenvironment of endometrial cancer

Tumour microenvironment (TME) is crucial to tumorigenesis. This study aimed to uncover the differences in immune phenotypes of TME in endometrial cancer (EC) using Uterine Corpus Endometrial Carcinoma (UCEC) cohort and explore the prognostic significance. We employed GVSA enrichment analysis to cluster The Cancer Genome Atlas (TCGA) EC samples into immune signature cluster modelling, evaluated immune cell profiling in UCEC cohort (n = 538) and defined four immune subtypes of EC. Next, we analysed the correlation between immune subtypes and clinical data including patient prognosis. Furthermore, we analysed the expression of immunomodulators and DNA methylation modification. The profiles of immune infiltration in TCGA UCEC cohort showed significant difference among four immune subtypes of EC. Among each immune subtype, natural killer T cells (NKT), dendritic cells (DCs) and CD8⁺T cells were significantly associated with EC patients survival. Each immune subtype exhibited specific molecular classification, immune cell characterization and immunomodulators expression. Moreover, the expression immunomodulators were significantly related to DNA methylation level. In conclusion, the identification of immune subtypes in EC tissues could reveal unique immune microenvironments in EC and predict the prognosis of EC patients.     

Identification of metabolism-associated molecular subtype in ovarian cancer

Ovarian cancer (OC) is the most lethal gynaecological cancer with genomic complexity and extensive heterogeneity. This study aimed to characterize the molecular features of OC based on the gene expression profile of 2752 previously characterized metabolism-relevant genes and provide new strategies to improve the clinical status of patients with OC. Finally, three molecular subtypes (C1, C2 and C3) were identified. The C2 subtype displayed the worst prognosis, upregulated immune-cell infiltration status and expression level of immune checkpoint genes, lower burden of copy number gains and losses and suboptimal response to targeted drug bevacizumab. The C1 subtype showed downregulated immune-cell infiltration status and expression level of immune checkpoint genes, the lowest incidence of BRCA mutation and optimal response to targeted drug bevacizumab. The C3 subtype had an intermediate immune status, the highest incidence of BRCA mutation and a secondary optimal response to bevacizumab. Gene signatures of C1 and C2 subtypes with an opposite expression level were mainly enriched in proteolysis and immune-related biological process. The C3 subtype was mainly enriched in the T cell-related biological process. The prognostic and immune status of subtypes were validated in the Gene Expression Omnibus (GEO) dataset, which was predicted with a 45-gene classifier. These findings might improve the understanding of the diversity and therapeutic strategies for OC.     

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.