Cardiovascular diseases (CVDs) impose a significant burden worldwide. Despite the elucidation of the etiology and underlying molecular mechanisms of CVDs by numerous studies and recent discovery of effective drugs, their morbidity, disability, and mortality are still high. Therefore, precise risk stratification and effective targeted therapies for CVDs are warranted. Recent improvements in single-cell RNA sequencing and spatial transcriptomics have improved our understanding of the mechanisms and cells involved in cardiovascular phylogeny and CVDs. Single-cell RNA sequencing can facilitate the study of the human heart at remarkably high resolution and cellular and molecular heterogeneity. However, this technique does not provide spatial information, which is essential for understanding homeostasis and disease. Spatial transcriptomics can elucidate intracellular interactions, transcription factor distribution, cell spatial localization, and molecular profiles of mRNA and identify cell populations causing the disease and their underlying mechanisms, including cell crosstalk. Herein, we introduce the main methods of RNA-seq and spatial transcriptomics analysis and highlight the latest advances in cardiovascular research. We conclude that single-cell RNA sequencing interprets disease progression in multiple dimensions, levels, perspectives, and dynamics by combining spatial and temporal characterization of the clinical phenome with multidisciplinary techniques such as spatial transcriptomics. This aligns with the dynamic evolution of CVDs (e.g., \"angina-myocardial infarction-heart failure\" in coronary artery disease). The study of pathways for disease onset and mechanisms (e.g., age, sex, comorbidities) in different patient subgroups should improve disease diagnosis and risk stratification. This can facilitate precise individualized treatment of CVDs.

BACKGROUND: At present, immunotherapy has become a powerful treatment for advanced gastric cancer (AGC), but not all patients can benefit from it. According to the latest research, the impact of B cell subpopulations on the immune microenvironment of gastric cancer (GC) is unknown. Exploring whether the interaction between B cells and tumor cells in GC affects the effectiveness of immunotherapy has attracted our interest.
METHODS: This study involved the re-analysis of single-cell RNA (scRNA) and spatial transcriptomics (ST) data from publicly available datasets. The focus was on investigating the subpopulations and differentiation trajectories of B cells in the gastric cancer (GC) tumor immune microenvironment (TIME). Spatial transcriptomics (ST) and multiple immunofluorescence (mIF) revealed a clear co-localization pattern between B cells and tumor cells. Multiple immunotherapy datasets were collected to identify unique immunotherapy biomarkers. The unique immunotherapeutic potential of targeting CCL28 was validated through a mouse gastric cancer model. In addition, flow cytometry revealed changes in the tumor immune microenvironment targeting CCL28.
RESULTS: The re-analysis of ST data from multiple cancer types revealed a co-localization pattern between B cells and tumor cells. A significant number of IgA plasma cells were identified in the GC TIME. Five different tumor-infiltrating B cell subpopulations and two unique B cell differentiation trajectories were characterized, along with seven GC-related states. By analyzing the communication between GC cells and B cells, it was further discovered that tumor cells can influence and recruit plasma cells through CCL28-CCR10 signaling. Additionally, there was a crosstalk between GC cells and B cells. Finally, we identified the LAMA/CD44 signaling axis as a potential prognostic marker for immunotherapy through a large amount of immunotherapy data. We also validated through various animal tumor models that targeting CCL28 can significantly promote CD8+T cell infiltration and function in the TME by regulating B cell and plasma cell functions, and has the ability to synergize immunotherapy.
CONCLUSIONS: The co-localization and crosstalk between GC cells and B cells significantly affect the efficacy of immunotherapy, and inhibiting the CCL28-CCR10 signal axis is a potential immunotherapy target for GC. Meanwhile, LAMA/CD44 pair may be a potential adverse indicator for immunotherapy and tumor prognosis.

Myasthenia gravis (MG) is etiologically associated with thymus abnormalities, but its pathology in the thymus remains unclear. In this study, we attempt to narrow down the features associated with MG using spatial transcriptome analysis of thymoma and thymic hyperplasia samples. We find that the majority of thymomas are constituted by the cortical region. However, the small medullary region is enlarged in seropositive thymomas and contains polygenic enrichment and MG-specific germinal center structures. Neuromuscular medullary thymic epithelial cells, previously identified as MG-specific autoantigen-producing cells, are enriched in the cortico-medullary junction. The medulla is characterized by a specific chemokine pattern and immune cell composition, including migratory dendritic cells and effector regulatory T cells. Similar germinal center structures and immune microenvironments are also observed in the thymic hyperplasia medulla. This study shows that the medulla and junction areas are linked to MG pathology and provides insights into future MG research.

UNASSIGNED: Microvascular invasion (MVI) stands as a pivotal pathological hallmark of hepatocellular carcinoma (HCC), closely linked to unfavorable prognosis, early recurrence, and metastatic progression. However, the precise mechanistic underpinnings governing its onset and advancement remain elusive.
UNASSIGNED: In this research, we downloaded bulk RNA-seq data from the TCGA and HCCDB repositories, single-cell RNA-seq data from the GEO database, and spatial transcriptomics data from the CNCB database. Leveraging the Scissor algorithm, we delineated prognosis-related cell subpopulations and discerned a distinct MVI-related malignant cell subtype. A comprehensive exploration of these malignant cell subpopulations was undertaken through pseudotime analysis and cell-cell communication scrutiny. Furthermore, we engineered a prognostic model grounded in MVI-related genes, employing 101 algorithm combinations integrated by 10 machine-learning algorithms on the TCGA training set. Rigorous evaluation ensued on internal testing sets and external validation sets, employing C-index, calibration curves, and decision curve analysis (DCA).
UNASSIGNED: Pseudotime analysis indicated that malignant cells, showing a positive correlation with MVI, were primarily concentrated in the early to middle stages of differentiation, correlating with an unfavorable prognosis. Importantly, these cells showed significant enrichment in the MYC pathway and were involved in extensive interactions with diverse cell types via the MIF signaling pathway. The association of malignant cells with the MVI phenotype was corroborated through validation in spatial transcriptomics data. The prognostic model we devised demonstrated exceptional sensitivity and specificity, surpassing the performance of most previously published models. Calibration curves and DCA underscored the clinical utility of this model.
UNASSIGNED: Through integrated multi-transcriptomics analysis, we delineated MVI-related malignant cells and elucidated their biological functions. This study provided novel insights for managing HCC, with the constructed prognostic model offering valuable support for clinical decision-making.

The cortical microenvironment surrounding malignant glioblastoma is a source of depolarizing crosstalk favoring hyperexcitability, tumor expansion, and immune evasion. Neosynaptogenesis, excess glutamate, and altered intrinsic membrane currents contribute to excitability dyshomeostasis, yet only half of the cases develop seizures, suggesting that tumor and host genomics, along with location, rather than mass effect, play a critical role. We analyzed the spatial contours and expression of 358 clinically validated human epilepsy genes in the human glioblastoma transcriptome compared to non-tumor adult and developing cortex datasets. Nearly half, including dosage-sensitive genes whose expression levels are securely linked to monogenic epilepsy, are strikingly enriched and aberrantly regulated at the leading edge, supporting a complex epistatic basis for peritumoral epileptogenesis. Surround hyperexcitability induced by complex patterns of proepileptic gene expression may explain the limited efficacy of narrowly targeted antiseizure medicines and the persistence of epilepsy following tumor resection and clarify why not all brain tumors provoke seizures.

The human adrenal gland is a complex endocrine tissue. Studies on adrenal renewal have been limited to animal models or human foetuses. Enhancing our understanding of adult human adrenal homeostasis is crucial for gaining insights into the pathogenesis of adrenal diseases, such as adrenocortical tumours. Here, we present a comprehensive cellular genomics analysis of the adult human normal adrenal gland, combining single-nuclei RNA sequencing and spatial transcriptome data to reconstruct adrenal gland homeostasis. As expected, we identified primary cells of the various zones of the adrenal cortex and medulla, but we also uncovered additional cell types. They constitute the adrenal microenvironment, including immune cells, mostly composed of a large population of M2 macrophages, and new cell populations, including different subpopulations of vascular-endothelial cells and cortical-neuroendocrine cells. Utilizing spatial transcriptome and pseudotime trajectory analysis, we support evidence of the centripetal dynamics of adrenocortical cell maintenance and the essential role played by Wnt/β-catenin, sonic hedgehog, and fibroblast growth factor pathways in the adult adrenocortical homeostasis. Furthermore, we compared single-nuclei transcriptional profiles obtained from six healthy adrenal glands and twelve adrenocortical adenomas. This analysis unveiled a notable heterogeneity in cell populations within the adenoma samples. In addition, we identified six distinct adenoma-specific clusters, each with varying distributions based on steroid profiles and tumour mutational status. Overall, our results provide novel insights into adrenal homeostasis and molecular mechanisms potentially underlying early adrenocortical tumorigenesis and/or autonomous steroid secretion. Our cell atlas represents a powerful resource to investigate other adrenal-related pathologies.

Bamboo with its remarkable growth rate and economic significance, offers an ideal system to investigate the molecular basis of organogenesis in rapidly growing plants, particular in monocots, where gene regulatory networks governing the maintenance and differentiation of shoot apical and intercalary meristems remain a subject of controversy. We employed both spatial and single-nucleus transcriptome sequencing on 10× platform to precisely dissect the gene functions in various tissues and early developmental stages of bamboo shoots. Our comprehensive analysis reveals distinct cell trajectories during shoot development, uncovering critical genes and pathways involved in procambium differentiation, intercalary meristem formation, and vascular tissue development. Spatial and temporal expression patterns of key regulatory genes, particularly those related to hormone signaling and lipid metabolism, strongly support the hypothesis that intercalary meristem origin from surrounded parenchyma cells. Specific gene expressions in intercalary meristem exhibit regular and dispersed distribution pattern, offering clues for understanding the intricate molecular mechanisms that drive the rapid growth of bamboo shoots. The single-nucleus and spatial transcriptome analysis reveal a comprehensive landscape of gene activity, enhancing the understanding of the molecular architecture of organogenesis and providing valuable resources for future genomic and genetic studies relying on identities of specific cell types.

BACKGROUND: Lung cancer is a highly aggressive and prevalent disease worldwide. By the time it is first diagnosed, distant metastases have usually already occurred. Among them, the prognosis of patients with brain metastasis from lung cancer is very poor. Therefore, it is particularly important to identify the evolutionary status of tumor cells during lung cancer brain metastases and discover the underlying mechanisms of lung cancer brain metastases.
METHODS: In this study, we analysed three types of data: single-cell RNA sequencing, bulk RNA sequencing, and spatial transcriptome. Firstly, we identified early metastatic epithelial cell clusters (EMEC) using CNV and trajectory analysis in scRNA-seq data. Secondly, we integrated scRNA-seq and spatial transcriptome data with the help of MIA (Multimodal intersection analysis) to explore the biological characteristics of EMEC. Finally, we used bulk RNA-seq data to validate the molecular characteristics of EMEC.
RESULTS: A total of 55,763 single cells were obtained and divided into 9 cell types. In brain metastasis, we found a significantly higher proportion of epithelial cells. In addition, we identified a specific subpopulation of epithelial cells, which was named as \"early metastatic epithelial cell clusters (EMEC)\". It is enriched in oxidative phosphorylation, coagulation, complement. Moreover, we also found that EMEC underwent cellular communication with other immune cells through ligand-receptor pairs such as MIF-(CD74 + CXCR4) and MIF-(CD74 + CD44). Next, we validated that EMEC were associated with poor clinical prognosis using three independent external datasets. Finally, spatial transcriptome analysis revealed specificity in the spatial distribution of EMEC, which shifted from the peripheral regions to the central regions of the tumour as the depth of tumor invasion progressed.
CONCLUSIONS: This study reveals the potential molecular mechanisms of lung cancer brain metastasis from both single-cell and spatial transcriptomic perspectives, providing biological insights and clinical reference value for detecting patients suffering from lung cancer brain metastasis.

UNASSIGNED: Hepatocellular carcinoma (HCC), representing more than 80% of primary liver cancer cases, lacks satisfactory etiology and diagnostic methods. This study aimed to elucidate the role of programmed cell death-associated genes (CDRGs) in HCC by constructing a diagnostic model using single-cell RNA sequencing (scRNA-seq) and RNA sequencing (RNA-seq) data.
UNASSIGNED: Six categories of CDRGs, including apoptosis, necroptosis, autophagy, pyroptosis, ferroptosis, and cuproptosis, were collected. RNA-seq data from blood-derived exosomes were sourced from the exoRBase database, RNA-seq data from cancer tissues from the TCGA database, and scRNA-seq data from the GEO database. Subsequently, we intersected the differentially expressed genes (DEGs) of the HCC cohort from exoRBase and TCGA databases with CDRGs, as well as DEGs obtained from single-cell datasets. Candidate biomarker genes were then screened using clinical indicators and a machine learning approach, resulting in the construction of a seven-gene diagnostic model for HCC. Additionally, scRNA-seq and spatial transcriptome sequencing (stRNA-seq) data of HCC from the Mendeley data portal were used to investigate the underlying mechanisms of these seven key genes and their association with immune checkpoint blockade (ICB) therapy. Finally, we validated the expression of key molecules in tissues and blood-derived exosomes through quantitative Polymerase Chain Reaction (qPCR) and immunohistochemistry experiments.
UNASSIGNED: Collectively, we obtained a total of 50 samples and 104,288 single cells. Following the meticulous screening, we established a seven-gene diagnostic model for HCC, demonstrating high diagnostic efficacy in both the exoRBase HCC cohort (training set: AUC = 1; testing set: AUC = 0.847) and TCGA HCC cohort (training set: AUC = 1; testing set: AUC = 0.976). Subsequent analysis revealed that HCC cluster 3 exhibited a higher stemness index and could serve as the starting point for the differentiation trajectory of HCC cells, also displaying more abundant interactions with other cell types in the microenvironment. Notably, key genes TRIB3 and NQO1 displayed elevated expression levels in HCC cells. Experimental validation further confirmed their elevated expression in both tumor tissues and blood-derived exosomes of cancer patients. Additionally, stRNA analysis not only substantiated these findings but also suggested that patients with high TRIB3 and NQO1 expression might respond more favorably to ICB therapy.
UNASSIGNED: The seven-gene diagnostic model demonstrated remarkable accuracy in HCC screening, with TRIB3 emerging as a promising diagnostic tool and therapeutic target for HCC.

Cuproptosis is a novel type to regulate cell death with copper-dependent manner, and has been reported to involve in the occurrence and development of various malignant tumors. However, the association between cuproptosis and the tumor microenvironment (TME) of clear cell renal cell carcinoma (ccRCC) remained unclear. To address this question, we integrated the single cell RNA sequencing (scRNA-seq) datasets of ccRCC across different stages, systematically examined the distinctive expression patterns of cuproptosis-related genes (CRGs) within the TME of ccRCC, and explored the crucial signatures using the spatial transcriptome sequencing (ST-seq) dataset. The cuproptosis activities reduced in cancer tissues along with the ccRCC development, and recovered after therapy. We identified HILPDA+ ccRCC1 subtype, characterized with hypoxia, as cuproptosis susceptible cells associated with a better prognosis. The main co-expression modules of HILPDA+ ccRCC1 subtype highlighted the role in anion transport, response to oxygen species and PD-L1-PD-1 pathway. Furthermore, the immunosuppressive cells might interact with HILPDA+ ccRCC1 subtype via HAVCR2-LGALS9, C3-C3AR1, HLA-A-CD8B and HLA-C-CD8A axises to shape the cuproptosis-related TME landscape. In summary, we anticipate that this study will offer valuable insights and potential strategies of cuproptosis for therapy of ccRCC.