Non-productive human TCR beta chain genes represent V-D-J diversity before selection upon function: insight into biased usage of TCRBD and TCRBJ genes and diversity of CDR3 region length

Composition and diversity analysis of the TCR CDR3 repertoire in patients with idiopathic orbital inflammation using high-throughput sequencing

BACKGROUND

Idiopathic orbital inflammation (IOI) is a nonspecific orbital inflammatory disease with the third highest prevalence among orbital diseases, and its pathogenesis is associated with T-cell-mediated immune responses. This study aimed to investigate the differences in T-cell receptor (TCR) expression between IOI patients and healthy subjects by high-throughput sequencing and to characterize TCR expression in patients with IOI and with respect to glucocorticoid response.

METHODS

A total of 19 subjects were enrolled in this study and were divided into the idiopathic orbital inflammation group (IOI group, n = 13) and the healthy control group (HC group, n = 6), and within the IOI group were further divided into the glucocorticoid therapy sensitive group (IOI(EF) group, n = 6) and the glucocorticoid therapy ineffective group (IOI(IN) group, n = 7) based on the degree of effectiveness to glucocorticoid therapy. High-throughput TCR sequencing was performed on peripheral blood mononuclear cells of IOI patients and healthy control individuals using 5' RACE technology combined with Unique Identifier (UID) digital tag correction technology. The TCR CDR3 region diversity, sharing patterns, and differential sequences between the IOI and HC groups, and between the IOI(EF) and IOI(IN) groups were analyzed.

RESULTS

It was found that the diversity of TCR CDR3 in the IOI group was significantly lower than that in the HC group, and the frequency of V gene use was significantly different between groups. The diversity of TCR CDR3 in patients in the IOI(EF) group was significantly lower than that in patients in the IOI(IN) group, and the frequency of V and J gene use was significantly different between the IOI(EF) group and the IOI(IN) group. Additionally, we found 133 nucleotide sequences shared in all IOI samples and screened two sequences with higher expression from them.

CONCLUSIONS

Our results suggested that abnormal clonal expansion of specific T-cells exists in IOI patients and that TCR diversity may had an impact on the prognosis of glucocorticoid-treated IOI. This study may contribute to a better understanding of the immune status of IOI and provide new insights for T-cell -associated IOI pathogenesis, diagnosis and treatment prediction.

T-cell tolerant fraction as a predictor of immune-related adverse events

BACKGROUND

Immune checkpoint inhibitor (ICI) therapies may cause unpredictable and potentially severe autoimmune toxicities termed immune-related adverse events (irAEs). Because T cells mediate ICI effects, T cell profiling may provide insight into the risk of irAEs. Here we evaluate a novel metric-the T-cell tolerant fraction-as a predictor of future irAEs.

METHODS

We examined T-cell receptor beta (TRB) locus sequencing from baseline pretreatment samples from an institutional registry and previously published studies. For each patient, we used TRB sequences to calculate the T-cell tolerant fraction, which was then assessed as a predictor of future irAEs (classified as Common Terminology Criteria for Adverse Event grade 0-1 vs grade ≥2). We then compared the tolerant fraction to TRB clonality and diversity. Finally, the tolerant fraction was assessed on (1) T cells enriched against napsin A, a potential autoantigen of irAEs; (2) thymic versus peripheral blood T cells; and (3) TRBs specific for various infections and autoimmune diseases.

RESULTS

A total of 77 patients with cancer (22 from an institutional registry and 55 from published studies) receiving ICI therapy (43 CTLA4, 19 PD1/PDL1, 15 combination CTLA4+PD1/PDL1) were included in the study. The tolerant fraction was significantly lower in cases with clinically significant irAEs (p<0.001) and had an area under the receiver operating curve (AUC) of 0.79. The tolerant fraction was lower for each ICI treatment category, reaching statistical significance for CTLA4 (p<0.001) and demonstrating non-significant trends for PD1/PDL1 (p=0.21) and combination ICI (p=0.18). The tolerant fraction for T cells enriched against napsin A was lower than other samples. The tolerant fraction was also lower in thymic versus peripheral blood samples, and lower in some (multiple sclerosis) but not other (type 1 diabetes) autoimmune diseases. In our study cohort, TRB clonality had an AUC of 0.62, and TRB diversity had an AUC of 0.60 for predicting irAEs.

CONCLUSIONS

Among patients receiving ICI, the baseline T-cell tolerant fraction may serve as a predictor of clinically significant irAEs.

The T cell receptor sequence influences the likelihood of T cell memory formation

T cell differentiation depends on activation through the T cell receptor (TCR), whose amino acid sequence varies cell to cell. Particular TCR amino acid sequences nearly guarantee Mucosal-Associated Invariant T (MAIT) and Natural Killer T (NKT) cell fates. To comprehensively define how TCR amino acids affects all T cell fates, we analyze the paired αβTCR sequence and transcriptome of 819,772 single cells. We find that hydrophobic CDR3 residues promote regulatory T cell transcriptional states in both the CD8 and CD4 lineages. Most strikingly, we find a set of TCR sequence features, concentrated in CDR2α, that promotes positive selection in the thymus as well as transition from naïve to memory in the periphery. Even among T cells that recognize the same antigen, these TCR sequence features help to explain which T cells form immunological memory, which is essential for effective pathogen response.

Two Step Selection for Bias in β Chain V-J Pairing

The β chain rearrangement in T cells is a two-step process where first Dβ and Jβ bind, and only then Vβ is joined to the complex. We here show that the frequency of human and mouse Vβ Jβ combinations deviates from the one expected based on each gene usage frequency. This bias is observed mainly in functional (F) rearrangements, but also slightly in non-functional (NF) rearrangements. Preferred Vβ Jβ combinations in F clones are shared between donors and samples, suggesting a common structural mechanism for these biases in addition to any host-specific antigen-induced peripheral selection. The sharing holds even in clones with Jβ1 that share the same Dβ1 gene. Vβ Jβ usage is correlated with the Molecular Weight and Isoelectric Point in F clones. The pairing is also observed in the Double Positive cells in mice thymocytes, suggesting that the selection leading to such a pairing occurs before thymic selection. These results suggest an additional structural checkpoint in the beta chain development prior to thymic selection during the T cell receptor expression. Understanding this structural selection is important for the distinction between normal and aberrant T cell development, and crucial for the design of engineered TCRs.

T cell receptor beta germline variability is revealed by inference from repertoire data

BACKGROUND

T and B cell receptor (TCR, BCR) repertoires constitute the foundation of adaptive immunity. Adaptive immune receptor repertoire sequencing (AIRR-seq) is a common approach to study immune system dynamics. Understanding the genetic factors influencing the composition and dynamics of these repertoires is of major scientific and clinical importance. The chromosomal loci encoding for the variable regions of TCRs and BCRs are challenging to decipher due to repetitive elements and undocumented structural variants.

METHODS

To confront this challenge, AIRR-seq-based methods have recently been developed for B cells, enabling genotype and haplotype inference and discovery of undocumented alleles. However, this approach relies on complete coverage of the receptors' variable regions, whereas most T cell studies sequence a small fraction of that region. Here, we adapted a B cell pipeline for undocumented alleles, genotype, and haplotype inference for full and partial AIRR-seq TCR data sets. The pipeline also deals with gene assignment ambiguities, which is especially important in the analysis of data sets of partial sequences.

RESULTS

From the full and partial AIRR-seq TCR data sets, we identified 39 undocumented polymorphisms in T cell receptor Beta V (TRBV) and 31 undocumented 5 ' UTR sequences. A subset of these inferences was also observed using independent genomic approaches. We found that a single nucleotide polymorphism differentiating between the two documented T cell receptor Beta D2 (TRBD2) alleles is strongly associated with dramatic changes in the expressed repertoire.

CONCLUSIONS

We reveal a rich picture of germline variability and demonstrate how a single nucleotide polymorphism dramatically affects the composition of the whole repertoire. Our findings provide a basis for annotation of TCR repertoires for future basic and clinical studies.

Reconstituting T cell receptor selection in-silico

Each T cell receptor (TCR) gene is created without regard for which substances (antigens) the receptor can recognize. T cell selection culls developing T cells when their TCRs (i) fail to recognize major histocompatibility complexes (MHCs) that act as antigen presenting platforms or (ii) recognize with high affinity self-antigens derived from healthy cells and tissue. While T cell selection has been thoroughly studied, little is known about which TCRs are retained or removed by this process. Therefore, we develop an approach using TCR gene sequencing and machine learning to identify patterns in TCR protein sequences influencing the outcome of T cell receptor selection. We verify the trained models classify TCRs from developing T cells as being before selection and TCRs from mature T cells as being after selection. Our approach may provide future avenues for studying the relationship between T cell selection and conditions like autoimmune diseases.

Selection influences naive CD8+ TCR-β repertoire sharing

Within each individual, the adaptive immune system generates a repertoire of cells expressing receptors capable of recognizing diverse potential pathogens. The theoretical diversity of the T-cell receptor (TCR) repertoire exceeds the actual size of the T-cell population in an individual by several orders of magnitude - making the observation of identical TCRs in different individuals extremely improbable if all receptors were equally likely. Despite this disparity between the theoretical and the realized diversity of the repertoire, these 'public' receptor sequences have been identified in autoimmune, cancer and pathogen interaction contexts. Biased generation processes explain the presence of public TCRs in the naive repertoire, but do not adequately explain the different abundances of these public TCRs. We investigate and characterize the distribution of genomic TCR-β sequences of naive CD8+ T cells from three genetically identical mice, comparing non-productive (non-functional sequences) and productive sequences. We find public TCR-β sequences at higher abundances compared with unshared sequences in the productive, but not in the non-productive, repertoire. We show that neutral processes such as recombination biases, codon degeneracy and generation probability do not fully account for these differences, and conclude that thymic or peripheral selection plays an important role in increasing the abundances of public TCR-β sequences.

High-throughput sequencing analysis of genes encoding the B-lymphocyte receptor heavy-chain CDR3 in renal and peripheral blood of IgA nephropathy

Aim: IgA nephropathy (IgAN) is one of the most common chronic glomerulonephritis. Its etiology and pathogenesis remain unclear. We thus explored the immune repertoire of the B-cell receptor (BCR) and the heavy-chain complementarity-determining region 3 (CDR3) in renal tissue and peripheral blood of IgAN patients. Method: Total RNAs extracted from renal tissues and peripheral blood of patients and peripheral blood of healthy controls (HCs) were analyzed via high-throughput multiplex PCR sequencing. We amplified and sequenced BCR heavy-chain CDR3 regions to explore repertoire diversity, V/J gene family distribution, CDR3 lengths, BCR heavy-chain variants, consistency between tissue and peripheral blood data, and clones ‘shared’ by these bodily compartments. Results: We identified the renal tissue and peripheral blood BCR heavy-chain CDR3 immune repertoires of 15 IgAN patients. Top1 could be more readily cloned from peripheral blood of patients than from controls (P<0.05), the average CDR3 length was significantly shorter in patients than in HCs (P<0.05), the variant frequency of the gene encoding the BCR heavy chain was higher in patients than in HCs (P<0.05), and the BCR variant frequency was highest in IgAN kidney tissue. Preliminary screening for ‘shared’ clones showed that, in at least 13 patients, the ‘ALYFHNSAY’, ‘ARWGPMYYYMDV’, ‘ARDQGALNA’, and ‘ARVDNPADF’ CDR3 sequences were evident in peripheral blood samples from patients, but not HCs. Conclusions: We found that the ‘ALYFHNSAY’, ‘ARWGPMYYYMDV’, ‘ARDQGALNA’, and ‘ARVDNPADF’ clonal sequences may be useful for noninvasive diagnosis and treatment planning in IgAN.

Clonotypic heterogeneity in cutaneous T-cell lymphoma (mycosis fungoides) revealed by comprehensive whole-exome sequencing

Mycosis fungoides (MF), the most common type of cutaneous T-cell lymphoma, is believed to represent a clonal expansion of a transformed skin-resident memory T cell. T-cell receptor (TCR) clonality (ie, identical sequences of rearranged TCRα, TCRβ, and TCRγ), the key premise of this hypothesis, has been difficult to document conclusively because malignant cells are not readily distinguishable from the tumor-infiltrating reactive lymphocytes that contribute to the TCR clonotypic repertoire of MF. Here, we have successfully adopted targeted whole-exome sequencing (WES) to identify the repertoire of rearranged TCR genes in tumor-enriched samples from patients with MF. Although some of the investigated MF biopsies had the expected frequency of monoclonal rearrangements of TCRγ corresponding to that of tumor cells, the majority of the samples presented multiple TCRγ, TCRα, and TCRβ clonotypes by WES. Our findings are compatible with the model in which the initial malignant transformation in MF does not occur in mature memory T cells but rather at the level of T-lymphocyte progenitors before TCRβ or TCRα rearrangements. We have also shown that WES can be combined with whole-transcriptome sequencing in the same sample, which enables comprehensive characterization of the TCR repertoire in relation to tumor content. WES/whole-transcriptome sequencing might be applicable to other types of T-cell lymphomas to determine clonal dominance and clonotypic heterogeneity in these malignancies.

MHC-II alleles shape the CDR3 repertoires of conventional and regulatory naïve CD4+ T cells

T cell maturation and activation depend upon T cell receptor (TCR) interactions with a wide variety of antigenic peptides displayed in a given major histocompatibility complex (MHC) context. Complementarity-determining region 3 (CDR3) is the most variable part of the TCRα and -β chains, which govern interactions with peptide-MHC complexes. However, it remains unclear how the CDR3 landscape is shaped by individual MHC context during thymic selection of naïve T cells. We established two mouse strains carrying distinct allelic variants of H2-A and analyzed thymic and peripheral production and TCR repertoires of naïve conventional CD4⁺ T (T_conv) and naïve regulatory CD4⁺ T (T_reg) cells. Compared with tuberculosis-resistant C57BL/6 (H2-A^b) mice, the tuberculosis-susceptible H2-A^j mice had fewer CD4⁺ T cells of both subsets in the thymus. In the periphery, this deficiency was only apparent for T_conv and was compensated for by peripheral reconstitution for T_reg We show that H2-A^j favors selection of a narrower and more convergent repertoire with more hydrophobic and strongly interacting amino acid residues in the middle of CDR3α and CDR3β, suggesting more stringent selection against a narrower peptide-MHC-II context. H2-A^j and H2-A^b mice have prominent reciprocal differences in CDR3α and CDR3β features, probably reflecting distinct modes of TCR fitting to MHC-II variants. These data reveal the mechanics and extent of how MHC-II shapes the naïve CD4⁺ T cell CDR3 landscape, which essentially defines adaptive response to infections and self-antigens.

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

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

Profiling of the TCRβ repertoire in non-model species using high-throughput sequencing

In recent years, immune repertoire profiling with high-throughput sequencing (HTS) has advanced our understanding of adaptive immunity. However, fast progress in the field applied mostly to human and mouse research, with only few studies devoted to other model vertebrates. We present the first in-depth characterization of the T-cell receptor (TCR) repertoire in a non-model mammal (bank vole, Myodes glareolus), widely used in ecological and evolutionary research. We used RNA from spleens, 5′RACE and HTS to describe V and J segments of TCRβ, qualitatively characterize preferential V–J segment usage and CDR3 length distribution. Overall orthology to murine genes was preserved, with 11 J and 37 V genes found in voles (although 3 V genes lacked a close orthologue). Further, we implemented unique molecular identifiers for quantitative analysis of CDR3 repertoire with stringent error correction. A conservative, lower bound estimation of the TCRβ repertoire was similar to that found for mice (1.7–2.3 × 10⁵ clonotypes). We hope that by providing an easy-to-follow molecular protocol and on-line bioinformatics tools that do not require reference sequences (AmpliTCR and AmpliCDR3), we will encourage HTS immune repertoire profiling in other non-model vertebrates, thus opening new research avenues in e.g. comparative immunology, ecology and evolutionary biology.

Compositional characteristics of human peripheral TRBV pseudogene rearrangements

The diversity of the T cell receptor (TCR) complementarity determining region 3 (CDR3) repertoire is the result of random combinations, insertions and deletions during recombination of the germline V, D and J gene fragments. During evolution, some human TCR beta chain variable (TRBV) pseudogenes have been retained. Many previous studies have focused on functional TRBV genes, while little attention has been given to TRBV pseudogenes. To describe the compositional characteristics of TRBV pseudogene rearrangements, we compared and analysed TRBV pseudogenes, TRBV open reading frames (ORFs) and functional TRBV genes via high-throughput sequencing of DNA obtained from the peripheral blood of 4 healthy volunteers and 4 patients. Our results revealed several differences in J and D gene usage. The V deletion distribution profile of the pseudogenes was significantly different from that of the ORFs and functional genes. In addition, arginine, lysine and cysteine were more frequently used in putative CDR3 pseudogene rearrangements, while functional rearrangements used more leucine. This study presents a comprehensive description of the compositional characteristics of peripheral TRBV pseudogene rearrangements, which will provide a reference for further research on TRBV pseudogenes.

Oligoclonal T-cell Receptor Repertoire in Colonic Biopsies of Patients with Microscopic Colitis and Ulcerative Colitis

BACKGROUND

Microscopic colitis (MC), comprising collagenous colitis (CC) and lymphocytic colitis (LC), is a type of variation of inflammatory bowel diseases. Local T-cell infiltration in the mucosa plays a major role in MC immunopathology.

METHODS

To understand diversity and clonality of infiltrating T cells, we analyzed the T-cell receptor beta (TCRβ) chains in colonic biopsies of MC, ulcerative colitis (UC), and their remission counterparts (CC/LC-HR [histological remission] or UC-R [remission]) compared with patients with noninflamed colons using next-generation sequencing.

RESULTS

Compared with controls and patients with CC, patients with LC had significantly lower diversity with significantly lower evenness and richness in TCRVβ-Jβ gene segments. Similarly, patients with LC-HR had lower diversity because of significantly lower TCRVβ-Jβ clone richness. Patients with UC and UC-R showed significantly higher diversity and richness. Univariate and multivariate analyses were performed to identify TCRVβ-Jβ gene segments differentiating disease types from controls or their remission counterparts. Patients with LC were discriminated from controls by 12 clones and from patients with CC by 8 clones. Neither univariate nor multivariate analyses showed significance for patients with CC or CC-HR compared with controls. Patients with UC and UC-R had 16 and 14 discriminating clones, respectively, compared with controls.

CONCLUSIONS

Altogether, patients with MC and UC showed an oligoclonal TCRβ distribution. TCRVβ-Jβ clone types and their diversity were distinctive between patients with CC and LC, as well as for patients with UC, suggesting different pathophysiological mechanisms according to disease type and stage. This study suggests that CC and LC are different entities because of differences in immunoregulatory responses, as mirrored by their T-cell repertoire.

Current status and recent advances of next generation sequencing techniques in immunological repertoire

To ward off a wide variety of pathogens, the human adaptive immune system harbors a vast array of T-cell receptors (TCRs) and B-cell receptors (BCRs), collectively referred to as the immune repertoire. High-throughput sequencing (HTS) of TCR/BCR genes allows in-depth molecular analysis of T/B-cell clones, providing an unprecedented level of detail when examining the T/B-cell repertoire of individuals. It can evaluate TCR/BCR complementarity-determining region 3 (CDR3) diversity and assess the clonal composition, including the size of the repertoire; similarities between repertoires; V(D)J segment use; nucleotide insertions and deletions; CDR3 lengths; and amino acid distributions along the CDR3s at sequence-level resolution. Deep sequencing of B-cell and T-cell repertoires offers the potential for a quantitative understanding of the adaptive immune system in healthy and disease states. Recently, paired sequencing strategies have also been developed, which can provide information about the identity of immune receptor pairs encoded by individual T or B lymphocytes. HTS technology provides a previously unimaginable amount of sequence data, accompanied, however, by numerous challenges associated with error correction and interpretation that remain to be solved. The review details some of the technologies and some of the recent achievements in this field.

Chromatin conformation governs T-cell receptor Jβ gene segment usage

T cells play fundamental roles in adaptive immunity, relying on a diverse repertoire of T-cell receptor (TCR) α and β chains. Diversity of the TCR β chain is generated in part by a random yet intrinsically biased combinatorial rearrangement of variable (Vβ), diversity (Dβ), and joining (Jβ) gene segments. The mechanisms that determine biases in gene segment use remain unclear. Here we show, using a high-throughput TCR sequencing approach, that a physical model of chromatin conformation at the DJβ genomic locus explains more than 80% of the biases in Jβ use that we measured in murine T cells. This model also predicts correctly how differences in intersegment genomic distances between humans and mice translate into differences in Jβ bias between TCR repertoires of these two species. As a consequence of these structural and other biases, TCR sequences are produced with different a priori frequencies, thus affecting their probability of becoming public TCRs that are shared among individuals. Surprisingly, we find that many more TCR sequences are shared among all five mice we studied than among only subgroups of three or four mice. We derive a necessary mathematical condition explaining this finding, which indicates that the TCR repertoire contains a core set of receptor sequences that are highly abundant among individuals, if their a priori probability of being produced by the recombination process is higher than a defined threshold. Our results provide evidence for an expanded role of chromatin conformation in VDJ rearrangement, from control of gene accessibility to precise determination of gene segment use.

Identification of errors introduced during high throughput sequencing of the T cell receptor repertoire

Background

Recent advances in massively parallel sequencing have increased the depth at which T cell receptor (TCR) repertoires can be probed by >3log10, allowing for saturation sequencing of immune repertoires. The resolution of this sequencing is dependent on its accuracy, and direct assessments of the errors formed during high throughput repertoire analyses are limited.

Results

We analyzed 3 monoclonal TCR from TCR transgenic, Rag^-/- mice using Illumina^® sequencing. A total of 27 sequencing reactions were performed for each TCR using a trifurcating design in which samples were divided into 3 at significant processing junctures. More than 20 million complementarity determining region (CDR) 3 sequences were analyzed. Filtering for lower quality sequences diminished but did not eliminate sequence errors, which occurred within 1-6% of sequences. Erroneous sequences were pre-dominantly of correct length and contained single nucleotide substitutions. Rates of specific substitutions varied dramatically in a position-dependent manner. Four substitutions, all purine-pyrimidine transversions, predominated. Solid phase amplification and sequencing rather than liquid sample amplification and preparation appeared to be the primary sources of error. Analysis of polyclonal repertoires demonstrated the impact of error accumulation on data parameters.

Conclusions

Caution is needed in interpreting repertoire data due to potential contamination with mis-sequence reads. However, a high association of errors with phred score, high relatedness of erroneous sequences with the parental sequence, dominance of specific nt substitutions, and skewed ratio of forward to reverse reads among erroneous sequences indicate approaches to filter erroneous sequences from repertoire data sets.

Comprehensive assessment of T-cell receptor beta-chain diversity in alphabeta T cells

The adaptive immune system uses several strategies to generate a repertoire of T- and B-cell antigen receptors with sufficient diversity to recognize the universe of potential pathogens. In alphabeta T cells, which primarily recognize peptide antigens presented by major histocompatibility complex molecules, most of this receptor diversity is contained within the third complementarity-determining region (CDR3) of the T-cell receptor (TCR) alpha and beta chains. Although it has been estimated that the adaptive immune system can generate up to 10(16) distinct alphabeta pairs, direct assessment of TCR CDR3 diversity has not proved amenable to standard capillary electrophoresis-based DNA sequencing. We developed a novel experimental and computational approach to measure TCR CDR3 diversity based on single-molecule DNA sequencing, and used this approach to determine the CDR3 sequence in millions of rearranged TCRbeta genes from T cells of 2 adults. We find that total TCRbeta receptor diversity is at least 4-fold higher than previous estimates, and the diversity in the subset of CD45RO(+) antigen-experienced alphabeta T cells is at least 10-fold higher than previous estimates. These methods should prove valuable for assessment of alphabeta T-cell repertoire diversity after hematopoietic cell transplantation, in states of congenital or acquired immunodeficiency, and during normal aging.

Profiling the T-cell receptor beta-chain repertoire by massively parallel sequencing

T-cell receptor (TCR) genomic loci undergo somatic V(D)J recombination, plus the addition/subtraction of nontemplated bases at recombination junctions, in order to generate the repertoire of structurally diverse T cells necessary for antigen recognition. TCR beta subunits can be unambiguously identified by their hypervariable CDR3 (Complement Determining Region 3) sequence. This is the site of V(D)J recombination encoding the principal site of antigen contact. The complexity and dynamics of the T-cell repertoire remain unknown because the potential repertoire size has made conventional sequence analysis intractable. Here, we use 5'-RACE, Illumina sequencing, and a novel short read assembly strategy to sample CDR3(beta) diversity in human T lymphocytes from peripheral blood. Assembly of 40.5 million short reads identified 33,664 distinct TCR(beta) clonotypes and provides precise measurements of CDR3(beta) length diversity, usage of nontemplated bases, sequence convergence, and preferences for TRBV (T-cell receptor beta variable gene) and TRBJ (T-cell receptor beta joining gene) gene usage and pairing. CDR3 length between conserved residues of TRBV and TRBJ ranged from 21 to 81 nucleotides (nt). TRBV gene usage ranged from 0.01% for TRBV17 to 24.6% for TRBV20-1. TRBJ gene usage ranged from 1.6% for TRBJ2-6 to 17.2% for TRBJ2-1. We identified 1573 examples of convergence where the same amino acid translation was specified by distinct CDR3(beta) nucleotide sequences. Direct sequence-based immunoprofiling will likely prove to be a useful tool for understanding repertoire dynamics in response to immune challenge, without a priori knowledge of antigen.

Evolutionarily conserved pattern of gene segment usage within the mammalian TCRbeta locus

Antigen receptor gene rearrangement is mediated by interactions between the VDJ recombinase and the recombination signal sequences that flank the antigen receptor gene segments. In this report I present phylogenetic analyses that suggest a remarkable evolutionary conservation of the recombination signal sequences flanking some of the orthologous T-cell receptor-beta locus gene segments between human and mouse. Comparison of published data on the usage of the same gene segments between human and mouse indicates similar conservation in the shape of the primary T-cell receptor-beta repertoire. I propose that interactions between the recombinase and its cognate recognition sequences play a hitherto underestimated role in the formation of the specific pattern of the primary, combinatorial antigen receptor repertoire and that this pattern appears to be conserved in diverse mammalian species. Generation of a conserved pattern of the primary T-cell receptor repertoire may be critical for efficient selection of immature T lymphocytes.

Assessing diversity: immune reconstitution and T-cell receptor BV spectratype analysis following stem cell transplantation

T-cell receptor (TCR) BV spectratyping provides information concerning immune reconstitution complementary to that derived from monoclonal antibody analysis of surface antigen expression. However, the most appropriate way to analyse and represent these data, the number of subfamilies that should be analysed, and the effects of CD4/8 ratio on observed diversity, remain unclear. A diversity scoring system was developed based on analysis of 11 cord blood (CB) and 12 normal adult BV spectratypes. CB subfamily spectratypes demonstrated minor deviations from a Gaussian pattern consistent with current knowledge about germline TCR rearrangements. Diversity scores were significantly lower in myeloma patients than normal adults (P < 0.001) and fell significantly following stem cell transplantation (P = 0.003). Subsequently increasing diversity was not detected by two previously described complexity scoring systems. The CD4/8 ratio was neither the major determinant of the absolute diversity score nor of the change in score for individual patients, suggesting that analysis of unselected mononuclear cells can provide information largely independent of CD4/8 ratios. There was a relatively low correlation between individual subfamily scores and overall diversity score. The novel and objective diversity scoring system described appears better able to document changes in spectratype patterns than previously described techniques and should aid the standardization of reporting.

T-cell receptor V? gene usage by T cells reactive with the tumor-rejection antigen SART-1 in oral squamous cell carcinoma

We recently described that the SART-1(690-698) peptide could induce HLA-A24-restricted cytotoxic T lymphocytes (CTLs), which recognize the SART-1(259) (+) tumor cells from peripheral blood mononuclear cells (PBMCs) of HLA-A24(+) cancer patients. In our study, in 5 of 14 HLA-A24(+) patients with oral squamous cell carcinomas (SCCs), CTLs could be induced with the SART-1(690-698) peptide from the PBMCs. In 2 of the patients from whom the highest CTL activities were induced, the T-cell receptor (TCR) Vbeta repertoire expressed by the SART-1(690-698)-specific CTLs was found to be restricted and multiple Vbeta families were predominantly expressed in each patient. Although the predominant Vbeta families were different between the 2 patients, Vbeta7 was highly and commonly predominant. The same predominant Vbeta families were also detected in the tumor-infiltrating lymphocytes (TILs) from each patient, and each Vbeta family contained one or more unique T-cell clonotypes. The unique T-cell clonotypes were found to be common between the TILs and SART-1(690-698)-specific CTLs from each patient, and especially 2 T-cell clonotypes with Vbeta7 were identical even in the 2 patients. One of the 2 T-cell clonotypes with Vbeta7 was detected in the TILs from 11 of 14 HLA-A24(+) patients and another was found in those from 8 of HLA-A24(+) patients, while none of 10 HLA-A24(-) patients demonstrated both T-cell clonotypes. These results strongly suggest that the T-cell clonotypes with Vbeta7 are major TCR Vbeta genes expressed by SART-1(690-698)-specific CTLs. Furthermore, autologous tumor cells from one of the HLA-A24(+) patients stimulated the PBMCs and regional lymph node cells (LNCs) to expand the same T-cell clonotypes as those in the SART-1(690-698)-specific CTLs. These results strongly suggest that the SART-1(690-698)-specific CTLs clearly accumulate in vivo, especially in the TILs, as a consequence of in situ antigenic stimulation by autologous tumor cells. The identification of the unique TCR Vbeta genes used by SART-1(259)-specific CTLs should help to improve the diagnosis of the specific immune response in patients with SART-1(259) (+) cancers, especially during anticancer immunotherapy.

Structural analysis of TCRalpha and beta chains from human T-Cell clones specific for small nuclear ribonucleoprotein polypeptides Sm-D, Sm-B and U1-70 kDa: TCR complementarity determining region 3 usage appears highly conserved

Systemic lupus erythematosus (SLE) and mixed connective tissue disease (MCTD) are systemic autoimmune diseases that are characterized by the presence of autoantibodies reactive with U small nuclear RNP (snRNP) autoantigens. Both B and T cells are important in the pathogenesis of the disease, and T- and B-cell immunity against snRNP polypeptides have been shown to be linked in vivo. Currently, several alternative hypotheses for the pathogenesis of these diseases have been proposed. These include loss of tolerance, modified self-antigens, molecular mimicry and nondirected immune activation. To help distinguish between the various models of disease pathogenesis, we have characterized the T-cell receptor (TCR) CDR3 from a large panel of well-characterized human T-cell clones and lines specific for individual snRNP polypeptides. The results presented here reveal highly restricted TCR usage across patients by the snRNP-reactive T cells based on the deduced amino acid sequence of the CDR3 loop. These data support the hypothesis that T-cell responses against self antigens in SLE and MCTD are antigen driven and that there are a limited number of T-cell epitopes present on the snRNP autoantigens.