High-throughput sequencing of the zebrafish antibody repertoire

Interwoven processes in fish development: microbial community succession and immune maturation

Fishes are hosts for many microorganisms that provide them with beneficial effects on growth, immune system development, nutrition and protection against pathogens. In order to avoid spreading of infectious diseases in aquaculture, prevention includes vaccinations and routine disinfection of eggs and equipment, while curative treatments consist in the administration of antibiotics. Vaccination processes can stress the fish and require substantial farmer's investment. Additionally, disinfection and antibiotics are not specific, and while they may be effective in the short term, they have major drawbacks in the long term. Indeed, they eliminate beneficial bacteria which are useful for the host and promote the raising of antibiotic resistance in beneficial, commensal but also in pathogenic bacterial strains. Numerous publications highlight the importance that plays the diversified microbial community colonizing fish (i.e., microbiota) in the development, health and ultimately survival of their host. This review targets the current knowledge on the bidirectional communication between the microbiota and the fish immune system during fish development. It explores the extent of this mutualistic relationship: on one hand, the effect that microbes exert on the immune system ontogeny of fishes, and on the other hand, the impact of critical steps in immune system development on the microbial recruitment and succession throughout their life. We will first describe the immune system and its ontogeny and gene expression steps in the immune system development of fishes. Secondly, the plurality of the microbiotas (depending on host organism, organ, and development stage) will be reviewed. Then, a description of the constant interactions between microbiota and immune system throughout the fish's life stages will be discussed. Healthy microbiotas allow immune system maturation and modulation of inflammation, both of which contribute to immune homeostasis. Thus, immune equilibrium is closely linked to microbiota stability and to the stages of microbial community succession during the host development. We will provide examples from several fish species and describe more extensively the mechanisms occurring in zebrafish model because immune system ontogeny is much more finely described for this species, thanks to the many existing zebrafish mutants which allow more precise investigations. We will conclude on how the conceptual framework associated to the research on the immune system will benefit from considering the relations between microbiota and immune system maturation. More precisely, the development of active tolerance of the microbiota from the earliest stages of life enables the sustainable establishment of a complex healthy microbial community in the adult host. Establishing a balanced host-microbiota interaction avoids triggering deleterious inflammation, and maintains immunological and microbiological homeostasis.

Single-cell RNA sequencing unveils the hidden powers of zebrafish kidney for generating both hematopoiesis and adaptive antiviral immunity

The vertebrate kidneys play two evolutionary conserved roles in waste excretion and osmoregulation. Besides, the kidney of fish is considered as a functional ortholog of mammalian bone marrow that serves as a hematopoietic hub for generating blood cell lineages and immunological responses. However, knowledge about the properties of kidney hematopoietic cells, and the functionality of the kidney in fish immune systems remains to be elucidated. To this end, our present study generated a comprehensive atlas with 59 hematopoietic stem/progenitor cell (HSPC) and immune-cells types from zebrafish kidneys via single-cell transcriptome profiling analysis. These populations included almost all known cells associated with innate and adaptive immunity, and displayed differential responses to viral infection, indicating their diverse functional roles in antiviral immunity. Remarkably, HSPCs were found to have extensive reactivities to viral infection, and the trained immunity can be effectively induced in certain HSPCs. In addition, the antigen-stimulated adaptive immunity can be fully generated in the kidney, suggesting the kidney acts as a secondary lymphoid organ. These results indicated that the fish kidney is a dual-functional entity with functionalities of both primary and secondary lymphoid organs. Our findings illustrated the unique features of fish immune systems, and highlighted the multifaced biology of kidneys in ancient vertebrates.

Immune aging in annual killifish

Turquoise killifish (Nothobranchius furzeri) evolved a naturally short lifespan of about six months and exhibit aging hallmarks that affect multiple organs. These hallmarks include protein aggregation, telomere shortening, cellular senescence, and systemic inflammation. Turquoise killifish possess the full spectrum of vertebrate-specific innate and adaptive immune system. However, during their recent evolutionary history, they lost subsets of mucosal-specific antibody isoforms that are present in other teleosts. As they age, the immune system of turquoise killifish undergoes dramatic cellular and systemic changes. These changes involve increased inflammation, reduced antibody diversity, an increased prevalence of pathogenic microbes in the intestine, and extensive DNA damage in immune progenitor cell clusters. Collectively, the wide array of age-related changes occurring in turquoise killifish suggest that, despite an evolutionary separation spanning hundreds of millions of years, teleosts and mammals share common features of immune system aging. Hence, the spontaneous aging observed in the killifish immune system offers an excellent opportunity for discovering fundamental and conserved aspects associated with immune system aging across vertebrates. Additionally, the species' naturally short lifespan of only a few months, along with its experimental accessibility, offers a robust platform for testing interventions to improve age-related dysfunctions in the whole organism and potentially inform the development of immune-based therapies for human aging-related diseases.

AntiRef: reference clusters of human antibody sequences

Motivation

Genetic biases in the human antibody repertoire result in publicly available antibody sequence datasets that contain many duplicate or highly similar sequences. Available datasets are further skewed by the predominance of studies focused on specific disease states, primarily cancer, autoimmunity, and a small number of infectious diseases that includes HIV, influenza, and SARS-CoV-2. These biases and redundancies are a barrier to rapid similarity searches and reduce the efficiency with which these datasets can be used to train statistical or machine-learning models. Identity-based clustering provides a solution; however, the extremely large size of available antibody sequence datasets makes such clustering operations computationally intensive and potentially out of reach for many scientists and researchers who would benefit from such data.

Results

Antibody Reference Clusters (AntiRef), which is modeled after UniRef, provides clustered datasets of filtered human antibody sequences. Due to the modular nature of recombined antibody genes, the clustering thresholds used by UniRef for general protein sequences are suboptimal for antibody clustering. Starting with an input dataset of ∼451M full-length, productive human antibody sequences, AntiRef provides reference datasets clustered at a range of antibody-optimized identity thresholds. AntiRef90 is one-third the size of the input dataset and less than half the size of the non-redundant AntiRef100.

Availability and implementation

AntiRef datasets are available on Zenodo (zenodo.org/record/7474336). All code used to generate AntiRef is available on GitHub (github.com/briney/antiref). The AntiRef versioning scheme (current version: v2022.12.14) refers to the date on which sequences were retrieved from OAS.

Characterization of RBD-specific cross-neutralizing antibodies responses against SARS-CoV-2 variants from COVID-19 convalescents

Introduction

The novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic has been posing a severe threat to global public health. Although broadly neutralizing antibodies have been used to prevent or treat corona virus disease 2019 (COVID-19), new emerging variants have been proven resistant to these antibodies.

Methods

In this study, we isolated receptor binding domain (RBD)-specific memory B cells using single-cell sorting method from two COVID-19 convalescents and expressed the antibody to test their neutralizing activity against diverse SARS-CoV-2 variants. Then, we resolved antibody-RBD complex structures of potent RBD-specific neutralizing antibodies by X-ray diffraction method. Finally, we analyzed the whole antibody repertoires of the two donors and studied the evolutionary pathway of potent neutralizing antibodies.

Results and discussion

We identified three potent RBD-specific neutralizing antibodies (1D7, 3G10 and 3C11) from two COVID-19 convalescents that neutralized authentic SARS-CoV-2 WH-1 and Delta variant, and one of them, 1D7, presented broadly neutralizing activity against WH-1, Beta, Gamma, Delta and Omicron authentic viruses. The resolved antibody-RBD complex structures of two antibodies, 3G10 and 3C11, indicate that both of them interact with the external subdomain of the RBD and that they belong to the RBD-1 and RBD-4 communities, respectively. From the antibody repertoire analysis, we found that the CDR3 frequencies of the light chain, which shared high degrees of amino acid identity with these three antibodies, were higher than those of the heavy chain. This research will contribute to the development of RBD-specific antibody-based drugs and immunogens against multiple variants.

Potential of TCR sequencing in graft-versus-host disease

Graft-versus-host disease (GvHD) remains one of the major complications following allogeneic haematopoietic stem cell transplantation (allo-HSCT). GvHD can occur in almost every tissue, with the skin, liver, and intestines being the mainly affected organs. T cells are implicated in initiating GvHD. T cells identify a broad range of antigens and mediate the immune response through receptors on their surfaces (T cell receptors, TCRs). The composition of TCRs within a T cell population defines the TCR repertoire of an individual, and this repertoire represents exposure to self and non-self proteins. Monitoring the changes in the TCR repertoire using TCR sequencing can provide an indication of the dynamics of a T cell population. Monitoring the frequency and specificities of specific TCR clonotypes longitudinally in different conditions and specimens (peripheral blood, GvHD-affected tissue samples) can provide insights into factors modulating immune reactions following allogeneic transplantation and will help to understand the underlying mechanisms mediating GvHD. This review provides insights into current studies of the TCR repertoire in GvHD and potential future clinical implications of TCR sequencing.

Molecular characterization and expression analysis of immunoglobulins (IgM and IgT) heavy chains in black rockfish (Sebastes schlegelii) that response to bacterial challenge

Sebastes schlegelii is a kind of fish with great economic values. Recently, with the continuous expansion of aquaculture scale and the continuous improvement of aquaculture density, outbreak of various diseases has caused huge economic losses to its aquaculture industry. Study of fish immune system can help to understand the mechanism of immune response to external pathogens and can promote the development of immune prevention and control methods. Immunoglobulins (Igs) are complex glycoproteins that appear to be unique to the vertebrates that can recognize a wide variety of pathogens and recruit immune cells and molecules to destroy pathogens, which are generated by a series of rearrangement and somatic mutations. We therefore studied the immunoglobulin genes of S. schlegelii in view of their important roles in resisting to external pathogen infections. In this study, the immunoglobulin heavy chain genes (sIgM, mIgM, sIgT, and mIgT) of S. schlegelii were successfully identified and cloned. Phylogenetic analysis showed that the IgM and IgT genes of S. schlegelii were clustered together with homologous genes of other species, indicating that they were highly conserved during the evolutionary process. Collinearity analysis showed that the immunoglobulin genes and their adjacent genes were aligned with zebrafish, Atlantic salmon and tilapia, which further confirmed the conserved immunoglobulin gene of teleost. Expression analysis of healthy tissues showed that the expression levels of sIgM, sIgT and mIgT were the highest in the skin, while mIgM was the highest in spleen. After different bacterial infection, IgM and IgT were significantly expressed in skin and gill, which may be because skin and gill are the first line of defense against the infection pathogens. Subcellular localization showed that the mIgT protein was expressed in both the cell membrane and cytoplasm. Meanwhile, recombinant protein of mIgT was obtained in vitro, which laid a foundation for subsequent protein function studies. These results provide a theoretical basis for understanding the immunity role of immunoglobulin in S. schlegelii.

Immunoglobulin VDJ repertoires reveal hallmarks of germinal centers in unique cell clusters isolated from zebrafish (Danio rerio) lymphoid tissues

DNA mutagenesis during antibody affinity maturation has potentially oncogenic or autoimmune outcomes if not tightly controlled as it is in mammalian germinal centers. Cold blooded vertebrates lack germinal centers, yet have a functional Ig gene mutator enzyme, Aicda. In fish there are clusters of Aicda⁺ cells encircled by pigmented 'melano-macrophages' and we test the hypothesis that these clusters are functionally analogous to germinal centers. Sequenced IgH VDJ repertoire libraries from individual isolated clusters showed evidence of B-cell clonal expansion and VDJ somatic hypermutation. Construction of Ig clonal lineage trees revealed that unlike surrounding lymphoid tissue, each cluster is dominated by a few B-cell VDJ clonotypes having hundreds of mutated variants. Recruitment of B-cells to the clusters appears to be ongoing, as there are additional Ig clones having smaller lineages. Finally, we show evidence for positive selection for replacement mutations in regions encoding the antigen contact loops, but not in the framework regions, consistent with functional antibody modification. Melano-macrophages appear to trap the Ag used for post-mutation B-cell selection, performing a role analogous to the follicular dendritic cells of mammalian germinal centers. These findings provide insights into the evolution of the affinity maturation process, the improvement of fish vaccines and possibly also the workings of atypical ectopic germinal centers generated in several human diseases.

Mutational patterns along different evolution paths of follicular lymphoma

Follicular lymphoma (FL) is an indolent disease, characterized by a median life expectancy of 18-20 years and by intermittent periods of relapse and remission. FL frequently transforms into the more aggressive diffuse large B cell lymphoma (t-FL). In previous studies, the analysis of immunoglobulin heavy chain variable region (IgHV) genes in sequential biopsies from the same patient revealed two different patterns of tumor clonal evolution: direct evolution, through acquisition of additional IgHV mutations over time, or divergent evolution, in which lymphoma clones from serial biopsies independently develop from a less-mutated common progenitor cell (CPC). Our goal in this study was to characterize the somatic hypermutation (SHM) patterns of IgHV genes in sequential FL samples from the same patients, and address the question of whether the mutation mechanisms (SHM targeting, DNA repair or both), or selection forces acting on the tumor clones, were different in FL samples compared to healthy control samples, or in late relapsed/transformed FL samples compared to earlier ones. Our analysis revealed differences in the distribution of mutations from each of the nucleotides when tumor and non-tumor clones were compared, while FL and transformed FL (t-FL) tumor clones displayed similar mutation distributions. Lineage tree measurements suggested that either initial clone affinity or selection thresholds were lower in FL samples compared to controls, but similar between FL and t-FL samples. Finally, we observed that both FL and t-FL tumor clones tend to accumulate larger numbers of potential N-glycosylation sites due to the introduction of new SHM. Taken together, these results suggest that transformation into t-FL, in contrast to initial FL development, is not associated with any major changes in DNA targeting or repair, or the selection threshold of the tumor clone.

NoisET: Noise Learning and Expansion Detection of T-Cell Receptors

High-throughput sequencing of T- and B-cell receptors makes it possible to track immune repertoires across time, in different tissues, in acute and chronic diseases and in healthy individuals. However, quantitative comparison between repertoires is confounded by variability in the read count of each receptor clonotype due to sampling, library preparation, and expression noise. We review methods for accounting for both biological and experimental noise and present an easy-to-use python package NoisET that implements and generalizes a previously developed Bayesian method. It can be used to learn experimental noise models for repertoire sequencing from replicates, and to detect responding clones following a stimulus. We test the package on different repertoire sequencing technologies and data sets. We review how such approaches have been used to identify responding clonotypes in vaccination and disease data. Availability: NoisET is freely available to use with source code at github.com/statbiophys/NoisET.

Phenotypic determinism and stochasticity in antibody repertoires of clonally expanded plasma cells

B cell clonal selection and expansion from a genetically diverse antibody repertoire guides the immune response to a target antigen. It remains unclear if clonal selection and expansion follow any deterministic rules or are stochastic with regards to phenotypic antibody properties such as antigen-binding, affinity, and epitope specificity. We perform the in-depth genotypic and phenotypic characterization of antibody repertoires following immunization in mice. We identify the degree to which clonal expansion is driven by antibody binding, affinity, and epitope specificity and as such may provide greater insight into vaccine-induced immunity.

The capacity of humoral B cell-mediated immunity to effectively respond to and protect against pathogenic infections is largely driven by the presence of a diverse repertoire of polyclonal antibodies in the serum, which are produced by plasma cells (PCs). Recent studies have started to reveal the balance between deterministic mechanisms and stochasticity of antibody repertoires on a genotypic level (i.e., clonal diversity, somatic hypermutation, and germline gene usage). However, it remains unclear if clonal selection and expansion of PCs follow any deterministic rules or are stochastic with regards to phenotypic antibody properties (i.e., antigen-binding, affinity, and epitope specificity). Here, we report on the in-depth genotypic and phenotypic characterization of clonally expanded PC antibody repertoires following protein immunization. We find that clonal expansion drives antigen specificity of the most expanded clones (top ∼10), whereas among the rest of the clonal repertoire antigen specificity is stochastic. Furthermore, we report both on a polyclonal repertoire and clonal lineage level that antibody-antigen binding affinity does not correlate with clonal expansion or somatic hypermutation. Last, we provide evidence for convergence toward targeting dominant epitopes despite clonal sequence diversity among the most expanded clones. Our results highlight the extent to which clonal expansion can be ascribed to antigen binding, affinity, and epitope specificity, and they have implications for the assessment of effective vaccines.

How Naive T-Cell Clone Counts Are Shaped By Heterogeneous Thymic Output and Homeostatic Proliferation

The specificity of T cells is that each T cell has only one T cell receptor (TCR). A T cell clone represents a collection of T cells with the same TCR sequence. Thus, the number of different T cell clones in an organism reflects the number of different T cell receptors (TCRs) that arise from recombination of the V(D)J gene segments during T cell development in the thymus. TCR diversity and more specifically, the clone abundance distribution, are important factors in immune functions. Specific recombination patterns occur more frequently than others while subsequent interactions between TCRs and self-antigens are known to trigger proliferation and sustain naive T cell survival. These processes are TCR-dependent, leading to clone-dependent thymic export and naive T cell proliferation rates. We describe the heterogeneous steady-state population of naive T cells (those that have not yet been antigenically triggered) by using a mean-field model of a regulated birth-death-immigration process. After accounting for random sampling, we investigate how TCR-dependent heterogeneities in immigration and proliferation rates affect the shape of clone abundance distributions (the number of different clones that are represented by a specific number of cells, or “clone counts”). By using reasonable physiological parameter values and fitting predicted clone counts to experimentally sampled clone abundances, we show that realistic levels of heterogeneity in immigration rates cause very little change to predicted clone-counts, but that modest heterogeneity in proliferation rates can generate the observed clone abundances. Our analysis provides constraints among physiological parameters that are necessary to yield predictions that qualitatively match the data. Assumptions of the model and potentially other important mechanistic factors are discussed.

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.

•

Ultra-deep sequencing of both IgM and IgG repertoires in chronic HIV-1 infection.

•

VDJ gene rearrangement patterns can be dramatically changed by HIV-1 infection.

•

Multiple mechanisms cause the high complexity of HIV-1-experienced antibodies.

•

Discovery of promising neutralizing HIV-1 antibodies from antibody repertoires.

Affinity maturation for an optimal balance between long-term immune coverage and short-term resource constraints

Humoral immunity relies on the mutation and selection of B cells to better recognize pathogens. This affinity maturation process produces cells with diverse recognition capabilities. Examining optimal immune strategies that maximize the long-term immune coverage at a minimal metabolic cost, we show when the immune system should mount a de novo response rather than rely on existing memory cells. Our theory recapitulates known modes of the B cell response, predicts the empirical form of the distribution of clone sizes, and rationalizes as a trade-off between metabolic and immune costs the antigenic imprinting effects that limit the efficacy of vaccines (original antigenic sin). Our predictions provide a framework to interpret experimental results that could be used to inform vaccination strategies.

In order to target threatening pathogens, the adaptive immune system performs a continuous reorganization of its lymphocyte repertoire. Following an immune challenge, the B cell repertoire can evolve cells of increased specificity for the encountered strain. This process of affinity maturation generates a memory pool whose diversity and size remain difficult to predict. We assume that the immune system follows a strategy that maximizes the long-term immune coverage and minimizes the short-term metabolic costs associated with affinity maturation. This strategy is defined as an optimal decision process on a finite dimensional phenotypic space, where a preexisting population of cells is sequentially challenged with a neutrally evolving strain. We show that the low specificity and high diversity of memory B cells—a key experimental result—can be explained as a strategy to protect against pathogens that evolve fast enough to escape highly potent but narrow memory. This plasticity of the repertoire drives the emergence of distinct regimes for the size and diversity of the memory pool, depending on the density of de novo responding cells and on the mutation rate of the strain. The model predicts power-law distributions of clonotype sizes observed in data and rationalizes antigenic imprinting as a strategy to minimize metabolic costs while keeping good immune protection against future strains.

Extensive age-dependent loss of antibody diversity in naturally short-lived turquoise killifish

Aging individuals exhibit a pervasive decline in adaptive immune function, with important implications for health and lifespan. Previous studies have found a pervasive loss of immune-repertoire diversity in human peripheral blood during aging; however, little is known about repertoire aging in other immune compartments, or in species other than humans. Here, we perform the first study of immune-repertoire aging in an emerging model of vertebrate aging, the African turquoise killifish (Nothobranchius furzeri). Despite their extremely short lifespans, these killifish exhibit complex and individualized heavy-chain repertoires, with a generative process capable of producing millions of distinct productive sequences. Whole-body killifish repertoires decline rapidly in within-individual diversity with age, while between-individual variability increases. Large, expanded B-cell clones exhibit far greater diversity loss with age than small clones, suggesting important differences in how age affects different B-cell populations. The immune repertoires of isolated intestinal samples exhibit especially dramatic age-related diversity loss, related to an elevated prevalence of expanded clones. Lower intestinal repertoire diversity was also associated with transcriptomic signatures of reduced B-cell activity, supporting a functional role for diversity changes in killifish immunosenescence. Our results highlight important differences in systemic vs. organ-specific aging dynamics in the adaptive immune system.

The Intra-Tumoral T Cell Receptor Repertoire: Steps Towards a Useful Clinical Biomarker

Adaptive immunity recognizes and responds to tumors, although they are part of the immunological "self." T cells, both CD4+ and CD8+, play a key role in the process, and the specific set of receptors which recognize tumor antigens therefore has the potential to provide prognostic biomarkers for tracking tumor growth after cancer therapy, including immunotherapy. Most published data on the T cell repertoire continue to rely on commercial proprietary methods, which often do not allow access to the raw data, and are difficult to validate. We describe an open-source protocol for amplifying, sequencing, and analyzing T cell receptors which is economical, robust, sensitive, and versatile. The key experimental step is the ligation of a single-stranded oligonucleotide to the 3' end of the T cell receptor cDNA, which allows easy amplification of all possible rearrangements using only a single set of primers per locus, while simultaneously introducing a unique molecular identifier to label each starting cDNA molecule. After sequencing, this molecular identifier can be used to correct both sequence errors and the effects of differential PCR amplification efficiency, thus producing a more accurate measure of the true T cell receptor frequency within the sample. Samples are then tagged with unique pairs of indices, facilitating robotic scale-up and significantly reducing cross-sample contamination from index hopping. This method has been applied to the analysis of tumor-infiltrating lymphocytes and matched peripheral blood samples from patients with a variety of solid tumors.

Individualized VDJ recombination predisposes the available Ig sequence space

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

Antibody upstream sequence diversity and its biological implications revealed by repertoire sequencing

The sequence upstream of the antibody variable region (antibody upstream sequence [AUS]) consists of a 5' untranslated region (5' UTR) and a preceding leader region. The sequence variations in AUS affect antibody engineering and PCR based antibody quantification and may also be implicated in mRNA transcription and translation. However, the diversity of AUSs remains elusive. Using 5' rapid amplification of cDNA ends and high-throughput antibody repertoire sequencing technique, we acquired full-length AUSs for human, rhesus macaque, cynomolgus macaque, mouse, and rat. We designed a bioinformatics pipeline and identified 3307 unique AUSs, corresponding to 3026 and 1457 unique sequences for 5' UTR and leader region, respectively. Comparative analysis indicated that 928 (63.69%) leader sequences are novel relative to those recorded in the international ImMunoGeneTics information system. Evolutionarily, leader sequences are more conserved than 5' UTR and seem to coevolve with their downstream V genes. Besides, single-nucleotide polymorphisms are position dependent for leader regions and may contribute to the functional reversal of the downstream V genes. Finally, the AUGs in AUSs were found to have little impact on gene expression. Taken together, our findings can facilitate primer design for capturing antibodies efficiently and provide a valuable resource for antibody engineering and molecule-level antibody studies.

Isolation of monoclonal antibodies from anti-synthetase syndrome patients and affinity maturation by recombination of independent somatic variants

The autoimmune disease known as Jo-1 positive anti-synthetase syndrome (ASS) is characterized by circulating antibody titers to histidyl-tRNA synthetase (HARS), which may play a role in modulating the non-canonical functions of HARS. Monoclonal antibodies to HARS were isolated by single-cell screening and sequencing from three Jo-1 positive ASS patients and shown to be of high affinity, covering diverse epitope space. The immune response was further characterized by repertoire sequencing from the most productive of the donor samples. In line with previous studies of autoimmune repertoires, these antibodies tended to have long complementarity-determining region H3 sequences with more positive-charged residues than average. Clones of interest were clustered into groups with related sequences, allowing us to observe different somatic mutations in related clones. We postulated that these had found alternate structural solutions for high affinity binding, but that mutations might be transferable between clones to further enhance binding affinity. Transfer of somatic mutations between antibodies within the same clonal group was able to enhance binding affinity in a number of cases, including beneficial transfer of a mutation from a lower affinity clone into one of higher affinity. Affinity enhancement was seen with mutation transfer both between related single-cell clones, and directly from related repertoire sequences. To our knowledge, this is the first demonstration of somatic hypermutation transfer from repertoire sequences to further mature in vivo derived antibodies, and represents an additional tool to aid in affinity maturation for the development of antibodies.

Short Read-Length Next Generation DNA Sequencing of Antibody CDR Combinations from Phage Selection Outputs

Phage display is commonly used to select target-binding antibody fragments from large libraries containing billions of unique antibody clones. In practice, selection outputs are often highly heterogenous, making it desirable to recover sequence information from the selected pool. Next Generation DNA Sequencing (NGS) enables the acquisition of sufficient sequencing reads to cover the pool diversity, however read-lengths are typically too short to capture paired antibody complementarity-determining regions (CDRs), which is needed to reconstruct target-binding antibody fragments. Here, we describe a simple in vitro protocol to bring the DNA encoding the antibody CDRs closer together. The final PCR product referred to as a "CDR strip" is suitable for short read-length NGS. In this method, phagemid ssDNA is recovered from antibody phage display biopanning and used as a template to create a heteroduplex with deletions between CDRs of interest. The shorter strand in the heteroduplex is preferentially PCR amplified to generate a CDR strip that is sequenced using NGS. We have also included a bioinformatics approach to analyze the CDR strip populations so that single antibody clones can be created from paired CDR sequences.

Cutting Edge: Neutralizing Public Antibody Responses Are an Ancient Form of Defense Conserved in Fish and Mammals

The repertoire of Abs is generated by genomic rearrangements during B cell differentiation. Although V(D)J rearrangements lead to repertoires mostly different between individuals, recent studies have shown that they contain a substantial fraction of overrepresented and shared "public" clones. We previously reported a strong public IgHμ clonotypic response against the rhabdovirus viral hemorrhagic septicemia virus in a teleost fish. In this study, we identified an IgL chain associated with this public response that allowed us to characterize its functionality. We show that this public Ab response has a potent neutralizing capacity that is typically associated with host protection during rhabdovirus infections. We also demonstrate that the public response is not restricted to a particular trout isogenic line but expressed in multiple genetic backgrounds and may be used as a marker of successful vaccination. Our work reveals that public B cell responses producing generic Abs constitute a mechanism of protection against infection conserved across vertebrates.

Human immune repertoire in hepatitis B virus infection

Hepatitis B virus (HBV) infection is a public health threat that affects 257 million people worldwide and can progress to liver cirrhosis, liver failure, and hepatocellular carcinoma. The HBV antigen- induced adaptive immune response plays an important role in HBV clearance. Immune repertoire sequencing (IRS) has been used to investigate the molecular mechanisms behind the immune system, find novel ways to treat HBV infection, and evaluate the genetic responses and immune characteristics of individuals infected by HBV or immunized by HBV vaccine. This review summarizes the human immune repertoire analysis methodology, and the application of the IRS in the prediction of HBV infection progression, treatment, and vaccination.

Do ectothermic vertebrates have a home in which to affinity mature their antibody responses?

There has been a longstanding question of whether affinity maturation occurs in ectotherms, and if it does, where in tissues this happens. Although cold-blooded vertebrates (ectotherms) lack histologically discernible germinal centers, they have a fully functional Ig gene mutator enzyme (activation-induced cytidine deaminase: AID or Aicda). Protein and Ig cDNA transcript analyses provide evidence that ectotherms can, under certain conditions, demonstrate antibody affinity maturation, and somatic hypermutation of their Ig genes during secondary immune responses. Here, we review the evidence for antibody affinity maturation and somatic hypermutation of Ig V(D)J exons. We argue that past evidence of long-term intact antigen retention, and recent studies of in situ expression of AID transcripts, point to fish melanomacrophage clusters as sites functionally analogous to a germinal center. Recent work in zebrafish provides a way forward to test these predictions through V(D)J repertoire analyses on isolated, intact melanomacrophage clusters. This work has implications not only for vaccine use in aquaculture, but also for antibody affinity maturation processes in all ectothermic vertebrates.

Diversity in immunogenomics: the value and the challenge

Immunogenomics studies have been largely limited to individuals of European ancestry, restricting the ability to identify variation in human adaptive immune responses across populations. Inclusion of a greater diversity of individuals in immunogenomics studies will substantially enhance our understanding of human immunology.

Biological controls for standardization and interpretation of adaptive immune receptor repertoire profiling

Use of adaptive immune receptor repertoire sequencing (AIRR-seq) has become widespread, providing new insights into the immune system with potential broad clinical and diagnostic applications. However, like many high-throughput technologies, it comes with several problems, and the AIRR Community was established to understand and help solve them. We, the AIRR Community’s Biological Resources Working Group, have surveyed scientists about the need for standards and controls in generating and annotating AIRR-seq data. Here, we review the current status of AIRR-seq, provide the results of our survey, and based on them, offer recommendations for developing AIRR-seq standards and controls, including future work.

High-Throughput Sequencing of Phage Display Libraries Reveals Parasitic Enrichment of Indel Mutants Caused by Amplification Bias

The combination of phage display technology with high-throughput sequencing enables in-depth analysis of library diversity and selection-driven dynamics. We applied short-read sequencing of the mutagenized region on focused display libraries of two homologous nucleic acid modification eraser proteins-AlkB and FTO-biopanned against methylated DNA. This revealed enriched genotypes with small indels and concomitant doubtful amino acid motifs within the FTO library. Nanopore sequencing of the entire display vector showed additional enrichment of large deletions overlooked by region-specific sequencing, and further impacted the interpretation of the obtained amino acid motifs. We could attribute enrichment of these corrupted clones to amplification bias due to arduous FTO display slowing down host cell growth as well as phage production. This amplification bias appeared to be stronger than affinity-based target selection. Recommendations are provided for proper sequence analysis of phage display data, which can improve motive discovery in libraries of proteins that are difficult to display.

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

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

High Frequency of Shared Clonotypes in Human T Cell Receptor Repertoires

The collection of T cell receptors (TCRs) generated by somatic recombination is large but unknown. We generate large TCR repertoire datasets as a resource to facilitate detailed studies of the role of TCR clonotypes and repertoires in health and disease. We estimate the size of individual human recombined and expressed TCRs by sequence analysis and determine the extent of sharing between individual repertoires. Our experiments reveal that each blood sample contains between 5 million and 21 million TCR clonotypes. Three individuals share 8% of TCRβ- or 11% of TCRα-chain clonotypes. Sorting by T cell phenotypes in four individuals shows that 5% of naive CD4+ and 3.5% of naive CD8+ subsets share their TCRβ clonotypes, whereas memory CD4+ and CD8+ subsets share 2.3% and 0.4% of their clonotypes, respectively. We identify the sequences of these shared TCR clonotypes that are of interest for studies of human T cell biology.

Soto et al. examine the extent to which five healthy adults share their T cell receptor (TCR) repertoire. Using sequencing and bioinformatics, they show a high prevalence of shared clonotypes even considering different T cell phenotypes. Possible functions for some clonotypes are inferred based on homology with TCRs in GenBank.

Applications of Machine and Deep Learning in Adaptive Immunity

Adaptive immunity is mediated by lymphocyte B and T cells, which respectively express a vast and diverse repertoire of B cell and T cell receptors and, in conjunction with peptide antigen presentation through major histocompatibility complexes (MHCs), can recognize and respond to pathogens and diseased cells. In recent years, advances in deep sequencing have led to a massive increase in the amount of adaptive immune receptor repertoire data; additionally, proteomics techniques have led to a wealth of data on peptide-MHC presentation. These large-scale data sets are now making it possible to train machine and deep learning models, which can be used to identify complex and high-dimensional patterns in immune repertoires. This article introduces adaptive immune repertoires and machine and deep learning related to biological sequence data and then summarizes the many applications in this field, which span from predicting the immunological status of a host to the antigen specificity of individual receptors and the engineering of immunotherapeutics.

High-throughput sequencing of the porcine antibody repertoire with or without PEDV infection: A proof-of-concept study

A host's adaptive immune system can protect against a wide variety of pathogens by producing diverse antibodies. The antibody repertoire is so vast that traditional low-throughput methods cannot fully sequence it. In this study, we developed a high-throughput sequencing (HTS) method for antibody repertoire assessment in swine, and tested it with or without porcine epidemic diarrhea virus (PEDV) infection. We isolated peripheral blood mononuclear cells from normal or PEDV-seropositive pigs and applied multiplex PCR to amplify the porcine B cell receptor heavy chain library, followed by HTS using the Illumina Miseq system to obtain full sequence information. The results from sequence analysis demonstrated that in normal conditions, several V gene segments were preferentially used, with IGHV1-4 and IGHV1S2 being the two most frequent. The IGHV usage in PEDV-seropositive pigs was not exactly the same as that of PEDV-seronegative pigs, with an increased usage of IGHV1-6. Our study provides an effective approach to comprehensively understand the overall porcine antibody repertoire, as well as to monitor broad antibody responses to viral challenge in pigs.

Restricted epitope specificity determined by variable region germline segment pairing in rodent antibody repertoires

Antibodies from B-cell clonal lineages share sequence and structural properties as well as epitope specificity. Clonally unrelated antibodies can similarly share sequence and specificity properties and are said to be convergent. Convergent antibody responses against several antigens have been described in humans and mice and include different classes of shared sequence features. In particular, some antigens and epitopes can induce convergent responses of clonally unrelated antibodies with restricted heavy (V_H) and light (V_L) chain variable region germline segment usage without similarity in the heavy chain third complementarity-determining region (CDR H3), a critical specificity determinant. Whether these V germline segment-restricted responses reflect a general epitope specificity restriction of antibodies with shared V_H/V_L pairing is not known. Here, we investigated this question by determining patterns of antigen binding competition between clonally unrelated antigen-specific rat antibodies from paired-chain deep sequencing datasets selected based solely on V_H/V_L pairing. We found that antibodies with shared V_H/V_L germline segment pairings but divergent CDR H3 sequences almost invariably have restricted epitope specificity indicated by shared binding competition patterns. This epitope restriction included 82 of 85 clonally unrelated antibodies with 13 different V_H/V_L pairings binding in 8 epitope groups in 2 antigens. The corollary that antibodies with shared V_H/V_L pairing and epitope-restricted binding can accommodate widely divergent CDR H3 sequences was confirmed by in vitro selection of variants of anti-human epidermal growth factor receptor 2 antibodies known to mediate critical antigen interactions through CDR H3. Our results show that restricted epitope specificity determined by V_H/V_L germline segment pairing is a general property of rodent antigen-specific antibodies.

Immune fingerprinting through repertoire similarity

Immune repertoires provide a unique fingerprint reflecting the immune history of individuals, with potential applications in precision medicine. However, the question of how personal that information is and how it can be used to identify individuals has not been explored. Here, we show that individuals can be uniquely identified from repertoires of just a few thousands lymphocytes. We present "Immprint," a classifier using an information-theoretic measure of repertoire similarity to distinguish pairs of repertoire samples coming from the same versus different individuals. Using published T-cell receptor repertoires and statistical modeling, we tested its ability to identify individuals with great accuracy, including identical twins, by computing false positive and false negative rates < 10-6 from samples composed of 10,000 T-cells. We verified through longitudinal datasets that the method is robust to acute infections and that the immune fingerprint is stable for at least three years. These results emphasize the private and personal nature of repertoire data.

Antigen receptor repertoires of one of the smallest known vertebrates

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

Antibody-guided structure-based vaccines

The vaccine field is pursuing diverse approaches to translate the molecular insights from analyses of effective antibodies and their targeted epitopes into immunogens capable of eliciting protective immune responses. Here we review current antibody-guided strategies including conformation-based, epitope-based, and lineage-based vaccine approaches, which are yielding promising vaccine candidates now being evaluated in clinical trials. We summarize directions being employed by the field, including the use of sequencing technologies to monitor and track developing immune responses for understanding and improving antibody-based immunity. We review opportunities and challenges to transform powerful new discoveries into safe and effective vaccines, which are encapsulated by vaccine efforts against a variety of pathogens including HIV-1, influenza A virus, malaria parasites, respiratory syncytial virus, and SARS-CoV-2. Overall, this review summarizes the extensive progress that has been made to realize antibody-guided structure-based vaccines, the considerable challenges faced, and the opportunities afforded by recently developed molecular approaches to vaccine development.

Learning the heterogeneous hypermutation landscape of immunoglobulins from high-throughput repertoire data

Somatic hypermutations of immunoglobulin (Ig) genes occurring during affinity maturation drive B-cell receptors’ ability to evolve strong binding to their antigenic targets. The landscape of these mutations is highly heterogeneous, with certain regions of the Ig gene being preferentially targeted. However, a rigorous quantification of this bias has been difficult because of phylogenetic correlations between sequences and the interference of selective forces. Here, we present an approach that corrects for these issues, and use it to learn a model of hypermutation preferences from a recently published large IgH repertoire dataset. The obtained model predicts mutation profiles accurately and in a reproducible way, including in the previously uncharacterized Complementarity Determining Region 3, revealing that both the sequence context of the mutation and its absolute position along the gene are important. In addition, we show that hypermutations occurring concomittantly along B-cell lineages tend to co-localize, suggesting a possible mechanism for accelerating affinity maturation.

Defining and Manipulating B Cell Immunodominance Hierarchies to Elicit Broadly Neutralizing Antibody Responses against Influenza Virus

The antibody repertoire possesses near-limitless diversity, enabling the adaptive immune system to accommodate essentially any antigen. However, this diversity explores the antigenic space unequally, allowing some pathogens like influenza virus to impose complex immunodominance hierarchies that distract antibody responses away from key sites of virus vulnerability. We developed a computational model of affinity maturation to map the patterns of immunodominance that evolve upon immunization with natural and engineered displays of hemagglutinin (HA), the influenza vaccine antigen. Based on this knowledge, we designed immunization protocols that subvert immune distraction and focus serum antibody responses upon a functionally conserved, but immunologically recessive, target of human broadly neutralizing antibodies. We tested in silico predictions by vaccinating transgenic mice in which antibody diversity was humanized to mirror clinically relevant humoral output. Collectively, our results demonstrate that complex patterns in antibody immunogenicity can be rationally defined and then manipulated to elicit engineered immunity.

Immune Literacy: Reading, Writing, and Editing Adaptive Immunity

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

PIRD: Pan Immune Repertoire Database

MOTIVATION

T and B cell receptors (TCRs and BCRs) play a pivotal role in the adaptive immune system by recognizing an enormous variety of external and internal antigens. Understanding these receptors is critical for exploring the process of immunoreaction and exploiting potential applications in immunotherapy and antibody drug design. Although a large number of samples have had their TCR and BCR repertoires sequenced using high-throughput sequencing in recent years, very few databases have been constructed to store these kinds of data. To resolve this issue, we developed a database.

RESULTS

We developed a database, the Pan Immune Repertoire Database (PIRD), located in China National GeneBank (CNGBdb), to collect and store annotated TCR and BCR sequencing data, including from Homo sapiens and other species. In addition to data storage, PIRD also provides functions of data visualization and interactive online analysis. Additionally, a manually curated database of TCRs and BCRs targeting known antigens (TBAdb) was also deposited in PIRD.

AVAILABILITY AND IMPLEMENTATION

PIRD can be freely accessed at https://db.cngb.org/pird.

Ultra-efficient sequencing of T Cell receptor repertoires reveals shared responses in muscle from patients with Myositis

BACKGROUND

Myositis, or idiopathic inflammatory myopathy (IIM), is a group disorders of unknown etiology characterized by the inflammation of skeletal muscle. The role of T cells and their antigenic targets in IIM initiation and progression is poorly understood. T cell receptor (TCR) repertoire sequencing is a powerful approach for characterizing complex T cell responses. However, current TCR sequencing methodologies are complex, expensive, or both, greatly limiting the scale of feasible studies.

METHODS

Here we present Framework Region 3 AmplifiKation sequencing ("FR3AK-seq"), a simplified multiplex PCR-based approach for the ultra-efficient and quantitative analysis of TCR complementarity determining region 3 (CDR3) repertoires. By using minimal primer sets targeting a conserved region immediately upstream of CDR3, undistorted amplicons are analyzed via short read, single-end sequencing. We also introduce the novel algorithm Inferring Sequences via Efficiency Projection and Primer Incorporation ("ISEPPI") for linking CDR3s to their associated variable genes.

FINDINGS

We find that FR3AK-seq is sensitive and quantitative, performing comparably to two different industry standards. FR3AK-seq and ISEPPI were used to efficiently and inexpensively characterize the T cell infiltrates of surgical muscle biopsies obtained from 145 patients with IIM and controls. A cluster of closely related TCRs was identified in samples from patients with sporadic inclusion body myositis (IBM).

INTERPRETATION

The ease and minimal cost of FR3AK-seq removes critical barriers to routine, large-scale TCR CDR3 repertoire analyses, thereby democratizing the quantitative assessment of human TCR repertoires in disease-relevant target tissues. Importantly, discovery of closely related TCRs in muscle from patients with IBM provides evidence for a shared antigen-driven T cell response in this disease of unknown pathogenesis.

FUNDING

This work was supported by NIH grant U24AI118633 and a Prostate Cancer Foundation Young Investigator Award.

FcRn Overexpression Expands Diversity of the Humoral Immune Response in bFcRn Transgenic Mice

The neonatal Fc receptor (FcRn) plays key roles in IgG and albumin homeostasis, maternal IgG transport, and antigen presentation of IgG-opsonized antigens. Previously, we reported that transgenic (Tg) mice that overexpress bovine FcRn (bFcRn) have augmented T-dependent humoral immune response with increased IgG protection, higher level of antigen-specific antibodies, greater number of antigen-specific B cells, and effective immune response even against weakly immunogenic epitopes. In this study we analyzed the diversity of the humoral immune response of bFcRn Tg mice, using a length distribution analysis (spectratyping) and next generation sequencing (NGS) of the immunoglobulin heavy chain variable regions. Our analysis showed that in response to immunization with ovalbumin or transfected cells that expressed a unique membrane protein, our Tg animals developed a more diverse plasma cell repertoire than controls, which manifested in greater numbers of different clones, and clusters with fewer highly expanded large clones, as identified by the variable region (CDR3) of the immunoglobulin heavy chain. The increased antibody diversity in Tg mice after immunization was observed at both IgM and IgG levels, indicating that the increased humoral immune diversity in Tg mice is due to a higher number of both activated, antigen-specific naïve and isotype switched B cells. We thus demonstrated that the BCR repertoire of the immunized bFcRn Tg animals is more diverse compared to wild type mice, which likely makes these Tg mice a better choice for monoclonal antibody production against challenging antigens, including the extracellular regions of cell membrane proteins.

PyIR: a scalable wrapper for processing billions of immunoglobulin and T cell receptor sequences using IgBLAST

BACKGROUND

Recent advances in DNA sequencing technologies have enabled significant leaps in capacity to generate large volumes of DNA sequence data, which has spurred a rapid growth in the use of bioinformatics as a means of interrogating antibody variable gene repertoires. Common tools used for annotation of antibody sequences are often limited in functionality, modularity and usability.

RESULTS

We have developed PyIR, a Python wrapper and library for IgBLAST, which offers a minimal setup CLI and API, FASTQ support, file chunking for large sequence files, JSON and Python dictionary output, and built-in sequence filtering.

CONCLUSIONS

PyIR offers improved processing speed over multithreaded IgBLAST (version 1.14) when spawning more than 16 processes on a single computer system. Its customizable filtering and data encapsulation allow it to be adapted to a wide range of computing environments. The API allows for IgBLAST to be used in customized bioinformatics workflows.

Sequence-dependent dynamics of synthetic and endogenous RSSs in V(D)J recombination

Developing lymphocytes of jawed vertebrates cleave and combine distinct gene segments to assemble antigen-receptor genes. This process called V(D)J recombination that involves the RAG recombinase binding and cutting recombination signal sequences (RSSs) composed of conserved heptamer and nonamer sequences flanking less well-conserved 12- or 23-bp spacers. Little quantitative information is known about the contributions of individual RSS positions over the course of the RAG-RSS interaction. We employ a single-molecule method known as tethered particle motion to track the formation, lifetime and cleavage of individual RAG-12RSS-23RSS paired complexes (PCs) for numerous synthetic and endogenous 12RSSs. We reveal that single-bp changes, including in the 12RSS spacer, can significantly and selectively alter PC formation or the probability of RAG-mediated cleavage in the PC. We find that some rarely used endogenous gene segments can be mapped directly to poor RAG binding on their adjacent 12RSSs. Finally, we find that while abrogating RSS nicking with Ca2+ leads to substantially shorter PC lifetimes, analysis of the complete lifetime distributions of any 12RSS even on this reduced system reveals that the process of exiting the PC involves unidentified molecular details whose involvement in RAG-RSS dynamics are crucial to quantitatively capture kinetics in V(D)J recombination.

Profiling immunoglobulin repertoires across multiple human tissues using RNA sequencing

Profiling immunoglobulin (Ig) receptor repertoires with specialized assays can be cost-ineffective and time-consuming. Here we report ImReP, a computational method for rapid and accurate profiling of the Ig repertoire, including the complementary-determining region 3 (CDR3), using regular RNA sequencing data such as those from 8,555 samples across 53 tissues types from 544 individuals in the Genotype-Tissue Expression (GTEx v6) project. Using ImReP and GTEx v6 data, we generate a collection of 3.6 million Ig sequences, termed the atlas of immunoglobulin repertoires (TAIR), across a broad range of tissue types that often do not have reported Ig repertoires information. Moreover, the flow of Ig clonotypes and inter-tissue repertoire similarities across immune-related tissues are also evaluated. In summary, TAIR is one of the largest collections of CDR3 sequences and tissue types, and should serve as an important resource for studying immunological diseases.

Information on immune receptor repertoire provides important insights on disease progression and therapy development, but can be expensive and time-consuming to obtain. Here the authors report ImReP, a computational method that can extract detailed immune repertoire information from existing tissue-specific RNA sequencing data.

OLGA: fast computation of generation probabilities of B- and T-cell receptor amino acid sequences and motifs

MOTIVATION

High-throughput sequencing of large immune repertoires has enabled the development of methods to predict the probability of generation by V(D)J recombination of T- and B-cell receptors of any specific nucleotide sequence. These generation probabilities are very non-homogeneous, ranging over 20 orders of magnitude in real repertoires. Since the function of a receptor really depends on its protein sequence, it is important to be able to predict this probability of generation at the amino acid level. However, brute-force summation over all the nucleotide sequences with the correct amino acid translation is computationally intractable. The purpose of this paper is to present a solution to this problem.

RESULTS

We use dynamic programming to construct an efficient and flexible algorithm, called OLGA (Optimized Likelihood estimate of immunoGlobulin Amino-acid sequences), for calculating the probability of generating a given CDR3 amino acid sequence or motif, with or without V/J restriction, as a result of V(D)J recombination in B or T cells. We apply it to databases of epitope-specific T-cell receptors to evaluate the probability that a typical human subject will possess T cells responsive to specific disease-associated epitopes. The model prediction shows an excellent agreement with published data. We suggest that OLGA may be a useful tool to guide vaccine design.

AVAILABILITY AND IMPLEMENTATION

Source code is available at https://github.com/zsethna/OLGA.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Somatic hypermutation analysis for improved identification of B cell clonal families from next-generation sequencing data

Adaptive immune receptor repertoire sequencing (AIRR-Seq) offers the possibility of identifying and tracking B cell clonal expansions during adaptive immune responses. Members of a B cell clone are descended from a common ancestor and share the same initial V(D)J rearrangement, but their B cell receptor (BCR) sequence may differ due to the accumulation of somatic hypermutations (SHMs). Clonal relationships are learned from AIRR-seq data by analyzing the BCR sequence, with the most common methods focused on the highly diverse junction region. However, clonally related cells often share SHMs which have been accumulated during affinity maturation. Here, we investigate whether shared SHMs in the V and J segments of the BCR can be leveraged along with the junction sequence to improve the ability to identify clonally related sequences. We develop independent distance functions that capture junction similarity and shared mutations, and combine these in a spectral clustering framework to infer the BCR clonal relationships. Using both simulated and experimental data, we show that this model improves both the sensitivity and specificity for identifying B cell clones. Source code for this method is freely available in the SCOPer (Spectral Clustering for clOne Partitioning) R package (version 0.2 or newer) in the Immcantation framework: www.immcantation.org under the AGPLv3 license.

Extreme genomic volatility characterizes the evolution of the immunoglobulin heavy chain locus in cyprinodontiform fishes

The evolution of the adaptive immune system has provided vertebrates with a uniquely sophisticated immune toolkit, enabling them to mount precise immune responses against a staggeringly diverse range of antigens. Like other vertebrates, teleost fishes possess a complex and functional adaptive immune system; however, our knowledge of the complex antigen-receptor genes underlying its functionality has been restricted to a small number of experimental and agricultural species, preventing systematic investigation into how these crucial gene loci evolve. Here, we analyse the genomic structure of the immunoglobulin heavy chain (IGH) gene loci in the cyprinodontiforms, a diverse and important group of teleosts present in many different habitats across the world. We reconstruct the complete IGH loci of the turquoise killifish (Nothobranchius furzeri) and the southern platyfish (Xiphophorus maculatus) and analyse their in vivo gene expression, revealing the presence of species-specific splice isoforms of transmembrane IGHM. We further characterize the IGH constant regions of 10 additional cyprinodontiform species, including guppy, Amazon molly, mummichog and mangrove killifish. Phylogenetic analysis of these constant regions suggests multiple independent rounds of duplication and deletion of the teleost-specific antibody class IGHZ in the cyprinodontiform lineage, demonstrating the extreme volatility of IGH evolution. Focusing on the cyprinodontiforms as a model taxon for comparative evolutionary immunology, this work provides novel genomic resources for studying adaptive immunity and sheds light on the evolutionary history of the adaptive immune system.

How repertoire data are changing antibody science

Antibodies are vital proteins of the immune system that recognize potentially harmful molecules and initiate their removal. Mammals can efficiently create vast numbers of antibodies with different sequences capable of binding to any antigen with high affinity and specificity. Because they can be developed to bind to many disease agents, antibodies can be used as therapeutics. In an organism, after antigen exposure, antibodies specific to that antigen are enriched through clonal selection, expansion, and somatic hypermutation. The antibodies present in an organism therefore report on its immune status, describe its innate ability to deal with harmful substances, and reveal how it has previously responded. Next-generation sequencing technologies are being increasingly used to query the antibody, or B-cell receptor (BCR), sequence repertoire, and the amount of BCR data in public repositories is growing. The Observed Antibody Space database, for example, currently contains over a billion sequences from 68 different studies. Repertoires are available that represent both the naive state (i.e. antigen-inexperienced) and that after immunization. This wealth of data has created opportunities to learn more about our immune system. In this review, we discuss the many ways in which BCR repertoire data have been or could be exploited. We highlight its utility for providing insights into how the naive immune repertoire is generated and how it responds to antigens. We also consider how structural information can be used to enhance these data and may lead to more accurate depictions of the sequence space and to applications in the discovery of new therapeutics.

Inter- and intraspecies comparison of phylogenetic fingerprints and sequence diversity of immunoglobulin variable genes

Protection and neutralization of a vast array of pathogens is accomplished by the tremendous diversity of the B cell receptor (BCR) repertoire. For jawed vertebrates, this diversity is initiated via the somatic recombination of immunoglobulin (Ig) germline elements. While it is clear that the number of these germline segments differs from species to species, the extent of cross-species sequence diversity remains largely uncharacterized. Here we use extensive computational and statistical methods to investigate the sequence diversity and evolutionary relationship between Ig variable (V), diversity (D), and joining (J) germline segments across nine commonly studied species ranging from zebrafish to human. Metrics such as guanine-cytosine (GC) content showed low redundancy across Ig germline genes within a given species. Other comparisons, including amino acid motifs, evolutionary selection, and sequence diversity, revealed species-specific properties. Additionally, we showed that the germline-encoded diversity differs across antibody (recombined V-D-J) repertoires of various B cell subsets. To facilitate future comparative immunogenomics analysis, we created VDJgermlines, an R package that contains the germline sequences from multiple species. Our study informs strategies for the humanization and engineering of therapeutic antibodies.

Complete characterization of the human immune cell transcriptome using accurate full-length cDNA sequencing

The human immune system relies on highly complex and diverse transcripts and the proteins they encode. These include transcripts encoding human leukocyte antigen (HLA) receptors as well as B cell and T cell receptors (BCR and TCR). Determining which alleles an individual possesses for each HLA gene (high-resolution HLA typing) is essential to establish donor–recipient compatibility in organ and bone marrow transplantations. In turn, the repertoires of millions of unique BCR and TCR transcripts in each individual carry a vast amount of health-relevant information. Both short-read RNA-seq-based HLA typing and BCR/TCR repertoire sequencing (AIRR-seq) currently rely on our incomplete knowledge of the genetic diversity at HLA and BCR/TCR loci. Here, we generated over 10,000,000 full-length cDNA sequences at a median accuracy of 97.9% using our nanopore sequencing-based Rolling Circle Amplification to Concatemeric Consensus (R2C2) protocol. We used this data set to (1) show that deep and accurate full-length cDNA sequencing can be used to provide isoform-level transcriptome analysis for more than 9000 loci, (2) generate accurate sequences of HLA alleles, and (3) extract detailed AIRR data for the analysis of the adaptive immune system. The HLA and AIRR analysis approaches we introduce here are untargeted and therefore do not require prior knowledge of the composition or genetic diversity of HLA and BCR/TCR loci.

Adaptive immune receptor repertoires, an overview of this exciting field

The adaptive immune response in jawed vertebrates relies on the huge diversity and specificity of the B cell and T cell antigen receptors, the immunoglobulins (IG) or antibodies and the T cell receptors (TR), respectively. The high level of diversity has represented a barrier to a comprehensive analysis of the adaptive immune response before the emergence of high-throughput sequencing (HTS) technologies. The size and complexity of HTS data requires the generation of novel computational and analytical approaches, which are transforming how the adaptive immune responses are deciphered to understand the clonal dynamics and properties of antigen-specific B and T cells in response to different kind of antigens. This exciting and rapidly evolving field is not only impacting human and clinical immunology but also comparative immunology. We are now closer to understanding the evolution of adaptive immune response in jawed vertebrates. This review provides an overview about classical and current strategies developed to assess the IG/TR diversity and their applications in basic and clinical immunology.