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.     

DNA polymorphism haplotypes of the human apolipoprotein APOA1-APOC3-APOA4 gene cluster

Summary The genes coding for apolipoproteins A1, C3, and A4 (APOA1, APOC3, APOA4) are closely linked and tandemly organized within a 15-kilobase (kb) DNA segment on the long arm of human chromosome 11. The nucleotide variability of a 61-kb DNA segment containing these genes and their flanking sequences was studied by restriction analysis of a sample of 18 unrelated Northern Europeans using seven different genomic DNA probes. Eleven restriction site polymorphisms located within this DNA segment were used for haplotype analysis of 129 Mediterranean and 67 American black chromosomes. Estimation of the extent of nonrandom association between these polymorphisms indicated considerable linkage disequilibrium within the APOA1-APOC3-APOA4 gene cluster. Several haplotypes arose by recombination, and the rate of recombination within this gene cluster was estimated to be at least 4 times greater than that expected based on uniform recombination. The polymorphism information content of each of these polymorphisms, taken individually, ranges between 0.053 and 0.375, while that of their haplotypes ranges between 0.858 and 0.862. Therefore, DNA polymorphism haplotypes in the APOA1-APOC3-APOA4 gene cluster constitute a highly informative genetic marker on the long arm of human chromosome 11.     

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 fine mapping of quantitative trait loci: A simulation study

Recently, the use of linkage disequilibrium (LD) to locate genes which affect quantitative traits (QTL) has received an increasing interest, but the plausibility of fine mapping using linkage disequilibrium techniques for QTL has not been well studied. The main objectives of this work were to (1) measure the extent and pattern of LD between a putative QTL and nearby markers in finite populations and (2) investigate the usefulness of LD in fine mapping QTL in simulated populations using a dense map of multiallelic or biallelic marker loci. The test of association between a marker and QTL and the power of the test were calculated based on single-marker regression analysis. The results show the presence of substantial linkage disequilibrium with closely linked marker loci after 100 to 200 generations of random mating. Although the power to test the association with a frequent QTL of large effect was satisfactory, the power was low for the QTL with a small effect and/or low frequency. More powerful, multi-locus methods may be required to map low frequent QTL with small genetic effects, as well as combining both linkage and linkage disequilibrium information. The results also showed that multiallelic markers are more useful than biallelic markers to detect linkage disequilibrium and association at an equal distance.     

Assignment of 128 genes localized on human chromosome 14q to the IMpRH map

To provide a gene-based comparative map and to examine a porcine genome assembly using bacterial artificial chromosome-based sequence, we have attempted to assign 128 genes localized on human chromosome 14q (HSA14q) to a porcine 7000-rad radiation hybrid (IMpRH) map. This study, together with earlier studies, has demonstrated the following. (i) 126 genes were incorporated into two SSC7 RH linkage groups by C artha G ene analysis. (ii) In the remaining two genes, TOX4 linked to TCRA located in SSC7 by two-point analysis, whereas SIP1 showed no significant linkage with any gene/marker registered in the IMpRH Web Server. (iii) In the two groups, the gene clusters located from 19.9 to 36.5 Mb on HSA14q11.2-q13.3 and from 64.0 to 104.3 Mb on HSA14q23-q32.33 respectively were assigned to SSC7q21-q26. (iv) Comparison of the gene order between the present RH map and the latest porcine sequence assembly revealed some inconsistencies, and a redundant arrangement of 16 genes in the sequence assembly.     

Haplotypes within genes of β-chemokines in 17q11 are associated with multiple sclerosis: a second phase study

We previously defined haplotypes of single nucleotide polymorphisms (SNP) with possible relevance to multiple sclerosis (MS) in 2 CC chemokine ligand (CCL) clusters in chromosome 17q11. The 17q11 region was also identified as a susceptibility locus by a meta-analysis of linkage studies. To confirm and refine the previous finding in a second, high resolution SNP scan in a new set of families. We genotyped 232 SNPs in 1369 individuals in 361 MS families. Transmission of marker alleles and haplotypes from unaffected parents to affected offspring was tested by using the pedigree disequilibrium test, the TRANSMIT 2.5 program, and the family and haplotype based association tests. Distribution of linkage disequilibrium (LD) was assessed by ldmax. In consensus with observations in the first scan, the present study identified haplotypes within CCL3 and CCL15 in the telomeric CCL cluster. There was also an overlap in the findings in the centromeric CCL cluster. Strong and extensive LD was detected both within the centromeric and telomeric CCL gene clusters. The present study replicates our previous findings and further suggests the existence of MS associated haplotypes within genes of CCL3 and CCL15. Haplotypes of interest are also present within the centromeric gene cluster (including CCL2, CCL7, CCL11, CCL8, and CCL13), but extensive LD prevents further refinement of these haplotypes by using the methods applied. Sequencing of the identified chromosomal segments and their flanking regions will be necessary to define specific variants with direct relevance to MS pathogenesis.     

Detecting high-resolution polymorphisms in human coding loci by combining PCR and single-strand conformation polymorphism (SSCP) analysis.

A strategy is described that allows the development of polymorphic genetic markers to be characterized in individual genes. Segments of the 3' untranslated regions are amplified, and polymorphisms are detected by digestion with frequently cutting enzymes and with the detection of single-stranded conformation polymorphisms. This allows these genes, or DNA segments, to be placed on the linkage maps of human chromosomes. Polymorphisms in two genes have been identified using this approach. A HaeIII polymorphism was detected in the KIT proto-oncogene, physically assigned to chromosome 4q11-12. This polymorphism is linked to other chromosome 4p markers and is in linkage disequilibrium with a HindIII polymorphism previously described at this locus. We have also identified in the insulin-like growth factor1 receptor gene (IGF1R) a 2-bp deletion that is present at a frequency of .25 in the Caucasian population. Pedigree analysis with this insertion/deletion polymorphism placed the IGF1R gene at the end of the current linkage map of chromosome 15q, consistent with the physical assignment of 15q2526. Thus, polymorphisms in specific genes can be used to related the physical, genetic, and comparative maps of mammalian genomes and to simplify the testing of candidate genes for human diseases.     

The predisposition to type 1 diabetes linked to the human leukocyte antigen complex includes at least one non-class II gene

The human leukocyte antigen (HLA) complex, encompassing 3.5 Mb of DNA from the centromeric HLA-DPB2 locus to the telomeric HLA-F locus on chromosome 6p21, encodes a major part of the genetic predisposition to develop type 1 diabetes, designated "IDDM1." A primary role for allelic variation of the class II HLA-DRB1, HLA-DQA1, and HLA-DQB1 loci has been established. However, studies of animals and humans have indicated that other, unmapped, major histocompatibility complex (MHC)-linked genes are participating in IDDM1. The strong linkage disequilibrium between genes in this complex makes mapping a difficult task. In the present paper, we report on the approach we have devised to circumvent the confounding effects of disequilibrium between class II alleles and alleles at other MHC loci. We have scanned 12 Mb of the MHC and flanking chromosome regions with microsatellite polymorphisms and analyzed the transmission of these marker alleles to diabetic probands from parents who were homozygous for the alleles of the HLA-DRB1, HLA-DQA1, and HLA-DQB1 genes. Our analysis, using three independent family sets, suggests the presence of an additional type I diabetes gene (or genes). This approach is useful for the analysis of other loci linked to common diseases, to verify if a candidate polymorphism can explain all of the association of a region or if the association is due to two or more loci in linkage disequilibrium with each other.     

LDB2000: sequence-based integrated maps of the human genome.

MOTIVATION: Integrated maps are useful for gene mapping and establishing the relationship between recombination and sequence. In this paper we describe algorithms and their implementation for constructing sequence-based integrated maps of the human chromosomes, which are presented in LDB2000, a web based resource. Gene mapping efforts are now focussing on linkage disequilibrium mapping and extension of the integrated map to represent the extent of linkage disequilibrium in different genomic regions would further increase the utility of these maps. RESULTS: Sequence-based integrated maps have been completed for chromosomes 21 and 22. These maps provide locations for genes and polymorphic markers in sequence and on genetic linkage, radiation hybrid and cytogenetic scales. Single nucleotide polymorphisms associated with genes in the maps are also included and their sequence locations indicated. Related locus information, such as aliases and expression information, can be searched on the WWW site.     

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.     

Comparison of the power and accuracy of biallelic and microsatellite markers in population-based gene-mapping methods.

Because of their great abundance and amenability to fully automated genotyping, single-nucleotide polymorphisms (SNPs) and simple insertion/deletion are emerging as a new generation of markers for positional cloning. Although the efficiency and cost associated with the markers are important in the mapping of human disease genes, the power to detect the linkage between the marker and the disease locus, as well as the accuracy of the estimation of the map location of the disease gene, dictate the selection of the markers. Both the power and the accuracy depend not only on the type of the markers but also on other factors, such as the age of the disease mutation, the magnitude of the genetic effect, the marker-allele distribution in the population, mutation rates of marker loci, the frequency of the disease allele, the recombination fraction, and the methods for mapping the human disease genes. In this article, we develop a mathematical framework and the analytical formulas for calculation of the power and the accuracy and investigate the impact that the aforementioned factors have on the power and the accuracy, by using two population-based gene-mapping methods-likelihood-based linkage-disequilibrium mapping and the transmission/disequilibrium test, for both biallelic SNPs and microsatellites. These studies provide not only guidance in selection of the markers and in the design of the sample scheme for positional cloning but also insight into the biological bases of the mapping of human disease genes.     

Sequence variation in the CHAT locus shows no association with late-onset Alzheimer's disease

There is substantial evidence for a susceptibility gene for late-onset Alzheimer's disease (AD) on chromosome 10. One of the characteristic features of AD is the degeneration and dysfunction of the cholinergic system. The genes encoding choline acetyltransferase (ChAT) and its vesicular transporter (VAChT), CHAT and SLC18A3 respectively, map to the linked region of chromosome 10 and are therefore both positional and obvious functional candidate genes for late-onset AD. We have screened both genes for sequence variants and investigated each for association with late-onset AD in up to 500 late-onset AD cases and 500 control DNAs collected in the UK. We detected a total of 17 sequence variants. Of these, 14 were in CHAT, comprising three non-synonymous variants (D7N in the S exon, A120T in exon 5 and L243F in exon 8), one synonymous change (H547H), nine single-nucleotide polymorphisms in intronic, untranslated or promoter regions, and a variable number of tandem repeats in intron 7. Three non-coding SNPs were detected in SLC18A3. None demonstrated any reproducible association with late-onset AD in our samples. Levels of linkage disequilibrium were generally low across the CHAT locus but two of the coding variants, D7N and A120T, proved to be in complete linkage disequilibrium.     

A primary linkage map of the human chromosome 11q22-23 region

We have constructed a genetic map of the human chromosomal region 11q22-23 by multipoint linkage analysis of 13 DNA polymorphisms that we have condensed into eight loci. An analysis for linkage disequilibrium between tightly linked probe/enzyme systems allows us to make specific recommendations for future DNA typing at these loci. The resulting sex-averaged multipoint map spans approximately 80 cM and differs considerably from previously reported genetic maps of this region. Our mathematically derived "most likely order" of the markers is compatible with physical mapping data using somatic cell hybrids. The known localizations of at least 14 functional genes and several disease loci to 11q22-23, including ataxia telangiectasia, make the mapping of this region especially relevant to studies of disease pathogenesis.     

A first high-density map of 981 biallelic markers on human chromosome 14

As the largest set of sequence variants, single-nucleotide polymorphisms (SNPs) constitute powerful assets for mapping genes and mutations related to common diseases and for pharmacogenetic studies. A major goal in human genetics is to establish a high-density map of the genome containing several hundred thousand SNPs. Here we assayed 3.7 Mb (154,397 bp in 24 alleles) of chromosome 14 expressed sequence tags (ESTs) and sequence-tagged sites, for sequence variation in DNA samples from 12 African individuals. We identified and mapped 480 biallelic markers (459 SNPs and 21 small insertions and deletions), equally distributed between EST and non-EST classes. Extensive research in public databases also yielded 604 chromosome 14 SNPs (dbSNPs), 520 of which could be mapped and 19 of which are common between CNG (i.e., identified at the Centre National de Génotypage) and dbSNP polymorphisms. We present a dense map of SNP variation of human chromosome 14 based on 981 nonredundant biallelic markers present among 1345 radiation hybrid mapped sequence objects. Next, bioinformatic tools allowed 945 significant sequence alignments to chromosome 14 contigs, giving the precise chromosome sequence position for 70% of the mapped sequences and SNPs. In addition, these tools also permitted the identification and mapping of 273 SNPs in 159 known genes. The availability of this SNP map will permit a wide range of genetic studies on a complete chromosome. The recognition of 45 genes with multiple SNPs, by allowing the construction of haplotypes, should facilitate pharmacogenetic studies in the corresponding regions.     

Statistical estimation and pedigree analysis of CCR2-CCR5 haplotypes

As more SNP marker data becomes available, researchers have used haplotypes of markers, rather than individual polymorphisms, for association analysis of candidate genes. In order to perform haplotype analysis in a population-based case-control study, haplotypes must be determined by estimation in the absence of family information or laboratory methods for establishing phase. Here, we test the accuracy of the Expectation-Maximization (EM) algorithm for estimating haplotype state and frequency in the CCR2-CCR5 gene region by comparison with haplotype state and frequency determined by pedigree analysis. To do this, we have characterized haplotypes comprising alleles at seven biallelic loci in the CCR2-CCR5 chemokine receptor gene region, a span of 20 kb on chromosome 3p21. Three-generation CEPH families (n=40), totaling 489 individuals, were genotyped by the 5'nuclease assay (TaqMan). Haplotype states and frequencies were compared in 103 grandparents who were assumed to have mated at random. Both pedigree analysis and the EM algorithm yielded the same small number of haplotypes for which linkage disequilibrium was nearly maximal. The haplotype frequencies generated by the two methods were nearly identical. These results suggest that the EM algorithm estimation of haplotype states, frequency, and linkage disequilibrium analysis will be an effective strategy in the CCR2-CCR5 gene region. For genetic epidemiology studies, CCR2-CCR5 allele and haplotype frequencies were determined in African-American (n=30), Hispanic (n=24) and European-American (n=34) populations.     

An alpha satellite DNA polymorphism specific for the centromeric region of chromosome 13

Alpha satellite DNA is composed of variants of a short consensus sequence that are repeated in tandem arrays in the centromeric heterochromatin of each human chromosome. To define centromeric markers for linkage studies, we screened human genomic DNA for restriction fragment length polymorphisms using a probe detecting alphoid sequences on chromosomes 13 and 21. We describe one such DNA polymorphism. Analysis of linkage of this DNA marker to other polymorphic markers in the CEPH pedigrees demonstrates linkage to markers on the proximal long arm of chromosome 13 and defines the centromeric end of the linkage map of this chromosome.     

Increased gene coverage and Alu frequency in large linkage disequilibrium blocks of the human genome

The human genome has linkage disequilibrium (LD) blocks, within which single-nucleotide polymorphisms show strong association with each other. We examined data from the International HapMap Project to define LD blocks and to detect DNA sequence features inside of them. We used permutation tests to determine the empirical significance of the association of LD blocks with genes and Alu repeats. Very large LD blocks (>200 kb) have significantly higher gene coverage and Alu frequency than the outcome obtained from permutation-based simulation, whereas there was no significant positive correlation between gene density and block size. We also observed a reduced frequency of Alu repeats at the gaps between large LD blocks, indicating that their enrichment in large LD blocks does not introduce recombination hotspots that would cause these gaps.     

Association of the widespread A149P hereditary fructose intolerance mutation with newly identified sequence polymorphisms in the aldolase B gene.

Hereditary fructose intolerance (HFI) is a potentially fatal autosomal recessive disease resulting from the catalytic deficiency of fructose 1-phosphate aldolase (aldolase B) in fructose-metabolizing tissues. The A149P mutation in exon 5 of the aldolase B gene, located on chromosome 9q21.3-q22.2, is widespread and the most common HFI mutation, accounting for 57% of HFI chromosomes. The possible origin of this mutation was studied by linkage to polymorphisms within the aldolase B gene. DNA fragments of the aldolase B gene containing the polymorphic marker loci from HFI patients homozygous for the A149P allele were amplified by PCR. Absolute linkage to a common PvuII RFLP allele was observed in 10 A149P homozygotes. In a more informative study, highly heterozygous polymorphisms were detected by direct sequence determination of a PCR-amplified aldolase B gene fragment. Two two-allele, single-base-pair polymorphisms, themselves in absolute linkage disequilibrium, in intron 8 (C at nucleotide 84 and A at nucleotide 105, or T at 84 and G at 105) of the aldolase B gene were identified. Mendelian segregation of these polymorphisms was confirmed in three families. Allele-specific oligonucleotide (ASO) hybridizations with probes for both sequence polymorphisms showed that 47% of 32 unrelated individuals were heterozygous at these loci; the calculated PIC value was .37. Finally, ASO hybridizations of PCR-amplified DNA from 15 HFI patients homozygous for the A149P allele with probes for these sequence polymorphisms revealed absolute linkage disequilibrium between the A149P mutation and the 84T/105G allele. These results are consistent with a single origin of the A149P allele and subsequent spread by genetic drift.