Print version ISSN 1415-4757
Genet. Mol. Biol. vol.22 n.3 São Paulo Sept. 1999
Intrapopulational genetic differentiation in Anopheles (N.) darlingi Root, 1926 (Diptera: Culicidae) in the amazon region
Joselita Maria Mendes dos Santos1, Jorge Arturo Lobo2, Wanderli Pedro Tadei1 and Eucleia P. Betioli Contel3
1Coordenação de Pesquisas em Ciências da Saúde, INPA, Alameda Cosme Ferreira, 1756, 69011-970 Manaus, AM, Brasil. Send correspondence to J.M.M.S. 2Escuela de Biologia, Universidad da Costa Rica, San José, Costa Rica. 3Departamento de Genética, Faculdade de Medicina de Ribeirão Preto, 14049-900 Ribeirão Preto, SP, Brasil.
Four natural Amazonian Anopheles darlingi populations were analyzed for 19 loci, 12 of which showed variation, i.e., 68.4% of polymorphic loci. The Ariquemes, Rondônia, population was the most variable, with a large number of alleles per locus (2.26 ± 0.27) and high polymorphism (P = 63.15). The highest values of observed and expected intralocus heterozygosity were observed in the Manaus, Amazonas, population (Ho = 0.432 ± 0.11; He = 0.375 ± 0.08), and the lowest in the Cachoeira Porteira, Pará, population (Ho = 0.236 ± 0.09; He = 0.290 ± 0.11). Wright's F-statistic revealed disequilibrium caused by an excess of homozygotes, as shown by the Fis > Fst values (Fis = 0.083 > 0.026) reflecting intrapopulational differentiation. The four populations studied were genetically similar, as indicated by distance and similarity values. Chromosomal variation of Amazon population also did not indicate geographical differentiation, and populations in the central region of the Amazon Basin showed high polymorphism in relation to the marginal populations, which were mainly monomorphic.
An increasing number of studies on biological and genetical variability of Anopheles darlingi have been published in the last few years due to the importance of this mosquito as the major vector of human malaria in Brazil, especially in the Amazon region. These studies have demonstrated plasticity in response to environmental changes, which might affect the transmission of this disease (Santos et al., 1981, 1985; Tadei and Santos, 1982; Tadei et al., 1982; Rosa-Freitas et al., 1992, 1995; Harbach et al., 1993).
Chromosome polymorphism studies have revealed significant genetic differentiation between populations of A. darlingi from northern and southern Brazil. High frequencies of inversion in heterozygotes were observed in populations from the Amazon region (Kreutzer et al., 1972; Tadei and Santos, 1982; Tadei et al., 1982).
Brazilian populations of A. darlingi revealed specific differences in hydrocarbon profiles from Ilha de Maracá, State of Roraima, and Jaru, State of Rondônia (Walsh, 1989, in Rosa-Freitas et al., 1995). Distinctive hydrocarbon patterns in this species were also found by Rosa-Freitas et al. (1992) in Dourado (São Paulo) and Costa Marques (Rondônia), although specimens from Juturnaíba were often misassigned. However, these authors did not find significant behavior or allozyme differences in these populations.
Mitochondrial DNA and morphological studies were carried out for six Brazilian populations of A. darlingi (Juturnaíba, Dourado, Costa Marques, Cuiabá, Manaus, and Capanema), but none showed morphological differences, and nucleotide divergence estimates were not high (Rosa-Freitas et al., 1995).
In the present study, four populations of A. darlingi were analyzed using enzymatic variation of 19 loci to understand the genetic structure of this species, which is the principal vector of human malaria in the Brazilian Amazon.
MATERIAL AND METHODS
Mosquitoes were obtained from natural populations from four locations in the Amazon region: Ariquemes (10o S, 63o W) and Porto Velho/Samuel (9o S, 63o W), State of Rondônia; Manaus (3o S, 60o W), State of Amazonas, and Cachoeira Porteira (1o S, 57o W), State of Pará (Figure 1).
Figure 1 - Collection sites of Anopheles darlingi populations in the Amazon region
Wild females laid eggs individually in the laboratory, and from each single brood an average of five individuals were used for each single enzymatic system (Table I). Electrophoretic techniques and enzyme recipes were those described in Steiner and Joslyn (1979) and Lima and Contel (1990). The gels were prepared as described by Santos et al. (1996). The following loci were analyzed: EST1, EST2, EST5, LAP1, LAP2, LAP5, IDH1, IDH2, ODH1, ODH2, ODH3, AO1, AO2, XDH1, GDH, G-6PDH, 6-PGDH, a-GPDH and PGM.
Allele frequencies were estimated directly from the data. In the case of null alleles, their frequencies were estimated by maximum likelihood, using the dominance model of inheritance. Proportion of polymorphic loci (P), heterozigosities and Wright's F-coefficients were estimated using the Biosys Program (Swofford and Selander, 1981). The fixation index within populations (Fis) was also calculated using ALPHA software of the Genioc Program collection (Cabello and Krieger, 1997), with option 4 (calculation of F between systems). A dendrogram was constructed using the UPGMA method (Nei, 1978).
Allelic frequencies and polymorphism
Twelve of the 19 loci studied presented genetic variability in the four populations: EST1, EST2, EST5, LAP5, IDH1, ODH1, ODH2, ODH3, AO1, AO2, a-PGDH, and PGM. LAP1, LAP2, IDH2, XDH1, GDH, G-6PDH and 6-PGDH were monomorphic in all populations. The locus 6-PGDH showed variation only in the Ariquemes population.
Table I shows the number of alleles per locus in each of the four populations. The Ariquemes population was the most polymorphic and presented the largest number of alleles per locus. The Manaus and Cachoeira Porteira populations showed the lowest degree of polymorphism and the lowest number of alleles per locus.
Most of the loci analyzed in the Ariquemes and Manaus populations were not in Hardy-Weinberg equilibrium (Table II). Porto Velho/Samuel and Cachoeira Porteira population loci were found to be in equilibrium, except for EST2 in the latter and IDH1 in both. However, EST2 homozygotes from the same population were more frequent than expected (EST2*F/EST2*F and EST2*S/EST2*S genotypes). The same was observed in Ariquemes, where homozygote classes for the same locus (EST2) presented significantly higher values than heterozygotes, and this pattern (excess of homozygotes) was not observed for all loci for which chi-square values were significant. In the Manaus population, the heterozygote class was the most frequent for the EST2 locus.
Estimate of heterozygosity and F coefficient of Wright
Levels of observed and expected intralocus heterozygosity varied within each population (Table III). The Manaus population showed the highest level of heterozygosity. This result contrasts with those obtained for number of alleles per locus and percentage of polymorphic loci, for which Manaus presented the lowest variability. The Porto Velho/Samuel population presented a wide range of observed heterozygosity and an amplitude of expected heterozygosity. Most of the loci analyzed also had excess homozygotes, with highly significant F values at the 1% level for loci EST2.
Despite the excess homozygotes observed in most of the loci analyzed, we also detected excess heterozygotes for locus EST5 in Ariquemes, ODH1 in Porto Velho/Samuel and EST2, IDH1 and PGM in Manaus. Of these, only EST2 in Manaus and EST5 in Ariquemes showed statistically significant (P < 0.05) excess heterozygosity.
The total F calculated by maximum likelihood (Table III) shows that there was significant disequilibrium in favor of excess homozygosity in the Ariquemes, Porto Velho/Samuel and Cachoeira Porteira populations. The mean F for Manaus was 0.01 ± 0.03, which does not differ significantly from zero. However, the standard error for the F estimates in the other populations ranged from 0.02 to 0.03, being much lower than the mean estimated values. The F coefficient was higher in Cachoeira Porteira and lower in Porto Velho/Samuel, possibly due to the influence of EST2 locus which did not show disequilibrium in Porto Velho/Samuel.
However, as can be seen from the values for each locus, mean F was significantly different between systems of the Ariquemes and Cachoeira Porteira populations, precisely those with the highest F values. This result was obtained from the chi-square values for heterogeneity (Table III).
Population genetic structure
Fis and Fst values varied considerably from locus to locus. Fis values were greater than Fst values for almost all loci, with the exception of EST5 and IDH1 (Table IV). Thus, we observed that gene frequency variability had an intrapopulational origin, and allele origin did not have much of an effect on the probability of sampling two different alleles in the set of populations.
Genetic distance and similarity data indicate that the four populations studied are genetically similar (Table V). The Ariquemes population is the closest to the Manaus population. Nevertheless, Ariquemes, Cachoeira Porteira and Porto Velho/Samuel have more in common genetically with each other than with Manaus (Figure II). The Manaus population is alone in a separate cluster.
Figure 2 - Dendrogram resulting from grouping populations on the basis of genetic similarity by the unweighted pairing group method with arithmetic mean (UPGMA), determined by Nei (1978). A, Ariquemes; B, Cachoeira Porteira; C, Porto Velho/Samuel, and D, Manaus.
DISCUSSION AND CONCLUSIONS
Levels of genetic variability
Allele frequency was not uniform for the different loci. Of the 19 loci analyzed, six presented variation. The observed deviation from genetic equilibrium was due in part to a larger number of observed than expected homozygotes, except for locus EST2 in the Manaus population.
Crouau-Roy (1988) pointed out that many studies have shown significant levels of heterozygote deficiency in natural populations, but usually these deficiencies are observed at some allozymic loci and/or only in some populations which have other loci in equilibrium. It has also been observed that different loci do not present the same heterozygote deficiency. These observations have been frequently made in mosquitoes (Santos et al., 1985; Van Driel et al., 1987; Hii et al., 1991) and other insects (Sturgeon and Mitton, 1986: Crouau-Roy, 1988). These deficiencies were attributed by these investigators to differential genotype selection, inbreeding, null alleles, founder effect (bottleneck) and even the population structure itself. The same results were observed by Scarpassa (1988) in a natural population of A. nuneztovari from Tucuruí for three nonspecific esterase loci. The author attributed these results to ecological alterations introduced into the area which probably affected the genotype frequencies, with different adaptive values for homozygote classes.
In the case of null alleles for locus ODH1 of Porto Velho/Samuel population, it can be seen that the locus is in Hardy-Weinberg equilibrium. The occurrence of null alleles was observed in loci A and B of EST-C of A. punctipennis (Narang and Kitzmiller, 1971) and in the aldehyde oxidase of A. nuneztovari (Narang et al., 1979).
On the basis of the data presented, though the alleles had different frequencies among populations, the loci analyzed showed a certain homogeneity, with no detection of diagnostic loci for the various populations.
The proportion of polymorphic loci in the 19 loci analyzed in A. darlingi indicates a high level of polymorphism even though the number of alleles per locus showed small variation. Narang et al. (1979) also found a high level of polymorphism (P = 0.632) for a population of A. darlingi from BR-174 (Manaus-Boa Vista). A relatively high degree of polymorphism has been also described in A. quadrimaculatus by Narang et al. (1989) in 50% of the 33 loci analyzed. However, Van Driel et al. (1987) found in strains of A. stephensi from Pakistan that only six loci were polymorphic of the 16 loci analyzed. This species also showed lower observed heterozygosity as compared with A. darlingi (Ho = 0.065 to 0.086/0.236 to 0.432). For 21 populations of different and incipient species of Anopheles, Cianchi et al. (1985) detected a wider range of both polymorphism and expected heterozygosity (0.31 £ P £ 0.78; 0.09 £ He £ 0.24). The large scope of genetic variability in Anopheles species has been closely related to the very structure of the population, such as differential mating, mutation, selection, and genetic drift.
Considering the values observed for the A. darlingi population studied, there was high heterozygosity with levels close to those for the genus Culex (0.290 £ He £ 0.375). The same is true for polymorphism (0.526 £ P £ 0.631). H and P values obtained in the present study contrast with those reported by Lanzaro et al. (1990) for the A. quadrimaculatus complex in which heterozygosity was much lower (0.10 £ H £ 0.18). However, polymorphism showed a greater amplitude of variation (0.32 £ P £ 0.60).
Genetic structure and population differentiation
The disequilibrium resulting from excess homozygotes was observed in strains of A. stephensi (Van Driel et al., 1987), in which Fis values were higher than Fst values for most of the loci analyzed. Mean fixation among the different loci was low (Fit = 0.107). However, locus Gpdh of A. stephensi and locus AO1 of A. darlingi had high fixation indexes (Fit = 0.831 and 0.663, respectively). The loci, especially Gpdh and A01, did not make an equal contribution in either species. For A. darlingi, the fixation index indicated the occurrence of some type of isolation in most of the loci analyzed, which interferred with random mating within each population.
Based on fixation index data, variability of allele frequencies in A. darlingi populations is judged to be of intrapopulation origin, resulting from differentiation within these populations. This differentiation was produced by excess homozygotes within each locus and was more relevant within than between populations. However, elevated genetic divergence between populations of C species of A. quadrimaculatus (Narang et al., 1990) was observed, as indicated by the Fst value of 0.219. Genetic divergence between populations can be observed when a certain degree of reproductive isolation has already become established in the populations, as described by Narang et al. (1990).
The genetic distance and similarity between the four A. darlingi populations indicate that these populations are genetically similar (0.976 £ I £ 0.990) although they are geographically isolated. The genetic distance of these populations was very small (0.010 £ D £ 0.024), with Ariquemes showing the smallest genetic distance (0.010) compared to Cachoeira Porteira. However, this difference was not significant, so that Ariquemes, Cachoeira Porteira and Porto Velho/Samuel were grouped into the same cluster of the dendrogram.
Comparing these data with those for other mosquito species shows that genetic similarity and distance are similar, except in interspecific cases. For example, for species of the A. marshalli complex, isoenzyme analysis showed little divergence, with a genetic similarity index of 0.803-0.972 (Lambert, 1983). These values are similar to some pairs of cryptic species in which the genetic distances were D = 0.13 between A. gambiae Giles and A. arabiensis Patton; D = 0.23 between A. melanoon Hackett and A. maculipennis Meigen, and D = 0.17 between species A and B of A. stephensi strains (Van Driel et al., 1987). These results are similar to the values obtained for the four populations of A. darlingi, although with a lower amplitude of variation (0.010-0.024). These data may be interpreted on the basis of the genetic structure of the populations themselves, in which matings may occur with preferential crosses. Thus, considering the present data as a whole, we note that the populations presented greater intrapopulational than interpopulational differentiation in a process that has extended to four Amazonian populations studied. Datas on chromosomal variation in different samples from Amazonia did not indicate geographically differentiated populations. Samples from the central region of Amazon Basin did, however, present high polymorphism in relation to the marginal populations, which were mainly monomorphic (Tadei and Santos, 1982; Tadei et al., 1982; Tadei, 1986).
The authors are grateful to Juracy de Freitas Maia for technical assistance. This research was supported in part by CNPq/ELETRONORTE. Publication supported by FAPESP.
Amostras de populações naturais de Anopheles darlingi foram analisadas para 19 locos, dos quais 12 apresentaram variação, com uma proporção de locos polimórficos de 68,40%. A população de Ariquemes mostrou maior variabilidade, com maior número de alelos por loco (2,26 ± 0,08) e maior polimorfismo (P = 63,15). Levando-se em conta os níveis de heterozigosidade intraloco observados e esperados, os maiores valores foram observados na população de Manaus (Ho = 0,432 ± 0,11; He = 0,375 ± 0,08) e os menores na população de Cachoeira Porteira (Ho = 0,236 ± 0,09; He = 0,290 ± 0,11). Quando analisados pelas estatísticas F de Wright, os resultados mostraram desequilíbrio decorrente de excesso de homozigotos, conforme mostram os valores de Fis > Fst (0,083 > 0,026), denotando uma certa diferenciação intrapopulacional. As quatro populações estudadas foram semelhantes geneticamente, conforme indicam os valores de distância e similaridade. As variações cromossômicas de populações da Amazônia não indicaram diferenciação geográfica, e na região central da bacia amazônica, estas foram altamente polimórficas em relação às marginais, que foram predominantemente monomórficas.
Cabello, P.H. and Krieger, H. (1997). Genioc. Sistema para Análises de Dados de Genética. FIOCRUZ, Rio de Janeiro. [ Links ]
Cianchi, R., Urbanelli, S., Villani, F., Sabatini, A. and Bullini, L. (1985). Electrophoretic studies in mosquitoes. Recent advances. Parassitologia 27: 157-167. [ Links ]
Crouau-Roy, B. (1988). Genetic structure of cave-dwelling beetle populations: significant deficiencies of heterozygotes. Heredity 60: 321-327. [ Links ]
Harbach, R.E., Roberts, D.R. and Manguin, S. (1993). Variation in the hindtarsal markings of Anopheles darlingi (Diptera: Culicidae) in Belize. Mosq. Syst. 25: 192-197. [ Links ]
Hii, J.L.K., Chew, M., Sang, V.Y., Munstermann, L.E., Tan, S.G., Panyim, S. and Yasothornsrikul, S. (1991). Population genetic analysis of host seeking and resting behaviors in the malaria vector, Anopheles balabacensis (Diptera: Culicidae). J. Med. Entomol. 28: 675-684. [ Links ]
Kreutzer, R.D., Kitzmiller, J.B. and Ferreira, E. (1972). Inversion polymorphism in the salivary gland chromosomes of Anopheles darlingi. Mosq. News 32: 555-556. [ Links ]
Lambert, M.D. (1983). A population genetical study of the African mosquito Anopheles marshallii (Theoband). Evolution 37: 484-495. [ Links ]
Lanzaro, G.C., Narang, S.K. and Seawright, J.A. (1990). Speciation in an Anopheline (Diptera: Culicidae) mosquito: Enzyme polymorphism and the genetic structure of populations. Ann. Ent. Soc. Am. 83: 578-585. [ Links ]
Lima, L.M.K.S. and Contel, E.P.B. (1990). Electrophoretic analysis of 12 proteins in natural populations of Spodoptera frugiperda (Lepidoptera: Noctuidae). Rev. Bras. Genet. 13: 711-729. [ Links ]
Narang, G., and Kitzmiller, J.B. (1971). Esterase polymorphism in a natural population of A. punctipennis. II. Analysis of the Est-C system. Can. J. Genet. Cytol. 13: 771-776.
Narang, S., Santos, J.M.M., Garcia, J.C., Cristakou, H.D. and Narang, N. (1979). Genética de populações de anofelinos IV. Estudos eletroforéticos das populações naturais de Anopheles nuneztovari e Anopheles darlingi. Correlação genética entre espécies. Acta Amazonica 9: 529-542. [ Links ]
Narang, S.K., Kaiser, P.E. and Seawright, J.A. (1989). Dichotomus electrophoretic taxonomic key for identification of sibling species A, B and C of the Anopheles quadrimaculatus complex (Diptera: Culicidae). J. Med. Ent. 26: 94-99. [ Links ]
Narang, S.K., Seawright, J.A. and Kaiser, P.E. (1990). Evidence for microgeographic genetic subdivision of Anopheles quadrimaculatus species C. J. Am. Mosq. Control Assoc. 6: 179-187. [ Links ]
Nei, M. (1978). Estimation of average heterozygosity and genetic distance from a small number of individuals. Genetics 89: 583-590. [ Links ]
Rosa-Freitas, M.G., Broomfield, G., Priestman, A., Milligan, P.J.M., Momen, H. and Molyneux, D.H. (1992). Cuticular hydrocarbons, isoenzymes and behavior of three populations of Anopheles darlingi from Brazil. J. Am. Mosq. Control Assoc. 8: 357-366. [ Links ]
Rosa-Freitas, M.G., Conn, J., Mitchell, S.E., Cockburn, A.F., Seawright, J.A. and Momen, H. (1995). Mitochondrial DNA and morphological analyses of Anopheles darlingi populations from Brazil (Diptera: Culicidae). Mosq. Syst. 27: 78-99. [ Links ]
Santos, J.M.M., Contel, E.P.B. and Kerr, W.E. (1981). Biologia de Anofelinos amazônicos. 1. Ciclo biológico, postura e estádios larvais de Anopheles darlingi Root, 1926 (Diptera: Culicidae) da Rodovia Manaus/Boa Vista. Acta Amazonica 11: 789-797. [ Links ]
Santos, J.M.M., Contel, E.P.B. and Kerr, W.E. (1985). Biology of Amazonian mosquitoes. III. Esterase isozymes in Anopheles darlingi. Acta Amazonica 15: 167-177. [ Links ]
Santos, J.M.M., Tadei, W.P. and Contel, E.P.B. (1996). Electrophoretic analysis of 11 enzymes in natural populations of Anopheles (N.) darlingi Root, 1926 (Diptera: Culicidae) in the Amazon region. Acta Amazonica 26: 97-114. [ Links ]
Scarpassa, V.M. (1988). Estudo do ciclo biológico e de isoenzimas na ontogênese de Anopheles (Nyssorrhynchus) nuñez-tovari Gabaldon, 1940 (Diptera: Culicidae). Master¢s thesis, INPA/FUA, Manaus. [ Links ]
Steiner, W.W.M. and Joslyn, D.J. (1979). Electrophoretic techniques for the genetic studies of mosquitoes. Mosq. News 39: 35-54. [ Links ]
Sturgeon, K.B. and Mitton, J.B. (1986). Allozyme and morphological differentiation of mountain pine beetles Dendroctonus pondersae Hopkins (Coleoptera: Scolytidae) associated with host tree. Evolution 40: 290-302. [ Links ]
Swofford, D.L. and Selander, R.B. (1981). Biosys-1: A Fortran program for the comprehensive analysis of electrophoretic data in population genetics and systematics. J. Hered. 72: 281-283. [ Links ]
Tadei, W.P. (1986). Estudos em espécies de Anopheles. Desempenho e rumos da Genética no Brasil. Rev. Bras. Genét. IX (Anais): 221-223. [ Links ]
Tadei, W.P. and Santos, J.M.M. (1982). Biologia de anofelinos amazônicos. VII. Estudo da variação de freqüência das inversões cromossômicas de Anopheles darlingi Root (Diptera: Culicidae). Acta Amazonica 12: 759-785. [ Links ]
Tadei, W.P., Santos, J.M.M. and Rabbani, M.G. (1982). Biologia de anofelinos amazônicos. V. Polimorfismo cromossômico de Anopheles darlingi Root (Diptera: Culicidae). Acta Amazonica 12: 353-369. [ Links ]
Van Driel, J.W., Sluiters, J.F. and Van Der Kaay, H.J. (1987). Allozyme variation in Anopheles stephensi Liston from Pakistan (Diptera: Culicidae). Biochem. Genet. 25: 789-802. [ Links ]
(Received August 14, 1997)