Linkage disequilibrium in the neurofibromatosis 1 (NF1) region: implications for gene mapping.

To test the usefulness of linkage disequilibrium for gene mapping, we compared physical distances and linkage disequilibrium among eight RFLPs in the neurofibromatosis 1 (NF1) region. Seven of the polymorphisms span most of the NF1 gene, while the remaining polymorphism lies approximately 70 kb 3' to a stop codon in exon 49. By using Centre d'Etude du Polymorphisme Humain (CEPH) kindreds, 91-110 unrelated parents were genotyped. A high degree of disequilibrium is maintained among the seven intragenic polymorphisms (r > .82, P < 10(-7)), even though they are separated by as much as 340 kb. The 3' polymorphism is only 68 kb distal to the next polymorphism, but disequilibrium between the 3' polymorphism and all others is comparatively low (magnitude of 4 < .33, P values .27-.001). This result was replicated in three sets of unrelated kindreds: the Utah CEPH families, the non-Utah CEPH families, and an independent set of NF1 families. Trigenic, quadrigenic, three-locus, and four-locus disequilibrium measures were also estimated. There was little evidence of higher-order linkage disequilibrium. As expected for a disease with multiple mutations, no disequilibrium was observed between the disease gene and any of the RFLPs. The observed pattern of high disequilibrium within the gene and a loss of disequilibrium 3' to the stop codon could have implications for gene mapping studies. These are discussed, and guidelines for linkage disequilibrium studies are suggested.     

Tight linkage between type III Gaucher's disease (Norrbottnian type) and a MspI polymorphism within the gene for human glucocerebrosidase

A MspI polymorphism was detected in the beta-glucocerebrosidase gene in 10 Swedish families affected by type III Gaucher's disease. The sizes of the polymorphic fragments were 1.70 and 1.75 kb and the disease was found to segregate with the 1.70-kb fragment in 32 meioses. Only the 1.75-kb fragment was detected in families with no history of Gaucher's disease. The results indicate that the mutation causing type III Gaucher's disease has occurred once within the Swedish population. The polymorphism is useful for carrier detection since biochemical tests sometimes give inconclusive results.     

RFLPs for Duchenne muscular dystrophy cDNA clones 9 and 10.

The complete 14-kb cDNA for the gene causing the X-linked recessive muscular dystrophy (MD) type Duchenne (DMD) and Becker (BMD) has recently been cloned and made available for deletion/duplication screening in patients. It detects 65 exon-containing nonpolymorphic HindIII fragments spread over a gene locus of about 2,000 kb. When the entire DMD cDNA is used, deletions/duplications can be found in about 65%-70% of affected patients, permitting direct carrier detection by densitometric scanning. But in cases where no deletion/duplication is detectable, RFLP analysis, specially favored within the gene, will be the method for carrier-status determination. Clones 9 and 10-1.2-kb and 0.7-kb fragments, respectively, of the 14-kb DMD cDNA--have been hybridized with human genomic DNA digested by nine different restriction enzymes. Five RFLPs, involving Asp700, PvuII, XbaI, and EcoRV sites, were detected, and Mendelian inheritance could be demonstrated. Since clones 9 and 10 are localized telomeric to the mutation-hot-spot region, their polymorphisms are thought to be very helpful as flanking markers for indirect carrier detection in families with a family history of DMD/BMD. Moreover, these RFLPs can be used for direct carrier detection or exclusion in families with patients showing a deletion/duplication in the region of p9 or p10.     

A method for developing high-density SNP maps and its application at the type 1 angiotensin II receptor (AGTR1) locus

Evaluating the potential genetic components of complex disease will likely be aided through the use of dense polymorphism maps. Previously, we reported evidence for linkage with diabetic nephropathy on chromosome 3q in a region encompassing the type 1 angiotensin II receptor (AGTR1) gene. To further investigate any role for this gene in disease onset, we set out to design a dense polymorphism map spanning the AGTR1 locus for the purpose of association studies. Toward this goal, we have developed a technique for rapid identification of polymorphisms in long stretches of genomic DNA. This approach uses long-range PCR, DNA pooling, and transposon-based DNA sequencing. Using this technique, we efficiently validated and genotyped 18 polymorphisms spanning the 60.5-kb AGTR1 locus. Our panel of polymorphisms has an average spacing of 3.2 kb and an average minor allele frequency of 24%.     

KpnI families of long, interspersed repetitive DNAs associated with the human β-globin gene cluster

KpnI families of long, interspersed repetitive DNAs are ubiquitous repetitive elements that occur in tens of thousands of copies in primate genomes. KpnI 1.2, 1.5 and two different KpnI 1.8-kb families were found within and flanking a 6.4-kb repeat beginning at 3 kb, 3' from the human β-globin gene. Thus, six different types of KpnI families have now been identified, and four of these are found next to each other in a specific 6.4-kb repeat. Clones of the distinct KpnI families were hybridized to clones of the 6.4-kb repeat and adjacent sequences encompassed within some 17.6 kb of DNA lying 3' to the β-globin gene cluster. The four KpnI families appear to make up the entire length of the 6.4-kb repeat. The linear order of the various cloned KpnI sequences in the repeat is 5'-pBK(1.8)26-pBK.(1.5)54-pBK(1.2)11-pBK(1.8)11-3'. KpnI 1.2-kb sequences were also detected downstream from the 6.4-kb repeat. As in the case of the KpnI 1.2 and 1.5-kb families, the two KpnI 1.8-kb sequence families described here each hybridized with about 15% of all plaques in two independently generated human genome libraries.     

Unraveling of the polymorphic C lambda 2-C lambda 3 amplification and the Ke+Oz- polymorphism in the human Ig lambda locus

Two polymorphisms of the human Ig(lambda) (IGL) locus have been described. The first polymorphism concerns a single, 2- or 3-fold amplification of 5.4 kb of DNA in the C(lambda)2-C(lambda)3 region. The second polymorphism is the Mcg(-)Ke(+)Oz(-) isotype, which has only been defined via serological analyses in Bence-Jones proteins of multiple myeloma patients and was assumed to be encoded by a polymorphic C(lambda)2 segment because of its high homology with the Mcg(-)Ke(-)Oz(-) C(lambda)2 isotype. It has been speculated that the Mcg(-)Ke(+)Oz(-) isotype might be encoded by a C(lambda) gene segment of the amplified C(lambda)2-C(lambda)3 region. We now unraveled both IGL gene polymorphisms. The amplification polymorphism appeared to result from a duplication, triplication, or quadruplication of a functional J-C(lambda)2 region and is likely to have originated from unequal crossing over of the J-C(lambda)2 and J-C(lambda)3 region via a 2.2-kb homologous repeat. The amplification polymorphism was found to result in the presence of one to five extra functional J-C(lambda)2 per genome regions, leading to decreased Ig(kappa):Ig(lambda) ratios on normal peripheral blood B cells. Via sequence analysis, we demonstrated that the Mcg(-)Ke(+)Oz(-) isotype is encoded by a polymorphic C(lambda)2 segment that differs from the normal C(lambda)2 gene segment at a single nucleotide position. This polymorphism was identified in only 1.5% (2 of 134) of individuals without J-C(lambda)2 amplification polymorphism and was not found in the J-C(lambda)2 amplification polymorphism of 44 individuals, indicating that the two IGL gene polymorphisms are not linked.     

Nonuniform recombination within the human beta-globin gene cluster.

Population genetic analysis of 15 restriction site polymorphisms demonstrates nonuniform recombination within the human beta-globin gene cluster. These DNA polymorphisms show two clusters of high nonrandom associations, one 5' and another 3' to the beta-globin structural gene, with no significant linkage disequilibrium between the two clusters. The 5'- and 3'-association clusters are 34.6 kilobases (kb) and 19.4 kb long, respectively, and are separated by 9.1 kb of DNA immediately 5' to the beta-globin gene. For each of these three DNA regions, we have observed a relationship between nonrandom associations and physical distance between the polymorphisms. However, this relationship differed for each of these regions. On the assumption that the effective population size (Ne) is 5,000-50,000, we estimate the total recombination rate to be 0.0017%-0.0002% in the 5' cluster, 0.0931%-0.0093% in the 3' cluster, and 0.2912%-0.0219% in the 9.1-kb region between them. The beta cluster thus shows nonuniformity in recombination. Moreover, the recombination rate in the 9.1-kb DNA segment is 3-30 times greater than expected and is thus a hot spot for meiotic recombination.     

The human fumarylacetoacetase gene: characterisation of restriction fragment length polymorphisms and identification of haplotypes in tyrosinemia type 1 and pseudodeficiency

Summary Deficiency of human fumarylacetoacetase (FAH) activity results in hereditary tyrosinemia type I. Using the restriction enzymes BglII, KpnI and StuI and a 1.3-kb cDNA probe for the FAH gene, we have found 6 restriction fragment length polymorphisms (RFLPs). These RFLPs were utilised in 3 tyrosinemia families in which one or both parents are carriers of both a tyrosinemia and a pseudodeficiency gene for FAH. Full information was obtained in two of these families. The polymorphisms identified 6 haplotypes. The haplotype distribution was significantly different in 32 unrelated tyrosinemia patients compared with a reference population of 100 individuals. The combined polymorphism information content was 0.77.     

Cloning of human heparan sulfate proteoglycan core protein, assignment of the gene (HSPG2) to 1p36.1----p35 and identification of a BamHI restriction fragment length polymorphism.

We have isolated a cDNA coding for the core protein of the large basement membrane heparan sulfate proteoglycan (HSPG) from a human fibrosarcoma cell (HT1080) library. The library was screened with a mouse cDNA probe and one clone obtained, with a 1.5-kb insert, was isolated and sequenced. The sequence contained an open reading frame coding for 507 amino acid residues with a 84% identity to the corresponding mouse sequence. This amino acid sequence contained several cysteine-rich internal repeats similar to those found in component chains of laminin. The HSPG cDNA clone was used to assign the gene (HSPG2) to the p36.1----p35 region of chromosome 1 using both somatic cell hybrid and in situ hybridization. In the study of the polymorphisms of the locus, a BamHI restriction fragment length polymorphism was identified in the gene. This polymorphism displayed bands of 23 and 12 kb with allele frequencies of 76 and 24%, respectively.     

Linkage disequilibrium predicts physical distance in the adenomatous polyposis coli region.

To test the reliability of linkage-disequilibrium analysis for gene mapping, we compared physical distance and linkage disequilibrium among seven polymorphisms in the adenomatous polyposis coli (APC) region on chromosome 5. Three of them lie within the APC gene, and two lie within the nearby MCC (mutated in colon cancer) gene. One polymorphism lies between the two genes, and one is likely to be 5' of MCC. Five of these polymorphisms are newly reported. All polymorphisms were typed in the CEPH kindreds, yielding 179-205 unrelated two-locus haplotypes. Linkage disequilibrium between each pair of polymorphisms is highly correlated with physical distance in this 550-kb region (correlation coefficient -.80, P < .006). This result is replicated in both the Utah and non-Utah CEPH kindreds. There is a tendency for greater disequilibrium among pairs of polymorphisms located within the same gene than among other pairs of polymorphisms. Trigenic, quadrigenic, three-locus, and four-locus disequilibrium measures were also estimated, but these measures revealed much less disequilibrium than did the two-locus disequilibrium measures. A review of 19 published disequilibrium studies, including this one, shows that linkage disequilibrium nearly always correlates significantly with physical distance in genomic regions > 50-60 kb but that it does not do so in smaller genomic regions. We show that this agrees with theoretical predictions. This finding helps to resolve controversies regarding the use of disequilibrium for inferring gene order. Disequilibrium mapping is unlikely to predict gene order correctly in regions < 50-60 kb in size but can often be applied successfully in regions of 50-500 kb or so in size. It is convenient that this is the range in which other mapping techniques, including chromosome walking and linkage mapping, become difficult.     

Nonrandom Distribution of Chloroplast Recombination Events in Chlamydomonas Reinhardtii: Evidence for a Hotspot and an Adjacent Cold Region

Intermolecular recombination of Chlamydomonas chloroplast genes has been analyzed in sexual crosses and following biolistic transformation. The pattern and position of specific exchange events within 15 kb of the 22-kb inverted repeat have been mapped with respect to known restriction fragment length polymorphism markers that distinguish the chloroplast genomes of the interfertile species Chlamydomonas reinhardtii and Chlamydomonas smithii. Recombinant progeny were selected from two- and three-factor crosses involving point mutations conferring herbicide (dr) and antibiotic resistance (er and spr) in the psbA, 23S and 16S ribosomal RNA genes, respectively. Exchange events were not randomly distributed over the 15-kb region, but were found to occur preferentially in a 0.7-kb sequence spanning the 3' end of the psbA gene and were much less common in an adjacent region of ca. 2.0 kb. These findings are corroborated by data showing that the dr mutation is unlinked genetically (3% recombination/kb) to the er and spr rRNA mutations, which are themselves linked and show ca. 1% recombination/kb. This discrepancy is significant since the dr-er and er-spr intervals are about the same length (ca. 7 kb). During chloroplast transformation, the 0.7-kb recombination hotspot also functions as a preferential site for exchange events leading to the integration of donor psbA gene sequences. The 0.7-kb hotspot region contains four classes of 18-37-bp direct repeats also found in other intergenic regions, but no open reading frame. Using deletion constructs in a chloroplast transformation assay, the hotspot was localized to a 500-bp region that lacks most of these repeats, which suggests that the repeats themselves are not responsible for the increased recombination frequency. Within this region, a 400-bp sequence is highly conserved between the chloroplast genomes of C. reinhardtii and C. smithii and includes several structural motifs characteristic of recombination hotspots in other systems.     

The V-J-intergenic region of the human kappa locus

The 23-kb region between the V kappa and J kappa gene clusters was investigated in some detail. The region was found to be free of V kappa genes or V kappa gene-like structures, confirming the previous supposition that the V kappa gene B3 is the J kappa proximal V kappa gene. The B3-J kappa distance of 23 kb was found to be the same in the DNAs of several individuals. A HindIII restriction fragment length polymorphism was detected within this region. A sequence of 533 bp located approximately in the middle of the region has a highly homologous counterpart (called homox) on another chromosome. The two sequences are 96% identical. Possible mechanisms for the generation of such a duplicate are discussed.     

Novel microsatellite markers and single nucleotide polymorphisms refine the tylosis with oesophageal cancer (TOC) minimal region on 17q25 to 42.5 kb: sequencing does not identify the causative gene

Tylosis (focal non-epidermolytic palmoplantar keratoderma) is associated with the early onset of squamous cell oesophageal cancer in three families. Linkage and haplotype analyses have previously mapped the tylosis with oesophageal cancer (TOC) locus to a 500-kb region on chromosome 17q25 that has also been implicated in sporadically occurring squamous cell oesophageal cancer. In the current study, 17 additional putative microsatellite markers were identified within this 500-kb region by using sequence data and seven of these were shown to be polymorphic in the UK and US families. In addition, our complete sequence analysis of the non-repetitive parts of the TOC minimal region identified 53 novel and six known single nucleotide polymorphisms (SNPs) in one or both of these families. Further fine mapping of the TOC disease locus by haplotype analysis of the seven polymorphic markers and 21 of the 59 SNPs allowed the reduction of the minimal region to 42.5 kb. One known and two putative genes are located within this region but none of these genes shows tylosis-specific mutations within their protein-coding regions. Alternative mechanisms of disease gene action must therefore be considered.     

Scanning the first part of the neurofibromatosis type 1 gene by RNA-SSCP: Identification of three novel mutations and of two new polymorphisms

Neurofibromatosis type 1 (NF1) of von Recklinghausen is a common autosomal dominant disorder, characterized by peripheral neurofibromas, café-au-lait spots and Lisch nodules of the iris. The high mutation rate at the NF1 locus results in a wide range of molecular abnormalities. We have scanned 14 different exons from the first part of the NF1 gene using the RNA-single strand conformation polymorphism (RNA-SSCP) method in a series of 40 NF1 patients. Three novel mutations, two nonsense and one missense, and two polymorphisms have been detected in familial cases. Genotype-phenotype correlations have been investigated, but no particular association has been detected. After this screening, the majority of NF1 chromosomes has not yet been characterized, confirming the difficulty in detecting the defect underlying NF1 in most families, even following extensive DNA analysis. Received: 4 August 1995 / Revised: 19 September 1995     

Polymorphic restriction sites in the horse (β-globin gene cluster

Horse DNA samples digested with PstI and probed with the rabbit beta 1 globin gene show three phenotypes determined by one fragment of variable length (about 5.1 or 3.3 kb). Family data demonstrate that these fragments segregate as Mendelian alleles. The frequencies of the two alleles are 0.66 for the 3.3-kb fragment and 0.34 for the 5.1-kb one. Another polymorphism has been detected with BamHI. Again three phenotypes determined by two alleles (fragments of 7.5 and 3.8 kb) have been observed. Allelic frequencies of the 7.5- and 3.8-kb fragments are 0.24 and 0.76 respectively. The two polymorphic sites are non-randomly associated.     

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.     

Direct detection of insertion/deletion polymorphisms in an autosomal region by analyzing high-density markers in individual spermatozoa

Direct polymerase chain reaction (PCR) detection of insertion/deletion (indel) polymorphisms requires sample homozygosity. For the indel polymorphisms that have the deletion allele with a relatively low frequency in the autosomal regions, direct PCR detection becomes difficult or impossible. The present study is, to our knowledge, the first designed to directly detect indel polymorphisms in a human autosomal region (i.e., the immunoglobulin V(H) region), through use of single haploid sperm cells as subjects. Unique marker sequences (n=32), spaced at approximately 5-kb intervals, were selected near the 3' end of the V(H) region. A two-round multiplex PCR protocol was used to amplify these sequences from single sperm samples from nine unrelated healthy donors. The parental haplotypes of the donors were determined by examining the presence or absence of these markers. Seven clustered markers in 6 of the 18 haplotypes were missing and likely represented a 35-40-kb indel polymorphism. The genotypes of the donors, with respect to this polymorphism, perfectly matched the expectation under Hardy-Weinberg equilibrium. Three V(H) gene segments, of which two are functional, are affected by this polymorphism. According to these results, >10% of individuals in the human population may not have these gene segments in their genome, and approximately 44% may have only one copy of these gene segments. The biological impact of this polymorphism would be very interesting to study. The approach used in the present study could be applied to understand the physical structure and diversity of all other autosomal regions.