Pseudocell Tracer-A method for inferring dynamic trajectories using scRNAseq and its application to B cells undergoing immunoglobulin class switch recombination

Conserved B cell signaling, activation, and differentiation in porcine jejunal and ileal Peyer's patches despite distinct immune landscapes

Peyer's patches (PPs) are B cell-rich sites of intestinal immune induction, yet PP-associated B cell signaling, activation, and differentiation are poorly defined. Single-cell and spatial transcriptomics were completed to study B cells from porcine jejunum and ileum containing PPs. Intestinal locations had distinct immune landscapes, including more follicular B cells in ileum and increased MHC-II-encoding gene expression in jejunal B cells. Despite distinct landscapes, conserved B cell dynamics were detected across intestinal locations, including B cell signaling to CD4+ macrophages that are putative phagocytic, cytotoxic, effector cells and deduced routes of B cell activation/differentiation, including resting B cells migrating into follicles to replicate/divide or differentiate into antibody-secreting cells residing in intestinal crypts. A six-biomarker panel recapitulated transcriptomics findings of B cell phenotypes, frequencies, and spatial locations via ex vivo and in situ staining. Findings convey conserved B cell dynamics across intestinal locations containing PPs, despite location-specific immune environments. Results establish a benchmark of B cell dynamics for understanding intestinal immune induction important to promoting gut/overall health.

Isotype-aware inference of B cell clonal lineage trees from single-cell sequencing data

Single-cell RNA sequencing (scRNA-seq) enables comprehensive characterization of the micro-evolutionary processes of B cells during an adaptive immune response, capturing features of somatic hypermutation (SHM) and class switch recombination (CSR). Existing phylogenetic approaches for reconstructing B cell evolution have primarily focused on the SHM process alone. Here, we present tree inference of B cell clonal lineages (TRIBAL), an algorithm designed to optimally reconstruct the evolutionary history of B cell clonal lineages undergoing both SHM and CSR from scRNA-seq data. Through simulations, we demonstrate that TRIBAL produces more comprehensive and accurate B cell lineage trees compared to existing methods. Using real-world datasets, TRIBAL successfully recapitulates expected biological trends in a model affinity maturation system while reconstructing evolutionary histories with more parsimonious class switching than state-of-the-art methods. Thus, TRIBAL significantly improves B cell lineage tracing, useful for modeling vaccine responses, disease progression, and the identification of therapeutic antibodies.

Temporal dynamics and genomic programming of plasma cell fates

Affinity-matured plasma cells (PCs) of varying lifespans are generated through a germinal center (GC) response. The developmental dynamics and genomic programs of antigen-specific PC precursors remain to be elucidated. Here, using a model antigen in mice, we demonstrate biphasic generation of PC precursors, with those generating long-lived bone marrow PCs preferentially produced in the late phase of GC response. Clonal tracing using single-cell RNA sequencing and B cell antigen receptor sequencing in spleen and bone marrow compartments, coupled with adoptive transfer experiments, reveals a new PC transition state that gives rise to functionally competent PC precursors. The latter undergo clonal expansion, dependent on inducible expression of TIGIT. We propose a model for the proliferation and programming of precursors of long-lived PCs, based on extended antigen encounters in the GC.

sciCSR infers B cell state transition and predicts class-switch recombination dynamics using single-cell transcriptomic data

Class-switch recombination (CSR) is an integral part of B cell maturation. Here we present sciCSR (pronounced 'scissor', single-cell inference of class-switch recombination), a computational pipeline that analyzes CSR events and dynamics of B cells from single-cell RNA sequencing (scRNA-seq) experiments. Validated on both simulated and real data, sciCSR re-analyzes scRNA-seq alignments to differentiate productive heavy-chain immunoglobulin transcripts from germline 'sterile' transcripts. From a snapshot of B cell scRNA-seq data, a Markov state model is built to infer the dynamics and direction of CSR. Applying sciCSR on severe acute respiratory syndrome coronavirus 2 vaccination time-course scRNA-seq data, we observe that sciCSR predicts, using data from an earlier time point in the collected time-course, the isotype distribution of B cell receptor repertoires of subsequent time points with high accuracy (cosine similarity ~0.9). Using processes specific to B cells, sciCSR identifies transitions that are often missed by conventional RNA velocity analyses and can reveal insights into the dynamics of B cell CSR during immune response.

TRIBAL: Tree Inference of B cell Clonal Lineages

B cells are a critical component of the adaptive immune system, responsible for producing antibodies that help protect the body from infections and foreign substances. Single cell RNA-sequencing (scRNA-seq) has allowed for both profiling of B cell receptor (BCR) sequences and gene expression. However, understanding the adaptive and evolutionary mechanisms of B cells in response to specific stimuli remains a significant challenge in the field of immunology. We introduce a new method, TRIBAL, which aims to infer the evolutionary history of clonally related B cells from scRNA-seq data. The key insight of TRIBAL is that inclusion of isotype data into the B cell lineage inference problem is valuable for reducing phylogenetic uncertainty that arises when only considering the receptor sequences. Consequently, the TRIBAL inferred B cell lineage trees jointly capture the somatic mutations introduced to the B cell receptor during affinity maturation and isotype transitions during class switch recombination. In addition, TRIBAL infers isotype transition probabilities that are valuable for gaining insight into the dynamics of class switching. Via in silico experiments, we demonstrate that TRIBAL infers isotype transition probabilities with the ability to distinguish between direct versus sequential switching in a B cell population. This results in more accurate B cell lineage trees and corresponding ancestral sequence and class switch reconstruction compared to competing methods. Using real-world scRNA-seq datasets, we show that TRIBAL recapitulates expected biological trends in a model affinity maturation system. Furthermore, the B cell lineage trees inferred by TRIBAL were equally plausible for the BCR sequences as those inferred by competing methods but yielded lower entropic partitions for the isotypes of the sequenced B cell. Thus, our method holds the potential to further advance our understanding of vaccine responses, disease progression, and the identification of therapeutic antibodies.

Temporal dynamics and genomic programming of plasma cell fates

Affinity-matured plasma cells (PCs) of varying lifespans are generated through a germinal center (GC) response. The developmental dynamics and genomic programs of antigen-specific PC precursors remain to be elucidated. Using a model antigen, we demonstrate biphasic generation of PC precursors, with those generating long-lived bone marrow PCs preferentially produced in the late phase of GC response. Clonal tracing using scRNA-seq+BCR-seq in spleen and bone marrow compartments, coupled with adoptive transfer experiments, reveal a novel PC transition state that gives rise to functionally competent PC precursors. The latter undergo clonal expansion, dependent on inducible expression of TIGIT. We propose a model for the proliferation and programming of precursors of long-lived PCs, based on extended antigen encounters followed by reduced antigen availability.

Cell fate dynamics and genomic programming of plasma cell precursors

This review is focused on the cellular dynamics and genomic programming of plasma cell (PC) precursors that arise during germinal center (GC) B cell responses in secondary lymphoid organs (SLOs) and give rise to PCs in the bone marrow. Considerable progress has been made in the phenotypic characterization of circulating and bone marrow PC precursors as well as their differentiated short-lived (SLPC) and long-lived (LLPC) counterparts, in the context of model antigen and vaccine responses. Importantly, it has been possible to infer the temporal dynamics of generation of PC precursors during a GC response. However, the nature of the PC precursors at their site of generation in SLOs, and their signaling and genomic states, remain to be elucidated. Our synthesis draws upon experimental studies conducted in murine models as well as in humans, the latter complemented with cell culture manipulations of PCs and their precursors. By integration of the studies in murine and human systems, which are being accelerated by new genomic methodologies, we highlight insights and hypotheses concerning the generation of PCs. This framework can be extended and explored from both fundamental and translational standpoints.