Single-cell trajectory detection uncovers progression and regulatory coordination in human B cell development

Proteomics appending a complementary dimension to precision oncotherapy

Recent advances in high-throughput proteomic profiling technologies have facilitated the precise quantification of numerous proteins across multiple specimens concurrently. Researchers have the opportunity to comprehensively analyze the molecular signatures in plentiful medical specimens or disease pattern cell lines. Along with advances in data analysis and integration, proteomics data could be efficiently consolidated and employed to recognize precise elementary molecular mechanisms and decode individual biomarkers, guiding the precision treatment of tumors. Herein, we review a broad array of proteomics technologies and the progress and methods for the integration of proteomics data and further discuss how to better merge proteomics in precision medicine and clinical settings.

Single-cell proteomics delineates murine systemic immune response to blast lung injury

Victims of explosive events frequently suffer from blast lung injuries. Immune system has been implicated in the pathogenesis of this disease. However, systemic immune responses underlying the progression and recovery of injury repair remain poorly understood. Here, we depict the systemic landscape of immune dysregulation during blast lung injury and uncover immune recovery patterns. Single-cell analyses reveal dramatic changes in neutrophils, macrophages, monocytes, dendritic cells, and eosinophils after a gas explosion, along with early involvement of CD4 T, CD8 T, and Th17 cells. We demonstrate that myeloid cells primarily exert functions during the acute phase, while the spleen serves as an alternative source of granulocytes. Granulopoiesis is initiated in the bone marrow at a later stage during blast lung injury recovery, rather than at the acute stage. These findings contribute to a better understanding of the pathogenesis and provide valuable insights for potential immune interventions in blast lung injury.

Mapping lineage-traced cells across time points with moslin

Simultaneous profiling of single-cell gene expression and lineage history holds enormous potential for studying cellular decision-making. Recent computational approaches combine both modalities into cellular trajectories; however, they cannot make use of all available lineage information in destructive time-series experiments. Here, we present moslin, a Gromov-Wasserstein-based model to couple cellular profiles across time points based on lineage and gene expression information. We validate our approach in simulations and demonstrate on Caenorhabditis elegans embryonic development how moslin predicts fate probabilities and putative decision driver genes. Finally, we use moslin to delineate lineage relationships among transiently activated fibroblast states during zebrafish heart regeneration.

Application of Mass Cytometry Platforms to Solid Organ Transplantation

Transplantation serves as the cornerstone of treatment for patients with end-stage organ disease. The prevalence of complications, such as allograft rejection, infection, and malignancies, underscores the need to dissect the complex interactions of the immune system at the single-cell level. In this review, we discuss studies using mass cytometry or cytometry by time-of-flight, a cutting-edge technology enabling the characterization of immune populations and cell-to-cell interactions in granular detail. We review the application of mass cytometry in human and experimental animal studies in the context of transplantation, uncovering invaluable contributions of the tool to understanding rejection and other transplant-related complications. We discuss recent innovations that have the potential to streamline and standardize mass cytometry workflows for application to multisite clinical trials. Additionally, we introduce imaging mass cytometry, a technique that couples the power of mass cytometry with spatial context, thereby mapping cellular interactions within tissue microenvironments. The synergistic integration of mass cytometry and imaging mass cytometry data with other omics data sets and high-dimensional data platforms to further define immune dynamics is discussed. In conclusion, mass cytometry technologies, when integrated with other tools and data, shed light on the intricate landscape of the immune response in transplantation. This approach holds significant potential for enhancing patient outcomes by advancing our understanding and facilitating the development of new diagnostics and therapeutics.

The pericardium forms as a distinct structure during heart formation

The heart integrates diverse cell lineages into a functional unit, including the pericardium, a mesothelial sac that supports heart movement, homeostasis, and immune responses. However, despite its critical roles, the developmental origins of the pericardium remain uncertain due to disparate models. Here, using live imaging, lineage tracking, and single-cell transcriptomics in zebrafish, we find the pericardium forms within the lateral plate mesoderm from dedicated anterior mesothelial progenitors and distinct from the classic heart field. Imaging of transgenic reporters in zebrafish documents lateral plate mesoderm cells that emerge lateral of the classic heart field and among a continuous mesothelial progenitor field. Single-cell transcriptomics and trajectories of hand2-expressing lateral plate mesoderm reveal distinct populations of mesothelial and cardiac precursors, including pericardial precursors that are distinct from the cardiomyocyte lineage. The mesothelial gene expression signature is conserved in mammals and carries over to post-natal development. Light sheet-based live-imaging and machine learning-supported cell tracking documents that during heart tube formation, pericardial precursors that reside at the anterior edge of the heart field migrate anteriorly and medially before fusing, enclosing the embryonic heart to form a single pericardial cavity. Pericardium formation proceeds even upon genetic disruption of heart tube formation, uncoupling the two structures. Canonical Wnt/β-catenin signaling modulates pericardial cell number, resulting in a stretched pericardial epithelium with reduced cell number upon canonical Wnt inhibition. We connect the pathological expression of secreted Wnt antagonists of the SFRP family found in pediatric dilated cardiomyopathy to increased pericardial stiffness: sFRP1 in the presence of increased catecholamines causes cardiomyocyte stiffness in neonatal rats as measured by atomic force microscopy. Altogether, our data integrate pericardium formation as an independent process into heart morphogenesis and connect disrupted pericardial tissue properties such as pericardial stiffness to pediatric cardiomyopathies.

Systemic deficits in lipid homeostasis promote aging-associated impairments in B cell progenitor development

Organismal aging has been associated with diverse metabolic and functional changes across tissues. Within the immune system, key features of physiological hematopoietic cell aging include increased fat deposition in the bone marrow, impaired hematopoietic stem and progenitor cell (HSPC) function, and a propensity towards myeloid differentiation. This shift in lineage bias can lead to pre-malignant bone marrow conditions such as clonal hematopoiesis of indeterminate potential (CHIP) or clonal cytopenias of undetermined significance (CCUS), frequently setting the stage for subsequent development of age-related cancers in myeloid or lymphoid lineages. At the systemic as well as sub-cellular level, human aging has also been associated with diverse lipid alterations, such as decreased phospholipid membrane fluidity that arises as a result of increased saturated fatty acid (FA) accumulation and a decay in n-3 polyunsaturated fatty acid (PUFA) species by the age of 80 years, however the extent to which impaired FA metabolism contributes to hematopoietic aging is less clear. Here, we performed comprehensive multi-omics analyses and uncovered a role for a key PUFA biosynthesis gene, ELOVL2 , in mouse and human immune cell aging. Whole transcriptome RNA-sequencing studies of bone marrow from aged Elovl2 mutant (enzyme-deficient) mice compared with age-matched controls revealed global down-regulation in lymphoid cell markers and expression of genes involved specifically in B cell development. Flow cytometric analyses of immune cell markers confirmed an aging-associated loss of B cell markers that was exacerbated in the bone marrow of Elovl2 mutant mice and unveiled CD79B, a vital molecular regulator of lymphoid progenitor development from the pro-B to pre-B cell stage, as a putative surface biomarker of accelerated immune aging. Complementary lipidomic studies extended these findings to reveal select alterations in lipid species in aged and Elovl2 mutant mouse bone marrow samples, suggesting significant changes in the biophysical properties of cellular membranes. Furthermore, single cell RNA-seq analysis of human HSPCs across the spectrum of human development and aging uncovered a rare subpopulation (<7%) of CD34 + HSPCs that expresses ELOVL2 in healthy adult bone marrow. This HSPC subset, along with CD79B -expressing lymphoid-committed cells, were almost completely absent in CD34 + cells isolated from elderly (>60 years old) bone marrow samples. Together, these findings uncover new roles for lipid metabolism enzymes in the molecular regulation of cellular aging and immune cell function in mouse and human hematopoiesis. In addition, because systemic loss of ELOVL2 enzymatic activity resulted in down-regulation of B cell genes that are also associated with lymphoproliferative neoplasms, this study sheds light on an intriguing metabolic pathway that could be leveraged in future studies as a novel therapeutic modality to target blood cancers or other age-related conditions involving the B cell lineage.

Inferring cell differentiation maps from lineage tracing data

During development, mulitpotent cells differentiate through a hierarchy of increasingly restricted progenitor cell types until they realize specialized cell types. A cell differentiation map describes this hierarchy, and inferring these maps is an active area of research spanning traditional single marker lineage studies to data-driven trajectory inference methods on single-cell RNA-seq data. Recent high-throughput lineage tracing technologies profile lineages and cell types at scale, but current methods to infer cell differentiation maps from these data rely on simple models with restrictive assumptions about the developmental process. We introduce a mathematical framework for cell differentiation maps based on the concept of potency, and develop an algorithm, Carta, that infers an optimal cell differentiation map from single-cell lineage tracing data. The key insight in Carta is to balance the trade-off between the complexity of the cell differentiation map and the number of unobserved cell type transitions on the lineage tree. We show that Carta more accurately infers cell differentiation maps on both simulated and real data compared to existing methods. In models of mammalian trunk development and mouse hematopoiesis, Carta identifies important features of development that are not revealed by other methods including convergent differentiation of specialized cell types, progenitor differentiation dynamics, and the refinement of routes of differentiation via new intermediate progenitors.

Tviblindi algorithm identifies branching developmental trajectories of human B-cell development and describes abnormalities in RAG-1 and WAS patients

Detailed knowledge of human B-cell development is crucial for the proper interpretation of inborn errors of immunity and malignant diseases. It is of interest to understand the kinetics of protein expression changes during development, but also to properly interpret the major and possibly alternative developmental trajectories. We have investigated human samples from healthy individuals with the aim of describing all B-cell developmental trajectories. We validated a 30-parameter mass cytometry panel and demonstrated the utility of "vaevictis" visualization of B-cell developmental stages. We used the trajectory inference tool "tviblindi" to exhaustively describe all trajectories leading to all developmental ends discovered in the data. Focusing on Natural Effector B cells, we demonstrated the dynamics of expression of nuclear factors (PAX-5, TdT, Ki-67, Bcl-2), cytokine and chemokine receptors (CD127, CXCR4, CXCR5) in relation to the canonical B-cell developmental stage markers. We observed branching of the memory development, where follicular memory formation was marked by CD73 expression. Lastly, we performed an analysis of two example cases of abnormal B-cell development caused by mutations in RAG-1 and Wiskott-Aldrich syndrome gene in patients with primary immunodeficiency. In conclusion, we developed, validated, and presented a comprehensive set of tools for the investigation of B-cell development in the bone marrow compartment.

Revealing gene function with statistical inference at single-cell resolution

Single-cell and spatial molecular profiling assays have shown large gains in sensitivity, resolution and throughput. Applying these technologies to specimens from human and model organisms promises to comprehensively catalogue cell types, reveal their lineage origins in development and discern their contributions to disease pathogenesis. Moreover, rapidly dropping costs have made well-controlled perturbation experiments and cohort studies widely accessible, illuminating mechanisms that give rise to phenotypes at the scale of the cell, the tissue and the whole organism. Interpreting the coming flood of single-cell data, much of which will be spatially resolved, will place a tremendous burden on existing computational pipelines. However, statistical concepts, models, tools and algorithms can be repurposed to solve problems now arising in genetic and molecular biology studies of development and disease. Here, I review how the questions that recent technological innovations promise to answer can be addressed by the major classes of statistical tools.

Deciphering bone marrow engraftment after allogeneic stem cell transplantation in humans using single-cell analyses

BACKGROUNDDonor cell engraftment is a prerequisite of successful allogeneic hematopoietic stem cell transplantation. Based on peripheral blood analyses, it is characterized by early myeloid recovery and T and B cell lymphopenia. However, cellular networks associated with bone marrow engraftment of allogeneic human cells have been poorly described.METHODSMass cytometry and CITE-Seq analyses were performed on bone marrow cells 3 months after transplantation in patients with acute myelogenous leukemia.RESULTSMass cytometric analyses in 26 patients and 20 healthy controls disclosed profound alterations in myeloid and B cell progenitors, with a shift toward terminal myeloid differentiation and decreased B cell progenitors. Unsupervised analysis separated recipients into 2 groups, one of them being driven by previous graft-versus-host disease (R2 patients). We then used single-cell CITE-Seq to decipher engraftment, which resolved 36 clusters, encompassing all bone marrow cellular components. Hematopoiesis in transplant recipients was sustained by committed myeloid and erythroid progenitors in a setting of monocyte-, NK cell-, and T cell-mediated inflammation. Gene expression revealed major pathways in transplant recipients, namely, TNF-α signaling via NF-κB and the IFN-γ response. The hallmark of allograft rejection was consistently found in clusters from transplant recipients, especially in R2 recipients.CONCLUSIONBone marrow cell engraftment of allogeneic donor cells is characterized by a state of emergency hematopoiesis in the setting of an allogeneic response driving inflammation.FUNDINGThis study was supported by the French National Cancer Institute (Institut National du Cancer; PLBIO19-239) and by an unrestricted research grant by Alexion Pharmaceuticals.

Single-cell multiomics analysis reveals cell/tissue-specific associations in bipolar disorder

This study investigates the cellular origin and tissue heterogeneity in bipolar disorder (BD) by integrating multiomics data. Four distinct datasets were employed, including single-cell RNA sequencing (scRNA-seq) data (embryonic and fetal brain, n = 8, 1,266 cells), BD Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq) data (adult brain, n = 210), BD bulk RNA-seq data (adult brain, n = 314), and BD genome-wide association study (GWAS) summary data (n = 413,466). The integration of scRNA-seq data with multiomics data relevant to BD was accomplished using the single-cell disease relevance score (scDRS) algorithm. We have identified a novel brain cell cluster named ADCY1, which exhibits distinct genetic characteristics. From a high-resolution genetic perspective, glial cells emerge as the primary cytopathology associated with BD. Specifically, astrocytes were significantly related to BD at the RNA-seq level, while microglia showed a strong association with BD across multiple panels, including the transcriptome-wide association study (TWAS), ATAC-seq, and RNA-seq. Additionally, oligodendrocyte precursor cells displayed a significant association with BD in both ATAC-seq and RNA-seq panel. Notably, our investigation of brain regions affected by BD revealed significant associations between BD and all three types of glial cells in the dorsolateral prefrontal cortex (DLPFC). Through comprehensive analyses, we identified several BD-associated genes, including CRMP1, SYT4, UCHL1, and ZBTB18. In conclusion, our findings suggest that glial cells, particularly in specific brain regions such as the DLPFC, may play a significant role in the pathogenesis of BD. The integration of multiomics data has provided valuable insights into the etiology of BD, shedding light on potential mechanisms underlying this complex psychiatric disorder.

scHiCyclePred: a deep learning framework for predicting cell cycle phases from single-cell Hi-C data using multi-scale interaction information

The emergence of single-cell Hi-C (scHi-C) technology has provided unprecedented opportunities for investigating the intricate relationship between cell cycle phases and the three-dimensional (3D) structure of chromatin. However, accurately predicting cell cycle phases based on scHi-C data remains a formidable challenge. Here, we present scHiCyclePred, a prediction model that integrates multiple feature sets to leverage scHi-C data for predicting cell cycle phases. scHiCyclePred extracts 3D chromatin structure features by incorporating multi-scale interaction information. The comparative analysis illustrates that scHiCyclePred surpasses existing methods such as Nagano_method and CIRCLET across various metrics including accuracy (ACC), F1 score, Precision, Recall, and balanced accuracy (BACC). In addition, we evaluate scHiCyclePred against the previously published CIRCLET using the dataset of complex tissues (Liu_dataset). Experimental results reveal significant improvements with scHiCyclePred exhibiting improvements of 0.39, 0.52, 0.52, and 0.39 over the CIRCLET in terms of ACC, F1 score, Precision, and Recall metrics, respectively. Furthermore, we conduct analyses on three-dimensional chromatin dynamics and gene features during the cell cycle, providing a more comprehensive understanding of cell cycle dynamics through chromatin structure. scHiCyclePred not only offers insights into cell biology but also holds promise for catalyzing breakthroughs in disease research. Access scHiCyclePred on GitHub at https:// github.com/HaoWuLab-Bioinformatics/ scHiCyclePred .

Single-cell omics: experimental workflow, data analyses and applications

Cells are the fundamental units of biological systems and exhibit unique development trajectories and molecular features. Our exploration of how the genomes orchestrate the formation and maintenance of each cell, and control the cellular phenotypes of various organismsis, is both captivating and intricate. Since the inception of the first single-cell RNA technology, technologies related to single-cell sequencing have experienced rapid advancements in recent years. These technologies have expanded horizontally to include single-cell genome, epigenome, proteome, and metabolome, while vertically, they have progressed to integrate multiple omics data and incorporate additional information such as spatial scRNA-seq and CRISPR screening. Single-cell omics represent a groundbreaking advancement in the biomedical field, offering profound insights into the understanding of complex diseases, including cancers. Here, we comprehensively summarize recent advances in single-cell omics technologies, with a specific focus on the methodology section. This overview aims to guide researchers in selecting appropriate methods for single-cell sequencing and related data analysis.

Automated cell type annotation and exploration of single-cell signaling dynamics using mass cytometry

Mass cytometry by time-of-flight (CyTOF) is an emerging technology allowing for in-depth characterization of cellular heterogeneity in cancer and other diseases. Unfortunately, high-dimensional analyses of CyTOF data remain quite demanding. Here, we deploy a bioinformatics framework that tackles two fundamental problems in CyTOF analyses namely (1) automated annotation of cell populations guided by a reference dataset and (2) systematic utilization of single-cell data for effective patient stratification. By applying this framework on several publicly available datasets, we demonstrate that the Scaffold approach achieves good trade-off between sensitivity and specificity for automated cell type annotation. Additionally, a case study focusing on a cohort of 43 leukemia patients reported salient interactions between signaling proteins that are sufficient to predict short-term survival at time of diagnosis using the XGBoost algorithm. Our work introduces an automated and versatile analysis framework for CyTOF data with many applications in future precision medicine projects.

Terminal deoxynucleotidyl transferase and CD84 identify human multi-potent lymphoid progenitors

Lymphoid specification in human hematopoietic progenitors is not fully understood. To better associate lymphoid identity with protein-level cell features, we conduct a highly multiplexed single-cell proteomic screen on human bone marrow progenitors. This screen identifies terminal deoxynucleotidyl transferase (TdT), a specialized DNA polymerase intrinsic to VDJ recombination, broadly expressed within CD34+ progenitors prior to B/T cell emergence. While these TdT+ cells coincide with granulocyte-monocyte progenitor (GMP) immunophenotype, their accessible chromatin regions show enrichment for lymphoid-associated transcription factor (TF) motifs. TdT expression on GMPs is inversely related to the SLAM family member CD84. Prospective isolation of CD84lo GMPs demonstrates robust lymphoid potentials ex vivo, while still retaining significant myeloid differentiation capacity, akin to LMPPs. This multi-omic study identifies human bone marrow lymphoid-primed progenitors, further defining the lympho-myeloid axis in human hematopoiesis.

Are We There Yet? Assessing the Readiness of Single-Cell Proteomics to Answer Biological Hypotheses

Single-cell analysis is an active area of research in many fields of biology. Measurements at single-cell resolution allow researchers to study diverse populations without losing biologically meaningful information to sample averages. Many technologies have been used to study single cells, including mass spectrometry-based single-cell proteomics (SCP). SCP has seen a lot of growth over the past couple of years through improvements in data acquisition and analysis, leading to greater proteomic depth. Because method development has been the main focus in SCP, biological applications have been sprinkled in only as proof-of-concept. However, SCP methods now provide significant coverage of the proteome and have been implemented in many laboratories. Thus, a primary question to address in our community is whether the current state of technology is ready for widespread adoption for biological inquiry. In this Perspective, we examine the potential for SCP in three thematic areas of biological investigation: cell annotation, developmental trajectories, and spatial mapping. We identify that the primary limitation of SCP is sample throughput. As proteome depth has been the primary target for method development to date, we advocate for a change in focus to facilitate measuring tens of thousands of single-cell proteomes to enable biological applications beyond proof-of-concept.

Quantifying cell-state densities in single-cell phenotypic landscapes using Mellon

Cell-state density characterizes the distribution of cells along phenotypic landscapes and is crucial for unraveling the mechanisms that drive diverse biological processes. Here, we present Mellon, an algorithm for estimation of cell-state densities from high-dimensional representations of single-cell data. We demonstrate Mellon's efficacy by dissecting the density landscape of differentiating systems, revealing a consistent pattern of high-density regions corresponding to major cell types intertwined with low-density, rare transitory states. We present evidence implicating enhancer priming and the activation of master regulators in emergence of these transitory states. Mellon offers the flexibility to perform temporal interpolation of time-series data, providing a detailed view of cell-state dynamics during developmental processes. Mellon facilitates density estimation across various single-cell data modalities, scaling linearly with the number of cells. Our work underscores the importance of cell-state density in understanding the differentiation processes, and the potential of Mellon to provide insights into mechanisms guiding biological trajectories.

Evaluation of proteome dynamics: Implications for statistical confidence in mass spectrometric determination

Single-cell proteomics is currently far less productive than other approaches. Still, the proteomic community is having trouble adapting to the limitation of having to examine fewer cells than they would like. Studies on a small number of cells should be carefully planned to maximize the chances of success in this situation. This study aims to determine how sample size and measurement speed (slope)/variation affect the accuracy of a protein proteome mass spectrometric determination. The determination accuracy was shown to increase, and the false positive rate was shown to decrease as the sample size increased from 7 to 100 cells and the measurement slope/variation (S/V) ratio increased from 1 to 6. Furthermore, it was discovered that the number of cells in the sample increased the accuracy of this estimate. Thus, for 100 cells, the measurement S/V ratio was typically estimated to be very close to the real-world value, with a standard deviation of 0.35. For sample sizes from 7 to 100 cells, this accuracy was seen when calculating the measurement S/V ratio. The findings can help researchers plan experiments for mass spectroscopic protein proteome determination and other research purposes.

Optimal sequencing budget allocation for trajectory reconstruction of single cells

BACKGROUND

Charting cellular trajectories over gene expression is key to understanding dynamic cellular processes and their underlying mechanisms. While advances in single-cell RNA-sequencing technologies and computational methods have pushed forward the recovery of such trajectories, trajectory inference remains a challenge due to the noisy, sparse, and high-dimensional nature of single-cell data. This challenge can be alleviated by increasing either the number of cells sampled along the trajectory (breadth) or the sequencing depth, i.e. the number of reads captured per cell (depth). Generally, these two factors are coupled due to an inherent breadth-depth tradeoff that arises when the sequencing budget is constrained due to financial or technical limitations.

RESULTS

Here we study the optimal allocation of a fixed sequencing budget to optimize the recovery of trajectory attributes. Empirical results reveal that reconstruction accuracy of internal cell structure in expression space scales with the logarithm of either the breadth or depth of sequencing. We additionally observe a power law relationship between the optimal number of sampled cells and the corresponding sequencing budget. For linear trajectories, non-monotonicity in trajectory reconstruction across the breadth-depth tradeoff can impact downstream inference, such as expression pattern analysis along the trajectory. We demonstrate these results for five single-cell RNA-sequencing datasets encompassing differentiation of embryonic stem cells, pancreatic beta cells, hepatoblast and multipotent hematopoietic cells, as well as induced reprogramming of embryonic fibroblasts into neurons. By addressing the challenges of single-cell data, our study offers insights into maximizing the efficiency of cellular trajectory analysis through strategic allocation of sequencing resources.

Single-cell morphodynamical trajectories enable prediction of gene expression accompanying cell state change

Extracellular signals induce changes to molecular programs that modulate multiple cellular phenotypes, including proliferation, motility, and differentiation status. The connection between dynamically adapting phenotypic states and the molecular programs that define them is not well understood. Here we develop data-driven models of single-cell phenotypic responses to extracellular stimuli by linking gene transcription levels to "morphodynamics" - changes in cell morphology and motility observable in time-lapse image data. We adopt a dynamics-first view of cell state by grouping single-cell trajectories into states with shared morphodynamic responses. The single-cell trajectories enable development of a first-of-its-kind computational approach to map live-cell dynamics to snapshot gene transcript levels, which we term MMIST, Molecular and Morphodynamics-Integrated Single-cell Trajectories. The key conceptual advance of MMIST is that cell behavior can be quantified based on dynamically defined states and that extracellular signals alter the overall distribution of cell states by altering rates of switching between states. We find a cell state landscape that is bound by epithelial and mesenchymal endpoints, with distinct sequences of epithelial to mesenchymal transition (EMT) and mesenchymal to epithelial transition (MET) intermediates. The analysis yields predictions for gene expression changes consistent with curated EMT gene sets and provides a prediction of thousands of RNA transcripts through extracellular signal-induced EMT and MET with near-continuous time resolution. The MMIST framework leverages true single-cell dynamical behavior to generate molecular-level omics inferences and is broadly applicable to other biological domains, time-lapse imaging approaches and molecular snapshot data.

Direct comparison of mass cytometry and single-cell RNA sequencing of human peripheral blood mononuclear cells

Single-cell methods offer a high-resolution approach for characterizing cell populations. Many studies rely on single-cell transcriptomics to draw conclusions regarding cell state and behavior, with the underlying assumption that transcriptomic readouts largely parallel their protein counterparts and subsequent activity. However, the relationship between transcriptomic and proteomic measurements is imprecise, and thus datasets that probe the extent of their concordance will be useful to refine such conclusions. Additionally, novel single-cell analysis tools often lack appropriate gold standard datasets for the purposes of assessment. Integrative (combining the two data modalities) and predictive (using one modality to improve results from the other) approaches in particular, would benefit from transcriptomic and proteomic data from the same sample of cells. For these reasons, we performed single-cell RNA sequencing, mass cytometry, and flow cytometry on a split-sample of human peripheral blood mononuclear cells. We directly compare the proportions of specific cell types resolved by each technique, and further describe the extent to which protein and mRNA measurements correlate within distinct cell types.

Establishing a conceptual framework for holistic cell states and state transitions

Cell states were traditionally defined by how they looked, where they were located, and what functions they performed. In this post-genomic era, the field is largely focused on a molecular view of cell state. Moving forward, we anticipate that the observables used to define cell states will evolve again as single-cell imaging and analytics are advancing at a breakneck pace via the collection of large-scale, systematic cell image datasets and the application of quantitative image-based data science methods. This is, therefore, a key moment in the arc of cell biological research to develop approaches that integrate the spatiotemporal observables of the physical structure and organization of the cell with molecular observables toward the concept of a holistic cell state. In this perspective, we propose a conceptual framework for holistic cell states and state transitions that is data-driven, practical, and useful to enable integrative analyses and modeling across many data types.

PhyloVelo enhances transcriptomic velocity field mapping using monotonically expressed genes

Single-cell RNA sequencing (scRNA-seq) is a powerful approach for studying cellular differentiation, but accurately tracking cell fate transitions can be challenging, especially in disease conditions. Here we introduce PhyloVelo, a computational framework that estimates the velocity of transcriptomic dynamics by using monotonically expressed genes (MEGs) or genes with expression patterns that either increase or decrease, but do not cycle, through phylogenetic time. Through integration of scRNA-seq data with lineage information, PhyloVelo identifies MEGs and reconstructs a transcriptomic velocity field. We validate PhyloVelo using simulated data and Caenorhabditis elegans ground truth data, successfully recovering linear, bifurcated and convergent differentiations. Applying PhyloVelo to seven lineage-traced scRNA-seq datasets, generated using CRISPR-Cas9 editing, lentiviral barcoding or immune repertoire profiling, demonstrates its high accuracy and robustness in inferring complex lineage trajectories while outperforming RNA velocity. Additionally, we discovered that MEGs across tissues and organisms share similar functions in translation and ribosome biogenesis.

Highly Multiplexing, Throughput and Efficient Single-Cell Protein Analysis with Digital Microfluidics

Proteins as crucial components of cells are responsible for the majority of cellular processes. Sensitive and efficient protein detection enables a more accurate and comprehensive investigation of cellular phenotypes and life activities. Here, a protein sequencing method with high multiplexing, high throughput, high cell utilization, and integration based on digital microfluidics (DMF-Protein-seq) is proposed, which transforms protein information into DNA sequencing readout via DNA-tagged antibodies and labels single cells with unique cell barcodes. In a 184-electrode DMF-Protein-seq system, ≈1800 cells are simultaneously detected per experimental run. The digital microfluidics device harnessing low-adsorbed hydrophobic surface and contaminants-isolated reaction space supports high cell utilization (>90%) and high mapping reads (>90%) with the input cells ranging from 140 to 2000. This system leverages split&pool strategy on the DMF chip for the first time to overcome DMF platform restriction in cell analysis throughput and replace the traditionally tedious bench-top combinatorial barcoding. With the benefits of high efficiency and sensitivity in protein analysis, the system offers great potential for cell classification and drug monitoring based on protein expression at the single-cell level.

The covariance environment defines cellular niches for spatial inference

A key challenge of analyzing data from high-resolution spatial profiling technologies is to suitably represent the features of cellular neighborhoods or niches. Here we introduce the covariance environment (COVET), a representation that leverages the gene-gene covariate structure across cells in the niche to capture the multivariate nature of cellular interactions within it. We define a principled optimal transport-based distance metric between COVET niches that scales to millions of cells. Using COVET to encode spatial context, we developed environmental variational inference (ENVI), a conditional variational autoencoder that jointly embeds spatial and single-cell RNA sequencing data into a latent space. ENVI includes two decoders: one to impute gene expression across the spatial modality and a second to project spatial information onto single-cell data. ENVI can confer spatial context to genomics data from single dissociated cells and outperforms alternatives for imputing gene expression on diverse spatial datasets.

An advanced comprehensive muti-cell-type-specific model for predicting anti-PD-1 therapeutic effect in melanoma

Rationale: Immune checkpoint inhibitors targeting the programmed cell death (PD)-1/PD-L1 pathway have promise in patients with advanced melanoma. However, drug resistance usually results in limited patient benefits. Recent single-cell RNA sequencing studies have elucidated that MM patients display distinctive transcriptional features of tumor cells, immune cells and interstitial cells, including loss of antigen presentation function of tumor cells, exhaustion of CD8+T and extracellular matrix secreted by fibroblasts to prevents immune infiltration, which leads to a poor response to immune checkpoint inhibitors (ICIs). However, cell subgroups beneficial to anti-tumor immunity and the model developed by them remain to be further identified. Methods: In this clinical study of neoadjuvant therapy with anti-PD-1 in advanced melanoma, tumor tissues were collected before and after treatment for single-nucleus sequencing, and the results were verified using multicolor immunofluorescence staining and public datasets. Results: This study describes four cell subgroups which are closely associated with the effectiveness of anti-PD-1 treatment. It also describes a cell-cell communication network, in which the interaction of the four cell subgroups contributes to anti-tumor immunity. Furthermore, we discuss a newly developed predictive model based on these four subgroups that holds significant potential for assessing the efficacy of anti-PD-1 treatment. Conclusions: These findings elucidate the primary mechanism of anti-PD-1 resistance and offer guidance for clinical drug administration for melanoma.

Autoimmunity-associated allele of tyrosine phosphatase gene PTPN22 enhances anti-viral immunity

The 1858C>T allele of the tyrosine phosphatase PTPN22 is present in 5-10% of the North American population and is strongly associated with numerous autoimmune diseases. Although research has been done to define how this allele potentiates autoimmunity, the influence PTPN22 and its pro-autoimmune allele has in anti-viral immunity remains poorly defined. Here, we use single cell RNA-sequencing and functional studies to interrogate the impact of this pro-autoimmune allele on anti-viral immunity during Lymphocytic Choriomeningitis Virus clone 13 (LCMV-cl13) infection. Mice homozygous for this allele (PEP-619WW) clear the LCMV-cl13 virus whereas wildtype (PEP-WT) mice cannot. This is associated with enhanced anti-viral CD4 T cell responses and a more immunostimulatory CD8α- cDC phenotype. Adoptive transfer studies demonstrated that PEP-619WW enhanced anti-viral CD4 T cell function through virus-specific CD4 T cell intrinsic and extrinsic mechanisms. Taken together, our data show that the pro-autoimmune allele of Ptpn22 drives a beneficial anti-viral immune response thereby preventing what is normally a chronic virus infection.

Resolving cellular dynamics using single-cell temporal transcriptomics

Cellular dynamics, the transition of a cell from one state to another, is central to understanding developmental processes and disease progression. Single-cell transcriptomics has been pushing the frontiers of cellular dynamics studies into a genome-wide and single-cell level. While most single-cell RNA sequencing approaches are disruptive and only provide a snapshot of cell states, the dynamics of a cell could be reconstructed by either exploiting temporal information hiding in the transcriptomics data or integrating additional information. In this review, we describe these approaches, highlighting their underlying principles, key assumptions, and the rationality to interpret the results as models. We also discuss the recently emerging nondisruptive live-cell transcriptomics methods, which are highly complementary to the computational models for their assumption-free nature.

The Shoot Apical Meristem: An Evolutionary Molding of Higher Plants

The shoot apical meristem (SAM) gives rise to the aerial structure of plants by producing lateral organs and other meristems. The SAM is responsible for plant developmental patterns, thus determining plant morphology and, consequently, many agronomic traits such as the number and size of fruits and flowers and kernel yield. Our current understanding of SAM morphology and regulation is based on studies conducted mainly on some angiosperms, including economically important crops such as maize (Zea mays) and rice (Oryza sativa), and the model species Arabidopsis (Arabidopsis thaliana). However, studies in other plant species from the gymnosperms are scant, making difficult comparative analyses that help us understand SAM regulation in diverse plant species. This limitation prevents deciphering the mechanisms by which evolution gave rise to the multiple plant structures within the plant kingdom and determines the conserved mechanisms involved in SAM maintenance and operation. This review aims to integrate and analyze the current knowledge of SAM evolution by combining the morphological and molecular information recently reported from the plant kingdom.

A novel mass cytometry protocol optimized for immunophenotyping of low-frequency antigen-specific T cells

Understanding antigen-specific T-cell responses, for example, following virus infections or allergen exposure, is of high relevance for the development of vaccines and therapeutics. We aimed on optimizing immunophenotyping of T cells after antigen stimulation by improving staining procedures for flow and mass cytometry. Our method can be used for primary cells of both mouse and human origin for the detection of low-frequency T-cell response using a dual-barcoding system for individual samples and conditions. First, live-cell barcoding was performed using anti-CD45 antibodies prior to an in vitro T-cell stimulation assay. Second, to discriminate between stimulation conditions and prevent cell loss, sample barcoding was combined with a commercial barcoding solution. This dual-barcoding approach is cell sparing and, therefore, particularly relevant for samples with low cell numbers. To further reduce cell loss and to increase debarcoding efficiency of multiplexed samples, we combined our dual-barcoding approach with a new centrifugation-free washing system by laminar flow (Curiox™). Finally, to demonstrate the benefits of our established protocol, we assayed virus-specific T-cell response in SARS-CoV-2-vaccinated and SARS-CoV-2-infected patients and compared with healthy non-exposed individuals by a high-parameter CyTOF analysis. We could reveal a heterogeneity of phenotypes among responding CD4, CD8, and gd-T cells following antigen-specific stimulations. Our protocol allows to assay antigen-specific responses of minute populations of T cells to virus-derived peptides, allergens, or other antigens from the same donor sample, in order to investigate qualitative and quantitative differences.

Time-resolved single-cell transcriptomics defines immune trajectories in glioblastoma

Deciphering the cell-state transitions underlying immune adaptation across time is fundamental for advancing biology. Empirical in vivo genomic technologies that capture cellular dynamics are currently lacking. We present Zman-seq, a single-cell technology recording transcriptomic dynamics across time by introducing time stamps into circulating immune cells, tracking them in tissues for days. Applying Zman-seq resolved cell-state and molecular trajectories of the dysfunctional immune microenvironment in glioblastoma. Within 24 hours of tumor infiltration, cytotoxic natural killer cells transitioned to a dysfunctional program regulated by TGFB1 signaling. Infiltrating monocytes differentiated into immunosuppressive macrophages, characterized by the upregulation of suppressive myeloid checkpoints Trem2, Il18bp, and Arg1, over 36 to 48 hours. Treatment with an antagonistic anti-TREM2 antibody reshaped the tumor microenvironment by redirecting the monocyte trajectory toward pro-inflammatory macrophages. Zman-seq is a broadly applicable technology, enabling empirical measurements of differentiation trajectories, which can enhance the development of more efficacious immunotherapies.

Stem cell-like reprogramming is required for leukemia-initiating activity in B-ALL

B cell acute lymphoblastic leukemia (B-ALL) is a multistep disease characterized by the hierarchical acquisition of genetic alterations. However, the question of how a primary oncogene reprograms stem cell-like properties in committed B cells and leads to a preneoplastic population remains unclear. Here, we used the PAX5::ELN oncogenic model to demonstrate a causal link between the differentiation blockade, the self-renewal, and the emergence of preleukemic stem cells (pre-LSCs). We show that PAX5::ELN disrupts the differentiation of preleukemic cells by enforcing the IL7r/JAK-STAT pathway. This disruption is associated with the induction of rare and quiescent pre-LSCs that sustain the leukemia-initiating activity, as assessed using the H2B-GFP model. Integration of transcriptomic and chromatin accessibility data reveals that those quiescent pre-LSCs lose B cell identity and reactivate an immature molecular program, reminiscent of human B-ALL chemo-resistant cells. Finally, our transcriptional regulatory network reveals the transcription factor EGR1 as a strong candidate to control quiescence/resistance of PAX5::ELN pre-LSCs as well as of blasts from human B-ALL.

A relay velocity model infers cell-dependent RNA velocity

RNA velocity provides an approach for inferring cellular state transitions from single-cell RNA sequencing (scRNA-seq) data. Conventional RNA velocity models infer universal kinetics from all cells in an scRNA-seq experiment, resulting in unpredictable performance in experiments with multi-stage and/or multi-lineage transition of cell states where the assumption of the same kinetic rates for all cells no longer holds. Here we present cellDancer, a scalable deep neural network that locally infers velocity for each cell from its neighbors and then relays a series of local velocities to provide single-cell resolution inference of velocity kinetics. In the simulation benchmark, cellDancer shows robust performance in multiple kinetic regimes, high dropout ratio datasets and sparse datasets. We show that cellDancer overcomes the limitations of existing RNA velocity models in modeling erythroid maturation and hippocampus development. Moreover, cellDancer provides cell-specific predictions of transcription, splicing and degradation rates, which we identify as potential indicators of cell fate in the mouse pancreas.

Single-cell expression predicts neuron-specific protein homeostasis networks

The protein homeostasis network keeps proteins in their correct shapes and avoids unwanted aggregation. In turn, the accumulation of aberrantly misfolded proteins has been directly associated with the onset of ageing-associated neurodegenerative diseases such as Alzheimer's and Parkinson's. However, a detailed and rational understanding of how protein homeostasis is achieved in health, and how it can be targeted for therapeutic intervention in diseases remains missing. Here, large-scale single-cell expression data from the Allen Brain Map are analysed to investigate the transcription regulation of the core protein homeostasis network across the human brain. Remarkably, distinct expression profiles suggest specialized protein homeostasis networks with systematic adaptations in excitatory neurons, inhibitory neurons and non-neuronal cells. Moreover, several chaperones and Ubiquitin ligases are found transcriptionally coregulated with genes important for synapse formation and maintenance, thus linking protein homeostasis to the regulation of neuronal function. Finally, evolutionary analyses highlight the conservation of an elevated interaction density in the chaperone network, suggesting that one of the most exciting aspects of chaperone action may yet be discovered in their collective action at the systems level. More generally, our work highlights the power of computational analyses for breaking down complexity and gaining complementary insights into fundamental biological problems.

Quantitative Modeling of Stemness in Single-Cell RNA Sequencing Data: A Nonlinear One-Class Support Vector Machine Method

Intratumoral heterogeneity and the presence of cancer stem cells are challenging issues in cancer therapy. An appropriate quantification of the stemness of individual cells for assessing the potential for self-renewal and differentiation from the cell of origin can define a measurement for quantifying different cell states, which is important in understanding the dynamics of cancer evolution, and might further provide possible targeted therapies aimed at tumor stem cells. Nevertheless, it is usually difficult to quantify the stemness of a cell based on molecular information associated with the cell. In this study, we proposed a stemness definition method with one-class Hadamard kernel support vector machine (OCHSVM) based on single-cell RNA sequencing (scRNA-seq) data. Applications of the proposed OCHSVM stemness are assessed by various data sets, including preimplantation embryo cells, induced pluripotent stem cells, or tumor cells. We further compared the OCHSVM model with state-of-the-art methods CytoTRACE, one-class logistic regression, or one-class SVM methods with different kernels. The computational results demonstrate that the OCHSVM method is more suitable for stemness identification using scRNA-seq data.

Gene module reconstruction elucidates cellular differentiation processes and the regulatory logic of specialized secretion

During differentiation, cells become structurally and functionally specialized, but comprehensive views of the underlying remodeling processes are elusive. Here, we leverage scRNA-seq developmental trajectories to reconstruct differentiation using two secretory tissues as a model system - the zebrafish notochord and hatching gland. First, we present an approach to integrate expression and functional similarities for gene module identification, revealing dozens of gene modules representing known and newly associated differentiation processes and their temporal ordering. Second, we focused on the unfolded protein response (UPR) transducer module to study how general versus cell-type specific secretory functions are regulated. By profiling loss- and gain-of-function embryos, we found that the UPR transcription factors creb3l1, creb3l2, and xbp1 are master regulators of a general secretion program. creb3l1/creb3l2 additionally activate an extracellular matrix secretion program, while xbp1 partners with bhlha15 to activate a gland-specific secretion program. Our study offers a multi-source integrated approach for functional gene module identification and illustrates how transcription factors confer general and specialized cellular functions.

Single cell proteomics characterization of bone marrow hematopoiesis with distinct Ras pathway lesions

Normal hematopoiesis requires constant prolific production of different blood cell lineages by multipotent hematopoietic stem cells (HSC). Stem- and progenitor- cells need to balance dormancy with proliferation. How genetic alterations impact frequency, lineage potential, and metabolism of HSC is largely unknown. Here, we compared induced expression of KRAS G12D or RasGRP1 to normal hematopoiesis. At low-resolution, both Ras pathway lesions result in skewing towards myeloid lineages. Single-cell resolution CyTOF proteomics unmasked an expansion of HSC- and progenitor- compartments for RasGRP1, contrasted by a depletion for KRAS G12D . SCENITH™ quantitates protein synthesis with single-cell precision and corroborated that immature cells display low metabolic SCENITH™ rates. Both RasGRP1 and KRAS G12D elevated mean SCENITH™ signals in immature cells. However, RasGRP1-overexpressing stem cells retain a metabolically quiescent cell-fraction, whereas this fraction diminishes for KRAS G12D . Our temporal single cell proteomics and metabolomics datasets provide a resource of mechanistic insights into altered hematopoiesis at single cell resolution.

SEACells infers transcriptional and epigenomic cellular states from single-cell genomics data

Metacells are cell groupings derived from single-cell sequencing data that represent highly granular, distinct cell states. Here we present single-cell aggregation of cell states (SEACells), an algorithm for identifying metacells that overcome the sparsity of single-cell data while retaining heterogeneity obscured by traditional cell clustering. SEACells outperforms existing algorithms in identifying comprehensive, compact and well-separated metacells in both RNA and assay for transposase-accessible chromatin (ATAC) modalities across datasets with discrete cell types and continuous trajectories. We demonstrate the use of SEACells to improve gene-peak associations, compute ATAC gene scores and infer the activities of critical regulators during differentiation. Metacell-level analysis scales to large datasets and is particularly well suited for patient cohorts, where per-patient aggregation provides more robust units for data integration. We use our metacells to reveal expression dynamics and gradual reconfiguration of the chromatin landscape during hematopoietic differentiation and to uniquely identify CD4 T cell differentiation and activation states associated with disease onset and severity in a Coronavirus Disease 2019 (COVID-19) patient cohort.

Harnessing CD3 diversity to optimize CAR T cells

Current US Food and Drug Administration-approved chimeric antigen receptor (CAR) T cells harbor the T cell receptor (TCR)-derived ζ chain as an intracellular activation domain in addition to costimulatory domains. The functionality in a CAR format of the other chains of the TCR complex, namely CD3δ, CD3ε and CD3γ, instead of ζ, remains unknown. In the present study, we have systematically engineered new CD3 CARs, each containing only one of the CD3 intracellular domains. We found that CARs containing CD3δ, CD3ε or CD3γ cytoplasmic tails outperformed the conventional ζ CAR T cells in vivo. Transcriptomic and proteomic analysis revealed differences in activation potential, metabolism and stimulation-induced T cell dysfunctionality that mechanistically explain the enhanced anti-tumor performance. Furthermore, dimerization of the CARs improved their overall functionality. Using these CARs as minimalistic and synthetic surrogate TCRs, we have identified the phosphatase SHP-1 as a new interaction partner of CD3δ that binds the CD3δ-ITAM on phosphorylation of its C-terminal tyrosine. SHP-1 attenuates and restrains activation signals and might thus prevent exhaustion and dysfunction. These new insights into T cell activation could promote the rational redesign of synthetic antigen receptors to improve cancer immunotherapy.

New perspectives on biology, disease progression, and therapy response of head and neck cancer gained from single cell RNA sequencing and spatial transcriptomics

Head and neck squamous cell carcinoma (HNSCC) is one of the most frequent cancers worldwide. The main risk factors are consumption of tobacco products and alcohol, as well as infection with human papilloma virus. Approved therapeutic options comprise surgery, radiation, chemotherapy, targeted therapy through epidermal growth factor receptor inhibition, and immunotherapy, but outcome has remained unsatisfactory due to recurrence rates of ~50% and the frequent occurrence of second primaries. The availability of the human genome sequence at the beginning of the millennium heralded the omics era, in which rapid technological progress has advanced our knowledge of the molecular biology of malignant diseases, including HNSCC, at an unprecedented pace. Initially, microarray-based methods, followed by approaches based on next-generation sequencing, were applied to study the genetics, epigenetics, and gene expression patterns of bulk tumors. More recently, the advent of single-cell RNA sequencing (scRNAseq) and spatial transcriptomics methods has facilitated the investigation of the heterogeneity between and within different cell populations in the tumor microenvironment (e.g., cancer cells, fibroblasts, immune cells, endothelial cells), led to the discovery of novel cell types, and advanced the discovery of cell-cell communication within tumors. This review provides an overview of scRNAseq, spatial transcriptomics, and the associated bioinformatics methods, and summarizes how their application has promoted our understanding of the emergence, composition, progression, and therapy responsiveness of, and intercellular signaling within, HNSCC.

A statistical framework for differential pseudotime analysis with multiple single-cell RNA-seq samples

Pseudotime analysis with single-cell RNA-sequencing (scRNA-seq) data has been widely used to study dynamic gene regulatory programs along continuous biological processes. While many methods have been developed to infer the pseudotemporal trajectories of cells within a biological sample, it remains a challenge to compare pseudotemporal patterns with multiple samples (or replicates) across different experimental conditions. Here, we introduce Lamian, a comprehensive and statistically-rigorous computational framework for differential multi-sample pseudotime analysis. Lamian can be used to identify changes in a biological process associated with sample covariates, such as different biological conditions while adjusting for batch effects, and to detect changes in gene expression, cell density, and topology of a pseudotemporal trajectory. Unlike existing methods that ignore sample variability, Lamian draws statistical inference after accounting for cross-sample variability and hence substantially reduces sample-specific false discoveries that are not generalizable to new samples. Using both real scRNA-seq and simulation data, including an analysis of differential immune response programs between COVID-19 patients with different disease severity levels, we demonstrate the advantages of Lamian in decoding cellular gene expression programs in continuous biological processes.

Mass cytometry as a tool in target validation and drug discovery

Mass cytometry provides highly multiparametric data at a single cell level, coupling the specificity and sensitivity of time-of-flight mass spectrometry with the single-cell throughput of flow cytometry. It offers great value in interrogating the potentially heterogenous impact that a drug may have on a biological system, allowing an investigator to capture not just changes in cell behavior, but how these changes may differ between cell subtypes. In this chapter, we review the technical details of the platform as well as its limitations, before describing our approach to planning and running a mass cytometry experiment. A series of method modules, spanning the staining process through to data cleaning, are described that are then combined to create three separate experiments. The first experiment illustrates a core process in mass cytometry: the validation and titration of a metal-conjugated antibody reporter. The second experiment explores the impact of a kinase inhibitor on cell cycle and apoptosis pathways of a human myeloma cell line. And the third experiment exploits the multiparametric capability of mass cytometry, by exploring the differential expression changes in a transcription factor upon drug treatment across the cellular compartments of a peripheral blood mononuclear cell sample.

Generation and molecular characterization of human pluripotent stem cell-derived pharyngeal foregut endoderm

Approaches to study human pharyngeal foregut endoderm-a developmental intermediate that is linked to various human syndromes involving pharynx development and organogenesis of tissues such as thymus, parathyroid, and thyroid-have been hampered by scarcity of tissue access and cellular models. We present an efficient stepwise differentiation method to generate human pharyngeal foregut endoderm from pluripotent stem cells. We determine dose and temporal requirements of signaling pathway engagement for optimized differentiation and characterize the differentiation products on cellular and integrated molecular level. We present a computational classification tool, "CellMatch," and transcriptomic classification of differentiation products on an integrated mouse scRNA-seq developmental roadmap confirms cellular maturation. Integrated transcriptomic and chromatin analyses infer differentiation stage-specific gene regulatory networks. Our work provides the method and integrated multiomic resource for the investigation of disease-relevant loci and gene regulatory networks and their role in developmental defects affecting the pharyngeal endoderm and its derivatives.

High-dimensional analysis reveals an immune atlas and novel neutrophil clusters in the lungs of model animals with Actinobacillus pleuropneumoniae-induced pneumonia

Due to the increase in bacterial resistance, improving the anti-infectious immunity of the host is rapidly becoming a new strategy for the prevention and treatment of bacterial pneumonia. However, the specific lung immune responses and key immune cell subsets involved in bacterial infection are obscure. Actinobacillus pleuropneumoniae (APP) can cause porcine pleuropneumonia, a highly contagious respiratory disease that has caused severe economic losses in the swine industry. Here, using high-dimensional mass cytometry, the major immune cell repertoire in the lungs of mice with APP infection was profiled. Various phenotypically distinct neutrophil subsets and Ly-6C+ inflammatory monocytes/macrophages accumulated post-infection. Moreover, a linear differentiation trajectory from inactivated to activated to apoptotic neutrophils corresponded with the stages of uninfected, onset, and recovery of APP infection. CD14+ neutrophils, which mainly increased in number during the recovery stage of infection, were revealed to have a stronger ability to produce cytokines, especially IL-10 and IL-21, than their CD14- counterparts. Importantly, MHC-II+ neutrophils with antigen-presenting cell features were identified, and their numbers increased in the lung after APP infection. Similar results were further confirmed in the lungs of piglets infected with APP and Klebsiella pneumoniae infection by using a single-cell RNA-seq technique. Additionally, a correlation analysis between cluster composition and the infection process yielded a dynamic and temporally associated immune landscape where key immune clusters, including previously unrecognized ones, marked various stages of infection. Thus, these results reveal the characteristics of key neutrophil clusters and provide a detailed understanding of the immune response to bacterial pneumonia.

GCSTI: A Single-Cell Pseudotemporal Trajectory Inference Method Based on Graph Compression

The single-cell pseudotemporal trajectory inference is an important way to explore the process of developmental changes within a cell. Due to the uneven rate of cell growth, changes in gene expression depend less on the time of data collection and more on a cell's "internal clock". To overcome the challenges of gene analysis, and replicate biological developmental processes, several strategies have been put forth. However, due to the size of single-cell datasets, locating relevant signposts usually necessitate clustering analysis or a sizable amount of priori information. To this end, we propose a novel single-cell pseudotemporal trajectory inference technique: GCSTI method, which is based on graph compression and doesn't rely on a priori knowledge or clustering procedures, can handle the trajectory inference problem for a large network in a stable and efficient manner. Additionally, we simultaneously improve the pseudotime defining method currently employed in this study in order to obtain more trustworthy and beneficial outcomes for trajectory inference. Finally, we validate the efficacy and stability of the GCSTI method using datasets from human skeletal muscle myogenic cells and four simulated datasets.

Dissecting and improving gene regulatory network inference using single-cell transcriptome data

Single-cell transcriptome data has been widely used to reconstruct gene regulatory networks (GRNs) controlling critical biological processes such as development and differentiation. Although a growing list of algorithms has been developed to infer GRNs using such data, achieving an inference accuracy consistently higher than random guessing has remained challenging. To address this, it is essential to delineate how the accuracy of regulatory inference is limited. Here, we systematically characterized factors limiting the accuracy of inferred GRNs and demonstrated that using pre-mRNA information can help improve regulatory inference compared to the typically used information (i.e., mature mRNA). Using kinetic modeling and simulated single-cell data sets, we showed that target genes' mature mRNA levels often fail to accurately report upstream regulatory activities because of gene-level and network-level factors, which can be improved by using pre-mRNA levels. We tested this finding on public single-cell RNA-seq data sets using intronic reads as proxies of pre-mRNA levels and can indeed achieve a higher inference accuracy compared to using exonic reads (corresponding to mature mRNAs). Using experimental data sets, we further validated findings from the simulated data sets and identified factors such as transcription factor activity dynamics influencing the accuracy of pre-mRNA-based inference. This work delineates the fundamental limitations of gene regulatory inference and helps improve GRN inference using single-cell RNA-seq data.

Evolving Approach to Clinical Cytometry for Immunodeficiencies and Other Immune Disorders

Primary immunodeficiencies were initially identified on the basis of recurrent, severe or unusual infections. Subsequently, it was noted that these diseases can also manifest with autoimmunity, autoinflammation, allergy, lymphoproliferation and malignancy, hence a conceptual change and their renaming as inborn errors of immunity. Ongoing advances in flow cytometry provide the opportunity to expand or modify the utility and scope of existing laboratory tests in this field to mirror this conceptual change. Here we have used the B cell subset, variably known as CD21low B cells, age-associated B cells and T-bet+ B cells, as an example to demonstrate this possibility.

Atonal homolog 7 (ATOH7) confers neuroprotection for photoreceptor cells in glaucoma via inhibition of the notch pathway

Single-cell RNA sequencing (scRNA-seq) technologies enable the profiling and analysis of the transcriptomes of single cells and hold promise for clarifying gene mechanisms at single-cell resolution. We based this study on scRNA-seq data to reveal glaucoma-related genes and downstream pathways with neuroprotection effects. The scRNA-seq datasets related to glaucoma of retinal tissue samples of human beings and Atonal Homolog 7 (ATOH7)-null mice were obtained from the GEO database. The 74 top marker genes and 20 cell clusters were obtained in human retinal tissue samples. The key gene ATOH7 was found after the intersection with genes from GeneCards data. In the ATOH7-null mouse retinal tissue samples, pseudotime inference demonstrated significant changes in cell differentiation. Moreover, mouse retinal photoreceptor cells (PRCs) were cultured and treated with lentivirus carrying oe-ATOH7 alone or in combination with Notch signaling pathway activator Jagged-1/FC, after which cell biological functions were determined. The involvement of ATOH7 in glaucoma was identified through regulating PRCs. Furthermore, ATOH7 conferred neuroprotection in PRCs in glaucoma by mediating the Notch signaling pathway. In vitro data confirmed that ATOH7 overexpression promoted the differentiation of PRCs and inhibited their apoptosis by suppressing the Notch signaling pathway. The evidence provided by our study highlighted the involvement of ATOH7 in the blockade of the Notch signaling pathway, resulting in the neuroprotection for PRCs in glaucoma.

Unravelling human hematopoietic progenitor cell diversity through association with intrinsic regulatory factors

Hematopoietic stem and progenitor cell (HSPC) transplantation is an essential therapy for hematological conditions, but finer definitions of human HSPC subsets with associated function could enable better tuning of grafts and more routine, lower-risk application. To deeply phenotype HSPCs, following a screen of 328 antigens, we quantified 41 surface proteins and functional regulators on millions of CD34+ and CD34- cells, spanning four primary human hematopoietic tissues: bone marrow, mobilized peripheral blood, cord blood, and fetal liver. We propose more granular definitions of HSPC subsets and provide new, detailed differentiation trajectories of erythroid and myeloid lineages. These aspects of our revised human hematopoietic model were validated with corresponding epigenetic analysis and in vitro clonal differentiation assays. Overall, we demonstrate the utility of using molecular regulators as surrogates for cellular identity and functional potential, providing a framework for description, prospective isolation, and cross-tissue comparison of HSPCs in humans.

Progressive plasticity during colorectal cancer metastasis

Metastasis is the principal cause of cancer death, yet we lack an understanding of metastatic cell states, their relationship to primary tumor states, and the mechanisms by which they transition. In a cohort of biospecimen trios from same-patient normal colon, primary and metastatic colorectal cancer, we show that while primary tumors largely adopt LGR5 + intestinal stem-like states, metastases display progressive plasticity. Loss of intestinal cell states is accompanied by reprogramming into a highly conserved fetal progenitor state, followed by non-canonical differentiation into divergent squamous and neuroendocrine-like states, which is exacerbated by chemotherapy and associated with poor patient survival. Using matched patient-derived organoids, we demonstrate that metastatic cancer cells exhibit greater cell-autonomous multilineage differentiation potential in response to microenvironment cues than their intestinal lineage-restricted primary tumor counterparts. We identify PROX1 as a stabilizer of intestinal lineage in the fetal progenitor state, whose downregulation licenses non-canonical reprogramming.