Search

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.

COVID-19 has caused severe threats to lives and damage to property worldwide. The immunopathology of the disease is of particular concern. Currently, researchers have used gene co-expression networks (GCNs) to deepen the study of molecular mechanisms of immune responses to COVID-19. However, most efforts have not fully explored dynamic changes of cell-type-specific molecular networks in the disease process. This study proposes a GCN construction pipeline named single-cell Disease Progression cellular module analysis (scDisProcema), which can trace dynamic changes of immune system response during disease progression using single-cell data. Here, scDisProcema considers changes in cell fate and expression patterns during disease development, identifying gene modules responsible for different immune cells. The hub genes are screened for each module by the specific expression level and the intercellular connectivity of modules. Based on functional items enriched by each gene module, we elucidate the biological processes of different cells involved in disease development and explain the molecular mechanisms underlying the process of cell depletion or proliferation caused by disease. Compared with traditional WGCNA methods, scDisProcema can make more convenient use of the heterogeneity information provided by scRNA-seq data and has great potential in exploring molecular changes during disease progression and organ development.

Spatially resolved transcriptomics provides genetic information in space toward elucidation of the spatial architecture in intact organs and the spatially resolved cell-cell communications mediating tissue homeostasis, development, and disease. To facilitate inference of spatially resolved cell-cell communications, we here present SpaTalk, which relies on a graph network and knowledge graph to model and score the ligand-receptor-target signaling network between spatially proximal cells by dissecting cell-type composition through a non-negative linear model and spatial mapping between single-cell transcriptomic and spatially resolved transcriptomic data. The benchmarked performance of SpaTalk on public single-cell spatial transcriptomic datasets is superior to that of existing inference methods. Then we apply SpaTalk to STARmap, Slide-seq, and 10X Visium data, revealing the in-depth communicative mechanisms underlying normal and disease tissues with spatial structure. SpaTalk can uncover spatially resolved cell-cell communications for single-cell and spot-based spatially resolved transcriptomic data universally, providing valuable insights into spatial inter-cellular tissue dynamics.

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 algorithm to distinguish spatially variable cell subclusters by reconstructing cells onto a pseudo-space with spatial transcriptome references (Visium, STARmap, Slide-seq, etc.). We demonstrated that scSpace can define biologically meaningful cell subpopulations neglected by single-cell RNA-seq or spatially resolved transcriptomics. The use of scSpace to uncover the spatial association within single-cell data, reproduced, the hierarchical distribution of cells in the brain cortex and liver lobules, and the regional variation of cells in heart ventricles and the intestinal villus. scSpace identified cell subclusters in intratelencephalic neurons, which were confirmed by their biomarkers. The application of scSpace in melanoma and Covid-19 exhibited a broad prospect in the discovery of spatial therapeutic markers.

Uncovering the tissue molecular architecture at single-cell resolution could help better understand organisms’ biological and pathological processes. However, bulk RNA-seq can only measure gene expression in cell mixtures, without revealing the transcriptional heterogeneity and spatial patterns of single cells. Herein, we introduce Bulk2Space (https://github.com/ZJUFanLab/bulk2space), a deep learning framework-based spatial deconvolution algorithm that can simultaneously disclose the spatial and cellular heterogeneity of bulk RNA-seq data using existing single-cell and spatial transcriptomics references. The use of bulk transcriptomics to validate Bulk2Space unveils, in particular, the spatial variance of immune cells in different tumor regions, the molecular and spatial heterogeneity of tissues during inflammation-induced tumorigenesis, and spatial patterns of novel genes in different cell types. Moreover, Bulk2Space is utilized to perform spatial deconvolution analysis on bulk transcriptome data from two different mouse brain regions derived from our in-house developed sequencing approach termed Spatial-seq. We have not only reconstructed the hierarchical structure of the mouse isocortex but also further annotated cell types that were not identified by original methods in the mouse hypothalamus.

Uncovering the tissue molecular architecture at single-cell resolution could help better understand organisms’ biological and pathological processes. However, bulk RNA-seq can only measure gene expression in cell mixtures, without revealing the transcriptional heterogeneity and spatial patterns of single cells. Herein, we introduce Bulk2Space (https://github.com/ZJUFanLab/bulk2space), a deep learning framework-based spatial deconvolution algorithm that can simultaneously disclose the spatial and cellular heterogeneity of bulk RNA-seq data using existing single-cell and spatial transcriptomics references. The use of bulk transcriptomics to validate Bulk2Space unveils, in particular, the spatial variance of immune cells in different tumor regions, the molecular and spatial heterogeneity of tissues during inflammation-induced tumorigenesis, and spatial patterns of novel genes in different cell types. Moreover, Bulk2Space is utilized to perform spatial deconvolution analysis on bulk transcriptome data from two different mouse brain regions derived from our in-house developed sequencing approach termed Spatial-seq. We have not only reconstructed the hierarchical structure of the mouse isocortex but also further annotated cell types that were not identified by original methods in the mouse hypothalamus.

Background Microbiome big data from population-scale cohorts holds the key to unleash the power of microbiomes to overcome critical challenges in disease control, treatment and precision medicine. However, variations introduced during data generation and processing limit the comparisons among independent studies in respect of interpretability. Although multiple databases have been constructed as platforms for data reuse, they are of limited value since only raw sequencing files are considered. Description Here, we present MetaGeneBank, a standardized database that provides details on sample collection and sequencing, and abundances of genes, microbiota and molecular functions for 4470 raw sequencing files (over 12 TB) collected from 16 studies covering over 10 types of diseases and 14 countries using a unified data-processing pipeline. The incorporation of tools that enable browsing and searching with descriptive attributes, gene sequences, microbiota and functions makes the database user-friendly. We found that the source of specimen contributes more than sequencing centers or platforms to the variations of microbiota. Special attention should be paid when re-analyzing sequencing files from different countries. Conclusions Collectively, MetaGeneBank provides a gateway to utilize the untapped potential of gut metagenomic data in helping fighting against human diseases. With the continuous updating of the database in terms of data volume, data types and sample types, MetaGeneBank would undoubtedly be the benchmarking database in the future in respect of data reuse, and would be valuable in translational science.

The heterogeneity of liver non-parenchymal cells (NPCs) is essential for liver structure and function. However, the current understanding of liver NPCs, especially in different liver diseases, remains incompletely elucidated. Here, a single-cell transcriptome atlas of 171,814 NPCs from healthy and 5 typical liver disease mouse models, including alcoholic liver disease, nonalcoholic steatohepatitis (NASH), drug-induced liver injury, cholestatic, and ischemia-reperfusion liver injury is constructed. The inter- and intra-group heterogeneity of 12 types (and numerous subtypes) of NPCs involving endothelial cells, hepatic stellate cells (HSCs), neutrophils, T cells, and mononuclear phagocytes (MPs) are summarized. A protective subtype of neutrophils characterized by Chil3high is validated and found significantly increasing only in drug-induced and cholestatic liver injury models. Transcriptional regulatory network analysis reveals disease-specific transcriptional reprogramming. Metabolic activity analysis indicates that fibrosis is accompanied by increases in glycolysis and retinol metabolism in activated HSCs and MPs. Moreover, we found that cell-cell interactions between cholangiocytes and immune cells contribute more to cholestatic liver fibrosis compared with NASH, while HSCs are more important for NASH fibrosis. Our atlas, together with an interactive website provides a systematic view of highly heterogeneous NPCs and a valuable resource to better understand pathological mechanisms underlying liver diseases.