A multigene deletion within the immunoglobulin heavy-chain region

The immunoglobulin heavy-chain genes are located in a cluster on chromosome 14. The simultaneous absence of the human IgG1, IgG2, IgG4, and IgA1 subclasses was previously reported in a healthy Tunisian Berber and was later shown to be due to a multigene deletion. We now describe a serological and molecular study of a different deletion observed in a healthy Tunisian. Blot hybridization analysis of the proband's DNA using gamma, epsilon, alpha, and mu switch probes showed that the deletion involves a large region of the immunoglobulin heavy-chain gene cluster: C psi epsilon, C alpha 1, C psi gamma, C gamma 2, and C gamma 4. Incidentally, we showed that the restriction enzyme EcoRI alone can be used with the alpha probe to differentiate A2m types. The deletion described, present in a person homozygous for GM-Am haplotypes (Gm1,17;..;5,14,11,13,10 A2m2), is consistent with previous location, by association analysis, of C psi gamma between C alpha 1 and C gamma 2. There is evidence to suggest that deletions may be more common in the Mediterranean region than in North American Caucasians.     

Nonuniform linkage disequilibrium within a 1,500-kb region of the human immunoglobulin heavy-chain complex.

We have characterized 10 VH polymorphic loci of the VH2, VH3, VH4, and VH5 families. Eight of 10 VH polymorphisms were found to be insertion/deletion polymorphisms, probably the result of nonhomologous recombination over the course of evolution of the current human VH repertoire. The 10 VH polymorphic loci were analyzed in 10 three-generation and 10 two-generation Canadian caucasoid families. Linkage disequilibrium (allelic association) was measured between pairs of VH polymorphic loci, and 12 significant associations were found. The degree of linkage disequilibrium measured between IGH polymorphic loci was then compared with the physical distance separating the loci. The physical distance between IGH polymorphic loci does not entirely determine the degree of linkage disequilibrium between polymorphic loci. Two regions, one in the VH region (between VH3f-2 and VH5-2 and one in the CH region (between C delta and C gamma 3), were found to have linkage disequilibrium values approximately 1/3,000 of that observed in other portions of the IGH region. The previous identification of recombinants in the C delta-to C gamma 3 region indicates that these areas of low linkage disequilibrium are consistent with the presence of recombination hot spots. The observed high amount of recombination in the subtelomeric portion of chromosome 14 therefore appears to be the result of specific hot spots for recombination, rather than a general increase in recombination in this region.     

Pulsed-field gel analysis of human immunoglobulin heavy-chain constant region gene deletions reveals the extent of unmapped regions within the locus

The human immunoglobulin heavy-chain constant region gene locus is organized in three main gene groups, the physical distances of which are unknown. Different types of gene deletions, originated by unequal crossingover, have been found to encompass one or more genes in the locus. We have analyzed some of these deletions by means of pulsed-field gel electrophoresis, which allows resolution of large DNA fragments. By identifying a fragment containing two of the main gene groups and by observing the size reduction of this fragment in subjects with deletions, we were able to estimate the distance between the two groups and better locate the pseudogene in-between.     

Complex protein binding within the mouse immunoglobulin heavy-chain enhancer.

We have begun to purify and characterize several proteins which bind to the mouse immunoglobulin heavy-chain enhancer to understand the molecular interactions important for enhancer activity. Three proteins which bind to different sites on the immunoglobulin heavy-chain enhancer have been chromatographically separated and partially purified. One protein binds a site which has not been reported previously and does not bind to other reported protein-binding sites on the immunoglobulin heavy-chain enhancer. Binding-site boundaries for the three partially purified proteins have been precisely mapped by methylation interference, DNase I footprinting, and orthophenanthroline/copper chemical nuclease footprinting. We have also characterized these three proteins with respect to dissociation rate constants.     

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.     

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     

Genetic analysis of eight loci tightly linked to neurofibromatosis 1

The genetic locus for neurofibromatosis 1 (NF1) has recently been mapped to the pericentromeric region of chromosome 17. We have genotyped eight previously identified RFLP probes on 50 NF1 families to determine the placement of the NF1 locus relative to the RFLP loci. Thirty-eight recombination events in the pericentromeric region were identified, eight involving crossovers between NF1 and loci on either chromosomal arm. Multipoint linkage analysis resulted in the unique placement of six loci at odds greater than 100:1 in the order of pter-A10-41-EW301-NF1-EW207-CRI-L581-CRI-L946 -qter. Owing to insufficient crossovers, three loci--D17Z1, EW206, and EW203--could not be uniquely localized. In this region female recombination rates were significantly higher than those of males. These data were part of a joint study aimed at the localization of both NF1 and tightly linked pericentromeric markers for chromosome 17.     

Gm typing by immunoglobulin heavy-chain gene RFLP analysis.

This study was undertaken to investigate a means of assigning Gm allotypes to Caucasians by RFLP analysis. A single immunoglobulin heavy-chain gamma-4 cDNA probe (HU gamma 4) was hybridized with genomic DNA digested separately with two restriction enzymes, TaqI and PvuII. Results showed excellent correlation (P less than .001) between serologically defined Gm allotypes G1m(1), G1m(2), G2m(23), and G1m;G3m (3;5,10) and RFLPs identified with the (HU gamma 4) probe. We conclude that it is now possible to define common Gm haplotypes in Caucasians by RFLP analysis. This method provides a useful adjunct to serological allotyping and indeed has several important advantages over traditional serology: it allows confident Gm assignment and the definition of homozygous and heterozygous Gm arrangements, is highly reproducible, and is readily executed in any molecular genetic laboratory.     

Nuclear factors binding to the human immunoglobulin heavy-chain gene enhancer.

The human immunoglobulin heavy chain (IgH) gene contains at least two tissue-specific regulatory regions, which are similar to the mouse IgH gene. One is the J-C enhancer and another is located in the 5' promoter region. Using an electrophoretic mobility shift assay and DNase I footprint, we have examined the interaction of factors in B cell nuclear extracts with the two regulatory regions of the human IgH gene. We have identified a nuclear factor in mouse B cell nuclear extracts which bound to specific sequence in the human IgH enhancer. This factor is apparently not present in mouse fibroblast nuclear extracts. We also found factor(s) which bound to the highly conserved octanucleotide sequence within the human IgH enhancer and 5' promoter regions.     

The immunoglobulin heavy-chain variable region in insulin-dependent diabetes mellitus: affected-sib-pair analysis and association studies.

We have analyzed immunoglobulin heavy-chain variable-region (VH) polymorphisms and genetic susceptibility to insulin-dependent diabetes mellitus (IDDM) by using a set of polymorphic loci that span approximately 1,000 kb of the VH region on chromosome 14q32. One hundred one Finnish families with at least two children affected with IDDM were studied. Conventional RFLPs determined by hybridization were used, since no microsatellite repeat markers have been available for this gene region. No evidence for linkage between the VH genes and IDDM could be obtained from haplotype-sharing analysis among the 133 diabetic sib pairs. The frequencies of various VH genotypes were also compared between 101 familial IDDM cases and 114 controls derived from the Finnish background population. The distribution of the genotypes at the VH2-B5 locus was significantly different between these groups (P=.004), the 3.4/3.4 genotype being less common in the IDDM cases. In addition, a different genotype distribution at the VH5-B2 locus was observed in the diabetic subjects (P = .022). When the IDDM cases were stratified by presence or absence of the high-risk HLA-DQB1*0302 allele, no differences in VH genotype frequencies were observed between the 0302-positive and 0302-negative cases. In the transmission test for linkage disequilibrium (TDT), no differences were found between the expected and observed frequencies of the transmitted alleles at the VH2-B5 or VH5-B2 locus. In conclusion, significant differences in VH genotype distributions were observed between the familial IDDM cases and the controls, but the observed associations could not be confirmed by the TDT. Haplotype sharing analysis provided no evidence for genetic linkage between the VH gene region and IDDM.     

Dispersed localization of D segments in the human immunoglobulin heavy-chain locus.

We have studied the organization of the human immunoglobulin heavy-chain genes by pulse field gel electrophoresis as well as by isolation of cosmid clones. The total length of the heavy-chain variable region locus was estimated to be approximately 3000 kb. We found that D segments including a recently isolated D5 segment were dispersed among VH segments. We identified a pseudo V segment 18 kb 3' to the D5 segment in isolated cosmid clones. A 300 kb fragment produced by MluI digestion contained VH, D, JH segments and the distance between VH and D was estimated to be approximately 240 kb. Overlapping cosmid clones containing the human D1, D2, D3, D4, JH, Cmu and C delta genes were isolated. Restriction maps of these regions indicated that the distance between D and JH is about 22 kb. A partial restriction map of the VH locus was constructed using the pulse field gel electrophoresis technique and deletion of VH segments in B cells.     

Regulation of DNA replication within the immunoglobulin heavy-chain locus during B cell commitment

The temporal order of replication of mammalian chromosomes appears to be linked to their functional organization, but the process that establishes and modifies this order during cell differentiation remains largely unknown. Here, we studied how the replication of the Igh locus initiates, progresses, and terminates in bone marrow pro-B cells undergoing B cell commitment. We show that many aspects of DNA replication can be quantitatively explained by a mechanism involving the stochastic firing of origins (across the S phase and the Igh locus) and extensive variations in their firing rate (along the locus). The firing rate of origins shows a high degree of coordination across Igh domains that span tens to hundreds of kilobases, a phenomenon not observed in simple eukaryotes. Differences in domain sizes and firing rates determine the temporal order of replication. During B cell commitment, the expression of the B-cell-specific factor Pax5 sharply alters the temporal order of replication by modifying the rate of origin firing within various Igh domains (particularly those containing Pax5 binding sites). We propose that, within the Igh C(H)-3'RR domain, Pax5 is responsible for both establishing and maintaining high rates of origin firing, mostly by controlling events downstream of the assembly of pre-replication complexes.     

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.     

Sequence dependence of chromosomal R-loops at the immunoglobulin heavy-chain Smu class switch region

The mechanism by which the cytidine deaminase activation-induced deaminase (AID) acts at immunoglobulin heavy-chain class switch regions during mammalian class switch recombination (CSR) remains unclear. R-loops have been proposed as a basis for this targeting. Here, we show that the difference between various forms of the Smu locus that can or cannot undergo CSR correlates well with the locations and detectability of R-loops. The Smu R-loops can initiate hundreds of base pairs upstream of the core repeat switch regions, and the area where the R-loops initiate corresponds to the zone where the AID mutation frequency begins to rise, despite a constant density of WRC sites in this region. The frequency of R-loops is 1 in 25 alleles, regardless of the presence of the core Smu repeats, again consistent with the initiation of most R-loops upstream of the core repeats. These findings explain the surprisingly high levels of residual CSR in B cells from mice lacking the core Smu repeats but the marked reduction in CSR in mice with deletions of the region upstream of the core Smu repeats. These studies also provide the first analysis of how R-loop formation in the eukaryotic chromosome depends on the DNA sequence.     

Functional analysis of defined mutations in the immunoglobulin heavy-chain enhancer in transgenic mice.

We have analyzed the effect of defined mutations in the mouse immunoglobulin heavy-chain enhancer after introduction into the germline of transgenic mice. We have tested a mutation of the enhancer octamer motif, a double mutation of the octamer motif and the microB-site, and a triple mutation in the microE2, microE3 and microE4-sites. All constructs are expressed in the spleen of transgenic mice. Furthermore, expression is exclusively detectable in lymphoid organs and not in several nonlymphoid tissues. Whereas mutations in the microE-sites have a more pronounced effect on transgene activity in thymocytes as compared to bone marrow and spleen cells, the octamer/microB double mutation shows significantly reduced expression levels only in B-cells. Finally, our results demonstrate that the intronic heavy-chain enhancer element does not contribute to the increase steady state levels of heavy-chain mRNA after stimulation of spleen cells with LPS.     

Sequencing of selected regions of the human immunoglobulin heavy-chain gene locus that completes the sequence from JH through the delta constant region.

Much of the nucleotide sequence between the start of the joining region and the end of the immunoglobulin heavy chain delta gene has already been determined. However, two gaps existed in potentially functionally important regions in this sequence: the region between the 3' end of the joining region and the heavy chain enhancer region and that between the enhancer and the mu constant region. We have determined the nucleotide sequences of these regions. The 734 bp between the joining and enhancer regions contained no additional joining regions. The 4525 bp region between the heavy chain enhancer and the mu constant region contains the mu switch region, which consists of pentameric repeats. Approximately 60% of these repeats are GGGCT and GAGCT. With the determination of these sequences, the entire region of the heavy chain locus starting upstream of the joining region to downstream of the last exon of the delta constant region (a total of more than 29 kb) has now been sequenced.     

Genetic analysis of the hypervariable region flanking the human insulin gene.

In order to study the mechanisms for the generation of length diversity within the 5' flanking region of the human insulin gene, we have isolated and sequenced a previously uncharacterized allele. This allele, of a size intermediate between those three already described in the literature, encompasses 1,156 base pairs (bp) and contains 81 reiterated tandem oligonucleotides of 14-15 bp each. Population analysis on 298 independently sampled individuals by Southern blotting of genomic DNA demonstrates that the polymorphic portion of the insulin 5' flanking region varies from 400 to more than 8,000 nucleotides, being encoded by from 30 to over 540 oligomeric repeats. Length variability 5' to the insulin gene is a result primarily of unequal crossing over, which generates an expansion or contraction in the number of tandem repeat units per chromosome. A similar mechanism probably accounts for nondispersed reiterated sequences at other loci in the human genome.     

Copy choice mechanism of immunoglobulin heavy-chain switch recombination.

The immunoglobulin heavy-chain switch is mediated by a recombination event between DNA switch regions associated with donor and recipient constant-region genes. We have determined that the mutations which can be found in some switch regions after recombination appear to arise on only one strand of DNA. This result suggests that switch recombination involves error-prone synthesis of one DNA strand and ligation of the other strand from preexisting DNA.