Polymorphisms of the immunoglobulin heavy-chain delta gene and association with other constant-region genes.

Polymorphisms have previously been reported for the C mu, C alpha, C epsilon, and C gamma genes of the immunoglobulin heavy-chain (IGH) gene cluster. Here we report polymorphisms of the IGH C delta gene region, observed using the enzymes ApaI, AvaII, TaqI, and XbaI. The TaqI and XbaI polymorphisms were used in an investigation of linkage disequilibrium throughout the cluster of constant-region genes. The TaqI polymorphism, located 5' to the C delta gene, is in linkage disequilibrium with a polymorphism of the C mu switch region. The XbaI polymorphism, which is in the vicinity of the C delta 2 exon, is not strongly associated with any other polymorphisms, including the TaqI polymorphism and the Gm polymorphism of C gamma 3. Although there is a high degree of association between most genes of the IGH region, there is a lack of association between C delta and C gamma 3, which may indicate a hot spot for recombination.     

Physical mapping of probes within 14q32, a subtelomeric region showing a high recombination frequency

The genetic linkage map of chromosome 14q32 contains 11 loci which span a distance of more than 60 cM. We have assigned 10 of these loci and the AKT1 proto-oncogene to segments of 14q32, using breakpoints derived from four independent chromosomal deletions or rearrangements. The most telomeric breakpoint was found in a proband (HSC 6) carrying a ring-14 chromosome. HSC 6 is monosomic for the distal part of 14q32, which contains the immunoglobulin heavy-chain locus (IGH), and random markers D14S20, D14S19, and D14S23. Two other chromosomal breakpoints, found in probands HSC 121 and HSC 981, could not be distinguished from each other using DNA probes, although the cytogenetic breakpoints appeared to be different at 14q32.32 and 14q32.31, respectively. The region between the breakpoints of HSC 6 and HSC 121 contains AKT1, D14S1, D14S17, and D14S16. The entire telomeric band 14q32 is assumed to contain about 10% of chromosome 14, or approximately 10 Mb. The 8 most telomeric loci, including D14S1, map to 14q32.32-qter, which measures only several megabases. However, these loci span a genetic distance of 23 cM. The high recombination frequency contrasts with the observation that two of the gamma genes in the IGH constant region show a high degree of linkage disequilibrium, though 180 kb apart. This finding suggests that a telomeric localization per se does not lead to a higher recombination frequency and favors the hypothesis that the higher recombination frequency at the telomeres may be due to specific "hot spots" for recombination.     

Linkage disequilibria between pairs of loci within a highly polymorphic region of chromosome 2Q.

A recently described region on chromosome 2q contains seven restriction fragment length polymorphisms (RFLPs) revealed by single-copy probes isolated from a 20-kilobase (kb) segment of a single cosmid insert. Analysis of six of these loci demonstrates modest amounts of linkage disequilibrium. This reflects the presence of a substantial number of different haplotypes in this chromosome region and indicates that the region could be used as one highly polymorphic locus. No consistent relationship is found between the amount of linkage disequilibrium and the physical distance between pairs of loci. For seven of the 10 pairs of diallelic loci studied, the observed disequilibrium can be attributed primarily to the absence of the minor haplotype from the population. These results suggest that, for small regions of the genome, factors such as mutation, genetic drift, and population admixture may have effects that outweight those of recombination. In addition, results are reviewed which show that estimates of linkage disequilibrium coefficients for tightly linked loci are very imprecise. Thus, the inference of gene order from linkage disequilibrium values must be regarded with caution.     

Linkage disequilibrium analysis of the human adenosine deaminase (ada) gene provides evidence for a lack of correlation between hot spots of equal and unequal homologous recombination

The linkage disequilibrium (LD) pattern within the adenosine deaminase (ADA) gene was analyzed by studying 13 polymorphic loci in 137 families from two European and three African populations. Evidence for the presence of a 12-kb meiotic crossover hot spot, spanning part of the first and the second intron and flanked by regions of reduced recombination activity, was obtained. Moreover, segregation analysis of 113 informative meioses revealed two recombination events that are internal or overlap the 12-kb region, thus suggesting a recombination rate for the hot-spot region about 50-fold higher than the mean rate across the human genome. Within the hot spot, a 144-bp palindromic sequence was also identified and its possible involvement in the recombination process is discussed. The 12-kb region characterized by the low degree of LD does not include the 3.2-kb region that is deleted, as a result of recurrent unequal homologous recombination between two Alu elements, in patients affected by autosomal severe combined immunodeficiency. This observation provides the first evidence for an absence of correlation between hot spots of equal and unequal homologous recombination.     

Localization and genetic linkage of the human immunoglobulin heavy chain genes and the creatine kinase brain (CKB) gene: Identification of a hot spot for recombination

The immunoglobulin heavy chain (IGH) gene cluster and the gene coding for the brain form of the enzyme creatine kinase (CKB) have previously been localized to chromosome 14, at 14q32.3 and 14q32, respectively. Here we report more precise regional localization of these genes by dosage studies using DNA from a child hemizygous for the region from 14q32.32 to 14qter. CKB and IGH are present in a single dose in the proband. Dosage studies in a second patient with a similar but smaller deletion due to a ring chromosome 14 show that CKB is proximal to the IGH cluster. An EcoRI restriction site polymorphism was found with probes for the CKB gene. Linkage analysis of family data indicates that CKB is closely linked to IGH. Linkage analysis also revealed unusually high recombination (beta = 3.2%) between the C delta and C gamma 3 genes of the IGH constant region, which are only 60 kb apart. This finding, in combination with a previous observation of linkage equilibrium in the region, suggests that the C delta-C gamma 3 region contains a recombination hot spot.     

The immunoglobulin heavy chain and disease association: application to pemphigus vulgaris

Genes involved in the immune response are generally encoded from a complex cluster of gene segments. Studies of the association of diseases with such genes require well-defined genetic markers throughout the selected region. A set of 15 polymorphic loci that span 1500 kb of the immunoglobulin heavy chain (IGH) complex, 8 in the variable (VH) region and 7 in the constant (CH) region, were selected for the study of disease association. We present a protocol for the use of multiple immunoglobulin heavy chain (IGH) polymorphisms for general application in disease association studies. No microsatellite repeat markers are available for this region. To demonstrate the applicability of this approach, we have examined these IGH polymorphisms in families with individuals affected with pemphigus vulgaris (PV), an autoimmune dermatologie disease. Allele frequencies in 12 patients with PV were compared with those found in their spouses, and with those in a white Canadian control population. A significant difference was found between PV patients and both control groups for the presence of the VH gene VH3f-R4, and possibly for the absence of VH3f-R3, suggesting the possibility of susceptibility factors in these regions. Examination of the frequencies of the IGH region C-haplotypes of PV patients indicated that, while the patients did not differ significantly from their spouses (²=1.79), both groups were found to differ significantly from the white Canadian control group (²=10.10), emphasizing the importance of matching the ethnic background of controls with that of the patient test group in disease association studies. Unexpectedly, two patients had large deletions of genes in the IGH constant region, which could play a role in the development of PV and require further investigation.     

Analysis of genetic variation reveals human immunoglobulin VH-region gene organization

We have investigated the extent of genetic variation and the number of germ-line heavy-chain-variable (VH) genes to obtain information on the organization and repertoire of the VH genes. Our studies revealed extensive genetic variation in this region, indicated by restriction-endonuclease site polymorphisms. Analysis of the distribution of selected polymorphic loci revealed evidence of linkage disequilibrium, particularly between VH2 and VH3 subclass loci, indicating that the subclasses are interdispersed in the human germ-line chromosome. Absolute correlation was detected between alleles of a VH2 locus and the alleles of three VH3 loci, evidence for an extra set of VH genes, which are present in 48% of the Caucasian population. A preliminary estimate of the number of VH genes, approximately 50, indicates a smaller number of VH genes than suggested by the amount of protein variation. The extensive genetic variation we have observed may be associated with genetic differences in the immune response and potentially with variable susceptibility to autoimmune disorders.     

Genetic analysis of eight linked polymorphisms within the human immunoglobulin heavy-chain region.

Genetic analyses of multiple restriction fragment length polymorphisms, revealed by a single DNA probe containing the switch region of the immunoglobulin constant heavy-chain (IgCH) mu gene, are presented here in detail. Five of the polymorphic loci segregate in complete linkage with IgCH allotypic markers, while one appears to be located at more than 10 centimorgans from the IgCH region. A study of over 100 random haplotypes typed at eight linked loci, including the Ig switch polymorphisms and the classical Gm-Am allotypes, allowed us to construct an evolutionary tree by which each haplotypic variant can be derived one from the other either by single-step mutation or by recombination. A few of the recombinant haplotypes appeared to carry large DNA duplications that could be explained by unequal crossing over; others might postulate gene-conversion events. Linkage disequilibria observed between the IgCH-linked loci were compared with expected ones. A heterogeneous distribution of recombination rates is clearly documented, a "hot" region of recombination being present between the gamma 2 and switch alpha 2 loci.     

Recombination and gene conversion in a 170-kb genomic region of Arabidopsis thaliana

Arabidopsis thaliana is a highly selfing plant that nevertheless appears to undergo substantial recombination. To reconcile its selfing habit with the observations of recombination, we have sampled the genetic diversity of A. thaliana at 14 loci of approximately 500 bp each, spread across 170 kb of genomic sequence centered on a QTL for resistance to herbivory. A total of 170 of the 6321 nucleotides surveyed were polymorphic, with 169 being biallelic. The mean silent genetic diversity (pi(s)) varied between 0.001 and 0.03. Pairwise linkage disequilibria between the polymorphisms were negatively correlated with distance, although this effect vanished when only pairs of polymorphisms with four haplotypes were included in the analysis. The absence of a consistent negative correlation between distance and linkage disequilibrium indicated that gene conversion might have played an important role in distributing genetic diversity throughout the region. We tested this by coalescent simulations and estimate that up to 90% of recombination is due to gene conversion.     

Linkage disequilibrium analyses and restriction mapping of four RFLPs at the proα2(I) collagen locus: lack of correlation between linkage disequilibrium and physical distance

Summary Restriction fragment length polymorphisms (RRLPs) located at short distances may demonstrate linkage disequilibrium. Under the assumption that the distances between the loci of the RFLPs are inversely related to the linkage disequilibria, gene order may be deduced. However, if the assumption is invalid, the results may be incorrect. We have studied four different DNA polymorphisms at the COLIA2 locus in 180 unrelated Norwegian individuals. Observed frequencies (presence/absence) for the different polymorphic sites were as follows: site A (EcoRI) 0.30/0.70, site B (MspI) 0.83/0.16, site C (StuI) 0.86/0.14, and site D (RsaI) 0.66/0.34. Of 16 possible haplotypes 12 were demonstrated, and 2 additional were deduced to be present. Restriction mapping of the four polymorphic sites gave the following order of the sites from the 5 to the 3 of the gene: A-D-B-C. Linkage disequilibrium was not found between the sites A and D; strong disequilibrium was found between sites A and C, and B and C; and less strong, between A and B, B and D, and C and D. Analysis of linkage disequilibrium coefficients between all pairs of loci demonstrated that there is no consistent relationship between linkage disequilibrium and physical distance (=-0.07). These results suggest that for a small region of the genome, factors such as deviating mutation rate and gene conversion may add significantly to rearrangements by recombination. Thus, a deduced gene order from linkage disequilibrium data has to be regarded with great caution.     

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.     

Prenatal diagnosis and linkage disequilibrium with cystic fibrosis for markers surrounding D7S8

Summary Three polymorphic DNA markers surrounding the D7S8 locus were tested for their usefulness in the diagnosis of cystic fibrosis (CF) by linkage analysis. The markers correspond to the loci D7S424 and D7S426. These polymorphisms were studied by centers in the U.S., the United Kingdom, the Netherlands, and Italy, using samples from populations throughout Europe and North America. The additional information provided by these probes increased the heterogeneity of the region from 50% to 58% and was essential for a completely informative diagnosis in one family. A very high degree of linkage disequilibrium was found between these markers, which span a distance of approximately 250kb. In addition, linkage disequilibrium with CF was noted. Significant heterogeneity of linkage disequilibrium was found among the populations, both for the marker-marker pairs and between the markers and CF.     

A high-resolution linkage map for the Z chromosome in chicken reveals hot spots for recombination

A comprehensive linkage map for chicken chromosome Z was constructed as the result of a large-scale screening of single nucleotide polymorphisms (SNPs). A total of 308 SNPs were assigned to Z based on the genotype distribution among 182 birds representing several populations. A linkage map comprising 210 markers and spanning 200.9 cM was established by analyzing a small Red junglefowl/White Leghorn intercross. There was excellent agreement between the linkage map for Z and a recently released assembly of the chicken genome (May 2006). Almost all SNPs assigned to chromosome Z in the present study are on Z in the new genome assembly. The remaining 12 loci are all found on unassigned contigs that can now be assigned to Z. The average recombination rate was estimated at 2.7 cM/Mb but there was a very uneven distribution of recombination events with both cold and hot spots of recombination. The existence of one of the major hot spots of recombination, located around position 39.4 Mb, was supported by the observed pattern of linkage disequilibrium. Thirteen markers from unassigned contigs were shown to be located on chromosome W. Three of these contigs included genes that have homologues on chromosome Z. The preliminary assignment of three more genes to the gene-poor W chromosome may be important for studies on the mechanism of sex determination and dosage compensation in birds.     

DNA haplotype analysis suggests linkage disequilibrium in the human insulin receptor gene

Summary Haplotypes of the insulin receptor gene were resolved in parents from Scandinavian nuclear families by studying the segregation of seven restriction fragment length polymorphisms (RFLPs). Of 97 unrelated parents, 41 had non-insulin-dependent diabetes mellitus (NIDDM). Considerable linkage disequilibrium in the region of the insulin receptor gene was found. Pairwise non-random associations were found between proximate RFLP sites, indicating the absence of recombinational hot spots between these sites. Thus, association studies between DNA polymorphisms at this locus and disease susceptibility genes could well be feasible in this population. Differences in the distribution of insulin receptor haplotypes were examined between NIDDM patients and healthy subjects. However, the differences observed were not statistically significant.     

Linkage disequilibrium detected between dystrophia myotonica and APOC2 locus in the Finnish population

Summary Three polymorphic loci APOC2, CKMM and p134C were used to haplotype 15 Finnish dystrophia myotonica (DM) families representing about one third of all DM patients in this isolated population. Compound APOC2 and CKMM haplotypes reveal linkage disequilibrium: 90% of DM chromosomes co-occur with the haplotypes that occur in 31% of normal chromosomes only. The same disequilibrium is present when only polymorphisms occurring at the APOC2 locus are used. Surprisingly, no statistically significant linkage disequilibrium was discovered at the CKMM locus alone. Of the meiotic events, 84% were informative when both APO2 and CKMM loci were used. When studied selectively, 60% of meiotic events were informative at the APOC2 locus, whereas CKMM alone resulted in 65% meiotic informativeness. The distal marker p134C was found to have an unfortunately low information content in our population.     

Polymorphism of the human Ig VH4 gene family.

In this study, we analyze the human VH4 gene family and find it to exhibit a level of polymorphism similar to that of the much larger VH3 family. A cloned VH4 probe detected an average of 10 hybridizing BgIII restriction fragments in genomic DNA derived from 75 unrelated individuals and a total of 15 distinct bands. Of these 15 restriction fragments, 12 were polymorphic, as demonstrated by band absence in some individuals. Oligonucleotide probes specific to CDR1 and CDR2 sequences of known VH4 genes detected limited numbers of bands and revealed sequence polymorphisms that correlated with several of the RFLP detected by the cloned probe. The prevalence of the individual polymorphic restriction fragments was highly variable, ranging from 1% to 97%, with a mean prevalence of 51%. These values resemble those previously observed among VH3 elements. Analysis of linkage disequilibrium suggests that most VH4 gene segments are in genetic equilibrium. These results indicate that the VH4 loci, like those of VH3, are dominated by relatively few, perhaps two to four, alleles/locus and further suggest that the haplotype organization of the human VH locus is very complex.     

Allelic recombination and linkage disequilibrium within Msp-1 of Plasmodium falciparum, the malignant human malaria parasite

The C-terminal, cysteine-rich 19kDa domain of merozoite surface protein-1 (MSP-1) of Plasmodium falciparum is a target of the host's humoral immunity and thus a malaria vaccine candidate. Although variation in the 19kDa domain is limited among parasite isolates, tertiary structure-dependent intramolecular associations between the 19kDa domain and other parts of MSP-1 are suggested to be involved in immune evasion by allowing competitive binding of protective and non-protective antibodies directed to their epitopes, which are conformationally in close proximity but separated at the primary structure. Since allelic recombination can account for the major variability of the Msp-1 gene, we examined whether linkage disequilibrium occurs between polymorphic loci in the 5'- and the 3'-region, the latter encoding the 19kDa domain. From 184 Thai field isolates, we selected 69 isolates with a single allelic type in six variable blocks of Msp-1 as determined by PCR-based allelic typing. All the isolates showed no evidence of recombination in blocks 6 to 16, whereas recombination was apparent in blocks 2 to 6. Sequencing of the 3'-region revealed two potential recombination sites in block 17. Strong linkage disequilibrium was seen between polymorphic loci in the 5'- and 3'-regions. The strength of this disequilibrium did not correlate with distance between loci. We discuss the possible role of epistatic selection on particular association types (haplotypes) of Msp-1.     

Tight linkage of pyruvate kinase (PKLR) and glucocerebrosidase (GBA) genes

Two polymorphisms, one in the liver-type pyruvate kinase gene (PKLR) and one in the glucocerebrosidase gene (GBA), both of which are on band q21 of chromosome 1, were found to be tightly linked. Each of three Gaucher disease mutations in 112 chromosomes studied was associated with a unique haplotype. With a conservative assumption about the length of time that the Gaucher disease mutation has been present in the Jewish population, we deduce that the genetic distance between these two loci is probably under 0.2 centimorgans. Four haplotypes are produced by these polymorphic loci, but two of these are relatively uncommon because the polymorphic sites are in linkage disequilibrium. Nonetheless these markers are potentially useful in the prenatal diagnosis of pyruvate kinase deficiency in families who have at least one affected child and may also be helpful in heterozygote detection in families with Gaucher disease where a specific mutation producing the disease in unknown.     

Carrier detection in haemophilia B using two further intragenic restriction fragment length polymorphisms.

A normal human population has been screened for the existence of further restriction fragment length polymorphisms (RFLPs) in the clotting factor IX gene in addition to the TaqI polymorphism already characterised (1,2). Two polymorphic loci were found, both within 6 Kb of the TaqI polymorphism within the body of the factor IX gene. One of the polymorphisms has been shown to be due to either the presence or absence of a particular recognition site for the restriction enzyme XmnI. The other, visualised as a difference in fragment pattern produced by digestion with either HinfI or DdeI, has two allelic forms differing by a 50 bp element of inserted DNA. Sequence analysis has shown the inserted element to be in a region of Z type DNA sequence, the insertion representing a duplication of flanking sequence on either side. The two polymorphisms are inherited in simple Mendelian fashion and have both been used to diagnose haemophilia B carrier status. It is estimated that the combined use of these polymorphisms in the factor IX gene, despite linkage disequilibrium between the 3 polymorphic loci, should enable carrier status to be determined in approximately 66% of all haemophilia B families.     

How can we distinguish between mutational "hot spots" and "old sites" in human mtDNA samples?

New research into variation in mutation rates across nucleotide positions in human mitochondrial DNA (mtDNA) calls into question population genetics models that assume a constant mutation rate for all sites in a sequence, particularly for hypervariable control region segments I and II. Related to this research is discovering the extent to which highly polymorphic sites are really mutational "hot spots" rather than "old" sites rooted early in the phylogenetic tree. This issue is addressed through the analysis of linkage disequilibrium patterns in the mtDNAs of 10 human populations. Hot spots can be expected to show little or no disequilibrium since they can be interpreted as having randomly expressed patterns. In fact, the results suggest that many highly polymorphic sites are not old sites, but instead are hot spots. Suspected hot spots are listed and compared with hypervariable sites given by Wakeley (1993) and Hasegawa et al. (1993).