Single-Cell RNA Sequencing-Based Computational Analysis to Describe Disease Heterogeneity

Highlighting roles of autophagy in human diseases: a perspective from single-cell RNA sequencing analyses

Autophagy, the lysosome-driven breakdown of intracellular components, is pivotal in regulating eukaryotic cellular processes and maintaining homeostasis, making it physiologically important even under normal conditions. Cellular mechanisms involving autophagy include the response to nutrient deprivation, intracellular quality control, early development, and cell differentiation. Despite its established health significance, the role of autophagy in cancer and other diseases remains complex and not fully understood. A comprehensive understanding of autophagy is crucial to facilitate the development of novel therapies and drugs that can protect and improve human health. High-throughput technologies, such as single-cell RNA sequencing (scRNA-seq), have enabled researchers to study transcriptional landscapes at single-cell resolution, significantly advancing our knowledge of autophagy pathways across diverse physiological and pathological contexts. This review discusses the latest advances in scRNA-seq for autophagy research and highlights its potential in the molecular characterization of various diseases.

Comprehensive review on single-cell RNA sequencing: A new frontier in Alzheimer's disease research

Alzheimer's disease (AD) is a multifaceted neurodegenerative condition marked by gradual cognitive deterioration and the loss of neurons. While conventional bulk RNA sequencing techniques have shed light on AD pathology, they frequently obscure the cellular diversity within brain tissues. The advent of single-cell RNA sequencing (scRNA-seq) has transformed our capability to analyze the cellular composition of AD, allowing for the detection of unique cell populations, rare cell types, and gene expression alterations at an individual cell level. This review examines the use of scRNA-seq in AD research, focusing on its contributions to understanding cellular diversity, disease progression, and potential therapeutic targets. We discuss key technological innovations, data analysis techniques, and challenges associated with scRNA-seq in studying AD. Furthermore, we highlight recent studies that have utilized scRNA-seq to identify novel biomarkers, uncover disease-associated pathways, and elucidate the role of non-neuronal cells, such as microglia and astrocytes, in AD pathogenesis. By providing a comprehensive overview of advancements in scRNA-seq for unraveling cellular heterogeneity in AD, this review highlights the transformative impact of scRNA-seq on our comprehension of disease mechanisms and the creation of targeted treatments.

Integrating multi-omics data of childhood asthma using a deep association model

Childhood asthma is one of the most common respiratory diseases with rising mortality and morbidity. The multi-omics data is providing a new chance to explore collaborative biomarkers and corresponding diagnostic models of childhood asthma. To capture the nonlinear association of multi-omics data and improve interpretability of diagnostic model, we proposed a novel deep association model (DAM) and corresponding efficient analysis framework. First, the Deep Subspace Reconstruction was used to fuse the omics data and diagnostic information, thereby correcting the distribution of the original omics data and reducing the influence of unnecessary data noises. Second, the Joint Deep Semi-Negative Matrix Factorization was applied to identify different latent sample patterns and extract biomarkers from different omics data levels. Third, our newly proposed Deep Orthogonal Canonical Correlation Analysis can rank features in the collaborative module, which are able to construct the diagnostic model considering nonlinear correlation between different omics data levels. Using DAM, we deeply analyzed the transcriptome and methylation data of childhood asthma. The effectiveness of DAM is verified from the perspectives of algorithm performance and biological significance on the independent test dataset, by ablation experiment and comparison with many baseline methods from clinical and biological studies. The DAM-induced diagnostic model can achieve a prediction AUC of 0.912, which is higher than that of many other alternative methods. Meanwhile, relevant pathways and biomarkers of childhood asthma are also recognized to be collectively altered on the gene expression and methylation levels. As an interpretable machine learning approach, DAM simultaneously considers the non-linear associations among samples and those among biological features, which should help explore interpretative biomarker candidates and efficient diagnostic models from multi-omics data analysis for human complex diseases.

P-CSN: single-cell RNA sequencing data analysis by partial cell-specific network

Although many single-cell computational methods proposed use gene expression as input, recent studies show that replacing 'unstable' gene expression with 'stable' gene-gene associations can greatly improve the performance of downstream analysis. To obtain accurate gene-gene associations, conditional cell-specific network method (c-CSN) filters out the indirect associations of cell-specific network method (CSN) based on the conditional independence of statistics. However, when there are strong connections in networks, the c-CSN suffers from false negative problem in network construction. To overcome this problem, a new partial cell-specific network method (p-CSN) based on the partial independence of statistics is proposed in this paper, which eliminates the singularity of the c-CSN by implicitly including direct associations among estimated variables. Based on the p-CSN, single-cell network entropy (scNEntropy) is further proposed to quantify cell state. The superiorities of our method are verified on several datasets. (i) Compared with traditional gene regulatory network construction methods, the p-CSN constructs partial cell-specific networks, namely, one cell to one network. (ii) When there are strong connections in networks, the p-CSN reduces the false negative probability of the c-CSN. (iii) The input of more accurate gene-gene associations further optimizes the performance of downstream analyses. (iv) The scNEntropy effectively quantifies cell state and reconstructs cell pseudo-time.

Noninvasive detection and interpretation of gastrointestinal diseases by collaborative serum metabolite and magnetically controlled capsule endoscopy

Gastrointestinal diseases are complex diseases that occur in the gastrointestinal tract. Common gastrointestinal diseases include chronic gastritis, peptic ulcers, inflammatory bowel disease, and gastrointestinal tumors. These diseases may manifest a long course, difficult treatment, and repeated attacks.

Gastroscopy and mucosal biopsy are the gold standard methods for diagnosing gastric and duodenal diseases, but they are invasive procedures and carry risks due to the necessity of sedation and anesthesia. Recently, several new approaches have been developed, including serological examination and magnetically controlled capsule endoscopy (MGCE). However, serological markers lack lesion information, while MGCE images lack molecular information.

This study proposes combining these two technologies in a collaborative noninvasive diagnostic scheme as an alternative to the standard procedures. We introduce an interpretable framework for the clinical diagnosis of gastrointestinal diseases. Based on collected blood samples and MGCE records of patients with gastrointestinal diseases and comparisons with normal individuals, we selected serum metabolite signatures by bioinformatic analysis, captured image embedding signatures by convolutional neural networks, and inferred the location-specific associations between these signatures.

Our study successfully identified five key metabolite signatures with functional relevance to gastrointestinal disease. The combined signatures achieved discrimination AUC of 0.88. Meanwhile, the image embedding signatures showed different levels of validation and testing accuracy ranging from 0.7 to 0.9 according to different locations in the gastrointestinal tract as explained by their specific associations with metabolite signatures. Overall, our work provides a new collaborative noninvasive identification pipeline and candidate metabolite biomarkers for image auxiliary diagnosis. This method should be valuable for the noninvasive detection and interpretation of gastrointestinal and other complex diseases.

Integrating temporal single-cell gene expression modalities for trajectory inference and disease prediction

BACKGROUND

Current methods for analyzing single-cell datasets have relied primarily on static gene expression measurements to characterize the molecular state of individual cells. However, capturing temporal changes in cell state is crucial for the interpretation of dynamic phenotypes such as the cell cycle, development, or disease progression. RNA velocity infers the direction and speed of transcriptional changes in individual cells, yet it is unclear how these temporal gene expression modalities may be leveraged for predictive modeling of cellular dynamics.

RESULTS

Here, we present the first task-oriented benchmarking study that investigates integration of temporal sequencing modalities for dynamic cell state prediction. We benchmark ten integration approaches on ten datasets spanning different biological contexts, sequencing technologies, and species. We find that integrated data more accurately infers biological trajectories and achieves increased performance on classifying cells according to perturbation and disease states. Furthermore, we show that simple concatenation of spliced and unspliced molecules performs consistently well on classification tasks and can be used over more memory intensive and computationally expensive methods.

CONCLUSIONS

This work illustrates how integrated temporal gene expression modalities may be leveraged for predicting cellular trajectories and sample-associated perturbation and disease phenotypes. Additionally, this study provides users with practical recommendations for task-specific integration of single-cell gene expression modalities.

Examining the Development of Chronic Thromboembolic Pulmonary Hypertension at the Single-Cell Level

BACKGROUND

The mechanism of chronic thromboembolic pulmonary hypertension (CTEPH) is known to be multifactorial but remains incompletely understood.

METHODS

In this study, single-cell RNA sequencing, which facilitates the identification of molecular profiles of samples on an individual cell level, was applied to investigate individual cell types in pulmonary endarterectomized tissues from 5 patients with CTEPH. The order of single-cell types was then traced along the developmental trajectory of CTEPH by trajectory inference analysis, and intercellular communication was characterized by analysis of ligand-receptor pairs between cell types. Finally, comprehensive bioinformatics tools were used to analyze possible functions of branch-specific cell types and the underlying mechanisms.

RESULTS

Eleven cell types were identified, with immune-related cell types (T cells, natural killer cells, macrophages, and mast cells) distributed in the left (early) branch of the pseudotime tree, cancer stem cells, and CRISPLD2+ cells as intermediate cell types, and classic disease-related cell types (fibroblasts, smooth muscle cells, myofibroblasts, and endothelial cells) in the right (later) branch. Ligand-receptor interactions revealed close communication between macrophages and disease-related cell types as well as between smooth muscle cells and fibroblasts or endothelial cells. Moreover, the ligands and receptors were significantly enriched in key pathways such as the PI3K/Akt signaling pathway. Furthermore, highly expressed genes specific to the undefined cell type were significantly enriched in important functions associated with regulation of endoplasmic reticulum stress.

CONCLUSIONS

This single-cell RNA sequencing analysis revealed the order of single cells along a developmental trajectory in CTEPH as well as close communication between different cell types in CTEPH pathogenesis.

Single-Cell Sequencing: Biological Insight and Potential Clinical Implications in Pediatric Leukemia

Single-cell sequencing (SCS) provides high-resolution insight into the genomic, epigenomic, and transcriptomic landscape of oncohematological malignancies including pediatric leukemia, the most common type of childhood cancer. Besides broadening our biological understanding of cellular heterogeneity, sub-clonal architecture, and regulatory network of tumor cell populations, SCS can offer clinically relevant, detailed characterization of distinct compartments affected by leukemia and identify therapeutically exploitable vulnerabilities. In this review, we provide an overview of SCS studies focused on the high-resolution genomic and transcriptomic scrutiny of pediatric leukemia. Our aim is to investigate and summarize how different layers of single-cell omics approaches can expectedly support clinical decision making in the future. Although the clinical management of pediatric leukemia underwent a spectacular improvement during the past decades, resistant disease is a major cause of therapy failure. Currently, only a small proportion of childhood leukemia patients benefit from genomics-driven therapy, as 15-20% of them meet the indication criteria of on-label targeted agents, and their overall response rate falls in a relatively wide range (40-85%). The in-depth scrutiny of various cell populations influencing the development, progression, and treatment resistance of different disease subtypes can potentially uncover a wider range of driver mechanisms for innovative therapeutic interventions.

Advances in single-cell sequencing: insights from organ transplantation

Single-cell RNA sequencing (scRNA-seq) is a comprehensive technical tool to analyze intracellular and intercellular interaction data by whole transcriptional profile analysis. Here, we describe the application in biomedical research, focusing on the immune system during organ transplantation and rejection. Unlike conventional transcriptome analysis, this method provides a full map of multiple cell populations in one specific tissue and presents a dynamic and transient unbiased method to explore the progression of allograft dysfunction, starting from the stress response to final graft failure. This promising sequencing technology remarkably improves individualized organ rejection treatment by identifying decisive cellular subgroups and cell-specific interactions.

Applying artificial intelligence in the microbiome for gastrointestinal diseases: A review

For a long time, gut bacteria have been recognized for their important roles in the occurrence and progression of gastrointestinal diseases like colorectal cancer, and the ever-increasing amounts of microbiome data combined with other high-quality clinical and imaging datasets are leading the study of gastrointestinal diseases into an era of biomedical big data. The "omics" technologies used for microbiome analysis continuously evolve, and the machine learning or artificial intelligence technologies are key to extract the relevant information from microbiome data. This review intends to provide a focused summary of recent research and applications of microbiome big data and to discuss the use of artificial intelligence to combat gastrointestinal diseases.

Reconstructing gene regulatory networks in single-cell transcriptomic data analysis

Gene regulatory networks play pivotal roles in our understanding of biological processes/mechanisms at the molecular level. Many studies have developed sample-specific or cell-type-specific gene regulatory networks from single-cell transcriptomic data based on a large amount of cell samples. Here, we review the state-of-the-art computational algorithms and describe various applications of gene regulatory networks in biological studies.

Understanding human gut diseases at single-cell resolution

Our understanding of gut functioning and pathophysiology has grown considerably in the past decades, and advancing technologies enable us to deepen this understanding. Single-cell RNA sequencing (scRNA-seq) has opened a new realm of cellular diversity and transcriptional variation in the human gut at a high, single-cell resolution. ScRNA-seq has pushed the science of the digestive system forward by characterizing the function of distinct cell types within complex intestinal cellular environments, by illuminating the heterogeneity within specific cell populations and by identifying novel cell types in the human gut that could contribute to a variety of intestinal diseases. In this review, we highlight recent discoveries made with scRNA-seq that significantly advance our understanding of the human gut both in health and across the spectrum of gut diseases, including inflammatory bowel disease, colorectal carcinoma and celiac disease.

The Application of Single-Cell RNA Sequencing in Vaccinology

Single-cell RNA sequencing allows highly detailed profiling of cellular immune responses from limited-volume samples, advancing prospects of a new era of systems immunology. The power of single-cell RNA sequencing offers various opportunities to decipher the immune response to infectious diseases and vaccines. Here, we describe the potential uses of single-cell RNA sequencing methods in prophylactic vaccine development, concentrating on infectious diseases including COVID-19. Using examples from several diseases, we review how single-cell RNA sequencing has been used to evaluate the immunological response to different vaccine platforms and regimens. By highlighting published and unpublished single-cell RNA sequencing studies relevant to vaccinology, we discuss some general considerations how the field could be enriched with the widespread adoption of this technology.

Computational Methods for Single-Cell RNA Sequencing

Single-cell RNA sequencing (scRNA-seq) has provided a high-dimensional catalog of millions of cells across species and diseases. These data have spurred the development of hundreds of computational tools to derive novel biological insights. Here, we outline the components of scRNA-seq analytical pipelines and the computational methods that underlie these steps. We describe available methods, highlight well-executed benchmarking studies, and identify opportunities for additional benchmarking studies and computational methods. As the biochemical approaches for single-cell omics advance, we propose coupled development of robust analytical pipelines suited for the challenges that new data present and principled selection of analytical methods that are suited for the biological questions to be addressed.

Comparison of Scanpy-based algorithms to remove the batch effect from single-cell RNA-seq data

With the development of single-cell RNA sequencing (scRNA-seq) technology, analysts need to integrate hundreds of thousands of cells with multiple experimental batches. It is becoming increasingly difficult for users to select the best integration methods to remove batch effects. Here, we compared the advantages and limitations of four commonly used Scanpy-based batch-correction methods using two representative and large-scale scRNA-seq datasets. We quantitatively evaluated batch-correction performance and efficiency. Furthermore, we discussed the performance differences among the evaluated methods at the algorithm level.

Essential gene expression pattern of head and neck squamous cell carcinoma revealed by tumor-specific expression rule based on single-cell RNA sequencing

Head and neck squamous cell carcinoma (HNSCC) has been widely reported and considered as one of the most threatening diseases to human health. Derived from complicated tissue subtypes, HNSCC has diverse symptoms and pathogenesis. They make the identification of the core carcinogenic factors of such diseases at the multi-cell level difficult. With the development of single-cell sequencing technologies, the effects of non-malignant cells on traditional bulk sequencing data can be eliminated directly. On the basis of fresh single-cell RNA-seq data, we set up a computational filtering strategy for tumor cell identification in an expression rule manner. This strategy can reveal the accurate expression distinction between tumor cells and adjacent tumor microenvironment, which are all supported by literature reports. Validated by several independent datasets, these rule genes can further group HNSCC patients with significant difference on survival risks. Thus, the establishment of our computational approach may not only provide an efficient tool to identify malignant cells in the tumor ecosystem but also deepen our understanding of tumor heterogeneity and tumorigenesis.

Temporal Dynamics and Heterogeneity of Cell Populations during Skeletal Muscle Regeneration

Mammalian skeletal muscle possesses a unique ability to regenerate, which is primarily mediated by a population of resident muscle stem cells (MuSCs) and requires a concerted response from other supporting cell populations. Previous targeted analysis has described the involvement of various specific populations in regeneration, but an unbiased and simultaneous evaluation of all cell populations has been limited. Therefore, we used single-cell RNA-sequencing to uncover gene expression signatures of over 53,000 individual cells during skeletal muscle regeneration. Cells clustered into 25 populations and subpopulations, including a subpopulation of immune gene enriched myoblasts (immunomyoblasts) and subpopulations of fibro-adipogenic progenitors. Our analyses also uncovered striking spatiotemporal dynamics in gene expression, population composition, and cell-cell interaction during muscle regeneration. These findings provide insights into the cellular and molecular underpinning of skeletal muscle regeneration.

Rainbow Kaposi's Sarcoma-Associated Herpesvirus Revealed Heterogenic Replication with Dynamic Gene Expression

Molecular mechanisms of Kaposi's sarcoma-associated herpesvirus (KSHV) reactivation have been studied primarily by measuring the total or average activity of an infected cell population, which often consists of a mixture of both nonresponding and reactivating cells that in turn contain KSHVs at various stages of replication. Studies on KSHV gene regulation at the individual cell level would allow us to better understand the basis for this heterogeneity, and new preventive measures could be developed based on findings from nonresponding cells exposed to reactivation stimuli. Here, we generated a recombinant reporter virus, which we named "Rainbow-KSHV," that encodes three fluorescence-tagged KSHV proteins (mBFP2-ORF6, mCardinal-ORF52, and mCherry-LANA). Rainbow-KSHV replicated similarly to a prototype reporter-KSHV, KSHVr.219, and wild-type BAC16 virus. Live imaging revealed unsynchronized initiation of reactivation and KSHV replication with diverse kinetics between individual cells. Cell fractionation revealed temporal gene regulation, in which early lytic gene expression was terminated in late protein-expressing cells. Finally, isolation of fluorescence-positive cells from nonresponders increased dynamic ranges of downstream experiments 10-fold. Thus, this study demonstrates a tool to examine heterogenic responses of KSHV reactivation for a deeper understanding of KSHV replication.IMPORTANCE Sensitivity and resolution of molecular analysis are often compromised by the use of techniques that measure the ensemble average of large cell populations. Having a research tool to nondestructively identify the KSHV replication stage in an infected cell would not only allow us to effectively isolate cells of interest from cell populations but also enable more precise sample selection for advanced single-cell analysis. We prepared a recombinant KSHV that can report on its replication stage in host cells by differential fluorescence emission. Consistent with previous host gene expression studies, our experiments reveal the highly heterogenic nature of KSHV replication/gene expression at individual cell levels. The utilization of a newly developed reporter-KSHV and initial characterization of KSHV replication in single cells are presented.

Clinical implications of intratumor heterogeneity: challenges and opportunities

In this review, we highlight the role of intratumoral heterogeneity, focusing on the clinical and biological ramifications this phenomenon poses. Intratumoral heterogeneity arises through complex genetic, epigenetic, and protein modifications that drive phenotypic selection in response to environmental pressures. Functionally, heterogeneity provides tumors with significant adaptability. This ranges from mutual beneficial cooperation between cells, which nurture features such as growth and metastasis, to the narrow escape and survival of clonal cell populations that have adapted to thrive under specific conditions such as hypoxia or chemotherapy. These dynamic intercellular interplays are guided by a Darwinian selection landscape between clonal tumor cell populations and the tumor microenvironment. Understanding the involved drivers and functional consequences of such tumor heterogeneity is challenging but also promises to provide novel insight needed to confront the problem of therapeutic resistance in tumors.