Stemness Related Genes Revealed by Network Analysis Associated With Tumor Immune Microenvironment and the Clinical Outcome in Lung Adenocarcinoma

NUSAP1 Promotes Immunity and Apoptosis by the SHCBP1/JAK2/STAT3 Phosphorylation Pathway to Induce Dendritic Cell Generation in Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) is the most common type of liver cancer and is associated with high morbidity and mortality rates. The aims of this study were to investigate the immune-promoting action of nucleolar and spindle-associated protein 1 (NUSAP1) and identify an immunotherapy target for HCC. The Cancer Genome Atlas (TCGA) was used to analyze interaction molecules and immune correlation. The interaction between NUSAP1 and SHC binding and spindle associated 1 (SHCBP1) was examined. The role of the SHCBP1/Janus kinase 2/signal transducer and activator of transcription 3 (SHCBP1/JAK2/STAT3) pathway in this process was explored. After co-culture with HCC cell lines, the differentiation of peripheral blood mononuclear cells (PBMCs) into dendritic cells (DC) was evaluated by measuring the expression of surface factors CD1a and CD86. Pathological tissues from 50 patients with HCC were collected to validate the results of cell experiments. The expression levels of CD1a and CD86 in tissues were also determined. The results show that NUSAP1 interacted with SHCBP1 and was positively correlated with DC. In HCC cell lines, an interaction was observed between NUSAP1 and SHCBP1. It was verified that NUSAP1 inhibited the JAK2/STAT3 phosphorylation pathway by blocking SHCBP1. After co-culture, the levels of CD1a and CD86 in PBMC were elevated. In the clinical specimens, CD1a and CD86 expression levels were significantly higher in the high-NUSAP1 group versus the low-NUSAP1 group. In Summary, NUSAP1 enhanced immunity by inhibiting the SHCBP1/JAK2/STAT3 phosphorylation pathway and promoted DC generation and HCC apoptosis. NUSAP1 may be a target of immunotherapy for HCC.

Role of cell division cycle-associated proteins in regulating cell cycle and promoting tumor progression

The cell division cycle-associated protein (CDCA) family is important in regulating cell division. High CDCA expression is significantly linked to tumor development. This review summarizes clinical and basic studies on CDCAs conducted in recent decades. Furthermore, it systematically introduces the molecular expression and function, key mechanisms, cell cycle regulation, and roles of CDCAs in tumor development, cell proliferation, drug resistance, invasion, and metastasis. Additionally, it presents the latest research on tumor diagnosis, prognosis, and treatment targeting CDCAs. These findings are pivotal for further in-depth studies on the role of CDCAs in promoting tumor development and provide theoretical support for their application as new anti-tumor targets.

Upregulation of sperm-associated antigen 5 expression in endometrial carcinoma was associated with poor prognosis and immune dysregulation, and promoted cell migration and invasion

Sperm-associated antigen 5 (SPAG5) regulates cancer cell invasion and is involved in the progression of many cancers. However, the role of SPAG5 in endometrial carcinoma (EC) is still unknown. The purpose of this study was to explore the role of SPAG5 in EC and its potential molecular mechanism. The UALCAN tool and cBioPortal were used to analyze the expression and alterations of SPAG5 in EC, respectively. OncoLnc was used for survival analysis. We analyzed the effects of SPAG5 on immune cell infiltration and the expression levels of immune checkpoints. We also overexpressed and knocked down SPAG5 in EC cells to explore the effect of SPAG5 regulation on migration, invasion, apoptosis, and the cell cycle of EC cells. We found that SPAG5 was overexpressed and the SPAG5 gene was often mutated in EC. High SPAG5 expression was significantly associated with poor overall survival in patients with EC. SPAG5 also affected the level of immune cell infiltration in the TIME and the expression of immune checkpoints lymphocyte activating 3 (LAG3) and T cell immunoreceptor with Ig and ITIM domains (TIGIT) in patients with EC. It may also be involved in the immunotherapy response in these patients. In vitro experiments showed that SPAG5 promotes cancer cell migration and invasion. In conclusion, this study lays the foundation for further understanding the molecular mechanisms of EC involving SPAG5 and contributes to diagnosing and managing this disease.

High-throughput transcriptome profiling indicates ribosomal RNAs to be associated with resistance to immunotherapy in non-small cell lung cancer (NSCLC)

BACKGROUND

Despite the impressive outcomes with immune checkpoint inhibitor (ICI) in non-small cell lung cancer (NSCLC), only a minority of the patients show long-term benefits from ICI. In this study, we used retrospective cohorts of ICI treated patients with NSCLC to discover and validate spatially resolved protein markers associated with resistance to programmed cell death protein-1 (PD-1) axis inhibition.

METHODS

Pretreatment samples from 56 patients with NSCLC treated with ICI were collected and analyzed in a tissue microarray (TMA) format in including four different tumor regions per patient using the GeoMx platform for spatially informed transcriptomics. 34 patients had assessable tissue with tumor compartment in all 4 TMA spots, 22 with leukocyte compartment and 12 with CD68 compartment. The patients' tissue that was not assessable in fourfold redundancy in each compartment was designated as the validation cohort; cytokeratin (CK) (N=22), leukocytes CD45 (N=31), macrophages, CD68 (N=43). The human whole transcriptome, represented by~18,000 individual genes assessed by oligonucleotide-tagged in situ hybridization, was sequenced on the NovaSeq platform to quantify the RNAs present in each region of interest.

RESULTS

54,000 gene variables were generated per case, from them 25,740 were analyzed after removing targets with expression lower than a prespecified frequency. Cox proportional-hazards model analysis was performed for overall and progression-free survival (OS, PFS, respectively). After identifying genes significantly associated with limited survival benefit (HR>1)/progression per spot per patient, we used the intersection of them across the four TMA spots per patient. This resulted in a list of 12 genes in the tumor-cell compartment (RPL13A, GNL3, FAM83A, CYBA, ACSL4, SLC25A6, EPAS1, RPL5, APOL1, HSPD1, RPS4Y1, ADI1). RPL13A, GNL3 in tumor-cell compartment were also significantly associated with OS and PFS, respectively, in the validation cohort (CK: HR, 2.48; p=0.02 and HR, 5.33; p=0.04). In CD45 compartment, secreted frizzled-related protein 2, was associated with OS in the discovery cohort but not in the validation cohort. Similarly, in the CD68 compartment ARHGAP and PNN interacting serine and arginine rich protein were significantly associated with PFS and OS, respectively, in the majority but not all four spots per patient.

CONCLUSION

This work highlights RPL13A and GNL3 as potential indicative biomarkers of resistance to PD-1 axis blockade that might help to improve precision immunotherapy strategies for lung cancer.

Unraveling the tumor immune microenvironment of lung adenocarcinoma using single-cell RNA sequencing

Tumor immune microenvironment (TIME) and its indications for lung cancer patient prognosis and therapeutic response have become new hotspots in cancer research in recent years. Tumor cells, immune cells, various regulatory factors, and their interactions in the TIME have been suggested to commonly influence lung cancer development and therapeutic outcome. The heterogeneity of TIME is composed of dynamic immune-related components, including various cancer cells, immune cells, cytokine/chemokine environments, cytotoxic activity, or immunosuppressive factors. The specific composition of cell subtypes may facilitate or hamper the response to immunotherapy and influence patient prognosis. Various markers have been found to stratify the patient prognosis or predict the therapeutic outcome. In this article, we systematically reviewed the recent advancement of TIME studies in lung adenocarcinoma (LUAD) using single-cell RNA sequencing (scRNA-seq) techniques, with specific focuses on the roles of TIME in LUAD development, TIME heterogeneity, indications of TIME in patient prognosis and therapeutic response during immunotherapy and drug resistance. The main findings in TIME heterogeneity and relevant markers or models for prognosis stratification and response prediction have been summarized. We hope that this review provides an overview of TIME status in LUAD and an inspiration for future development of strategies and biomarkers in LUAD treatment.

Targeted immunotherapy to cancer stem cells: A novel strategy of anticancer immunotherapy

Cancer is a major disease that endangers human health. Cancer drug resistance and relapse are the two main causes contributing to cancer treatment failure. Cancer stem cells (CSCs) are a small fraction of tumor cells that are responsible for tumorigenesis, metastasis, relapse, and resistance to conventional anticancer therapies. Therefore, CSCs are considered to be the root of cancer recurrence, metastasis, and drug resistance. Novel anticancer strategies need to face this new challenge and explore their efficacy against CSCs. Recently, immunotherapy has made rapid advances in cancer treatment, and its potential against CSCs is also an interesting area of research. Meanwhile, immunotherapy strategies are novel therapeutic modalities with promising results in targeting CSCs. In this review, we summarize the targeting of CSCs by various immunotherapy strategies such as monoclonal antibodies(mAb), tumor vaccines, immune checkpoint inhibitors, and chimeric antigen receptor-T cells(CAR-T) in pre-clinical and clinical studies. This review provides new insights into the application of these immunotherapeutic approaches to potential anti-tumor therapies in the future.

KIF2C: An important factor involved in signaling pathways, immune infiltration, and DNA damage repair in tumorigenesis

BACKGROUNDS

Poorly regulated mitosis and chromosomal instability are common characteristics in malignant tumor cells. Kinesin family member 2 C (KIF2C), also known as mitotic centromere-associated kinesin (MCAK) is an essential component during mitotic regulation. In recent years, KIF2C was shown to be dysregulated in several tumors and was involved in many aspects of tumor self-regulation. Research on KIF2C may be a new direction and target for anti-tumor therapy.

OBJECT

The article aims at reviewing current literatures and summarizing the research status of KIF2C in malignant tumors as well as the oncogenic signaling pathways associated with KIF2C and its role in immune infiltration.

RESULT

In this review, we summarize the KIF2C mechanisms and signaling pathways in different malignant tumors, and briefly describe its involvement in pathways related to classical chemotherapeutic drug resistance, such as MEK/ERK, mTOR, Wnt/β-catenin, P53 and TGF-β1/Smad pathways. KIF2C upregulation was shown to promote tumor cell migration, invasion, chemotherapy resistance and inhibit DNA damage repair. It was also highly correlated with microRNAs, and CD4 +T cell and CD8 +T cell tumor immune infiltration.

CONCLUSION

This review shows that KIF2C may function as a new anticancer drug target with great potential for malignant tumor treatment and the mitigation of chemotherapy resistance.

Oncogenic and immunological roles of RACGAP1 in pan-cancer and its potential value in nasopharyngeal carcinoma

A particular GTPase-activating protein called RACGAP1 is involved in apoptosis, proliferation, invasion, metastasis, and drug resistance in a variety of malignancies. Nevertheless, the role of RACGAP1 in pan-cancer was less studied, and its value of the expression and prognostic of nasopharyngeal carcinoma (NPC) has not been explored. Hence, the goal of this study was to investigate the oncogenic and immunological roles of RACGAP1 in various cancers and its potential value in NPC. We comprehensively analyzed RACGAP1 expression, prognostic value, function, methylation levels, relationship with immune cells, immune infiltration, and immunotherapy response in pan-cancer utilizing multiple databases. The results discovered that RACGAP1 expression was elevated in most cancers and suggested poor prognosis, which could be related to the involvement of RACGAP1 in various cancer-related pathways such as the cell cycle and correlated with RACGAP1 methylation levels, immune cell infiltration and reaction to immunotherapy, and chemoresistance. RACGAP1 could inhibit anti-tumor immunity and immunotherapy responses by fostering immune cell infiltration and cytotoxic T lymphocyte dysfunction. Significantly, we validated that RACGAP1 mRNA and protein were highly expressed in NPC. The Gene Expression Omnibus database revealed that elevated RACGAP1 expression was associated with shorter PFS in patients with NPC, and RACGAP1 potentially influenced cell cycle progression, DNA replication, metabolism, and immune-related pathways, resulting in the recurrence and metastasis of NPC. This study indicated that RACGAP1 could be a potential biomarker in pan-cancer and NPC.

Low-density lipoprotein receptor promotes crosstalk between cell stemness and tumor immune microenvironment in breast cancer: a large data-based multi-omics study

BACKGROUND

Tumor cells with stemness in breast cancer might facilitate the immune microenvironment's suppression process and led to anti-tumor immune effects. The primary objective of this study was to identify potential targets to disrupt the communication between cancer cell stemness and the immune microenvironment.

METHODS

In this study, we initially isolated tumor cells with varying degrees of stemness using a spheroid formation assay. Subsequently, we employed RNA-seq and proteomic analyses to identify genes associated with stemness through gene trend analysis. These stemness-related genes were then subjected to pan-cancer analysis to elucidate their functional roles in a broader spectrum of cancer types. RNA-seq data of 3132 patients with breast cancer with clinical data were obtained from public databases. Using the identified stemness genes, we constructed two distinct stemness subtypes, denoted as C1 and C2. We subsequently conducted a comprehensive analysis of the differences between these subtypes using pathway enrichment methodology and immune infiltration algorithms. Furthermore, we identified key immune-related stemness genes by employing lasso regression analysis and a Cox survival regression model. We conducted in vitro experiments to ascertain the regulatory impact of the key gene on cell stemness. Additionally, we utilized immune infiltration analysis and pan-cancer analysis to delineate the functions attributed to this key gene. Lastly, single-cell RNA sequencing (scRNA-seq) was employed to conduct a more comprehensive examination of the key gene's role within the microenvironment.

RESULTS

In our study, we initially identified a set of 65 stemness-related genes in breast cancer cells displaying varying stemness capabilities. Subsequently, through survival analysis, we pinpointed 41 of these stemness genes that held prognostic significance. We observed that the C2 subtype exhibited a higher stemness capacity compared to the C1 subtype and displayed a more aggressive malignancy profile. Further analysis using Lasso-Cox algorithm identified LDLR as a pivotal immune-related stemness gene. It became evident that LDLR played a crucial role in shaping the immune microenvironment. In vitro experiments demonstrated that LDLR regulated the cell stemness of breast cancer. Immune infiltration analysis and pan-cancer analysis determined that LDLR inhibited the proliferation of immune cells and might promote tumor cell progression. Lastly, in our scRNA-seq analysis, we discovered that LDLR exhibited associations with stemness marker genes within breast cancer tissues. Moreover, LDLR demonstrated higher expression levels in tumor cells compared to immune cells, further emphasizing its relevance in the context of breast cancer.

CONCLUSION

LDLR is an important immune stemness gene that regulates cell stemness and enhances the crosstalk between breast cancer cancer cell stemness and tumor immune microenvironment.

Stemness genes and miR-1247-3p expression associate with clinicopathological parameters and prognosis in lung adenocarcinoma

Lung cancer makes up one-fourth of all cancer-related mortality with the highest mortality rate among all cancers. Despite recent scientific advancements in cancer therapeutics, the 5-year survival rate of lung adenocarcinoma (LUAD) cancer patients remains below 15 percent. It has been suggested that the high mortality rate of LUAD is linked to the acquisition of progenitor-like cells with stem-like characteristics that assist the whole tumor in regulating immune cell infiltration. To examine this hypothesis further, this study mined several databases to explore the presence of stemness-related genes and miRNAs in LUAD cancers. We examine their association with immune and accessory cell infiltration rates and patient survival. We found 3 stem cell-related genes, ORC1L, KIF20A, and DLGAP5, present in LUAD that also correlate with changes in immune infiltration rates and reduced patient survival rates. Additionally, the modulation in myeloid-derived suppressor cell (MDSC) infiltration and miRNA hsa-mir-1247-3p mediated targeting of tumor suppressor SLC24A4 and oncogenes RAB3B and HJURP appears to primarily regulate LUAD patient survival. Given these findings, hsa-mir-1247-3p and/or its associated gene targets may offer a promising avenue to enhance patient survivability.

Combination of downregulating FEN1 and PD-1 blockade enhances antitumor activity of CD8+ T cells against HNSCC cells in vitro

BACKGROUND

Programmed cell death ligand 1 (PD-L1) and human leukocyte antigen/major histocompatibility complex (HLA/MHC) are two main kinds of immunophenotypes affecting the susceptibility to anti-PD therapy. Our previous study found that down-regulation of flap endonuclease-1 (FEN1) could not only inhibit PD-L1 expression, but also upregulate HLA expression in head and neck squamous cell carcinoma (HNSCC). We aimed to clarify whether downregulating FEN1 cloud enhance the response to PD-1 blockade, and possible mechanisms in HNSCC in vitro.

METHODS

Differential expression of FEN1 in HNSCC tumor and normal tissues were explored in the TIMER and TISIDB datasets. A HNSCC cells/CD8+ T cells co-culture model was established. HNSCC cell cycle and apoptosis were recorded by flow cytometry. Immune activity markers of granzyme A, granzyme B, and PRF1 expressed in the CD8+ T cells, and IFN-γ, IL-2, and TNF-α secreted in the supernatants were detected by western blot, ELISA, respectively.

RESULTS

FEN1 was highly expressed in HNSCC and associated with low immune infiltration. Downregulating FEN1 could induce HLA class I expression, and inhibit PD-L1 expression in HNSCC cells. Functionally, FEN1 knockdown enhanced the response to αPD-1 mAb by mediating G2/M phase arrest, apoptosis of HNSCC cells. Mechanistically, targeting FEN1 synergized with αPD-1 mAb could reinforce the antitumor response of CD8+ T cells against HNSCC cells, as indicated by increasing granzyme A, granzyme B, and PRF1 expressions, and promoting IFN-γ, IL-2, and TNF-α secretions.

CONCLUSION

These findings might offer a potential combined strategy for patients resistant to anti-PD therapy via combining FEN1 knockdown and PD-1 blockade.

SPAG5 Expression Predicts Poor Prognosis and is Associated With Adverse Immune Infiltration in Lung Adenocarcinomas

Background

Sperm-associated antigen 5 (SPAG5) has been identified as a novel driver oncogene involved in multiple cancers; however, its role in lung adenocarcinoma (LUAD) needs further investigation. Our study aims to elucidate the potential significance of SPAG5 in LUAD prognosis and its implications for the efficacy of immunotherapy.

Methods

In this study, we used bioinformatics analysis and tissue microarray (TMA) staining to examine the potential role of SPAG5 in LUAD survival and response to immunotherapy. We used the Oncomine, TIMER2.0, Gene Expression Profiling Interactive Analysis (GEPIA), Sangerbox, PredicScan, and Kaplan-Meier Plotter databases to examine the expression and prognostic role of SPAG5 in the LUAD of The Cancer Genome Atlas (TCGA), Gene Expression Omnibus (GEO), and other databases. We also used Cancer Single-cell State Atlas (CancerSEA) and Tumor Immune Estimation Resource (TIMER2.0) to analyze the association of SPAG5 with malignant phenotype and tumor immune microenvironment. Furthermore, Immune Cell Abundance Identifier (ImmuCellAI) analysis of TCGA sequencing data was used to predict the role of SPAG5 in determining the response to immune checkpoint blockade (ICB) treatment in LUAD. Co-expression analysis of programmed death-ligand 1 (PD-L1) and SPAG5 was performed using LUAD TMA immunohistochemistry (IHC) analysis.

Results

Our findings indicate that SPAG5 is overexpressed in LUAD and is positively correlated with advanced clinical stage, poor overall survival, relapse-free survival, and progression-free survival outcomes. SPAG5 may be involved in regulating the cell cycle, proliferation, invasion, DNA damage and repair, and tumor immunosuppression. Furthermore, TMA IHC analysis showed a positive correlation between PD-L1 expression in LUAD and SPAG5 which suggests that SPAG5 may serve as a potential predictor of response to ICB therapy in LUAD.

Conclusions

Our results highlight the role of SAPG5 in promoting a tumor malignancy phenotype and immunosuppression in LUAD and suggest that SPAG5 may serve as a potential response marker for ICB therapy.

CCNB1 is a novel prognostic biomarker and promotes proliferation, migration and invasion in Wilms tumor

BACKGROUND

Wilms tumour (WT) is a mixed type of embryonal tumour that usually occurs in early childhood. However, our knowledge of the pathogenesis or progression mechanism of WT is inadequate, and there is a scarcity of beneficial therapeutic strategies.

METHODS

High-throughput RNA sequencing was employed in this study to identify differentially expressed genes (DEGs) in clinical tumor samples and matching normal tissues. The STRING database was utilized to build a protein-protein interaction (PPI) network, and the Cytohubba method was used to identify the top 10 highly related HUB genes. Then, the key genes were further screened by univariate COX survival analysis. Subsequently, the XCELL algorithm was used to evaluate the tumour immune infiltration. RT-PCR, WB, and IF were used to verify the expression level of key genes in clinical tissues and tumour cell lines. Finally, the function of the key gene was further verified by loss-of-function experiments.

RESULTS

We initially screened 1612 DEGs, of which 1030 were up-regulated and 582 were down-regulated. The GO and KEGG enrichment analysis suggested these genes were associated with 'cell cycle', 'DNA replication'. Subsequently, we identified 10 key HUB genes, among them CCNB1 was strongly related to WT patients' overall survival. Multiple survival analyses showed that CCNB1 was an independent indicator of WT prognosis. Thus, we constructed a nomogram of CCNB1 combined with other clinical indicators. Single gene GSEA and immune infiltration analysis revealed that CCNB1 was associated with the degree of infiltration or activation status of multiple immune cells. TIDE analysis indicated that this gene was correlated with multiple key immune checkpoint molecules and TIDE scores. Finally, we validated the differential expression level of CCNB1 in an external gene set, the pan-cancer, clinical samples, and cell lines. CCNB1 silencing significantly inhibited the proliferation, migration, and invasive capabilities of WIT-49 cells, also, promoted apoptosis, and in turn induced G2 phase cell cycle arrest in loss-of-function assays.

CONCLUSION

Our study suggests that CCNB1 is closely related to WT progression and prognosis, and serves as a potential target.

AURKA, TOP2A and MELK are the key genes identified by WGCNA for the pathogenesis of lung adenocarcinoma

The comprehensive analysis of single or multiple microarray datasets is currently available in Gene Expression Omnibus (GEO) databases, with several studies having identified genes strongly associated with the development of lung adenocarcinoma (LUAD). However, the mechanisms of LUAD development remain largely unknown and has not yet been systematically studied; thus, further studies are required in this field. In the present study, weighted gene co-expression network analysis (WGCNA) was used for the evaluation of key genes with potential high risk of LUAD, and to provide more reliable evidence concerning its pathogenesis. The GSE140797 dataset from the high-throughput GEO database was downloaded and was first analyzed using the Limma package in the R language in order to determine the differentially expressed genes. The dataset was then analyzed using the WGCNA package to analyze the co-expressed genes, and the modular genes with the highest correlation with the clinical phenotype were identified. Subsequently, the pathogenic genes shared in common between the result of the two analyses were imported into the STRING database for protein-protein interaction network analysis. The hub genes were screened out using Cytoscape, and then The Cancer Genome Atlas analysis, receiver operating characteristic analysis and survival analysis were subsequently performed. Finally, the key genes were evaluated using reverse transcription-quantitative PCR and western blot analysis. Bioinformatics analysis of the GSE140797 dataset revealed eight key genes: AURKA, BUB1, CCNB1, CDK1, MELK, NUSAP1, TOP2A and PBK. Finally, the AURKA, TOP2A and MELK genes were evaluated in samples from patients with lung cancer using WGCNA and RT-qPCR, western blot analysis experiments, providing basis for further research on the mechanisms of LUAD development and targeted therapy.

KIF2C accelerates the development of non-small cell lung cancer and is suppressed by miR-186-3p via the AKT-GSK3β-β-catenin pathway

This study aimed to explore how kinesin family member 2C (KIF2C) influences the progression of non-small cell lung cancer (NSCLC). The levels of KIF2C and microRNA-186-3p (miR-186-3p) were examined by quantitative real-time polymerase chain reaction (qRT-PCR). Through the utilization of cell counting kit-8 (CCK-8) assay, colony formation assay, wound closure assay, and Transwell assay, NSCLC cell proliferation, migration, and invasion were identified, respectively. NSCLC cell apoptosis was assessed using the TUNEL assay and flow cytometry (FCM) assay. Luciferase reporter analysis was used to investigate the relationship between KIF2C and miR-186-3p. Western blot assays were conducted to investigate the influence of KIF2C on the AKT-GSK3β-β-catenin pathway. The results showed that KIF2C was up-regulated in NSCLC cells, which predicted poor prognosis. KIF2C overexpression promoted the proliferation, migration, and invasion of NSCLC cells as well as inhibited NSCLC cell apoptosis. KIF2C was as a key target of miR-186-3p. High expression of KIF2C, meanwhile, increased the levels of β-catenin, p-GSK-3β and phosphorylated protein kinase B (p-AKT). KIF2C downregulation and miR-186-3p upregulation reversed these outcomes. As an oncogenic factor, KIF2C is negatively regulated by miR-186-3p and participates in the progression of NSCLC through the AKT-GSK3β-β-catenin pathway.

Enhancement of immune surveillance in breast cancer by targeting hypoxic tumor endothelium: Can it be an immunological switch point?

Breast cancer ranks second among the causes of cancer-related deaths in women. In spite of the recent advances achieved in the diagnosis and treatment of breast cancer, further study is required to overcome the risk of cancer resistance to treatment and thereby improve the prognosis of individuals with advanced-stage breast cancer. The existence of a hypoxic microenvironment is a well-known event in the development of mutagenesis and rapid proliferation of cancer cells. Tumor cells, purposefully cause local hypoxia in order to induce angiogenesis and growth factors that promote tumor growth and metastatic characteristics, while healthy tissue surrounding the tumor suffers damage or mutate. It has been found that these settings with low oxygen levels cause immunosuppression and a lack of immune surveillance by reducing the activation and recruitment of tumor infiltrating leukocytes (TILs). The immune system is further suppressed by hypoxic tumor endothelium through a variety of ways, which creates an immunosuppressive milieu in the tumor microenvironment. Non responsiveness of tumor endothelium to inflammatory signals or endothelial anergy exclude effector T cells from the tumor milieu. Expression of endothelial specific antigens and immunoinhibitory molecules like Programmed death ligand 1,2 (PDL-1, 2) and T cell immunoglobulin and mucin-domain containing-3 (TIM-3) by tumor endothelium adds fuel to the fire by inhibiting T lymphocytes while promoting regulatory T cells. The hypoxic microenvironment in turn recruits Myeloid Derived Suppressor Cells (MDSCs), Tumor Associated Macrophages (TAMs) and T regulatory cells (Treg). The structure and function of newly generated blood vessels within tumors, on the other hand, are aberrant, lacking the specific organization of normal tissue vasculature. Vascular normalisation may work for a variety of tumour types and show to be an advantageous complement to immunotherapy for improving tumour access. By enhancing immune response in the hypoxic tumor microenvironment, via immune-herbal therapeutic and immune-nutraceuticals based approaches that leverage immunological evasion of tumor, will be briefly reviewed in this article. Whether these tactics may be the game changer for emerging immunological switch point to attenuate the breast cancer growth and prevent metastatic cell division, is the key concern of the current study.

NUSAP1 Binds ILF2 to Modulate R-Loop Accumulation and DNA Damage in Prostate Cancer

Increased expression of NUSAP1 has been identified as a robust prognostic biomarker in prostate cancer and other malignancies. We have previously shown that NUSAP1 is positively regulated by E2F1 and promotes cancer invasion and metastasis. To further understand the biological function of NUSAP1, we used affinity purification and mass spectrometry proteomic analysis to identify NUSAP1 interactors. We identified 85 unique proteins in the NUSAP1 interactome, including ILF2, DHX9, and other RNA-binding proteins. Using proteomic approaches, we uncovered a function for NUSAP1 in maintaining R-loops and in DNA damage response through its interaction with ILF2. Co-immunoprecipitation and colocalization using confocal microscopy verified the interactions of NUSAP1 with ILF2 and DHX9, and RNA/DNA hybrids. We showed that the microtubule and charged helical domains of NUSAP1 were necessary for the protein-protein interactions. Depletion of ILF2 alone further increased camptothecin-induced R-loop accumulation and DNA damage, and NUSAP1 depletion abolished this effect. In human prostate adenocarcinoma, NUSAP1 and ILF2 mRNA expression levels are positively correlated, elevated, and associated with poor clinical outcomes. Our study identifies a novel role for NUSAP1 in regulating R-loop formation and accumulation in response to DNA damage through its interactions with ILF2 and hence provides a potential therapeutic target.

Transmission of Exosomal TPX2 Promotes Metastasis and Resistance of NSCLC Cells to Docetaxel

Background

Lung cancer, most of which is non-small cell lung cancer (NSCLC), is the most common tumor in the world, and drug resistance, as a major problem in clinical treatment, has attracted extensive attention. However, the role and mechanism of Targeting protein for Xenopus kinesin-like protein 2 (TPX2), which is highly expressed in NSCLC, is still unclear.

Methods

Bioinformatics analysis was used to analyze the relationship between TPX2 and the clinicopathological features of NSCLC. Stable TPX2 overexpression cell lines with were constructed by lentivirus infection, and the effect of TPX2 on proliferation, migration, invasion and chemoresistance to docetaxel was characterized by the CCK8, wound healing, transwell, colony formation assay and FACS. An in vivo lung homing mouse model was used to further confirmed the role of TPX2 on metastasis. Exosomes were extracted by differential centrifugation from the culture supernatant, and their functions were investigated by co-culture with tumor cells. Gene expression was detected via Western blot and real time PCR (RT-qPCR).

Results

Overexpression of TPX2 was related to the poor prognosis of NSCLC. Promoted migration, invasion and metastasis, and reduced the sensitivity of NSCLC cells to docetaxel. The abundance of TPX2 can be packaged in vesicles and transported to other cells. In addition, overexpression of TPX2 induced the accumulation of β-catenin and C-myc.

Conclusion

Our findings indicated that intercellular transfer of exosomal TPX2 triggered metastasis and resistance against to docetaxel in lung cancer cells, through activating downstream WNT/β-catenin signaling pathway.

The Epitope Basis of Embryonic Stem Cell-Induced Antitumor Immunity against Bladder Cancer

Embryonic stem cells (ESCs) and induced pluripotent stem cells (iPSCs) share many cellular and molecular features with cancer cells. Taking advantage of these similarities, stem cells are effective vaccines against cancers in animal models. However, the molecular basis is not well understood, which hinders the development of effective cancer vaccines. Here, prophylactic and therapeutic bladder cancer vaccines composed of allogeneic ESCs and CpG with or without granulocyte macrophage colony stimulating factor are tested. The ESC-based cancer vaccines are able to induce specific antitumor immunity including stimulating cytotoxic CD8⁺ T cells and memory CD4⁺ T cells, reducing myeloid-derived suppressor cells, and preventing bladder cancer growth in mouse models. Furthermore, several genes that are overexpressed in both ESCs and tumors are identified. An epitope-based vaccine designed with shared overexpressed proteins induces specific antitumor immunity and reduces bladder cancer growth. Functional epitopes underlying the action of stem cell-based vaccines against bladder cancer are identified and it is confirmed that ESC-based anticancer vaccines have great potential. A systematic approach is provided here to developing novel effective epitope-based cancer vaccines in the future.

NUSAP1, a novel stemness-related protein, promotes early recurrence of hepatocellular carcinoma

Early recurrence (within 2 years after resection) is the primary cause of poor outcomes among hepatocellular carcinoma (HCC) patients, and liver cancer stem cells are the main contributors to postsurgical HCC recurrence. Nucleolar and spindle-associated protein 1 (NUSAP1) has been reported to be involved in tumor progression. We investigated the function and clinical value of NUSAP1 in early recurrence of HCC. Data from public datasets and our cohort were used to assess the association between NUSAP1 expression and early HCC recurrence. Gain- and loss-of-function experiments were carried out in vivo and in vitro. The predictive effect of NUSAP1 on early HCC recurrence was further evaluated by a validation cohort. We found that elevated NUSAP1 expression in HCC specimens was correlated with poor outcome, especially in cases with postoperative early recurrence. Functional studies indicated that NUSAP1 significantly promotes HCC progression. A postsurgical recurrence murine model further revealed that upregulated NUSAP1 dramatically increased the likelihood of HCC early recurrence. RNA sequencing data revealed that the gene sets of cancer stemness and the signal transducer and activator of transcription 3 (STAT3) pathway were enriched by NUSAP1 overexpression. Mechanistically, NUSAP1 enhanced cancer stemness through stimulating STAT3 nuclear translocation and activation through receptor of activated protein C kinase 1 (RACK1). In a validation cohort with 112 HCC patients, NUSAP1 effectively predicted HCC early recurrence. Our results indicated that NUSAP1 promotes early recurrence of HCC by sustaining cancer stemness and could serve as a valuable predictive indicator for postsurgical intervention in HCC patients.

MAGE-A3 regulates tumor stemness in gastric cancer through the PI3K/AKT pathway

Gastric cancer remains a malignant disease of the digestive tract with high mortality and morbidity worldwide. However, due to its complex pathological mechanisms and lack of effective clinical therapies, the survival rate of patients after receiving treatment is not satisfactory. A increasing number of studies have focused on cancer stem cells and their regulatory properties. In this study, we first constructed a co-expression network based on the WGCNA algorithm to identify modules with different degrees of association with tumor stemness indices. After selecting the most positively correlated modules of the stemness index, we performed a consensus clustering analysis on gastric cancer samples and constructed the co-expression network again. We then selected the modules of interest and applied univariate COX regression analysis to the genes in this module for preliminary screening. The results of the screening were then used in LASSO regression analysis to construct a risk prognostic model and subsequently a sixteen-gene model was obtained. Finally, after verifying the accuracy of the module and screening for risk genes, we identified MAGE-A3 as the final study subject. We then performed in vivo and in vitro experiments to verify its effect on tumor stemness and tumour proliferation. Our data supports that MAGE-A3 is a tumor stemness regulator and a potent prognostic biomarker which can help the prediction and treatment of gastric cancer patients.

Silencing of CPSF7 inhibits the proliferation, migration, and invasion of lung adenocarcinoma cells by blocking the AKT/mTOR signaling pathway

Cleavage and polyadenylation specific factor 7 (CPSF7) is an important participator in the cleavage and polyadenylation of pre-mRNAs. This study aims to uncover the function and underlying mechanism of CPSF7 in lung adenocarcinoma (LUAD). CPSF7 expression in LUAD cells was measured using real time-quantitative polymerase chain reaction and Western blotting. Our results showed that CPSF7 expression was upregulated in LUAD cell lines (A549, H1299, and HCC827). To explore the function of CPSF7 on LUAD, CPSF7 was silenced by the si-CPSF7 transfection and overexpressed by the oe-CPSF7 transfection in A549 cells. Cell proliferation was measured using cell counting kit-8 and colony formation assays. Cell migration and invasion were measured by wound healing and Transwell assays, respectively. Our data revealed that CPSF7 silencing inhibited the viability, colony formation, migration, and invasion of LUAD cells. On the contrary, CPSF7 overexpression enhanced the malignant characteristics of LUAD cells. Additionally, expression of AKT/mTOR pathway-related proteins was detected using Western blotting. CPSF7 silencing blocked the AKT/mTOR signaling pathway. The intervention of SC79 (an activator of the AKT/mTOR pathway) weakened the antitumor effects of CPSF7 silencing in LUAD cells. Silencing of CPSF7 inhibits the malignant characteristics of LUAD cells by blocking the AKT/mTOR signaling pathway.

Development and Optimization of a Prognostic Model Associated with Stemness Genes in Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) is one of the most lethal cancers worldwide, which is associated with a variety of risk factors. Cancer stem cells are self-renewal cells, which can promote the occurrence and metastasis of tumors and enhance the drug resistance of tumor treatment. This study aimed to develop a stemness score model to assess the prognosis of hepatocellular carcinoma (HCC) patients for the optimization of treatment. The single-cell sequencing data GSE149614 was downloaded from the GEO database. Then, we compared the gene expression of hepatic stem cells and other hepatocytes in tumor samples to screen differentially expressed genes related to stemness. R package "clusterProfiler" was used to explore the potential function of stemness-related genes. We then constructed a prognostic model using LASSO regression analysis based on the TCGA and GSE14520 cohorts. The associations of stemness score with clinical features, drug sensitivity, gene mutation, and tumor immune microenvironment were further explored. R package "rms" was used to construct the nomogram model. A total of 18 stemness-related genes were enrolled to construct the prognosis model. Kaplan-Meier analysis proved the good performance of the stemness score model at predicting overall survival (OS) of HCC patients. The stemness score was closely associated with clinical features, drug sensitivity, and tumor immune microenvironment of HCC. The infiltration level of CD8+ T cells was lower, and tumor-associated macrophages were higher in patients with high-stemness score, indicating an immunosuppressive microenvironment. Our study established an 18 stemness-related gene model that reliably predicts OS in HCC. The findings may help clarify the biological characteristics and progression of HCC and help the future diagnosis and therapy of HCC.

Integrative Analysis of Bulk RNA-Seq and Single-Cell RNA-Seq Unveils the Characteristics of the Immune Microenvironment and Prognosis Signature in Prostate Cancer

The immune microenvironment is a culmination of the collaborative effort of immune cells and is important in cancer development. The underlying mechanisms of the tumor immune microenvironment in regulating prostate cancer (PRAD) are unclear. In the current study, 144 natural killer cell-related genes were identified using differential expression, single-sample gene set enrichment analysis, and weighted gene coexpression network analysis. Furthermore, VCL, ACTA2, MYL9, MYLK, MYH11, TPM1, ACTG2, TAGLN, and FLNC were selected as hub genes via the protein-protein interaction network. Based on the expression patterns of the hub genes, endothelial, epithelial, and tissue stem cells were identified as key cell subpopulations, which could regulate PRAD via immune response, extracellular signaling, and protein formation. Moreover, 27 genes were identified as prognostic signatures and used to construct the risk score model. Receiver operating characteristic curves revealed the good performance of the risk score model in both the training and testing datasets. Different chemotherapeutic responses were observed between the low- and high-risk groups. Additionally, a nomogram based on the risk score and other clinical features was established to predict the 1-, 3-, and 5-year progression-free interval of patients with PRAD. This study provides novel insights into the molecular mechanisms of the immune microenvironment and its role in the pathogenesis of PARD. The identification of key cell subpopulations has a potential therapeutic and prognostic use in PRAD.

Drug-Tolerant Persister Cells in Cancer Therapy Resistance

One of the current stumbling blocks in our fight against cancer is the development of acquired resistance to therapy, which is attributable to approximately 90% of cancer-related deaths. Undercutting this process during treatment could significantly improve cancer management. In many cases, drug resistance is mediated by a drug-tolerant persister (DTP) cell subpopulation present in tumors, often referred to as persister cells. This review provides a summary of currently known persister cell subpopulations and approaches to target them. A specific DTP cell subpopulation with elevated levels of aldehyde dehydrogenase (ALDH) activity has stem cell-like characteristics and a high level of plasticity, enabling them to switch rapidly between high and low ALDH activity. Further studies are required to fully elucidate the functions of ALDH-high DTP cells, how they withstand drug concentrations that kill other cells, and how they rapidly adapt under levels of high cellular stress and eventually lead to more aggressive, recurrent, and drug-resistant cancer. Furthermore, this review addresses the processes used by the ALDH-high persister cell subpopulation to enable cancer progression, the ALDH isoforms important in these processes, interactions of ALDH-high DTPs with the tumor microenvironment, and approaches to therapeutically modulate this subpopulation in order to more effectively manage cancer.

CDC20 is a novel biomarker for improved clinical predictions in epithelial ovarian cancer

Epithelial ovarian cancer (EOC), a common tumor of the female reproductive system, ranks first in fatalities among gynecological malignancies. Most patients find tumors at late stage and have extremely poor prognoses, which necessitates improvements in early detection. This study applied bioinformatic methods to identify potential biomarkers of EOC. First, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed on differentially expressed genes (DEGs) and hub genes, and a protein-protein interaction (PPI) network was constructed. The network of hub genes was analyzed using GeneMANIA, and an analysis of biological processes was constructed with BINGO. Lastly, hub genes were analyzed for EOC-related oncology using the Oncomine and TCGA databases, and the cBioPortal online platform. Overall, cell division cycle 20 (CDC20) was identified as a key gene in EOC. Short hairpin RNA (shRNA) was used to silence CDC20 to explore its effects on EOC cell proliferation, apoptosis and SRY-related HMG-box 2 (SOX2) expression. DEGs were enriched in pathways related to cell cycle signaling, cancer, progesterone-mediated oocyte maturation, Wnt signaling and P53 signaling. Analysis revealed high expression of CDC20 in EOC tissues and a correlation with histology and tumor grade. CDC20 levels are highest in serous adenocarcinoma, when compared to ovarian clear cell carcinoma, ovarian endometrioid carcinoma and mucinous adenocarcinoma. High CDC20 expression within the tumor is associated with poor EOC prognosis. After silencing CDC20, EOC cell proliferation and migration decreased, apoptosis increased, and SOX2 expression decreased. In conclusion, CDC20 is likely a key biomarker of EOC and may act as an upstream regulator of SOX2 to mediate the SOX2 signaling in the progression of EOC. Future application of CDC20 analysis to early detection may improve prognosis, and it has the potential to be a therapeutic target.

Immune- and Stemness-Related Genes Revealed by Comprehensive Analysis and Validation for Cancer Immunity and Prognosis and Its Nomogram in Lung Adenocarcinoma

Objective

Lung adenocarcinoma (LUAD) is a familiar lung cancer with a very poor prognosis. This study investigated the immune- and stemness-related genes to develop model related with cancer immunity and prognosis in LUAD.

Method

The Cancer Genome Atlas (TCGA) was utilized for obtaining original transcriptome data and clinical information. Differential expression, prognostic value, and correlation with clinic parameter of mRNA stemness index (mRNAsi) were conducted in LUAD. Significant mRNAsi-related module and hub genes were screened using weighted gene coexpression network analysis (WGCNA). Meanwhile, immune-related differential genes (IRGs) were screened in LUAD. Stem cell index and immune-related differential genes (SC-IRGs) were screened and further developed to construct prognosis-related model and nomogram. Comprehensive analysis of hub genes and subgroups, involving enrichment in the subgroup [gene set enrichment analysis (GSEA)], gene mutation, genetic correlation, gene expression, immune, tumor mutation burden (TMB), and drug sensitivity, used bioinformatics and reverse transcription polymerase chain reaction (RT-PCR) for verification.

Results

Through difference analysis, mRNAsi of LUAD group was markedly higher than that of normal group. Clinical parameters (age, gender, and T staging) were ascertained to be highly relevant to mRNAsi. MEturquoise and MEblue were found to be the most significant modules (including positive and negative correlations) related to mRNAsi via WGCNA. The functions and pathways of the two mRNAsi-related modules were mainly enriched in tumorigenesis, development, and metastasis. Combining stem cell index–related differential genes and immune-related differential genes, 30 prognosis-related SC-IRGs were screened via Cox regression analysis. Then, 16 prognosis-related SC-IRGs were screened to construct a LASSO regression model at last. In addition, the model was successfully validated by using TCGA-LUAD and GSE68465, whereas c-index and the calibration curves were utilized to demonstrate the clinical value of our nomogram. Following the validation of the model, GSEA, immune cell correlation, TMB, clinical relevance, etc., have found significant difference in high- and low-risk groups, and 16-gene expression of the SC-IRG model also was tested by RT-PCR. ADRB2, ANGPTL4, BDNF, CBLC, CX3CR1, and IL3RA were found markedly different expression between the tumor and normal group.

Conclusion

The SC-IRG model and the prognostic nomogram could accurately predict LUAD survival. Our study used mRNAsi combined with immunity that may lay a foundation for the future research studies in LUAD.

Immune Landscape and Classification in Lung Adenocarcinoma Based on a Novel Cell Cycle Checkpoints Related Signature for Predicting Prognosis and Therapeutic Response

Lung adenocarcinoma (LUAD) is one of the most common malignancies with the highest mortality globally, and it has a poor prognosis. Cell cycle checkpoints play a central role in the entire system of monitoring cell cycle processes, by regulating the signalling pathway of the cell cycle. Cell cycle checkpoints related genes (CCCRGs) have potential utility in predicting survival, and response to immunotherapies and chemotherapies. To examine this, based on CCCRGs, we identified two lung adenocarcinoma subtypes, called cluster1 and cluster2, by consensus clustering. Enrichment analysis revealed significant discrepancies between the two subtypes in gene sets associated with cell cycle activation and tumor progression. In addition, based on Least Absolute Shrinkage and Selection Operator (LASSO) Cox regression, we have developed and validated a cell cycle checkpoints-related risk signature to predict prognosis, tumour immune microenvironment: (TIME), immunotherapy and chemotherapy responses for lung adenocarcinoma patients. Results from calibration plot, decision curve analysis (DCA), and time-dependent receiver operating characteristic curve (ROC) revealed that combining age, gender, pathological stages, and risk score in lung adenocarcinoma patients allowed for a more accurate and predictive nomogram. The area under curve for lung adenocarcinoma patients with 1-, 3-, 5-, and 10-year overall survival was: 0.74, 0.73, 0.75, and 0.81, respectively. Taken together, our proposed 4-CCCRG signature can serve as a clinically useful indicator to help predict patients outcomes, and could provide important guidance for immunotherapies and chemotherapies decision for lung adenocarcinoma patients.

Chromobox 4 (CBX4) promotes tumor progression and stemness via activating CDC20 in gastric cancer

BACKGROUND

The Chromobox homolog 4 (CBX4) has been found to be overexpressed in multiple malignancies. However, the associations between CBX4 and gastric cancer (GC) have remained unclear. This study aimed to determine the biological roles of CBX4 in GC and identify effective therapeutic targets.

METHODS

The 3-(4,5-dimethylthiazol-2-yl) (MTT) assays were used to screen CBX family members. Differential analysis was utilized to evaluate the CBX4 levels. Kaplan-Meier analysis was used to perform prognostic analysis. Western blotting assay, quantitative polymerase chain reaction (qPCR) assay and immunohistochemistry (IHC) were used to assess CBX4 expressions. Colony formation assay, Cell Counting Kit-8 (CCK-8) assay, and Transwell assay were used to assess progression features of cells. The tail vein injection model was utilized to determine the metastatic efficacy of GC cells. Tumor sphere formation assay was used to assess tumor stemness maintenance ability. Chromatin immunoprecipitation (ChIP)-qPCR assay was used to evaluate the associations between CBX4 and CDC20. A subcutaneous tumor model was used to assess the in vivo growth ability of GC.

RESULTS

The MTT assay revealed that only CBX4 inhibition could lead to notable restriction of GC growth, as compared to others. Differential analysis suggested that CBX4 was upregulated in tumor samples relative to normal tissues. Less favorable overall survival (OS) outcomes were noticed in GC patients with high CBX4 in comparison to those with low CBX4. High CBX4 could notably enhance cell proliferation capacity, migration ability, and in vivo metastatic efficacy. Gene set enrichment analysis (GSEA) indicated the relationships between CBX4 and GC stemness, and CBX4 overexpression could remarkably elevate self-renewal ability of GC cells. In addition, CBX4 could mainly promote CDC20 messenger RNA (mRNA) levels, and targeting CBX4 suppressed the relative CDC20 levels. The ChIP-qPCR assay further demonstrated that CBX4 coordinated with H3K4me3 to bind at the CDC20 promoter region. Additionally, CBX4 depended on CDC20 to drive GC growth. Lastly, downregulated CBX4 could notably inhibit the growth of GC in vivo.

CONCLUSIONS

This study highlights the oncogenic roles of CBX4 in GC. CBX4 activates CDC20 to maintain stemness features of GC, thereby creating therapeutic vulnerabilities in the treatment of GC.

An immune gene signature to predict prognosis and immunotherapeutic response in lung adenocarcinoma

Lung adenocarcinoma is one of the most common malignant tumors worldwide. The purpose of this study was to construct a stable immune gene signature for prediction of prognosis (IGSPP) and response to immune checkpoint inhibitors (ICIs) therapy in LUAD patients. Five genes were screened by weighted gene coexpression network analysis, Cox regression and LASSO regression analyses and were used to construct the IGSPP. The survival rate of the IGSPP low-risk group was higher than that of the IGSPP high-risk group. Multivariate Cox regression analysis showed that IGSPP could be used as an independent prognostic factor for the overall survival of LUAD patients. IGSPP genes were enriched in cell cycle pathways. IGSPP gene mutation rates were higher in the high-risk group. CD4 memory-activated T cells, M0 and M1 macrophages had higher infiltration abundance in the high-risk group, which was associated with poor overall survival. In contrast, the abundance of resting CD4 memory T cells, monocytes, resting dendritic cells and resting mast cells associated with a better prognosis was higher in the low-risk group. TIDE scores and the expressions of different immune checkpoints showed that patients in the high-risk IGSPP group benefited more from ICIs treatment. In short, an IGSPP of LUAD was constructed and characterized. It could be used to predict the prognosis and benefits of ICIs treatment in LUAD patients.

The Predictive Competing Endogenous RNA Regulatory Networks and Potential Prognostic and Immunological Roles of Cyclin A2 in Pan-Cancer Analysis

Although accumulating evidence has verified the relationship between CCNA2 and cancers, no pan-cancer analysis about the function and the upstream molecular mechanism of CCNA2 is available. For the first time, we analyzed potential oncogenic roles of CCNA2 in 33 cancer types via The Cancer Genome Atlas (TCGA) database. Overexpression of CCNA2 is widespread in almost all cancer types, and it is related to poor prognosis and advanced pathological stages in most cases. Moreover, we conducted upstream miRNAs and lncRNAs of CCNA2 to establish upstream regulatory networks in kidney renal clear cell carcinoma (LINC00997/miR-27b-3p/CCNA2), liver hepatocellular carcinoma (SNHG16, GUSBP11, FGD5-AS1, LINC00630, CD27-AS1, LINC00997/miR-22-3p/CCNA2, miR-29b-3p/CCNA2, miR-29c-3p/CCNA2, and miR-204-5p/CCNA2), and lung adenocarcinoma (miRNA-218-5p/CCNA2 and miR-204-5p/CCNA2) by expression analysis, survival analysis, and correlation analysis. The CCNA2 expression is positively correlated with Th2 cell infiltration and negatively correlated with CD4+ central memory and effector memory T-cell infiltration in different cancer types. Furthermore, CCNA2 is positively associated with expressions of immune checkpoints (CD274, CTLA4, HAVCR2, LAG3, PDCD1, and TIGIT) in most cancer types. Our first CCNA2 pan-cancer study contributes to understanding the prognostic and immunological roles and potential upstream molecular mechanisms of CCNA2 in different cancers.

Potential molecular target screening and bioinformatics analysis of cholangiocarcinoma based on GEO database

BACKGROUND

Cholangiocarcinoma is a highly malignant tumor with a poor prognosis. Targeted therapy is important for the treatment of cholangiocarcinoma, and it is therefore of great clinical importance to identify novel molecular targets for targeted therapy of this malignancy.

AIM

To identify potential molecular targets for the treatment of cholangiocarcinoma and identify the key genes involved in cholangiocarcinoma by bioinformatics analysis.

METHODS

We downloaded two sets of cholangiocarcinoma expression profile data from GEO database. GEO2R online analysis tool was used to screen differentially expressed genes in cholangiocarcinoma tumor tissues and normal tissues, and we performed GO enrichment analysis, KEGG pathway analysis, and protein interaction network for differentially expressed genes. We used Cytoscape software to calculate key genes. The GEPIA database was used to verify the expression of hub genes in cholangiocarcinoma tissues.

RESULTS

A total of 158 differentially expressed genes were identified. GO enrichment analysis showed that these differential genes were mainly involved in the cellular response to zinc ion, negative regulation of growth, cell adhesion, metabolic process, and protein homotetramerization. They were enriched in exosomes, extracellular spaces, elastic fibers, and organelle membranes. The main molecular functions are related to heparin binding, cysteine-type endopeptidase inhibitor activity, protein homodimerization activity, receptor binding, and pyridoxal phosphate binding. KEGG pathway analysis showed that differential genes are mainly involved in processes such as mineral absorption, carbon and propanoate metabolism, PPAR signaling pathway, and fatty acid degradation. A protein interaction network diagram was constructed based on the String database, and the CytoHubba plug-in of the Cytoscape software was used to calculate the key genes. The key genes were all up-regulated ones. GEPIA analysis verified that the expression of key genes in cholangiocarcinoma tissues was significantly higher than that in normal tissues.

CONCLUSION

In this study, eight key genes related to cholangiocarcinoma were identified, including NUSAP1, TOP2A, RAD51AP1, MCM4, KIAA0101, CDCA5, TYMS, and ZWINT. These genes provide new ideas for in-depth study of the targeted therapy of cholangiocarcinoma, and are expected to become new molecular therapeutic targets.

Identification of a ZC3H12D-regulated competing endogenous RNA network for prognosis of lung adenocarcinoma at single-cell level

Background

To identify hub genes from the competing endogenous RNA (ceRNA) network of lung adenocarcinoma (LUAD) and to explore their potential functions on prognosis of patients from a single-cell perspective.

Methods

We performed RNA-sequencing of LUAD to construct ceRNA regulatory network, integrating with public databases to identify the vital pathways related to patients’ prognosis and to reveal the expression level of hub genes under different conditions, the functional enrichment of co-expressed genes and their potential immune-related mechanisms.

Results

ZC3H12D-hsa-miR-4443-ENST00000630242 axis was found to be related with LUAD. Lower ZC3H12D expression was significantly associated with shorter overall survival (OS) of patients (HR = 2.007, P < 0.05), and its expression was higher in early-stage patients, including T1 (P < 0.05) and N0 (P < 0.05). Additionally, ZC3H12D expression was higher in immune cells displayed by single-cell RNA-sequencing data, especially in Treg cells of lung cancer and CD8 T cells, B cells and CD4 T cells of LUAD. The functional enrichment analysis showed that the co-expressed genes mainly played a role in lymphocyte activation and cytokine-cytokine receptor interaction. In addition, ZC3H12D was associated with multiple immune cells and immune molecules, including immune checkpoints CTLA4, CD96 and TIGIT.

Conclusion

ZC3H12D-hsa-miR-4443-ENST00000630242 ceRNA network was identified in LUAD. ZC3H12D could affect prognosis of patients by regulating mRNA, miRNA, lncRNA, immune cells and immune molecules. Therefore, it may serve as a vital predictive marker and could be regarded as a potential therapeutic target for LUAD in the future.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-021-08992-1.

Novel evidence of obesity paradox in esophageal adenocarcinoma: perspective on genes that uncouple adiposity from dismal outcomes

Background: Obesity is a strong risk factor for esophageal adenocarcinoma (EAC). Nevertheless, not all the patients with EAC are obesity, and a substantial proportion of obesity patients don't suffer from poor prognoses. The mechanisms behind the “obesity paradox” that uncouple obesity from dismal outcomes in EAC are unclear. This study aimed to explore the “obesity-guarding” genes (OGG) profiles and their prognostic values in patients with EAC.

Methods: Gene expression data and clinical information of patients with EAC were downloaded from The Cancer Genome Atlas (TCGA) database. Enrichment analysis was used to explore the OGG functions and pathways. Cox regression analysis and nomogram model were performed to investigate the OGG prognostic values for overall survival (OS). In addition, relations between OGG and immune cells were assessed by the “CIBERSORT” algorithm and the Tumor IMmune Estimation Resource (TIMER) tool. Finally, the results were experimentally validated in real-world study.

Results: A total of 69 OGG were retrieved, and 17 significantly differentially expressed genes (SDEG) were identified between normal and EAC tissues. Enrichment analysis showed the OGG were enriched in the mitochondrion-related and various receptor pathways. Univariate Cox regression results showed that the MCM6, ATXN2 and CSK were significantly associated with OS (P=0.036, 0.039, 0.046, respectively). Multivariate Cox regression analysis showed MCM6 and CSK were independent prognostic genes for OS (P=0.025, 0.041, respectively). Nomogram demonstrated that the OGG had good predictive abilities for the 1-, 2-, and 3-year OS. Immunity analysis demonstrated that OGG were significantly associated with immune cells (P <0.05). In addition, clinical correlation analysis revealed that the OGG had significant relations with clinical parameters (P <0.05). The experiment results confirmed that the SDEG were significantly different between normal and EAC tissues (P <0.05).

Conclusions: We identified the OGG expression profiles that may uncouple obesity from poor survival in patients with EAC. They have prognostic values in predicting patients' OS, and may be exploited for prognostic biomarkers.

Single-Cell RNA Sequencing in Lung Cancer: Revealing Phenotype Shaping of Stromal Cells in the Microenvironment

The lung tumor microenvironment, which is composed of heterogeneous cell populations, plays an important role in the progression of lung cancer and is closely related to therapeutic efficacy. Increasing evidence has shown that stromal components play a key role in regulating tumor invasion, metastasis and drug resistance. Therefore, a better understanding of stromal components in the tumor microenvironment is helpful for the diagnosis and treatment of lung cancer. Rapid advances in technology have brought our understanding of disease into the genetic era, and single-cell RNA sequencing has enabled us to describe gene expression profiles with unprecedented resolution, enabling quantitative analysis of gene expression at the single-cell level to reveal the correlations among heterogeneity, signaling pathways, drug resistance and microenvironment molding in lung cancer, which is important for the treatment of this disease. In this paper, several common single-cell RNA sequencing methods and their advantages and disadvantages are briefly introduced to provide a reference for selection of suitable methods. Furthermore, we review the latest progress of single-cell RNA sequencing in the study of stromal cells in the lung tumor microenvironment.

Classification of lung adenocarcinoma based on stemness scores in bulk and single cell transcriptomes

•

This study explores tumor stemness based on both bulk tumor and single cell transcriptome datasets.

•

High-stemness tumors are less responsive to immunotherapy and targeted therapy compared with low-stemness tumors.

•

Many high-stemness cells are at the beginning of the cell evolution trajectory, while most low-stemness cells are in the terminal or later phase.

•

The correlations of tumor stemness with intratumor heterogeneity and tumor immunity were in the opposite direction between bulk tumors and single cells.

Tumor stemness is associated with tumor progression and therapy resistance. The recent advances in sequencing, genomics, and computational technologies have facilitated investigation into the tumor stemness cell-like characteristics. We identified subtypes of lung adenocarcinoma (LUAD) in bulk tumors or single cells based on the enrichment scores of 12 stemness signatures by clustering analysis of their transcriptomic profiles. Three stemness subtypes of LUAD were identified: St-H, St-M, and St-L, having high, medium, and low stemness signatures, respectively, consistently in six different datasets. Among the three subtypes, St-H was the most enriched in epithelial-mesenchymal transition, invasion, and metastasis signaling, genomically instable, irresponsive to immunotherapies and targeted therapies, and hence had the worst prognosis. We observed that intratumor heterogeneity was significantly higher in high-stemness than in low-stemness bulk tumors, but significantly lower in high-stemness than in low-stemness single cancer cells. Moreover, tumor immunity was stronger in high-stemness than in low-stemness cancer cells, but weaker in high-stemness than in low-stemness bulk tumors. These differences between bulk tumors and single cancer cells could be attributed to the non-tumor cells in bulk tumors that confounded the results of correlation analysis. Furthermore, pseudotime analysis showed that many St-H cells were at the beginning of the cell evolution trajectory, compared to most St-L cells in the terminal or later phase, suggesting that many low-stemness cells are originated from high-stemness cells. The stemness-based classification of LUAD may provide novel insights into the tumor biology as well as precise clinical management of this disease.

System Analysis of ROS-Related Genes in the Prognosis, Immune Infiltration, and Drug Sensitivity in Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) is an aggressive malignant tumor with a poor prognosis. Reactive oxygen species (ROS) play an important role in tumors; however, the role of ROS-related genes is still unclear in HCC. Therefore, we analyzed the role of ROS-related genes in HCC via bioinformatics methods. Firstly, a prognosis model was constructed using LASSO Cox regression and multivariate analyses. We also investigated the potential function of the ROS-related genes and the correlation with immune infiltration, tumor stemness, and drug sensitivity. ICGC database was used for validation. Secondly, we further analyzed the role of 11 ROS-related genes in HCC. As a member of ROS gene family, the role of STK25 has remained unclear in HCC. We explored the biological function of STK25 using in vitro experiments. The present study was the first to construct a ROS-related prognostic model in HCC. The correlation of ROS-related genes with immune infiltration, tumor stemness, and drug sensitivity was dissected. Furthermore, we demonstrated that STK25 knockdown could increase the proliferation, migration, and invasion capacity of HCC cells.

Nuclear accumulation of KPNA2 impacts radioresistance through positive regulation of the PLSCR1-STAT1 loop in lung adenocarcinoma

Lung adenocarcinoma (ADC) is the predominant histological type of lung cancer, and radiotherapy is one of the current therapeutic strategies for lung cancer treatment. Unfortunately, biological complexity and cancer heterogeneity contribute to radioresistance development. Karyopherin α2 (KPNA2) is a member of the importin α family that mediates the nucleocytoplasmic transport of cargo proteins. KPNA2 overexpression is observed across cancer tissues of diverse origins. However, the role of KPNA2 in lung cancer radioresistance is unclear. Herein, we demonstrated that high expression of KPNA2 is positively correlated with radioresistance and cancer stem cell (CSC) properties in lung ADC cells. Radioresistant cells exhibited nuclear accumulation of KPNA2 and its cargos (OCT4 and c‐MYC). Additionally, KPNA2 knockdown regulated CSC‐related gene expression in radioresistant cells. Next‐generation sequencing and bioinformatic analysis revealed that STAT1 activation and nuclear phospholipid scramblase 1 (PLSCR1) are involved in KPNA2‐mediated radioresistance. Endogenous PLSCR1 interacting with KPNA2 and PLSCR1 knockdown suppressed the radioresistance induced by KPNA2 expression. Both STAT1 and PLSCR1 were found to be positively correlated with dysregulated KPNA2 in radioresistant cells and ADC tissues. We further demonstrated a potential positive feedback loop between PLSCR1 and STAT1 in radioresistant cells, and this PLSCR1‐STAT1 loop modulates CSC characteristics. In addition, AKT1 knockdown attenuated the nuclear accumulation of KPNA2 in radioresistant lung cancer cells. Our results collectively support a mechanistic understanding of a novel role for KPNA2 in promoting radioresistance in lung ADC cells.

CCNB1 promotes the development of hepatocellular carcinoma by mediating DNA replication in the cell cycle

In our studies, cyclin B1 (CCNB1) mRNA and protein were overexpressed in hepatocellular carcinoma (HCC) tissues compared with non-HCC tissues. Moreover, CCNB1 was overexpressed in the serum of HCC patients. The expression of CCNB1 was associated with several crucial clinicopathologic characteristics, and the HCC patients with overexpressed CCNB1 had worse overall survival outcomes. In the screening of interactional genes, a total of 266 upregulated co-expression genes, which were positively associated with CCNB1, were selected from the datasets, and 67 downregulated co-expression genes, which were negatively associated with CCNB1, were identified. The key genes might be functionally enriched in DNA replication and the cell cycle pathways. CDC20, CCNA2, PLK1, and FTCD were selected for further research because they were highly connected in the protein-protein interaction networks. Upregulated CDC20, CCNA2, and PLK1 and downregulated FTCD might result in undesirable overall survival outcomes for HCC patients. The univariate Cox analysis results showed that CDC20 and PLK1 might be two independent risk factors, while FTCD might be protective in HCC. Therefore, CCNB1 may participate in the cell cycle of HCC by regulating DNA replication, and CCNB1 may provide a direction for the diagnosis of early-stage HCC and targeted HCC therapy.

A Novel Cancer Stemness-Related Signature for Predicting Prognosis in Patients with Colon Adenocarcinoma

Objective

To explore the cancer stemness features and develop a novel cancer stemness-related prognostic signature for colon adenocarcinoma (COAD).

Methods

We downloaded the mRNA expression data and clinical data of COAD from TCGA database and GEO database. Stemness index, mRNAsi, was utilized to investigate cancer stemness features. Weighted gene coexpression network analysis (WGCNA) was used to identify cancer stemness-related genes. Univariate and multivariate Cox regression analyses were applied to construct a prognostic risk cancer stemness-related signature. We then performed internal and external validation. The relationship between cancer stemness and COAD immune microenvironment was investigated.

Results

COAD patients with higher mRNAsi score or EREG-mRNAsi score have significant longer overall survival (OS). We identified 483 differently expressed genes (DEGs) between the high and low mRNAsi score groups. We developed a cancer stemness-related signature using fifteen genes (including RAB31, COL6A3, COL5A2, CCDC80, ADAM12, VGLL3, ECM2, POSTN, DPYSL3, PCDH7, CRISPLD2, COLEC12, NRP2, ISLR, and CCDC8) for prognosis prediction of COAD. Low-risk score was associated with significantly preferable OS in comparison with high-risk score, and the area under the ROC curve (AUC) for OS prediction was 0.705. The prognostic signature was an independent predictor for OS of COAD. Macrophages, mast cells, and T helper cells were the vital infiltration immune cells, and APC costimulation and type II IFN response were the vital immune pathways in COAD.

Conclusions

We developed and validated a novel cancer stemness-related prognostic signature for COAD, which would contribute to understanding of molecular mechanism in COAD.

Integrated Analysis of Stemness-Related LncRNAs Helps Predict the Immunotherapy Responsiveness of Gastric Cancer Patients

The immune microenvironment plays a critical role in tumor biology. As a critical feature of cancers, stemness is acknowledged as a contributor to the development of drug resistance in gastric cancers (GCs). Long non-coding RNAs (lncRNAs) have been revealed to participate in this process. In this study, we aimed to develop a stemness-related lncRNA signature (SRLncSig) with guiding significance for immunotherapy. Three cohorts (TCGA, Zhongshan, and IMvigor210) were enrolled for analysis. A list of stemness-related lncRNAs (SRlncRNAs) was collected by co-expression strategy under the threshold of coefficient value >0.35 and p-value < 0.05. Cox and Lasso regression analysis was further applied to find out the SRlncRNAs with prognosis-predictive value to establish the SRLncSig in the TCGA cohort. IPS and TIDE algorithms were further applied to predict the efficacy of SRLncSig in TCGA and Zhongshan cohorts. IMvigor210 was composed of patients with clinical outcomes of immunotherapy. The results indicated that SRLncSig not only was confirmed as an independent risk factor for GCs but also identified as a robust indicator for immunotherapy. The patient with a lower SRLncSig score was more likely to benefit from immunotherapy, and the results were highly consistent in three cohorts. In conclusion, our study not only could clarify the correlations between stemness and immunotherapy in GC patients but also provided a model to guide the applications of immunotherapy in clinical practice.

Stemness-related LncRNA pair signature for predicting therapy response in gastric cancer

Objective

As a critical feature of cancers, stemness is acknowledged as a contributor to the development of drug resistance in gastric cancer (GC). LncRNAs have been revealed to participate in this process. In this study, we tried to develop a stemness-related lncRNA pair signature as guidance for clinical decisions.

Methods

The analysis was initiated by collecting stemness-related lncRNAs in TCGA cohort. The differentially expressed stemness-related lncRNAs between normal and tumor tissues in GC patients from TCGA datasets were further collected to establish the signature based on Lasso and Cox regression analyses. The predictive efficacy of the signature for chemotherapy and immunotherapy was also tested. The practicality of this signature was also validated by Zhongshan cohort.

Results

A 13-DEsrlncRNA pair-based signature was established. The cutoff point acquired by the AIC algorithm divided the TCGA cohort into high and low risk groups. We found that the low-risk group presented with better survival (Kaplan-Meier analysis, p < 0.001). Cox regression analyse was also conducted to confirm the signature as an independent risk factor for GC {p < 0.001, HR = 1.300, 95% CI (1.231–1.373)]}. As for the practicality of this signature, the IC50 of cytotoxic chemotherapeutics was significantly higher in the high-risk group. The low-risk group also presented with higher immunophenoscore (IPS) in both the “CTLA4+ PD1+” (Mann-Whitney U test, p = 0.019) and “CTLA4- PD1+” (Mann-Whitney U test, p = 0.013) groups, indicating higher sensitivity to immunotherapy. The efficacy of the signature was also validated by Zhongshan cohort.

Conclusions

This study could not only provide a stemness-related lncRNA signature for survival prediction in GC patients but also established a model with predictive potentials for GC patients’ sensitivity to chemotherapy and immunotherapy.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-021-08798-1.

Mechanisms of Progression and Heterogeneity in Multiple Nodules of Lung Adenocarcinoma

Lung cancer remains the leading cause of cancer-related death worldwide. Lung adenocarcinoma (LUAD) is thought to be caused by precursor lesions of atypical adenoma-like hyperplasia and may have extensive in situ growth before infiltration. To explore the relevant factors in heterogeneity and evolution of lung adenocarcinoma subtypes, the authors perform single-cell RNA sequencing (scRNA-seq) on tumor and normal tissue from five multiple nodules' LUAD patients and conduct a thorough gene expression profiling of cancer cells and cells in their microenvironment at single-cell level. This study gives a deep understanding of heterogeneity and evolution in early glandular neoplasia of the lung. This dataset leads to discovery of the changes in the immune microenvironment during the development of LUAD, and the development process from adenocarcinoma in situ (AIS) to invasive adenocarcinoma (IAC). This work sheds light on the direction of early tumor development and whether they are homologous.

[Progress in Single-cell RNA Sequencing of Lung Adenocarcinoma]

肺腺癌（lung adenocarcinoma, LUAD）是临床上肺癌最常见的亚型，是癌症相关死亡最主要的原因之一。过去十几年中，随着薄层计算机断层扫描（computed tomography, CT）广泛用于常规肺癌筛查，影像学上表现为小结节的LUAD发病率显著增高，其发生发展机制复杂，个体预后差异显著。尽管近年来针对LUAD的靶向和免疫疗法取得了重大进展，但肿瘤细胞的耐药性始终未得到有效解决，从而限制了患者获益。随着人类基因组计划的完成，以测序为基本手段的基因组学及转录组学进入临床和科研人员的视野。单细胞测序作为近年来受到高度关注的新型测序手段，与二代测序相比，其能在单细胞水平上对细胞群体进行特异性分析，揭示出每种细胞类型独特的变化，在单细胞水平上对许多异质基质细胞和癌细胞进行较精准地评估，从而揭示了分子成分的复杂性以及与非恶性组织中相应成分的区别。综上，通过单细胞测序深入了解LUAD发生发展机制和肿瘤微环境（tumor microenvironment, TME）的异质性及其耐药性形成机制，从而发现新的治疗靶点是临床医生和基础科学家迫切的需求。本文综合论述了单细胞测序在LUAD中的具体应用和研究进展。

Weighted Gene Co-expression Network Analysis Identifies CALD1 as a Biomarker Related to M2 Macrophages Infiltration in Stage III and IV Mismatch Repair-Proficient Colorectal Carcinoma

Immunotherapy has achieved efficacy for advanced colorectal cancer (CRC) patients with a mismatch-repair-deficient (dMMR) subtype. However, little immunotherapy efficacy was observed in patients with the mismatch repair-proficient (pMMR) subtype, and hence, identifying new immune therapeutic targets is imperative for those patients. In this study, transcriptome data of stage III/IV CRC patients were retrieved from the Gene Expression Omnibus database. The CIBERSORT algorithm was used to quantify immune cellular compositions, and the results revealed that M2 macrophage fractions were higher in pMMR patients as compared with those with the dMMR subtype; moreover, pMMR patients with higher M2 macrophage fractions experienced shorter overall survival (OS). Subsequently, weighted gene co-expression network analysis and protein–protein interaction network analysis identified six hub genes related to M2 macrophage infiltrations in pMMR CRC patients: CALD1, COL6A1, COL1A2, TIMP3, DCN, and SPARC. Univariate and multivariate Cox regression analyses then determined CALD1 as the independent prognostic biomarker for OS. CALD1 was upregulated specifically the in CMS4 CRC subtype, and single-sample Gene Set Enrichment Analysis (ssGSEA) revealed that CALD1 was significantly correlated with angiogenesis and TGF-β signaling gene sets enrichment scores in stage III/IV pMMR CRC samples. The Estimation of STromal and Immune cells in MAlignant Tumor tissues using Expression data (ESTIMATE) algorithm and correlation analysis revealed that CALD1 was significantly associated with multiple immune and stromal components in a tumor microenvironment. In addition, GSEA demonstrated that high expression of CALD1 was significantly correlated with antigen processing and presentation, chemokine signaling, leukocyte transendothelial migration, vascular smooth muscle contraction, cytokine–cytokine receptor interaction, cell adhesion molecules, focal adhesion, MAPK, and TGF-beta signaling pathways. Furthermore, the proliferation, invasion, and migration abilities of cancer cells were suppressed after reducing CALD1 expression in CRC cell lines. Taken together, multiple bioinformatics analyses and cell-level assays demonstrated that CALD1 could serve as a prognostic biomarker and a prospective therapeutic target for stage III/IV pMMR CRCs.

Prognostic Gene Expression, Stemness and Immune Microenvironment in Pediatric Tumors

Simple Summary

Tumors in children and young adults are rare and diagnostically distinct from those occurring in older patients. They frequently arise from developing cells, resembling stem cells, which may explain some of the clinical and biologic differences observed. The aim of this retrospective transcriptome study was to investigate the prognostic landscape, immune tumor microenvironment (TME) and stemness in a cohort of 4068 transcriptomes of such tumors. We find that patients’ prognosis correlates with distinct gene expression patterns similar to adult tumor types. Stemness defined by a computational stemness score (mRNAsi) correlates with clinical and molecular parameters that is distinct for each tumor type. In Wilms tumors that recapitulate normal kidney development microscopically, stemness correlates with distinct patterns of immune cell infiltration by transcriptome analysis and by cell localization in tumor tissue.

Abstract

Pediatric tumors frequently arise from embryonal cells, often displaying a stem cell-like (“small round blue”) morphology in tissue sections. Because recently “stemness” has been associated with a poor immune response in tumors, we investigated the association of prognostic gene expression, stemness and the immune microenvironment systematically using transcriptomes of 4068 tumors occurring mostly at the pediatric and young adult age. While the prognostic landscape of gene expression (PRECOG) and infiltrating immune cell types (CIBERSORT) is similar to that of tumor entities occurring mainly in adults, the patterns are distinct for each diagnostic entity. A high stemness score (mRNAsi) correlates with clinical and morphologic subtype in Wilms tumors, neuroblastomas, synovial sarcomas, atypical teratoid rhabdoid tumors and germ cell tumors. In neuroblastomas, a high mRNAsi is associated with shortened overall survival. In Wilms tumors a high mRNAsi correlates with blastemal morphology, whereas tumors with predominant epithelial or stromal differentiation have a low mRNAsi and a high percentage of M2 type macrophages. This could be validated in Wilms tumor tissue (n = 78). Here, blastemal areas are low in M2 macrophage infiltrates, while nearby stromal differentiated areas contain abundant M2 macrophages, suggesting local microanatomic regulation of the immune response.