High-throughput sequencing of the zebrafish antibody repertoire

Massively Parallel Sequencing of Peritoneal and Splenic B Cell Repertoires Highlights Unique Properties of B-1 Cell Antibodies

B-1 cells are a unique subset of B cells that are positively selected for expressing autoreactive BCRs. We isolated RNA from peritoneal (B-1a, B-1b, B-2) and splenic (B-1a, marginal zone, follicular) B cells from C57BL/6 mice and used 5'-RACE to amplify the IgH V region using massively parallel sequencing. By analyzing 379,000 functional transcripts, we demonstrate that B-1a cells use a distinct and restricted repertoire. All B-1 cell subsets, especially peritoneal B-1a cells, had a high proportion of sequences without N additions, suggesting predominantly prenatal development. Their transcripts differed markedly and uniquely contained V_H11 and V_H12 genes, which were rearranged only with a restricted selection of D and J genes, unlike other V genes. Compared to peritoneal B-1a, the peritoneal B-1b repertoire was larger, had little overlap with B-1a, and most sequences contained N additions. Similarly, the splenic B-1a repertoire differed from peritoneal B-1a sequences, having more unique sequences and more frequent N additions, suggesting influx of B-1a cells into the spleen from nonperitoneal sites. Two CDR3s, previously described as Abs to bromelain-treated RBCs, comprised 43% of peritoneal B-1a sequences. We show that a single-chain variable fragment designed after the most prevalent B-1a sequence bound oxidation-specific epitopes such as the phosphocholine of oxidized phospholipids. In summary, we provide the IgH V region library of six murine B cell subsets, including, to our knowledge for the first time, a comparison between B-1a and B-1b cells, and we highlight qualities of B-1 cell Abs that indicate unique selection processes.

Herpesvirus Infection Induces both Specific and Heterologous Antiviral Antibodies in Carp

IgM antibody diversity induced by viral infection in teleost fish sera remains largely unexplored despite several studies performed on their transcript counterparts in lymphoid organs. Here, IgM binding to microarrays containing ~20,000 human proteins was used to study sera from carp (Cyprinus carpio) populations having high titers of viral neutralization in vitro after surviving an experimental infection with cyprinid herpes virus 3 (CyHV-3). The range of diversity of the induced antibodies was unexpectedly high, showing CyHV-3 infection-dependent, non-specific IgM-binding activity of a ~20-fold wider variety than that found in sera from healthy carp (natural antibodies) with no anti-CyHV-3 neutralization titers. An inverse correlation between the IgM-binding levels in healthy versus infection-survivor/healthy ratios suggests that an infection-dependent feed back-like mechanism may control such clonal expansion. Surprisingly, among the infection-expanded levels, not only specific anti-frgII_CyHV-3 and anti-CyHV-3 IgM-binding antibodies but also antibodies recognizing recombinant fragment epitopes from heterologous fish rhabdoviruses were detected in infection-survivor carp sera. Some alternative explanations for these findings in lower vertebrates are discussed.

Perspectives on theory at the interface of physics and biology

Theoretical physics is the search for simple and universal mathematical descriptions of the natural world. In contrast, much of modern biology is an exploration of the complexity and diversity of life. For many, this contrast is prima facie evidence that theory, in the sense that physicists use the word, is impossible in a biological context. For others, this contrast serves to highlight a grand challenge. I am an optimist, and believe (along with many colleagues) that the time is ripe for the emergence of a more unified theoretical physics of biological systems, building on successes in thinking about particular phenomena. In this essay I try to explain the reasons for my optimism, through a combination of historical and modern examples.

Distorted antibody repertoire developed in the absence of pre-B cell receptor formation

The pre-B cell receptor (pre-BCR), consisting of the μ heavy chain (μHC) and the surrogate light chain (SLC, Vpre-B and λ5), plays important roles during B cell development. The formation of the pre-BCR, which enables the nascent immunoglobulin HC to associate with the SLC, is considered a prerequisite for B cell development. However, a significant number of peripheral mature (leaky) B cells exist in SLC-deficient mice. These leaky B cells develop in the absence of pre-BCR and do not undergo the pre-BCR checkpoint. The antibody repertoires of leaky B cells thus reflect the absence of pre-BCR function. To investigate how the absence of the pre-BCR is circumvented by these leaky-B cells and examine the effect of the pre-BCR checkpoint on the antibody system, we analyzed the antibody repertoires of λ5-deficient (λ5^-/-) mice using next-generation sequencing. In λ5^-/- mice, spleen B cells displayed different patterns of VDJ-usage, relative to those in wild-type (WT) mice. Moreover, leaky B cells were neither derived from unusual B2 cells, characterized by particular LC gene rearrangements in the absence of pre-BCR signaling, nor from B1 cells, originating from different B cell progenitors. Analysis of the CDR-H3 amino acid sequences of μ-chain repertoires revealed that certain bone marrow B cells with particular CDR-H3 profiles undergo clonal expansion in λ5^-/- mice. Part of these CDR-H3s contain arginine(s) in the middle of the CDR-H3 loop in λ5^-/- mice, whereas few arginine(s) exist in this middle loop in WT CDR-H3s in the absence of clonal expansion. This CDR-H3 feature in λ5^-/- mice presumably reflects the role of the pre-BCR in autoantibody regulation, since arginine(s) are often found in the antigen-binding site of autoantibodies. Here, we present a unique viewpoint on the role of pre-BCR, by assessing the whole antibody repertoire formed in SLC-deficient mice.

Reproducibility and Reuse of Adaptive Immune Receptor Repertoire Data

High-throughput sequencing (HTS) of immunoglobulin (B-cell receptor, antibody) and T-cell receptor repertoires has increased dramatically since the technique was introduced in 2009 (1–3). This experimental approach explores the maturation of the adaptive immune system and its response to antigens, pathogens, and disease conditions in exquisite detail. It holds significant promise for diagnostic and therapy-guiding applications. New technology often spreads rapidly, sometimes more rapidly than the understanding of how to make the products of that technology reliable, reproducible, or usable by others. As complex technologies have developed, scientific communities have come together to adopt common standards, protocols, and policies for generating and sharing data sets, such as the MIAME protocols developed for microarray experiments. The Adaptive Immune Receptor Repertoire (AIRR) Community formed in 2015 to address similar issues for HTS data of immune repertoires. The purpose of this perspective is to provide an overview of the AIRR Community’s founding principles and present the progress that the AIRR Community has made in developing standards of practice and data sharing protocols. Finally, and most important, we invite all interested parties to join this effort to facilitate sharing and use of these powerful data sets (join@airr-community.org).

Integrating high-throughput screening and sequencing for monoclonal antibody discovery and engineering

Monoclonal antibody discovery and engineering is a field that has traditionally been dominated by high-throughput screening platforms (e.g. hybridomas and surface display). In recent years the emergence of high-throughput sequencing has made it possible to obtain large-scale information on antibody repertoire diversity. Additionally, it has now become more routine to perform high-throughput sequencing on antibody repertoires to also directly discover antibodies. In this review, we provide an overview of the progress in this field to date and show how high-throughput screening and sequencing are converging to deliver powerful new workflows for monoclonal antibody discovery and engineering.

Immune Repertoire after Immunization As Seen by Next-Generation Sequencing and Proteomics

The immune system produces a diverse repertoire of immunoglobulins in response to foreign antigens. During B-cell development, VDJ recombination and somatic mutations generate diversity, whereas selection processes remove it. Using both proteomic and NGS approaches, we characterized the immune repertoires in groups of rats after immunization with purified antigens. Proteomics and NGS data on the repertoire are in qualitative agreement, but did show quantitative differences that may relate to differences between the biological niches that were sampled for these approaches. Both methods contributed complementary information in the characterization of the immune repertoire. It was found that the immune repertoires resulting from each antigen had many similarities that allowed samples to cluster together, and that mutated immunoglobulin peptides were shared among animals with a response to the same antigen significantly more than for different antigens. However, the number of shared sequences decreased in a log-linear fashion relative to the number of animals that share them, which may affect future applications. A phylogenetic analysis on the NGS reads showed that reads from different individuals immunized with the same antigen populated distinct branches of the phylogram, an indication that the repertoire had converged. Also, similar mutation patterns were found in branches of the phylogenetic tree that were associated with antigen-specific immunoglobulins through proteomics data. Thus, data from different analysis methods and different experimental platforms show that the immunoglobulin repertoires of immunized animals have overlapping and converging features. With additional research, this may enable interesting applications in biotechnology and clinical diagnostics.

Quantitative Characterization of the T Cell Receptor Repertoire of Naïve and Memory Subsets Using an Integrated Experimental and Computational Pipeline Which Is Robust, Economical, and Versatile

The T cell receptor (TCR) repertoire can provide a personalized biomarker for infectious and non-infectious diseases. We describe a protocol for amplifying, sequencing, and analyzing TCRs which is robust, sensitive, and versatile. The key experimental step is ligation of a single-stranded oligonucleotide to the 3' end of the TCR cDNA. This allows amplification of all possible rearrangements using a single set of primers per locus. It also introduces a unique molecular identifier to label each starting cDNA molecule. This molecular identifier is used to correct for sequence errors and for effects of differential PCR amplification efficiency, thus producing more accurate measures of the true TCR frequency within the sample. This integrated experimental and computational pipeline is applied to the analysis of human memory and naive subpopulations, and results in consistent measures of diversity and inequality. After error correction, the distribution of TCR sequence abundance in all subpopulations followed a power law over a wide range of values. The power law exponent differed between naïve and memory populations, but was consistent between individuals. The integrated experimental and analysis pipeline we describe is appropriate to studies of T cell responses in a broad range of physiological and pathological contexts.

Reconstructing Antibody Repertoires from Error-Prone Immunosequencing Reads

Transforming error-prone immunosequencing datasets into Ab repertoires is a fundamental problem in immunogenomics, and a prerequisite for studies of immune responses. Although various repertoire reconstruction algorithms were released in the last 3 y, it remains unclear how to benchmark them and how to assess the accuracy of the reconstructed repertoires. We describe an accurate IgReC algorithm for constructing Ab repertoires from high-throughput immunosequencing datasets and a new framework for assessing the quality of reconstructed repertoires. Surprisingly, Ab repertoires constructed by IgReC from barcoded immunosequencing datasets in the blind mode (without using information about unique molecular identifiers) improved upon the repertoires constructed by the state-of-the-art tools that use barcoding. This finding suggests that IgReC may alleviate the need to generate repertoires using the barcoding technology (the workhorse of current immunogenomics efforts) because our computational approach to error correction of immunosequencing data is nearly as powerful as the experimental approach based on barcoding.

Describing the diversity of Ag specific receptors in vertebrates: Contribution of repertoire deep sequencing

During the last decades, gene and cDNA cloning identified TCR and Ig genes across vertebrates; genome sequencing of TCR and Ig loci in many species revealed the different organizations selected during evolution under the pressure of generating diverse repertoires of Ag receptors. By detecting clonotypes over a wide range of frequency, deep sequencing of Ig and TCR transcripts provides a new way to compare the structure of expressed repertoires in species of various sizes, at different stages of development, with different physiologies, and displaying multiple adaptations to the environment. In this review, we provide a short overview of the technologies currently used to produce global description of immune repertoires, describe how they have already been used in comparative immunology, and we discuss the future potential of such approaches. The development of these methodologies in new species holds promise for new discoveries concerning particular adaptations. As an example, understanding the development of adaptive immunity across metamorphosis in frogs has been made possible by such approaches. Repertoire sequencing is now widely used, not only in basic research but also in the context of immunotherapy and vaccination. Analysis of fish responses to pathogens and vaccines has already benefited from these methods. Finally, we also discuss potential advances based on repertoire sequencing of multigene families of immune sensors and effectors in invertebrates.

Omics in fish mucosal immunity

The mucosal immune system of fish is a complex network of immune cells and molecules that are constantly surveilling the environment and protecting the host from infection. A number of "omics" tools are now available and utilized to understand the complexity of mucosal immune systems in non-traditional animal models. This review summarizes recent advances in the implementation of "omics" tools pertaining to the four mucosa-associated lymphoid tissues in teleosts. Genomics, transcriptomics, proteomics, and "omics" in microbiome research require interdisciplinary collaboration and careful experimental design. The data-rich datasets generated are proving really useful at discovering new innate immune players in fish mucosal secretions, identifying novel markers of specific mucosal immune responses, unraveling the diversity of the B and T cell repertoires and characterizing the diversity of the microbial communities present in teleost mucosal surfaces. Bioinformatics, data analysis and storage platforms should be developed to facilitate rapid processing of large datasets, especially when mammalian tools such as bioinformatics analysis software are not available in fishes.

Accurate immune repertoire sequencing reveals malaria infection driven antibody lineage diversification in young children

Accurately measuring antibody repertoire sequence composition in a small amount of blood is challenging yet important for understanding repertoire responses to infection and vaccination. We develop molecular identifier clustering-based immune repertoire sequencing (MIDCIRS) and use it to study age-related antibody repertoire development and diversification before and during acute malaria in infants (< 12 months old) and toddlers (12-47 months old) with 4-8 ml of blood. Here, we show this accurate and high-coverage repertoire-sequencing method can use as few as 1000 naive B cells. Unexpectedly, we discover high levels of somatic hypermutation in infants as young as 3 months old. Antibody clonal lineage analysis reveals that somatic hypermutation levels are increased in both infants and toddlers upon infection, and memory B cells isolated from individuals who previously experienced malaria continue to induce somatic hypermutations upon malaria rechallenge. These results highlight the potential of antibody repertoire diversification in infants and toddlers.Somatic hypermutation of antibodies can occur in infants but are difficult to track. Here the authors present a new method called MIDCIRS for deep quantitative repertoire sequencing with few cells, and show infants as young as 3 months can expand antibody lineage complexity in response to malaria infection.

Systems immunology: just getting started

Systems-biology approaches in immunology take various forms, but here we review strategies for measuring a broad swath of immunological functions as a means of discovering previously unknown relationships and phenomena and as a powerful way of understanding the immune system as a whole. This approach has rejuvenated the field of vaccine development and has fostered hope that new ways will be found to combat infectious diseases that have proven refractory to classical approaches. Systems immunology also presents an important new strategy for understanding human immunity directly, taking advantage of the many ways the immune system of humans can be manipulated.

Structure-based cross-docking analysis of antibody-antigen interactions

Antibody–antigen interactions are critical to our immune response, and understanding the structure-based biophysical determinants for their binding specificity and affinity is of fundamental importance. We present a computational structure-based cross-docking study to test the identification of native antibody–antigen interaction pairs among cognate and non-cognate complexes. We picked a dataset of 17 antibody–antigen complexes of which 11 have both bound and unbound structures available, and we generated a representative ensemble of cognate and non-cognate complexes. Using the Rosetta interface score as a classifier, the cognate pair was the top-ranked model in 80% (14/17) of the antigen targets using bound monomer structures in docking, 35% (6/17) when using unbound, and 12% (2/17) when using the homology-modeled backbones to generate the complexes. Increasing rigid-body diversity of the models using RosettaDock’s local dock routine lowers the discrimination accuracy with the cognate antibody–antigen pair ranking in bound and unbound models but recovers additional top-ranked cognate complexes when using homology models. The study is the first structure-based cross-docking attempt aimed at distinguishing antibody–antigen binders from non-binders and demonstrates the challenges to address for the methods to be widely applicable to supplement high-throughput experimental antibody sequencing workflows.

Next-generation sequencing in chronic lymphocytic leukemia: recent findings and new horizons

The rapid progress in next-generation sequencing technologies has significantly contributed to our knowledge of the genetic events associated with the development, progression and treatment resistance of chronic lymphocytic leukemia patients. Together with the discovery of new driver mutations, next-generation sequencing has revealed an immense degree of both intra- and inter-tumor heterogeneity and enabled us to describe marked clonal evolution. Advances in immunogenetics may be implemented to detect minimal residual disease more sensitively and to track clonal B cell populations, their dynamics and molecular characteristics. The interpretation of these aspects is indispensable to thoroughly examine the genetic background of chronic lymphocytic leukemia. We review and discuss the recent results provided by the different next-generation sequencing techniques used in studying the chronic lymphocytic leukemia genome, as well as future perspectives in the methodologies and applications.

Molecular and genetic inflammation networks in major human diseases

It has been well-recognized that inflammation alongside tissue repair and damage maintaining tissue homeostasis determines the initiation and progression of complex diseases. Albeit with the accomplishment of having captured the most critical inflammation-involved molecules, genetic susceptibilities, epigenetic factors, and environmental factors, our schemata on the role of inflammation in complex diseases remain largely patchy, in part due to the success of reductionism in terms of research methodology per se. Omics data alongside the advances in data integration technologies have enabled reconstruction of molecular and genetic inflammation networks which shed light on the underlying pathophysiology of complex diseases or clinical conditions. Given the proven beneficial role of anti-inflammation in coronary heart disease as well as other complex diseases and immunotherapy as a revolutionary transition in oncology, it becomes timely to review our current understanding of the molecular and genetic inflammation networks underlying major human diseases. In this review, we first briefly discuss the complexity of infectious diseases and then highlight recently uncovered molecular and genetic inflammation networks in other major human diseases including obesity, type II diabetes, coronary heart disease, late onset Alzheimer's disease, Parkinson's disease, and sporadic cancer. The commonality and specificity of these molecular networks are addressed in the context of genetics based on genome-wide association study (GWAS). The double-sword role of inflammation, such as how the aberrant type 1 and/or type 2 immunity leads to chronic and severe clinical conditions, remains open in terms of the inflammasome and the core inflammatome network features. Increasingly available large Omics and clinical data in tandem with systems biology approaches have offered an exciting yet challenging opportunity toward reconstruction of more comprehensive and dynamic molecular and genetic inflammation networks, which hold great promise in transiting network snapshots to video-style multi-scale interplays of disease mechanisms, in turn leading to effective clinical intervention.

Ig-seq: Deep sequencing of the variable region of Atlantic salmon IgM heavy chain transcripts

Immunoglobulin M plays a key role in systemic protection of Atlantic salmon against pathogens. Until recent, studies have focused on antigen-specific antibodies and little is known about the IgM repertoire: its size, developmental changes and responses to antigens. We report the development of deep sequencing protocol to characterize the repertoire of IgM heavy chain variable region. Its structure and changes were examined at the early stages of life and after infection with virus of cardiac myopathy. Clonotypes are identified by the V and J gene segments and amino acid sequences of CDR3, which determine the contribution of the heavy chain to the antigen binding properties. A major fraction of transcripts are functional while the rest are either sterile (transcribed from noncoding parts of Ig loci) or include stop codons. Despite marked difference in frequencies of combinations of V and J genes, the size of repertoire is large. The IgM diversity steadily increases after hatch followed with temporal reduction during smoltification and recovery after seawater transfer. Most clonotypes are present only in one fish. However multiple transcripts in uninfected fish are produced exclusively from a small fraction of shared clonotypes. While only 4.7% of clonotypes are detected in three and more fish, they comprise 35% of transcripts. Increased frequencies of most abundant clonotypes were detected in the head kidney and blood at ten weeks after viral infection and all were shared. Occurrence of the same clonotypes in multiple individuals can be explained with either their simple structure or exposure to common antigens. Complexity of CDR3 assessed by contents of non complementary nucleotides is slightly lower in shared clonotypes but difference is small. High nucleotide diversity of CDR3 with identical amino acid sequences suggests selection.

Novel Approaches to Analyze Immunoglobulin Repertoires

Analysis of immunoglobulin (Ig) repertoires aims to comprehend Ig diversity with the goal of predicting humoral immune responses in the context of infection, vaccination, autoimmunity, and malignancies. The first next-generation sequencing (NGS) analyses of bulk B cell populations dramatically advanced sampling depth over previous low-throughput single-cell-based protocols, albeit at the expense of accuracy and loss of chain-pairing information. In recent years the field has substantially differentiated, with bulk analyses becoming more accurate while single-cell approaches have gained in throughput. Additionally, new platforms striving to combine high throughput and chain pairing have been developed as well as various computational tools for analysis. Here we review the developments of the past 4-5 years and discuss the open challenges.

Assessment of antibody library diversity through next generation sequencing and technical error compensation

Antibody libraries are important resources to derive antibodies to be used for a wide range of applications, from structural and functional studies to intracellular protein interference studies to developing new diagnostics and therapeutics. Whatever the goal, the key parameter for an antibody library is its complexity (also known as diversity), i.e. the number of distinct elements in the collection, which directly reflects the probability of finding in the library an antibody against a given antigen, of sufficiently high affinity. Quantitative evaluation of antibody library complexity and quality has been for a long time inadequately addressed, due to the high similarity and length of the sequences of the library. Complexity was usually inferred by the transformation efficiency and tested either by fingerprinting and/or sequencing of a few hundred random library elements. Inferring complexity from such a small sampling is, however, very rudimental and gives limited information about the real diversity, because complexity does not scale linearly with sample size. Next-generation sequencing (NGS) has opened new ways to tackle the antibody library complexity quality assessment. However, much remains to be done to fully exploit the potential of NGS for the quantitative analysis of antibody repertoires and to overcome current limitations. To obtain a more reliable antibody library complexity estimate here we show a new, PCR-free, NGS approach to sequence antibody libraries on Illumina platform, coupled to a new bioinformatic analysis and software (Diversity Estimator of Antibody Library, DEAL) that allows to reliably estimate the complexity, taking in consideration the sequencing error.

Deciphering antigen-responding antibody repertoires by using next-generation sequencing and confirming them through antibody-gene synthesis

Vast diversity and high specificity of antigen recognition by antibodies are hallmarks of the acquired immune system. Although the molecular mechanisms that yield the extremely large antibody repertoires are precisely understood, comprehensive description of the global antibody repertoire generated in individual bodies has been hindered by the lack of powerful measures. To obtain holistic understanding of the antibody-repertoire space, we used next-generation sequencing (NGS) to analyze the deep profiles of naive and antigen-responding repertoires of the IgM, IgG1, and IgG2c classes formed in individual mice. The overall landscapes of naive IgM repertoires were almost the same for each mouse, whereas those of IgG1 and IgG2c differed considerably among naive individuals. Next, we immunized mice with a model antigen, nitrophenol (NP)-hapten linked to chicken γ-globulin (CGG) carrier, and compared the antigen-responding repertoires in individual mice. To extract the complete antigen response, we developed an intelligible method for detecting common components of antigen-responding repertoires. The major responding antibodies were IGHV1-72/IGHD1-1/IGHJ2 for NP-hapten and IGHV9-3/IGHD3-1/IGHJ2 for CGG-carrier protein. The antigen-binding specificities of the identified antibodies were confirmed through ELISA after antibody-gene synthesis and expression of the corresponding NGS reads. Thus, we deciphered antigen-responding antibody repertoires by inclusively analyzing the antibody-repertoire space generated in individual bodies by using NGS, which avoided inadvertent omission of key antibody repertoires.

Immune engineering: from systems immunology to engineering immunity

The smallpox vaccine represents the earliest attempt in engineering immunity. The recent success of chimeric antigen receptor T cells (CAR-T cells) in cancer once again demonstrates the clinical potential of immune engineering. Inspired by this success, diverse approaches have been used to boost various aspects of immunity: engineering dendritic cells (DCs), natural killer (NK) cells, T cells, antibodies, cytokines, small peptides, and others. With recent development of various high-throughput technologies (of which engineers, especially biomedical engineers/bioengineers contributed significantly), such as immune repertoire sequencing, and analytical methods, a systems level of understanding immunity (or the lack of it) beyond model animals has provided critical insights into the human immune system. This review focuses on recent progressed made in systems biology and the engineering of adaptive immunity.

The role of adaptive immunity as an ecological filter on the gut microbiota in zebrafish

All animals live in intimate association with communities of microbes, collectively referred to as their microbiota. Certain host traits can influence which microbial taxa comprise the microbiota. One potentially important trait in vertebrate animals is the adaptive immune system, which has been hypothesized to act as an ecological filter, promoting the presence of some microbial taxa over others. Here we surveyed the intestinal microbiota of 68 wild-type zebrafish, with functional adaptive immunity, and 61 rag1^- zebrafish, lacking functional B- and T-cell receptors, to test the role of adaptive immunity as an ecological filter on the intestinal microbiota. In addition, we tested the robustness of adaptive immunity's filtering effects to host-host interaction by comparing the microbiota of fish populations segregated by genotype to those containing both genotypes. The presence of adaptive immunity individualized the gut microbiota and decreased the contributions of neutral processes to gut microbiota assembly. Although mixing genotypes led to increased phylogenetic diversity in each, there was no significant effect of adaptive immunity on gut microbiota composition in either housing condition. Interestingly, the most robust effect on microbiota composition was co-housing within a tank. In all, these results suggest that adaptive immunity has a role as an ecological filter of the zebrafish gut microbiota, but it can be overwhelmed by other factors, including transmission of microbes among hosts.

The expressed TCRβ CDR3 repertoire is dominated by conserved DNA sequences in channel catfish

We analyzed by high-throughput sequencing T cell receptor beta CDR3 repertoires expressed by αβ T cells in outbred channel catfish before and after an immunizing infection with the parasitic protozoan Ichthyophthirius multifiliis. We compared CDR3 repertoires in caudal fin before infection and at three weeks after infection, and in skin, PBL, spleen and head kidney at seven and twenty-one weeks after infection. Public clonotypes with the same CDR3 amino acid sequence were expressed by αβ T cells that underwent clonal expansion following development of immunity. These clonally expanded αβ T cells were primarily located in spleen and skin, which is a site of infection. Although multiple DNA sequences were expected to code for each public clonotype, each public clonotype was predominately coded by an identical CDR3 DNA sequence in combination with the same J gene in all fish. The processes underlying this shared use of CDR3 DNA sequences are not clear.

Insights into immune system development and function from mouse T-cell repertoires

The ability of the adaptive immune system to respond to arbitrary pathogens stems from the broad diversity of immune cell surface receptors. This diversity originates in a stochastic DNA editing process (VDJ recombination) that acts on the surface receptor gene each time a new immune cell is created from a stem cell. By analyzing T-cell receptor (TCR) sequence repertoires taken from the blood and thymus of mice of different ages, we quantify the changes in the VDJ recombination process that occur from embryo to young adult. We find a rapid increase with age in the number of random insertions and a dramatic increase in diversity. Because the blood accumulates thymic output over time, blood repertoires are mixtures of different statistical recombination processes, and we unravel the mixture statistics to obtain a picture of the time evolution of the early immune system. Sequence repertoire analysis also allows us to detect the statistical impact of selection on the output of the VDJ recombination process. The effects we find are nearly identical between thymus and blood, suggesting that our analysis mainly detects selection for proper folding of the TCR receptor protein. We further find that selection is weaker in laboratory mice than in humans and it does not affect the diversity of the repertoire.

Hierarchical Clustering Can Identify B Cell Clones with High Confidence in Ig Repertoire Sequencing Data

Adaptive immunity is driven by the expansion, somatic hypermutation, and selection of B cell clones. Each clone is the progeny of a single B cell responding to Ag, with diversified Ig receptors. These receptors can now be profiled on a large scale by next-generation sequencing. Such data provide a window into the microevolutionary dynamics that drive successful immune responses and the dysregulation that occurs with aging or disease. Clonal relationships are not directly measured, but they must be computationally inferred from these sequencing data. Although several hierarchical clustering-based methods have been proposed, they vary in distance and linkage methods and have not yet been rigorously compared. In this study, we use a combination of human experimental and simulated data to characterize the performance of hierarchical clustering-based methods for partitioning sequences into clones. We find that single linkage clustering has high performance, with specificity, sensitivity, and positive predictive value all >99%, whereas other linkages result in a significant loss of sensitivity. Surprisingly, distance metrics that incorporate the biases of somatic hypermutation do not outperform simple Hamming distance. Although errors were more likely in sequences with short junctions, using the entire dataset to choose a single distance threshold for clustering is near optimal. Our results suggest that hierarchical clustering using single linkage with Hamming distance identifies clones with high confidence and provides a fully automated method for clonal grouping. The performance estimates we develop provide important context to interpret clonal analysis of repertoire sequencing data and allow for rigorous testing of other clonal grouping algorithms.

SpotLight Proteomics: uncovering the hidden blood proteome improves diagnostic power of proteomics

The human blood proteome is frequently assessed by protein abundance profiling using a combination of liquid chromatography and tandem mass spectrometry (LC-MS/MS). In traditional sequence database search, many good-quality MS/MS data remain unassigned. Here we uncover the hidden part of the blood proteome via novel SpotLight approach. This method combines de novo MS/MS sequencing of enriched antibodies and co-extracted proteins with subsequent label-free quantification of new and known peptides in both enriched and unfractionated samples. In a pilot study on differentiating early stages of Alzheimer’s disease (AD) from Dementia with Lewy Bodies (DLB), on peptide level the hidden proteome contributed almost as much information to patient stratification as the apparent proteome. Intriguingly, many of the new peptide sequences are attributable to antibody variable regions, and are potentially indicative of disease etiology. When the hidden and apparent proteomes are combined, the accuracy of differentiating AD (n = 97) and DLB (n = 47) increased from ≈85% to ≈95%. The low added burden of SpotLight proteome analysis makes it attractive for use in clinical settings.

Single amino acid fingerprinting of the human antibody repertoire with high density peptide arrays

The antibody species that patrol in a patient's blood are an invaluable part of the immune system. While most of them shield us from life-threatening infections, some of them do harm in autoimmune diseases. If we knew exactly all the antigens that elicited all the antibody species within a group of patients, we could learn which ones correlate with immune protection, are irrelevant, or do harm. Here, we demonstrate an approach to this question: First, we use a plethora of phage-displayed peptides to identify many different serum antibody binding peptides. Next, we synthesize identified peptides in the array format and rescreen the serum used for phage panning to validate antibody binding peptides. Finally, we systematically vary the sequence of validated antibody binding peptides to identify those amino acids within the peptides that are crucial for binding "their" antibody species. The resulting immune fingerprints can then be used to trace them back to potential antigens. We investigated the serum of an individual in this pipeline, which led to the identification of 73 antibody fingerprints. Some fingerprints could be traced back to their most likely antigen, for example the immunodominant capsid protein VP1 of enteroviruses, most likely elicited by the ubiquitous poliovirus vaccination. Thus, with our approach, it is possible, to pinpoint those antibody species that correlate with a certain antigen, without any pre-information. This can help to unravel hitherto enigmatic diseases.

Transcriptomic analysis of common carp anterior kidney during Cyprinid herpesvirus 3 infection: Immunoglobulin repertoire and homologue functional divergence

Cyprinid herpesvirus 3 (CyHV-3) infects koi and common carp and causes widespread mortalities. While the virus is a significant concern for aquaculture operations in many countries, in Australia the virus may be a useful biocontrol agent for pest carp. However, carp immune responses to CyHV-3, and the molecular mechanisms underpinning resistance, are not well understood. Here we used RNA-Seq on carp during different phases of CyHV-3 infection to detect the gene expression dynamics of both host and virus simultaneously. During acute CyHV-3 infection, the carp host modified the expression of genes involved in various immune systems and detoxification pathways. Moreover, the activated pathways were skewed toward humoral immune responses, which may have been influenced by the virus itself. Many immune-related genes were duplicated in the carp genome, and often these were expressed differently across the infection phases. Of particular interest were two interleukin-10 homologues that were not expressed synchronously, suggesting neo- or sub-functionalization. The carp immunoglobulin repertoire significantly diversified during active CyHV-3 infection, which was followed by the selection of high-affinity B-cells. This is indicative of a developing adaptive immune response, and is the first attempt to use RNA-Seq to understand this process in fish during a viral infection.

Characterization of the diversity of T cell receptor γδ complementary determinant region 3 in human peripheral blood by Immune Repertoire Sequencing

γδ T cells function as sentinels in early host response to infections and malignancies. Although γδ T cells are regarded as innate immune cells and recognize antigens in a non-MHC restricted manner, they possess a huge diversity of complementary determinant region 3 (CDR3) of T cell receptor (TCR) generated by the rearrangement of germ-line gene V- (D) -J-C fragments. However, the detailed characteristics of the TCRγδ CDR3 repertoire remain unclear. A comprehensive analysis would answer fundamental questions about the diversity of the TCRγδ CDR3 repertoire and elucidate the mechanism underlying γδ T cell recognition of pathogens and tumor antigens. In this study, we used Immune Repertoire Sequencing (IR-SEQ) to analyze the diversity of TCRγδ CDR3 repertoires from 30 healthy donors. The results show that IR-SEQ had sufficient repeatability to analyze the TCRγδ CDR3 repertoire. The diversity of TCRγδ CDR3 repertoire is quite dispersed and individually different. The TCR δ chain (TRD) repertoire displayed more diversity and less sharing among individuals compared with TCR γ chain (TRG). To our knowledge, this is the first study to use IR-SEQ to characterize the repertoire of TCRγδ CDR3 in human peripheral blood γδ T cells by using IR-SEQ. Our findings provide a basic understanding of the diversity of TCRγδ repertoire in the physiological condition, which provides a clue to the underlying mechanism of γδ T cell recognition of pathogens and tumor antigens.

Comparative analysis of the feline immunoglobulin repertoire

Next-Generation Sequencing combined with bioinformatics is a powerful tool for analyzing the large number of DNA sequences present in the expressed antibody repertoire and these data sets can be used to advance a number of research areas including antibody discovery and engineering. The accurate measurement of the immune repertoire sequence composition, diversity and abundance is important for understanding the repertoire response in infections, vaccinations and cancer immunology and could also be useful for elucidating novel molecular targets. In this study 4 individual domestic cats (Felis catus) were subjected to antibody repertoire sequencing with total number of sequences generated 1079863 for VH for IgG, 1050824 VH for IgM, 569518 for VK and 450195 for VL. Our analysis suggests that a similar VDJ expression patterns exists across all cats. Similar to the canine repertoire, the feline repertoire is dominated by a single subgroup, namely VH3. The antibody paratope of felines showed similar amino acid variation when compared to human, mouse and canine counterparts. All animals show a similarly skewed VH CDR-H3 profile and, when compared to canine, human and mouse, distinct differences are observed. Our study represents the first attempt to characterize sequence diversity in the expressed feline antibody repertoire and this demonstrates the utility of using NGS to elucidate entire antibody repertoires from individual animals. These data provide significant insight into understanding the feline immune system function.

B cell repertoires in HLA-sensitized kidney transplant candidates undergoing desensitization therapy

Background

Kidney transplantation is the most effective treatment for end-stage renal disease. Sensitization refers to pre-existing antibodies against human leukocyte antigen (HLA) protein and remains a major barrier to successful transplantation. Despite implementation of desensitization strategies, many candidates fail to respond. Our objective was to determine whether measuring B cell repertoires could differentiate candidates that respond to desensitization therapy.

Methods

We developed an assay based on high-throughput DNA sequencing of the variable domain of the heavy chain of immunoglobulin genes to measure changes in B cell repertoires in 19 highly HLA-sensitized kidney transplant candidates undergoing desensitization and 7 controls with low to moderate HLA sensitization levels. Responders to desensitization had a decrease of 5% points or greater in cumulated calculated panel reactive antibody (cPRA) levels, and non-responders had no decrease in cPRA.

Results

Dominant B cell clones were not observed in highly sensitized candidates, suggesting that the B cells responsible for sensitization are either not present in peripheral blood or present at comparable levels to other circulating B cells. Candidates that responded to desensitization therapy had pre-treatment repertoires composed of a larger fraction of class-switched (IgG and IgA) isotypes compared to non-responding candidates. After B cell depleting therapy, the proportion of switched isotypes increased and the mutation frequencies of the remaining non-switched isotypes (IgM and IgD) increased in both responders and non-responders, perhaps representing a shift in the repertoire towards memory B cells or plasmablasts. Conversely, after transplantation, non-switched isotypes with fewer mutations increased, suggesting a shift in the repertoire towards naïve B cells.

Conclusions

Relative abundance of different B cell isotypes is strongly perturbed by desensitization therapy and transplantation, potentially reflecting changes in the relative abundance of memory and naïve B cell compartments. Candidates that responded to therapy experienced similar changes to those that did not respond. Further studies are required to understand differences between these two groups of highly sensitized kidney transplant candidates.

Electronic supplementary material

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

Direct measurement of B-cell receptor repertoire's composition and variation in systemic lupus erythematosus

Systemic lupus erythematosus (SLE) is a systemic autoimmune disease that is known to be associated with polyclonal B-cell hyper-reactivity. B-cell receptor (BCR) has a central role in B-cell development, activation, survival and apoptosis, and thus is a critical component of the regulation of both protective and autoreactive B cells. In this study, we applied multiplex PCR and Illumina high-throughput sequencing to study the composition and variation of the BCRs in peripheral blood mononuclear cells from SLE patients and healthy donors (NC). We found that SLE group displayed significantly shorter CDR3 average length (14.86±0.76aa vs 15.70±0.43aa), more arginine percentage of CDR3 amino acids (7.57±0.20% vs 7.32±0.19%) and poorer immunological diversity than the healthy ones. CDR3 sequence YGMDV present in all SLE samples may provide more information in generating more effective B-cell targeted diagnosis/therapies strategies.

Can the Immune System Perform a t-Test?

The self-nonself discrimination hypothesis remains a landmark concept in immunology. It proposes that tolerance breaks down in the presence of nonself antigens. In strike contrast, in statistics, occurrence of nonself elements in a sample (i.e., outliers) is not obligatory to violate the null hypothesis. Very often, what is crucial is the combination of (self) elements in a sample. The two views on how to detect a change seem challengingly different and it could seem difficult to conceive how immunological cellular interactions could trigger responses with a precision comparable to some statistical tests. Here it is shown that frustrated cellular interactions reconcile the two views within a plausible immunological setting. It is proposed that the adaptive immune system can be promptly activated either when nonself ligands are detected or self-ligands occur in abnormal combinations. In particular we show that cellular populations behaving in this way could perform location statistical tests, with performances comparable to t or KS tests, or even more general data mining tests such as support vector machines or random forests. In more general terms, this work claims that plausible immunological models should provide accurate detection mechanisms for host protection and, furthermore, that investigation on mechanisms leading to improved detection in “in silico” models can help unveil how the real immune system works.

Reproducibility and Reuse of Adaptive Immune Receptor Repertoire Data

High-throughput sequencing (HTS) of immunoglobulin (B-cell receptor, antibody) and T-cell receptor repertoires has increased dramatically since the technique was introduced in 2009 (1-3). This experimental approach explores the maturation of the adaptive immune system and its response to antigens, pathogens, and disease conditions in exquisite detail. It holds significant promise for diagnostic and therapy-guiding applications. New technology often spreads rapidly, sometimes more rapidly than the understanding of how to make the products of that technology reliable, reproducible, or usable by others. As complex technologies have developed, scientific communities have come together to adopt common standards, protocols, and policies for generating and sharing data sets, such as the MIAME protocols developed for microarray experiments. The Adaptive Immune Receptor Repertoire (AIRR) Community formed in 2015 to address similar issues for HTS data of immune repertoires. The purpose of this perspective is to provide an overview of the AIRR Community's founding principles and present the progress that the AIRR Community has made in developing standards of practice and data sharing protocols. Finally, and most important, we invite all interested parties to join this effort to facilitate sharing and use of these powerful data sets (join@airr-community.org).

Measuring the sequence-affinity landscape of antibodies with massively parallel titration curves

Despite the central role that antibodies play in the adaptive immune system and in biotechnology, much remains unknown about the quantitative relationship between an antibody’s amino acid sequence and its antigen binding affinity. Here we describe a new experimental approach, called Tite-Seq, that is capable of measuring binding titration curves and corresponding affinities for thousands of variant antibodies in parallel. The measurement of titration curves eliminates the confounding effects of antibody expression and stability that arise in standard deep mutational scanning assays. We demonstrate Tite-Seq on the CDR1H and CDR3H regions of a well-studied scFv antibody. Our data shed light on the structural basis for antigen binding affinity and suggests a role for secondary CDR loops in establishing antibody stability. Tite-Seq fills a large gap in the ability to measure critical aspects of the adaptive immune system, and can be readily used for studying sequence-affinity landscapes in other protein systems.

DOI:http://dx.doi.org/10.7554/eLife.23156.001

Antibodies are proteins produced by cells of the immune system to tag or neutralize potential threats to the body, such as foreign substances and disease-causing microbes. Antibodies do this by binding to target molecules called antigens. An antibody’s ability to bind to an antigen depends on the sequence of amino acids – the building blocks of proteins – that make up the antibody. Through a process that randomizes this sequence of amino acids, the immune system generates a vast pool of antibodies that are able to target almost any foreign antigen that exists in nature.

Currently, little is understood about how the sequence of amino acids in an antibody determines how strongly that antibody binds to its antigen target – a property referred to as the antibody’s binding affinity. Answering this fundamental question requires techniques that can measure the affinities of many different antibodies at the same time. However, previous high-throughput methods have been unable to provide quantitative measurements of binding affinities. These kinds of measurements are difficult because an antibody’s amino acid sequence governs more than just binding affinity: it also affects how easy it is to produce that antibody, and what fraction of antibody molecules work properly.

Adams et al. now describe a new method, named “Tite-Seq”, that overcomes these issues. First, thousands of different antibodies are displayed on the surface of yeast cells, with each cell carrying a single kind of antibody. These cells are then incubated with fluorescently labeled antigen at a wide range of different concentrations. Next, the yeast cells are sorted based on how brightly they glow; brighter cells have more antigen bound to them, and so it is possible to calculate how much of the antigen is bound to each kind of antibody at each concentration. Plotting these data provides a “binding curve” for each antibody, which is then used to read off the antibody’s binding affinity in a way that is not affected by the factors that have plagued other high-throughput methods.

Tite-Seq is thus able to measure the binding affinities for thousands of different antibodies at the same time. This will potentially allow researchers to address many fundamental and yet unanswered questions about how the immune system works. Tite-Seq can also be used to measure how amino acid sequence affects the binding affinity of proteins other than antibodies.

DOI:http://dx.doi.org/10.7554/eLife.23156.002

T Cell Repertoire Diversity Is Decreased in Type 1 Diabetes Patients

Type 1 diabetes mellitus (T1D) is an immune-mediated disease. The autoreactive T cells in T1D patients attack and destroy their own pancreatic cells. In order to systematically investigate the potential autoreactive T cell receptors (TCRs), we used a high-throughput immune repertoire sequencing technique to profile the spectrum of TCRs in individual T1D patients and controls. We sequenced the T cell repertoire of nine T1D patients, four type 2 diabetes (T2D) patients, and six nondiabetic controls. The diversity of the T cell repertoire in T1D patients was significantly decreased in comparison with T2D patients (P=7.0E-08 for CD4⁺ T cells, P=1.4E-04 for CD8⁺ T cells) and nondiabetic controls (P=2.7E-09 for CD4⁺ T cells, P=7.6E-06 for CD8⁺ T cells). Moreover, T1D patients had significantly more highly-expanded T cell clones than T2D patients (P=5.2E-06 for CD4⁺ T cells, P=1.9E-07 for CD8⁺ T cells) and nondiabetic controls (P=1.7E-07 for CD4⁺ T cells, P=3.3E-03 for CD8⁺ T cells). Furthermore, we identified a group of highly-expanded T cell receptor clones that are shared by more than two T1D patients. Although further validation in larger cohorts is needed, our data suggest that T cell receptor diversity measurements may become a valuable tool in investigating diabetes, such as using the diversity as an index to distinguish different types of diabetes.

Single-cell barcoding and sequencing using droplet microfluidics

Single-cell RNA sequencing has recently emerged as a powerful tool for mapping cellular heterogeneity in diseased and healthy tissues, yet high-throughput methods are needed for capturing the unbiased diversity of cells. Droplet microfluidics is among the most promising candidates for capturing and processing thousands of individual cells for whole-transcriptome or genomic analysis in a massively parallel manner with minimal reagent use. We recently established a method called inDrops, which has the capability to index >15,000 cells in an hour. A suspension of cells is first encapsulated into nanoliter droplets with hydrogel beads (HBs) bearing barcoding DNA primers. Cells are then lysed and mRNA is barcoded (indexed) by a reverse transcription (RT) reaction. Here we provide details for (i) establishing an inDrops platform (1 d); (ii) performing hydrogel bead synthesis (4 d); (iii) encapsulating and barcoding cells (1 d); and (iv) RNA-seq library preparation (2 d). inDrops is a robust and scalable platform, and it is unique in its ability to capture and profile >75% of cells in even very small samples, on a scale of thousands or tens of thousands of cells.

An elaborate landscape of the human antibody repertoire against enterovirus 71 infection is revealed by phage display screening and deep sequencing

Enterovirus 71 (EV71) causes outbreaks of hand, foot and mouth disease (HFMD), primarily in the Asia-Pacific area, that are often associated with complications of severe to fatal neurological symptoms. There are currently no anti-viral therapies or vaccines available for the treatment of EV71 infection. Illustrating human antibody responses neutralizing EV71 infection could potentially provide valuable information for the development of effective therapies and vaccines. Here, we constructed a comprehensive phage display library based on peripheral blood of eight EV71-infected donors and identified 27 EV71-specific human antibodies, of which four have neutralizing activity in in vitro experiments. Deep sequencing analysis of the antibody heavy chains at the transcript level of another three independent EV71-infected donors and three controls demonstrates that heavy chains of the EV71-specific antibodies are conserved among EV71-infected individuals but absent in controls, suggesting convergent evolution of human antibodies against EV71.

Accurate Structure Prediction of CDR H3 Loops Enabled by a Novel Structure-Based C-Terminal Constraint

Ab structure prediction has made great strides, but accurately modeling CDR H3 loops remains elusive. Unlike the other five CDR loops, CDR H3 does not adopt canonical conformations and must be modeled de novo. During Antibody Modeling Assessment II, we found that biasing simulations toward kinked conformations enables generating low-root mean square deviation models (Weitzner et al. 2014. Proteins 82: 1611-1623), and since then, we have presented new geometric parameters defining the kink conformation (Weitzner et al. 2015. Structure 23: 302-311). In this study, we use these parameters to develop a new biasing constraint. When applied to a benchmark set of high-quality CDR H3 loops, the average minimum root mean square deviation sampled is 0.93 Å, compared with 1.34 Å without the constraint. We then test the performance of the constrained de novo method for homology modeling and rigid-body docking and present the results for 1) the Antibody Modeling Assessment II targets, 2) the 2009 RosettaAntibody benchmark set, and 3) the high-quality set.

Advanced Methodologies in High-Throughput Sequencing of Immune Repertoires

In recent years, major efforts have been made to develop sophisticated experimental and bioinformatic workflows for sequencing adaptive immune repertoires. The immunological insight gained has been applied to fields as varied as lymphocyte biology, immunodiagnostics, vaccines, cancer immunotherapy, and antibody engineering. In this review, we provide a detailed overview of these advanced methodologies, focusing specifically on strategies to reduce sequencing errors and bias and to achieve high-throughput pairing of variable regions (e.g., heavy-light or alpha-beta chains). In addition, we highlight recent technologies for single-cell transcriptome sequencing that can be integrated with immune repertoires. Finally, we provide a perspective on advanced immune repertoire sequencing and its ability to impact basic immunology, biopharmaceutical drug discovery and development, and cancer immunotherapy.

Minimal Residual Disease Detection and Evolved IGH Clones Analysis in Acute B Lymphoblastic Leukemia Using IGH Deep Sequencing

Acute B lymphoblastic leukemia (B-ALL) is one of the most common types of childhood cancer worldwide and chemotherapy is the main treatment approach. Despite good response rates to chemotherapy regiments, many patients eventually relapse and minimal residual disease (MRD) is the leading risk factor for relapse. The evolution of leukemic clones during disease development and treatment may have clinical significance. In this study, we performed immunoglobulin heavy chain (IGH) repertoire high throughput sequencing (HTS) on the diagnostic and post-treatment samples of 51 pediatric B-ALL patients. We identified leukemic IGH clones in 92.2% of the diagnostic samples and nearly half of the patients were polyclonal. About one-third of the leukemic clones have correct open reading frame in the complementarity determining region 3 (CDR3) of IGH, which demonstrates that the leukemic B cells were in the early developmental stage. We also demonstrated the higher sensitivity of HTS in MRD detection and investigated the clinical value of using peripheral blood in MRD detection and monitoring the clonal IGH evolution. In addition, we found leukemic clones were extensively undergoing continuous clonal IGH evolution by variable gene replacement. Dynamic frequency change and newly emerged evolved IGH clones were identified upon the pressure of chemotherapy. In summary, we confirmed the high sensitivity and universal applicability of HTS in MRD detection. We also reported the ubiquitous evolved IGH clones in B-ALL samples and their response to chemotherapy during treatment.

Comparative Analysis of Immune Repertoires between Bactrian Camel's Conventional and Heavy-Chain Antibodies

Compared to classical antibodies, camel heavy chain antibodies (HCAbs) are smaller in size due to lack of the light chain and the first constant domain of the heavy chain (CH1 region). The variable regions of HCAbs (VHHs) are more soluble and stable than that of conventional antibodies (VHs). Even with such simple structure, they are still functional in antigen binding. Although HCAbs have been extensively investigated over the past two decades, most efforts have been based upon low throughput sequence analysis, and there are only limited reports trying to analyze and describe the complete immune repertoire (IR) of camel HCAbs. Here we leveraged the high-throughput data generated by Next Generation Sequencing (NGS) of the variable domains of the antibody heavy chains from three Bactrian camels to conduct in-depth comparative analyses of the immunoglobulin repertoire. These include analyses of the complementary determining region 3 (CDR3) length and distribution, mutation rate, antibody characteristic amino acids, the distribution of the cysteine (Cys) codons, and the non-classical VHHs. We found that there is higher diversity in the CDR2 than in the other sub-regions, and there is a higher mutation rate in the VHHs than in the VHs (P < 0.05). In addition to substitutions at amino acid (AA) residue positions NO.49/50/52 between VH and VHH clones, we also observed other substitutions at the positions NO.40/54/57/96/101 that could lead to additional structural alterations. We also found that VH-derived VHH clones, referred to as non-classical VHH clones in this study, accounted for about 8% of all clones. Further, only 5%-10% clones had the Trp > Arg AA substitution at the first position of framework 4 for all types of clones. We present, for the first time, a relatively complete picture of the Bactrian camel antibody immune repertoire, including conventional antibody (Ab) and HCAbs, using PCR and in silico analysis based on high-throughput NGS data.

Deep sequencing in library selection projects: what insight does it bring?

High throughput sequencing is poised to change all aspects of the way antibodies and other binders are discovered and engineered. Millions of available sequence reads provide an unprecedented sampling depth able to guide the design and construction of effective, high quality naïve libraries containing tens of billions of unique molecules. Furthermore, during selections, high throughput sequencing enables quantitative tracing of enriched clones and position-specific guidance to amino acid variation under positive selection during antibody engineering. Successful application of the technologies relies on specific PCR reagent design, correct sequencing platform selection, and effective use of computational tools and statistical measures to remove error, identify antibodies, estimate diversity, and extract signatures of selection from the clone down to individual structural positions. Here we review these considerations and discuss some of the remaining challenges to the widespread adoption of the technology.

High-Throughput Sequencing-Based Immune Repertoire Study during Infectious Disease

The selectivity of the adaptive immune response is based on the enormous diversity of T and B cell antigen-specific receptors. The immune repertoire, the collection of T and B cells with functional diversity in the circulatory system at any given time, is dynamic and reflects the essence of immune selectivity. In this article, we review the recent advances in immune repertoire study of infectious diseases, which were achieved by traditional techniques and high-throughput sequencing (HTS) techniques. HTS techniques enable the determination of complementary regions of lymphocyte receptors with unprecedented efficiency and scale. This progress in methodology enhances the understanding of immunologic changes during pathogen challenge and also provides a basis for further development of novel diagnostic markers, immunotherapies, and vaccines.