Single-cell RNA-seq reveals the invasive trajectory and molecular cascades underlying glioblastoma progression

Turning attention to tumor-host interface and focus on the peritumoral heterogeneity of glioblastoma

Approximately 90% of glioblastoma recurrences occur in the peritumoral brain zone (PBZ), while the spatial heterogeneity of the PBZ is not well studied. In this study, two PBZ tissues and one tumor tissue sample are obtained from each patient via preoperative imaging. We assess the microenvironment and the characteristics of infiltrating immune/tumor cells using various techniques. Our data indicate there are one or more regions with higher cerebral blood flow in PBZ, which we collectively name the "higher cerebral blood flow interface" (HBI). The HBI exhibited more neovascularization than the "lower cerebral blood flow interfaces" (LBI). The HBI tend to have increased infiltration of macrophages and T lymphocytes infiltration compared with that in LBI. There are more tumor cells in the HBI than in LBI, with substantial differences in the gene expression profiles of these tumor cells. HBI may be the key area of PBZ-targeting therapy after surgical resection.

Spatial Transcriptomics in a Case of Follicular Thyroid Carcinoma Reveals Clone-Specific Dysregulation of Genes Regulating Extracellular Matrix in the Invading Front

Follicular thyroid carcinoma (FTC) is recognized by its ability to invade the tumor capsule and blood vessels, although the exact molecular signals orchestrating this phenotype remain elusive. In this study, the spatial transcriptional landscape of an FTC is detailed with comparisons between the invasive front and histologically indolent central core tumor areas. The Visium spatial gene expression platform allowed us to interrogate and visualize the whole transcriptome in 2D across formalin-fixated paraffin-embedded (FFPE) tissue sections. Four different 6 × 6 mm areas of an FTC were scrutinized, including regions with capsular and vascular invasion, capsule-near area without invasion, and a central core area of the tumor. Following successful capturing and sequencing, several expressional clusters were identified with regional variation. Most notably, invasive tumor cell clusters were significantly over-expressing genes associated with pathways interacting with the extracellular matrix (ECM) remodeling and epithelial-to-mesenchymal transition (EMT). Subsets of these genes (POSTN and DPYSL3) were additionally validated using immunohistochemistry in an independent cohort of follicular thyroid tumors showing a clear gradient pattern from the core to the periphery of the tumor. Moreover, the reconstruction of the evolutionary tree identified the invasive clones as late events in follicular thyroid tumorigenesis. To our knowledge, this is one of the first 2D global transcriptional mappings of FTC using this platform to date. Invasive FTC clones develop in a stepwise fashion and display significant dysregulation of genes associated with the ECM and EMT - thus highlighting important molecular crosstalk for further investigations.

uPAR (PLAUR) Marks Two Intra-Tumoral Subtypes of Glioblastoma: Insights from Single-Cell RNA Sequencing

Urokinase plasminogen activator receptor (uPAR) encoded by the PLAUR gene is known as a clinical marker for cell invasiveness in glioblastoma multiforme (GBM). It is additionally implicated in various processes, including angiogenesis and inflammation within the tumor microenvironment. However, there has not been a comprehensive study that depicts the overall functions and molecular cooperators of PLAUR with respect to intra-tumoral subtypes of GBM. Using single-cell RNA sequencing data from 37 GBM patients, we identified PLAUR as a marker gene for two distinct subtypes in GBM. One subtype is featured by inflammatory activities and the other subtype is marked by ECM remodeling processes. Using the whole-transcriptome data from single cells, we are able to uncover the molecular cooperators of PLAUR for both subtypes without presuming biological pathways. Two protein networks comprise the molecular context of PLAUR, with each of the two subtypes characterized by a different dominant network. We concluded that targeting PLAUR directly influences the mechanisms represented by these two protein networks, regardless of the subtype of the targeted cell.

Spatial metabolic heterogeneity of oligodendrogliomas at single-cell resolution

Oligodendrogliomas are a type of rare and incurable gliomas whose metabolic profiles have yet to be fully examined. The present study examined the spatial differences in metabolic landscapes underlying oligodendrogliomas and should provide unique insights into the metabolic characteristics of these uncommon tumors. Single-cell RNA-sequencing expression profiles from 4044 oligodendroglioma cells derived from tumors resected from four locations frontal, temporal, parietal, and frontotemporoinsular) and in which 1p/19q co-deletion and IDH1 or IDH2 mutations were confirmed were computationally analyzed through a robust workflow to elucidate relative differences in metabolic pathway activities among the different locations. Dimensionality reduction using metabolic expression profiles exhibited clustering corresponding to each location subgroup. From the 80 metabolic pathways examined, over 70 pathways had significantly different activity scores between location subgroups. Further analysis of metabolic heterogeneity suggests that mitochondrial oxidative phosphorylation accounts for considerable metabolic variation within the same locations. Steroid and fatty acid metabolism pathways were also found to be major contributors to heterogeneity. Oligodendrogliomas display distinct spatial metabolic differences in addition to intra-location metabolic heterogeneity.

Quantitative Evaluation of Stem-like Markers of Human Glioblastoma Using Single-Cell RNA Sequencing Datasets

Targeting glioblastoma (GBM) stem-like cells (GSCs) is a common interest in both the laboratory investigation and clinical treatment of GBM. Most of the currently applied GBM stem-like markers lack validation and comparison with common standards regarding their efficiency and feasibility in various targeting methods. Using single-cell RNA sequencing datasets from 37 GBM patients, we obtained a large pool of 2173 GBM stem-like marker candidates. To evaluate and select these candidates quantitatively, we characterized the efficiency of the candidate markers in targeting the GBM stem-like cells by their frequencies and significance of being the stem-like cluster markers. This was followed by further selection based on either their differential expression in GBM stem-like cells compared with normal brain cells or their relative expression level compared with other expressed genes. The cellular location of the translated protein was also considered. Different combinations of selection criteria highlight different markers for different application scenarios. By comparing the commonly used GSCs marker CD133 (PROM1) with markers selected by our method regarding their universality, significance, and abundance, we revealed the limitations of CD133 as a GBM stem-like marker. Overall, we propose BCAN, PTPRZ1, SOX4, etc. for laboratory-based assays with samples free of normal cells. For in vivo targeting applications that require high efficiency in targeting the stem-like subtype, the ability to distinguish GSCs from normal brain cells, and a high expression level, we recommend the intracellular marker TUBB3 and the surface markers PTPRS and GPR56.

A Novel Molecular Classification Method for Glioblastoma Based on Tumor Cell Differentiation Trajectories

The latest 2021 WHO classification redefines glioblastoma (GBM) as the hierarchical reporting standard by eliminating glioblastoma, IDH-mutant and only retaining the tumor entity of "glioblastoma, IDH-wild type." Knowing that subclassification of tumors based on molecular features is supposed to facilitate the therapeutic choice and increase the response rate in cancer patients, it is necessary to carry out molecular classification of the newly defined GBM. Although differentiation trajectory inference based on single-cell sequencing (scRNA-seq) data holds great promise for identifying cell heterogeneity, it has not been used in the study of GBM molecular classification. Single-cell transcriptome sequencing data from 10 GBM samples were used to identify molecular classification based on differentiation trajectories. The expressions of identified features were validated by public bulk RNA-sequencing data. Clinical feasibility of the classification system was examined in tissue samples by immunohistochemical (IHC) staining and immunofluorescence, and their clinical significance was investigated in public cohorts and clinical samples with complete clinical follow-up information. By analyzing scRNA-seq data of 10 GBM samples, four differentiation trajectories from the glioblastoma stem cell-like (GSCL) cluster were identified, based on which malignant cells were classified into five characteristic subclusters. Each cluster exhibited different potential drug sensitivities, pathways, functions, and transcriptional modules. The classification model was further examined in TCGA and CGGA datasets. According to the different abundance of five characteristic cell clusters, the patients were classified into five groups which we named Ac-G, Class-G, Neo-G, Opc-G, and Undiff-G groups. It was found that the Undiff-G group exhibited the worst overall survival (OS) in both TCGA and CGGA cohorts. In addition, the classification model was verified by IHC staining in 137 GBM samples to further clarify the difference in OS between the five groups. Furthermore, the novel biomarkers of glioblastoma stem cells (GSCs) were also described. In summary, we identified five classifications of GBM and found that they exhibited distinct drug sensitivities and different prognoses, suggesting that the new grouping system may be able to provide important prognostic information and have certain guiding significance for the treatment of GBM, and identified the GSCL cluster in GBM tissues and described its characteristic program, which may help develop new potential therapeutic targets for GSCs in GBM.

RUNX1/CD44 axis regulates the proliferation, migration, and immunotherapy of gliomas: A single-cell sequencing analysis

Background

Glioma is one of the most common, primary, and lethal adult brain tumors because of its extreme aggressiveness and poor prognosis. Several recent studies relevant to the immune function of CD44, a transmembrane glycoprotein as a significant hyaluronic acid receptor, have achieved great success, revealing the critical role of CD44 in immune infiltration in gliomas. The overexpression of CD44 has been verified to correlate with cancer aggressiveness and migration, while the clinical and immune features of CD44 expression have not yet been thoroughly characterized in gliomas.

Methods

Molecular and clinical data of glioma collected from publicly available genomic databases were analyzed.

Results

CD44 was up-expressed in malignant gliomas, notably in the 1p/19q non-codeletion cases, isocitrate dehydrogenase (IDH) wild-type, and mesenchymal subtypes in GBM samples. CD44 expression level strongly correlates with stromal and immune cells, mainly infiltrating the glioma microenvironment by single-cell sequencing analysis. Meanwhile, CD44 can be a promising biomarker in predicting immunotherapy responses and mediating the expression of PD-L1. Finally, RUNX1/CD44 axis could promote the proliferation and migration of gliomas.

Conclusions

Therefore, CD44 was responsible for glioma growth and progression. It could potentially lead to a novel target for glioma immunotherapy or a prognostic biomarker.

Molecular targeted therapy: A new avenue in glioblastoma treatment

Glioblastoma, also referred to as glioblastoma multiforme (GBM), is grade IV astrocytoma characterized by being fast-growing and the most aggressive brain tumor. In adults, it is the most prevalent type of malignant brain tumor. Despite the advancements in both diagnosis tools and therapeutic treatments, GBM is still associated with poor survival rate without any statistically significant improvement in the past three decades. Patient's genome signature is one of the key factors causing the development of this tumor, in addition to previous radiation exposure and other environmental factors. Researchers have identified genomic and subsequent molecular alterations affecting core pathways that trigger the malignant phenotype of this tumor. Targeting intrinsically altered molecules and pathways is seen as a novel avenue in GBM treatment. The present review shed light on signaling pathways and intrinsically altered molecules implicated in GBM development. It discussed the main challenges impeding successful GBM treatment, such as the blood brain barrier and tumor microenvironment (TME), the plasticity and heterogeneity of both GBM and TME and the glioblastoma stem cells. The present review also presented current advancements in GBM molecular targeted therapy in clinical trials. Profound and comprehensive understanding of molecular participants opens doors for innovative, more targeted and personalized GBM therapeutic modalities.

Spatial analysis of the glioblastoma proteome reveals specific molecular signatures and markers of survival

Molecular heterogeneity is a key feature of glioblastoma that impedes patient stratification and leads to large discrepancies in mean patient survival. Here, we analyze a cohort of 96 glioblastoma patients with survival ranging from a few months to over 4 years. 46 tumors are analyzed by mass spectrometry-based spatially-resolved proteomics guided by mass spectrometry imaging. Integration of protein expression and clinical information highlights three molecular groups associated with immune, neurogenesis, and tumorigenesis signatures with high intra-tumoral heterogeneity. Furthermore, a set of proteins originating from reference and alternative ORFs is found to be statistically significant based on patient survival times. Among these proteins, a 5-protein signature is associated with survival. The expression of these 5 proteins is validated by immunofluorescence on an additional cohort of 50 patients. Overall, our work characterizes distinct molecular regions within glioblastoma tissues based on protein expression, which may help guide glioblastoma prognosis and improve current glioblastoma classification.

Identification of hub genes and transcription factor regulatory network for heart failure using RNA-seq data and robust rank aggregation analysis

Background

Heart failure (HF) is the end stage of various cardiovascular diseases with a high mortality rate. Novel diagnostic and therapeutic biomarkers for HF are urgently required. Our research aims to identify HF-related hub genes and regulatory networks using bioinformatics and validation assays.

Methods

Using four RNA-seq datasets in the Gene Expression Omnibus (GEO) database, we screened differentially expressed genes (DEGs) of HF using Removal of Unwanted Variation from RNA-seq data (RUVSeq) and the robust rank aggregation (RRA) method. Then, hub genes were recognized using the STRING database and Cytoscape software with cytoHubba plug-in. Furthermore, reliable hub genes were validated by the GEO microarray datasets and quantitative reverse transcription polymerase chain reaction (qRT-PCR) using heart tissues from patients with HF and non-failing donors (NFDs). In addition, R packages “clusterProfiler” and “GSVA” were utilized for enrichment analysis. Moreover, the transcription factor (TF)–DEG regulatory network was constructed by Cytoscape and verified in a microarray dataset.

Results

A total of 201 robust DEGs were identified in patients with HF and NFDs. STRING and Cytoscape analysis recognized six hub genes, among which ASPN, COL1A1, and FMOD were confirmed as reliable hub genes through microarray datasets and qRT-PCR validation. Functional analysis showed that the DEGs and hub genes were enriched in T-cell-mediated immune response and myocardial glucose metabolism, which were closely associated with myocardial fibrosis. In addition, the TF–DEG regulatory network was constructed, and 13 significant TF–DEG pairs were finally identified.

Conclusion

Our study integrated different RNA-seq datasets using RUVSeq and the RRA method and identified ASPN, COL1A1, and FMOD as potential diagnostic biomarkers for HF. The results provide new insights into the underlying mechanisms and effective treatments of HF.

Bioinformatic profiling identifies the glutaminase to be a potential novel cuproptosis-related biomarker for glioma

Glioma is the most common tumour of the central nervous system, with a poor prognosis and an increasing trend of incidence in recent years; it is also beginning to affect younger age groups more. Added to this, cuproptosis is a new form of cell death. Indeed, when a certain amount of copper accumulates in a cell, it affects specific mitochondrial metabolic enzymes in that cell and leads to cell death–a phenomenon known as cuproptosis. In this study, we applied bioinformatics analysis, and, according to the results of the study analysis and Gene Ontology (GO), as well as the Kyoto Encyclopedia of Genes and Genomes KyotoEncyclopediaofGenesandGenomes, the glutaminase (GLS) genes affect the prognosis and tumour mutation of glioma patients through cuproptosis. Interestingly, however, GLS is not involved in the immune escape of glioma. Glutaminase genes are a class of glucose metabolism-related genes that are involved in the tricarboxylic acid cycle of cells. At the same time, the expression of the glutaminase gene was positively correlated with the degree of immune cell infiltration and the expression of various immune cell markers, and thus affected the prognosis of glioma patients. Therefore, we believe that the cuproptosis-related glutaminase gene can be an important factor in determining the prognosis of glioma patients.

Nanotherapeutic treatment of the invasive glioblastoma tumor microenvironment

Glioblastoma (GBM) is the most common malignant adult brain cancer with no curative treatment strategy. A significant hurdle in GBM treatment is effective therapeutic delivery to the brain-invading tumor cells that remain following surgery within functioning brain regions. Developing therapies that can either directly target these brain-invading tumor cells or act on other cell types and molecular processes supporting tumor cell invasion and recurrence are essential steps in advancing new treatments in the clinic. This review highlights some of the drug delivery strategies and nanotherapeutic technologies that are designed to target brain-invading GBM cells or non-neoplastic, invasion-supporting cells residing within the GBM tumor microenvironment.

Clinical Characterization and Prognostic Value of TPM4 and Its Correlation with Epithelial–Mesenchymal Transition in Glioma

Tropomyosin 4 (TPM4) has been reported as an oncogenic gene across different malignancies. However, the role of TPM4 in glioma remains unclear. This study aimed to determine the clinical characterization and prognostic value of TPM4 in gliomas. Transcriptome expression and clinical information were collected from the CGGA and TCGA datasets, which included 998 glioma patients. ScRNA-seq data were obtained from CGGA. R software was utilized for statistical analyses. There was a positive correlation between TPM4 and WHO grades. IDH-wildtype and mesenchymal subtype gliomas were accompanied by TPM4 upregulation. GO and GSEA analysis suggested that TPM4 was profoundly associated with epithelial-to-mesenchymal transition (EMT). Subsequent GSVA revealed a robust correlation between TPM4 and three signaling pathways of EMT (hypoxia, TGF-β, PI3K/AKT). Furthermore, TPM4 showed a synergistic effect with mesenchymal biomarkers, particularly with N-cadherin, Slug, Snail, TWIST1, and vimentin. ScRNA-seq analysis suggested that higher TPM4 was mainly attributed to tumor cells and macrophages and associated with tumor cell progression and macrophage polarization. Finally, high TPM4 was significantly associated with unfavorable outcomes. In conclusion, our findings indicate that TPM4 is significantly correlated with more malignant characteristics of gliomas, potentially through involvement in EMT. TPM4 could predict worse survival for patients with glioma.

Pan-Cancer Analysis of the Immunological Role of PDIA5: A Potential Target for Immunotherapy

The aberrant protein disulfide isomerase A5 (PDIA5) expression was relevant to the poor prognosis of patients with human cancers. However, its relationship with the epigenetic and genetic alterations and its effect on tumor immunity is still lacking. In the present study, we comprehensively analyzed the immune infiltration role of PDIA5 in human cancers based on large-scale bioinformatics analyses and in vitro experiments. Obvious DNA methylation and moderate alteration frequency of PDIA5 were observed in human cancers. The expression level of PDIA5 was significantly correlated with infiltrated immune cells, immune pathways, and other immune signatures. We found that cancer cells and macrophages exhibited high PDIA5 expression in human cancers using the single-cell RNA sequencing analysis. We also demonstrated the interaction between PDIA5 and immune cells in glioblastoma multiforme (GBM). Multiplex immunofluorescence staining showed the upregulated expression level of PDIA5 and the increased number of M2 macrophage markers-CD163 positive cells in pan-cancer samples. Notably, PDIA5 silencing resulted in upregulated expression of PD-L1 and SPP1 in U251 cells. Silencing of PDIA5 in hepG2 cells, U251 cells, and PC3 cells contributed to a decline in their ability of proliferation, clone formation, and invasion and inhibited the migration of cocultured M2 macrophages. Additionally, PDIA5 also displayed predictive value in the immunotherapy response of both murine and human cancer cohorts. Overall, our findings indicated that PDIA5 might be a potential target for immunotherapies in cancers.

Prognostic Significance of Pseudotime from Texture Parameters of FDG PET/CT in Locally Advanced Non-Small-Cell Lung Cancer with Tri-Modality Therapy

Simple Summary

Although texture parameters of F-18 fluorodeoxyglucose positron emission tomography/computed tomography images were known to associate tumor biology and clinical features, the types and implications of parameters are too various and complicated. To overcome the limitation of texture parameter, we attempted to produce a new simplified parameter from texture parameters of F-18 fluorodeoxyglucose positron emission tomography/computed tomography images in lung cancer patients using pseudotime analysis. Pseudotime analysis is a recently developed method to explore changes in cell or tissue characteristics based on transcriptomic expression. It is the first study to apply pseudotime analysis into radiomics dataset other than transcriptomics data. Herein, we demonstrated that pseudotime can be successfully estimated from texture parameters. In the aspect of prognostic prediction, pseudotime was an independent prognostic factor for overall survival in contrast to conventional parameters such as metabolic tumor volume and total lesion glycolysis. This study showed possibility of integrating various texture parameters into single parameter which reflects disease progression status. Pseudotime, as a concrete value of disease progression, is expected to be used in clinical field to evaluate disease and predict prognosis.

Abstract

Texture analysis provides image parameters from F-18 fluorodeoxyglucose positron emission tomography/computed tomography (FDG PET/CT). Although some parameters are associated with tumor biology and clinical features, the types and implications of these parameters are complicated. We applied pseudotime analysis, which has recently been used to estimate changes in individual sample characteristics, to texture parameters from FDG PET/CT images of locally advanced non-small-cell lung cancer (NSCLC) patients undergoing neoadjuvant concurrent chemoradiation therapy (CCRT) followed by surgery. Our subjects were 303 NSCLC patients who underwent pretherapeutic FDG PET/CT and tri-modality therapy. Texture parameters of the primary tumor were calculated from FDG PET/CT images acquired before neoadjuvant CCRT. Pseudotime analysis was performed using the PhenoPath tool. Clinicopathologic features including survival data were collected and survival analysis was performed to compare the prognostic significances of pseudotime parameters with those of conventional PET parameters. Pseudotime was successfully estimated from texture parameters. Normalized co-occurrence homogeneity, normalized co-occurrence inverse difference moment, and black–white symmetry showed positive correlations with pseudotime, short run emphasis, normalized co-occurrence dissimilarity, and short zone emphasis negative correlation. The maximum standardized uptake value (SUV) and mean SUV were not associated with overall survival. Pseudotime, metabolic tumor volume (MTV), and total lesion glycolysis (TLG) showed significant associations with overall survival. In contrast to MTV and TLG, pseudotime was an independent prognostic factor for overall survival. Various metabolic texture parameters can be integrated into a single parameter using pseudotime analysis. Pseudotime of the primary tumor, estimated from FDG PET/CT images, better predicts overall survival in locally advanced NSCLC patients treated with tri-modality therapy than conventional PET parameters.

Intercellular Communication-Related Molecular Subtypes and a Gene Signature Identified by the Single-Cell RNA Sequencing Combined with a Transcriptomic Analysis

Background. The tumor microenvironment (TME) of lung adenocarcinoma (LUAD) comprise various cell types that communicate with each other through ligand-receptor interactions. This study focused on the identification of cell types in LUAD by single-cell RNA sequencing (scRNA-seq) data and screening of intercellular communication-related genes. Methods. The Gene Expression Omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/geo) provided the RNA-seq data of LUAD patients in the GSE149655, GSE31210, and GSE72094 datasets. Quality control of the scRNA-seq data in GSE149655 was performed by the Seurat package (http://seurat.r-forge.r-project.org) for identifying highly variable genes for principal component analysis (PCA) and cell clustering. The CellPhoneDB (http://www.cellphonedb.org) was used for filtering intercellular communication-related ligand-receptor pairs. According to ligand and receptor expressions, LUAD samples were clustered using ConsensusClusterPlus (https://www.bioconductor.org/packages/release/bioc/html/ConsensusClusterPlus). Additionally, the identification of prognosis-related ligand and receptor genes was conducted along with the development of a risk prediction model by the least absolute shrinkage and selection operator (LASSO) Cox regression analysis. Results. This study identified twelve cell types in 8170 cells of LUAD tissues along with 219 ligand and receptor genes. LUAD was classified into three different molecular subtypes, among which cluster 3 (C3) had the longest overall survival (OS) time and cluster (C1) had the shortest OS time. In comparison with the other two molecular subtypes, it was observed that C1 had a higher rate of somatic mutations and lower levels of infiltrating immune cells and immune scores. Ten genes were screened from the total ligand and receptor genes to construct a risk model, which showed a strong prediction power in the prognosis of patients with LUAD. Conclusion. The results of this study revealed cell types specific to LUAD, which were classified into different molecular subtypes according to intercellular communication-related genes. A novel prognostic risk model was developed in this study, providing new insights into prognostic assessment models for LUAD.

Single-cell transcriptome highlights a multilayer regulatory network on an invasive trajectory within colorectal cancer progression

PURPOSE

Colorectal cancer (CRC) is one of the most common and fatal gastrointestinal malignancies, in which cancer stem cells (CSCs) were identified to enable tumor heterogeneity and initiate tumor formation. However, the process from CSCs to invasion cells is unconfirmed.

METHODS

Several bioinformatics methods, including clustering, pseudotime analysis, gene set variation analysis and gene ontology enrichment, were used to construct a path of gradual transformation of CSCs to invasive cells, called "stem-to-invasion path". A large amount of signaling interactions were collated to build the multilayer regulatory network. Kaplan-Meier curve and time-dependent ROC method were applied to reveal prognostic values.

RESULTS

We validated the heterogeneity of cells in the tumor microenvironment and revealed the presence of malignant epithelial cells with high invasive potential within primary colonic carcinomas. Next, the "stem-to-invasion path" was identified through constructing a branching trajectory with cancer cells arranged in order. A multilayer regulatory network considered as the vital factor involved in acquiring invasion characteristics underlying the path was built to elucidate the interactions between tumor cell and tumor-associated microenvironment. Then we further identified a novel combinatorial biomarker that can be used to assess the prognosis for CRC patients, and validated its predictive robustness on the independent dataset.

CONCLUSION

Our work provides new insights into the acquisition of invasive potential in primary tumor cells, as well as potential therapeutic targets for CRC invasiveness, which may be useful for the cancer research and clinical treatment.

TNFSF13 Is a Novel Onco-Inflammatory Marker and Correlates With Immune Infiltration in Gliomas

Existing therapeutic strategies for gliomas are restricted; hence, exploration for novel diagnostic indicator and treatment is essential. Here, we performed bioinformatic analyses for TNFSF13 (also known as APRIL), a proliferation-inducing ligand of the tumor necrosis factor (TNF) superfamily, aiming to assess its potential for predicting glioma patient’s prognosis and targeted therapy. TNFSF13 expression was upregulated in the increase of tumor grades based on Xiangya cohort. In high TNFSF13 gliomas, somatic mutation was proved to correlate with amplification of EGFR and deletion of CDKN2A; while mutation of IDH1 was more frequently observed in low TNFSF13 group. We also confirmed the positive correlation between TNFSF13 and infiltrating immune and stromal cells in glioma microenvironment. Further, TNFSF13 was found to be involved in immunosuppression via diverse immunoregulation pathways and was associated with other immune checkpoints and inflammation. Single-cell sequencing revealed an abundant expression of TNFSF13 in neoplastic cells and M2 macrophages, which TNFSF13 might potentially regulate the cell communication via IL-8, C3, and CD44. Lastly, TNFSF13 mediated the activities of transcription factors including FOXO3, MEIS2, and IRF8. Our analyses demonstrated the relevance between TNFSF13 and glioma progress and indicated the potential of TNFSF13 as a novel diagnostic onco-inflammatory biomarker and immunotherapy target of gliomas.

Expression and crucial role of long non-coding RNA FGD5-AS1 in human cancers

Long non-coding RNAs (lncRNAs) are transcribed by RNA polymerase II and are longer than 200 nucleotides. Several studies have revealed that lncRNAs are important regulators of cancer progression. The lncRNA FGD5-AS1, first identified in 2018, has emerged as a crucial regulator of processes related to carcinogenesis. The expression levels of FGD5-AS1 are known to be significantly up-regulated in a variety of human cancers. Moreover, FGD5-AS1 expression closely correlates with clinical features and poor prognosis and its expression has been shown to attenuate cell proliferation, cell migration, cell invasiveness, drug resistance, and the epithelial-mesenchymal transition through several pathways. Here, we provide an overview of the role of FGD5-AS1 in various cancers and discuss its potential clinical utility in tumor progression. In addition, we used a gene expression profiling interactive analysis dataset to explore associations between FGD5-AS1 pan-cancer expressions and prognoses.

Uncovering Spatiotemporal Heterogeneity of High-Grade Gliomas: From Disease Biology to Therapeutic Implications

Glioblastomas (GBM) are the most common and aggressive tumors of the central nervous system. Rapid tumor growth and diffuse infiltration into healthy brain tissue, along with high intratumoral heterogeneity, challenge therapeutic efficacy and prognosis. A better understanding of spatiotemporal tumor heterogeneity at the histological, cellular, molecular, and dynamic levels would accelerate the development of novel treatments for this devastating brain cancer. Histologically, GBM is characterized by nuclear atypia, cellular pleomorphism, necrosis, microvascular proliferation, and pseudopalisades. At the cellular level, the glioma microenvironment comprises a heterogeneous landscape of cell populations, including tumor cells, non-transformed/reactive glial and neural cells, immune cells, mesenchymal cells, and stem cells, which support tumor growth and invasion through complex network crosstalk. Genomic and transcriptomic analyses of gliomas have revealed significant inter and intratumoral heterogeneity and insights into their molecular pathogenesis. Moreover, recent evidence suggests that diverse dynamics of collective motion patterns exist in glioma tumors, which correlate with histological features. We hypothesize that glioma heterogeneity is not stochastic, but rather arises from organized and dynamic attributes, which favor glioma malignancy and influences treatment regimens. This review highlights the importance of an integrative approach of glioma histopathological features, single-cell and spatially resolved transcriptomic and cellular dynamics to understand tumor heterogeneity and maximize therapeutic effects.

Single-cell RNA sequencing in cancer: Applications, advances, and emerging challenges

Cancer has become one of the greatest threats to human health, and new technologies are urgently needed to further clarify the mechanisms of cancer so that better detection and treatment strategies can be developed. At present, extensive genomic analysis and testing of clinical specimens shape the insights into carcinoma. Nevertheless, carcinoma of humans is a complex ecosystem of cells, including carcinoma cells and immunity-related and stroma-related subsets, with accurate characteristics obscured by extensive genome-related approaches. A growing body of research shows that sequencing of single-cell RNA (scRNA-seq) is emerging to be an effective way for dissecting human tumor tissue at single-cell resolution, presenting one prominent way for explaining carcinoma biology. This review summarizes the research progress of scRNA-seq in the field of tumors, focusing on the application of scRNA-seq in tumor circulating cells, tumor stem cells, tumor drug resistance, the tumor microenvironment, and so on, which provides a new perspective for tumor research.

A comprehensive prognostic signature for glioblastoma patients based on transcriptomics and single cell sequencing

PURPOSE

Glioblastoma (GBM) is the most common and deadly brain tumor. We aimed to reveal potential prognostic GBM marker genes, elaborate their functions, and build an effective a prognostic model for GBM patients.

METHODS

Through data mining of The Cancer Genome Atlas (TCGA) and the Chinese Glioma Genome Atlas (CGGA), we screened for significantly differentially expressed genes (DEGs) to calculate risk scores for individual patients. Published data of somatic mutation and copy number variation profiles were analyzed for distinct genomic alterations associated with risk scores. In addition, single-cell sequencing was used to explore the biological functions of the identified prognostic marker genes. By combining risk scores and other clinical features, we built a comprehensive prognostic GBM model.

RESULTS

Seven DEGs (CLEC5A, HOXC6, HOXA5, CCL2, GPRASP1, BSCL2 and PTX3) were identified as being prognostic for GBM. Expression of these genes was confirmed in different GBM cell lines using real-time PCR. Risk scores calculated from the seven DEGs revealed prognostic value irrespective of other clinical factors, including IDH mutation status, and were negatively correlated with TP53 expression. The prognostic genes were found to be associated with tumor proliferation and progression based on pseudo-time analysis in neoplastic cells. A final prognostic model was developed and validated with a good performance, especially in geriatric GBM patients.

CONCLUSIONS

Using genetic profiles, age, IDH mutation status, and chemotherapy and radiotherapy, we constructed a comprehensive prognostic model for GBM patients. The model has a good performance, especially in geriatric GBM patients.

Musashi1 Contribution to Glioblastoma Development via Regulation of a Network of DNA Replication, Cell Cycle and Division Genes

RNA-binding proteins (RBPs) function as master regulators of gene expression. Alterations in their levels are often observed in tumors with numerous oncogenic RBPs identified in recent years. Musashi1 (Msi1) is an RBP and stem cell gene that controls the balance between self-renewal and differentiation. High Msi1 levels have been observed in multiple tumors including glioblastoma and are often associated with poor patient outcomes and tumor growth. A comprehensive genomic analysis identified a network of cell cycle/division and DNA replication genes and established these processes as Msi1's core regulatory functions in glioblastoma. Msi1 controls this gene network via two mechanisms: direct interaction and indirect regulation mediated by the transcription factors E2F2 and E2F8. Moreover, glioblastoma lines with Msi1 knockout (KO) displayed increased sensitivity to cell cycle and DNA replication inhibitors. Our results suggest that a drug combination strategy (Msi1 + cell cycle/DNA replication inhibitors) could be a viable route to treat glioblastoma.

Involvement of the Catecholamine Pathway in Glioblastoma Development

Glioblastoma (GBM) is the most aggressive tumor of the central nervous system (CNS). The standard of care improves the overall survival of patients only by a few months. Explorations of new therapeutic targets related to molecular properties of the tumor are under way. Even though neurotransmitters and their receptors normally function as mediators of interneuronal communication, growing data suggest that these molecules are also involved in modulating the development and growth of GBM by acting on neuronal and glioblastoma stem cells. In our previous DNA CpG methylation studies, gene ontology analyses revealed the involvement of the monoamine pathway in sequential GBM. In this follow-up study, we quantitated the expression levels of four selected catecholamine pathway markers (alpha 1D adrenergic receptor-ADRA1D; adrenergic beta receptor kinase 1 or G protein-coupled receptor kinase 2-ADRBK1/GRK2; dopamine receptor D2-DRD2; and synaptic vesicle monoamine transporter-SLC18A2) by immunohistochemistry, and compared the histological scores with the methylation levels within the promoters + genes of these markers in 21 pairs of sequential GBM and in controls. Subsequently, we also determined the promoter and gene methylation levels of the same markers in an independent database cohort of sequential GBM pairs. These analyses revealed partial inverse correlations between the catecholamine protein expression and promoter + gene methylation levels, when the tumor and control samples were compared. However, we found no differences in the promoter + gene methylation levels of these markers in either our own or in the database primary-recurrent GBM pairs, despite the higher protein expression of all markers in the primary samples. This observation suggests that regulation of catecholamine expression is only partially related to CpG methylation within the promoter + gene regions, and additional mechanisms may also influence the expression of these markers in progressive GBM. These analyses underscore the involvement of certain catecholamine pathway markers in GBM development and suggest that these molecules mediating or modulating tumor growth merit further exploration.

PDIA5 is Correlated With Immune Infiltration and Predicts Poor Prognosis in Gliomas

Gliomas are the most common and lethal primary malignant tumor of the brain. Routine treatment including surgical resection, chemotherapy, and radiotherapy produced limited therapeutic effect, while immunotherapy targeting the glioma microenvironment has offered a novel therapeutic option. PDIA5 protein is the member of PDI family, which is highly expressed in glioma and participates in glioma progression. Based on large-scale bioinformatics analysis, we discovered that PDIA5 expression level is upregulated in aggressive gliomas, with high PDIA5 expression predicting poor clinical outcomes. We also observed positive correlation between PDIA5 and immune infiltrating cells, immune related pathways, inflammatory activities, and other immune checkpoint members. Patients with high PDIA5 high-expression benefited from immunotherapies. Additionally, immunohistochemistry revealed that PDIA5 and macrophage biomarker CD68 were upregulated in high-grade gliomas, and patients with low PDIA5 level experienced favorable outcomes among 33 glioma patients. Single cell RNA sequencing exhibited that PDIA5 was in high level presenting in neoplastic cells and macrophages. Cell transfection and co-culture of glioma cells and macrophages revealed that PDIA5 in tumor cells mediated macrophages exhausting. Altogether, our findings indicate that PDIA5 overexpression is associated with immune infiltration in gliomas, and may be a promising therapeutic target for glioma immunotherapy.

Dissecting the Invasion-Associated Long Non-coding RNAs Using Single-Cell RNA-Seq Data of Glioblastoma

Glioblastoma (GBM) is characterized by rapid and lethal infiltration of brain tissue, which is the primary cause of treatment failure and deaths for GBM. Therefore, understanding the molecular mechanisms of tumor cell invasion is crucial for the treatment of GBM. In this study, we dissected the single-cell RNA-seq data of 3345 cells from four patients and identified dysregulated genes including long non-coding RNAs (lncRNAs), which were involved in the development and progression of GBM. Based on co-expression network analysis, we identified a module (M1) that significantly overlapped with the largest number of dysregulated genes and was confirmed to be associated with GBM invasion by integrating EMT signature, experiment-validated invasive marker and pseudotime trajectory analysis. Further, we denoted invasion-associated lncRNAs which showed significant correlations with M1 and revealed their gradually increased expression levels along the tumor cell invasion trajectory, such as VIM-AS1, WWTR1-AS1, and NEAT1. We also observed the contribution of higher expression of these lncRNAs to poorer survival of GBM patients. These results were mostly recaptured in another validation data of 7930 single cells from 28 GBM patients. Our findings identified lncRNAs that played critical roles in regulating or controlling cell invasion and migration of GBM and provided new insights into the molecular mechanisms underlying GBM invasion as well as potential targets for the treatment of GBM.

The essential role of long non-coding RNA GAS5 in glioma: interaction with microRNAs, chemosensitivity and potential as a biomarker

Glioma is a malignant brain tumor with a generally poor prognosis. Dysregulation of a long non-coding RNA, GAS5, has been detected in numerous cancers, including glioma. Previous studies have suggested that GAS5 plays a significant functional role in glioma, affecting proliferation, metastasis, invasion, and apoptosis. In this review, we describe the roles and mechanisms of GAS5 in glioma. GAS5 may be a biomarker for diagnosis and prognosis, and even a potential target for glioma treatment, and therefore warrants further investigation.

Toward uncharted territory of cellular heterogeneity: advances and applications of single-cell RNA-seq

Among single-cell analysis technologies, single-cell RNA-seq (scRNA-seq) has been one of the front runners in technical inventions. Since its induction, scRNA-seq has been well received and undergone many fast-paced technical improvements in cDNA synthesis and amplification, processing and alignment of next generation sequencing reads, differentially expressed gene calling, cell clustering, subpopulation identification, and developmental trajectory prediction. scRNA-seq has been exponentially applied to study global transcriptional profiles in all cell types in humans and animal models, healthy or with diseases, including cancer. Accumulative novel subtypes and rare subpopulations have been discovered as potential underlying mechanisms of stochasticity, differentiation, proliferation, tumorigenesis, and aging. scRNA-seq has gradually revealed the uncharted territory of cellular heterogeneity in transcriptomes and developed novel therapeutic approaches for biomedical applications. This review of the advancement of scRNA-seq methods provides an exploratory guide of the quickly evolving technical landscape and insights of focused features and strengths in each prominent area of progress.

Towards multi-omics characterization of tumor heterogeneity: a comprehensive review of statistical and machine learning approaches

The multi-omics molecular characterization of cancer opened a new horizon for our understanding of cancer biology and therapeutic strategies. However, a tumor biopsy comprises diverse types of cells limited not only to cancerous cells but also to tumor microenvironmental cells and adjacent normal cells. This heterogeneity is a major confounding factor that hampers a robust and reproducible bioinformatic analysis for biomarker identification using multi-omics profiles. Besides, the heterogeneity itself has been recognized over the years for its significant prognostic values in some cancer types, thus offering another promising avenue for therapeutic intervention. A number of computational approaches to unravel such heterogeneity from high-throughput molecular profiles of a tumor sample have been proposed, but most of them rely on the data from an individual omics layer. Since the heterogeneity of cells is widely distributed across multi-omics layers, methods based on an individual layer can only partially characterize the heterogeneous admixture of cells. To help facilitate further development of the methodologies that synchronously account for several multi-omics profiles, we wrote a comprehensive review of diverse approaches to characterize tumor heterogeneity based on three different omics layers: genome, epigenome and transcriptome. As a result, this review can be useful for the analysis of multi-omics profiles produced by many large-scale consortia. Contact:sunkim.bioinfo@snu.ac.kr.

Advances in histone deacetylase inhibitors in targeting glioblastoma stem cells

Glioblastoma multiforme (GBM) is a lethal grade IV glioma (WHO classification) and widely prevalent primary brain tumor in adults. GBM tumors harbor cellular heterogeneity with the presence of a small subpopulation of tumor cells, described as GBM cancer stem cells (CSCs) that pose resistance to standard anticancer regimens and eventually mediate aggressive relapse or intractable progressive GBM. Existing conventional anticancer therapies for GBM do not target GBM stem cells and are mostly palliative; therefore, exploration of new strategies to target stem cells of GBM has to be prioritized for the development of effective GBM therapy. Recent developments in the understanding of GBM pathophysiology demonstrated dysregulation of epigenetic mechanisms along with the genetic changes in GBM CSCs. Altered expression/activity of key epigenetic regulators, especially histone deacetylases (HDACs) in GBM stem cells has been associated with poor prognosis; inhibiting the activity of HDACs using histone deacetylase inhibitors (HDACi) has been promising as mono-therapeutic in targeting GBM and in sensitizing GBM stem cells to an existing anticancer regimen. Here, we review the development of pan/selective HDACi as potential anticancer agents in targeting the stem cells of glioblastoma as a mono or combination therapy.

A Single Cell but Many Different Transcripts: A Journey into the World of Long Non-Coding RNAs

In late 2012 it was evidenced that most of the human genome is transcribed but only a small percentage of the transcripts are translated. This observation supported the importance of non-coding RNAs and it was confirmed in several organisms. The most abundant non-translated transcripts are long non-coding RNAs (lncRNAs). In contrast to protein-coding RNAs, they show a more cell-specific expression. To understand the function of lncRNAs, it is fundamental to investigate in which cells they are preferentially expressed and to detect their subcellular localization. Recent improvements of techniques that localize single RNA molecules in tissues like single-cell RNA sequencing and fluorescence amplification methods have given a considerable boost in the knowledge of the lncRNA functions. In recent years, single-cell transcription variability was associated with non-coding RNA expression, revealing this class of RNAs as important transcripts in the cell lineage specification. The purpose of this review is to collect updated information about lncRNA classification and new findings on their function derived from single-cell analysis. We also retained useful for all researchers to describe the methods available for single-cell analysis and the databases collecting single-cell and lncRNA data. Tables are included to schematize, describe, and compare exposed concepts.