TIRTL-seq: Deep, quantitative, and affordable paired TCR repertoire sequencing

ɑ/β T cells are key players in adaptive immunity. The specificity of T cells is determined by the sequences of the hypervariable T cell receptor (TCR) ɑ and β chains. Although bulk TCR sequencing offers a cost-effective approach for in-depth TCR repertoire profiling, it does not provide chain pairings, which are essential for determining T cell specificity. In contrast, single-cell TCR sequencing technologies produce paired chain data, but are limited in throughput to thousands of cells and are cost-prohibitive for cohort-scale studies. Here, we present TIRTL-seq (Throughput-Intensive Rapid TCR Library sequencing), a novel approach that generates ready-to-sequence TCR libraries from live cells in less than 7 hours. The protocol is optimized for use with non-contact liquid handlers in an automation-friendly 384-well plate format. Reaction volume miniaturization reduces library preparation costs to <$0.50 per well. The core principle of TIRTL-seq is the parallel generation of hundreds of libraries providing multiple biological replicates from a single sample that allows precise inference of both frequencies of individual clones and TCR chain pairings from well-occurrence patterns. We demonstrate scalability of our approach up to 1 million unique paired αβTCR clonotypes corresponding to over 30 million T cells per sample at a cost of less than $2000. For a sample of 10 million cells the cost is ~$200. We benchmarked TIRTL-seq against state-of-the-art 5'RACE bulk TCR-seq and 10x Genomics Chromium technologies on longitudinal samples. We show that TIRTL-seq is able to quantitatively identify expanding and contracting clonotypes between timepoints while providing accurate TCR chain pairings, including distinct temporal dynamics of SARS-CoV-2-specific and EBV-specific CD8+ T cell responses after infection. While clonal expansion was followed by sharp contraction for SARS-CoV-2 specific TCRs, EBV-specific TCRs remained stable once established. The sequences of both ɑ and β TCR chains are essential for determining T cell specificity. As the field moves towards greater applications in diagnostics and immunotherapy that rely on TCR specificity, we anticipate that our scalable paired TCR sequencing methodology will be instrumental for collecting large paired-chain datasets and ultimately extracting therapeutically relevant information from the TCR repertoire.

Influenza vaccination stimulates maturation of the human T follicular helper cell response

The differentiation and specificity of human CD4+ T follicular helper cells (TFH cells) after influenza vaccination have been poorly defined. Here we profiled blood and draining lymph node (LN) samples from human volunteers for over 2 years after two influenza vaccines were administered 1 year apart to define the evolution of the CD4+ TFH cell response. The first vaccination induced an increase in the frequency of circulating TFH (cTFH) and LN TFH cells at week 1 postvaccination. This increase was transient for cTFH cells, whereas the LN TFH cells further expanded during week 2 and remained elevated in frequency for at least 3 months. We observed several distinct subsets of TFH cells in the LN, including pre-TFH cells, memory TFH cells, germinal center (GC) TFH cells and interleukin-10+ TFH cell subsets beginning at baseline and at all time points postvaccination. The shift toward a GC TFH cell phenotype occurred with faster kinetics after the second vaccine compared to the first vaccine. We identified several influenza-specific TFH cell clonal lineages, including multiple responses targeting internal influenza virus proteins, and found that each TFH cell state was attainable within a clonal lineage. Thus, human TFH cells form a durable and dynamic multitissue network.

Benchmarking and integrating human B-cell receptor genomic and antibody proteomic profiling

Immunoglobulins (Ig), which exist either as B-cell receptors (BCR) on the surface of B cells or as antibodies when secreted, play a key role in the recognition and response to antigenic threats. The capability to jointly characterize the BCR and antibody repertoire is crucial for understanding human adaptive immunity. From peripheral blood, bulk BCR sequencing (bulkBCR-seq) currently provides the highest sampling depth, single-cell BCR sequencing (scBCR-seq) allows for paired chain characterization, and antibody peptide sequencing by tandem mass spectrometry (Ab-seq) provides information on the composition of secreted antibodies in the serum. Yet, it has not been benchmarked to what extent the datasets generated by these three technologies overlap and complement each other. To address this question, we isolated peripheral blood B cells from healthy human donors and sequenced BCRs at bulk and single-cell levels, in addition to utilizing publicly available sequencing data. Integrated analysis was performed on these datasets, resolved by replicates and across individuals. Simultaneously, serum antibodies were isolated, digested with multiple proteases, and analyzed with Ab-seq. Systems immunology analysis showed high concordance in repertoire features between bulk and scBCR-seq within individuals, especially when replicates were utilized. In addition, Ab-seq identified clonotype-specific peptides using both bulk and scBCR-seq library references, demonstrating the feasibility of combining scBCR-seq and Ab-seq for reconstructing paired-chain Ig sequences from the serum antibody repertoire. Collectively, our work serves as a proof-of-principle for combining bulk sequencing, single-cell sequencing, and mass spectrometry as complementary methods towards capturing humoral immunity in its entirety.

Secondary bone marrow graft loss after third-party virus-specific T cell infusion: Case report of a rare complication

Virus-specific T cells (VST) from partially-HLA matched donors have been effective for treatment of refractory viral infections in immunocompromised patients in prior studies with a good safety profile, but rare adverse events have been described. Here we describe a unique and severe adverse event of VST therapy in an infant with severe combined immunodeficiency, who receives, as part of a clinical trial (NCT03475212), third party VSTs for treating cytomegalovirus viremia following bone marrow transplantation. At one-month post-VST infusion, rejection of graft and reversal of chimerism is observed, as is an expansion of T cells exclusively from the VST donor. Single-cell gene expression and T cell receptor profiling demonstrate a narrow repertoire of predominantly activated CD4+ T cells in the recipient at the time of rejection, with the repertoire overlapping more with that of peripheral blood from VST donor than the infused VST product. This case thus demonstrates a rare but serious side effect of VST therapy.

Spatiotemporal development of the human T follicular helper cell response to Influenza vaccination

We profiled blood and draining lymph node (LN) samples from human volunteers after influenza vaccination over two years to define evolution in the T follicular helper cell (TFH) response. We show LN TFH cells expanded in a clonal-manner during the first two weeks after vaccination and persisted within the LN for up to six months. LN and circulating TFH (cTFH) clonotypes overlapped but had distinct kinetics. LN TFH cell phenotypes were heterogeneous and mutable, first differentiating into pre-TFH during the month after vaccination before maturing into GC and IL-10+ TFH cells. TFH expansion, upregulation of glucose metabolism, and redifferentiation into GC TFH cells occurred with faster kinetics after re-vaccination in the second year. We identified several influenza-specific TFH clonal lineages, including multiple responses targeting internal influenza proteins, and show each TFH state is attainable within a lineage. This study demonstrates that human TFH cells form a durable and dynamic multi-tissue network.

Unique features of the TCR repertoire of reactivated memory T cells in the experimental mouse tumor model

T cell engineering with T cell receptors (TCR) specific to tumor antigens has become a breakthrough towards personalized cancer adoptive cell immunotherapy. However, the search for therapeutic TCRs is often challenging, and effective strategies are strongly required for the identification and enrichment of tumor-specific T cells that express TCRs with superior functional characteristics. Using an experimental mouse tumor model, we studied sequential changes in TCR repertoire features of T cells involved in the primary and secondary immune responses to allogeneic tumor antigens. In-depth bioinformatics analysis of TCR repertoires showed differences in reactivated memory T cells compared to primarily activated effectors. After cognate antigen re-encounter, memory cells were enriched with clonotypes that express α-chain TCR with high potential cross-reactivity and enhanced strength of interaction with both MHC and docked peptides. Our findings suggest that functionally true memory T cells could be a better source of therapeutic TCRs for adoptive cell therapy. No marked changes were observed in the physicochemical characteristics of TCRβ in reactivated memory clonotypes, indicative of the dominant role of TCRα in the secondary allogeneic immune response. The results of this study could further contribute to the development of TCR-modified T cell products based on the phenomenon of TCR chain centricity.

PTEN directs developmental and metabolic signaling for innate-like T cell fate and tissue homeostasis

Phosphatase and tensin homologue (PTEN) is frequently mutated in human cancer, but its roles in lymphopoiesis and tissue homeostasis remain poorly defined. Here we show that PTEN orchestrates a two-step developmental process linking antigen receptor and IL-23-Stat3 signalling to type-17 innate-like T cell generation. Loss of PTEN leads to pronounced accumulation of mature IL-17-producing innate-like T cells in the thymus. IL-23 is essential for their accumulation, and ablation of IL-23 or IL-17 signalling rectifies the reduced survival of female PTEN-haploinsufficient mice that model human patients with PTEN mutations. Single-cell transcriptome and network analyses revealed the dynamic regulation of PTEN, mTOR and metabolic activities that accompanied type-17 cell programming. Furthermore, deletion of mTORC1 or mTORC2 blocks PTEN loss-driven type-17 cell accumulation, and this is further shaped by the Foxo1 and Stat3 pathways. Collectively, our study establishes developmental and metabolic signalling networks underpinning type-17 cell fate decisions and their functional effects at coordinating PTEN-dependent tissue homeostasis.

TCR repertoire profiling revealed antigen-driven CD8+ T cell clonal groups shared in synovial fluid of patients with spondyloarthritis

Spondyloarthritis (SpA) comprises a number of inflammatory rheumatic diseases with overlapping clinical manifestations. Strong association with several HLA-I alleles and T cell infiltration into an inflamed joint suggest involvement of T cells in SpA pathogenesis. In this study, we performed high-throughput T cell repertoire profiling of synovial fluid (SF) and peripheral blood (PB) samples collected from a large cohort of SpA patients. We showed that synovial fluid is enriched with expanded T cell clones that are shared between patients with similar HLA genotypes and persist during recurrent synovitis. Using an algorithm for identification of TCRs involved in immune response we discovered several antigen-driven CD8+ clonal groups associated with risk HLA-B*27 or HLA-B*38 alleles. We further show that these clonal groups were enriched in SF and had higher frequency in PB of SpA patients vs healthy donors, implying their relevance to SpA pathogenesis. Several of the groups were shared among patients with different SpAs that suggests a common immunopathological mechanism of the diseases. In summary, our results provide evidence for the role of specific CD8+ T cell clones in pathogenesis of SpA.

Common Trajectories of Highly Effective CD19-Specific CAR T Cells Identified by Endogenous T-cell Receptor Lineages

Current chimeric antigen receptor-modified (CAR) T-cell products are evaluated in bulk, without assessing functional heterogeneity. We therefore generated a comprehensive single-cell gene expression and T-cell receptor (TCR) sequencing data set using pre- and postinfusion CD19-CAR T cells from blood and bone marrow samples of pediatric patients with B-cell acute lymphoblastic leukemia. We identified cytotoxic postinfusion cells with identical TCRs to a subset of preinfusion CAR T cells. These effector precursor cells exhibited a unique transcriptional profile compared with other preinfusion cells, corresponding to an unexpected surface phenotype (TIGIT+, CD62Llo, CD27-). Upon stimulation, these cells showed functional superiority and decreased expression of the exhaustion-associated transcription factor TOX. Collectively, these results demonstrate diverse effector potentials within preinfusion CAR T-cell products, which can be exploited for therapeutic applications. Furthermore, we provide an integrative experimental and analytic framework for elucidating the mechanisms underlying effector development in CAR T-cell products.

SIGNIFICANCE

Utilizing clonal trajectories to define transcriptional potential, we find a unique signature of CAR T-cell effector precursors present in preinfusion cell products. Functional assessment of cells with this signature indicated early effector potential and resistance to exhaustion, consistent with postinfusion cellular patterns observed in patients. This article is highlighted in the In This Issue feature, p. 2007.

Clonotypic IgH Response against Systemic Viral infection in Pronephros and Spleen of a Teleost Fish

Upon infection, B lymphocytes develop clonal responses. In teleost fish, which lack lymph nodes, the kinetics and location of B cell responses remain poorly characterized. Fish pronephros is the site of B cell differentiation and the main niche for persistence of plasma cells. In this study, we undertook the analysis of the rainbow trout IgHμ repertoire in this critical tissue for humoral adaptive immunity after primary immunization and boost with a rhabdovirus, the viral hemorrhagic septicemia virus (VHSV). We used a barcoded 5' RACE-cDNA sequencing approach to characterize modifications of the IgHμ repertoire, including VH usage in expressed V(D)J rearrangements, clonal diversity, and clonotype sharing between individual fish and treatments. In the pronephros, our approach quantified the clonotype frequency across the whole IgH repertoire (i.e., with all VH), measuring the frequency of Ag-responding clonotypes. Viral infection led to extensive modifications of the pronephros B cell repertoire, implicating several VH subgroups after primary infection. In contrast, only modest changes in repertoire persisted 5 mo later, including VHSV-specific public expansions. The IgM public response implicating IgHV1-18 and JH5, previously described in spleen, was confirmed in pronephros in all infected fish, strongly correlated to the response. However, the distribution of top clonotypes showed that pronephros and spleen B cells constitute distinct compartments with different IgH repertoires. Unexpectedly, after boost, the frequency of anti-VHSV clonotypes decreased both in pronephros and spleen, raising questions about B cell circulation. A better monitoring of B cell response kinetics in lymphoid tissues will be an essential step to understand B memory and plasmocyte formation mechanisms in fish.

Combining genotypes and T cell receptor distributions to infer genetic loci determining V(D)J recombination probabilities

Every T cell receptor (TCR) repertoire is shaped by a complex probabilistic tangle of genetically determined biases and immune exposures. T cells combine a random V(D)J recombination process with a selection process to generate highly diverse and functional TCRs. The extent to which an individual’s genetic background is associated with their resulting TCR repertoire diversity has yet to be fully explored. Using a previously published repertoire sequencing dataset paired with high-resolution genome-wide genotyping from a large human cohort, we infer specific genetic loci associated with V(D)J recombination probabilities using genome-wide association inference. We show that V(D)J gene usage profiles are associated with variation in the TCRB locus and, specifically for the functional TCR repertoire, variation in the major histocompatibility complex locus. Further, we identify specific variations in the genes encoding the Artemis protein and the TdT protein to be associated with biasing junctional nucleotide deletion and N-insertion, respectively. These results refine our understanding of genetically-determined TCR repertoire biases by confirming and extending previous studies on the genetic determinants of V(D)J gene usage and providing the first examples of trans genetic variants which are associated with modifying junctional diversity. Together, these insights lay the groundwork for further explorations into how immune responses vary between individuals.

SARS-CoV-2 mRNA vaccination elicits a robust and persistent T follicular helper cell response in humans

SARS-CoV-2 mRNA vaccines induce robust anti-spike (S) antibody and CD4+ T cell responses. It is not yet clear whether vaccine-induced follicular helper CD4+ T (TFH) cell responses contribute to this outstanding immunogenicity. Using fine-needle aspiration of draining axillary lymph nodes from individuals who received the BNT162b2 mRNA vaccine, we evaluated the T cell receptor sequences and phenotype of lymph node TFH. Mining of the responding TFH T cell receptor repertoire revealed a strikingly immunodominant HLA-DPB1∗04-restricted response to S167-180 in individuals with this allele, which is among the most common HLA alleles in humans. Paired blood and lymph node specimens show that while circulating S-specific TFH cells peak one week after the second immunization, S-specific TFH persist at nearly constant frequencies for at least six months. Collectively, our results underscore the key role that robust TFH cell responses play in establishing long-term immunity by this efficacious human vaccine.

Analysis of T-Cell Receptor Repertoire in Transplantation: Fingerprint of T Cell-mediated Alloresponse

T cells play a key role in determining allograft function by mediating allogeneic immune responses to cause rejection, and recent work pointed their role in mediating tolerance in transplantation. The unique T-cell receptor (TCR) expressed on the surface of each T cell determines the antigen specificity of the cell and can be the specific fingerprint for identifying and monitoring. Next-generation sequencing (NGS) techniques provide powerful tools for deep and high-throughput TCR profiling, and facilitate to depict the entire T cell repertoire profile and trace antigen-specific T cells in circulation and local tissues. Tailing T cell transcriptomes and TCR sequences at the single cell level provides a full landscape of alloreactive T-cell clones development and biofunction in alloresponse. Here, we review the recent advances in TCR sequencing techniques and computational tools, as well as the recent discovery in overall TCR profile and antigen-specific T cells tracking in transplantation. We further discuss the challenges and potential of using TCR sequencing-based assays to profile alloreactive TCR repertoire as the fingerprint for immune monitoring and prediction of rejection and tolerance.

ImmunoDataAnalyzer: a bioinformatics pipeline for processing barcoded and UMI tagged immunological NGS data

Background

Next-generation sequencing (NGS) is nowadays the most used high-throughput technology for DNA sequencing. Among others NGS enables the in-depth analysis of immune repertoires. Research in the field of T cell receptor (TCR) and immunoglobulin (IG) repertoires aids in understanding immunological diseases. A main objective is the analysis of the V(D)J recombination defining the structure and specificity of the immune repertoire. Accurate processing, evaluation and visualization of immune repertoire NGS data is important for better understanding immune responses and immunological behavior.

Results

ImmunoDataAnalyzer (IMDA) is a pipeline we have developed for automatizing the analysis of immunological NGS data. IMDA unites the functionality from carefully selected immune repertoire analysis software tools and covers the whole spectrum from initial quality control up to the comparison of multiple immune repertoires. It provides methods for automated pre-processing of barcoded and UMI tagged immune repertoire NGS data, facilitates the assembly of clonotypes and calculates key figures for describing the immune repertoire. These include commonly used clonality and diversity measures, as well as indicators for V(D)J gene segment usage and between sample similarity. IMDA reports all relevant information in a compact summary containing visualizations, calculations, and sample details, all of which serve for a more detailed overview. IMDA further generates an output file including key figures for all samples, designed to serve as input for machine learning frameworks to find models for differentiating between specific traits of samples.

Conclusions

IMDA constructs TCR and IG repertoire data from raw NGS reads and facilitates descriptive data analysis and comparison of immune repertoires. The IMDA workflow focus on quality control and ease of use for non-computer scientists. The provided output directly facilitates the interpretation of input data and includes information about clonality, diversity, clonotype overlap as well as similarity, and V(D)J gene segment usage. IMDA further supports the detection of sample swaps and cross-sample contamination that potentially occurred during sample preparation. In summary, IMDA reduces the effort usually required for immune repertoire data analysis by providing an automated workflow for processing raw NGS data into immune repertoires and subsequent analysis. The implementation is open-source and available on https://bioinformatics.fh-hagenberg.at/immunoanalyzer/.

Dual UMIs and Dual Barcodes With Minimal PCR Amplification Removes Artifacts and Acquires Accurate Antibody Repertoire

Antibody repertoire sequencing (Rep-seq) has been widely used to reveal repertoire dynamics and to interrogate antibodies of interest at single nucleotide-level resolution. However, polymerase chain reaction (PCR) amplification introduces extensive artifacts including chimeras and nucleotide errors, leading to false discovery of antibodies and incorrect assessment of somatic hypermutations (SHMs) which subsequently mislead downstream investigations. Here, a novel approach named DUMPArts, which improves the accuracy of antibody repertoires by labeling each sample with dual barcodes and each molecule with dual unique molecular identifiers (UMIs) via minimal PCR amplification to remove artifacts, is developed. Tested by ultra-deep Rep-seq data, DUMPArts removed inter-sample chimeras, which cause artifactual shared clones and constitute approximately 15% of reads in the library, as well as intra-sample chimeras with erroneous SHMs and constituting approximately 20% of the reads, and corrected base errors and amplification biases by consensus building. The removal of these artifacts will provide an accurate assessment of antibody repertoires and benefit related studies, especially mAb discovery and antibody-guided vaccine design.

Tools for fundamental analysis functions of TCR repertoires: a systematic comparison

The full set of T cell receptors (TCRs) in an individual is known as his or her TCR repertoire. Defining TCR repertoires under physiological conditions and in response to a disease or vaccine may lead to a better understanding of adaptive immunity and thus has great biological and clinical value. In the past decade, several high-throughput sequencing-based tools have been developed to assign TCRs to germline genes and to extract complementarity-determining region 3 (CDR3) sequences using different algorithms. Although these tools claim to be able to perform the full range of fundamental TCR repertoire analyses, there is no clear consensus of which tool is best suited to particular projects. Here, we present a systematic analysis of 12 available TCR repertoire analysis tools using simulated data, with an emphasis on fundamental analysis functions. Our results shed light on the detailed functions of TCR repertoire analysis tools and may therefore help researchers in the field to choose the right tools for their particular experimental design.

Deciphering Repertoire of B16 Melanoma Reactive TCRs by Immunization, In Vitro Restimulation and Sequencing of IFNγ-Secreting T Cells

Recent advances in cancer immunotherapy have great promise for the treatment of solid tumors. One of the key limiting factors that hamper the decoding of physiological responses to these therapies is the inability to distinguish between specific and nonspecific responses. The identification of tumor-specific lymphocytes is also the most challenging step in cancer cell therapies such as adoptive cell transfer and T cell receptor (TCR) cloning. Here, we have elaborated a protocol for the identification of tumor-specific T lymphocytes and the deciphering of their repertoires. B16 melanoma engraftment following anti-PD1 checkpoint therapy provides better antitumor immunity compared to repetitive immunization with heat-shocked tumor cells. We have also revealed that the most error-prone part of dendritic cell (DC) generation, i.e., their maturation step, can be omitted if DCs are cultured at a sufficiently high density. Using this optimized protocol, we have achieved a robust IFNγ response to B16F0 antigens, but only within CD4+ T helper cells. A comparison of the repertoires of IFNγ-positive and -negative cells shows a prominent enrichment of certain clones with putative tumor specificity among the IFNγ+ fraction. In summary, our optimized protocol and the data provided here will aid in the acquisition of broad statistical data and the creation of a meaningful database of B16-specific TCRs.

TCRpair: prediction of functional pairing between HLA-A*02:01-restricted T cell receptor α and β chains

SUMMARY

The ability of a T cell to recognize foreign peptides is defined by a single α and a single β hypervariable complementarity determining region (CDR3), which together form the T cell receptor (TCR) heterodimer. In ∼30%-35% of T cells, two α chains are expressed at the mRNA level but only one α chain is part of the functional TCR. This effect can also be observed for β chains, although it is less common. The identification of functional α/β chain pairs is instrumental in high-throughput characterization of therapeutic TCRs. TCRpair is the first method that predicts whether an α and β chain pair forms a functional, HLA-A*02:01 specific TCR without requiring the sequence of a recognized peptide. By taking additional amino acids flanking the CDR3 regions into account, TCRpair achieves an AUC of 0.71.

AVAILABILITY

TCRpair is implemented in Python using TensorFlow 2.0 and is freely available at https://www.github.com/amoesch/TCRpair.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A comprehensive model based on temporal dynamics of peripheral T cell repertoire for predicting post-treatment distant metastasis of nasopharyngeal carcinoma

Many nasopharyngeal carcinoma (NPC) patients develop distant metastases after treatment, leading to poor outcomes. To date, there are no peripheral biomarkers suitable for all NPC patients to predict distant metastasis. Hence, we purposed to develop a noninvasive comprehensive model for predicting post-treatment distant metastasis of all NPC. Since T-cell receptor β chain (TCRB) repertoire has achieved prognostic prediction in many cancers, the clinical characteristics and parameters of TCRB repertoire of 71 cases of peripheral blood samples (pairwise pre-treatment and post-treatment samples from 40 NPC patients who without (nM, n = 21) or with (M, n = 19) post-treatment distant metastasis) were collected. The least absolute shrinkage and selection operator algorithm was used to construct a distant metastasis prediction model. In terms of TCRB repertoire parameters, the diversity of TCRB repertoire was significantly decreased in M group after treatment but not in nM group. Ascending TCRB diversity and higher similarity between pre- and post-treatment samples showed better distant metastasis-free survival (DMFS). The similarity still had robust DMFS prediction in patients with reduced TCRB diversity. More importantly, the 5-factor comprehensive model consisting of basic clinical characteristics and TCRB repertoire indices showed a higher prognostic accuracy than any one individual factor in DMFS predicting. In conclusion, treatment had different effects on the composition of TCRB repertoire in patients without and with post-treatment distant metastasis. The dynamics of TCRB diversity, the similarity of TCRB repertoires, and combinations of these factors with basic clinical characteristics could serve as noninvasive DMFS predictors for all NPC patients.

Foxp3 enhancers synergize to maximize regulatory T cell suppressive capacity

T reg cells bearing a diverse antigen receptor repertoire suppress pathogenic T cells and maintain immune homeostasis during their long lifespan. How their robust function is determined genetically remains elusive. Here, we investigate the regulatory space of the cis-regulatory elements of T reg lineage-specifying factor Foxp3. Foxp3 enhancers are known as distinct readers of environmental cues controlling T reg cell induction or lineage stability. However, their single deficiencies cause mild, if any, immune dysregulation, leaving the key transcriptional mechanisms determining Foxp3 expression and thereby T reg cell suppressive capacity uncertain. We examined the collective activities of Foxp3 enhancers and found that they coordinate to maximize T reg cell induction, Foxp3 expression level, or lineage stability through distinct modes and that ablation of synergistic enhancers leads to lethal autoimmunity in young mice. Thus, the induction and maintenance of a diverse, stable T reg cell repertoire rely on combinatorial Foxp3 enhancers, suggesting broad, stage-specific, synergistic activities of cell-intrinsic factors and cell-extrinsic cues in determining T reg cell suppressive capacity.

High-throughput and single-cell T cell receptor sequencing technologies

T cells express T cell receptors (TCRs) composed of somatically recombined TCRα and TCRβ chains, which mediate recognition of major histocompatibility complex (MHC)-antigen complexes and drive the antigen-specific adaptive immune response to pathogens and cancer. The TCR repertoire in each individual is highly diverse, which allows for recognition of a wide array of foreign antigens, but also presents a challenge in analyzing this response using conventional methods. Recent studies have developed high-throughput sequencing technologies to identify TCR sequences, analyze their antigen specificities using experimental and computational tools, and pair TCRs with transcriptional and epigenetic cell state phenotypes in single cells. In this Review, we highlight these technological advances and describe how they have been applied to discover fundamental insights into T cell-mediated immunity.

UMIc: A Preprocessing Method for UMI Deduplication and Reads Correction

A recent refinement in high-throughput sequencing involves the incorporation of unique molecular identifiers (UMIs), which are random oligonucleotide barcodes, on the library preparation steps. A UMI adds a unique identity to different DNA/RNA input molecules through polymerase chain reaction (PCR) amplification, thus reducing bias of this step. Here, we propose an alignment free framework serving as a preprocessing step of fastq files, called UMIc, for deduplication and correction of reads building consensus sequences from each UMI. Our approach takes into account the frequency and the Phred quality of nucleotides and the distances between the UMIs and the actual sequences. We have tested the tool using different scenarios of UMI-tagged library data, having in mind the aspect of a wide application. UMIc is an open-source tool implemented in R and is freely available from https://github.com/BiodataAnalysisGroup/UMIc.

Generation of TCRα-transduced T cells for adoptive transfer therapy of salmonellosis in mice

Adoptive transfer therapy has great potential to treat diseases such as cancer as well as autoimmune and infectious diseases. Identification of chain-centric T cell receptors (TCRs) with the dominant-active antigen-specific α-chains (TCRα) can significantly improve the efficacy of adoptive cell therapy while reducing time, labor, and costs of generation of TCR-modified antigen-specific T cells. This protocol describes how to generate salmonella-specific TCRα-modified mouse T cells by retroviral transduction and evaluate their functional activity in vivo in the mouse model of salmonellosis.

For complete details on the use and execution of this protocol, please refer to Kalinina et al. (2020).

•

Adoptive T cell therapy has great potential for treatment of infections

•

Protocol describes generating T cells transduced with dominant pathogen-specific TCRα

•

Protocol describes usage of specific TCRα-transduced T cells in mice salmonellosis

Distinct organization of adaptive immunity in the long-lived rodent Spalax galili

A balanced immune response is a cornerstone of healthy aging. Here, we uncover distinctive features of the long-lived blind mole-rat (Spalax spp.) adaptive immune system, relative to humans and mice. The T-cell repertoire remains diverse throughout the Spalax lifespan, suggesting a paucity of large long-lived clones of effector-memory T cells. Expression of master transcription factors of T-cell differentiation, as well as checkpoint and cytotoxicity genes, remains low as Spalax ages. The thymus shrinks as in mice and humans, while interleukin-7 and interleukin-7 receptor expression remains high, potentially reflecting the sustained homeostasis of naive T cells. With aging, immunoglobulin hypermutation level does not increase and the immunoglobulin-M repertoire remains diverse, suggesting shorter B-cell memory and sustained homeostasis of innate-like B cells. The Spalax adaptive immune system thus appears biased towards sustained functional and receptor diversity over specialized, long-lived effector-memory clones-a unique organizational strategy that potentially underlies this animal's extraordinary longevity and healthy aging.

Functionally specialized human CD4+ T-cell subsets express physicochemically distinct TCRs

The organizational integrity of the adaptive immune system is determined by functionally discrete subsets of CD4⁺ T cells, but it has remained unclear to what extent lineage choice is influenced by clonotypically expressed T-cell receptors (TCRs). To address this issue, we used a high-throughput approach to profile the αβ TCR repertoires of human naive and effector/memory CD4⁺ T-cell subsets, irrespective of antigen specificity. Highly conserved physicochemical and recombinatorial features were encoded on a subset-specific basis in the effector/memory compartment. Clonal tracking further identified forbidden and permitted transition pathways, mapping effector/memory subsets related by interconversion or ontogeny. Public sequences were largely confined to particular effector/memory subsets, including regulatory T cells (Tregs), which also displayed hardwired repertoire features in the naive compartment. Accordingly, these cumulative repertoire portraits establish a link between clonotype fate decisions in the complex world of CD4⁺ T cells and the intrinsic properties of somatically rearranged TCRs.

MAUI-seq: Metabarcoding using amplicons with unique molecular identifiers to improve error correction

Sequencing and PCR errors are a major challenge when characterizing genetic diversity using high-throughput amplicon sequencing (HTAS). We have developed a multiplexed HTAS method, MAUI-seq, which uses unique molecular identifiers (UMIs) to improve error correction by exploiting variation among sequences associated with a single UMI. Erroneous sequences are recognized because, across the data set, they are over-represented among the minor sequences associated with UMIs. We show that two main advantages of this approach are efficient elimination of chimeric and other erroneous reads, outperforming dada2 and unoise3, and the ability to confidently recognize genuine alleles that are present at low abundance or resemble chimeras. The method provides sensitive and flexible profiling of diversity and is readily adaptable to most HTAS applications, including microbial 16S rRNA profiling and metabarcoding of environmental DNA.

T-cell tracking, safety, and effect of low-dose donor memory T-cell infusions after αβ T cell-depleted hematopoietic stem cell transplantation

The delayed recovery of adaptive immunity underlies transplant-related mortality (TRM) after αβ T cell-depleted hematopoietic stem cell transplantation (HSCT). We tested the use of low-dose memory donor lymphocyte infusions (mDLIs) after engraftment of αβ T cell-depleted grafts.A cohort of 131 pediatric patients (median age 9 years) were grafted with αβ T cell-depleted products from either haplo (n = 79) or unrelated donors (n = 52). After engraftment, patients received mDLIs prepared by CD45RA depletion. Cell dose was escalated monthly from 25 × 10³ to 100 × 10³/kg (haplo) and from 100 × 10³ to 300 × 10³ /kg (MUD). In a subcohort of 16 patients, T-cell receptor (TCR) repertoire profiling with deep sequencing was used to track T-cell clones and to evaluate the contribution of mDLI to the immune repertoire.In total, 343 mDLIs were administered. The cumulative incidence (CI) of grades II and III de novo acute graft-versus-host disease (aGVHD) was 5% and 2%, respectively, and the CI of chronic graft-versus-host disease was 7%. Half of the patients with undetectable CMV-specific T cells before mDLI recovered CMV-specific T cells. TCR repertoire profiling confirmed that mDLI-derived T cells significantly contribute to the TCR repertoire up to 1 year after HSCT and include persistent, CMV-specific T-cell clones.

TCRs with segment TRAV9-2 or a CDR3 histidine are overrepresented among nickel-specific CD4+ T cells

BACKGROUND

Nickel is the most frequent cause of T cell-mediated allergic contact dermatitis worldwide. In vitro, CD4+ T cells from all donors respond to nickel but the involved αβ T cell receptor (TCR) repertoire has not been comprehensively analyzed.

METHODS

We introduce CD154 (CD40L) upregulation as a fast, unbiased, and quantitative method to detect nickel-specific CD4+ T cells ex vivo in blood of clinically characterized allergic and non allergic donors. Naïve (CCR7+ CD45RA+) and memory (not naïve) CD154+ CD4+ T cells were analyzed by flow cytometry after 5 hours of stimulation with 200 µmol/L NiSO₄ ., TCR α- and β-chains of sorted nickel-specific and control cells were studied by high-throughput sequencing.

RESULTS

Stimulation of PBMCs with NiSO₄ induced CD154 expression on ~0.1% (mean) of naïve and memory CD4+ T cells. In allergic donors with recent positive patch test, memory frequencies further increased ~13-fold and were associated with markers of in vivo activation. CD154 expression was TCR-mediated since single clones could be specifically restimulated. Among nickel-specific CD4+ T cells of allergic and non allergic donors, TCRs expressing the α-chain segment TRAV9-2 or a histidine in their α- or β-chain complementarity determining region 3 (CDR3) were highly overrepresented.

CONCLUSIONS

Induced CD154 expression represents a reliable method to study nickel-specific CD4+ T cells. TCRs with particular features respond in all donors, while strongly increased blood frequencies indicate nickel allergy for some donors. Our approach may be extended to other contact allergens for the further development of diagnostic and predictive in vitro tests.

Immune Literacy: Reading, Writing, and Editing Adaptive Immunity

Advances in reading, writing, and editing DNA are providing unprecedented insights into the complexity of immunological systems. This combination of systems and synthetic biology methods is enabling the quantitative and precise understanding of molecular recognition in adaptive immunity, thus providing a framework for reprogramming immune responses for translational medicine. In this review, we will highlight state-of-the-art methods such as immune repertoire sequencing, immunoinformatics, and immunogenomic engineering and their application toward adaptive immunity. We showcase novel and interdisciplinary approaches that have the promise of transforming the design and breadth of molecular and cellular immunotherapies.

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.

Measuring Intratumoral Heterogeneity of Immune Repertoires

There is considerable clinical and fundamental value in measuring the clonal heterogeneity of T and B cell expansions in tumors and tumor-associated lymphoid structures—along with the associated heterogeneity of the tumor neoantigen landscape—but such analyses remain challenging to perform. Here, we propose a straightforward approach to analyze the heterogeneity of immune repertoires between different tissue sections in a quantitative and controlled way, based on a beta-binomial noise model trained on control replicates obtained at the level of single-cell suspensions. This approach allows to identify local clonal expansions with high accuracy. We reveal in situ proliferation of clonal T cells in a mouse model of melanoma, and analyze heterogeneity of immunoglobulin repertoires between sections of a metastatically-infiltrated lymph node in human melanoma and primary human colon tumor. On the latter example, we demonstrate the importance of training the noise model on datasets with depth and content that is comparable to the samples being studied. Altogether, we describe here the crucial basic instrumentarium needed to facilitate proper experimental setup planning in the rapidly evolving field of intratumoral immune repertoires, from the wet lab to bioinformatics analysis.

TCR Convergence in Individuals Treated With Immune Checkpoint Inhibition for Cancer

Tumor antigen-driven selection may expand T cells having T cell receptors (TCRs) of shared antigen specificity but different amino acid or nucleotide sequence in a process known as TCR convergence. Substitution sequencing errors introduced by TCRβ (TCRB) repertoire sequencing may create artifacts resembling TCR convergence. Given the anticipated differences in substitution error rates across different next-generation sequencing platforms, the choice of platform could be consequential. To test this, we performed TCRB sequencing on the same peripheral blood mononuclear cells (PBMC) from individuals with cancer receiving anti-CTLA-4 or anti-PD-1 using an Illumina-based approach (Sequenta) and an Ion Torrent-based approach (Oncomine TCRB-LR). While both approaches found similar TCR diversity, clonality, and clonal overlap, we found that Illumina-based sequencing resulted in higher TCR convergence than with the Ion Torrent approach. To build upon this initial observation we conducted a systematic comparison of Illumina-based TCRB sequencing assays, including those employing molecular barcodes, with the Oncomine assay, revealing differences in the frequency of convergent events, purportedly artifactual rearrangements, and sensitivity of detection. Finally, we applied the Ion Torrent-based approach to evaluate clonality and convergence in a cohort of individuals receiving anti-CTLA-4 blockade for cancer. We found that clonality and convergence independently predicted response and could be combined to improve the accuracy of a logistic regression classifier. These results demonstrate the importance of the sequencing platform in assessing TCRB convergence.

Deep sequencing of B cell receptor repertoire

Immune repertoire is a collection of enormously diverse adaptive immune cells within an individual. As the repertoire shapes and represents immunological conditions, identification of clones and characterization of diversity are critical for understanding how to protect ourselves against various illness such as infectious diseases and cancers. Over the past several years, fast growing technologies for high throughput sequencing have facilitated rapid advancement of repertoire research, enabling us to observe the diversity of repertoire at an unprecedented level. Here, we focus on B cell receptor (BCR) repertoire and review approaches to B cell isolation and sequencing library construction. These experiments should be carefully designed according to BCR regions to be interrogated, such as heavy chain full length, complementarity determining regions, and isotypes. We also highlight preprocessing steps to remove sequencing and PCR errors with unique molecular index and bioinformatics techniques. Due to the nature of massive sequence variation in BCR, caution is warranted when interpreting repertoire diversity from error-prone sequencing data. Furthermore, we provide a summary of statistical frameworks and bioinformatics tools for clonal evolution and diversity. Finally, we discuss limitations of current BCR-seq technologies and future perspectives on advances in repertoire sequencing.

Nasal Vaccination Drives Modifications of Nasal and Systemic Antibody Repertoires in Rainbow Trout

Bony fish represent the most basal vertebrate branch with a dedicated mucosal immune system, which comprises immunologically heterogeneous microenvironments armed with innate and adaptive components. In rainbow trout (Oncorhynchus mykiss), a nasopharynx-associated lymphoid tissue (NALT) was recently described as a diffuse network of myeloid and lymphoid cells located in the olfactory organ of fish. Several studies have demonstrated high levels of protection conferred by nasal vaccines against viral and bacterial pathogens; however, the mechanisms underlying the observed protection are not well understood. We applied 5'RACE and a deep sequencing-based approach to investigate the clonal structure of the systemic and mucosal rainbow trout B cell repertoire. The analysis of Ig repertoire in control trout suggests different structures of IgM and IgT spleen and NALT repertoires, with restricted repertoire diversity in NALT. Nasal and injection vaccination with a bacterial vaccine revealed unique dynamics of IgM and IgT repertoires at systemic and mucosal sites and the remarkable ability of nasal vaccines to induce spleen Ig responses. Our findings provide an important immunological basis for the effectiveness of nasal vaccination in fish and other vertebrate animals and will help the design of future nasal vaccination strategies.

RNase H-dependent PCR-enabled T-cell receptor sequencing for highly specific and efficient targeted sequencing of T-cell receptor mRNA for single-cell and repertoire analysis

RNase H-dependent PCR-enabled T-cell receptor sequencing (rhTCRseq) can be used to determine paired alpha/beta T-cell receptor (TCR) clonotypes in single cells or perform alpha and beta TCR repertoire analysis in bulk RNA samples. With the enhanced specificity of RNase H-dependent PCR (rhPCR), it achieves TCR-specific amplification and addition of dual-index barcodes in a single PCR step. For single cells, the protocol includes sorting of single cells into plates, generation of cDNA libraries, a TCR-specific amplification step, a second PCR on pooled sample to generate a sequencing library, and sequencing. In the bulk method, sorting and cDNA library steps are replaced with a reverse-transcriptase (RT) reaction that adds a unique molecular identifier (UMI) to each cDNA molecule to improve the accuracy of repertoire-frequency measurements. Compared to other methods for TCR sequencing, rhTCRseq has a streamlined workflow and the ability to analyze single cells in 384-well plates. Compared to TCR reconstruction from single-cell transcriptome sequencing data, it improves the success rate for obtaining paired alpha/beta information and ensures recovery of complete complementarity-determining region 3 (CDR3) sequences, a prerequisite for cloning/expression of discovered TCRs. Although it has lower throughput than droplet-based methods, rhTCRseq is well-suited to analysis of small sorted populations, especially when analysis of 96 or 384 single cells is sufficient to identify predominant T-cell clones. For single cells, sorting typically requires 2-4 h and can be performed days, or even months, before library construction and data processing, which takes ~4 d; the bulk RNA protocol takes ~3 d.

The Pipeline Repertoire for Ig-Seq Analysis

With the advent of high-throughput sequencing of immunoglobulin genes (Ig-Seq), the understanding of antibody repertoires and their dynamics among individuals and populations has become an exciting area of research. There is an increasing number of computational tools that aid in every step of the immune repertoire characterization. However, since not all tools function identically, every pipeline has its unique rationale and capabilities, creating a rich blend of useful features that may appear intimidating for newcomer laboratories with the desire to plunge into immune repertoire analysis to expand and improve their research; hence, all pipeline strengths and differences may not seem evident. In this review we provide a practical and organized list of the current set of computational tools, focusing on their most attractive features and differences in order to carry out the characterization of antibody repertoires so that the reader better decides a strategic approach for the experimental design, and computational pathways for the analyses of immune repertoires.

Sequencing the peripheral blood B and T cell repertoire - Quantifying robustness and limitations

The adaptive immune response generates a large repertoire of T cells with T-cell receptors (TCR_alpha and TCR_beta) and B cells with immunoglobulins (Ig). The repertoire changes in response to antigen stimulation both through amplification of specific cells (clonal expansion) as well as somatic hypermutation of immunoglobulins. Alterations of the immune repertoire have been observed in response to acute disease, such as external pathogens, or chronic diseases, such as autoimmunity and cancer. Here we establish experimental and analytical protocols for quantifying the peripheral blood of healthy human individuals by profiling the immune repertoire for the Complementarity determining region 3 (CDR3) of the variable regions of TCR_beta (CDRβ3) and the IgG heavy chain (CDRH1, CDRH2, CDRH3). The results demonstrate that 40 ml of blood are sufficient to reliably capture the 10,000 most common TCR_beta and 1000 most common IgG and determine their relative frequency in the circulation. We conclude that by using an accessible sample size of human PBMC one is able to robustly monitor alterations in the immune repertoire.

Rare Event Detection Using Error-corrected DNA and RNA Sequencing

Conventional next-generation sequencing techniques (NGS) have allowed for immense genomic characterization for over a decade. Specifically, NGS has been used to analyze the spectrum of clonal mutations in malignancy. Though far more efficient than traditional Sanger methods, NGS struggles with identifying rare clonal and subclonal mutations due to its high error rate of ~0.5-2.0%. Thus, standard NGS has a limit of detection for mutations that are >0.02 variant allele fraction (VAF). While the clinical significance for mutations this rare in patients without known disease remains unclear, patients treated for leukemia have significantly improved outcomes when residual disease is <0.0001 by flow cytometry. In order to mitigate this artefactual background of NGS, numerous methods have been developed. Here we describe a method for Error-corrected DNA and RNA Sequencing (ECS), which involves tagging individual molecules with both a 16 bp random index for error-correction and an 8 bp patient-specific index for multiplexing. Our method can detect and track clonal mutations at variant allele fractions (VAFs) two orders of magnitude lower than the detection limit of NGS and as rare as 0.0001 VAF.

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

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

High-throughput sequencing of the T-cell receptor repertoire: pitfalls and opportunities

T-cell specificity is determined by the T-cell receptor, a heterodimeric protein coded for by an extremely diverse set of genes produced by imprecise somatic gene recombination. Massively parallel high-throughput sequencing allows millions of different T-cell receptor genes to be characterized from a single sample of blood or tissue. However, the extraordinary heterogeneity of the immune repertoire poses significant challenges for subsequent analysis of the data. We outline the major steps in processing of repertoire data, considering low-level processing of raw sequence files and high-level algorithms, which seek to extract biological or pathological information. The latest generation of bioinformatics tools allows millions of DNA sequences to be accurately and rapidly assigned to their respective variable V and J gene segments, and to reconstruct an almost error-free representation of the non-templated additions and deletions that occur. High-level processing can measure the diversity of the repertoire in different samples, quantify V and J usage and identify private and public T-cell receptors. Finally, we discuss the major challenge of linking T-cell receptor sequence to function, and specifically to antigen recognition. Sophisticated machine learning algorithms are being developed that can combine the paradoxical degeneracy and cross-reactivity of individual T-cell receptors with the specificity of the overall T-cell immune response. Computational analysis will provide the key to unlock the potential of the T-cell receptor repertoire to give insight into the fundamental biology of the adaptive immune system and to provide powerful biomarkers of disease.

High-throughput sequencing of murine immunoglobulin heavy chain repertoires using single side unique molecular identifiers on an Ion Torrent PGM

With the advent of high-throughput sequencing (HTS), profiling immunoglobulin (IG) repertoires has become an essential part of immunological research. Advances in sequencing technology enable the IonTorrent Personal Genome Machine (PGM) to cover the full-length of IG mRNA transcripts. Nucleotide insertions and deletions (indels) are the dominant errors of the PGM sequencing platform and can critically influence IG repertoire assessments. Here, we present a PGM-tailored IG repertoire sequencing approach combining error correction through unique molecular identifier (UID) barcoding and indel detection through ImMunoGeneTics (IMGT), the most commonly used sequence alignment database for IG sequences. Using artificially falsified sequences for benchmarking, we found that IMGT's underlying algorithms efficiently detect 98% of the introduced indels. Undetected indels are either located at the end of the sequences or produce masked frameshifts with an insertion and deletion in close proximity. The complementary determining regions 3 (CDR3s) are returned correct for up to 3 insertions or 3 deletions through conservative culling. We further show, that our PGM-tailored unique molecular identifiers result in highly accurate HTS data if combined with the presented processing strategy. In this regard, considering sequences with at least two copies from datasets with UID families of minimum 3 reads result in correct sequences with over 99% confidence. Finally, we show that the protocol can readily be used to generate homogenous datasets for bulk sequencing of murine bone marrow samples. Taken together, this approach will help to establish benchtop-scale sequencing of IG heavy chain transcripts in the field of IG repertoire research.

RepSeq Data Representativeness and Robustness Assessment by Shannon Entropy

High-throughput sequencing (HTS) has the potential to decipher the diversity of T cell repertoires and their dynamics during immune responses. Applied to T cell subsets such as T effector and T regulatory cells, it should help identify novel biomarkers of diseases. However, given the extreme diversity of TCR repertoires, understanding how the sequencing conditions, including cell numbers, biological and technical sampling and sequencing depth, impact the experimental outcome is critical to proper use of these data. Here, we assessed the representativeness and robustness of TCR repertoire diversity assessment according to experimental conditions. By comparative analyses of experimental datasets and computer simulations, we found that (i) for small samples, the number of clonotypes recovered is often higher than the number of cells per sample, even after removing the singletons; (ii) high-sequencing depth for small samples alters the clonotype distributions, which can be corrected by filtering the datasets using Shannon entropy as a threshold; and (iii) a single sequencing run at high depth does not ensure a good coverage of the clonotype richness in highly polyclonal populations, which can be better covered using multiple sequencing. Altogether, our results warrant better understanding and awareness of the limitation of TCR diversity analyses by HTS and justify the development of novel computational tools for improved modeling of the highly complex nature of TCR repertoires.

Analyzing Immunoglobulin Repertoires

Somatic assembly of T cell receptor and B cell receptor (BCR) genes produces a vast diversity of lymphocyte antigen recognition capacity. The advent of efficient high-throughput sequencing of lymphocyte antigen receptor genes has recently generated unprecedented opportunities for exploration of adaptive immune responses. With these opportunities have come significant challenges in understanding the analysis techniques that most accurately reflect underlying biological phenomena. In this regard, sample preparation and sequence analysis techniques, which have largely been borrowed and adapted from other fields, continue to evolve. Here, we review current methods and challenges of library preparation, sequencing and statistical analysis of lymphocyte receptor repertoire studies. We discuss the general steps in the process of immune repertoire generation including sample preparation, platforms available for sequencing, processing of sequencing data, measurable features of the immune repertoire, and the statistical tools that can be used for analysis and interpretation of the data. Because BCR analysis harbors additional complexities, such as immunoglobulin (Ig) (i.e., antibody) gene somatic hypermutation and class switch recombination, the emphasis of this review is on Ig/BCR sequence analysis.

History, applications, and challenges of immune repertoire research

The diversity of T and B cells in terms of their receptor sequences is huge in the vertebrate's immune system and provides broad protection against the vast diversity of pathogens. Immune repertoire is defined as the sum of T cell receptors and B cell receptors (also named immunoglobulin) that makes the organism's adaptive immune system. Before the emergence of high-throughput sequencing, the studies on immune repertoire were limited by the underdeveloped methodologies, since it was impossible to capture the whole picture by the low-throughput tools. The massive paralleled sequencing technology suits perfectly the researches on immune repertoire. In this article, we review the history of immune repertoire studies, in terms of technologies and research applications. Particularly, we discuss several aspects of challenges in this field and highlight the efforts to develop potential solutions, in the era of high-throughput sequencing of the immune repertoire.

Immune Repertoire Sequencing Using Molecular Identifiers Enables Accurate Clonality Discovery and Clone Size Quantification

Unique molecular identifiers (MIDs) have been demonstrated to effectively improve immune repertoire sequencing (IR-seq) accuracy, especially to identify somatic hypermutations in antibody repertoire sequencing. However, evaluating the sensitivity to detect rare T cells and the degree of clonal expansion in IR-seq has been difficult due to the lack of knowledge of T cell receptor (TCR) RNA molecule copy number and a generalized approach to estimate T cell clone size from TCR RNA molecule quantification. This limited the application of TCR repertoire sequencing (TCR-seq) in clinical settings, such as detecting minimal residual disease in lymphoid malignancies after treatment, evaluating effectiveness of vaccination and assessing degree of infection. Here, we describe using an MID Clustering-based IR-Seq (MIDCIRS) method to quantitatively study TCR RNA molecule copy number and clonality in T cells. First, we demonstrated the necessity of performing MID sub-clustering to eliminate erroneous sequences. Further, we showed that MIDCIRS enables a sensitive detection of a single cell in as many as one million naïve T cells and an accurate estimation of the degree of T cell clonal expression. The demonstrated accuracy, sensitivity, and wide dynamic range of MIDCIRS TCR-seq provide foundations for future applications in both basic research and clinical settings.

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.

Novel Method for High-Throughput Full-Length IGHV-D-J Sequencing of the Immune Repertoire from Bulk B-Cells with Single-Cell Resolution

Efficient and accurate high-throughput DNA sequencing of the adaptive immune receptor repertoire (AIRR) is necessary to study immune diversity in healthy subjects and disease-related conditions. The high complexity and diversity of the AIRR coupled with the limited amount of starting material, which can compromise identification of the full biological diversity makes such sequencing particularly challenging. AIRR sequencing protocols often fail to fully capture the sampled AIRR diversity, especially for samples containing restricted numbers of B lymphocytes. Here, we describe a library preparation method for immunoglobulin sequencing that results in an exhaustive full-length repertoire where virtually every sampled B-cell is sequenced. This maximizes the likelihood of identifying and quantifying the entire IGHV-D-J repertoire of a sample, including the detection of rearrangements present in only one cell in the starting population. The methodology establishes the importance of circumventing genetic material dilution in the preamplification phases and incorporates the use of certain described concepts: (1) balancing the starting material amount and depth of sequencing, (2) avoiding IGHV gene-specific amplification, and (3) using Unique Molecular Identifier. Together, this methodology is highly efficient, in particular for detecting rare rearrangements in the sampled population and when only a limited amount of starting material is available.

Overview of methodologies for T-cell receptor repertoire analysis

Background

The T-cell receptor (TCR), located on the surface of T cells, is responsible for the recognition of the antigen-major histocompatibility complex, leading to the initiation of an inflammatory response. Analysing the TCR repertoire may help to gain a better understanding of the immune system features and of the aetiology and progression of diseases, in particular those with unknown antigenic triggers. The extreme diversity of the TCR repertoire represents a major analytical challenge; this has led to the development of specialized methods which aim to characterize the TCR repertoire in-depth. Currently, next generation sequencing based technologies are most widely employed for the high-throughput analysis of the immune cell repertoire.

Results

Here, we report on the latest methodological advancements in the field by describing and comparing the available tools; from the choice of the starting material and library preparation method, to the sequencing technologies and data analysis.

Finally, we provide a practical example and our own experience by reporting some exemplary results from a small internal benchmark study, where current approaches from the literature and the market are employed and compared.

Conclusions

Several valid methods for clonotype identification and TCR repertoire analysis exist, however, a gold standard method for the field has not yet been identified. Depending on the purpose of the scientific study, some approaches may be more suitable than others. Finally, due to possible method specific biases, scientists must be careful when comparing results obtained using different methods.

Electronic supplementary material

The online version of this article (doi:10.1186/s12896-017-0379-9) contains supplementary material, which is available to authorized users.

Application of nonsense-mediated primer exclusion (NOPE) for preparation of unique molecular barcoded libraries

Background

Recently we proposed efficient method to exclude undesirable primers at any stage of amplification reaction, here termed NOPE (NOnsense-mediated Primer Exclusion). According to this method, added oligonucleotide overlapping with the 3′-end of unwanted amplification primer (NOPE oligo) simultaneously provides a template for its elongation. This elongation disrupts specificity of unwanted primer, preventing its further participation in PCR. The suggested approach allows to rationally manage the course of PCR reactions in order to facilitate analysis of complex DNA mixtures as well as to perform multistage PCR bypassing intermediate purification steps.

Results

Here we apply NOPE method to DNA library preparation for the high-throughput sequencing (HTS) with the PCR-based introduction of unique molecular identifiers (UMI). We show that NOPE oligo efficiently neutralizes UMI-containing oligonucleotides after introduction of UMI into sample DNA molecules, thus allowing to proceed with further amplification steps without purification and associated loss of starting material. At the same time, NOPE oligo does not affect the efficiency of target PCR amplification.

Conclusion

We describe a simple, robust and cheap modification of UMI-labeled HTS libraries preparation procedure, that allows to bypass purification step and thus to preserve starting material which may be limited, e.g. circulating tumor DNA, circulating fetal DNA, or small amounts of isolated cells of interest. Furthermore, demonstrated simplicity and robustness of NOPE method should make it popular in various PCR protocols.