Longitudinal immunosequencing in healthy people reveals persistent T cell receptors rich in highly public receptors.

BACKGROUND: The adaptive immune system maintains a diversity of T cells capable of recognizing a broad array of antigens. Each T cell's specificity for antigens is determined by its T cell receptors (TCRs), which together across all T cells form a repertoire of millions of unique receptors in each individual. Although many studies have examined how TCR repertoires change in response to disease or drugs, few have explored the temporal dynamics of the TCR repertoire in healthy individuals. RESULTS: Here we report immunosequencing of TCR ß chains (TCRß) from the blood of three healthy individuals at eight time points over one year. TCRß repertoires of all peripheral-blood T cells and sorted memory T cells clustered clearly by individual, systematically demonstrating that TCRß repertoires are specific to individuals across time. This individuality was absent from TCRßs from naive T cells, suggesting that the differences resulted from an individual's antigen exposure history, not genetic background. Many characteristics of the TCRß repertoire (e.g., diversity, clonality) were stable across time, although we found evidence of T cell expansion dynamics even within healthy individuals. We further identified a subset of "persistent" TCRßs present across all time points. These receptors were rich in clonal and highly public receptors and may play a key role in immune system maintenance. CONCLUSIONS: Our results highlight the importance of longitudinal sampling of the immune system, providing a much-needed baseline for TCRß dynamics in healthy individuals. Such a baseline will improve interpretation of changes in the TCRß repertoire during disease or treatment.

High-throughput sequencing of T-cell receptors reveals a homogeneous repertoire of tumour-infiltrating lymphocytes in ovarian cancer.

The cellular adaptive immune system mounts a response to many solid tumours mediated by tumour-infiltrating T lymphocytes (TILs). Basic measurements of these TILs, including total count, show promise as prognostic markers for a variety of cancers, including ovarian and colorectal. In addition, recent therapeutic advances are thought to exploit this immune response to effectively fight melanoma, with promising studies showing efficacy in additional cancers. However, many of the basic properties of TILs are poorly understood, including specificity, clonality, and spatial heterogeneity of the T-cell response. We utilize deep sequencing of rearranged T-cell receptor beta (TCRB) genes to characterize the basic properties of TILs in ovarian carcinoma. Due to somatic rearrangement during T-cell development, the TCR beta chain sequence serves as a molecular tag for each T-cell clone. Using these sequence tags, we assess similarities and differences between infiltrating T cells in discretely sampled sections of large tumours and compare to T cells from peripheral blood. Within the limits of sensitivity of our assay, the TIL repertoires show strong similarity throughout each tumour and are distinct from the circulating T-cell repertoire. We conclude that the cellular adaptive immune response within ovarian carcinomas is spatially homogeneous and distinct from the T-cell compartment of peripheral blood.

Tumor-infiltrating lymphocytes in colorectal tumors display a diversity of T cell receptor sequences that differ from the T cells in adjacent mucosal tissue.

Tumors from colorectal cancer (CRC) are generally immunogenic and commonly infiltrated with T lymphocytes. However, the details of the adaptive immune reaction to these tumors are poorly understood. We have accrued both colon tumor samples and adjacent healthy mucosal samples from 15 CRC patients to study lymphocytes infiltrating these tissues. We apply a method for detailed sequencing of T-cell receptor (TCR) sequences from tumor-infiltrating lymphocytes (TILs) in CRC tumors at high throughput to probe T-cell clones in comparison with the TCRs from adjacent healthy mucosal tissue. In parallel, we captured TIL counts using standard immunohistochemistry. The variation in diversity of the TIL repertoire was far wider than the variation of T-cell clones in the healthy mucosa, and the oligoclonality was higher on average in the tumors. However, the diversity of the T-cell repertoire in both CRC tumors and healthy mucosa was on average 100-fold lower than in peripheral blood. Using the TCR sequences to identify and track clones between mucosal and tumor samples, we determined that the immune response in the tumor is different than in the adjacent mucosal tissue, and the number of shared clones is not dependent on distance between the samples. Together, these data imply that CRC tumors induce a specific adaptive immune response, but that this response differs widely in strength and breadth between patients.

A Public Database of Memory and Naive B-Cell Receptor Sequences.

The vast diversity of B-cell receptors (BCR) and secreted antibodies enables the recognition of, and response to, a wide range of epitopes, but this diversity has also limited our understanding of humoral immunity. We present a public database of more than 37 million unique BCR sequences from three healthy adult donors that is many fold deeper than any existing resource, together with a set of online tools designed to facilitate the visualization and analysis of the annotated data. We estimate the clonal diversity of the naive and memory B-cell repertoires of healthy individuals, and provide a set of examples that illustrate the utility of the database, including several views of the basic properties of immunoglobulin heavy chain sequences, such as rearrangement length, subunit usage, and somatic hypermutation positions and dynamics.

T-cell receptor sequencing reveals the clonal diversity and overlap of colonic effector and FOXP3+ T cells in ulcerative colitis.

BACKGROUND: FOXP3 regulatory T cell prevent inflammation but are paradoxically increased in ulcerative colitis (UC). Local T-cell activation has been hypothesized to account for increased FOXP3 expression in colon lamina propria (LP) T cells. METHODS: To see if human FOXP3 LP T cells are an activated fraction of otherwise FOXP3 effector T cells and explore their clonal diversity in health and disease, we deep sequenced clonally unique T-cell receptor hypervariable regions of FOXP3 and FOXP3CD4 T-cell subpopulations from inflamed versus noninflamed colon LP or mesenteric lymph nodes of patients with or without UC. RESULTS: The clonal diversity of each LP T-cell population was not different between patients with versus without UC. Repertoire overlap was only seen between a minority of FOXP3 and FOXP3 cells, including recently activated CD38 cells and Th17-like CD161 effector T cells, but this repertoire overlap was not different between patients with versus without UC and was no larger than the overlap between Helios and Helios FOXP3 cells. CONCLUSIONS: Thus, at steady state, only a minority of FOXP3, and particularly Helios, T cells share a T-cell receptor sequence with FOXP3 effector populations in the colon LP, even in UC, revealing distinct clonal origins for LP regulatory T cell and effector T cells in humans.

Annotation of pseudogenic gene segments by massively parallel sequencing of rearranged lymphocyte receptor loci.

BACKGROUND: The adaptive immune system generates a remarkable range of antigen-specific T-cell receptors (TCRs), allowing the recognition of a diverse set of antigens. Most of this diversity is encoded in the complementarity determining region 3 (CDR3) of the ß chain of the αß TCR, which is generated by somatic recombination of noncontiguous variable (V), diversity (D), and joining (J) gene segments. Deletion and non-templated insertion of nucleotides at the D-J and V-DJ junctions further increases diversity. Many of these gene segments are annotated as non-functional owing to defects in their primary sequence, the absence of motifs necessary for rearrangement, or chromosomal locations outside the TCR locus. METHODS: We sought to utilize a novel method, based on high-throughput sequencing of rearranged TCR genes in a large cohort of individuals, to evaluate the use of functional and non-functional alleles. We amplified and sequenced genomic DNA from the peripheral blood of 587 healthy volunteers using a multiplexed polymerase chain reaction assay that targets the variable region of the rearranged TCRß locus, and we determined the presence and the proportion of productive rearrangements for each TCRß V gene segment in each individual. We then used this information to annotate the functional status of TCRß V gene segments in this cohort. RESULTS: For most TCRß V gene segments, our method agrees with previously reported functional annotations. However, we identified novel non-functional alleles for several gene segments, some of which were used exclusively in our cohort to the detriment of reported functional alleles. We also saw that some gene segments reported to have both functional and non-functional alleles consistently behaved in our cohort as either functional or non-functional, suggesting that some reported alleles were not present in the population studied. CONCLUSIONS: In this proof-of-principle study, we used high-throughput sequencing of the TCRß locus of a large cohort of healthy volunteers to evaluate the use of functional and non-functional alleles of individual TCRß V gene segments. With some modifications, our method has the potential to be extended to gene segments in the α, γ, and δ TCR loci, as well as the genes encoding for B-cell receptor chains.

High-throughput pairing of T cell receptor α and ß sequences.

The T cell receptor (TCR) protein is a heterodimer composed of an α chain and a ß chain. TCR genes undergo somatic DNA rearrangements to generate the diversity of T cell binding specificities needed for effective immunity. Recently, high-throughput immunosequencing methods have been developed to profile the TCR α (TCRA) and TCR ß (TCRB) repertoires. However, these methods cannot determine which TCRA and TCRB chains combine to form a specific TCR, which is essential for many functional and therapeutic applications. We describe and validate a method called pairSEQ, which can leverage the diversity of TCR sequences to accurately pair hundreds of thousands of TCRA and TCRB sequences in a single experiment. Our TCR pairing method uses standard laboratory consumables and equipment without the need for single-cell technologies. We show that pairSEQ can be applied to T cells from both blood and solid tissues, such as tumors.

Using synthetic templates to design an unbiased multiplex PCR assay.

T and B cell receptor loci undergo combinatorial rearrangement, generating a diverse immune receptor repertoire, which is vital for recognition of potential antigens. Here we use a multiplex PCR with a mixture of primers targeting the rearranged variable and joining segments to capture receptor diversity. Differential hybridization kinetics can introduce significant amplification biases that alter the composition of sequence libraries prepared by multiplex PCR. Using a synthetic immune receptor repertoire, we identify and minimize such biases and computationally remove residual bias after sequencing. We apply this method to a multiplex T cell receptor gamma sequencing assay. To demonstrate accuracy in a biological setting, we apply the method to monitor minimal residual disease in acute lymphoblastic leukaemia patients. A similar methodology can be extended to any adaptive immune locus.

Deep sequencing of the human TCRγ and TCRß repertoires suggests that TCRß rearranges after αß and γδ T cell commitment.

T lymphocytes respond to a broad array of pathogens with the combinatorial diversity of the T cell receptor (TCR). This adaptive response is possible because of the unique structure of the TCR, which is composed of two chains, either αß or γδ, that undergo genetic rearrangement in the thymus. αß and γδ T cells are functionally distinct within the host but are derived from a common multipotent precursor. The canonical model for T cell lineage commitment assumes that the γ, δ, and ß chains rearrange before αß or γδ T cell commitment. To test the standard model in humans, we used high-throughput sequencing to catalog millions of TCRγ and TCRß chains from peripheral blood αß and γδ T cells from three unrelated individuals. Almost all sampled αß and γδ T cells had rearranged TCRγ sequences. Although sampled αß T cells had a diverse repertoire of rearranged TCRß chains, less than 4% of γδ T cells in peripheral blood had a rearranged TCRß chain. Our data suggest that TCRγ rearranges in all T lymphocytes, consistent with TCRγ rearranging before T cell lineage commitment. However, rearrangement of the TCRß locus appears to be restricted after T cell precursors commit to the αß T cell lineage. Indeed, in T cell leukemias and lymphomas, TCRγ is almost always rearranged and TCRß is only rearranged in a subset of cancers. Because high-throughput sequencing of TCRs is translated into the clinic for monitoring minimal residual for leukemia/lymphoma, our data suggest the sequencing target should be TCRγ.

The major clades of Loasaceae: phylogenetic analysis using the plastid matK and trnL-trnF regions.

Phylogenetic analyses of Loasaceae that apply DNA sequence data from the plastid trnL-trnF region and matK gene in both maximum-parsimony and maximum-likelihood searches are presented. The results place subfamily Loasoideae as the sister of a subfamily Gronovioideae-Mentzelia clade. Schismocarpus is the sister of the Loasoideae-Gronovioideae-Mentzelia clade. The Schismocarpus-Loasoideae-Gronovioideae-Mentzelia clade is the sister of Eucnide. Several clades in Loasoideae receive strong support, providing insights on generic circumscription problems. Within Mentzelia, several major clades receive strong support, which clarifies relationships among previously circumscribed sections. Prior taxonomic and phylogenetic hypotheses are modeled using topology constraints in parsimony and likelihood analyses; tree lengths and likelihoods, respectively, are compared from constrained and unconstrained analyses to evaluate the relative support for various hypotheses. We use the Shimodaira-Hasegawa (SH) test to establish the significance of the differences between constrained and unconstrained topologies. The SH test rejects topologies based on hypotheses for (1) the placement of gronovioids as the sister of the rest of Loasaceae, (2) the monophyly of subfamily Mentzelioideae as well as Gronovioideae and Loasoideae, (3) the monophyly of Loasa sensu lato as circumscribed by Urban and Gilg, and (4) the monophyly of Mentzelia torreyi and Mentzelia sect. Bartonia.