Immunology in the Era of Single-Cell Technologies.

Immune cells are characterized by diversity, specificity, plasticity, and adaptability-properties that enable them to contribute to homeostasis and respond specifically and dynamically to the many threats encountered by the body. Single-cell technologies, including the assessment of transcriptomics, genomics, and proteomics at the level of individual cells, are ideally suited to studying these properties of immune cells. In this review we discuss the benefits of adopting single-cell approaches in studying underappreciated qualities of immune cells and highlight examples where these technologies have been critical to advancing our understanding of the immune system in health and disease.

A single-embryo, single-cell time-resolved model for mouse gastrulation.

Mouse embryonic development is a canonical model system for studying mammalian cell fate acquisition. Recently, single-cell atlases comprehensively charted embryonic transcriptional landscapes, yet inference of the coordinated dynamics of cells over such atlases remains challenging. Here, we introduce a temporal model for mouse gastrulation, consisting of data from 153 individually sampled embryos spanning 36 h of molecular diversification. Using algorithms and precise timing, we infer differentiation flows and lineage specification dynamics over the embryonic transcriptional manifold. Rapid transcriptional bifurcations characterize the commitment of early specialized node and blood cells. However, for most lineages, we observe combinatorial multi-furcation dynamics rather than hierarchical transcriptional transitions. In the mesoderm, dozens of transcription factors combinatorially regulate multifurcations, as we exemplify using time-matched chimeric embryos of Foxc1/Foxc2 mutants. Our study rejects the notion of differentiation being governed by a series of binary choices, providing an alternative quantitative model for cell fate acquisition.

Longitudinal Multi-omics Analyses Identify Responses of Megakaryocytes, Erythroid Cells, and Plasmablasts as Hallmarks of Severe COVID-19.

Temporal resolution of cellular features associated with a severe COVID-19 disease trajectory is needed for understanding skewed immune responses and defining predictors of outcome. Here, we performed a longitudinal multi-omics study using a two-center cohort of 14 patients. We analyzed the bulk transcriptome, bulk DNA methylome, and single-cell transcriptome (>358,000 cells, including BCR profiles) of peripheral blood samples harvested from up to 5 time points. Validation was performed in two independent cohorts of COVID-19 patients. Severe COVID-19 was characterized by an increase of proliferating, metabolically hyperactive plasmablasts. Coinciding with critical illness, we also identified an expansion of interferon-activated circulating megakaryocytes and increased erythropoiesis with features of hypoxic signaling. Megakaryocyte- and erythroid-cell-derived co-expression modules were predictive of fatal disease outcome. The study demonstrates broad cellular effects of SARS-CoV-2 infection beyond adaptive immune cells and provides an entry point toward developing biomarkers and targeted treatments of patients with COVID-19.

The Biogenesis of SRP RNA Is Modulated by an RNA Folding Intermediate Attained during Transcription.

The signal recognition particle (SRP), responsible for co-translational protein targeting and delivery to cellular membranes, depends on the native long-hairpin fold of its RNA to confer functionality. Since RNA initiates folding during its synthesis, we used high-resolution optical tweezers to follow in real time the co-transcriptional folding of SRP RNA. Surprisingly, SRP RNA folding is robust to transcription rate changes and the presence or absence of its 5'-precursor sequence. The folding pathway also reveals the obligatory attainment of a non-native hairpin intermediate (H1) that eventually rearranges into the native fold. Furthermore, H1 provides a structural platform alternative to the native fold for RNase P to bind and mature SRP RNA co-transcriptionally. Delays in attaining the final native fold are detrimental to the cell, altogether showing that a co-transcriptional folding pathway underpins the proper biogenesis of function-essential SRP RNA.

Changes in above- versus belowground biomass distribution in permafrost regions in response to climate warming.

Permafrost regions contain approximately half of the carbon stored in land ecosystems and have warmed at least twice as much as any other biome. This warming has influenced vegetation activity, leading to changes in plant composition, physiology, and biomass storage in aboveground and belowground components, ultimately impacting ecosystem carbon balance. Yet, little is known about the causes and magnitude of long-term changes in the above- to belowground biomass ratio of plants (η). Here, we analyzed η values using 3,013 plots and 26,337 species-specific measurements across eight sites on the Tibetan Plateau from 1995 to 2021. Our analysis revealed distinct temporal trends in η for three vegetation types: a 17% increase in alpine wetlands, and a decrease of 26% and 48% in alpine meadows and alpine steppes, respectively. These trends were primarily driven by temperature-induced growth preferences rather than shifts in plant species composition. Our findings indicate that in wetter ecosystems, climate warming promotes aboveground plant growth, while in drier ecosystems, such as alpine meadows and alpine steppes, plants allocate more biomass belowground. Furthermore, we observed a threefold strengthening of the warming effect on η over the past 27 y. Soil moisture was found to modulate the sensitivity of η to soil temperature in alpine meadows and alpine steppes, but not in alpine wetlands. Our results contribute to a better understanding of the processes driving the response of biomass distribution to climate warming, which is crucial for predicting the future carbon trajectory of permafrost ecosystems and climate feedback.

Age-related behavioral resilience in smartphone touchscreen interaction dynamics.

We experience a life that is full of ups and downs. The ability to bounce back after adverse life events such as the loss of a loved one or serious illness declines with age, and such isolated events can even trigger accelerated aging. How humans respond to common day-to-day perturbations is less clear. Here, we infer the aging status from smartphone behavior by using a decision tree regression model trained to accurately estimate the chronological age based on the dynamics of touchscreen interactions. Individuals (N = 280, 21 to 87 y of age) expressed smartphone behavior that appeared younger on certain days and older on other days through the observation period that lasted up to ~4 y. We captured the essence of these fluctuations by leveraging the mathematical concept of critical transitions and tipping points in complex systems. In most individuals, we find one or more alternative stable aging states separated by tipping points. The older the individual, the lower the resilience to forces that push the behavior across the tipping point into an older state. Traditional accounts of aging based on sparse longitudinal data spanning decades suggest a gradual behavioral decline with age. Taken together with our current results, we propose that the gradual age-related changes are interleaved with more complex dynamics at shorter timescales where the same individual may navigate distinct behavioral aging states from one day to the next. Real-world behavioral data modeled as a complex system can transform how we view and study aging.

Dissection and integration of bursty transcriptional dynamics for complex systems.

RNA velocity estimation is a potentially powerful tool to reveal the directionality of transcriptional changes in single-cell RNA-sequencing data, but it lacks accuracy, absent advanced metabolic labeling techniques. We developed an approach, TopicVelo, that disentangles simultaneous, yet distinct, dynamics by using a probabilistic topic model, a highly interpretable form of latent space factorization, to infer cells and genes associated with individual processes, thereby capturing cellular pluripotency or multifaceted functionality. Focusing on process-associated cells and genes enables accurate estimation of process-specific velocities via a master equation for a transcriptional burst model accounting for intrinsic stochasticity. The method obtains a global transition matrix by leveraging cell topic weights to integrate process-specific signals. In challenging systems, this method accurately recovers complex transitions and terminal states, while our use of first-passage time analysis provides insights into transient transitions. These results expand the limits of RNA velocity, empowering future studies of cell fate and functional responses.

scTOP: physics-inspired order parameters for cellular identification and visualization.

Advances in single-cell RNA sequencing provide an unprecedented window into cellular identity. The abundance of data requires new theoretical and computational frameworks to analyze the dynamics of differentiation and integrate knowledge from cell atlases. We present 'single-cell Type Order Parameters' (scTOP): a statistical, physics-inspired approach for quantifying cell identity given a reference basis of cell types. scTOP can accurately classify cells, visualize developmental trajectories and assess the fidelity of engineered cells. Importantly, scTOP does this without feature selection, statistical fitting or dimensional reduction (e.g. uniform manifold approximation and projection, principle components analysis, etc.). We illustrate the power of scTOP using human and mouse datasets. By reanalyzing mouse lung data, we characterize a transient hybrid alveolar type 1/alveolar type 2 cell population. Visualizations of lineage tracing hematopoiesis data using scTOP confirm that a single clone can give rise to multiple mature cell types. We assess the transcriptional similarity between endogenous and donor-derived cells in the context of murine pulmonary cell transplantation. Our results suggest that physics-inspired order parameters can be an important tool for understanding differentiation and characterizing engineered cells. scTOP is available as an easy-to-use Python package.

Defining developmental trajectories of prosensory cells in human inner ear organoids at single-cell resolution.

The inner ear sensory epithelia contain mechanosensitive hair cells and supporting cells. Both cell types arise from SOX2-expressing prosensory cells, but the mechanisms underlying the diversification of these cell lineages remain unclear. To determine the transcriptional trajectory of prosensory cells, we established a SOX2-2A-ntdTomato human embryonic stem cell line using CRISPR/Cas9, and performed single-cell RNA-sequencing analyses with SOX2-positive cells isolated from inner ear organoids at various time points between differentiation days 20 and 60. Our pseudotime analysis suggests that vestibular type II hair cells arise primarily from supporting cells, rather than bi-fated prosensory cells in organoids. Moreover, ion channel- and ion-transporter-related gene sets were enriched in supporting cells versus prosensory cells, whereas Wnt signaling-related gene sets were enriched in hair cells versus supporting cells. These findings provide valuable insights into how prosensory cells give rise to hair cells and supporting cells during human inner ear development, and may provide a clue to promote hair cell regeneration from resident supporting cells in individuals with hearing loss or balance disorders.

Power in Numbers: Single-Cell RNA-Seq Strategies to Dissect Complex Tissues.

The growing scale and declining cost of single-cell RNA-sequencing (RNA-seq) now permit a repetition of cell sampling that increases the power to detect rare cell states, reconstruct developmental trajectories, and measure phenotype in new terms such as cellular variance. The characterization of anatomy and developmental dynamics has not had an equivalent breakthrough since groundbreaking advances in live fluorescent microscopy. The new resolution obtained by single-cell RNA-seq is a boon to genetics because the novel description of phenotype offers the opportunity to refine gene function and dissect pleiotropy. In addition, the recent pairing of high-throughput genetic perturbation with single-cell RNA-seq has made practical a scale of genetic screening not previously possible.

Characterizing efficient feature selection for single-cell expression analysis.

Unsupervised feature selection is a critical step for efficient and accurate analysis of single-cell RNA-seq data. Previous benchmarks used two different criteria to compare feature selection methods: (i) proportion of ground-truth marker genes included in the selected features and (ii) accuracy of cell clustering using ground-truth cell types. Here, we systematically compare the performance of 11 feature selection methods for both criteria. We first demonstrate the discordance between these criteria and suggest using the latter. We then compare the distribution of selected genes in their means between feature selection methods. We show that lowly expressed genes exhibit seriously high coefficients of variation and are mostly excluded by high-performance methods. In particular, high-deviation- and high-expression-based methods outperform the widely used in Seurat package in clustering cells and data visualization. We further show they also enable a clear separation of the same cell type from different tissues as well as accurate estimation of cell trajectories.

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.

scCRT: a contrastive-based dimensionality reduction model for scRNA-seq trajectory inference.

Trajectory inference is a crucial task in single-cell RNA-sequencing downstream analysis, which can reveal the dynamic processes of biological development, including cell differentiation. Dimensionality reduction is an important step in the trajectory inference process. However, most existing trajectory methods rely on cell features derived from traditional dimensionality reduction methods, such as principal component analysis and uniform manifold approximation and projection. These methods are not specifically designed for trajectory inference and fail to fully leverage prior information from upstream analysis, limiting their performance. Here, we introduce scCRT, a novel dimensionality reduction model for trajectory inference. In order to utilize prior information to learn accurate cells representation, scCRT integrates two feature learning components: a cell-level pairwise module and a cluster-level contrastive module. The cell-level module focuses on learning accurate cell representations in a reduced-dimensionality space while maintaining the cell-cell positional relationships in the original space. The cluster-level contrastive module uses prior cell state information to aggregate similar cells, preventing excessive dispersion in the low-dimensional space. Experimental findings from 54 real and 81 synthetic datasets, totaling 135 datasets, highlighted the superior performance of scCRT compared with commonly used trajectory inference methods. Additionally, an ablation study revealed that both cell-level and cluster-level modules enhance the model's ability to learn accurate cell features, facilitating cell lineage inference. The source code of scCRT is available at https://github.com/yuchen21-web/scCRT-for-scRNA-seq.

Longitudinal mapping of the development of cortical thickness and surface area in rhesus macaques during the first three years.

Studying dynamic spatiotemporal patterns of early brain development in macaque monkeys is critical for understanding the cortical organization and evolution in humans, given the phylogenetic closeness between humans and macaques. However, due to huge challenges in the analysis of early brain Magnetic Resonance Imaging (MRI) data typically with extremely low contrast and dynamic imaging appearances, our knowledge of the early macaque cortical development remains scarce. To fill this critical gap, this paper characterizes the early developmental patterns of cortical thickness and surface area in rhesus macaques by leveraging advanced computing tools tailored for early developing brains based on a densely sampled longitudinal dataset with 140 rhesus macaque MRI scans seamlessly covering from birth to 36 mo of age. The average cortical thickness exhibits an inverted U-shaped trajectory with peak thickness at around 4.3 mo of age, which is remarkably in line with the age of peak thickness at 14 mo in humans, considering the around 3:1 age ratio of human to macaque. The total cortical surface area in macaques increases monotonically but with relatively lower expansions than in humans. The spatial distributions of thicker and thinner regions are quite consistent during development, with gyri having a thicker cortex than sulci. By 4 mo of age, over 81% of cortical vertices have reached their peaks in thickness, except for the insula and medial temporal cortices, while most cortical vertices keep expanding in surface area, except for the occipital cortex. These findings provide important insights into early brain development and evolution in primates.

Mapping human brain charts cross-sectionally and longitudinally.

Brain scans acquired across large, age-diverse cohorts have facilitated recent progress in establishing normative brain aging charts. Here, we ask the critical question of whether cross-sectional estimates of age-related brain trajectories resemble those directly measured from longitudinal data. We show that age-related brain changes inferred from cross-sectionally mapped brain charts can substantially underestimate actual changes measured longitudinally. We further find that brain aging trajectories vary markedly between individuals and are difficult to predict with population-level age trends estimated cross-sectionally. Prediction errors relate modestly to neuroimaging confounds and lifestyle factors. Our findings provide explicit evidence for the importance of longitudinal measurements in ascertaining brain development and aging trajectories.

Single-nuclei transcriptome analysis of the shoot apex vascular system differentiation in Populus.

Differentiation of stem cells in the plant apex gives rise to aerial tissues and organs. Presently, we lack a lineage map of the shoot apex cells in woody perennials - a crucial gap considering their role in determining primary and secondary growth. Here, we used single-nuclei RNA-sequencing to determine cell type-specific transcriptomes of the Populus vegetative shoot apex. We identified highly heterogeneous cell populations clustered into seven broad groups represented by 18 transcriptionally distinct cell clusters. Next, we established the developmental trajectories of the epidermis, leaf mesophyll and vascular tissue. Motivated by the high similarities between Populus and Arabidopsis cell population in the vegetative apex, we applied a pipeline for interspecific single-cell gene expression data integration. We contrasted the developmental trajectories of primary phloem and xylem formation in both species, establishing the first comparison of vascular development between a model annual herbaceous and a woody perennial plant species. Our results offer a valuable resource for investigating the principles underlying cell division and differentiation conserved between herbaceous and perennial species while also allowing us to examine species-specific differences at single-cell resolution.

Tracing human life trajectory using gut microbial communities by context-aware deep learning.

The gut microbial communities are highly plastic throughout life, and the human gut microbial communities show spatial-temporal dynamic patterns at different life stages. However, the underlying association between gut microbial communities and time-related factors remains unclear. The lack of context-awareness, insufficient data, and the existence of batch effect are the three major issues, making the life trajection of the host based on gut microbial communities problematic. Here, we used a novel computational approach (microDELTA, microbial-based deep life trajectory) to track longitudinal human gut microbial communities' alterations, which employs transfer learning for context-aware mining of gut microbial community dynamics at different life stages. Using an infant cohort, we demonstrated that microDELTA outperformed Neural Network for accurately predicting the age of infant with different delivery mode, especially for newborn infants of vaginal delivery with the area under the receiver operating characteristic curve of microDELTA and Neural Network at 0.811 and 0.436, respectively. In this context, we have discovered the influence of delivery mode on infant gut microbial communities. Along the human lifespan, we also applied microDELTA to a Chinese traveler cohort, a Hadza hunter-gatherer cohort and an elderly cohort. Results revealed the association between long-term dietary shifts during travel and adult gut microbial communities, the seasonal cycling of gut microbial communities for the Hadza hunter-gatherers, and the distinctive microbial pattern of elderly gut microbial communities. In summary, microDELTA can largely solve the issues in tracing the life trajectory of the human microbial communities and generate accurate and flexible models for a broad spectrum of microbial-based longitudinal researches.

Context-dependent gene regulatory network reveals regulation dynamics and cell trajectories using unspliced transcripts.

Gene regulatory networks govern complex gene expression programs in various biological phenomena, including embryonic development, cell fate decisions and oncogenesis. Single-cell techniques are increasingly being used to study gene expression, providing higher resolution than traditional approaches. However, inferring a comprehensive gene regulatory network across different cell types remains a challenge. Here, we propose to construct context-dependent gene regulatory networks (CDGRNs) from single-cell RNA sequencing data utilizing both spliced and unspliced transcript expression levels. A gene regulatory network is decomposed into subnetworks corresponding to different transcriptomic contexts. Each subnetwork comprises the consensus active regulation pairs of transcription factors and their target genes shared by a group of cells, inferred by a Gaussian mixture model. We find that the union of gene regulation pairs in all contexts is sufficient to reconstruct differentiation trajectories. Functions specific to the cell cycle, cell differentiation or tissue-specific functions are enriched throughout the developmental process in each context. Surprisingly, we also observe that the network entropy of CDGRNs decreases along differentiation trajectories, indicating directionality in differentiation. Overall, CDGRN allows us to establish the connection between gene regulation at the molecular level and cell differentiation at the macroscopic level.

Screening single-cell trajectories via continuity assessments for cell transition potential.

Advances in single-cell sequencing and data analysis have made it possible to infer biological trajectories spanning heterogeneous cell populations based on transcriptome variation. These trajectories yield a wealth of novel insights into dynamic processes such as development and differentiation. However, trajectory analysis relies on an assumption of trajectory continuity, and experimental limitations preclude some real-world scenarios from meeting this condition. The current lack of assessment metrics makes it difficult to ascertain if/when a given trajectory deviates from continuity, and what impact such a divergence would have on inference accuracy is unclear. By analyzing simulated breaks introduced into in silico and real single-cell data, we found that discontinuity caused precipitous drops in the accuracy of trajectory inference. We then generate a simple scoring algorithm for assessing trajectory continuity, and found that continuity assessments in real-world cases of intestinal stem cell development and CD8 + T cells differentiation efficiently identifies trajectories consistent with empirical knowledge. This assessment approach can also be used in cases where a priori knowledge is lacking to screen a pool of inferred lineages for their adherence to presumed continuity, and serve as a means for weighing higher likelihood trajectories for validation via empirical studies, as exemplified by our case studies in psoriatic arthritis and acute kidney injury. This tool is freely available through github at qingshanni/scEGRET.

scMaSigPro: Differential Expression Analysis along Single-Cell Trajectories.

MOTIVATION: Understanding the dynamics of gene expression across different cellular states is crucial for discerning the mechanisms underneath cellular differentiation. Genes that exhibit variation in mean expression as a function of Pseudotime and between branching trajectories are expected to govern cell fate decisions. We introduce scMaSigPro, a method for the identification of differential gene expression patterns along Pseudotime and branching paths simultaneously. RESULTS: We assessed the performance of scMaSigPro using synthetic and public datasets. Our evaluation shows that scMaSigPro outperforms existing methods in controlling the False Positive Rate and is computationally efficient. AVAILABILITY AND IMPLEMENTATION: scMaSigPro is available as a free R package (version 4.0 or higher) under the GPL(≥2) license on GitHub at 'github.com/BioBam/scMaSigPro' and archived with version 0.03 on Zenodo at 'zenodo.org/records/12568922'.