Benchmarking computational methods for B-cell receptor reconstruction from single-cell RNA-seq data

The landscape of cell lineage tracing

Cell fate changes play a crucial role in the processes of natural development, disease progression, and the efficacy of therapeutic interventions. The definition of the various types of cell fate changes, including cell expansion, differentiation, transdifferentiation, dedifferentiation, reprogramming, and state transitions, represents a complex and evolving field of research known as cell lineage tracing. This review will systematically introduce the research history and progress in this field, which can be broadly divided into two parts: prospective tracing and retrospective tracing. The initial section encompasses an array of methodologies pertaining to isotope labeling, transient fluorescent tracers, non-fluorescent transient tracers, non-fluorescent genetic markers, fluorescent protein, genetic marker delivery, genetic recombination, exogenous DNA barcodes, CRISPR-Cas9 mediated DNA barcodes, and base editor-mediated DNA barcodes. The second part of the review covers genetic mosaicism, genomic DNA alteration, TCR/BCR, DNA methylation, and mitochondrial DNA mutation. In the final section, we will address the principal challenges and prospective avenues of enquiry in the field of cell lineage tracing, with a particular focus on the sequencing techniques and mathematical models pertinent to single-cell genetic lineage tracing, and the value of pursuing a more comprehensive investigation at both the spatial and temporal levels in the study of cell lineage tracing.

bcRflow: a Nextflow pipeline for characterizing B cell receptor repertoires from non-targeted transcriptomic data

B cells play a critical role in the adaptive recognition of foreign antigens through diverse receptor generation. While targeted immune sequencing methods are commonly used to profile B cell receptors (BCRs), they have limitations in cost and tissue availability. Analyzing B cell receptor profiling from non-targeted transcriptomics data is a promising alternative, but a systematic pipeline integrating tools for accurate immune repertoire extraction is lacking. Here, we present bcRflow, a Nextflow pipeline designed to characterize BCR repertoires from non-targeted transcriptomics data, with functional modules for alignment, processing, and visualization. bcRflow is a comprehensive, reproducible, and scalable pipeline that can run on high-performance computing clusters, cloud-based computing resources like Amazon Web Services (AWS), the Open OnDemand framework, or even local desktops. bcRflow utilizes institutional configurations provided by nf-core to ensure maximum portability and accessibility. To demonstrate the functionality of the bcRflow pipeline, we analyzed a public dataset of bulk transcriptomic samples from COVID-19 patients and healthy controls. We have shown that bcRflow streamlines the analysis of BCR repertoires from non-targeted transcriptomics data, providing valuable insights into the B cell immune response for biological and clinical research. bcRflow is available at https://github.com/Bioinformatics-Core-at-Childrens/bcRflow.

Single-cell immune repertoire analysis

Single-cell T cell and B cell antigen receptor-sequencing data analysis can potentially perform in-depth assessments of adaptive immune cells that inform on understanding immune cell development to tracking clonal expansion in disease and therapy. However, it has been extremely challenging to analyze and interpret T cells and B cells and their adaptive immune receptor repertoires at the single-cell level due to not only the complexity of the data but also the underlying biology. In this Review, we delve into the computational breakthroughs that have transformed the analysis of single-cell T cell and B cell antigen receptor-sequencing data.

Leveraging multi-omics data to empower quantitative systems pharmacology in immuno-oncology

Understanding the intricate interactions of cancer cells with the tumor microenvironment (TME) is a pre-requisite for the optimization of immunotherapy. Mechanistic models such as quantitative systems pharmacology (QSP) provide insights into the TME dynamics and predict the efficacy of immunotherapy in virtual patient populations/digital twins but require vast amounts of multimodal data for parameterization. Large-scale datasets characterizing the TME are available due to recent advances in bioinformatics for multi-omics data. Here, we discuss the perspectives of leveraging omics-derived bioinformatics estimates to inform QSP models and circumvent the challenges of model calibration and validation in immuno-oncology.

Systematic evaluation of B-cell clonal family inference approaches

The reconstruction of clonal families (CFs) in B-cell receptor (BCR) repertoire analysis is a crucial step to understand the adaptive immune system and how it responds to antigens. The BCR repertoire of an individual is formed throughout life and is diverse due to several factors such as gene recombination and somatic hypermutation. The use of Adaptive Immune Receptor Repertoire sequencing (AIRR-seq) using next generation sequencing enabled the generation of full BCR repertoires that also include rare CFs. The reconstruction of CFs from AIRR-seq data is challenging and several approaches have been developed to solve this problem. Currently, most methods use the heavy chain (HC) only, as it is more variable than the light chain (LC). CF reconstruction options include the definition of appropriate sequence similarity measures, the use of shared mutations among sequences, and the possibility of reconstruction without preliminary clustering based on V- and J-gene annotation. In this study, we aimed to systematically evaluate different approaches for CF reconstruction and to determine their impact on various outcome measures such as the number of CFs derived, the size of the CFs, and the accuracy of the reconstruction. The methods were compared to each other and to a method that groups sequences based on identical junction sequences and another method that only determines subclones. We found that after accounting for data set variability, in particular sequencing depth and mutation load, the reconstruction approach has an impact on part of the outcome measures, including the number of CFs. Simulations indicate that unique junctions and subclones should not be used as substitutes for CF and that more complex methods do not outperform simpler methods. Also, we conclude that different approaches differ in their ability to correctly reconstruct CFs when not considering the LC and to identify shared CFs. The results showed the effect of different approaches on the reconstruction of CFs and highlighted the importance of choosing an appropriate method.

IGH repertoire analysis at scale: deciphering the complexity of B cell infiltration and migration in esophageal squamous cell carcinoma

Tumor-infiltrating B-lineage cells have become predictors of prognosis and immunotherapy responses in various cancers. However, limited knowledge about their infiltration and migration patterns has hindered the understanding of their anti-tumor functions. Here, we examined the immunoglobulin heavy chain (IGH) repertoires in 496 multi-regional tumor, 107 normal tissue, and 48 metastatic lymph node samples obtained from 107 patients with esophageal squamous cell carcinoma (ESCC). Our study revealed higher IgG-type B-lineage cells infiltration in tumors than in healthy tissue, which was associated with improved patient outcomes. Genes such as ACTN1, COL6A5, and pathways like focal adhesion, which shapes the physical structure of tumors, could affect B-lineage cell infiltration. Notably, the IGH sequence was used as an identity-tag to monitor B cell migration, and their infiltration schema within the tumor were depicted based on our multi-regional tumor specimens. This analysis revealed an escalation in B cell clones overlapped between metastatic lymph nodes and tumors. Therefore, the Lymph Node Activation Index was defined, which could predict the outcomes of patients with lymph node metastasis. This research introduces a novel framework for probing B cell infiltration and migration within the tumor microenvironment using large-scale transcriptome data, while simultaneously providing fresh perspectives on B cell immunology within ESCC.

Best practices for single-cell analysis across modalities

Recent advances in single-cell technologies have enabled high-throughput molecular profiling of cells across modalities and locations. Single-cell transcriptomics data can now be complemented by chromatin accessibility, surface protein expression, adaptive immune receptor repertoire profiling and spatial information. The increasing availability of single-cell data across modalities has motivated the development of novel computational methods to help analysts derive biological insights. As the field grows, it becomes increasingly difficult to navigate the vast landscape of tools and analysis steps. Here, we summarize independent benchmarking studies of unimodal and multimodal single-cell analysis across modalities to suggest comprehensive best-practice workflows for the most common analysis steps. Where independent benchmarks are not available, we review and contrast popular methods. Our article serves as an entry point for novices in the field of single-cell (multi-)omic analysis and guides advanced users to the most recent best practices.

Antigen-specificity measurements are the key to understanding T cell responses

Antigen-specific T cells play a central role in the adaptive immune response and come in a wide range of phenotypes. T cell receptors (TCRs) mediate the antigen-specificities found in T cells. Importantly, high-throughput TCR sequencing provides a fingerprint which allows tracking of specific T cells and their clonal expansion in response to particular antigens. As a result, many studies have leveraged TCR sequencing in an attempt to elucidate the role of antigen-specific T cells in various contexts. Here, we discuss the published approaches to studying antigen-specific T cells and their specific TCR repertoire. Further, we discuss how these methods have been applied to study the TCR repertoire in various diseases in order to characterize the antigen-specific T cells involved in the immune control of disease.