Search

Endothelial cells (ECs) play an important role in tumor progression. Currently, the main target of anti-angiogenic therapy is the vascular endothelial growth factor (VEGF) pathway. Some patients do benefit from anti-VEGF/VEGFR therapy; however, a large number of patients do not have response or acquire drug resistance after treatment. Moreover, anti-VEGF/VEGFR therapy may lead to nephrotoxicity and cardiovascular-related side effects due to its action on normal ECs. Therefore, it is necessary to identify targets that are specific to tumor ECs and could be applied to various cancer types. We integrated single-cell RNA sequencing data from six cancer types and constructed a multi-cancer EC atlas to decode the characteristic of tumor ECs. We found that tip-like ECs mainly exist in tumor tissues but barely exist in normal tissues. Tip-like ECs are involved in the promotion of tumor angiogenesis and inhibition on anti-tumor immune responses. Moreover, tumor cells, myeloid cells, and pericytes are the main sources of pro-angiogenic factors. High proportion of tip-like ECs is associated with poor prognosis in multiple cancer types. We also identified that prostate-specific membrane antigen (PSMA) is a specific marker for tip-like ECs in all the cancer types we studied. In summary, we demonstrate that tip-like ECs are the main differential EC subcluster between tumors and normal tissues. Tip-like ECs may promote tumor progression through promoting angiogenesis while inhibiting anti-tumor immune responses. PSMA was a specific marker for tip-like ECs, which could be used as a potential target for the diagnosis and treatment of non-prostate cancers., Competing Interests: The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (© 2022 The Authors.)

Accumulating evidence has recognized that cancer-associated fibroblasts (CAFs) are major players in the desmoplastic stroma of ovarian cancer, modulating tumor progression and therapeutic response. However, it is unclear regarding the signatures of CAFs could be utilized to predict the clinical outcomes of ovarian cancer patients. To fill in this gap, we explored the intratumoral compartment of ovarian cancer by analyzing the single-cell RNA-sequencing (scRNA-seq) datasets of ovarian carcinoma samples, and identified two distinct CAFs (tumor-promoting CAF_c1 subtype and myofibroblasts-like CAF_c2 subtype). The clinical significance of CAF subtypes was further validated in The Cancer Genomics Atlas (TCGA) database and other independent immunotherapy response datasets, and the results revealed that the patients with a higher expression of CAF_c1 signatures had a worse prognosis and showed a tendency of resistance to immunotherapy. This work uncovered the signatures of the CAF_c1 subtype that could serve as a novel prognostic indicator and predictive marker for immunotherapy., Competing Interests: The authors declare the following financial interests/personal relationships which may be considered as potential competing interests: Yixuan Huang and Xinlei Zhang are employees of Beijing Cloudna Technology Co., Ltd., and the other authors declare no competing financial interests., (© 2022 The Author(s).)

Tumor necrosis factor-α-inducible protein 8-like 2 (TIPE2) is encoded by TNFAIP8L2 and is a newly identified negative regulator of natural and acquired immunity that plays a critical function in maintaining immune homeostasis. Recently, CAR-NK immune cell therapy has been a focus of major research efforts as a novel cancer therapeutic strategy. TIPE2 is a potential checkpoint molecule for immune cell maturation and antitumor immunity that could be used as a novel NK cell-based immunotherapeutic approach. In this study, we explored the expression of TNFAIP8L2 across various tumor types and found that TNFAIP8L2 was highly expressed in most tumor types and correlated with prognosis. Survival analysis showed that TNFAIP8L2 expression was predictive of improved survival in cervical-squamous-cell-carcinoma (CESC), sarcoma (SARC) and skin-cutaneous-melanoma (SKCM). Conversely, TNFAIP8L2 expression predicted poorer survival in acute myeloid leukemia (LAML), lower-grade-glioma (LGG), kidney-renal-clear-cell-carcinoma (KIRC) and uveal-melanoma (UVM). Analysis of stemness features and immune cell infiltration indicated that TNFAIP8L2 was significantly associated with cancer stem cell index and increased macrophage and dendritic cell infiltration. Our data suggest that TNFAIP8L2 may be a novel immune checkpoint biomarker across different tumor types, particularly in LAML, LGG, KIRC and UVM, and may have further utility as a potential target for immunotherapy., Competing Interests: The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (© 2022 The Author(s).)

Centrosome and spindle pole-associated protein ( CSPP1 ) is a centrosome and microtubule-binding protein that plays a role in cell cycle-dependent cytoskeleton organization and cilia formation. Previous studies have suggested that CSPP1 plays a role in tumorigenesis; however, no pan-cancer analysis has been performed. This study systematically investigates the expression of CSPP1 and its potential clinical outcomes associated with diagnosis, prognosis, and therapy. CSPP1 is widely present in tissues and cells and its aberrant expression serves as a diagnostic biomarker for cancer. CSPP1 dysregulation is driven by multi-dimensional mechanisms involving genetic alterations, DNA methylation, and miRNAs. Phosphorylation of CSPP1 at specific sites may play a role in tumorigenesis. In addition, CSPP1 correlates with clinical features and outcomes in multiple cancers. Take brain low-grade gliomas (LGG) with a poor prognosis as an example, functional enrichment analysis implies that CSPP1 may play a role in ferroptosis and tumor microenvironment (TME), including regulating epithelial-mesenchymal transition, stromal response, and immune response. Further analysis confirms that CSPP1 dysregulates ferroptosis in LGG and other cancers, making it possible for ferroptosis-based drugs to be used in the treatment of these cancers. Importantly, CSPP1 -associated tumors are infiltrated in different TMEs, rendering immune checkpoint blockade therapy beneficial for these cancer patients. Our study is the first to demonstrate that CSPP1 is a potential diagnostic and prognostic biomarker associated with ferroptosis and TME, providing a new target for drug therapy and immunotherapy in specific cancers., Competing Interests: The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (© 2022 The Authors. Published by Elsevier B.V. on behalf of Research Network of Computational and Structural Biotechnology.)

Shc SH2-domain binding protein 1 (SHCBP1), a protein specific binding to SH2 domain of Src homolog and collagen homolog (Shc), takes part in the regulation of various signal transduction pathways, which has been reported to be associated with tumorigenesis and progression. However, the pathological mechanisms are not completely investigated. Thus, this study aimed to comprehensively elucidate the potential functions of SHCBP1 in multiple cancer types. The comprehensive analyses for SHCBP1 in various tumors, including gene expression, diagnosis, prognosis, immune-related features, genetic alteration, and function enrichment, were conducted based on multiple databases and analysis tools. SHCBP1 was upregulated in most types of cancers. The results of qRT-PCR had confirmed that SHCBP1 mRNA was significantly upregulated in lung adenocarcinoma (LUAD) and liver hepatocellular carcinoma (LIHC) cell lines. Based on the receiver operating characteristic (ROC) and survival analysis, SHCBP1 was considered as a potential diagnostic and prognostic biomarker. Furthermore, SHCBP1 expression was linked with tumor immunity and immunosuppressive microenvironment according to the correlation analysis of SHCBP1 expression with immune cells infiltration, immune checkpoint genes, and immune-related genes (MHC genes, chemokines, and chemokines receptors). Moreover, SHCBP1 expression correlated with tumor mutational burden (TMB), microsatellite instability (MSI), and neoantigens. The feature of SHCBP1 mutational landscape in pan-cancer was identified. Finally, we focused on investigating the clinical significance and the potential biological role of SHCBP1 in LUAD. Our study comprehensively uncovered that SHCBP1 could be identified as an immune-related biomarker for cancer diagnosis and prognosis, and a potential therapeutic target for tumor immunotherapy., Competing Interests: The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (© 2022 The Author(s).)

Background: Cox proportional hazard regression (CPH) model relies on the proportional hazard (PH) assumption: the hazard of variables is independent of time. CPH has been widely used to identify prognostic markers of the transcriptome. However, the comprehensive investigation on PH assumption in transcriptomic data has lacked., Results: The whole transcriptomic data of the 9,056 patients from 32 cohorts of The Cancer Genome Atlas and the 3 lung cancer cohorts from Gene Expression Omnibus were collected to construct CPH model for each gene separately for fitting the overall survival. An average of 8.5% gene CPH models violated the PH assumption in TCGA pan-cancer cohorts. In the gene interaction networks, both hub and non-hub genes in CPH models were likely to have non-proportional hazards. Violations of PH assumption for the same gene models were not consistent in 5 non-small cell lung cancer datasets (all kappa coefficients < 0.2), indicating that the non-proportionality of gene CPH models depended on the datasets. Furthermore, the introduction of log(t) or sqrt(t) time-functions into CPH improved the performance of gene models on overall survival fitting in most tumors. The time-dependent CPH changed the significance of log hazard ratio of the 31.9% gene variables., Conclusions: Our analysis resulted that non-proportional hazards should not be ignored in transcriptomic data. Introducing time interaction term ameliorated performance and interpretability of non-proportional hazards of transcriptome data in CPH., (© 2022 The Author(s).)

The SARS-CoV-2 is constantly mutating, and the new coronavirus such as Omicron has spread to many countries around the world. Anexelekto (AXL) is a transmembrane protein with biological functions such as promoting cell growth, migration, aggregation, metastasis and adhesion, and plays an important role in cancers and coronavirus disease 2019 (COVID-19). Unlike angiotensin-converting enzyme 2 (ACE2), AXL was highly expressed in respiratory system cells. In this study, we verified the AXL expression in cancer and normal tissues and found AXL expression was strongly correlated with cancer prognosis, tumor mutation burden (TMB), the microsatellite instability (MSI) in most tumor types. Immune infiltration analysis also demonstrated that there was an inextricable link between AXL expression and immune scores in cancer patients, especially in BLCA, BRCA and CESC. The NK-cells, plasmacytoid dendritic cells, myeloid dendritic cells, as one of the important components of the tumor microenvironment, were highly expressed AXL. In addition, AXL-related tumor neoantigens were identified and might provide the novel potential targets for tumor vaccines or SARS-Cov-2 vaccines research in cancer patients., Competing Interests: The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (© 2022 The Author(s).)

Cholangiocarcinomas (CCAs) are tumors that arise from the cholangiocytes. Although some genes have been shown with important roles in pathological process, interactions or cross-talks among different RNAs are important to understand the detailed molecular mechanisms in cancer development, especially discussing cross-talks among isomiRs and other RNAs. Herein, to characterize crucial genes in CCA, the protein expression profile was performed to survey potential crucial mRNAs and related non-coding RNAs (ncRNAs) in mRNA-ncRNA network, mainly including miRNAs/isomiRs and lncRNAs. Deregulated mRNAs were firstly obtained if consistent expression patterns were found at protein and mRNA levels, and related miRNAs/isomiRs were screened according to regulatory relationships. Diverse isomiRs from a given miRNA locus also contributed to interactions between the small RNAs and target mRNAs, and miRNAs were further used to survey related lncRNAs to expand the interactions. Thus, several groups of RNAs were constructed as candidate competitive endogenous RNA (ceRNA) networks. Finally, we found that RAB11FIP1:miR-101-3p:MIR3142HG may be a potential ceRNA network, and the interactions among them may be more complex due to variety of isomiRs. Simultaneously, RAB11FIP1 and miR-194-5p were also detected other related lncRNAs (FBXL19-AS1, SNHG1 and PVT1) that may be crucial in coding-non-coding RNA regulatory network. Our results show that diverse isomiRs with sequence and expression heterogeneities contribute to ceRNA regulatory network that may have crucial roles in CCA, which will expand our understanding of interactions among diverse RNAs and their contributions in cancer development., Competing Interests: The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (© 2021 The Author(s).)

Background: Gliomas are one of the most common types of primary tumors in central nervous system. Previous studies have found that macrophages actively participate in tumor growth., Methods: Weighted gene co-expression network analysis was used to identify meaningful macrophage-related gene genes for clustering. Pamr, SVM, and neural network were applied for validating clustering results. Somatic mutation and methylation were used for defining the features of identified clusters. Differentially expressed genes (DEGs) between the stratified groups after performing elastic regression and principal component analyses were used for the construction of MScores. The expression of macrophage-specific genes were evaluated in tumor microenvironment based on single cell sequencing analysis. A total of 2365 samples from 15 glioma datasets and 5842 pan-cancer samples were used for external validation of MScore., Results: Macrophages were identified to be negatively associated with the survival of glioma patients. Twenty-six macrophage-specific DEGs obtained by elastic regression and PCA were highly expressed in macrophages at single-cell level. The prognostic value of MScores in glioma was validated by the active proinflammatory and metabolic profile of infiltrating microenvironment and response to immunotherapies of samples with this signature. MScores managed to stratify patient survival probabilities in 15 external glioma datasets and pan-cancer datasets, which predicted worse survival outcome. Sequencing data and immunohistochemistry of Xiangya glioma cohort confirmed the prognostic value of MScores. A prognostic model based on MScores demonstrated high accuracy rate., Conclusion: Our findings strongly support a modulatory role of macrophages, especially M2 macrophages in glioma progression and warrants further experimental studies., Competing Interests: The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (© 2021 The Authors. Published by Elsevier B.V. on behalf of Research Network of Computational and Structural Biotechnology.)

Synthetic lethality is thought to play an important role in anticancer therapies. Herein, to understand the potential distributions and relationships between synthetic lethal interactions between genes, especially for pairs deriving from different sources, we performed an integrative analysis of genes at multiple molecular levels. Based on inter-species phylogenetic conservation of synthetic lethal interactions, gene pairs from yeast and humans were analyzed; a total of 37,588 candidate gene pairs containing 7,816 genes were collected. Of these, 49.74% of genes had 2-10 interactions, 22.93% were involved in hallmarks of cancer, and 21.61% were identified as core essential genes. Many genes were shown to have important biological roles via functional enrichment analysis, and 65 were identified as potentially crucial in the pathophysiology of cancer. Gene pairs with dysregulated expression patterns had higher prognostic values. Further screening based on mutation and expression levels showed that remaining gene pairs were mainly derived from human predicted or validated pairs, while most predicted pairs from yeast were filtered from analysis. Genes with synthetic lethality were further analyzed with their interactive microRNAs (miRNAs) at the isomiR level which have been widely studied as negatively regulatory molecules. The miRNA-mRNA interaction network revealed that many synthetic lethal genes contributed to the cell cycle (seven of 12 genes), cancer pathways (five of 12 genes), oocyte meiosis, the p53 signaling pathway, and hallmarks of cancer. Our study contributes to the understanding of synthetic lethal interactions and promotes the application of genetic interactions in further cancer precision medicine., Competing Interests: The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (© 2020 The Author(s).)

Classical microRNA (miRNA) has been so far believed as a single sequence, but it indeed contains multiple miRNA isoforms (isomiR) with various sequences and expression patterns. It is not clear whether these diverse isomiRs have potential relationships and whether they contribute to miRNA:mRNA interactions. Here, we aimed to reveal the potential evolutionary and functional relationships of multiple isomiRs based on let-7 and miR-10 gene families that are prone to clustering together on chromosomes. Multiple isomiRs within gene families showed similar functions to their canonical miRNAs, indicating selection of the predominant sequence. IsomiRs containing novel seed regions showed increased/decreased biological function depending on whether they had more/less specific target mRNAs than their annotated seed. Few gene ontology(GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were shared among the target genes of the annotated seeds and the novel seeds. Various let-7 isomiRs with novel seed regions may cause opposing drug responses despite the fact that they are generated from the same miRNA locus and have highly similar sequences. IsomiRs, especially the dominant isomiRs with shifted seeds, may disturb the coding-non-coding RNA regulatory network. These findings provide insight into the multiple isomiRs and isomiR-mediated control of gene expression in the pathogenesis of cancer., Competing Interests: The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (© 2020 The Author(s).)

Black/African-American (B/AA) breast cancer patients tend to have more aggressive tumor biology compared to White/Caucasians. In this study, a variety of breast tumor molecular expression profiles of patients derived from the two racial groupings were investigated. Breast invasive carcinoma sample data (RNASeq version 2, Reverse Phase Protein Array, mutation, and miRSeq data) from the Cancer Genome Atlas were examined. The results affirm that B/AA patients are more likely than Caucasian patients to harbor the aggressive basal-like or the poor prognosis-associated HER2-enriched molecular subtypes of breast cancer. There is also a higher incidence of the triple-negative breast cancer (TNBC) among B/AA patients than the general population, a fact reflected in the mutation patterns of genes such as PIK3CA and TP53 . Furthermore, an immortalization signature gene set, is enriched in samples from B/AA patients. Among stage III patients, TERT, DRAP1, and PQBP1, all members of the immortalization gene signature set, are among master-regulators with increased activity in B/AA patients. Master-regulators driving differences in expression profiles between the two groups include immortalization markers, senescence markers, and immune response and redox gene products. Differences in expression, between B/AA and Caucasian patients, of RB1 , hsa-let-7a , E2F1 , c-MYC , TERT , and other biomolecules appear to cooperate to enhance entry into the S-phase of the cell cycle in B/AA patients. Higher expression of miR-221 , an oncomiR that facilitates entry into the cell cycle S-phase, is regulated by c-MYC , which is expressed more in breast cancer samples from B/AA patients. Furthermore, the cell migration- and invasion-promoting miRNA, miR-135b , has increased relative expression in B/AA patients. Knock down of the immortalization marker TERT inhibited triple-negative breast cancer cell lines (MDA-MB-231 and MDA-MB-468) cell viability and decreased expression of TERT, MYC and WNT11. For those patients with available survival data, prognosis of stage II patients 50 years of age or younger at diagnosis, was distinctly poorer in B/AA patients. Also associated with this subset of B/AA patients are missense mutations in HUWE1 and PTEN expression loss. Relative to Caucasian non-responders to endocrine therapy, B/AA non-responders show suppressed expression of a signature gene set on which biological processes including signaling by interleukins , circadian clock , regulation of lipid metabolism by PPARα , FOXO-mediated transcription , and regulation of TP53 degradation are over-represented. Thus, we identify molecular expression patterns suggesting diminished response to oxidative stress, changes in regulation of tumor suppressors/facilitators, and enhanced immortalization in B/AA patients are likely important in defining the more aggressive molecular tumor phenotype reported in B/AA patients., Competing Interests: The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (© 2020 The Author(s).)

Microsatellite instability (MSI) is a genomic property of the cancers with defective DNA mismatch repair and is a useful marker for cancer diagnosis and treatment in diverse cancer types. In particular, MSI has been associated with the active immune checkpoint blockade therapy response in cancer. Most of computational methods for predicting MSI are based on DNA sequencing data and a few are based on mRNA expression data. Using the RNA-Seq pan-cancer datasets for three cancer cohorts (colon, gastric, and endometrial cancers) from The Cancer Genome Atlas (TCGA) program, we developed an algorithm (PreMSIm) for predicting MSI from the expression profiling of a 15-gene panel in cancer. We demonstrated that PreMSIm had high prediction performance in predicting MSI in most cases using both RNA-Seq and microarray gene expression datasets. Moreover, PreMSIm displayed superior or comparable performance versus other DNA or mRNA-based methods. We conclude that PreMSIm has the potential to provide an alternative approach for identifying MSI in cancer., Competing Interests: The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (© 2020 The Authors.)