Note on the prior probability of autosomal linkage
J. Ott
DOI: https://doi.org/10.1111/j.1469-1809.1976.tb00148.x
1976-05-01
Annals of Human Genetics
Abstract:In a recent paper, Elston & Lange (1975) have confirmed a calculation by Renwick that the prior probability that two random loci should be on the same chromosome is equal to 0.054. In addition, they have extended these calculations to a number m of markers computing the prior probability that a trait T is within a distance d (call this ‘close’) to at least one of the m markers Mi, i = 1, . . . , m. This note has two purposes, (i) to give an intuitive explanation for a phenomenon unexplained by the authors and (ii) to show that by a simple device, the authors’ calculations can be very well approximated. The events C, = ‘ T close to M,’ and C, = ‘ T close to M,’ with i 4 j are not independent which is not easy to see intuitively and is not explained by the authors. For a single chromosome of length 1 , this dependency still persists except when d B 1 . It can easily be shown to be most pronounced at d = 112 and to decrease symmetrically around this point with increasing distance Id 1/21. For the case of a single chromosome, the intuitive explanation of this dependency goes as follows. T is more likely to be close to a marker M, when Mi is not located within a distance d of the ends of the chromosome. So, knowledge of the fact that T is close to Mi tells us that T is probably located away from the endpoints of the chromosome where it is also more likely to be close to another marker Mi. As a matter of fact, P(C, and Cj) > P(C,) P(C,) (the difference between these two expressions was taken above as a measure of dependency). Elston & Lange have taken these dependencies properly into account for all chromosomes which resulted in the quite complicated calculations presented. However, these dependencies may not be too important numerically. The assumption that the events ‘ T close to Mi’ for all markers Mi are independent permits an easy calculation of the prior probabilities of linkage for any number of markers. It is interesting to see that under this simplifying assumption, the results deviate very little from the correct ones obtained by Elston & Lange. I f p is the prior probability for T to be close to a single marker, then, under independence, the prior probability that T be close to at least one of the m markers is simply 1 ( 1 -p)”. Assuming, e.g. p = 1/46 corresponding to a recombination fraction of 40 %, 1 (45/46)20 = 0.36 as compared to the correct value of 0.37. Similarly, with p = 0.054 corresponding to a distance equalling the length of the largest human chromosome, 1 0 ~ 9 4 6 ~ ~ = 0.67 (correct value: 0.64). In practice, the most interesting range of distances is for recombination fractions smaller than 40%. For these, a consideration of the arguments in the second paragraph, extended to 22 chromosomes of unequal lengths, shows that the error becomes smaller with decreasing distance. For the case of two traits considered later in Elston & Lange’s work, we may again assume that closeness of one trait to the nearest marker is independent of whether or not the other trait is close to at least one marker. This independence assumption leads to results even closer to the correct ones than before. For a fixed distance d , let p,, be the authors’ prior probability for one trait and