Distribution of nonrandom associations between pairs of protein loci along the third chromosome of Drosophila melanogaster

Abstract

The within-chromosome distribution of gametic disequilibrium (GD) between protein loci, and the underlying evolutionary factors of this distribution, are still largely unknown. Here, we report a detailed study of GD between a large number of protein loci (15) spanning 87% of the total length of the third chromosome of Drosophila melanogaster in a large sample of haplotypes (600) drawn from a single natural population. We used a sign-based GD estimation method recently developed for multiallelic systems, which considerably increases both the statistical power and the accuracy of estimation of the intensity of GD. We found that strong GD between pairs of protein loci was widespread throughout the chromosome. In total, 22% of both the pairs of alleles and pairs of loci were in significant GD, with mean intensities (as measured by D' coefficients) of 0.43 and 0.31, respectively. In addition, strong GD often occurs between loci that are far apart. By way of illustration, 32% of the allele pairs in significant GD occurred within pairs of loci separated by effective frequencies of recombination (EFRs) of 15-20 cM, the mean D' value being 0.49. These observations are in sharp contrast with previous studies showing that GD between protein loci is rarely found in natural populations of outcrossing species, even between very closely linked loci. Interestingly, we found that most instances of significant interallelic GD (68%) involved functionally related protein loci. Specifically, GD was markedly more frequent between protein loci related by the functions of hormonal control, molybdenum control, antioxidant defense system, and reproduction than between loci without known functional relationship, which is indicative of epistatic selection. Furthermore, long-distance GD between functionally related loci (mean EFR 9 cM) suggests that epistatic interactions must be very strong along the chromosome. This evidence is hardly compatible with the neutral theory and has far-reaching implications for understanding the multilocus architecture of the functional genome. Our findings also suggest that GD may be a useful tool for discovering networks of functionally interacting proteins.