Andrew M. Collins

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 individuals 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.

iHMMune-align

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/.

A Temporal Model of Human IgE and IgG Antibody Function

The diversity of the human antibody repertoire that is generated by V(D)J gene rearrangement is extended by nine constant region genes that give antibodies their complex array of effector functions. The application of high throughput sequencing to the study of V(D)J gene rearrangements has led to significant recent advances in our understanding of the antigen-binding repertoire. In contrast, our understanding of antibody function has changed little, and mystery still surrounds the existence of four distinctive IgG subclasses. Recent observations from murine models and from human studies of VDJ somatic point mutations suggest that the timing of emergence of cells from the germinal center may vary as a consequence of class switching. This should lead to predictable differences in affinity between isotypes. These differences, and varying abilities of the isotypes to fix complement and bind FcRs, could help coordinate the humoral defenses over the time course of a response. We therefore propose a Temporal Model of human IgE and IgG function in which early emergence of IgE sensitizes sentinel mast cells while switching to IgG3 recruits FcγR-mediated functions to the early response. IgG1 then emerges as the major effector of antigen clearance, and subsequently IgG2 competes with IgG1 to produce immune complexes that slow the inflammatory drive. Persisting antigen may finally stimulate high affinity IgG4 that outcompetes other isotypes and can terminate IgG1/FcγR-mediated activation via the inhibitory FcγRIIB. In this way, IgG antibodies of different subclasses, at different concentrations and with sometimes opposing functions deliver cohesive, protective immune function.

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.

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.

Human immunoglobulin classes and subclasses show variability in VDJ gene mutation levels

Somatic point mutations provide glimpses into B‐cell histories, and mutation numbers generally correlate with antibody affinity. We recently proposed a model of human isotype function, based in part on mutation analysis, in which the dominant pathway of isotype switching involves B cells moving sequentially through the four immunoglobulin (Ig) G subclasses. This should result in predictable differences in affinity between isotypes, and this helps explain how different isotypes work together. The model built on analysis of rearranged immunoglobulin heavy chain sequences amplified from Papua New Guinean villagers, which showed highly significant differences in the mean number of V‐REGION mutations in sequences, associated with the different IgG subclasses. To determine whether this relationship between mutation levels and isotypes is a more general phenomenon, the present study was conducted in healthy, urban residents of Sydney, Australia. VDJ sequences were generated from eight individuals, using 454 pyrosequencing, from cells expressing all isotypes except IgD and IgE. This resulted in 35 118 unique, productive VDJ sequences for the study. The data confirm that VDJ genes associated with progressively more 3′ Ig heavy chain gamma (IGHG) constant region genes show increasing levels of point mutation. Mean V‐REGION mutations in IgA1 and IgA2 sequences were similar. Patterns of mutations also differed between isotypes. Despite their association with T‐independent responses, IgG2 sequences showed significantly more mutational evidence of antigen selection than other IgG isotypes. Antigen selection was also significantly higher in IgA2 than in IgA1 sequences, raising the possibility of a preferential switch pathway from IGHG2 to IGHA2.

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.

Clustering-based identification of clonally-related immunoglobulin gene sequence sets

BackgroundClonal expansion of B lymphocytes coupled with somatic mutation and antigen selection allow the mammalian humoral immune system to generate highly specific immunoglobulins (IG) or antibodies against invading bacteria, viruses and toxins. The availability of high-throughput DNA sequencing methods is providing new avenues for studying this clonal expansion and identifying the factors guiding the generation of antibodies. The identification of groups of rearranged immunoglobulin gene sequences descended from the same rearrangement (clonally-related sets) in very large sets of sequences is facilitated by the availability of immunoglobulin gene sequence alignment and partitioning software that can accurately predict component germline gene, but has required painstaking visual inspection and analysis of sequences.ResultsWe have developed and implemented an algorithm for identifying sets of clonally-related sequences in large human immunoglobulin heavy chain gene variable region sequence sets. The program processes sequences that have been partitioned using iHMMune-align, and uses pairwise comparisons of CDR3 sequences and similarity in IGHV and IGHJ germline gene assignments to construct a distance matrix. Agglomerative hierarchical clustering is then used to identify likely groups of clonally-related sequences. The program is available for download from http://www.cse.unsw.edu.au/~ihmmune/ClonalRelate/ClonalRelate.zip.ConclusionsThe method was evaluated on several benchmark datasets and provided a more accurate and considerably faster identification of clonally-related immunoglobulin gene sequences than visual inspection by domain experts.

Analysis of IgE Antibodies from a Patient with Atopic Dermatitis: Biased V Gene Usage and Evidence for Polyreactive IgE Heavy Chain Complementarity-Determining Region 3

To better understand V gene usage, specificity, and clonal origins of IgE Abs in allergic reactions, we have constructed a combinatorial Ab library from the mRNA of an adult patient with atopic dermatitis. Sequence analysis of random clones revealed that 33% of clones used the IGHV6-1 H chain V gene segment, the only member of the VH6 gene family. IGHV6-1 is rarely used in the expressed adult repertoire; however, it is associated with fetal derived Abs. Features of the VH6 rearrangements included short complementarity-determining region 3, frequent use of IGHD7-27 D gene, and little nucleotide addition at the D-J junction. There was also a low level of mutation compared with VH1, VH3, and VH4 rearrangements. The library was expressed as phage-Fab fusions, and specific phage selected by panning on the egg allergen ovomucoid. Upon expression as soluble IgE Fabs, 12 clones demonstrated binding to ovomucoid, skim milk, and BSA by ELISA. Nucleotide sequencing demonstrated that the IGHV6-1 V gene segment encoded each of the 12 multiply reactive IgE Fabs. A cyclic peptide was designed from the complementarity-determining region 3 of several of these clones. The cyclic peptide bound both self and nonself Ags, including ovomucoid, human IgG, tetanus toxoid, and human and bovine von Willebrand factor. These results suggest that some IgE Abs may bind more than one Ag, which would have important implications for understanding the multiple sensitivities seen in conditions such as atopic dermatitis.

IL-16 Regulation of Human Mast Cells/Basophils and Their Susceptibility to HIV-1

AIDS patients often contain HIV-1-infected mast cells (MCs)/basophils in their peripheral blood, and in vivo-differentiated MCs/basophils have been isolated from the blood of asthma patients that are HIV-1 susceptible ex vivo due to their surface expression of CD4 and varied chemokine receptors. Because IL-16 is a ligand for CD4 and/or an undefined CD4-associated protein, the ability of this multifunctional cytokine to regulate the development of human MCs/basophils from nongranulated progenitors residing in cord or peripheral blood was evaluated. After 3 wk of culture in the presence of c-kit ligand, IL-16 induced the progenitors residing in the blood of normal individuals to increase their expression of chymase and tryptase about 20-fold. As assessed immunohistochemically, >80% of these tryptase+ and/or chymase+ cells expressed CD4. The resulting cells responded to IL-16 in an in vitro chemotaxis assay, and this biologic response could be blocked by anti-IL-16 and anti-CD4 Abs as well as by a competitive peptide inhibitor corresponding to a sequence in the C-terminal domain of IL-16. The additional finding that IL-16 induces calcium mobilization in the HMC-1 cell line indicates that IL-16 acts directly on MCs and their committed progenitors. IL-16-treated MCs/basophils also are less susceptible to infection by an M/R5-tropic strain of HIV-1. Thus, IL-16 regulates MCs/basophils at a number of levels, including their vulnerability to retroviral infection.