A novel unconventional T cell population enriched in Crohn's disease.

OBJECTIVE: One of the current hypotheses to explain the proinflammatory immune response in IBD is a dysregulated T cell reaction to yet unknown intestinal antigens. As such, it may be possible to identify disease-associated T cell clonotypes by analysing the peripheral and intestinal T-cell receptor (TCR) repertoire of patients with IBD and controls. DESIGN: We performed bulk TCR repertoire profiling of both the TCR alpha and beta chains using high-throughput sequencing in peripheral blood samples of a total of 244 patients with IBD and healthy controls as well as from matched blood and intestinal tissue of 59 patients with IBD and disease controls. We further characterised specific T cell clonotypes via single-cell RNAseq. RESULTS: We identified a group of clonotypes, characterised by semi-invariant TCR alpha chains, to be significantly enriched in the blood of patients with Crohn's disease (CD) and particularly expanded in the CD8+ T cell population. Single-cell RNAseq data showed an innate-like phenotype of these cells, with a comparable gene expression to unconventional T cells such as mucosal associated invariant T and natural killer T (NKT) cells, but with distinct TCRs. CONCLUSIONS: We identified and characterised a subpopulation of unconventional Crohn-associated invariant T (CAIT) cells. Multiple evidence suggests these cells to be part of the NKT type II population. The potential implications of this population for CD or a subset thereof remain to be elucidated, and the immunophenotype and antigen reactivity of CAIT cells need further investigations in future studies.

The energy sensor AMPK orchestrates metabolic and translational adaptation in expanding T helper cells.

The importance of cellular metabolic adaptation in inducing robust T cell responses is well established. However, the mechanism by which T cells link information regarding nutrient supply to clonal expansion and effector function is still enigmatic. Herein, we report that the metabolic sensor adenosine monophosphate-activated protein kinase (AMPK) is a critical link between cellular energy demand and translational activity and, thus, orchestrates optimal expansion of T cells in vivo. AMPK deficiency did not affect T cell fate decision, activation, or T effector cell generation; however, the magnitude of T cell responses in murine in vivo models of T cell activation was markedly reduced. This impairment was global, as all T helper cell subsets were similarly sensitive to loss of AMPK which resulted in reduced T cell accumulation in peripheral organs and reduced disease severity in pathophysiologically as diverse models as T cell transfer colitis and allergic airway inflammation. T cell receptor repertoire analysis confirmed similar clonotype frequencies in different lymphoid organs, thereby supporting the concept of a quantitative impairment in clonal expansion rather than a skewed qualitative immune response. In line with these findings, in-depth metabolic analysis revealed a decrease in T cell oxidative metabolism, and gene set enrichment analysis indicated a major reduction in ribosomal biogenesis and mRNA translation in AMPK-deficient T cells. We, thus, provide evidence that through its interference with these delicate processes, AMPK orchestrates the quantitative, but not the qualitative, manifestation of primary T cell responses in vivo.

Local fitness landscape of the green fluorescent protein.

Fitness landscapes depict how genotypes manifest at the phenotypic level and form the basis of our understanding of many areas of biology, yet their properties remain elusive. Previous studies have analysed specific genes, often using their function as a proxy for fitness, experimentally assessing the effect on function of single mutations and their combinations in a specific sequence or in different sequences. However, systematic high-throughput studies of the local fitness landscape of an entire protein have not yet been reported. Here we visualize an extensive region of the local fitness landscape of the green fluorescent protein from Aequorea victoria (avGFP) by measuring the native function (fluorescence) of tens of thousands of derivative genotypes of avGFP. We show that the fitness landscape of avGFP is narrow, with 3/4 of the derivatives with a single mutation showing reduced fluorescence and half of the derivatives with four mutations being completely non-fluorescent. The narrowness is enhanced by epistasis, which was detected in up to 30% of genotypes with multiple mutations and mostly occurred through the cumulative effect of slightly deleterious mutations causing a threshold-like decrease in protein stability and a concomitant loss of fluorescence. A model of orthologous sequence divergence spanning hundreds of millions of years predicted the extent of epistasis in our data, indicating congruence between the fitness landscape properties at the local and global scales. The characterization of the local fitness landscape of avGFP has important implications for several fields including molecular evolution, population genetics and protein design.

The Absence of Retroelement Activity Is Characteristic for Childhood Acute Leukemias and Adult Acute Lymphoblastic Leukemia.

Retroelements (RE) have been proposed as important players in cancerogenesis. Different cancer types are characterized by a different level of tumor-specific RE insertions. In previous studies, small cohorts of hematological malignancies, such as acute myeloid leukemia, multiple myeloma, and chronic lymphocytic leukemia have been characterized by a low level of RE insertional activity. Acute lymphoblastic leukemia (ALL) in adults and childhood acute leukemias have not been studied in this context. We performed a search for new RE insertions (Alu and L1) in 44 childhood ALL, 14 childhood acute myeloid leukemia, and 14 adult ALL samples using a highly sensitive NGS-based approach. First, we evaluated the method sensitivity revealing the 1% detection threshold for the proportion of cells with specific RE insertion. Following this result, we did not identify new tumor-specific RE insertions in the tested cohort of acute leukemia samples at the established level of sensitivity. Additionally, we analyzed the transcription levels of active L1 copies and found them increased. Thus, the increased transcription of active L1 copies is not sufficient for overt elevation of L1 retrotranspositional activity in leukemia.

Lineage Conversion in Pediatric B-Cell Precursor Acute Leukemia under Blinatumomab Therapy.

We report incidence and deep molecular characteristics of lineage switch in 182 pediatric patients affected by B-cell precursor acute lymphoblastic leukemia (BCP-ALL), who were treated with blinatumomab. We documented six cases of lineage switch that occurred after or during blinatumomab exposure. Therefore, lineage conversion was found in 17.4% of all resistance cases (4/27) and 3.2% of relapses (2/63). Half of patients switched completely from BCP-ALL to CD19-negative acute myeloid leukemia, others retained CD19-positive B-blasts and acquired an additional CD19-negative blast population: myeloid or unclassifiable. Five patients had KMT2A gene rearrangements; one had TCF3::ZNF384 translocation. The presented cases showed consistency of gene rearrangements and fusion transcripts across initially diagnosed leukemia and lineage switch. In two of six patients, the clonal architecture assessed by IG/TR gene rearrangements was stable, while in others, loss of clones or gain of new clones was noted. KMT2A-r patients demonstrated very few additional mutations, while in the TCF3::ZNF384 case, lineage switch was accompanied by a large set of additional mutations. The immunophenotype of an existing leukemia sometimes changes via different mechanisms and with different additional molecular changes. Careful investigation of all BM compartments together with all molecular -minimal residual disease studies can lead to reliable identification of lineage switch.

Precise tracking of vaccine-responding T cell clones reveals convergent and personalized response in identical twins.

T cell receptor (TCR) repertoire data contain information about infections that could be used in disease diagnostics and vaccine development, but extracting that information remains a major challenge. Here we developed a statistical framework to detect TCR clone proliferation and contraction from longitudinal repertoire data. We applied this framework to data from three pairs of identical twins immunized with the yellow fever vaccine. We identified 600 to 1,700 responding TCRs in each donor and validated them using three independent assays. While the responding TCRs were mostly private, albeit with higher overlap between twins, they could be well-predicted using a classifier based on sequence similarity. Our method can also be applied to samples obtained postinfection, making it suitable for systematic discovery of new infection-specific TCRs in the clinic.

High-throughput sequencing of T-cell receptor alpha chain clonal rearrangements at the DNA level in lymphoid malignancies.

Rearrangements of T- and B-cell receptor (TCR and BCR) genes are useful markers for clonality assessment as well as for minimal residual disease (MRD) monitoring during the treatment of haematological malignancies. Currently, rearrangements of three out of four TCR and all BCR loci are used for this purpose. The fourth TCR gene, TRA, has not been used so far due to the lack of a method for its rearrangement detection in genomic DNA. Here we propose the first high-throughput sequencing based method for the identification of clonal TRA gene rearrangements at the DNA level. The method is based on target amplification of the rearranged TRA locus using an advanced multiplex polymerase chain reaction system and high-throughput sequencing, and has been tested on DNA samples from peripheral blood of healthy donors. Combinations of all functional V- and J-segments were detected, indicating the high sensitivity of the method. Additionally, we identified clonal TRA rearrangements in 57 out of 112 tested DNA samples of patients with various T-lineage lymphoproliferative disorders. The method fills the existing gap in utilizing the TRA gene for a wide range of studies, including clonality assessment, MRD monitoring and clonal evolution analysis in different lymphoid malignancies.

T-cell receptor and B-cell receptor repertoire profiling in adaptive immunity.

B-cell receptors and T-cell receptors are the key molecules responsible for specific antigen recognition in adaptive immunity. The huge diversity of immune receptor repertoires constrained their comprehensive studies in the past. More recently, however, high-throughput sequencing based techniques have revolutionized the field of immune receptor repertoire profiling enabling new insights into the development and function of the adaptive immune system. In this review we describe current methods for immune receptor profiling and software tools used for repertoire reconstruction from raw sequencing data. We also provide examples of how immune repertoire profiling can be used to study adaptive immunity in disease and in the course of organ and bone marrow transplantation.

Comparative analysis of murine T-cell receptor repertoires.

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

CD8+ T cells with characteristic T cell receptor beta motif are detected in blood and expanded in synovial fluid of ankylosing spondylitis patients.

Objective: The risk of AS is associated with genomic variants related to antigen presentation and specific cytokine signalling pathways, suggesting the involvement of cellular immunity in disease initiation/progression. The aim of the present study was to explore the repertoire of TCR sequences in healthy donors and AS patients to uncover AS-linked TCR variants. Methods: Using quantitative molecular-barcoded 5'-RACE, we performed deep TCR ß repertoire profiling of peripheral blood (PB) and SF samples for 25 AS patients and 108 healthy donors. AS-linked TCR variants were identified using a new computational approach that relies on a probabilistic model of the VDJ rearrangement process. Results: Using the donor-agnostic probabilistic model, we reveal a TCR ß motif characteristic for PB of AS patients, represented by eight highly homologous amino acid sequence variants. Some of these variants were previously reported in SF and PB of patients with ReA and in PB of AS patients. We demonstrate that identified AS-linked clones have a CD8+ phenotype, present at relatively low frequencies in PB, and are significantly enriched in matched SF samples of AS patients. Conclusion: Our results suggest the involvement of a particular antigen-specific subset of CD8+ T cells in AS pathogenesis, confirming and expanding earlier findings. The high similarity of the clonotypes with the ones found in ReA implies common mechanisms for the initiation of the diseases.

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

The diversity of T-cell receptors recognizing foreign pathogens is generated through a highly stochastic recombination process, making the independent production of the same sequence rare. Yet unrelated individuals do share receptors, which together constitute a "public" repertoire of abundant clonotypes. The TCR repertoire is initially formed prenatally, when the enzyme inserting random nucleotides is downregulated, producing a limited diversity subset. By statistically analyzing deep sequencing T-cell repertoire data from twins, unrelated individuals of various ages, and cord blood, we show that T-cell clones generated before birth persist and maintain high abundances in adult organisms for decades, slowly decaying with age. Our results suggest that large, low-diversity public clones are created during pre-natal life, and survive over long periods, providing the basis of the public repertoire.

Dynamics of Individual T Cell Repertoires: From Cord Blood to Centenarians.

The diversity, architecture, and dynamics of the TCR repertoire largely determine our ability to effectively withstand infections and malignancies with minimal mistargeting of immune responses. In this study, we have employed deep TCRß repertoire sequencing with normalization based on unique molecular identifiers to explore the long-term dynamics of T cell immunity. We demonstrate remarkable stability of repertoire, where approximately half of all T cells in peripheral blood are represented by clones that persist and generally preserve their frequencies for 3 y. We further characterize the extremes of lifelong TCR repertoire evolution, analyzing samples ranging from umbilical cord blood to centenarian peripheral blood. We show that the fetal TCR repertoire, albeit structurally maintained within regulated borders due to the lower numbers of randomly added nucleotides, is not limited with respect to observed functional diversity. We reveal decreased efficiency of nonsense-mediated mRNA decay in umbilical cord blood, which may reflect specific regulatory mechanisms in development. Furthermore, we demonstrate that human TCR repertoires are functionally more similar at birth but diverge during life, and we track the lifelong behavior of CMV- and EBV-specific T cell clonotypes. Finally, we reveal gender differences in dynamics of TCR diversity constriction, which come to naught in the oldest age. Based on our data, we propose a more general explanation for the previous observations on the relationships between longevity and immunity.

Towards error-free profiling of immune repertoires.

Deep profiling of antibody and T cell-receptor repertoires by means of high-throughput sequencing has become an attractive approach for adaptive immunity studies, but its power is substantially compromised by the accumulation of PCR and sequencing errors. Here we report MIGEC (molecular identifier groups-based error correction), a strategy for high-throughput sequencing data analysis. MIGEC allows for nearly absolute error correction while fully preserving the natural diversity of complex immune repertoires.

Distinctive properties of identical twins' TCR repertoires revealed by high-throughput sequencing.

Adaptive immunity in humans is provided by hypervariable Ig-like molecules on the surface of B and T cells. The final set of these molecules in each organism is formed under the influence of two forces: individual genetic traits and the environment, which includes the diverse spectra of alien and self-antigens. Here we assess the impact of individual genetic factors on the formation of the adaptive immunity by analyzing the T-cell receptor (TCR) repertoires of three pairs of monozygous twins by next-generation sequencing. Surprisingly, we found that an overlap between the TCR repertoires of monozygous twins is similar to an overlap between the TCR repertoires of nonrelated individuals. However, the number of identical complementary determining region 3 sequences in two individuals is significantly increased for twin pairs in the fraction of highly abundant TCR molecules, which is enriched by the antigen-experienced T cells. We found that the initial recruitment of particular TCR V genes for recombination and subsequent selection in the thymus is strictly determined by individual genetic factors. J genes of TCRs are selected randomly for recombination; however, the subsequent selection in the thymus gives preference to some α but not ß J segments. These findings provide a deeper insight into the mechanism of TCR repertoire generation.

Age-related decrease in TCR repertoire diversity measured with deep and normalized sequence profiling.

The decrease of TCR diversity with aging has never been studied by direct methods. In this study, we combined high-throughput Illumina sequencing with unique cDNA molecular identifier technology to achieve deep and precisely normalized profiling of TCR ß repertoires in 39 healthy donors aged 6-90 y. We demonstrate that TCR ß diversity per 10(6) T cells decreases roughly linearly with age, with significant reduction already apparent by age 40. The percentage of naive T cells showed a strong correlation with measured TCR diversity and decreased linearly up to age 70. Remarkably, the oldest group (average age 82 y) was characterized by a higher percentage of naive CD4(+) T cells, lower abundance of expanded clones, and increased TCR diversity compared with the previous age group (average age 62 y), suggesting the influence of age selection and association of these three related parameters with longevity. Interestingly, cross-analysis of individual TCR ß repertoires revealed a set >10,000 of the most representative public TCR ß clonotypes, whose abundance among the top 100,000 clones correlated with TCR diversity and decreased with aging.

tcR: an R package for T cell receptor repertoire advanced data analysis.

BACKGROUND: The Immunoglobulins (IG) and the T cell receptors (TR) play the key role in antigen recognition during the adaptive immune response. Recent progress in next-generation sequencing technologies has provided an opportunity for the deep T cell receptor repertoire profiling. However, a specialised software is required for the rational analysis of massive data generated by next-generation sequencing. RESULTS: Here we introduce tcR, a new R package, representing a platform for the advanced analysis of T cell receptor repertoires, which includes diversity measures, shared T cell receptor sequences identification, gene usage statistics computation and other widely used methods. The tool has proven its utility in recent research studies. CONCLUSIONS: tcR is an R package for the advanced analysis of T cell receptor repertoires after primary TR sequences extraction from raw sequencing reads. The stable version can be directly installed from The Comprehensive R Archive Network ( http://cran.r-project.org/mirrors.html ). The source code and development version are available at tcR GitHub ( http://imminfo.github.io/tcr/ ) along with the full documentation and typical usage examples.

Pairing of T-cell receptor chains via emulsion PCR.

Our ability to analyze adaptive immunity and engineer its activity has long been constrained by our limited ability to identify native pairs of heavy-light antibody chains and alpha-beta T-cell receptor (TCR) chains--both of which comprise coupled "halves of a key", collectively capable of recognizing specific antigens. Here, we report a cell-based emulsion RT-PCR approach that allows the selective fusion of the native pairs of amplified TCR alpha and beta chain genes for complex samples. A new type of PCR suppression technique was developed that makes it possible to amplify the fused library with minimal noise for subsequent analysis by high-throughput paired-end Illumina sequencing. With this technique, single analysis of a complex blood sample allows identification of multiple native TCR chain pairs. This approach may be extended to identify native antibody chain pairs and, more generally, pairs of mRNA molecules that are coexpressed in the same living cells.

Two-dimensional high-throughput on-cell screening of immunoglobulins against broad antigen repertoires.

Identifying high-affinity antibodies in human serum is challenging due to extremely low number of circulating B cells specific to the desired antigens. Delays caused by a lack of information on the immunogenic proteins of viral origin hamper the development of therapeutic antibodies. We propose an efficient approach allowing for enrichment of high-affinity antibodies against pathogen proteins with simultaneous epitope mapping, even in the absence of structural information about the pathogenic immunogens. To screen therapeutic antibodies from blood of recovered donors, only pathogen transcriptome is required to design an antigen polypeptide library, representing pathogen proteins, exposed on the bacteriophage surface. We developed a two-dimensional screening approach enriching lentiviral immunoglobulin libraries from the convalescent or vaccinated donors against bacteriophage library expressing the overlapping set of polypeptides covering the spike protein of SARS-CoV-2. This platform is suitable for pathogen-specific immunoglobulin enrichment and allows high-throughput selection of therapeutic human antibodies.

Next generation sequencing for TCR repertoire profiling: platform-specific features and correction algorithms.

The TCR repertoire is a mirror of the human immune system that reflects processes caused by infections, cancer, autoimmunity, and aging. Next generation sequencing (NGS) is becoming a powerful tool for deep TCR profiling; yet, questions abound regarding the methodological approaches for sample preparation and correct data interpretation. Accumulated PCR and sequencing errors along with library preparation bottlenecks and uneven PCR efficiencies lead to information loss, biased quantification, and generation of huge artificial TCR diversity. Here, we compare Illumina, 454, and Ion Torrent platforms for individual TCR profiling, evaluate the rate and character of errors, and propose advanced platform-specific algorithms to correct massive sequencing data. These developments are applicable to a wide variety of next generation sequencing applications. We demonstrate that advanced correction allows the removal of the majority of artificial TCR diversity with concomitant rescue of most of the sequencing information. Thus, this correction enhances the accuracy of clonotype identification and quantification as well as overall TCR diversity measurements.