Persisting fetal clonotypes influence the structure and overlap of adult human T cell receptor repertoires

Exploring the pre-immune landscape of antigen-specific T cells

Background

Adaptive immune responses to newly encountered pathogens depend on the mobilization of antigen-specific clonotypes from a vastly diverse pool of naive T cells. Using recent advances in immune repertoire sequencing technologies, models of the immune receptor rearrangement process, and a database of annotated T cell receptor (TCR) sequences with known specificities, we explored the baseline frequencies of T cells specific for defined human leukocyte antigen (HLA) class I-restricted epitopes in healthy individuals.

Methods

We used a database of TCR sequences with known antigen specificities and a probabilistic TCR rearrangement model to estimate the baseline frequencies of TCRs specific to distinct antigens epitopespecificT-cells. We verified our estimates using a publicly available collection of TCR repertoires from healthy individuals. We also interrogated a database of immunogenic and non-immunogenic peptides is used to link baseline T-cell frequencies with epitope immunogenicity.

Results

Our findings revealed a high degree of variability in the prevalence of T cells specific for different antigens that could be explained by the physicochemical properties of the corresponding HLA class I-bound peptides. The occurrence of certain rearrangements was influenced by ancestry and HLA class I restriction, and umbilical cord blood samples contained higher frequencies of common pathogen-specific TCRs. We also identified a quantitative link between specific T cell frequencies and the immunogenicity of cognate epitopes presented by defined HLA class I molecules.

Conclusions

Our results suggest that the population frequencies of specific T cells are strikingly non-uniform across epitopes that are known to elicit immune responses. This inference leads to a new definition of epitope immunogenicity based on specific TCR frequencies, which can be estimated with a high degree of accuracy in silico, thereby providing a novel framework to integrate computational and experimental genomics with basic and translational research efforts in the field of T cell immunology.

Electronic supplementary material

The online version of this article (10.1186/s13073-018-0577-7) contains supplementary material, which is available to authorized users.

The Changing Landscape of Naive T Cell Receptor Repertoire With Human Aging

Human aging is associated with a profound loss of thymus productivity, yet naïve T lymphocytes still maintain their numbers by division in the periphery for many years. The extent of such proliferation may depend on the cytokine environment, including IL-7 and T-cell receptor (TCR) “tonic” signaling mediated by self pMHCs recognition. Additionally, intrinsic properties of distinct subpopulations of naïve T cells could influence the overall dynamics of aging-related changes within the naïve T cell compartment. Here, we investigated the differences in the architecture of TCR beta repertoires for naïve CD4, naïve CD8, naïve CD4⁺CD25⁻CD31⁺ (enriched with recent thymic emigrants, RTE), and mature naïve CD4⁺CD25⁻CD31⁻ peripheral blood subsets between young and middle-age/old healthy individuals. In addition to observing the accumulation of clonal expansions (as was shown previously), we reveal several notable changes in the characteristics of T cell repertoire. We observed significant decrease of CDR3 length, NDN insert, and number of non-template added N nucleotides within TCR beta CDR3 with aging, together with a prominent change of physicochemical properties of the central part of CDR3 loop. These changes were similar across CD4, CD8, RTE-enriched, and mature CD4 subsets of naïve T cells, with minimal or no difference observed between the latter two subsets for individuals of the same age group. We also observed an increase in “publicity” (fraction of shared clonotypes) of CD4, but not CD8 naïve T cell repertoires. We propose several explanations for these phenomena built upon previous studies of naïve T-cell homeostasis, and call for further studies of the mechanisms causing the observed changes and of consequences of these changes in respect of the possible holes formed in the landscape of naïve T cell TCR repertoire.

Predicting the spectrum of TCR repertoire sharing with a data-driven model of recombination

Despite the extreme diversity of T-cell repertoires, many identical T-cell receptor (TCR) sequences are found in a large number of individual mice and humans. These widely shared sequences, often referred to as "public," have been suggested to be over-represented due to their potential immune functionality or their ease of generation by V(D)J recombination. Here, we show that even for large cohorts, the observed degree of sharing of TCR sequences between individuals is well predicted by a model accounting for the known quantitative statistical biases in the generation process, together with a simple model of thymic selection. Whether a sequence is shared by many individuals is predicted to depend on the number of queried individuals and the sampling depth, as well as on the sequence itself, in agreement with the data. We introduce the degree of publicness conditional on the queried cohort size and the size of the sampled repertoires. Based on these observations, we propose a public/private sequence classifier, "PUBLIC" (Public Universal Binary Likelihood Inference Classifier), based on the generation probability, which performs very well even for small cohort sizes.

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.

Method for identification of condition-associated public antigen receptor sequences

Diverse repertoires of hypervariable immunoglobulin receptors (TCR and BCR) recognize antigens in the adaptive immune system. The development of immunoglobulin receptor repertoire sequencing methods makes it possible to perform repertoire-wide disease association studies of antigen receptor sequences. We developed a statistical framework for associating receptors to disease from only a small cohort of patients, with no need for a control cohort. Our method successfully identifies previously validated Cytomegalovirus and type one diabetes responsive TCRβ sequences .

Quantitative profiling reveals minor changes of T cell receptor repertoire in response to subunit inactivated influenza vaccine

Vaccination against influenza is widely used to protect against seasonal flu epidemic although its effectiveness is debated. Here we performed deep quantitative T cell receptor repertoire profiling in peripheral blood of a healthy volunteer in response to trivalent subunit influenza vaccine. We did not observe significant rebuilding of peripheral blood T cell receptors composition in response to vaccination. However, we found several clonotypes in memory T cell fraction that were undetectable before the vaccination and had a maximum concentration at day 45 after vaccine administration. These cells were found in lower concentration in the course of repertoire monitoring for two years period. Our observation suggests a potential for recruitment of only a limited number of new T cells after each seasonal influenza vaccination.

High-throughput immune repertoire analysis with IGoR

High-throughput immune repertoire sequencing is promising to lead to new statistical diagnostic tools for medicine and biology. Successful implementations of these methods require a correct characterization, analysis, and interpretation of these data sets. We present IGoR (Inference and Generation Of Repertoires)-a comprehensive tool that takes B or T cell receptor sequence reads and quantitatively characterizes the statistics of receptor generation from both cDNA and gDNA. It probabilistically annotates sequences and its modular structure can be used to investigate models of increasing biological complexity for different organisms. For B cells, IGoR returns the hypermutation statistics, which we use to reveal co-localization of hypermutations along the sequence. We demonstrate that IGoR outperforms existing tools in accuracy and estimate the sample sizes needed for reliable repertoire characterization.

[IGoR: a tool for learning and simulating the random generation of antigen receptors]

Antigen receptors, which form the base of the adaptive immune system, are created stochastically by a DNA editing process called V(D)J recombination. As high-throughput sequencing enables to study the repertoire of these receptors, it is now possible to learn the probabilistic laws of this random process, and to use them to analyse receptors of interest, generate synthetic repertoires to create controls, or aid the identification of receptors that are specific to diseases, with possible applications for medical diagnostics. This article describes how these tasks can be performed using the IGoR software, which can learn statistical models from data, annotate existing sequences, or generate new synthetic ones with the same laws as the recombination process.

Public Clonotypes and Convergent Recombination Characterize the Naïve CD8+ T-Cell Receptor Repertoire of Extremely Preterm Neonates

Respiratory support improvements have aided survival of premature neonates, but infection susceptibility remains a predominant problem. We previously reported that neonatal mice have a rapidly evolving T-cell receptor (TCR) repertoire that impairs CD8⁺ T cell immunity. To understand the impact of prematurity on the human CD8⁺ TCR repertoire, we performed next-generation sequencing of the complementarity-determining region 3 (CDR3) from the rearranged TCR variable beta (Vβ) in sorted, naïve CD8⁺ T cells from extremely preterm neonates (23–27 weeks gestation), term neonates (37–41 weeks gestation), children (16–56 months), and adults (25–50 years old). Strikingly, preterm neonates had an increased frequency of public clonotypes shared between unrelated individuals. Public clonotypes identified in preterm infants were encoded by germline gene sequences, and some of these clonotypes persisted into adulthood. The preterm neonatal naïve CD8⁺ TCR repertoire exhibited convergent recombination, characterized by different nucleotide sequences encoding the same amino acid CDR3 sequence. As determined by Pielou’s evenness and iChao1 metrics, extremely preterm neonates have less clonality, and a much lower bound for the number of unique TCR within an individual preterm neonate, which indicates a less rich and diverse repertoire, as compared to term neonates, children, and adults. This suggests that T cell selection in the preterm neonate may be less stringent or different. Our analysis is the first to compare the TCR repertoire of naïve CD8⁺ T cells between viable preterm neonates and term neonates. We find preterm neonates have a repertoire immaturity which potentially contributes to their increased infection susceptibility. A developmentally regulated, evenly distributed repertoire in preterm neonates may lead to the inclusion of public TCR CDR3β sequences that overlap between unrelated individuals in the preterm repertoire.

Comparative analysis of murine T-cell receptor repertoires

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

Establishment and stability of the latent HIV-1 DNA reservoir

HIV-1 infection cannot be cured because the virus persists as integrated proviral DNA in long-lived cells despite years of suppressive antiretroviral therapy (ART). In a previous paper (Zanini et al, 2015) we documented HIV-1 evolution in 10 untreated patients. Here we characterize establishment, turnover, and evolution of viral DNA reservoirs in the same patients after 3–18 years of suppressive ART. A median of 14% (range 0–42%) of the DNA sequences were defective due to G-to-A hypermutation. Remaining DNA sequences showed no evidence of evolution over years of suppressive ART. Most sequences from the DNA reservoirs were very similar to viruses actively replicating in plasma (RNA sequences) shortly before start of ART. The results do not support persistent HIV-1 replication as a mechanism to maintain the HIV-1 reservoir during suppressive therapy. Rather, the data indicate that DNA variants are turning over as long as patients are untreated and that suppressive ART halts this turnover.

DOI:http://dx.doi.org/10.7554/eLife.18889.001