Native N-glycome profiling of single cells and ng-level blood isolates using label-free capillary electrophoresis-mass spectrometry.

The development of reliable single-cell dispensers and substantial sensitivity improvement in mass spectrometry made proteomic profiling of individual cells achievable. Yet, there are no established methods for single-cell glycome analysis due to the inability to amplify glycans and sample losses associated with sample processing and glycan labeling. In this work, we present an integrated platform coupling online in-capillary sample processing with high-sensitivity label-free capillary electrophoresis-mass spectrometry for N-glycan profiling of single mammalian cells. Direct and unbiased quantitative characterization of single-cell surface N-glycomes are demonstrated for HeLa and U87 cells, with the detection of up to 100 N-glycans per single cell. Interestingly, N-glycome alterations are unequivocally detected at the single-cell level in HeLa and U87 cells stimulated with lipopolysaccharide. The developed workflow is also applied to the profiling of ng-level amounts (5-500 ng) of blood-derived protein, extracellular vesicle, and total plasma isolates, resulting in over 170, 220, and 370 quantitated N-glycans, respectively.

Diagnostic leukapheresis reveals distinct phenotypes of NSCLC circulating tumor cells.

BACKGROUND: Circulating tumor cells (CTCs) hold immense promise for unraveling tumor heterogeneity and understanding treatment resistance. However, conventional methods, especially in cancers like non-small cell lung cancer (NSCLC), often yield low CTC numbers, hindering comprehensive analyses. This study addresses this limitation by employing diagnostic leukapheresis (DLA) to cancer patients, enabling the screening of larger blood volumes. To leverage DLA's full potential, this study introduces a novel approach for CTC enrichment from DLAs. METHODS: DLA was applied to six advanced stage NSCLC patients. For an unbiased CTC enrichment, a two-step approach based on negative depletion of hematopoietic cells was used. Single-cell (sc) whole-transcriptome sequencing was performed, and CTCs were identified based on gene signatures and inferred copy number variations. RESULTS: Remarkably, this innovative approach led to the identification of unprecedented 3,363 CTC transcriptomes. The extensive heterogeneity among CTCs was unveiled, highlighting distinct phenotypes related to the epithelial-mesenchymal transition (EMT) axis, stemness, immune responsiveness, and metabolism. Comparison with sc transcriptomes from primary NSCLC cells revealed that CTCs encapsulate the heterogeneity of their primary counterparts while maintaining unique CTC-specific phenotypes. CONCLUSIONS: In conclusion, this study pioneers a transformative method for enriching CTCs from DLA, resulting in a substantial increase in CTC numbers. This allowed the creation of the first-ever single-cell whole transcriptome in-depth characterization of the heterogeneity of over 3,300 NSCLC-CTCs. The findings not only confirm the diagnostic value of CTCs in monitoring tumor heterogeneity but also propose a CTC-specific signature that can be exploited for targeted CTC-directed therapies in the future. This comprehensive approach signifies a major leap forward, positioning CTCs as a key player in advancing our understanding of cancer dynamics and paving the way for tailored therapeutic interventions.

Comparative single-cell multiplex immunophenotyping of therapy-naive patients with rheumatoid arthritis, systemic sclerosis, and systemic lupus erythematosus shed light on disease-specific composition of the peripheral immune system.

Introduction: Systemic autoimmune diseases (SADs) are a significant burden on the healthcare system. Understanding the complexity of the peripheral immunophenotype in SADs may facilitate the differential diagnosis and identification of potential therapeutic targets. Methods: Single-cell mass cytometric immunophenotyping was performed on peripheral blood mononuclear cells (PBMCs) from healthy controls (HCs) and therapy-naive patients with rheumatoid arthritis (RA), progressive systemic sclerosis (SSc), and systemic lupus erythematosus (SLE). Immunophenotyping was performed on 15,387,165 CD45+ live single cells from 52 participants (13 cases/group), using an antibody panel to detect 34 markers. Results: Using the t-SNE (t-distributed stochastic neighbor embedding) algorithm, the following 17 main immune cell types were determined: CD4+/CD57- T cells, CD4+/CD57+ T cells, CD8+/CD161- T cells, CD8+/CD161+/CD28+ T cells, CD8dim T cells, CD3+/CD4-/CD8- T cells, TCRγ/δ T cells, CD4+ NKT cells, CD8+ NKT cells, classic NK cells, CD56dim/CD98dim cells, B cells, plasmablasts, monocytes, CD11cdim/CD172dim cells, myeloid dendritic cells (mDCs), and plasmacytoid dendritic cells (pDCs). Seven of the 17 main cell types exhibited statistically significant frequencies in the investigated groups. The expression levels of the 34 markers in the main populations were compared between HCs and SADs. In summary, 59 scatter plots showed significant differences in the expression intensities between at least two groups. Next, each immune cell population was divided into subpopulations (metaclusters) using the FlowSOM (self-organizing map) algorithm. Finally, 121 metaclusters (MCs) of the 10 main immune cell populations were found to have significant differences to classify diseases. The single-cell T-cell heterogeneity represented 64MCs based on the expression of 34 markers, and the frequency of 23 MCs differed significantly between at least twoconditions. The CD3- non-T-cell compartment contained 57 MCs with 17 MCs differentiating at least two investigated groups. In summary, we are the first to demonstrate the complexity of the immunophenotype of 34 markers over 15 million single cells in HCs vs. therapy-naive patients with RA, SSc, and SLE. Disease specific population frequencies or expression patterns of peripheral immune cells provide a single-cell data resource to the scientific community.

ARTseq-FISH reveals position-dependent differences in gene expression of micropatterned mESCs.

Differences in gene-expression profiles between individual cells can give rise to distinct cell fate decisions. Yet how localisation on a micropattern impacts initial changes in mRNA, protein, and phosphoprotein abundance remains unclear. To identify the effect of cellular position on gene expression, we developed a scalable antibody and mRNA targeting sequential fluorescence in situ hybridisation (ARTseq-FISH) method capable of simultaneously profiling mRNAs, proteins, and phosphoproteins in single cells. We studied 67 (phospho-)protein and mRNA targets in individual mouse embryonic stem cells (mESCs) cultured on circular micropatterns. ARTseq-FISH reveals relative changes in both abundance and localisation of mRNAs and (phospho-)proteins during the first 48 hours of exit from pluripotency. We confirm these changes by conventional immunofluorescence and time-lapse microscopy. Chemical labelling, immunofluorescence, and single-cell time-lapse microscopy further show that cells closer to the edge of the micropattern exhibit increased proliferation compared to cells at the centre. Together these data suggest that while gene expression is still highly heterogeneous position-dependent differences in mRNA and protein levels emerge as early as 12 hours after LIF withdrawal.

CASi: A framework for cross-timepoint analysis of single-cell RNA sequencing data.

Single-cell RNA sequencing (scRNA-seq) technology has been widely used to study the differences in gene expression at the single cell level, providing insights into the research of cell development, differentiation, and functional heterogeneity. Various pipelines and workflows of scRNA-seq analysis have been developed but few considered multi-timepoint data specifically. In this study, we develop CASi, a comprehensive framework for analyzing multiple timepoints' scRNA-seq data, which provides users with: (1) cross-timepoint cell annotation, (2) detection of potentially novel cell types emerged over time, (3) visualization of cell population evolution, and (4) identification of temporal differentially expressed genes (tDEGs). Through comprehensive simulation studies and applications to a real multi-timepoint single cell dataset, we demonstrate the robust and favorable performance of the proposal versus existing methods serving similar purposes.

Biclustering analysis on tree-shaped time-series single cell gene expression data of Caenorhabditis elegans.

BACKGROUND: In recent years, gene clustering analysis has become a widely used tool for studying gene functions, efficiently categorizing genes with similar expression patterns to aid in identifying gene functions. Caenorhabditis elegans is commonly used in embryonic research due to its consistent cell lineage from fertilized egg to adulthood. Biologists use 4D confocal imaging to observe gene expression dynamics at the single-cell level. However, on one hand, the observed tree-shaped time-series datasets have characteristics such as non-pairwise data points between different individuals. On the other hand, the influence of cell type heterogeneity should also be considered during clustering, aiming to obtain more biologically significant clustering results. RESULTS: A biclustering model is proposed for tree-shaped single-cell gene expression data of Caenorhabditis elegans. Detailedly, a tree-shaped piecewise polynomial function is first employed to fit non-pairwise gene expression time series data. Then, four factors are considered in the objective function, including Pearson correlation coefficients capturing gene correlations, p-values from the Kolmogorov-Smirnov test measuring the similarity between cells, as well as gene expression size and bicluster overlapping size. After that, Genetic Algorithm is utilized to optimize the function. CONCLUSION: The results on the small-scale dataset analysis validate the feasibility and effectiveness of our model and are superior to existing classical biclustering models. Besides, gene enrichment analysis is employed to assess the results on the complete real dataset analysis, confirming that the discovered biclustering results hold significant biological relevance.

Data-driven selection of analysis decisions in single-cell RNA-seq trajectory inference.

Single-cell RNA sequencing (scRNA-seq) experiments have become instrumental in developmental and differentiation studies, enabling the profiling of cells at a single or multiple time-points to uncover subtle variations in expression profiles reflecting underlying biological processes. Benchmarking studies have compared many of the computational methods used to reconstruct cellular dynamics; however, researchers still encounter challenges in their analysis due to uncertainty with respect to selecting the most appropriate methods and parameters. Even among universal data processing steps used by trajectory inference methods such as feature selection and dimension reduction, trajectory methods' performances are highly dataset-specific. To address these challenges, we developed Escort, a novel framework for evaluating a dataset's suitability for trajectory inference and quantifying trajectory properties influenced by analysis decisions. Escort evaluates the suitability of trajectory analysis and the combined effects of processing choices using trajectory-specific metrics. Escort navigates single-cell trajectory analysis through these data-driven assessments, reducing uncertainty and much of the decision burden inherent to trajectory inference analyses. Escort is implemented in an accessible R package and R/Shiny application, providing researchers with the necessary tools to make informed decisions during trajectory analysis and enabling new insights into dynamic biological processes at single-cell resolution.

Single-Cell RNA Sequencing Reveals Cardiac Fibroblast-Specific Transcriptomic Changes in Dilated Cardiomyopathy.

Dilated cardiomyopathy (DCM) is the most common cause of heart failure, with a complex aetiology involving multiple cell types. We aimed to detect cell-specific transcriptomic alterations in DCM through analysis that leveraged recent advancements in single-cell analytical tools. Single-cell RNA sequencing (scRNA-seq) data from human DCM cardiac tissue were subjected to an updated bioinformatic workflow in which unsupervised clustering was paired with reference label transfer to more comprehensively annotate the dataset. Differential gene expression was detected primarily in the cardiac fibroblast population. Bulk RNA sequencing was performed on an independent cohort of human cardiac tissue and compared with scRNA-seq gene alterations to generate a stratified list of higher-confidence, fibroblast-specific expression candidates for further validation. Concordant gene dysregulation was confirmed in TGFß-induced fibroblasts. Functional assessment of gene candidates showed that AEBP1 may play a significant role in fibroblast activation. This unbiased approach enabled improved resolution of cardiac cell-type-specific transcriptomic alterations in DCM.

Identification and validation of a prognostic anoikis-related gene signature in papillary thyroid carcinoma by integrated analysis of single-cell and bulk RNA-sequencing.

Papillary thyroid carcinoma (PTC) prognosis may be deteriorated due to the metastases, and anoikis palys an essential role in the tumor metastasis. However, the potential effect of anoikis-related genes on the prognosis of PTC was unclear. The mRNA and clinical information were obtained from the cancer genome atlas database. Hub genes were identified and risk model was constructed using Cox regression analysis. Kaplan-Meier (K-M) curve was applied for the survival analysis. Immune infiltration and immune therapy response were calculated using CIBERSORT and TIDE. The identification of cell types and cell interaction was performed by Seurat, SingleR and CellChat packages. GO, KEGG, and GSVA were applied for the enrichment analysis. Protein-protein interaction network was constructed in STRING and Cytoscape. Drug sensitivity was assessed in GSCA. Based on bulk RNA data, we identified 4 anoikis-related risk signatures, which were oncogenes, and constructed a risk model. The enrichment analysis found high risk group was enriched in some immune-related pathways. High risk group had higher infiltration of Tregs, higher TIDE score and lower levels of monocytes and CD8 T cells. Based on scRNA data, we found that 4 hub genes were mainly expressed in monocytes and macrophages, and they interacted with T cells. Hub genes were significantly related to immune escape-related genes. Drug sensitivity analysis suggested that cyclin dependent kinase inhibitor 2A may be a better chemotherapy target. We constructed a risk model which could effectively and steadily predict the prognosis of PTC. We inferred that the immune escape may be involved in the development of PTC.

Single-cell sequencing of tumour infiltrating T cells efficiently identifies tumour-specific T cell receptors based on the T cell activation score.

Adoptively transferred T cell receptor-engineered T cells are a promising cancer treatment strategy, and the identification of tumour-specific TCRs is essential. Previous studies reported that tumour-reactive T cells and TCRs could be isolated based on the expression of activation markers. However, since T cells with different cell states could not respond uniformly to activation but show a heterogeneous expression profile of activation and effector molecules, isolation of tumour-reactive T cells based on single activation or effector molecules could result in the absence of tumour-reactive T cells; thus, combinations of multiple activation and effector molecules could improve the efficiency of isolating tumour-specific TCRs. We enrolled two patients with lung adenocarcinoma and obtained their tumour infiltrating lymphocytes (TILs) and autologous tumour cells (ATCs). TILs were cocultured with the corresponding ATCs for 12 h and subjected to single-cell RNA sequencing. First, we identified three TCRs with the highest expression levels of IFNG and TNFRSF9 mRNA for each patient, yet only the top one or two recognized the corresponding ATCs in each patient. Next, we defined the activation score based on normalized expression levels of IFNG, IL2, TNF, IL2RA, CD69, TNFRSF9, GZMB, GZMA, GZMK, and PRF1 mRNA for each T cell and then identified three TCRs with the highest activation score for each patient. We found that all three TCRs in each patient could specifically identify corresponding ATCs. In conclusion, we established an efficient approach to isolate tumour-reactive TCRs based on combinations of multiple activation and effector molecules through single-cell RNA sequencing.

Multiplexed bulk and single-cell RNA-seq hybrid enables cost-efficient disease modeling with chimeric organoids.

Disease modeling with isogenic Induced Pluripotent Stem Cell (iPSC)-differentiated organoids serves as a powerful technique for studying disease mechanisms. Multiplexed coculture is crucial to mitigate batch effects when studying the genetic effects of disease-causing variants in differentiated iPSCs or organoids, and demultiplexing at the single-cell level can be conveniently achieved by assessing natural genetic barcodes. Here, to enable cost-efficient time-series experimental designs via multiplexed bulk and single-cell RNA-seq of hybrids, we introduce a computational method in our Vireo Suite, Vireo-bulk, to effectively deconvolve pooled bulk RNA-seq data by genotype reference, and thereby quantify donor abundance over the course of differentiation and identify differentially expressed genes among donors. Furthermore, with multiplexed scRNA-seq and bulk RNA-seq, we demonstrate the usefulness and necessity of a pooled design to reveal donor iPSC line heterogeneity during macrophage cell differentiation and to model rare WT1 mutation-driven kidney disease with chimeric organoids. Our work provides an experimental and analytic pipeline for dissecting disease mechanisms with chimeric organoids.

Modification of Seurat v4 for the Development of a Phase Assignment Tool Able to Distinguish between G2 and Mitotic Cells.

Single-cell RNA sequencing (scRNAseq) is a rapidly advancing field enabling the characterisation of heterogeneous gene expression profiles within a population. The cell cycle phase is a major contributor to gene expression variance between cells and computational analysis tools have been developed to assign cell cycle phases to cells within scRNAseq datasets. Whilst these tools can be extremely useful, all have the drawback that they classify cells as only G1, S or G2/M. Existing discrete cell phase assignment tools are unable to differentiate between G2 and M and continuous-phase-assignment tools are unable to identify a region corresponding specifically to mitosis in a pseudo-timeline for continuous assignment along the cell cycle. In this study, bulk RNA sequencing was used to identify differentially expressed genes between mitotic and interphase cells isolated based on phospho-histone H3 expression using fluorescence-activated cell sorting. These gene lists were used to develop a methodology which can distinguish G2 and M phase cells in scRNAseq datasets. The phase assignment tools present in Seurat were modified to allow for cell cycle phase assignment of all stages of the cell cycle to identify a mitotic-specific cell population.

Single-Cell Transcriptional Profile Construction of Rat Pituitary Glands before and after Sexual Maturation and Identification of Novel Marker Spp1 in Gonadotropes.

The mammalian pituitary gland drives highly conserved physiological processes such as somatic cell growth, pubertal transformation, fertility, and metabolism by secreting a variety of hormones. Recently, single-cell transcriptomics techniques have been used in pituitary gland research. However, more studies have focused on adult pituitary gland tissues from different species or different sexes, and no research has yet resolved cellular differences in pituitary gland tissue before and after sexual maturation. Here, we identified a total of 15 cell clusters and constructed single-cell transcriptional profiles of rats before and after sexual maturation. Furthermore, focusing on the gonadotrope cluster, 106 genes were found to be differentially expressed before and after sexual maturation. It was verified that Spp1, which is specifically expressed in gonadotrope cells, could serve as a novel marker for this cell cluster and has a promotional effect on the synthesis and secretion of follicle-stimulating hormone. The results provide a new resource for further resolving the regulatory mechanism of pituitary gland development and pituitary hormone synthesis and secretion.

Pancreatic ß-Cell Identity Change through the Lens of Single-Cell Omics Research.

The main hallmark in the development of both type 1 and type 2 diabetes is a decline in functional ß-cell mass. This decline is predominantly attributed to ß-cell death, although recent findings suggest that the loss of ß-cell identity may also contribute to ß-cell dysfunction. This phenomenon is characterized by a reduced expression of key markers associated with ß-cell identity. This review delves into the insights gained from single-cell omics research specifically focused on ß-cell identity. It highlights how single-cell omics based studies have uncovered an unexpected level of heterogeneity among ß-cells and have facilitated the identification of distinct ß-cell subpopulations through the discovery of cell surface markers, transcriptional regulators, the upregulation of stress-related genes, and alterations in chromatin activity. Furthermore, specific subsets of ß-cells have been identified in diabetes, such as displaying an immature, dedifferentiated gene signature, expressing significantly lower insulin mRNA levels, and expressing increased ß-cell precursor markers. Additionally, single-cell omics has increased insight into the detrimental effects of diabetes-associated conditions, including endoplasmic reticulum stress, oxidative stress, and inflammation, on ß-cell identity. Lastly, this review outlines the factors that may influence the identification of ß-cell subpopulations when designing and performing a single-cell omics experiment.

Revealing Genetic Dynamics: scRNA-seq Unravels Modifications in Human PDL Cells across In Vivo and In Vitro Environments.

The periodontal ligament (PDL) is a highly specialized fibrous tissue comprising heterogeneous cell populations of an intricate nature. These complexities, along with challenges due to cell culture, impede a comprehensive understanding of periodontal pathophysiology. This study aims to address this gap, employing single-cell RNA sequencing (scRNA-seq) technology to analyze the genetic intricacies of PDL both in vivo and in vitro. Primary human PDL samples (n = 7) were split for direct in vivo analysis and cell culture under serum-containing and serum-free conditions. Cell hashing and sorting, scRNA-seq library preparation using the 10x Genomics protocol, and Illumina sequencing were conducted. Primary analysis was performed using Cellranger, with downstream analysis via the R packages Seurat and SCORPIUS. Seven distinct PDL cell clusters were identified comprising different cellular subsets, each characterized by unique genetic profiles, with some showing donor-specific patterns in representation and distribution. Formation of these cellular clusters was influenced by culture conditions, particularly serum presence. Furthermore, certain cell populations were found to be inherent to the PDL tissue, while others exhibited variability across donors. This study elucidates specific genes and cell clusters within the PDL, revealing both inherent and context-driven subpopulations. The impact of culture conditions-notably the presence of serum-on cell cluster formation highlights the critical need for refining culture protocols, as comprehending these influences can drive the creation of superior culture systems vital for advancing research in PDL biology and regenerative therapies. These discoveries not only deepen our comprehension of PDL biology but also open avenues for future investigations into uncovering underlying mechanisms.

Coupling imaging mass cytometry with Alcian blue histochemical staining for a single-slide approach.

Imaging mass cytometry (IMC) is a metal mass spectrometry-based method allowing highly multiplex immunophenotyping of cells within tissue samples. However, some limitations of IMC are its 1-µm resolution and its time and costs of analysis limiting respectively the detailed histopathological analysis of IMC-produced images and its application to small selected tissue regions of interest (ROI) of one to few square millimeters. Coupling on a single-tissue section, IMC and histopathological analyses could permit a better selection of the ROI for IMC analysis as well as co-analysis of immunophenotyping and histopathological data until the single-cell level. The development of this method is the aim of the present study in which we point to the feasibility of applying the IMC process to tissue sections previously Alcian blue-stained and digitalized before IMC tissue destructive analyses. This method could help to improve the process of IMC in terms of ROI selection, time of analysis, and the confrontation between histopathological and immunophenotypic data of cells.

The future of rapid and automated single-cell data analysis using reference mapping.

As the number of single-cell datasets continues to grow rapidly, workflows that map new data to well-curated reference atlases offer enormous promise for the biological community. In this perspective, we discuss key computational challenges and opportunities for single-cell reference-mapping algorithms. We discuss how mapping algorithms will enable the integration of diverse datasets across disease states, molecular modalities, genetic perturbations, and diverse species and will eventually replace manual and laborious unsupervised clustering pipelines.

Cell metabolism: Functional and phenotypic single cell approaches.

Several metabolic pathways are essential for the physiological regulation of immune cells, but their dysregulation can cause immune dysfunction. Hypermetabolic and hypometabolic states represent deviations in the magnitude and flexibility of effector cells in different contexts, for example in autoimmunity, infections or cancer. To study immunometabolism, most methods focus on bulk populations and rely on in vitro activation assays. Nowadays, thanks to the development of single-cell technologies, including multiparameter flow cytometry, mass cytometry, RNA cytometry, among others, the metabolic state of individual immune cells can be measured in a variety of samples obtained in basic, translational and clinical studies. Here, we provide an overview of different single-cell approaches that are employed to investigate both mitochondrial functions and cell dependence from mitochondria metabolism. Moreover, besides the description of the appropriate experimental settings, we discuss the strengths and weaknesses of different approaches with the aim to suggest how to study cell metabolism in the settings of interest.

Spectral analysis and sorting of microbial organisms using a spectral sorter.

This chapter discusses the problems related to the application of conventional flow cytometers to microbiology. To address some of those limitations, the concept of spectral flow cytometry is introduced and the advantages over conventional flow cytometry for bacterial sorting are presented. We demonstrate by using ThermoFisher's Bigfoot spectral sorter where the spectral signatures of different stains for staining bacteria are demonstrated with an example of performing unmixing on spectral datasets. In addition to the Bigfoot's spectral analysis, the special biosafety features of this instrument are discussed. Utilizing these biosafety features, the sorting and patterning at the single cell level is optimized using non-pathogenic bacteria. Finally, the chapter is concluded by presenting a novel, label free, non-destructive, and rapid phenotypic method called Elastic Light Scattering (ELS) technology for identification of the patterned bacterial cells based on their unique colony scatter patterns.

Harmonized analysis of PBMC by mass cytometry.

Mass cytometry permits the high dimensional analysis of cellular systems at single-cell resolution with high throughput in various areas of biomedical research. Here, we provide a state-of-the-art protocol for the analysis of human peripheral blood mononuclear cells (PBMC) by mass cytometry. We focus on the implementation of measures promoting the harmonization of large and complex studies to aid robustness and reproducibility of immune phenotyping data.