Linkage Disequilibrium Diversity Research Articles

Estimating linkage disequilibria

Estimating linkage disequilibria

Pronounced inter- and intrachromosomal variation in linkage disequilibrium across the zebra finch genome

The extent of nonrandom association of alleles at two or more loci, termed linkage disequilibrium (LD), can reveal much about population demography, selection, and recombination rate, and is a key consideration when designing association mapping studies. Here, we describe a genome-wide analysis of LD in the zebra finch (Taeniopygia guttata) using 838 single nucleotide polymorphisms and present LD maps for all assembled chromosomes. We found that LD declined with physical distance approximately five times faster on the microchromosomes compared to macrochromosomes. The distribution of LD across individual macrochromosomes also varied in a distinct pattern. In the center of the macrochromosomes there were large blocks of markers, sometimes spanning tens of mega bases, in strong LD whereas on the ends of macrochromosomes LD declined more rapidly. Regions of high LD were not simply the result of suppressed recombination around the centromere and this pattern has not been observed previously in other taxa. We also found evidence that this pattern of LD has remained stable across many generations. The variability in LD between and within chromosomes has important implications for genome wide association studies in birds and for our understanding of the distribution of recombination events and the processes that govern them.

Aspects of observing and claiming allele flips in association studies

Significant allele flipping, where associations for the same disease occur at opposite alleles of the same bi-allelic locus, is increasing. But when is a significant allele flip genuine? We address the statistical issues of claiming and observing genuine allele flips in actual samples. We show that unless an allele flip is genuine, the probability of observing a significant allele flip in samples ascertained similarly from a common population is negligible. We derive expressions for the expected values of commonly used measures of association, which confirm previous findings that the underlying mechanism of a genuine allele flip is variation in the haplotype frequencies and show further how this variation interacts with variation in the genetic effects to impact allele flipping. We show that for association testing at proxy SNPs, common in genome-wide association studies, variation in haplotype frequencies must coincide with a reversal in the sign of linkage disequilibrium (LD) to trigger genuine allele flips. Using HapMap data and r, rather than r(2), to highlight previously unobserved effects, we show that unless genetic effects are large, variation in LD is unlikely to cause genuine allele flips in samples drawn from the same population. However, as populations diverge, it is an increasingly viable cause of a genuine allele flip for sufficiently large genetic effect and/or sample sizes. We conclude that evidence of variation in local patterns of LD, ancestral composition of study samples, and environmental exposures between study populations can provide compelling practical evidence in defense of a genuine allele flip.

SCAN: SNP and copy number annotation

Genome-wide association studies (GWAS) generate relationships between hundreds of thousands of single nucleotide polymorphisms (SNPs) and complex phenotypes. The contribution of the traditionally overlooked copy number variations (CNVs) to complex traits is also being actively studied. To facilitate the interpretation of the data and the designing of follow-up experimental validations, we have developed a database that enables the sensible prioritization of these variants by combining several approaches, involving not only publicly available physical and functional annotations but also multilocus linkage disequilibrium (LD) annotations as well as annotations of expression quantitative trait loci (eQTLs). For each SNP, the SCAN database provides: (i) summary information from eQTL mapping of HapMap SNPs to gene expression (evaluated by the Affymetrix exon array) in the full set of HapMap CEU (Caucasians from UT, USA) and YRI (Yoruba people from Ibadan, Nigeria) samples; (ii) LD information, in the case of a HapMap SNP, including what genes have variation in strong LD (pairwise or multilocus LD) with the variant and how well the SNP is covered by different high-throughput platforms; (iii) summary information available from public databases (e.g. physical and functional annotations); and (iv) summary information from other GWAS. For each gene, SCAN provides annotations on: (i) eQTLs for the gene (both local and distant SNPs) and (ii) the coverage of all variants in the HapMap at that gene on each high-throughput platform. For each genomic region, SCAN provides annotations on: (i) physical and functional annotations of all SNPs, genes and known CNVs within the region and (ii) all genes regulated by the eQTLs within the region. http://www.scandb.org. Supplementary data are available at Bioinformatics online.

Genome-wide comparisons of variation in linkage disequilibrium

Current genome-wide surveys of common diseases and complex traits fundamentally aim to detect indirect associations where the single nucleotide polymorphisms (SNPs) carrying the association signals are not biologically active but are in linkage disequilibrium (LD) with some unknown functional polymorphisms. Reproducing any novel discoveries from these genome-wide scans in independent studies is now a prerequisite for the putative findings to be accepted. Significant differences in patterns of LD between populations can affect the portability of phenotypic associations when the replication effort or meta-analyses are attempted in populations that are distinct from the original population in which the genome-wide study is performed. Here, we introduce a novel method for genome-wide analyses of LD variations between populations that allow the identification of candidate regions with different patterns of LD. The evidence of LD variation provided by the introduced method correlated with the degree of differences in the frequencies of the most common haplotype across the populations. Identified regions also resulted in greater variation in the success of replication attempts compared with random regions in the genome. A separate permutation strategy introduced for assessing LD variation in the absence of genome-wide data also correctly identified the expected variation in LD patterns in two well-established regions undergoing strong population-specific evolutionary pressure. Importantly, this method addresses whether a failure to reproduce a disease association in a disparate population is due to underlying differences in LD structure with an unknown functional polymorphism, which is vital in the current climate of replicating and fine-mapping established findings from genome-wide association studies.

Linkage disequilibrium in the North American Holstein population

Linkage disequilibrium was estimated using 7119 single nucleotide polymorphism markers across the genome and 200 animals from the North American Holstein cattle population. The analysis of maternally inherited haplotypes revealed strong linkage disequilibrium (r(2) > 0.8) in genomic regions of approximately 50 kb or less. While linkage disequilibrium decays as a function of genomic distance, genomic regions within genes showed greater linkage disequilibrium and greater variation in linkage disequilibrium compared with intergenic regions. Identification of haplotype blocks could characterize the most common haplotypes. Although maximum haplotype block size was over 1 Mb, mean block size was 26-113 kb by various definitions, which was larger than that observed in humans ( approximately 10 kb). Effective population size of the dairy cattle population was estimated from linkage disequilibrium between single nucleotide polymorphism marker pairs in various haplotype ranges. Rapid reduction of effective population size of dairy cattle was inferred from linkage disequilibrium in recent generations. This result implies a loss of genetic diversity because of the high rate of inbreeding and high selection intensity in dairy cattle. The pattern observed in this study indicated linkage disequilibrium in the current dairy cattle population could be exploited to refine mapping resolution. Changes in effective population size during past generations imply a necessity of plans to maintain polymorphism in the Holstein population.

ATOM: a powerful gene-based association test by combining optimally weighted markers

Large-scale candidate-gene and genome-wide association studies genotype multiple SNPs within or surrounding a gene, including both tag and functional SNPs. The immense amount of data generated in these studies poses new challenges to analysis. One particularly challenging yet important question is how to best use all genetic information to test whether a gene or a region is associated with the trait of interest. Here we propose a powerful gene-based Association Test by combining Optimally Weighted Markers (ATOM) within a genomic region. Due to variation in linkage disequilibrium, different markers often associate with the trait of interest at different levels. To appropriately apportion their contributions, we assign a weight to each marker that is proportional to the amount of information it captures about the trait locus. We analytically derive the optimal weights for both quantitative and binary traits, and describe a procedure for estimating the weights from a reference database such as the HapMap. Compared with existing approaches, our method has several distinct advantages, including (i) the ability to borrow information from an external database to increase power, (ii) the theoretical derivation of optimal marker weights and (iii) the scalability to simultaneous analysis of all SNPs in candidate genes and pathways. Through extensive simulations and analysis of the FTO gene in our ongoing genome-wide association study on childhood obesity, we demonstrate that ATOM increases the power to detect genetic association as compared with several commonly used multi-marker association tests.

Linkage disequilibrium, genetic association mapping and gene localization in crop plants

DNA-based molecular markers have been extensively utilized for a variety of studies in both plant and animal systems. One of the major uses of these markers is the construction of genome-wide molecular maps and the genetic analysis of simple and complex traits. However, these studies are generally based on linkage analysis in mapping populations, thus placing serious limitations in using molecular markers for genetic analysis in a variety of plant populations. Therefore, alternative approach has been suggested, linkage disequilibrium-based association analysis which detects and locates quantitative trait loci (QTL) by the strength of the correlation between a trait and a marker. Although association analysis has already been used for studies on genetics of complex traits in humans, its use in plants has newly started. In the present review, we describe what is known about variation in linkage disequilibrium (LD) and summarize published results on association studies in crop plant species. We give a list of different factors affecting LD, and discuss the current issues of LD research in plants. Later, we also describe the various uses of LD in crop plants research and summarize the present status of LD researches in different plant genomes. Finally, future key issues about the application of these studies on the localization of genes in these crop plants have been also discussed.

Linkage Disequilibrium Under Skewed Offspring Distribution Among Individuals in a Population

Correlations in coalescence times between two loci are derived under selectively neutral population models in which the offspring of an individual can number on the order of the population size. The correlations depend on the rates of recombination and random drift and are shown to be functions of the parameters controlling the size and frequency of these large reproduction events. Since a prediction of linkage disequilibrium can be written in terms of correlations in coalescence times, it follows that the prediction of linkage disequilibrium is a function not only of the rate of recombination but also of the reproduction parameters. Low linkage disequilibrium is predicted if the offspring of a single individual frequently replace almost the entire population. However, high linkage disequilibrium can be predicted if the offspring of a single individual replace an intermediate fraction of the population. In some cases the model reproduces the standard Wright-Fisher predictions. Contrary to common intuition, high linkage disequilibrium can be predicted despite frequent recombination, and low linkage disequilibrium under infrequent recombination. Simulations support the analytical results but show that the variance of linkage disequilibrium is very large.

Accommodating Linkage Disequilibrium in Genetic-Association Analyses via Ridge Regression

Large-scale genetic-association studies that take advantage of an extremely dense set of genetic markers have begun to produce very compelling statistical associations between multiple makers exhibiting strong linkage disequilibrium (LD) in a single genomic region and a phenotype of interest. However, the ultimate biological or "functional" significance of these multiple associations has been difficult to discern. In fact, the LD relationships between not only the markers found to be associated with the phenotype but also potential functionally or causally relevant genetic variations that reside near those markers have been exploited in such studies. Unfortunately, LD, especially strong LD, between variations at neighboring loci can make it difficult to distinguish the functionally relevant variations from nonfunctional variations. Although there are (rare) situations in which it is impossible to determine the independent phenotypic effects of variations in LD, there are strategies for accommodating LD between variations at different loci, and they can be used to tease out their independent effects on a phenotype. These strategies make it possible to differentiate potentially causative from noncausative variations. We describe one such approach involving ridge regression. We showcase the method by using both simulated and real data. Our results suggest that ridge regression and related techniques have the potential to distinguish causative from noncausative variations in association studies.

GABRB2 Association with Schizophrenia: Commonalities and Differences Between Ethnic Groups and Clinical Subtypes

Single nucleotide polymorphisms (SNPs) and haplotypes in intron 8 of type A gamma-aminobutyric acid (GABA(A)) receptor beta2 subunit gene (GABRB2) were initially found to be associated with schizophrenia in Chinese. This finding was subjected to cross-validation in this study with Japanese (JP) and German Caucasian (GE) subjects. Single nucleotide polymorphisms discovery and genotyping were carried out through resequencing of a 1839 base pair (bp) region in GABRB2. Tagging SNPs (tSNPs) were selected based on linkage disequilibrium (LD), combinations of which were analyzed with Bonferroni correction and permutation for disease association. Random resampling was applied to generate size- and gender-balanced cases and control subjects. Out of the 17 SNPs (9.2/kilobase [kb]) revealed, 6 were population-specific. Population variations in LD were observable, and at least two low LD points were identified in both populations. Although disease association at single SNP level was only shown in GE, strong association was demonstrated in both JP (p = .0002 - .0191) and GE (p = .0033 - .0410) subjects, centering on haplotypes containing rs1816072 and rs1816071. Among different clinical subtypes, the most significant association was exhibited by systematic schizophrenia. Cross-population validation of GABRB2 association with schizophrenia has been obtained with JP and GE subjects, with the genotype-disease correlations being strongest in systematic schizophrenia, the most severe subtype of the disease.

Evaluation of an algorithm of tagging SNPs selection by linkage disequilibrium

Background: Single nucleotide polymorphisms (SNPs) are the most abundant kind of genetic polymorphism in the human genome. They are important in both genetic research and genetic testing in a clinical setting, such as in the area of pharmacogenetics. In order to improve efficiency, tagging SNPs (tagSNPs) are selected in genes of interest to represent other co-related SNPs in linkage disequilibrium (LD) with the tagSNPs. Various algorithms have been proposed to identify a subset of single nucleotide polymorphisms as tagSNPs. Most algorithms of tagSNPs selection are haplotype-based, in which the spatial relationship between SNPs is considered. Currently, a more efficient cluster-based algorithm is proposed which clusters SNPs solely by a LD parameter, such as r 2. Here, we evaluated the sample distribution of r 2 and its effect on the cluster-based tagSNPs selection. Design and methods: The genotype data of 198 individual within a 500-kb region on 5q31 was used to evaluate the sample distribution of r 2 and its effect on the cluster-based tagSNPs selection. Results: It was found that the degree of variation of LD depends on the LD structure of genes. Conclusion: As a cluster-based tagSNPs selection algorithm does not take into account the spatial position of SNPs, a more stringent r 2 threshold is required to achieve more reliable tagSNPs selection.

Sequence features in regions of weak and strong linkage disequilibrium

We use genotype data generated by the International HapMap Project to dissect the relationship between sequence features and the degree of linkage disequilibrium in the genome. We show that variation in linkage disequilibrium is broadly similar across populations and examine sequence landscape in regions of strong and weak disequilibrium. Linkage disequilibrium is generally low within approximately 15 Mb of the telomeres of each chromosome and noticeably elevated in large, duplicated regions of the genome as well as within approximately 5 Mb of centromeres and other heterochromatic regions. At a broad scale (100-1000 kb resolution), our results show that regions of strong linkage disequilibrium are typically GC poor and have reduced polymorphism. In addition, these regions are enriched for LINE repeats, but have fewer SINE, DNA, and simple repeats than the rest of the genome. At a fine scale, we examine the sequence composition of "hotspots" for the rapid breakdown of linkage disequilibrium and show that they are enriched in SINEs, in simple repeats, and in sequences that are conserved between species. Regions of high and low linkage disequilibrium (the top and bottom quartiles of the genome) have a higher density of genes and coding bases than the rest of the genome. Closer examination of the data shows that whereas some types of genes (including genes involved in immune response and sensory perception) are typically located in regions of low linkage disequilibrium, other genes (including those involved in DNA and RNA metabolism, response to DNA damage, and the cell cycle) are preferentially located in regions of strong linkage disequilibrium. Our results provide a detailed analysis of the relationship between sequence features and linkage disequilibrium and suggest an evolutionary justification for the heterogeneity in linkage disequilibrium in the genome.

Platelet glycoprotein I(b)alpha and integrin alpha2 beta1 polymorphisms: gene frequencies and linkage disequilibrium in a population diversity panel.

Genetic variants in the GP1BA and ITGA2 genes have been proposed as potential modifiers for arterial vascular disease and bleeding disorders. Since ancestry may play an important role in the prevalence of these variants, we sought to determine their allele frequency and linkage disequilibrium in a collection of 1064 DNA samples from 51 ethnic groups. We studied haplotypes of ITGA2 defined by single nucleotide substitutions at positions -52, 807, and 1648, and GP1BA variants defined by sequence changes in positions -5 (Kozak), 1018 (T145M, HPA-2) and 1285 (VNTR A, B, C and D). Frequency of haplotypes of ITGA2 showed considerable variation across the different groups, with a higher prevalence of the haplotype -52C or T/807C/1648A observed in African compared with caucasian and Asian populations. The haplotypes 52C/807T/1648A and -52T/807T/1648A were not observed in caucasians or South Americans. While relative frequencies of the GP1BA Kozak alleles were comparable across groups, the methionine allele (HPA-2b) showed a higher frequency in Africa (0.26) than in the other groups. We also observed a high prevalence of the VNTR B allele in the African and Israeli populations. Haplotype analysis revealed incomplete linkage disequilibrium between the HPA-2 and VNTR alleles. Incorporation of GP1BA variants into the set of SNPs already genotyped by the HapMap project disrupted the pre-existing haplotype block. These data provide a valuable resource for optimal selection of variants best tailored for association studies of vascular disease or bleeding disorders when examining individuals of different ancestral origins.

Linkage disequilibrium patterns vary substantially among populations

A major initiative to create a global human haplotype map has recently been launched as a tool to improve the efficiency of disease gene mapping. The 'HapMap' project will study common variants in depth in four (and to a lesser degree in up to 12) populations to catalogue haplotypes that are expected to be common to all populations. A hope of the 'HapMap' project is that much of the genome occurs in regions of limited diversity such that only a few of the SNPs in each region will capture the diversity and be relevant around the world. In order to explore the implications of studying only a limited number of populations, we have analyzed linkage disequilibrium (LD) patterns of three 175-320 kb genomic regions in 16 diverse populations with an emphasis on African and European populations. Analyses of these three genomic regions provide empiric demonstration of marked differences in frequencies of the same few haplotypes, resulting in differences in the amount of LD and very different sets of haplotype frequencies. These results highlight the distinction between the statistical concept of LD and the biological reality of haplotypes and their frequencies. The significant quantitative and qualitative variation in LD among populations, even for populations within a geographic region, emphasizes the importance of studying diverse populations in the HapMap project to assure broad applicability of the results.

Serotonin gene polymorphisms and bipolar I disorder: Focus on the serotonin transporter

The pathogenesis of bipolar disorder may involve, at least in part, aberrations in serotonergic neurotransmission. Hence, serotonergic genes are attractive targets for association studies of bipolar disorder. We have reviewed the literature in this field. It is difficult to synthesize results as only one polymorphism per gene was typically investigated in relatively small samples. Nevertheless, suggestive associations are available for the 5HT2A receptor and the serotonin transporter genes. With the availability of extensive polymorphism data and high throughput genotyping techniques, comprehensive evaluation of these genes using adequately powered samples is warranted. We also report on our investigations of the serotonin transporter, SLC6A4 (17q11.1‐q12). An insertion/deletion polymorphism (5HTTLPR) in the promoter region of this gene has been investigated intensively. However, the results have been inconsistent. We reasoned that other polymorphism/s may contribute to the associations and the inconsistencies may be due to variations in linkage disequilibrium (LD) patterns between samples. Therefore, we conducted LD analyses, as well as association and linkage using 12 polymorphisms, including 5HTTLPR. We evaluated two samples. The first sample consisted of 135 US Caucasian nuclear families having a proband with bipolar I disorder (BDI, DSM IV criteria) and available parents. For case‐control analyses, the patients from these families were compared with cord blood samples from local Caucasian live births (n = 182). Our second, independent sample was recruited through the Systematic Treatment Enhancement Program for Bipolar Disorder (STEP‐BD, 545 cases, 548 controls). No significant associations were detected at the individual polymorphism or haplotype level using the case‐control or family‐based analyses. Our analyses do not support association between SLC6A4 and BDI families. Further studies using sub‐groups of BDI are worthwhile.

Linkage disequilibrium in young genetically isolated Dutch population

The design and feasibility of genetic studies of complex diseases are critically dependent on the extent and distribution of linkage disequilibrium (LD) across the genome and between different populations. We have examined genomewide and region-specific LD in a young genetically isolated population identified in the Netherlands by genotyping approximately 800 Short Tandem Repeat markers distributed genomewide across 58 individuals. Several regions were analyzed further using a denser marker map. The permutation-corrected measure of LD was used for analysis. A significant (P<0.0004) relation between LD and genetic distance on a genomewide scale was found. Distance explained 4% of the total LD variation. For fine-mapping data, distance accounted for a larger proportion of LD variation (up to 39%). A notable similarity in the genomewide distribution of LD was revealed between this population and other young genetically isolated populations from Micronesia and Costa Rica. Our study population and experiment was simulated in silico to confirm our knowledge of the history of the population. High agreement was observed between results of analysis of simulated and empirical data. We conclude that our population shows a high level of LD similar to that demonstrated previously in other young genetic isolates. In Europe, there may be a large number of young genetically isolated populations that are similar in history to ours. In these populations, a similar degree of LD is expected and thus they may be effectively used for linkage or LD mapping.

Genetic demography of Antioquia (Colombia) and the Central Valley of Costa Rica.

We report a comparative genetic characterization of two population isolates with parallel demographic histories: the Central Valley of Costa Rica (CVCR) and Antioquia (in northwest Colombia). The analysis of mtDNA, Y-chromosome and autosomal polymorphisms shows that Antioquia and the CVCR are genetically very similar, indicating that closely related parental populations founded these two isolates. In both populations, the male ancestry is predominantly European, whereas the female ancestry is mostly Amerind. In agreement with their isolation, the Amerindian mtDNA diversity of Antioquia and the CVCR is typical of ethnically-defined native populations and is markedly lower than in other Latin American populations. A comparison of linkage disequilibrium (LD) at 18 marker pairs in Antioquia and the CVCR shows that markers in LD in both populations are located at short genetic distances (<approximately 1 cM), whereas markers separated by greater distances are in LD only in the CVCR. This difference probably reflects stochastic variation of LD at the limited number of genome regions compared. The genetic similarity of the populations from Antioquia and the CVCR together with differences in LD between them should be exploitable for the identification and fine mapping of shared disease-related gene variants.

Population variation in linkage disequilibrium across the COMT gene considering promoter region and coding region variation.

Catechol-O-methyl transferase (COMT) catalyzes the first step in one of the major pathways in the degradation of catecholamines. The COMT gene on chromosome 22 has been considered a candidate gene for many neuropsychiatric disorders, in part because an exon 4 single nucleotide polymorphism (SNP) in COMT causes an amino acid substitution associated with significantly altered enzyme activity. This functional variant, detected as an NlaIII restriction site polymorphism (RSP), is polymorphic in populations from around the world. A four-site haplotype spanning 28 kb effectively encompasses COMT. This haplotype is comprised of two novel polymorphisms [a tetranucleotide short tandem repeat polymorphism (STRP) in intron 1 and a HindIII RSP at the 5' end of COMT], the NlaIII site, and another previously published site - a BglI RSP at the 3' end of the gene. Overall linkage disequilibrium (LD) for this haplotype is strong and significant in 32 population samples from around the world. Conditional probabilities indicate that, in spite of moderate to strong disequilibrium in most non-African populations, the NlaIII site, although often used for prediction, would not always be a reliable predictor of allelic variation at the other sites. Because other functional variation might exist, especially regulatory variation, these findings indicate that haplotypes would be more effective indicators of possible involvement of COMT in disease etiology.

Mating signal variation and bimodality in a mosaic hybrid zone between Chorthippus grasshopper species.

The grasshoppers Chorthippus brunneus and Chorthippus jacobsi are easily distinguished by male calling song and the number of stridulatory pegs on the hind femur, and form a mosaic hybrid zone in northern Spain. In this paper, we fit a two-dimensional cline to variation in male calling song characters, which are of particular interest as they are likely to be involved in mate choice by females. As with variation in peg number, local habitat makes only a small contribution in explaining deviations in mean song score from clinal expectations. However, the fitted width of the cline for song characters is significantly narrower than for peg number, suggesting that mating signals may be associated with reduced hybrid fitness in the field and that recombination rates are sufficient to allow clines for different characters to diverge in width. Despite this, estimates for the overall elevation in linkage disequilibrium at the zone center, based on covariance between peg and song characters, reveal a substantial overrepresentation of parental genotypes at the cline center relative to the expectations of a tension zone of similar width. Examination of covariance at individual sites reveals that this inflated estimate of linkage disequilibrium is caused by several sites where the distribution of phenotypes is effectively bimodal. This substantial variation in linkage disequilibrium at the cline center could result from local variation in the strength of assortative mating or selection against hybrids, or may reflect the long-distance colonization of empty habitat from outside the hybrid zone, which would continually create new contacts between parental genotypes at the cline center. Hybrid zones like this, in which strong linkage disequilibrium occurs in some situations but not in others, are of particular relevance to speciation research and allow investigation of the spread of combinations of alleles through different genetic and ecological backgrounds.