Population genetic structure of Plasmodium falciparum in the two main African vectors, Anopheles gambiae and Anopheles funestus

We investigated patterns of genetic diversity of Plasmodium falciparum associated with its two main African vectors: Anopheles gambiae and Anopheles funestus. We dissected 10,296 wild-caught mosquitoes from three tropical sites, two in Cameroon (Simbock and Tibati, separated by 320 km) and one in Kenya (Rota, >2,000 km from the other two sites). We assayed seven microsatellite loci in 746 oocysts from 183 infected mosquito guts. Genetic polymorphism was very high in parasites isolated from both vector species. The expected heterozygosity (HE) was 0.79 in both species; the observed heterozygosities (HO) were 0.32 in A. funestus and 0.42 in A. gambiae, indicating considerable inbreeding within both vector species. Mean selfing (s) between genetically identical gametes was s = 0.33. Differences in the rate of inbreeding were statistically insignificant among sites and between the two vector species. As expected, because of the high rate of inbreeding, linkage disequilibrium was very high; it was significant for all 21 loci pairs in A. gambiae and for 15 of 21 pairs in A. funestus, although only two pairwise comparisons were between loci on the same chromosome. Overall, the genetic population structure of P. falciparum, as evaluated by F statistics, was predominantly clonal rather than panmictic, a population structure that facilitates the spread of antimalarial drug and vaccine resistance and thus may impair the effectiveness of malaria control efforts.

Differences in the biology and ecology of these two main vectors might entail a differential impact on the genetic composition of the P. falciparum populations they harbor and transmit.We compared patterns of genetic diversity of P. falciparum isolated from midguts of wild-collected A. gambiae and A. funestus.The parasite population genetic structure in its mosquito vectors is of primary relevance to understanding the disease's epidemiology and its evolution (including predicting and monitoring the spread of drug resistance alleles) and to devising effective means for its control.
Most population genetic investigations of P. falciparum have focused on the haploid blood stages of the parasite in its human host (12)(13)(14)(15)(16).In contrast, we investigated the parasite's genetic diversity with oocysts dissected from mosquito guts, which allows the observation of the products of meiosis and, thus, the direct estimation of the population structure of the parasite.A few earlier studies considered oocysts' genetic diversity but targeted a small number of antigen-coding loci subjected to strong host immune selection (17)(18), thus giving a possibly biased view of P. falciparum population structure.Using seven presumably neutral microsatellite loci, we have recently demonstrated that P. falciparum population structure is predominantly clonal in A. gambiae from Kenya (19).We now extend this approach to parasite populations from the two major African vector species, A. gambiae s.s. and A. funestus, collected in Central (Cameroon) and East Africa (Kenya) across successive years.
The number of oocysts per infected mosquito (Table 1) did not differ significantly between the two species: There were 6.8 oocysts per midgut (C.I. ϭ 5.6-8.6)for A. gambiae and 5.9 (C.I. ϭ 4.5-7.9)for A. funestus, either overall or within each site (bootstrap twosample t test, P ϭ 0.4).The frequency distribution of oocysts per mosquito was similar for the two species (stochastic equality test, P ϭ 0.2).The index of discrepancy was 0.99 for both species; the variance-to-mean ratios were 17.7 and 17.6 (Table 1).These values are characteristic of the highly aggregated oocyst distributions generally observed in anopheline mosquitoes, with a high proportion of the parasites concentrated in a small fraction of hosts (19)(20)(21).
All loci were highly polymorphic, showing an average (ϮSE) of 15.9 (Ϯ2.2) alleles per locus, ranging from 7 in Pfg377 to 23 in TA109.The mean allelic richness of parasite subpopulations within mosquitoes across sites and samples was identical in the two vectors (1.79 Ϯ 0.05 for A. gambiae and 1.79 Ϯ 0.04 for A. funestus) and did not differ significantly within each site between the two species (Wilcoxon two-sample test, P ϭ 0.7).
The genotypic frequencies suggested considerable inbreeding in the P. falciparum populations of either vector.The total observed and expected frequencies of heterozygotes were, respectively, H O ϭ 0.42 Ϯ 0.03 and H E ϭ 0.79 Ϯ 0.05 in A. gambiae and H O ϭ 0.32 Ϯ 0.03 and H E ϭ 0.79 Ϯ 0.04 in A. funestus.The two species did not differ in the mean expected heterozygosities of their parasite populations (SI Table 3) within each site or within the total sample (Wilcoxon two-sample tests, P Ͼ 0.3).
The discrepancy between observed and expected heterozygosities remained high and significant for each locus and within each mosquito species when we considered only mosquito guts with 10 or more oocysts per gut (H O ϭ 0. Population Structure.We used F statistics to explore the pattern and distribution of heterozygote deficits.For the whole population of oocysts distributed in their respective mosquito hosts (subpopulations), the total heterozygote deficit F IT , which measures the overall departure from panmixia, was F IT ϭ 0.53 (C.I. ϭ 0.49-0.56,P Ͻ 10 Ϫ4 ).For oocyst subpopulations from the two species separately, F IT ϭ 0.47 (C.I. ϭ 0.43-0.51,P Ͻ 10 Ϫ4 ) for A. gambiae and F IT ϭ 0.60 (C.I. ϭ 0.57-0.64,P Ͻ 10 Ϫ4 ) for A. funestus.This difference was due to the fraction of the total allelic variance of P. falciparum oocysts that was distributed among mosquitoes in the two vector species (see section titled F ST between mosquitoes).The total heterozygote deficit F IT can be partitioned into two components: F IS , which measures the deficit due to the deviation from panmixia within a mosquito gut, i.e., the deviation from random union of gametes within each oocyst relative to gametes from different oocysts within the same mosquito; and F ST , which measures the heterozygote deficit due to the partitioning of the total allelic variance among populations.This latter measure can be further split into a component that measures heterozygote deficit due to differences between subpopulations (oocysts within mosquito guts) within a sample (F ST between mosquitoes) and one that measures heterozygote deficit due to the three additional parameters that contribute to population structure: sampling date, site, and Anopheles species.FIS within mosquitoes.F IS measures departure from panmixia due to the nonrandom association of alleles within oocysts in mosquitoes (23).F IS was calculated for the 112 mosquitoes (subpopulations) with at least two oocysts.The average F IS was 0.18 (C.I. ϭ 0.14-0.23,P Ͻ 10 Ϫ4 ) for A. gambiae (68 subpopulations) and 0.23 (C.I. ϭ 0.17-0.28,P Ͻ 10 Ϫ4 ) for A. funestus (44 subpopulations).Single-locus F IS estimates were also positive and significant (Fig. 1A) and did not differ between species when all loci were combined (Wilcoxon two-sample test, P ϭ 0.3) or treated separately (see C.I. 95% of Fig. 1 A).

HIERFSTAT (see Materials and Methods).
We have also used FSTAT version 2.9.4 to assess differentiation among sites because the values obtained can be compared with published data on genetic differentiation of oocysts.The estimates of differentiation between sites ranged from 0.02 to 0.12, indicating weak but significant spatial differentiation in all cases (Table 2).
FST between sampling dates.Simbock was sampled in 2002, 2003, and 2004; Rota was sampled in 2002 and 2003.The allele frequency distributions were similar at different times (for all HIERFSTAT tests, P Ͼ 0.05).We therefore combined the samples from different times at these two sites.Genetic differentiation between Anopheles species.We used canonical correspondence analysis (24) to determine the relative contribution of the vector species to the global genetic structure of the P. falciparum populations.A graphic representation of the results shows largely overlapping distributions (centroids and ellipses of the 95% C.I. in Fig. 2) of the genetic variability of P. falciparum between vector species and among sites.Thus, the vector species did not contribute to the genetic population structure of the parasite.
A Monte Carlo permutation test on the first four canonical correspondence analysis axes that combines species and sites is statistically significant (P ϭ 0.001).This suggests some degree of geographic differentiation (see F ST between sites; Table 2), although not between species.This interpretation was corroborated by FSTAT version 2.9.4 analysis, performing 20 random samplings of two sets of oocyst subpopulations within each species at each site.The interspecies F ST values obtained were all included inside the range of intraspecies F ST values.
Linkage Disequilibrium.Linkage disequilibrium was significant for each of the 21 pairs of loci in A. gambiae, with R ϭ 0.60-0.69(P Ͻ 10 Ϫ3 with Bonferroni correction), and for 15 of 21 pairs in A. funestus, R ϭ 0.53-0.79(P Ͻ 10 Ϫ3 ).The difference between the two species is likely due to the smaller number of oocyst subpopulations (i.e., infected mosquitoes with more than one oocyst per gut) in A. funestus (44)      The seven microsatellite loci are distributed among five different chromosomes; therefore, only 2 of 21 pairwise comparisons involved loci on the same chromosome (19).

Discussion
Our results address three issues with important epidemiological consequences: the degree of inbreeding and the extent to which the population structure of P. falciparum is closer to clonal rather to panmictic; the degree of genetic differentiation among geographically separated populations of P. falciparum; and the distribution of the Plasmodium oocyst pool within its major mosquito vector species, A. gambiae and A. funestus, which may be indicative of a differential effectiveness in malaria transmission between the two species.
The ''clonal theory'' of parasitic protozoa (25) proposed that the population structure of many parasitic protozoa, and in particular malaria, appears to be a result of clonal propagation rather than of randomly mating parasites, although the relative importance of both reproductive mechanisms can vary within and among parasitic species in different regions and epidemiological settings.Plasmodium was considered a likely exception because this idea seems to be in conflict with the biology of the species: The transmission of the parasite from a human host to another necessarily requires a sexual stage in the mosquito vector and, thus, the formation of gametes, fertilization, and meiosis.However, although the sexual reproduction in the mosquito host is obligatory, it does not necessarily result in recombination between genetically distinct gametes.If the combining gametes both derive from the same clonal population of haploid ancestors, fertilization will occur between genetically identical male and female gametocytes, and thus physiological sexuality does not necessarily imply genetic sexuality.The issue remains of the extent to which fertilization will imply genetically identical or different gametes, in natural populations of P. falciparum sampled from different vector species, and/or in different epidemiological settings.
The extent of cross-fertilization and panmixia relative to clonal reproduction in P. falciparum has indeed been the subject of considerable debate during the last decade (13,19,(26)(27)(28)(29)(30).The current consensus is that this parasite species displays a range of population structures associated with the levels of transmission intensity (13).Regions of high P. falciparum transmission intensity correlate with high prevalence (31)(32) and high genetic diversity of the parasite (13,15,33), which result in a high proportion of (genetically) mixed infections in humans.Increased genetic diversity in the human host increases the probability that the mosquito will obtain genetically distinct gametocytes with its blood meal.This, in turn, increases the probability of parasite outcrossing within the mosquito gut (34)(35).Conversely, higher levels of inbreeding and clonal population structure are expected in regions with lower transmission rates (13,36).
We collected our mosquitoes in three different African regions of high-transmission intensity, with an annual entomologic inoculation rate ranging between 100 and 300 infective bites per person per year (2,37).Accordingly, we observed high genetic diversity of P. falciparum in each of the two vector species: H E ϭ 0.79 Ϯ 0.05 for A. gambiae and H E ϭ 0.79 Ϯ 0.04 for A. funestus.These values are similar to those obtained in a recent study using the same microsatellite loci to explore the level of genetic diversity in P. falciparum oocyst stages in Kenya (H E ϭ 0.76-0.88;ref. 19).Our results are also consistent with several investigations conducted on haploid blood stages in highly endemic African populations (12,13,15,38,39).
In agreement with the clonal theory and despite the high endemicity, we found strong and highly significant heterozygote deficits at the subpopulation (mosquito) level: In both A. gambiae and A. funestus, the observed heterozygosities were only approximately one-half of the expected values under the hypothesis of random union of the gametes.This led to estimates of the self-fertilization index (s) of 0.30 and 0.37, respectively, for A. gambiae and A. funestus, i.e., that approximately one gametocyte of every three or four is involved in selfing.These selfing values are similar to those previously reported for the oocyst populations of A. gambiae in Kenya: s ϭ 0.25 (F IS ϭ 0.15) (19).
Read et al. (40) estimated from observations of gametocyte sex ratios in the peripheral blood of patients from Madang in Papua New Guinea that 62% of all P. falciparum matings are between genetically identical gametes.Based on evidence from oocyst genotyping using antigen-coding loci, overall inbreeding coefficients were estimated to be 0.90 in Papua New Guinea (18) but 0.34 in Tanzania (17).This difference was attributed to the lower levels of malaria transmission in Papua New Guinea as compared with Tanzania (18).However, in an earlier study (19), we observed that selfing in P. falciparum is not a local strategy that might occur in some but not other localities or occasions, but rather that it is a generalized life history trait of P. falciparum, which occurs even in high-transmission regions.Moreover, the observation made earlier (19) that there is in P. falciparum a considerable proportion of additional inbreeding due to nonrandom distribution of oocyst genotypes among mosquito guts (average F ST ϭ 0.35 for 11 P. falciparum populations from Kenya infecting A. gambiae) was confirmed in the present study for the two vectors, A. gambiae and A. funestus (average F ST ϭ 0.35 and 0.50, respectively).
The evidence of considerable inbreeding derived from heterozygote deficits was confirmed by the high level of linkage disequilibrium between loci observed in all three localities and in either vector species.High linkage disequilibrium can be due to the high level of inbreeding arising both from selfing and from the consequent high relatedness among oocysts within each mosquito, so that there is little opportunity for recombination.Similar results have been obtained with parasitized blood samples (13,17,19,39,41) from different African regions and different years or seasons displaying a very wide variety of transmission intensities.The conclusion is that the population structure of P. falciparum is predominantly clonal and that this is the case for parasites obtained from both of the two major vectors of P. falciparum in Africa, A. gambiae and A. funestus.Thus, P. falciparum cannot be considered a panmictic species even in African regions with the highest endemicity.This has important epidemiological consequences: nonrandom mating and ''clonality'' in highly endemic regions can facilitate the spread of antimalarial drug resistance (42) and thus rapidly impair the effectiveness of measures for malaria control if monitoring lags behind.
How much genetic differentiation occurs among different regions or among vector species also has considerable epidemiological consequences, one of the most important being the understanding of the spread of the parasite resistance to drugs and vaccine (42).The degree of geographic differentiation between populations separated by Ϸ320 km (Simbock to Tibati) or Ͼ2,000 km (Cameroon to Kenya) is significant but quite low, F ST ϭ 0.02 to 0.12, by one measure (FSTAT version 2.9.4.) and not significant, F ST ϭ 0.004 to 0.02, by a different statistical method (HIERFSTAT; see Results).Bogreau et al. (39) have observed in blood samples a substantial degree of differentiation between urban and rural sites at the scale of the African continent (F ST ϭ 0.17-0.24).However, microsatellite investigations of blood stages, conducted in highly endemic regions, are consistent with those reported in the present paper, showing low levels of geographic differentiation even over considerable distances.Geographic differentiation is weak among four Sudan sites separated by 420 km (F ST ϭ 0.006-0.11;ref. 12).Similarly, the genetic differentiation among microsatellite loci from three different countries is very low (F ST ϭ 0.04; ref. 27).From human bloodstages, Anderson et al. (13) noted that the genetic variation of P. falciparum at the African scale is distributed mostly within localities and not at all or very little among them (F ST ϭ 0.003-0.012),despite distances Ͼ2,000 km.We found the same feature of the genetic variation distribution in Africa at the oocyst level by sampling the parasite within the mosquito definitive host.
Finally, our study addresses the role of the mosquito vector in shaping P. falciparum population structure.Approximately 10 Anopheles species are vectors of P. falciparum in Africa (3,4).The most important ones are A. gambiae and A. funestus because they are widely distributed throughout the continent, are highly anthropophilic and endophilic, and exhibit high survival rates (6, 7).In Kenya, the incidence of infected mosquitoes was not significantly different between the two vector species A. gambiae and A. funestus, but in the two Cameroon populations, the proportion of infected mosquitoes was significantly higher in A. funestus than in A. gambiae and as it was for all samples combined (3.1% vs. 1.7%,Table 1).
For large-scale malaria vector control, it is critical to understand whether ecological and other biological differences between the two most important vectors affect the population structure of the parasite with possible consequences on their efficiency in malaria transmission.One obvious major concern is the possibility for the P. falciparum parasite to evolve and adapt differentially to the different Anopheles species present in a malaria transmission area.If this is the case, implementation of a strategy aiming at neutralizing one vector species by introgression of a transgenic resistance trait may buffer the advantage created by the introduction of the transgene because it may liberate a niche for a rare parasite strain adapted to a different vector species to exploit, with major consequences on the efficiency of vector control strategies.Despite its potentially great significance, we know of no study, other than the present one, that has addressed this issue.Our results demonstrate a lack of evidence for any significant influence of malaria vector species diversity on the population structure of the transmitted pathogen, both A. gambiae and A. funestus being equally prone to be infected by any P. falciparum genotype.The population structure of P. falciparum is thus similar for the two species, with no segregation of particular P. falciparum genotypes within any of the two major African vector species.

Materials and Methods
Sampling Sites.Mosquitoes were collected in two sites from two ecoclimatic areas of Cameroon, Central Africa, Simbock (lat 3°51ЈN, long 11°30ЈE) and Tibati (lat 6°28ЈN, long 12°37ЈE), and from Rota (lat 0°06.3ЈS,long 34°40.3ЈE) in Western Kenya, East Africa.Simbock is a suburb of Yaounde ´, the capital of Cameroon, within the equatorial climate zone.Annual rainfall is Ϸ1,500 mm spread throughout the year.There are two rainy seasons, one from September to November and a shorter one from March to June.Tibati is in a rural setting Ϸ320 km north of Simbock, within the highland tropical climatic area of the Adamawa Plateau, with annual rainfall Ͼ1,500 mm and mean annual temperature of Ϸ22°C.Vegetation is of the Sudan-Guinean type.The single rainy season extends from March to October (43).Rota is in a rural setting, on the marshes of Lake Victoria in Western Kenya, with a tropical climate and annual rainfall Ͼ1,000 mm.In all three highly endemic malaria localities, transmission is perennial, with inoculation rates of Ϸ200 (Tibati), Ϸ300 (Simbock), and 100-300 (Rota) infective bites per human per year (2,3,37).In Cameroon, mosquitoes were captured between 2000 and 0600, indoors and outdoors, before feeding, from 10 houses Ϸ1 km or less apart in Simbock.In Tibati, mosquitoes were captured from seven houses, Ͻ5 km from each other.Specimens from Rota were aspirated from the walls of 10 huts in the morning and reared for 1 week in the laboratory.A. gambiae s.s. and A. funestus females were found in all three sites surveyed (along with other P. falciparum-transmitting Anopheles species; ref. 3).Once identified by means of morphological keys (6,7), the mosquitoes were dissected in a saline solution on the morning of their collection in Cameroon and 1 week after collection in Rota to allow for the development of oocysts in infected females.Dissected midguts were checked for the presence of oocysts under a compound microscope.Oocystinfected guts were stored individually, along with the remainder of the mosquito body, in 1.5-ml microtubes with 80% ethanol.Samples were transported to the Montpellier laboratory and stored at 4°C.Infected females belonging to the A. gambiae complex and A. funestus group were identified to species by PCR diagnostic assays (44,45).
Oocyst Dissection.Oocysts were dissected from each gut after gradual rehydration of the specimens in distilled water by using a Leica (Deerfield, IL) inverted DMIRB microscope.Individual oocysts were teased from the basal lamina of the midgut epithelium and transferred to individual 0.5-ml tubes with 10 l of distilled water.Tips and dissecting tools were washed in Javel water and 70% ethanol and changed regularly to avoid contamination.Parasite DNA from each oocyst was extracted by using the DNeasy Tissue Kit of Qiagen (Valencia, CA).
Vector Infection Rate and Parasite Distribution.Rates of mosquito infection, number of oocysts per mosquito, mean intensity, indices of parasite distribution within mosquitoes (variance-to-mean ratio and index of discrepancy; ref. 46), and comparisons between Anopheles species were computed by using the QUANTITATIVE PARASITOLOGY 3.0 computer program (47).
Genetic Analysis.Genetic variability was quantified as allelic richness by using FSTAT version 2.9.4 (48,49), as well as observed (H O ) and expected heterozygosity under panmixia (H E ), by using the unbiased estimator of Nei (50), correcting for small samples.Comparisons between populations were performed by using nonparametric tests on two or more samples (Wilcoxon two-sample and Kruskal-Wallis tests) by using the R software (51,52).
We used FSTAT version 2.9.4 (49) to calculate F statistics.For a two-level population structure (individual and subpopulation levels), the total F IT , measuring the inbreeding coefficient of an individual (I) relative to the total (T), can be partitioned into F IS , inbreeding coefficient of an individual (I) relative to the subpopulation (S), which measures departure from Hardy-Weinberg expectations under the hypothesis of panmixia within each subpopulation; and F ST , which measures genetic differentiation between subpopulations (S) in the total population (T): The data set was organized into three levels: (i) oocysts within the mosquito (subpopulation or the mosquito level), (ii) oocysts infecting mosquitoes within a site (population or site level), and (iii) all oocysts infecting mosquitoes in the whole geographic area sampled (metapopulation or geographic area level).
F IS was estimated and tested by randomizing alleles between oocysts within each of the 112 oocyst subpopulations (112 infected mosquitoes with more than one oocyst).F ST between subpopulations (the fraction of the total allelic variance due to clustering of mosquitoes within geographical sites) was estimated and tested by randomizing oocyst genotypes between subpopulations (mosquitoes) within each population (site).Because FSTAT version 2.9.4 is designed for a two-level model and cannot take into account the three-level structure of our data set, differentiation between sites was estimated in two ways.First, the fixation index F ST was tested by using FSTAT version 2.9.4 on all oocysts sampled in each locality by randomizing oocyst genotypes between populations relative to total.Second, because the highly significant differentiation observed between subpopulations strongly affects population differentiation computed between oocyst populations from different vectors, sites, or years (53), the fixation index was also estimated by using HIERFSTAT version 0.04-2 (54), which computes F statistics from any number of hierarchical levels, enabling us to take into account the highly significant effect of the mosquito level on parasite population structure to estimate spatial differentiation between populations.Significance of the hierarchical F statistics was assessed by 1,000 permutations of subpopulations (mosquitoes) among populations (sites).The crossed factor ''date of sampling'' was also estimated, and its significance was assessed by using HIERFSTAT version 0.04-2.
The Anopheles species cannot be considered as a hierarchical level because mosquitoes from different species can feed on the same human host.To measure and illustrate levels of genetic differentiation between P. falciparum harbored by the two vectors in the three sites, a canonical correspondence analysis was first carried out by using CANOCO software.Canonical correspondence analysis searches for multivariate relationships between two data sets (e.g., genetic data and species environmental data) (55).
Only oocysts for which all seven microsatellite loci could be scored were considered.The significance of the canonical axes was tested with a Monte Carlo permutation test (24), which allows estimation of the 95% C.I.s of the centroids of each population.A second analysis used FSTAT version 2.9.4., performing 20 random samplings of two equal sets of oocyst subpopulations within each Anopheles species in each site, to estimate and assess the significance of 20 F ST pair values for each species.These intraspecies F ST values were compared with the interspecies F ST values of the oocyst subpopulations of the two species in each site.
Linkage disequilibrium between pairs of loci was measured with the correlation coefficient R generalized for the multiallelic case.Statistical significance was assessed by a permutation test: Genotypes at two loci are associated at random a number of times, and the log-likelihood G test statistic is recalculated on the randomized data.We applied a Bonferroni correction to the P values (P value times the number of tests), as implemented in FSTAT version 2.9.4.

Fig. 1 .
Fig. 1.Inbreeding coefficient statistics for P. falciparum oocysts collected from A. gambiae (G) and A. funestus (F) mosquitoes in Simbock and Tibati (Cameroon) and in Rota (Kenya).(A) F IS measures the deficiency of heterozygotes resulting from nonrandom union of gametocytes within individual mosquito guts; bars indicate the 95% C.I. for each F IS value.(B) FST measures the deficiency of heterozygotes resulting from nonrandom distribution of oocyst genotypes among mosquito guts.L1 to L7 are the following seven microsatellite loci, respectively: POLYa, TA60, ARA2, Pfg377, PfPK2, TA87, and TA109 (ref.19).

FST
and 95% C.I. (15,000 bootstrap simulations over loci) are shown.Departure of FST from zero was tested by using the Bonferroni-corrected P value as implemented in FSTAT version 2.9.4.*Statistically significant (P Ͻ 0.001).

Fig. 2 .
Fig. 2. Results of the canonical correspondence analysis.Relative contribution of the variables ''site'' and ''Anopheles species'' to the genetic structure of P. falciparum populations infecting the two Anopheles species in the three sites.Oocyst populations were projected on the first two axes of the canonical analyses.Centroids (dots) of each parasite population are surrounded by the 95% C.I.s (ellipses).

Mosquito
Collections and Field Processing.Anopheles mosquitoes were collected in Simbock in November 2002, March 2003, and November 2004; in Tibati in September 2003; and in Rota from February to May in 2002 and 2003 (two replicate samples).

Table 1 . Distribution of P. falciparum oocysts in A. gambiae and A. funestus in three African sites and different years
7988 ͉ www.pnas.org͞cgi͞doi͞10.1073͞pnas.0702715104Annanet al.
than in A. gambiae (68).In Simbock, 20 of 21 pairs display significant linkage disequilibrium in 2002, 19 of 21 pairs in 2003, and 19 of 21 pairs in 2004.In Rota, all 21 pairs were in linkage disequilibrium in 2002 and 2003.Thus, linkage disequilibrium is stable over years in each site, corroborating that the ''clonal'' population structure persists over the generations.