Search

SummaryIn the era where transcriptome profiling moves towards single-cell and spatial resolutions, the traditional co-expression analysis lacks the power to fully utilize such rich information to unravel spatial gene associations. Here we present a Python package called Spatial Enrichment Analysis of Gene Associations using L-index (SEAGAL) to detect and visualize spatial gene correlations at both single-gene and gene-set levels. Our package takes spatial transcriptomics data sets with gene expression and the aligned spatial coordinates as input. It allows for analyzing and visualizing spatial correlations at both single-gene and gene-set levels. The output could be visualized as volcano plots and heatmaps with a few lines of code, thus providing an easy-yet-comprehensive tool for mining spatial gene associations.Availability and ImplementationThe Python package SEAGAL can be installed using pip:https://pypi.org/project/seagal/. The source code and step-by-step tutorials are available at:https://github.com/linhuawang/SEAGAL.Contactlinhuaw@bcm.edu

MotivationWith the development in single-cell multi-omics sequencing technology and data integration algorithms, we have entered the single-cell multi-omics era. Current multi-omics analysis algorithms failed to systematically dissect the heterogeneity within the datasets when inferring cis-regulatory events. Thus, there is a need for cis-regulatory element inferring algorithms that considers the cellular heterogeneity.ResultsHere, we propose scGREAT, a single-cell multi-omics regulatory state analysis Python package with a rapid graph-based correlation measurementL. The graph-based correlation method assigns each cell a localLindex, pinpointing specific cell groups of certain regulatory states. Such single-cell resolved regulatory state information enables the heterogeneity analysis equipped in the package. Applying scGREAT to the 10X Multiome PBMC dataset, we demonstrated how it could help subcluster cell types, infer regulation-based pseudo-time trajectory, discover feature modules, and find cluster-specific regulatory gene-peak pairs. Besides, we showed that global L index, which is the average of all local L values, is a better replacement for Pearson’s r in ruling out confounding regulatory relationships that are not of research interests.Availabilityhttps://github.com/ChaozhongLiu/scGREAT

Sequencing-based spatial transcriptomics (ST) techniques have been groundbreaking in dissecting cell-cell communications within tissues by profiling positional gene expression. However, the most widely used ST technique, Visium Spatial Gene Expression by 10x Genomics (Visium), does not provide single-cell resolution, making it difficult to profile cell type-level information. Many reference-based deconvolution methods have been developed to increase its resolution, but the platform and batch effects between the reference and ST data compromise their accuracy. Here, we propose a new approach,Region-based cellSorting (ReSort), that generates a pseudo-internal-reference to reduce these platform effects. By simulating ST datasets under various scenarios, we demonstrate that ReSort significantly improves the accuracy of six state-of-the-art reference-based deconvolution methods. Moreover, applying ReSort to a mouse breast cancer tumor bearing both epithelial and mesenchymal clones identifies the spatial differences of immune cells between the clones, providing important insights for understanding the relationship between epithelial-mesenchymal transition and immune infiltration in breast cancer.

With the growing availability of Alzheimer’s disease (AD) transcriptomic data, several studies have nominated new therapeutic targets. However, a major challenge is accounting for latent (hidden) factors which affect the discovery of therapeutic targets. Using unsupervised machine learning, we identified a latent factor in brain tissue, and we validated the factor in AD and normal samples, across multiple studies, and different brain tissues. Moreover, significant metabolic differences were observed due to the latent factor. The latent factor was found to reflect cell-type heterogeneity in the brain and after adjusting for it, we were able to identify new biological pathways. The changes observed at both transcriptomic and metabolomic levels support the importance of identifying any latent factors before pursuing downstream analysis to accurately identify biomarkers.

Background Single-cell omics technology is rapidly developing to measure the epigenome, genome, and transcriptome across a range of cell types. However, it is still challenging to integrate omics data from different modalities. Here, we propose a variation of the Siamese neural network framework called MinNet, which is trained to integrate multi-omics data on the single-cell resolution by using graph-based contrastive loss. Results By training the model and testing it on several benchmark datasets, we showed its accuracy and generalizability in integrating scRNA-seq with scATAC-seq, and scRNA-seq with epitope data. Further evaluation demonstrated our model's unique ability to remove the batch effect, a common problem in actual practice. To show how the integration impacts downstream analysis, we established model-based smoothing and cis-regulatory element-inferring method and validated it with external pcHi-C evidence. Finally, we applied the framework to a COVID-19 dataset to bolster the original work with integration-based analysis, showing its necessity in single-cell multi-omics research. Conclusions MinNet is a novel deep-learning framework for single-cell multi-omics sequencing data integration. It ranked top among other methods in benchmarking and is especially suitable for integrating datasets with batch and biological variances. With the single-cell resolution integration results, analysis of the interplay between genome and transcriptome can be done to help researchers understand their data and question.

We are pleased to introduce a first-of-its-kind algorithm that combines in-silico region detection and spatial gene-expression imputation. Spatial transcriptomics by 10X Visium (ST) is a new technology used to dissect gene and cell spatial organization. Analyzing this new type of data has two main challenges: automatically annotating the major tissue regions and excessive zero values of gene-expression due to high dropout rates. We developed a computational tool—MIST—that addresses both challenges by automatically identifying tissue regions and estimating missing gene-expression values for individual tissue regions. We validated MIST detected regions across multiple datasets using manual annotation on the histological staining images as references. We also demonstrated that MIST can accurately recover ST’s missing values through hold-out experiments. Furthermore, we showed that MIST could identify subtle intra-tissue heterogeneity and recover spatial gene-gene interaction signals. We therefore strongly encourage using MIST prior to downstream ST analysis because it provides unbiased region annotations and enables accurately de-noised spatial gene-expression profiles.

As more datasets, tools, workflows, APIs, and other digital resources are produced by the research community, it is becoming increasingly difficult to harmonize and organize these efforts for maximal synergistic integrated utilization. The Findable, Accessible, Interoperable, and Reusable (FAIR) guiding principles have prompted many stakeholders to consider strategies for tackling this challenge by making these digital resources follow common standards and best practices so that they can become more integrated and organized. Faced with the question of how to make digital resources more FAIR, it has become imperative to measure what it means to be FAIR. The diversity of resources, communities, and stakeholders have different goals and use cases and this makes assessment of FAIRness particularly challenging. To begin resolving this challenge, the FAIRshake toolkit was developed to enable the establishment of community-driven FAIR metrics and rubrics paired with manual, semi- and fully-automated FAIR assessment capabilities. The FAIRshake toolkit contains a database that lists registered digital resources, with their associated metrics, rubrics, and assessments. The FAIRshake toolkit also has a browser extension and a bookmarklet that enables viewing and submitting assessments from any website. The FAIR assessment results are visualized as an insignia that can be viewed on the FAIRshake website, or embedded within hosting websites. Using FAIRshake, a variety of bioinformatics tools, datasets listed on dbGaP, APIs registered in SmartAPI, workflows in Dockstore, and other biomedical digital resources were manually and automatically assessed for FAIRness. In each case, the assessments revealed room for improvement, which prompted enhancements that significantly upgraded FAIRness scores of several digital resources.