Linkage Disequilibrium in Natural and Experimental Populations of Drosophila Melanogaster

We have studied linkage disequilibrium in Drosophila melanogaster in two samples from a wild population and in four large laboratory populations derived from the wild samples. We have assayed four polymorphic enzyme loci, fairly closely linked in the third chromosome: Sod Est-6, Pgm, and Odh. The assay method used allows us to identify the allele associations separately in each of the two homologous chromosomes from each male sampled. We have detected significant linkage disequilibrium between two loci in 16.7% of the cases in the wild samples and in 27.8% of the cases in the experimental populations, considerably more than would be expected by chance alone. We have also found three-locus disequilibria in more instances than would be expected by chance. Some disequilibria present in the wild samples disappear in the experimental populations derived from them, but new ones appear over the generations. The effective population sizes required to generate the observed disequilibria by randomness range from 40 to more than 60,000 individuals in the natural population, depending on which locus pair is considered, and from 100 to more than 60,000 in the experimental populations. These population sizes are unrealistic; the fact that different locus-pairs yield disparate estimates within the same population argues against the likelihood that the disequilibria may have arisen as a consequence of population bottlenecks. Migration, or population mixing, cannot be excluded as the process generating the disequilibria in the wild samples, but can in the experimental populations. We conclude that linkage disequilibrium in these populations is most likely due to natural selection acting on the allozymes, or on loci very tightly linked to them.     

Simulation of the effect of population size on the evolution of the recombination fraction

A study, by means of computer simulation, has been performed on the evolution of recombination rate modifier genes in a system with three diallelic loci (A, B and C). The locus C, selectively neutral, is responsible for the modification of the recombination fraction between the major loci (A and B) which are subjected to selection. Two models have been analysed, the modifier allele being recessive in one of them, and codominant in the other, with infinite and finite populations. Distinct initial genic frequencies of the major loci and different selection coefficients have been utilised. We have found that the frequency of the allele which favours recombination increases in finite populations, and decreases slightly in infinite populations. These results are consistent with previous theory; presumably, selection favours alleles reducing recombination between epistatically interacting loci in a infinite population, since this reduces the breakup of advantageous combinations of alleles. However, in finite populations, selection favours the breakup of the random linkage disequilibria which are produced by random drift.     

Estimation of multilocus linkage disequilibria in diploid populations with dominant markers

Analysis of population structure and organization with DNA-based markers can provide important information regarding the history and evolution of a species. Linkage disequilibrium (LD) analysis based on allelic associations between different loci is emerging as a viable tool to unravel the genetic basis of population differentiation. In this article, we derive the EM algorithm to obtain the maximum-likelihood estimates of the linkage disequilibria between dominant markers, to study the patterns of genetic diversity for a diploid species. The algorithm was expanded to estimate and test linkage disequilibria of different orders among three dominant markers and can be technically extended to manipulate an arbitrary number of dominant markers. The feasibility of the proposed algorithm is validated by an example of population genetic studies of hickory trees, native to southeastern China, using dominant random amplified polymorphic DNA markers. Extensive simulation studies were performed to investigate the statistical properties of this algorithm. The precision of the estimates of linkage disequilibrium between dominant markers was compared with that between codominant markers. Results from simulation studies suggest that three-locus LD analysis displays increased power of LD detection relative to two-locus LD analysis. This algorithm is useful for studying the pattern and amount of genetic variation within and among populations.     

Linkage Disequilibrium in Natural Populations of DROSOPHILA MELANOGASTER

Two large, stable populations (Texas and Japan) of Drosophila melanogaster were surveyed at 21 allozyme loci on the second and third chromosomes and for chromosomal gene arrangements on those two chromosomes. Over 220 independent gametes were sampled from each population. The types and frequencies of the surveyed genetic variation are similar to those observed previously and suggest only slight differentiation among geographically distant populations. Linkage disequilibrium among linked allozymes loci is only slightly, if at all, detectable with these sample sizes. Linkage disequilibrium between linked inversions and allozymes loci is common especially when located in the same arm. These disequilibria appear to be in the same direction for most comparisons in the two population samples. This result is interpreted as evidence of similar selective environments (ecological and genetic) in the two populations. It is also noted that the direction of these linkage disequilibria appears to be oriented with respect to the gene frequencies at the component loci.     

Migration, selection and histocompatibility

The extreme polymorphism of the HLA system allows a very precise definition of populations by gene frequencies. The existence of several linked loci (HLA-A, B, C, D, DR) induces the possibility of detection of preferential chromosomes in populations (particular associations of alleles or linkage disequilibria). These linkage disequilibria indicate that the conditions required for Hardy-Weinberg equilibrium are not met. The linkage disequilibria can be used as indicators for the evaluation of genetic dift, founder effect and inbreeding in small populations or for the detection of selection and migration in larger populations.     

Genetic differentiation of neutral markers and quantitative traits in predominantly selfing metapopulations: confronting theory and experiments with Arabidopsis thaliana

The comparison of the genetic differentiation of quantitative traits (QST) and molecular markers (FST) can inform on the strength and spatial heterogeneity of selection in natural populations, provided that markers behave neutrally. However, selection may influence the behaviour of markers in selfing species with strong linkage disequilibria among loci, therefore invalidating this test of detection of selection. We address this issue by monitoring the genetic differentiation of five microsatellite loci (FST) and nine quantitative traits (QST) in experimental metapopulations of the predominantly selfing species Arabidopsis thaliana, that evolved during eight generations. Metapopulations differed with respect to population size and selection heterogeneity. In large populations, the genetic differentiation of neutral microsatellites was much larger under heterogeneous selection than under uniform selection. Using simulations, we show that this influence of selection heterogeneity on FST can be attributable to initial linkage disequilibria among loci, creating stronger genetic differentiation of QTL than expected under a simple additive model with no initial linkage. We found no significant differences between FST and QST regardless of selection heterogeneity, despite a demonstrated effect of selection on QST values. Additional data are required to validate the role of mating system and linkage disequilibria in the joint evolution of neutral and selected genetic differentiation, but our results suggest that FST/QST comparisons can be conservative tests to detect selection in selfing species.     

The strength of the selection barrier between populations

Genetic differences among populations exposed to selection form barriers against genetic exchange by mortality among hybrids. The strength of such a selection barrier, with which one (recipient) population reacts against immigration from another (donor) population, may be measured as the cumulative mean fitness of hybrids and their descendants relative to the fitness of the recipient population. Previous work analysed a case of weak selection with pairwise epistatic interactions by assuming small genetic distance between two populations in contact. The present study allows large genetic difference between the donor and recipient populations and considers weak multilocus selection with arbitrary epistatic interactions between two or more linked loci. An approximate analytical expression for the barrier strength is obtained as an expansion in which the strength of selection plays the role of a small parameter. It is shown that allele frequencies and gametic linkage disequilibria contribute in different ways to the strength of the selection barrier.     

Selection and genetic drift of polymorphisms within the merozoite surface protein-1 gene of Plasmodium falciparum

Intragenic recombination in the merozoite surface protein-1 gene (Msp-1) of Plasmodium falciparum is a major mechanism for allelic variation among natural parasite populations. The frequency of recombination depends on the intensity of transmission in the vector mosquito. In the present study, linkage disequilibrium between polymorphic 'loci' in the 5'- and 3'-regions of Msp-1 was examined in parasite populations from Brazilian Amazon and southern Vietnam and compared with that in a Thai population previously reported. The R2 test identified clusters of linkage disequilibria between the 5'- and 3'-regions, which are different among the three populations. However, the overall strength of linkage disequilibria was stronger in Brazil, a hypoendemic area, than in Vietnam and Thailand, mesoendemic areas, suggesting that linkage disequilibrium in Msp-1 inversely correlates with the intensity of transmission. To investigate possible mechanisms for linkage disequilibrium in Msp-1, we applied the Fst index, which measures the inter-population variance in allele frequency, to 'loci' in Msp-1 among the three populations. The Fst test identified two distinct regions with respect to inter-population allele frequency in Msp-1: one for highly divergent 'loci' in the 5'-region and the other for non-divergent 'loci' in the 3'-region. These results suggest that genetic drift is not the sole mechanism for linkage disequilibrium, but selection operates on 'loci' in the 3'-region in hypo- and mesoendemic areas of malaria.     

Distinguishing recombination and intragenic gene conversion by linkage disequilibrium patterns

Deterministic theory suggests that reciprocal recombination and intragenic, interallelic conversion have different effects on the linkage disequilibrium between a pair of genetic markers. Under a model of reciprocal recombination, the decay rate of linkage disequilibrium depends on the distance between the two markers, while under conversion the decay rate is independent of this distance, provided that conversion tracts are short. A population genetic three-locus model provides a function Q of two-locus linkage disequilibria. Viewed as a random variable, Q is the basis for a test of the relative impact of conversion and recombination. This test requires haplotype frequency data of a sufficiently variable three-locus system. One of the few examples currently available is data from the Human Leukocyte Antigen (HLA) class I genes of three Amerindian populations. We find that conversion may have played a dominant role in shaping haplotype patterns over short stretches of DNA, whereas reciprocal recombination may have played a greater role over longer stretches of DNA. However, in order to draw firm conclusions more independent data are necessary.     

Comparative Effects of Pollen and Seed Migration on the Cytonuclear Structure of Plant Populations. I. Maternal Cytoplasmic Inheritance

We explicitly solve and analyze a series of deterministic continent-island models to delimit the effects of pollen and seed migration on cytonuclear frequencies and disequilibria in random-mating, mixed-mating and self-fertilized populations. Given the critical assumption of maternal cytoplasmic inheritance, five major findings are (i) nonzero cytonuclear disequilibria will be maintained in the island population if and only if at least some migration occurs each generation through seeds with nonrandom cytonuclear associations; (ii) immigrant seeds with no cytonuclear disequilibria can strongly affect the genetic structure of the island population by generating significant and long-lasting transient associations; (iii) with all else being equal, substantially greater admixture disequilibria are generally found with higher rates of seed migration into, or higher levels of self-fertilization within, the island population (with the possible exception of the heterozygote disequilibrium); (iv) pollen migration can either enhance or reduce the cytonuclear disequilibria caused by seed migration, or that due to mixed-mating in the absence of seed migration, but the effect is usually small and appears primarily to make a noticeable difference in predominantly outcrossing populations; and (v) pollen migration alone cannot generate even transient disequilibria de novo in populations with completely random associations. This same basic behavior is exhibited as long as there is some random outcrossing in the island population. Self-fertilized populations represent a special case, however, in that they are necessarily closed to pollen migration, and nonzero disequilibria can be maintained even in the absence of seed migration. All of these general results hold whether the population is censused as adults or as seeds, but the ability to detect nonrandom cytonuclear associations can depend strongly on the life stage censused in populations with a significant level of random outcrossing. We suggest how these models might be used for the estimation of seed and pollen migration.     

Selection for recombination in structured populations

In finite populations, linkage disequilibria generated by the interaction of drift and directional selection (Hill-Robertson effect) can select for sex and recombination, even in the absence of epistasis. Previous models of this process predict very little advantage to recombination in large panmictic populations. In this article we demonstrate that substantial levels of linkage disequilibria can accumulate by drift in the presence of selection in populations of any size, provided that the population is subdivided. We quantify (i) the linkage disequilibrium produced by the interaction of drift and selection during the selective sweep of beneficial alleles at two loci in a subdivided population and (ii) the selection for recombination generated by these disequilibria. We show that, in a population subdivided into n demes of large size N, both the disequilibrium and the selection for recombination are equivalent to that expected in a single population of a size intermediate between the size of each deme (N) and the total size (nN), depending on the rate of migration among demes, m. We also show by simulations that, with small demes, the selection for recombination is stronger than both that expected in an unstructured population (m = 1 - 1/n) and that expected in a set of isolated demes (m = 0). Indeed, migration maintains polymorphisms that would otherwise be lost rapidly from small demes, while population structure maintains enough local stochasticity to generate linkage disequilibria. These effects are also strong enough to overcome the twofold cost of sex under strong selection when sex is initially rare. Overall, our results show that the stochastic theories of the evolution of sex apply to a much broader range of conditions than previously expected.     

Further analysis of complex allozyme polymorphisms in a barley population

New theory has recently been developed for two-locus models. In the light of this theory, an earlier analysis of esterase allozyme data from an experimental barley population has been modified to take proper account of initial gametic phase (linkage) disequilibria. The results show that the directions in which two-locus genotypic frequencies deviated from products of one-locus frequencies in this population followed those predicted by neutral descent theory. The observed departures were, however, much larger in size than predicted by the new descent measure theory, indicating that selection is operating in the population.     

Epistasis between deleterious mutations and the evolution of recombination

Epistasis and the evolution of recombination are closely intertwined: epistasis generates linkage disequilibria (i.e. statistical associations between alleles), whereas recombination breaks them up. The mutational deterministic hypothesis (MDH) states that high recombination rates are maintained because the breaking up of linkage disequilibria generated by negative epistasis enables more efficient purging of deleterious mutations. However, recent theoretical and experimental work challenges the MDH. Experimental evidence suggests that negative epistasis, required by the MDH, is relatively uncommon. On the theoretical side, population genetic models suggest that, compared with the combined effects of drift and selection, epistasis generates a negligible amount of linkage disequilibria. Here, we assess these criticisms and discuss to what extent they invalidate the MDH as an explanation for the evolution of recombination.     

Dynamics of correlated genetic systems. II. Simulation studies of chromosomal segments under selection

The dynamics of chromosomal segments under selection are investigated by comparing experimental data to simulations of simple models of selection. The simulations assume 93 loci distributed evenly along an entire chromosome. The two issues addressed in this paper concern rates of decay of linkage disequilibria for chromosomes under selection and rates of gene frequency change after perturbation of gametic frequencies to states near the edge of the gametic frequency simplex. The findings are: (1) If reasonable values of inbreeding depression are assumed, linkage disequilibria decays to zero but at a rate nearly twice that expected from neutral theory. Experimental results also show accelerated decay rates. The acceleration of decay seems to be a simple consequence of the increased heterozygosity produced by selection. It is, therefore, argued that massive linkage disequilibria, of the kind found by Franklin and Lewontin (1970) in their simulations, are unlikely to characterize the genetic structure of natural populations of random mating organisms. (2) It is possible to distinguish between two time-honored models of multilocus selection, known as the symmetric overdominant and classical models, on the basis of gene frequency change near the edge of the gametic frequency simplex assuming linkage disequilibria is intense. (3) Examination of experimental data from perturbation experiments shows that neither of these elementary models adequately account for observed rates of gene frequency change, although the symmetric overdominant model does provide the better fit. Instead the experimental data suggest a markedly nonuniform distribution of selective effects along the chromosome. The data also suggest that these selective effects combine in markedly nonadditive ways in determining joint fitness.     

RESPONSE TO SELECTION IN PARTIALLY SELF-FERTILIZING POPULATIONS. I. SELECTION ON A SINGLE TRAIT

Self-fertilization is a common form of reproduction in plants and it has important implications for quantitative trait evolution. Here, I present a model of selection on quantitative traits that can accommodate any level of self-fertilization. The “structured linear model” (SLM) predicts the evolution of the mean phenotype as a function of three distinct quantities: the mean additive genetic value, the directional dominance, and the mean inbreeding coefficient. Stochastic simulations of truncation selection demonstrate the accuracy of the SLM in predicting changes in the mean and variance of a quantitative trait over the full range of selfing rates. They also illustrate how complex interactions between selection and mating system determine the population distribution of inbreeding coefficients and also the amount of linkage disequilibrium. Changes in the genetic variance due to linkage disequilibria, which are commonly referred to as the “Bulmer effect,” are greatly magnified by selfing. This complicates the relationship between selfing rate and response to selection. Like the random mating theory, the parameters of the SLM can be estimated from phenotypic data.     

Gametic associations between inversion and allozyme polymorphisms in Drosophila buzzatii.

Gametic disequilibria between second chromosome polymorphic arrangements and seven linked allozyme loci were estimated in seven populations of Drosophila buzzatii from Argentina. Significant and consistent associations across populations were detected for Est-1, Est-2, Aldox, and XDH: Phenograms based on Nei's genetic distance showed that chromosomes carrying the 2ST arrangement were more similar to each other, irrespective of the population from which they were extracted, than to chromosomes carrying the derived 2J and 2JZ3.Restriction of recombination in heterokaryotypes seems to be the best explanation for the significant linkage disequilibria between inversions and the loci located inside the rearranged segments, for example, Est-1 and Aldox, or close to the break points, for example, Est-2. However, epistatic interactions between Xdh, which is outside the inversions and not near the break points, and loci tightly linked to the inversions, is the most likely explanation for the association between Xdh and chromosomal arrangements. Some of the associations detected in endemic Argentinean populations are coincident with data obtained in colonizing populations of the Old World and Australia. Thus historical processes that took place in the original area of the species' distribution can account for these linkage disequilibria in colonized populations of D. buzzatii.     

Population structures and the role of genetic exchange in the zoonotic pathogen Cryptosporidium parvum

Apicomplexan protozoan parasites include some of the most globally important human and animal pathogens, all of which have obligatory sexual cycles in their definitive hosts. Despite their importance and the relevance of understanding the population genetic structure and role of genetic exchange in generating diversity, population genetic analysis has largely been restricted to Plasmodium spp. and Toxoplasma gondii. These species show a considerable diversity of population structure suggesting different strategies for transmission and survival in mammalian hosts. We have undertaken a population genetic analysis of a further apicomplexan species (Cryptosporidium parvum) to extend our understanding of the diversity of genetic structures and test whether it has a clonal population structure. Nothing is known about the population structure of this parasite. We have analyzed 180 parasite isolates from both humans and cattle derived from a single discrete geographical area, using three minisatellite and four microsatellite markers that define 38 multilocus genotypes. Analysis of linkage disequilibria between pairs of loci combined with measures of genetic distance and similarity provides evidence that the sample comprises four genetically isolated populations. One group of human isolates consists primarily of two closely related multilocus genotypes (clonal), while the major subtypes of a second group, common to both humans and animals, show a panmictic population structure. The data provide an important step in understanding the role of genetic exchange in these parasites, which is an essential prerequisite for determining the value of multilocus genotyping for the analysis of sources of human infection as well as future molecular epidemiological studies.     

Sex ratio, reproductive mode and genetic diversity in Triops cancriformis

1. Aquatic invertebrates display a wide array of alternative reproductive modes from apomixis to hermaphroditism and cyclical parthenogenesis. These have important effects on genetic diversity and population structure. Populations of the 'living fossil' Triops cancriformis display a range of sex ratios, and various reproductive modes are thought to underlie this variation. Using sex ratio information and histological analyses European populations have been inferred to be gonochoric (with separate males and females), selfing hermaphroditic and androdioecious, a rare reproductive mode in which selfing hermaphrodites coexist with variable proportions of males. In addition, some populations have been described as meiotic parthenogens.
2. Here we use population genetic analysis using microsatellite loci in populations with a range of sex ratios including a gonochoric population, and marker segregation patterns in heterozygote individuals reared in isolation, to clarify the reproductive mode in this species.
3. Our data show that populations in general have very low levels of genetic diversity. Non-gonochoric populations show lower genetic diversity, more heterozygote deficiencies, higher inbreeding coefficients and stronger linkage disequilibria than the gonochoric population. The maintenance of some heterozygosity in populations is consistent with some male influence in T. cancriformis populations, as would be expected from an androdioecious reproductive system. Results of marker segregation in eggs produced in isolation from non-gonochoric populations indicate that meiosis occurs and are consistent with two reproductive modes: selfing hermaphroditism and a type of ameiotic parthenogenesis.
4. Overall, our data indicate that androdioecy and selfing hermaphroditism are the most likely reproductive modes of non-gonochoric European Triops populations. Triops populations are strongly structured, suggesting high genetic drift and low levels of gene flow.     

The meiotic consequences of side-arm bridge separation

Lewontin and Krakauer's test (Genetics 74: 175–195, 1973) of the spatial heterogeneity of F values has been applied to the data from four natural populations of Drosophila subobscura. It was found that genes behaving as selected were the ones that are found in linkage disequilibria with inversions within which they are located. What was really tested was then the behavior of these inversions. It seems desirable when this test is applied to have an assurance that the genes tested are not in linkage disequilibria with inversions.     

Small effective population sizes in a widespread selfing species, Lymnaea truncatula (Gastropoda: Pulmonata)

We present here a spatial and temporal population genetic survey of a common freshwater snail, also a predominantly selfing species, Lymnaea truncatula. The rate of genetic diversity loss was quantified by estimating the effective size (Ne) of the snail populations, using two different methods. A temporal survey allowed estimation of a variance effective size of the populations, and a spatial survey allowed the estimation of an inbreeding effective size, from two-locus identity disequilibria estimates. Both methods were consistent and provided low Ne values. Drift due to (i) high amounts of selfing and (ii) fluctuations in population sizes because of temporary habitats, and also selection coupled to genome-wide linkage disequilibria, could explain such reductions in Ne. The loss of genetic diversity appears to be counterbalanced only very partially by low apparent rates of gene flow.