Search

Background Lung adenocarcinoma (LUAD), as the most common subtype of lung cancer, is the leading cause of cancer deaths in the world. The accumulation of driver gene mutations enables cancer cells to gradually acquire growth advantage. Therefore, it is important to understand the functions and interactions of driver gene mutations in cancer progression. Methods We obtained gene mutation data and gene expression profile of 506 LUAD tumors from The Cancer Genome Atlas (TCGA). The subtypes of tumors with driver gene mutations were identified by consensus cluster analysis. Results We found 21 significantly mutually exclusive pairs consisting of 20 genes among 506 LUAD patients. Because of the increased transcriptomic heterogeneity of mutations, we identified subtypes among tumors with non‐silent mutations in driver genes. There were 494 mutually exclusive pairs found among driver gene mutations within different subtypes. Furthermore, we identified functions of mutually exclusive pairs based on the hypothesis of functional redundancy of mutual exclusivity. These mutually exclusive pairs were significantly enriched in nuclear division and humoral immune response, which played crucial roles in cancer initiation and progression. We also found 79 mutually exclusive triples among subtypes of tumors with driver gene mutations, which were key roles in cell motility and cellular chemical homeostasis. In addition, two mutually exclusive triples and one mutually exclusive triple were associated with the overall survival and disease‐specific survival of LUAD patients, respectively. Conclusions We revealed novel mutual exclusivity and generated a comprehensive functional landscape of driver gene mutations, which could offer a new perspective to understand the mechanisms of cancer development and identify potential biomarkers for LUAD therapy., We revealed novel mutual exclusivity and generated a comprehensive functional landscape of driver gene mutations, which could offer a new perspective to understand the mechanisms of cancer development and identify biomarkers for LUAD therapy.

Simple Summary The aim of our study was to evaluate the predictive performance of transcriptomic biomarkers to immune response. The study collected 22 transcriptomic biomarkers and constructed multiple benchmark datasets to evaluate their predictive performance of immune checkpoint blockade (ICB) response in pre-treatment patients with distinct ICB agents in diverse cancers. We found “Immune-checkpoint molecule” biomarkers PD-L1, PD-L2, CTLA-4 and IMPRES and the “Effector molecule” biomarker CYT showed significant associations with ICB response and clinical outcomes. These immune-checkpoint biomarkers and another immune effector IFN-gamma presented predictive ability in melanoma, urothelial cancer and clear cell renal-cell cancer. Interestingly, for anti-PD-1 therapy and anti-CTLA-4 therapy, the top-performing response biomarkers were usually mutually exclusive even though in the same biomarker category and most of biomarkers with outstanding predictive power were observed in patients with combined anti-PD-1 and anti-CTLA-4 therapy. Abstract Background: Immune checkpoint blockade (ICB) therapy has yielded successful clinical responses in treatment of a minority of patients in certain cancer types. Substantial efforts were made to establish biomarkers for predicting responsiveness to ICB. However, the systematic assessment of these ICB response biomarkers remains insufficient. Methods: We collected 22 transcriptome-based biomarkers for ICB response and constructed multiple benchmark datasets to evaluate the associations with clinical response, predictive performance, and clinical efficacy of them in pre-treatment patients with distinct ICB agents in diverse cancers. Results: Overall, “Immune-checkpoint molecule” biomarkers PD-L1, PD-L2, CTLA-4 and IMPRES and the “Effector molecule” biomarker CYT showed significant associations with ICB response and clinical outcomes. These immune-checkpoint biomarkers and another immune effector IFN-gamma presented predictive ability in melanoma, urothelial cancer (UC) and clear cell renal-cell cancer (ccRCC). In non-small cell lung cancer (NSCLC), only PD-L2 and CTLA-4 showed preferable correlation with clinical response. Under different ICB therapies, the top-performing biomarkers were usually mutually exclusive in patients with anti-PD-1 and anti-CTLA-4 therapy, and most of biomarkers presented outstanding predictive power in patients with combined anti-PD-1 and anti-CTLA-4 therapy. Conclusions: Our results show these biomarkers had different performance in predicting ICB response across distinct ICB agents in diverse cancers.

Background Glioblastoma multiforme (GBM) is a highly, malignant tumor of the primary central nervous system. Patients diagnosed with this type of tumor have a poor prognosis. Lymphocyte activation plays important roles in the development of cancers and its therapeutic treatments. Objective We sought to identify an efficient lymphocyte activation-associated gene signature that could predict the progression and prognosis of GBM. Methods We used univariate Cox proportional hazards regression and stepwise regression algorithm to develop a lymphocyte activation-associated gene signature in the training dataset (TCGA, n = 525). Then, the signature was validated in two datasets, including GSE16011 (n = 150) and GSE13041 (n = 191) using the Kaplan Meier method. Univariate and multivariate Cox proportional hazards regression models were used to adjust for clinicopathological factors. Results We identified a lymphocyte activation-associated gene signature (TCF3, IGFBP2, TYRO3 and NOD2) in the training dataset and classified the patients into high-risk and low-risk groups with significant differences in overall survival (median survival 15.33 months vs 12.57 months, HR = 1.55, 95% CI [1.28–1.87], log-rank test P Conclusions In summary, our results suggested that the lymphocyte activation-associated gene signature is a promising factor for the survival of patients, which is helpful for the prognosis of GBM patients.

Breast cancer (BRCA) is the most common cancer and a major cause of death in women. Long non-coding RNAs (lncRNAs) are emerging as key regulators and have been implicated in carcinogenesis and prognosis. In this study, we aimed to develop a lncRNA signature of BRCA patients to improve risk stratification. In the training cohort (GSE21653, n = 232), 17 lncRNAs were identified by univariate Cox proportional hazards regression, which were significantly associated with patients’ survival. The least absolute shrinkage and selection operator-penalized Cox proportional hazards regression analysis was used to identify a six-lncRNA signature. According to the median of the signature risk score, patients were divided into a high-risk group and a low-risk group with significant disease-free survival differences in the training cohort. A similar phenomenon was observed in validation cohorts (GSE42568, n = 101; GSE20711, n = 87). The six-lncRNA signature remained as independent prognostic factors after adjusting for clinical factors in these two cohorts. Furthermore, this signature significantly predicted the survival of grade III patients and estrogen receptor-positive patients. Furthermore, in another cohort (GSE19615, n = 115), the low-risk patients that were treated with tamoxifen therapy had longer disease-free survival than those who underwent no therapy. Overall, the six-lncRNA signature can be a potential prognostic tool used to predict disease-free survival of patients and to predict the benefits of tamoxifen treatment in BRCA, which will be helpful in guiding individualized treatments for BRCA patients.

Long non-coding RNAs (lncRNAs) are a major class of non-coding RNAs, and the functional deregulations of lncRNAs have been shown to be associated with the development and progression of BC. In this work, we conduct an integrative analysis on five re-annotated lncRNA expression datasets from the Gene Expression Omnibus (GEO) which included a total of 891 BC samples. We identified a five-lncRNA signature that was significantly associated with DFS in the training cohort of 327 patients. We found the five-lncRNA signature could effectively stratify patients in the training dataset into high- and low-risk groups with significantly different DFS (p = 3.29 × 10-5 , log-rank test). The five-lncRNA signature was effectively validated in four independent cohorts, and prognostic analysis results showed that the five-lncRNA signature was independent of clinical prognostic factors, such as BC subtypes and adjuvant treatments. Furthermore, GSEA suggested that the five-lncRNA signature was involved in BC metastasis-related pathways. Our findings indicate that these five lncRNAs may be implicated in BC pathogenesis, and further, these lncRNAs may potentially serve as novel candidate biomarkers for the identification of BC patients at high risk for tumor recurrence.

Our knowledge of lncRNA is very limited and discovering novel disease-related long non-coding RNA (lncRNA) has been a major research challenge in cancer studies. In this work, we developed an LncRNA Network-based Prioritization approach, named "LncNetP" based on the competing endogenous RNA (ceRNA) and disease phenotype association assumptions. Through application to 11 cancer types with 3089 common lncRNA and miRNA samples from the Cancer Genome Atlas (TCGA), our approach yielded an average area under the ROC curve (AUC) of 83.87%, with the highest AUC (95.22%) for renal cell carcinoma, by the leave-one-out cross validation strategy. Moreover, we demonstrated the excellent performance of our approach by evaluating the influencing factors including disease phenotype associations, known disease lncRNAs and the numbers of cancer types. Comparisons with previous methods further suggested the integrative importance of our approach. Taking hepatocellular carcinoma (LIHC) as a case study, we predicted four candidate lncRNA genes, RHPN1-AS1, AC007389.1, LINC01116 and BMS1P20 that may serve as novel disease risk factors for disease diagnosis and prognosis. In summary, our lncRNA prioritization strategy can efficiently identify disease-related lncRNAs and help researchers better understand the important roles of lncRNAs in human cancers.

Glioblastoma (GBM) is the most common and aggressive primary brain tumor, in which GBM stem cells (GSCs) were identified to contribute to aggressive phenotypes and poor prognosis. Yet, how GSCs progress to invasive cells remains largely unexplored. Here, we revealed the cell subpopulations with distinct functional status and the existence of cells with high invasive potential within heterogeneous primary GBM tumors. We reconstructed a branched trajectory by pseudotemporal ordering of single tumor cells, in which the root showed GSC‐like phenotype while the end displayed high invasive activity. Thus, we further determined a path along which GSCs gradually transformed to invasive cells, called the ‘stem‐to‐invasion path’. Along this path, cells showed incremental expression of GBM invasion‐associated signatures and diminishing expression of GBM stem cell markers. These findings were validated in an independent single‐cell data set of GBM. Through analyzing the molecular cascades underlying the path, we identify crucial factors controlling the attainment of invasive potential of tumor cells, including transcription factors and long noncoding RNAs. Our work provides novel insights into GBM progression, especially the attainment of invasive potential in primary tumor cells, and supports the cancer stem cell model, with valuable implications for GBM therapy., Glioblastoma (GBM) is the most common and aggressive primary brain tumor with poor prognosis. Here, through reconstructing a branched trajectory by pseudotemporal ordering of single tumor cells, we determine the ‘stem‐to‐invasion path’, along which cells show incremental invasive potential and diminishing expression of GBM stem cell markers. We further identify crucial factors controlling the attainment of invasive potential of tumor cells.

Long non-coding RNAs (lncRNAs) are a major class of non-coding RNAs, and the functional deregulations of lncRNAs have been shown to be associated with the development and progression of BC. In this work, we conduct an integrative analysis on five re-annotated lncRNA expression datasets from the Gene Expression Omnibus (GEO) which included a total of 891 BC samples. We identified a five-lncRNA signature that was significantly associated with DFS in the training cohort of 327 patients. We found the five-lncRNA signature could effectively stratify patients in the training dataset into high- and low-risk groups with significantly different DFS (p = 3.29×10(−5), log-rank test). The five-lncRNA signature was effectively validated in four independent cohorts, and prognostic analysis results showed that the five-lncRNA signature was independent of clinical prognostic factors, such as BC subtypes and adjuvant treatments. Furthermore, GSEA suggested that the five-lncRNA signature was involved in BC metastasis-related pathways. Our findings indicate that these five lncRNAs may be implicated in BC pathogenesis, and further, these lncRNAs may potentially serve as novel candidate biomarkers for the identification of BC patients at high risk for tumor recurrence.