Nomenclature of the human immunoglobulin heavy (IGH) genes

De novo antibody discovery in human blood from full-length single B cell transcriptomics and matching haplotyped-resolved germline assemblies

Immunoglobulin (IGH, IGK, IGL) loci in the human genome are highly polymorphic regions that encode the building blocks of the light and heavy chain IG proteins that dimerize to form antibodies. The processes of V(D)J recombination and somatic hypermutation in B cells are responsible for creating an enormous reservoir of highly specific antibodies capable of binding a vast array of possible antigens. However, the antibody repertoire is fundamentally limited by the set of variable (V), diversity (D), and joining (J) alleles present in the germline IG loci. To better understand how the germline IG haplotypes contribute to the expressed antibody repertoire, we combined genome sequencing of the germline IG loci with single-cell transcriptome sequencing of B cells from the same donor. Sequencing and assembly of the germline IG loci captured the IGH locus in a single fully-phased contig where the maternal and paternal contributions to the germline V, D, and J repertoire can be fully resolved. The B cells were collected following a measles, mumps, and rubella (MMR) vaccination, resulting in a population of cells that were activated in response to this specific immune challenge. Single-cell, full-length transcriptome sequencing of these B cells resulted in whole transcriptome characterization of each cell, as well as highly-accurate consensus sequences for the somatically rearranged and hypermutated light and heavy chain IG transcripts. A subset of antibodies synthesized based on their consensus heavy and light chain transcript sequences demonstrated binding to measles antigens and neutralization of measles live virus.

AIRR-C IG Reference Sets: curated sets of immunoglobulin heavy and light chain germline genes

Introduction

Analysis of an individual's immunoglobulin (IG) gene repertoire requires the use of high-quality germline gene reference sets. When sets only contain alleles supported by strong evidence, AIRR sequencing (AIRR-seq) data analysis is more accurate and studies of the evolution of IG genes, their allelic variants and the expressed immune repertoire is therefore facilitated.

Methods

The Adaptive Immune Receptor Repertoire Community (AIRR-C) IG Reference Sets have been developed by including only human IG heavy and light chain alleles that have been confirmed by evidence from multiple high-quality sources. To further improve AIRR-seq analysis, some alleles have been extended to deal with short 3' or 5' truncations that can lead them to be overlooked by alignment utilities. To avoid other challenges for analysis programs, exact paralogs (e.g. IGHV1-69*01 and IGHV1-69D*01) are only represented once in each set, though alternative sequence names are noted in accompanying metadata.

Results and discussion

The Reference Sets include less than half the previously recognised IG alleles (e.g. just 198 IGHV sequences), and also include a number of novel alleles: 8 IGHV alleles, 2 IGKV alleles and 5 IGLV alleles. Despite their smaller sizes, erroneous calls were eliminated, and excellent coverage was achieved when a set of repertoires comprising over 4 million V(D)J rearrangements from 99 individuals were analyzed using the Sets. The version-tracked AIRR-C IG Reference Sets are freely available at the OGRDB website (https://ogrdb.airr-community.org/germline_sets/Human) and will be regularly updated to include newly observed and previously reported sequences that can be confirmed by new high-quality data.

Antibody Sequence and Structure Analyses Using IMGT®: 30 Years of Immunoinformatics

IMGT^®, the international ImMunoGeneTics information system^®, http://www.imgt.org , the global reference in immunogenetics and immunoinformatics, was created in 1989 by Marie-Paule Lefranc (Université de Montpellier and CNRS) to manage the huge diversity of the antigen receptors, immunoglobulins (IG) or antibodies, and T cell receptors (TR) of the adaptive immune responses. The founding of IMGT^® marked the advent of immunoinformatics, which emerged at the interface between immunogenetics and bioinformatics. IMGT^® standardized analysis of the IG, TR, and major histocompatibility (MH) genes and proteins bridges the gap between sequences and three-dimensional (3D) structures, for all jawed vertebrates from fish to humans. This is achieved through the IMGT Scientific chart rules, based on the IMGT-ONTOLOGY axioms, and primarily CLASSIFICATION (IMGT gene and allele nomenclature) and NUMEROTATION (IMGT unique numbering and IMGT Colliers de Perles). IMGT^® comprises seven databases (IMGT/LIGM-DB for nucleotide sequences, IMGT/GENE-DB for genes and alleles, etc.), 17 tools (IMGT/V-QUEST, IMGT/JunctionAnalysis, IMGT/HighV-QUEST for NGS, etc.), and more than 20,000 Web resources. In this chapter, the focus is on the tools for amino acid sequences per domain (IMGT/DomainGapAlign and IMGT/Collier-de-Perles), and on the databases for receptors (IMGT/2Dstructure-DB and IMGT/3D-structure-DB) described per receptor, chain, and domain and, for 3D, with contact analysis, paratope, and epitope. The IMGT/mAb-DB is the query interface for monoclonal antibodies (mAb), fusion proteins for immune applications (FPIA), composite proteins for clinical applications (CPCA), and related proteins of interest (RPI) with links to IMGT^® 2D and 3D databases and to the World Health Organization (WHO) International Nonproprietary Names (INN) program lists. The chapter includes the human IG allotypes and antibody engineered variants for effector properties used in the description of therapeutical mAb.

IMGT ® Nomenclature of Engineered IGHG Variants Involved in Antibody Effector Properties and Formats

The constant region of the immunoglobulin (IG) or antibody heavy gamma chain is frequently engineered to modify the effector properties of the therapeutic monoclonal antibodies. These variants are classified in regards to their effects on effector functions, antibody-dependent cytotoxicity (ADCC), antibody-dependent phagocytosis (ADCP), complement-dependent cytotoxicity (CDC) enhancement or reduction, B cell inhibition by the coengagement of antigen and FcγR on the same cell, on half-life increase, and/or on structure such as prevention of IgG4 half-IG exchange, hexamerisation, knobs-into-holes and the heteropairing H-H of bispecific antibodies, absence of disulfide bridge inter H-L, absence of glycosylation site, and site-specific drug attachment engineered cysteine. The IMGT engineered variant identifier is comprised of the species and gene name (and eventually allele), the letter 'v' followed by a number (assigned chronologically), and for each concerned domain (e.g, CH1, h, CH2 and CH3), the novel AA (single letter abbreviation) and IMGT position according to the IMGT unique numbering for the C-domain and between parentheses, the Eu numbering. IMGT engineered variants are described with detailed amino acid changes, visualized in motifs based on the IMGT numbering bridging genes, sequences, and structures for higher order description.

Cross-reactive antibodies targeting surface-exposed non-structural protein 1 (NS1) of dengue virus-infected cells recognize epitopes on the spaghetti loop of the β-ladder domain

Non-structural protein 1 (NS1) is a glycoprotein component of dengue virus (DENV) that is essential for viral replication, infection and immune evasion. Immunization with NS1 has been shown to elicit antibody-mediated immune responses which protect mice against DENV infections. Here, we obtained peripheral blood mononuclear cells from human subjects with secondary dengue infections, which were used to construct a dengue immune phage library displaying single-chain variable fragments. Phage selective for DENV NS1 were obtained by biopanning. Twenty-one monoclonal antibodies (mAbs) against DENV NS1 were generated from the selected phage and characterized in detail. We found most anti-NS1 mAbs used IGHV1 heavy chain antibody genes. The mAbs were classified into strongly and weakly-reactive groups based on their binding to NS1 expressed in dengue virus 2 (DENV2)-infected cells. Antibody binding experiments with recombinant NS1 proteins revealed that the mAbs recognize conformational epitopes on the β-ladder domain (amino acid residues 178–273) of DENV NS1. Epitope mapping studies on alanine-substituted NS1 proteins identified distinct but overlapping epitopes. Protruding amino acids distributed around the spaghetti loop are required for the binding of the strongly-reactive mAbs, whereas the recognition residues of the weakly-reactive mAbs are likely to be located in inaccessible sites facing toward the cell membrane. This information could guide the design of an NS1 epitope-based vaccine that targets cross-reactive conserved epitopes on cell surface-associated DENV NS1.

Efficient human-like antibody repertoire and hybridoma production in trans-chromosomic mice carrying megabase-sized human immunoglobulin loci

Trans-chromosomic (Tc) mice carrying mini-chromosomes with megabase-sized human immunoglobulin (Ig) loci have contributed to the development of fully human therapeutic monoclonal antibodies, but mitotic instability of human mini-chromosomes in mice may limit the efficiency of hybridoma production. Here, we establish human antibody-producing Tc mice (TC-mAb mice) that stably maintain a mouse-derived, engineered chromosome containing the entire human Ig heavy and kappa chain loci in a mouse Ig-knockout background. Comprehensive, high-throughput DNA sequencing shows that the human Ig repertoire, including variable gene usage, is well recapitulated in TC-mAb mice. Despite slightly altered B cell development and a delayed immune response, TC-mAb mice have more subsets of antigen-specific plasmablast and plasma cells than wild-type mice, leading to efficient hybridoma production. Our results thus suggest that TC-mAb mice offer a valuable platform for obtaining fully human therapeutic antibodies, and a useful model for elucidating the regulation of human Ig repertoire formation.

Trans-chromosomic (Tc) mice have helped the development of therapeutic antibodies, but chromosome instability limits its application. Here the authors develop a new line of Tc mice with full human Ig heavy and kappa loci integrated into the mouse artificial chromosome for stable passage, and confirm efficient generation of B cell responses and specific antibodies.

A novel digital PCR-based method to quantify (switched) B cells reveals the extent of allelic involvement in different recombination processes in the IGH locus

B cells fulfill an important role in the adaptive immunity. Upon activation and immunoglobulin (IG) class switching, these cells function in the humoral immunity compartment as plasma cells. For clinical applications, it can be important to quantify (switched) B cells accurately in a variety of body fluids and tissues of benign, inflammatory and malignant origin. For decades, flow cytometry and immunohistochemistry (IHC) have been the preferred methods for quantification. Although these methods are widely used, both depend on the accessibility of B cell epitopes and therefore require intact (fixed) cells. Whenever samples are low in quantity and/or quality, accurate quantification can be difficult. By shifting the focus from epitopes to DNA markers, quantification of B cells remains achievable. During differentiation and maturation, B cells are subjected to programmed genetic recombination processes like VDJ rearrangements and class switch recombination (CSR), which result in deletion of specific sequences of the IGH locus. These cell type-specific DNA "scars" (loss of sequences) in IG genes can be exploited as B cell markers in digital PCR (dPCR) based quantification methods. Here, we describe a novel, specific and sensitive digital PCR-based method to quantify mature and switched B cells in DNA specimens of benign and (copy number unstable) malignant origin. We compared this novel way of B cell quantitation with flow cytometric and immunohistochemical methods. Through cross-validation with flow cytometric sorted B cell subpopulations, we gained quantitative insights into allelic involvement in different recombination processes in the IGH locus. Our newly developed method is accurate and independent of the cellular context, offering new possibilities for quantification, even for (limited) small samples like liquid biopsies.

IMGT® Biocuration and Analysis of the Rhesus Monkey IG Loci

The adaptive immune system, along with the innate immune system, are the two main biological processes that protect an organism from pathogens. The adaptive immune system is characterized by the specificity and extreme diversity of its antigen receptors. These antigen receptors are the immunoglobulins (IG) or antibodies of the B cells and the T cell receptors (TR) of the T cells. The IG are proteins that have a dual role in immunity: they recognize antigens and trigger elimination mechanisms, to rid the body of foreign cells. The synthesis of the immunoglobulin heavy and light chains requires gene rearrangements at the DNA level in the IGH, IGK, and IGL loci. The rhesus monkey (Macaca mulatta) is one of the most widely used nonhuman primate species in biomedical research. In this manuscript, we provide a thorough analysis of the three IG loci of the Mmul_10 assembly of rhesus monkey, integrating IMGT previously existing data. Detailed characterization of IG genes includes their localization and position in the loci, the determination of the allele functionality, and the description of the regulatory elements of their promoters as well as the sequences of the conventional recombination signals (RS). This complete annotation of the genomic IG loci of Mmul_10 assembly and the highly detailed IG gene characterization could be used as a model, in additional rhesus monkey assemblies, for the analysis of the IG allelic polymorphism and structural variation, which have been described in rhesus monkeys.

DNA End Joining: G0-ing to the Core

Humans have evolved a series of DNA double-strand break (DSB) repair pathways to efficiently and accurately rejoin nascently formed pairs of double-stranded DNA ends (DSEs). In G0/G1-phase cells, non-homologous end joining (NHEJ) and alternative end joining (A-EJ) operate to support covalent rejoining of DSEs. While NHEJ is predominantly utilized and collaborates extensively with the DNA damage response (DDR) to support pairing of DSEs, much less is known about A-EJ collaboration with DDR factors when NHEJ is absent. Non-cycling lymphocyte progenitor cells use NHEJ to complete V(D)J recombination of antigen receptor genes, initiated by the RAG1/2 endonuclease which holds its pair of targeted DSBs in a synapse until each specified pair of DSEs is handed off to the NHEJ DSB sensor complex, Ku. Similar to designer endonuclease DSBs, the absence of Ku allows for A-EJ to access RAG1/2 DSEs but with random pairing to complete their repair. Here, we describe recent insights into the major phases of DSB end joining, with an emphasis on synapsis and tethering mechanisms, and bring together new and old concepts of NHEJ vs. A-EJ and on RAG2-mediated repair pathway choice.

Immunoglobulins or Antibodies: IMGT® Bridging Genes, Structures and Functions

IMGT®, the international ImMunoGeneTics® information system founded in 1989 by Marie-Paule Lefranc (Université de Montpellier and CNRS), marked the advent of immunoinformatics, a new science at the interface between immunogenetics and bioinformatics. For the first time, the immunoglobulin (IG) or antibody and T cell receptor (TR) genes were officially recognized as 'genes' as well as were conventional genes. This major breakthrough has allowed the entry, in genomic databases, of the IG and TR variable (V), diversity (D) and joining (J) genes and alleles of Homo sapiens and of other jawed vertebrate species, based on the CLASSIFICATION axiom. The second major breakthrough has been the IMGT unique numbering and the IMGT Collier de Perles for the V and constant (C) domains of the IG and TR and other proteins of the IG superfamily (IgSF), based on the NUMEROTATION axiom. IMGT-ONTOLOGY axioms and concepts bridge genes, sequences, structures and functions, between biological and computational spheres in the IMGT® system (Web resources, databases and tools). They provide the IMGT Scientific chart rules to identify, to describe and to analyse the IG complex molecular data, the huge diversity of repertoires, the genetic (alleles, allotypes, CNV) polymorphisms, the IG dual function (paratope/epitope, effector properties), the antibody humanization and engineering.

Genotyping and Copy Number Analysis of Immunoglobin Heavy Chain Variable Genes Using Long Reads

One of the remaining challenges to describing an individual's genetic variation lies in the highly heterogeneous and complex genomic regions that impede the use of classical reference-guided mapping and assembly approaches. Once such region is the Immunoglobulin heavy chain locus (IGH), which is critical for the development of antibodies and the adaptive immune system. We describe ImmunoTyper, the first PacBio-based genotyping and copy number calling tool specifically designed for IGH V genes (IGHV). We demonstrate that ImmunoTyper's multi-stage clustering and combinatorial optimization approach represents the most comprehensive IGHV genotyping approach published to date, through validation using gold-standard IGH reference sequence. This preliminary work establishes the feasibility of fine-grained genotype and copy number analysis using error-prone long reads in complex multi-gene loci and opens the door for in-depth investigation into IGHV heterogeneity using accessible and increasingly common whole-genome sequence.

Predicting B cell receptor substitution profiles using public repertoire data

B cells develop high affinity receptors during the course of affinity maturation, a cyclic process of mutation and selection. At the end of affinity maturation, a number of cells sharing the same ancestor (i.e. in the same “clonal family”) are released from the germinal center; their amino acid frequency profile reflects the allowed and disallowed substitutions at each position. These clonal-family-specific frequency profiles, called “substitution profiles”, are useful for studying the course of affinity maturation as well as for antibody engineering purposes. However, most often only a single sequence is recovered from each clonal family in a sequencing experiment, making it impossible to construct a clonal-family-specific substitution profile. Given the public release of many high-quality large B cell receptor datasets, one may ask whether it is possible to use such data in a prediction model for clonal-family-specific substitution profiles. In this paper, we present the method “Substitution Profiles Using Related Families” (SPURF), a penalized tensor regression framework that integrates information from a rich assemblage of datasets to predict the clonal-family-specific substitution profile for any single input sequence. Using this framework, we show that substitution profiles from similar clonal families can be leveraged together with simulated substitution profiles and germline gene sequence information to improve prediction. We fit this model on a large public dataset and validate the robustness of our approach on two external datasets. Furthermore, we provide a command-line tool in an open-source software package (https://github.com/krdav/SPURF) implementing these ideas and providing easy prediction using our pre-fit models.

Antibody engineering can be greatly informed by knowledge about the underlying affinity maturation process. As such this can be probed by sequencing, but unfortunately, in practice often only one member of the clonal family is sequenced, making it difficult to determine a set of possible amino acid mutations that would retain the original antibody antigen binding affinity. We overcome this data sparsity by developing a statistical learning approach that leverages vast information about amino acid preferences available in public immune system repertoire data. We use a penalized regression approach to devise a flexible statistical model that integrates multiple sources of information into a coherent prediction framework and validate our prediction algorithm using subsampling and held out data.

Different Somatic Hypermutation Levels among Antibody Subclasses Disclosed by a New Next-Generation Sequencing-Based Antibody Repertoire Analysis

A diverse antibody repertoire is primarily generated by the rearrangement of V, D, and J genes and subsequent somatic hypermutation (SHM). Class-switch recombination (CSR) produces various isotypes and subclasses with different functional properties. Although antibody isotypes and subclasses are considered to be produced by both direct and sequential CSR, it is still not fully understood how SHMs accumulate during the process in which antibody subclasses are generated. Here, we developed a new next-generation sequencing (NGS)-based antibody repertoire analysis capable of identifying all antibody isotype and subclass genes and used it to examine the peripheral blood mononuclear cells of 12 healthy individuals. Using a total of 5,480,040 sequences, we compared percentage frequency of variable (V), junctional (J) sequence, and a combination of V and J, diversity, length, and amino acid compositions of CDR3, SHM, and shared clones in the IgM, IgD, IgG3, IgG1, IgG2, IgG4, IgA1, IgE, and IgA2 genes. The usage and diversity were similar among the immunoglobulin (Ig) subclasses. Clonally related sequences sharing identical V, D, J, and CDR3 amino acid sequences were frequently found within multiple Ig subclasses, especially between IgG1 and IgG2 or IgA1 and IgA2. SHM occurred most frequently in IgG4, while IgG3 genes were the least mutated among all IgG subclasses. The shared clones had almost the same SHM levels among Ig subclasses, while subclass-specific clones had different levels of SHM dependent on the genomic location. Given the sequential CSR, these results suggest that CSR occurs sequentially over multiple subclasses in the order corresponding to the genomic location of IGHCs, but CSR is likely to occur more quickly than SHMs accumulate within Ig genes under physiological conditions. NGS-based antibody repertoire analysis should provide critical information on how various antibodies are generated in the immune system.

Dysregulation of B Cell Repertoire Formation in Myasthenia Gravis Patients Revealed through Deep Sequencing

Myasthenia gravis (MG) is a prototypical B cell-mediated autoimmune disease affecting 20-50 people per 100,000. The majority of patients fall into two clinically distinguishable types based on whether they produce autoantibodies targeting the acetylcholine receptor (AChR-MG) or muscle specific kinase (MuSK-MG). The autoantibodies are pathogenic, but whether their generation is associated with broader defects in the B cell repertoire is unknown. To address this question, we performed deep sequencing of the BCR repertoire of AChR-MG, MuSK-MG, and healthy subjects to generate ∼518,000 unique V_H and V_L sequences from sorted naive and memory B cell populations. AChR-MG and MuSK-MG subjects displayed distinct gene segment usage biases in both V_H and V_L sequences within the naive and memory compartments. The memory compartment of AChR-MG was further characterized by reduced positive selection of somatic mutations in the V_H CDR and altered V_H CDR3 physicochemical properties. The V_L repertoire of MuSK-MG was specifically characterized by reduced V-J segment distance in recombined sequences, suggesting diminished V_L receptor editing during B cell development. Our results identify large-scale abnormalities in both the naive and memory B cell repertoires. Particular abnormalities were unique to either AChR-MG or MuSK-MG, indicating that the repertoires reflect the distinct properties of the subtypes. These repertoire abnormalities are consistent with previously observed defects in B cell tolerance checkpoints in MG, thereby offering additional insight regarding the impact of tolerance defects on peripheral autoimmune repertoires. These collective findings point toward a deformed B cell repertoire as a fundamental component of MG.

Envelope-specific B-cell populations in African green monkeys chronically infected with simian immunodeficiency virus

African green monkeys (AGMs) are natural primate hosts of simian immunodeficiency virus (SIV). Interestingly, features of the envelope-specific antibody responses in SIV-infected AGMs are distinct from that of HIV-infected humans and SIV-infected rhesus monkeys, including gp120-focused responses and rapid development of autologous neutralization. Yet, the lack of genetic tools to evaluate B-cell lineages hinders potential use of this unique non-human primate model for HIV vaccine development. Here we define features of the AGM Ig loci and compare the proportion of Env-specific memory B-cell populations to that of HIV-infected humans and SIV-infected rhesus monkeys. AGMs appear to have a higher proportion of Env-specific memory B cells that are mainly gp120 directed. Furthermore, AGM gp120-specific monoclonal antibodies display robust antibody-dependent cellular cytotoxicity and CD4-dependent virion capture activity. Our results support the use of AGMs to model induction of functional gp120-specific antibodies by HIV vaccine strategies.

Infection of African green monkeys with simian immunodeficiency virus is a potential model for HIV vaccine development. Here, Zhang et al. catalogue the immunoglobulin loci present in the genome of these animals, and experimentally study their B-cell response to the viral envelope protein.

Maximum-Entropy Models of Sequenced Immune Repertoires Predict Antigen-Antibody Affinity

The immune system has developed a number of distinct complex mechanisms to shape and control the antibody repertoire. One of these mechanisms, the affinity maturation process, works in an evolutionary-like fashion: after binding to a foreign molecule, the antibody-producing B-cells exhibit a high-frequency mutation rate in the genome region that codes for the antibody active site. Eventually, cells that produce antibodies with higher affinity for their cognate antigen are selected and clonally expanded. Here, we propose a new statistical approach based on maximum entropy modeling in which a scoring function related to the binding affinity of antibodies against a specific antigen is inferred from a sample of sequences of the immune repertoire of an individual. We use our inference strategy to infer a statistical model on a data set obtained by sequencing a fairly large portion of the immune repertoire of an HIV-1 infected patient. The Pearson correlation coefficient between our scoring function and the IC50 neutralization titer measured on 30 different antibodies of known sequence is as high as 0.77 (p-value 10-6), outperforming other sequence- and structure-based models.

Harnessing the protective potential of HIV-1 neutralizing antibodies

Recent biological, structural, and technical advances are converging within the HIV-1 vaccine field to harness the power of antibodies for prevention and therapy. Numerous monoclonal antibodies with broad neutralizing activity against diverse HIV-1 isolates have now been identified, revealing at least five sites of vulnerability on the envelope (Env) glycoproteins. While there are practical and technological barriers blocking a clear path from broadly neutralizing antibodies (bNAb) to a protective vaccine, this is not a dead end. Scientists are revisiting old approaches with new technology, cutting new trails through unexplored territory, and paving new roads in the hopes of preventing HIV-1 infection. Other promising avenues to capitalize on the power of bNAbs are also being pursued, such as passive antibody immunotherapy and gene therapy approaches. Moreover, non-neutralizing antibodies have inhibitory activities that could have protective potential, alone or in combination with bNAbs. With a new generation of bNAbs, and a clinical trial that associated antibodies with reduced acquisition, the field is closer than ever to developing strategies to use antibodies against HIV-1.

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.

V genes in primates from whole genome sequencing data

The adaptive immune system uses V genes for antigen recognition. However, the evolutionary diversification and selection processes within and across species and orders remain poorly understood. Here, we studied the amino acid (AA) sequences obtained from the translated in-frame V exons of immunoglobulins (IG) and T cell receptors (TR) from 16 primate species whose genomes have been sequenced. Multi-species comparative analysis supports the hypothesis that V genes in the IG loci undergo birth/death processes, thereby permitting rapid adaptability over evolutionary time. We also show that multiple cladistic groupings exist in the TRA (35 clades) and TRB (25 clades) V gene loci and that each primate species typically contributes at least one V gene to each of these clades. The results demonstrate that IG V genes and TR V genes have quite different evolutionary pathways; multiple duplications can explain the IG loci results, while coevolutionary pressures can explain the phylogenetic results of the TR V gene loci. Our results suggest that there exist evolutionary relationships between V gene clades in the TRA and TRB loci. Due to the long-standing preservation of these clades, such genes may have specific and necessary roles for the viability of a species.

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.

Determination of allergen specificity by heavy chains in grass pollen allergen-specific IgE antibodies

BACKGROUND

Affinity and clonality of allergen-specific IgE antibodies are important determinants for the magnitude of IgE-mediated allergic inflammation.

OBJECTIVE

We sought to analyze the contribution of heavy and light chains of human allergen-specific IgE antibodies for allergen specificity and to test whether promiscuous pairing of heavy and light chains with different allergen specificity allows binding and might affect affinity.

METHODS

Ten IgE Fabs specific for 3 non-cross-reactive major timothy grass pollen allergens (Phl p 1, Phl p 2, and Phl p 5) obtained by means of combinatorial cloning from patients with grass pollen allergy were used to construct stable recombinant single chain variable fragments (ScFvs) representing the original Fabs and shuffled ScFvs in which heavy chains were recombined with light chains from IgE Fabs with specificity for other allergens by using the pCANTAB 5 E expression system. Possible ancestor genes for the heavy chain and light chain variable region-encoding genes were determined by using sequence comparison with the ImMunoGeneTics database, and their chromosomal locations were determined. Recombinant ScFvs were tested for allergen specificity and epitope recognition by means of direct and sandwich ELISA, and affinity by using surface plasmon resonance experiments.

RESULTS

The shuffling experiments demonstrate that promiscuous pairing of heavy and light chains is possible and maintains allergen specificity, which is mainly determined by the heavy chains. ScFvs consisting of different heavy and light chains exhibited different affinities and even epitope specificity for the corresponding allergen.

CONCLUSION

Our results indicate that allergen specificity of allergen-specific IgE is mainly determined by the heavy chains. Different heavy and light chain pairings in allergen-specific IgE antibodies affect affinity and epitope specificity and thus might influence clinical reactivity to allergens.

The bovine genomic DNA sequence data reveal three IGHV subgroups, only one of which is functionally expressed

A comprehensive analysis of cattle shotgun sequencing data reveals 36 immunoglobulin heavy chain variable genes. The previously described bovine subgroup IGHV1 contains 10 functional genes with a conserved promoter including the consensus octamer and several other transcription factor binding sites, intact exons and matching cDNA sequences. Subgroups IGHV2 and IGHV3 consist entirely of pseudogenes. Thus, the bovine germline IGHV repertoire is very limited. The IGHV genes are distributed in mammalian clans I and II, while no clan III genes were detected. Clan-specific PCR of genomic DNA from cattle, sheep, Eurasian elk, white-tailed deer, pig and dolphin indicates highly dynamic evolution of IGHV gene usage within Cetartiodactyla. The bovine germline IGHV repertoire was probably generated by recent duplications of an IGHV1-IGHV2 homology unit. Immunoglobulin heavy chain genes are largely incorrectly assembled in the current cattle genome versions Btau_4.2 and UMD_3.1. FISH experiments confirm an IGHV locus close to terminus of BTA21.

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.

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

The existence of many highly similar genes in the lymphocyte receptor gene loci makes them difficult to investigate, and the determination of phased "haplotypes" has been particularly problematic. However, V(D)J gene rearrangements provide an opportunity to infer the association of Ig genes along the chromosomes. The chromosomal distribution of H chain genes in an Ig genotype can be inferred through analysis of VDJ rearrangements in individuals who are heterozygous at points within the IGH locus. We analyzed VDJ rearrangements from 44 individuals for whom sufficient unique rearrangements were available to allow comprehensive genotyping. Nine individuals were identified who were heterozygous at the IGHJ6 locus and for whom sufficient suitable VDJ rearrangements were available to allow comprehensive haplotyping. Each of the 18 resulting IGHV│IGHD│IGHJ haplotypes was unique. Apparent deletion polymorphisms were seen that involved as many as four contiguous, functional IGHV genes. Two deletion polymorphisms involving multiple contiguous IGHD genes were also inferred. Three previously unidentified gene duplications were detected, where two sequences recognized as allelic variants of a single gene were both inferred to be on a single chromosome. Phased genomic data brings clarity to the study of the contribution of each gene to the available repertoire of rearranged VDJ genes. Analysis of rearrangement frequencies suggests that particular genes may have substantially different yet predictable propensities for rearrangement within different haplotypes. Together with data highlighting the extent of haplotypic variation within the population, this suggests that there may be substantial variability in the available Ab repertoires of different individuals.

Clinical significance of productive immunoglobulin heavy chain gene rearrangements in childhood acute lymphoblastic leukemia

We analyzed the CDR3 region of 80 children with B-cell acute lymphoblastic leukemia (B-ALL) using the ImMunoGeneTics Information System and JOINSOLVER. In total, 108 IGH@ rearrangements were analyzed. Most of them (75.3%) were non-productive. IGHV@ segments proximal to IGHD-IGHJ@ were preferentially rearranged (45.3%). Increased utilization of IGHV3 segments IGHV3-13 (11.3%) and IGHV3-15 (9.3%), IGHD3 (30.5%), and IGHJ4 (34%) was noted. In pro-B ALL more frequent were IGHV3-11 (33.3%) and IGHV6-1 (33.3%), IGHD2-21 (50%), IGHJ4 (50%), and IGHJ6 (50%) segments. Shorter CDR3 length was observed in IGHV@6, IGHD7, and IGHJ1 segments, whereas increased CDR3 length was related to IGHV3, IGHD2, and IGHJ4 segments. Increased risk of relapse was found in patients with productive sequences. Specifically, the relapse-free survival rate at 5 years in patients with productive sequences at diagnosis was 75% (standard error [SE] ±9%), whereas in patients with non-productive sequences it was 97% (SE ±1.92%) (p-value =0.0264). Monoclonality and oligoclonality were identified in 81.2% and 18.75% cases at diagnosis, respectively. Sequence analysis revealed IGHV@ to IGHDJ joining only in 6.6% cases with oligoclonality. The majority (75%) of relapsed patients had monoclonal IGH@ rearrangements. The preferential utilization of IGHV@ segments proximal to IGHDJ depended on their location on the IGHV@ locus. Molecular mechanisms occurring during IGH@ rearrangement might play an essential role in childhood ALL prognosis. In our study, the productivity of the rearranged sequences at diagnosis proved to be a significant prognostic factor.

A preliminary analysis of the immunoglobulin genes in the African elephant (Loxodonta africana)

The genomic organization of the IgH (Immunoglobulin heavy chain), Igκ (Immunoglobulin kappa chain), and Igλ (Immunoglobulin lambda chain) loci in the African elephant (Loxodonta africana) was annotated using available genome data. The elephant IgH locus on scaffold 57 spans over 2,974 kb, and consists of at least 112 V(H) gene segments, 87 D(H) gene segments (the largest number in mammals examined so far), six J(H) gene segments, a single μ, a δ remnant, and eight γ genes (α and ε genes are missing, most likely due to sequence gaps). The Igκ locus, found on three scaffolds (202, 50 and 86), contains a total of 153 V(κ) gene segments, three J(κ) segments, and a single C(κ) gene. Two different transcriptional orientations were determined for these V(κ) gene segments. In contrast, the Igλ locus on scaffold 68 includes 15 V(λ) gene segments, all with the same transcriptional polarity as the downstream J(λ)-C(λ) cluster. These data suggest that the elephant immunoglobulin gene repertoire is highly diverse and complex. Our results provide insights into the immunoglobulin genes in a placental mammal that is evolutionarily distant from humans, mice, and domestic animals.

Tracing antigen signatures in the human IgE repertoire

Allergen recognition by IgE antibodies is a key event in allergic inflammation.

In this study, the IgE IGHV repertoires of individuals with allergy to the major birch pollen allergen, Bet v 1, were analyzed over a four years period of allergen exposure by RT-PCR and sequencing of cDNA.

Approximately half of the IgE transcripts represented non-redundant sequences, which belonged to seventeen different IGHV genes. Most variable regions contained somatic mutations but also non-mutated sequences were identified. There was no evidence for relevant increases of somatic mutations over time of allergen exposure. Highly similar IgE variable regions were found after four years of allergen exposure in the same and in genetically non-related individuals.

Our results indicate that allergens select and shape a limited number of similar IgE variable regions in the human IgE repertoire.

Individual variation in the germline Ig gene repertoire inferred from variable region gene rearrangements

Individual variation in the Ig germline gene repertoire leads to individual differences in the combinatorial diversity of the Ab repertoire, but the study of such variation has been problematic. The application of high-throughput DNA sequencing to the study of rearranged Ig genes now makes this possible. The sequencing of thousands of VDJ rearrangements from an individual, either from genomic DNA or expressed mRNA, should allow their germline IGHV, IGHD, and IGHJ repertoires to be inferred. In addition, where previously mere glimpses of diversity could be gained from sequencing studies, new large data sets should allow the rearrangement frequency of different genes and alleles to be seen with clarity. We analyzed the DNA of 108,210 human IgH chain rearrangements from 12 individuals and determined their individual IGH genotypes. The number of reportedly functional IGHV genes and allelic variants ranged from 45 to 60, principally because of variable levels of gene heterozygosity, and included 14 previously unreported IGHV polymorphisms. New polymorphisms of the IGHD3-16 and IGHJ6 genes were also seen. At heterozygous loci, remarkably different rearrangement frequencies were seen for the various IGHV alleles, and these frequencies were consistent between individuals. The specific alleles that make up an individual's Ig genotype may therefore be critical in shaping the combinatorial repertoire. The extent of genotypic variation between individuals is highlighted by an individual with aplastic anemia who appears to lack six contiguous IGHD genes on both chromosomes. These deletions significantly alter the potential expressed IGH repertoire, and possibly immune function, in this individual.

Chronic lymphocytic leukaemia: an immunobiology approach

B cell chronic lymphocytic leukaemia (CLL) is the most common adult leukaemia that follows an extremely variable clinical course. Several important prognostic parameters defining pathogenic and clinical subgroups of CLL have been identified and validated recently. The biological significance of immunoglobulin (Ig) heavy chain variable region gene (IgHV) mutational status and associated ZAP-70 over-expression, CD38 and chromosomal aberrations have enabled to identify patients at high risk for early disease progression and inferior survival. Moreover, studies of the B cell antigen receptor (BCR) structure and receptor signalling have been most helpful in revealing some new aspects of the biology of this disease. In particular, the analysis of IG genes has revealed that the expressed IgHV/IgKV/IgLV gene repertoires of CLL cells differ from those of normal B cells. A further unique feature of the CLL IG repertoire is the existence of subsets of cases with "stereotyped" BCRs. Accumulating molecular and phenotypic data support the notion that CLL development and evolution is not a simple scholastic event and strongly indicates a role for antigen in driving the cell of origin for at least some subsets of CLL cases.

Gross deletions involving IGHM, BTK, or Artemis: a model for genomic lesions mediated by transposable elements

Most genetic disruptions underlying human disease are microlesions, whereas gross lesions are rare with gross deletions being most frequently found (6%). Similar observations have been made in primary immunodeficiency genes, such as BTK, but for unknown reasons the IGHM and DCLRE1C (Artemis) gene defects frequently represent gross deletions ( approximately 60%). We characterized the gross deletion breakpoints in IGHM-, BTK-, and Artemis-deficient patients. The IGHM deletion breakpoints did not show involvement of recombination signal sequences or immunoglobulin switch regions. Instead, five IGHM, eight BTK, and five unique Artemis breakpoints were located in or near sequences derived from transposable elements (TE). The breakpoints of four out of five disrupted Artemis alleles were located in highly homologous regions, similar to Ig subclass deficiencies and Vh deletion polymorphisms. Nevertheless, these observations suggest a role for TEs in mediating gross deletions. The identified gross deletion breakpoints were mostly located in TE subclasses that were specifically overrepresented in the involved gene as compared to the average in the human genome. This concerned both long (LINE1) and short (Alu, MIR) interspersed elements, as well as LTR retrotransposons (ERV). Furthermore, a high total TE content (>40%) was associated with an increased frequency of gross deletions. Both findings were further investigated and confirmed in a total set of 20 genes disrupted in human disease. Thus, to our knowledge for the first time, we provide evidence that a high TE content, irrespective of the type of element, results in the increased incidence of gross deletions as gene disruption underlying human disease.

Evolutionary dynamics of the immunoglobulin heavy chain variable region genes in vertebrates

Immunoglobulin heavy chains are polypeptides encoded by four genes: variable (IGHV), joining (IGHJ), diversity (IGHD), and constant (IGHC) region genes. The number of IGHV genes varies from species to species. To understand the evolution of the IGHV multigene family, we identified and analyzed the IGHV sequences from 16 vertebrate species. The results show that the numbers of functional and nonfunctional IGHV genes among different species are positively correlated. The number of IGHV genes is relatively stable in teleosts, but the intragenomic sequence variation is generally higher in teleosts than in tetrapods. The IGHV genes in tetrapods can be classified into three phylogenetic clans (I, II, and III). The clan III and/or II genes are relatively abundant, whereas clan I genes exist in small numbers or are absent in most species. The genomic organization of clan I, II, and III IGHV genes varies considerably among species, but the entire IGHV locus seems to be conserved in the subtelomeric or near-centromeric region of chromosome. The presence or absence of specific IGHV clan members and the lineage-specific expansion and contraction of IGHV genes indicate that the IGHV locus continues to evolve in a species-specific manner. Our results suggest that the evolution of IGHV multigene family is more complex than previously thought and that several factors may act synergistically for the development of antibody repertoire.

Many human immunoglobulin heavy-chain IGHV gene polymorphisms have been reported in error

The identification of the genes that make up rearranged immunoglobulin genes is critical to many studies. For example, the enumeration of mutations in immunoglobulin genes is important for the prognosis of chronic lymphocytic leukemia, and this requires the accurate identification of the germline genes from which a particular sequence is derived. The immunoglobulin heavy-chain variable (IGHV) gene repertoire is generally considered to be highly polymorphic. In this report, we describe a bioinformatic analysis of germline and rearranged immunoglobulin gene sequences which casts doubt on the existence of a substantial proportion of reported germline polymorphisms. We report a five-level classification system for IGHV genes, which indicates the likelihood that the genes have been reported accurately. The classification scheme also reflects the likelihood that germline genes could be incorrectly identified in mutated VDJ rearrangements, because of similarities to other alleles. Of the 226 IGHV alleles that have previously been reported, our analysis suggests that 104 of these alleles almost certainly include sequence errors, and should be removed from the available repertoire. The analysis also highlights the presence of common mismatches, with respect to the germline, in many rearranged heavy-chain sequences, suggesting the existence of twelve previously unreported alleles. Sequencing of IGHV genes from six individuals in this study confirmed the existence of three of these alleles, which we designate IGHV3-49*04, IGHV3-49*05 and IGHV4-39*07. We therefore present a revised repertoire of expressed IGHV genes, which should substantially improve the accuracy of immunoglobulin gene analysis.

iHMMune-align: hidden Markov model-based alignment and identification of germline genes in rearranged immunoglobulin gene sequences

MOTIVATION

Immunoglobulin heavy chain (IGH) genes in mature B lymphocytes are the result of recombination of IGHV, IGHD and IGHJ germline genes, followed by somatic mutation. The correct identification of the germline genes that make up a variable VH domain is essential to our understanding of the process of antibody diversity generation as well as to clinical investigations of some leukaemias and lymphomas.

RESULTS

We have developed iHMMune-align, an alignment program that uses a hidden Markov model (HMM) to model the processes involved in human IGH gene rearrangement and maturation. The performance of iHMMune-align was compared to that of other immunoglobulin gene alignment utilities using both clonally related and randomly selected IGH sequences. This evaluation suggests that iHMMune-align provides a more accurate identification of component germline genes than other currently available IGH gene characterization programs.

AVAILABILITY

iHMMune-align cross-platform Java executable and web interface are freely available to academic users and can be accessed at http://www.emi.unsw.edu.au/~ihmmune/.

Use of IGHJ and IGHD gene mutations in analysis of immunoglobulin sequences for the prognosis of chronic lymphocytic leukemia

The level of somatic point mutation in immunoglobulin genes is an important prognostic indicator for patients with chronic lymphocytic leukemia (CLL). Mutation analysis presently focuses solely upon the heavy chain IGHV gene, however mutation is a stochastic process that also targets IGHD and IGHJ genes. Here, we evaluate the completeness and reliability of the reported IGHJ gene repertoire, and demonstrate the likely consequences of the inclusion of IGHD and IGHJ mutations in CLL analysis, using a dataset of 607 sequences. Inclusion of these mutations would lead to the re-classification of many sequences, which should significantly improve the prognostic value of mutation analysis.

The analysis and quantification of a clonal B cell response in a hyperimmunized anti-D donor

Healthy volunteers are hyperimmunized with RhD-positive red cells in order to obtain plasma containing high titres of anti-D immunoglobulin, which is used for the prevention of haemolytic disease of the fetus and newborn. We analysed the anti-D immune response in a donor who had been hyperimmunized for 7 years and who showed declining anti-D titres despite re-immunization. A phage display library representing the complete immunorepertoire and a second library representing the IGHV3 superspecies family genes (IGHV3s) repertoire in the donor were constructed and analysed. A clonal Ig-gene rearrangement was quantified in the peripheral blood by limiting dilution polymerase chain reaction (PCR) All RhD-binding phages from both libraries, except one, had heavy chains with IGH-VDJ rearrangements of the same clonal origin, but with different patterns of somatic mutations and joined with different light chains. Limiting dilution PCR performed on mRNA and genomic DNA showed a frequency of 1 clonal B cell in 2000 IgG1/3-positive B cells. We show the presence of clonally related RhD-specific B cells in a hyperimmunized anti-D donor who had declining anti-D titres and who was unresponsive to re-immunization. Furthermore, we found a high frequency of clonal B cells. These results contribute to the understanding of the immune response against RhD in hyperimmunized anti-D donors.

Immunoglobulin heavy- and light-chain repertoire in splenic marginal zone lymphoma

The considerable heterogeneity in morphology, immunophenotype, genotype, and clinical behavior of splenic marginal zone lymphoma (SMZL) hinders firm conclusions on the origin and differentiation stage of the neoplastic cells. Immunoglobulin (IG) gene usage and somatic mutation patterns were studied in a series of 43 SMZL cases. Clonal IGHV-D-J rearrangements were amplified in 42/43 cases (4 cases carried double rearrangements). Among IGHV-D-J rearrangements, IGHV3 and IGHV4 subgroup genes were used with the highest frequency. Nineteen IGHV genes were unmutated (> 98% homology to the closest germline IGHV gene), whereas 27/46 were mutated. Clonal IGKV-J and IGLV-J gene rearrangements were amplified in 36/43 cases, including 31 IGKV-J (8/31 in lambda light-chain expressing cases) and 12 IGLV-J rearrangements; 9/31 IGKV and 6/12 IGLV sequences were mutated. IGKV-J and IGLV-J rearrangements used 14 IGKV and 9 IGLV different germline genes. Significant evidence for positive selection by classical T-dependent antigen was found in only 5/27 IGHV and 6/15 IGKV+IGLV mutated genes. These results provide evidence for the diverse B-cell subpopulations residing in the SMZ, which could represent physiologic equivalents of distinct SMZL subtypes. Furthermore, they indicate that in SMZL, as in other B cell malignancies, a complementarity imprint of antigen selection might be witnessed either by IGHV, IGKV, or IGLV rearranged sequences.

Molecular biology of lymphoma in the microarray era

This review will focus on the molecular biology of lymphoproliferative disorders with emphasis on lymphomas. The spectrum of known recurrent gene rearrangements found in lymphomas will be outlined and their relevance to diagnosis and subclassification of disease will be discussed. Finally, a survey of the current trends in gene expression profiling of lymphomas by microarray technology will be presented with reference to implications for diagnosis, classification, prognosis and treatment.

Reconsidering the human immunoglobulin heavy-chain locus:

We have used a bioinformatics approach to evaluate the completeness and functionality of the reported human immunoglobulin heavy-chain IGHD gene repertoire. Using the hidden Markov-model-based iHMMune-align program, 1,080 relatively unmutated heavy-chain sequences were aligned against the reported repertoire. These alignments were compared with alignments to 1,639 more highly mutated sequences. Comparisons of the frequencies of gene utilization in the two databases, and analysis of features of aligned IGHD gene segments, including their length, the frequency with which they appear to mutate, and the frequency with which specific mutations were seen, were used to determine the reliability of alignments to the less commonly seen IGHD genes. Analysis demonstrates that IGHD4-23 and IGHD5-24, which have been reported to be open reading frames of uncertain functionality, are represented in the expressed gene repertoire; however, the functionality of IGHD6-25 must be questioned. Sequence similarities make the unequivocal identification of members of the IGHD1 gene family problematic, although all genes except IGHD1-14*01 appear to be functional. On the other hand, reported allelic variants of IGHD2-2 and of the IGHD3 gene family appear to be nonfunctional, very rare, or nonexistent. Analysis also suggests that the reported repertoire is relatively complete, although one new putative polymorphism (IGHD3-10*p03) was identified. This study therefore confirms a surprising lack of diversity in the available IGHD gene repertoire, and restriction of the germline sequence databases to the functional set described here will substantially improve the accuracy of IGHD gene alignments and therefore the accuracy of analysis of the V-D-J junction.

Identification of two new alleles, IGHV3-23*04 and IGHJ6*04, and the complete sequence of the IGHV3-h pseudogene in the human immunoglobulin locus and their prevalences in Danish Caucasians

More than 100 variable (V), 27 diversity (D), and six joining (J) genes are encoded in the human heavy chain locus, and many of these genes exists in different allelic forms. The number of genes and the allelic differences help to create diversity in the immunoglobulin receptors, a key feature of the adaptive immune system. We here report the identification of two novel and seemingly functional alleles of human heavy chain genes. The variable IGHV3-23*04 allele is found with an allele frequency of 0.21 amongst Danish Caucasians, whereas the novel joining IGHJ6*04 allele is rare (allele frequency 0.02). We also report the full sequence of IGHV3-h. The gene exists in two allelic forms but is only found in 58% of the Danish Caucasians studied. The methionine translation initiation codon is mutated, ATG-->AAG, and we therefore propose that the gene is a pseudogene incapable of being translated.

Immunoglobulin light chain repertoire in chronic lymphocytic leukemia

Immunoglobulin kappa (IGK) and immunoglobulin lambda (IGL) light chain repertoire was analyzed in 276 chronic lymphocytic leukemia (CLL) cases and compared with the relevant repertoires from normal, autoreactive, and neoplastic cells. Twenty-one functional IGKV genes were used in IGKV-J rearrangements of 179 kappa-CLL cases; the most frequent genes were IGKV3-20(A27), IGKV1-39/1D-39(O2/O12), IGKV1-5(L12), IGKV4-1(B3), and IGKV2-30(A17); 90 (50.3%) of 179 IGK sequences were mutated (similarity < 98%). Twenty functional IGLV genes were used in IGLV-J rearrangements of 97 lambda-CLL cases; the most frequent genes were IGLV3-21(VL2-14), IGLV2-8(VL1-2), and IGLV2-14(VL1-4); 44 of 97 IGL sequences (45.4%) were mutated. Subsets with "CLL-biased" homologous complementarity-determining region 3 (CDR3) were identified: (1) IGKV2-30-IGKJ2, 7 sequences with homologous kappa CDR3 (KCDR3), 5 of 7 associated with homologous IGHV4-34 heavy chains; (2) IGKV1-39/1D-39-IGKJ1/4, 4 unmutated sequences with homologous KCDR3, 2 of 4 associated with homologous IGHV4-39 heavy chains; (3) IGKV1-5-IGKJ1/3, 4 sequences with homologous KCDR3, 2 of 4 associated with unmutated nonhomologous IGHV4-39 heavy chains; (4) IGLV1-44-IGLJ2/3, 2 sequences with homologous lambda CDR3 (LCDR3), associated with homologous IGHV4-b heavy chains; and (5) IGLV3-21-IGLJ2/3, 9 sequences with homologous LCDR3, 3 of 9 associated with homologous IGHV3-21 heavy chains. The existence of subsets that comprise given IGKV-J/IGLV-J domains associated with IGHV-D-J domains that display homologous CDR3 provides further evidence for the role of antigen in CLL pathogenesis.

Production of recombinant Ig molecules from antigen-selected single B cells and restricted usage of Ig-gene segments by anti-D antibodies

The Ig-genes of the heavy chains in anti-D-specific hybridomas and Fab/scFv-fragments selected from phage-display libraries are restricted to a group of closely related genes (IGHV3s genes). We analyzed the Ig-gene repertoire in anti-D-specific B cells of two hyperimmunized donors using a completely different method. Single B cells were cultured for 10 days in an EL4.B5 culture system. mRNA from anti-D-producing B cells was reverse transcribed into cDNA. Heavy- and light-chain gene rearrangements were amplified by PCR reactions, sequenced and cloned into a pNUT-vector system, thereby allowing the production of complete IgG and IgM. Eleven anti-D-specific B-cell clones were isolated and analyzed. Eight of these clones (including IgM-producing clones) had IGHV3s genes. We demonstrated that functional anti-D-specific IgM (4 clones) and IgG (2 clones) was produced. Using a new method, we analyzed the IGHV gene repertoire of anti-D-specific B cells of hyperimmunized donors and showed that it is indeed restricted. Moreover, we found a high frequency (1:100 and 1:500) of anti-D-specific B cells in the peripheral B cells of hyperimmunized donors. We suggest that this approach could be applied for the selection of human mAbs from immunized donors and for the analysis of Ig-gene repertoires at the single-B cell level.

Dramatic variation of the vomeronasal pheromone receptor gene repertoire among five orders of placental and marsupial mammals

Pheromones are chemicals emitted and sensed by conspecifics to elicit social and sexual responses and are perceived in terrestrial vertebrates primarily by the vomeronasal organ (VNO). Pheromone receptors in the mammalian VNO are encoded by the V1R and V2R gene superfamilies. The V1R superfamily contains 187 and 102 putatively functional genes in the mouse and rat, respectively. To investigate whether this large repertoire size is typical among mammals with functional VNOs, we here describe the V1R repertoires of dog, cow, and opossum based on their draft genome sequences. The dog and cow have only 8 and 32 intact V1R genes, respectively. Thus, the intact V1R repertoire size varies by at least 23-fold among placental mammals with functional VNOs. To our knowledge, this size ratio represents the greatest among-species variation in gene family size of all mammalian gene families. Phylogenetic analysis of placental V1R genes suggests multiple losses of ancestral genes in carnivores and artiodactyls and gains of many new genes by gene duplication in rodents, manifesting massive gene births and deaths. We also identify 49 intact opossum V1R genes and discover independent expansions of the repertoire in placentals and marsupials. We further show a concordance between the V1R repertoire size and the complexity of VNO morphology, suggesting that the latter could indicate the sophistication of pheromone communications within species. In sum, our results demonstrate tremendous diversity and rapid evolution of mammalian V1R gene inventories and caution the generalization of VNO biology from rodents to all mammals.

Assembly and analysis of the mouse immunoglobulin kappa gene sequence

The mechanisms regulating V gene usage leading to the immunoglobulin (Ig) repertoire have been of interest for many years but are only partially defined. To gain insight into these processes, we have assembled the nucleotide sequence of the Mus musculus Igkappa locus using data recently made available from genome-wide sequencing efforts. We found the locus to be 3.21 Mb in length and mapped all known functional, pseudo- and relic V gene segments onto the sequence, along with known regulatory elements. We corrected errors in former gene assignments, positions and orientations and identified a novel Vkappa4 gene segment. This assembly allowed the establishment of a unified nomenclature for the V genes based on their relative positions similar to the nomenclature system adopted for the human Ig loci. The 5' boundary of the locus is defined by the presence of the tumor-associated calcium-signal transducer-2 gene located 19 kb upstream of Vkappa24-140, the most distal V gene. No non- Vkappa genes were found in the sequence of the locus. Detailed analysis of the sequences 0.5 kb upstream, within, and 0.5 kb downstream of each potentially functional V gene revealed interesting patterns of statistically significant clustering of transcription factor consensus binding sites, generally specific to a particular family. We found E boxes were clustered not only in promoter regions, but also nearby recombination signal sequences. Family members of Vkappa4/5 genes exhibit a conserved pattern of octamer sites in their downstream regions, as well as Ebf sites in their introns, and Lef-1 sites in their upstream regions. We discuss potential functional implications of these findings in the context of possible combinatorial mechanisms for targeting V genes for rearrangement. The assembled sequence and its analyses are available as a resource to the scientific community.

Exonuclease activity and P nucleotide addition in the generation of the expressed immunoglobulin repertoire

Background

Immunoglobulin rearrangement involves random and imprecise processes that act to both create and constrain diversity. Two such processes are the loss of nucleotides through the action of unknown exonuclease(s) and the addition of P nucleotides. The study of such processes has been compromised by difficulties in reliably aligning immunoglobulin genes and in the partitioning of nucleotides between segment ends, and between N and P nucleotides.

Results

A dataset of 294 human IgM sequences was created and partitioned with the aid of a probabilistic model. Non-random removal of nucleotides is seen between the three IGH gene types with the IGHV gene averaging removals of 1.2 nucleotides compared to 4.7 for the other gene ends (p < 0.001). Individual IGHV, IGHD and IGHJ gene subgroups also display statistical differences in the level of nucleotide loss. For example, within the IGHJ group, IGHJ3 has average removals of 1.3 nucleotides compared to 6.4 nucleotides for IGHJ6 genes (p < 0.002). Analysis of putative P nucleotides within the IgM and pooled datasets revealed only a single putative P nucleotide motif (GTT at the 3' D-REGION end) to occur at a frequency significantly higher then would be expected from random N nucleotide addition.

Conclusions

The loss of nucleotides due to the action of exonucleases is not random, but is influenced by the nucleotide composition of the genes. P nucleotides do not make a significant contribution to diversity of immunoglobulin sequences. Although palindromic sequences are present in 10% of immunologlobulin rearrangements, most of the 'palindromic' nucleotides are likely to have been inserted into the junction during the process of N nucleotide addition. P nucleotides can only be stated with confidence to contribute to diversity of less than 1% of sequences. Any attempt to identify P nucleotides in immunoglobulins is therefore likely to introduce errors into the partitioning of such sequences.

IMGT-ONTOLOGY and IMGT databases, tools and Web resources for immunogenetics and immunoinformatics

The international ImMunoGeneTics information system (IMGT; http://imgt.cines.fr), is a high quality integrated information system specialized in immunoglobulins (IG), T cell receptors (TR), major histocompatibility complex (MHC), and related proteins of the immune system (RPI) of human and other vertebrates, created in 1989, by the Laboratoire d'ImmunoGénétique Moléculaire (LIGM; Université Montpellier II and CNRS) at Montpellier, France. IMGT provides a common access to standardized data which include nucleotide and protein sequences, oligonucleotide primers, gene maps, genetic polymorphisms, specificities, 2D and 3D structures. IMGT consists of several sequence databases (IMGT/LIGM-DB, IMGT/MHC-DB, IMGT/PRIMER-DB), one genome database (IMGT/GENE-DB) and one 3D structure database (IMGT/3Dstructure-DB), interactive tools for sequence analysis (IMGT/V-QUEST, IMGT/JunctionAnalysis, IMGT/PhyloGene, IMGT/Allele-Align), for genome analysis (IMGT/GeneSearch, IMGT/GeneView, IMGT/LocusView) and for 3D structure analysis (IMGT/StructuralQuery), and Web resources ("IMGT Marie-Paule page") comprising 8000 HTML pages. IMGT other accesses include SRS, FTP, search by BLAST, etc. By its high quality and its easy data distribution, IMGT has important implications in medical research (repertoire in autoimmune diseases, AIDS, leukemias, lymphomas, myelomas), veterinary research, genome diversity and genome evolution studies of the adaptive immune responses, biotechnology related to antibody engineering (single chain Fragment variable (scFv), phage displays, combinatorial libraries) and therapeutical approaches (grafts, immunotherapy). IMGT is freely available at http://imgt.cines.fr.