Ultrasensitive allele inference from immune repertoire sequencing data with MiXCR.

Allelic variability in the adaptive immune receptor loci, which harbor the gene segments that encode B cell and T cell receptors (BCR/TCR), has been shown to be of critical importance for immune responses to pathogens and vaccines. In recent years, B cell and T cell receptor repertoire sequencing (Rep-Seq) has become widespread in immunology research making it the most readily available source of information about allelic diversity in immunoglobulin (IG) and T cell receptor (TR) loci in different populations. Here we present a novel algorithm for extra-sensitive and specific variable (V) and joining (J) gene allele inference and genotyping allowing reconstruction of individual high-quality gene segment libraries. The approach can be applied for inferring allelic variants from peripheral blood lymphocyte BCR and TCR repertoire sequencing data, including hypermutated isotype-switched BCR sequences, thus allowing high-throughput genotyping and novel allele discovery from a wide variety of existing datasets. The developed algorithm is a part of the MiXCR software ( https://mixcr.com ) and can be incorporated into any pipeline utilizing upstream processing with MiXCR. We demonstrate the accuracy of this approach using Rep-Seq paired with long-read genomic sequencing data, comparing it to a widely used algorithm, TIgGER. We applied the algorithm to a large set of IG heavy chain (IGH) Rep-Seq data from 450 donors of ancestrally diverse population groups, and to the largest reported full-length TCR alpha and beta chain (TRA; TRB) Rep-Seq dataset, representing 134 individuals. This allowed us to assess the genetic diversity of genes within the IGH, TRA and TRB loci in different populations and demonstrate the connection between antibody repertoire gene usage and the number of allelic variants present in the population. Finally we established a database of allelic variants of V and J genes inferred from Rep-Seq data and their population frequencies with free public access at https://vdj.online .

Memory persistence and differentiation into antibody-secreting cells accompanied by positive selection in longitudinal BCR repertoires.

The stability and plasticity of B cell-mediated immune memory ensures the ability to respond to the repeated challenges. We have analyzed the longitudinal dynamics of immunoglobulin heavy chain repertoires from memory B cells, plasmablasts, and plasma cells from the peripheral blood of generally healthy volunteers. We reveal a high degree of clonal persistence in individual memory B cell subsets, with inter-individual convergence in memory and antibody-secreting cells (ASCs). ASC clonotypes demonstrate clonal relatedness to memory B cells, and are transient in peripheral blood. We identify two clusters of expanded clonal lineages with differing prevalence of memory B cells, isotypes, and persistence. Phylogenetic analysis revealed signs of reactivation of persisting memory B cell-enriched clonal lineages, accompanied by new rounds of affinity maturation during proliferation and differentiation into ASCs. Negative selection contributes to both persisting and reactivated lineages, preserving the functionality and specificity of B cell receptors (BCRs) to protect against current and future pathogens.

Deconvolution of B cell receptor repertoire in multiple sclerosis patients revealed a delay in tBreg maturation.

Background: B lymphocytes play a pivotal regulatory role in the development of the immune response. It was previously shown that deficiency in B regulatory cells (Bregs) or a decrease in their anti-inflammatory activity can lead to immunological dysfunctions. However, the exact mechanisms of Bregs development and functioning are only partially resolved. For instance, only a little is known about the structure of their B cell receptor (BCR) repertoires in autoimmune disorders, including multiple sclerosis (MS), a severe neuroinflammatory disease with a yet unknown etiology. Here, we elucidate specific properties of B regulatory cells in MS. Methods: We performed a prospective study of the transitional Breg (tBreg) subpopulations with the CD19+CD24highCD38high phenotype from MS patients and healthy donors by (i) measuring their content during two diverging courses of relapsing-remitting MS: benign multiple sclerosis (BMS) and highly active multiple sclerosis (HAMS); (ii) analyzing BCR repertoires of circulating B cells by high-throughput sequencing; and (iii) measuring the percentage of CD27+ cells in tBregs. Results: The tBregs from HAMS patients carry the heavy chain with a lower amount of hypermutations than tBregs from healthy donors. The percentage of transitional CD24highCD38high B cells is elevated, whereas the frequency of differentiated CD27+ cells in this transitional B cell subset was decreased in the MS patients as compared with healthy donors. Conclusions: Impaired maturation of regulatory B cells is associated with MS progression.

Surface NKG2C Identifies Differentiated αßT-Cell Clones Expanded in Peripheral Blood.

T cells that express CD56 in peripheral blood of healthy humans represent a heterogeneous and poorly studied subset. In this work, we analyzed this subset for NKG2C expression. In both CD56+ and CD56- subsets most of the NKG2C+ T cells had a phenotype of highly differentiated CD8+ TEMRA cells. The CD56+NKG2C+ T cells also expressed a number of NK cell receptors, such as NKG2D, CD16, KIR2DL2/DL3, and maturation marker CD57 more often than the CD56-NKG2C+CD3+ cells. TCR ß-chain repertoire of the CD3+CD56+NKG2C+ cell fraction was limited by the prevalence of one or several clonotypes which can be found within the most abundant clonotypes in total or CD8+ T cell fraction TCRß repertoire. Thus, NKG2C expression in highly differentiated CD56+ T cells was associated with the most expanded αß T cell clones. NKG2C+ T cells produced almost no IFN-γ in response to stimulation with HCMV pp65-derived peptides. This may be partially due to the high content of CD45RA+CD57+ cells in the fraction. CD3+NKG2C+ cells showed signs of activation, and the frequency of this T-cell subset in HCMV-positive individuals was positively correlated with the frequency of NKG2C+ NK cells that may imply a coordinated in a certain extent development of the NKG2C+ T and NK cell subsets under HCMV infection.

Exploring the pre-immune landscape of antigen-specific T cells.

BACKGROUND: Adaptive immune responses to newly encountered pathogens depend on the mobilization of antigen-specific clonotypes from a vastly diverse pool of naive T cells. Using recent advances in immune repertoire sequencing technologies, models of the immune receptor rearrangement process, and a database of annotated T cell receptor (TCR) sequences with known specificities, we explored the baseline frequencies of T cells specific for defined human leukocyte antigen (HLA) class I-restricted epitopes in healthy individuals. METHODS: We used a database of TCR sequences with known antigen specificities and a probabilistic TCR rearrangement model to estimate the baseline frequencies of TCRs specific to distinct antigens epitopespecificT-cells. We verified our estimates using a publicly available collection of TCR repertoires from healthy individuals. We also interrogated a database of immunogenic and non-immunogenic peptides is used to link baseline T-cell frequencies with epitope immunogenicity. RESULTS: Our findings revealed a high degree of variability in the prevalence of T cells specific for different antigens that could be explained by the physicochemical properties of the corresponding HLA class I-bound peptides. The occurrence of certain rearrangements was influenced by ancestry and HLA class I restriction, and umbilical cord blood samples contained higher frequencies of common pathogen-specific TCRs. We also identified a quantitative link between specific T cell frequencies and the immunogenicity of cognate epitopes presented by defined HLA class I molecules. CONCLUSIONS: Our results suggest that the population frequencies of specific T cells are strikingly non-uniform across epitopes that are known to elicit immune responses. This inference leads to a new definition of epitope immunogenicity based on specific TCR frequencies, which can be estimated with a high degree of accuracy in silico, thereby providing a novel framework to integrate computational and experimental genomics with basic and translational research efforts in the field of T cell immunology.