Localization of a small genomic region associated with elevated ACE

Defining the relationship between multiple polymorphisms in a small genomic region and an underlying quantitative trait locus (QTL) represents a major challenge in human genetics. Pedigree analyses have shown that angiotensin I-converting enzyme (ACE) levels are influenced by a QTL located within or close to the ACE gene and most likely resides in the 3' region of this locus. We genotyped seven polymorphisms spanning 13 kb in the 3' end of ACE in 159 Afro-Caribbean subjects to evaluate the linkage disequilibrium between these sites and to narrow the genomic region associated with an elevated ACE level using a cladistic analysis. The linkage disequilibrium measurement D' and a haplotype tree revealed three distinct haplotype segments, presumably because of recombination. The value of the linkage disequilibrium parameter p(excess) was highest for site 22982, which is located in the middle segment. A series of nested, cladistic analyses confirmed that the other two regions are unlikely to be the ACE-linked QTL and that the variant resides in the middle region. Analyses of the same polymorphisms in 98 unrelated Europeans in the Monitoring Trends and Determinants in Cardiovascular Diseases (MONICA) study resulted in fewer haplotypes than were observed among the Afro-Caribbean subjects, suggesting that populations with greater genetic diversity may be especially informative for fine-scale mapping.     

A locus on chromosome 11p with multiple restriction site polymorphisms.

We have discovered and characterized a new polymorphic locus on chromosome 11p, D11S12, defined by an arbitrary genomic DNA segment cloned in the plasmid pADJ762. Four different polymorphic restriction sites with minor allele frequencies greater than 5% are revealed by Southern hybridization of this probe and its derivatives to digests of human DNAs. These include two MspI sites, a TaqI site, and a BclI site. The frequencies of the common haplotypes at this locus have been determined in a Utah population. Significant linkage disequilibrium has been demonstrated to exist between some pairs of polymorphic sites. A molecular map of this region has been determined, and the polymorphic sites have been localized. Comparison of physical separation with degree of linkage disequilibrium reveals an interesting case where an MspI site and a TaqI site that are separated by 6.8 kilobases (kb) show a greater degree of disequilibrium with each other than they do with two polymorphic sites located between them. One of the two interior sites is a BclI site that is approximately 0.2 kb away from the TaqI site but shows the same degree of disequilibrium with the TaqI site as with the MspI site 6.7 kb away. Although there is significant linkage disequilibrium at this locus, there are four major haplotypes with frequencies of 5% or greater, and the polymorphic information content (PIC) of this locus is .64.     

Protein variation in Adh and Adh-related in Drosophila pseudoobscura. Linkage disequilibrium between single nucleotide polymorphisms and protein alleles

A 3.5-kb segment of the alcohol dehydrogenase (Adh) region that includes the Adh and Adh-related genes was sequenced in 139 Drosophila pseudoobscura strains collected from 13 populations. The Adh gene encodes four protein alleles and rejects a neutral model of protein evolution with the McDonald-Kreitman test, although the number of segregating synonymous sites is too high to conclude that adaptive selection has operated. The Adh-related gene encodes 18 protein haplotypes and fails to reject an equilibrium neutral model. The populations fail to show significant geographic differentiation of the Adh-related haplotypes. Eight of 404 single nucleotide polymorphisms (SNPs) in the Adh region were in significant linkage disequilibrium with three ADHR protein alleles. Coalescent simulations with and without recombination were used to derive the expected levels of significant linkage disequilibrium between SNPs and 18 protein haplotypes. Maximum levels of linkage disequilibrium are expected for protein alleles at moderate frequencies. In coalescent models without recombination, linkage disequilibrium decays between SNPs and high frequency haplotypes because common alleles mutate to haplotypes that are rare or that reach moderate frequency. The implication of this study is that linkage disequilibrium mapping has the highest probability of success with disease-causing alleles at frequencies of 10%.     

Molecular Analysis of an Allozyme Cline: Alcohol Dehydrogenase in Drosophila Melanogaster on the East Coast of North America

Clines may either be selectively maintained or be the by-product of nonadaptive processes related to population structure and history. Drosophila melanogaster populations on the east coast of North America show a latitudinal cline in the frequencies of two common electrophoretically distinguishable alleles at the alcohol dehydrogenase locus (Adh), designated Adh-S and Adh-F. This cline may either be adaptive or an artifact of a possible recent dual founding of North American D. melanogaster populations in which frequencies of Adh alleles differed between founder populations. By means of a high resolution restriction-mapping technique, we studied the distribution of 113 haplotypes derived from 44 polymorphic DNA markers within the Adh region in 1533 individuals from 25 populations throughout the cline. We found significant clinal differentiation at the polymorphisms determining the mobility-difference causing amino acid replacement between Adh-F and Adh-S alleles. Hitchhiking was limited, despite extensive linkage disequilibrium, and other sites did not vary clinally. Such a pattern of differentiation implies that selection is responsible for the cline. To investigate whether selection acts only on the Adh-F/S site, we performed a ``selective equivalence'''' test under the assumption that all variability within the specified allelic class is selectively neutral. This revealed selective equivalence among Adh-S-bearing haplotypes, whose frequencies showed no differentiation throughout the cline, implying high levels of frequency-homogenizing gene flow. Geographical heterogeneity among Adh-F-bearing haplotypes implied the action of selection on one or more additional variants in linkage disequilibrium with Adh-F. In a further study of a subset of the data (n = 1076 from 18 populations), we found a combined insertion/deletion polymorphism, designated &1, located in the 5'' adult intron and in linkage disequilibrium with Adh-F, to show more marked clinal variation than Adh-F/S. Although the unequivocal identification of the precise target(s) of selection requires further study, we suggest that clinal selection may be acting epistatically on the Adh-F/S and &1 polymorphisms.     

Recombination disequilibrium in ideal and natural populations

Following Hardy-Weinberg disequilibrium (HWD) occurring at a single locus and linkage disequilibrium (LD) between two loci in generations, we here proposed the third genetic disequilibrium in a population: recombination disequilibrium (RD). RD is a measurement of crossover interference among multiple loci in a random mating population. In natural populations besides recombination interference, RD may also be due to selection, mutation, gene conversion, drift and/or migration. Therefore, similarly to LD, RD will also reflect the history of natural selection and mutation. In breeding populations, RD purely results from recombination interference and hence can be used to build or evaluate and correct a linkage map. Practical examples from F₂, testcross and human populations indeed demonstrate that RD is useful for measuring recombination interference between two short intervals and evaluating linkage maps. As with LD, RD will be important for studying genetic mapping, association of haplotypes with disease, plant breading and population history.     

A cladistic analysis of phenotypic associations with haplotypes inferred from restriction endonuclease mapping. I. Basic theory and an analysis of alcohol dehydrogenase activity in Drosophila

Because some genes have been cloned that have a known biochemical or physiological function, genetic variation can be measured in a population at loci that may directly influence a phenotype of interest. With this measured genotype approach, specific alleles or haplotypes in the probed DNA region can be assigned phenotypic effects. In this paper we address several problems encountered in implementing the measured genotype approach with restriction site data. A number of analytical problems arise in part as a consequence of the linkage disequilibrium that is commonly encountered when dealing with small DNA regions: 1) different restriction site polymorphisms are not statistically independent, 2) the sites being measured are not likely to be the direct cause of the associated phenotypic effects, 3) haplotype classes may be phenotypically heterogeneous, and 4) the sites that are most strongly associated with phenotypic effects are not necessarily the most closely linked to the actual genetic cause of the effects. When recombination and gene conversion are rare, the primary cause of linkage disequilibrium is history (mutational origin, genetic drift, hitchhiking, etc.). We deal with historical association directly by producing a cladogram that partially reconstructs the evolutionary history of the present-day haplotype variability. The cladogram defines a nested analysis of variance that simultaneously detects phenotypic effects, localizes the effects within the cladogram, and identifies haplotypes that are potentially heterogeneous in their phenotypic associations. The power of this approach is illustrated by an analysis of the associations between alcohol dehydrogenase (ADH) activity and restriction site variability in a 13-kb fragment surrounding the ADH locus in Drosophila melanogaster.     

Haplotypes at ATM identify coding-sequence variation and indicate a region of extensive linkage disequilibrium

Genetic variation in the human population may lead to functional variants of genes that contribute to risk for common chronic diseases such as cancer. In an effort to detect such possible predisposing variants, we constructed haplotypes for a candidate gene and tested their efficacy in association studies. We developed haplotypes consisting of 14 biallelic neutral-sequence variants that span 142 kb of the ATM locus. ATM is the gene responsible for the autosomal recessive disease ataxia-telangiectasia (AT). These ATM noncoding single-nucleotide polymorphisms (SNPs) were genotyped in nine CEPH families (89 individuals) and in 260 DNA samples from four different ethnic origins. Analysis of these data with an expectation-maximization algorithm revealed 22 haplotypes at this locus, with three major haplotypes having frequencies > or = .10. Tests for recombination and linkage disequilibrium (LD) show reduced recombination and extensive LD at the ATM locus, in all four ethnic groups studied. The most striking example was found in the study population of European ancestry, in which no evidence for recombination could be discerned. The potential of ATM haplotypes for detection of genetic variants through association studies was tested by analysis of 84 individuals carrying one of three ATM coding SNPs. Each coding SNP was detected by association with an ATM haplotype. We demonstrate that association studies with haplotypes for candidate genes have significant potential for the detection of genetic backgrounds that contribute to disease.     

Disequilibrium likelihoods for fine-scale mapping of a rare allele.

Genetic linkage studies based on pedigree data have limited resolution, because of the relatively small number of segregations. Disequilibrium mapping, which uses population associations to infer the location of a disease mutation, provides one possible strategy for narrowing the candidate region. The coalescent process provides a model for the ancestry of a sample of disease alleles, and recombination events between disease locus and marker may be placed on this ancestral phylogeny. These events define the recombinant classes, the sets of sampled disease copies descending from the meiosis at which a given recombination occurred. We show how Monte Carlo generation of the recombinant classes leads to a linkage likelihood for fine-scale mapping from disease haplotypes. We compare single-marker disequilibrium mapping with interval-disequilibrium mapping and discuss how the approach may be extended to multipoint-disequilibrium mapping. The method and its properties are illustrated with an example of simulated data, constructed to be typical of fine-scale mapping of a rare disease in the Japanese population. The method can take into account known features of population history, such as changing patterns of population growth.     

Linkage disequilibrium and recombination rate estimates in the self-incompatibility region of Arabidopsis lyrata

Genetic diversity is unusually high at loci in the S-locus region of the self-incompatible species of the flowering plant, Arabidopsis lyrata, not just in the S loci themselves, but also at two nearby loci. In a previous study of a single natural population from Iceland, we attributed this elevated polymorphism to linkage disequilibrium (LD) between variants at loci close to the S locus and the S alleles, which are maintained in the population by balancing selection. With the four S-flanking loci whose diversity we previously studied, we could not determine the extent of the region linked to the S loci in which neutral sites are affected. We also could not exclude the possibility of a population bottleneck, or of admixture, as causes of the LD. We have now studied four more distant loci flanking the S-locus region, and more populations, and we analyze the results using a theoretical model of the effect of balancing selection on diversity at linked neutral sites within and between different functional S-allelic classes. In the model, diversity is a function of the number of selectively maintained alleles and the recombination distances from the selectively maintained sites. We use the model to estimate the number of different functional S alleles, their turnover rate, and recombination rates between the S-locus region and other loci. Our estimates suggest that there is a small region of very low recombination surrounding the S-locus region.     

Geographic and haplotype structure of candidate type 2 diabetes susceptibility variants at the calpain-10 locus

Recently, a positional cloning study proposed that haplotypes at the calpain-10 locus (CAPN10) are associated with increased risk of type 2 diabetes, or non-insulin-dependent diabetes mellitus, in Mexican Americans, Finns, and Germans. To inform the interpretation of the original mapping results and to look for evidence for the action of natural selection on CAPN10, we undertook a population-based genotyping survey of the candidate susceptibility variants. First, we genotyped sites 43, 19, and 63 (the haplotype-defining variants previously proposed) and four closely linked SNPs, in 561 individuals from 11 populations from five continents, and we examined the linkage disequilibrium among them. We then examined the ancestral state of these sites by sequencing orthologous portions of CAPN10 in chimpanzee and orangutan (the identity of sites 43 and 19 was further investigated in a limited sample of other great apes and Old World and New World monkeys). Our survey suggests larger-than-expected differences in the distribution of CAPN10 susceptibility variants between African and non-African populations, with common, derived haplotypes in European and Asian samples (including one of two proposed risk haplotypes) being rare or absent in African samples. These results suggest a history of positive natural selection at the locus, resulting in significant geographic differences in polymorphism frequencies. The relationship of these differences to disease risk is discussed.     

Characterization of recombination in the HLA class II region.

Studies of linkage disequilibrium across the HLA class II region have been useful in predicting where recombination is most likely to occur. The strong associations between genes within the 85-kb region from DQB1 to DRB1 are consistent with low frequency of recombination in this segment of DNA. Conversely, a lack of association between alleles of TAP1 and TAP2 (approximately 15 kb) has been observed, suggesting that recombination occurs here with relatively high frequency. Much of the HLA class II region has now been sequenced, providing the tools to undertake detailed analysis of recombination. Twenty-seven families containing one or two recombinant chromosomes within the 500-kb interval between the DPB1 and DRB1 genes were used to determine patterns of recombination across this region. SSCP analysis and microsatellite typing yielded identification of 127 novel polymorphic markers distributed throughout the class II region, allowing refinement of the site of crossover in 30 class II recombinant chromosomes. The three regions where recombination was observed most frequently are as follows: the 45-kb interval between HLA-DNA and RING3 (11 cases), the 50-kb interval between DQB3 and DQB1 (6 cases), and an 8.8-kb segment of the TAP2 gene (3 cases). Six of the 10 remaining recombinants await further characterization, pending identification of additional informative markers, while four recombinants were localized to other intervals (outliers). Analysis of association between markers flanking HLA-DNA to RING3 (45 kb), as well as TAP1 to TAP2 (15 kb), by use of independent CEPH haplotypes indicated little or no linkage disequilibrium, supporting the familial recombination data. A notable sequence motif located within a region associated with increased rates of recombination consisted of a (TGGA)12 tandem repeat within the TAP2 gene.     

An excess of nonsynonymous polymorphism and extensive haplotype structure at the PtABI1B locus in European aspen (Populus tremula): a case of balancing selection in an obligately outcrossing plant?

Here, we describe an unusually pronounced haplotype structure at the PtABI1B locus in the obligately outcrossing tree Populus tremula. Both nucleotide diversity and divergence at PtABI1B was low compared to other P. tremula genes suggesting that the gene is located in a region with a low mutation rate. Despite this, PtABI1B shows a very marked excess of nonsynonymous polymorphisms across the entire coding region and linkage disequilibrium (LD) extending across the entire PtABI1B region of approximately 2.6 kb. Such extensive LD is normally not seen in P. tremula. The extensive LD at PtABI1B is caused by the presence of two distinct haplotypes. The haplotype structure is not caused by a lack of recombination in the region, because evidence of recombination can be detected. In addition, several statistical tests strongly reject neutrality for the PtABI1B region, suggesting that the unusual haplotype structure could be actively maintained by balancing selection.     

Selecting a maximally informative set of single-nucleotide polymorphisms for association analyses using linkage disequilibrium

Common genetic polymorphisms may explain a portion of the heritable risk for common diseases. Within candidate genes, the number of common polymorphisms is finite, but direct assay of all existing common polymorphism is inefficient, because genotypes at many of these sites are strongly correlated. Thus, it is not necessary to assay all common variants if the patterns of allelic association between common variants can be described. We have developed an algorithm to select the maximally informative set of common single-nucleotide polymorphisms (tagSNPs) to assay in candidate-gene association studies, such that all known common polymorphisms either are directly assayed or exceed a threshold level of association with a tagSNP. The algorithm is based on the r(2) linkage disequilibrium (LD) statistic, because r(2) is directly related to statistical power to detect disease associations with unassayed sites. We show that, at a relatively stringent r(2) threshold (r2>0.8), the LD-selected tagSNPs resolve >80% of all haplotypes across a set of 100 candidate genes, regardless of recombination, and tag specific haplotypes and clades of related haplotypes in nonrecombinant regions. Thus, if the patterns of common variation are described for a candidate gene, analysis of the tagSNP set can comprehensively interrogate for main effects from common functional variation. We demonstrate that, although common variation tends to be shared between populations, tagSNPs should be selected separately for populations with different ancestries.     

Genetic linkage and transmission disequilibrium of marker haplotypes at chromosome 1q41 in human systemic lupus erythematosus

Systemic lupus erythematosus (SLE) is a chronic autoimmune disease characterized by the production of autoantibodies to a wide range of self-antigens. Recent genome screens have implicated numerous chromosomal regions as potential SLE susceptibility loci. Among these, the 1q41 locus is of particular interest, because evidence for linkage has been found in several independent SLE family collections. Additionally, the 1q41 locus appears to be syntenic with a susceptibility interval identified in the NZM2410 mouse model for SLE. Here, we report the results of genotyping of 11 microsatellite markers within the 1q41 region in 210 SLE sibpair and 122 SLE trio families. These data confirm the modest evidence for linkage at 1q41 in our family collection (LOD = 1.21 at marker D1S2616). Evidence for significant linkage disequilibrium in this interval was also found. Multiple markers in the region exhibit transmission disequilibrium, with the peak single marker multiallelic linkage disequilibrium noted at D1S490 (pedigree disequilibrium test [PDT] global P value = 0.0091). Two- and three-marker haplotypes from the 1q41 region similarly showed strong transmission distortion in the collection of 332 SLE families. The finding of linkage together with significant transmission disequilibrium provides strong evidence for a susceptibility locus at 1q41 in human SLE.     

Direct Visualization of the Highly Polymorphic RNU2 Locus in Proximity to the BRCA1 Gene

Although the breast cancer susceptibility gene BRCA1 is one of the most extensively characterized genetic loci, much less is known about its upstream variable number tandem repeat element, the RNU2 locus. RNU2 encodes the U2 small nuclear RNA, an essential splicing element, but this locus is missing from the human genome assembly due to the inherent difficulty in the assembly of repetitive sequences. To fill the gap between RNU2 and BRCA1, we have reconstructed the physical map of this region by re-examining genomic clone sequences of public databases, which allowed us to precisely localize the RNU2 array 124 kb telomeric to BRCA1. We measured by performing FISH analyses on combed DNA for the first time the exact number of repeats carried by each of the two alleles in 41 individuals and found a range of 6-82 copies and a level of heterozygosity of 98%. The precise localisation of the RNU2 locus in the genome reference assembly and the implementation of a new technical tool to study it will make the detailed exploration of this locus possible. This recently neglected macrosatellite could be valuable for evaluating the potential role of structural variations in disease due to its location next to a major cancer susceptibility gene.     

Polymorphic DNA haplotypes at the human phenylalanine hydroxylase locus and their relationship with phenylketonuria

Summary Eight polymorphic restriction enzyme sites at the phenylalanine hydroxylase (PAH) locus were analyzed from the parental chromosomes in 33 Danish nuclear families with at least one phenylketonuric (PKU) child. Determination of haplotypes of 66 normal chromosomes and 66 chromosomes bearing mutant allele (S) demonstrated that there are at least two haplotypes which occur predominantly on PKU chromosomes and rarely otherwise. Overall, the relative frequencies of the various haplotypes are significantly different on PKU-and normal-allele bearing chromosomes, even though there is no predominantly occurring unique haplotype which can characterize the PKU chromosomes. In addition, no significant association (linkage disequilibrium) between any single polymorphic site and the mutant allele (s) was observed. The results suggest that either the phenylketonuric mutation was very ancient so that the polymorphic sites and the mutation have reached linkage equilibrium or the mutant allele (s) are the results of multiple mutations in the phenylalanine hydroxylase gene in man. Furthermore, a crude relationship between standardized linkage disequilibria and physical map distances of the polymorphic sites indicates that there is no apparent recombination hot-spot in the human phenylalanine hydroxylase gene, since the recombination rate within the locus apears to be uniform and likely to be occurring at a rate similar to that within the HLA gene cluster. The limitations of this later analysis are discussed in view of the sampling errors of disequilibrium measure used, and the potential untility of the PAH haplotypes for prenatal diagnosis and detection of PKU carriers is established.     

The BRCA1 Ashkenazi founder mutations occur on common haplotypes and are not highly correlated with anonymous single nucleotide polymorphisms likely to be used in genome-wide case-control association studies

Background

We studied linkage disequilibrium (LD) patterns at the BRCA1 locus, a susceptibility gene for breast and ovarian cancer, using a dense set of 114 single nucleotide polymorphisms in 5 population groups. We focused on Ashkenazi Jews in whom there are known founder mutations, to address the question of whether we would have been able to identify the 185delAG mutation in a case-control association study (should one have been done) using anonymous genetic markers. This mutation is present in approximately 1% of the general Ashkenazi population and 4% of Ashkenazi breast cancer cases. We evaluated LD using pairwise and haplotype-based methods, and assessed correlation of SNPs with the founder mutations using Pearson's correlation coefficient.

Results

BRCA1 is characterized by very high linkage disequilibrium in all populations spanning several hundred kilobases. Overall, haplotype blocks and pair-wise LD bins were highly correlated, with lower LD in African versus non-African populations. The 185delAG and 5382insC founder mutations occur on the two most common haplotypes among Ashkenazim. Because these mutations are rare, even though they are in strong LD with many other SNPs in the region as measured by D-prime, there were no strong associations when assessed by Pearson's correlation coefficient, r (maximum of 0.04 for the 185delAG).

Conclusion

Since the required sample size is related to the inverse of r, this suggests that it would have been difficult to map BRCA1 in an Ashkenazi case-unrelated control association study using anonymous markers that were linked to the founder mutations.     

Genetic polymorphisms at the human liver/islet glucose transporter (GLUT2) gene locus in Caucasian and West Indian subjects with type 2 (non-insulin-dependent) diabetes mellitus.

The liver/islet glucose transporter (GLUT2) is mainly expressed in the hepatocytes of the liver and the beta-cells of the pancreatic islets and a defect in this transporter could lead to diabetic phenotypes, such as relative hypoinsulinaemia and reduced uptake and metabolism of glucose in the liver. DNA from unrelated individuals was digested with the restriction endonucleases Bgl-I and Taq-I followed by blotting and hybridisation with a 32P-labelled GLUT2 cDNA which revealed three restriction fragment length polymorphisms (RFLPs) (B1, T1 and T2) in both Caucasian and West Indian populations. Linkage analysis between these variant sites demonstrated that the alleles of these polymorphisms were in strong linkage disequilibrium. Disease association of genetic variants at the GLUT2 locus with type 2 diabetes was examined with these RFLPs in both Caucasian (n = 54) and West Indian (n = 46) populations with type 2 diabetes. There were no significant differences in the frequency of alleles, genotypes or haplotypes between diabetic patients and non-diabetic controls. However, there were significant differences in the allele frequencies of all these three polymorphisms between Caucasian and West Indian populations.     

Molecular Genetic Variation in the Centromeric Region of the X Chromosome in Three Drosophila Ananassae Populations. I. Contrasts between the Vermilion and Forked Loci

We have surveyed three natural populations of Drosophila ananassae for restriction map variation at the forked (f) and vermilion (v) loci, using 6-cutter restriction enzymes. Both loci are located in the centromeric region of the X chromosome. Two major conclusions can be drawn from the data. First, we found strong evidence for population subdivision, i.e., significant differences in the frequency distributions of polymorphisms and/or haplotypes between the Burma, India, and Brazil populations. Secondly, the pattern of DNA sequence variation between the two loci is unexpectedly different. The level of nucleotide variation in the v locus region is reduced (relative to f), especially in the Burma population. Furthermore, in contrast to v, we found no insertions/deletions larger than 700 bp and no significant linkage disequilibrium at f. The genetic differentiation among subpopulations can readily be attributed to restricted migration as the predominant evolutionary force. According to population genetics theory, the differences in DNA polymorphisms between the two loci are in qualitative agreement with the hypothesis that recombination is reduced in the v locus region ("centromere effect") but not at f. In order to test this hypothesis directly, we determined the cytogenetic positions of several loci in the centromeric region by in situ hybridization and found by comparison with the genetic map that recombination at v is indeed very low, much lower than at f.