Systematic investigation of the prognostic impact of clonal status of somatic mutations across multiple cancer types.

Tumors are genetically heterogeneous and many mutations are actually present in subclonal populations. The clonal status of mutations is valuable for accurate prognosis, clinical management. The aim of this study was to identify the clonal status of somatic mutations and systematically evaluate their prognostic values across various cancer types. We totally identified 227 clonal and 432 subclonal mutations contributed to prognosis and demonstrated the importance of clonal status in improving mutation-related clinical guidance. We further developed a customized multi-step approach to identify gene-specific prognostic patterns of clonal status at pan-cancer level and found some cancer-specific prognostic patterns. The 'subclonal-dependent risk' subpattern was one of the most common subpatterns, it usually accompanied by high genomic in-stability and high extent of intra-tumor heterogeneity and could be used to improve the accuracy of prognostic analysis. Our results revealed the importance of clonal status, especially subclonal mutation in clinical survival.

IDH clonal heterogeneity segregates a subgroup of non-1p/19q codeleted gliomas with unfavourable clinical outcome.

AIMS: Diffuse gliomas (DGs) are classified into three major molecular subgroups following the revised World Health Organisation (WHO) classification criteria based on their IDH mutation and 1p/19q codeletion status. However, substantial biological heterogeneity and differences in the clinical course are apparent within each subgroup, which remain to be resolved. We sought to assess the clonal status of somatic mutations and explore whether additional molecular subgroups exist within DG. METHODS: A computational framework that integrates the variant allele frequency, local copy number and tumour purity was used to infer the clonality of somatic mutations in 876 DGs from The Cancer Genome Atlas (TCGA). We performed an unsupervised cluster analysis to identify molecular subgroups and characterised their clinical and biological significance. RESULTS: DGs showed widespread genetic intratumoural heterogeneity (ITH), with nearly all driver genes harbouring subclonal mutations, even for known glioma initiating event IDH1 (17.1%). Gliomas with subclonal IDH mutation and without 1p/19q codeletion showed shorter overall and disease-specific survival, higher ITH and exhibited differences in genomic patterns, transcript levels and proliferative potential, when compared with IDH clonal mutation and no 1p/19q codeletion gliomas. We defined a refined stratification system based on the current WHO glioma molecular classification, which showed close correlations with patients' clinical outcomes. CONCLUSIONS: For the first time, we integrated the clonal status of somatic mutations into cancer genomic classification and highlighted the necessity of considering IDH clonal architectures in glioma precision stratification.

CellMarker: a manually curated resource of cell markers in human and mouse.

One of the most fundamental questions in biology is what types of cells form different tissues and organs in a functionally coordinated fashion. Larger-scale single-cell sequencing and biology experiment studies are now rapidly opening up new ways to track this question by revealing substantial cell markers for distinguishing different cell types in tissues. Here, we developed the CellMarker database (http://biocc.hrbmu.edu.cn/CellMarker/ or http://bio-bigdata.hrbmu.edu.cn/CellMarker/), aiming to provide a comprehensive and accurate resource of cell markers for various cell types in tissues of human and mouse. By manually curating over 100 000 published papers, 4124 entries including the cell marker information, tissue type, cell type, cancer information and source, were recorded. At last, 13 605 cell markers of 467 cell types in 158 human tissues/sub-tissues and 9148 cell makers of 389 cell types in 81 mouse tissues/sub-tissues were collected and deposited in CellMarker. CellMarker provides a user-friendly interface for browsing, searching and downloading markers of diverse cell types of different tissues. Furthermore, a summarized marker prevalence in each cell type is graphically and intuitively presented through a vivid statistical graph. We believe that CellMarker is a comprehensive and valuable resource for cell researches in precisely identifying and characterizing cells, especially at the single-cell level.

SEECancer: a resource for somatic events in evolution of cancer genome.

Cancer cells progressively evolve from a premalignant to a malignant state, which is driven by accumulating somatic alterations that confer normal cells a fitness advantage. Improvements in high-throughput sequencing techniques have led to an increase in construction of tumor phylogenetics and identification of somatic driver events that specifically occurred in different tumor progression stages. Here, we developed the SEECancer database (http://biocc.hrbmu.edu.cn/SEECancer), which aims to present the comprehensive cancer evolutionary stage-specific somatic events (including early-specific, late-specific, relapse-specific, metastasis-specific, drug-resistant and drug-induced genomic events) and their temporal orders. By manually curating over 10 000 published articles, 1231 evolutionary stage-specific genomic events and 5772 temporal orders involving 82 human cancers and 23 tissue origins were collected and deposited in the SEECancer database. Each entry contains the somatic event, evolutionary stage, cancer type, detection approach and relevant evidence. SEECancer provides a user-friendly interface for browsing, searching and downloading evolutionary stage-specific somatic events and temporal relationships in various cancers. With increasing attention on cancer genome evolution, the necessary information in SEECancer will facilitate understanding of cancer etiology and development of evolutionary therapeutics, and help clinicians to discover biomarkers for monitoring tumor progression.

Revealing clonality and subclonality of driver genes for clinical survival benefits in breast cancer.

PURPOSE: Genomic studies have revealed that genomic aberrations play important roles in the progression of this disease. The aim of this study was to evaluate the associations between clinical survival outcomes of the clonality and subclonality status of driver genes in breast cancer. METHODS: We performed an integrated analysis to infer the clonal status of 55 driver genes in breast cancer data from TCGA. We used the chi-squared test to assess the relations between clonality of driver gene mutations and clinicopathological factors. The Kaplan-Meier method was performed for the visualization and the differences between survival curves were calculated by log-rank test. Univariate and multivariate Cox proportional hazards regression models were used to adjust for clinicopathological factors. RESULTS: We identified a high proportion of clonal mutations in these driver genes. Among them, there were 17 genes showing significant associations between their clonality and multiple clinicopathologic factors. Performing survival analysis on BRCA patients with clonal or subclonal driver gene mutations, we found that clonal ERBB2, FOXA1, and KMT2C mutations and subclonal GATA3 and RB1 mutations predicted shorter overall survival compared with those with wild type. Furthermore, clonal ERBB2 and FOXA1 mutations and subclonal GATA3 and RB1 mutations independently predicted for shorter overall survival after adjusting for clinicopathological factors. By longitudinal analysis, the clonality of ERBB2, FOXA1, GATA3, and RB1 significantly predicted patients' outcome within some specific BRCA tumor stages and histological subtypes. CONCLUSIONS: In summary, these clonal or subclonal mutations of driver genes have implications for diagnosis, prognosis, and treatment with BRCA patients.

Prognostic impact of a lymphocyte activation-associated gene signature in GBM based on transcriptome analysis.

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 < 0.001). This signature showed similar prognostic values in the other two datasets. Further, univariate and multivariate Cox proportional hazards regression models analysis indicated that the signature was an independent prognostic factor for GBM patients. Moreover, we determined that there were differences in lymphocyte activity between the high- and low-risk groups of GBM patients among all datasets. Furthermore, the lymphocyte activation-associated gene signature could significantly predict the survival of patients with certain features, including IDH-wildtype patients and patients undergoing radiotherapy. In addition, the signature may also improve the prognostic power of age. 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.

Identification of a Six-lncRNA Signature With Prognostic Value for Breast Cancer 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.

Sex difference of mutation clonality in diffuse glioma evolution.

BACKGROUND: Sex differences in glioma incidence and outcome have been previously reported but remain poorly understood. Many sex differences that affect the cancer risk were thought to be associated with cancer evolution. METHODS: In this study, we used an integrated framework to infer the timing and clonal status of mutations in ~600 diffuse gliomas from The Cancer Genome Atlas (TCGA) including glioblastomas (GBMs) and low-grade gliomas (LGGs), and investigated the sex difference of mutation clonality. RESULTS: We observed higher overall and subclonal mutation burden in female patients with different grades of gliomas, which could be largely explained by the mutations of the X chromosome. Some well-established drivers were identified showing sex-biased clonality, such as CDH18 and ATRX. Focusing on glioma subtypes, we further found a higher subclonal mutation burden in females than males in the majority of glioma subtypes, and observed opposite clonal tendency of several drivers between male and female patients in a specific subtype. Moreover, analysis of clinically actionable genes revealed that mutations in genes of the mitogen-activated protein kinase (MAPK) signaling pathway were more likely to be clonal in female patients with GBM, whereas mutations in genes involved in the receptor tyrosine kinase signaling pathway were more likely to be clonal in male patients with LGG. CONCLUSIONS: The patients with diffuse glioma showed sex-biased mutation clonality (eg, different subclonal mutation number and different clonal tendency of cancer genes), highlighting the need to consider sex as an important variable for improving glioma therapy and clinical care.