Comment on 'rigorous benchmarking of T cell receptor repertoire profiling methods for cancer RNA sequencing'.

Transcriptome sequencing has become common in cancer research, resulting in the generation of a substantial volume of RNA sequencing (RNA-Seq) data. The ability to extract immune repertoires from these data is crucial for obtaining information on infiltrating T- and B-lymphocyte clones when dedicated amplicon T-cell/B-cell receptors sequencing (TCR-Seq/BCR-Seq) methods are unavailable. In response to this demand, several dedicated computational methods have been developed, including MiXCR, TRUST and ImRep. In the recent publication in Briefings in Bioinformatics, Peng et al. have conducted an intensive, systematic comparison of the three previously mentioned tools. Although their effort is commendable, we do have a few constructive critiques regarding technical elements of their analysis.

Targeted depletion of TRBV9+ T cells as immunotherapy in a patient with ankylosing spondylitis.

Autoimmunity is intrinsically driven by memory T and B cell clones inappropriately targeted at self-antigens. Selective depletion or suppression of self-reactive T cells remains a holy grail of autoimmune therapy, but disease-associated T cell receptors (TCRs) and cognate antigenic epitopes remained elusive. A TRBV9-containing CD8+ TCR motif was recently associated with the pathogenesis of ankylosing spondylitis, psoriatic arthritis and acute anterior uveitis, and cognate HLA-B*27-presented epitopes were identified. Following successful testing in nonhuman primate models, here we report human TRBV9+ T cell elimination in ankylosing spondylitis. The patient achieved remission within 3 months and ceased anti-TNF therapy after 5 years of continuous use. Complete remission has now persisted for 4 years, with three doses of anti-TRBV9 administered per year. We also observed a profound improvement in spinal mobility metrics and the Bath Ankylosing Spondylitis Metrology Index (BASMI). This represents a possibly curative therapy of an autoimmune disease via selective depletion of a TRBV-defined group of T cells. The anti-TRBV9 therapy could potentially be applicable to other HLA-B*27-associated spondyloarthropathies. Such targeted elimination of the underlying cause of the disease without systemic immunosuppression could offer a new generation of safe and efficient therapies for autoimmunity.

Defective peripheral B cell selection in common variable immune deficiency patients with autoimmune manifestations.

Common variable immune deficiency (CVID) is a heterogeneous disorder characterized by recurrent infections, low levels of serum immunoglobulins, and impaired vaccine responses. Autoimmune manifestations are common, but B cell central and peripheral selection mechanisms in CVID are incompletely understood. Here, we find that receptor editing, a measure of central tolerance, is increased in transitional B cells from CVID patients and that these cells have a higher immunoglobulin κ:λ ratio in CVID patients with autoimmune manifestations than in those with infection only. Contrariwise, the selection pressure in the germinal center on CD27bright memory B cells is decreased in CVID patients with autoimmune manifestations. Finally, functionally, T cell-dependent activation showed that naive B cells in CVID patients are badly equipped for activation and induction of mismatch repair genes. We conclude that central tolerance is functional whereas peripheral selection is defective in CVID patients with autoimmune manifestations, which could underpin the development of autoimmunity.

Different classes of genomic inserts contribute to human antibody diversity.

Recombination of antibody genes in B cells can involve distant genomic loci and contribute a foreign antigen-binding element to form hybrid antibodies with broad reactivity for Plasmodium falciparum. So far, antibodies containing the extracellular domain of the LAIR1 and LILRB1 receptors represent unique examples of cross-chromosomal antibody diversification. Here, we devise a technique to profile non-VDJ elements from distant genes in antibody transcripts. Independent of the preexposure of donors to malaria parasites, non-VDJ inserts were detected in 80% of individuals at frequencies of 1 in 104 to 105 B cells. We detected insertions in heavy, but not in light chain or T cell receptor transcripts. We classify the insertions into four types depending on the insert origin and destination: 1) mitochondrial and 2) nuclear DNA inserts integrated at VDJ junctions; 3) inserts originating from telomere proximal genes; and 4) fragile sites incorporated between J-to-constant junctions. The latter class of inserts was exclusively found in memory and in in vitro activated B cells, while all other classes were already detected in naïve B cells. More than 10% of inserts preserved the reading frame, including transcripts with signs of antigen-driven affinity maturation. Collectively, our study unravels a mechanism of antibody diversification that is layered on the classical V(D)J and switch recombination.

Benchmarking of T cell receptor repertoire profiling methods reveals large systematic biases.

Monitoring the T cell receptor (TCR) repertoire in health and disease can provide key insights into adaptive immune responses, but the accuracy of current TCR sequencing (TCRseq) methods is unclear. In this study, we systematically compared the results of nine commercial and academic TCRseq methods, including six rapid amplification of complementary DNA ends (RACE)-polymerase chain reaction (PCR) and three multiplex-PCR approaches, when applied to the same T cell sample. We found marked differences in accuracy and intra- and inter-method reproducibility for T cell receptor α (TRA) and T cell receptor ß (TRB) TCR chains. Most methods showed a lower ability to capture TRA than TRB diversity. Low RNA input generated non-representative repertoires. Results from the 5' RACE-PCR methods were consistent among themselves but differed from the RNA-based multiplex-PCR results. Using an in silico meta-repertoire generated from 108 replicates, we found that one genomic DNA-based method and two non-unique molecular identifier (UMI) RNA-based methods were more sensitive than UMI methods in detecting rare clonotypes, despite the better clonotype quantification accuracy of the latter.

Functionally specialized human CD4+ T-cell subsets express physicochemically distinct TCRs.

The organizational integrity of the adaptive immune system is determined by functionally discrete subsets of CD4+ T cells, but it has remained unclear to what extent lineage choice is influenced by clonotypically expressed T-cell receptors (TCRs). To address this issue, we used a high-throughput approach to profile the αß TCR repertoires of human naive and effector/memory CD4+ T-cell subsets, irrespective of antigen specificity. Highly conserved physicochemical and recombinatorial features were encoded on a subset-specific basis in the effector/memory compartment. Clonal tracking further identified forbidden and permitted transition pathways, mapping effector/memory subsets related by interconversion or ontogeny. Public sequences were largely confined to particular effector/memory subsets, including regulatory T cells (Tregs), which also displayed hardwired repertoire features in the naive compartment. Accordingly, these cumulative repertoire portraits establish a link between clonotype fate decisions in the complex world of CD4+ T cells and the intrinsic properties of somatically rearranged TCRs.

Two subsets of stem-like CD8+ memory T cell progenitors with distinct fate commitments in humans.

T cell memory relies on the generation of antigen-specific progenitors with stem-like properties. However, the identity of these progenitors has remained unclear, precluding a full understanding of the differentiation trajectories that underpin the heterogeneity of antigen-experienced T cells. We used a systematic approach guided by single-cell RNA-sequencing data to map the organizational structure of the human CD8+ memory T cell pool under physiological conditions. We identified two previously unrecognized subsets of clonally, epigenetically, functionally, phenotypically and transcriptionally distinct stem-like CD8+ memory T cells. Progenitors lacking the inhibitory receptors programmed death-1 (PD-1) and T cell immunoreceptor with Ig and ITIM domains (TIGIT) were committed to a functional lineage, whereas progenitors expressing PD-1 and TIGIT were committed to a dysfunctional, exhausted-like lineage. Collectively, these data reveal the existence of parallel differentiation programs in the human CD8+ memory T cell pool, with potentially broad implications for the development of immunotherapies and vaccines.

The Interplay between CD27dull and CD27bright B Cells Ensures the Flexibility, Stability, and Resilience of Human B Cell Memory.

Memory B cells (MBCs) epitomize the adaptation of the immune system to the environment. We identify two MBC subsets in peripheral blood, CD27dull and CD27bright MBCs, whose frequency changes with age. Heavy chain variable region (VH) usage, somatic mutation frequency replacement-to-silent ratio, and CDR3 property changes, reflecting consecutive selection of highly antigen-specific, low cross-reactive antibody variants, all demonstrate that CD27dull and CD27bright MBCs represent sequential MBC developmental stages, and stringent antigen-driven pressure selects CD27dull into the CD27bright MBC pool. Dynamics of human MBCs are exploited in pregnancy, when 50% of maternal MBCs are lost and CD27dull MBCs transit to the more differentiated CD27bright stage. In the postpartum period, the maternal MBC pool is replenished by the expansion of persistent CD27dull clones. Thus, the stability and flexibility of human B cell memory is ensured by CD27dull MBCs that expand and differentiate in response to change.

CD4+ T Follicular Helper Cells in Human Tonsils and Blood Are Clonally Convergent but Divergent from Non-Tfh CD4+ Cells.

T follicular helper (Tfh) cells are fundamental for B cell selection and antibody maturation in germinal centers. Circulating Tfh (cTfh) cells constitute a minor proportion of the CD4+ T cells in peripheral blood, but their clonotypic relationship to Tfh populations resident in lymph nodes and the extent to which they differ from non-Tfh CD4+ cells have been unclear. Using donor-matched blood and tonsil samples, we investigate T cell receptor (TCR) sharing between tonsillar Tfh cells and peripheral Tfh and non-Tfh cell populations. TCR transcript sequencing reveals considerable clonal overlap between peripheral and tonsillar Tfh cell subsets as well as a clear distinction between Tfh and non-Tfh cells. Furthermore, influenza-specific cTfh cell clones derived from blood can be found in the repertoire of tonsillar Tfh cells. Therefore, human blood samples can be used to gain insight into the specificity of Tfh responses occurring in lymphoid tissues, provided that cTfh subsets are studied.

CXCR3 Identifies Human Naive CD8+ T Cells with Enhanced Effector Differentiation Potential.

In mice, the ability of naive T (TN) cells to mount an effector response correlates with TCR sensitivity for self-derived Ags, which can be quantified indirectly by measuring surface expression levels of CD5. Equivalent findings have not been reported previously in humans. We identified two discrete subsets of human CD8+ TN cells, defined by the absence or presence of the chemokine receptor CXCR3. The more abundant CXCR3+ TN cell subset displayed an effector-like transcriptional profile and expressed TCRs with physicochemical characteristics indicative of enhanced interactions with peptide-HLA class I Ags. Moreover, CXCR3+ TN cells frequently produced IL-2 and TNF in response to nonspecific activation directly ex vivo and differentiated readily into Ag-specific effector cells in vitro. Comparative analyses further revealed that human CXCR3+ TN cells were transcriptionally equivalent to murine CXCR3+ TN cells, which expressed high levels of CD5. These findings provide support for the notion that effector differentiation is shaped by heterogeneity in the preimmune repertoire of human CD8+ T cells.

CD49b defines functionally mature Treg cells that survey skin and vascular tissues.

Regulatory T (Treg) cells prevent autoimmunity by limiting immune responses and inflammation in the secondary lymphoid organs and nonlymphoid tissues. While unique subsets of Treg cells have been described in some nonlymphoid tissues, their relationship to Treg cells in secondary lymphoid organs and circulation remains unclear. Furthermore, it is possible that Treg cells from similar tissue types share largely similar properties. We have identified a short-lived effector Treg cell subset that expresses the α2 integrin, CD49b, and exhibits a unique tissue distribution, being abundant in peripheral blood, vasculature, skin, and skin-draining lymph nodes, but uncommon in the intestines and in viscera-draining lymph nodes. CD49b+ Treg cells, which display superior functionality revealed by in vitro and in vivo assays, appear to develop after multiple rounds of cell division and TCR-dependent activation. Accordingly, single-cell RNA-seq analysis placed these cells at the apex of the Treg developmental trajectory. These results shed light on the identity and development of a functionally potent subset of mature effector Treg cells that recirculate through and survey peripheral tissues.

Comparative Analysis of B-Cell Receptor Repertoires Induced by Live Yellow Fever Vaccine in Young and Middle-Age Donors.

Age-related changes can significantly alter the state of adaptive immune system and often lead to attenuated response to novel pathogens and vaccination. In present study we employed 5'RACE UMI-based full length and nearly error-free immunoglobulin profiling to compare plasma cell antibody repertoires in young (19-26 years) and middle-age (45-58 years) individuals vaccinated with a live yellow fever vaccine, modeling a newly encountered pathogen. Our analysis has revealed age-related differences in the responding antibody repertoire ranging from distinct IGH CDR3 repertoire properties to differences in somatic hypermutation intensity and efficiency and antibody lineage tree structure. Overall, our findings suggest that younger individuals respond with a more diverse antibody repertoire and employ a more efficient somatic hypermutation process than elder individuals in response to a newly encountered pathogen.

Comparative analysis of murine T-cell receptor repertoires.

For understanding the rules and laws of adaptive immunity, high-throughput profiling of T-cell receptor (TCR) repertoires becomes a powerful tool. The structure of TCR repertoires is instructive even before the antigen specificity of each particular receptor becomes available. It embodies information about the thymic and peripheral selection of T cells; the readiness of an adaptive immunity to withstand new challenges; the character, magnitude and memory of immune responses; and the aetiological and functional proximity of T-cell subsets. Here, we describe our current analytical approaches for the comparative analysis of murine TCR repertoires, and show several examples of how these approaches can be applied for particular experimental settings. We analyse the efficiency of different metrics used for estimation of repertoire diversity, repertoire overlap, V-gene and J-gene segments usage similarity, and amino acid composition of CDR3. We discuss basic differences of these metrics and their advantages and limitations in different experimental models, and we provide guidelines for choosing an efficient way to lead a comparative analysis of TCR repertoires. Applied to the various known and newly developed mouse models, such analysis should allow us to disentangle multiple sophisticated puzzles in adaptive immunity.

Dynamics of Individual T Cell Repertoires: From Cord Blood to Centenarians.

The diversity, architecture, and dynamics of the TCR repertoire largely determine our ability to effectively withstand infections and malignancies with minimal mistargeting of immune responses. In this study, we have employed deep TCRß repertoire sequencing with normalization based on unique molecular identifiers to explore the long-term dynamics of T cell immunity. We demonstrate remarkable stability of repertoire, where approximately half of all T cells in peripheral blood are represented by clones that persist and generally preserve their frequencies for 3 y. We further characterize the extremes of lifelong TCR repertoire evolution, analyzing samples ranging from umbilical cord blood to centenarian peripheral blood. We show that the fetal TCR repertoire, albeit structurally maintained within regulated borders due to the lower numbers of randomly added nucleotides, is not limited with respect to observed functional diversity. We reveal decreased efficiency of nonsense-mediated mRNA decay in umbilical cord blood, which may reflect specific regulatory mechanisms in development. Furthermore, we demonstrate that human TCR repertoires are functionally more similar at birth but diverge during life, and we track the lifelong behavior of CMV- and EBV-specific T cell clonotypes. Finally, we reveal gender differences in dynamics of TCR diversity constriction, which come to naught in the oldest age. Based on our data, we propose a more general explanation for the previous observations on the relationships between longevity and immunity.

Quantitative profiling of immune repertoires for minor lymphocyte counts using unique molecular identifiers.

Emerging high-throughput sequencing methods for the analyses of complex structure of TCR and BCR repertoires give a powerful impulse to adaptive immunity studies. However, there are still essential technical obstacles for performing a truly quantitative analysis. Specifically, it remains challenging to obtain comprehensive information on the clonal composition of small lymphocyte populations, such as Ag-specific, functional, or tissue-resident cell subsets isolated by sorting, microdissection, or fine needle aspirates. In this study, we report a robust approach based on unique molecular identifiers that allows profiling Ag receptors for several hundred to thousand lymphocytes while preserving qualitative and quantitative information on clonal composition of the sample. We also describe several general features regarding the data analysis with unique molecular identifiers that are critical for accurate counting of starting molecules in high-throughput sequencing applications.