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. Electronic supplementary materials The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

The immunoglobulin <Emphasis Type="Italic">IGHD</Emphasis> gene locus in C57BL/6 mice

Four immunoglobulin heavy chain diversity (IGHD) gene subgroups (DFL16, DSP2, DQ52, and DST4) have been identified previously in BALB/c mice. Although the locations of most IGHD genes have been established based on restriction map and Southern blot analysis, a complete mouse IGHD gene locus map at the sequence level is still not available. In addition, a previous restriction fragment length polymorphism study suggested that significant difference in the IGHD gene locus exists between C57BL/6 and BALB/c mice. The author has now analyzed the C57BL/6 mouse genomic data and established a complete map of the IGHD gene locus. All four IGHD subgroups previously identified in BALB/c mice were found to be present in C57BL/6 mice. However, unlike the BALB/c mice, which have at least 13 IGHD genes, the C57BL/6 genome contains only ten IGHD genes, which include one DFL16, six DSP2, one DQ52, and two DST4 genes. There are also differences in the coding regions of the DST4 and DQ52 genes between the two mouse strains.     

Structure of the catfish IGH locus: analysis of the region including the single functional IGHM gene

The catfish IGH locus is large (∼1 Mb) and complex, having undergone multiple internal duplications and transpositions. To define the structure of the locus that contains the single expressed IGHM gene, two overlapping bacterial-artificial-chromosome (BAC) clones spanning the most 3′ end of the channel catfish immunoglobulin heavy (IGH) chain locus have been completely sequenced. The analyses created a contig of 257,153 bp containing 55 VH, 6 D, 12 JH genes and the IGH constant region genes encoding the functional secreted and membrane forms of IgM and the membrane form of IgD. This analysis revealed three major features. First, no C-region genes were found aside from the previously described IGHM1 and IGHD1, with the latter gene being the most 3′ C-region gene of the catfish IGH locus. There was no evidence in the region sequenced for genes that could encode an Ig class similar to the IgZ/IgT described in zebrafish, trout and pufferfish. Second, there are a high number of VH pseudogenes, 28 out of 55 (51%). In contrast, the entire zebrafish IGH locus has 40 functional VH genes and eight pseudogenes (17%). Third, an internal duplication of a 52.4-kb block of VH genes has occurred. These observations suggest that the IGH locus of teleost fish varies significantly from species to species in the diversity of C-region genes as well as the numbers of genes encoding V regions.     

Linkage mapping of the locus responsible for male pseudohermaphroditism (mp) on rat chromosome 7.

TF is a mutant rat strain showing male pseudohermaphroditism controlled by an autosomal single recessive gene (mp). The affected rats show lack of Leydig cells and androgen deficiency. In this study, we performed linkage analysis using F(2) progeny of crosses between TF and BN strains to determine the chromosomal localization of the mp locus. The mp locus was mapped in a 4 cM region of the distal region of rat chromosome 7 between D7Rat3 and D7Rat115 or D7Rat94. Comparison of the linkage map with corresponding regions of the published rat genome sequence revealed several candidate genes for the mp mutation, including the Dhh, Tegt, Gdp3, and Amhr2 genes.     

Accuracy and coverage assessment of Oryctolagus cuniculus (rabbit) genes encoding immunoglobulins in the whole genome sequence assembly (OryCun2.0) and localization of the IGH locus to chromosome 20

We report on the analyses of genes encoding immunoglobulin heavy and light chains in the rabbit 6.51× whole genome assembly. This OryCun2.0 assembly confirms previous mapping of the duplicated IGK1 and IGK2 loci to chromosome 2 and the IGL lambda light chain locus to chromosome 21. The most frequently rearranged and expressed IGHV1 that is closest to IG DH and IGHJ genes encodes rabbit VHa allotypes. The partially inbred Thorbecke strain rabbit used for whole-genome sequencing was homozygous at the IGK but heterozygous with the IGHV1a1 allele in one of 79 IGHV-containing unplaced scaffolds and IGHV1a2, IGHM, IGHG, and IGHE sequences in another. Some IGKV, IGLV, and IGHA genes are also in other unplaced scaffolds. By fluorescence in situ hybridization, we assigned the previously unmapped IGH locus to the q-telomeric region of rabbit chromosome 20. An approximately 3-Mb segment of human chromosome 14 including IGH genes predicted to map to this telomeric region based on synteny analysis could not be located on assembled chromosome 20. Unplaced scaffold chrUn0053 contains some of the genes that comparative mapping predicts to be missing. We identified discrepancies between previous targeted studies and the OryCun2.0 assembly and some new BAC clones with IGH sequences that can guide other studies to further sequence and improve the OryCun2.0 assembly. Complete knowledge of gene sequences encoding variable regions of rabbit heavy, kappa, and lambda chains will lead to better understanding of how and why rabbits produce antibodies of high specificity and affinity through gene conversion and somatic hypermutation.     

Detailed chromosomal and radiation hybrid mapping in the proximal part of rat Chromosome 10 and gene order comparison with mouse and human

The rat provides valuable and sometimes unique models of human complex diseases. To fully exploit the rat models in biomedical research, it is important to have access to detailed knowledge of the rat genome organization as well as its relation to the human genome. Rat Chromosome 10 (RNO10) harbors several important cancer-related genes. Deletions in the proximal part of RNO10 were repeatedly found in a rat model for endometrial cancer. To identify functional and positional candidate genes in the affected region, we used radiation hybrid (RH) mapping and single- and dual-color fluorescence in situ hybridization (FISH) techniques to construct a detailed chromosomal map of the proximal part of RNO10. The regional localization of 14 genes, most of them cancer-related (Grin2a, Gspt1, Crebbp, Gfer, Tsc2, Tpsb1, Il9r, Il4, Irf1, Csf2, Sparc, Tp53, Thra1, Gh1), and of five microsatellite markers (D10Mit10, D10Rat42, D10Rat50, D10Rat72, and D10Rat165) was determined on RNO10. For a fifteenth gene, Ppm1b, which had previously been assigned to RNO10, the map position was corrected to RNO6q12-q13.     

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.     

Duplications and deletions of Vh genes in inbred strains of mice

evolution of variable region (Vh) gene family copy number and polymorphism was investigated by the analysis of the immunoglobulin heavy chain variable region (Igh-V) locus in 74 inbred strains and substrains of mice. Several strains were found to have slight differences from Igh-V haplotypes previously identified, usually of a single Vh gene family. These results indicate that the evolution of copy number in the mouse Igh-V locus proceeds largely by the accumulation of incremental changes, reflecting the clustered organization of the mouse Igh-V locus. We have found no evidence of very large or frequent duplication or deletion events indicative of rapid expansion or contraction processes. The existence of one or more particularly large Vh gene families most likely reflects random copy number variation, rather than selection for the amplification of their members. The identification of strains with recombinant Vh gene arrays demonstrates that recombination, both within and between haplotypes, appears to be the predominant mechanism generating the high restriction fragment length polymorphism in the Igh-V locus.Abbreviations used in this paper Igh-V immunoglobulin heavy chain variable region locus - Vh heavy chain variable region gene - Dh heavy chain diversity region gene - V immunoglobulin kappa light chain variable region gene - V T-cell receptor beta chain variable region gene     

Content and organization of the human Ig VH locus: definition of three new VH families and linkage to the Ig CH locus.

We present a detailed analysis of the content and organization of the human immunoglobulin VH locus. Human VH genes representing five distinct families were isolated, including novel members belonging to two out of three of the known VH gene families (VH1 and VH3) as well as members of three new families (VH4, VH5, and VH6). We report the nucleotide sequence of 21 novel human VH genes, many of which belong to the three new VH gene families. In addition, we provide a preliminary analysis of the organization of these gene segments over the full extent of the locus. We find that the five multi-segment families (VH1-5) have members interspersed over nearly the full 1500-2000 kb of the VH locus, and estimate that the entire heavy chain locus covers 2500 kb or less. Finally, we provide the first report of the physical linkage of the variable and constant loci of a human Ig gene family by demonstrating that the most proximal known human VH segments lie within 100 kb of the constant region locus.     

Genetic analysis of the LEW.1AR1-iddm rat: an animal model for spontaneous diabetes mellitus

The LEW.1AR1-iddm/Ztm rat is a new animal model of type 1 diabetes mellitus, which shows an autosomal recessive mode of inheritance for the diabetes-inducing gene. The aim of this study was to define predisposing loci of the diabetic syndrome by linkage analysis using microsatellite markers. A backcross population of 218 rats (BN × LEW.1AR1-iddm) × LEW.1AR1-iddm was analyzed using 157 polymorphic microsatellite markers covering the entire genome. Three genomic regions showed a significant linkage to the diabetic syndrome. The first susceptibility locus on rat Chromosome (RNO) 1 (LOD score 4.13) mapped to the region 1q51–55, which codes for potential candidate genes like Ins1 and Nkx2-3. The second susceptibility locus was also localized on RNO1 in the centromeric region 1p11 (LOD score 2.7) encompassing the Sod2 gene. The third quantitative trait loci (LOD score 2.97) was located on RNO20 within the major histocompatibility complex region. Comparative mapping revealed that the homologous regions in the human genome contain the IDDM loci 1, 5, 8, and 17. The identification of diabetes susceptibility regions of the genetically uniform LEW.1AR1-iddm rat strain will pave the way toward a detailed characterization of the loci conferring diabetes development as well as their functional relevance for the pathogenesis of type 1 diabetes mellitus.     

Homoeologous set of NBS-LRR genes located at leaf and stripe rust resistance loci on short arms of chromosome 1 of wheat

Homoeologous group 1 chromosomes of wheat contain important genes that confer resistance to leaf, stem and stripe rusts, powdery mildew and Russian wheat aphid. A disease resistance gene analog encoding nucleotide binding site-leucine rich repeat (NBS-LRR), designated RgaYr10, was previously identified at the stripe rust resistant locus, Yr10, located on chromosome 1BS distal to the storage protein, Gli-B1 locus. RgaYr10 identified gene members in the homoeologous region of chromosome 1DS cosegregating with the leaf rust resistance gene, Lr21, which originally was transferred from a diploid D genome progenitor. Four RgaYr10 gene members were isolated from chromosome 1DS and compared to two gene members previously isolated from the chromosome 1BS homeologue. NBS-LRR genes tightly linked to stripe rust resistance gene Yr10 on chromosome 1BS were closely related in sequence and structure to NBS-LRR genes tightly linked to leaf rust resistance gene Lr21 located within the homoeologous region on chromosome 1DS. The level of sequence homology was similar between NBS-LRR genes that were isolated from different genomes as compared to genes from the same genome. Electronic Publication     

Sucrose Phosphate Synthase Genes in Plants Belong to Three Different Families

We present phylogenetic analyses to demonstrate that there are three families of sucrose phosphate synthase (SPS) genes present in higher plants. Two data sets were examined, one consisting of full-length proteins and a second larger set that covered a highly conserved region including the 14-3-3 binding region and the UDPGlu active site. Analysis of both datasets showed a well supported separation of known genes into three families, designated A, B, and C. The genomic sequences of Arabidopsis thaliana include a member in each family: two genes on chromosome 5 belong to Family A, one gene on chromosome 1 to Family B, and one gene on chromosome 4 to Family C. Each of three Citrus genes belong to one of the three families. Intron/exon organization of the four Arabidopsis genes differed according to phylogenetic analysis, with members of the same family from different species having similar genomic organization of their SPS genes. The two Family A genes on Arabidopsis chromosome 5 appear to be due to a recent duplication. Analysis of published literature and ESTs indicated that functional differentiation of the families was not obvious, although B family members appear not to be expressed in roots. B family genes were cloned from two Actinidia species and southern analysis indicated the presence of a single gene family, which contrasts to the multiple members of Family A in Actinidia. Only two family C genes have been reported to date. Received: 17 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.     

On the genomics of immunoglobulins in the gray,short-tailed opossum <Emphasis Type="Italic">Monodelphis domestica</Emphasis>

Annotated maps of the IGH, IGK, and IGL loci in the gray, short-tailed opossum Monodelphis domestica were generated from analyses of the available whole genome sequence for this species. Analyses of their content and organization confirmed a number of previous conclusions based on characterization of complementary DNAs encoding opossum immunoglobulin heavy and light chains and limited genomic analysis, including (a) the predominance of a single immunoglobulin heavy chain variable region (IGHV) subgroup and clan, (b) the presence of a single immunoglobulin (Ig)G subclass, (c) the apparent absence of an IgD, and (d) the general organization and V gene complexity of the IGK and IGL light chain loci. In addition, several unexpected discoveries were made including the presence of a partial V to D, germline-joined IGHV segment, the first germline-joined Ig V gene to be found in a mammal. In addition was the presence of a larger number of IGKV subgroups than had been previously identified. With this report, annotated maps of the major histocompatibility complex, T-cell receptor, and immunoglobulin loci have been completed for M. domestica, the only non-eutherian mammalian species for which this has been accomplished, strengthening the utility of this species as a model organism. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

MouseV k gene classification by nucleic acid sequence similarity

Analyses of immunoglobulin (Ig) variable (V) region gene usage in the immune response, estimates ofV gene germline complexity, and other nucleic acid hybridization-based studies depend on the extent to which such genes are related (i. e., sequence similarity) and their organization in gene families. While mouseIgh heavy chainV region (V _H ) gene families are relatively well-established, a corresponding systematic classification ofIgk light chainV region (V _k ) genes has not been reported. The present analysis, in the course of which we reviewed the known extent of theV _k germline gene repertoire andV _k gene usage in a variety of responses to foreign and self antigens, provides a classification of mouseV _k genes in gene families composed of members with >80% overall nucleic acid sequence similarity. This classification differed in several aspects from that ofV _H genes: only someV _k gene families were as clearly separated (by >25% sequence dissimilarity) as typicalV _H gene families; mostV _k gene families were closely related and, in several instances, members from different families were very similar (>80%) over large sequence portions; frequently, classification by nucleic acid sequence similarity diverged from existing classifications based on amino-terminal protein sequence similarity. Our data have implications forV _k gene analyses by nucleic acid hybridization and describe potentially important differences in sequence organization betweenV _H andV _k genes.     

Two large Arabidopsis thaliana gene families are homologous to the Brassica gene superfamily that encodes pollen coat proteins and the male component of the self-incompatibility response

The male component of the self-incompatibility response in Brassica has recently been shown to be encoded by the S locus cysteine-rich gene (SCR). SCR is related, at the sequence level, to the pollen coat protein (PCP) gene family whose members encode small, cysteine-rich proteins located in the proteo-lipidic surface layer (tryphine) of Brassica pollen grains. Here we show that the Arabidopsis genome includes two large gene families with homology to SCR and to the PCP gene family, respectively. These genes are poorly predicted by gene-identification algorithms and, with few exceptions, have been missed in previous annotations. Based on sequence comparison and an analysis of the expression patterns of several members of each family, we discuss the possible functions of these genes. In particular, we consider the possibility that SCR-related genes in Arabidopsis may encode ligands for the S gene family of receptor-like kinases in this species.     

Immunoglobulin light-chain genes in the rhesus macaque I: kappa light-chain germline sequences for subgroups IGKV1, IGKV and IGKV3

The combined processes of immunoglobulin (IG) gene rearrangement and somatic hypermutation allow for the creation of an extremely diverse antibody repertoire. Knowledge of the germline sequence of the IG genes is required so that hypermutation and the affinity matured humoral response can be properly studied. Variable region genes can be arranged into subgroups; in humans, there are 11 IGLV subgroups and six IGKV subgroups. The rhesus macaque (Macaca mulatta) is a relevant non-human primate model for human immunological systems. A number of macaque IGHV, IGHD and IGHJ genes have already been reported, but only one light-chain germline gene has been published so far. Here we report the isolation of new macaque IGKV genes by polymerase chain reaction (PCR) amplification from macaque genomic DNA using primers based on the human sequences. Twenty-eight IGKV1, 22 IGKV2 and 12 IGKV3 germline genes for the macaque were found, the open reading frames of which exhibit high homology to their human counterparts (>96, >99 and >96%, respectively).Supplementary material is available in the online version of this article at .Nucleotide sequence data reported are available in the DDBJ/EMBL/GenBank databases under accession numbers AY963709-AY963773.Authors W.A. Howard and J.M. Bible contributed equally to this work.     

Genomic structure, mapping, and expression analysis of the mammalian Lunatic, Manic, and Radical fringe genes

The three members of the mammalian fringe gene family, Manic fringe (Mfng), Radical fringe (Rfng), and Lunatic fringe (Lfng), were identified on the basis of their similarity to Drosophila fringe (fng) and their participation in the evolutionarily conserved Notch receptor signaling pathway. Fringe genes encode pioneer secretory proteins with weak similarity to glycosyltransferases. Both expression patterns and functional studies support an important role for Fringe genes in patterning during embryonic development and an association with cellular transformation. We have now further characterized the expression and determined the chromosomal localization and genomic structure of the mouse Mfng, Rfng, and Lfng genes; the genomic structure and conceptual open reading frame of the human RFNG gene; and the refined chromosomal localization of the three human fringe genes. The mouse Fringe genes are expressed in the embryo and in adult tissues. The mouse and human Fringe family members map to three different chromosomes in regions of conserved synteny: Mfng maps to mouse Chr 15, and MFNG maps to human Chr 22q13.1 in the region of two cancer-associated loci; Lfng maps to mouse Chr 5, and LFNG maps to human Chr 7p22; Rfng maps to mouse Chr 11, and RFNG maps to human Chr 17q25 in the minimal region for a familial psoriasis susceptibility locus. Characterization of the genomic loci of the Fringe gene family members reveals a conserved genomic organization of 8 exons. Comparative analysis of mammalian Fringe genomic organization suggests that the first exon is evolutionarily labile and that the Fringe genes have a genomic structure distinct from those of previously characterized glycosyltransferases. Received: 19 February 1999 / Accepted: 22 February 1999     

Genomic screening by 454 pyrosequencing identifies a new human IGHV gene and sixteen other new IGHV allelic variants

Complete and accurate knowledge of the genes and allelic variants of the human immunoglobulin gene loci is critical for studies of B cell repertoire development and somatic point mutation, but evidence from studies of VDJ rearrangements suggests that our knowledge of the available immunoglobulin gene repertoire is far from complete. The reported repertoire has changed little over the last 15 years. This is, in part, a consequence of the inefficiencies involved in searching for new members of large, multigenic gene families by cloning and sequencing. The advent of high-throughput sequencing provides a new avenue by which the germline repertoire can be explored. In this report, we describe pyrosequencing studies of the heavy chain IGHV1, IGHV3 and IGHV4 gene subgroups in ten Papua New Guineans. Thousands of 454 reads aligned with complete identity to 51 previously reported functional IGHV genes and allelic variants. A new gene, IGHV3-NL1*01, was identified, which differs from the nearest previously reported gene by 15 nucleotides. Sixteen new IGHV alleles were also identified, 15 of which varied from previously reported functional IGHV genes by between one and four nucleotides, while one sequence appears to be a functional variant of the pseudogene IGHV3-25. BLAST searches suggest that at least six of these new genes are carried within the relatively well-studied populations of North America, Europe or Asia. This study substantially expands the known immunoglobulin gene repertoire and demonstrates that genetic variation of immunoglobulin genes can now be efficiently explored in different human populations using high-throughput pyrosequencing.