Single cell RNA-sequencing: A powerful yet still challenging technology to study cellular heterogeneity

Integrated analysis of single-cell and bulk RNA-sequencing identifies a signature based on drug response genes to predict prognosis and therapeutic response in ovarian cancer

Ovarian cancer represents a severe gynecological malignancy with a dire prognosis, underscoring the imperative need for dependable biomarkers that can accurately predict drug response and guide therapeutic choices. In this study, we harnessed online single-cell RNA sequencing (scRNAseq) and bulk RNA sequencing (RNAseq) datasets, applying the Scissor algorithm to identify cells responsive to paclitaxel. From these cells, we derived a gene signature, subsequently used to construct a prognostic model that demonstrated high sensitivity and specificity in predicting patient outcomes. Moreover, we conducted pathway and functional enrichment analyses to uncover potential molecular mechanisms driving the prognostic gene signature. This study illustrates the critical role of scRNAseq and bulk RNAseq in developing precise prognostic models for ovarian cancer, potentially transforming clinical decision-making.

Revealing the molecular landscape of human placenta: a systematic review and meta-analysis of single-cell RNA sequencing studies

BACKGROUND

With increasing significance of developmental programming effects associated with placental dysfunction, more investigations are devoted to improving the characterization and understanding of placental signatures in health and disease. The placenta is a transitory but dynamic organ adapting to the shifting demands of fetal development and available resources of the maternal supply throughout pregnancy. Trophoblasts (cytotrophoblasts, syncytiotrophoblasts, and extravillous trophoblasts) are placental-specific cell types responsible for the main placental exchanges and adaptations. Transcriptomic studies with single-cell resolution have led to advances in understanding the placenta's role in health and disease. These studies, however, often show discrepancies in characterization of the different placental cell types.

OBJECTIVE AND RATIONALE

We aim to review the knowledge regarding placental structure and function gained from the use of single-cell RNA sequencing (scRNAseq), followed by comparing cell-type-specific genes, highlighting their similarities and differences. Moreover, we intend to identify consensus marker genes for the various trophoblast cell types across studies. Finally, we will discuss the contributions and potential applications of scRNAseq in studying pregnancy-related diseases.

SEARCH METHODS

We conducted a comprehensive systematic literature review to identify different cell types and their functions at the human maternal-fetal interface, focusing on all original scRNAseq studies on placentas published before March 2023 and published reviews (total of 28 studies identified) using PubMed search. Our approach involved curating cell types and subtypes that had previously been defined using scRNAseq and comparing the genes used as markers or identified as potential new markers. Next, we reanalyzed expression matrices from the six available scRNAseq raw datasets with cell annotations (four from first trimester and two at term), using Wilcoxon rank-sum tests to compare gene expression among studies and annotate trophoblast cell markers in both first trimester and term placentas. Furthermore, we integrated scRNAseq raw data available from 18 healthy first trimester and nine term placentas, and performed clustering and differential gene expression analysis. We further compared markers obtained with the analysis of annotated and raw datasets with the literature to obtain a common signature gene list for major placental cell types.

OUTCOMES

Variations in the sampling site, gestational age, fetal sex, and subsequent sequencing and analysis methods were observed between the studies. Although their proportions varied, the three trophoblast types were consistently identified across all scRNAseq studies, unlike other non-trophoblast cell types. Notably, no marker genes were shared by all studies for any of the investigated cell types. Moreover, most of the newly defined markers in one study were not observed in other studies. These discrepancies were confirmed by our analysis on trophoblast cell types, where hundreds of potential marker genes were identified in each study but with little overlap across studies. From 35 461 and 23 378 cells of high quality in the first trimester and term placentas, respectively, we obtained major placental cell types, including perivascular cells that previously had not been identified in the first trimester. Importantly, our meta-analysis provides marker genes for major placental cell types based on our extensive curation.

WIDER IMPLICATIONS

This review and meta-analysis emphasizes the need for establishing a consensus for annotating placental cell types from scRNAseq data. The marker genes identified here can be deployed for defining human placental cell types, thereby facilitating and improving the reproducibility of trophoblast cell annotation.

Single-Cell Sequencing Technology in Ruminant Livestock: Challenges and Opportunities

Advancements in single-cell sequencing have transformed the genomics field by allowing researchers to delve into the intricate cellular heterogeneity within tissues at greater resolution. While single-cell omics are more widely applied in model organisms and humans, their use in livestock species is just beginning. Studies in cattle, sheep, and goats have already leveraged single-cell and single-nuclei RNA-seq as well as single-cell and single-nuclei ATAC-seq to delineate cellular diversity in tissues, track changes in cell populations and gene expression over developmental stages, and characterize immune cell populations important for disease resistance and resilience. Although challenges exist for the use of this technology in ruminant livestock, such as the precise annotation of unique cell populations and spatial resolution of cells within a tissue, there is vast potential to enhance our understanding of the cellular and molecular mechanisms underpinning traits essential for healthy and productive livestock. This review intends to highlight the insights gained from published single-cell omics studies in cattle, sheep, and goats, particularly those with publicly accessible data. Further, this manuscript will discuss the challenges and opportunities of this technology in ruminant livestock and how it may contribute to enhanced profitability and sustainability of animal agriculture in the future.

Single-Cell Omics for Transcriptome CHaracterization (SCOTCH): isoform-level characterization of gene expression through long-read single-cell RNA sequencing

The advent of long-read single-cell transcriptome sequencing (lr-scRNA-Seq) represents a significant leap forward in single-cell genomics. With the recent introduction of R10 flowcells by Oxford Nanopore, we propose that previous computational methods designed to handle high sequencing error rates are no longer relevant, and that the prevailing approach using short reads to compile "barcode space" (candidate barcode list) to de-multiplex long reads are no longer necessary. Instead, computational methods should now shift focus on harnessing the unique benefits of long reads to analyze transcriptome complexity. In this context, we introduce a comprehensive suite of computational methods named Single-Cell Omics for Transcriptome CHaracterization (SCOTCH). Our method is compatible with the single-cell library preparation platform from both 10X Genomics and Parse Biosciences, facilitating the analysis of special cell populations, such as neurons, hepatocytes and developing cardiomyocytes. We specifically re-formulated the transcript mapping problem with a compatibility matrix and addressed the multiple-mapping issue using probabilistic inference, which allows the discovery of novel isoforms as well as the detection of differential isoform usage between cell populations. We evaluated SCOTCH through analysis of real data across different combinations of single-cell libraries and sequencing technologies (10X + Illumina, Parse + Illumina, 10X + Nanopore_R9, 10X + Nanopore_R10, Parse + Nanopore_R10), and showed its ability to infer novel biological insights on cell type-specific isoform expression. These datasets enhance the availability of publicly available data for continued development of computational approaches. In summary, SCOTCH allows extraction of more biological insights from the new advancements in single-cell library construction and sequencing technologies, facilitating the examination of transcriptome complexity at the single-cell level.

Single-cell and spatial RNA sequencing reveal the spatiotemporal trajectories of fruit senescence

The senescence of fruit is a complex physiological process, with various cell types within the pericarp, making it highly challenging to elucidate their individual roles in fruit senescence. In this study, a single-cell expression atlas of the pericarp of pitaya (Hylocereus undatus) is constructed, revealing exocarp and mesocarp cells undergoing the most significant changes during the fruit senescence process. Pseudotime analysis establishes cellular differentiation and gene expression trajectories during senescence. Early-stage oxidative stress imbalance is followed by the activation of resistance in exocarp cells, subsequently senescence-associated proteins accumulate in the mesocarp cells at late-stage senescence. The central role of the early response factor HuCMB1 is unveiled in the senescence regulatory network. This study provides a spatiotemporal perspective for a deeper understanding of the dynamic senescence process in plants.

Chronic arsenic exposure induces malignant transformation of human HaCaT cells through both deterministic and stochastic changes in transcriptome expression

Biological processes are inherently stochastic, i.e., are partially driven by hard to predict random probabilistic processes. Carcinogenesis is driven both by stochastic and deterministic (predictable non-random) changes. However, very few studies systematically examine the contribution of stochastic events leading to cancer development. In differential gene expression studies, the established data analysis paradigms incentivize expression changes that are uniformly different across the experimental versus control groups, introducing preferential inclusion of deterministic changes at the expense of stochastic processes that might also play a crucial role in the process of carcinogenesis. In this study, we applied simple computational techniques to quantify: (i) The impact of chronic arsenic (iAs) exposure as well as passaging time on stochastic gene expression and (ii) Which genes were expressed deterministically and which were expressed stochastically at each of the three stages of cancer development. Using biological coefficient of variation as an empirical measure of stochasticity we demonstrate that chronic iAs exposure consistently suppressed passaging related stochastic gene expression at multiple time points tested, selecting for a homogenous cell population that undergo transformation. Employing multiple balanced removal of outlier data, we show that chronic iAs exposure induced deterministic and stochastic changes in the expression of unique set of genes, that populate largely unique biological pathways. Together, our data unequivocally demonstrate that both deterministic and stochastic changes in transcriptome-wide expression are critical in driving biological processes, pathways and networks towards clonal selection, carcinogenesis, and tumor heterogeneity.

Characterization of tumor-associated endothelial cells and the development of a prognostic model in pancreatic ductal adenocarcinoma

Pancreatic ductal adenocarcinoma (PDAC) is a highly aggressive malignancy characterized by a complex tumor microenvironment. Angiogenesis is of paramount importance in the proliferation and metastasis of PDAC. However, currently, there are no well-defined biomarkers available to guide the prognosis and treatment of PDAC. In this study, we investigated the interactions between tumor-associated endothelial cells (TAECs) and tumor cells in PDAC, and identified a specific subset of TAECs characterized by high expression of COL4A1. COL4A1+ endothelial cells interact with tumor cells through the COLLAGEN signaling pathway to promote tumor cell proliferation, migration, and invasion. We also observed activation of HOXD9 in COL4A1+ endothelial cells. Based on these findings, we developed a prognostic model called TaEMS, which accurately predicts patient prognosis. TaEMS identified high-risk patients enriched in cell cycle-related pathways and low-risk patients enriched in focal adhesions, smooth muscle regulation, and immune pathways. Moreover, high-risk patients displayed a reduced level of immune cell infiltration, indicating the presence of a "cold tumor" phenotype. Overall, our study uncovered an intricate crosstalk between TAECs and tumor cells in PDAC, emphasizing the role of HOXD9 and highlighting the potential of TaEMS as a prognostic biomarker for precise therapies.

Unveiling the immune system aging in single-cell resolution

This commentary investigates the findings presented in the article by Yang et al. in 2023, published in the Journal of Leukocyte Biology. This commentary first summarizes the spatial-temporal dynamics of regulatory T cells derived from mice (Tabula Muris Senis) of different ages (3, 18, and 24 mo) at different anatomic niches like lymph nodes and bone marrow. We also reported possible combinations of receptor-ligand interactions among T follicular regulatory cells, T follicular helper cells, and germinal center B cells, such as the calmodulin/Fas axis and PSGL-1/L-selectin axis. Then, we have elaborated on the significance of understanding aging regulatory T cells and offered some possible future research directions for Yang et al., contributing to a critical analysis of their recent study. Building on these foundations, further investigations and studies can be conducted to delve deeper into the mechanisms by which regulatory T cells influence health upon aging, potentially unveiling novel therapeutic targets to ameliorate age-related pathogenicity.

The Simple prEservatioN of Single cElls method for cryopreservation enables the generation of single-cell immune profiles from whole blood

Introduction

Current multistep methods utilized for preparing and cryopreserving single-cell suspensions from blood samples for single-cell RNA sequencing (scRNA-seq) are time-consuming, requiring trained personnel and special equipment, so limiting their clinical adoption. We developed a method, Simple prEservatioN of Single cElls (SENSE), for single-step cryopreservation of whole blood (WB) along with granulocyte depletion during single-cell assay, to generate high quality single-cell profiles (SCP).

Methods

WB was cryopreserved using the SENSE method and peripheral blood mononuclear cells (PBMCs) were isolated and cryopreserved using the traditional density-gradient method (PBMC method) from the same blood sample (n=6). The SCPs obtained from both methods were processed using a similar pipeline and quality control parameters. Further, entropy calculation, differential gene expression, and cellular communication analysis were performed to compare cell types and subtypes from both methods.

Results

Highly viable (86.3 ± 1.51%) single-cell suspensions (22,353 cells) were obtained from the six WB samples cryopreserved using the SENSE method. In-depth characterization of the scRNA-seq datasets from the samples processed with the SENSE method yielded high-quality profiles of lymphoid and myeloid cell types which were in concordance with the profiles obtained with classical multistep PBMC method processed samples. Additionally, the SENSE method cryopreserved samples exhibited significantly higher T-cell enrichment, enabling deeper characterization of T-cell subtypes. Overall, the SENSE and PBMC methods processed samples exhibited transcriptional, and cellular communication network level similarities across cell types with no batch effect except in myeloid lineage cells.

Discussion

Comparative analysis of scRNA-seq datasets obtained with the two cryopreservation methods i.e., SENSE and PBMC methods, yielded similar cellular and molecular profiles, confirming the suitability of the former method's incorporation in clinics/labs for cryopreserving and obtaining high-quality single-cells for conducting critical translational research.

Computational single cell oncology: state of the art

Single cell computational analysis has emerged as a powerful tool in the field of oncology, enabling researchers to decipher the complex cellular heterogeneity that characterizes cancer. By leveraging computational algorithms and bioinformatics approaches, this methodology provides insights into the underlying genetic, epigenetic and transcriptomic variations among individual cancer cells. In this paper, we present a comprehensive overview of single cell computational analysis in oncology, discussing the key computational techniques employed for data processing, analysis, and interpretation. We explore the challenges associated with single cell data, including data quality control, normalization, dimensionality reduction, clustering, and trajectory inference. Furthermore, we highlight the applications of single cell computational analysis, including the identification of novel cell states, the characterization of tumor subtypes, the discovery of biomarkers, and the prediction of therapy response. Finally, we address the future directions and potential advancements in the field, including the development of machine learning and deep learning approaches for single cell analysis. Overall, this paper aims to provide a roadmap for researchers interested in leveraging computational methods to unlock the full potential of single cell analysis in understanding cancer biology with the goal of advancing precision oncology. For this purpose, we also include a notebook that instructs on how to apply the recommended tools in the Preprocessing and Quality Control section.

Macrophage heterogeneity in the single-cell era: facts and artifacts

In this spotlight, we review technical issues that compromise single-cell analysis of tissue macrophages, including limited and unrepresentative yields, fragmentation and generation of remnants, and activation during tissue disaggregation. These issues may lead to a misleading definition of subpopulations of macrophages and the expression of macrophage-specific transcripts by unrelated cells. Recognition of the technical limitations of single-cell approaches is required in order to map the full spectrum of tissue-resident macrophage heterogeneity and assess its biological significance.

Methods for cell-type annotation on scRNA-seq data: A recent overview

The evolution of single-cell technology is ongoing, continually generating massive amounts of data that reveal many mysteries surrounding intricate diseases. However, their drawbacks continue to constrain us. Among these, annotating cell types in single-cell gene expressions pose a substantial challenge, despite the myriad of tools at our disposal. The rapid growth in data, resources, and tools has consequently brought about significant alterations in this area over the years. In our study, we spotlight all note-worthy cell type annotation techniques developed over the past four years. We provide an overview of the latest trends in this field, showcasing the most advanced methods in taxonomy. Our research underscores the demand for additional tools that incorporate a biological context and also predicts that the rising trend of graph neural network approaches will likely lead this research field in the coming years.

Investigation of Tumor Heterogeneity Using Integrated Single-Cell RNA Sequence Analysis to Focus on Genes Related to Breast Cancer-, EMT-, CSC-, and Metastasis-Related Markers in Patients with HER2-Positive Breast Cancer

Human epidermal growth factor receptor 2 (HER2) protein, which is characterized by the amplification of ERBB2, is a molecular target for HER2-overexpressing breast cancer. Many targeted HER2 strategies have been well developed thus far. Furthermore, intratumoral heterogeneity in HER2 cases has been observed with immunohistochemical staining and has been considered one of the reasons for drug resistance. Therefore, we conducted an integrated analysis of the breast cancer single-cell gene expression data for HER2-positive breast cancer cases from both scRNA-seq data from public datasets and data from our cohort and compared them with those for luminal breast cancer datasets. In our results, heterogeneous distribution of the expression of breast cancer-related genes (ESR1, PGR, ERBB2, and MKI67) was observed. Various gene expression levels differed at the single-cell level between the ERBB2-high group and ERBB2-low group. Moreover, molecular functions and ERBB2 expression levels differed between estrogen receptor (ER)-positive and ER-negative HER2 cases. Additionally, the gene expression levels of typical breast cancer-, CSC-, EMT-, and metastasis-related markers were also different across each patient. These results suggest that diversity in gene expression could occur not only in the presence of ERBB2 expression and ER status but also in the molecular characteristics of each patient.

Impact of dual-baculovirus infection on the Sf9 insect cell transcriptome during rAAV production using single-cell RNA-seq

The insect cell-baculovirus expression vector system (IC-BEVS) has shown to be a powerful platform to produce complex biopharmaceutical products, such as recombinant proteins and virus-like particles. More recently, IC-BEVS has also been used as an alternative to produce recombinant adeno-associated virus (rAAV). However, little is known about the variability of insect cell populations and the potential effect of heterogeneity (e.g., stochastic infection process and differences in infection kinetics) on product titer and/or quality. In this study, transcriptomics analysis of Sf9 insect cells during the production of rAAV of serotype 2 (rAAV2) using a low multiplicity of infection, dual-baculovirus system was performed via single-cell RNA-sequencing (scRNA-seq). Before infection, the principal source of variability in Sf9 insect cells was associated with the cell cycle. Over the course of infection, an increase in transcriptional heterogeneity was detected, which was linked to the expression of baculovirus genes as well as to differences in rAAV transgenes (rep, cap and gfp) expression. Noteworthy, at 24 h post-infection, only 29.4% of cells enclosed all three necessary rAAV transgenes to produce packed rAAV2 particles, indicating limitations of the dual-baculovirus system. In addition, the trajectory analysis herein performed highlighted that biological processes such as protein folding, metabolic processes, translation, and stress response have been significantly altered upon infection. Overall, this work reports the first application of scRNA-seq to the IC-BEVS and highlights significant variations in individual cells within the population, providing insight into the rational cell and process engineering toward improved rAAV2 production in IC-BEVS.

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.

Advances in single-cell RNA sequencing and its applications in cancer research

Cancers are a group of heterogeneous diseases characterized by the acquisition of functional capabilities during the transition from a normal to a neoplastic state. Powerful experimental and computational tools can be applied to elucidate the mechanisms of occurrence, progression, metastasis, and drug resistance; however, challenges remain. Bulk RNA sequencing techniques only reflect the average gene expression in a sample, making it difficult to understand tumor heterogeneity and the tumor microenvironment. The emergence and development of single-cell RNA sequencing (scRNA-seq) technologies have provided opportunities to understand subtle changes in tumor biology by identifying distinct cell subpopulations, dissecting the tumor microenvironment, and characterizing cellular genomic mutations. Recently, scRNA-seq technology has been increasingly used in cancer studies to explore tumor heterogeneity and the tumor microenvironment, which has increased the understanding of tumorigenesis and evolution. This review summarizes the basic processes and development of scRNA-seq technologies and their increasing applications in cancer research and clinical practice.

Construction of stemness gene score by bulk and single-cell transcriptome to characterize the prognosis of breast cancer

Breast cancer (BC) is a heterogeneous disease characterized by significant differences in prognosis and therapy response. Numerous prognostic tools have been developed for breast cancer. Usually these tools are based on bulk RNA-sequencing (RNA-Seq) and ignore tumor heterogeneity. Consequently, the goal of this study was to construct a single-cell level tool for predicting the prognosis of BC patients. In this study, we constructed a stemness-risk gene score (SGS) model based on single-sample gene set enrichment analysis (ssGSEA). Patients were divided into two groups based on the median SGS. Patients with a high SGS scores had a significantly worse prognosis than those with a low SGS, and these groups exhibited differences in several tumor characteristics, such as immune infiltration, gene mutations, and copy number variants. Our results indicate that the SGS is a reliable tool for predicting prognosis and response to immunotherapy in BC patients.

The impact of survivorship bias in glioblastoma research

Despite advances in the therapy of Central Nervous System (CNS) malignancies, treatment of glioblastoma (GB) poses significant challenges due to GB resistance and high recurrence rates following post-operative radio-chemotherapy. The majority of prognostic and predictive GB biomarkers are currently developed using tumour samples obtained through surgical interventions. However, the selection criteria adopted by different neurosurgeons to determine which cases are suitable for surgery make operated patients not representative of all GB cases. Particularly, geriatric and frail individuals are excluded from surgical consideration in some cancer centers. Such selection generates a survival (or selection) bias that introduces limitations, rendering the patients or data chosen for downstream analyses not representative of the entire community. In this review, we discuss the implication of survivorship bias on current and novel biomarkers for patient selection, stratification, therapy, and outcome analyses.

Deciphering the single-cell transcriptome network in keloids with intra-lesional injection of triamcinolone acetonide combined with 5-fluorouracil

Objectives

Keloid is a highly aggressive fibrotic disease resulting from excessive extracellular matrix deposition after dermal injury. Intra-lesional injection of triamcinolone acetonide (TAC) in combination with 5-fluorouracil (5-FU) is a commonly used pharmacological regimen and long-term repeated injections can achieve sustained inhibition of keloid proliferation. However, the molecular mechanisms underlying the inhibitory effect on keloids remain insufficiently investigated.

Methods and materials

This study performed single-cell RNA sequencing analysis of keloids treated with TAC+5-FU injections, keloids, and skins to explore patterns of gene expression regulation and cellular reprogramming.

Results

The results revealed that TAC+5-FU interrupted the differentiation trajectory of fibroblasts toward pro-fibrotic subtypes and induced keloid atrophy possibly by inhibiting the FGF signaling pathway in intercellular communication. It also stimulated partial fibroblasts to develop the potential for self-replication and multidirectional differentiation, which may be a possible cellular source of keloid recurrence. T cell dynamics demonstrated elevated expression of secretory globulin family members, which may be possible immunotherapeutic targets. Schwann cell populations achieved functional changes by increasing the proportion of apoptotic or senescence-associated cell populations and reducing cell clusters that promote epidermal development and fibroblast proliferation.

Conclusions

Our findings elucidated the molecular and cellular reprogramming of keloids by intra-lesional injection of TAC+5-FU, which will provide new insights to understand the mechanism of action and therapeutic targets.

RevGel-seq: instrument-free single-cell RNA sequencing using a reversible hydrogel for cell-specific barcoding

Progress in sample preparation for scRNA-seq is reported based on RevGel-seq, a reversible-hydrogel technology optimized for samples of fresh cells. Complexes of one cell paired with one barcoded bead are stabilized by a chemical linker and dispersed in a hydrogel in the liquid state. Upon gelation on ice the complexes are immobilized and physically separated without requiring nanowells or droplets. Cell lysis is triggered by detergent diffusion, and RNA molecules are captured on the adjacent barcoded beads for further processing with reverse transcription and preparation for cDNA sequencing. As a proof of concept, analysis of PBMC using RevGel-seq achieves results similar to microfluidic-based technologies when using the same original sample and the same data analysis software. In addition, a clinically relevant application of RevGel-seq is presented for pancreatic islet cells. Furthermore, characterizations carried out on cardiomyocytes demonstrate that the hydrogel technology readily accommodates very large cells. Standard analyses are in the 10,000-input cell range with the current gelation device, in order to satisfy common requirements for single-cell research. A convenient stopping point after two hours has been established by freezing at the cell lysis step, with full preservation of gene expression profiles. Overall, our results show that RevGel-seq represents an accessible and efficient instrument-free alternative, enabling flexibility in terms of experimental design and timing of sample processing, while providing broad coverage of cell types.

Complex Analysis of Single-Cell RNA Sequencing Data

Single-cell RNA sequencing (scRNA-seq) is a revolutionary tool for studying the physiology of normal and pathologically altered tissues. This approach provides information about molecular features (gene expression, mutations, chromatin accessibility, etc.) of cells, opens up the possibility to analyze the trajectories/phylogeny of cell differentiation and cell-cell interactions, and helps in discovery of new cell types and previously unexplored processes. From a clinical point of view, scRNA-seq facilitates deeper and more detailed analysis of molecular mechanisms of diseases and serves as a basis for the development of new preventive, diagnostic, and therapeutic strategies. The review describes different approaches to the analysis of scRNA-seq data, discusses the advantages and disadvantages of bioinformatics tools, provides recommendations and examples of their successful use, and suggests potential directions for improvement. We also emphasize the need for creating new protocols, including multiomics ones, for the preparation of DNA/RNA libraries of single cells with the purpose of more complete understanding of individual cells.