Cell-cell communications shape tumor microenvironment and predict clinical outcomes in clear cell renal carcinoma.

BACKGROUND: Cell-cell communications of various cell populations within tumor microenvironment play an essential role in primary tumor growth, metastasis evolution, and immune escape. Nevertheless, comprehensive investigation of cell-cell communications in the ccRCC (Clear cell renal carcinoma) microenvironment and how this interplay affects prognosis still remains limited. METHODS: Intercellular communications were characterized by single-cell data. Firstly, we employed "CellChat" package to characterize intercellular communications across all types of cells in microenvironment in VHL mutated and non-mutated samples from 8 patients, respectively. And pseudotime trajectory analyses were performed with monocle analyses. Finally clinical prognosis and immunotherapy efficacy with different landscapes of intercellular interplay are evaluated by TCGA-KIRC and immunotherapy cohort. RESULTS: Firstly, the VHL phenotype may be related to the intercellular communication landscape. And trajectory analysis reveals the potential relationship of cell-cell communication molecules with T cells and Myeloid cells differentiation. Furthermore, those molecules also correlate with the infiltration of T cells and Myeloid cells. A tumor cluster with highly expressed ligands was defined by quantitative analysis and transcription factor enrichment analysis, which was identified to be pivotal for intercellular communications in tumor microenvironment. Finally, bulk data indicates bulk that different clusters with different intercellular communications have significant predictive value for prognosis and distinguished immunotherapy efficiency. CONCLUSIONS: The intercellular communication landscapes of VHL wild and VHL mutant ccRCC vary. Intercellular communications within the tumor microenvironment also influence T cell and myeloid cell development and infiltration, as well as predict clinical prognosis and immunotherapy efficacy in ccRCC.

Pan-cancer dissection of vasculogenic mimicry characteristic to provide potential therapeutic targets.

Introduction: Vasculogenic mimicry (VM) represents a novel form of tumor angiogenesis that is associated with tumor invasiveness and drug resistance. However, the VM landscape across cancer types remains poorly understood. In this study, we elucidate the characterizations of VM across cancers based on multi-omics data and provide potential targeted therapeutic strategies. Methods: Multi-omics data from The Cancer Genome Atlas was used to conduct comprehensive analyses of the characteristics of VM related genes (VRGs) across cancer types. Pan-cancer vasculogenic mimicry score was established to provide a depiction of the VM landscape across cancer types. The correlation between VM and cancer phenotypes was conducted to explore potential regulatory mechanisms of VM. We further systematically examined the relationship between VM and both tumor immunity and tumor microenvironment (TME). In addition, cell communication analysis based on single-cell transcriptome data was used to investigate the interactions between VM cells and TME. Finally, transcriptional and drug response data from the Genomics of Drug Sensitivity in Cancer database were utilized to identify potential therapeutic targets and drugs. The impact of VM on immunotherapy was also further clarified. Results: Our study revealed that VRGs were dysregulated in tumor and regulated by multiple mechanisms. Then, VM level was found to be heterogeneous among different tumors and correlated with tumor invasiveness, metastatic potential, malignancy, and prognosis. VM was found to be strongly associated with epithelial-mesenchymal transition (EMT). Further analyses revealed cancer-associated fibroblasts can promote EMT and VM formation. Furthermore, the immune-suppressive state is associated with a microenvironment characterized by high levels of VM. VM score can be used as an indicator to predict the effect of immunotherapy. Finally, seven potential drugs targeting VM were identified. Conclusion: In conclusion, we elucidate the characteristics and key regulatory mechanisms of VM across various cancer types, underscoring the pivotal role of CAFs in VM. VM was further found to be associated with the immunosuppressive TME. We also provide clues for the research of drugs targeting VM. Our study provides an initial overview and reference point for future research on VM, opening up new avenues for therapeutic intervention.

Single-cell multi-omics analysis reveals dysfunctional Wnt signaling of spermatogonia in non-obstructive azoospermia.

Background: Non-obstructive azoospermia (NOA) is the most severe type that leads to 1% of male infertility. Wnt signaling governs normal sperm maturation. However, the role of Wnt signaling in spermatogonia in NOA has incompletely been uncovered, and upstream molecules regulating Wnt signaling remain unclear. Methods: Bulk RNA sequencing (RNA-seq) of NOA was used to identify the hub gene module in NOA utilizing weighted gene co-expression network analyses (WGCNAs). Single-cell RNA sequencing (scRNA-seq) of NOA was employed to explore dysfunctional signaling pathways in the specific cell type with gene sets of signaling pathways. Single-cell regulatory network inference and clustering (pySCENIC) for Python analysis was applied to speculate putative transcription factors in spermatogonia. Moreover, single-cell assay for transposase-accessible chromatin sequencing (scATAC-seq) determined the regulated genes of these transcription factors. Finally, spatial transcriptomic data were used to analyze cell type and Wnt signaling spatial distribution. Results: The Wnt signaling pathway was demonstrated to be enriched in the hub gene module of NOA by bulk RNA-seq. Then, scRNA-seq data revealed the downregulated activity and dysfunction of Wnt signaling of spermatogonia in NOA samples. Conjoint analyses of the pySCENIC algorithm and scATAC-seq data indicated that three transcription factors (CTCF, AR, and ARNTL) were related to the activities of Wnt signaling in NOA. Eventually, spatial expression localization of Wnt signaling was identified to be in accordance with the distribution patterns of spermatogonia, Sertoli cells, and Leydig cells. Conclusion: In conclusion, we identified that downregulated Wnt signaling of spermatogonia in NOA and three transcription factors (CTCF, AR, and ARNTL) may be involved in this dysfunctional Wnt signaling. These findings provide new mechanisms for NOA and new therapeutic targets for NOA patients.

A novel classification based on non-apoptosis cell death predicts clinical outcomes and immunotherapy response of clear renal cell carcinoma.

BACKGROUND: Non-apoptosis cell death could be a secondary consequence of the immune response, which profoundly influences tumor microenvironment (TME), escaping from chemotherapy/immunotherapy-induced apoptosis resistance effects. Whereas, systemic analysis of non-apoptosis regulated cell death associated with TME and clinical outcomes remains unveiled. METHODS: Our kidney clear carcinoma (KIRC) samples from The Cancer Genome Atlas (TCGA) were stratified into three clusters based on the activity of autophagic cell death, ferroptosis, pyroptosis and necroptosis. Clinical prognosis, TME landscape, biological functions and somatic mutation frequency were compared among the clusters. Additionally, to identify a gene signature highly correlated with clinical prognosis, a risk score model was constructed, and the clinical prognosis, immune infiltration, somatic mutation and biological pathways of risk score subgroups were investigated. RESULTS: Our non-apoptosis cell death clusters are robustly predictive of immunotherapy responses. Patients in Cluster B are the most sensitive to immune checkpoint blockades-depended immunotherapy. Our risk score model was also verified as a promising biomarker for clinical prognosis and immunotherapy efficiency. Where, the High-risk score group was more sensitive to immunotherapy. CONCLUSIONS: The novel non-apoptosis cell death-based classification and risk score model could predict the outcome of immunotherapy, and highly associate with immune infiltration. These findings may provide a novel strategy to aid in identificatin of biomarkers and selecting personalized therapeutic strategies.