The inference of phased haplotypes for the immunoglobulin H chain V region gene loci by analysis of VDJ gene rearrangements

High-Throughput DNA Sequencing Analysis of Antibody Repertoires

New high-throughput DNA sequencing (HTS) technologies developed in the past decade have begun to be applied to the study of the complex gene rearrangements that encode human antibodies. This article first reviews the genetic features of Ig loci and the HTS technologies that have been applied to human repertoire studies, then discusses key choices for experimental design and data analysis in these experiments and the insights gained in immunological and infectious disease studies with the use of these approaches.

Genomic variation in the porcine immunoglobulin lambda variable region

Production of a vast antibody repertoire is essential for the protection against pathogens. Variable region germline complexity contributes to repertoire diversity and is a standard feature of mammalian immunoglobulin loci, but functional V region genes are limited in swine. For example, the porcine lambda light chain locus is composed of 23 variable (V) genes and 4 joining (J) genes, but only 10 or 11 V and 2 J genes are functional. Allelic variation in V and J may increase overall diversity within a population, yet lead to repertoire holes in individuals lacking key alleles. Previous studies focused on heavy chain genetic variation, thus light chain allelic diversity is not known. We characterized allelic variation of the porcine immunoglobulin lambda variable (IGLV) region genes. All intact IGLV genes in 81 pigs were amplified, sequenced, and analyzed to determine their allelic variation and functionality. We observed mutational variation across the entire length of the IGLV genes, in both framework and complementarity determining regions (CDRs). Three recombination hotspot motifs were also identified suggesting that non-allelic homologous recombination is an evolutionarily alternative mechanism for generating germline antibody diversity. Functional alleles were greatest in the most highly expressed families, IGLV3 and IGLV8. At the population level, allelic variation appears to help maintain the potential for broad antibody repertoire diversity in spite of reduced gene segment choices and limited germline sequence modification. The trade-off may be a reduction in repertoire diversity within individuals that could result in an increased variation in immunity to infectious disease and response to vaccination.

DJ Pairing during VDJ Recombination Shows Positional Biases That Vary among Individuals with Differing IGHD Locus Immunogenotypes

Human IgH diversity is influenced by biases in the pairing of IGHD and IGHJ genes, but these biases have not been described in detail. We used high-throughput sequencing of VDJ rearrangements to explore DJ pairing biases in 29 individuals. It was possible to infer three contrasting IGHD-IGHJ haplotypes in nine of these individuals, and two of these haplotypes include deletion polymorphisms involving multiple contiguous IGHD genes. Therefore, we were able to explore how the underlying genetic makeup of the H chain locus influences the formation of particular DJ pairs. Analysis of nonproductive rearrangements demonstrates that 3' IGHD genes tend to pair preferentially with 5' IGHJ genes, whereas 5' IGHD genes pair preferentially with 3' IGHJ genes; the relationship between IGHD gene pairing frequencies and IGHD gene position is a near linear one for each IGHJ gene. However, striking differences are seen in individuals who carry deletion polymorphisms in the D locus. The absence of different blocks of IGHD genes leads to increases in the utilization frequencies of just a handful of genes, and these genes have no clear positional relationships to the deleted genes. This suggests that pairing frequencies may be influenced by additional complex positional relationships that perhaps arise from chromatin structure. In contrast to IGHD gene usage, IGHJ gene usage is unaffected by the IGHD gene-deletion polymorphisms. Such an outcome would be expected if the recombinase complex associates with an IGHJ gene before associating with an IGHD gene partner.

Practical guidelines for B-cell receptor repertoire sequencing analysis

High-throughput sequencing of B-cell immunoglobulin repertoires is increasingly being applied to gain insights into the adaptive immune response in healthy individuals and in those with a wide range of diseases. Recent applications include the study of autoimmunity, infection, allergy, cancer and aging. As sequencing technologies continue to improve, these repertoire sequencing experiments are producing ever larger datasets, with tens- to hundreds-of-millions of sequences. These data require specialized bioinformatics pipelines to be analyzed effectively. Numerous methods and tools have been developed to handle different steps of the analysis, and integrated software suites have recently been made available. However, the field has yet to converge on a standard pipeline for data processing and analysis. Common file formats for data sharing are also lacking. Here we provide a set of practical guidelines for B-cell receptor repertoire sequencing analysis, starting from raw sequencing reads and proceeding through pre-processing, determination of population structure, and analysis of repertoire properties. These include methods for unique molecular identifiers and sequencing error correction, V(D)J assignment and detection of novel alleles, clonal assignment, lineage tree construction, somatic hypermutation modeling, selection analysis, and analysis of stereotyped or convergent responses. The guidelines presented here highlight the major steps involved in the analysis of B-cell repertoire sequencing data, along with recommendations on how to avoid common pitfalls.

Soma-to-germline feedback is implied by the extreme polymorphism at IGHV relative to MHC: The manifest polymorphism of the MHC appears greatly exceeded at Immunoglobulin loci, suggesting antigen-selected somatic V mutants penetrate Weismann's Barrier

Soma-to-germline feedback is forbidden under the neo-Darwinian paradigm. Nevertheless, there is a growing realization it occurs frequently in immunoglobulin (Ig) variable (V) region genes. This is a surprising development. It arises from a most unlikely source in light of the exposure of co-author EJS to the haplotype data of RL Dawkins and others on the polymorphism of the Major Histocompatibility Complex, which is generally assumed to be the most polymorphic region in the genome (spanning ∼4 Mb). The comparison between the magnitude of MHC polymorphism with estimates for the human heavy chain immunoglobulin V locus (spanning ∼1 Mb), suggests IGHV could be many orders of magnitude more polymorphic than the MHC. This conclusion needs airing in the literature as it implies generational churn and soma-to-germline gene feedback. Pedigree-based experimental strategies to resolve the IGHV issue are outlined.

Rheumatoid arthritis response to treatment across IgG1 allotype - anti-TNF incompatibility: a case-only study

Introduction

We have hypothesized that incompatibility between the G1m genotype of the patient and the G1m1 and G1m17 allotypes carried by infliximab (INX) and adalimumab (ADM) could decrease the efficacy of these anti-tumor necrosis factor (anti-TNF) antibodies in the treatment of rheumatoid arthritis (RA).

Methods

The G1m genotypes were analyzed in three collections of patients with RA totaling 1037 subjects. The first, used for discovery, comprised 215 Spanish patients. The second and third were successively used for replication. They included 429 British and Greek patients and 393 Spanish and British patients, respectively. Two outcomes were considered: change in the Disease Activity Score in 28 joint (ΔDAS28) and the European League Against Rheumatism (EULAR) response criteria.

Results

An association between less response to INX and incompatibility of the G1m1,17 allotype was found in the discovery collection at 6 months of treatment (P = 0.03). This association was confirmed in the replications (P = 0.02 and 0.08, respectively) leading to a global association (P = 0.001) that involved a mean difference in ΔDAS28 of 0.4 units between compatible and incompatible patients (2.3 ± 1.5 in compatible patients vs. 1.9 ± 1.5 in incompatible patients) and an increase in responders and decrease in non-responders according to the EULAR criteria (P = 0.03). A similar association was suggested for patients treated with ADM in the discovery collection, but it was not supported by replication.

Conclusions

Our results suggest that G1m1,17 allotypes are associated with response to INX and could aid improved therapeutic targeting in RA.

Electronic supplementary material

The online version of this article (doi:10.1186/s13075-015-0571-z) contains supplementary material, which is available to authorized users.

Automated analysis of high-throughput B-cell sequencing data reveals a high frequency of novel immunoglobulin V gene segment alleles

Individual variation in germline and expressed B-cell immunoglobulin (Ig) repertoires has been associated with aging, disease susceptibility, and differential response to infection and vaccination. Repertoire properties can now be studied at large-scale through next-generation sequencing of rearranged Ig genes. Accurate analysis of these repertoire-sequencing (Rep-Seq) data requires identifying the germline variable (V), diversity (D), and joining (J) gene segments used by each Ig sequence. Current V(D)J assignment methods work by aligning sequences to a database of known germline V(D)J segment alleles. However, existing databases are likely to be incomplete and novel polymorphisms are hard to differentiate from the frequent occurrence of somatic hypermutations in Ig sequences. Here we develop a Tool for Ig Genotype Elucidation via Rep-Seq (TIgGER). TIgGER analyzes mutation patterns in Rep-Seq data to identify novel V segment alleles, and also constructs a personalized germline database containing the specific set of alleles carried by a subject. This information is then used to improve the initial V segment assignments from existing tools, like IMGT/HighV-QUEST. The application of TIgGER to Rep-Seq data from seven subjects identified 11 novel V segment alleles, including at least one in every subject examined. These novel alleles constituted 13% of the total number of unique alleles in these subjects, and impacted 3% of V(D)J segment assignments. These results reinforce the highly polymorphic nature of human Ig V genes, and suggest that many novel alleles remain to be discovered. The integration of TIgGER into Rep-Seq processing pipelines will increase the accuracy of V segment assignments, thus improving B-cell repertoire analyses.

Sequencing of the human IG light chain loci from a hydatidiform mole BAC library reveals locus-specific signatures of genetic diversity

Germline variation at immunoglobulin (IG) loci is critical for pathogen-mediated immunity, but establishing complete haplotype sequences in these regions has been problematic because of complex sequence architecture and diploid source DNA. We sequenced BAC clones from the effectively haploid human hydatidiform mole cell line, CHM1htert, across the light chain IG loci, kappa (IGK) and lambda (IGL), creating single haplotype representations of these regions. The IGL haplotype generated here is 1.25 Mb of contiguous sequence, including four novel IGLV alleles, one novel IGLC allele, and an 11.9-kb insertion. The CH17 IGK haplotype consists of two 644 kb proximal and 466 kb distal contigs separated by a large gap of unknown size; these assemblies added 49 kb of unique sequence extending into this gap. Our analysis also resulted in the characterization of seven novel IGKV alleles and a 16.7-kb region exhibiting signatures of interlocus sequence exchange between distal and proximal IGKV gene clusters. Genetic diversity in IGK/IGL was compared with that of the IG heavy chain (IGH) locus within the same haploid genome, revealing threefold (IGK) and sixfold (IGL) higher diversity in the IGH locus, potentially associated with increased levels of segmental duplication and the telomeric location of IGH.

pRESTO: a toolkit for processing high-throughput sequencing raw reads of lymphocyte receptor repertoires

Driven by dramatic technological improvements, large-scale characterization of lymphocyte receptor repertoires via high-throughput sequencing is now feasible. Although promising, the high germline and somatic diversity, especially of B-cell immunoglobulin repertoires, presents challenges for analysis requiring the development of specialized computational pipelines. We developed the REpertoire Sequencing TOolkit (pRESTO) for processing reads from high-throughput lymphocyte receptor studies. pRESTO processes raw sequences to produce error-corrected, sorted and annotated sequence sets, along with a wealth of metrics at each step. The toolkit supports multiplexed primer pools, single- or paired-end reads and emerging technologies that use single-molecule identifiers. pRESTO has been tested on data generated from Roche and Illumina platforms. It has a built-in capacity to parallelize the work between available processors and is able to efficiently process millions of sequences generated by typical high-throughput projects.

AVAILABILITY AND IMPLEMENTATION

pRESTO is freely available for academic use. The software package and detailed tutorials may be downloaded from http://clip.med.yale.edu/presto.

The promise and challenge of high-throughput sequencing of the antibody repertoire

Georgiou and colleagues discuss rapidly evolving methods for high-throughput sequencing of the antibody repertoire, and how the resulting data may be applied to answer basic and translational research questions.

Efforts to determine the antibody repertoire encoded by B cells in the blood or lymphoid organs using high-throughput DNA sequencing technologies have been advancing at an extremely rapid pace and are transforming our understanding of humoral immune responses. Information gained from high-throughput DNA sequencing of immunoglobulin genes (Ig-seq) can be applied to detect B-cell malignancies with high sensitivity, to discover antibodies specific for antigens of interest, to guide vaccine development and to understand autoimmunity. Rapid progress in the development of experimental protocols and informatics analysis tools is helping to reduce sequencing artifacts, to achieve more precise quantification of clonal diversity and to extract the most pertinent biological information. That said, broader application of Ig-seq, especially in clinical settings, will require the development of a standardized experimental design framework that will enable the sharing and meta-analysis of sequencing data generated by different laboratories.

Convergent antibody signatures in human dengue

Dengue is the most prevalent mosquito-borne viral disease in humans, and the lack of early prognostics, vaccines, and therapeutics contributes to immense disease burden. To identify patterns that could be used for sequence-based monitoring of the antibody response to dengue, we examined antibody heavy-chain gene rearrangements in longitudinal peripheral blood samples from 60 dengue patients. Comparing signatures between acute dengue, postrecovery, and healthy samples, we found increased expansion of B cell clones in acute dengue patients, with higher overall clonality in secondary infection. Additionally, we observed consistent antibody sequence features in acute dengue in the highly variable major antigen-binding determinant, complementarity-determining region 3 (CDR3), with specific CDR3 sequences highly enriched in acute samples compared to postrecovery, healthy, or non-dengue samples. Dengue thus provides a striking example of a human viral infection where convergent immune signatures can be identified in multiple individuals. Such signatures could facilitate surveillance of immunological memory in communities.

Deep sequencing and Circos analyses of antibody libraries reveal antigen-driven selection of Ig VH genes during HIV-1 infection

The vast diversity of antibody repertoires is largely attributed to heavy chain (V(H)) recombination of variable (V), diversity (D) and joining (J) gene segments. We used 454 sequencing information of the variable domains of the antibody heavy chain repertoires from neonates, normal adults and an HIV-1-infected individual, to analyze, with Circos software, the VDJ pairing patterns at birth, adulthood and a time-dependent response to HIV-1 infection. Our comparative analyses of the Ig VDJ repertoires from these libraries indicated that, from birth to adulthood, VDJ recombination patterns remain the same with some slight changes, whereas some V(H) families are selected and preferentially expressed after long-term infection with HIV-1. We also demonstrated that the immune system responds to HIV-1 chronic infection by selectively expanding certain HV families in an attempt to combat infection. Our findings may have implications for understanding immune responses in pathology as well as for development of new therapeutics and vaccines.

The human IgE repertoire

IgE is a key mediator in allergic diseases. However, in strong contrast to other antibody isotypes, many details of the composition of the human IgE repertoire are poorly defined. The low levels of human IgE in the circulation and the rarity of IgE-producing B cells are important reasons for this lack of knowledge. In this review, we summarize the current knowledge on these repertoires both in terms of their complexity and activity, i.e. knowledge which despite the difficulties encountered when studying the molecular details of human IgE has been acquired in recent years. We also take a look at likely future developments, for instance through improvements in sequencing technology and methodology that allow the isolation of additional allergen-specific human antibodies mimicking IgE, as this certainly will support our understanding of human IgE in the context of human disease in the years to come.

The shape of the lymphocyte receptor repertoire: lessons from the B cell receptor

Both the B cell receptor (BCR) and the T cell receptor (TCR) repertoires are generated through essentially identical processes of V(D)J recombination, exonuclease trimming of germline genes, and the random addition of non-template encoded nucleotides. The naïve TCR repertoire is constrained by thymic selection, and TCR repertoire studies have therefore focused strongly on the diversity of MHC-binding complementarity determining region (CDR) CDR3. The process of somatic point mutations has given B cell studies a major focus on variable (IGHV, IGLV, and IGKV) genes. This in turn has influenced how both the naïve and memory BCR repertoires have been studied. Diversity (D) genes are also more easily identified in BCR VDJ rearrangements than in TCR VDJ rearrangements, and this has allowed the processes and elements that contribute to the incredible diversity of the immunoglobulin heavy chain CDR3 to be analyzed in detail. This diversity can be contrasted with that of the light chain where a small number of polypeptide sequences dominate the repertoire. Biases in the use of different germline genes, in gene processing, and in the addition of non-template encoded nucleotides appear to be intrinsic to the recombination process, imparting "shape" to the repertoire of rearranged genes as a result of differences spanning many orders of magnitude in the probabilities that different BCRs will be generated. This may function to increase the precursor frequency of naïve B cells with important specificities, and the likely emergence of such B cell lineages upon antigen exposure is discussed with reference to public and private T cell clonotypes.

Complete haplotype sequence of the human immunoglobulin heavy-chain variable, diversity, and joining genes and characterization of allelic and copy-number variation

The immunoglobulin heavy-chain locus (IGH) encodes variable (IGHV), diversity (IGHD), joining (IGHJ), and constant (IGHC) genes and is responsible for antibody heavy-chain biosynthesis, which is vital to the adaptive immune response. Programmed V-(D)-J somatic rearrangement and the complex duplicated nature of the locus have impeded attempts to reconcile its genomic organization based on traditional B-lymphocyte derived genetic material. As a result, sequence descriptions of germline variation within IGHV are lacking, haplotype inference using traditional linkage disequilibrium methods has been difficult, and the human genome reference assembly is missing several expressed IGHV genes. By using a hydatidiform mole BAC clone resource, we present the most complete haplotype of IGHV, IGHD, and IGHJ gene regions derived from a single chromosome, representing an alternate assembly of ∼1 Mbp of high-quality finished sequence. From this we add 101 kbp of previously uncharacterized sequence, including functional IGHV genes, and characterize four large germline copy-number variants (CNVs). In addition to this germline reference, we identify and characterize eight CNV-containing haplotypes from a panel of nine diploid genomes of diverse ethnic origin, discovering previously unmapped IGHV genes and an additional 121 kbp of insertion sequence. We genotype four of these CNVs by using PCR in 425 individuals from nine human populations. We find that all four are highly polymorphic and show considerable evidence of stratification (Fst = 0.3-0.5), with the greatest differences observed between African and Asian populations. These CNVs exhibit weak linkage disequilibrium with SNPs from two commercial arrays in most of the populations tested.

Chromatin conformation governs T-cell receptor Jβ gene segment usage

T cells play fundamental roles in adaptive immunity, relying on a diverse repertoire of T-cell receptor (TCR) α and β chains. Diversity of the TCR β chain is generated in part by a random yet intrinsically biased combinatorial rearrangement of variable (Vβ), diversity (Dβ), and joining (Jβ) gene segments. The mechanisms that determine biases in gene segment use remain unclear. Here we show, using a high-throughput TCR sequencing approach, that a physical model of chromatin conformation at the DJβ genomic locus explains more than 80% of the biases in Jβ use that we measured in murine T cells. This model also predicts correctly how differences in intersegment genomic distances between humans and mice translate into differences in Jβ bias between TCR repertoires of these two species. As a consequence of these structural and other biases, TCR sequences are produced with different a priori frequencies, thus affecting their probability of becoming public TCRs that are shared among individuals. Surprisingly, we find that many more TCR sequences are shared among all five mice we studied than among only subgroups of three or four mice. We derive a necessary mathematical condition explaining this finding, which indicates that the TCR repertoire contains a core set of receptor sequences that are highly abundant among individuals, if their a priori probability of being produced by the recombination process is higher than a defined threshold. Our results provide evidence for an expanded role of chromatin conformation in VDJ rearrangement, from control of gene accessibility to precise determination of gene segment use.

The immunoglobulin heavy chain locus: genetic variation, missing data, and implications for human disease

The immunoglobulin (IG) loci consist of repeated and highly homologous sets of genes of different types, variable (V), diversity (D) and junction (J), that rearrange in developing B cells to produce an individual's highly variable repertoire of expressed antibodies, designed to bind to a vast array of pathogens. This repeated structure makes these loci susceptible to a high frequency of insertion and deletion events through evolutionary time, and also makes them difficult to characterize at the genomic level or assay with high-throughput techniques. Given the central role of antibodies in the adaptive immune system, it is not surprising that early candidate gene approaches showed that germline polymorphisms in these regions correlated with susceptibility to both infectious and autoimmune diseases. However, more recent studies, particularly those using high-throughput genome-wide arrays, have failed to implicate these loci in disease. In this review of the IG heavy chain variable gene cluster (IGHV), we examine how poorly we understand the distribution of haplotype variation in this genomic region, and we argue that this lack of information may mask candidate loci in the IGHV gene cluster as causative factors for infectious and autoimmune diseases.