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

Human SARS-CoV-2 challenge uncovers local and systemic response dynamics

The COVID-19 pandemic is an ongoing global health threat, yet our understanding of the dynamics of early cellular responses to this disease remains limited1. Here in our SARS-CoV-2 human challenge study, we used single-cell multi-omics profiling of nasopharyngeal swabs and blood to temporally resolve abortive, transient and sustained infections in seronegative individuals challenged with pre-Alpha SARS-CoV-2. Our analyses revealed rapid changes in cell-type proportions and dozens of highly dynamic cellular response states in epithelial and immune cells associated with specific time points and infection status. We observed that the interferon response in blood preceded the nasopharyngeal response. Moreover, nasopharyngeal immune infiltration occurred early in samples from individuals with only transient infection and later in samples from individuals with sustained infection. High expression of HLA-DQA2 before inoculation was associated with preventing sustained infection. Ciliated cells showed multiple immune responses and were most permissive for viral replication, whereas nasopharyngeal T cells and macrophages were infected non-productively. We resolved 54 T cell states, including acutely activated T cells that clonally expanded while carrying convergent SARS-CoV-2 motifs. Our new computational pipeline Cell2TCR identifies activated antigen-responding T cells based on a gene expression signature and clusters these into clonotype groups and motifs. Overall, our detailed time series data can serve as a Rosetta stone for epithelial and immune cell responses and reveals early dynamic responses associated with protection against infection.

FUME-TCRseq Enables Sensitive and Accurate Sequencing of the T-cell Receptor from Limited Input of Degraded RNA

Genomic analysis of the T-cell receptor (TCR) reveals the strength, breadth, and clonal dynamics of the adaptive immune response to pathogens or cancer. The diversity of the TCR repertoire, however, means that sequencing is technically challenging, particularly for samples with low-quality, degraded nucleic acids. Here, we developed and validated FUME-TCRseq, a robust and sensitive RNA-based TCR sequencing methodology that is suitable for formalin-fixed paraffin-embedded samples and low amounts of input material. FUME-TCRseq incorporates unique molecular identifiers into each molecule of cDNA, allowing correction for sequencing errors and PCR bias. Using RNA extracted from colorectal and head and neck cancers to benchmark the accuracy and sensitivity of FUME-TCRseq against existing methods demonstrated excellent concordance between the datasets. Furthermore, FUME-TCRseq detected more clonotypes than a commercial RNA-based alternative, with shorter library preparation time and significantly lower cost. The high sensitivity and the ability to sequence RNA of poor quality and limited amount enabled quantitative analysis of small numbers of cells from archival tissue sections, which is not possible with other methods. Spatially resolved FUME-TCRseq analysis of colorectal cancers using macrodissected archival samples revealed the shifting T-cell landscapes at the transition to an invasive phenotype and between tumor subclones containing distinct driver alterations. In summary, FUME-TCRseq represents an accurate, sensitive, and low-cost tool for the characterization of T-cell repertoires, particularly in samples with low-quality RNA that have not been accessible using existing methodology.

SIGNIFICANCE

FUME-TCRseq is a TCR sequencing methodology that supports sensitive and spatially resolved detection of TCR clones in archival clinical specimens, which can facilitate longitudinal tracking of immune responses through disease course and treatment.

T-cell receptor diversity and allergen sensitivity in childhood asthma and atopic dermatitis

BACKGROUND

Childhood allergies of asthma and atopic dermatitis (AD) involve an overactive T-cell immune response triggered by allergens. However, the impact of T-cell receptor (TCR) repertoires on allergen sensitization and their role in mediating different phenotypes of asthma and AD in early childhood remains unclear.

METHODS

A total of 78 children, comprising 26 with asthma alone, 26 with AD alone, and 26 healthy controls (HC), were enrolled. TCR repertoire profiles were determined using a unique molecular identifier system for next-generation sequencing. Integrative analyses of their associations with allergen-specific IgE levels and allergies were performed.

RESULTS

The diversity in TCR alpha variable region (TRAV) genes of TCR repertoires and complementarity determining region 3 (CDR3) clonality in TRAV/TRBV (beta) genes were significantly higher in children with AD compared with those with asthma and HC (p < .05). Compared with HC, the expression of TRAV13-1 and TRAV4 genes was significantly higher in both asthma and AD (p < .05), with a significant positive correlation with mite-specific IgE levels (p < .01). In contrast, TRBV7-9 gene expression was significantly lower in both asthma and AD (p < .01), with this gene showing a significant negative correlation with mite-specific IgE levels (p < .01). Furthermore, significantly higher TRAV8-3 gene expression, positively correlated with food-specific IgE levels, was found in children with AD compared with those with asthma (p < .05).

CONCLUSION

Integrated TCR repertoires analysis provides clinical insights into the diverse TCR genes linked to antigen specificity, offering potential for precision immunotherapy in childhood allergies.

Enhancing comparative T cell receptor repertoire analysis in small biological samples through pooling homologous cell samples from multiple mice

Accurate characterization and comparison of T cell receptor (TCR) repertoires from small biological samples present significant challenges. The main challenge is the low material input, which compromises the quality of bulk sequencing and hinders the recovery of sufficient TCR sequences for robust analyses. We aimed to address this limitation by implementing a strategic approach to pool homologous biological samples. Our findings demonstrate that such pooling indeed enhances the TCR repertoire coverage, particularly for cell subsets of constrained sizes, and enables accurate comparisons of TCR repertoires at different levels of complexity across T cell subsets with different sizes. This methodology holds promise for advancing our understanding of T cell repertoires in scenarios where sample size constraints are a prevailing concern.

Personalized neoantigen-based T cell therapy triggers cytotoxic lymphocytes expressing polyclonal TCR against metastatic ovarian cancer

Neoantigen-reactive cytotoxic T lymphocytes play a vital role in precise cancer cell elimination. In this study, we demonstrate the effectiveness of personalized neoantigen-based T cell therapy in inducing tumor regression in two patients suffering from heavily-burdened metastatic ovarian cancer. Our approach involved the development of a robust pipeline for ex vivo expansion of neoantigen-reactive T lymphocytes. Neoantigen peptides were designed and synthesized based on the somatic mutations of the tumors and their predicted HLA binding affinities. These peptides were then presented to T lymphocytes through co-culture with neoantigen-loaded dendritic cells for ex vivo expansion. Subsequent to cell therapy, both patients exhibited significant reductions in tumor marker levels and experienced substantial tumor regression. One patient achieved repeated cancer regression through infusions of T cell products generated from newly identified neoantigens. Transcriptomic analyses revealed a remarkable increase in neoantigen-reactive cytotoxic lymphocytes in the peripheral blood of the patients following cell therapy. These cytotoxic T lymphocytes expressed polyclonal T cell receptors (TCR) against neoantigens, along with abundant cytotoxic proteins and pro-inflammatory cytokines. The efficacy of neoantigen targeting was significantly associated with the immunogenicity and TCR polyclonality. Notably, the neoantigen-specific TCR clonotypes persisted in the peripheral blood after cell therapy. Our findings indicate that personalized neoantigen-based T cell therapy triggers cytotoxic lymphocytes expressing polyclonal TCR against ovarian cancer, suggesting its promising potential in cancer immunotherapy.

High Epstein-Barr virus capsid antigen IgG level associates with the carriership of CD8+ T cell somatic mutations in the STAT3 SH2 domain

High carrier prevalence of STAT3 SH2 domain somatic mutations was recently discovered in CD8+ T cells. We found these low-allele-fraction clones in 26% of donors, without difference between multiple sclerosis (MS) patients and controls. Here we tested whether anti-viral antibodies associate with the carriership of these mutant clones. We compared antibody responses against common viruses in mutation carriers vs. non-carriers. Plasma samples of 152 donors (92 MS patients, 60 controls) were analyzed for antibodies against cytomegalovirus (CMV), Epstein-Barr virus (EBV), human herpesvirus-6A and parvovirus B19. The mutation carrier status associated with EBV VCA IgG level (p = 0.005) and remained significant after logistic regression (p = 0.036). This association was contributed similarly by MS patients and controls. These results suggest that EBV contributes to the generation or growth of these clones. The pathogenic role of the STAT3 mutant clones in MS is presently unclear, but their detailed characterization warrants further study.

Early activation of inflammatory pathways in UBA1-mutated hematopoietic stem and progenitor cells in VEXAS

VEXAS (vacuoles, E1 enzyme, X-linked, autoinflammatory, somatic) syndrome is a pleiotropic, severe autoinflammatory disease caused by somatic mutations in the ubiquitin-like modifier activating enzyme 1 (UBA1) gene. To elucidate VEXAS pathophysiology, we performed transcriptome sequencing of single bone marrow mononuclear cells and hematopoietic stem and progenitor cells (HSPCs) from VEXAS patients. HSPCs are biased toward myeloid (granulocytic) differentiation, and against lymphoid differentiation in VEXAS. Activation of multiple inflammatory pathways (interferons and tumor necrosis factor alpha) occurs ontogenically early in primitive hematopoietic cells and particularly in the myeloid lineage in VEXAS, and inflammation is prominent in UBA1-mutated cells. Dysregulation in protein degradation likely leads to higher stress response in VEXAS HSPCs, which positively correlates with inflammation. TCR usage is restricted and there are increased cytotoxicity and IFN-γ signaling in T cells. In VEXAS syndrome, both aberrant inflammation and myeloid predominance appear intrinsic to hematopoietic stem cells mutated in UBA1.

BERTrand-peptide:TCR binding prediction using Bidirectional Encoder Representations from Transformers augmented with random TCR pairing

MOTIVATION

The advent of T-cell receptor (TCR) sequencing experiments allowed for a significant increase in the amount of peptide:TCR binding data available and a number of machine-learning models appeared in recent years. High-quality prediction models for a fixed epitope sequence are feasible, provided enough known binding TCR sequences are available. However, their performance drops significantly for previously unseen peptides.

RESULTS

We prepare the dataset of known peptide:TCR binders and augment it with negative decoys created using healthy donors' T-cell repertoires. We employ deep learning methods commonly applied in Natural Language Processing to train part a peptide:TCR binding model with a degree of cross-peptide generalization (0.69 AUROC). We demonstrate that BERTrand outperforms the published methods when evaluated on peptide sequences not used during model training.

AVAILABILITY AND IMPLEMENTATION

The datasets and the code for model training are available at https://github.com/SFGLab/bertrand.

Large clones of pre-existing T cells drive early immunity against SARS-COV-2 and LCMV infection

T cell responses precede antibody and may provide early control of infection. We analyzed the clonal basis of this rapid response following SARS-COV-2 infection. We applied T cell receptor (TCR) sequencing to define the trajectories of individual T cell clones immediately. In SARS-COV-2 PCR+ individuals, a wave of TCRs strongly but transiently expand, frequently peaking the same week as the first positive PCR test. These expanding TCR CDR3s were enriched for sequences functionally annotated as SARS-COV-2 specific. Epitopes recognized by the expanding TCRs were highly conserved between SARS-COV-2 strains but not with circulating human coronaviruses. Many expanding CDR3s were present at high frequency in pre-pandemic repertoires. Early response TCRs specific for lymphocytic choriomeningitis virus epitopes were also found at high frequency in the preinfection naive repertoire. High-frequency naive precursors may allow the T cell response to respond rapidly during the crucial early phases of acute viral infection.

Viral infection reveals hidden sharing of TCR CDR3 sequences between individuals

The T cell receptor is generated by a process of random and imprecise somatic recombination. The number of possible T cell receptors which this process can produce is enormous, greatly exceeding the number of T cells in an individual. Thus, the likelihood of identical TCRs being observed in multiple individuals (public TCRs) might be expected to be very low. Nevertheless such public TCRs have often been reported. In this study we explore the extent of TCR publicity in the context of acute resolving Lymphocytic choriomeningitis virus (LCMV) infection in mice. We show that the repertoire of effector T cells following LCMV infection contains a population of highly shared TCR sequences. This subset of TCRs has a distribution of naive precursor frequencies, generation probabilities, and physico-chemical CDR3 properties which lie between those of classic public TCRs, which are observed in uninfected repertoires, and the dominant private TCR repertoire. We have named this set of sequences "hidden public" TCRs, since they are only revealed following infection. A similar repertoire of hidden public TCRs can be observed in humans after a first exposure to SARS-COV-2. The presence of hidden public TCRs which rapidly expand following viral infection may therefore be a general feature of adaptive immunity, identifying an additional level of inter-individual sharing in the TCR repertoire which may form an important component of the effector and memory response.

Quantitative annotations of T-Cell repertoire specificity

The specificity of a T-cell receptor (TCR) repertoire determines personalized immune capacity. Existing methods have modeled the qualitative aspects of TCR specificity, while the quantitative aspects remained unaddressed. We developed a package, TCRanno, to quantify the specificity of TCR repertoires. We created deep-learning-based, epitope-aware vector embeddings to infer individual TCR specificity. Then we aggregated clonotype frequencies of TCRs to obtain a quantitative profile of repertoire specificity at epitope, antigen and organism levels. Applying TCRanno to 4195 TCR repertoires revealed quantitative changes in repertoire specificity upon infections, autoimmunity and cancers. Specifically, TCRanno found cytomegalovirus-specific TCRs in seronegative healthy individuals, supporting the possibility of abortive infections. TCRanno discovered age-accumulated fraction of severe acute respiratory syndrome coronavirus 2 specific TCRs in pre-pandemic samples, which may explain the aggressive symptoms and age-related severity of coronavirus disease 2019. TCRanno also identified the encounter of Hepatitis B antigens as a potential trigger of systemic lupus erythematosus. TCRanno annotations showed capability in distinguishing TCR repertoires of healthy and cancers including melanoma, lung and breast cancers. TCRanno also demonstrated usefulness to single-cell TCRseq+gene expression data analyses by isolating T-cells with the specificity of interest.

Correction of T-Cell Repertoire and Autoimmune Diabetes in NOD Mice by Non-myeloablative T-Cell Depleted Allogeneic HSCT

The induction of partial tolerance toward pancreatic autoantigens in the treatment of type 1 diabetes mellitus (T1DM) can be attained by autologous hematopoietic stem cell transplantation (HSCT). However, most patients treated by autologous HSCT eventually relapse. Furthermore, allogeneic HSCT which could potentially provide a durable non-autoimmune T-cell receptor (TCR) repertoire is associated with a substantial risk for transplant-related mortality. We have previously demonstrated an effective approach for attaining engraftment without graft versus host disease (GVHD) of allogeneic T-cell depleted HSCT, following non-myeloablative conditioning, using donor-derived anti-3rd party central memory CD8 veto T cells (Tcm). In the present study, we investigated the ability of this relatively safe transplant modality to eliminate autoimmune T-cell clones in the NOD mouse model which spontaneously develop T1DM. Our results demonstrate that using this approach, marked durable chimerism is attained, without any transplant-related mortality, and with a very high rate of diabetes prevention. TCR sequencing of transplanted mice showed profound changes in the T-cell repertoire and decrease in the prevalence of specific autoimmune T-cell clones directed against pancreatic antigens. This approach could be considered as strategy to treat people destined to develop T1DM but with residual beta cell function, or as a platform for prevention of beta cell destruction after transplantation of allogenic beta cells.

Naive T cells inhibit the outgrowth of intractable antigen-activated memory T cells: implications for T-cell immunotherapy

BACKGROUND

The wider application of T cells targeting viral tumor-antigens via their native receptors is hampered by the failure to expand potent tumor-specific T cells from patients. Here, we examine reasons for and solutions to this failure, taking as our model the preparation of Epstein-Barr virus (EBV)-specific T cells (EBVSTs) for the treatment of EBV-positive lymphoma. EBVSTs could not be manufactured from almost one-third of patients, either because they failed to expand, or they expanded, but lacked EBV specificity. We identified an underlying cause of this problem and established a clinically feasible approach to overcome it.

METHODS

CD45RO+CD45RA- memory compartment residing antigen-specific T cells were enriched by depleting CD45RA positive (+) peripheral blood mononuclear cells (PBMCs) that include naïve T cells, among other subsets, prior to EBV antigen stimulation. We then compared the phenotype, specificity, function and T-cell receptor (TCR) Vβ repertoire of EBVSTs expanded from unfractionated whole (W)-PBMCs and CD45RA-depleted (RAD)-PBMCs on day 16. To identify the CD45RA component that inhibited EBVST outgrowth, isolated CD45RA+ subsets were added back to RAD-PBMCs followed by expansion and characterization. The in vivo potency of W-EBVSTs and RAD-EBVSTs was compared in a murine xenograft model of autologous EBV+ lymphoma.

RESULTS

Depletion of CD45RA+ PBMCs before antigen stimulation increased EBVST expansion, antigen-specificity and potency in vitro and in vivo. TCR sequencing revealed a selective outgrowth in RAD-EBVSTs of clonotypes that expanded poorly in W-EBVSTs. Inhibition of antigen-stimulated T cells by CD45RA+ PBMCs could be reproduced only by the naïve T-cell fraction, while CD45RA+ regulatory T cells, natural killer cells, stem cell memory and effector memory subsets lacked inhibitory activity. Crucially, CD45RA depletion of PBMCs from patients with lymphoma enabled the outgrowth of EBVSTs that failed to expand from W-PBMCs. This enhanced specificity extended to T cells specific for other viruses.

CONCLUSION

Our findings suggest that naïve T cells inhibit the outgrowth of antigen-stimulated memory T cells, highlighting the profound effects of intra-T-cell subset interactions. Having overcome our inability to generate EBVSTs from many patients with lymphoma, we have introduced CD45RA depletion into three clinical trials: NCT01555892 and NCT04288726 using autologous and allogeneic EBVSTs to treat lymphoma and NCT04013802 using multivirus-specific T cells to treat viral infections after hematopoietic stem cell transplantation.

The use of non-functional clonotypes as a natural calibrator for quantitative bias correction in adaptive immune receptor repertoire profiling

High-throughput sequencing of adaptive immune receptor repertoires is a valuable tool for receiving insights in adaptive immunity studies. Several powerful TCR/BCR repertoire reconstruction and analysis methods have been developed in the past decade. However, detecting and correcting the discrepancy between real and experimentally observed lymphocyte clone frequencies are still challenging. Here, we discovered a hallmark anomaly in the ratio between read count and clone count-based frequencies of non-functional clonotypes in multiplex PCR-based immune repertoires. Calculating this anomaly, we formulated a quantitative measure of V- and J-genes frequency bias driven by multiplex PCR during library preparation called Over Amplification Rate (OAR). Based on the OAR concept, we developed an original software for multiplex PCR-specific bias evaluation and correction named iROAR: immune Repertoire Over Amplification Removal (https://github.com/smiranast/iROAR). The iROAR algorithm was successfully tested on previously published TCR repertoires obtained using both 5' RACE (Rapid Amplification of cDNA Ends)-based and multiplex PCR-based approaches and compared with a biological spike-in-based method for PCR bias evaluation. The developed approach can increase the accuracy and consistency of repertoires reconstructed by different methods making them more applicable for comparative analysis.

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.

Into the multi-omics era: Progress of T cells profiling in the context of solid organ transplantation

T cells are the common type of lymphocyte to mediate allograft rejection, remaining long-term allograft survival impeditive. However, the heterogeneity of T cells, in terms of differentiation and activation status, the effector function, and highly diverse T cell receptors (TCRs) have thus precluded us from tracking these T cells and thereby comprehending their fate in recipients due to the limitations of traditional detection approaches. Recently, with the widespread development of single-cell techniques, the identification and characterization of T cells have been performed at single-cell resolution, which has contributed to a deeper comprehension of T cell heterogeneity by relevant detections in a single cell - such as gene expression, DNA methylation, chromatin accessibility, surface proteins, and TCR. Although these approaches can provide valuable insights into an individual cell independently, a comprehensive understanding can be obtained when applied joint analysis. Multi-omics techniques have been implemented in characterizing T cells in health and disease, including transplantation. This review focuses on the thesis, challenges, and advances in these technologies and highlights their application to the study of alloreactive T cells to improve the understanding of T cell heterogeneity in solid organ transplantation.

Substantial uneven proliferation of CD4+ T cells during recovery from acute HIV infection is sufficient to explain the observed expanded clones in the HIV reservoir

The HIV reservoir is a population of 1-10 million anatomically dispersed, latently infected memory CD4+ T cells in which HIV DNA is quiescently integrated into human chromosomal DNA. When antiretroviral therapy (ART) is stopped and HIV replication initiates in one of these cells, systemic viral spread resumes, rekindling progression to AIDS. Therefore, HIV latency prevents cure. The detection of many populations of identical HIV sequences at unique integration sites implicates CD4+ T cell proliferation as the critical driver of reservoir sustainment after a prolonged period of effective ART. Initial reservoir formation occurs during the first week of primary infection usually before ART is started. While empirical data indicates that both de novo infection and cellular proliferation generate latently infected cells during early untreated infection, it is not known which of these mechanisms is predominant. We developed a mathematical model that recapitulates the profound depletion and brisk recovery of CD4+ T cells, reservoir creation, and viral load trajectory during primary HIV infection. We extended the model to stochastically simulate individual HIV reservoir clones. This model predicts the first detection of HIV infected clones approximately 5 weeks after infection as has recently been shown in vivo and suggests that substantial, uneven proliferation among clones during the recovery from CD4+ lymphopenia is the most plausible explanation for the observed clonal reservoir distribution during the first year of infection.

Phenomenology and dynamics of competitive ecosystems beyond the niche-neutral regimes

Structure, composition, and stability of ecological populations are shaped by the inter- and intraspecies interactions within their communities. It remains to be fully understood how the interplay of these interactions with other factors, such as immigration, controls the structure, the diversity, and the long-term stability of ecological systems in the presence of noise and fluctuations. We address this problem using a minimal model of interacting multispecies ecological communities that incorporates competition, immigration, and demographic noise. We find that a complete phase diagram exhibits rich behavior with multiple regimes that go beyond the classical "niche" and "neutral" regimes, extending and modifying the "rare biosphere" or "niche-like" dichotomy. In particular, we observe regimes that cannot be characterized as either niche or neutral where a multimodal species abundance distribution is observed. We characterize the transitions between the different regimes and show how these arise from the underlying kinetics of the species turnover, extinction, and invasion. Our model serves as a minimal null model of noisy competitive ecological systems, against which more complex models that include factors such as mutations and environmental noise can be compared.

Spatial maps of T cell receptors and transcriptomes reveal distinct immune niches and interactions in the adaptive immune response

T cells mediate antigen-specific immune responses to disease through the specificity and diversity of their clonotypic T cell receptors (TCRs). Determining the spatial distributions of T cell clonotypes in tissues is essential to understanding T cell behavior, but spatial sequencing methods remain unable to profile the TCR repertoire. Here, we developed Slide-TCR-seq, a 10-μm-resolution method, to sequence whole transcriptomes and TCRs within intact tissues. We confirmed the ability of Slide-TCR-seq to map the characteristic locations of T cells and their receptors in mouse spleen. In human lymphoid germinal centers, we identified spatially distinct TCR repertoires. Profiling T cells in renal cell carcinoma and melanoma specimens revealed heterogeneous immune responses: T cell states and infiltration differed intra- and inter-clonally, and adjacent tumor and immune cells exhibited distinct gene expression. Altogether, our method yields insights into the spatial relationships between clonality, neighboring cell types, and gene expression that drive T cell responses.

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.

Longitudinal analysis reveals age-related changes in the T cell receptor repertoire of human T cell subsets

A diverse T cell receptor (TCR) repertoire is essential for protection against a variety of pathogens, and TCR repertoire size is believed to decline with age. However, the precise size of human TCR repertoires, in both total and subsets of T cells, as well as their changes with age, are not fully characterized. We conducted a longitudinal analysis of the human blood TCRα and TCRβ repertoire of CD4⁺ and CD8⁺ T cell subsets using a unique molecular identifier–based (UMI-based) RNA-seq method. Thorough analysis of 1.9 × 10⁸ T cells yielded the lower estimate of TCR repertoire richness in an adult at 3.8 × 10⁸. Alterations of the TCR repertoire with age were observed in all 4 subsets of T cells. The greatest reduction was observed in naive CD8⁺ T cells, while the greatest clonal expansion was in memory CD8⁺ T cells, and the highest increased retention of TCR sequences was in memory CD8⁺ T cells. Our results demonstrated that age-related TCR repertoire attrition is subset specific and more profound for CD8⁺ than CD4⁺ T cells, suggesting that aging has a more profound effect on cytotoxic as opposed to helper T cell functions. This may explain the increased susceptibility of older adults to novel infections.

How Naive T-Cell Clone Counts Are Shaped By Heterogeneous Thymic Output and Homeostatic Proliferation

The specificity of T cells is that each T cell has only one T cell receptor (TCR). A T cell clone represents a collection of T cells with the same TCR sequence. Thus, the number of different T cell clones in an organism reflects the number of different T cell receptors (TCRs) that arise from recombination of the V(D)J gene segments during T cell development in the thymus. TCR diversity and more specifically, the clone abundance distribution, are important factors in immune functions. Specific recombination patterns occur more frequently than others while subsequent interactions between TCRs and self-antigens are known to trigger proliferation and sustain naive T cell survival. These processes are TCR-dependent, leading to clone-dependent thymic export and naive T cell proliferation rates. We describe the heterogeneous steady-state population of naive T cells (those that have not yet been antigenically triggered) by using a mean-field model of a regulated birth-death-immigration process. After accounting for random sampling, we investigate how TCR-dependent heterogeneities in immigration and proliferation rates affect the shape of clone abundance distributions (the number of different clones that are represented by a specific number of cells, or “clone counts”). By using reasonable physiological parameter values and fitting predicted clone counts to experimentally sampled clone abundances, we show that realistic levels of heterogeneity in immigration rates cause very little change to predicted clone-counts, but that modest heterogeneity in proliferation rates can generate the observed clone abundances. Our analysis provides constraints among physiological parameters that are necessary to yield predictions that qualitatively match the data. Assumptions of the model and potentially other important mechanistic factors are discussed.

Rapid synchronous type 1 IFN and virus-specific T cell responses characterize first wave non-severe SARS-CoV-2 infections

Effective control of SARS-CoV-2 infection on primary exposure may reveal correlates of protective immunity to future variants, but we lack insights into immune responses before or at the time virus is first detected. We use blood transcriptomics, multiparameter flow cytometry, and T cell receptor (TCR) sequencing spanning the time of incident non-severe infection in unvaccinated virus-naive individuals to identify rapid type 1 interferon (IFN) responses common to other acute respiratory viruses and cell proliferation responses that discriminate SARS-CoV-2 from other viruses. These peak by the time the virus is first detected and sometimes precede virus detection. Cell proliferation is most evident in CD8 T cells and associated with specific expansion of SARS-CoV-2-reactive TCRs, in contrast to virus-specific antibodies, which lag by 1-2 weeks. Our data support a protective role for early type 1 IFN and CD8 T cell responses, with implications for development of universal T cell vaccines.

TCRpower: quantifying the detection power of T-cell receptor sequencing with a novel computational pipeline calibrated by spike-in sequences

T-cell receptor (TCR) sequencing has enabled the development of innovative diagnostic tests for cancers, autoimmune diseases and other applications. However, the rarity of many T-cell clonotypes presents a detection challenge, which may lead to misdiagnosis if diagnostically relevant TCRs remain undetected. To address this issue, we developed TCRpower, a novel computational pipeline for quantifying the statistical detection power of TCR sequencing methods. TCRpower calculates the probability of detecting a TCR sequence as a function of several key parameters: in-vivo TCR frequency, T-cell sample count, read sequencing depth and read cutoff. To calibrate TCRpower, we selected unique TCRs of 45 T-cell clones (TCCs) as spike-in TCRs. We sequenced the spike-in TCRs from TCCs, together with TCRs from peripheral blood, using a 5' RACE protocol. The 45 spike-in TCRs covered a wide range of sample frequencies, ranging from 5 per 100 to 1 per 1 million. The resulting spike-in TCR read counts and ground truth frequencies allowed us to calibrate TCRpower. In our TCR sequencing data, we observed a consistent linear relationship between sample and sequencing read frequencies. We were also able to reliably detect spike-in TCRs with frequencies as low as one per million. By implementing an optimized read cutoff, we eliminated most of the falsely detected sequences in our data (TCR α-chain 99.0% and TCR β-chain 92.4%), thereby improving diagnostic specificity. TCRpower is publicly available and can be used to optimize future TCR sequencing experiments, and thereby enable reliable detection of disease-relevant TCRs for diagnostic applications.

The importance of taking ART appropriately in children and adolescents with HIV-1 to reach the highest capacity of immune function later in life

Current antiretroviral therapy (ART) guidelines recommend treating all children with HIV-1 infection. This has changed from the broader use of ART to treat children to improve morbidity and minimise mortality. However, prior to current recommendations, not everyone with HIV-1 received timely treatment. What happens to the paediatric immune system when HIV-1 replication is not appropriately supressed remains unclear. 11 samples from adolescents with HIV-1 on ART and uninfected controls in the UK, aged 12-25 years, were examined; overall, adolescents with CD4⁺ counts > 500/μl and a viral load < 50 copies/ml were compared with adolescents with CD4⁺ counts < 500/μl and a viral load > 50 copies/ml at time of sampling. Measurements of thymic output were combined with high throughput next generation sequencing and bioinformatics to systematically organize CD4⁺ and CD8⁺ T cell receptor (TCR) repertoires. TCR repertoire diversity, clonal expansions, TCR sequence sharing, and formation of TCR clusters in HIV-1 infected adolescents with successful HIV-1 suppression were compared to adolescents with ineffective HIV-1 suppression. Thymic output and CD4⁺ T cell numbers were decreased in HIV-1 infected adolescents with poor HIV-1 suppression. A strong homeostatic TCR response, driven by the decreased CD4⁺ T cell compartment and reduced thymic output was observed in the virally uncontrolled HIV-1-infected adolescents. Formation of abundant robust TCR clusters and structurally related TCRs were found in the adolescents with effective HIV-1 suppression. Numerous CD4⁺ T cell numbers in the virally controlled adolescents emphasize the importance of high thymic output and formation of robust TCR clusters in the maintenance of HIV-1 suppression. While the profound capacity for immune recovery in children may allow better opportunity to deal with immunological stress, when ART is taken appropriately, this study demonstrates new insights into the unique paediatric immune system and the immunological changes when HIV-1 replication is ongoing.

OUP accepted manuscript

The study and manipulation of T cell receptors (TCRs) is central to multiple fields across basic and translational immunology research. Produced by V(D)J recombination, TCRs are often only recorded in the literature and data repositories as a combination of their V and J gene symbols, plus their hypervariable CDR3 amino acid sequence. However, numerous applications require full-length coding nucleotide sequences. Here we present Stitchr, a software tool developed to specifically address this limitation. Given minimal V/J/CDR3 information, Stitchr produces complete coding sequences representing a fully spliced TCR cDNA. Due to its modular design, Stitchr can be used for TCR engineering using either published germline or novel/modified variable and constant region sequences. Sequences produced by Stitchr were validated by synthesizing and transducing TCR sequences into Jurkat cells, recapitulating the expected antigen specificity of the parental TCR. Using a companion script, Thimble, we demonstrate that Stitchr can process a million TCRs in under ten minutes using a standard desktop personal computer. By systematizing the production and modification of TCR sequences, we propose that Stitchr will increase the speed, repeatability, and reproducibility of TCR research. Stitchr is available on GitHub.

The Intra-Tumoral T Cell Receptor Repertoire: Steps Towards a Useful Clinical Biomarker

Adaptive immunity recognizes and responds to tumors, although they are part of the immunological "self." T cells, both CD4+ and CD8+, play a key role in the process, and the specific set of receptors which recognize tumor antigens therefore has the potential to provide prognostic biomarkers for tracking tumor growth after cancer therapy, including immunotherapy. Most published data on the T cell repertoire continue to rely on commercial proprietary methods, which often do not allow access to the raw data, and are difficult to validate. We describe an open-source protocol for amplifying, sequencing, and analyzing T cell receptors which is economical, robust, sensitive, and versatile. The key experimental step is the ligation of a single-stranded oligonucleotide to the 3' end of the T cell receptor cDNA, which allows easy amplification of all possible rearrangements using only a single set of primers per locus, while simultaneously introducing a unique molecular identifier to label each starting cDNA molecule. After sequencing, this molecular identifier can be used to correct both sequence errors and the effects of differential PCR amplification efficiency, thus producing a more accurate measure of the true T cell receptor frequency within the sample. Samples are then tagged with unique pairs of indices, facilitating robotic scale-up and significantly reducing cross-sample contamination from index hopping. This method has been applied to the analysis of tumor-infiltrating lymphocytes and matched peripheral blood samples from patients with a variety of solid tumors.

CD8+ cell somatic mutations in multiple sclerosis patients and controls-Enrichment of mutations in STAT3 and other genes implicated in hematological malignancies

Somatic mutations have a central role in cancer but their role in other diseases such as common autoimmune disorders is not clear. Previously we and others have demonstrated that especially CD8+ T cells in blood can harbor persistent somatic mutations in some patients with multiple sclerosis (MS) and rheumatoid arthritis. Here we concentrated on CD8+ cells in more detail and tested (i) how commonly somatic mutations are detectable, (ii) does the overall mutation load differ between MS patients and controls, and (iii) do the somatic mutations accumulate non-randomly in certain genes? We separated peripheral blood CD8+ cells from newly diagnosed relapsing MS patients (n = 21) as well as matched controls (n = 21) and performed next-generation sequencing of the CD8+ cells' DNA, limiting our search to a custom panel of 2524 immunity and cancer related genes, which enabled us to obtain a median sequencing depth of over 2000x. We discovered nonsynonymous somatic mutations in all MS patients' and controls' CD8+ cell DNA samples, with no significant difference in number between the groups (p = 0.60), at a median allelic fraction of 0.5% (range 0.2-8.6%). The mutations showed statistically significant clustering especially to the STAT3 gene, and also enrichment to the SMARCA2, DNMT3A, SOCS1 and PPP3CA genes. Known activating STAT3 mutations were found both in MS patients and controls and overall 1/5 of the mutations were previously described cancer mutations. The detected clustering suggests a selection advantage of the mutated CD8+ clones and calls for further research on possible phenotypic effects.

Metabolite and thymocyte development defects in ADA-SCID mice receiving enzyme replacement therapy

Deficiency of adenosine deaminase (ADA, EC3.5.4.4), a housekeeping enzyme intrinsic to the purine salvage pathway, leads to severe combined immunodeficiency (SCID) both in humans and mice. Lack of ADA results in the intracellular accumulation of toxic metabolites which have effects on T cell development and function. While untreated ADA-SCID is a fatal disorder, there are different therapeutic options available to restore ADA activity and reconstitute a functioning immune system, including enzyme replacement therapy (ERT). Administration of ERT in the form of pegylated bovine ADA (PEG-ADA) has proved a life-saving though non-curative treatment for ADA-SCID patients. However, in many patients treated with PEG-ADA, there is suboptimal immune recovery with low T and B cell numbers. Here, we show reduced thymus cellularity in ADA-SCID mice despite weekly PEG-ADA treatment. This was associated with lack of effective adenosine (Ado) detoxification in the thymus. We also show that thymocyte development in ADA-deficient thymi is arrested at the DN3-to-DN4 stage transition with thymocytes undergoing dATP-induced apoptosis rather than defective TCRβ rearrangement or β-selection. Our studies demonstrate at a detailed level that exogenous once-a-week enzyme replacement does not fully correct intra-thymic metabolic or immunological abnormalities associated with ADA deficiency.

Determinants of anti-PD-1 response and resistance in clear cell renal cell carcinoma

ADAPTeR is a prospective, phase II study of nivolumab (anti-PD-1) in 15 treatment-naive patients (115 multiregion tumor samples) with metastatic clear cell renal cell carcinoma (ccRCC) aiming to understand the mechanism underpinning therapeutic response. Genomic analyses show no correlation between tumor molecular features and response, whereas ccRCC-specific human endogenous retrovirus expression indirectly correlates with clinical response. T cell receptor (TCR) analysis reveals a significantly higher number of expanded TCR clones pre-treatment in responders suggesting pre-existing immunity. Maintenance of highly similar clusters of TCRs post-treatment predict response, suggesting ongoing antigen engagement and survival of families of T cells likely recognizing the same antigens. In responders, nivolumab-bound CD8+ T cells are expanded and express GZMK/B. Our data suggest nivolumab drives both maintenance and replacement of previously expanded T cell clones, but only maintenance correlates with response. We hypothesize that maintenance and boosting of a pre-existing response is a key element of anti-PD-1 mode of action.

High Throughput Human T Cell Receptor Sequencing: A New Window Into Repertoire Establishment and Alloreactivity

Recent advances in high throughput sequencing (HTS) of T cell receptors (TCRs) and in transcriptomic analysis, particularly at the single cell level, have opened the door to a new level of understanding of human immunology and immune-related diseases. In this article, we discuss the use of HTS of TCRs to discern the factors controlling human T cell repertoire development and how this approach can be used in combination with human immune system (HIS) mouse models to understand human repertoire selection in an unprecedented manner. An exceptionally high proportion of human T cells has alloreactive potential, which can best be understood as a consequence of the processes governing thymic selection. High throughput TCR sequencing has allowed assessment of the development, magnitude and nature of the human alloresponse at a new level and has provided a tool for tracking the fate of pre-transplant-defined donor- and host-reactive TCRs following transplantation. New insights into human allograft rejection and tolerance obtained with this method in combination with single cell transcriptional analyses are reviewed here.

Postpartum breast cancer has a distinct molecular profile that predicts poor outcomes

Young women's breast cancer (YWBC) has poor prognosis and known interactions with parity. Women diagnosed within 5-10 years of childbirth, defined as postpartum breast cancer (PPBC), have poorer prognosis compared to age, stage, and biologic subtype-matched nulliparous patients. Genomic differences that explain this poor prognosis remain unknown. In this study, using RNA expression data from clinically matched estrogen receptor positive (ER+) cases (n = 16), we observe that ER+ YWBC can be differentiated based on a postpartum or nulliparous diagnosis. The gene expression signatures of PPBC are consistent with increased cell cycle, T-cell activation and reduced estrogen receptor and TP53 signaling. When applied to a large YWBC cohort, these signatures for ER+ PPBC associate with significantly reduced 15-year survival rates in high compared to low expressing cases. Cumulatively these results provide evidence that PPBC is a unique entity within YWBC with poor prognostic phenotypes.

Characteristics of T-cell receptor repertoire of stem cell-like memory CD4+ T cells

Stem cell-like memory T cells (Tscm) combine phenotypes of naïve and memory. However, it remains unclear how T cell receptor (TCR) characteristics contribute to heterogeneity in Tscm and other memory T cells. We compared the TCR-beta (TRB) repertoire characteristics of CD4+ Tscm with those of naïve and other CD4+ memory (Tm) in 16 human subjects. Compared with Tm, Tscm had an increased diversity across all stretches of TRB repertoire structure, a skewed gene usage, and a shorter length distribution of CDR3 region. These distinctions between Tscm and Tm were enlarged in top1000 abundant clonotypes. Furthermore, top1000 clonotypes in Tscm were more public than those in Tm and grouped in more clusters, implying more epitope types recognized by top1000 clonotypes in Tscm. Importantly, self-reactive clonotypes were public and enriched in Tscm rather than Tm, of type one diabetes patients. Therefore, this study highlights the unique features of Tscm different from those of other memory subsets and provides clues to understand the physiological and pathological functions of Tscm.

Immune Equilibrium Depends on the Interaction Between Recognition and Presentation Landscapes

In this review, we described the structure and organization of antigen-recognizing repertoires of B and T cells from the standpoint of modern immunology. We summarized the latest advances in bioinformatics analysis of sequencing data from T and B cell repertoires and also presented contemporary ideas about the mechanisms of clonal diversity formation at different stages of organism development. At the same time, we focused on the importance of the allelic variants of the HLA genes and spectra of presented antigens for the formation of T-cell receptors (TCR) landscapes. The main idea of this review is that immune equilibrium and proper functioning of immunity are highly dependent on the interaction between the recognition and the presentation landscapes of antigens. Certain changes in these landscapes can occur during life, which can affect the protective function of adaptive immunity. We described some mechanisms associated with these changes, for example, the conversion of effector cells into regulatory cells and vice versa due to the trans-differentiation or bystander effect, changes in the clonal organization of the general TCR repertoire due to homeostatic proliferation or aging, and the background for the altered presentation of some antigens due to SNP mutations of MHC, or the alteration of the presenting antigens due to post-translational modifications. The authors suggest that such alterations can lead to an increase in the risk of the development of oncological and autoimmune diseases and influence the sensitivity of the organism to different infectious agents.

Plasticity of the Immune System in Children Following Treatment Interruption in HIV-1 Infection

It is intriguing that, unlike adults with HIV-1, children with HIV-1 reach a greater CD4+ T cell recovery following planned treatment cessation. The reasons for the better outcomes in children remain unknown but may be related to increased thymic output and diversity of T cell receptor repertoires. HIV-1 infected children from the PENTA 11 trial tolerated planned treatment interruption without adverse long-term clinical, virological, or immunological consequences, once antiretroviral therapy was re-introduced. This contrasts to treatment interruption trials of HIV-1 infected adults, who had rapid changes in T cells and slow recovery when antiretroviral therapy was restarted. How children can develop such effective immune responses to planned treatment interruption may be critical for future studies. PENTA 11 was a randomized, phase II trial of planned treatment interruptions in HIV-1-infected children (ISRCTN 36694210). In this sub-study, eight patients in long-term follow-up were chosen with CD4+ count>500/ml, viral load <50c/ml at baseline: four patients on treatment interruption and four on continuous treatment. Together with measurements of thymic output, we used high-throughput next generation sequencing and bioinformatics to systematically organize memory CD8+ and naïve CD4+ T cell receptors according to diversity, clonal expansions, sequence sharing, antigen specificity, and T cell receptor similarities following treatment interruption compared to continuous treatment. We observed an increase in thymic output following treatment interruption compared to continuous treatment. This was accompanied by an increase in T cell receptor clonal expansions, increased T cell receptor sharing, and higher sequence similarities between patients, suggesting a more focused T cell receptor repertoire. The low numbers of patients included is a limitation and the data should be interpreted with caution. Nonetheless, the high levels of thymic output and the high diversity of the T cell receptor repertoire in children may be sufficient to reconstitute the T cell immune repertoire and reverse the impact of interruption of antiretroviral therapy. Importantly, the effective T cell receptor repertoires following treatment interruption may inform novel therapeutic strategies in children infected with HIV-1.

Autoencoder based local T cell repertoire density can be used to classify samples and T cell receptors

Recent advances in T cell repertoire (TCR) sequencing allow for the characterization of repertoire properties, as well as the frequency and sharing of specific TCR. However, there is no efficient measure for the local density of a given TCR. TCRs are often described either through their Complementary Determining region 3 (CDR3) sequences, or theirV/J usage, or their clone size. We here show that the local repertoire density can be estimated using a combined representation of these components through distance conserving autoencoders and Kernel Density Estimates (KDE). We present ELATE-an Encoder-based LocAl Tcr dEnsity and show that the resulting density of a sample can be used as a novel measure to study repertoire properties. The cross-density between two samples can be used as a similarity matrix to fully characterize samples from the same host. Finally, the same projection in combination with machine learning algorithms can be used to predict TCR-peptide binding through the local density of known TCRs binding a specific target.

Sex-Based Dimorphism of Anticancer Immune Response and Molecular Mechanisms of Immune Evasion

Purpose

We previously demonstrated that sex influences response to immune-checkpoint inhibitors.

Here we investigate sex-based differences in the molecular mechanisms of anticancer immune-response and immune evasion in patients with NSCLC.

Experimental Design

We analyzed a) transcriptome-data of 2575 early-stage NSCLCs from 7 different datasets; b) 327 tumor-samples extensively characterized at the molecular level from the TRACERx lung study; c) two independent cohorts of respectively 329 and 391 patients with advanced NSCLC treated with anti-PD1/anti-PDL1 drugs.

Results

As compared with men, the tumor microenvironment (TME) of women was significantly enriched for a number of innate and adaptive immune cell-types, including specific T-cell subpopulations.

NSCLCs of men and women exploited different mechanisms of immune evasion.

The TME of females was characterized by significantly greater T-cell dysfunction status, higher expression of inhibitory immune-checkpoint molecules and higher abundance of immune-suppressive cells, including Cancer Associated Fibroblasts, MDSCs and Regulatory T-cells.

By contrast, the TME of males was significantly enriched for a T-cells excluded phenotype.

We reported data supporting impaired neoantigens presentation to immune system in tumors of men, as molecular mechanism explaining the findings observed.

Finally, in line with our results, we showed significant sex-based differences in the association between TMB and outcome of patients with advanced NSCLC treated with anti-PD1/PDL1 drugs.

Conclusions

We demonstrated meaningful sex-based differences of anticancer immune response and immune evasion mechanisms, that may be exploited to improve immunotherapy efficacy for both women and men.

Decombinator V4: an improved AIRR compliant-software package for T-cell receptor sequence annotation?

MOTIVATION

Analysis of the T-cell receptor repertoire is rapidly entering the general toolbox used by researchers interested in cellular immunity. The annotation of T-cell receptors (TCRs) from raw sequence data poses specific challenges, which arise from the fact that TCRs are not germline encoded, and because of the stochastic nature of the generating process.

RESULTS

In this study, we report the release of Decombinator V4, a tool for the accurate and fast annotation of large sets of TCR sequences. Decombinator was one of the early Python software packages released to analyse the rapidly increasing flow of T-cell receptor repertoire sequence data. The Decombinator package now provides Python 3 compatibility, incorporates improved sequencing error and PCR bias correction algorithms, and provides output which conforms to the international standards proposed by the Adaptive Immune Receptor Repertoire Community.

AVAILABILITY AND IMPLEMENTATION

The entire Decombinator suite is freely available at: https://github.com/innate2adaptive/Decombinator.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Global analysis of shared T cell specificities in human non-small cell lung cancer enables HLA inference and antigen discovery

To identify disease-relevant T cell receptors (TCRs) with shared antigen specificity, we analyzed 778,938 TCRβ chain sequences from 178 non-small cell lung cancer patients using the GLIPH2 (grouping of lymphocyte interactions with paratope hotspots 2) algorithm. We identified over 66,000 shared specificity groups, of which 435 were clonally expanded and enriched in tumors compared to adjacent lung. The antigenic epitopes of one such tumor-enriched specificity group were identified using a yeast peptide-HLA A^∗02:01 display library. These included a peptide from the epithelial protein TMEM161A, which is overexpressed in tumors and cross-reactive epitopes from Epstein-Barr virus and E. coli. Our findings suggest that this cross-reactivity may underlie the presence of virus-specific T cells in tumor infiltrates and that pathogen cross-reactivity may be a feature of multiple cancers. The approach and analytical pipelines generated in this work, as well as the specificity groups defined here, present a resource for understanding the T cell response in cancer.

Persistent T Cell Repertoire Perturbation and T Cell Activation in HIV After Long Term Treatment

Objective

In people living with HIV (PLHIV), we sought to test the hypothesis that long term anti-retroviral therapy restores the normal T cell repertoire, and investigate the functional relationship of residual repertoire abnormalities to persistent immune system dysregulation.

Methods

We conducted a case-control study in PLHIV and HIV-negative volunteers, of circulating T cell receptor repertoires and whole blood transcriptomes by RNA sequencing, complemented by metadata from routinely collected health care records.

Results

T cell receptor sequencing revealed persistent abnormalities in the clonal T cell repertoire of PLHIV, characterized by reduced repertoire diversity and oligoclonal T cell expansion correlated with elevated CD8 T cell counts. We found no evidence that these expansions were driven by cytomegalovirus or another common antigen. Increased frequency of long CDR3 sequences and reduced frequency of public sequences among the expanded clones implicated abnormal thymic selection as a contributing factor. These abnormalities in the repertoire correlated with systems level evidence of persistent T cell activation in genome-wide blood transcriptomes.

Conclusions

The diversity of T cell receptor repertoires in PLHIV on long term anti-retroviral therapy remains significantly depleted, and skewed by idiosyncratic clones, partly attributable to altered thymic output and associated with T cell mediated chronic immune activation. Further investigation of thymic function and the antigenic drivers of T cell clonal selection in PLHIV are critical to efforts to fully re-establish normal immune function.

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

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

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.

Immune fingerprinting through repertoire similarity

Immune repertoires provide a unique fingerprint reflecting the immune history of individuals, with potential applications in precision medicine. However, the question of how personal that information is and how it can be used to identify individuals has not been explored. Here, we show that individuals can be uniquely identified from repertoires of just a few thousands lymphocytes. We present "Immprint," a classifier using an information-theoretic measure of repertoire similarity to distinguish pairs of repertoire samples coming from the same versus different individuals. Using published T-cell receptor repertoires and statistical modeling, we tested its ability to identify individuals with great accuracy, including identical twins, by computing false positive and false negative rates < 10-6 from samples composed of 10,000 T-cells. We verified through longitudinal datasets that the method is robust to acute infections and that the immune fingerprint is stable for at least three years. These results emphasize the private and personal nature of repertoire data.

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.

CRISPR-Mediated Base Conversion Allows Discriminatory Depletion of Endogenous T Cell Receptors for Enhanced Synthetic Immunity

Emerging base editing technology exploits CRISPR RNA-guided DNA modification effects for highly specific C > T conversion, which has been used to efficiently disrupt gene expression. These tools can enhance synthetic T cell immunity by restricting specificity, addressing histocompatibility leukocyte antigen (HLA) barriers, and promoting persistence. We report lentiviral delivery of a hepatitis B-virus (HBV)-specific recombinant T cell receptor (rTCR) and a linked CRISPR single-guide RNA for simultaneous disruption of endogenous TCRs (eTCRs) when combined with transient cytosine deamination. Discriminatory depletion of eTCR and coupled expression of rTCR resulted in enrichment of HBV-specific populations from 55% (SEM, ±2.4%) to 95% (SEM, ±0.5%). Intensity of rTCR expression increased 1.8- to 2.9-fold compared to that in cells retaining their competing eTCR, and increased cytokine production and killing of HBV antigen-expressing hepatoma cells in a 3D microfluidic model were exhibited. Molecular signatures confirmed that seamless conversion of C > T (G > A) had created a premature stop codon in TCR beta constant 1/2 loci, with no notable activity at predicted off-target sites. Thus, targeted disruption of eTCR by cytosine deamination and discriminatory enrichment of antigen-specific T cells offers the prospect of enhanced, more specific T cell therapies against HBV-associated hepatocellular carcinoma (HCC) as well as other viral and tumor antigens.

Combined TCR Repertoire Profiles and Blood Cell Phenotypes Predict Melanoma Patient Response to Personalized Neoantigen Therapy plus Anti-PD-1

T cells use highly diverse receptors (TCRs) to identify tumor cells presenting neoantigens arising from genetic mutations and establish anti-tumor activity. Immunotherapy harnessing neoantigen-specific T cells to target tumors has emerged as a promising clinical approach. To assess whether a comprehensive peripheral mononuclear blood cell analysis predicts responses to a personalized neoantigen cancer vaccine combined with anti-PD-1 therapy, we characterize the TCR repertoires and T and B cell frequencies in 21 patients with metastatic melanoma who received this regimen. TCR-α/β-chain sequencing reveals that prolonged progression-free survival (PFS) is strongly associated with increased clonal baseline TCR repertoires and longitudinal repertoire stability. Furthermore, the frequencies of antigen-experienced T and B cells in the peripheral blood correlate with repertoire characteristics. Analysis of these baseline immune features enables prediction of PFS following treatment. This method offers a pragmatic clinical approach to assess patients’ immune state and to direct therapeutic decision making.

Pre-treatment blood-based factors predict response to immunotherapy

TCR repertoire clonality and stability associate with improved clinical outcomes

Baseline T and B cell memory phenotypes associate with improved clinical outcomes

Combined baseline TCR repertoire and PBMC phenotypes predict immunotherapy response

Second-Strand Synthesis-Based Massively Parallel scRNA-Seq Reveals Cellular States and Molecular Features of Human Inflammatory Skin Pathologies

High-throughput single-cell RNA-sequencing (scRNA-seq) methodologies enable characterization of complex biological samples by increasing the number of cells that can be profiled contemporaneously. Nevertheless, these approaches recover less information per cell than low-throughput strategies. To accurately report the expression of key phenotypic features of cells, scRNA-seq platforms are needed that are both high fidelity and high throughput. To address this need, we created Seq-Well S3 ("Second-Strand Synthesis"), a massively parallel scRNA-seq protocol that uses a randomly primed second-strand synthesis to recover complementary DNA (cDNA) molecules that were successfully reverse transcribed but to which a second oligonucleotide handle, necessary for subsequent whole transcriptome amplification, was not appended due to inefficient template switching. Seq-Well S3 increased the efficiency of transcript capture and gene detection compared with that of previous iterations by up to 10- and 5-fold, respectively. We used Seq-Well S3 to chart the transcriptional landscape of five human inflammatory skin diseases, thus providing a resource for the further study of human skin inflammation.

Ultra-efficient sequencing of T Cell receptor repertoires reveals shared responses in muscle from patients with Myositis

BACKGROUND

Myositis, or idiopathic inflammatory myopathy (IIM), is a group disorders of unknown etiology characterized by the inflammation of skeletal muscle. The role of T cells and their antigenic targets in IIM initiation and progression is poorly understood. T cell receptor (TCR) repertoire sequencing is a powerful approach for characterizing complex T cell responses. However, current TCR sequencing methodologies are complex, expensive, or both, greatly limiting the scale of feasible studies.

METHODS

Here we present Framework Region 3 AmplifiKation sequencing ("FR3AK-seq"), a simplified multiplex PCR-based approach for the ultra-efficient and quantitative analysis of TCR complementarity determining region 3 (CDR3) repertoires. By using minimal primer sets targeting a conserved region immediately upstream of CDR3, undistorted amplicons are analyzed via short read, single-end sequencing. We also introduce the novel algorithm Inferring Sequences via Efficiency Projection and Primer Incorporation ("ISEPPI") for linking CDR3s to their associated variable genes.

FINDINGS

We find that FR3AK-seq is sensitive and quantitative, performing comparably to two different industry standards. FR3AK-seq and ISEPPI were used to efficiently and inexpensively characterize the T cell infiltrates of surgical muscle biopsies obtained from 145 patients with IIM and controls. A cluster of closely related TCRs was identified in samples from patients with sporadic inclusion body myositis (IBM).

INTERPRETATION

The ease and minimal cost of FR3AK-seq removes critical barriers to routine, large-scale TCR CDR3 repertoire analyses, thereby democratizing the quantitative assessment of human TCR repertoires in disease-relevant target tissues. Importantly, discovery of closely related TCRs in muscle from patients with IBM provides evidence for a shared antigen-driven T cell response in this disease of unknown pathogenesis.

FUNDING

This work was supported by NIH grant U24AI118633 and a Prostate Cancer Foundation Young Investigator Award.

The TCR Repertoire Reconstitution in Multiple Sclerosis: Comparing One-Shot and Continuous Immunosuppressive Therapies

Natalizumab (NTZ) and autologous hematopoietic stem cell transplantation (AHSCT) are two successful treatments for relapsing-remitting multiple sclerosis (RRMS), an autoimmune T-cell-driven disorder affecting the central nervous system that is characterized by relapses interspersed with periods of complete or partial recovery. Both RRMS treatments have been documented to impact T-cell subpopulations and the T-cell receptor (TCR) repertoire in terms of clone frequency, but, so far, the link between T-cell naive and memory populations, autoimmunity, and treatment outcome has not yet been established hindering insight into the post-treatment TCR landscape of MS patients. To address this important knowledge gap, we tracked peripheral T-cell subpopulations (naïve and memory CD4+ and CD8+) across 15 RRMS patients before and after two years of continuous treatment (NTZ) and a single treatment course (AHSCT) by high-throughput TCRß sequencing. We found that the two MS treatments left treatment-specific multidimensional traces in patient TCRß repertoire dynamics with respect to clonal expansion, clonal diversity and repertoire architecture. Comparing MS TCR sequences with published datasets suggested that the majority of public TCRs belonged to virus-associated sequences. In summary, applying multi-dimensional computational immunology to a TCRß dataset of treated MS patients, we show that qualitative changes of TCRß repertoires encode treatment-specific information that may be relevant for future clinical trials monitoring and personalized MS follow-up, diagnosis and treatment regimes. Natalizumab (NTZ) and autologous hematopoietic stem cell transplantation (AHSCT) are two successful treatments for relapsing-remitting multiple sclerosis (RRMS), an autoimmune T-cell-driven disorder affecting the central nervous system that is characterized by relapses interspersed with periods of complete or partial recovery. Both RRMS treatments have been documented to impact T-cell subpopulations and the T-cell receptor (TCR) repertoire in terms of clone frequency, but, so far, the link between T-cell naive and memory populations, autoimmunity, and treatment outcome has not yet been established hindering insight into the posttreatment TCR landscape of MS patients. To address this important knowledge gap, we tracked peripheral T-cell subpopulations (naive and memory CD4+ and CD8+) across 15 RRMS patients before and after 2 years of continuous treatment (NTZ) and a single treatment course (AHSCT) by high-throughput TCRβ sequencing. We found that the two MS treatments left treatment-specific multidimensional traces in patient TCRβ repertoire dynamics with respect to clonal expansion, clonal diversity, and repertoire architecture. Comparing MS TCR sequences with published datasets suggested that the majority of public TCRs belonged to virus-associated sequences. In summary, applying multidimensional computational immunology to a TCRβ dataset of treated MS patients, we show that qualitative changes of TCRβ repertoires encode treatment-specific information that may be relevant for future clinical trials monitoring and personalized MS follow-up, diagnosis, and treatment regimens.

Molecular T-Cell Repertoire Analysis as Source of Prognostic and Predictive Biomarkers for Checkpoint Blockade Immunotherapy

The T cells are key players of the response to checkpoint blockade immunotherapy (CBI) and monitoring the strength and specificity of antitumor T-cell reactivity remains a crucial but elusive component of precision immunotherapy. The entire assembly of T-cell receptor (TCR) sequences accounts for antigen specificity and strength of the T-cell immune response. The TCR repertoire hence represents a “footprint” of the conditions faced by T cells that dynamically evolves according to the challenges that arise for the immune system, such as tumor neo-antigenic load. Hence, TCR repertoire analysis is becoming increasingly important to comprehensively understand the nature of a successful antitumor T-cell response, and to improve the success and safety of current CBI.