Tracking global changes induced in the CD4 T-cell receptor repertoire by immunization with a complex antigen using short stretches of CDR3 protein sequence

Principles and therapeutic applications of adaptive immunity

Adaptive immunity provides protection against infectious and malignant diseases. These effects are mediated by lymphocytes that sense and respond with targeted precision to perturbations induced by pathogens and tissue damage. Here, we review key principles underlying adaptive immunity orchestrated by distinct T cell and B cell populations and their extensions to disease therapies. We discuss the intracellular and intercellular processes shaping antigen specificity and recognition in immune activation and lymphocyte functions in mediating effector and memory responses. We also describe how lymphocytes balance protective immunity against autoimmunity and immunopathology, including during immune tolerance, response to chronic antigen stimulation, and adaptation to non-lymphoid tissues in coordinating tissue immunity and homeostasis. Finally, we discuss extracellular signals and cell-intrinsic programs underpinning adaptive immunity and conclude by summarizing key advances in vaccination and engineering adaptive immune responses for therapeutic interventions. A deeper understanding of these principles holds promise for uncovering new means to improve human health.

Influence of pre-B cell receptor deficiency on the immunoglobulin repertoires in peripheral blood B cells before and after immunization

During B cell development, pre-B cell receptor (pre-BCR), comprising the immunoglobulin heavy chain (HC) and surrogate light chain (SLC), plays a crucial role. The expression of pre-BCR serves as a certification of HC quality, confirming its ability to associate with the SLC and light chain (LC). In mice lacking SLC, the absence of this quality control mechanism leads to a distorted repertoire of HCs in the spleen and bone marrow. In this study, we conducted a comparative analysis of the immunoglobulin gene repertoire in peripheral blood cells of both wild-type mice and pre-BCR-deficient mice. Our findings reveal differences not only in the μ HC repertoire but also in the α HC and κ LC repertoires of the pre-BCR-deficient mice. These results suggest that the pre-BCR-mediated quality check of HC influences the selection of class-switched HC and LC repertoires. To further explore the impact of pre-BCR deficiency, we immunized these mice with thymus-dependent antigens and compared the antigen-responding repertoires. Our observations indicate that the affinity maturation pathways remain consistent between wild-type mice and pre-BCR-deficient mice, albeit with variations in the degree of maturation.

Quantitative systems-based prediction of antimicrobial resistance evolution

Predicting evolution is a fundamental problem in biology with practical implications for treating antimicrobial resistance, which is a complex system-level phenomenon. In this perspective article, we explore the limits of predicting antimicrobial resistance evolution, quantitatively define the predictability and repeatability of microevolutionary processes, and speculate on how these quantities vary across temporal, biological, and complexity scales. The opportunities and challenges for predicting antimicrobial resistance in the context of systems biology are also discussed. Based on recent research, we conclude that the evolution of antimicrobial resistance can be predicted using a systems biology approach integrating quantitative models with multiscale data from microbial evolution experiments.

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.

Signatures of T cell immunity revealed using sequence similarity with TCRDivER algorithm

Changes in the T cell receptor (TCR) repertoires have become important markers for monitoring disease or therapy progression. With the rise of immunotherapy usage in cancer, infectious and autoimmune disease, accurate assessment and comparison of the "state" of the TCR repertoire has become paramount. One important driver of change within the repertoire is T cell proliferation following immunisation. A way of monitoring this is by investigating large clones of individual T cells believed to bind epitopes connected to the disease. However, as a single target can be bound by many different TCRs, monitoring individual clones cannot fully account for T cell cross-reactivity. Moreover, T cells responding to the same target often exhibit higher sequence similarity, which highlights the importance of accounting for TCR similarity within the repertoire. This complexity of binding relationships between a TCR and its target convolutes comparison of immune responses between individuals or comparisons of TCR repertoires at different timepoints. Here we propose TCRDivER algorithm (T cell Receptor Diversity Estimates for Repertoires), a global method of T cell repertoire comparison using diversity profiles sensitive to both clone size and sequence similarity. This approach allowed for distinction between spleen TCR repertoires of immunised and non-immunised mice, showing the need for including both facets of repertoire changes simultaneously. The analysis revealed biologically interpretable relationships between sequence similarity and clonality. These aid in understanding differences and separation of repertoires stemming from different biological context. With the rise of availability of sequencing data we expect our tool to find broad usage in clinical and research applications.

TCR sequencing: applications in immuno-oncology research

•T-cell receptor (TCR) interaction with major histocompatibility complex-antigen complexes leads to antitumour responses.•TCR sequencing analysis allows characterisation of T cells that recognise tumour neoantigens.•T-cell clonal revival and clonal replacement potentially underpin immunotherapy responses.

Quantifying changes in the T cell receptor repertoire during thymic development

One of the feats of adaptive immunity is its ability to recognize foreign pathogens while sparing the self. During maturation in the thymus, T cells are selected through the binding properties of their antigen-specific T-cell receptor (TCR), through the elimination of both weakly (positive selection) and strongly (negative selection) self-reactive receptors. However, the impact of thymic selection on the TCR repertoire is poorly understood. Here, we use transgenic Nur77-mice expressing a T-cell activation reporter to study the repertoires of thymic T cells at various stages of their development, including cells that do not pass selection. We combine high-throughput repertoire sequencing with statistical inference techniques to characterize the selection of the TCR in these distinct subsets. We find small but significant differences in the TCR repertoire parameters between the maturation stages, which recapitulate known differentiation pathways leading to the CD4⁺ and CD8⁺ subtypes. These differences can be simulated by simple models of selection acting linearly on the sequence features. We find no evidence of specific sequences or sequence motifs or features that are suppressed by negative selection. These results favour a collective or statistical model for T-cell self non-self discrimination, where negative selection biases the repertoire away from self recognition, rather than ensuring lack of self-reactivity at the single-cell level.

Dynamic kernel matching for non-conforming data: A case study of T cell receptor datasets

Most statistical classifiers are designed to find patterns in data where numbers fit into rows and columns, like in a spreadsheet, but many kinds of data do not conform to this structure. To uncover patterns in non-conforming data, we describe an approach for modifying established statistical classifiers to handle non-conforming data, which we call dynamic kernel matching (DKM). As examples of non-conforming data, we consider (i) a dataset of T-cell receptor (TCR) sequences labelled by disease antigen and (ii) a dataset of sequenced TCR repertoires labelled by patient cytomegalovirus (CMV) serostatus, anticipating that both datasets contain signatures for diagnosing disease. We successfully fit statistical classifiers augmented with DKM to both datasets and report the performance on holdout data using standard metrics and metrics allowing for indeterminant diagnoses. Finally, we identify the patterns used by our statistical classifiers to generate predictions and show that these patterns agree with observations from experimental studies.

simAIRR: simulation of adaptive immune repertoires with realistic receptor sequence sharing for benchmarking of immune state prediction methods

BACKGROUND

Machine learning (ML) has gained significant attention for classifying immune states in adaptive immune receptor repertoires (AIRRs) to support the advancement of immunodiagnostics and therapeutics. Simulated data are crucial for the rigorous benchmarking of AIRR-ML methods. Existing approaches to generating synthetic benchmarking datasets result in the generation of naive repertoires missing the key feature of many shared receptor sequences (selected for common antigens) found in antigen-experienced repertoires.

RESULTS

We demonstrate that a common approach to generating simulated AIRR benchmark datasets can introduce biases, which may be exploited for undesired shortcut learning by certain ML methods. To mitigate undesirable access to true signals in simulated AIRR datasets, we devised a simulation strategy (simAIRR) that constructs antigen-experienced-like repertoires with a realistic overlap of receptor sequences. simAIRR can be used for constructing AIRR-level benchmarks based on a range of assumptions (or experimental data sources) for what constitutes receptor-level immune signals. This includes the possibility of making or not making any prior assumptions regarding the similarity or commonality of immune state-associated sequences that will be used as true signals. We demonstrate the real-world realism of our proposed simulation approach by showing that basic ML strategies perform similarly on simAIRR-generated and real-world experimental AIRR datasets.

CONCLUSIONS

This study sheds light on the potential shortcut learning opportunities for ML methods that can arise with the state-of-the-art way of simulating AIRR datasets. simAIRR is available as a Python package: https://github.com/KanduriC/simAIRR.

Comparing T cell receptor repertoires using optimal transport

The complexity of entire T cell receptor (TCR) repertoires makes their comparison a difficult but important task. Current methods of TCR repertoire comparison can incur a high loss of distributional information by considering overly simplistic sequence- or repertoire-level characteristics. Optimal transport methods form a suitable approach for such comparison given some distance or metric between values in the sample space, with appealing theoretical and computational properties. In this paper we introduce a nonparametric approach to comparing empirical TCR repertoires that applies the Sinkhorn distance, a fast, contemporary optimal transport method, and a recently-created distance between TCRs called TCRdist. We show that our methods identify meaningful differences between samples from distinct TCR distributions for several case studies, and compete with more complicated methods despite minimal modeling assumptions and a simpler pipeline.

Reference-based comparison of adaptive immune receptor repertoires

B and T cell receptor (immune) repertoires can represent an individual's immune history. While current repertoire analysis methods aim to discriminate between health and disease states, they are typically based on only a limited number of parameters. Here, we introduce immuneREF: a quantitative multidimensional measure of adaptive immune repertoire (and transcriptome) similarity that allows interpretation of immune repertoire variation by relying on both repertoire features and cross-referencing of simulated and experimental datasets. To quantify immune repertoire similarity landscapes across health and disease, we applied immuneREF to >2,400 datasets from individuals with varying immune states (healthy, [autoimmune] disease, and infection). We discovered, in contrast to the current paradigm, that blood-derived immune repertoires of healthy and diseased individuals are highly similar for certain immune states, suggesting that repertoire changes to immune perturbations are less pronounced than previously thought. In conclusion, immuneREF enables the population-wide study of adaptive immune response similarity across immune states.

Disorder and the Neural Representation of Complex Odors

Animals smelling in the real world use a small number of receptors to sense a vast number of natural molecular mixtures, and proceed to learn arbitrary associations between odors and valences. Here, we propose how the architecture of olfactory circuits leverages disorder, diffuse sensing and redundancy in representation to meet these immense complementary challenges. First, the diffuse and disordered binding of receptors to many molecules compresses a vast but sparsely-structured odor space into a small receptor space, yielding an odor code that preserves similarity in a precise sense. Introducing any order/structure in the sensing degrades similarity preservation. Next, lateral interactions further reduce the correlation present in the low-dimensional receptor code. Finally, expansive disordered projections from the periphery to the central brain reconfigure the densely packed information into a high-dimensional representation, which contains multiple redundant subsets from which downstream neurons can learn flexible associations and valences. Moreover, introducing any order in the expansive projections degrades the ability to recall the learned associations in the presence of noise. We test our theory empirically using data from Drosophila. Our theory suggests that the neural processing of sparse but high-dimensional olfactory information differs from the other senses in its fundamental use of disorder.

A hierarchy of selection pressures determines the organization of the T cell receptor repertoire

We systematically examine the receptor repertoire in T cell subsets in young, adult, and LCMV-infected mice. Somatic recombination generates diversity, resulting in the limited overlap between nucleotide sequences of different repertoires even within the same individual. However, statistical features of the repertoire, quantified by the V gene and CDR3 k-mer frequency distributions, are highly conserved. A hierarchy of immunological processes drives the evolution of this structure. Intra-thymic divergence of CD4+ and CD8+ lineages imposes subtle but dominant differences observed across repertoires of all subpopulations in both young and adult mice. Differentiation from naive through memory to effector phenotype imposes an additional gradient of repertoire diversification, which is further influenced by age in a complex and lineage-dependent manner. The distinct repertoire of CD4+ regulatory T cells is more similar to naive cells in young mice and to effectors in adults. Finally, we describe divergent (naive and memory) and convergent (CD8+ effector) evolution of the repertoire following acute infection with LCMV. This study presents a quantitative framework that captures the structure of the repertoire in terms of its fundamental statistical properties and describes how this structure evolves as individual T cells differentiate, migrate and mature in response to antigen exposure.

Comparative Study of Repertoire Classification Methods Reveals Data Efficiency of k-mer Feature Extraction

The repertoire of T cell receptors encodes various types of immunological information. Machine learning is indispensable for decoding such information from repertoire datasets measured by next-generation sequencing (NGS). In particular, the classification of repertoires is the most basic task, which is relevant for a variety of scientific and clinical problems. Supported by the recent appearance of large datasets, efficient but data-expensive methods have been proposed. However, it is unclear whether they can work efficiently when the available sample size is severely restricted as in practical situations. In this study, we demonstrate that their performances can be impaired substantially below critical sample sizes. To complement this drawback, we propose MotifBoost, which exploits the information of short k-mer motifs of TCRs. MotifBoost can perform the classification as efficiently as a deep learning method on large datasets while providing more stable and reliable results on small datasets. We tested MotifBoost on the four small datasets which consist of various conditions such as Cytomegalovirus (CMV), HIV, α-chain, β-chain and it consistently preserved the stability. We also clarify that the robustness of MotifBoost can be attributed to the efficiency of k-mer motifs as representation features of repertoires. Finally, by comparing the predictions of these methods, we show that the whole sequence identity and sequence motifs encode partially different information and that a combination of such complementary information is necessary for further development of repertoire analysis.

Profiling the baseline performance and limits of machine learning models for adaptive immune receptor repertoire classification

BACKGROUND

Machine learning (ML) methodology development for the classification of immune states in adaptive immune receptor repertoires (AIRRs) has seen a recent surge of interest. However, so far, there does not exist a systematic evaluation of scenarios where classical ML methods (such as penalized logistic regression) already perform adequately for AIRR classification. This hinders investigative reorientation to those scenarios where method development of more sophisticated ML approaches may be required.

RESULTS

To identify those scenarios where a baseline ML method is able to perform well for AIRR classification, we generated a collection of synthetic AIRR benchmark data sets encompassing a wide range of data set architecture-associated and immune state-associated sequence patterns (signal) complexity. We trained ≈1,700 ML models with varying assumptions regarding immune signal on ≈1,000 data sets with a total of ≈250,000 AIRRs containing ≈46 billion TCRβ CDR3 amino acid sequences, thereby surpassing the sample sizes of current state-of-the-art AIRR-ML setups by two orders of magnitude. We found that L1-penalized logistic regression achieved high prediction accuracy even when the immune signal occurs only in 1 out of 50,000 AIR sequences.

CONCLUSIONS

We provide a reference benchmark to guide new AIRR-ML classification methodology by (i) identifying those scenarios characterized by immune signal and data set complexity, where baseline methods already achieve high prediction accuracy, and (ii) facilitating realistic expectations of the performance of AIRR-ML models given training data set properties and assumptions. Our study serves as a template for defining specialized AIRR benchmark data sets for comprehensive benchmarking of AIRR-ML methods.

Tertiary lymphoid structure and decreased CD8+ T cell infiltration in minimally invasive adenocarcinoma

Knowledge of the tumor microenvironment (TME) in patients with early lung cancer, especially in comparison with the matched adjacent tissues, remains lacking. To characterize TME of early-stage lung adenocarcinoma, we performed RNA-seq profiling on 58 pairs of minimally invasive adenocarcinoma (MIA) tumors and matched adjacent normal tissues. MIA tumors exhibited an adaptive TME characterized by high CD4⁺ T cell infiltration, high B-cell activation, and low CD8⁺ T cell infiltration. The high expression of markers for B cells, activated CD4⁺ T cells, and follicular helper T (Tfh) cells in bulk MIA samples and three independent single-cell RNA-seq datasets implied tertiary lymphoid structures (TLS) formation. Multiplex immunohistochemistry staining validated TLS formation and revealed an enrichment of follicular regulatory T cells (Tfr) in TLS follicles, which may explain the lower CD8⁺ T cell infiltration and attenuated anti-tumor immunity in MIA. Our study demonstrates how integrating transcriptome and pathology characterize TME and elucidate potential mechanisms of tumor immune evasion.

•

Higher infiltration and activation of B and CD4⁺ T cell characterize MIA tumors

•

MIA tumors are infiltrated with lower CD8⁺ T cells than normal tissues

•

TLS constituted by B, CD4⁺ T cells, and CD35⁺ FDC is validated in MIA tumors

•

Decreased CD8⁺ T is associated with Tfr-mediated immunosuppression in MIA tumor

Application of T cell receptor (TCR) repertoire analysis for the advancement of cancer immunotherapy

T cell receptor (TCR) sequencing has emerged as a powerful new technology in analysis of the host-tumour interaction. The advances in NextGen sequencing technologies, coupled with powerful novel bioinformatic tools, allow quantitative and reproducible characterisation of repertoires from tumour and blood samples from an increasing number of patients with a variety of solid cancers. In this review, we consider how global metrics such as T cell clonality and diversity can be extracted from these repertoires and used to give insight into the mechanism of action of immune checkpoint blockade. Furthermore, we explore how the analysis of TCR overlap between repertories can help define spatial and temporal heterogeneity of the anti-tumoural immune response. Finally, we review how analysis of TCR sequence and structure, either of individual TCRs or from sets of related TCRs can be used to annotate the antigenic specificity, with important implications for the development of personalised adoptive cellular immunotherapies.

SARS-CoV-2 antigen exposure history shapes phenotypes and specificity of memory CD8 T cells

Although mRNA vaccine efficacy against severe COVID-19 remains high, variant emergence and breakthrough infections have changed vaccine policy to include booster immunizations. However, the effect of diverse and repeated antigen exposures on SARS-CoV-2 memory T cells is poorly understood. Here, we utilize DNA-barcoded MHC-multimers combined with scRNAseq and scTCRseq to capture the ex vivo profile of SARS-CoV-2-responsive T cells within a cohort of individuals with one, two, or three antigen exposures, including vaccination, primary infection, and breakthrough infection. We found that the order of exposure determined the relative distribution between spike- and non-spike-specific responses, with vaccination after infection leading to further expansion of spike-specific T cells and differentiation to a CCR7-CD45RA+ effector phenotype. In contrast, individuals experiencing a breakthrough infection mount vigorous non-spike-specific responses. In-depth analysis of over 4,000 epitope-specific T cell receptor sequences demonstrates that all types of exposures elicit diverse repertoires characterized by shared, dominant TCR motifs, with no evidence for repertoire narrowing from repeated exposure. Our findings suggest that breakthrough infections diversify the T cell memory repertoire and that current vaccination protocols continue to expand and differentiate spike-specific memory responses.

Grouping T-Cell Antigen Receptors by Specificity

Grouping TCRs on the similarity of CDR3 sequences could effectively cluster them by specificity. Three versions of the GLIPH algorithm are described briefly here, with instructions to use GLIPH algorithms to cluster TCRs by specificity.

Predicting T Cell Receptor Antigen Specificity From Structural Features Derived From Homology Models of Receptor-Peptide-Major Histocompatibility Complexes

The physical interaction between the T cell receptor (TCR) and its cognate antigen causes T cells to activate and participate in the immune response. Understanding this physical interaction is important in predicting TCR binding to a target epitope, as well as potential cross-reactivity. Here, we propose a way of collecting informative features of the binding interface from homology models of T cell receptor-peptide-major histocompatibility complex (TCR-pMHC) complexes. The information collected from these structures is sufficient to discriminate binding from non-binding TCR-pMHC pairs in multiple independent datasets. The classifier is limited by the number of crystal structures available for the homology modelling and by the size of the training set. However, the classifier shows comparable performance to sequence-based classifiers requiring much larger training sets.

Receptor repertoires of murine follicular T helper cells reveal a high clonal overlap in separate lymph nodes in autoimmunity

Follicular T helper cells (Tfh) are a specialized subset of CD4 effector T cells that are crucial for germinal center (GC) reactions and for selecting B cells to undergo affinity maturation. Despite this central role for humoral immunity, only few data exist about their clonal distribution when multiple lymphoid organs are exposed to the same antigen (Ag) as it is the case in autoimmunity. Here, we used an autoantibody-mediated disease model of the skin and injected one auto-Ag into the two footpads of the same mouse and analyzed the T cell receptor (TCR)β sequences of Tfh located in GCs of both contralateral draining lymph nodes. We found that over 90% of the dominant GC-Tfh clonotypes were shared in both lymph nodes but only transiently. The initially dominant Tfh clonotypes especially declined after establishment of chronic disease while GC reaction and autoimmune disease continued. Our data demonstrates a dynamic behavior of Tfh clonotypes under autoimmune conditions and emphasizes the importance of the time point for distinguishing auto-Ag-specific Tfh clonotypes from potential bystander activated ones.

Machine Learning in Discovery of New Antivirals and Optimization of Viral Infections Therapy

Nowadays, computational approaches play an important role in the design of new drug-like compounds and optimization of pharmacotherapeutic treatment of diseases. The emerging growth of viral infections, including those caused by the Human Immunodeficiency Virus (HIV), Ebola virus, recently detected coronavirus, and some others, leads to many newly infected people with a high risk of death or severe complications. A huge amount of chemical, biological, clinical data is at the disposal of the researchers. Therefore, there are many opportunities to find the relationships between the particular features of chemical data and the antiviral activity of biologically active compounds based on machine learning approaches. Biological and clinical data can also be used for building models to predict relationships between viral genotype and drug resistance, which might help determine the clinical outcome of treatment. In the current study, we consider machine-learning approaches in the antiviral research carried out during the past decade. We overview in detail the application of machine-learning methods for the design of new potential antiviral agents and vaccines, drug resistance prediction, and analysis of virus-host interactions. Our review also covers the perspectives of using the machine-learning approaches for antiviral research, including Dengue, Ebola viruses, Influenza A, Human Immunodeficiency Virus, coronaviruses, and some others.

Clustering based approach for population level identification of condition-associated T-cell receptor β-chain CDR3 sequences

BACKGROUND

Deep immune receptor sequencing, RepSeq, provides unprecedented opportunities for identifying and studying condition-associated T-cell clonotypes, represented by T-cell receptor (TCR) CDR3 sequences. However, due to the immense diversity of the immune repertoire, identification of condition relevant TCR CDR3s from total repertoires has mostly been limited to either "public" CDR3 sequences or to comparisons of CDR3 frequencies observed in a single individual. A methodology for the identification of condition-associated TCR CDR3s by direct population level comparison of RepSeq samples is currently lacking.

RESULTS

We present a method for direct population level comparison of RepSeq samples using immune repertoire sub-units (or sub-repertoires) that are shared across individuals. The method first performs unsupervised clustering of CDR3s within each sample. It then finds matching clusters across samples, called immune sub-repertoires, and performs statistical differential abundance testing at the level of the identified sub-repertoires. It finally ranks CDR3s in differentially abundant sub-repertoires for relevance to the condition. We applied the method on total TCR CDR3β RepSeq datasets of celiac disease patients, as well as on public datasets of yellow fever vaccination. The method successfully identified celiac disease associated CDR3β sequences, as evidenced by considerable agreement of TRBV-gene and positional amino acid usage patterns in the detected CDR3β sequences with previously known CDR3βs specific to gluten in celiac disease. It also successfully recovered significantly high numbers of previously known CDR3β sequences relevant to each condition than would be expected by chance.

CONCLUSION

We conclude that immune sub-repertoires of similar immuno-genomic features shared across unrelated individuals can serve as viable units of immune repertoire comparison, serving as proxy for identification of condition-associated CDR3s.

Machine Learning Analysis of Naïve B-Cell Receptor Repertoires Stratifies Celiac Disease Patients and Controls

Celiac disease (CeD) is a common autoimmune disorder caused by an abnormal immune response to dietary gluten proteins. The disease has high heritability. HLA is the major susceptibility factor, and the HLA effect is mediated via presentation of deamidated gluten peptides by disease-associated HLA-DQ variants to CD4+ T cells. In addition to gluten-specific CD4+ T cells the patients have antibodies to transglutaminase 2 (autoantigen) and deamidated gluten peptides. These disease-specific antibodies recognize defined epitopes and they display common usage of specific heavy and light chains across patients. Interactions between T cells and B cells are likely central in the pathogenesis, but how the repertoires of naïve T and B cells relate to the pathogenic effector cells is unexplored. To this end, we applied machine learning classification models to naïve B cell receptor (BCR) repertoires from CeD patients and healthy controls. Strikingly, we obtained a promising classification performance with an F1 score of 85%. Clusters of heavy and light chain sequences were inferred and used as features for the model, and signatures associated with the disease were then characterized. These signatures included amino acid (AA) 3-mers with distinct bio-physiochemical characteristics and enriched V and J genes. We found that CeD-associated clusters can be identified and that common motifs can be characterized from naïve BCR repertoires. The results may indicate a genetic influence by BCR encoding genes in CeD. Analysis of naïve BCRs as presented here may become an important part of assessing the risk of individuals to develop CeD. Our model demonstrates the potential of using BCR repertoires and in particular, naïve BCR repertoires, as disease susceptibility markers.

The clonal structure and dynamics of the human T cell response to an organic chemical hapten

Diphenylcyclopropenone (DPC) is an organic chemical hapten which induces allergic contact dermatitis and is used in the treatment of warts, melanoma, and alopecia areata. This therapeutic setting therefore provided an opportunity to study T cell receptor (TCR) repertoire changes in response to hapten sensitization in humans. Repeated exposure to DPC induced highly dynamic transient expansions of a polyclonal diverse T cell population. The number of TCRs expanded early after sensitization varies between individuals and predicts the magnitude of the allergic reaction. The expanded TCRs show preferential TCR V and J gene usage and consist of clusters of TCRs with similar sequences, two characteristic features of antigen-driven responses. The expanded TCRs share subtle sequence motifs that can be captured using a dynamic Bayesian network. These observations suggest the response to DPC is mediated by a polyclonal population of T cells recognizing a small number of dominant antigens.

Mapping the functional landscape of T cell receptor repertoires by single-T cell transcriptomics

Many experimental and bioinformatics approaches have been developed to characterize the human T cell receptor (TCR) repertoire. However, the unknown functional relevance of TCR profiling hinders unbiased interpretation of the biology of T cells. To address this inadequacy, we developed tessa, a tool to integrate TCRs with gene expression of T cells to estimate the effect that TCRs confer on the phenotypes of T cells. Tessa leveraged techniques combining single-cell RNA-sequencing with TCR sequencing. We validated tessa and showed its superiority over existing approaches that investigate only the TCR sequences. With tessa, we demonstrated that TCR similarity constrains the phenotypes of T cells to be similar and dictates a gradient in antigen targeting efficiency of T cell clonotypes with convergent TCRs. We showed this constraint could predict a functional dichotomization of T cells postimmunotherapy treatment and is weakened in tumor contexts.

Early life imprints the hierarchy of T cell clone sizes

The adaptive immune system responds to pathogens by selecting clones of cells with specific receptors. While clonal selection in response to particular antigens has been studied in detail, it is unknown how a lifetime of exposures to many antigens collectively shape the immune repertoire. Here, using mathematical modeling and statistical analyses of T cell receptor sequencing data, we develop a quantitative theory of human T cell dynamics compatible with the statistical laws of repertoire organization. We find that clonal expansions during a perinatal time window leave a long-lasting imprint on the human T cell repertoire, which is only slowly reshaped by fluctuating clonal selection during adult life. Our work provides a mechanism for how early clonal dynamics imprint the hierarchy of T cell clone sizes with implications for pathogen defense and autoimmunity.

Population variability in the generation and selection of T-cell repertoires

The diversity of T-cell receptor (TCR) repertoires is achieved by a combination of two intrinsically stochastic steps: random receptor generation by VDJ recombination, and selection based on the recognition of random self-peptides presented on the major histocompatibility complex. These processes lead to a large receptor variability within and between individuals. However, the characterization of the variability is hampered by the limited size of the sampled repertoires. We introduce a new software tool SONIA to facilitate inference of individual-specific computational models for the generation and selection of the TCR beta chain (TRB) from sequenced repertoires of 651 individuals, separating and quantifying the variability of the two processes of generation and selection in the population. We find not only that most of the variability is driven by the VDJ generation process, but there is a large degree of consistency between individuals with the inter-individual variance of repertoires being about ∼2% of the intra-individual variance. Known viral-specific TCRs follow the same generation and selection statistics as all TCRs.

A Conserved TCRβ Signature Dominates a Highly Polyclonal T-Cell Expansion During the Acute Phase of a Murine Malaria Infection

CD4+ αβ T-cells are key mediators of the immune response to a first Plasmodium infection, undergoing extensive activation and splenic expansion during the acute phase of an infection. However, the clonality and clonal composition of this expansion has not previously been described. Using a comparative infection model, we sequenced the splenic CD4+ T-cell receptor repertoires generated over the time-course of a Plasmodium chabaudi infection. We show through repeat replicate experiments, single-cell RNA-seq, and analyses of independent RNA-seq data, that following a first infection - within a highly polyclonal expansion - T-effector repertoires are consistently dominated by TRBV3 gene usage. Clustering by sequence similarity, we find the same dominant clonal signature is expanded across replicates in the acute phase of an infection, revealing a conserved pathogen-specific T-cell response that is consistently a hallmark of a first infection, but not expanded upon re-challenge. Determining the host or parasite factors driving this conserved response may uncover novel immune targets for malaria therapeutic purposes.

Predicting Cross-Reactivity and Antigen Specificity of T Cell Receptors

Adaptive immune recognition is mediated by specific interactions between heterodimeric T cell receptors (TCRs) and their cognate peptide-MHC (pMHC) ligands, and the methods to accurately predict TCR:pMHC interaction would have profound clinical, therapeutic and pharmaceutical applications. Herein, we review recent developments in predicting cross-reactivity and antigen specificity of TCR recognition. We discuss current experimental and computational approaches to investigate cross-reactivity and antigen-specificity of TCRs and highlight how integrating kinetic, biophysical and structural features may offer valuable insights in modeling immunogenicity. We further underscore the close inter-relationship of these two interconnected notions and the need to investigate each in the light of the other for a better understanding of T cell responsiveness for the effective clinical applications.

T cell receptor sequence clustering and antigen specificity

There has been increasing interest in the role of T cells and their involvement in cancer, autoimmune and infectious diseases. However, the nature of T cell receptor (TCR) epitope recognition at a repertoire level is not yet fully understood. Due to technological advances a plethora of TCR sequences from a variety of disease and treatment settings has become readily available. Current efforts in TCR specificity analysis focus on identifying characteristics in immune repertoires which can explain or predict disease outcome or progression, or can be used to monitor the efficacy of disease therapy. In this context, clustering of TCRs by sequence to reflect biological similarity, and especially to reflect antigen specificity have become of paramount importance. We review the main TCR sequence clustering methods and the different similarity measures they use, and discuss their performance and possible improvement. We aim to provide guidance for non-specialists who wish to use TCR repertoire sequencing for disease tracking, patient stratification or therapy prediction, and to provide a starting point for those aiming to develop novel techniques for TCR annotation through clustering.

Methods for sequence and structural analysis of B and T cell receptor repertoires

B cell receptors (BCRs) and T cell receptors (TCRs) make up an essential network of defense molecules that, collectively, can distinguish self from non-self and facilitate destruction of antigen-bearing cells such as pathogens or tumors. The analysis of BCR and TCR repertoires plays an important role in both basic immunology as well as in biotechnology. Because the repertoires are highly diverse, specialized software methods are needed to extract meaningful information from BCR and TCR sequence data. Here, we review recent developments in bioinformatics tools for analysis of BCR and TCR repertoires, with an emphasis on those that incorporate structural features. After describing the recent sequencing technologies for immune receptor repertoires, we survey structural modeling methods for BCR and TCRs, along with methods for clustering such models. We review downstream analyses, including BCR and TCR epitope prediction, antibody-antigen docking and TCR-peptide-MHC Modeling. We also briefly discuss molecular dynamics in this context.

On the viability of unsupervised T-cell receptor sequence clustering for epitope preference

MOTIVATION

The T-cell receptor (TCR) is responsible for recognizing epitopes presented on cell surfaces. Linking TCR sequences to their ability to target specific epitopes is currently an unsolved problem, yet one of great interest. Indeed, it is currently unknown how dissimilar TCR sequences can be before they no longer bind the same epitope. This question is confounded by the fact that there are many ways to define the similarity between two TCR sequences. Here we investigate both issues in the context of TCR sequence unsupervised clustering.

RESULTS

We provide an overview of the performance of various distance metrics on two large independent datasets with 412 and 2835 TCR sequences respectively. Our results confirm the presence of structural distinct TCR groups that target identical epitopes. In addition, we put forward several recommendations to perform unsupervised T-cell receptor sequence clustering.

AVAILABILITY AND IMPLEMENTATION

Source code implemented in Python 3 available at https://github.com/pmeysman/TCRclusteringPaper.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

The Diagnostic, Prognostic, and Therapeutic Potential of Adaptive Immune Receptor Repertoire Profiling in Cancer

Lymphocytes play a critical role in antitumor immune responses. They are directly targeted by some therapies, and the composition and spatial organization of intratumor T-cell populations is prognostic in some cancer types. A better understanding of lymphocyte population dynamics over the course of disease and in response to therapy is urgently needed to guide therapy decisions and to develop new therapy targets. Deep sequencing of the repertoire of antigen receptor-encoding genes expressed in a lymphocyte population has become a widely used approach for profiling the population's immune status. Lymphocyte antigen receptor repertoire deep sequencing data can be used to assess the clonal richness and diversity of lymphocyte populations; to track clone members over time, between tissues, and across lymphocyte subsets; to detect clonal expansion; and to detect the recruitment of new clones into a tissue. Repertoire sequencing is thus a critical complement to other methods of lymphocyte and immune profiling in cancer. This review describes the current state of knowledge based on repertoire sequencing studies conducted on human cancer patients, with a focus on studies of the T-cell receptor beta chain locus. The review then outlines important questions left unanswered and suggests future directions for the field.

Identifying and Tracking Low-Frequency Virus-Specific TCR Clonotypes Using High-Throughput Sequencing

Tracking antigen-specific T cell responses over time within individuals is difficult because of lack of knowledge of antigen-specific TCR sequences, limitations in sample size, and assay sensitivities. We hypothesized that analyses of high-throughput sequencing of TCR clonotypes could provide functional readouts of individuals' immunological histories. Using high-throughput TCR sequencing, we develop a database of TCRβ sequences from large cohorts of mice before (naive) and after smallpox vaccination. We computationally identify 315 vaccine-associated TCR sequences (VATS) that are used to train a diagnostic classifier that distinguishes naive from vaccinated samples in mice up to 9 months post-vaccination with >99% accuracy. We determine that the VATS library contains virus-responsive TCRs by in vitro expansion assays and virus-specific tetramer sorting. These data outline a platform for advancing our capabilities to identify pathogen-specific TCR sequences, which can be used to identify and quantitate low-frequency pathogen-specific TCR sequences in circulation over time with exceptional sensitivity.

sumrep: A Summary Statistic Framework for Immune Receptor Repertoire Comparison and Model Validation

The adaptive immune system generates an incredible diversity of antigen receptors for B and T cells to keep dangerous pathogens at bay. The DNA sequences coding for these receptors arise by a complex recombination process followed by a series of productivity-based filters, as well as affinity maturation for B cells, giving considerable diversity to the circulating pool of receptor sequences. Although these datasets hold considerable promise for medical and public health applications, the complex structure of the resulting adaptive immune receptor repertoire sequencing (AIRR-seq) datasets makes analysis difficult. In this paper we introduce sumrep, an R package that efficiently performs a wide variety of repertoire summaries and comparisons, and show how sumrep can be used to perform model validation. We find that summaries vary in their ability to differentiate between datasets, although many are able to distinguish between covariates such as donor, timepoint, and cell type for BCR and TCR repertoires. We show that deletion and insertion lengths resulting from V(D)J recombination tend to be more discriminative characterizations of a repertoire than summaries that describe the amino acid composition of the CDR3 region. We also find that state-of-the-art generative models excel at recapitulating gene usage and recombination statistics in a given experimental repertoire, but struggle to capture many physiochemical properties of real repertoires.

Deep sequencing of B cell receptor repertoire

Immune repertoire is a collection of enormously diverse adaptive immune cells within an individual. As the repertoire shapes and represents immunological conditions, identification of clones and characterization of diversity are critical for understanding how to protect ourselves against various illness such as infectious diseases and cancers. Over the past several years, fast growing technologies for high throughput sequencing have facilitated rapid advancement of repertoire research, enabling us to observe the diversity of repertoire at an unprecedented level. Here, we focus on B cell receptor (BCR) repertoire and review approaches to B cell isolation and sequencing library construction. These experiments should be carefully designed according to BCR regions to be interrogated, such as heavy chain full length, complementarity determining regions, and isotypes. We also highlight preprocessing steps to remove sequencing and PCR errors with unique molecular index and bioinformatics techniques. Due to the nature of massive sequence variation in BCR, caution is warranted when interpreting repertoire diversity from error-prone sequencing data. Furthermore, we provide a summary of statistical frameworks and bioinformatics tools for clonal evolution and diversity. Finally, we discuss limitations of current BCR-seq technologies and future perspectives on advances in repertoire sequencing.

Recent progress in the analysis of αβT cell and B cell receptor repertoires

T cell receptors (TCRs) and B cell receptors (BCRs) are vertebrate evolution's best answer to the threat of microbial pathogens that can evolve much faster than ourselves. These antigen receptors are generated during T cell or B cell development by combinatorial rearrangement of germline genome V, D and J gene segments, and with junctional residues capable of enormous diversity. For decades the complexity of these receptor repertoires has limited their analysis, but advances in DNA sequencing technology and an array of complementary tools have now made their study much more tractable, filling a major gap in our ability to understand immunology as a system. Here, we summarize the recent approaches and discoveries that are enabling these advances, with some suggestions as to what may lie ahead.

Capturing the differences between humoral immunity in the normal and tumor environments from repertoire-seq of B-cell receptors using supervised machine learning

BACKGROUND

The recent success of immunotherapy in treating tumors has attracted increasing interest in research related to the adaptive immune system in the tumor microenvironment. Recent advances in next-generation sequencing technology enabled the sequencing of whole T-cell receptors (TCRs) and B-cell receptors (BCRs)/immunoglobulins (Igs) in the tumor microenvironment. Since BCRs/Igs in tumor tissues have high affinities for tumor-specific antigens, the patterns of their amino acid sequences and other sequence-independent features such as the number of somatic hypermutations (SHMs) may differ between the normal and tumor microenvironments. However, given the high diversity of BCRs/Igs and the rarity of recurrent sequences among individuals, it is far more difficult to capture such differences in BCR/Ig sequences than in TCR sequences. The aim of this study was to explore the possibility of discriminating BCRs/Igs in tumor and in normal tissues, by capturing these differences using supervised machine learning methods applied to RNA sequences of BCRs/Igs.

RESULTS

RNA sequences of BCRs/Igs were obtained from matched normal and tumor specimens from 90 gastric cancer patients. BCR/Ig-features obtained in Rep-Seq were used to classify individual BCR/Ig sequences into normal or tumor classes. Different machine learning models using various features were constructed as well as gradient boosting machine (GBM) classifier combining these models. The results demonstrated that BCR/Ig sequences between normal and tumor microenvironments exhibit their differences. Next, by using a GBM trained to classify individual BCR/Ig sequences, we tried to classify sets of BCR/Ig sequences into normal or tumor classes. As a result, an area under the curve (AUC) value of 0.826 was achieved, suggesting that BCR/Ig repertoires have distinct sequence-level features in normal and tumor tissues.

CONCLUSIONS

To the best of our knowledge, this is the first study to show that BCR/Ig sequences derived from tumor and normal tissues have globally distinct patterns, and that these tissues can be effectively differentiated using BCR/Ig repertoires.

Checkpoint blockade immunotherapy enhances the frequency and effector function of murine tumor-infiltrating T cells but does not alter TCRβ diversity

Checkpoint blockade immunotherapy is now a first-line treatment option for patients with melanoma. Despite achieving objective responses in about half of patients, the exact immune mechanisms elicited and those required for therapeutic success have not been clearly identified. Insight into these mechanisms is key for improving outcomes in a broader range of cancer patients. We used a murine melanoma model to track responses by different subsets of tumor-infiltrating lymphocytes (TIL) during checkpoint blockade immunotherapy. Tumors from treated mice had increased frequencies of both CD4⁺ and CD8⁺ T cells, which also showed evidence of functional reinvigoration and elevated effector cytokine production after immunotherapy. We predicted that increased T cell numbers and function within tumors reflected either infiltration by new T cells or clonal expansion by a few high-affinity tumor-reactive T cells. To address this, we compared TIL diversity before and after immunotherapy by sequencing the complementarity determining region 3 (CDR3) of all T cell receptor beta (TCRβ) genes. While checkpoint blockade effectively slowed tumor progression and increased T cell frequencies, the diversity of intratumoral T cells remained stable. This was true when analyzing total T cells and when focusing on smaller subsets of effector CD4⁺ and CD8⁺ TIL as well as regulatory T cells. Our study suggests that checkpoint blockade immunotherapy does not broaden the T cell repertoire within murine melanoma tumors, but rather expands existing T cell populations and enhances effector capabilities.

Using T Cell Receptor Repertoires to Understand the Principles of Adaptive Immune Recognition

Adaptive immune recognition is mediated by antigen receptors on B and T cells generated by somatic recombination during lineage development. The high level of diversity resulting from this process posed technical limitations that previously limited the comprehensive analysis of adaptive immune recognition. Advances over the last ten years have produced data and approaches allowing insights into how T cells develop, evolutionary signatures of recombination and selection, and the features of T cell receptors that mediate epitope-specific binding and T cell activation. The size and complexity of these data have necessitated the generation of novel computational and analytical approaches, which are transforming how T cell immunology is conducted. Here we review the development and application of novel biological, theoretical, and computational methods for understanding T cell recognition and discuss the potential for improved models of receptor:antigen interactions.

Biophysicochemical Motifs in T-cell Receptor Sequences Distinguish Repertoires from Tumor-Infiltrating Lymphocyte and Adjacent Healthy Tissue

Immune repertoire deep sequencing allows comprehensive characterization of antigen receptor-encoding genes in a lymphocyte population. We hypothesized that this method could enable a novel approach to diagnose disease by identifying antigen receptor sequence patterns associated with clinical phenotypes. In this study, we developed statistical classifiers of T-cell receptor (TCR) repertoires that distinguish tumor tissue from patient-matched healthy tissue of the same organ. The basis of both classifiers was a biophysicochemical motif in the complementarity determining region 3 (CDR3) of TCRβ chains. To develop each classifier, we extracted 4-mers from every TCRβ CDR3 and represented each 4-mer using biophysicochemical features of its amino acid sequence combined with quantification of 4-mer (or receptor) abundance. This representation was scored using a logistic regression model. Unlike typical logistic regression, the classifier is fitted and validated under the requirement that at least 1 positively labeled 4-mer appears in every tumor repertoire and no positively labeled 4-mers appear in healthy tissue repertoires. We applied our method to publicly available data in which tumor and adjacent healthy tissue were collected from each patient. Using a patient-holdout cross-validation, our method achieved classification accuracy of 93% and 94% for colorectal and breast cancer, respectively. The parameter values for each classifier revealed distinct biophysicochemical properties for tumor-associated 4-mers within each cancer type. We propose that such motifs might be used to develop novel immune-based cancer screening assays. SIGNIFICANCE: This study presents a novel computational approach to identify T-cell repertoire differences between normal and tumor tissue.See related commentary by Zoete and Coukos, p. 1299.

The expanding role of systems immunology in decoding the T cell receptor repertoire

T cells play a crucial role in the immune system's defense against many infectious diseases, including persistent infections for which no effective vaccines currently exist. The T cell component of the adaptive immune system is highly complex involving a constantly evolving landscape of various inter-related T cell populations. These T cell populations are characterized by their phenotypic and functional properties as well as the collection, or repertoire, of T cell receptors (TCR) that mediate T cell recognition of antigenic peptides derived from pathogens. Understanding the various processes and factors that impact the development and evolution of the broader T cell repertoire available to recognize and respond to pathogens and the characteristics of antigen-experienced T cell repertoires associated with effective immune control of pathogens is critical to the rational design of T cell-based vaccines and therapies. In this article we discuss, using examples of recent research, the promise that systems immunology approaches, involving quantitative analysis and mathematical and computational modeling of immunological data, hold for decoding the complex TCR repertoire system in the current era of advancing technologies.

Human T cell receptor occurrence patterns encode immune history, genetic background, and receptor specificity

The T cell receptor (TCR) repertoire encodes immune exposure history through the dynamic formation of immunological memory. Statistical analysis of repertoire sequencing data has the potential to decode disease associations from large cohorts with measured phenotypes. However, the repertoire perturbation induced by a given immunological challenge is conditioned on genetic background via major histocompatibility complex (MHC) polymorphism. We explore associations between MHC alleles, immune exposures, and shared TCRs in a large human cohort. Using a previously published repertoire sequencing dataset augmented with high-resolution MHC genotyping, our analysis reveals rich structure: striking imprints of common pathogens, clusters of co-occurring TCRs that may represent markers of shared immune exposures, and substantial variations in TCR-MHC association strength across MHC loci. Guided by atomic contacts in solved TCR:peptide-MHC structures, we identify sequence covariation between TCR and MHC. These insights and our analysis framework lay the groundwork for further explorations into TCR diversity.

High-throughput sequencing of the T-cell receptor repertoire: pitfalls and opportunities

T-cell specificity is determined by the T-cell receptor, a heterodimeric protein coded for by an extremely diverse set of genes produced by imprecise somatic gene recombination. Massively parallel high-throughput sequencing allows millions of different T-cell receptor genes to be characterized from a single sample of blood or tissue. However, the extraordinary heterogeneity of the immune repertoire poses significant challenges for subsequent analysis of the data. We outline the major steps in processing of repertoire data, considering low-level processing of raw sequence files and high-level algorithms, which seek to extract biological or pathological information. The latest generation of bioinformatics tools allows millions of DNA sequences to be accurately and rapidly assigned to their respective variable V and J gene segments, and to reconstruct an almost error-free representation of the non-templated additions and deletions that occur. High-level processing can measure the diversity of the repertoire in different samples, quantify V and J usage and identify private and public T-cell receptors. Finally, we discuss the major challenge of linking T-cell receptor sequence to function, and specifically to antigen recognition. Sophisticated machine learning algorithms are being developed that can combine the paradoxical degeneracy and cross-reactivity of individual T-cell receptors with the specificity of the overall T-cell immune response. Computational analysis will provide the key to unlock the potential of the T-cell receptor repertoire to give insight into the fundamental biology of the adaptive immune system and to provide powerful biomarkers of disease.

The Bayesian optimist's guide to adaptive immune receptor repertoire analysis

Probabilistic modeling is fundamental to the statistical analysis of complex data. In addition to forming a coherent description of the data-generating process, probabilistic models enable parameter inference about given datasets. This procedure is well developed in the Bayesian perspective, in which one infers probability distributions describing to what extent various possible parameters agree with the data. In this paper, we motivate and review probabilistic modeling for adaptive immune receptor repertoire data then describe progress and prospects for future work, from germline haplotyping to adaptive immune system deployment across tissues. The relevant quantities in immune sequence analysis include not only continuous parameters such as gene use frequency but also discrete objects such as B-cell clusters and lineages. Throughout this review, we unravel the many opportunities for probabilistic modeling in adaptive immune receptor analysis, including settings for which the Bayesian approach holds substantial promise (especially if one is optimistic about new computational methods). From our perspective, the greatest prospects for progress in probabilistic modeling for repertoires concern ancestral sequence estimation for B-cell receptor lineages, including uncertainty from germline genotype, rearrangement, and lineage development.

CD154 Costimulation Shifts the Local T-Cell Receptor Repertoire Not Only During Thymic Selection but Also During Peripheral T-Dependent Humoral Immune Responses

CD154 is a transmembrane cytokine expressed transiently on activated CD4 T cells upon T-cell receptor (TCR) stimulation that interacts with CD40 on antigen-presenting cells. The signaling via CD154:CD40 is essential for B-cell maturation and germinal center formation and also for the final differentiation of CD4 T cells during T-dependent humoral immune responses. Recent data demonstrate that CD154 is critically involved in the selection of T-cell clones during the negative selection process in the thymus. Whether CD154 signaling influences the TCR repertoire during peripheral T-dependent humoral immune responses has not yet been elucidated. To find out, we used CD154-deficient mice and assessed the global TCRβ repertoire in T-cell zones (TCZ) of spleens by high-throughput sequencing after induction of a Th2 response to the multiepitopic antigen sheep red blood cells. Qualitative and quantitative comparison of the splenic TCZ-specific TCRβ repertoires revealed that CD154 deficiency shifts the distribution of Vβ-Jβ genes after antigen exposure. This data led to the conclusion that costimulation via CD154:CD40 during the interaction of T cells with CD40-matured B cells contributes to the recruitment of T-cell clones into the immune response and thereby shapes the peripheral TCR repertoire.

RepSeq Data Representativeness and Robustness Assessment by Shannon Entropy

High-throughput sequencing (HTS) has the potential to decipher the diversity of T cell repertoires and their dynamics during immune responses. Applied to T cell subsets such as T effector and T regulatory cells, it should help identify novel biomarkers of diseases. However, given the extreme diversity of TCR repertoires, understanding how the sequencing conditions, including cell numbers, biological and technical sampling and sequencing depth, impact the experimental outcome is critical to proper use of these data. Here, we assessed the representativeness and robustness of TCR repertoire diversity assessment according to experimental conditions. By comparative analyses of experimental datasets and computer simulations, we found that (i) for small samples, the number of clonotypes recovered is often higher than the number of cells per sample, even after removing the singletons; (ii) high-sequencing depth for small samples alters the clonotype distributions, which can be corrected by filtering the datasets using Shannon entropy as a threshold; and (iii) a single sequencing run at high depth does not ensure a good coverage of the clonotype richness in highly polyclonal populations, which can be better covered using multiple sequencing. Altogether, our results warrant better understanding and awareness of the limitation of TCR diversity analyses by HTS and justify the development of novel computational tools for improved modeling of the highly complex nature of TCR repertoires.