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.     

A landscape for the dynamics of an immune response

Dynamic shaping of the antibody repertoire is essential for effective immunity. We describe here a novel approach for clarifying how the antibody repertoire is shaped over time for development of a specific immune response. We obtained over 500 immunoglobulin G1 clones harboring VH186.2 from the spleen, bone marrow, and microdissected individual germinal centers of (4-hydroxy-3-nitrophenyl) acetyl-immunized C57BL/6 mice at various time points postimmunization. Statistical analyses provided an index for defining clonal diversity and cluster analyses gave us a three-dimensional landscape with which clone distance was visualized with the expression level of antibodies. This landscape approach facilitated our understanding of the dynamics shaping the actual antibody repertoire, in which pre-existing naturally occurring population persisted and provided a significant impact upon the repertoire. To the established model for describing production of the antibody-forming cells, we were able to append an indispensable issue in considering the maturation of humoral immune response.     

The diversity of T cell receptors specific for self MHC gene products

Cytolytic and helper T cells exhibit, in addition to their specificity for foreign antigen, a restriction specificity for self MHC gene products. The present study was designed to assess the degree of diversity within the repertoire of receptors that are involved in T cell recognition of self MHC gene products. For this purpose, we generated a series of murine cytolytic T lymphocyte (CTL) clones specific for a hapten antigen and restricted to the self MHC gene product H-2Kb. An analysis of the hapten fine specificity of these clones by using hapten analogues revealed the presence of substantial diversity within the repertoire of CTL receptors specific for the hapten. The degree of diversity within the repertoire of self H-2 recognition structures on these clones was assessed by testing clones on panels of syngeneic, congenic H-2K disparate, and H-2Kb mutant target cells bearing varying amounts of antigen. A striking degree of heterogeneity in H-2K recognition was found among these H-2Kb restricted CTL. We estimate that there are probably a minimum of 65 different patterns of H-2K recognition among these clones. Our results suggest a high degree of diversity exists within the repertoire of self MHC recognition structures on antigen-specific T cells restricted to a single self MHC gene product.     

Resource competition determines selection of B cell repertoires

Previous experiments with mouse chimeras demonstrated that cellular competition for antigen-specific survival signals plays a crucial role in the maintenance of the naive B cell repertoire. Transgenic (Tg) B cell populations in these chimeras have a shortened lifespan and poor competitive abilities as compared to more diverse non-Tg populations in the same mice. We develop a mathematical model to investigate the mechanism of B cell competition. The model allows for various B cell clones, generated in the bone marrow, to go into the peripheral circulation, where they compete specifically for various ligands providing survival signals. In the model we also find the observed poor competitive abilities of the Tg repertoire. Investigating the nature of the competition in the model, we find that most of the competition is "intraspecific" occurring largely within the clone of truly Tg B cells, and within the repertoire of leaky Tg and non-Tg B cells. This is confirmed by analysing a simplified version of the model, which only allows for intraspecific competition, and resembles a simple ecological model with density-dependent death. The fact that our model accounts for the data, casts doubt on a previous interpretation of the same data arguing that more diverse repertoires outcompete repertoires of lower diversity. Here, we conclude that most of the data can be explained with intraspecific competition, and formulate an experimental prediction that allows one to distinguish between the previous interpretation of inter-specific competition between repertoires, and the current interpretation of intraspecific competition.     

MAJOR HISTOCOMPATIBILITY COMPLEX I POLYMORPHISM IN A CETACEAN: THE GRAY WHALE (ESCHRICHTIUS ROBUSTUS)

Vertebrate Major Histocompatibility Complex-I ( Mbc -I) proteins bind and display self and foreign peptides on the cell surface. Mbc -I polymorphism is considered critical for eliciting immune responses to a diversity of antigens, and thus, for the health and conservation of a species. Based on restriction fragment length polymorphism it was concluded that whales generally have low Mbc -I polymorphism. This was attributed to the weak pathogenic pressure in aquatic habitats. Since gray whales' habits might favor their encounter with diverse pathogens, resulting in selection pressures for Mbc variability, we searched for functional polymorphism of gray whale's Mbc -I exon 2 sequences that encode the α-1 protein domain of the antigen-binding region. We sequenced twelve Mbc -I exon 2 clones from each of seven gray whales. Most obtained sequences were similar to functional artiodactyl sequences, and their reconstructed phylogeny consistently showed three whale Mbc -I exon 2 sequence clusters (A, B, and C). Sequences within a given cluster and from distinct clusters, showed a greater number of non-synonymous than synonymous substitutions, that commonly shifted the physiochemical properties of the involved residue, suggesting the existence of a diverse repertoire of Mbc -I antigen recognition sites in gray whales, and that selection maintains gray whale's Mbc-I allelic diversity.     

Enhanced CHO Clone Screening: Application of Targeted Locus Amplification and Next‐Generation Sequencing Technologies for Cell Line Development

Early analytical clone screening is important during Chinese hamster ovary (CHO) cell line development of biotherapeutic proteins to select a clonally derived cell line with most favorable stability and product quality. Sensitive sequence confirmation methods using mass spectrometry have limitations in throughput and turnaround time. Next‐generation sequencing (NGS) technologies emerged as alternatives for CHO clone analytics. We report an efficient NGS workflow applying the targeted locus amplification (TLA) strategy for genomic screening of antibody expressing CHO clones. In contrast to previously reported RNA sequencing approaches, TLA allows for targeted sequencing of genomic integrated transgenic DNA without prior locus information, robust detection of single‐nucleotide variants (SNVs) and transgenic rearrangements. During clone selection, TLA/NGS revealed CHO clones with high‐level SNVs within the antibody gene and we report in another case the utility of TLA/NGS to identify rearrangements at transgenic DNA level. We also determined detection limits for SNVs calling and the potential to identify clone contaminations by TLA/NGS. TLA/NGS also allows to identify genetically identical clones. In summary, we demonstrate that TLA/NGS is a robust screening method useful for routine clone analytics during cell line development with the potential to process up to 24 CHO clones in less than 7 workdays.     

Low level sequence variant analysis of recombinant proteins: an optimized approach

Sequence variants in recombinant biopharmaceuticals may have a relevant and unpredictable impact on clinical safety and efficacy. Hence, their sensitive analysis is important throughout bioprocess development. The two stage analytical approach presented here provides a quick multi clone comparison of candidate production cell lines as a first stage, followed by an in-depth analysis including identification and quantitation of aberrant sequence variants of selected clones as a second stage. We show that the differential analysis is a suitable tool for sensitive and fast batch to batch comparison of recombinant proteins. The optimized approach allows for detection of not only single amino acid substitutions in unmodified peptides, but also substitutions in posttranslational modified peptides such as glycopeptides, for detection of truncated or elongated sequence variants as well as double amino acid substitutions or substitution with amino acid structural isomers within one peptide. In two case studies we were able to detect sequence variants of different origin down to a sub percentage level. One of the sequence variants (Thr → Asn) could be correlated to a cytosine to adenine substitution at DNA (desoxyribonucleic acid) level. In the second case we were able to correlate the sub percentage substitution (Phe → Tyr) to amino acid limitation in the chemically defined fermentation medium.     

MHC class I gene conversion mutations alter the CD8 T cell repertoire

MHC class I molecules are highly polymorphic within populations. This diversity is thought to be the result of selective maintenance of new class I alleles formed by gene conversion. It has been proposed that rare alleles are maintained by their ability to confer resistance to common pathogens. Investigation has focused on differences in the presentation of foreign Ags by class I alleles, but the majority of peptides presented by class I molecules are self peptides used in shaping the naive T cell repertoire. We propose that the key substrate for the natural selection of class I gene conversion variants is the diversity in immune potential formed by new alleles. We show that T cells compete with each other for niches in the thymus and spleen during development, and that competition between different clones is dramatically affected by class I mutations. We also show that peripheral naive T cells proliferate preferentially in the presence of the class I variant that directed T cell development. The data argue that class I gene conversion mutations dramatically affect both the development and the maintenance of the naive CD8 T cell repertoire.     

Trade-offs in antibody repertoires to complex antigens

Pathogens vary in their antigenic complexity. While some pathogens such as measles present a few relatively invariant targets to the immune system, others such as malaria display considerable antigenic diversity. How the immune response copes in the presence of multiple antigens, and whether a trade-off exists between the breadth and efficacy of antibody (Ab)-mediated immune responses, are unsolved problems. We present a theoretical model of affinity maturation of B-cell receptors (BCRs) during a primary infection and examine how variation in the number of accessible antigenic sites alters the Ab repertoire. Naive B cells with randomly generated receptor sequences initiate the germinal centre (GC) reaction. The binding affinity of a BCR to an antigen is quantified via a genotype–phenotype map, based on a random energy landscape, that combines local and distant interactions between residues. In the presence of numerous antigens or epitopes, B-cell clones with different specificities compete for stimulation during rounds of mutation within GCs. We find that the availability of many epitopes reduces the affinity and relative breadth of the Ab repertoire. Despite the stochasticity of somatic hypermutation, patterns of immunodominance are strongly shaped by chance selection of naive B cells with specificities for particular epitopes. Our model provides a mechanistic basis for the diversity of Ab repertoires and the evolutionary advantage of antigenically complex pathogens.     

Immune repertoire mining for rapid affinity optimization of mouse monoclonal antibodies

Traditional hybridoma and B cell cloning antibody discovery platforms have inherent limits in immune repertoire sampling depth. One consequence is that monoclonal antibody (mAb) leads often lack the necessary affinity for therapeutic applications, thus requiring labor-intensive and time-consuming affinity in vitro engineering optimization steps. Here, we show that high-affinity variants of mouse-derived mAbs can be rapidly obtained by testing of somatic sequence variants obtained by deep sequencing of antibody variable regions in immune repertories from immunized mice, even with a relatively sparse sampling of sequence variants from large sequence datasets. Affinity improvements can be achieved for mAbs with a wide range of affinities. The optimized antibody variants derived from immune repertoire mining have no detectable in vitro off-target binding and have in vivo clearance comparable to the parental mAbs, essential properties in therapeutic antibody leads. As generation of antibody variants in vitro is replaced by mining of variants generated in vivo, the procedure can be applied to rapidly identify affinity-optimized mAb variants.     

Stochastic niche structure and diversity maintenance in the T cell repertoire

The reliability of the immune response to pathogenic challenge depends critically on the size and diversity of the T cell repertoire. We study naïve T cell repertoire diversity maintenance by a stochastic model that incorporates the concept of competition between T cells for survival stimuli emanating from self-antigen presenting cells (APCs). In the mean field approximation we show that clonotype extinction is certain and compute mean extinction times. We introduce the concept of mean niche overlap and show that clones with a mean niche overlap greater than one have a short repertoire lifespan. This selection differential induces minimal recognition commonality between T cell receptors (TCRs) resulting in a diverse T cell repertoire.     

Ageing of the B-cell repertoire

Older people are more susceptible to infection, less responsive to vaccination and have a more inflammatory immune environment. Using spectratype analysis, we have previously shown that the B-cell repertoire of older people shows evidence of inappropriate clonal expansions in the absence of challenge, and that this loss of B-cell diversity correlates with poor health. Studies on response to vaccination, using both spectratyping and high-throughput sequencing of the repertoire, indicate that older responses to challenge are lacking in magnitude and/or delayed significantly. Also that some of the biologically significant differences may be in different classes of antibody. We have also previously shown that normal young B-cell repertoires can vary between different phenotypic subsets of B cells. In this paper, we present an analysis of immunoglobulin repertoire in different subclasses of antibody in five different populations of B cell, and show how the repertoire in these different groups changes with age. Although some age-related repertoire differences occur in naive cells, before exogenous antigen exposure, we see indications that there is a general dysregulation of the selective forces that shape memory B-cell populations in older people.     

Characterization of human IgM and IgG repertoires in individuals with chronic HIV-1 infection

Advancements in high-throughput sequencing (HTS) of antibody repertoires (Ig-Seq) have unprecedentedly improved our ability to characterize the antibody repertoires on a large scale. However, currently, only a few studies explored the influence of chronic HIV-1 infection on human antibody repertoires and many of them reached contradictory conclusions, possibly limited by inadequate sequencing depth and throughput. To better understand how HIV-1 infection would impact humoral immune system, in this study, we systematically analyzed the differences between the IgM (HIV-IgM) and IgG (HIV-IgG) heavy chain repertoires of HIV-1 infected patients, as well as between antibody repertoires of HIV-1 patients and healthy donors (HH). Notably, the public unique clones accounted for only a negligible proportion between the HIV-IgM and HIV-IgG repertoires libraries, and the diversity of unique clones in HIV-IgG remarkably reduced. In aspect of somatic mutation rates of CDR1 and CDR2, the HIV-IgG repertoire was higher than HIV-IgM. Besides, the average length of CDR3 region in HIV-IgM was significant longer than that in the HH repertoire, presumably caused by the great number of novel VDJ rearrangement patterns, especially a massive use of IGHJ6. Moreover, some of the B cell clonotypes had numerous clones, and somatic variants were detected within the clonotype lineage in HIV-IgG, indicating HIV-1 neutralizing activities. The in-depth characterization of HIV-IgG and HIV-IgM repertoires enriches our knowledge in the profound effect of HIV-1 infection on human antibody repertoires and may have practical value for the discovery of therapeutic antibodies.     

The Organization of Song Repertoires in Song Sparrows: Themes and Variations

Although songtypes are generally considered to be important functional units in birdsong, they have not been well-characterized in terms of within- and between-songtype variation. We analyzed the song repertoires of 12 adult male song sparrows (Melospiza melodia) from a population in New York. We identified minimal units of production (MUPs), and calculated the degree to which song variants within a bird's repertoire were similar to each other based on an analysis of MUP sharing. Using statistical techniques from numerical taxonomy, we assessed how song variants clustered into songtypes, and we derived quantitative measures of within-songtype and between-songtype similarity. We found that birds produced a limited number of songtypes, but constantly produced new song variants within the framework of these types. Most song variants were produced very rarely; over 43% of song variants were produced only once. Repertoires differed in the degree of between-songtype similarity and in the number of songtypes defined, but there was comparatively less variation in within-songtype similarity. Between-songtype similarity and repertoire size were positively correlated. We argue that song sparrow songtypes are probabilistic units of song production, and discuss the functional and evolutionary implications of having vocal motor patterns organized in this way.     

Restricted utilization of germ-line VH genes and diversity of D regions in rabbit splenic Ig mRNA

Previous studies have suggested that the majority of rabbit germ-line VH genes encode molecules that are rarely found in serum or secretory Ig. To examine the repertoire of expressed VH genes, we prepared a cDNA library from splenic mRNA of an alpha 1/alpha 1 rabbit and isolated 10 complete VH-encoding cDNA clones. None of the cDNA clones hybridized to an oligomer that had hybridized to more than 50% of cloned germ-line VH genes. These data indicate that only a subset of germ-line VH genes are used in functional VDJ rearrangements. DNA sequence analysis demonstrated that the 10 cDNA clones contained highly similar VH regions, further suggesting that the repertoire of utilized VH genes is limited. In contrast, the D regions of each of the 10 clones exhibited little similarity to one another, suggesting that the rabbit has a large D region repertoire. We propose that the apparent lack of diversity within the VH segment of VDJ rearrangements is offset by extensive D region diversity.     

MHC gene copy number variation in Tasmanian devils: implications for the spread of a contagious cancer

Tasmanian devils face extinction owing to the emergence of a contagious cancer. Devil facial tumour disease (DFTD) is a clonal cancer spread owing to a lack of major histocompatibility complex (MHC) barriers in Tasmanian devil populations. We present a comprehensive screen of MHC diversity in devils and identify 25 MHC types and 53 novel sequences, but conclude that overall levels of MHC diversity at the sequence level are low. The majority of MHC Class I variation can be explained by allelic copy number variation with two to seven sequence variants identified per individual. MHC sequences are divided into two distinct groups based on sequence similarity. DFTD cells and most devils have sequences from both groups. Twenty per cent of individuals have a restricted MHC repertoire and contain only group I or only group II sequences. Counterintuitively, we postulate that the immune system of individuals with a restricted MHC repertoire may recognize foreign MHC antigens on the surface of the DFTD cell. The implication of these results for management of DFTD and this endangered species are discussed.     

An isozyme study of clone diversity and relative importance of sexual and vegetative recruitment in the grass Brachypodium pinnatum

We asked whether differences in abundance and seed prodtiction of Brachypodium pinnatum after 16 yr of four different experimental land use regimes were reflected in differences in vegetative dispersal distance, clone diversity, clone area, and the proportions of sexual and vegetative recruitment. Mean vegetative dispersal distance was 5.5 mm yr'. Electrophoresis of 5 polymorphic isozyme loci of 20 tillers sampled at defined positions in each of twelve 1 × 6 m sampling areas (3 per treatment) revealed considerable clonal diversity. Per sampling area we found on average 9.98 enzyme phenotypes (clones), mean Simpson index was 0.825. and mean Shannon index 0.801. The mean ratio of sexual vs vegetative recruitment was about 1:32000. Despite this low ratio, clonal diversity within the population of" B. pinnatum was higher than reported for other clonal plant populations, possibly because of its high ramet densities. Moan clone area was 5.73 m²-. i.e. mean clone radius was 1.35 m. None of the 10 pairwise correlations between abundance and seed production on the one hand, and number of clones per plot sample, plot Simpson index, plot Shannon index, ratio of vegetative vs sexual recruitment, and clone area on the other, was significant. Mean clone radius was 245 times larger than the mean distance of yearly vegetative dispersal which suggests old ages and low turnover rates of clones. The time scale of the inert response of clonal diversity of B. pinnatum to changes in land use appears to largely exceed the experimental period of 16 yr.     

Somatic hypermutation and selection of B cells in thymic germinal centers responding to acetylcholine receptor in myasthenia gravis

The muscle weakness in myasthenia gravis (MG) is mediated by autoantibodies against the nicotinic acetylcholine receptor (AChR) at the neuromuscular junction. Production of these pathogenic autoantibodies is believed to be associated with germinal centers (GC) and anti-AChR-secreting plasma cells in the hyperplastic thymus of patients with early onset MG (EOMG). Here, we describe the repertoire of rearranged heavy chain V genes and their clonal origins in GC from a typical EOMG patient. Three hundred fifteen rearranged Ig V(H) genes were amplified, cloned, and sequenced from sections of four thymic GC containing AChR-specific B cells. We found that thymic GC contain a remarkably heterogeneous population of B cells. Both naive and circulating memory B cells undergo Ag-driven clonal proliferation, somatic hypermutation, and selection. Numerous B cell clones were present, with no individual clone dominating the response. Comparisons of B cell clonal sequences from different GC and known anti-AChR Abs from other patients showed convergent mutations in the complementarity determining regions. These results are consistent with AChR driving an ongoing GC response in the thymus of EOMG patients. This is the first detailed analysis of B cell clones in human GC responding to a defined protein Ag, and the response we observed may reflect the effects of chronic stimulation by autoantigen.     

Size estimate of the alpha beta TCR repertoire of naive mouse splenocytes

The diversity of the T cell repertoire of mature T splenocytes is generated, in the thymus, by pairing of alpha and beta variable domains of the alpha beta TCR and by the rearrangements of various gene segments encoding these domains. In the periphery, it results from competition between various T cell subpopulations including recent thymic migrants and long-lived T cells. Quantitative data on the actual size of the T cell repertoire are lacking. Using PCR methods and extensive sequencing, we have measured for the first time the size of the TCR-alpha beta repertoire of naive mouse T splenocytes. There are 5-8 x 105 different nucleotide sequences of BV chains in the whole spleen of young adult mice. We have also determined the size of the BV repertoire in a subpopulation of AV2+ T splenocytes, which allows us to provide a minimum estimate of the alpha beta repertoire. We find that the mouse spleen harbors about 2 x 106 clones of about 10 cells each. This figure, although orders of magnitude smaller than the maximum theoretical diversity (estimated up to 1015), is still large enough to maintain a high functional diversity.