The Heterogeneous Cellular States of Glioblastoma Stem Cells Revealed by Single Cell Analysis.

Glioblastoma stem cells (GSCs) contributed to the progression, treatment resistance, and relapse of glioblastoma (GBM). However, current researches on GSCs were performed usually outside the human tumor microenvironment, ignoring the importance of the cellular states of primary GSCs. In this study, we leveraged single-cell transcriptome sequencing data of 6 independent GBM cohorts from public databases, and combined lineage and stemness features to identify primary GSCs. We dissected the cell states of GSCs and correlated them with the clinical outcomes of patients. As a result, we constructed a cellular hierarchy where GSCs resided at the center. In addition, we identified and characterized 2 different and recurrent GSCs subpopulations: proliferative GSCs (pGSCs) and quiescent GSCs (qGSCs). The pGSCs showed high cell cycle activity, indicating rapid cell division, while qGSCs showed a quiescent state. Then we traced the processes of tumor development by pseudo-time analysis and tumor phylogeny, and found that GSCs accumulated throughout the whole tumor development period. During the process, pGSCs mainly contributed to the early stage and qGSCs were enriched in the later stage. Finally, we constructed an 8-gene prognostic signature reflecting pGSCs activity and found that patients whose tumors were enriched for the pGSC signature had poor clinical outcomes. Our study highlights the primary GSCs heterogeneity and its correlation to tumor development and clinical outcomes, providing the potential targets for GBM treatment.

Single-Cell Transcriptomic Analysis Reveals the Crosstalk Propensity Between the Tumor Intermediate State and the CD8+ T Exhausted State to be Associated with Clinical Benefits in Melanoma.

Heterogeneous crosstalk between tumor cells and CD8+ T cells leads to substantial variation in clinical benefits from immunotherapy in melanoma. Due to spatial distribution and functional state heterogeneity, it is still unknown whether there is a crosstalk propensity between tumor cells and CD8+ T cells in melanoma, and how this crosstalk propensity affects the clinical outcome of patients. Using public single-cell transcriptome data, extensive heterogeneous functional states and ligand-receptor interactions of tumor cells and CD8+ T cells were revealed in melanoma. Furthermore, based on the association between cell-cell communication intensity and cell state activity in a single cell, we identified a crosstalk propensity between the tumor intermediate state and the CD8+ T exhausted state. This crosstalk propensity was further verified by pseudo-spatial proximity, spatial co-location, and the intra/intercellular signal transduction network. At the sample level, the tumor intermediate state and the CD8+ T exhausted state synergistically indicated better prognosis and both reduced in immunotherapy-resistant samples. The risk groups defined based on these two cell states could comprehensively reflect tumor genomic mutations and anti-tumor immunity information. The low-risk group had a higher BRAF mutation fraction as well as stronger antitumor immune response. Our findings highlighted the crosstalk propensity between the tumor intermediate state and the CD8+ T exhausted state, which may serve as a reference to guide the development of diagnostic biomarkers for risk stratification and therapeutic targets for new therapeutic strategies.

Revealing the contribution of somatic gene mutations to shaping tumor immune microenvironment.

Interaction between tumor cells and immune cells determined highly heterogeneous microenvironments across patients, leading to substantial variation in clinical benefits from immunotherapy. Somatic gene mutations were found not only to elicit adaptive immunity but also to influence the composition of tumor immune microenvironment and various processes of antitumor immunity. However, due to an incomplete view of associations between gene mutations and immunophenotypes, how tumor cells shape the immune microenvironment and further determine the clinical benefit of immunotherapy is still unclear. To address this, we proposed a computational approach, inference of mutation effect on immunophenotype by integrated gene set enrichment analysis (MEIGSEA), for tracing back the genomic factor responsible for differences in immunophenotypes. MEIGSEA was demonstrated to accurately identify the previous confirmed immune-associated gene mutations, and systematic evaluation in simulation data further supported its performance. We used MEIGSEA to investigate the influence of driver gene mutations on the infiltration of 22 immune cell types across 19 cancers from The Cancer Genome Atlas. The top associated gene mutations with infiltration of CD8 T cells, such as CASP8, KRAS and EGFR, also showed extensive impact on other immune components; meanwhile, immune effector cells shared critical gene mutations that collaboratively contribute to shaping distinct tumor immune microenvironment. Furthermore, we highlighted the predictive capacity of gene mutations that are positively associated with CD8 T cells for the clinical benefit of immunotherapy. Taken together, we present a computational framework to help illustrate the potential of somatic gene mutations in shaping the tumor immune microenvironment.

Applicability of Anticancer Drugs for the Triple-Negative Breast Cancer Based on Homologous Recombination Repair Deficiency.

Triple-negative breast cancer (TNBC) is a highly aggressive disease with historically poor outcomes, primarily due to the lack of effective targeted therapies. Here, we established a drug sensitivity prediction model based on the homologous recombination deficiency (HRD) using 83 TNBC patients from TCGA. Through analyzing the effect of HRD status on response efficacy of anticancer drugs and elucidating its related mechanisms of action, we found rucaparib (PARP inhibitor) and doxorubicin (anthracycline) sensitive in HR-deficient patients, while paclitaxel sensitive in the HR-proficient. Further, we identified a HRD signature based on gene expression data and constructed a transcriptomic HRD score, for analyzing the functional association between anticancer drug perturbation and HRD. The results revealed that CHIR99021 (GSK3 inhibitor) and doxorubicin have similar expression perturbation patterns with HRD, and talazoparib (PARP inhibitor) could kill tumor cells by reversing the HRD activity. Genomic characteristics indicated that doxorubicin inhibited tumor cells growth by hindering the process of DNA damage repair, while the resistance of cisplatin was related to the activation of angiogenesis and epithelial-mesenchymal transition. The negative correlation of HRD signature score could interpret the association of doxorubicin pIC50 with worse chemotherapy response and shorter survival of TNBC patients. In summary, these findings explain the applicability of anticancer drugs in TNBC and underscore the importance of HRD in promoting personalized treatment development.

Analysis of Mutations and Dysregulated Pathways Unravels Carcinogenic Effect and Clinical Actionability of Mutational Processes.

Somatic mutations accumulate over time in cancer cells as a consequence of mutational processes. However, the role of mutational processes in carcinogenesis remains poorly understood. Here, we infer the causal relationship between mutational processes and somatic mutations in 5,828 samples spanning 34 cancer subtypes. We found most mutational processes cause abundant recurrent mutations in cancer genes, while exceptionally ultraviolet exposure and altered activity of the error-prone polymerase bring a large number of recurrent non-driver mutations. Furthermore, some mutations are specifically induced by a certain mutational process, such as IDH1 p.R132H which is mainly caused by spontaneous deamination of 5-methylcytosine. At the pathway level, clock-like mutational processes extensively trigger mutations to dysregulate cancer signal transduction pathways. In addition, APOBEC mutational process destroys DNA double-strand break repair pathway, and bladder cancer patients with high APOBEC activity, though with homologous recombination proficient, show a significantly longer overall survival with platinum regimens. These findings help to understand how mutational processes act on the genome to promote carcinogenesis, and further, presents novel insights for cancer prevention and treatment, as our results showing, APOBEC mutagenesis and HRD synergistically contributed to the clinical benefits of platinum-based treatment.

Dynamic regulatory networks of T cell trajectory dissect transcriptional control of T cell state transition.

T cells exhibit heterogeneous functional states, which correlate with responsiveness to immune checkpoint blockade and prognosis of tumor patients. However, the molecular regulatory mechanisms underlying the dynamic process of T cell state transition remain largely unknown. Based on single-cell transcriptome data of T cells in non-small cell lung cancer, we combined cell states and pseudo-times to propose a pipeline to construct dynamic regulatory networks for dissecting the process of T cell dysfunction. Candidate regulators at different stages were revealed in the process of tumor-infiltrating T cell dysfunction. Through comparing dynamic networks across the T cell state transition, we revealed frequent regulatory interaction rewiring and further refined critical regulators mediating each state transition. Several known regulators were identified, including TCF7, EOMES, ID2, and TOX. Notably, one of the critical regulators, TSC22D3, was frequently identified in the state transitions from the intermediate state to the pre-dysfunction and dysfunction state, exerting diverse roles in each state transition by regulatory interaction rewiring. Moreover, higher expression of TSC22D3 was associated with the clinical outcome of tumor patients. Our study embedded transcription factors (TFs) within the temporal dynamic networks, providing a comprehensive view of dynamic regulatory mechanisms controlling the process of T cell state transition.

Single-Cell Transcriptomic Analysis Reveals a Tumor-Reactive T Cell Signature Associated With Clinical Outcome and Immunotherapy Response In Melanoma.

The infiltration of tumor-reactive T cells in the tumor site is associated with better survival and immunotherapy response. However, tumor-reactive T cells were often represented by the infiltration of total CD8+ T cells, which was confounded by the presence of bystander T cells. To identify tumor-reactive T cells at the cancer lesion, we performed integration analyses of three scRNA-seq data sets of T cells in melanoma. Extensive heterogeneous functional states of T cells were revealed in the tumor microenvironment. Among these states, we identified a subset of tumor-reactive T cells which specifically expressed tumor-reactive markers and T cell activation signature, and were strongly enriched for larger T cell receptor (TCR) clones. We further identified and validated a tumor-reactive T cell signature (TRS) to evaluate the tumor reactivity of T cells in tumor patients. Patients with high TRS scores have strong immune activity and high mutation burden in the TCGA-SKCM cohort. We also demonstrated a significant association of the TRS with the clinical outcomes of melanoma patients, with higher TRS scores representing better survival, which was validated in four external independent cohorts. Furthermore, the TRS scores exhibited greater performance on prognosis prediction than infiltration levels of CD8+ T cells and previously published prognosis-related signatures. Finally, we observed the capability of TRS to predict immunotherapy response in melanoma. Together, based on integrated analysis of single-cell RNA-sequencing, we developed and validated a tumor-reactive-related signature that demonstrated significant association with clinical outcomes and response to immunotherapy.

Combined homologous recombination repair deficiency and immune activation analysis for predicting intensified responses of anthracycline, cyclophosphamide and taxane chemotherapy in triple-negative breast cancer.

BACKGROUND: Triple-negative breast cancer (TNBC) is a clinically aggressive disease with abundant variants that cause homologous recombination repair deficiency (HRD). Whether TNBC patients with HRD are sensitive to anthracycline, cyclophosphamide and taxane (ACT), and whether the combination of HRD and tumour immunity can improve the recognition of ACT responders are still unknown. METHODS: Data from 83 TNBC patients in The Cancer Genome Atlas (TCGA) was used as a discovery cohort to analyse the association between HRD and ACT chemotherapy benefits. The combined effects of HRD and immune activation on ACT chemotherapy were explored at both the genome and the transcriptome levels. Independent cohorts from the Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) and Gene Expression Omnibus (GEO) were adopted to validate our findings. RESULTS: HRD was associated with a longer ACT chemotherapy failure-free interval (FFI) with a hazard ratio of 0.16 (P = 0.004) and improved patient prognosis (P = 0.0063). By analysing both HRD status and ACT response, we identified patients with a distinct TNBC subtype (ACT-S&HR-P) that showed higher tumour lymphocyte infiltration, IFN-γ activity and NK cell levels. Patients with ACT-S&HR-P had significantly elevated immune inhibitor levels and presented immune activation associated with the increased activities of both innate immune cells and adaptive immune cells, which suggested treatment with immune checkpoint blockade as an option for this subtype. Our analysis revealed that the combination of HRD and immune activation enhanced the efficiency of identifying responders to ACT chemotherapy (AUC = 0.91, P = 1.06e-04) and synergistically contributed to the clinical benefits of TNBC patients. A transcriptional HRD signature of ACT response-related prognostic factors was identified and independently validated to be significantly associated with improved survival in the GEO cohort (P = 0.0038) and the METABRIC dataset (P < 0.0001). CONCLUSIONS: These findings highlight that HR deficiency prolongs FFI and predicts intensified responses in TNBC patients by combining HRD and immune activation, which provides a molecular basis for identifying ACT responders.

Revealing the subtyping of non-small cell lung cancer based on genomic evolutionary patterns by multi-region sequencing.

Accurately classifying patients with non-small cell lung cancer (NSCLC) from the perspective of tumor evolution has not been systematically studied to date. Here, we reconstructed phylogenetic relationships of somatic mutations in 100 early NSCLC patients (327 lesions) through reanalyzing the TRACERx data. Based on the genomic evolutionary patterns presented on the phylogenetic trees, we grouped NSCLC patients into three evolutionary subtypes. The phylogenetic trees among three subtypes exhibited distinct branching structures, with one subtype representing branched evolution and another reflecting the early accumulation of genomic variation. However, in the evolutionary pattern of the third subtype, some mutations experienced selective sweeps and were gradually replaced by multiple newly formed subclonal populations. The subtype patients with poor prognosis had higher intra-tumor heterogeneity and subclonal diversity. We combined genomic heterogeneity with clinical phenotypes analysis and found that subclonal expansion results in the progression and deterioration of the tumor. The molecular mechanisms of subtype-specific Early Driver Feature (EDF) genes differed across the evolutionary subtypes, reflecting the characteristics of the subtype itself. In summary, our study provided new insights on the stratification of NSCLC patients based on genomic evolution that can be valuable for us to understand the development of pulmonary tumor profoundly.

Prioritizing Gene Cascading Paths to Model Colorectal Cancer Through Engineered Organoids.

Engineered organoids by sequential introduction of key mutations could help modeling the dynamic cancer progression. However, it remains difficult to determine gene paths which were sufficient to capture cancer behaviors and to broadly explain cancer mechanisms. Here, as a case study of colorectal cancer (CRC), functional and dynamic characterizations of five types of engineered organoids with different mutation combinations of five driver genes (APC, SMAD4, KRAS, TP53, and PIK3CA) showed that sequential introductions of all five driver mutations could induce enhanced activation of more hallmark signatures, tending to cancer. Comparative analysis of engineered organoids and corresponding CRC tissues revealed sequential introduction of key mutations could continually shorten the biological distance from engineered organoids to CRC tissues. Nevertheless, there still existed substantial biological gaps between the engineered organoid even with five key mutations and CRC samples. Thus, we proposed an integrative strategy to prioritize gene cascading paths for shrinking biological gaps between engineered organoids and CRC tissues. Our results not only recapitulated the well-known adenoma-carcinoma sequence model (e.g., AKST-organoid with driver mutations in APC, KRAS, SMAD4, and TP53), but also provided potential paths for delineating alternative pathogenesis underlying CRC populations (e.g., A-organoid with APC mutation). Our strategy also can be applied to both organoids with more mutations and other cancers, which can improve and innovate mechanism across cancer patients for drug design and cancer therapy.

Identifying Cancer Driver lncRNAs Bridged by Functional Effectors through Integrating Multi-omics Data in Human Cancers.

The accumulation of somatic driver mutations in the human genome enables cells to gradually acquire a growth advantage and contributes to tumor development. Great efforts on protein-coding cancer drivers have yielded fruitful discoveries and clinical applications. However, investigations on cancer drivers in non-coding regions, especially long non-coding RNAs (lncRNAs), are extremely scarce due to the limitation of functional understanding. Thus, to identify driver lncRNAs integrating multi-omics data in human cancers, we proposed a computational framework, DriverLncNet, which dissected the functional impact of somatic copy number alteration (CNA) of lncRNAs on regulatory networks and captured key functional effectors in dys-regulatory networks. Applying it to 5 cancer types from The Cancer Genome Atlas (TCGA), we portrayed the landscape of 117 driver lncRNAs and revealed their associated cancer hallmarks through their functional effectors. Moreover, lncRNA RP11-571M6.8 was detected to be highly associated with immunotherapeutic targets (PD-1, PD-L1, and CTLA-4) and regulatory T cell infiltration level and their markers (IL2RA and FCGR2B) in glioblastoma multiforme, highlighting its immunosuppressive function. Meanwhile, a high expression of RP11-1020A11.1 in bladder carcinoma was predictive of poor survival independent of clinical characteristics, and CTD-2256P15.2 in lung adenocarcinoma responded to the sensitivity of methyl ethyl ketone (MEK) inhibitors. In summary, this study provided a framework to decipher the mechanisms of tumorigenesis from driver lncRNA level, established a new landscape of driver lncRNAs in human cancers, and offered potential clinical implications for precision oncology.

Discovering Rare Genes Contributing to Cancer Stemness and Invasive Potential by GBM Single-Cell Transcriptional Analysis.

Single-cell RNA sequencing presents the sophisticated delineation of cell transcriptomes in many cancer types and highlights the tumor heterogeneity at higher resolution, which provides a new chance to explore the molecular mechanism in a minority of cells. In this study, we utilized publicly available single-cell RNA-seq data to discover and comprehensively dissect rare genes existing in few glioblastoma (GBM) cells. Moreover, we designed a framework to systematically identify 51 rare protein-coding genes (PCGs) and 47 rare long non-coding RNAs (lncRNAs) in GBM. Patients with high expression levels of rare genes like CYB5R2 and TPPP3 had worse overall survival and disease-free survival, implying their potential implication in GBM progression and prognosis. We found that these rare genes tended to be specifically expressed in GBM cancer stem cells, which emphasized their ability to characterize stem-like cancer cells and implied their contribution to GBM growth. Furthermore, rare genes were enriched in a 17-cell subset, which was located in an individual branch of the pseudotime trajectory of cancer progression and exhibited high cell cycle activity and invasive potential. Our study captures the rare genes highly expressed in few cells, deepens our understanding of special states during GBM tumorigenesis and progression such as cancer stemness and invasion, and proposes potential targets for cancer therapy.

CancerSEA: a cancer single-cell state atlas.

High functional heterogeneity of cancer cells poses a major challenge for cancer research. Single-cell sequencing technology provides an unprecedented opportunity to decipher diverse functional states of cancer cells at single-cell resolution, and cancer scRNA-seq datasets have been largely accumulated. This emphasizes the urgent need to build a dedicated resource to decode the functional states of cancer single cells. Here, we developed CancerSEA (http://biocc.hrbmu.edu.cn/CancerSEA/ or http://202.97.205.69/CancerSEA/), the first dedicated database that aims to comprehensively explore distinct functional states of cancer cells at the single-cell level. CancerSEA portrays a cancer single-cell functional state atlas, involving 14 functional states (including stemness, invasion, metastasis, proliferation, EMT, angiogenesis, apoptosis, cell cycle, differentiation, DNA damage, DNA repair, hypoxia, inflammation and quiescence) of 41 900 cancer single cells from 25 cancer types. It allows querying which functional states are associated with the gene (or gene list) of interest in different cancers. CancerSEA also provides functional state-associated PCG/lncRNA repertoires across all cancers, in specific cancers, and in individual cancer single-cell datasets. In summary, CancerSEA provides a user-friendly interface for comprehensively searching, browsing, visualizing and downloading functional state activity profiles of tens of thousands of cancer single cells and the corresponding PCGs/lncRNAs expression profiles.

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.

A pan-cancer atlas of cancer hallmark-associated candidate driver lncRNAs.

Substantial cancer genome sequencing efforts have discovered many important driver genes contributing to tumorigenesis. However, very little is known about the genetic alterations of long non-coding RNAs (lncRNAs) in cancer. Thus, there is a need for systematic surveys of driver lncRNAs. Through integrative analysis of 5918 tumors across 11 cancer types, we revealed that lncRNAs have undergone dramatic genomic alterations, many of which are mutually exclusive with well-known cancer genes. Using the hypothesis of functional redundancy of mutual exclusivity, we developed a computational framework to identify driver lncRNAs associated with different cancer hallmarks. Applying it to pan-cancer data, we identified 378 candidate driver lncRNAs whose genomic features highly resemble the known cancer driver genes (e.g. high conservation and early replication). We further validated the candidate driver lncRNAs involved in 'Tissue Invasion and Metastasis' in lung adenocarcinoma and breast cancer, and also highlighted their potential roles in improving clinical outcomes. In summary, we have generated a comprehensive landscape of cancer candidate driver lncRNAs that could act as a starting point for future functional explorations, as well as the identification of biomarkers and lncRNA-based target therapy.

Capturing functional long non-coding RNAs through integrating large-scale causal relations from gene perturbation experiments.

Characterizing functions of long noncoding RNAs (lncRNAs) remains a major challenge, mostly due to the lack of lncRNA-involved regulatory relationships. A wide array of genome-wide expression profiles generated by gene perturbation have been widely used to capture causal links between perturbed genes and response genes. Through annotating >600 gene perturbation profiles, over 354,000 causal relationships between perturbed genes and lncRNAs were identified. This large-scale resource of causal relations inspired us to develop a novel computational approach LnCAR for inferring lncRNAs' functions, which showed a higher accuracy than the co-expression based approach. By application of LnCAR to the cancer hallmark processes, we identified 38 lncRNAs involved in distinct carcinogenic processes. The "activating invasion & metastasis" related lncRNAs were strongly associated with metastatic progression in various cancer types and could act as a predictor of cancer metastasis. Meanwhile, the "evading immune destruction" related lncRNAs showed significant associations with immune infiltration of various immune cells and, importantly, can predict response to anti-PD-1 immunotherapy, suggesting their potential roles as biomarkers for immune therapy. Taken together, our approach provides a novel way to systematically reveal functions of lncRNAs, which will be helpful for further experimental exploration and clinical translational research of lncRNAs.

TIP: A Web Server for Resolving Tumor Immunophenotype Profiling.

: Systematically tracking the tumor immunophenotype is required to understand the mechanisms of cancer immunity and improve clinical benefit of cancer immunotherapy. However, progress in current research is hindered by the lack of comprehensive immune activity resources and easy-to-use tools for biologists, clinicians, and researchers to conveniently evaluate immune activity during the "cancer-immunity cycle." We developed a user-friendly one-stop shop web tool called TIP to comprehensively resolve tumor immunophenotype. TIP has the capability to rapidly analyze and intuitively visualize the activity of anticancer immunity and the extent of tumor-infiltrating immune cells across the seven-step cancer-immunity cycle. Also, we precalculated the pan-cancer immunophenotype for 11,373 samples from 33 The Cancer Genome Atlas human cancers that allow users to obtain and compare immunophenotype of pan-cancer samples. We expect TIP to be useful in a large number of emerging cancer immunity studies and development of effective immunotherapy biomarkers. TIP is freely available for use at http://biocc.hrbmu.edu.cn/TIP/. SIGNIFICANCE: TIP is a one-stop shop platform that can help biologists, clinicians, and researchers conveniently evaluate anticancer immune activity with their own gene expression data.See related commentary by Hirano, p. 6536.

Optimizing prognosis-related key miRNA-target interactions responsible for cancer metastasis.

Increasing evidence suggests that the abnormality of microRNAs (miRNAs) and their downstream targets is frequently implicated in the pathogenesis of human cancers, however, the clinical benefit of causal miRNA-target interactions has been seldom studied. Here, we proposed a computational method to optimize prognosis-related key miRNA-target interactions by combining transcriptome and clinical data from thousands of TCGA tumors across 16 cancer types. We obtained a total of 1,956 prognosis-related key miRNA-target interactions between 112 miRNAs and 1,443 their targets. Interestingly, these key target genes are specifically involved in tumor progression-related functions, such as 'cell adhesion' and 'cell migration'. Furthermore, they are most significantly correlated with 'tissue invasion and metastasis', a hallmark of metastasis, in ten distinct types of cancer through the hallmark analysis. These results implicated that the prognosis-related key miRNA-target interactions were highly associated with cancer metastasis. Finally, we observed that the combination of these key miRNA-target interactions allowed to distinguish patients with good prognosis from those with poor prognosis both in most TCGA cancer types and independent validation sets, highlighting their roles in cancer metastasis. We provided a user-friendly database named miRNATarget (freely available at http://biocc.hrbmu.edu.cn/miRNATar/), which provides an overview of the prognosis-related key miRNA-target interactions across 16 cancer types.