rCASC: reproducible classification analysis of single-cell sequencing data

Single-cell transcriptomics analysis of bullous pemphigoid unveils immune-stromal crosstalk in type 2 inflammatory disease

Bullous pemphigoid (BP) is a type 2 inflammation- and immunity-driven skin disease, yet a comprehensive understanding of the immune landscape, particularly immune-stromal crosstalk in BP, remains elusive. Herein, using single-cell RNA sequencing (scRNA-seq) and in vitro functional analyzes, we pinpoint Th2 cells, dendritic cells (DCs), and fibroblasts as crucial cell populations. The IL13-IL13RA1 ligand-receptor pair is identified as the most significant mediator of immune-stromal crosstalk in BP. Notably, fibroblasts and DCs expressing IL13RA1 respond to IL13-secreting Th2 cells, thereby amplifying Th2 cell-mediated cascade responses, which occurs through the specific upregulation of PLA2G2A in fibroblasts and CCL17 in myeloid cells, creating a positive feedback loop integral to immune-stromal crosstalk. Furthermore, PLA2G2A and CCL17 contribute to an increased titer of pathogenic anti-BP180-NC16A autoantibodies in BP patients. Our work provides a comprehensive insight into BP pathogenesis and shows a mechanism governing immune-stromal interactions, providing potential avenues for future therapeutic research.

A single cell RNAseq benchmark experiment embedding "controlled" cancer heterogeneity

Single-cell RNA sequencing (scRNA-seq) has emerged as a vital tool in tumour research, enabling the exploration of molecular complexities at the individual cell level. It offers new technical possibilities for advancing tumour research with the potential to yield significant breakthroughs. However, deciphering meaningful insights from scRNA-seq data poses challenges, particularly in cell annotation and tumour subpopulation identification. Efficient algorithms are therefore needed to unravel the intricate biological processes of cancer. To address these challenges, benchmarking datasets are essential to validate bioinformatics methodologies for analysing single-cell omics in oncology. Here, we present a 10XGenomics scRNA-seq experiment, providing a controlled heterogeneous environment using lung cancer cell lines characterised by the expression of seven different driver genes (EGFR, ALK, MET, ERBB2, KRAS, BRAF, ROS1), leading to partially overlapping functional pathways. Our dataset provides a comprehensive framework for the development and validation of methodologies for analysing cancer heterogeneity by means of scRNA-seq.

p140Cap inhibits β-Catenin in the breast cancer stem cell compartment instructing a protective anti-tumor immune response

The p140Cap adaptor protein is a tumor suppressor in breast cancer associated with a favorable prognosis. Here we highlight a function of p140Cap in orchestrating local and systemic tumor-extrinsic events that eventually result in inhibition of the polymorphonuclear myeloid-derived suppressor cell function in creating an immunosuppressive tumor-promoting environment in the primary tumor, and premetastatic niches at distant sites. Integrative transcriptomic and preclinical studies unravel that p140Cap controls an epistatic axis where, through the upstream inhibition of β-Catenin, it restricts tumorigenicity and self-renewal of tumor-initiating cells limiting the release of the inflammatory cytokine G-CSF, required for polymorphonuclear myeloid-derived suppressor cells to exert their local and systemic tumor conducive function. Mechanistically, p140Cap inhibition of β-Catenin depends on its ability to localize in and stabilize the β-Catenin destruction complex, promoting enhanced β-Catenin inactivation. Clinical studies in women show that low p140Cap expression correlates with reduced presence of tumor-infiltrating lymphocytes and more aggressive tumor types in a large cohort of real-life female breast cancer patients, highlighting the potential of p140Cap as a biomarker for therapeutic intervention targeting the β-Catenin/ Tumor-initiating cells /G-CSF/ polymorphonuclear myeloid-derived suppressor cell axis to restore an efficient anti-tumor immune response.

Single-Cell RNAseq Complexity Reduction

An important step in single-cell RNAseq data analysis is the preparation of the single cell transcription data for cell sub-population partitioning. In this chapter, we describe how to perform complexity reduction for 3' end single-cell RNAseq transcriptomics data.

Using "Galaxy-rCASC": A Public Galaxy Instance for Single-Cell RNA-Seq Data Analysis

rCASC is a modular workflow providing an integrated environment for single-cell RNA-seq (scRNA-Seq) data analysis exploiting Docker containers to achieve functional and computational reproducibility. It was initially developed as an R package usable also through a Java GUI. However, the Java frontend cannot be employed when running rCASC on a remote server, a typical setup due to the significant computational resources commonly needed to analyze scRNA-Seq data.To allow the use of rCASC through a graphical user interface on the client side and to harness the many advantages provided by the Galaxy platform, we have made rCASC available as a Galaxy set of tools, also providing a dedicated public instance of Galaxy named "Galaxy-rCASC." To integrate rCASC into Galaxy, all its functions, originally implemented as a set of Docker containers to maximize reproducibility, have been extensively reworked to become independent from the R package functions that launch them in the original implementation. Furthermore, suitable Galaxy wrappers have been developed for most functions of rCASC. We provide a detailed reference document to the use of Galaxy-rCASC with insights and explanations on the platform functionalities, parameters, and output while guiding the reader through the typical rCASC analysis workflow of a scRNA-Seq dataset.

Guidance on Processing the 10x Genomics Single Cell Gene Expression Assay

The demand for technologies that allow the study of gene expression at single cell resolution continues to increase. One such assay was launched in 2016 by the US-based company 10x Genomics Inc. Utilizing the power of the single cell on a large scale (Zheng et al. Nat Commun 8:14049, 2017)-capturing thousands of cells at once-has shaped life sciences ever since and allowed researchers to discover new insights within their respective fields of study such as oncology, neurobiology, and immunology (among others). Obtaining high-data quality is the key to being able to make these meaningful discoveries, which in turn is directly linked to the quality of the initial cell (or nuclei) suspension that is used to load the 10x Genomics Chromium Single Cell Gene Expression assay. A successful workflow relies on a cell suspension which is fully dissociated, extremely clean, and of high viability. While the workflow itself has been detailed elsewhere (De Simone et al. Methods Mol Biol 1979:87-110, 2019), in this chapter we will focus on the importance of the quality of the initial cell suspension, as well as common mistakes that can occur while running a Single Cell Gene Expression assay. The descriptions of these tips and tricks refer to the current version of the 10x Genomics User Guide (Chromium Single Cell 3' Reagent Kits User Guide (v3.1 Chemistry Dual Index). https://support10xgenomicscom/single-cell-gene-expression/index/doc/user-guide-chromium-single-cell-3-reagent-kits-user-guide-v31-chemistry-dual-index) which can be downloaded from the Support section on the 10x Genomics website (10x Genomics website. https://www10xgenomicscom). These documents and user guides are continuously improved and updated; hence, it is important to regularly check the company's website for the most recent version.

Functional-Feature-Based Data Reduction Using Sparsely Connected Autoencoders

Single-cell RNA sequencing (scRNA-seq) allows for the creation of large collections of individual cells transcriptome. Unsupervised clustering is an essential element for the analysis of these data, and it represents the initial step for the identification of different cell types to investigate the cell subpopulation structure of a biological sample. However, it is possible that the clustering aggregation features do not perfectly match the underlying biology since scRNA-seq data are characterized by high noise. In this chapter, we describe a functional feature-driven data reduction approach, which could provide a better link among cell clusters and their underlying cell biology.

Single-Cell RNAseq Data QC and Preprocessing

The first step in single-cell RNAseq data analysis is the evaluation of the overall quality of the cell transcriptome and the preparation of the single-cell transcription data for clustering. In this chapter, we describe one of the possible approaches to perform single-cell data preprocessing for 3' end single-cell RNAseq transcriptomics data.

Single-Cell RNAseq Clustering

Single-cell RNA sequencing (scRNA-seq) allows the creation of large collections of individual cells transcriptome. Unsupervised clustering is an essential element for the analysis of these data, and it represents the initial step for the identification of different cell types to investigate the cell subpopulation organization of a sample. In this chapter, we describe how to approach the clustering of single-cell RNAseq transcriptomics data using various clustering tools, and we provide some information on the limitations affecting the clustering procedure.

Tissue RNA Integrity in Visium Spatial Protocol (Fresh Frozen Samples)

The transcriptome of a tissue can be acquired both by single-cell RNAseq (scRNA-seq) and by spatial transcriptomics (ST). The dissociation step, which is mandatory in scRNA-seq methods, might lead to the loss of fragile cells and of spatial information, thus limiting the acquisition of the tissue cellular organization. Spatial transcriptomics methods moderate the above-mentioned issues and provide single-cell transcripts detection over an intact fresh frozen tissue section. Visium platform, commercialized from 10× Genomics, provides a whole transcriptome spatial transcriptomics platform, which does not require dedicated instruments, other than those available in any pathology laboratory. In spatial transcriptomics, proper tissue handling is mandatory to preserve the morphological quality of the tissue sections and the integrity of mRNA transcripts. Proper tissue handling is critical for downstream library preparation and sequencing performance. In this chapter, we describe the most critical steps of Visium protocol on fresh frozen tissues and we provide indications on how to interpret the data obtained from the quality control analysis recommended during the workflow.

Bringing Cell Subpopulation Discovery on a Cloud-HPC Using rCASC and StreamFlow

The idea behind novel single-cell RNA sequencing (scRNA-seq) pipelines is to isolate single cells through microfluidic approaches and generate sequencing libraries in which the transcripts are tagged to track their cell of origin. Modern scRNA-seq platforms are capable of analyzing up to many thousands of cells in each run. Then, combined with massive high-throughput sequencing producing billions of reads, scRNA-seq allows the assessment of fundamental biological properties of cell populations and biological systems at unprecedented resolution.In this chapter, we describe how cell subpopulation discovery algorithms, integrated into rCASC, could be efficiently executed on cloud-HPC infrastructure. To achieve this task, we focus on the StreamFlow framework which provides container-native runtime support for scientific workflows in cloud/HPC environments.

Identifying Gene Markers Associated with Cell Subpopulations

An important point of the analysis of a single-cell RNA experiment is the identification of the key elements, i.e., genes, characterizing each cell subpopulation cluster. In this chapter, we describe the use of sparsely connected autoencoder, as a tool to convert single-cell clusters in pseudo-RNAseq experiments to be used as input for differential expression analysis, and the use of COMET, as a tool to depict cluster-specific gene markers.

Stardust: improving spatial transcriptomics data analysis through space-aware modularity optimization-based clustering

Background

Spatial transcriptomics (ST) combines stained tissue images with spatially resolved high-throughput RNA sequencing. The spatial transcriptomic analysis includes challenging tasks like clustering, where a partition among data points (spots) is defined by means of a similarity measure. Improving clustering results is a key factor as clustering affects subsequent downstream analysis. State-of-the-art approaches group data by taking into account transcriptional similarity and some by exploiting spatial information as well. However, it is not yet clear how much the spatial information combined with transcriptomics improves the clustering result.

Results

We propose a new clustering method, Stardust, that easily exploits the combination of space and transcriptomic information in the clustering procedure through a manual or fully automatic tuning of algorithm parameters. Moreover, a parameter-free version of the method is also provided where the spatial contribution depends dynamically on the expression distances distribution in the space. We evaluated the proposed methods results by analyzing ST data sets available on the 10x Genomics website and comparing clustering performances with state-of-the-art approaches by measuring the spots' stability in the clusters and their biological coherence. Stability is defined by the tendency of each point to remain clustered with the same neighbors when perturbations are applied.

Conclusions

Stardust is an easy-to-use methodology allowing to define how much spatial information should influence clustering on different tissues and achieving more stable results than state-of-the-art approaches.

Deep learning-based advances and applications for single-cell RNA-sequencing data analysis

The rapid development of single-cell RNA-sequencing (scRNA-seq) technology has raised significant computational and analytical challenges. The application of deep learning to scRNA-seq data analysis is rapidly evolving and can overcome the unique challenges in upstream (quality control and normalization) and downstream (cell-, gene- and pathway-level) analysis of scRNA-seq data. In the present study, recent advances and applications of deep learning-based methods, together with specific tools for scRNA-seq data analysis, were summarized. Moreover, the future perspectives and challenges of deep-learning techniques regarding the appropriate analysis and interpretation of scRNA-seq data were investigated. The present study aimed to provide evidence supporting the biomedical application of deep learning-based tools and may aid biologists and bioinformaticians in navigating this exciting and fast-moving area.

Sparsely Connected Autoencoders: A Multi-Purpose Tool for Single Cell omics Analysis

Background: Biological processes are based on complex networks of cells and molecules. Single cell multi-omics is a new tool aiming to provide new incites in the complex network of events controlling the functionality of the cell. Methods: Since single cell technologies provide many sample measurements, they are the ideal environment for the application of Deep Learning and Machine Learning approaches. An autoencoder is composed of an encoder and a decoder sub-model. An autoencoder is a very powerful tool in data compression and noise removal. However, the decoder model remains a black box from which is impossible to depict the contribution of the single input elements. We have recently developed a new class of autoencoders, called Sparsely Connected Autoencoders (SCA), which have the advantage of providing a controlled association among the input layer and the decoder module. This new architecture has the benefit that the decoder model is not a black box anymore and can be used to depict new biologically interesting features from single cell data. Results: Here, we show that SCA hidden layer can grab new information usually hidden in single cell data, like providing clustering on meta-features difficult, i.e. transcription factors expression, or not technically not possible, i.e. miRNA expression, to depict in single cell RNAseq data. Furthermore, SCA representation of cell clusters has the advantage of simulating a conventional bulk RNAseq, which is a data transformation allowing the identification of similarity among independent experiments. Conclusions: In our opinion, SCA represents the bioinformatics version of a universal “Swiss-knife” for the extraction of hidden knowledgeable features from single cell omics data.

Laniakea@ReCaS: exploring the potential of customisable Galaxy on-demand instances as a cloud-based service

BACKGROUND

Improving the availability and usability of data and analytical tools is a critical precondition for further advancing modern biological and biomedical research. For instance, one of the many ramifications of the COVID-19 global pandemic has been to make even more evident the importance of having bioinformatics tools and data readily actionable by researchers through convenient access points and supported by adequate IT infrastructures. One of the most successful efforts in improving the availability and usability of bioinformatics tools and data is represented by the Galaxy workflow manager and its thriving community. In 2020 we introduced Laniakea, a software platform conceived to streamline the configuration and deployment of "on-demand" Galaxy instances over the cloud. By facilitating the set-up and configuration of Galaxy web servers, Laniakea provides researchers with a powerful and highly customisable platform for executing complex bioinformatics analyses. The system can be accessed through a dedicated and user-friendly web interface that allows the Galaxy web server's initial configuration and deployment.

RESULTS

"Laniakea@ReCaS", the first instance of a Laniakea-based service, is managed by ELIXIR-IT and was officially launched in February 2020, after about one year of development and testing that involved several users. Researchers can request access to Laniakea@ReCaS through an open-ended call for use-cases. Ten project proposals have been accepted since then, totalling 18 Galaxy on-demand virtual servers that employ ~ 100 CPUs, ~ 250 GB of RAM and ~ 5 TB of storage and serve several different communities and purposes. Herein, we present eight use cases demonstrating the versatility of the platform.

CONCLUSIONS

During this first year of activity, the Laniakea-based service emerged as a flexible platform that facilitated the rapid development of bioinformatics tools, the efficient delivery of training activities, and the provision of public bioinformatics services in different settings, including food safety and clinical research. Laniakea@ReCaS provides a proof of concept of how enabling access to appropriate, reliable IT resources and ready-to-use bioinformatics tools can considerably streamline researchers' work.

Single-cell biology to decode the immune cellular composition of kidney inflammation

Single-cell biology is transforming the ability of researchers to understand cellular signaling and identity across medical and biological disciplines. Especially for immune-mediated diseases, a single-cell look at immune cell subtypes, signaling, and activity might yield fundamental insights into the disease etiology, mechanisms, and potential therapeutic interventions. In this review, we highlight recent advances in the field of single-cell RNA profiling and their application to understand renal function in health and disease. With a focus on the immune system, in particular on T cells, we propose some key directions of understanding renal inflammation using single-cell approaches. We detail the benefits and shortcomings of the various technological approaches outlined and give advice on potential pitfalls and challenges in experimental setup and computational analysis. Finally, we conclude with a brief outlook into a promising future for single-cell technologies to elucidate kidney function.

Sparsely-connected autoencoder (SCA) for single cell RNAseq data mining

Single-cell RNA sequencing (scRNAseq) is an essential tool to investigate cellular heterogeneity. Thus, it would be of great interest being able to disclose biological information belonging to cell subpopulations, which can be defined by clustering analysis of scRNAseq data. In this manuscript, we report a tool that we developed for the functional mining of single cell clusters based on Sparsely-Connected Autoencoder (SCA). This tool allows uncovering hidden features associated with scRNAseq data. We implemented two new metrics, QCC (Quality Control of Cluster) and QCM (Quality Control of Model), which allow quantifying the ability of SCA to reconstruct valuable cell clusters and to evaluate the quality of the neural network achievements, respectively. Our data indicate that SCA encoded space, derived by different experimentally validated data (TF targets, miRNA targets, Kinase targets, and cancer-related immune signatures), can be used to grasp single cell cluster-specific functional features. In our implementation, SCA efficacy comes from its ability to reconstruct only specific clusters, thus indicating only those clusters where the SCA encoding space is a key element for cells aggregation. SCA analysis is implemented as module in rCASC framework and it is supported by a GUI to simplify it usage for biologists and medical personnel.

Regenerative Reprogramming of the Intestinal Stem Cell State via Hippo Signaling Suppresses Metastatic Colorectal Cancer

Although the Hippo transcriptional coactivator YAP is considered oncogenic in many tissues, its roles in intestinal homeostasis and colorectal cancer (CRC) remain controversial. Here, we demonstrate that the Hippo kinases LATS1/2 and MST1/2, which inhibit YAP activity, are required for maintaining Wnt signaling and canonical stem cell function. Hippo inhibition induces a distinct epithelial cell state marked by low Wnt signaling, a wound-healing response, and transcription factor Klf6 expression. Notably, loss of LATS1/2 or overexpression of YAP is sufficient to reprogram Lgr5+ cancer stem cells to this state and thereby suppress tumor growth in organoids, patient-derived xenografts, and mouse models of primary and metastatic CRC. Finally, we demonstrate that genetic deletion of YAP and its paralog TAZ promotes the growth of these tumors. Collectively, our results establish the role of YAP as a tumor suppressor in the adult colon and implicate Hippo kinases as therapeutic vulnerabilities in colorectal malignancies.

The intersectional genetics landscape for humans

BACKGROUND

The human body is made up of hundreds-perhaps thousands-of cell types and states, most of which are currently inaccessible genetically. Intersectional genetic approaches can increase the number of genetically accessible cells, but the scope and safety of these approaches have not been systematically assessed. A typical intersectional method acts like an "AND" logic gate by converting the input of 2 or more active, yet unspecific, regulatory elements (REs) into a single cell type specific synthetic output.

RESULTS

Here, we systematically assessed the intersectional genetics landscape of the human genome using a subset of cells from a large RE usage atlas (Functional ANnoTation Of the Mammalian genome 5 consortium, FANTOM5) obtained by cap analysis of gene expression sequencing (CAGE-seq). We developed the heuristics and algorithms to retrieve and quality-rank "AND" gate intersections. Of the 154 primary cell types surveyed, >90% can be distinguished from each other with as few as 3 to 4 active REs, with quantifiable safety and robustness. We call these minimal intersections of active REs with cell-type diagnostic potential "versatile entry codes" (VEnCodes). Each of the 158 cancer cell types surveyed could also be distinguished from the healthy primary cell types with small VEnCodes, most of which were robust to intra- and interindividual variation. Methods for the cross-validation of CAGE-seq-derived VEnCodes and for the extraction of VEnCodes from pooled single-cell sequencing data are also presented.

CONCLUSIONS

Our work provides a systematic view of the intersectional genetics landscape in humans and demonstrates the potential of these approaches for future gene delivery technologies.

rCASC: reproducible classification analysis of single-cell sequencing data

Background

Single-cell RNA sequencing is essential for investigating cellular heterogeneity and highlighting cell subpopulation-specific signatures. Single-cell sequencing applications have spread from conventional RNA sequencing to epigenomics, e.g., ATAC-seq. Many related algorithms and tools have been developed, but few computational workflows provide analysis flexibility while also achieving functional (i.e., information about the data and the tools used are saved as metadata) and computational reproducibility (i.e., a real image of the computational environment used to generate the data is stored) through a user-friendly environment.

Findings

rCASC is a modular workflow providing an integrated analysis environment (from count generation to cell subpopulation identification) exploiting Docker containerization to achieve both functional and computational reproducibility in data analysis. Hence, rCASC provides preprocessing tools to remove low-quality cells and/or specific bias, e.g., cell cycle. Subpopulation discovery can instead be achieved using different clustering techniques based on different distance metrics. Cluster quality is then estimated through the new metric "cell stability score" (CSS), which describes the stability of a cell in a cluster as a consequence of a perturbation induced by removing a random set of cells from the cell population. CSS provides better cluster robustness information than the silhouette metric. Moreover, rCASC's tools can identify cluster-specific gene signatures.

Conclusions

rCASC is a modular workflow with new features that could help researchers define cell subpopulations and detect subpopulation-specific markers. It uses Docker for ease of installation and to achieve a computation-reproducible analysis. A Java GUI is provided to welcome users without computational skills in R.