Comparison of visualization tools for single-cell RNAseq data

Exploring cellular changes in ruptured human quadriceps tendons at single-cell resolution

Tendon ruptures in humans have often been studied during the chronic phase of injury, particularly in the context of rotator cuff disease. However, the early response to acute tendon ruptures remains less investigated. Quadriceps tendons, which require prompt surgical treatment, offer a model to investigate this early response. Therefore, this study aimed to explore the early cellular changes in ruptured compared to healthy human quadriceps tendons. Quadriceps tendon samples were collected from patients undergoing fracture repair (healthy) or tendon repair surgery (collected 7-8 days post-injury). Nuclei were isolated for single-nucleus RNA sequencing, and comprehensive transcriptomic analysis was conducted. The transcriptomes of 12,808 nuclei (7268 from healthy and 5540 from ruptured quadriceps tendons) were profiled, revealing 12 major cell types and several cell subtypes and states. Rupture samples showed increased expression of genes related to extracellular matrix organisation and cell cycle signalling, and a decrease in expression of genes in lipid metabolism pathways. These changes were predominantly driven by gene expression changes in the fibroblast, vascular endothelial cell (VEC), mural cell, and macrophage populations: fibroblasts shift to an activated phenotype upon rupture and there is an increase in the proportion of capillary and dividing VECs. A diverse immune environment was observed, with a shift from homeostatic to activated macrophages following rupture. Cell-cell interactions increased in number and diversity in rupture, and primarily involved fibroblast and VEC populations. Collectively, this transcriptomic analysis suggests that fibroblasts and endothelial cells are key orchestrators of the early injury response within ruptured quadriceps tendon. KEY POINTS: Tendon ruptures in humans have regularly been studied during the chronic phase of injury, but less is known about the early injury response after acute tendon ruptures. This study explored the early cellular changes in ruptured compared to healthy human quadriceps tendons at single-cell resolution. Fibroblasts and endothelial cells seem to be the key orchestrators of the early injury response within ruptured quadriceps tendon. Therefore, these cell types are obvious targets for interventions to enhance tendon healing. Overall, this study highlights that the development of more effective therapeutic options for tendon injury requires better understanding of the cellular, extracellular, and mechanical landscape of tendon tissue.

Integration of multi-omics and single-cell transcriptome reveals mitochondrial outer membrane protein-2 (MTX-2) as a prognostic biomarker and characterizes ubiquinone metabolism in lung adenocarcinoma

Lung adenocarcinoma (LUAD) remains to be one of the most prevalent and highly invasive forms of cancer. Mitochondrial outer membrane protein-2 or Metaxin-2 (MTX2), a key regulator of mitochondrial function, has been linked to cellular bioenergetics and stress response mechanisms. However, its roles in the progression and prognosis of LUAD remain largely unexplored. This study, employed a multi-omics approach, integrating transcriptomic and clinical patient data from public databases, to evaluate the expression and prognostic relevance of MTX2 in LUAD. Single-cell RNA sequencing was utilized to further explore MTX2's role in immune infiltration and interactions within the tumor microenvironment. Additionally, we validated these findings through a series of molecular biology and functional assays. Our results demonstrated that MTX2 expression was higher in LUAD tissues compared to normal lung tissues. Elevated MTX2 levels were significantly associated with poorer overall survival in LUAD patients. Functional analyses revealed that MTX2 regulates mitochondrial bioenergetics and facilitates tumor cell proliferation. Additionally, MTX2 expression was associated with increased immune cell infiltration. A pathway analysis identified cell metabolic and tumor growth pathways regulated by MTX2, supporting its role in tumor progression. Our research identifies MTX2 as a promising prognostic biomarker and therapeutic target for LUAD. Increased expression of MTX2 promotes tumor growth by altering metabolic pathways and modulating the immune response, underscoring its potential as a new target for LUAD treatment.

Atf3 controls transitioning in female mitochondrial cardiomyopathy as identified by spatial and single-cell transcriptomics

Oxidative phosphorylation defects result in now intractable mitochondrial diseases (MD) with cardiac involvement markedly affecting prognosis. The mechanisms underlying the transition from compensation to dysfunction in response to metabolic deficiency remain unclear. Here, we used spatially resolved transcriptomics and single-nucleus RNA sequencing (snRNA-seq) on the heart of a patient with mitochondrial cardiomyopathy (MCM), combined with an MCM mouse model with cardiac-specific Ndufs6 knockdown (FS6KD). Cardiomyocytes demonstrated the most heterogeneous expression landscape among cell types caused by metabolic perturbation, and pseudotime trajectory analysis revealed dynamic cellular states transitioning from compensation to severe compromise. This progression coincided with the transient up-regulation of a transcription factor, ATF3. Genetic ablation of Atf3 in FS6KD corroborated its pivotal role, effectively delaying cardiomyopathy progression in a female-specific manner. Our findings highlight a fate-determining role of ATF3 in female MCM progression and that the latest transcriptomic analysis will help decipher the mechanisms underlying MD progression.

MEF2C controls segment-specific gene regulatory networks that direct heart tube morphogenesis

The gene regulatory networks (GRNs) that control early heart formation are beginning to be understood, but lineage-specific GRNs remain largely undefined. We investigated networks controlled by the vital transcription factor MEF2C, with a time course of single-nucleus RNA- and ATAC-sequencing in wild-type and Mef2c -null embryos. We identified a "posteriorized" cardiac gene signature and chromatin landscape in the absence of MEF2C. Integrating our multiomics data in a deep learning-based model, we constructed developmental trajectories for each of the outflow tract, ventricular, and inflow tract segments, and alterations of these in Mef2c -null embryos. We computationally identified segment-specific MEF2C-dependent enhancers, with activity in the developing zebrafish heart. Finally, using inferred GRNs we discovered that the Mef2c -null heart malformations are partly driven by increased activity of the nuclear hormone receptor NR2F2. Our results delineate lineage-specific GRNs in the early heart tube and provide a generalizable framework for dissecting transcriptional networks governing developmental processes.

Infusing structural assumptions into dimensionality reduction for single-cell RNA sequencing data to identify small gene sets

Dimensionality reduction greatly facilitates the exploration of cellular heterogeneity in single-cell RNA sequencing data. While most of such approaches are data-driven, it can be useful to incorporate biologically plausible assumptions about the underlying structure or the experimental design. We propose the boosting autoencoder (BAE) approach, which combines the advantages of unsupervised deep learning for dimensionality reduction and boosting for formalizing assumptions. Specifically, our approach selects small sets of genes that explain latent dimensions. As illustrative applications, we explore the diversity of neural cell identities and temporal patterns of embryonic development.

paraCell: a novel software tool for the interactive analysis and visualization of standard and dual host-parasite single-cell RNA-seq data

Advances in sequencing technology have led to a dramatic increase in the number of single-cell transcriptomic datasets. In the field of parasitology, these datasets typically describe the gene expression patterns of a given parasite species at the single-cell level under experimental conditions, in specific hosts or tissues, or at different life cycle stages. However, while this wealth of available data represents a significant resource, analysing these datasets often requires expert computational skills, preventing a considerable proportion of the parasitology community from meaningfully integrating existing single-cell data into their work. Here, we present paraCell, a novel software tool that allows the user to visualize and analyse pre-loaded single-cell data without requiring any programming ability. The source code is free to allow remote installation. On our web server, we demonstrated how to visualize and re-analyse published Plasmodium and Trypanosoma datasets. We have also generated Toxoplasma-mouse and Theileria-cow scRNA-seq datasets to highlight the functionality of paraCell for pathogen-host interaction. The analysis of the data highlights the impact of the host interferon-γ response and gene expression profiles associated with disease susceptibility by these intracellular parasites, respectively.

Metabolic Adaptations Rewire CD4 T Cells in a Subset-Specific Manner in Human Critical Illness with and without Sepsis

Host immunity in sepsis has features of hyperinflammation together with progressive immunosuppression, particularly among CD4 T cells, that can predispose to secondary infections and ineffectual organ recovery. Metabolic and immunologic dysfunction are archetypal findings in critically ill patients with sepsis, but whether these factors are mechanistically linked remains incompletely defined. We characterized functional metabolic properties of human CD4 T cells from critically ill patients with and without sepsis and healthy adults. CD4 T cells in critical illness showed increased subset-specific metabolic plasticity, with regulatory T cells (Tregs) acquiring glycolytic capacity that stabilized suppressive markers FOXP3 and TIGIT and correlated with clinical illness severity. Single-cell transcriptomics identified differential kynurenine metabolism in Tregs, which was validated ex vivo as a mechanism of Treg glycolytic adaptation and suppressive rewiring. These findings underscore immunometabolic dysfunction as a driver of CD4 T cell remodeling in sepsis and suggest therapeutic avenues to restore an effective immune response.

Broad repression of DNA repair genes in senescent cells identified by integration of transcriptomic data

Cellular senescence plays a significant role in tissue aging. Senescent cells, which resist apoptosis while remaining metabolically active, generate endogenous DNA-damaging agents, primarily reactive oxygen species. Efficient DNA repair is therefore crucial in these cells, especially when they undergo senescence escape, resuming DNA replication and cellular proliferation. To investigate whether senescent cell transcriptomes reflect adequate DNA repair capacity, we conducted a comprehensive meta-analysis of 60 transcriptomic datasets comparing senescent to proliferating cells. Our analysis revealed a striking downregulation of genes encoding essential components across DNA repair pathways in senescent cells. This includes pathways active in different cell cycle phases such as nucleotide excision repair, base excision repair, nonhomologous end joining and homologous recombination repair of double-strand breaks, mismatch repair and interstrand crosslink repair. The downregulation observed suggests a significant accumulation of DNA lesions. Experimental monitoring of DNA repair readouts in cells that underwent radiation-induced senescence supported this conclusion. This phenomenon was consistent across various senescence triggers and was also observed in primary cell lines from aging individuals. These findings highlight the potential of senescent cells as 'ticking bombs' in aging-related diseases and tumors recurring following therapy-induced senescence.

CZI Cell Science Program, Abdulla S, Aevermann B, Assis P, Badajoz S, Bell SM, Bezzi E, Cakir B, Chaffer J, Chambers S, Cherry J, Chi T, Chien J, Dorman L, Garcia-Nieto P, Gloria N, Hastie M, Hegeman D, Hilton J, Huang T, Infeld A, Istrate AM, Jelic I, Katsuya K, Kim YJ, Liang K, Lin M, Lombardo M, Marshall B, Martin B, McDade F, Megill C, Patel N, Predeus A, Raymor B, Robatmili B, Rogers D, Rutherford E, Sadgat D, Shin A, Small C, Smith T, Sridharan P, Tarashansky A, Tavares N, Thomas H, Tolopko A, Urisko M, Yan J, Yeretssian G, Zamanian J, Mani A, Cool J, Carr A. CZ CELLxGENE Discover: a single-cell data platform for scalable exploration, analysis and modeling of aggregated data

Hundreds of millions of single cells have been analyzed using high-throughput transcriptomic methods. The cumulative knowledge within these datasets provides an exciting opportunity for unlocking insights into health and disease at the level of single cells. Meta-analyses that span diverse datasets building on recent advances in large language models and other machine-learning approaches pose exciting new directions to model and extract insight from single-cell data. Despite the promise of these and emerging analytical tools for analyzing large amounts of data, the sheer number of datasets, data models and accessibility remains a challenge. Here, we present CZ CELLxGENE Discover (cellxgene.cziscience.com), a data platform that provides curated and interoperable single-cell data. Available via a free-to-use online data portal, CZ CELLxGENE hosts a growing corpus of community-contributed data of over 93 million unique cells. Curated, standardized and associated with consistent cell-level metadata, this collection of single-cell transcriptomic data is the largest of its kind and growing rapidly via community contributions. A suite of tools and features enables accessibility and reusability of the data via both computational and visual interfaces to allow researchers to explore individual datasets, perform cross-corpus analysis, and run meta-analyses of tens of millions of cells across studies and tissues at the resolution of single cells.

Single-cell omics: experimental workflow, data analyses and applications

Cells are the fundamental units of biological systems and exhibit unique development trajectories and molecular features. Our exploration of how the genomes orchestrate the formation and maintenance of each cell, and control the cellular phenotypes of various organismsis, is both captivating and intricate. Since the inception of the first single-cell RNA technology, technologies related to single-cell sequencing have experienced rapid advancements in recent years. These technologies have expanded horizontally to include single-cell genome, epigenome, proteome, and metabolome, while vertically, they have progressed to integrate multiple omics data and incorporate additional information such as spatial scRNA-seq and CRISPR screening. Single-cell omics represent a groundbreaking advancement in the biomedical field, offering profound insights into the understanding of complex diseases, including cancers. Here, we comprehensively summarize recent advances in single-cell omics technologies, with a specific focus on the methodology section. This overview aims to guide researchers in selecting appropriate methods for single-cell sequencing and related data analysis.

Dissection of Gene Expression at the Single-Cell Level: scRNA-seq

Sequencing approaches that allowed for a better resolution of the transcriptome have been a major goal in the transcriptomics field since the development of RNA-seq techniques. While RNA-seq provides gene expression data from one entire sample in bulk, single-cell analysis allows for a better characterization of gene expression associated to specific cell types. Single-cell RNA-seq (scRNA-seq) is a reliable technique to unravel transcriptomic features of the tissues of interest dissociated at a single-cell level. The main feature of the single-cell technique is its ability to generate barcoded individual cells that allow for tracking the origin of thousands to millions of transcripts and reveal new cell types associated to diseases and different cell types and states. In this chapter, we discuss how scRNA-seq has become the gold standard to deepen the understanding of the gene expression with single-cell resolution.

Protocol for analyzing functional gene module perturbation during the progression of diseases using a single-cell Bayesian biclustering framework

The pathogenesis of complex diseases involves intricate gene regulation across cell types, necessitating a comprehensive analysis approach. Here, we present a protocol for analyzing functional gene module (FGM) perturbation during the progression of diseases using a single-cell Bayesian biclustering (scBC) framework. We describe steps for setting up the scBC workspace, preparing and exploring input data, training the model, and reconstructing the data matrix. We then detail procedures for Bayesian biclustering, exploring biclustering results, and uncovering pathway perturbations. For complete details on the use and execution of this protocol, please refer to Gong et al.1.

easySCF: a tool for enhancing interoperability between R and Python for efficient single-cell data analysis

SUMMARY

This study introduces easySCF, a tool designed to enhance the interoperability of single-cell data between the two major bioinformatics platforms, R and Python. By supporting seamless data exchange, easySCF improves the efficiency and accuracy of single-cell data analysis.

AVAILABILITY AND IMPLEMENTATION

easySCF utilizes a unified data format (.h5 format) to facilitate data transfer between R and Python platforms. The tool has been evaluated for data processing speed, memory efficiency, and disk usage, as well as its capability to handle large-scale single-cell datasets. easySCF is available as an open-source package, with implementation details and documentation accessible at https://github.com/xleizi/easySCF.

Joint representation and visualization of derailed cell states with Decipher

Biological insights often depend on comparing conditions such as disease and health, yet we lack effective computational tools for integrating single-cell genomics data across conditions or characterizing transitions from normal to deviant cell states. Here, we present Decipher, a deep generative model that characterizes derailed cell-state trajectories. Decipher jointly models and visualizes gene expression and cell state from normal and perturbed single-cell RNA-seq data, revealing shared and disrupted dynamics. We demonstrate its superior performance across diverse contexts, including in pancreatitis with oncogene mutation, acute myeloid leukemia, and gastric cancer.

CIRCLE-An interactive visual exploration tool for single cell RNA-Seq data

sCIRCLE (single-Cell Interactive Real-time Computer visualization for Low-dimensional Exploration) is a tool for exploratory analysis of single cell RNA-seq (scRNA-seq) data sets, with a focus on bacterial scRNA-seq. The software takes an information design perspective to re-envision visually and interactively exploring low dimensional representations of scRNA-Seq data. Users can project cells in various 3D and 2D spaces and interactively query and paint cells using rich metadata sets reporting on cell cluster, gene function, and gene expression. As a standalone application it contains, among other features, options for dimensionality reduction, navigation and interaction with data in 3d and 2d space, gene filtering, fold change and metacell computation as well as various capabilities for visualization, data import and export.

Single nucleus and spatial transcriptomic profiling of healthy human hamstring tendon

The molecular and cellular basis of health in human tendons remains poorly understood. Among human tendons, hamstring tendon has markedly low pathology and can provide a prototypic healthy tendon reference. The aim of this study was to determine the transcriptomes and location of all cell types in healthy hamstring tendon. Using single nucleus RNA sequencing, we profiled the transcriptomes of 10 533 nuclei from four healthy donors and identified 12 distinct cell types. We confirmed the presence of two fibroblast cell types, endothelial cells, mural cells, and immune cells, and identified cell types previously unreported in tendons, including different skeletal muscle cell types, satellite cells, adipocytes, and undefined nervous system cells. The location of these cell types within tendon was defined using spatial transcriptomics and imaging, and potential transcriptional networks and cell-cell interactions were analyzed. We demonstrate that fibroblasts have the highest number of potential cell-cell interactions in our dataset, are present throughout the tendon, and play an important role in the production and organization of extracellular matrix, thus confirming their role as key regulators of hamstring tendon homeostasis. Overall, our findings underscore the complexity of the cellular networks that underpin healthy human tendon function and the central role of fibroblasts as key regulators of hamstring tendon tissue homeostasis.

scX: A user-friendly tool for scRNA-seq exploration

Single-cell RNA sequencing (scRNA-seq) has transformed our ability to explore biological systems. Nevertheless, proficient expertise is essential for handling and interpreting the data. In this paper, we present scX, an R package built on the Shiny framework that streamlines the analysis, exploration, and visualization of single-cell experiments. With an interactive graphic interface, implemented as a web application, scX provides easy access to key scRNAseq analyses, including marker identification, gene expression profiling, and differential gene expression analysis. Additionally, scX seamlessly integrates with commonly used single-cell Seurat and Single-CellExperiment R objects, resulting in efficient processing and visualization of varied datasets. Overall, scX serves as a valuable and user-friendly tool for effortless exploration and sharing of single-cell data, simplifying some of the complexities inherent in scRNAseq analysis.

Statistical method scDEED for detecting dubious 2D single-cell embeddings and optimizing t-SNE and UMAP hyperparameters

Two-dimensional (2D) embedding methods are crucial for single-cell data visualization. Popular methods such as t-distributed stochastic neighbor embedding (t-SNE) and uniform manifold approximation and projection (UMAP) are commonly used for visualizing cell clusters; however, it is well known that t-SNE and UMAP's 2D embeddings might not reliably inform the similarities among cell clusters. Motivated by this challenge, we present a statistical method, scDEED, for detecting dubious cell embeddings output by a 2D-embedding method. By calculating a reliability score for every cell embedding based on the similarity between the cell's 2D-embedding neighbors and pre-embedding neighbors, scDEED identifies the cell embeddings with low reliability scores as dubious and those with high reliability scores as trustworthy. Moreover, by minimizing the number of dubious cell embeddings, scDEED provides intuitive guidance for optimizing the hyperparameters of an embedding method. We show the effectiveness of scDEED on multiple datasets for detecting dubious cell embeddings and optimizing the hyperparameters of t-SNE and UMAP.

Multiple Cell Lineages Give Rise to a Cell Type

It has long been assumed that a specific cell type arises following stepwise specification of cells corresponding to the branching of cell lineages. However, accumulating evidence indicates that multiple and even remote cell lineages can lead to the development of the same cells. Four examples giving different yet new insights will be discussed: skeletal muscle development from precursors with distinct initial histories of transcriptional regulation, lens cell development from remote lineages yet sharing basic transcription factors, blood cell development under intersectional pathways, and neural tissue development from cardiac precursors through the manipulation of just one component of epigenetic regulation. These examples provide flexible and nondogmatic perspectives on developmental cell regulation, fundamentally revising the old model relying on cell lineages.

scDEED: a statistical method for detecting dubious 2D single-cell embeddings and optimizing t-SNE and UMAP hyperparameters

Two-dimensional (2D) embedding methods are crucial for single-cell data visualization. Popular methods such as t-SNE and UMAP are commonly used for visualizing cell clusters; however, it is well known that t-SNE and UMAP's 2D embedding might not reliably inform the similarities among cell clusters. Motivated by this challenge, we developed a statistical method, scDEED, for detecting dubious cell embeddings output by any 2D-embedding method. By calculating a reliability score for every cell embedding, scDEED identifies the cell embeddings with low reliability scores as dubious and those with high reliability scores as trustworthy. Moreover, by minimizing the number of dubious cell embeddings, scDEED provides intuitive guidance for optimizing the hyperparameters of an embedding method. Applied to multiple scRNA-seq datasets, scDEED demonstrates its effectiveness for detecting dubious cell embeddings and optimizing the hyperparameters of t-SNE and UMAP.

A Roadmap for the Human Gut Cell Atlas

The number of studies investigating the human gastrointestinal tract using various single-cell profiling methods has increased substantially in the past few years. Although this increase provides a unique opportunity for the generation of the first comprehensive Human Gut Cell Atlas (HGCA), there remains a range of major challenges ahead. Above all, the ultimate success will largely depend on a structured and coordinated approach that aligns global efforts undertaken by a large number of research groups. In this Roadmap, we discuss a comprehensive forward-thinking direction for the generation of the HGCA on behalf of the Gut Biological Network of the Human Cell Atlas. Based on the consensus opinion of experts from across the globe, we outline the main requirements for the first complete HGCA by summarizing existing data sets and highlighting anatomical regions and/or tissues with limited coverage. We provide recommendations for future studies and discuss key methodologies and the importance of integrating the healthy gut atlas with related diseases and gut organoids. Importantly, we critically overview the computational tools available and provide recommendations to overcome key challenges.

Single cell genomics as a transformative approach for aquaculture research and innovation

Single cell genomics encompasses a suite of rapidly maturing technologies that measure the molecular profiles of individual cells within target samples. These approaches provide a large up-step in biological information compared to long-established 'bulk' methods that profile the average molecular profiles of all cells in a sample, and have led to transformative advances in understanding of cellular biology, particularly in humans and model organisms. The application of single cell genomics is fast expanding to non-model taxa, including aquaculture species, where numerous research applications are underway with many more envisaged. In this review, we highlight the potential transformative applications of single cell genomics in aquaculture research, considering barriers and potential solutions to the broad uptake of these technologies. Focusing on single cell transcriptomics, we outline considerations for experimental design, including the essential requirement to obtain high quality cells/nuclei for sequencing in ectothermic aquatic species. We further outline data analysis and bioinformatics considerations, tailored to studies with the under-characterized genomes of aquaculture species, where our knowledge of cellular heterogeneity and cell marker genes is immature. Overall, this review offers a useful source of knowledge for researchers aiming to apply single cell genomics to address biological challenges faced by the global aquaculture sector though an improved understanding of cell biology.

Transitioning single-cell genomics into the clinic

The use of genomics is firmly established in clinical practice, resulting in innovations across a wide range of disciplines such as genetic screening, rare disease diagnosis and molecularly guided therapy choice. This new field of genomic medicine has led to improvements in patient outcomes. However, most clinical applications of genomics rely on information generated from bulk approaches, which do not directly capture the genomic variation that underlies cellular heterogeneity. With the advent of single-cell technologies, research is rapidly uncovering how genomic data at cellular resolution can be used to understand disease pathology and mechanisms. Both DNA-based and RNA-based single-cell technologies have the potential to improve existing clinical applications and open new application spaces for genomics in clinical practice, with oncology, immunology and haematology poised for initial adoption. However, challenges in translating cellular genomics from research to a clinical setting must first be overcome.

scViewer: An Interactive Single-Cell Gene Expression Visualization Tool

Single-cell RNA sequencing (scRNA-seq) is an attractive technology for researchers to gain valuable insights into the cellular processes and cell type diversity present in all tissues. The data generated by the scRNA-seq experiment are high-dimensional and complex in nature. Several tools are now available to analyze the raw scRNA-seq data from public databases; however, simple and easy-to-explore single-cell gene expression visualization tools focusing on differential expression and co-expression are lacking. Here, we present scViewer, an interactive graphical user interface (GUI) R/Shiny application designed to facilitate the visualization of scRNA-seq gene expression data. With the processed Seurat RDS object as input, scViewer utilizes several statistical approaches to provide detailed information on the loaded scRNA-seq experiment and generates publication-ready plots. The major functionalities of scViewer include exploring cell-type-specific gene expression, co-expression analysis of two genes, and differential expression analysis with different biological conditions considering both cell-level and subject-level variations using negative binomial mixed modeling. We utilized a publicly available dataset (brain cells from a study of Alzheimer's disease to demonstrate the utility of our tool. scViewer can be downloaded from GitHub as a Shiny app with local installation. Overall, scViewer is a user-friendly application that will allow researchers to visualize and interpret the scRNA-seq data efficiently for multi-condition comparison by performing gene-level differential expression and co-expression analysis on the fly. Considering the functionalities of this Shiny app, scViewer can be a great resource for collaboration between bioinformaticians and wet lab scientists for faster data visualizations.

scRNASequest: an ecosystem of scRNA-seq analysis, visualization, and publishing

BACKGROUND

Single-cell RNA sequencing is a state-of-the-art technology to understand gene expression in complex tissues. With the growing amount of data being generated, the standardization and automation of data analysis are critical to generating hypotheses and discovering biological insights.

RESULTS

Here, we present scRNASequest, a semi-automated single-cell RNA-seq (scRNA-seq) data analysis workflow which allows (1) preprocessing from raw UMI count data, (2) harmonization by one or multiple methods, (3) reference-dataset-based cell type label transfer and embedding projection, (4) multi-sample, multi-condition single-cell level differential gene expression analysis, and (5) seamless integration with cellxgene VIP for visualization and with CellDepot for data hosting and sharing by generating compatible h5ad files.

CONCLUSIONS

We developed scRNASequest, an end-to-end pipeline for single-cell RNA-seq data analysis, visualization, and publishing. The source code under MIT open-source license is provided at https://github.com/interactivereport/scRNASequest . We also prepared a bookdown tutorial for the installation and detailed usage of the pipeline: https://interactivereport.github.io/scRNAsequest/tutorial/docs/ . Users have the option to run it on a local computer with a Linux/Unix system including MacOS, or interact with SGE/Slurm schedulers on high-performance computing (HPC) clusters.

siVAE: interpretable deep generative models for single-cell transcriptomes

Neural networks such as variational autoencoders (VAE) perform dimensionality reduction for the visualization and analysis of genomic data, but are limited in their interpretability: it is unknown which data features are represented by each embedding dimension. We present siVAE, a VAE that is interpretable by design, thereby enhancing downstream analysis tasks. Through interpretation, siVAE also identifies gene modules and hubs without explicit gene network inference. We use siVAE to identify gene modules whose connectivity is associated with diverse phenotypes such as iPSC neuronal differentiation efficiency and dementia, showcasing the wide applicability of interpretable generative models for genomic data analysis.

Polyphony: an Interactive Transfer Learning Framework for Single-Cell Data Analysis

Reference-based cell-type annotation can significantly reduce time and effort in single-cell analysis by transferring labels from a previously-annotated dataset to a new dataset. However, label transfer by end-to-end computational methods is challenging due to the entanglement of technical (e.g., from different sequencing batches or techniques) and biological (e.g., from different cellular microenvironments) variations, only the first of which must be removed. To address this issue, we propose Polyphony, an interactive transfer learning (ITL) framework, to complement biologists' knowledge with advanced computational methods. Polyphony is motivated and guided by domain experts' needs for a controllable, interactive, and algorithm-assisted annotation process, identified through interviews with seven biologists. We introduce anchors, i.e., analogous cell populations across datasets, as a paradigm to explain the computational process and collect user feedback for model improvement. We further design a set of visualizations and interactions to empower users to add, delete, or modify anchors, resulting in refined cell type annotations. The effectiveness of this approach is demonstrated through quantitative experiments, two hypothetical use cases, and interviews with two biologists. The results show that our anchor-based ITL method takes advantage of both human and machine intelligence in annotating massive single-cell datasets.

scAnalyzeR: A Comprehensive Software Package With Graphical User Interface for Single-Cell RNA Sequencing Analysis and its Application on Liver Cancer

Introduction: The application of single-cell RNA sequencing to delineate tissue heterogeneity and complexity has become increasingly popular. Given its tremendous resolution and high-dimensional capacity for in-depth investigation, single-cell RNA sequencing offers an unprecedented research power. Although some popular software packages are available for single-cell RNA sequencing data analysis and visualization, it is still a big challenge for their usage, as they provide only a command-line interface and require significant level of bioinformatics skills. Methods: We have developed scAnalyzeR, which is a single-cell RNA sequencing analysis pipeline with an interactive and user-friendly graphical interface for analyzing and visualizing single-cell RNA sequencing data. It accepts single-cell RNA sequencing data from various technology platforms and different model organisms (human and mouse) and allows flexibility in input file format. It provides functionalities for data preprocessing, quality control, basic summary statistics, dimension reduction, unsupervised clustering, differential gene expression, gene set enrichment analysis, correlation analysis, pseudotime cell trajectory inference, and various visualization plots. It also provides default parameters for easy usage and allows a wide range of flexibility and optimization by accepting user-defined options. It has been developed as a docker image that can be run in any docker-supported environment including Linux, Mac, and Windows, without installing any dependencies. Results: We compared the performance of scAnalyzeR with 2 other graphical tools that are popular for analyzing single-cell RNA sequencing data. The comparison was based on the comprehensiveness of functionalities, ease of usage and flexibility, and execution time. In general, scAnalyzeR outperformed the other tested counterparts in various aspects, demonstrating its superior overall performance. To illustrate the usefulness of scAnalyzeR in cancer research, we have analyzed the in-house liver cancer single-cell RNA sequencing dataset. Liver cancer tumor cells were revealed to have multiple subpopulations with distinctive gene expression signatures. Conclusion: scAnalyzeR has comprehensive functionalities and demonstrated usability. We anticipate more functionalities to be adopted in the future development.

SingleCAnalyzer: Interactive Analysis of Single Cell RNA-Seq Data on the Cloud

Single-cell RNA sequencing (scRNA-Seq) enables researchers to quantify the transcriptomes of individual cells. The capacity of researchers to perform this type of analysis has allowed researchers to undertake new scientific goals. The usefulness of scRNA-Seq has depended on the development of new computational biology methods, which have been designed to meeting challenges associated with scRNA-Seq analysis. However, the proper application of these computational methods requires extensive bioinformatics expertise. Otherwise, it is often difficult to obtain reliable and reproducible results. We have developed SingleCAnalyzer, a cloud platform that provides a means to perform full scRNA-Seq analysis from FASTQ within an easy-to-use and self-exploratory web interface. Its analysis pipeline includes the demultiplexing and alignment of FASTQ files, read trimming, sample quality control, feature selection, empty droplets detection, dimensional reduction, cellular type prediction, unsupervised clustering of cells, pseudotime/trajectory analysis, expression comparisons between groups, functional enrichment of differentially expressed genes and gene set expression analysis. Results are presented with interactive graphs, which provide exploratory and analytical features. SingleCAnalyzer is freely available at https://singleCAnalyzer.eu.

SCANNER: a web platform for annotation, visualization and sharing of single cell RNA-seq data

In recent years, efficient scRNA-seq methods have been developed, enabling the transcriptome profiling of single cells massively in parallel. Meanwhile, its high dimensionality and complexity bring challenges to the data analysis and require extensive collaborations between biologists and bioinformaticians and/or biostatisticians. The communication between these two units demands a platform for easy data sharing and exploration. Here we developed Single-Cell Transcriptomics Annotated Viewer (SCANNER), as a public web resource for the scientific community, for sharing and analyzing scRNA-seq data in a collaborative manner. It is easy-to-use without requiring special software or extensive coding skills. Moreover, it equipped a real-time database for secure data management and enables an efficient investigation of the activation of gene sets on a single-cell basis. Currently, SCANNER hosts a database of 19 types of cancers and COVID-19, as well as healthy samples from lungs of smokers and non-smokers, human brain cells and peripheral blood mononuclear cells (PBMC). The database will be frequently updated with datasets from new studies. Using SCANNER, we identified a larger proportion of cancer-associated fibroblasts cells and more active fibroblast growth-related genes in melanoma tissues in female patients compared to male patients. Moreover, we found ACE2 is mainly expressed in lung pneumocytes, secretory cells and ciliated cells and differentially expressed in lungs of smokers and never smokers.

Single-cell sequencing reveals lineage-specific dynamic genetic regulation of gene expression during human cardiomyocyte differentiation

Dynamic and temporally specific gene regulatory changes may underlie unexplained genetic associations with complex disease. During a dynamic process such as cellular differentiation, the overall cell type composition of a tissue (or an in vitro culture) and the gene regulatory profile of each cell can both experience significant changes over time. To identify these dynamic effects in high resolution, we collected single-cell RNA-sequencing data over a differentiation time course from induced pluripotent stem cells to cardiomyocytes, sampled at 7 unique time points in 19 human cell lines. We employed a flexible approach to map dynamic eQTLs whose effects vary significantly over the course of bifurcating differentiation trajectories, including many whose effects are specific to one of these two lineages. Our study design allowed us to distinguish true dynamic eQTLs affecting a specific cell lineage from expression changes driven by potentially non-genetic differences between cell lines such as cell composition. Additionally, we used the cell type profiles learned from single-cell data to deconvolve and re-analyze data from matched bulk RNA-seq samples. Using this approach, we were able to identify a large number of novel dynamic eQTLs in single cell data while also attributing dynamic effects in bulk to a particular lineage. Overall, we found that using single cell data to uncover dynamic eQTLs can provide new insight into the gene regulatory changes that occur among heterogeneous cell types during cardiomyocyte differentiation.

Immune system cells from COVID-19 patients display compromised mitochondrial-nuclear expression co-regulation and rewiring toward glycolysis

Mitochondria are pivotal for bioenergetics, as well as in cellular response to viral infections. Nevertheless, their role in COVID-19 was largely overlooked. Here, we analyzed available bulk RNA-seq datasets from COVID-19 patients and corresponding healthy controls (three blood datasets, N = 48 healthy, 119 patients; two respiratory tract datasets, N = 157 healthy, 524 patients). We found significantly reduced mtDNA gene expression in blood, but not in respiratory tract samples from patients. Next, analysis of eight single-cells RNA-seq datasets from peripheral blood mononuclear cells, nasopharyngeal samples, and Bronchoalveolar lavage fluid (N = 1,192,243 cells), revealed significantly reduced mtDNA gene expression especially in immune system cells from patients. This is associated with elevated expression of nuclear DNA-encoded OXPHOS subunits, suggesting compromised mitochondrial-nuclear co-regulation. This, together with elevated expression of ROS-response genes and glycolysis enzymes in patients, suggest rewiring toward glycolysis, thus generating beneficial conditions for SARS-CoV-2 replication. Our findings underline the centrality of mitochondrial dysfunction in COVID-19.

•

mtDNA gene expression is downregulated in COVID-19 blood, but not in respiratory tract

•

Decreased mtDNA gene expression disrupts mito-nuclear coordination

•

mtDNA is downregulated and rewired toward glycolysis especially in immune system cells

•

Mitochondrial dysfunction is central to the etiology of COVID19

UCSC Cell Browser: visualize your single-cell data

SUMMARY

As the use of single-cell technologies has grown, so has the need for tools to explore these large, complicated datasets. The UCSC Cell Browser is a tool that allows scientists to visualize gene expression and metadata annotation distribution throughout a single-cell dataset or multiple datasets.

AVAILABILITY AND IMPLEMENTATION

We provide the UCSC Cell Browser as a free website where scientists can explore a growing collection of single-cell datasets and a freely available python package for scientists to create stable, self-contained visualizations for their own single-cell datasets. Learn more at https://cells.ucsc.edu.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

RNA-Scoop: interactive visualization of transcripts in single-cell transcriptomes

Recent advances in single-cell RNA sequencing technologies have made detection of transcripts in single cells possible. The level of resolution provided by these technologies can be used to study changes in transcript usage across cell populations and help investigate new biology. Here, we introduce RNA-Scoop, an interactive cell cluster and transcriptome visualization tool to analyze transcript usage across cell categories and clusters. The tool allows users to examine differential transcript expression across clusters and investigate how usage of specific transcript expression mechanisms varies across cell groups.

CellexalVR: A virtual reality platform to visualize and analyze single-cell omics data

Single-cell RNAseq is a routinely used method to explore heterogeneity within cell populations. Data from these experiments are often visualized using dimension reduction methods such as UMAP and tSNE, where each cell is projected in two or three dimensional space. Three-dimensional projections can be more informative for larger and complex datasets because they are less prone to merging and flattening similar cell-types/clusters together. However, visualizing and cross-comparing 3D projections using current software on conventional flat-screen displays is far from optimal as they are still essentially 2D, and lack meaningful interaction between the user and the data. Here we present CellexalVR (www.cellexalvr.med.lu.se), a feature-rich, fully interactive virtual reality environment for the visualization and analysis of single-cell experiments that allows researchers to intuitively and collaboratively gain an understanding of their data.

•

Single-cell experiments are often visualized when embedded into three dimensions

•

CellexalVR is a virtual reality environment to visualize all data simultaneously

•

Teams can analyze single-cell experiments together in VR regardless of location

The Dual Role of Innate Lymphoid and Natural Killer Cells in Cancer. from Phenotype to Single-Cell Transcriptomics, Functions and Clinical Uses

Simple Summary

Innate lymphoid cells (ILCs), a family of innate immune cells including natural killers (NKs), play a multitude of roles in first-line cancer control, in escape from immunity and in cancer progression. In this review, we summarize preclinical and clinical data on ILCs and NK cells concerning their phenotype, function and clinical applications in cellular therapy trials. We also describe how single-cell transcriptome sequencing has been used and forecast how it will be used to better understand ILC and NK involvement in cancer control and progression as well as their therapeutic potential.

Abstract

The role of innate lymphoid cells (ILCs), including natural killer (NK) cells, is pivotal in inflammatory modulation and cancer. Natural killer cell activity and count have been demonstrated to be regulated by the expression of activating and inhibitory receptors together with and as a consequence of different stimuli. The great majority of NK cell populations have an anti-tumor activity due to their cytotoxicity, and for this reason have been used for cellular therapies in cancer patients. On the other hand, the recently classified helper ILCs are fundamentally involved in inflammation and they can be either helpful or harmful in cancer development and progression. Tissue niche seems to play an important role in modulating ILC function and conversion, as observed at the transcriptional level. In the past, these cell populations have been classified by the presence of specific cellular receptor markers; more recently, due to the advent of single-cell RNA sequencing (scRNA-seq), it has been possible to also explore them at the transcriptomic level. In this article we review studies on ILC (and NK cell) classification, function and their involvement in cancer. We also summarize the potential application of NK cells in cancer therapy and give an overview of the most recent studies involving ILCs and NKs at scRNA-seq, focusing on cancer. Finally, we provide a resource for those who wish to start single-cell transcriptomic analysis on the context of these innate lymphoid cell populations.

Examining heterogeneity of stromal cells in tumor microenvironment based on pan-cancer single-cell RNA sequencing data

Tumor tissues contain both tumor and non-tumor cells, which include infiltrated immune cells and stromal cells, collectively called the tumor microenvironment (TME). Single-cell RNA sequencing (scRNAseq) enables the examination of heterogeneity of tumor cells and TME. In this review, we examined scRNAseq datasets for multiple cancer types and evaluated the heterogeneity of major cell type composition in different cancer types. We further showed that endothelial cells and fibroblasts/myofibroblasts in different cancer types can be classified into common subtypes, and the subtype composition is clearly associated with cancer characteristic and therapy response.

CovidExpress: an interactive portal for intuitive investigation on SARS-CoV-2 related transcriptomes

Infection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in humans could cause coronavirus disease 2019 (COVID-19). Since its first discovery in Dec 2019, SARS-CoV-2 has become a global pandemic and caused 3.3 million direct/indirect deaths (2021 May). Amongst the scientific community's response to COVID-19, data sharing has emerged as an essential aspect of the combat against SARS-CoV-2. Despite the ever-growing studies about SARS-CoV-2 and COVID-19, to date, only a few databases were curated to enable access to gene expression data. Furthermore, these databases curated only a small set of data and do not provide easy access for investigators without computational skills to perform analyses. To fill this gap and advance open-access to the growing gene expression data on this deadly virus, we collected about 1,500 human bulk RNA-seq datasets from publicly available resources, developed a database and visualization tool, named CovidExpress (https://stjudecab.github.io/covidexpress). This open access database will allow research investigators to examine the gene expression in various tissues, cell lines, and their response to SARS-CoV-2 under different experimental conditions, accelerating the understanding of the etiology of this disease to inform the drug and vaccine development. Our integrative analysis of this big dataset highlights a set of commonly regulated genes in SARS-CoV-2 infected lung and Rhinovirus infected nasal tissues, including OASL that were under-studied in COVID-19 related reports. Our results also suggested a potential FURIN positive feedback loop that might explain the evolutional advantage of SARS-CoV-2.

Single-Cell Transcriptome Analysis as a Promising Tool to Study Pluripotent Stem Cell Reprogramming

Cells are the basic units of all organisms and are involved in all vital activities, such as proliferation, differentiation, senescence, and apoptosis. A human body consists of more than 30 trillion cells generated through repeated division and differentiation from a single-cell fertilized egg in a highly organized programmatic fashion. Since the recent formation of the Human Cell Atlas consortium, establishing the Human Cell Atlas at the single-cell level has been an ongoing activity with the goal of understanding the mechanisms underlying diseases and vital cellular activities at the level of the single cell. In particular, transcriptome analysis of embryonic stem cells at the single-cell level is of great importance, as these cells are responsible for determining cell fate. Here, we review single-cell analysis techniques that have been actively used in recent years, introduce the single-cell analysis studies currently in progress in pluripotent stem cells and reprogramming, and forecast future studies.

Understanding the Adult Mammalian Heart at Single-Cell RNA-Seq Resolution

During the last decade, extensive efforts have been made to comprehend cardiac cell genetic and functional diversity. Such knowledge allows for the definition of the cardiac cellular interactome as a reasonable strategy to increase our understanding of the normal and pathologic heart. Previous experimental approaches including cell lineage tracing, flow cytometry, and bulk RNA-Seq have often tackled the analysis of cardiac cell diversity as based on the assumption that cell types can be identified by the expression of a single gene. More recently, however, the emergence of single-cell RNA-Seq technology has led us to explore the diversity of individual cells, enabling the cardiovascular research community to redefine cardiac cell subpopulations and identify relevant ones, and even novel cell types, through their cell-specific transcriptomic signatures in an unbiased manner. These findings are changing our understanding of cell composition and in consequence the identification of potential therapeutic targets for different cardiac diseases. In this review, we provide an overview of the continuously changing cardiac cellular landscape, traveling from the pre-single-cell RNA-Seq times to the single cell-RNA-Seq revolution, and discuss the utilities and limitations of this technology.

The current landscape of single-cell transcriptomics for cancer immunotherapy

Immunotherapies such as immune checkpoint blockade and adoptive cell transfer have revolutionized cancer treatment, but further progress is hindered by our limited understanding of tumor resistance mechanisms. Emerging technologies now enable the study of tumors at the single-cell level, providing unprecedented high-resolution insights into the genetic makeup of the tumor microenvironment and immune system that bulk genomics cannot fully capture. Here, we highlight the recent key findings of the use of single-cell RNA sequencing to deconvolute heterogeneous tumors and immune populations during immunotherapy. Single-cell RNA sequencing has identified new crucial factors and cellular subpopulations that either promote tumor progression or leave tumors vulnerable to immunotherapy. We anticipate that the strategic use of single-cell analytics will promote the development of the next generation of successful, rationally designed immunotherapeutics.

Dimensionality Reduction of Single-Cell RNA-Seq Data

Dimensionality reduction is a crucial step in essentially every single-cell RNA-sequencing (scRNA-seq) analysis. In this chapter, we describe the typical dimensionality reduction workflow that is used for scRNA-seq datasets, specifically highlighting the roles of principal component analysis, t-distributed stochastic neighborhood embedding, and uniform manifold approximation and projection in this setting. We particularly emphasize efficient computation; the software implementations used in this chapter can scale to datasets with millions of cells.

Toward a Consensus View of Mammalian Adipocyte Stem and Progenitor Cell Heterogeneity

White adipose tissue (WAT) is a cellularly heterogeneous endocrine organ that not only serves as an energy reservoir, but also actively participates in metabolic homeostasis. Among the main constituents of adipose tissue are adipocytes, which arise from adipose stem and progenitor cells (ASPCs). While it is well known that these ASPCs reside in the stromal vascular fraction (SVF) of adipose tissue, their molecular heterogeneity and functional diversity is still poorly understood. Driven by the resolving power of single-cell transcriptomics, several recent studies provided new insights into the cellular complexity of ASPCs among different mammalian fat depots. In this review, we present current knowledge on ASPCs, their population structure, hierarchy, fat depot-specific nature, function, and regulatory mechanisms, and discuss not only the similarities, but also the differences between mouse and human ASPC biology.