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.     

<Emphasis Type="Italic">IGHV1</Emphasis>,<Emphasis Type="Italic"> IGHV5</Emphasis> and<Emphasis Type="Italic"> IGHV7</Emphasis> subgroup genes in the Rhesus macaque

The diversity of the antibody response is achieved, in part, by rearrangement of different immunoglobulin (Ig) genes. The Ig heavy chain is made up of a variable region (IGHV), a diversity region (IGHD) and a joining region (IGHJ). Human germline IGHV genes have been grouped into seven multigene subgroups. Size and usage of these subgroups is not equal, the IGHV3 subgroup is the most commonly used (36%), followed by IGHV1/7 (26%), then IGHV4, IGHV5, IGHV2, IGHV6 (15%, 12%, 4%, 3% respectively). The rhesus macaque (Macaca mulatta) is a useful non-human primate model for studies of infection and the database of germline Ig genes for the macaque is gradually growing to become a useful tool in the study of B-cell responses. The proportions of IGHV subgroup usage in the macaque are similar to those in man. Representatives from IGHV3 and IGHV4 subgroups for the macaque have been published, as have germline sequences of the IGHD and IGHJ genes. However, to date there have been no sequences published from the second largest IGHV subgroup, IGHV1. We report the isolation and sequencing of a genomic fragment containing an IGHV1 gene from the macaque. Polymerase chain reaction (PCR) primers designed from this sequence enabled us to amplify and sequence 25 new IGHV1 germline genes. We also isolated two IGHV7 genes, using the same primers, and two IGHV5 genes, using human IGHV5 primers.     

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.     

The mouse antibody heavy chain repertoire is germline-focused and highly variable between inbred strains

The human and mouse antibody repertoires are formed by identical processes, but like all small animals, mice only have sufficient lymphocytes to express a small part of the potential antibody repertoire. In this study, we determined how the heavy chain repertoires of two mouse strains are generated. Analysis of IgM- and IgG-associated VDJ rearrangements generated by high-throughput sequencing confirmed the presence of 99 functional immunoglobulin heavy chain variable (IGHV) genes in the C57BL/6 genome, and inferred the presence of 164 IGHV genes in the BALB/c genome. Remarkably, only five IGHV sequences were common to both strains. Compared with humans, little N nucleotide addition was seen in the junctions of mouse VDJ genes. Germline human IgG-associated IGHV genes are rare, but many murine IgG-associated IGHV genes were unmutated. Together these results suggest that the expressed mouse repertoire is more germline-focused than the human repertoire. The apparently divergent germline repertoires of the mouse strains are discussed with reference to reports that inbred mouse strains carry blocks of genes derived from each of the three subspecies of the house mouse. We hypothesize that the germline genes of BALB/c and C57BL/6 mice may originally have evolved to generate distinct germline-focused antibody repertoires in the different mouse subspecies.     

Characterization of the expressed immunoglobulin IGHV repertoire in the New World marmoset Callithrix jacchus

The common marmoset (Callithrix jacchus jacchus) is a member of the Callithrichinae, a family of outbred New World primates with limited MHC polymorphisms and a propensity to develop spontaneous or experimentally induced autoimmunity. C. jacchus marmosets are susceptible to experimental allergic encephalomyelitis (EAE), and spontaneously develop autoimmune colitis and thyroiditis. Such disease models approximate the complexity of human autoimmune disorders, and allow an investigation of the respective roles of T-cell and antibody responses to self-antigens in outbred species. A key issue for further definition of the pathogenic antibody responses in human autoimmunity is to understand the diversity of the immunoglobulin repertoire in primate models. Here, we characterized the expressed immunoglobulin IGHV repertoire of the C. jacchus marmoset. Six IGHV subgroups were identified which show a high degree of sequence similarity to their human IGHV counterparts (IGHV1, IGHV3, IGHV4, IGHV5, IGHV6, and IGHV7). As in the expressed human IGHV repertoire, the framework regions are more conserved when compared to the complementarity-determining regions (CDRs), with the greatest degree of variability located in CDR3. Predicted structural features are highly conserved between C. jacchus and human IGHV. This information now provides a framework for studies of the antigen-specific repertoire of pathogenic antibodies in EAE and other immune-mediated diseases.     

Single-nucleotide polymorphism analysis by pyrosequencing

There is a growing demand for high-throughput methods for analysis of single-nucleotide polymorphic (SNP) positions. Here, we have evaluated a novel sequencing approach, pyrosequencing, for such purposes. Pyrosequencing is a sequencing-by-synthesis method in which a cascade of enzymatic reactions yields detectable light, which is proportional to incorporated nucleotides. One feature of typing SNPs with pyrosequencing is that each allelic variant will give a unique sequence compared to the two other variants. These variants can easily be distinguished by a pattern recognition software. The software displays the allelic alternatives and allows for direct comparison with the pyrosequencing raw data. For optimal determination of SNPs, various protocols of nucleotide dispensing order were investigated. Here, we demonstrate that typing of SNPs can efficiently be performed by pyrosequencing using an automated system for parallel analysis of 96 samples in approximately 5 min, suitable for large-scale screening and typing of SNPs.     

Combined Patterns of IGHV Repertoire and Cytogenetic/Molecular Alterations in Monoclonal B Lymphocytosis versus Chronic Lymphocytic Leukemia

Background

Chronic lymphocytic leukemia (CLL)-like monoclonal B lymphocytosis (MBL) with (MBL^hi) or without (MBL^lo) absolute B-lymphocytosis precedes most CLL cases,the specific determinants for malignant progression remaining unknown.

Methodology/Principal Findings

For this purpose, simultaneous iFISH and molecular analysis of well-established cytogenetic alterations of chromosomes 11, 12, 13, 14 and 17 together with the pattern of rearrangement of the IGHV genes were performed in CLL-like cells from MBL and CLL cases. Our results based on 78 CLL-like MBL and 117 CLL clones from 166 subjects living in the same geographical area, show the existence of three major groups of clones with distinct but partially overlapping patterns of IGHV gene usage, IGHV mutational status and cytogenetic alterations. These included a group enriched in MBL^lo clones expressing specific IGHV subgroups (e.g. VH3-23) with no or isolated good-prognosis cytogenetic alterations, a second group which mainly consisted of clinical MBL^hi and advanced stage CLL with a skewed but different CLL-associated IGHV gene repertoire (e.g. VH1-69), frequently associated with complex karyotypes and poor-prognosis cytogenetic alterations, and a third group of clones with intermediate features, with prevalence of mutated IGHV genes, and higher numbers of del(13q)⁺ clonal B-cells.

Conclusions/Significance

These findings suggest that the specific IGHV repertoire and IGHV mutational status of CLL-like B-cell clones may modulate the type of cytogenetic alterations acquired, their rate of acquisition and/or potentially also their clinical consequences. Further long-term follow-up studies investigating the IGHV gene repertoire of MBL^lo clones in distinct geographic areas and microenvironments are required to confirm our findings and shed light on the potential role of some antigen-binding BCR specificities contributing to clonal evolution.     

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.     

Biased use of the IGHV4 family and evidence for antigen selection in Chlamydophila psittaci-negative ocular adnexal extranodal marginal zone lymphomas

Extranodal marginal zone lymphomas (EMZL) are the most common lymphomas in the ocular adnexa. The etiology and potential role for antigenic stimulation in these lymphomas are still controversial. We have examined IGHV gene usage and mutations in 67 Chlamydophila psittaci-negative ocular adnexal EMZL. Clonal IGHV gene sequences were identified in 43 tumors originating from the orbit (19), conjunctivae (18) and lacrimal gland (6). Forty four potentially functional clonal IGHV gene sequences were detected with overrepresentation of the IGHV4 family and IGHV4-34 gene. All but 3 sequences were mutated with the average percent homology to the germ line of 93.5±6.1. Multinomial model and Focused binomial test demonstrated evidence for positive and/or negative antigen selection in 59% of the potentially functional IGHV genes. Intraclonal variation was detected in 8 of 11 tumor specimens. Overall our findings demonstrate that C. psittaci-negative ocular adnexal EMZL exhibit biased usage of IGHV families and genes with evidence for intraclonal heterogeneity and antigen selection in multiple tumors, implicating B-cell receptor-mediated antigen stimulation in the pathogenesis of these lymphomas.     

The human antibody response to porcine xenoantigens is encoded by IGHV3-11 and IGHV3-74 IgVH germline progenitors.

Preformed and induced Ab responses present a major immunological barrier to the use of pig organs for human xenotransplantation. We generated IgM and IgG gene libraries established from lymphocytes of patients treated with a bioartificial liver (BAL) containing pig hepatocytes and used these libraries to identify IgVH genes that encode human Ab responses to pig xenoantigens. Genes encoded by the VH3 family are increased in expression in patients following BAL treatment. cDNA libraries representing the VH3 gene family were generated, and the relative frequency of expression of genes used to encode the Ab response was determined at days 0, 10, and 21. Ig genes derived from the IGHV3-11 and IGHV3-74 germline progenitors increase in frequency post-BAL. The IGHV3-11 gene encodes 12% of VH3 cDNA clones expressed as IgM Abs at day 0 and 32.4-39.0% of cDNA clones encoding IgM Abs in two patients at day 10. IGHV3-11 and IGHV3-74 genes encoding IgM Abs in these patients are expressed without evidence of somatic mutation. By day 21, an isotype switch occurs and IGHV3-11 IgVH progenitors encode IgG Abs that demonstrate somatic mutation. We cloned these genes into a phagemid vector, expressed these clones as single-chain Abs, and demonstrated that the IGHV3-11 gene encodes Abs with the ability to bind to the gal alpha (1,3) gal epitope. Our results demonstrate that the xenoantibody response in humans is encoded by IgVH genes restricted to IGHV3-11 and IGHV3-74 germline progenitors. IgM Abs are expressed in germline configuration and IgG Abs demonstrate somatic mutations by day 21.     

Uneven Distribution of Repetitive Trinucleotide Motifs in Human Immunoglobulin Heavy Variable Genes

Insertions and deletions of entire codons have recently been discovered as a mechanism by which B cells, in addition to conventional base substitution, evolve the antibodies produced by their immunoglobulin genes. These events frequently seem to involve repetitive sequence motifs in the antibody-encoding genes, and it has been suggested that they occur through polymerase slippage. In order to better understand the process of codon deletion, we have analyzed the human immunoglobulin heavy variable (IGHV) germline gene repertoire for the presence of trinucleotide repeats. Such repeats would ensure that the reading frame is maintained in the case of a deletional event, as slippage over multiples of three bases would be favored. We demonstrate here that IGHV genes specifically carry repetitive trinucleotide motifs in the complementarity-determining regions (CDR) 1 and 2, thus making these parts of the genes that encode highly flexible structures particularly prone to functional deletions. We propose that the human IGHV repertoire carries inherent motifs that allow an antibody response to develop efficiently by targeting codon deletion events to the parts of the molecule that are likely to be able to harbor such modifications. Received: 10 April 2001 / Accepted: 27 August 2001     

Intraclonal cell expansion and selection driven by B cell receptor in chronic lymphocytic leukemia

The mutational status of the immunoglobulin heavy-chain variable region (IGHV) genes utilized by chronic lymphocytic leukemia (CLL) clones defines two disease subgroups. Patients with unmutated IGHV have a more aggressive disease and a worse outcome than patients with cells having somatic IGHV gene mutations. Moreover, up to 30% of the unmutated CLL clones exhibit very similar or identical B cell receptors (BcR), often encoded by the same IG genes. These "stereotyped" BcRs have been classified into defined subsets. The presence of an IGHV gene somatic mutation and the utilization of a skewed gene repertoire compared with normal B cells together with the expression of stereotyped receptors by unmutated CLL clones may indicate stimulation/selection by antigenic epitopes. This antigenic stimulation may occur prior to or during neoplastic transformation, but it is unknown whether this stimulation/selection continues after leukemogenesis has ceased. In this study, we focused on seven CLL cases with stereotyped BcR Subset #8 found among a cohort of 700 patients; in six, the cells expressed IgG and utilized IGHV4-39 and IGKV1-39/IGKV1D-39 genes, as reported for Subset #8 BcR. One case exhibited special features, including expression of IgM or IgG by different subclones consequent to an isotype switch, allelic inclusion at the IGH locus in the IgM-expressing cells and a particular pattern of cytogenetic lesions. Collectively, the data indicate a process of antigenic stimulation/selection of the fully transformed CLL cells leading to the expansion of the Subset #8 IgG-bearing subclone.     

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.Electronic Supplementary Material Supplementary material is available for this article at     

抗人AFP抗体噬菌体呈现文库的构建筛选及基因序列测定

从人ＡＦＰ免疫小鼠脾细胞ｍＲＮＡ中扩增出全套抗体Ｖ区基因并随机拼接为ＳｃＦｖ基因,构建全套ＳｃＦｖ基因噬菌体呈现文库,经两轮ｐａｎｎｉｎｇ筛选富集,一次从９４个单个重组噬菌体克隆中筛选到１４个具有ＡＦＰ结合活性的克隆,测定两个阳性克隆中ＳｃＦｖ基因的核苷酸序列,获得一个Ｖ＿Ｈ基因和两个Ｖ＿Ｋ基因,其基因序列分别与鼠ＩｇＶ＿ＨＪ５５８、Ｖ＿ｋ－ＯＸ１和Ｖ＿ｋ４／５家族同源性最高;推导出的氨基酸序列中均含有抗体Ｖ区特征性的两个恒定的半胱氨酸残基、具有明确的三个ＣＤＲ和四个ＦＲ序列,表明这三个基因均系新发现的功能性鼠抗体Ｖ区基因序列。     

Organization of the variable region of the immunoglobulin heavy-chain gene locus of the rat

We have mapped and annotated the variable region of the immunoglobulin heavy (IGH) gene locus of the Brown Norway (BN) rat (assembly V3.4; Rat Genomic Sequence Consortium). In addition to known variable region genes, we found 12 novel previously unidentified functional IGHV genes and 1 novel functional IGHD gene. In total, the variable region of the rat IGH locus is composed of at least 353 unique IGHV genes, 21 IGHD genes, and 5 IGHJ genes, of which 131, 14, and 4 are potentially functional genes, respectively. Of all species studied so far, the rat seems to have the highest number of functional IGHV genes in the genome. Rat IGHV genes can be classified into 13 IGHV families based on nucleotide sequence identity. The variable region of the BN rat spans a total length of approximately 4.9 Mb and is organized in a typical translocon organization. Like the mouse, members of the various IGHV gene families are more or less grouped together on the genome, albeit some members of IGHV gene families are found intermingled with each other. In the rat, the largest IGHV gene families are IGHV1, IGHV2, and IGHV5. The overall conclusion is that the genomic organization of the variable region of the rat IGH locus is strikingly similar to that of the mouse, illustrating the close evolutionary relationship between these two species.     

Method enabling fast partial sequencing of cDNA clones

Pyrosequencing is a nonelectrophoretic single-tube DNA sequencing method that takes advantage of cooperativity between four enzymes to monitor DNA synthesis. To investigate the feasibility of the recently developed technique for tag sequencing, 64 colonies of a selected cDNA library from human were sequenced by both pyrosequencing and Sanger DNA sequencing. To determine the needed length for finding a unique DNA sequence, 100 sequence tags from human were retrieved from the database and different lengths from each sequence were randomly analyzed. An homology search based on 20 and 30 nucleotides produced 97 and 98% unique hits, respectively. An homology search based on 100 nucleotides could identify all searched genes. Pyrosequencing was employed to produce sequence data for 30 nucleotides. A similar search using BLAST revealed 16 different genes. Forty-six percent of the sequences shared homology with one gene at different positions. Two of the 64 clones had unique sequences. The search results from pyrosequencing were in 100% agreement with conventional DNA sequencing methods. The possibility of using a fully automated pyrosequencer machine for future high-throughput tag sequencing is discussed.     

Systematic inference of copy-number genotypes from personal genome sequencing data reveals extensive olfactory receptor gene content diversity

Copy-number variations (CNVs) are widespread in the human genome, but comprehensive assignments of integer locus copy-numbers (i.e., copy-number genotypes) that, for example, enable discrimination of homozygous from heterozygous CNVs, have remained challenging. Here we present CopySeq, a novel computational approach with an underlying statistical framework that analyzes the depth-of-coverage of high-throughput DNA sequencing reads, and can incorporate paired-end and breakpoint junction analysis based CNV-analysis approaches, to infer locus copy-number genotypes. We benchmarked CopySeq by genotyping 500 chromosome 1 CNV regions in 150 personal genomes sequenced at low-coverage. The assessed copy-number genotypes were highly concordant with our performed qPCR experiments (Pearson correlation coefficient 0.94), and with the published results of two microarray platforms (95-99% concordance). We further demonstrated the utility of CopySeq for analyzing gene regions enriched for segmental duplications by comprehensively inferring copy-number genotypes in the CNV-enriched >800 olfactory receptor (OR) human gene and pseudogene loci. CopySeq revealed that OR loci display an extensive range of locus copy-numbers across individuals, with zero to two copies in some OR loci, and two to nine copies in others. Among genetic variants affecting OR loci we identified deleterious variants including CNVs and SNPs affecting ~15% and ~20% of the human OR gene repertoire, respectively, implying that genetic variants with a possible impact on smell perception are widespread. Finally, we found that for several OR loci the reference genome appears to represent a minor-frequency variant, implying a necessary revision of the OR repertoire for future functional studies. CopySeq can ascertain genomic structural variation in specific gene families as well as at a genome-wide scale, where it may enable the quantitative evaluation of CNVs in genome-wide association studies involving high-throughput sequencing.