||This article may require cleanup to meet Wikipedia's quality standards. (September 2009)|
In population genetics, linkage disequilibrium is the non-random association of alleles at two or more loci, that descend from single, ancestral chromosomes. Linkage disequilibrium is wholly a measurement of proximal genomic space. It is necessary to refer to this as gametic phase disequilibrium or simply gametic disequilibrium because it is described through DNA recombination. In other words, linkage disequilibrium is the occurrence of some combinations of alleles or genetic markers in a population more often or less often than would be expected from a random formation of haplotypes from alleles based on their frequencies. It is a second order phenomenon derived from linkage, which is the presence of two or more loci on a chromosome with limited recombination between them. The amount of linkage disequilibrium depends on the difference between observed allelic frequencies and those expected from a homogenous, randomly distributed model. Populations where combinations of alleles or genotypes can be found in the expected proportions are said to be in linkage equilibrium.
The level of linkage disequilibrium is influenced by a number of factors, including genetic linkage, selection, the rate of recombination, the rate of mutation, genetic drift, non-random mating, and population structure. A limiting example of the effect of rate of recombination may be seen in some organisms (such as bacteria) that reproduce asexually and hence exhibit no recombination to break down the linkage disequilibrium. An example of the effect of population structure is the phenomenon of Finnish disease heritage, which is attributed to a population bottleneck.
Consider the haplotypes for two loci A and B with two alleles each—a two-locus, two-allele model. Then the following table defines the frequencies of each combination:
Note that these are relative frequencies. One can use the above frequencies to determine the frequency of each of the alleles:
If the two loci and the alleles are independent from each other, then one can express the observation as " is found and is found". The table above lists the frequencies for , , and for, , hence the frequency of is , and according to the rules of elementary statistics .
In the genetic literature the phrase "two alleles are in LD" usually means that D ≠ 0. Contrariwise, "linkage equilibrium" means D = 0.
The following table illustrates the relationship between the haplotype frequencies and allele frequencies and D.
is easy to calculate with, but has the disadvantage of depending on the frequencies of the alleles. This is evident since frequencies are between 0 and 1. If any locus has an allele frequency 0 or 1 no disequilibrium can be observed. When the allelic frequencies are 0.5, the disequilibrium is maximal. Lewontin suggested normalising D by dividing it by the theoretical maximum for the observed allele frequencies.
Another measure of LD which is an alternative to is the correlation coefficient between pairs of loci, expressed as
This is also adjusted to the loci having different allele frequencies.
In summary, linkage disequilibrium reflects the difference between the expected haplotype frequencies under the assumption of independence, and observed haplotype frequencies. A value of 0 for indicates that the examined loci are in fact independent of one another, while a value of 1 demonstrates complete dependency.
Role of recombination
In the absence of evolutionary forces other than random mating, Mendelian segregation, random chromosomal assortment, and chromosomal crossover (i.e. in the absence of natural selection, inbreeding, and genetic drift), the linkage disequilibrium measure converges to zero along the time axis at a rate depending on the magnitude of the recombination rate between the two loci.
Using the notation above, , we can demonstrate this convergence to zero as follows. In the next generation, , the frequency of the haplotype , becomes
This follows because a fraction of the haplotypes in the offspring have not recombined, and are thus copies of a random haplotype in their parents. A fraction of those are . A fraction have recombined these two loci. If the parents result from random mating, the probability of the copy at locus having allele is and the probability of the copy at locus having allele is , and as these copies are initially in the two different gametes that formed the diploid genotype, these are independent events so that the probabilities can be multiplied.
This formula can be rewritten as
where at the -th generation is designated as . Thus we have
If , then so that converges to zero.
If at some time we observe linkage disequilibrium, it will disappear in the future due to recombination. However, the smaller the distance between the two loci, the smaller will be the rate of convergence of to zero.
Example: Human Leukocyte Antigen (HLA) alleles
HLA constitutes a group of cell surface antigens as MHC of humans. Because HLA genes are located at adjacent loci on the particular region of a chromosome and presumed to exhibit epistasis with each other or with other genes, a sizable fraction of alleles are in linkage disequilibrium.
|No. of individuals|
expression () frequency of antigen :
expression () frequency of antigen :
frequency of gene :
Denoting the '―' alleles at antigen i to be 'x,' and at antigen j to be 'y,' the observed frequency of haplotype xy is
and the estimated frequency of haplotype xy is
Then LD measure is expressed as
Standard errors are obtained as follows:
exceeds 2 in its absolute value, the magnitude of is large statistically significantly. For data in Table 1 it is 20.9, thus existence of statistically significant LD between A1 and B8 in the population is admitted.
|HLA-A alleles i||HLA-B alleles j|
Table 2 shows some of the combinations of HLA-A and B alleles where significant LD was observed among pan-europeans.
Vogel and Motulsky (1997) argued how long would it take that linkage disequilibrium between loci of HLA-A and B disappeared. Recombination between loci of HLA-A and B was considered to be of the order of magnitude 0.008. We will argue similarly to Vogel and Motulsky below. In case LD measure was observed to be 0.003 in Pan-europeans in the list of Mittal it is mostly non-significant. If had reduced from 0.07 to 0.003 under recombination effect as shown by , then . Suppose a generation took 25 years, this means 10,000 years. The time span seems rather short in the history of humans. Thus observed linkage disequilibrium between HLA-A and B loci might indicate some sort of interactive selection.
The presence of linkage disequilibrium between an HLA locus and a presumed major gene of disease susceptibility corresponds to any of the following phenomena:
- Relative risk for the person having a specific HLA allele to become suffered from a particular disease is greater than 1.
- The HLA antigen frequency among patients exceeds more than that among a healthy population. This is evaluated by value to exceed 0.
- 2x2 association table of patients and healthy controls with HLA alleles shows a significant deviation from the equilibrium state deduced from the marginal frequencies.
(1) Relative risk
Relative risk of an HLA allele for a disease is approximated by the odds ratio in the 2x2 association table of the allele with the disease. Table 3 shows association of HLA-B27 with ankylosing spondylitis among a Dutch population. Relative risk of this allele is approximated by
Woolf's method is applied to see if there is statistical significance. Let
follows the chi-square distribution with . In the data of Table 3, the significant association exists at the 0.1% level. Haldane's modification applies to the case when either of is zero, where replace and with
|Disease||HLA allele||Relative risk (%)||FAD (%)||FAP (%)|
|Spondylitis in inflammatory bowel disease||B27||10||50||8||0.46|
|Systemic lupus erythematosus||DR3||3||45||20||0.31|
|Diabetes mellitus type 1||DR4||6||75||30||0.64|
In Table 4, some examples of association between HLA alleles and diseases are presented.
(1a) Allele frequency excess among patients over controls
Even high relative risks between HLA alleles and the diseases were observed, only the magnitude of relative risk would not be able to determine the strength of association. value is expressed by
where and are HLA allele frequencies among patients and healthy populations, respectively. In Table 4, column was added in this quotation. Putting aside 2 diseases with high relative risks both of which are also with high values, among other diseases, juvenile diabetes mellitus (type 1) has a strong association with DR4 even with a low relative risk.
(2) Discrepancies from expected values from marginal frequencies in 2x2 association table of HLA alleles and disease
This can be confirmed by test calculating
A comparison of different measures of LD is provided by Devlin & Risch 
- PLINK - whole genome association analysis toolset, which can calculate LD among other things
- LdCompare— open-source software for calculating LD.
- SNP and Variation Suite- commercial software with interactive LD plot.
- GOLD - Graphical Overview of Linkage Disequilibrium
- TASSEL -software to evaluate linkage disequilibrium, traits associations, and evolutionary patterns
- Hardy-Weinberg principle
- Genetic linkage
- Genealogical DNA test
- Tag SNP
- Association Mapping
- Family based QTL mapping
- Reich, DE (2001). "Linkage disequilibrium in the human genome". Nature (Nautre Publishing Group) 411: 199–204. doi:10.1038/35075590. Retrieved 12 July 2013.
- Falconer, DS; Mackay, TFC (1996). Introduction to Quantitative Genetics (4th ed.). Harlow, Essex, UK: Addison Wesley Longman. ISBN 0-582-24302-5.
- Robbins, R.B. (1 July 1918). "Some applications of mathematics to breeding problems III". Genetics 3 (4): 375–389. PMC 1200443. PMID 17245911.
- R.C. Lewontin and K. Kojima (1960). "The evolutionary dynamics of complex polymorphisms". Evolution 14 (4): 458–472. doi:10.2307/2405995. ISSN 0014-3820. JSTOR 2405995.
- Lewontin, R. C. (1964). "The interaction of selection and linkage. I. General considerations; heterotic models". Genetics 49 (1): 49–67. PMC 1210557. PMID 17248194.
- Svejgaard A, Hauge M, Jersild C, Plaz P, Ryder LP, Staub Nielsen L, Thomsen M (1979). The HLA System: An Introductory Survey, 2nd ed. Basel; London; Chichester: Karger; Distributed by Wiley, ISBN 3805530498(pbk).
- Vogel F, Motulsky AG (1997). Human Genetics : Problems and Approaches, 3rd ed.Berlin; London: Springer, ISBN 3-540-60290-9.
- Mittal KK, Hasegawa T, Ting A, Mickey MR, Terasaki PI (1973). "Genetic variation in the HL-A system between Ainus, Japanese, and Caucasians," In Dausset J, Colombani J, eds. Histocompatibility Testing, 1972, pp. 187-195, Copenhagen: Munksgaard, ISBN 87-16-01101-5.
- Yasuda N, Tsuji K (1975). "A counting method of maximum likelihood for estimating haplotype frequency in the HL-A system." Jinrui Idengaku Zasshi 20(1): 1-15, PMID 1237691.
- Mittal KK (1976). "The HLA polymorphism and susceptibility to disease." Vox Sang 31: 161-173, PMID 969389.
- Gregersen PK (2009). "Genetics of rheumatic diseases," InFirestein GS, Budd RC, Harris ED Jr, McInnes IB, Ruddy S, Sergent JS, eds. (2009). Kelley's Textbook of Rheumatology, pp. 305-321, Philadelphia, PA: Saunders/Elsevier, ISBN 978-1-4160-3285-4.
- Bengtsson BO, Thomson G (1981). "Measuring the strength of associations between HLA antigens and diseases." Tissue Antigens18(5): 356-363, PMID 7344182.
- Nijenhuis LE (1977). "Genetic considerations on association between HLA and disease." Hum Genet38(2): 175-182, PMID 908564.
- Woolf B (1955). "On estimating the relation between blood group and disease." Ann Hum Genet 19(4): 251-253, PMID 14388528.
- Haldane JB (1956). "The estimation and significance of the logarithm of a ratio of frequencies." Ann Hum Genet20(4): 309-311, PMID 13314400.
- Sokal RR, Rohlf FJ (1981). Biometry: The Principles and Practice of Statistics in Biological Research. Oxford: W.H. Freeman, ISBN 0-7167-1254-7.
- Devlin B., Risch N. (1995). "A Comparison of Linkage Disequilibrium Measures for Fine-Scale Mapping". Genomics 29 (2): 311–322. doi:10.1006/geno.1995.9003. PMID 8666377.
- Hao K., Di X., Cawley S. (2007). "LdCompare: rapid computation of single- and multiple-marker r2 and genetic coverage". Bioinformatics 23 (2): 252–254. doi:10.1093/bioinformatics/btl574. PMID 17148510.
- Hedrick, Philip W. (2005). Genetics of Populations (3rd ed.). Sudbury, Boston, Toronto, London, Singapore: Jones and Bartlett Publishers. ISBN 0-7637-4772-6.
- Bibliography: Linkage Disequilibrium Analysis : a bibliography of more than one thousand articles on Linkage disequilibrium published since 1918.