The Journal of Heredity 2002:93(4)
© 2002 The American Genetic Association 93:249-253
Genetic Differentiation of Anopheles gambiae s.s. Populations in Mali, West Africa, Using Microsatellite Loci
From the Department of Organismic Biology, Ecology and Evolution, University of California, Los Angeles, CA 90095-1606 (Carnahan, Taylor, and Diuk-Wasser); School of Medicine, Epidemiology and Public Health, Yale University, New Haven, CT 06520 (Zheng); Malaria Research and Training Center, Faculté de Médecine, de Pharmacie et d'Odonto-Stomatologie, Bamako, B.P. 1805, Mali (Touré and Dolo); School of Hygiene and Public Health, Johns Hopkins University, Baltimore, MD 21205 (Norris); and Department of Pathology, University of Texas Medical Branch, Galveston, TX 77555-0609 (Lanzaro).
Address correspondence to Charles Taylor, UCLA Department of Organismic Biology, Ecology and Evolution, 621 Charles Young Dr. S, Los Angeles, CA 90095-1606, or e-mail: taylor{at}biology.ucla.edu.
| Abstract |
|---|
|
|
|---|
Anopheles gambiae sensu stricto is a principal vector of malaria through much of sub-Saharan Africa, where this disease is a major cause of morbidity and mortality in human populations. Accordingly, population sizes and gene flow in this species have received special attention, as these parameters are important in attempts to control malaria by impacting its mosquito vector. Past measures of genetic differentiation have sometimes yielded conflicting results, in some cases suggesting that gene flow is extensive over vast distances (6000 km) and is disrupted only by major geological disturbances and/or barriers. Using microsatellite DNA loci from populations in Mali, West Africa, we measured genetic differentiation over uniform habitats favorable to the species across distances ranging from 62 to 536 km. Gene flow was strongly correlated with distance (r2 = 0.77), with no major differences among chromosomes. We conclude that in this part of Africa, at least, genetic differentiation for microsatellite DNA loci is consistent with traditional models of isolation by distance.
| Introduction |
|---|
|
|
|---|
Anopheles gambiae sensu stricto, a principal vector of malaria, is widely distributed in Africa south of the Sahara Desert to northern South Africa. Throughout the region it is highly polymorphic for gene arrangements on the second chromosome (Bryan et al. 1982; Coluzzi et al. 1979). The genetic structure of this species is complex, with several recognized "chromosomal forms" that differ ecologically. These chromosomal forms are of uncertain taxonomic status; for the purposes of this article we will simply treat them as distinct units that potentially interbreed within themselves. In Mali, the site of our study, three chromosomal forms are widely distributed: Bamako, Savanna, and Mopti. These forms commonly occur in sympatry and their relative abundance varies greatly from location to location and season to season. Correlation studies have established that the composition by chromosomal form varies with climate (Touré et al. 1998b), the 2Rb arrangement is closely tied to the amount of vegetation (Thomson et al. 1997), and the canonical correlation between all gene arrangement frequencies and climate, especially dry-cold season temperature and precipitation, is very high (r > 0.88) (Massumi et al. 2000).
Geographic variation for other genetic markers has been less clear. In particular, for mitochondrial DNA (mtDNA) Besansky et al. (1997) found little or no population subdivision across considerable distances in East Africa and reported that there was no evidence for isolation by distance there. Likewise, Lehmann et al. (1996, 1997) found no differences in microsatellite DNA frequencies between populations of the Savanna form on the east and west coasts of Africa, thousands of kilometers apart. In that study, the one barrier to gene flow they did find separated populations on opposite sides of the Rift Valley complex in East Africa, a pattern also observed using mtDNA markers (Lehmann et al. 2000). At the same time Lanzaro et al. (1998) reported that populations in Mali separated by uniformly favorable habitat showed differences at microsatellite DNA loci comparable to those reported by Lehmann et al. (2000) for the Rift Valley.
Extensive karyotype data from Mali, West Africa, has shown populations to be geographically structured for frequencies of certain chromosomal arrangements (Touré et al. 1998b). While these data do not directly show populations of the same form to be isolated by distance, they do support the notion that such arrangements vary both seasonally and geographically. Unfortunately these chromosomal inversions are known to be nonneutral and therefore may not be adequate markers to test whether gene flow is restricted across geographic distance.
The purpose of this study was to inquire further into the question of geographic population genetic structure for microsatellite DNA loci. We looked at populations of An. gambiae s.s. in Mali, West Africa, separated by a range of distances of continuously suitable habitat, and examined whether these distances were correlated with differences in microsatellite allele frequencies.
| Methods |
|---|
|
|
|---|
An. arabiensis and An. gambiae female mosquitoes were collected from human habitations in six villages in Mali, West Africa: Banambani (12°48'N, 08°03'W), Mopti (14°31'N, 04°16'W), Pimperena (11°28'N, 05°42'W), Selinkenyi (11°42'N, 08°17'W), Soulouba (11°47'N, 07°43'W), and Tene (13°24'N, 04°35'W). Their locations are shown in Figure 1. The six villages are located along the Niger River valley and the geographic distance between them ranges from 62 to 536 km. Descriptions of the plant communities of this region may be found in White (1983) and Touré et al. (1998b). The vegetation throughout this region is quite uniform, as shown by the land cover map derived from remotely sensed data by the University of Maryland Global Land Cover Facility (Figure 1). It is based on NASA/NOAA Pathfinder land datasets (AVHRR sensor) with a spatial resolution of 1 km, collected from 1981 to 1994. Land cover was inferred using a supervised classification procedure from the Landsat Multispectral Scanner System, Landsat Thematic Mapper, and Linear Imaging Self-Scanning Sensors (Hansen et al. 2000). The area included in our samples was almost continuously covered by grassland, shrubland, or woodland. These are all known to support high densities of An. gambiae (Touré et al. 1998b). The dispersal distance of this species is several hundred meters to several kilometers per day (Touré et al. 1998a), so any significant barriers to dispersal should be recognizable from the 1 km resolution of this figure. Because no such barriers are evident from the figure, we infer that there were no significant impediments to dispersal in these populations.
|
Female mosquitoes were captured daily in the morning using mouth aspirators. They were temporarily stored in small, netted cardboard containers that were held in the shade and kept humidified until they were transported to the laboratory at the Malaria Research and Training Center in Bamako.
An. gambiae s.s. populations at Tene and Mopti are known to contain only the Mopti form (Touré et al. 1998b), whereas populations at the other sites are typically mixed, and therefore each specimen collected from these sites required cytological examination to determine the form to which it belonged. Chromosome preparations were made by extracting ovaries from the abdomen of each sample, following established protocols (Coluzzi 1968; Hunt 1973). Individuals from Mopti and Tene were identified to species by diagnostic polymerase chain reaction (PCR) (Scott et al. 1993). For the remaining locations, species identification and karyotype scoring was done using a phase-contrast microscope and classified according to Coluzzi et al. (1979). The rare individuals that were classified as hybrids by karyotype were excluded from this analysis. Following cytological examination, the carcass of each specimen was placed in a correspondingly numbered 1.5 ml tube containing silica gel for DNA extraction. Individual female mosquitoes were removed from the tubes in which they had been stored and DNA was extracted from the carcass using a standard extraction protocol (Post et al. 1993). The microsatellite primer sequences and PCR conditions employed were as described by Zheng et al. (1996). Loci were selected because previous studies had found them to be highly polymorphic in An. gambiae populations from Mali (Lanzaro et al. 1995, 1998). Microsatellites were amplified by PCR using fluorescent-labeled primers. The PCR products were fractionated on high-resolution polyacrylamide gels and scanned using an ABI-377 automated sequencer and GeneScan software (Applied Biosystems). Genotyping of individual mosquitoes was then performed using the Genotyper software package (Applied Biosystems).
A total of 11 populations, defined by chromosomal form at a given location, were compared. The number of microsatellite loci used in each comparison ranged from 5 to 23 and the number of individuals in each population ranged from 4 to 190 (Table 1). The loci used (by form and location) were as follows: Savanna form, Pimperena-Banamabani, AGXH7, AGXH8, AGXH25, AG2H26, AGXH38, AG2H79, AG3H83, AG3H88, AG2H95, AG3H119, AG2H125, AG3H128, AG2H143, AG2H156, AG3H158, AG2H175, AGXH293, AG3H577, AG2H603, AG2H637; Savanna form, Selinkenyi-Soulouba, AGXH7, AGXH36, AGXH53, AG2H85, AG3H93, AG2H117, AG2H135, AGXH145, AG3H158, AG2H187, AG2H197, AG2H290, AG2H417, AGXH503, AG3H525, AG3H544, AG3H555, AG2H603, AGXH678, AGXH711, AG3H750, AG3H776, AG2H786, AG3H817; Bamako form, Selinkenyi-Banambani, AGXH7, AGXH8, AGXH25, AG2H26, AG2H79, AG3H88, AG2H95, AG3H119, AG2H125, AG3H128, AG2H143, AG2H156, AG3H158, AG2H175, AGXH293, AG3H577, AG2H603, AG2H637; Bamako form, Selinkenyi-Soulouba, AGXH7, AGXH36, AGXH53, AG2H85, AG3H93, AG2H117, AG2H135, AGXH145, AG3H158, AG2H187, AG2H197, AG2H290, AG2H417, AGXH503, AG3H525, AG3H544, AG3H555, AG2H603, AGXH678, AGXH711, AG3H750, AG3H776, AG2H786, AG3H817; Mopti form, Selinkenyi-Banambani, AGXH7, AGXH8, AGXH25, AG2H26, AG2H79, AG3H88, AG2H95, AG3H119, AG2H125, AG3H128, AG2H143, AG2H156, AG3H158, AG2H175, AGXH293, AG3H577, AG2H603, AG2H637; Mopti form, Tene-Mopti, AGXH36, AGXH53, AG3H93, AG2H143, AGXH145, AG3H577, AGXH678, AGXH711, AG3H750; Mopti form, Tene-Selinkenyi, AGXH36, AGXH53, AG3H93, AGXH145, AGXH678, AGXH711, AG3H750; and Mopti form, Mopti-Selinkenyi, AGXH36, AGXH53, AG3H93, AGXH145, AGXH678, AGXH711, AG3H750.
|
The distribution of forms varies from location to location and estimates of gene flow would be confounded by analyses between both location and form. To factor out this effect, all estimates of genetic distance were made within a single chromosomal form. Unfortunately not all loci were the same for each sample, so a complete genetic distance matrix between all populations was impossible to obtain. In total we could make eight pairwise comparisons between the six locations listed above within the three chromosomal forms sampled.
Data were analyzed using the Arlequin software package developed by Schneider et al. (2000). Pairwise genetic differences between populations were determined from weighted estimates of FST and RST over all loci, using the number of different alleles or the sum of squared size difference, respectively. As expected, FST had smaller confidence limits than RST, so we used the former as the measure of differentiation. Values of Nem derived from FST were used to correlate with distance, where Ne is the effective population size and m is the fraction of the population replaced by immigrants. In some models of gene flow, Nem can be interpreted as an unbiased estimate of the number of individuals moving between populations, and when correlated with geographic distance, can be used to test isolation by distance among populations. The confidence limits of FST estimates were determined by a permutation approach where, p is equal to the number of permutations leading to an equal or larger value of FST. We assumed that the populations are more or less in equilibrium and that mutation rate is negligible compared to migration, thus calculating the number of migrants (Nem) from estimates of FST using the relation FST = 1/(Nem + 1) (Crow and Kimura 1970).
Strict applicability of this estimate of Nem depends on many assumptions, especially that the populations are at equilibrium with respect to drift and migration. Additional assumptions or shortcomings have been discussed by Goldstein and Pollack (1997), Whitlock and McCauley (1999), and in this particular context by Taylor et al. (2001). While this estimate for Nem derived from FST may not be strictly interpretable as the number of migrants, there is a monotonic relation between the two measures, so it still offers a heuristically useful way to compare levels of gene flow. So long as the assumptions are borne in mind, Nem derived from FST should still be appropriate for relative comparison of genetic differentiation as it is employed in this study.
| Results |
|---|
|
|
|---|
The loci in this study were highly polymorphic, with an average of 15.5 alleles per locus and an average gene diversity of 0.74. Significance tests of departures from HardyWeinberg expectations were performed for each locus using the method described by Guo and Thomson (1992). This method is equivalent to Fisher's exact test using a contingency table of arbitrary size and computing the probability of observing a table with the same marginal counts under the null hypothesis of no association. There were few significant deviations (P > .001) from HardyWeinberg expectations at each of the loci: 5 of the 13 populations sampled had no significant deviations at any of the loci, and all but 1 of the remaining 8 populations had no more than two loci that deviated. We attribute this to the presence of null alleles, as described in Lanzaro et al. (1995, 1998). The one exception, Pimperena-Savanna, had a total of nine loci that showed a significant departure from HardyWeinberg. The reason for this exception is unclear. The majority of the loci away from HardyWeinberg equilibrium were located on chromosome 2, which is known to be under strong selection (Touré et al. 1998b), but some loci on both chromosomes 3 and X were also affected.
Some of the populations, such as Soulouba-Bamako, had small sample sizes. They were included in this study only if they were offset by a large number of loci sampled (more than 20). Conversely, samples with small numbers of usable loci were used only if they were offset by a larger number of individuals.
The pairwise estimates of FST and Nem are shown with the geographic distance between each location in Table 1. Estimates of log(Nem) are plotted against geographic distance for all loci together and for the loci on each chromosome separately (Figure 2). The correlation between genetic and geographic distance was determined using least-squares regression. The significance test here is one-tailed because we do not expect a positive correlation between Nem and geographic distance. For all loci, there was a significant (P < .005) inverse correlation with r2 = 0.77 (Figure 2A). Parallel analyses of loci from each of the three individual chromosomes gave similar results. Correlations on both the third (r2 = 0.55) and X (r2 = 0.68) chromosomes were statistically significant (P < .05); correlations on the second chromosome were not statistically significant (P > .2), but were of similar magnitude to the others.
|
| Discussion |
|---|
|
|
|---|
These results demonstrate that different levels of genetic differentiation exist between populations of An. gambiae in West Africa. A strong correlation exists between genetic and geographic distance. This correlation is significant for all microsatellite loci across chromosomes and individually on two out of three chromosomes. Nonsignificant correlation for loci on the second chromosome may be due to the small number of comparisons or because of selection for inversions, documented by Touré et al. (1998b). All else being equal, selection would be confounded with the effects of isolation by distance, giving an overall lower correlation between differentiation and distance. Multilocus values of Nem ranged between 1.21 and 64.43, with the number of migrants between populations closely associated with their distance apart.
Recent mark-release-recapture (MRR) experiments have been used to estimate gene flow among populations in Mali by determining the amount of intervillage movement (Taylor et al. 2001). For the period from 1993 to 1998, movement between Banambani and its nearest neighboring village, Donéguébougou, was 0.039. The distance between the two neighboring villages is about 7 km. If we assume that the intervillage rate is analogous to the interpopulation rate, then the estimates for Nem were 35.9, 59.4, 75.0 for the Bamako, Savanna, and Mopti forms, respectively, quite in line with what was observed in this study.
These results contrast with several recent studies in both East and West Africa. Lehmann et al. (1996) found only small differences in microsatellite allele frequencies between populations of An. gambiae in Kenya and Senegal, a distance of approximately 6000 km. The average value of Nem was 2.5 using RST and 7.7 using FST for microsatellite loci. These results differ markedly from the results of this study, though the distances in our study were not so large. One possible reason the authors suggest for this low differentiation is range constraints on the evolution of the loci. If the number of repeats is bounded then, according to the stepwise mutation model, estimates of genetic differentiation using FST and RST in this way may be problematic for populations separated for long periods of time (Feldman et al. 1997; Goldstein and Pollack 1997; Nauta and Weissing 1996). So one possible explanation for the discrepancy is that the populations in Mali have not been separated by a large enough distance or long enough time for the range constraints to have an effect.
In East Africa, only the Rift Valley complex appears to be a barrier to gene flow for both microsatellite and mtDNA loci. Lehmann et al. (1997, 2000) found differences across the Rift Valley complex, while differences among populations on either side were much lower. Microsatellite values for FST and RST across the Rift Valley complex were 0.104 and 0.032, respectively, with an average distance between sites of approximately 700 km. On either side of the Rift Valley complex, both values of FST and RST were less than 0.01 with geographic distances up to 220 km. In our study, similar values (0.167, 0.140, and 0.034) of FST to those found across the Rift Valley complex were found for populations separated by 200500 km with no obvious barriers. Our results do not discount the notion that barriers exist. However, given that we observed such a strong correlation between geographic and genetic distance across uniformly favorable habitat, it is evident that barriers are not the only cause for a reduction in gene flow in this species. The situation in Mali, West Africa, is perhaps markedly different than that in East Africa.
| Acknowledgments |
|---|
This article was supported by a grant from the National Institutes of Health.
| Footnotes |
|---|
Corresponding Editor: Ross MacIntyre
Received November 11, 2000
Accepted November 26, 2001
| References |
|---|
|
|
|---|
-
Besansky NJ, Lehmann T, Fahey TG, Fontenille D, Braack LE, Hawley WA, and Collins FH, 1997. Patterns of mitochondrial variation within and between African malaria vectors, Anopheles gambiae and An. arabiensis, suggest extensive gene flow. Genetics 147:18171828.[Abstract]
Bryan JH, Di Deco MA, Petrarca V, and Coluzzi M, 1982 Inversion polymorphism and incipient speciation in Anopheles gambiae s. str. in the Gambia, West Africa. Genetica 59:167176.[CrossRef]
Coluzzi M, 1968. Cromosomi politenici delle cellule nutrici ovariche nel complesso gambiae del genere Anopheles. Parassitologia 10:179185.
Coluzzi M, Sabatini A, Petrarca V, and DiDeco MA, 1979. Chromosomal differentiation and adaptation to human environments in the Anopheles gambiae complex. Trans R Soc Trop Med Hyg 73:483497.[CrossRef][ISI][Medline]
Crow JF, Kimura M, 1970. An introduction to population genetics theory. New York: Harper & Row.
Feldman MW, Bergman A, Pollock DD, and Goldstein DB, 1997. Microsatellite genetic distances with range constraints: analytic description and problems of estimation. Genetics 145:207216.[Abstract]
Goldstein DB and Pollock DD, 1997. Launching microsatellites: a review of mutation processes and methods of phylogenetic inference. J Hered 88:335342.
Guo S and Thompson E, 1992. Performing the exact test of HardyWeinberg proportion for multiple alleles. Biometrics 48:361372.[CrossRef][ISI][Medline]
Hansen M, DeFries R, Townshend JRG, and Sohlberg R, 2000. Global land cover classification at 1 km resolution using a decision tree classifier. Int J Remote Sens 21:13311365.[CrossRef]
Hunt RH, 1973. A cytological technique for the study of Anopheles gambiae complex. Parassitologia 15:137139.[Medline]
Lanzaro GC, Touré YT, Carnahan J, Zheng L, Dolo G, Traoré S, Petrarca V, Vernick KD, and Taylor CE, 1998. Complexities in the genetic structure of Anopheles gambiae populations in West Africa as revealed by microsatellite DNA analysis. Proc Natl Acad Sci USA 95:1426014265.
Lanzaro GC, Zheng L, Touré YT, Traore SF, Kafatos FC, and Vernick KD, 1995. Microsatellite DNA and isozyme variability in a West African population of Anopheles gambiae. Insect Mol Biol 4:105112.[ISI][Medline]
Lehmann T, Besansky NJ, Hawley WA, Kamau L, Fahey TG, and Collins FH, 1997. Microgeographic structure of Anopheles gambiae in western Kenya based on mtDNA and microsatellite loci. Mol Ecol 6:243253.[CrossRef][Medline]
Lehmann T, Blackston CR, Besansky NJ, Escalante AA, Collins FH, and Hawley WA, 2000. The Rift Valley complex as a barrier to gene flow for Anopheles gambiae in Kenya: the mtDNA perspective. J Hered 91:165168.
Lehmann T, Hawley WA, Kamau L, Fontenille D, Simard F, and Collins FH, 1996. Genetic differentiation of Anopheles gambiae populations from East and West Africa: comparison of microsatellite and allozyme loci. Heredity 77:192208.
Massumi PP, Diuk-Wasser M, Touré YT, and Taylor CE, 2000. Canonical correlations between gene arrangement frequencies of Anopheles gambiae s.s. and weather variables in Mali, West Africa. Presented at the 49th Annual Meeting of the American Society of Tropical Medicine and Hygiene, Houston, TX, October 29November 2, 2000.
Nauta MJ and Weissing FJ, 1996. Constraints on allele size at microsatellite loci: implications for genetic differentiation. Genetics 143:10211032.[Abstract]
Post RJ, Flook PK, and Millest AL, 1993. Methods for the preservation if insects for DNA studies. Biochem Syst Ecol 21:8599.[CrossRef]
Schneider S, Roessli D, and Excoffier L, 2000. Arlequin: a software for population genetics data analysis, version 2.000. Geneva, Switzerland: Genetics and Biometry Lab, University of Geneva.
Scott JA, Brogdon WG, and Collins FH, 1993. Identification of single specimens of the Anopheles gambiae complex by the polymerase chain reaction. Am J Trop Med Hyg 49:520529.
Taylor CE, Touré YT, Carnahan J, Norris DE, Dolo G, Traoré SF, Edillo FE, and Lanzaro GC, 2001. Gene flow among populations of the malaria vector, Anopheles gambiae, in Mali, West Africa. Genetics 157:743750.
Thomson MC, Connor SJ, Milligan PJM, and Flasse SP, 1997. Mapping malaria risk in Africa: what can satellite data contribute? Parasitol Today 13:313318.[CrossRef][ISI][Medline]
Touré YT, Dolo G, Petrarca V, Traoré SF, Bouré M, Dao M, Carnahan J, and Taylor CE, 1998a. Mark-release-recapture experiments with Anopheles gambiae s.l. in Banambani, Mali, to determine population size and structure. Med Vet Entomol 12:7483.[CrossRef][ISI][Medline]
Touré YT, Petrarca V, Traoré SF, Coulibaly A, Maiga HM, Sankaré O, Sow M, Di Deco MA, and Coluzzi M, 1998b. The distribution and inversion polymorphism of chromosomally recognized taxa of Anopheles gambiae complex in Mali, West Africa. Parassitologia 40:477511.[Medline]
Wang R, Kafatos FC, and Zheng L, 1999. Microsatellite markers and genotyping procedures for Anopheles gambiae. Parasitol Today 15:3337.[CrossRef][ISI][Medline]
White F, 1983. The vegetation of Africa. A descriptive memoir to accompany the UNESCO/AETFAT/INSO vegetation map of Africa. Geneva, Switzerland: UNESCO.
Whitlock MC and McCauley DE, 1999. Indirect measures of gene flow and migration: FST
1/(4Nm + 1). Heredity 82:117125.
This article has been cited by other articles:
![]() |
N. C. Manoukis, J. R. Powell, M. B. Toure, A. Sacko, F. E. Edillo, M. B. Coulibaly, S. F. Traore, C. E. Taylor, and N. J. Besansky A test of the chromosomal theory of ecotypic speciation in Anopheles gambiae PNAS, February 26, 2008; 105(8): 2940 - 2945. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. J. Kent, S. Mharakurwa, and D. E. Norris Spatial and Temporal Genetic Structure of Anopheles arabiensis in Southern Zambia over Consecutive Wet and Drought Years Am J Trop Med Hyg, August 1, 2007; 77(2): 316 - 323. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. W. Chambers, J. K. Meece, J. A. McGowan, D. D. Lovin, R. R. Hemme, D. D. Chadee, K. McAbee, S. E. Brown, D. L. Knudson, and D. W. Severson Microsatellite Isolation and Linkage Group Identification in the Yellow Fever Mosquito Aedes aegypti J. Hered., May 1, 2007; 98(3): 202 - 210. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. E. CONN, J. H. VINEIS, J. P. BOLLBACK, D. Y. ONYABE, R. C. WILKERSON, and M. M. POVOA POPULATION STRUCTURE OF THE MALARIA VECTOR ANOPHELES DARLINGI IN A MALARIA-ENDEMIC REGION OF EASTERN AMAZONIAN BRAZIL. Am J Trop Med Hyg, May 1, 2006; 74(5): 798 - 806. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. MOLINA-CRUZ, A. M. P. DE MERIDA, K. MILLS, F. RODRIGUEZ, C. SCHOUA, M. M. YURRITA, E. MOLINA, M. PALMIERI, and W. C. BLACK IV GENE FLOW AMONG ANOPHELES ALBIMANUS POPULATIONS IN CENTRAL AMERICA, SOUTH AMERICA, AND THE CARIBBEAN ASSESSED BY MICROSATELLITES AND MITOCHONDRIAL DNA Am J Trop Med Hyg, September 1, 2004; 71(3): 350 - 359. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||




