- Open Access
Genetic structure of Anopheles gambiae populations on islands in northwestern Lake Victoria, Uganda
Malaria Journal volume 4, Article number: 59 (2005)
Alternative means of malaria control are urgently needed. Evaluating the effectiveness of measures that involve genetic manipulation of vector populations will be facilitated by identifying small, genetically isolated vector populations. The study was designed to use variation in microsatellite markers to look at genetic structure across four Lake Victoria islands and two surrounding mainland populations and for evidence of any restriction to free gene flow.
Four Islands (from 20–50 km apart) and two surrounding mainland populations (96 km apart) were studied. Samples of indoor resting adult mosquitoes, collected over two consecutive years, were genotyped at microsatellite loci distributed broadly throughout the genome and analysed for genetic structure, effective migration (Nem) and effective population size (Ne).
Ne estimates showed island populations to consist of smaller demes compared to the mainland ones. Most populations were significantly differentiated geographically, and from one year to the other. Average geographic pair-wise F ST ranged from 0.014–0.105 and several pairs of populations had Ne m < 3. The loci showed broad heterogeneity at capturing or estimating population differences.
These island populations are significantly genetically differentiated. Differences reoccurred over the study period, between the two mainland populations and between each other. This appears to be the product of their separation by water, dynamics of small populations and local adaptation. With further characterisation these islands could become possible sites for applying measures evaluating effectiveness of control by genetic manipulation.
Malaria kills over a million people annually, most from sub-Saharan Africa . Additionally, malaria mortality is on the rise, largely because of the emergence over the past two decades of widespread Plasmodium resistance to affordable antimalarial drugs . Control approaches such as insecticide impregnated bed nets are also being challenged by the emergence of insecticide resistance in Anopheles gambiae and Anopheles funestus, the two primary malaria vectors in sub-Saharan Africa [3, 4].
An alternative malaria control strategy being investigated in a number of laboratories is to genetically modify the vectorial capacity of vector populations by driving a genetic construct into the natural population. Genes that influence blood meal host selection, mosquito longevity, or Plasmodium survival have all been considered in genetic control, but most work has mainly focused on the identification of target genes that could modify the mosquito's ability to support Plasmodium sporogonic development [5–10]. The overall genetic control strategy depends not only on the identification and isolation of target genes but also on the development of effective transformation and drive systems and the development of potential field testing sites with vector populations that have been well characterized from the perspective of population biology and genetics. Although major advances are evident in genome resource development [11–13], target gene discoveries [14–18] and in genetic tool development [18, 19], less progress has been made in characterizing vector populations in potential field trial sites.
Studies of population genetic structure are vital to any vector-targeted control measure, especially where A. gambiae is one of the vectors . This species has a distribution that covers almost all of sub-Saharan Africa and genetic differentiation across populations of A. gambiae in Africa is complex. Microsatellite-, allozyme-and mitochondria-based studies have suggested extensive gene flow between populations in Senegal and western Kenya, a geographical distance of 6,000 km [21, 22]. In contrast, analyses of frequencies of paracentric chromosomal inversions and ribosomal DNA markers have revealed high levels of population structure within sympatric populations of A. gambiae in West Africa [23–26]and high differentiation has been observed within Kenya across distances of 700 km traversing the Rift Valley . In addition, A. gambiae island populations in Sao Tome  and A. arabiensis from the islands of Madagascar, Mauritius and Reunion have also shown extensive differentiation . It is not clear if the lack of extensive differentiation among A. gambiae populations across wide geographical distances (Senegal and Kenya) is due to high rates of gene flow among large populations or shared ancestral polymorphisms from a recent population expansion event . Physical barriers such as large areas of water and the Rift Valley are implicated in some instances where populations are highly differentiated, but chromosome inversion and molecular data also show clear evidence of pre-mating barriers producing reproductive isolation among sympatric populations [31, 32].
This study is focused on population structure of A. gambiae on islands in Lake Victoria, a part of Africa where A. gambiae populations are generally thought to consist exclusively of the Savanna chromosomal form and the S molecular form. The purpose of this study was to use variation in microsatellite markers to investigate the genetic structure of populations of A. gambiae s.s on several islands in northern Lake Victoria with a view to determining whether geographic separation of these islands (from 20–50 km) was associated with any evidence suggesting restriction to gene flow. Chen and others  in a study very similar to this in objective, design and geographic area, looked at A. gambiae populations on islands 2.5–21 km apart in eastern Lake Victoria. Genetic structuring among the island populations and between islands and surrounding mainland populations was still detectable, though low. Isolated populations are potentially useful sites for studies to evaluate the potential impact of malaria control measures that involve genetic manipulation of natural vector populations.
Study sites and field collections
The study area lies in Uganda, sub-Saharan Africa where malaria is endemic. Six A. gambiae populations involving one from each of four islands in northern half of Lake Victoria and two from the surrounding southern Uganda mainland were studied (Figure 1). The two inland populations consisted of Entebbe (EB), a peninsular jutting into Lake Victoria and Wamala (WL) located by the shores of a small inland lake 96 km away from Entebbe. The four islands are Nsadzi (NZ), Bugala (BL), Sserinya (SY) and Bukasa (BK). The islands are remote, but variable in size and ease of accessibility from the mainland and each other. Bugala, the largest, can be accessed from mainland by small boats and ferry whereas Nsadzi is the smallest and accessed only by boat. Bukasa lies farthest from the mainland sites. Apart from Entebbe, which is mainly residential, the rest of the locations are inhabited with people living by traditional farming subsisted with a little of fishing. Indoor resting adults from each population were captured at two or three separate villages by insecticide spraying done between 6 and 7 am. Populations were sampled as year (yr) 1 and as year 2 collections within a period of one to two years. Year 1 collections were made between November 2001 and February 2002. Year 2 collection, a replicate effort was performed between December 2002 and May 2003.
A. gambiae sensu lato (s.l.) of both sexes morphologically identified from other anophelines, based on a species identification key , were preserved in 80% alcohol and sent to the Center for Tropical Disease Research and Training-, University of Notre Dame (USA), for molecular identification and further analysis.
Molecular species identification and marker genotyping
Genomic DNA was extracted from single mosquitoes by kit, using procedures in Mukwaya et al. , or, for yr 2, in 96 well-plates (Wizard SV-96 Genomic DNA Purification System, Promega) processed with a Biomek FX workstation (Beckman Coulter). Molecular species identification (PCR) was according to Scott et al . Individuals that either did not amplify or gave incorrect sized product were excluded from subsequent genotyping analysis. This left 32 (Entebbe), 20 (Wamala), 33 (Bukasa), 32 (Sserinya), 36 (Bugala) and 36 (Nsadzi) individuals for the yr 1 sample set and 43 (Entebbe), 45 (Wamala), 47 (Bukasa), 20 (Sserinya), 92 (Bugala), 47 (Nsadzi) for yr 2. Yr 1 individuals were genotyped for variation at 17 microsatellite loci. Yr 2 replicate effort consisted of a 10 loci subset of the 17. Most loci used have been described elsewhere [21, 37, 38]. For those not previously described additional details are provided (see additional files 1, 2, 3).
Genotyping PCR was as follows: each 25 μl reaction contained 3.75 ng genomic DNA, 125 mM KCl, 25 mM Tris-HCl, PH 8.3, variable concentrations of MgCl2, 0.2 mM dNTP (Invitrogen, Carlsbad, CA) 0.011 mM each of either Fam, Tet and Hex or Blue, Green and Black Beckman coulter dye tagged forward primer and unlabeled reverse primer (Gibco/Brl, Gaithersburg, Md or Proligo LLC, Boulder, Co or Invitrogen) and 0.25 μl of home-made Taq DNA Polymerase. Yr 1 amplification was from GeneAmp 9600, whereas GeneAmp 9700 thermocycler (Applied Biosystems) was used for yr 2. The cycling program consisted of one cycle at 96°C, 5 minutes; thirty-five cycles of 94°C, 30 seconds; 55°C or optimal, 20 seconds; 72°C, 30 seconds; and one cycle of 72°C, 5 minutes. The Fam-Tet-Hex labeled PCR products constituted five of the 17 yr 1 loci set and were fragment size scored on the ABI 377 automatic sequencer using default settings of the genotyper software (Applied Biosystems). The remaining 12 loci of the data set were pool-plexed (two groups of six loci each) and genotyped using dye-labelled chemistry on the CEQ 8000 Beckman-Coulter capillary array genetic analysis system. Yr 2 were also pool-plexed into two groups (one of four and other of 6 loci) and similarly genotyped. A pool comprised; 1 μl product of each of 6 PCR reactions, 0.5 μl of a 400 bp size standard (Beckman-Coulter) and 30 μl SLS buffer (Beckman-Coulter). Both genotypers generate output fragment/allele sizes that are of within system reproducible non-integer lengths. Sizing of the outputs into integer length format useable by input files of the various genetic analysis programs is necessary. All Beckman-Coulter run samples were sized by binning, an automated process that relies on prior knowledge of the spectrum range of most possible apparent sizes for the generation of nominal fragment length sizes, with CEQ8000 software. This created an allele list that was used repeatedly to identify alleles whenever a locus was run under the same conditions. Sized alleles were manually inspected for correctness. Proper use of the binning option is described in the CEQ 8000 Genetic Analysis System User's Guide (Beckman-Coulter PN 608315).
Within population deviations from Hardy-Weinberg (HW) expectations at each locus were tested by exact tests using an online (web) version of GENEPOP an update of version 1.2  and also by ARLEQUIN . Input files for both programs were conversions from the program Microsatellite Analyser (MSA) . Conformity to Hardy-Weinberg expectations [H0 = of random union of gametes] was tested using the probability test. The possibility that heterozygosity deficiency may be the cause for departure from expectations was determined by setting the GENEPOP option [H1 = heterozygote deficiency]. To identify and correct genotyping errors in the data set the program MICRO-CHECKER  was used. Wherever presence of null alleles was suggested the data set adjustment procedure was accordingly applied to correct allele and genotype frequencies. The null-allele-adjusted data set was then used to explore the effect of null alleles on differentiation values resulting from the analysis. Linkage disequilibria, tests for independence between loci pairs, were done with web GENEPOP. Significance came from probability tests generated using Markov chain method at default parameter settings. Assessments of the six population deme sizes were achieved through estimations of effective population size (Ne) calculated from genetic data using the program MLNE . The single isolated population option was used. Ne calculations by hand were performed to verify the MLNE results. Equations used in the hand calculations have been adequately described [28, 44, 45]. Essentially, current Ne, an estimate based on temporal variation in allele frequencies from one sampling time to another, was calculated across the ten shared yr 1 and yr 2 loci. The allele frequencies for both data sets were from MSA basic descriptive statistics outputs. The allele frequency change variance estimator Fc was chosen over Fa because it is less affected by the presence of an allele at time t but not time 0, and over Fk for its superior Ne estimation when > 3 alleles per locus are present. Fc was calculated according to Nei and Tajima  and was weighted for multiple loci using equation (8) in Tajima and Nei ;Waples  before substituting it into equation (11) in Waples  to get Ne. Twelve generations per year was adopted for t in equation 11 above. The presence of genetic differences across populations was determined from three measures of genetic variability; genic differentiation that tests for allelic distribution and genotypic differentiation for genotypic distribution; both done with GENEPOP. The third measure looked for variation in frequencies of observed heterozygosity among populations. This was done with the Friedman test from the statistical program package SPSS. The measures described only show the presence or absence of differences. For magnitude of differences or population structure three indices of differentiation were performed; multi-loci population pair-wise Wright's F-statistics (F ST); R ST  an index that differs from F ST mainly in assumption for model of microsatellite evolution; and Nm an index of migration rate. Pair-wise F ST were generated using MSA, R ST were got from the program ARLEQUIN  and Nem were estimated from formula; F ST = 1/1+4Nem adopted from equation 5.17 . To further evaluate structure results, population pair-wise yr 1 and yr 2 F ST distributions were compared using paired t-test and the Wilcoxon-signed rank test from the SPSS package. Isolation by distance as the model explaining the observed population structure was tested by regression of Pair-wise Population F ST/(1 - F ST) against natural logarithm (ln) of pair wise geographical distances (Spearman Rank Correlation Test). The procedure was carried out online as computed in GENEPOP. Significance of the correlation coefficient was from Mantel tests. The geographical distances used were straight-line measurements between map points.
Population composition, HW proportions and independence of loci
Molecular species identification  showed all samples that generated a PCR product, except some from Bukasa, to be A. gambiae. In Bukasa, all year one (yr 1) samples were A. gambiae, while the yr 2 collection was composed of about 80% A. gambiae and 20% Anopheles arabiensis. Within population Hardy-Weinberg (HW) equilibrium tests (Ho = random union of gametes, H1 = heterozygote deficit) found eight of 17 yr 1 loci in HW equilibrium across all populations. H544 was the only locus out of equilibrium in every population. The equilibrium status of the remaining 8 loci varied in a population dependent manner (see additional files 1, 2, 3). The Wamala population had the fewest loci departing from HW equilibrium, with only1 of the 17 with a heterozygote deficit. Bugala had the highest levels of departure from HW equilibrium, with six of 17 out of equilibrium. Yr 2 exhibited some deviations from equilibrium as well with significantly positive F is values in 17 of 60 tests. These HW deviations in both data sets indicated heterozygote deficiencies. MICRO-CHECKER, a program that statistically discerns out HW equilibrium errors resulting from null alleles from those by inbreeding or Wahlund effects based on distinctive allele class distribution signatures that each error carries , attributed all observed loci heterozygote deficiencies to null alleles. Linkage disequilibrium (LD) tests for loci pairings across the six populations were overall insignificant (P > 0.05) except in three out of 136 (2%) pairings for yr 1. The three loci pairs that showed non-random association were H93 vs 29C1, H117 vs H544 and H117 vs MBP1B. All loci pairings used in yr 2 showed random association (LD tests P > 0.05).
Population genetic variability and differentiation
The loci were highly polymorphic in all populations as seen from number of alleles and heterozygosities (additional files 1, 2, 3). Although there were no significant across population differences in mean observed heterozygosities (Ho) in both years (Friedman test: χ2 0.05,5,17 = 5.662, P = 0.340 for yr 1; yr 2 was similar), differences in allele composition and manner of pairing were evident from the highly significant genic and genotypic differentiation all P <<0.0001. Genic and genotypic differentiation tests are for allelic and genotypic distributions across populations, with the null hypothesis being (H0 = distribution identical across populations).
The effective population size (Ne), which is the size of an ideal population that behaves, with respect to allele fluctuations, like the observed real population, was estimated from the program MLNE. . The Ne estimates showed differences in deme sizes between island and mainland populations (Table 1). The islands consisted of much smaller A. gambiae population sizes compared to mainland. Hand calculated Ne estimates (not shown) corroborated the MLNE values.
Degrees of genetic differentiation and population structure
Multilocus yr 1 F ST comparisons between population pairs revealed significant differentiation (Table 2). The across years population comparisons revealed substantial subdivision, except for the two mainland sites, in that comparisons of a particular location a certain year to itself another year were no lesser differentiated than those to different locations another year (Table 3). Likewise within yr 1 versus within yr 2 population pair comparisons comprised numerous instances of F ST variation in magnitudes (Table 4), even though statistically the yr 2 F ST distributions couldn't be shown to significantly differ from those of yr 1 (P = 0.119, Wilcoxon signed ranks test). The yr 1 F ST distribution from a survey across the 10 loci used in year 2 (Table 4) was not significantly different from the distribution calculated using all 17 loci (t 0.05,14 = 0.05, p = 0.961, paired t-test). MICRO-CHECKER null allele adjusted data sets, when re-analysed for F ST gave similar levels of population differentiation as the unadjusted ones. Global F ST differentiation across combined all yr 1 loci among the four islands (F ST = 0.042, P <0.001) was comparable to that between island and mainland populations (F ST = 0.044, P <0.001) and only a little lower than was observed between the two mainland populations (F ST = 0.054, P <0.001) (see Table 5). The study loci were broadly spread across the genome and varied in their ability to capture inter population differences. Three adjacent study loci, MBP1A, MBP1B and 22C1, on the left arm of chromosome 2 starkly stood out from the others at capturing extreme population genetic differentiation values, all across except between island and mainland comparisons (Table 5). These three loci lie in the 2La inversion at the proximal end and around its breakpoint neighborhood. When those three and H79 on 2R, the other inversion spanning locus, were excluded from the analysis, the between mainland population differences and the among islands differences substantially dropped leaving the between mainland and island and comparisons involving Bukasa as the remaining appreciable differentiations (Table 6). Moreover, H79, MBP1A, MBP1B and 22C1 alone account for nearly all the drop in F ST values observed when all null allele associated loci were excluded from the analysis (additional file 4).
Estimates across all the 17 yr 1 loci of the effective migration (Nem) showed the existence of structuring with varying degrees of restriction to gene flow between population pairs (Table 7). Geographical distance as the main factor explaining differentiation patterns was found to be insufficient. The observed population structure was not compatible with the isolation by distance model when regression between F ST/(1 - F ST) versus ln distance was evaluated (Mantle test; P = 0.787), in that there was little correlation between geographical distance and degree of differentiation (Fig 2).
The studied samples consisted of indoor resting, insecticide spray-catch specimens. Although there have been occasional indicators from other studies of A. gambiae that certain genotypes are associated with different resting behaviors [50–52], overall the A. gambiae populations in East Africa are panmictic, even taking into account different resting behaviors . So it can be taken that indoor sampling was adequately representative.
Neutrality from selection and genetic independence of loci used in genetic studies are required prior to analysing genetic variation at multiple microsatellite loci for population structure. Three pairings involving 5 loci in this study showed nonrandom association. All loci used in the study have known chromosome map locations (Figure 3). It is likely that H93 and 29C1 are unlinked because they are one chromosomal subdivision apart and located at the telomeric end, a region of chromosome where recombination is less restricted. However, in the islands population study by Chen et al  linkage disequilibria was also found among some of their loci pairs so it is plausible H93 and 29C1 linkage disequilibria could be quite incomplete through hitch-hiking to a nearby gene under selection. There is no direct genetic evidence to support this though. H117 and MBP1B although situated on the same chromosome arm, the two loci are far apart and sit in different chromosomal environments. H117 sits on telomeric end whereas MBP1B is located in an inversion and for standard arrangement more than six divisions upstream (Figure 3). Therefore, little possibility for linkage is expected, be it in the standard or inverted arrangement. H117 and H544, the last of non-freely associating pairs, map to different chromosomes and hence are not in the same linkage group so they are more likely to be unlinked. Finally, three instances of significance out of 136 tests (~2%), as is the case for this data set, are not above the range expected by chance alone at α = 0.05.
Deviations from HW were registered at certain loci. Deviation from frequencies expected from HW is not uncommon and while a potential indicator for selection at a locus  it is considered unlikely in most of the loci, as majority of them (15 of the 17, see additional file 1) have previously been used without evidence of selection. Moreover, departure from HW can arise from a variety of other causes including presence of null alleles [55, 56], hidden sub-structure and inbreeding in a population . These collections were made from more than one village so patchy distribution within each population could, if present, affect the equilibrium. Little is actually known about breeding behavior, deme sizes and distribution in these populations. Although some slight inbreeding has recently been suggested for natural A. gambiae populations in East Africa , which if present could account for the deviations, the expected associated inbreeding signature of genome-wide departures from HW equilibrium was not found. Inbreeding in these samples being the cause of non-equilibrium was ruled out due to lack of such genome-wide departures from HW equilibrium in any of the populations. The observed HW equilibrium departures were locus specific. Moreover, earlier studies on other East African populations found random mating [21, 51]. The HW deviations were attributed to null alleles by the MICRO-CHECKER program. This program statistically discerns out HW equilibrium errors resulting from null alleles from those by inbreeding or Wahlund effects as each carries a distinct allele class distribution signature. In fact such, locus-specific, HW deviation patterns resulting from null alleles have been previously encountered by other investigators. Donnelly et al  found null alleles responsible for 5 of 6 loci HW deviations while studying structure in A. arabiensis, whereas Lehmann et al  had 4 of 9 loci showing some instances of non-equilibrium in their A. gambiae study. The impact of null alleles in the data on the analysis was negligible based on the fact that re-analysis of adjusted data sets returned similar differentiation values.
Populations, other than the mainland ones, were significantly differentiated across the years to the extent that they were substantially different even from themselves, from one year to another. This is evidence, in these populations, for demographic instability probably emanating from seasonal changes and is indicative of small population sizes on the islands. In spite of overall differentiation across the years; the within yr 1 F ST distribution when arrayed against the within yr 2 F ST distribution did not statistically significantly differ according to the Wilcoxon-signed ranks test perhaps because of some population pair differences that were exactly recaptured a year later. The yr 1 vs yr 2 irregularity of appearance of A. arabiensis, in Bukasa samples is probably a sampling-time artifact that probably caught them out of synch. Yr 1 one samples were collected in the months of November through February, a period that falls in the dry season, whereas the yr 2 collections spanned through a dry and wet season (see methods). Bukasa yr 1 samples were collected November/December 2001, while yr 2 got collected during months of April/May 2003. These samples, while spatially true replicates, were not replicates temporally. Relative frequencies of various members of the A. gambiae complex are known to fluctuate with season and geographical location .
Effective population size (Ne) comparisons across populations are usually not factored into structure analysis due to lack of reliable direct methods of estimates [27, 44]. The study generated indirect Ne estimates show that the islands on the whole have lower deme sizes compared to the mainland. The island Ne's were in the hundreds, whereas mainland effective populations sizes were in the thousands, a result that is consistent with the conclusion arrived at earlier that small population sizes exist on these islands. In contrast, the western Kenya island study  inferred a large effective population size, in both, the islands and mainland, based on their comparable degrees of polymorphism in terms of average number of alleles and levels of observed heterozygosity. However, Nes inferred that way are only qualitative and do not take into account actual allele constitution or make up the way changes in individual allele frequencies in the method of Ne calculation used this study does. Therefore, the present study's Nes because of their being quantitative are more exact. A previous study on A. gambiae population size in Kenya  corroborates the large mainland Ne estimates.
Within population genetic diversity was high both on the islands and the mainland considering heterozygosity levels and the number of alleles seen (additional files 1, 2, 3). Across population differentiation, with respect to allele frequencies and genotype constitution, was high in all cases. The level of genetic differentiation among islands and mainland populations was considerable according to multi-loci pair-wise F ST (Table 2). F ST and R ST both estimate the amount of differentiation but each suits different scenarios. F ST assumes infinite allele mutational (IAM) model while R ST assumes and requires strict adherence to a step-wise mutation (SMM) model for microsatellite evolution . Of the repeat motif classes in the marker sets used only the tri-nucleotide (3 bp) repeat loci satisfactorily conformed to the SMM with regard to generating products consistent with a series predictable from the repeat motif inside a constant flanking sequence; because several alleles among the dinucleotide loci appeared to be separated by only one nucleotide which leads to inconsistencies and mis-scoring. Therefore, F ST values were regarded as the more robust ones. Low but significant genetic structure was found among the island population (F ST = 0.019) and between island and mainland populations (F ST = 0.003) situated from 3–20 km apart in the Western Kenya-Lake Victoria study . These Ugandan island populations situated 20–50 km apart are more differentiated (Table 5) than those in the Western Kenya Lake Victoria island study, perhaps due to the longer separation distances involved. Across 17 loci, the observed levels of differentiation among the island populations did not much differ from those seen between islands to mainland or between the two mainland populations. However, this effect was not identical genome-wide in that all loci did not capture it to the same extent. They greatly varied in their ability to capture inter-population differences. Among those loci that captured significant group differences (Table 5) it is apparent that each had its own independent differentiation rate across the populations. The loci in the inversions particularly the three involved with 2La extremely differentiate the populations. Excluding them from the analysis substantially drops most inter island differences and the inter mainland pair difference although island to mainland difference and Bukasa differences are less affected (Table 6). Although the effect of inversions on gene flow in A. gambiae is unknown the above result points to possible role of inversion situated loci in driving population differentiation. In fact, 2La and some 3R inversions have shown clines with aridity [23, 59, 60] and association with particular resting behaviors[51, 61], such that genes within them are probably involved in environmental adaptations. Site ecological differences are evident across these populations. The islands are mostly forested and covered in rush green natural vegetation. The inland Lake Wamala population lies in a lesser naturally-vegetated, drier, wooded grassland-like region with farm crops. The mainland Entebbe area is somewhat intermediate; a peninsular extending from a forested mainland on one end and becoming less vegetated heading towards the lake. While it is possible in light of the above that some of the observed variation between the populations is shaped by differential adaptation and small population size effects, the rest, at mutation equilibrium, is then accounted for by restrictions to gene flow. This gene flow restriction is not likely to arise from chromosomal form diversity because populations in this region are thought to consist of only the savanna form. It is possibly arising from barriers to dispersal.
The indices of effective migration (Nem) indicate that gene flow is indeed substantially though not completely restricted, between many pairs (Table 7). It is strongly evident that the nature of the barrier responsible for the observed population structure has less to do with sheer geographical distance (Figure 2), than with water separation: Entebbe peninsular is geographically farther from the inland Wamala population than from any island population, however, it is less isolated genetically from Wamala than from any of the islands. The distances separating these populations (see Figure 1) are beyond both the normal 1 km A. gambiae flight range  and 7 km wind assisted flight range . While it is not absolutely inconceivable that wind could be a factor in this, mosquito dispersal between these populations is more likely to be man assisted. However, conclusions about effective migration levels derived from Nem values should be interpreted with care for several reasons: Foremost, Nem were indirectly estimated from F ST. The relationship between Nem and F ST is non linear so any errors in F ST are magnified in Nem. Secondly, although an F ST gives a measure of relative amount of differentiation between a population pair it is still confounded by time in sense that the derived Nem is based on structure that has been generated over many generations so cannot distinguish recurrent from ancestral gene flow. Actually it is advised that all indirectly calculated migration rates be viewed cautiously . This study had scope to primarily study differentiation and not to measure present active migrations or actual dispersals (Nm) between populations and so the Nem are only portrayals of gene flow rates in terms of effective migration in light of the observed levels of differentiation. To get actual dispersal or migration levels would require use of direct methods of acquisition such as capture-recapture. The cost of these direct methods has become affordable in recent years .
It was found that these island populations in North Western Lake Victoria region are substantially differentiated from the mainland and some of each other. It also is that this differentiation is strongly shaped by physical barriers to dispersal or gene flow, processes associated with small population sizes and possibly also by ecological adaptation because the levels of differentiation found contrast starkly with what has mostly been reported for A. gambiae populations around the continent. Most of the F ST were much higher than (F ST = 0.014) expected for mainland populations at similar range of separation distances [64–66]. The differentiations in several instances were more like those seen across the Rift Valley (mean F ST = 0.104, ); island populations in Sao Tome  and A. arabiensis amongst Madagascar, Mauritius and Reunion (F ST 0.08 - 0.215, ) that involved barriers to gene flow. Although not in complete genetic isolation since only gene flow from Nem levels of 2 and less could allow this , they are some of the most differentiated A. gambiae populations among those studied to date. This high differentiation and smaller population size confers to them some practical importance in fight against malaria because completely, or in their absence, even nearly isolated small vector populations could be used as field sites for evaluating impact of malaria control measures including those using genetic manipulations. However, before they are adopted for this role extensive additional studies must be carried out. There is need to establish for example, what the exact nature of the barrier is. Is it just water or is there more to it like some other yet unknown physical aspect? In this way, potential ways of its compromise could be monitored during duration of trials. It would be interesting to figure out the origin of the observed differentiation. It makes a huge difference to understand whether this is recurrent or historical gene flow. Among the recurrent processes involved it is crucial to know the relative significance of the factors at play. For instance, is the differentiation primarily driven by extinctions on islands followed by re-colonization from elsewhere or just drift fluctuation followed by recovery from extensive births without significant immigrants impacts. These pertinent studies could be done with use of markers that have lower mutation rates to microsatellites and are able to look farther back into the past and incorporating the findings with those from direct measures of present day migration rates. It is still intriguing that there is substantial differentiation amongst these populations in spite of possible passive mosquito dispersal (by human activity) across the barrier through ferry or boat traffic (Fig. 1). This could mean that passive dispersal, though commonly implicated, might not be as effective as widely thought. The role and extent of passive mosquito dispersal in natural conditions need to be empirically determined.
These lake islands are significantly genetically differentiated from the two mainland populations. Several of them are also differentiated from one another. The genetic differences are real for they reappeared in yr 2. These genetic differentiations are possibly the product of several factors: the islands physical separation across water, effects of their small population size and local ecological adaptation. Although the relative contribution of each differentiating factor is yet to be quantified, when done these islands could become candidate sites for measures evaluating effectiveness of control by genetic manipulation. Lastly, this study adds to the body of data that has found substantial structure among A. gambiae populations across physical barriers.
Marshall E: Malaria. A renewed assault on an old and deadly foe. Science. 2000, 290: 428-430. 10.1126/science.290.5491.428.
Wellems TE: Plasmodium chloroquine resistance and the search for a replacement antimalarial drug. Science. 2002, 298: 124-126. 10.1126/science.1078167.
Hemingway J, Field L, Vontas J: An overview of insecticide resistance. Science. 2002, 298: 96-97. 10.1126/science.1078052.
Weill M, Malcolm C, Chandre F, Mogensen K, Berthomieu A, Marquine M, Raymond M: The unique mutation in ace-1 giving high insecticide resistance is easily detectable in mosquito vectors. Insect Mol Biol. 2004, 13: 1-7. 10.1111/j.1365-2583.2004.00452.x.
James AA: Mosquito molecular genetics: the hands that feed bite back. Science. 1992, 257: 37-38.
James AA: Blocking malaria parasite invasion of mosquito salivary glands. J Exp Biol. 2003, 206: 3817-3821. 10.1242/jeb.00616.
Collins FH, Besansky NJ: Vector biology and the control of malaria in Africa. Science. 1994, 264: 1874-1875.
Crampton JM, Warren A, Lycett GJ, Hughes MA, Comley IP, Eggleston P: Genetic manipulation of insect vectors as a strategy for the control of vector-borne disease. Ann Trop Med Parasitol. 1994, 88: 3-12.
Alphey L, Beard CB, Billingsley P, Coetzee M, Crisanti A, Curtis C, Eggleston P, Godfray C, Hemingway J, Jacobs-Lorena M, James AA, Kafatos FC, Mukwaya LG, Paton M, Powell JR, Schneider W, Scott TW, Sina B, Sinden R, Sinkins S, Spielman A, Toure Y, Collins FH: Malaria control with genetically manipulated insect vectors. Science. 2002, 298: 119-121. 10.1126/science.1078278.
Nirmala X, James AA: Engineering Plasmodium-refractory phenotypes in mosquitoes. Trends Parasitol. 2003, 19: 384-387. 10.1016/S1471-4922(03)00188-0.
Holt RA, Subramanian GM, Halpern A, Sutton GG, Charlab R, Nusskern DR, Wincker P, Clark AG, Ribeiro JM, Wides R, Salzberg SL, Loftus B, Yandell M, Majoros WH, Rusch DB, Lai Z, Kraft CL, Abril JF, Anthouard V, Arensburger P, Atkinson PW, Baden H, de Berardinis V, Baldwin D, Benes V, Biedler J, Blass C, Bolanos R, Boscus D, Barnstead M, Cai S, Center A, Chaturverdi K, Christophides GK, Chrystal MA, Clamp M, Cravchik A, Curwen V, Dana A, Delcher A, Dew I, Evans CA, Flanigan M, Grundschober-Freimoser A, Friedli L, Gu Z, Guan P, Guigo R, Hillenmeyer ME, Hladun SL, Hogan JR, Hong YS, Hoover J, Jaillon O, Ke Z, Kodira C, Kokoza E, Koutsos A, Letunic I, Levitsky A, Liang Y, Lin JJ, Lobo NF, Lopez JR, Malek JA, McIntosh TC, Meister S, Miller J, Mobarry C, Mongin E, Murphy SD, O'Brochta DA, Pfannkoch C, Qi R, Regier MA, Remington K, Shao H, Sharakhova MV, Sitter CD, Shetty J, Smith TJ, Strong R, Sun J, Thomasova D, Ton LQ, Topalis P, Tu Z, Unger MF, Walenz B, Wang A, Wang J, Wang M, Wang X, Woodford KJ, Wortman JR, Wu M, Yao A, Zdobnov EM, Zhang H, Zhao Q, Zhao S, Zhu SC, Zhimulev I, Coluzzi M, della Torre A, Roth CW, Louis C, Kalush F, Mural RJ, Myers EW, Adams MD, Smith HO, Broder S, Gardner MJ, Fraser CM, Birney E, Bork P, Brey PT, Venter JC, Weissenbach J, Kafatos FC, Collins FH, Hoffman SL: The genome sequence of the malaria mosquito Anopheles gambiae. Science. 2002, 298: 129-149. 10.1126/science.1076181.
Mongin E, Louis C, Holt RA, Birney E, Collins FH: The Anopheles gambiae genome: an update. Trends Parasitol. 2004, 20: 49-52. 10.1016/j.pt.2003.11.003.
Hill CA, Kafatos FC, Stansfield SK, Collins FH: Arthropod-borne diseases: vector control in the genomics era. Nat Rev Microbiol. 2005, 3: 262-268. 10.1038/nrmicro1101.
Collins FH, Zheng L, Paskewitz SM, Kafatos FC: Progress in the map-based cloning of the Anopheles gambiae genes responsible for the encapsulation of malarial parasites. Ann Trop Med Parasitol. 1997, 91: 517-521. 10.1080/00034989760888.
Zheng L, Cornel AJ, Wang R, Erfle H, Voss H, Ansorge W, Kafatos FC, Collins FH: Quantitative trait loci for refractoriness of Anopheles gambiae to Plasmodium cynomolg i B. Science. 1997, 276: 425-428. 10.1126/science.276.5311.425.
Collins FH, Saunders RD, Kafatos FC, Roth C, Ke Z, Wang X, Dymbrowski K, Ton L, Hogan J: Genetics in the study of mosquito susceptibility to Plasmodium. Parassitologia. 1999, 41: 163-168.
de Lara Capurro M, Coleman J, Beerntsen BT, Myles KM, Olson KE, Rocha E, Krettli AU, James AA: Virus-expressed, recombinant single-chain antibody blocks sporozoite infection of salivary glands in Plasmodium gallinaceum-infected Aedes aegypti. Am J Trop Med Hyg. 2000, 62: 427-433.
Ito J, Ghosh A, Moreira LA, Wimmer EA, Jacobs-Lorena M: Transgenic anopheline mosquitoes impaired in transmission of a malaria parasite. Nature. 2002, 417: 452-455. 10.1038/417452a.
Catteruccia F, Nolan T, Loukeris TG, Blass C, Savakis C, Kafatos FC, Crisanti A: Stable germline transformation of the malaria mosquito Anopheles stephensi. Nature. 2000, 405: 959-962. 10.1038/35016096.
Collins FH, Kamau L, Ranson HA, Vulule JM: Molecular entomology and prospects for malaria control. Bull World Health Organ. 2000, 78: 1412-1423.
Lehmann T, Hawley WA, Kamau L, Fontenille D, Simard F, Collins FH: Genetic differentiation of Anopheles gambiae populations from East and west Africa: comparison of microsatellite and allozyme loci. Heredity. 1996, 77: 192-200.
Besansky NJ, Lehmann T, Fahey GT, Fontenille D, Braack LE, Hawley WA, Collins FH: Patterns of mitochondrial variation within and between African malaria vectors, Anopheles gambiae and An. arabiensis, suggest extensive gene flow. Genetics. 1997, 147: 1817-1828.
Toure YT, Petrarca V, Traore SF, Coulibaly A, Maiga HM, Sankare O, Sow M, Di Deco MA, Coluzzi M: The distribution and inversion polymorphism of chromosomally recognized taxa of the Anopheles gambiae complex in Mali, West Africa. Parassitologia. 1998, 40: 477-511.
Gentile G, Slotman M, Ketmaier V, Powell JR, Caccone A: Attempts to molecularly distinguish cryptic taxa in Anopheles gambiae s.s. Insect Mol Biol. 2001, 10: 25-32. 10.1046/j.1365-2583.2001.00237.x.
Mukabayire O, Caridi J, Wang X, Toure YT, Coluzzi M, Besansky NJ: Patterns of DNA sequence variation in chromosomally recognized taxa of Anopheles gambiae: evidence from rDNA and single-copy loci. Insect Mol Biol. 2001, 10: 33-46. 10.1046/j.1365-2583.2001.00238.x.
Favia G, Lanfrancotti A, Spanos L, Siden-Kiamos I, Louis C: Molecular characterization of ribosomal DNA polymorphisms discriminating among chromosomal forms of Anopheles gambiae s.s. Insect Mol Biol. 2001, 10: 19-23. 10.1046/j.1365-2583.2001.00236.x.
Lehmann T, Hawley WA, Grebert H, Danga M, Atieli F, Collins FH: The Rift Valley complex as a barrier to gene flow for Anopheles gambiae in Kenya. J Hered. 1999, 90: 613-621. 10.1093/jhered/90.6.613.
Pinto J, Donnelly MJ, Sousa CA, Malta-Vacas J, Gil V, Ferreira C, Petrarca V, do Rosario VE, Charlwood JD: An island within an island: genetic differentiation of Anopheles gambiae in Sao Tome, West Africa, and its relevance to malaria vector control. Heredity. 2003, 91: 407-414. 10.1038/sj.hdy.6800348.
Simard F, Fontenille D, Lehmann T, Girod R, Brutus L, Gopaul R, Dournon C, Collins FH: High amounts of genetic differentiation between populations of the malaria vector Anopheles arabiensis from West Africa and eastern outer islands. Am J Trop Med Hyg. 1999, 60: 1000-1009.
Donnelly MJ, Licht MC, Lehmann T: Evidence for recent population expansion in the evolutionary history of the malaria vectors Anopheles arabiensis and Anopheles gambiae. Mol Biol Evol. 2001, 18: 1353-1364.
Tripet F, Toure YT, Taylor CE, Norris DE, Dolo G, Lanzaro GC: DNA analysis of transferred sperm reveals significant levels of gene flow between molecular forms of Anopheles gambiae. Mol Ecol. 2001, 10: 1725-1732. 10.1046/j.0962-1083.2001.01301.x.
Tripet F, Toure YT, Dolo G, Lanzaro GC: Frequency of multiple inseminations in field-collected Anopheles gambiae females revealed by DNA analysis of transferred sperm. Am J Trop Med Hyg. 2003, 68: 1-5.
Chen H, Minakawa N, Beier J, Yan G: Population genetic structure of Anopheles gambiae mosquitoes on Lake Victoria islands, west Kenya. Malar J. 2004, 3: 48-10.1186/1475-2875-3-48. [http://www.malariajournal.com/content/3/1/48]
Gillies MT, De Meillon B: The Anophelinae of Africa south of the Sahara. 1968, Johannesburg: Publications of the South African Institute for Medical Research No 54, 2
Mukwaya LG, Kayondo JK, Crabtree MB, Savage HM, Biggerstaff BJ, Miller BR: Genetic differentiation in the yellow fever virus vector, Aedes simpsoni complex, in Africa: sequence variation in the ribosomal DNA internal transcribed spacers of anthropophilic and non-anthropophilic populations. Insect Mol Biol. 2000, 9: 85-91. 10.1046/j.1365-2583.2000.00161.x.
Scott JA, Brogdon WG, Collins FH: Identification of single specimens of the Anopheles gambiae complex by the polymerase chain reaction. Am J Trop Med Hyg. 1993, 49: 520-529.
Wang R, Kafatos FC, Zheng L: Microsatellite markers and genotyping procedures for Anopheles gambiae. Parasitol Today. 1999, 15: 33-37. 10.1016/S0169-4758(98)01360-X.
Zheng L, Benedict MQ, Cornel AJ, Collins FH, Kafatos FC: An integrated genetic map of the African human malaria vector mosquito, Anopheles gambiae. Genetics. 1996, 143: 941-952.
Raymond M, Rousset F: GENEPOP (Version 1.2): Population genetics software for exact tests and ecumenism. J Heredity. 1995, 86: 248-249.
Schneider S, Roessli D, Excoffier L: ARLEQUIN, Version 2.000: A software for population genetics data analysis. 2000, Produced by the Genetics and Biometry Laboratory, University of Geneva, Switzerland
Dieringer D, Schlotterer C: Microsatellite analyzer (MSA): a platform independent analysis tool for large microsatellite data sets. Mol Ecol Notes. 2003, 3: 167-169. 10.1046/j.1471-8286.2003.00351.x.
Van Oosterhout C, Hutchinson WF, Wills DPM, Shipley P: MICROCHECKER: Software for identifying and correcting genotyping errors in microsatellite data. Mol Ecol Notes. 2004, 4: 535-538. 10.1111/j.1471-8286.2004.00684.x.
Wang J, Whitlock MC: Estimating effective population size and migration rates from genetic samples over space and time. Genetics. 2003, 163: 429-446.
Waples RS: A generalized approach for estimating effective population size from temporal changes in allele frequency. Genetics. 1989, 121: 379-391.
Lehmann T, Hawley WA, Grebert H, Collins FH: The effective population size of Anopheles gambiae in Kenya: implications for population structure. Mol Biol Evol. 1998, 15: 264-276.
Nei M, Tajima F: Genetic drift and estimation of effective population size. Genetics. 1981, 98: 625-640.
Tajima F, Nei M: Note on genetic drift and estimation of effective population size. Genetics. 1984, 106: 569-574.
Slatkin M: A measure of population subdivision based on microsatellite allele frequencies. Genetics. 1995, 139: 457-462.
Hartl DL, Clark AG: Principles of Population Genetics. 1997, Sunderland, MA: Sinauer Associates, Inc
Colluzi M, Sabatini A, Petrarca V, Di Deco MA: Chromosomal differentiation and adaptation to human environments in the Anopheles gambiae complex. Trans R Soc Trop Med Hyg. 1979, 73: 483-497. 10.1016/0035-9203(79)90036-1.
Petrarca V, Beier JC: Intraspecific chromosomal polymorphism in the Anopheles gambiae complex as a factor affecting malaria transmission in the Kisumu area of Kenya. Am J Trop Med Hyg. 1992, 46: 229-237.
Mnzava AE, Rwegoshora RT, Wilkes TJ, Tanner M, Curtis CF: Anopheles arabiensis and An. gambiae chromosomal inversion polymorphism, feeding and resting behaviour in relation to insecticide house-spraying in Tanzania. Med Vet Entomol. 1995, 9: 316-324.
Smits A, Roelants P, Van Bortel W, Coosemans M: Enzyme polymorphisms in the Anopheles gambiae (Diptera:Culicidae) complex related to feeding and resting behavior in the Imbo Valley, Burundi. J Med Entomol. 1996, 33: 545-553.
Garcia De Leon FJ, Chikhi L, Bonhomme F: Microsatellite polymorphism and population subdivision in natural populations of European Sea Bass (Dicentrarchus labrax Linnaeus, 1758). Mol Ecol. 1997, 6: 51-62. 10.1046/j.1365-294X.1997.t01-1-00151.x.
Callen DF, Thompson AD, Shen Y, Phillips HA, Richards RI, Mulley JC, Sutherland GR: Incidence and origin of "null" alleles in the (AC)n microsatellite markers. Am J Hum Genet. 1993, 52: 922-927.
Paetkau D, Strobeck C: The molecular basis and evolutionary history of a microsatellite null allele in bears. Mol Ecol. 1995, 4: 519-520.
Donnelly MJ, Cuamba N, Charlwood JD, Collins FH, Townson H: Population structure in the malaria vector, Anopheles arabiensis Patton, in East Africa. Heredity. 1999, 83: 408-417. 10.1038/sj.hdy.6885930.
Lehmann T, Licht M, Gimnig JE, Hightower A, Vulule JM, Hawley WA: Spatial and temporal variation in kinship among Anopheles gambiae (Diptera: Culicidae) mosquitoes. J Med Entomol. 2003, 40: 421-429.
Coluzzi M, Sabatini A, Petrarca V, DiDeco MA: Chromosomal differentiation and adpatation to human environments in the Anopheles gambiae complex. Trans R Soc Trop Med Hyg. 1979, 73: 483-497. 10.1016/0035-9203(79)90036-1.
Coluzzi M, Petrarca V, Di Deco MA: Chromosomal inversion intergradation and incipient speciation in Anopheles gambiae. Boll Zool. 1985, 52: 45-63.
Coluzzi M, Sabatini A, Petrarca V, Di Deco MA: Behavioural divergences between mosquitoes with different inversion karyotypes in polymorphic populations of the Anopheles gambiae complex. Nature. 1977, 266: 832-833. 10.1038/266832a0.
Costantini C, Li SG, Della Torre A, Sagnon N, Coluzzi M, Taylor CE: Density, survival and dispersal of Anopheles gambiae complex mosquitoes in a west African Sudan savanna village. Med Vet Entomol. 1996, 10: 203-219.
Whitlock MC, McCauley DE: Indirect measures of gene flow and migration: FST not equal to 1/(4Nm + 1). Heredity. 1999, 82: 117-125. 10.1038/sj.hdy.6884960.
Kamau L, Lehmann T, Hawley WA, Orago AS, Collins FH: Microgeographic genetic differentiation of Anopheles gambiae mosquitoes from Asembo Bay, western Kenya: a comparison with Kilifi in coastal Kenya. Am J Trop Med Hyg. 1998, 58: 64-69.
Lehmann T, Besansky NJ, Hawley WA, Fahey TG, Kamau L, Collins FH: Microgeographic structure of Anopheles gambiae in western Kenya based on mtDNA and microsatellite loci. Mol Ecol. 1997, 6: 243-253. 10.1046/j.1365-294X.1997.00177.x.
Lehmann T, Licht M, Elissa N, Maega BT, Chimumbwa JM, Watsenga FT, Wondji CS, Simard F, Hawley WA: Population structure of Anopheles gambiae in Africa. J Hered. 2003, 94: 133-147. 10.1093/jhered/esg024.
We acknowledge: Dr. Brandon J. Hacket at Purdue University, Dr. Martin J. Donnelly of Liverpool School of Tropical Medicine, Dr. Liangbiao Zheng at Yale University School of Medicine, Mr Fred Ssenfuka and entomology field collections team of Uganda Virus Research Institute (UVRI), Julie Niedibalski, Marcia Kern, Maria Unger, James Hogan, Jigar Patel from Center for Tropical Disease Research and Training, University of Notre Dame and the Evolutionary Discussion Group (EDG) at the University of Notre Dame for their insightful suggestions. This work was supported by funds from the Uganda Ministry of Health to LGM and NIH grant #PO1AI45123 to FHC.
JKK carried out study design, sample processing, data acquisition, analysis and interpretation and manuscript preparation. LGM conceived of the study. AS substantially participated in data analysis. APM was key in data collection techniques and analysis. MBC greatly helped draft the manuscript. NJB helped with marker selection resources. FHC participated in the design of the study and substantially helped draft the manuscript.
Electronic supplementary material
About this article
Cite this article
Kayondo, J.K., Mukwaya, L.G., Stump, A. et al. Genetic structure of Anopheles gambiae populations on islands in northwestern Lake Victoria, Uganda. Malar J 4, 59 (2005) doi:10.1186/1475-2875-4-59
- Gene Flow
- Null Allele
- Effective Population Size
- Island Population
- Small Population Size