De novo analysis of bulk RNA-seq data at spatially resolved single-cell resolution

Single-cell and spatial omics: exploring hypothalamic heterogeneity

Elucidating the complex dynamic cellular organization in the hypothalamus is critical for understanding its role in coordinating fundamental body functions. Over the past decade, single-cell and spatial omics technologies have significantly evolved, overcoming initial technical challenges in capturing and analyzing individual cells. These high-throughput omics technologies now offer a remarkable opportunity to comprehend the complex spatiotemporal patterns of transcriptional diversity and cell-type characteristics across the entire hypothalamus. Current single-cell and single-nucleus RNA sequencing methods comprehensively quantify gene expression by exploring distinct phenotypes across various subregions of the hypothalamus. However, single-cell/single-nucleus RNA sequencing requires isolating the cell/nuclei from the tissue, potentially resulting in the loss of spatial information concerning neuronal networks. Spatial transcriptomics methods, by bypassing the cell dissociation, can elucidate the intricate spatial organization of neural networks through their imaging and sequencing technologies. In this review, we highlight the applicative value of single-cell and spatial transcriptomics in exploring the complex molecular-genetic diversity of hypothalamic cell types, driven by recent high-throughput achievements.

The molecular subtypes of autoimmune diseases

Autoimmune diseases (ADs) are characterized by their complexity and a wide range of clinical differences. Despite patients presenting with similar symptoms and disease patterns, their reactions to treatments may vary. The current approach of personalized medicine, which relies on molecular data, is seen as an effective method to address the variability in these diseases. This review examined the pathologic classification of ADs, such as multiple sclerosis and lupus nephritis, over time. Acknowledging the limitations inherent in pathologic classification, the focus shifted to molecular classification to achieve a deeper insight into disease heterogeneity. The study outlined the established methods and findings from the molecular classification of ADs, categorizing systemic lupus erythematosus (SLE) into four subtypes, inflammatory bowel disease (IBD) into two, rheumatoid arthritis (RA) into three, and multiple sclerosis (MS) into a single subtype. It was observed that the high inflammation subtype of IBD, the RA inflammation subtype, and the MS "inflammation & EGF" subtype share similarities. These subtypes all display a consistent pattern of inflammation that is primarily driven by the activation of the JAK-STAT pathway, with the effective drugs being those that target this signaling pathway. Additionally, by identifying markers that are uniquely associated with the various subtypes within the same disease, the study was able to describe the differences between subtypes in detail. The findings are expected to contribute to the development of personalized treatment plans for patients and establish a strong basis for tailored approaches to treating autoimmune diseases.

Comprehensive multi-omics profiling identifies novel molecular subtypes of pancreatic ductal adenocarcinoma

Pancreatic cancer, a highly fatal malignancy, is predicted to rank as the second leading cause of cancer-related death in the next decade. This highlights the urgent need for new insights into personalized diagnosis and treatment. Although molecular subtypes of pancreatic cancer were well established in genomics and transcriptomics, few known molecular classifications are translated to guide clinical strategies and require a paradigm shift. Notably, chronically developing and continuously improving high-throughput technologies and systems serve as an important driving force to further portray the molecular landscape of pancreatic cancer in terms of epigenomics, proteomics, metabonomics, and metagenomics. Therefore, a more comprehensive understanding of molecular classifications at multiple levels using an integrated multi-omics approach holds great promise to exploit more potential therapeutic options. In this review, we recapitulated the molecular spectrum from different omics levels, discussed various subtypes on multi-omics means to move one step forward towards bench-to-beside translation of pancreatic cancer with clinical impact, and proposed some methodological and scientific challenges in store.

Integrated analysis of single-cell RNA-seq and bulk RNA-seq reveals immune suppression subtypes and establishes a novel signature for determining the prognosis in lung adenocarcinoma

BACKGROUND

Lung adenocarcinoma (LUAD) is the most common histological type of lung cancer with lower survival rates. Recent advancements in targeted therapies and immunotherapies targeting immune checkpoints have achieved remarkable success, there is still a large percentage of LUAD that lacks available therapeutic options. Due to tumor heterogeneity, the diagnosis and treatment of LUAD are challenging. Exploring the biology of LUAD and identifying new biomarker and therapeutic targets options are essential.

METHOD

We performed single-cell RNA sequencing (scRNA-seq) of 6 paired primary and adjacent LUAD tissues, and integrative omics analysis of the scRNA-seq, bulk RNA-seq and whole-exome sequencing data revealed molecular subtype characteristics. Our experimental results confirm that CDC25C gene can serve as a potential marker for poor prognosis in LUAD.

RESULTS

We investigated aberrant gene expression in diverse cell types in LUAD via the scRNA-seq data. Moreover, multi-omics clustering revealed four subgroups defined by transcriptional profile and molecular subtype 4 (MS4) with poor survival probability, and immune cell infiltration signatures revealed that MS4 tended to be the immunosuppressive subtype. Our study revealed that the CDC25C gene can be a distinct prognostic biomarker that indicates immune infiltration levels and response to immunotherapy in LUAD patients. Our experimental results concluded that CDC25C expression affects lung cancer cell invasion and migration, might play a key role in regulating Epithelial-Mesenchymal Transition (EMT) pathways.

CONCLUSIONS

Our multi-omics result revealed a comprehensive set of molecular attributes associated with prognosis-related genes in LUAD at the cellular and tissue level. Identification of a subtype of immunosuppressive TME and prognostic signature for LUAD. We identified the cell cycle regulation gene CDC25C affects lung cancer cell invasion and migration, which can be used as a potential biomarker for LUAD.

Next-generation spatial transcriptomics: unleashing the power to gear up translational oncology

The growing advances in spatial transcriptomics (ST) stand as the new frontier bringing unprecedented influences in the realm of translational oncology. This has triggered systemic experimental design, analytical scope, and depth alongside with thorough bioinformatics approaches being constantly developed in the last few years. However, harnessing the power of spatial biology and streamlining an array of ST tools to achieve designated research goals are fundamental and require real-world experiences. We present a systemic review by updating the technical scope of ST across different principal basis in a timeline manner hinting on the generally adopted ST techniques used within the community. We also review the current progress of bioinformatic tools and propose in a pipelined workflow with a toolbox available for ST data exploration. With particular interests in tumor microenvironment where ST is being broadly utilized, we summarize the up-to-date progress made via ST-based technologies by narrating studies categorized into either mechanistic elucidation or biomarker profiling (translational oncology) across multiple cancer types and their ways of deploying the research through ST. This updated review offers as a guidance with forward-looking viewpoints endorsed by many high-resolution ST tools being utilized to disentangle biological questions that may lead to clinical significance in the future.

Single-cell RNA-sequencing reveals a profound immune cell response in human cytomegalovirus-infected humanized mice

Human cytomegalovirus (HCMV) is a common herpesvirus that persistently infects a large portion of the world's population. Despite the robust host immune response, HCMV is able to replicate, evade host defenses, and establish latency throughout the lifespan by developing multiple immunomodulatory strategies, making the studies on the interaction between HCMV infection and host response particularly important. HCMV has a strict host specificity that specifically infects humans. Therefore, most of the in vivo researches of HCMV rely on clinical samples. Fortunately, the establishment of humanized mouse models allows for convenient in-lab animal experiments involving HCMV infection. Single-cell RNA sequencing enables the study of the relationship between viral and host gene expressions at the single-cell level within host cells. In this study, we assessed the gene expression alterations of PBMCs at the single-cell level within HCMV-infected humanized mice, which sheds light onto the virus-host interactions in the context of HCMV infection of humanized mice and provides a valuable dataset for the related researches.

OmicVerse: a framework for bridging and deepening insights across bulk and single-cell sequencing

Single-cell sequencing is frequently affected by "omission" due to limitations in sequencing throughput, yet bulk RNA-seq may contain these ostensibly "omitted" cells. Here, we introduce the single cell trajectory blending from Bulk RNA-seq (BulkTrajBlend) algorithm, a component of the OmicVerse suite that leverages a Beta-Variational AutoEncoder for data deconvolution and graph neural networks for the discovery of overlapping communities. This approach effectively interpolates and restores the continuity of "omitted" cells within single-cell RNA sequencing datasets. Furthermore, OmicVerse provides an extensive toolkit for both bulk and single cell RNA-seq analysis, offering seamless access to diverse methodologies, streamlining computational processes, fostering exquisite data visualization, and facilitating the extraction of significant biological insights to advance scientific research.

scRank infers drug-responsive cell types from untreated scRNA-seq data using a target-perturbed gene regulatory network

Cells respond divergently to drugs due to the heterogeneity among cell populations. Thus, it is crucial to identify drug-responsive cell populations in order to accurately elucidate the mechanism of drug action, which is still a great challenge. Here, we address this problem with scRank, which employs a target-perturbed gene regulatory network to rank drug-responsive cell populations via in silico drug perturbations using untreated single-cell transcriptomic data. We benchmark scRank on simulated and real datasets, which shows the superior performance of scRank over existing methods. When applied to medulloblastoma and major depressive disorder datasets, scRank identifies drug-responsive cell types that are consistent with the literature. Moreover, scRank accurately uncovers the macrophage subpopulation responsive to tanshinone IIA and its potential targets in myocardial infarction, with experimental validation. In conclusion, scRank enables the inference of drug-responsive cell types using untreated single-cell data, thus providing insights into the cellular-level impacts of therapeutic interventions.

Empowering artificial intelligence in characterizing the human primary pacemaker of the heart at single cell resolution

The sinus node (SN) serves as the primary pacemaker of the heart and is the first component of the cardiac conduction system. Due to its anatomical properties and sample scarcity, the cellular composition of the human SN has been historically challenging to study. Here, we employed a novel deep learning deconvolution method, namely Bulk2space, to characterise the cellular heterogeneity of the human SN using existing single-cell datasets of non-human species. As a proof of principle, we used Bulk2Space to profile the cells of the bulk human right atrium using publicly available mouse scRNA-Seq data as a reference. 18 human cell populations were identified, with cardiac myocytes being the most abundant. Each identified cell population correlated to its published experimental counterpart. Subsequently, we applied the deconvolution to the bulk transcriptome of the human SN and identified 11 cell populations, including a population of pacemaker cardiomyocytes expressing pacemaking ion channels (HCN1, HCN4, CACNA1D) and transcription factors (SHOX2 and TBX3). The connective tissue of the SN was characterised by adipocyte and fibroblast populations, as well as key immune cells. Our work unravelled the unique single cell composition of the human SN by leveraging the power of a novel machine learning method.

Simulating multiple variability in spatially resolved transcriptomics with scCube

A pressing challenge in spatially resolved transcriptomics (SRT) is to benchmark the computational methods. A widely-used approach involves utilizing simulated data. However, biases exist in terms of the currently available simulated SRT data, which seriously affects the accuracy of method evaluation and validation. Herein, we present scCube ( https://github.com/ZJUFanLab/scCube ), a Python package for independent, reproducible, and technology-diverse simulation of SRT data. scCube not only enables the preservation of spatial expression patterns of genes in reference-based simulations, but also generates simulated data with different spatial variability (covering the spatial pattern type, the resolution, the spot arrangement, the targeted gene type, and the tissue slice dimension, etc.) in reference-free simulations. We comprehensively benchmark scCube with existing single-cell or SRT simulators, and demonstrate the utility of scCube in benchmarking spot deconvolution, gene imputation, and resolution enhancement methods in detail through three applications.

spatialHeatmap: visualizing spatial bulk and single-cell assays in anatomical images

Visualizing spatial assay data in anatomical images is vital for understanding biological processes in cell, tissue, and organ organizations. Technologies requiring this functionality include traditional one-at-a-time assays, and bulk and single-cell omics experiments, including RNA-seq and proteomics. The spatialHeatmap software provides a series of powerful new methods for these needs, and allows users to work with adequately formatted anatomical images from public collections or custom images. It colors the spatial features (e.g. tissues) annotated in the images according to the measured or predicted abundance levels of biomolecules (e.g. mRNAs) using a color key. This core functionality of the package is called a spatial heatmap plot. Single-cell data can be co-visualized in composite plots that combine spatial heatmaps with embedding plots of high-dimensional data. The resulting spatial context information is essential for gaining insights into the tissue-level organization of single-cell data, or vice versa. Additional core functionalities include the automated identification of biomolecules with spatially selective abundance patterns and clusters of biomolecules sharing similar abundance profiles. To appeal to both non-expert and computational users, spatialHeatmap provides a graphical and a command-line interface, respectively. It is distributed as a free, open-source Bioconductor package (https://bioconductor.org/packages/spatialHeatmap) that users can install on personal computers, shared servers, or cloud systems.

Streamlining spatial omics data analysis with Pysodb

Advances in spatial omics technologies have improved the understanding of cellular organization in tissues, leading to the generation of complex and heterogeneous data and prompting the development of specialized tools for managing, loading and visualizing spatial omics data. The Spatial Omics Database (SODB) was established to offer a unified format for data storage and interactive visualization modules. Here we detail the use of Pysodb, a Python-based tool designed to enable the efficient exploration and loading of spatial datasets from SODB within a Python environment. We present seven case studies using Pysodb, detailing the interaction with various computational methods, ensuring reproducibility of experimental data and facilitating the integration of new data and alternative applications in SODB. The approach offers a reference for method developers by outlining label and metadata availability in representative spatial data that can be loaded by Pysodb. The tool is supplemented by a website ( https://protocols-pysodb.readthedocs.io/ ) with detailed information for benchmarking analysis, and allows method developers to focus on computational models by facilitating data processing. This protocol is designed for researchers with limited experience in computational biology. Depending on the dataset complexity, the protocol typically requires ~12 h to complete.

Metabolic heterogeneity in clear cell renal cell carcinoma revealed by single-cell RNA sequencing and spatial transcriptomics

BACKGROUND

Clear cell renal cell carcinoma is a prototypical tumor characterized by metabolic reprogramming, which extends beyond tumor cells to encompass diverse cell types within the tumor microenvironment. Nonetheless, current research on metabolic reprogramming in renal cell carcinoma mostly focuses on either tumor cells alone or conducts analyses of all cells within the tumor microenvironment as a mixture, thereby failing to precisely identify metabolic changes in different cell types within the tumor microenvironment.

METHODS

Gathering 9 major single-cell RNA sequencing databases of clear cell renal cell carcinoma, encompassing 195 samples. Spatial transcriptomics data were selected to conduct metabolic activity analysis with spatial localization. Developing scMet program to convert RNA-seq data into scRNA-seq data for downstream analysis.

RESULTS

Diverse cellular entities within the tumor microenvironment exhibit distinct infiltration preferences across varying histological grades and tissue origins. Higher-grade tumors manifest pronounced immunosuppressive traits. The identification of tumor cells in the RNA splicing state reveals an association between the enrichment of this particular cellular population and an unfavorable prognostic outcome. The energy metabolism of CD8+ T cells is pivotal not only for their cytotoxic effector functions but also as a marker of impending cellular exhaustion. Sphingolipid metabolism evinces a correlation with diverse macrophage-specific traits, particularly M2 polarization. The tumor epicenter is characterized by heightened metabolic activity, prominently marked by elevated tricarboxylic acid cycle and glycolysis while the pericapsular milieu showcases a conspicuous enrichment of attributes associated with vasculogenesis, inflammatory responses, and epithelial-mesenchymal transition. The scMet facilitates the transformation of RNA sequencing datasets sourced from TCGA into scRNA sequencing data, maintaining a substantial degree of correlation.

CONCLUSIONS

The tumor microenvironment of clear cell renal cell carcinoma demonstrates significant metabolic heterogeneity across various cell types and spatial dimensions. scMet exhibits a notable capability to transform RNA sequencing data into scRNA sequencing data with a high degree of correlation.

Targeting Histamine and Histamine Receptors for Memory Regulation: An Emotional Perspective

Histamine has long been accepted as a pro-cognitive agent. However, lines of evidence have suggested that the roles of histamine in learning and memory processes are much more complex than previously thought. When explained by the spatial perspectives, there are many contradictory results. However, using emotional memory perspectives, we suspect that the histaminergic system may interplay with stress, reward inhibition, and attention to modulate emotional memory formation. The functional diversity of histamine makes it a viable target for clinical management of neuropsychiatric disorders. Here, we update the current knowledge about the functions of histamine in emotional memory and summarize the underlying molecular and neural circuit mechanisms. Finally, we review the main clinical studies about the impacts of histamine-related compounds on memory and discuss insights into future research on the roles of histamine in emotional memory. Despite the recent progress in histamine research, the histaminergic emotional memory circuits are poorly understood, and it is also worth verifying the functions of histamine receptors in a more spatiotemporally specific manner.

Single-cell analysis reveals the immune heterogeneity and interactions in lungs undergoing hepatic ischemia-reperfusion

Hepatic ischemia-reperfusion IR (HIR) is an unavoidable pathophysiological process during liver transplantation, resulting in systematic sterile inflammation and remote organ injury. Acute lung injury (ALI) is a serious complication after liver transplantation with high postoperative morbidity and mortality. However, the underlying mechanism is still unclear. To assess the phenotype and plasticity of various cell types in the lung tissue microenvironment after HIR at the single-cell level, single-cell RNA sequencing (scRNA-seq) was performed using the lungs from HIR-induced mice. In our results, we identified 23 cell types in the lungs after HIR and found that this highly complex ecosystem was formed by subpopulations of bone marrow-derived cells that signaled each other and mediated inflammatory responses in different states and different intervals. We described the unique transcriptional profiles of lung cell clusters and discovered two novel cell subtypes (Tspo+Endothelial cells and Vcan+ monocytes), as well as the endothelial cell-immune cell and immune cell-T cell clusters interactome. In addition, we found that S100 calcium binding protein (S100a8/a9), specifically and highly expressed in immune cell clusters of lung tissues and exhibited detrimental effects. Finally, the cellular landscape of the lung tissues after HIR was established, highlighting the heterogeneity and cellular interactions between major immune cells in HIR-induced lungs. Our findings provided new insights into the mechanisms of HIR-induced ALI and offered potential therapeutic target to prevent ALI after liver transplantation.

Multiscale modeling of drug resistance in glioblastoma with gene mutations and angiogenesis

Drug resistance is a prominent impediment to the efficacy of targeted therapies across various cancer types, including glioblastoma (GBM). However, comprehending the intricate intracellular and extracellular mechanisms underlying drug resistance remains elusive. Empirical investigations have elucidated that genetic aberrations, such as gene mutations, along with microenvironmental adaptation, notably angiogenesis, act as pivotal drivers of tumor progression and drug resistance. Nonetheless, mathematical models frequently compartmentalize these factors in isolation. In this study, we present a multiscale agent-based model of GBM, encompassing cellular dynamics, intricate signaling pathways, gene mutations, angiogenesis, and therapeutic interventions. This integrative framework facilitates an exploration of the interplay between genetic mutations and the vascular microenvironment in shaping the dynamic evolution of tumors during treatment with tyrosine kinase inhibitor. Our simulations unveil that mutations influencing the migration and proliferation of tumor cells expedite the emergence of phenotype heterogeneity, thereby exacerbating tumor invasion under both treated and untreated conditions. Moreover, angiogenesis proximate to the tumor fosters a protumoral milieu, augmenting mutation-induced drug resistance by increasing the survival rate of tumor cells. Collectively, our findings underscore the dual roles of intrinsic genetic mutations and extrinsic microenvironmental adaptations in steering tumor growth and drug resistance. Finally, we substantiate our model predictions concerning the impact of gene mutations and angiogenesis on the responsiveness of targeted therapies by integrating single-cell RNA-seq, spatial transcriptomics, bulk RNA-seq, and clinical data from GBM patients. The multidimensional approach enhances our understanding of the complexities governing drug resistance in glioma and offers insights into potential therapeutic strategies.

Laser capture microdissection for biomedical research: towards high-throughput, multi-omics, and single-cell resolution

Spatial omics technologies have become powerful methods to provide valuable insights into cells and tissues within a complex context, significantly enhancing our understanding of the intricate and multifaceted biological system. With an increasing focus on spatial heterogeneity, there is a growing need for unbiased, spatially resolved omics technologies. Laser capture microdissection (LCM) is a cutting-edge method for acquiring spatial information that can quickly collect regions of interest (ROIs) from heterogeneous tissues, with resolutions ranging from single cells to cell populations. Thus, LCM has been widely used for studying the cellular and molecular mechanisms of diseases. This review focuses on the differences among four types of commonly used LCM technologies and their applications in omics and disease research. Key attributes of application cases are also highlighted, such as throughput and spatial resolution. In addition, we comprehensively discuss the existing challenges and the great potential of LCM in biomedical research, disease diagnosis, and targeted therapy from the perspective of high-throughput, multi-omics, and single-cell resolution.

Whole genome deconvolution unveils Alzheimer's resilient epigenetic signature

Assay for Transposase Accessible Chromatin by sequencing (ATAC-seq) accurately depicts the chromatin regulatory state and altered mechanisms guiding gene expression in disease. However, bulk sequencing entangles information from different cell types and obscures cellular heterogeneity. To address this, we developed Cellformer, a deep learning method that deconvolutes bulk ATAC-seq into cell type-specific expression across the whole genome. Cellformer enables cost-effective cell type-specific open chromatin profiling in large cohorts. Applied to 191 bulk samples from 3 brain regions, Cellformer identifies cell type-specific gene regulatory mechanisms involved in resilience to Alzheimer's disease, an uncommon group of cognitively healthy individuals that harbor a high pathological load of Alzheimer's disease. Cell type-resolved chromatin profiling unveils cell type-specific pathways and nominates potential epigenetic mediators underlying resilience that may illuminate therapeutic opportunities to limit the cognitive impact of the disease. Cellformer is freely available to facilitate future investigations using high-throughput bulk ATAC-seq data.

Promise of spatially resolved omics for tumor research

Tumors are spatially heterogeneous tissues that comprise numerous cell types with intricate structures. By interacting with the microenvironment, tumor cells undergo dynamic changes in gene expression and metabolism, resulting in spatiotemporal variations in their capacity for proliferation and metastasis. In recent years, the rapid development of histological techniques has enabled efficient and high-throughput biomolecule analysis. By preserving location information while obtaining a large number of gene and molecular data, spatially resolved metabolomics (SRM) and spatially resolved transcriptomics (SRT) approaches can offer new ideas and reliable tools for the in-depth study of tumors. This review provides a comprehensive introduction and summary of the fundamental principles and research methods used for SRM and SRT techniques, as well as a review of their applications in cancer-related fields.

Integrative multi-omics and systems bioinformatics in translational neuroscience: A data mining perspective

Bioinformatic analysis of large and complex omics datasets has become increasingly useful in modern day biology by providing a great depth of information, with its application to neuroscience termed neuroinformatics. Data mining of omics datasets has enabled the generation of new hypotheses based on differentially regulated biological molecules associated with disease mechanisms, which can be tested experimentally for improved diagnostic and therapeutic targeting of neurodegenerative diseases. Importantly, integrating multi-omics data using a systems bioinformatics approach will advance the understanding of the layered and interactive network of biological regulation that exchanges systemic knowledge to facilitate the development of a comprehensive human brain profile. In this review, we first summarize data mining studies utilizing datasets from the individual type of omics analysis, including epigenetics/epigenomics, transcriptomics, proteomics, metabolomics, lipidomics, and spatial omics, pertaining to Alzheimer's disease, Parkinson's disease, and multiple sclerosis. We then discuss multi-omics integration approaches, including independent biological integration and unsupervised integration methods, for more intuitive and informative interpretation of the biological data obtained across different omics layers. We further assess studies that integrate multi-omics in data mining which provide convoluted biological insights and offer proof-of-concept proposition towards systems bioinformatics in the reconstruction of brain networks. Finally, we recommend a combination of high dimensional bioinformatics analysis with experimental validation to achieve translational neuroscience applications including biomarker discovery, therapeutic development, and elucidation of disease mechanisms. We conclude by providing future perspectives and opportunities in applying integrative multi-omics and systems bioinformatics to achieve precision phenotyping of neurodegenerative diseases and towards personalized medicine.

Dissecting the brain with spatially resolved multi-omics

Recent studies have highlighted spatially resolved multi-omics technologies, including spatial genomics, transcriptomics, proteomics, and metabolomics, as powerful tools to decipher the spatial heterogeneity of the brain. Here, we focus on two major approaches in spatial transcriptomics (next-generation sequencing-based technologies and image-based technologies), and mass spectrometry imaging technologies used in spatial proteomics and spatial metabolomics. Furthermore, we discuss their applications in neuroscience, including building the brain atlas, uncovering gene expression patterns of neurons for special behaviors, deciphering the molecular basis of neuronal communication, and providing a more comprehensive explanation of the molecular mechanisms underlying central nervous system disorders. However, further efforts are still needed toward the integrative application of multi-omics technologies, including the real-time spatial multi-omics analysis in living cells, the detailed gene profile in a whole-brain view, and the combination of functional verification.

Reconstruction of the cell pseudo-space from single-cell RNA sequencing data with scSpace

Tissues are highly complicated with spatial heterogeneity in gene expression. However, the cutting-edge single-cell RNA-seq technology eliminates the spatial information of individual cells, which contributes to the characterization of cell identities. Herein, we propose single-cell spatial position associated co-embeddings (scSpace), an integrative method to identify spatially variable cell subpopulations by reconstructing cells onto a pseudo-space with spatial transcriptome references (Visium, STARmap, Slide-seq, etc.). We benchmark scSpace with both simulated and biological datasets, and demonstrate that scSpace can accurately and robustly identify spatially variated cell subpopulations. When employed to reconstruct the spatial architectures of complex tissue such as the brain cortex, the small intestinal villus, the liver lobule, the kidney, the embryonic heart, and others, scSpace shows promising performance on revealing the pairwise cellular spatial association within single-cell data. The application of scSpace in melanoma and COVID-19 exhibits a broad prospect in the discovery of spatial therapeutic markers.

Aged G Protein-Coupled Receptor Kinase 3 (Grk3)-Deficient Mice Exhibit Enhanced Osteoclastogenesis and Develop Bone Lesions Analogous to Human Paget's Disease of Bone

Paget's Disease of Bone (PDB) is a metabolic bone disease that is characterized by dysregulated osteoclast function leading to focal abnormalities of bone remodeling. It can lead to pain, fracture, and bone deformity. G protein-coupled receptor kinase 3 (GRK3) is an important negative regulator of G protein-coupled receptor (GPCR) signaling. GRK3 is known to regulate GPCR function in osteoblasts and preosteoblasts, but its regulatory function in osteoclasts is not well defined. Here, we report that Grk3 expression increases during osteoclast differentiation in both human and mouse primary cells and established cell lines. We also show that aged mice deficient in Grk3 develop bone lesions similar to those seen in human PDB and other Paget's Disease mouse models. We show that a deficiency in Grk3 expression enhances osteoclastogenesis in vitro and proliferation of hematopoietic osteoclast precursors in vivo but does not affect the osteoclast-mediated bone resorption function or cellular senescence pathway. Notably, we also observe decreased Grk3 expression in peripheral blood mononuclear cells of patients with PDB compared with age- and gender-matched healthy controls. Our data suggest that GRK3 has relevance to the regulation of osteoclast differentiation and that it may have relevance to the pathogenesis of PDB and other metabolic bone diseases associated with osteoclast activation.

A guidebook of spatial transcriptomic technologies, data resources and analysis approaches

Advances in transcriptomic technologies have deepened our understanding of the cellular gene expression programs of multicellular organisms and provided a theoretical basis for disease diagnosis and therapy. However, both bulk and single-cell RNA sequencing approaches lose the spatial context of cells within the tissue microenvironment, and the development of spatial transcriptomics has made overall bias-free access to both transcriptional information and spatial information possible. Here, we elaborate development of spatial transcriptomic technologies to help researchers select the best-suited technology for their goals and integrate the vast amounts of data to facilitate data accessibility and availability. Then, we marshal various computational approaches to analyze spatial transcriptomic data for various purposes and describe the spatial multimodal omics and its potential for application in tumor tissue. Finally, we provide a detailed discussion and outlook of the spatial transcriptomic technologies, data resources and analysis approaches to guide current and future research on spatial transcriptomics.