The Impact of Heterogeneity on Single-Cell Sequencing

Single-cell RNA-seq reveals the effects of the FecB mutation on the transcriptome profile in ovine cumulus cells

Genetic variations in the ovine ovulation rate, which are associated with the FecB mutation, provide useful models by which to explore the mechanisms regulating the development of mammalian antral follicles. In order to study the effects of the FecB mutation on cumulus cell differentiation, preovulatory follicles were aspirated and cumulus cells were isolated from three FecB genotypes (homozygous, heterozygous and wild type) of Small Tail Han (STH) sheep superstimulated with FSH. Transcriptome information from tens of thousands of cumulus cells was determined with the 10 × Genomics single-cell RNA-seq technology. Under the superovulation treatment, the observed number of preovulatory follicles in the ovaries of FecB carriers was still significantly higher than that in the wild-type (P < 0.05). The expression patterns of cumulus cells differed between FecB carriers and wild-type ewes. The screened cumulus cells could also be further divided into different cell clusters, and the differentiation states and fates of each group of cumulus cells also remained different, which supports the notion that heterogeneity in gene expression is prevalent in single cells. The oxidative phosphorylation pathway was significantly enriched in differentially expressed genes among the cell differentiation branch nodes of cumulus cells and among the differentially expressed genes of cumulus cells from the three genotypes. Combined with the important role of oxidative phosphorylation in the maturation of COCs, we suggest that the oxidative phosphorylation pathway of cumulus cells plays a crucial role in the differentiation process of cumulus cells and the mutation effect of the FecB gene.

Advances in single-cell long-read sequencing technologies

With an increase in accuracy and throughput of long-read sequencing technologies, they are rapidly being assimilated into the single-cell sequencing pipelines. For transcriptome sequencing, these techniques provide RNA isoform-level information in addition to the gene expression profiles. Long-read sequencing technologies not only help in uncovering complex patterns of cell-type specific splicing, but also offer unprecedented insights into the origin of cellular complexity and thus potentially new avenues for drug development. Additionally, single-cell long-read DNA sequencing enables high-quality assemblies, structural variant detection, haplotype phasing, resolving high-complexity regions, and characterization of epigenetic modifications. Given that significant progress has primarily occurred in single-cell RNA isoform sequencing (scRiso-seq), this review will delve into these advancements in depth and highlight the practical considerations and operational challenges, particularly pertaining to downstream analysis. We also aim to offer a concise introduction to complementary technologies for single-cell sequencing of the genome, epigenome and epitranscriptome. We conclude by identifying certain key areas of innovation that may drive these technologies further and foster more widespread application in biomedical science.

A spatial transcriptome map of the developing maize ear

A comprehensive understanding of inflorescence development is crucial for crop genetic improvement, as inflorescence meristems give rise to reproductive organs and determine grain yield. However, dissecting inflorescence development at the cellular level has been challenging owing to a lack of specific marker genes to distinguish among cell types, particularly in different types of meristems that are vital for organ formation. In this study, we used spatial enhanced resolution omics-sequencing (Stereo-seq) to construct a precise spatial transcriptome map of the developing maize ear primordium, identifying 12 cell types, including 4 newly defined cell types found mainly in the inflorescence meristem. By extracting the meristem components for detailed clustering, we identified three subtypes of meristem and validated two MADS-box genes that were specifically expressed at the apex of determinate meristems and involved in stem cell determinacy. Furthermore, by integrating single-cell RNA transcriptomes, we identified a series of spatially specific networks and hub genes that may provide new insights into the formation of different tissues. In summary, this study provides a valuable resource for research on cereal inflorescence development, offering new clues for yield improvement.

Molecular classification and intratumoral heterogeneity of gastric adenocarcinoma

Gastric cancers frequently harbor striking histological complexity and diversity between lesions as well as within single lesions, known as inter- and intratumoral heterogeneity, respectively. The latest World Health Organization Classification of Tumors designated more than 30 histological subtypes for gastric epithelial tumors, assigning 12 subtypes for gastric adenocarcinoma (GAD). Meanwhile, recent advances in genome-wide analyses have provided molecular aspects to the histological classification of GAD, and consequently revealed different molecular traits underlying these histological subtypes. Moreover, accumulating knowledge of comprehensive molecular profiles has led to establishing molecular classifications of GAD, which are often associated with clinical biomarkers for therapeutics and prognosis. However, most of our knowledge of GAD molecular profiles is based on inter-tumoral heterogeneity, and the molecular profiles underlying intratumoral heterogeneity are yet to be determined. In this review, recently established molecular classifications of GAD are introduced in the aspect of pathological diagnosis and are discussed in the context of intratumoral heterogeneity.

Powerful and accurate detection of temporal gene expression patterns from multi-sample multi-stage single-cell transcriptomics data with TDEseq

We present a non-parametric statistical method called TDEseq that takes full advantage of smoothing splines basis functions to account for the dependence of multiple time points in scRNA-seq studies, and uses hierarchical structure linear additive mixed models to model the correlated cells within an individual. As a result, TDEseq demonstrates powerful performance in identifying four potential temporal expression patterns within a specific cell type. Extensive simulation studies and the analysis of four published scRNA-seq datasets show that TDEseq can produce well-calibrated p-values and up to 20% power gain over the existing methods for detecting temporal gene expression patterns.

Compound Activity Prediction with Dose-Dependent Transcriptomic Profiles and Deep Learning

Predicting compound activity in assays is a long-standing challenge in drug discovery. Computational models based on compound-induced gene expression signatures from a single profiling assay have shown promise toward predicting compound activity in other, seemingly unrelated, assays. Applications of such models include predicting mechanisms-of-action (MoA) for phenotypic hits, identifying off-target activities, and identifying polypharmacologies. Here, we introduce transcriptomics-to-activity transformer (TAT) models that leverage gene expression profiles observed over compound treatment at multiple concentrations to predict the compound activity in other biochemical or cellular assays. We built TAT models based on gene expression data from a RASL-seq assay to predict the activity of 2692 compounds in 262 dose-response assays. We obtained useful models for 51% of the assays, as determined through a realistic held-out set. Prospectively, we experimentally validated the activity predictions of a TAT model in a malaria inhibition assay. With a 63% hit rate, TAT successfully identified several submicromolar malaria inhibitors. Our results thus demonstrate the potential of transcriptomic responses over compound concentration and the TAT modeling framework as a cost-efficient way to identify the bioactivities of promising compounds across many assays.

Raman Flow Cytometry and Its Biomedical Applications

Raman flow cytometry (RFC) uniquely integrates the "label-free" capability of Raman spectroscopy with the "high-throughput" attribute of traditional flow cytometry (FCM), offering exceptional performance in cell characterization and sorting. Unlike conventional FCM, RFC stands out for its elimination of the dependency on fluorescent labels, thereby reducing interference with the natural state of cells. Furthermore, it significantly enhances the detection information, providing a more comprehensive chemical fingerprint of cells. This review thoroughly discusses the fundamental principles and technological advantages of RFC and elaborates on its various applications in the biomedical field, from identifying and characterizing cancer cells for in vivo cancer detection and surveillance to sorting stem cells, paving the way for cell therapy, and identifying metabolic products of microbial cells, enabling the differentiation of microbial subgroups. Moreover, we delve into the current challenges and future directions regarding the improvement in sensitivity and throughput. This holds significant implications for the field of cell analysis, especially for the advancement of metabolomics.

Antigen presenting cells in cancer immunity and mediation of immune checkpoint blockade

Antigen-presenting cells (APCs) are pivotal mediators of immune responses. Their role has increasingly been spotlighted in the realm of cancer immunology, particularly as our understanding of immunotherapy continues to evolve and improve. There is growing evidence that these cells play a non-trivial role in cancer immunity and have roles dependent on surface markers, growth factors, transcription factors, and their surrounding environment. The main dendritic cell (DC) subsets found in cancer are conventional DCs (cDC1 and cDC2), monocyte-derived DCs (moDC), plasmacytoid DCs (pDC), and mature and regulatory DCs (mregDC). The notable subsets of monocytes and macrophages include classical and non-classical monocytes, macrophages, which demonstrate a continuum from a pro-inflammatory (M1) phenotype to an anti-inflammatory (M2) phenotype, and tumor-associated macrophages (TAMs). Despite their classification in the same cell type, each subset may take on an immune-activating or immunosuppressive phenotype, shaped by factors in the tumor microenvironment (TME). In this review, we introduce the role of DCs, monocytes, and macrophages and recent studies investigating them in the cancer immunity context. Additionally, we review how certain characteristics such as abundance, surface markers, and indirect or direct signaling pathways of DCs and macrophages may influence tumor response to immune checkpoint blockade (ICB) therapy. We also highlight existing knowledge gaps regarding the precise contributions of different myeloid cell subsets in influencing the response to ICB therapy. These findings provide a summary of our current understanding of myeloid cells in mediating cancer immunity and ICB and offer insight into alternative or combination therapies that may enhance the success of ICB in cancers.

Advanced optical imaging for the rational design of nanomedicines

Despite the enormous potential of nanomedicines to shape the future of medicine, their clinical translation remains suboptimal. Translational challenges are present in every step of the development pipeline, from a lack of understanding of patient heterogeneity to insufficient insights on nanoparticle properties and their impact on material-cell interactions. Here, we discuss how the adoption of advanced optical microscopy techniques, such as super-resolution optical microscopies, correlative techniques, and high-content modalities, could aid the rational design of nanocarriers, by characterizing the cell, the nanomaterial, and their interaction with unprecedented spatial and/or temporal detail. In this nanomedicine arena, we will discuss how the implementation of these techniques, with their versatility and specificity, can yield high volumes of multi-parametric data; and how machine learning can aid the rapid advances in microscopy: from image acquisition to data interpretation.

The many paths to artemisinin resistance in Plasmodium falciparum

Emerging resistance against artemisinin (ART) poses a major challenge in controlling malaria. Parasites with mutations in PfKelch13, the major marker for ART resistance, are known to reduce hemoglobin endocytosis, induce unfolded protein response (UPR), elevate phosphatidylinositol-3-phosphate (PI3P) levels, and stimulate autophagy. Nonetheless, PfKelch13-independent resistance is also reported, indicating extensive complementation by reconfiguration in the parasite metabolome and transcriptome. These findings implicate that there may not be a single 'universal identifier' of ART resistance. This review sheds light on the molecular, transcriptional, and metabolic pathways associated with ART resistance, while also highlighting the interplay between cellular heterogeneity, environmental stress, and ART sensitivity.

Estimating genome-wide DNA methylation heterogeneity with methylation patterns

BACKGROUND

In a heterogeneous population of cells, individual cells can behave differently and respond variably to the environment. This cellular diversity can be assessed by measuring DNA methylation patterns. The loci with variable methylation patterns are informative of cellular heterogeneity and may serve as biomarkers of diseases and developmental progression. Cell-to-cell methylation heterogeneity can be evaluated through single-cell methylomes or computational techniques for pooled cells. However, the feasibility and performance of these approaches to precisely estimate methylation heterogeneity require further assessment.

RESULTS

Here, we proposed model-based methods adopted from a mathematical framework originally from biodiversity, to estimate genome-wide DNA methylation heterogeneity. We evaluated the performance of our models and the existing methods with feature comparison, and tested on both synthetic datasets and real data. Overall, our methods have demonstrated advantages over others because of their better correlation with the actual heterogeneity. We also demonstrated that methylation heterogeneity offers an additional layer of biological information distinct from the conventional methylation level. In the case studies, we showed that distinct profiles of methylation heterogeneity in CG and non-CG methylation can predict the regulatory roles between genomic elements in Arabidopsis. This opens up a new direction for plant epigenomics. Finally, we demonstrated that our score might be able to identify loci in human cancer samples as putative biomarkers for early cancer detection.

CONCLUSIONS

We adopted the mathematical framework from biodiversity into three model-based methods for analyzing genome-wide DNA methylation heterogeneity to monitor cellular heterogeneity. Our methods, namely MeH, have been implemented, evaluated with existing methods, and are open to the research community.

Quantitative biochemical phenotypic heterogeneity of senescent macrophage at a single cell level by Synchrotron Radiation Fourier Transform Infrared Microspectroscopy

Macrophage senescence plays an important role in pathophysiological process of age-related diseases such as atherosclerosis, chronic obstructive pulmonary disease (COPD), pulmonary fibrosis, and lung cancer. After macrophage senescence, the biochemical phenotypes related to biological functions showed great heterogeneity. However, the biochemical phenotype and phenotypic heterogeneity of senescent macrophage has not been fully understood. Exploring the phenotype of biochemical substances in senescent macrophage will be helpful for understanding the function of senescent macrophage and finding out the potential mechanism between immune macrophage senescence and age-related diseases. In this study, we employed SR-FTIR microspectroscopy to detect the biochemical phenotype and phenotypic heterogeneity of single macrophage. The whole infrared spectra of senescent macrophages shifted, indicating biochemical substance changes within senescent macrophages. PCA and intercellular Euclidean distance statistical analysis based on specific spectra regions revealed dynamic changes of lipids and proteins during macrophage senescence. This proved that SR-FTIR microspectroscopy is an effective tool to detect the single cell biochemical phenotype transformation and phenotypic heterogeneity during macrophage senescence. It is of great significance to provide an evaluation method or clue for the study of cellular functions related to intracellular biochemical substances.

Advances in novel strategies for isolation, characterization, and analysis of CTCs and ctDNA

Over the past decade, the detection and analysis of liquid biopsy biomarkers such as circulating tumor cells (CTCs) and circulating tumor DNA (ctDNA) have advanced significantly. They have received recognition for their clinical usefulness in detecting cancer at an early stage, monitoring disease, and evaluating treatment response. The emergence of liquid biopsy has been a helpful development, as it offers a minimally invasive, rapid, real-time monitoring, and possible alternative to traditional tissue biopsies. In resource-limited settings, the ideal platform for liquid biopsy should not only extract more CTCs or ctDNA from a minimal sample volume but also accurately represent the molecular heterogeneity of the patient's disease. This review covers novel strategies and advancements in CTC and ctDNA-based liquid biopsy platforms, including microfluidic applications and comprehensive analysis of molecular complexity. We discuss these systems' operational principles and performance efficiencies, as well as future opportunities and challenges for their implementation in clinical settings. In addition, we emphasize the importance of integrated platforms that incorporate machine learning and artificial intelligence in accurate liquid biopsy detection systems, which can greatly improve cancer management and enable precision diagnostics.

Quantitative biochemical phenotypic heterogeneity of macrophages after myelin debris phagocytosis at a single cell level by synchrotron radiation fourier transform infrared microspectroscopy

During the immuno-inflammatory pathophysiological process of spinal cord injury, traumatic brain injury, and ischemic stroke, macrophages play an important role in phagocytizing and clearing degenerated myelin debris. After phagocytizing myelin debris, the biochemical phenotypes related to the biological function of macrophages show vast heterogeneity; however, it is not fully understood. Detecting biochemical changes after myelin debris phagocytosis by macrophages at a single-cell level is helpful to characterize phenotypic and functional heterogeneity. In this study, based on the cell model of myelin debris phagocytosis by macrophages in vitro, the biochemical changes in macrophages were investigated using Synchrotron radiation-based Fourier transform infrared (SR-FTIR) microspectroscopy. Infrared spectrum fluctuations, principal component analysis, and cell-to-cell Euclidean distance statistical analysis of specific spectrum regions revealed dynamic and significant changes in proteins and lipids within macrophages after myelin debris phagocytosis. Thus, SR-FTIR microspectroscopy is a powerful identification toolkit for exploring biochemical phenotype heterogeneity transformation that may be of great importance to providing an evaluation strategy for studying cell functions related to cellular substance distribution and metabolism.

A posterior probability based Bayesian method for single-cell RNA-seq data imputation

Single-cell RNA-sequencing (scRNA-seq) data suffer from a lot of zeros. Such dropout events impede the downstream data analyses. We propose BayesImpute to infer and impute dropouts from the scRNA-seq data. Using the expression rate and coefficient of variation of the genes within the cell subpopulation, BayesImpute first determines likely dropouts, and then constructs the posterior distribution for each gene and uses the posterior mean to impute dropout values. Some simulated and real experiments show that BayesImpute can effectively identify dropout events and reduce the introduction of false positive signals. Additionally, BayesImpute successfully recovers the true expression levels of missing values, restores the gene-to-gene and cell-to-cell correlation coefficient, and maintains the biological information in bulk RNA-seq data. Furthermore, BayesImpute boosts the clustering and visualization of cell subpopulations and improves the identification of differentially expressed genes. We further demonstrate that, in comparison to other statistical-based imputation methods, BayesImpute is scalable and fast with minimal memory usage.

Regulatory role of Chitinase 3-like 1 gene in papillary thyroid carcinoma proved by integration analyses of single-cell sequencing with cohort and experimental validations

Papillary thyroid carcinoma (PTC) is one of the most common thyroid carcinomas. The gross extrathyroidal extension and extensive metastases of PTC lead to high rates of recurrence and poor clinical outcomes. However, the mechanisms underlying PTC development are poorly understood. In this study, using single-cell RNA sequencing, the transcriptome profiles of two PTC patients were addressed, including PTC1 with low malignancy and good prognosis and PTC2 with high malignancy and poor prognosis. We found that epithelial subcluster Epi02 was the most associated with the malignant development of PTC cells, with which the fold change of Chitinase 3-like 1 (CHI3L1) is on the top of the differentially expressed genes between PTC1 and PTC2 (P < 0.001). However CHI3L1 is rarely investigated in PTC as far. We then studied its role in PTC with a series of experiments. Firstly, qRT-PCR analysis of 14 PTC patients showed that the expression of CHI3L1 was positively correlated with malignancy. In addition, overexpression or silencing of CHI3L1 in TPC-1 cells, a PTC cell line, cultured in vitro showed that the proliferation, invasion, and metastasis of the cells were promoted or alleviated by CHI3L1. Further, immunohistochemistry analysis of 110 PTC cases revealed a significant relationship between CHI3L1 protein expression and PTC progression, especially the T (P < 0.001), N (P < 0.001), M stages (P = 0.007) and gross ETE (P < 0.001). Together, our results prove that CHI3L1 is a positive regulator of malignant development of PTC, and it promotes proliferation, invasion, and metastasis of PTC cells. Our study improves understanding of the molecular mechanisms underlying the progression of PTC and provides new insights for the clinical diagnosis and treatment of PTC.

There's more to it: uncovering genomewide DNA methylation heterogeneity

Tweetable abstract Monitoring changes in methylation heterogeneity can be powerful in detecting disease progression early. This editorial highlights the importance of profiling methylation heterogeneity and identifies existing measures and research gaps.

Single nuclei transcriptomics delineates complex immune and kidney cell interactions contributing to kidney allograft fibrosis

Chronic allograft dysfunction (CAD), characterized histologically by interstitial fibrosis and tubular atrophy, is the major cause of kidney allograft loss. Here, using single nuclei RNA sequencing and transcriptome analysis, we identified the origin, functional heterogeneity, and regulation of fibrosis-forming cells in kidney allografts with CAD. A robust technique was used to isolate individual nuclei from kidney allograft biopsies and successfully profiled 23,980 nuclei from five kidney transplant recipients with CAD and 17,913 nuclei from three patients with normal allograft function. Our analysis revealed two distinct states of fibrosis in CAD; low and high extracellular matrix (ECM) with distinct kidney cell subclusters, immune cell types, and transcriptional profiles. Imaging mass cytometry analysis confirmed increased ECM deposition at the protein level. Proximal tubular cells transitioned to an injured mixed tubular (MT1) phenotype comprised of activated fibroblasts and myofibroblast markers, generated provisional ECM which recruited inflammatory cells, and served as the main driver of fibrosis. MT1 cells in the high ECM state achieved replicative repair evidenced by dedifferentiation and nephrogenic transcriptional signatures. MT1 in the low ECM state showed decreased apoptosis, decreased cycling tubular cells, and severe metabolic dysfunction, limiting the potential for repair. Activated B, T and plasma cells were increased in the high ECM state, while macrophage subtypes were increased in the low ECM state. Intercellular communication between kidney parenchymal cells and donor-derived macrophages, detected several years post-transplantation, played a key role in injury propagation. Thus, our study identified novel molecular targets for interventions aimed to ameliorate or prevent allograft fibrogenesis in kidney transplant recipients.

The CellPhe toolkit for cell phenotyping using time-lapse imaging and pattern recognition

With phenotypic heterogeneity in whole cell populations widely recognised, the demand for quantitative and temporal analysis approaches to characterise single cell morphology and dynamics has increased. We present CellPhe, a pattern recognition toolkit for the unbiased characterisation of cellular phenotypes within time-lapse videos. CellPhe imports tracking information from multiple segmentation and tracking algorithms to provide automated cell phenotyping from different imaging modalities, including fluorescence. To maximise data quality for downstream analysis, our toolkit includes automated recognition and removal of erroneous cell boundaries induced by inaccurate tracking and segmentation. We provide an extensive list of features extracted from individual cell time series, with custom feature selection to identify variables that provide greatest discrimination for the analysis in question. Using ensemble classification for accurate prediction of cellular phenotype and clustering algorithms for the characterisation of heterogeneous subsets, we validate and prove adaptability using different cell types and experimental conditions.

Deciphering endothelial heterogeneity in health and disease at single-cell resolution: progress and perspectives

Endothelial cells (ECs) constitute the inner lining of vascular beds in mammals and are crucial for homeostatic regulation of blood vessel physiology, but also play a key role in pathogenesis of many diseases, thereby representing realistic therapeutic targets. However, it has become evident that ECs are heterogeneous, encompassing several subtypes with distinct functions, which makes EC targeting and modulation in diseases challenging. The rise of the new single-cell era has led to an emergence of studies aimed at interrogating transcriptome diversity along the vascular tree, and has revolutionized our understanding of EC heterogeneity from both a physiological and pathophysiological context. Here, we discuss recent landmark studies aimed at teasing apart the heterogeneous nature of ECs. We cover driving (epi)genetic, transcriptomic, and metabolic forces underlying EC heterogeneity in health and disease, as well as current strategies used to combat disease-enriched EC phenotypes, and propose strategies to transcend largely descriptive heterogeneity towards prioritization and functional validation of therapeutically targetable drivers of EC diversity. Lastly, we provide an overview of the most recent advances and hurdles in single EC OMICs.

Dissecting insect cell heterogeneity during influenza VLP production using single-cell transcriptomics

The insect cell-baculovirus expression vector system (IC-BEVS) has been widely used to produce recombinant protein at high titers, including complex virus-like particles (VPLs). However, cell-to-cell variability upon infection is yet one of the least understood phenomena in virology, and little is known about its impact on production of therapeutic proteins. This study aimed at dissecting insect cell population heterogeneity during production of influenza VLPs in IC-BEVS using single-cell RNA-seq (scRNA-seq). High Five cell population was shown to be heterogeneous even before infection, with cell cycle being one of the factors contributing for this variation. In addition, infected insect cells were clustered according to the timing and level of baculovirus genes expression, with each cluster reporting similar influenza VLPs transgenes (i.e., hemagglutinin and M1) transcript counts. Trajectory analysis enabled to track infection progression throughout pseudotime. Specific pathways such as translation machinery, protein folding, sorting and degradation, endocytosis and energy metabolism were identified as being those which vary the most during insect cell infection and production of Influenza VLPs. Overall, this study lays the ground for the application of scRNA-seq in IC-BEVS processes to isolate relevant biological mechanisms during recombinant protein expression towards its further optimization.

Photothermal plasmonic microballs for non-contact single-cell calcium ionophore delivery in heterogeneous cells

Hierarchical plasmonic nanostructures comprising gold nanorod (AuNR)-covered microballs via syringe-injection reduction show good potential for selective single-cell calcium ionophore (A23187) delivery and apoptosis induction in heterogenous cancer cells.

Monocyte/Macrophage Heterogeneity during Skin Wound Healing in Mice

Monocytes (Mos)/macrophages (Mϕs) orchestrate biological processes critical for efficient skin wound healing. However, current understanding of skin wound Mo/Mϕ heterogeneity is limited by traditional experimental approaches such as flow cytometry and immunohistochemistry. Therefore, we sought to more fully explore Mo/Mϕ heterogeneity and associated state transitions during the course of excisional skin wound healing in mice using single-cell RNA sequencing. The live CD45+CD11b+Ly6G- cells were isolated from skin wounds of C57BL/6 mice on days 3, 6, and 10 postinjury and captured using the 10x Genomics Chromium platform. A total of 2813 high-quality cells were embedded into a uniform manifold approximation and projection space, and eight clusters of distinctive cell populations were identified. Cluster dissimilarity and differentially expressed gene analysis categorized those clusters into three groups: early-stage/proinflammatory, late-stage/prohealing, and Ag-presenting phenotypes. Signature gene and Gene Ontology analysis of each cluster provided clues about the different functions of the Mo/Mϕ subsets, including inflammation, chemotaxis, biosynthesis, angiogenesis, proliferation, and cell death. Quantitative PCR assays validated characteristics of early- versus late-stage Mos/Mϕs inferred from our single-cell RNA sequencing dataset. Additionally, cell trajectory analysis by pseudotime and RNA velocity and adoptive transfer experiments indicated state transitions between early- and late-state Mos/Mϕs as healing progressed. Finally, we show that the chemokine Ccl7, which was a signature gene for early-stage Mos/Mϕs, preferentially induced the accumulation of proinflammatory Ly6C+F4/80lo/- Mos/Mϕs in mouse skin wounds. In summary, our data demonstrate the complexity of Mo/Mϕ phenotypes, their dynamic behavior, and diverse functions during normal skin wound healing.

scATACpipe: A nextflow pipeline for comprehensive and reproducible analyses of single cell ATAC-seq data

Single cell ATAC-seq (scATAC-seq) has become the most widely used method for profiling open chromatin landscape of heterogeneous cell populations at a single-cell resolution. Although numerous software tools and pipelines have been developed, an easy-to-use, scalable, reproducible, and comprehensive pipeline for scATAC-seq data analyses is still lacking. To fill this gap, we developed scATACpipe, a Nextflow pipeline, for performing comprehensive analyses of scATAC-seq data including extensive quality assessment, preprocessing, dimension reduction, clustering, peak calling, differential accessibility inference, integration with scRNA-seq data, transcription factor activity and footprinting analysis, co-accessibility inference, and cell trajectory prediction. scATACpipe enables users to perform the end-to-end analysis of scATAC-seq data with three sub-workflow options for preprocessing that leverage 10x Genomics Cell Ranger ATAC software, the ultra-fast Chromap procedures, and a set of custom scripts implementing current best practices for scATAC-seq data preprocessing. The pipeline extends the R package ArchR for downstream analysis with added support to any eukaryotic species with an annotated reference genome. Importantly, scATACpipe generates an all-in-one HTML report for the entire analysis and outputs cluster-specific BAM, BED, and BigWig files for visualization in a genome browser. scATACpipe eliminates the need for users to chain different tools together and facilitates reproducible and comprehensive analyses of scATAC-seq data from raw reads to various biological insights with minimal changes of configuration settings for different computing environments or species. By applying it to public datasets, we illustrated the utility, flexibility, versatility, and reliability of our pipeline, and demonstrated that our scATACpipe outperforms other workflows.

ATAXIC: An Algorithm to Quantify Transcriptomic Perturbation Heterogeneity in Single Cancer Cells

The single-cell RNA sequencing (scRNA-seq) has recently been widely utilized to quantify transcriptomic profiles in single cells of bulk tumors. The transcriptomic profiles in single cells facilitate the investigation of intratumor heterogeneity that is unlikely confounded by the nontumor components. We proposed an algorithm (ATAXIC) to quantify the heterogeneity of transcriptomic perturbations (TPs) in single cancer cells. ATAXIC calculated the TP heterogeneity level of a single cell based on the standard deviations of the absolute z-scored gene expression values for tens of thousands of genes, reflecting the asynchronous degree of transcriptomic alterations relative to the central (mean) tendency. By analyzing scRNA-seq datasets for eight cancer types, we revealed that ATAXIC scores were likely to correlate positively with the enrichment scores of various proliferation and oncogenic signatures, DNA damage repair, treatment resistance, and unfavorable phenotypes and outcomes in cancer. The ATAXIC scores varied among different cancer types, with lung cancer and melanoma having the lowest average scores and clear cell renal cell carcinoma having the highest average scores. The low TP heterogeneity in lung cancer and melanoma could bestow relatively higher response rates to immune checkpoint inhibitors on both cancer types. In conclusion, ATAXIC is a useful algorithm to quantify the TP heterogeneity in single cancer cells, as well as providing new insights into tumor biology.

Cellular heterogeneity in disease severity and clinical outcome: Granular understanding of immune response is key

During an infectious disease progression, it is crucial to understand the cellular heterogeneity underlying the differential immune response landscape that will augment the precise information of the disease severity modulators, leading to differential clinical outcome. Patients with COVID-19 display a complex yet regulated immune profile with a heterogeneous array of clinical manifestation that delineates disease severity sub-phenotypes and worst clinical outcomes. Therefore, it is necessary to elucidate/understand/enumerate the role of cellular heterogeneity during COVID-19 disease to understand the underlying immunological mechanisms regulating the disease severity. This article aims to comprehend the current findings regarding dysregulation and impairment of immune response in COVID-19 disease severity sub-phenotypes and relate them to a wide array of heterogeneous populations of immune cells. On the basis of the findings, it suggests a possible functional correlation between cellular heterogeneity and the COVID-19 disease severity. It highlights the plausible modulators of age, gender, comorbidities, and hosts' genetics that may be considered relevant in regulating the host response and subsequently the COVID-19 disease severity. Finally, it aims to highlight challenges in COVID-19 disease that can be achieved by the application of single-cell genomics, which may aid in delineating the heterogeneity with more granular understanding. This will augment our future pandemic preparedness with possibility to identify the subset of patients with increased diseased severity.

Immunosenescence, aging and successful aging

Aging induces a series of immune related changes, which is called immunosenescence, playing important roles in many age-related diseases, especially neurodegenerative diseases, tumors, cardiovascular diseases, autoimmune diseases and coronavirus disease 2019(COVID-19). However, the mechanism of immunosenescence, the association with aging and successful aging, and the effects on diseases are not revealed obviously. In order to provide theoretical basis for preventing or controlling diseases effectively and achieve successful aging, we conducted the review and found that changes of aging-related phenotypes, deterioration of immune organ function and alterations of immune cell subsets participated in the process of immunosenescence, which had great effects on the occurrence and development of age-related diseases.

Associating pathways with diseases using single-cell expression profiles and making inferences about potential drugs

Finding direct dependencies between genetic pathways and diseases has been the target of multiple studies as it has many applications. However, due to cellular heterogeneity and limitations of the number of samples for bulk expression profiles, such studies have faced hurdles in the past. Here, we propose a method to perform single-cell expression-based inference of association between pathway, disease and cell-type (sci-PDC), which can help to understand their cause and effect and guide precision therapy. Our approach highlighted reliable relationships between a few diseases and pathways. Using the example of diabetes, we have demonstrated how sci-PDC helps in tracking variation of association between pathways and diseases with changes in age and species. The variation in pathways-disease associations in mice and humans revealed critical facts about the suitability of the mouse model for a few pathways in the context of diabetes. The coherence between results from our method and previous reports, including information about the drug target pathways, highlights its reliability for multidimensional utility.

Modeling Susceptibility to Cardiotoxicity in Cancer Therapy Using Human iPSC-Derived Cardiac Cells and Systems Biology

The development of human-induced pluripotent stem cell-derived cardiac cell types has created a new paradigm in assessing drug-induced cardiotoxicity. Advances in genomics and epigenomics have also implicated several genomic loci and biological pathways that may contribute to susceptibility to cancer therapies. In this review, we first provide a brief overview of the cardiotoxicity associated with chemotherapy. We then provide a detailed summary of systems biology approaches being applied to elucidate potential molecular mechanisms involved in cardiotoxicity. Finally, we discuss combining systems biology approaches with iPSC technology to help discover molecular mechanisms associated with cardiotoxicity.

Timelapse viability assay to detect division and death of primary multiple myeloma cells in response to drug treatments with single cell resolution

Heterogeneity among cancer cells and in the tumor microenvironment (TME) is thought to be a significant contributor to the heterogeneity of clinical therapy response observed between patients and can evolve over time. A primary example of this is multiple myeloma (MM), a generally incurable cancer where such heterogeneity contributes to the persistent evolution of drug resistance. However, there is a paucity of functional assays for studying this heterogeneity in patient samples or for assessing the influence of the patient TME on therapy response. Indeed, the population-averaged data provided by traditional drug response assays and the large number of cells required for screening remain significant hurdles to advancement. To address these hurdles, we developed a suite of accessible technologies for quantifying functional drug response to a panel of therapies in ex vivo three-dimensional culture using small quantities of a patient's own cancer and TME components. This suite includes tools for label-free single-cell identification and quantification of both cell division and death events with a standard brightfield microscope, an open-source software package for objective image analysis and feasible data management of multi-day timelapse experiments, and a new approach to fluorescent detection of cell death that is compatible with long-term imaging of primary cells. These new tools and capabilities are used to enable sensitive, objective, functional characterization of primary MM cell therapy response in the presence of TME components, laying the foundation for future studies and efforts to enable predictive assessment drug efficacy for individual patients.

Pervasive sequence-level variation in the transcriptome of Plasmodium falciparum

Single-nucleotide variations (SNVs) in RNA, arising from co- and post-transcriptional phenomena including transcription errors and RNA-editing, are well studied in a range of organisms. In the malaria parasite Plasmodium falciparum, stage-specific and non-specific gene-expression variations accompany the parasite's array of developmental and morphological phenotypes over the course of its complex life cycle. However, the extent, rate and effect of sequence-level variation in the parasite's transcriptome are unknown. Here, we report the presence of pervasive, non-specific SNVs in the P. falciparum transcriptome. SNV rates for a gene were correlated to gene length (r[Formula: see text]0.65-0.7) but not to the AT-content of that gene. Global SNV rates for the P. falciparum lines we used, and for publicly available P. vivax and P. falciparum clinical isolate datasets, were of the order of 10-3 per base, ∼10× higher than rates we calculated for bacterial datasets. These variations may reflect an intrinsic transcriptional error rate in the parasite, and RNA editing may be responsible for a subset of them. This seemingly characteristic property of the parasite may have implications for clinical outcomes and the basic biology and evolution of P. falciparum and parasite biology more broadly. We anticipate that our study will prompt further investigations into the exact sources, consequences and possible adaptive roles of these SNVs.

Integrated analysis of an in vivo model of intra-nasal exposure to instilled air pollutants reveals cell-type specific responses in the placenta

The placenta is a heterogeneous organ whose development involves complex interactions of trophoblasts with decidual, vascular, and immune cells at the fetal-maternal interface. It maintains a critical balance between maternal and fetal homeostasis. Placental dysfunction can lead to adverse pregnancy outcomes including intra-uterine growth restriction, pre-eclampsia, or pre-term birth. Exposure to environmental pollutants contributes to the development of placental abnormalities, with poorly understood molecular underpinning. Here we used a mouse (C57BL/6) model of environmental pollutant exposure by administration of a particulate matter (SRM1649b at 300 μg/day/mouse) suspension intra-nasally beginning 2 months before conception and during gestation, in comparison to saline-exposed controls. Placental transcriptomes, at day 19 of gestation, were determined using bulk RNA-seq from whole placentas of exposed (n = 4) and control (n = 4) animals and scRNAseq of three distinct placental layers, followed by flow cytometry analysis of the placental immune cell landscape. Our results indicate a reduction in vascular placental cells, especially cells responsible for structural integrity, and increase in trophoblast proliferation in animals exposed to particulate matter. Pollution-induced inflammation was also evident, especially in the decidual layer. These data indicate that environmental exposure to air pollutants triggers changes in the placental cellular composition, mediating adverse pregnancy outcomes.

The role of serum amyloid A1 in the adipogenic differentiation of human adipose-derived stem cells basing on single-cell RNA sequencing analysis

Background

Adipose-derived stem cells (ASCs) are obtained from a variety of sources in vivo where they present in large quantities. These cells are suitable for use in autologous transplantation and the construction of tissue-engineered adipose tissue. Studies have shown that ASCs differentiation is in a high degree of heterogeneity, yet the molecular basis including key regulators of differentiation remains to clarify.

Methods

We performed single-cell RNA sequencing and bioinformatics analysis on both undifferentiated (ASC-GM group) and adipogenically differentiated human ASCs (ASC-AD group, ASCs were cultured in adipogenic inducing medium for 1 week). And then, we verified the results of serum amyloid A1 (SAA1) with western blotting, immunofluorescence staining, oil red O staining. After these experiments, we down-regulated the expression of serum amyloid A1 (SAA1) gene to verify the adipogenic differentiation ability of ASCs.

Results

In single-cell RNA sequence analyzing, we obtained 4415 cells in the ASC-GM group and 4634 cells in the ASC-AD group. The integrated sample cells could be divided into 11 subgroups (0–10 cluster). The cells in cluster 0, 2, 5 were came from ASC-GM group and the cells in cluster 1, 3, 7 came from ASC-AD group. The cells of cluster 4 and 6 came from both ASC-GM and ASC-AD groups. Fatty acid binding protein 4, fatty acid binding protein 5, complement factor D, fatty acid desaturase 1, and insulin like growth factor binding protein 5 were high expressed in category 1 and 7. Regulation of inflammatory response is the rank 1 biological processes. And cellular responses to external stimuli, negative regulation of defense response and acute inflammatory response are included in top 20 biological processes. Based on the MCODE results, we found that SAA1, C-C Motif Chemokine Ligand 5 (CCL5), and Annexin A1 (ANXA1) significantly highly expressed during adipogenic differentiation. Western blot and immunofluorescent staining results showed that SAA1 increased during adipogenesis. And the area of ORO positive staining in siSAA1 cells was significantly lower than in the siControl (negative control) cells.

Conclusions

Our results also indicated that our adipogenic induction was successful, and there was great heterogeneity in the adipogenic differentiation of ASCs. SAA1 with the regulation of inflammatory response were involved in adipogenesis of ASCs based on single-cell RNA sequencing analysis. The data obtained will help to elucidate the intrinsic mechanism of heterogeneity in the differentiation process of stem cells, thus, guiding the regulation of self-renewal and differentiation of adult stem cells.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13287-022-02873-5.

Deciphering Tumour Heterogeneity: From Tissue to Liquid Biopsy

Simple Summary

Most malignant tumours are highly heterogeneous at molecular and phenotypic levels. Tumour variability poses challenges for the management of patients, as it arises between patients and even evolves in space and time within a single patient. Currently, treatment-decision making usually relies on the molecular characteristics of a limited tumour tissue sample at the time of diagnosis or disease progression but does not take into account the complexity of the bulk tumours and their constant evolution over time. In this review, we explore the extent of tumour heterogeneity and report the mechanisms that promote and sustain this diversity in cancers. We summarise the clinical strikes of tumour diversity in the management of patients with cancer. Finally, we discuss the current material and technological approaches that are relevant to adequately appreciate tumour heterogeneity.

Abstract

Human solid malignancies harbour a heterogeneous set of cells with distinct genotypes and phenotypes. This heterogeneity is installed at multiple levels. A biological diversity is commonly observed between tumours from different patients (inter-tumour heterogeneity) and cannot be fully captured by the current consensus molecular classifications for specific cancers. To extend the complexity in cancer, there are substantial differences from cell to cell within an individual tumour (intra-tumour heterogeneity, ITH) and the features of cancer cells evolve in space and time. Currently, treatment-decision making usually relies on the molecular characteristics of a limited tumour tissue sample at the time of diagnosis or disease progression but does not take into account the complexity of the bulk tumours and their constant evolution over time. In this review, we explore the extent of tumour heterogeneity with an emphasis on ITH and report the mechanisms that promote and sustain this diversity in cancers. We summarise the clinical strikes of ITH in the management of patients with cancer. Finally, we discuss the current material and technological approaches that are relevant to adequately appreciate ITH.

EMT-mediated regulation of CXCL1/5 for resistance to anti-EGFR therapy in colorectal cancer

The emergence of RAS/RAF mutant clone is the main feature of EGFR inhibitor resistance in KRAS wild-type colon cancer. However, its molecular mechanism is thought to be multifactorial, mainly due to cellular heterogeneity. In order to better understand the resistance mechanism in a single clone level, we successfully isolated nine cells with cetuximab-resistant (CR) clonality from in vitro system. All CR cells harbored either KRAS or BRAF mutations. Characteristically, these cells showed a higher EMT (Epithelial to mesenchymal transition) signature, showing increased EMT markers such as SNAI2. Moreover, the expression level of CXCL1/5, a secreted protein, was significantly higher in CR cells compared to the parental cells. In these CR cells, CXCL1/5 expression was coordinately regulated by SNAI2/NFKB and transactivated EGFR through CXCR/MMPI/EGF axis via autocrine singling. We also observed that combined cetuximab/MEK inhibitor not only showed growth inhibition but also reduced the secreted amounts of CXCL1/5. We further found that serum CXCL1/5 level was positively correlated with the presence of RAS/RAF mutation in colon cancer patients during cetuximab therapy, suggesting its role as a biomarker. These data indicated that the application of serum CXCL1/5 could be a potential serologic biomarker for predicting resistance to EGFR therapy in colorectal cancer.

Single-Cell Sequencing: Current Applications in Precision Onco-Genomics and Cancer Therapeutics

Simple Summary

Single-cell sequencing technologies are growing, advancing, and supporting new opportunities to better understand cancer. A variety of technologies are available that analyze the human transcriptome, genome, epigenome, and proteome, enabling integrated datasets. As a result, these integrated datasets contribute to new mechanistic insights and areas with therapeutic potential. This review summarizes the various single-cell sequencing techniques and provides examples of recent high-impact findings from the utilization of these technologies. Additionally, the translational relevance of these technologies and their use in clinical trials is described, along with the future potential for novel findings using these innovative methods.

Abstract

Single-cell sequencing encompasses a variety of technologies that evaluate cells at the genomic, transcriptomic, epigenomic, and proteomic levels. Each of these levels can be split into additional techniques that enable specific and optimized sequencing for a specialized purpose. At the transcriptomic level, single-cell sequencing has been used to understand immune-malignant cell networks, as well as differences between primary versus metastatic tumors. At the genomic and epigenomic levels, single-cell sequencing technology has been used to study genetic mutations involved in tumor evolution or the reprogramming of regulatory elements present in metastasized disease, respectively. Lastly, at the proteomic level, single-cell sequencing has been used to identify biomarkers important for predicting patient prognosis, as well as biomarkers essential for evaluating optimal treatment strategies. Integrated databases and atlases, as a result of large sequencing experiments, provide a vast array of information that can be applied to various studies and accessed by researchers to further answer scientific questions. This review summarizes recent, high-impact literature covering these aspects, as well as single-cell sequencing in the translational setting. Specifically, we review the potential that single-cell sequencing has in the clinic and its implementation in current clinical studies.

Microfluidics for single cell analysis

Cells have several internal molecules that are present in low amounts and any fluctuation in its number drives a change in cell behavior. These molecules present inside the cells are continuously fluctuating, thus producing noises in the intrinsic environment and thereby directly affecting the cellular behavior. Single-cell analysis using microfluidics is an important tool for monitoring cell behavior by analyzing internal molecules. Several gene circuits have been designed for this purpose that are labeled with fluorescence encoding genes for monitoring cell dynamics and behavior. We discuss herewith designed and fabricated microfluidics devices that are used for trapping and tracking cells under controlled environmental conditions. This chapter highlights microfluidics chip for monitoring cells to promote their basic understanding.

Surfactin cyclic lipopeptides change the plasma membrane composition and lateral organization in mammalian cells

The specific structure and composition of the cell plasma membrane (PM) is crucial for many cellular processes and can be targeted by various substances with potential medical applications. In this context, biosurfactants (BS) constitute a promising group of natural compounds that possess several biological functions, including anticancer activity. Despite the efficiency of BS, their mode of action had never been elucidated before. Here, we demonstrate the influence of cyclic lipopeptide surfactin (SU) on the PM of CHO-K1 cells. Both FLIM and svFCS experiments show that even a low concentration of SU causes significant changes in the membrane fluidity and dynamic molecular organization. Further, we demonstrate that SU causes a relevant dose-dependent reduction of cellular cholesterol by extracting it from the PM. Finally, we show that CHO-25RA cells characterized by increased cholesterol levels are more sensitive to SU treatment than CHO-K1 cells. We propose that sterols organizing the PM raft nanodomains, constitute a potential target for SU and other biosurfactants. In our opinion, the anticancer activity of biosurfactants is directly related with the higher cholesterol content found in many cancer cells.

Integrated Single-Cell Bioinformatics Analysis Reveals Intrinsic and Extrinsic Biological Characteristics of Hematopoietic Stem Cell Aging

Hematopoietic stem cell (HSC) aging, which is accompanied by loss of self-renewal capacity, myeloid-biased differentiation and increased risks of hematopoietic malignancies, is an important focus in stem cell research. However, the mechanisms underlying HSC aging have not been fully elucidated. In the present study, we integrated 3 independent single-cell transcriptome datasets of HSCs together and identified Stat3 and Ifngr1 as two markers of apoptosis-biased and inflammatory aged HSCs. Besides, common differentially expressed genes (DEGs) between young and aged HSCs were identified and further validated by quantitative RT-PCR. Functional enrichment analysis revealed that these DEGs were predominantly involved in the cell cycle and the tumor necrosis factor (TNF) signaling pathway. We further found that the Skp2-induced signaling pathway (Skp2→Cip1→CycA/CDK2→DP-1) contributed to a rapid transition through G1 phase in aged HSCs. In addition, analysis of the extrinsic alterations on HSC aging revealed the increased expression levels of inflammatory genes in bone marrow microenvironment. Colony formation unit assays showed that inflammatory cytokines promoted cellular senescence and that blockade of inflammatory pathway markedly rejuvenated aged HSC functions and increased B cell output. Collectively, our study elucidated the biological characteristics of HSC aging, and the genes and pathways we identified could be potential biomarkers and targets for the identification and rejuvenation of aged HSCs.

Single-Cell RNA Sequencing Reveals Cellular Heterogeneity and Stage Transition under Temperature Stress in Synchronized Plasmodium falciparum Cells

The malaria parasite has a complex life cycle exhibiting phenotypic and morphogenic variations in two different hosts by existing in heterogeneous developmental states. To investigate this cellular heterogeneity of the parasite within the human host, we performed single-cell RNA sequencing of synchronized Plasmodium cells under control and temperature treatment conditions. Using the Malaria Cell Atlas (https://www.sanger.ac.uk/science/tools/mca) as a guide, we identified 9 subtypes of the parasite distributed across known intraerythrocytic stages. Interestingly, temperature treatment results in the upregulation of the AP2-G gene, the master regulator of sexual development in a small subpopulation of the parasites. Moreover, we identified a heterogeneous stress-responsive subpopulation (clusters 5, 6, and 7 [∼10% of the total population]) that exhibits upregulation of stress response pathways under normal growth conditions. We also developed an online exploratory tool that will provide new insights into gene function under normal and temperature stress conditions. Thus, our study reveals important insights into cell-to-cell heterogeneity in the parasite population under temperature treatment that will be instrumental toward a mechanistic understanding of cellular adaptation and population dynamics in Plasmodium falciparum. IMPORTANCE The malaria parasite has a complex life cycle exhibiting phenotypic variations in two different hosts accompanied by cell-to-cell variability that is important for stress tolerance, immune evasion, and drug resistance. To investigate cellular heterogeneity determined by gene expression, we performed single-cell RNA sequencing (scRNA-seq) of about 12,000 synchronized Plasmodium cells under physiologically relevant normal (37°C) and temperature stress (40°C) conditions phenocopying the cyclic bouts of fever experienced during malarial infection. In this study, we found that parasites exhibit transcriptional heterogeneity in an otherwise morphologically synchronized culture. Also, a subset of parasites is continually committed to gametocytogenesis and stress-responsive pathways. These observations have important implications for understanding the mechanisms of drug resistance generation and vaccine development against the malaria parasite.

Single-cell transcriptomic sequencing analyses of cell heterogeneity during osteogenesis of human adipose-derived mesenchymal stem cells

Mesenchymal stem cells (MSCs) are known for their multilineage differentiation potential with immune-modulatory properties. The molecular underpinnings of differentiation remain largely undefined. In this study, we investigated the cellular and molecular features of chemically induced osteogenesis from MSC isolated from human adipose tissue (human adipose MSCs, hAMSCs) using single-cell RNA-sequencing (scRNA-seq). We found that a near complete differentiation of osteogenic clusters from hAMSCs under a directional induction. Both groups of cells are heterogeneous, and some of the hAMSCs cells are intrinsically prepared for osteogenesis, while variant OS clusters seems in cooperation with a due division of the general function. We identified a set of genes related to cell stress response highly expressed during the differentiation. We also characterized a series of transitional transcriptional waves throughout the process from hAMSCs to osteoblast and specified the unique gene networks and epigenetic status as key markers of osteogenesis.

A likelihood-based deconvolution of bulk gene expression data using single-cell references

Direct comparison of bulk gene expression profiles is complicated by distinct cell type mixtures in each sample which obscure whether observed differences are actually due to changes in expression levels themselves or simply due to differing cell type compositions. Single-cell technology has made it possible to measure gene expression in individual cells, achieving higher resolution at the expense of increased noise. If carefully incorporated, such single-cell data can be used to deconvolve bulk samples to yield accurate estimates of the true cell type proportions, thus enabling one to disentangle the effects of differential expression and cell type mixtures. Here, we propose a generative model and a likelihood-based inference method that uses asymptotic statistical theory and a novel optimization procedure to perform deconvolution of bulk RNA-seq data to produce accurate cell type proportion estimates. We demonstrate the effectiveness of our method, called RNA-Sieve, across a diverse array of scenarios involving real data and discuss extensions made uniquely possible by our probabilistic framework, including a demonstration of well-calibrated confidence intervals.

Label-free image-encoded microfluidic cell sorter with a scanning Bessel beam

The microfluidic-based, label-free image-guided cell sorter offers a low-cost, high information content, and disposable solution that overcomes many limitations in conventional cell sorters. However, flow confinement for most microfluidic devices is generally only one-dimensional using sheath flow. As a result, the equilibrium distribution of cells spreads beyond the focal plane of commonly used Gaussian laser excitation beams, resulting in a large number of blurred images that hinder subsequent cell sorting based on cell image features. To address this issue, we present a Bessel-Gaussian beam image-guided cell sorter with an ultra-long depth of focus, enabling focused images of >85% of passing cells. This system features label-free sorting capabilities based on features extracted from the output temporal waveform of a photomultiplier tube (PMT) detector. For the sorting of polystyrene beads, SKNO1 leukemia cells, and Scenedesmus green algae, our results indicate a sorting purity of 97%, 97%, and 98%, respectively, showing that the temporal waveforms from the PMT outputs have strong correlations with cell image features. These correlations are also confirmed by off-line reconstructed cell images from a temporal-spatial transformation algorithm tailored to the scanning Bessel-Gaussian beam.

Agent-based modeling reveals benefits of heterogeneous and stochastic cell populations during cGAS-mediated IFNβ production

MOTIVATION

The cGAS pathway is a component of the innate immune system responsible for the detection of pathogenic DNA and upregulation of interferon beta (IFNβ). Experimental evidence shows that IFNβ signaling occurs in highly heterogeneous cells and is stochastic in nature; however, the benefits of these attributes remain unclear. To investigate how stochasticity and heterogeneity affect IFNβ production, an agent-based model is developed to simulate both DNA transfection and viral infection.

RESULTS

We show that heterogeneity can enhance IFNβ responses during infection. Furthermore, by varying the degree of IFNβ stochasticity, we find that only a percentage of cells (20-30%) need to respond during infection. Going beyond this range provides no additional protection against cell death or reduction of viral load. Overall, these simulations suggest that heterogeneity and stochasticity are important for moderating immune potency while minimizing cell death during infection.

AVAILABILITY AND IMPLEMENTATION

Model repository is available at: https://github.com/ImmuSystems-Lab/AgentBasedModel-cGASPathway.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Machine learning for perturbational single-cell omics

Cell biology is fundamentally limited in its ability to collect complete data on cellular phenotypes and the wide range of responses to perturbation. Areas such as computer vision and speech recognition have addressed this problem of characterizing unseen or unlabeled conditions with the combined advances of big data, deep learning, and computing resources in the past 5 years. Similarly, recent advances in machine learning approaches enabled by single-cell data start to address prediction tasks in perturbation response modeling. We first define objectives in learning perturbation response in single-cell omics; survey existing approaches, resources, and datasets (https://github.com/theislab/sc-pert); and discuss how a perturbation atlas can enable deep learning models to construct an informative perturbation latent space. We then examine future avenues toward more powerful and explainable modeling using deep neural networks, which enable the integration of disparate information sources and an understanding of heterogeneous, complex, and unseen systems.

A comprehensive survey of regulatory network inference methods using single cell RNA sequencing data

Gene regulatory network is a complicated set of interactions between genetic materials, which dictates how cells develop in living organisms and react to their surrounding environment. Robust comprehension of these interactions would help explain how cells function as well as predict their reactions to external factors. This knowledge can benefit both developmental biology and clinical research such as drug development or epidemiology research. Recently, the rapid advance of single-cell sequencing technologies, which pushed the limit of transcriptomic profiling to the individual cell level, opens up an entirely new area for regulatory network research. To exploit this new abundant source of data and take advantage of data in single-cell resolution, a number of computational methods have been proposed to uncover the interactions hidden by the averaging process in standard bulk sequencing. In this article, we review 15 such network inference methods developed for single-cell data. We discuss their underlying assumptions, inference techniques, usability, and pros and cons. In an extensive analysis using simulation, we also assess the methods' performance, sensitivity to dropout and time complexity. The main objective of this survey is to assist not only life scientists in selecting suitable methods for their data and analysis purposes but also computational scientists in developing new methods by highlighting outstanding challenges in the field that remain to be addressed in the future development.

High-Throughput Raman Flow Cytometry and Beyond

Flow cytometry is a powerful tool with applications in diverse fields such as microbiology, immunology, virology, cancer biology, stem cell biology, and metabolic engineering. It rapidly counts and characterizes large heterogeneous populations of cells in suspension (e.g., blood cells, stem cells, cancer cells, and microorganisms) and dissociated solid tissues (e.g., lymph nodes, spleen, and solid tumors) with typical throughputs of 1,000-100,000 events per second (eps). By measuring cell size, cell granularity, and the expression of cell surface and intracellular molecules, it provides systematic insights into biological processes. Flow cytometers may also include cell sorting capabilities to enable subsequent additional analysis of the sorted sample (e.g., electron microscopy and DNA/RNA sequencing), cloning, and directed evolution. Unfortunately, traditional flow cytometry has several critical limitations as it mainly relies on fluorescent labeling for cellular phenotyping, which is an indirect measure of intracellular molecules and surface antigens. Furthermore, it often requires time-consuming preparation protocols and is incompatible with cell therapy. To overcome these difficulties, a different type of flow cytometry based on direct measurements of intracellular molecules by Raman spectroscopy, or "Raman flow cytometry" for short, has emerged. Raman flow cytometry obtains a chemical fingerprint of the cell in a nondestructive manner, allowing for single-cell metabolic phenotyping. However, its slow signal acquisition due to the weak light-molecule interaction of spontaneous Raman scattering prevents the throughput necessary to interrogate large cell populations in reasonable time frames, resulting in throughputs of about 1 eps. The remedy to this throughput limit lies in coherent Raman scattering methods such as stimulated Raman scattering (SRS) and coherent anti-Stokes Raman scattering (CARS), which offer a significantly enhanced light-sample interaction and hence enable high-throughput Raman flow cytometry, Raman imaging flow cytometry, and even Raman image-activated cell sorting (RIACS). In this Account, we outline recent advances, technical challenges, and emerging opportunities of coherent Raman flow cytometry. First, we review the principles of various types of SRS and CARS and introduce several techniques of coherent Raman flow cytometry such as CARS, multiplex CARS, Fourier-transform CARS, SRS, SRS imaging flow cytometry, and RIACS. Next, we discuss a unique set of applications enabled by coherent Raman flow cytometry, from microbiology and lipid biology to cancer detection and cell therapy. Finally, we describe future opportunities and challenges of coherent Raman flow cytometry including increasing sensitivity and throughput, integration with droplet microfluidics, utilizing machine learning techniques, or achieving in vivo flow cytometry. This Account summarizes the growing field of high-throughput Raman flow cytometry and the bright future it can bring.

Spatially Resolved Mass Spectrometry at the Single Cell: Recent Innovations in Proteomics and Metabolomics

Biological systems are composed of heterogeneous populations of cells that intercommunicate to form a functional living tissue. Biological function varies greatly across populations of cells, as each single cell has a unique transcriptome, proteome, and metabolome that translates to functional differences within single species and across kingdoms. Over the past decade, substantial advancements in our ability to characterize omic profiles on a single cell level have occurred, including in multiple spectroscopic and mass spectrometry (MS)-based techniques. Of these technologies, spatially resolved mass spectrometry approaches, including mass spectrometry imaging (MSI), have shown the most progress for single cell proteomics and metabolomics. For example, reporter-based methods using heavy metal tags have allowed for targeted MS investigation of the proteome at the subcellular level, and development of technologies such as laser ablation electrospray ionization mass spectrometry (LAESI-MS) now mean that dynamic metabolomics can be performed in situ. In this Perspective, we showcase advancements in single cell spatial metabolomics and proteomics over the past decade and highlight important aspects related to high-throughput screening, data analysis, and more which are vital to the success of achieving proteomic and metabolomic profiling at the single cell scale. Finally, using this broad literature summary, we provide a perspective on how the next decade may unfold in the area of single cell MS-based proteomics and metabolomics.