Estimates of Linkage Disequilibrium and the Recombination Parameter Determined from Segregating Nucleotide Sites in the Alcohol Dehydrogenase Region of Drosophila Pseudoobscura

The alcohol dehydrogenase (Adh) region of Drosophila pseudoobscura, which includes the two genes Adh and Adh-Dup, was used to examine the pattern and organization of linkage disequilibrium among pairs of segregating nucleotide sites. A collection of 99 strains from the geographic range of D. pseudoobscura were nucleotide-sequenced with polymerase chain reaction-mediated techniques. All pairs of the 359 polymorphic sites in the 3.5-kb Adh region were tested for significant linkage disequilibrium with Fisher's exact test. Of the 74,278 pairwise comparisons of segregating sites, 127 were in significant linkage disequilibrium at the 5% level. The distribution of five linkage disequilibrium estimators D(ij), D(2), r(ij), r(2) and Dij were compared to theoretical distributions. The observed distributions of D(ij), D(2), r(ij) and r(2) were consistent with the theoretical distribution given an infinite sites model. The observed distribution of Dij differed from the theoretical distribution because of an excess of values at -1 and 1. No spatial pattern was observed in the linkage disequilibrium pattern in the Adh region except for two clusters of sites nonrandomly associated in the adult intron and intron 2 of Adh. The magnitude of linkage disequilibrium decreases significantly as nucleotide distance increases, or a distance effect. Adh-Dup had a larger estimate of the recombination parameter, 4Nc, than Adh, where N is the effective population size and c is the recombination rate. A comparison of the mutation and recombination parameters shows that 7-17 recombination events occur for each mutation event. The heterogeneous estimates of the recombination parameter and the inverse relationship between linkage disequilibrium and nucleotide distance are no longer significant when the two clusters of Adh intron sites are excluded from analyses. The most likely explanation for the two clusters of linkage disequilibria is epistatic selection between sites in the cluster to maintain pre-mRNA secondary structure.     

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.     

Significant linkage disequilibrium between the Huntington disease gene and the loci D4S10 and D4S95 in the Dutch population.

Significant linkage disequilibrium has been found between the Huntington disease (HD) gene and DNA markers located around D4S95 and D4S98. The linkage-disequilibrium studies favor the proximal location of the HD gene, in contrast to the conflicting results of recombination analyses. We have analyzed 45 Dutch HD families with 19 DNA markers and have calculated the strength of linkage disequilibrium. Highly significant linkage disequilibrium has been detected with D4S95, consistent with the studies in other populations. In contrast with most other studies, however, the area of linkage disequilibrium extends from D4S10 proximally to D4S95, covering 1,100 kb. These results confirm that the HD gene most likely maps near D4S95.     

Linkage disequilibrium in natural populations of Trypanosoma cruzi (flagellate), the agent of Chagas' disease

We have studied linkage disequilibrium in natural populations of Trypanosoma cruzi, the agent of Chagas' disease, by analyzing (i) a set of 524 stocks from the whole geographical range of the parasite, characterized at four gene loci coding for enzymes; (ii) a subsample of 121 stocks characterized at 12 enzyme loci; and (iii) a subset of 386 stocks from six locations in Bolivia, characterized by four enzyme loci. Our results show that the linkage disequilibrium reaches the maximum possible value, given the observed allelic frequencies, for almost all the locus pairs. This result is most consistent with the hypothesis that genetic recombination is absent or very rare in T. cruzi natural populations. Partition of the linkage disequilibrium variance for the six Bolivian populations shows that both inter- and intrapopulation components are substantial and that the relationships among the components are D2IS less than D2ST, and D'2IS less than D'2ST. These inequalities are interpreted as the result of an interplay between genetic drift, rare or absent mating, and clonal selection in generating linkage disequilibrium in T. cruzi populations.     

Haplotypes and linkage disequilibrium at the phenylalanine hydroxylase locus, PAH, in a global representation of populations

Because defects in the phenylalanine hydroxylase gene (PAH) cause phenylketonuria (PKU), PAH was studied for normal polymorphisms and linkage disequilibrium soon after the gene was cloned. Studies in the 1980s concentrated on European populations in which PKU was common and showed that haplotype-frequency variation exists between some regions of the world. In European populations, linkage disequilibrium generally was found not to exist between RFLPs at opposite ends of the gene but was found to exist among the RFLPs clustered at each end. We have now undertaken the first global survey of normal variation and disequilibrium across the PAH gene. Four well-mapped single-nucleotide polymorphisms (SNPs) spanning approximately 75 kb, two near each end of the gene, were selected to allow linkage disequilibrium across most of the gene to be examined. These SNPs were studied as PCR-RFLP markers in samples of, on average, 50 individuals for each of 29 populations, including, for the first time, multiple populations from Africa and from the Americas. All four sites are polymorphic in all 29 populations. Although all but 5 of the 16 possible haplotypes reach frequencies >5% somewhere in the world, no haplotype was seen in all populations. Overall linkage disequilibrium is highly significant in all populations, but disequilibrium between the opposite ends is significant only in Native American populations and in one African population. This study demonstrates that the physical extent of linkage disequilibrium can differ substantially among populations from different regions of the world, because of both ancient genetic drift in the ancestor common to a large regional group of modern populations and recent genetic drift affecting individual populations.     

Linkage disequilibrium and recent selection at three immunity receptor loci in Drosophila simulans

Immune system genes in a California population sample of Drosophila simulans were shown to bear several hallmarks of the effects of past directional selection. One potential effect of directional selection is an increase in linkage disequilibrium among the polymorphic sites that are linked to the site under selection. In this study, we focus on three D. simulans immunity loci, Hmu, Sr-CI/Sr-CIII, and Tehao, for which the polymorphic sites are in nearly perfect linkage disequilibrium, an unusual finding even with respect to other immunity genes sampled from the same lines. The most likely explanation for this finding is that, at each locus, two divergent alleles have been selected to intermediate frequencies in the recent past. The extent to which the linkage disequilibrium extends to the flanks of each of the immunity genes is minimal, suggesting that the favored mutations actually occurred within the immunity genes themselves. Furthermore, the excess linkage disequilibrium found in the California population is not found in an African D. simulans population sample and may be a result of novel pathogen-mediated selection pressures encountered during establishment of non-African populations.     

Multiple linked β and α globin genes in Atlantic cod: A PCR based strategy of genomic exploration

Allozyme variation in Atlantic cod hemoglobins shows various signs of natural selection. We report a genomic exploration of globin genes in this non-model organism. Applying a PCR based strategy with a strict criterion of phylogenetically informative sites we estimate the number of linked β and α globin genes. We estimate PCR error rate by PCR of cloned DNA and recloning and by analysis of singleton variable sites among clones. Based on the error rate we exclude variable sites so that the remaining variation meets successively stricter criteria of doubleton and triplet variable site. Applying these criteria we find ten clusters of linked β/α globin genes in the genome of Atlantic cod. Six variable amino acid changes in both genes were found in linkage disequilibrium with silent nucleotide substitutions. A phylogenetic tree, based on our strictly phylogenetically informative sites among 57 clones from 19 individuals, is split into two major branches by an amino acid change in a β gene. This change is supported by extensive linkage disequilibrium between the amino acid change and numerous other phylogenetically informative silent nucleotide sites. The different gene sets in the genome may represent different loci encoding different globins and/or allelic variation at some loci.     

Molecular population genetics of the alcohol dehydrogenase locus in the Hawaiian drosophilid D. mimica

Sequence variation among 10 alleles of the alcohol dehydrogenase (Adh) gene of the Hawaiian drosophilid D. mimica was analyzed with reference to the evolutionary history of the Hawaiian subgroup as well as to levels and patterns of polymorphism of the Adh gene in continental drosophilid species. The Adh gene of D. mimica is less polymorphic than that of other drosophilid species, and no replacement substitutions were found. Statistical analyses of the Adh alleles suggested the action of balancing selection and revealed significant linkage disequilibrium among three of the variable sites. The effective population size was estimated to be only slightly smaller than that of continental species and, surprisingly, on the same order of magnitude as the actual size.      

Amino Acid Covariation in a Functionally Important Human Immunodeficiency Virus Type 1 Protein Region Is Associated With Population Subdivision

The frequently reported amino acid covariation of the highly polymorphic human immunodeficiency virus type 1 (HIV-1) exterior envelope glycoprotein V3 region has been assumed to reflect fitness epistasis between residues. However, nonrandom association of amino acids, or linkage disequilibrium, has many possible causes, including population subdivision. If the amino acids at a set of sequence sites differ in frequencies between subpopulations, then analysis of the whole population may reveal linkage disequilibrium even if it does not exist in any subpopulation. HIV-1 has a complex population structure, and the effects of this structure on linkage disequilibrium were investigated by estimating within- and among-subpopulation components of variance in linkage disequilibrium. The amino acid covariation previously reported is explained by differences in amino acid frequencies among virus subpopulations in different patients and by nonsystematic disequilibrium among patients. Disequilibrium within patients appears to be entirely due to differences in amino acid frequencies among sampling time points and among chemokine coreceptor usage phenotypes of virus particles, but not source tissues. Positive selection explains differences in allele frequencies among time points and phenotypes, indicating that these differences are adaptive rather than due to genetic drift. However, the absence of a correlation between linkage disequilibrium and phenotype suggests that fitness epistasis is an unlikely cause of disequilibrium. Indeed, when population structure is removed by analyzing sequences from a single time point and phenotype, no disequilibrium is detectable within patients. These results caution against interpreting amino acid covariation and coevolution as evidence for fitness epistasis.     

Search for a founder mutation in idiopathic focal dystonia from Northern Germany.

Both the discovery of the DYT1 gene on chromosome 9q34 in autosomal dominant early-onset torsion dystonia and the detection of linkage for one form of adult-onset focal dystonia to chromosome 18p (DYT7) in a family from northern Germany provide the opportunity to further investigate genetic factors in the focal dystonias. Additionally, reports of linkage disequilibrium between several chromosome 18 markers and focal dystonia, both in sporadic patients from northern Germany and in members of affected families from central Europe suggest the existence of a founder mutation underlying focal dystonia in this population. To evaluate the role of these loci in focal dystonia, we tested 85 patients from northern Germany who had primary focal dystonia, both for the GAG deletion in the DYT1 gene on chromosome 9q34 and for linkage disequilibrium at the chromosome 18p markers D18S1105, D18S1098, D18S481, and D18S54. None of these patients had the GAG deletion in the DYT1 gene. Furthermore, Hardy-Weinberg analysis of markers on 18p in our patient population and in 85 control subjects from the same region did not support linkage disequilibrium. Taken together, these results suggest that most cases of focal dystonia in patients of northern German or central European origin are due neither to the GAG deletion in DYT1 nor to a proposed founder mutation on chromosome 18p but must be caused by other genetic or environmental factors.     

A method for estimating the mutation,gene conversion and recombination parameters in small multigene families

A simple two-locus gene conversion model is considered to investigate the amounts of DNA variation and linkage disequilibrium in small multigene families. The exact solutions for the expectations and variances of the amounts of variation within and between two loci are obtained. It is shown that gene conversion increases the amount of variation within each locus and decreases the amount of variation between two loci. The expectation and variance of the amount of linkage disequilibrium are also obtained. Gene conversion generates positive linkage disequilibrium and the degree of linkage disequilibrium decreases as the recombination rate is increased. Using the theoretical results, a method for estimating the mutation, gene conversion, and recombination parameters is developed and applied to the data of the Amy multigene family in Drosophila melanogaster. The gene conversion rate is estimated to be approximately 60-165 times higher than the mutation rate for synonymous sites.     

Population structure and its effect on haplotype diversity and linkage disequilibrium surrounding the xa5 locus of rice (Oryza sativa L.)

To assess the usefulness of linkage disequilibrium mapping in an autogamous, domesticated species, we have characterized linkage disequilibrium in the candidate region for xa5, a recessive gene conferring race-specific resistance to bacterial blight in rice. This trait and locus have good mapping information, a tractable phenotype, and available sequence data, but no cloned gene. We sampled 13 short segments from the 70-kb candidate region in 114 accessions of Oryza sativa. Five additional segments were sequenced from the adjacent 45-kb region in resistant accessions to estimate the distance at which linkage disequilibrium decays. The data show significant linkage disequilibrium between sites 100 kb apart. The presence of the xa5 resistant reaction in two ecotypes and in accessions with different haplotypes in the candidate region may indicate multiple origins or genetic heterogeneity for resistance. In addition, genetic differentiation between ecotypes emphasizes the need for controlling for population structure in the design of linkage disequilibrium studies in rice.     

Developments of the Price equation and natural selection under uncertainty

Many approaches to the study of adaptation, following Darwin, centre on the number of offspring of individuals. Population genetics theory makes clear that predicting gene frequency changes requires more detailed knowledge, for example of linkage and linkage disequilibrium and mating systems. Because gene frequency changes underlie adaptation, this can lead to a suspicion that approaches ignoring these sophistications are approximate or tentative or wrong. Stochastic environments and sexual selection are two topics in which there are widespread views that focusing on number of offspring of individuals is not enough, and that proper treatments require the introduction of further details, namely variability in offspring number and linkage disequilibrium, respectively. However, the bulk of empirical research on adaptation and a great deal of theoretical work continue to employ these approaches. Here, a new theoretical development arising from the Price equation provides a formal justification in very general circumstances for focusing on the arithmetic average of the relative number of offspring of individuals.     

Refined mapping of the CSNU3 gene to a 1.8-Mb region on chromosome 19q13.1 using historical recombinants in Libyan Jewish cystinuria patients.

Cystinuria is a genetic disease manifested by the development of kidney stones. In some patients, the disease is caused by mutations in the SLC3A1 gene located on chromosome 2p. In others, the disease is caused by a gene that maps to chromosome 19q, but has not yet been cloned. Cystinuria is very common among Jews of Libyan ancestry living in Israel. Previously we have shown that the disease-causing gene in Libyan Jews maps to an 8-cM interval on chromosome 19q between the markers D19S409 and D19S208. Several markers from chromosome 19q showed strong linkage disequilibrium, and a specific haplotype was found in more than half of the carrier chromosomes. In this study we have analyzed Libyan Jewish cystinuria families with eight markers from within the interval containing the gene. Seven of these markers showed significant linkage disequilibrium. A common haplotype was found in 16 of the 17 carrier chromosomes. Analysis of historical recombinants placed the gene in a 1.8-Mb interval between the markers D19S430 and D19S874. Two segments of the historical carrier chromosome used to calculate the mutation's age revealed that the disease-causing mutation was introduced into this population 7-16 generations ago.     

Nucleotide Variation and Conservation at the Dpp Locus, a Gene Controlling Early Development in Drosophila

A study of polymorphism and species divergence of the dpp gene of Drosophila has been made. Eighteen lines from a population of D. melanogaster were sequenced for 5200 bp of the Hin region of the gene, coding for the dpp polypeptide. A comparison was made with sequence from D. simulans. Ninety-six silent polymorphisms and three amino acid replacement polymorphisms were found. The overall silent polymorphism (0.0247) is low, but haplotype diversity (0.0066 for effectively silent sites and 0.0054 for all sites) is in the range found for enzyme loci. Amino acid variation is absent in the N-terminal signal peptide, the C-terminal TGF-β peptide and in the N-terminal half of the pro-protein region. At the nucleotide level there is strong conservation in the middle half of the large intron and in the 3' untranslated sequence of the last exon. The 3' untranslated conservation, which is perfect for 110 bp among all the divergent species, is unexplained. There is strong positive linkage disequilibrium among polymorphic sites, with stretches of apparent gene conversion among originally divergent sequences. The population apparently is a migration mixture of divergent clades.     

Linkage disequilibrium mapping of the bolting gene in sea beet using AFLP markers

The possibility of using linkage disequilibrium mapping in natural plant populations was assessed. In studying linkage disequilibrium among 137 mapped AFLP markers in four populations of sea beet (Beta vulgaris ssp. maritima (L.) Arcang.) it was shown that tightly linked loci could be detected by screening for associations. It was hypothesized that the short distances spanned by linkage disequilibrium enable markers that are very tightly linked to a target gene to be identified. The hypothesis was tested by whole-genome screening of AFLP markers for association with the gene for the annual growth habit, the B gene, in a sample of 106 sea beets. Despite the dominant nature of AFLP, two markers showing significant linkage disequilibrium with the B gene were detected. The results indicate the potential use of linkage disequilibrium for gene mapping in natural plant populations.     

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.     

Vitiligo: complex segregation and linkage disequilibrium analyses with respect to microsatellite loci spanning the HLA

Familial clustering and linkage disequilibrium studies suggest that genetic factors predispose to vitiligo, although a clear transmission pattern and cosegregation of vitiligo with specific mutations have not been demonstrated. We collected pedigree data on vitiligo from a set of 56 multigeneration families belonging to the Paisa community from Antioquia, Colombia, with the goal of applying the unified model of complex segregation and linkage disequilibrium analyses to test the hypotheses of the existence of a major gene predisposing to vitiligo and that allelic or haplotype polymorphisms of microsatellite loci at 6p21.3-21.4 spanning HLA (D6S276, D6S265, D6S273, and D6S291) are associated with this predisposition. Minimum sibship sample size to discriminate dominant and recessive inheritance models was largely accomplished. Between the 15 models of complex segregation used, the one that best fitted the data was that of a major dominant gene and the existence of strong environmental effects acting on the recessive genotype. The penetrance and risk estimations discriminated two sets of vitiligo patients: those with early onset of vitiligo cosegregating with a dominant mode of inheritance without environmental effects, and those with late onset of vitiligo cosegregating with the recessive genotype and being influenced by environmental effects. After establishing the normal distribution of allelic frequencies and performing multiple comparisons correction, the linkage disequilibrium analysis suggested that a major genetic factor could be located at 6p21.3-21.4, because we detected significant case-control differences for allele 122 at D6S265 ( Pc=0.0264) and significant linkage disequilibrium between loci D6S276 and D6S273 in the cases but not in the controls. We cannot explain these results as a consequence of evolutionary forces or as genetic stratification acting differentially on cases and controls, because there was neither deviation from the Hardy-Weinberg expectations nor genetic subdivision between cases and controls, as theta; (non-biased F(ST)) was not significantly different from 0.     

Linkage disequilibrium under genetic hitchhiking in finite populations

The model of genetic hitchhiking predicts a reduction in sequence diversity at a neutral locus closely linked to a beneficial allele. In addition, it has been shown that the same process results in a specific pattern of correlations (linkage disequilibrium) between neutral polymorphisms along the chromosome at the time of fixation of the beneficial allele. During the hitchhiking event, linkage disequilibrium on either side of the beneficial allele is built up whereas it is destroyed across the selected site. We derive explicit formulas for the expectation of the covariance measure D and standardized linkage disequilibrium sigma 2D between a pair of polymorphic sites. For our analysis we use the approximation of a star-like genealogy at the selected site. The resulting expressions are approximately correct in the limit of large selection coefficients. Using simulations we show that the resulting pattern of linkage disequilibrium is quickly-i.e., in <0.1N generations-destroyed after the fixation of the beneficial allele for moderately distant neutral loci, where N is the diploid population size.     

Association of seven polymorphisms of the D2 dopamine receptor gene with brain receptor-binding characteristics

Association of alleles at the Taql A, Taql B, intron 6, Taql D, exon 7, exon 8, and promoter-141C sites of the D2 dopamine receptor gene with D2 dopamine receptor binding characteristics in the caudate nucleus of Caucasian alcoholic and nonalcoholic subjects was determined. For the Taql D, exon 7, exon 8, and promoter-141C sites there were no significant allelic differences in Bmax (number of binding sites) or Kd (binding affinity) of the D2 dopamine receptors. However, subjects having the minor alleles at the Taql A, Taql B, and intron 6 sites had significantly lower Bmax than subjects not having them. None of these three polymorphisms had any significant effect on Kd. Highly significant linkage disequilibria were observed among the Taql A, Taql B, and intron 6 polymorphic sites, but linkage disequilibria between these three sites and each of the Taql D, exon 7, exon 8, and promoter-141C sites were of lesser or of no significance. Taken together, these findings suggest that the Taql A, Taql B, and intron 6 polymorphisms, but not the Taql D, exon 7, exon 8, and promoter-141C polymorphisms, are in linkage disequilibrium with a functional allelic variant that affects D2 dopamine receptor expression.