Individual heritable differences result in unique cell lymphocyte receptor repertoires of naïve and antigen-experienced cells

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.

Quantitative trait loci mapping provides insights into the genetic regulation of dendritic cell numbers in mouse tissues

To explore the influence of genetics on homeostatic regulation of dendritic cell (DC) numbers, we present a screen of DCs and their progenitors in lymphoid and non-lymphoid tissues in Collaborative Cross (CC) and Diversity Outbred (DO) mice. We report 30 and 71 loci with logarithm of the odds (LOD) scores >8.18 and ranging from 6.67 to 8.19, respectively. The analysis reveals the highly polygenic and pleiotropic architecture of this complex trait, including many of the previously identified genetic regulators of DC development and maturation. Two SNPs in genes potentially underlying variation in DC homeostasis, a splice variant in Gramd4 (rs235532740) and a missense variant in Orai3 (rs216659754), are confirmed by gene editing using CRISPR-Cas9. Gramd4 is a central regulator of DC homeostasis that impacts the entire DC lineage, and Orai3 regulates cDC2 numbers in tissues. Overall, the data reveal a large number of candidate genes regulating DC homeostasis in vivo.

Premature skewing of T cell receptor clonality and delayed memory expansion in HIV-exposed infants

While preventing vertical HIV transmission has been very successful, HIV-exposed uninfected infants (iHEU) experience an elevated risk to infections compared to HIV-unexposed and uninfected infants (iHUU). Here we present a longitudinal multimodal analysis of infant immune ontogeny that highlights the impact of HIV/ARV exposure. Using mass cytometry, we show alterations in T cell memory differentiation between iHEU and iHUU being significant from week 15 of life. The altered memory T cell differentiation in iHEU was preceded by lower TCR Vβ clonotypic diversity and linked to TCR clonal depletion within the naïve T cell compartment. Compared to iHUU, iHEU had elevated CD56loCD16loPerforin+CD38+CD45RA+FcεRIγ+ NK cells at 1 month postpartum and whose abundance pre-vaccination were predictive of vaccine-induced pertussis and rotavirus antibody responses post 3 months of life. Collectively, HIV/ARV exposure disrupted the trajectory of innate and adaptive immunity from birth which may underlie relative vulnerability to infections in iHEU.

Monoclonal antibodies derived from B cells in subjects with cystic fibrosis reduce Pseudomonas aeruginosa burden in mice

Pseudomonas aeruginosa (PA) is an opportunistic, frequently multidrug-resistant pathogen that can cause severe infections in hospitalized patients. Antibodies against the PA virulence factor, PcrV, protect from death and disease in a variety of animal models. However, clinical trials of PcrV-binding antibody-based products have thus far failed to demonstrate benefit. Prior candidates were derivations of antibodies identified using protein-immunized animal systems and required extensive engineering to optimize binding and/or reduce immunogenicity. Of note, PA infections are common in people with cystic fibrosis (pwCF), who are generally believed to mount normal adaptive immune responses. Here we utilized a tetramer reagent to detect and isolate PcrV-specific B cells in pwCF and, via single-cell sorting and paired-chain sequencing, identified the B cell receptor (BCR) variable region sequences that confer PcrV-specificity. We derived multiple high affinity anti-PcrV monoclonal antibodies (mAbs) from PcrV-specific B cells across 3 donors, including mAbs that exhibit potent anti-PA activity in a murine pneumonia model. This robust strategy for mAb discovery expands what is known about PA-specific B cells in pwCF and yields novel mAbs with potential for future clinical use.

Identification of apolipoprotein B-reactive CDR3 motifs allows tracking of atherosclerosis-related memory CD4+T cells in multiple donors

Introduction

Atherosclerosis is a major pathological condition that underlies many cardiovascular diseases (CVDs). Its etiology involves breach of tolerance to self, leading to clonal expansion of autoreactive apolipoprotein B (APOB)-reactive CD4+T cells that correlates with clinical CVD. The T-cell receptor (TCR) sequences that mediate activation of APOB-specific CD4+T cells are unknown.

Methods

In a previous study, we had profiled the hypervariable complementarity determining region 3 (CDR3) of CD4+T cells that respond to six immunodominant APOB epitopes in most donors. Here, we comprehensively analyze this dataset of 149,065 APOB-reactive and 199,211 non-reactive control CDR3s from six human leukocyte antigen-typed donors.

Results

We identified 672 highly expanded (frequency threshold > 1.39E-03) clones that were significantly enriched in the APOB-reactive group as compared to the controls (log10 odds ratio ≥1, Fisher's test p < 0.01). Analysis of 114,755 naïve, 91,001 central memory (TCM) and 29,839 effector memory (TEM) CDR3 sequences from the same donors revealed that APOB+ clones can be traced to the complex repertoire of unenriched blood T cells. The fraction of APOB+ clones that overlapped with memory CDR3s ranged from 2.2% to 46% (average 16.4%). This was significantly higher than their overlap with the naïve pool, which ranged from 0.7% to 2% (average 1.36%). CDR3 motif analysis with the machine learning-based in-silico tool, GLIPHs (grouping of lymphocyte interactions by paratope hotspots), identified 532 APOB+ motifs. Analysis of naïve and memory CDR3 sequences with GLIPH revealed that ~40% (209 of 532) of these APOB+ motifs were enriched in the memory pool. Network analysis with Cytoscape revealed extensive sharing of the memory-affiliated APOB+ motifs across multiple donors. We identified six motifs that were present in TCM and TEM CDR3 sequences from >80% of the donors and were highly enriched in the APOB-reactive TCR repertoire.

Discussion

The identified APOB-reactive expanded CD4+T cell clones and conserved motifs can be used to annotate and track human atherosclerosis-related autoreactive CD4+T cells and measure their clonal expansion.

The multilevel extensive diversity across the cynomolgus macaque captured by ultra-deep adaptive immune receptor repertoire sequencing

The extent to which AIRRs differ among and within individuals remains elusive. Via ultra-deep repertoire sequencing of 22 and 25 tissues in three cynomolgus macaques, respectively, we identified 84 and 114 novel IGHV and TRBV alleles, confirming 72 (85.71%) and 100 (87.72%) of them. The heterogeneous V gene usage patterns were influenced, in turn, by genetics, isotype (for BCRs only), tissue group, and tissue. A higher proportion of intragroup shared clones in the intestinal tissues than those in other tissues suggests a close intra-intestinal adaptive immunity network. Significantly higher mutation burdens in the public clones and the inter-tissue shared IgM and IgD clones indicate that they might target the shared antigens. This study reveals the extensive heterogeneity of the AIRRs at various levels and has broad fundamental and clinical implications. The data generated here will serve as an invaluable resource for future studies on adaptive immunity in health and diseases.

Characterizing adjuvants' effects at murine immunoglobulin repertoire level

Generating large-scale, high-fidelity sequencing data is challenging and, furthermore, not much has been done to characterize adjuvants' effects at the repertoire level. Thus, we introduced an IgSeq pipeline that standardized library prep protocols and data analysis functions for accurate repertoire profiling. We then studied systemically effects of CpG and Alum on the Ig heavy chain repertoire using the ovalbumin (OVA) murine model. Ig repertoires of different tissues (spleen and bone marrow) and isotypes (IgG and IgM) were examined and compared in IGHV mutation, gene usage, CDR3 length, clonal diversity, and clonal selection. We found Ig repertoires of different compartments exhibited distinguishable profiles at the non-immunized steady state, and distinctions became more pronounced upon adjuvanted immunizations. Notably, Alum and CpG effects exhibited different tissue- and isotype-preferences. The former led to increased diversity of abundant clones in bone marrow, and the latter promoted the selection of IgG clones in both tissues.

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

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

Genetic variation in the immunoglobulin heavy chain locus shapes the human antibody repertoire

Variation in the antibody response has been linked to differential outcomes in disease, and suboptimal vaccine and therapeutic responsiveness, the determinants of which have not been fully elucidated. Countering models that presume antibodies are generated largely by stochastic processes, we demonstrate that polymorphisms within the immunoglobulin heavy chain locus (IGH) impact the naive and antigen-experienced antibody repertoire, indicating that genetics predisposes individuals to mount qualitatively and quantitatively different antibody responses. We pair recently developed long-read genomic sequencing methods with antibody repertoire profiling to comprehensively resolve IGH genetic variation, including novel structural variants, single nucleotide variants, and genes and alleles. We show that IGH germline variants determine the presence and frequency of antibody genes in the expressed repertoire, including those enriched in functional elements linked to V(D)J recombination, and overlapping disease-associated variants. These results illuminate the power of leveraging IGH genetics to better understand the regulation, function, and dynamics of the antibody response in disease.

Application of germline antibody features to vaccine development, antibody discovery, antibody optimization and disease diagnosis

Although the efficacy and commercial success of vaccines and therapeutic antibodies have been tremendous, designing and discovering new drug candidates remains a labor-, time- and cost-intensive endeavor with high risks. The main challenges of vaccine development are inducing a strong immune response in broad populations and providing effective prevention against a group of highly variable pathogens. Meanwhile, antibody discovery faces several great obstacles, especially the blindness in antibody screening and the unpredictability of the developability and druggability of antibody drugs. These challenges are largely due to poorly understanding of germline antibodies and the antibody responses to pathogen invasions. Thanks to the recent developments in high-throughput sequencing and structural biology, we have gained insight into the germline immunoglobulin (Ig) genes and germline antibodies and then the germline antibody features associated with antigens and disease manifestation. In this review, we firstly outline the broad associations between germline antibodies and antigens. Moreover, we comprehensively review the recent applications of antigen-specific germline antibody features, physicochemical properties-associated germline antibody features, and disease manifestation-associated germline antibody features on vaccine development, antibody discovery, antibody optimization, and disease diagnosis. Lastly, we discuss the bottlenecks and perspectives of current and potential applications of germline antibody features in the biotechnology field.

Archaic humans have contributed to large-scale variation in modern human T cell receptor genes

Human T cell receptors (TCRs) are critical for mediating immune responses to pathogens and tumors and regulating self-antigen recognition. Yet, variations in the genes encoding TCRs remain insufficiently defined. Detailed analysis of expressed TCR alpha, beta, gamma, and delta genes in 45 donors from four human populations-African, East Asian, South Asian, and European-revealed 175 additional TCR variable and junctional alleles. Most of these contained coding changes and were present at widely differing frequencies in the populations, a finding confirmed using DNA samples from the 1000 Genomes Project. Importantly, we identified three Neanderthal-derived, introgressed TCR regions including a highly divergent TRGV4 variant, which mediated altered butyrophilin-like molecule 3 (BTNL3) ligand reactivity and was frequent in all modern Eurasian population groups. Our results demonstrate remarkable variation in TCR genes in both individuals and populations, providing a strong incentive for including allelic variation in studies of TCR function in human biology.

The evolutionary and functional significance of germline immunoglobulin gene variation

The recombination between immunoglobulin (IG) gene segments determines an individual's naïve antibody repertoire and, consequently, (auto)antigen recognition. Emerging evidence suggests that mammalian IG germline variation impacts humoral immune responses associated with vaccination, infection, and autoimmunity - from the molecular level of epitope specificity, up to profound changes in the architecture of antibody repertoires. These links between IG germline variants and immunophenotype raise the question on the evolutionary causes and consequences of diversity within IG loci. We discuss why the extreme diversity in IG loci remains a mystery, why resolving this is important for the design of more effective vaccines and therapeutics, and how recent evidence from multiple lines of inquiry may help us do so.

Exploration of shared features of B cell receptor and T cell receptor repertoires reveals distinct clonotype clusters

Although B cells and T cells are integral players of the adaptive immune system and act in co-dependent ways to orchestrate immune responses, existing methods to study the immune repertoire have largely focused on separate analyses of B cell receptor (BCR) and T cell receptor (TCR) repertoires. Based on our hypothesis that the shared history of immune exposures and the shared cellular machinery for recombination result in similarities between BCR and TCR repertoires in an individual, we examine any commonalities and interrelationships between BCR and TCR repertoires. We find that the BCR and TCR repertoires have covarying clonal architecture and diversity, and that the pattern of correlations appears to be altered in immune-mediated diseases. Furthermore, hierarchical clustering of public B and T cell clonotypes in both health and disease based on correlation of clonal proportion revealed distinct clusters of B and T cell clonotypes that exhibit increased sequence similarity, share motifs, and have distinct amino acid characteristics. Our findings point to common principles governing memory formation, recombination, and clonal expansion to antigens in B and T cells within an individual. A significant proportion of public BCR and TCR repertoire can be clustered into nonoverlapping and correlated clusters, suggesting a novel way of grouping B and T cell clonotypes.

Evolution and modulation of antigen-specific T cell responses in melanoma patients

Analyzing antigen-specific T cell responses at scale has been challenging. Here, we analyze three types of T cell receptor (TCR) repertoire data (antigen-specific TCRs, TCR-repertoire, and single-cell RNA + TCRαβ-sequencing data) from 515 patients with primary or metastatic melanoma and compare it to 783 healthy controls. Although melanoma-associated antigen (MAA) -specific TCRs are restricted to individuals, they share sequence similarities that allow us to build classifiers for predicting anti-MAA T cells. The frequency of anti-MAA T cells distinguishes melanoma patients from healthy and predicts metastatic recurrence from primary melanoma. Anti-MAA T cells have stem-like properties and frequent interactions with regulatory T cells and tumor cells via Galectin9-TIM3 and PVR-TIGIT -axes, respectively. In the responding patients, the number of expanded anti-MAA clones are higher after the anti-PD1(+anti-CTLA4) therapy and the exhaustion phenotype is rescued. Our systems immunology approach paves the way for understanding antigen-specific responses in human disorders.

Previous studies have characterized the diversity and dynamics of the T cell receptor (TCR) repertoire in patients with solid cancer. Here, by analyzing TCR repertoire data from multiple datasets, the authors report that melanoma-associated antigen-specific TCRs can be used to separate metastatic melanoma patients from healthy controls and to follow anti-tumor responses in patients treated with immunotherapy.

AIRRscape: An interactive tool for exploring B-cell receptor repertoires and antibody responses

The sequencing of antibody repertoires of B-cells at increasing coverage and depth has led to the identification of vast numbers of immunoglobulin heavy and light chains. However, the size and complexity of these Adaptive Immune Receptor Repertoire sequencing (AIRR-seq) datasets makes it difficult to perform exploratory analyses. To aid in data exploration, we have developed AIRRscape, an R Shiny-based interactive web browser application that enables B-cell receptor (BCR) and antibody feature discovery through comparisons among multiple repertoires. Using AIRR-seq data as input, AIRRscape starts by aggregating and sorting repertoires into interactive and explorable bins of germline V-gene, germline J-gene, and CDR3 length, providing a high-level view of the entire repertoire. Interesting subsets of repertoires can be quickly identified and selected, and then network topologies of CDR3 motifs can be generated for further exploration. Here we demonstrate AIRRscape using patient BCR repertoires and sequences of published monoclonal antibodies to investigate patterns of humoral immunity to three viral pathogens: SARS-CoV-2, HIV-1, and DENV (dengue virus). AIRRscape reveals convergent antibody sequences among datasets for all three pathogens, although HIV-1 antibody datasets display limited convergence and idiosyncratic responses. We have made AIRRscape available as a web-based Shiny application, along with code on GitHub to encourage its open development and use by immuno-informaticians, virologists, immunologists, vaccine developers, and other scientists that are interested in exploring and comparing multiple immune receptor repertoires.

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.

Haematopoietic Stem Cell Transplantation Results in Extensive Remodelling of the Clonal T Cell Repertoire in Multiple Sclerosis

Autologous haematopoietic stem cell transplantation (AHSCT) is a vital therapeutic option for patients with highly active multiple sclerosis (MS). Rates of remission suggest AHSCT is the most effective form of immunotherapy in controlling the disease. Despite an evolving understanding of the biology of immune reconstitution following AHSCT, the mechanism by which AHSCT enables sustained disease remission beyond the period of lymphopenia remains to be elucidated. Auto-reactive T cells are considered central to MS pathogenesis. Here, we analyse T cell reconstitution for 36 months following AHSCT in a cohort of highly active MS patients. Through longitudinal analysis of sorted naïve and memory T cell clones, we establish that AHSCT induces profound changes in the dominant T cell landscape of both CD4+ and CD8+ memory T cell clones. Lymphopenia induced homeostatic proliferation is followed by clonal attrition; with only 19% of dominant CD4 (p <0.025) and 13% of dominant CD8 (p <0.005) clones from the pre-transplant repertoire detected at 36 months. Recovery of a thymically-derived CD4 naïve T cell repertoire occurs at 12 months and is ongoing at 36 months, however diversity of the naïve populations is not increased from baseline suggesting the principal mechanism of durable remission from MS after AHSCT relates to depletion of putative auto-reactive clones. In a cohort of MS patients expressing the MS risk allele HLA DRB1*15:01, public clones are probed as potential biomarkers of disease. AHSCT appears to induce sustained periods of disease remission with dynamic changes in the clonal T cell repertoire out to 36 months post-transplant.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Adaptive Immune Receptor Repertoire (AIRR) Community Guide to TR and IG Gene Annotation

High-throughput sequencing of adaptive immune receptor repertoires (AIRR, i.e., IG and TR) has revolutionized the ability to carry out large-scale experiments to study the adaptive immune response. Since the method was first introduced in 2009, AIRR sequencing (AIRR-Seq) has been applied to survey the immune state of individuals, identify antigen-specific or immune-state-associated signatures of immune responses, study the development of the antibody immune response, and guide the development of vaccines and antibody therapies. Recent advancements in the technology include sequencing at the single-cell level and in parallel with gene expression, which allows the introduction of multi-omics approaches to understand in detail the adaptive immune response. Analyzing AIRR-seq data can prove challenging even with high-quality sequencing, in part due to the many steps involved and the need to parameterize each step. In this chapter, we outline key factors to consider when preprocessing raw AIRR-Seq data and annotating the genetic origins of the rearranged receptors. We also highlight a number of common difficulties with common AIRR-seq data processing and provide strategies to address them.

Individualized VDJ recombination predisposes the available Ig sequence space

The process of recombination between variable (V), diversity (D), and joining (J) immunoglobulin (Ig) gene segments determines an individual's naive Ig repertoire and, consequently, (auto)antigen recognition. VDJ recombination follows probabilistic rules that can be modeled statistically. So far, it remains unknown whether VDJ recombination rules differ between individuals. If these rules differed, identical (auto)antigen-specific Ig sequences would be generated with individual-specific probabilities, signifying that the available Ig sequence space is individual specific. We devised a sensitivity-tested distance measure that enables inter-individual comparison of VDJ recombination models. We discovered, accounting for several sources of noise as well as allelic variation in Ig sequencing data, that not only unrelated individuals but also human monozygotic twins and even inbred mice possess statistically distinguishable immunoglobulin recombination models. This suggests that, in addition to genetic, there is also nongenetic modulation of VDJ recombination. We demonstrate that population-wide individualized VDJ recombination can result in orders of magnitude of difference in the probability to generate (auto)antigen-specific Ig sequences. Our findings have implications for immune receptor-based individualized medicine approaches relevant to vaccination, infection, and autoimmunity.

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

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

Germline polymorphisms and alternative splicing of human immunoglobulin light chain genes

Inference of germline polymorphisms in immunoglobulin genes from B cell receptor repertoires is complicated by somatic hypermutations, sequencing/PCR errors, and by varying length of reference alleles. The light chain inference is particularly challenging owing to large gene duplications and absence of D genes. We analyzed the light chain cDNA sequences from naïve B cell receptor repertoires from 100 individuals. We optimized light chain allele inference by tweaking parameters of the TIgGER functions, extending the germline reference sequences, and establishing mismatch frequency patterns at polymorphic positions to filter out false-positive candidates. We identified 48 previously unreported variants of light chain variable genes. We selected 14 variants for validation and successfully validated 11 by Sanger sequencing. Clustering of light chain 5'UTR, L-PART1, and L-PART2 revealed partial intron retention in 11 kappa and 9 lambda V alleles. Our results provide insight into germline variation in human light chain immunoglobulin loci.

Different B cell subpopulations show distinct patterns in their IgH repertoire metrics

Several human B cell subpopulations are recognised in the peripheral blood, which play distinct roles in the humoral immune response. These cells undergo developmental and maturational changes involving VDJ recombination, somatic hypermutation and class switch recombination, altogether shaping their immunoglobulin heavy chain (IgH) repertoire. Here, we sequenced the IgH repertoire of naïve, marginal zone, switched and plasma cells from 10 healthy adults along with matched unsorted and in silico separated CD19⁺ bulk B cells. Using advanced bioinformatic analysis and machine learning, we show that sorted B cell subpopulations are characterised by distinct repertoire characteristics on both the individual sequence and the repertoire level. Sorted subpopulations shared similar repertoire characteristics with their corresponding in silico separated subsets. Furthermore, certain IgH repertoire characteristics correlated with the position of the constant region on the IgH locus. Overall, this study provides unprecedented insight over mechanisms of B cell repertoire control in peripherally circulating B cell subpopulations.

Hiding in plain sight: structure and sequence analysis reveals the importance of the antibody DE loop for antibody-antigen binding

Antibody variable domains contain “complementarity-determining regions” (CDRs), the loops that form the antigen binding site. CDRs1-3 are recognized as the canonical CDRs. However, a fourth loop sits adjacent to CDR1 and CDR2 and joins the D and E strands on the antibody v-type fold. This “DE loop” is usually treated as a framework region, even though mutations in the loop affect the conformation of the CDRs and residues in the DE loop occasionally contact antigen. We analyzed the length, structure, and sequence features of all DE loops in the Protein Data Bank (PDB), as well as millions of sequences from HIV-1 infected and naïve patients. We refer to the DE loop as H4 and L4 in the heavy and light chains, respectively. Clustering the backbone conformations of the most common length of L4 (6 residues) reveals four conformations: two κ-only clusters, one λ-only cluster, and one mixed κ/λ cluster. Most H4 loops are length-8 and exist primarily in one conformation; a secondary conformation represents a small fraction of H4-8 structures. H4 sequence variability exceeds that of the antibody framework in naïve human high-throughput sequences, and both L4 and H4 sequence variability from λ and heavy germline sequences exceed that of germline framework regions. Finally, we identified dozens of structures in the PDB with insertions in the DE loop, all related to broadly neutralizing HIV-1 antibodies (bNabs), as well as antibody sequences from high-throughput sequencing studies of HIV-infected individuals, illuminating a possible role in humoral immunity to HIV-1.

Landscape of T-cell repertoires with public COVID-19-associated T-cell receptors in pre-pandemic risk cohorts

Objectives

T cells have an essential role in the antiviral defence. Public T-cell receptor (TCR) clonotypes are expanded in a substantial proportion of COVID-19 patients. We set out to exploit their potential use as read-out for COVID-19 T-cell immune responses.

Methods

We searched for COVID-19-associated T-cell clones with public TCRs, as defined by identical complementarity-determining region 3 (CDR3) beta chain amino acid sequence that can be reproducibly detected in the blood of COVID-19 patients. Of the different clonotype identification algorithms used in this study, deep sequencing of brain tissue of five patients with fatal COVID-19 delivered 68 TCR clonotypes with superior representation across 140 immune repertoires of unrelated COVID-19 patients.

Results

Mining of immune repertoires from subjects not previously exposed to the virus showed that these clonotypes can be found in almost 20% of pre-pandemic immune repertoires of healthy subjects, with lower representation in repertoires from risk groups like individuals above the age of 60 years or patients with cancer.

Conclusion

Together, our data show that at least a proportion of the SARS-CoV-2 T-cell response is mediated by public TCRs that are present in repertoires of unexposed individuals. The lower representation of these clones in repertoires of risk groups or failure to expand such clones may contribute to more unfavorable clinical COVID-19 courses.

Diversity and shared T-cell receptor repertoire analysis in esophageal squamous cell carcinoma

The tumor immune response is dependent on the interaction between tumor cells and the T-cell subset expressing the T-cell receptor (TCR) repertoire that infiltrates into the tumor microenvironment. The present study explored the diversity and shared TCR repertoires expressed on the surface of locoregional T cells and identified the T lymphocyte subsets infiltrating into esophageal squamous cell carcinoma (ESCC), in order to provide insight into the efficiency of immunotherapy and the development of a novel immune-oriented therapeutic strategy. A total of 53 patients with ESCC were enrolled in the present study, and immunohistochemical analysis of CD3, CD8, CD45RO, FOXP3, CD274, HLA class I and AE1/AE3 was performed. Digital pathological assessment was performed to evaluate the expression level of each marker. The clinicopathological significance of the immuno relation high (IR-Hi) group was assessed. Adaptor ligation PCR and next-generation sequencing were performed to explore the diversity of the TCR repertoire and to investigate the shared TCR repertoire in the IR-Hi group. Repertoire dissimilarity index (RDI) analysis was performed to assess the diversity of TCR, and the existence of shared TCRα and TCRβ was also investigated. Further stratification was performed according to the expression of markers of different T-cell subsets. Patients were stratified into IR-Hi and immuno relation low (IR-Lo) groups. Cancer-specific survival and recurrence-free survival rates were significantly improved in the IR-Hi group compared with in the IT-Lo group. The diversity of the TCR repertoire was significantly higher in the IR-Hi group. TCR repertoire analysis revealed 27 combinations of TCRα and 23 combinations of TCRβ VJ regions that were shared among the IR-Hi group. The IR-Hi group was divided into three clusters. Overall, the current findings revealed that the IR-Hi group maintained the diversity of TCR, and a portion of the IR-Hi cases held the T cells with shared TCR repertoires, implying recognition of shared antigens. The prognosis of patients with ESCC was affected by the existence of immune response cells and may possibly be stratified by the T-cell subsets.

Igh Locus Polymorphism May Dictate Topological Chromatin Conformation and V Gene Usage in the Ig Repertoire

Vast repertoires of unique antigen receptors are created in developing B and T lymphocytes. The antigen receptor loci contain many variable (V), diversity (D) and joining (J) gene segments that are arrayed across very large genomic expanses and are joined to form variable-region exons of expressed immunoglobulins and T cell receptors. This process creates the potential for an organism to respond to large numbers of different pathogens. Here, we consider the possibility that genetic polymorphisms with alterations in a vast array of regulatory elements in the immunoglobulin heavy chain (IgH) locus lead to changes in locus topology and impact immune-repertoire formation.

Profiling of T Cell Repertoire in SARS-CoV-2-Infected COVID-19 Patients Between Mild Disease and Pneumonia

The coronavirus disease 2019 (COVID-19) pandemic, caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has become a public health emergency. The most common symptoms of COVID-19 are fever, cough, and fatigue. While most patients with COVID-19 present with mild illness, some patients develop pneumonia, an important risk factor for mortality, at early stage of viral infection, putting these patients at increased risk of death. So far, little has been known about differences in the T cell repertoires between COVID-19 patients with and without pneumonia during SARS-CoV-2 infection. Herein, we aimed to investigate T cell receptor (TCR) repertoire profiles and patient-specific SARS-CoV-2-associated TCR clusters between COVID-19 patients with mild disease (no sign of pneumonia) and pneumonia. The TCR sequencing was conducted to characterize the peripheral TCR repertoire profile and diversity. The TCR clustering and CDR3 annotation were exploited to further discover groups of patient-specific TCR clonotypes with potential SARS-CoV-2 antigen specificities. Our study indicated a slight decrease in the TCR repertoire diversity and a skewed CDR3 length usage in patients with pneumonia compared to those with mild disease. The SARS-CoV-2-associated TCR clusters enriched in patients with mild disease exhibited significantly higher TCR generation probabilities and most of which were highly shared among patients, compared with those from pneumonia patients. Importantly, using similarity network-based clustering followed by the sequence conservation analysis, we found different patterns of CDR3 sequence motifs between mild disease- and pneumonia-specific SARS-CoV-2-associated public TCR clusters. Our results showed that characteristics of overall TCR repertoire and SARS-CoV-2-associated TCR clusters/clonotypes were divergent between COVID-19 patients with mild disease and patients with pneumonia. These findings provide important insights into the correlation between the TCR repertoire and disease severity in COVID-19 patients.

Large-scale analysis of 2,152 Ig-seq datasets reveals key features of B cell biology and the antibody repertoire

Antibody repertoire sequencing enables researchers to acquire millions of B cell receptors and investigate these molecules at the single-nucleotide level. This power and resolution in studying humoral responses have led to its wide applications. However, most of these studies were conducted with a limited number of samples. Given the extraordinary diversity, assessment of these key features with a large sample set is demanded. Thus, we collect and systematically analyze 2,152 high-quality heavy-chain antibody repertoires. Our study reveals that 52 core variable genes universally contribute to more than 99% of each individual's repertoire; a distal interspersed preferences characterize V gene recombination; the number of public clones between two repertoires follows a linear model, and the positive selection dominates at RGYW motif in somatic hypermutations. Thus, this population-level analysis resolves some critical features of the antibody repertoire and may have significant value to the large cadre of scientists.

Selection influences naive CD8+ TCR-β repertoire sharing

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

Machine Learning Analysis of Naïve B-Cell Receptor Repertoires Stratifies Celiac Disease Patients and Controls

Celiac disease (CeD) is a common autoimmune disorder caused by an abnormal immune response to dietary gluten proteins. The disease has high heritability. HLA is the major susceptibility factor, and the HLA effect is mediated via presentation of deamidated gluten peptides by disease-associated HLA-DQ variants to CD4+ T cells. In addition to gluten-specific CD4+ T cells the patients have antibodies to transglutaminase 2 (autoantigen) and deamidated gluten peptides. These disease-specific antibodies recognize defined epitopes and they display common usage of specific heavy and light chains across patients. Interactions between T cells and B cells are likely central in the pathogenesis, but how the repertoires of naïve T and B cells relate to the pathogenic effector cells is unexplored. To this end, we applied machine learning classification models to naïve B cell receptor (BCR) repertoires from CeD patients and healthy controls. Strikingly, we obtained a promising classification performance with an F1 score of 85%. Clusters of heavy and light chain sequences were inferred and used as features for the model, and signatures associated with the disease were then characterized. These signatures included amino acid (AA) 3-mers with distinct bio-physiochemical characteristics and enriched V and J genes. We found that CeD-associated clusters can be identified and that common motifs can be characterized from naïve BCR repertoires. The results may indicate a genetic influence by BCR encoding genes in CeD. Analysis of naïve BCRs as presented here may become an important part of assessing the risk of individuals to develop CeD. Our model demonstrates the potential of using BCR repertoires and in particular, naïve BCR repertoires, as disease susceptibility markers.

Generation of self-reactive, shared T-cell receptor α chains in the human thymus

The T-cell receptor (TCR) repertoire is generated in a semistochastic process of gene recombination and pairing of TCRα to TCRβ chains with the estimated total TCR diversity of >10⁸. Despite this high diversity, similar or identical TCR chains are found to recur in immune responses. Here, we analyzed the thymic generation of TCR sequences previously associated with recognition of self- and nonself-antigens, represented by sequences associated with autoimmune diabetes and HIV, respectively. Unexpectedly, in the CD4⁺ compartment TCRα chains associated with the recognition of self-antigens were generated in significantly higher numbers than TCRα chains associated with the recognition of nonself-antigens. The analysis of the circulating repertoire further showed that these chains are not lost in negative selection nor predominantly converted to the regulatory T-cell lineage. The high abundance of self-reactive TCRα chains in multiple individuals suggests that the human thymus has a predilection to generate self-reactive TCRα chains independently of the HLA-type and that the individual risk of autoimmunity may be modulated by the TCRβ repertoire associated with these chains.

VJ Segment Usage of TCR-Beta Repertoire in Monozygotic Cystic Fibrosis Twins

Sixteen monozygotic cystic fibrosis (CF) twin pairs of whom 14 pairs were homozygous for the most common p.Phe508del CFTR mutation were selected from the European Cystic Fibrosis Twin and Sibling Study Cohort. The monozygotic twins were examined in their T cell receptor (TCR) repertoire in peripheral blood by amplicon sequencing of the CDR3 variable region of the ß-chain. The recruitment of TCR J and V genes for recombination and selection in the thymus showed a strong genetic influence in the CF twin cohort as indicated by the shortest Jensen-Shannon distance to the twin individual. Exceptions were the clinically most discordant and/or most severely affected twin pairs where clonal expansion probably caused by recurrent pulmonary infections overshadowed the impact of the identical genomic blueprint. In general the Simpson clonality was low indicating that the population of TCRß clonotypes of the CF twins was dominated by the naïve T-cell repertoire. Intrapair sharing of clonotypes was significantly more frequent among monozygotic CF twins than among pairs of unrelated CF patients. Complete nucleotide sequence identity was observed in about 0.11% of CDR3 sequences which partially should represent persisting fetal clones derived from the same progenitor T cells. Complete amino acid sequence identity was noted in 0.59% of clonotypes. Of the nearly 40,000 frequent amino acid clonotypes shared by at least two twin siblings 99.8% were already known within the immuneACCESS database and only 73 had yet not been detected indicating that the CDR3ß repertoire of CF children and adolescents does not carry a disease-specific signature but rather shares public clones with that of the non-CF community. Clonotypes shared within twin pairs and between unrelated CF siblings were highly abundant among healthy non-CF people, less represented in individuals with infectious disease and uncommon in patients with cancer. This subset of shared CF clonotypes defines CDR3 amino acid sequences that are more common in health than in disease.

Immunoglobulin germline gene variation and its impact on human disease

Immunoglobulins (Ig) play an important role in the immune system both when expressed as antigen receptors on the cell surface of B cells and as antibodies secreted into extracellular fluids. The advent of high-throughput sequencing methods has enabled the investigation of human Ig repertoires at unprecedented depth. This has led to the discovery of many previously unreported germline Ig alleles. Moreover, it is becoming clear that convergent and stereotypic antibody responses are common where different individuals recognise defined antigenic epitopes with the use of the same Ig V genes. Thus, germline V gene variation is increasingly being linked to the differential capacity of generating an effective immune response, which might lead to varying disease susceptibility. Here, we review recent evidence of how germline variation in Ig genes impacts the Ig repertoire and its subsequent effects on the adaptive immune response in vaccination, infection, and autoimmunity.

Antigen receptor repertoires of one of the smallest known vertebrates

The rules underlying the structure of antigen receptor repertoires are not yet fully defined, despite their enormous importance for the understanding of adaptive immunity. With current technology, the large antigen receptor repertoires of mice and humans cannot be comprehensively studied. To circumvent the problems associated with incomplete sampling, we have studied the immunogenetic features of one of the smallest known vertebrates, the cyprinid fish Paedocypris sp. "Singkep" ("minifish"). Despite its small size, minifish has the key genetic facilities characterizing the principal vertebrate lymphocyte lineages. As described for mammals, the frequency distributions of immunoglobulin and T cell receptor clonotypes exhibit the features of fractal systems, demonstrating that self-similarity is a fundamental property of antigen receptor repertoires of vertebrates, irrespective of body size. Hence, minifish achieve immunocompetence via a few thousand lymphocytes organized in robust scale-free networks, thereby ensuring immune reactivity even when cells are lost or clone sizes fluctuate during immune responses.

The adaptive immune receptor repertoire community as a model for FAIR stewardship of big immunology data

Systems biology involves network-oriented, computational approaches to modeling biological systems through analysis of big biological data. To contribute maximally to scientific progress, big biological data should be FAIR: findable, accessible, interoperable, and reusable. Here, we describe high-throughput sequencing data that characterize the vast diversity of B- and T-cell clones comprising the adaptive immune receptor repertoire (AIRR-seq data) and its contribution to our understanding of COVID-19 (coronavirus disease 19). We describe the accomplishments of the AIRR community, a grass-roots network of interdisciplinary laboratory scientists, bioinformaticians, and policy wonks, in creating and publishing standards, software and repositories for AIRR-seq data based on the FAIR principles.

Position-Dependent Differential Targeting of Somatic Hypermutation

Somatic hypermutation (SHM) generates much of the Ab diversity necessary for affinity maturation and effective humoral immunity. The activation-induced cytidine deaminase-induced DNA lesions and error-prone repair that underlie SHM are known to exhibit intrinsic biases when targeting the Ig sequences. Computational models for SHM targeting often model the targeting probability of a nucleotide in a motif-based fashion, assuming that the same DNA motif is equally likely to be targeted regardless of its position along the Ig sequence. The validity of this assumption, however, has not been rigorously studied in vivo. In this study, by analyzing a large collection of 956,157 human Ig sequences while controlling for the confounding influence of selection, we show that the likelihood of a DNA 5-mer motif being targeted by SHM is not the same at different positions in the same Ig sequence. We found position-dependent differential SHM targeting for about three quarters of the 38 and 269 unique motifs from more than half of the 292 and 1912 motif-allele pairs analyzed using productive and nonproductive Ig sequences, respectively. The direction of the differential SHM targeting was largely conserved across individuals with no allele-specific effect within an IgH variable gene family, but was not consistent with general decay of SHM targeting with increasing distance from the transcription start site. However, SHM targeting did correlate positively with the mutability of the wider sequence neighborhood surrounding the motif. These findings provide insights and future directions for computational efforts toward modeling SHM.

Analyzing the Mycobacterium tuberculosis immune response by T-cell receptor clustering with GLIPH2 and genome-wide antigen screening

CD4+ T cells are critical to fighting pathogens, but a comprehensive analysis of human T-cell specificities is hindered by the diversity of HLA alleles (>20,000) and the complexity of many pathogen genomes. We previously described GLIPH, an algorithm to cluster T-cell receptors (TCRs) that recognize the same epitope and to predict their HLA restriction, but this method loses efficiency and accuracy when >10,000 TCRs are analyzed. Here we describe an improved algorithm, GLIPH2, that can process millions of TCR sequences. We used GLIPH2 to analyze 19,044 unique TCRβ sequences from 58 individuals latently infected with Mycobacterium tuberculosis (Mtb) and to group them according to their specificity. To identify the epitopes targeted by clusters of Mtb-specific T cells, we carried out a screen of 3,724 distinct proteins covering 95% of Mtb protein-coding genes using artificial antigen-presenting cells (aAPCs) and reporter T cells. We found that at least five PPE (Pro-Pro-Glu) proteins are targets for T-cell recognition in Mtb.

Intra-Vκ Cluster Recombination Shapes the Ig Kappa Locus Repertoire

During V(D)J recombination, RAG proteins introduce DNA double-strand breaks (DSBs) at recombination signal sequences (RSSs) that contain either 12- or 23-nt spacer regions. Coordinated 12/23 cleavage predicts that DSBs at variable (V) gene segments should equal the level of breakage at joining (J) segments. Contrary to this, here we report abundant RAG-dependent DSBs at multiple Vκ gene segments independent of V-J rearrangement. We find that a large fraction of Vκ gene segments are flanked not only by a bone-fide 12 spacer but also an overlapping, 23-spacer flipped RSS. These compatible pairs of RSSs mediate recombination and deletion inside the Vκ cluster even in the complete absence of Jκ gene segments and support a V(D)J recombination center (RC) independent of the conventional Jκ-centered RC. We propose an improved model of Vκ-Jκ repertoire formation by incorporating these surprisingly frequent, evolutionarily conserved intra-Vκ cluster recombination events.

A Novel Framework for Characterizing Genomic Haplotype Diversity in the Human Immunoglobulin Heavy Chain Locus

An incomplete ascertainment of genetic variation within the highly polymorphic immunoglobulin heavy chain locus (IGH) has hindered our ability to define genetic factors that influence antibody-mediated processes. Due to locus complexity, standard high-throughput approaches have failed to accurately and comprehensively capture IGH polymorphism. As a result, the locus has only been fully characterized two times, severely limiting our knowledge of human IGH diversity. Here, we combine targeted long-read sequencing with a novel bioinformatics tool, IGenotyper, to fully characterize IGH variation in a haplotype-specific manner. We apply this approach to eight human samples, including a haploid cell line and two mother-father-child trios, and demonstrate the ability to generate high-quality assemblies (>98% complete and >99% accurate), genotypes, and gene annotations, identifying 2 novel structural variants and 15 novel IGH alleles. We show multiplexing allows for scaling of the approach without impacting data quality, and that our genotype call sets are more accurate than short-read (>35% increase in true positives and >97% decrease in false-positives) and array/imputation-based datasets. This framework establishes a desperately needed foundation for leveraging IG genomic data to study population-level variation in antibody-mediated immunity, critical for bettering our understanding of disease risk, and responses to vaccines and therapeutics.

The ADC API: A Web API for the Programmatic Query of the AIRR Data Commons

The Adaptive Immune Receptor Repertoire (AIRR) Community is a research-driven group that is establishing a clear set of community-accepted data and metadata standards; standards-based reference implementation tools; and policies and practices for infrastructure to support the deposit, curation, storage, and use of high-throughput sequencing data from B-cell and T-cell receptor repertoires (AIRR-seq data). The AIRR Data Commons is a distributed system of data repositories that utilizes a common data model, a common query language, and common interoperability formats for storage, query, and downloading of AIRR-seq data. Here is described the principal technical standards for the AIRR Data Commons consisting of the AIRR Data Model for repertoires and rearrangements, the AIRR Data Commons (ADC) API for programmatic query of data repositories, a reference implementation for ADC API services, and tools for querying and validating data repositories that support the ADC API. AIRR-seq data repositories can become part of the AIRR Data Commons by implementing the data model and API. The AIRR Data Commons allows AIRR-seq data to be reused for novel analyses and empowers researchers to discover new biological insights about the adaptive immune system.

Immunotyping of clinically divergent p.Phe508del homozygous monozygous cystic fibrosis twins

Blood of the three clinically most concordant and most discordant p.Phe508del homozygous monozygous twin pairs of the European Cystic Fibrosis Twin and Sibling Study was examined in two postzygotic attributes that generate diversity between monozygous twins, i.e. the repertoire of the CDR3 region of the T-cell receptor ß chains and the DNA methylation at 450,000 genomic CpG sites. Methylation patterns in peripheral blood of twins changed at selected cell-type-independent positions and the immune cells of the twins showed individual profiles of the T cell receptor repertoire reflecting the plasticity of the immune system of genetically identical humans with cystic fibrosis to cope with the environment.

Determinants governing T cell receptor α/β-chain pairing in repertoire formation of identical twins

The T cell repertoire in each individual includes T cell receptors (TCRs) of enormous sequence diversity through the pairing of diverse TCR α- and β-chains, each generated by somatic recombination of paralogous gene segments. Whether the TCR repertoire contributes to susceptibility to infectious or autoimmune diseases in concert with disease-associated major histocompatibility complex (MHC) polymorphisms is unknown. Due to a lack in high-throughput technologies to sequence TCR α-β pairs, current studies on whether the TCR repertoire is shaped by host genetics have so far relied only on single-chain analysis. Using a high-throughput single T cell sequencing technology, we obtained the largest paired TCRαβ dataset so far, comprising 965,523 clonotypes from 15 healthy individuals including 6 monozygotic twin pairs. Public TCR α- and, to a lesser extent, TCR β-chain sequences were common in all individuals. In contrast, sharing of entirely identical TCRαβ amino acid sequences was very infrequent in unrelated individuals, but highly increased in twins, in particular in CD4 memory T cells. Based on nucleotide sequence identity, a subset of these shared clonotypes appeared to be the progeny of T cells that had been generated during fetal development and had persisted for more than 50 y. Additional shared TCRαβ in twins were encoded by different nucleotide sequences, implying that genetic determinants impose structural constraints on thymic selection that favor the selection of TCR α-β pairs with entire sequence identities.

sumrep: A Summary Statistic Framework for Immune Receptor Repertoire Comparison and Model Validation

The adaptive immune system generates an incredible diversity of antigen receptors for B and T cells to keep dangerous pathogens at bay. The DNA sequences coding for these receptors arise by a complex recombination process followed by a series of productivity-based filters, as well as affinity maturation for B cells, giving considerable diversity to the circulating pool of receptor sequences. Although these datasets hold considerable promise for medical and public health applications, the complex structure of the resulting adaptive immune receptor repertoire sequencing (AIRR-seq) datasets makes analysis difficult. In this paper we introduce sumrep, an R package that efficiently performs a wide variety of repertoire summaries and comparisons, and show how sumrep can be used to perform model validation. We find that summaries vary in their ability to differentiate between datasets, although many are able to distinguish between covariates such as donor, timepoint, and cell type for BCR and TCR repertoires. We show that deletion and insertion lengths resulting from V(D)J recombination tend to be more discriminative characterizations of a repertoire than summaries that describe the amino acid composition of the CDR3 region. We also find that state-of-the-art generative models excel at recapitulating gene usage and recombination statistics in a given experimental repertoire, but struggle to capture many physiochemical properties of real repertoires.

The FluPRINT dataset, a multidimensional analysis of the influenza vaccine imprint on the immune system

Machine learning has the potential to identify novel biological factors underlying successful antibody responses to influenza vaccines. The first attempts have revealed a high level of complexity in establishing influenza immunity, and many different cellular and molecular components are involved. Of note is that the previously identified correlates of protection fail to account for the majority of individual responses across different age groups and influenza seasons. Challenges remain from the small sample sizes in most studies and from often limited data sets, such as transcriptomic data. Here we report the creation of a unified database, FluPRINT, to enable large-scale studies exploring the cellular and molecular underpinnings of successful antibody responses to influenza vaccines. Over 3,000 parameters were considered, including serological responses to influenza strains, serum cytokines, cell phenotypes, and cytokine stimulations. FluPRINT, facilitates the application of machine learning algorithms for data mining. The data are publicly available and represent a resource to uncover new markers and mechanisms that are important for influenza vaccine immunogenicity.

Measurement(s)	immune response trait
Technology Type(s)	digital curation
Factor Type(s)	gender • race • visit_age • bmi • flu_vaccination_history • statin_use • influenza_infection_history • influenza_hospitalization • cmv_status • ebv_status
Sample Characteristic - Organism	Homo sapiens

Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.9902447

cAb-Rep: A Database of Curated Antibody Repertoires for Exploring Antibody Diversity and Predicting Antibody Prevalence

The diversity of B cell receptors provides a basis for recognizing numerous pathogens. Antibody repertoire sequencing has revealed relationships between B cell receptor sequences, their diversity, and their function in infection, vaccination, and disease. However, many repertoire datasets have been deposited without annotation or quality control, limiting their utility. To accelerate investigations of B cell immunoglobulin sequence repertoires and to facilitate development of algorithms for their analysis, we constructed a comprehensive public database of curated human B cell immunoglobulin sequence repertoires, cAb-Rep (https://cab-rep.c2b2.columbia.edu), which currently includes 306 immunoglobulin repertoires from 121 human donors, who were healthy, vaccinated, or had autoimmune disease. The database contains a total of 267.9 million V(D)J heavy chain and 72.9 million VJ light chain transcripts. These transcripts are full-length or near full-length, have been annotated with gene origin, antibody isotype, somatic hypermutations, and other biological characteristics, and are stored in FASTA format to facilitate their direct use by most current repertoire-analysis programs. We describe a website to search cAb-Rep for similar antibodies along with methods for analysis of the prevalence of antibodies with specific genetic signatures, for estimation of reproducibility of somatic hypermutation patterns of interest, and for delineating frequencies of somatically introduced N-glycosylation. cAb-Rep should be useful for investigating attributes of B cell sequence repertoires, for understanding characteristics of affinity maturation, and for identifying potential barriers to the elicitation of effective neutralizing antibodies in infection or by vaccination.

A TRAV26-1-encoded recognition motif focuses the biased T cell response in celiac disease

The semi-public T-cell response towards the gluten epitope DQ2.5-glia-α2 uses a prototypic TCR encoded by the germline segments TRAV26-1 and TRBV7-2. Through mutagenesis experiments, we show that a TRAV26-1encoded recognition motif contacts the MHC β-chain and the TCR CDR3β loop underpinning this conserved T-cell response restricted to the prototypic TCRs.