Estimating T-cell repertoire diversity: limitations of classical estimators and a new approach

Deep sequencing of the T cell receptor β repertoire reveals signature patterns and clonal drift in atherosclerotic plaques and patients

The T cell receptor (TCR) β repertoire directly reflects the status of T cell function. Meanwhile, the immune/inflammatory responses regulated by T cells are the critical determinants of atherosclerosis development. However, due to technical limitations, the composition and molecular characteristics of the TCR repertoire in atherosclerotic patients have not been fully elucidated. In the present study, we use powerful immune repertoire sequencing technology to study this issue. Results show that the utilization of V and/or J genes and the diversity of TCRβ repertoire in atherosclerotic plaques are significantly reduced compared to those in the peripheral blood of normal subjects and atherosclerotic patients. The frequencies of the common T cell clones with certain lengths of the complement determining region 3 regions are notably different among all groups. The high-frequency common clones are also increased in the atherosclerotic plaques compared to that in the other two groups. The expansion of several T cell clonotypes (V29-1J2-1, V20-1J1-6, V6-3J2-7 and V11-2J2-2) is validated in atherosclerotic patients. In short, this study reveals that the diversity of TCR β repertoire significantly decreases in atherosclerotic plaques, probably because of the reduced utilization of VJ genes and marked expansion of some T cell subclones. It provides the basis for understanding the roles of T lymphocytes in the pathogenesis of atherosclerosis.

Quantification of Inter-Sample Differences in T-Cell Receptor Repertoires Using Sequence-Based Information

Inter-sample comparisons of T-cell receptor (TCR) repertoires are crucial for gaining a better understanding of the immunological states determined by different collections of T cells from different donor sites, cell types, and genetic and pathological backgrounds. For quantitative comparison, most previous studies utilized conventional methods in ecology, which focus on TCR sequences that overlap between pairwise samples. Some recent studies attempted another approach that is categorized into Poisson abundance models using the abundance distribution of observed TCR sequences. However, these methods ignore the details of the measured sequences and are consequently unable to identify sub-repertoires that might have important contributions to the observed inter-sample differences. Moreover, the sparsity of sequence data due to the huge diversity of repertoires hampers the performance of these methods, especially when few overlapping sequences exist. In this paper, we propose a new approach for REpertoire COmparison in Low Dimensions (RECOLD) based on TCR sequence information, which can estimate the low-dimensional structure by embedding the pairwise sequence dissimilarities in high-dimensional sequence space. The inter-sample differences between repertoires are then quantified by information-theoretic measures among the distributions of data estimated in the embedded space. Using datasets of mouse and human TCR repertoires, we demonstrate that RECOLD can accurately identify the inter-sample hierarchical structures, which have a good correspondence with our intuitive understanding about sample conditions. Moreover, for the dataset of transgenic mice that have strong restrictions on the diversity of their repertoires, our estimated inter-sample structure was consistent with the structure estimated by previous methods based on abundance or overlapping sequence information. For the dataset of human healthy donors and Sézary syndrome patients, our method also showed robust estimation performance even under the condition of high sparsity in TCR sequences, while previous studies failed to estimate the structure. In addition, we identified the sequences that contribute to the pairwise-sample differences between the repertoires with the different genetic backgrounds of mice. Such identification of the sequences contributing to variation in immune cell repertoires may provide substantial insight for the development of new immunotherapies and vaccines.

Local bone-marrow exposure: how to interpret the data on stable chromosome aberrations in circulating lymphocytes? (some comments on the use of FISH method for dose reconstruction for Techa riverside Residents)

The method of fluorescence in situ hybridization (FISH) applied to peripheral blood T lymphocytes is used for retrospective dose estimation, and the results obtained from the analysis of stable chromosomal aberrations are usually interpreted as a dose accumulated in the red bone marrow (RBM). However, after local internal exposure of the RBM, doses derived from FISH were found to be lower than those derived from direct measurements of radionuclides accumulated in the bodies of exposed persons. These results were obtained for people residing near the Techa River contaminated by ^89,90Sr (beta-emitters) in 1949-1956 (Chelyabinsk Oblast, Russia). A new analysis has been performed of the combined results of FISH studies (n = 178) undertaken during 1994-2012 for persons living on the Techa Riverside. Analysis confirms the lower slope of the translocation yield per Gy (8.0 ± 0.7 × 10^-3) for Techa residents in comparison with FISH data for donors with external exposures (11.6 ± 1.6 × 10^-3, Tawn et al., Radiat Res 184(3):296-303, 2015). It was suggested that some portion of T cells remained unexposed, because they represented the descendants of T cell progenitors, which had migrated to the thymus before the start of ^89,90Sr intakes. To clarify this problem, the dynamics of T-cell Genera (TG), combining all descendants of specific T-cell progenitor reaching the thymus, was considered. Rates of TGs produced by RBM over different age periods of human life were estimated with the use of the mathematic model of T-cell homeostasis (Bains, Mathematical modeling of T-cell homeostasis. A thesis submitted for the degree of Doctor of Philosophy of the University College London. http://discovery.ucl.ac.uk/20159/1/20159.pdf , 2010). The rate of TG loss during the lifetime was assumed to be very small in comparison with production rate. The recirculation of mature T lymphocytes in contaminated RBM was taken into account. According to our model estimates, at the time of blood sampling, the fraction of exposed T lymphocytes (whose progenitors were irradiated) ranged from 20 to 80% depending on the donors' age at the start of exposure to ^89,90Sr. Dose to T lymphocytes, estimated from FISH studies, should be about 0.6-0.9 of RBM dose for residents of the upper Techa region and about 0.4-0.8 in the middle Techa region. Our results could explain the lower value of translocation yield per Gy obtained for Techa residents. The approaches for further model improvement and validation are discussed in this paper.

Quantitative Characterization of the T Cell Receptor Repertoire of Naïve and Memory Subsets Using an Integrated Experimental and Computational Pipeline Which Is Robust, Economical, and Versatile

The T cell receptor (TCR) repertoire can provide a personalized biomarker for infectious and non-infectious diseases. We describe a protocol for amplifying, sequencing, and analyzing TCRs which is robust, sensitive, and versatile. The key experimental step is ligation of a single-stranded oligonucleotide to the 3' end of the TCR cDNA. This allows amplification of all possible rearrangements using a single set of primers per locus. It also introduces a unique molecular identifier to label each starting cDNA molecule. This molecular identifier is used to correct for sequence errors and for effects of differential PCR amplification efficiency, thus producing more accurate measures of the true TCR frequency within the sample. This integrated experimental and computational pipeline is applied to the analysis of human memory and naive subpopulations, and results in consistent measures of diversity and inequality. After error correction, the distribution of TCR sequence abundance in all subpopulations followed a power law over a wide range of values. The power law exponent differed between naïve and memory populations, but was consistent between individuals. The integrated experimental and analysis pipeline we describe is appropriate to studies of T cell responses in a broad range of physiological and pathological contexts.

Successful and Maladaptive T Cell Aging

Throughout life, the T cell system adapts to shifting resources and demands, resulting in a fundamentally restructured immune system in older individuals. Here we review the cellular and molecular features of an aged immune system and discuss the trade-offs inherent to these adaptive mechanisms. Processes include homeostatic proliferation that maintains compartment size at the expense of partial loss in stemness and incomplete differentiation and the activation of negative regulatory programs, which constrain effector T cell expansion and prevent increasing oligoclonality but also interfere with memory cell generation. We propose that immune failure occurs when adaptive strategies developed by the aging T cell system fail and also discuss how, in some settings, the programs associated with T cell aging culminates in a maladaptive response that directly contributes to chronic inflammatory disease.

The impact of HTLV-1 on the cellular genome

Human T-lymphotropic virus type-1 (HTLV-1) is the causative agent of adult T-cell leukaemia/lymphoma (ATL), an aggressive CD4+ T-cell malignancy. The mechanisms of leukaemogenesis in ATL are incompletely understood. Insertional mutagenesis has not previously been thought to contribute to the pathogenesis of ATL. However, the recent discovery that HTLV-1 binds the key chromatin architectural protein CTCF raises the hypothesis that HTLV-1 deregulates host gene expression by causing abnormal chromatin looping, bringing the strong HTLV-1 promoter-enhancer near to host genes that lie up to 2Mb from the integrated provirus. Here we review current opinion on the mechanisms of oncogenesis in ATL, with particular emphasis on the local and distant impact of HTLV-1 on the structure and expression of the host genome.

Overview of methodologies for T-cell receptor repertoire analysis

Background

The T-cell receptor (TCR), located on the surface of T cells, is responsible for the recognition of the antigen-major histocompatibility complex, leading to the initiation of an inflammatory response. Analysing the TCR repertoire may help to gain a better understanding of the immune system features and of the aetiology and progression of diseases, in particular those with unknown antigenic triggers. The extreme diversity of the TCR repertoire represents a major analytical challenge; this has led to the development of specialized methods which aim to characterize the TCR repertoire in-depth. Currently, next generation sequencing based technologies are most widely employed for the high-throughput analysis of the immune cell repertoire.

Results

Here, we report on the latest methodological advancements in the field by describing and comparing the available tools; from the choice of the starting material and library preparation method, to the sequencing technologies and data analysis.

Finally, we provide a practical example and our own experience by reporting some exemplary results from a small internal benchmark study, where current approaches from the literature and the market are employed and compared.

Conclusions

Several valid methods for clonotype identification and TCR repertoire analysis exist, however, a gold standard method for the field has not yet been identified. Depending on the purpose of the scientific study, some approaches may be more suitable than others. Finally, due to possible method specific biases, scientists must be careful when comparing results obtained using different methods.

Electronic supplementary material

The online version of this article (doi:10.1186/s12896-017-0379-9) contains supplementary material, which is available to authorized users.

Pattern formation in a nonlocal mathematical model for the multiple roles of the TGF-β pathway in tumour dynamics

The growth and invasion of cancer cells are very complex processes, which can be regulated by the cross-talk between various signalling pathways, or by single signalling pathways that can control multiple aspects of cell behaviour. TGF-β is one of the most investigated signalling pathways in oncology, since it can regulate multiple aspects of cell behaviour: cell proliferation and apoptosis, cell-cell adhesion and epithelial-to-mesenchimal transition via loss of cell adhesion. In this study, we use a mathematical modelling approach to investigate the complex roles of TGF-β signalling pathways on the inhibition and growth of tumours, as well as on the epithelial-to-mesenchimal transition involved in the metastasis of tumour cells. We show that the nonlocal mathematical model derived here to describe repulsive and adhesive cell-cell interactions can explain the formation of new tumour cell aggregations at positions in space that are further away from the main aggregation. Moreover, we show that the increase in cell-cell adhesion leads to fewer but larger aggregations, and the increase in TGF-β molecules - whose late-stage effect is to decrease cell adhesion - leads to many small cellular aggregations. Finally, we perform a sensitivity analysis on some parameters associated with TGF-β dynamics, and use it to investigate the relation between the tumour size and its metastatic spread.

Sequence and Structural Analyses Reveal Distinct and Highly Diverse Human CD8+ TCR Repertoires to Immunodominant Viral Antigens

A diverse T cell receptor (TCR) repertoire is essential for controlling viral infections. However, information about TCR repertoires to defined viral antigens is limited. We performed a comprehensive analysis of CD8⁺ TCR repertoires for two dominant viral epitopes: pp65_495-503 (NLV) of cytomegalovirus and M1_58-66 (GIL) of influenza A virus. The highly individualized repertoires (87-5,533 α or β clonotypes per subject) comprised thousands of unique TCRα and TCRβ sequences and dozens of distinct complementary determining region (CDR)3α and CDR3β motifs. However, diversity is effectively restricted by preferential V-J combinations, CDR3 lengths, and CDR3α/CDR3β pairings. Structures of two GIL-specific TCRs bound to GIL-HLA-A2 provided a potential explanation for the lower diversity of GIL-specific versus NLV-specific repertoires. These anti-viral TCRs occupied up to 3.4% of the CD8⁺ TCRβ repertoire, ensuring broad T cell responses to single epitopes. Our portrait of two anti-viral TCR repertoires may inform the development of predictors of immune protection.

Methodologic Considerations in the Application of Next-Generation Sequencing of Human TRB Repertoires for Clinical Use

Next-generation sequencing (NGS) of immune receptors has become a standard tool to assess minimal residual disease (MRD) in patients treated for lymphoid malignancy, and it is being used to study the T-cell repertoire in many clinical settings. To better understanding the potential clinical utility and limitations of this application outside of MRD, we developed a BIOMED-2 primer-based NGS method and characterized its performance in controls and patients with graft-versus-host disease (GVHD) after allogeneic hematopoietic transplant. For controls and patients with GVHD, replicate sequencing of the same T-cell receptor β (TRB) libraries was highly reproducible. Higher variability was observed in sequencing of different TRB libraries made from the same DNA stock. Variability was increased in patients with GVHD compared with controls; patients with GVHD also had lower diversity than controls. In the T-cell repertoire of a healthy person, approximately 99.6% of the CDR3 clones were in low abundance, with frequency <10^-3. A single library could identify >93% of the clones with frequency ≥10^-3 in the repertoire. Sequencing in duplicate increased the average detection rate to >97%. This work demonstrates that NGS reliably and robustly characterizes TRB populations in healthy individuals and patients with GVHD with frequency ≥10^-3 and provides a methodologic framework for applying NGS immune repertoire methods to clinical testing applications beyond MRD.

On the organization of human T-cell receptor loci: log-periodic distribution of T-cell receptor gene segments

The human T-cell repertoire is complex and is generated by the rearrangement of variable (V), diversity (D) and joining (J) segments on the T-cell receptor (TCR) loci. The T-cell repertoire demonstrates self-similarity in terms clonal frequencies when defined by V, D and J gene segment usage; therefore to determine whether the structural ordering of these gene segments on the TCR loci contributes to the observed clonal frequencies, the TCR loci were examined for self-similarity and periodicity in terms of gene segment organization. Logarithmic transformation of numeric sequence order demonstrated that the V and J gene segments for both T-cell receptor α (TRA) and β (TRB) loci are arranged in a self-similar manner when the spacing between the adjacent segments was considered as a function of the size of the neighbouring gene segment, with an average fractal dimension of approximately 1.5. Accounting for the gene segments occurring on helical DNA molecules with a logarithmic distribution, sine and cosine functions of the log-transformed angular coordinates of the start and stop nucleotides of successive TCR gene segments showed an ordered progression from the 5' to the 3' end of the locus, supporting a log-periodic organization. T-cell clonal frequency estimates, based on V and J segment usage, from normal stem cell donors were plotted against the V and J segment on TRB locus and demonstrated a periodic distribution. We hypothesize that this quasi-periodic variation in gene-segment representation in the T-cell clonal repertoire may be influenced by the location of the gene segments on the periodic-logarithmically scaled TCR loci. Interactions between the two strands of DNA in the double helix may influence the probability of gene segment usage by means of either constructive or destructive interference resulting from the superposition of the two helices.

Mathematical Models for Immunology: Current State of the Art and Future Research Directions

The advances in genetics and biochemistry that have taken place over the last 10 years led to significant advances in experimental and clinical immunology. In turn, this has led to the development of new mathematical models to investigate qualitatively and quantitatively various open questions in immunology. In this study we present a review of some research areas in mathematical immunology that evolved over the last 10 years. To this end, we take a step-by-step approach in discussing a range of models derived to study the dynamics of both the innate and immune responses at the molecular, cellular and tissue scales. To emphasise the use of mathematics in modelling in this area, we also review some of the mathematical tools used to investigate these models. Finally, we discuss some future trends in both experimental immunology and mathematical immunology for the upcoming years.

Single-cell TCRseq: paired recovery of entire T-cell alpha and beta chain transcripts in T-cell receptors from single-cell RNAseq

Accurate characterization of the repertoire of the T-cell receptor (TCR) alpha and beta chains is critical to understanding adaptive immunity. Such characterization has many applications across such fields as vaccine development and response, clone-tracking in cancer, and immunotherapy. Here we present a new methodology called single-cell TCRseq (scTCRseq) for the identification and assembly of full-length rearranged V(D)J T-cell receptor sequences from paired-end single-cell RNA sequencing reads. The method allows accurate identification of the V(D)J rearrangements for each individual T-cell and has the novel ability to recover paired alpha and beta segments. Source code is available at https://github.com/ElementoLab/scTCRseq.

Host Resistance and Immune Aging

Human immune system aging results in impaired responses to pathogens or vaccines. In the innate immune system, which mediates the earliest pro-inflammatory responses to immunologic challenge, processes ranging from Toll-like Receptor function to Neutrophil Extracellular Trap formation are generally diminished in older adults. Dysregulated, enhanced basal inflammation with age reflecting activation by endogenous damage-associated ligands contributes to impaired innate immune responses. In the adaptive immune system, T and B cell subsets and function alter with age. The control of cytomegalovirus infection, particularly in the T lineage, plays a dominant role in the differentiation and diversity of the T cell compartment.

Tracking TCRβ Sequence Clonotype Expansions during Antiviral Therapy Using High-Throughput Sequencing of the Hypervariable Region

To maintain a persistent infection viruses such as hepatitis C virus (HCV) employ a range of mechanisms that subvert protective T cell responses. The suppression of antigen-specific T cell responses by HCV hinders efforts to profile T cell responses during chronic infection and antiviral therapy. Conventional methods of detecting antigen-specific T cells utilize either antigen stimulation (e.g., ELISpot, proliferation assays, cytokine production) or antigen-loaded tetramer staining. This limits the ability to profile T cell responses during chronic infection due to suppressed effector function and the requirement for prior knowledge of antigenic viral peptide sequences. Recently, high-throughput sequencing (HTS) technologies have been developed for the analysis of T cell repertoires. In the present study, we have assessed the feasibility of HTS of the TCRβ complementarity determining region (CDR)3 to track T cell expansions in an antigen-independent manner. Using sequential blood samples from HCV-infected individuals undergoing antiviral therapy, we were able to measure the population frequencies of >35,000 TCRβ sequence clonotypes in each individual over the course of 12 weeks. TRBV/TRBJ gene segment usage varied markedly between individuals but remained relatively constant within individuals across the course of therapy. Despite this stable TRBV/TRBJ gene segment usage, a number of TCRβ sequence clonotypes showed dramatic changes in read frequency. These changes could not be linked to therapy outcomes in the present study; however, the TCRβ CDR3 sequences with the largest fold changes did include sequences with identical TRBV/TRBJ gene segment usage and high junction region homology to previously published CDR3 sequences from HCV-specific T cells targeting the HLA-B*0801-restricted ¹³⁹⁵HSKKKCDEL¹⁴⁰³ and HLA-A*0101-restricted ¹⁴³⁵ATDALMTGY¹⁴⁴³ epitopes. The pipeline developed in this proof of concept study provides a platform for the design of future experiments to accurately address the question of whether T cell responses contribute to SVR upon antiviral therapy. This pipeline represents a novel technique to analyze T cell dynamics in situations where conventional antigen-dependent methods are limited due to suppression of T cell functions and highly diverse antigenic sequences.

High-throughput sequencing of immune repertoires in multiple sclerosis

T cells and B cells are crucial in the initiation and maintenance of multiple sclerosis (MS), and the activation of these cells is believed to be mediated through specific recognition of antigens by the T‐ and B‐cell receptors. The antigen receptors are highly polymorphic due to recombination (T‐ and B‐cell receptors) and mutation (B‐cell receptors) of the encoding genes, which can therefore be used as fingerprints to track individual T‐ and B‐cell clones. Such studies can shed light on mechanisms driving the immune responses and provide new insights into the pathogenesis. Here, we summarize studies that have explored the T‐ and B‐cell receptor repertoires using earlier methodological approaches, and we focus on how high‐throughput sequencing has provided new knowledge by surveying the immune repertoires in MS in even greater detail and with unprecedented depth.

High-throughput sequencing insights into T-cell receptor repertoire diversity in aging

Decline in T-cell generation leading to T-cell receptor repertoire contraction is a cornerstone of immune system aging, and consequent disorders. High-throughput sequencing enables in-depth immune repertoire characterization, but blood samples are too small to capture its total diversity. New computational models could enable accurate estimation of this diversity.

In-Depth Assessment of Within-Individual and Inter-Individual Variation in the B Cell Receptor Repertoire

High-throughput sequencing of the B cell receptor (BCR) repertoire can provide rapid characterization of the B cell response in a wide variety of applications in health, after vaccination and in infectious, inflammatory and immune-driven disease, and is starting to yield clinical applications. However, the interpretation of repertoire data is compromised by a lack of studies to assess the intra and inter-individual variation in the BCR repertoire over time in healthy individuals. We applied a standardized isotype-specific BCR repertoire deep sequencing protocol to a single highly sampled participant, and then evaluated the method in 9 further participants to comprehensively describe such variation. We assessed total repertoire metrics of mutation, diversity, VJ gene usage and isotype subclass usage as well as tracking specific BCR sequence clusters. There was good assay reproducibility (both in PCR amplification and biological replicates), but we detected striking fluctuations in the repertoire over time that we hypothesize may be due to subclinical immune activation. Repertoire properties were unique for each individual, which could partly be explained by a decrease in IgG2 with age, and genetic differences at the immunoglobulin locus. There was a small repertoire of public clusters (0.5, 0.3, and 1.4% of total IgA, IgG, and IgM clusters, respectively), which was enriched for expanded clusters containing sequences with suspected specificity toward antigens that should have been historically encountered by all participants through prior immunization or infection. We thus provide baseline BCR repertoire information that can be used to inform future study design, and aid in interpretation of results from these studies. Furthermore, our results indicate that BCR repertoire studies could be used to track changes in the public repertoire in and between populations that might relate to population immunity against infectious diseases, and identify the characteristics of inflammatory and immunological diseases.