Skip to main content

Seasonal genetic partitioning in the neotropical malaria vector, Anopheles darlingi

Abstract

Background

Anopheles darlingi is the main malaria mosquito vector in the Amazonia region. In spite of being considered a riverine, forest-dwelling species, this mosquito is becoming more abundant in peri-urban areas, increasing malaria risk. This has been associated with human-driven environmental changes such as deforestation.

Methods

Microsatellites were used to characterize A. darlingi from seven localities along the Madeira River, Rondônia (Brazil), collected in the early and late periods of the rainy season.

Results

Two genetically distinct subpopulations were detected: one (subpopulation A) was associated with the late rainfall period and seems to be ecologically closer to the typical forest A. darlingi; the other (subpopulation B) was associated with the early rainfall period and is probably more adapted to drier conditions by exploiting permanent anthropogenic breeding sites. Results suggest also a pattern of asymmetric introgression, with more subpopulation A alleles introgressed into subpopulation B. Both subpopulations (and admixed mosquitoes) presented similar malaria infection rates, highlighting the potential for perennial malaria transmission in the region.

Conclusions

The co-occurrence of two genetically distinct subpopulations of A. darlingi adapted to different periods of rainfall may promote a more perennial transmission of malaria throughout the year. These findings, in a context of strong environmental impact due to deforestation and dam construction, have serious implications for malaria epidemiology and control in the Amazonian region.

Background

Anopheles darlingi is a major malaria mosquito vector in the Americas. It is of particular importance in the Amazon basin where most malaria cases occur in South America[1, 2]. This species has a widespread geographic distribution (from Mexico to Argentina) is highly anthropophilic and is susceptible to human Plasmodium malaria parasites. Anopheles darlingi is typically described as a lowland, riverine, forest-dwelling species[3]. However, mosquito populations of this species display great behavioural plasticity and morphological, biological and genetic diversity[1, 36]. Moreover, this vector has been found at high densities in areas affected by deforestation and anthropogenic environmental changes[7].

Due to its intrinsic diversity the question of A. darlingi being a species complex has been raised by many authors. Charlwood[8] found chromosomal, morphological and behavioural differences between A. darlingi from central Amazon and further north (Venezuela and Guyana). However, several isozyme and DNA-based studies reported levels of genetic differentiation consistent with a monotypic species throughout the species' range[911]. Although mtDNA history revealed some genetic divergence between South and Central American populations[1214], complete mtDNA genomes of A. darlingi did not provide support for speciation within the taxon[6]. Nuclear microsatellite markers revealed restricted contemporary gene flow among A. darlingi populations over broad geographical ranges (continental scale) but mainly due to isolation by distance[1517].

Precipitation is considered to be the main ecological variable influencing A. darlingi distribution and abundance[3]. The annual density of A. darlingi populations is seasonal and dependent upon the availability of larval habitats[1]. In riverine areas, seasonality varies between less populated rural sites, with high mosquito densities after the peak of rainfall[18], and urban and suburban sites, where two density peaks may occur in the transition from the dry to the rainy season and vice versa[19, 20]. In drier inland areas (further from the influence of rivers) with human disturbance, a peak of mosquito density may occur in the dry season, associated with the presence of artificial permanent breeding sites, as dams and artificial lakes[18, 21].

Variable density peaks may be the result of seasonal fluctuations of the same population or the consequence of the presence of two differentially adapted subpopulations of the same species. In both cases, there is a growing concern that the increased anthropogenic influence in riverine, forested habitats of A. darlingi may result in changes in the seasonality patterns of this mosquito with implications for malaria transmission. Such is the case of the region of Porto Velho, in Rondônia state, Brazil. This city lies along the Madeira River, an important tributary of the Amazona River. A major environmental disturbance is expected in this region due to the construction of two hydroelectric power dams along the river, with a planned total flooded area of 380 sq km, which will likely impact local mosquito populations.

In the present study, microsatellite loci were used to analyse genetic variation of Amazonian A. darlingi from the region of Porto Velho (RO). Results suggest a season-dependent genetic partition, with two distinct mosquito subpopulations that display similar Plasmodium infection rates.

Methods

Mosquito collections

Porto Velho (ca. 436,000 inhabitants) is the capital city of the state of Rondônia (Brazil), which lies along the shores of Madeira River in the western part of the Amazonian forest, close to the border with Bolivia. Most of the human population lives on the right bank of the Madeira River, while the left bank is sparsely populated with small farms and settlements created by the National Institute of Colonization and Agrarian Reform. The state is one of the most affected by deforestation in Brazil[22]. The region has a tropical monsoon climate (Koppen classification, Am) with a dry season from May to September and a rainy season from October to April[23]. Anopheles darlingi abundance is seasonal and follows the rainfall regime. During the rainy season, human-biting mosquito densities reach a peak in March-April, with over 30 bites per human per hour, whereas in the dry season biting rates fall to nearly zero[19]. In urban and suburban areas of the region, an additional density peak is recorded between the end of the dry season and the beginning of the next rainy season (October-November) and malaria transmission continues throughout the year[19]. Mean annual rainfall and average temperature values between 1961 and 1990 were obtained for the region from Instituto Nacional de Meteorologia (INMET) (Figure 1A).

Figure 1
figure 1

Porto Velho mean annual rainfall and average temperature values between 1961 and 1990 (A) and Map of Brazil, showing de state of Rondonia and the collection sites along the Madeira River (B). Rainfall and temperature data was obtained from Instituto Nacional de Meteorologia (INMET). The Landsat image was obtained from http://glovis.usgs.gov/ (U. S. Geological Survey). 1: Jaci Parana; 2: Amazonas; 3: Teotonio; 4: Bate Estaca; 5: Santo Antonio; 6: Engenho Velho; 7: Vila Candelaria.

Mosquito sampling took place in seven localities along the Madeira River near Porto Velho (Figure 1B; Table 1): Vila Candelária, a riverside district of the city, less than 2 km from the city centre; the suburban riverine localities of Bate Estaca, Engenho Velho, Santo Antônio, Teotônio, and Amazonas, located 5 to 20 km southwest of the city centre; and, Jaci Paraná, a more inland locality located in a deforested area, ca. 75 km southwest of Porto Velho.

Table 1 Characterization of the sampling sites of A. darlingi in Madeira River, Rondônia state, Brazil (2007)

Collections took place between 2007 and 2008 during two different periods of the rainy season: February and March (hereafter termed late rainfall period) corresponding to the end of the rainy season and when A. darlingi abundance reaches its peak in the region[19]; October and December (hereafter termed early rainfall period) corresponding to the end of the dry season and the beginning of the rainy season, when mosquito abundance begins to increase. Mosquitoes were captured by outdoor landing catches with hand-held aspirators. Six 6-hour (18.00 to 00.00) collections were performed at each site.

Collected mosquitoes were morphologically identified to species using the identification keys of the subgenus of Nyssorhynchus[24]. Anopheles darlingi specimens were preserved individually in tubes with 70% isopropyl alcohol.

Microsatellites genotyping

DNA was extracted from each mosquito (whole body) using 5% Chelex solution (BioRad), according to the manufacturer's recommendations and eluted in 200 μl of double-distilled water.

Ten microsatellite loci were genotyped (see Additional file1). Eight of the primer pairs were originally described by Conn et al.[25]. Due to amplification problems, the primers of locus ADC107 were redesigned. Furthermore, primers for two additional loci (ADMP2 and ADMP9) were obtained from the GenBank sequences reported by Calvo et al.[26], both from 3’UTR of unknown transcripts. Each locus was amplified by PCR using fluorescently labelled (FAM, NED, or HEX; Applied Biosystems) reverse primers. PCR reactions were conducted in a 20-μl PCR mix containing 1× PCR GoTaq® Flexi Buffer (Promega), 1.5 mM MgCl2, 200 μM dNTPs, 0.5 μM of each primer, 0.5 U of Taq polymerase and 1 μl of DNA template. Cycling conditions included an initial denaturation at 94°C for 5 min, followed by 30 cycles at 94°C for 30 sec, annealing (52-57°C, Additional file2) for 30 sec, elongation at 72°C for 35 sec; and a final extension step of 5–10 min at 72°C. Amplified fragments were separated by capillary electrophoresis in an automatic sequencer (ABI3730, Applied Biosystems) and sizes scored using the software GeneMarker (SoftGenetics) and GeneMapper (Applied Biosystems).

Plasmodium detection

The detection of Plasmodium sp. DNA in individual A. darlingi mosquitoes was performed by real-time PCR using a universal primer (P1 forward primer [5’ -ACGATCAGATACCGTCGTAATCTT-3’]) and the species-specific nucleotide sequences of the 18S rRNA genes of Plasmodium falciparum (5’-CAATCTAAAAGTCACCTCGAAAGATG-3’) and Plasmodium vivax (5’-CAATCTAAGAATAAACTCCGAAGAGAAA-3’)[27]. Each 15-μl reaction mix contained 2 μl of sample DNA, 7.5 μl of 1× DNA SYBR Green reagent, and 0.5 μM concentrations of each primer. The PCR conditions consisted of 50°C for 2 min, followed by 95°C for 10 min and 40 cycles each with 15 sec at 95°C and 1 min at 60°C. Positive samples, when difference from blank Ct was greater than 15 cycles, were re-analysed for confirmation. All positive samples were confirmed.

Microsatellites data analysis

Estimates of expected heterozygosity (He) allele richness (Rs) and private allele richness (Ps) were calculated in FSTAT v 2.9.3.2[28] and HP-RARE[29]. Hardy-Weinberg (HW) equilibrium tests were performed with the software ARLEQUIN 3.1.1[30]. Linkage disequilibrium (LD) between pairs of loci was tested using the log likelihood ratio statistic available in GENEPOP[31, 32]. The software MICRO-CHECKER[33] was used to detect the presence of null alleles.

The genetic structure of the A. darlingi populations was assessed by computing locus-specific and pairwise F ST estimates[34] and by analysis of molecular variance (AMOVA)[35], both available in Arlequin. Permutation-based statistical tests were performed with 5,000 permutations. Bayesian clustering analysis with STRUCTURE 2.3.4[36] was also performed, without prior information on sampling groups, using the admixture model and assuming correlated allele frequencies among populations (λ was set at 1, default value). Twenty independent runs with 105 burn-in steps and 105 iterations were done for each value of K (K varied from 1 to 13). The number of genetic clusters (K) was inferred using the ΔK statistic[37] with the software STRUCTURE HARVESTER[38]. The information from the outputs of each K (20 runs) was compiled by the Greedy method implemented in CLUMPP software[39]. Individual assignment to clusters was performed using a probability threshold (Tq) of 0.80 following the recommendations of Vaha and Primmer[40].

The Bayesian method implemented in NewHybrids[41] was used to estimate the posterior probability of individuals falling into six genotypic classes: pure-bred, F1, F2 and backcrosses. The approach of uniform priors was used. Results were based on the average of five independent runs of 105 iterations. Individual assignment to classes was based on a posterior probability threshold, Tq > 0.5.

In order to detect recent population perturbations, Wilcoxon signed-ranks tests of heterozygosity were performed with BOTTLENECK 1.2.02[42]. Calculations were done under the stepwise mutation model (SMM) and a two-phase model (TPM) with 5% of indels larger than one repeat[43]. Analyses were based on 1,000 replications. Wilcoxon signed-ranks tests to compare means and Chi-square or exact Fisher’s tests on contingency tables were performed with IBM SPSS Statistics v. 22 for Windows (IBM Corp). Whenever multiple tests were performed, the nominal significance level (α = 0.05) was adjusted by the sequential Bonferroni procedure[44].

Results

A total of 440 A. darlingi were analysed in this study. Of these, 158 were collected in the early rainfall period whereas 282 were collected in the late rainfall period (Table 1).

Genetic diversity

All ten microsatellite loci were polymorphic and a total of 212 alleles were detected across all loci and samples. The number of alleles per locus ranged from five (ADSP2) to 44 (ADC01). Exact tests of LD showed one significant association for the pair ADC29-ADC110 in the sample of Amazonas, early rainfall period (P <0.001; adjusted significance level for 45 pairwise tests per population). Significant departures from H-W equilibrium were detected at loci ADC107n (all samples) and ADC29 (in 11 samples) suggesting the presence of null alleles (Table 2). These loci were thus excluded from subsequent analyses.

Table 2 Microsatellite variation for A. darlingi in Porto Velho (RO)

Estimates of expected heterozygosity (He) and allelic richness (Rs) per locality and sampling period are presented in Table 2. The mean Rs over loci ranged from 9 to 10. The mean He over loci ranged from 0.142 in Jaci Paraná to 0.838 in Santo Antônio (both early rainfall periods). The largest difference in He between sampling periods was found in the most inland sample of Jaci Paraná, with an estimate of 0.142 in the early rainfall and 0.794 in the late rainfall.

Genetic structure

Bayesian clustering analysis implemented in STRUCTURE revealed two genetically distinct subpopulations (best K = 2), hereafter denoted as subpopulations A and B (Figure 2). A total of 133 individuals showed admixed ancestry between the two subpopulations (i e, 0.2 < qi < 0.8). When these individuals were excluded and subpopulations A and B were analysed separately, each consisted of a single genetic cluster (best K = 1). There was no particular association between these subpopulations and any of the localities sampled (Figure 2A and Additional file3). However, each one was associated with the period of the collection. Subpopulation A grouped 173 individuals, 171 (99%) of which were collected in the late rainfall period. Subpopulation B grouped 134 individuals, 110 (82%) of which were collected in the early rainfall period and 24 (18%) in the late rainfall (Figure 2B). The seasonal pattern of alternation between subpopulations A and B from late to early rainfall periods was consistent in all riverine localities, with a predominance of a single subpopulation in both months of collection at each period. In the inland Jaci Paraná, however, subpopulation A predominated over B only in the first month of collection (February), whereas in April, B was already the most abundant (Figure 2A).

Figure 2
figure 2

Bayesian clustering analysis by STRUCTURE [[36]] of A. darlingi. The multilocus genotype of each individual is represented by a column with two colours showing the genetic proportion assigned to each cluster (K = 2). Blue: Lineage A; Red: Lineage B. Individuals were ordered by locality and sampling period (ER: early rain; LR: late rain).

The pattern of genetic structuring between sampling periods (i e, early and late rainfall) was confirmed by a significant estimate of global F ST over loci (0.023, P < 0.0001). When all individuals were included significant F ST was obtained in two (P2 and P9) out of the eight loci analysed (Table 3). The overall F ST estimate more than doubled when admixed individuals (i e, those with qi < 0.80 for both subpopulations, from STRUCTURE) were removed from the analysis (0.062, P <0.0001), and six out of eight loci showed a significant F ST (Table 3). Pairwise F ST estimates among sampling sites within subpopulations (without admixed individuals) ranged from 0 to 0.036 in subpopulation A (with three out of 21 significant comparisons, see Additional file2), while in subpopulation B F ST values varied between 0.002 and 0.032 and were all non-significant (Additional file2).

Table 3 Estimates of F ST per locus in Anopheles darlingi populations with and without admixed individuals

Furthermore, the component of variance between sampling periods was higher than that among sampling sites within each subpopulation when admixed individuals were excluded, showing 2.87% of genetic variation between sampling periods, 14 times higher than the non-significant percentage of variation among populations within subpopulation (Additional file4).

The analysis with NewHybrids showed a homogeneity of subpopulation A, associated with the late rainfall sampling period, with 96% of the individuals being assigned as purebred A and 4% as hybrids. In subpopulation B, associated with the early rainfall period, 46% individuals were classified as purebred B, 6% as purebred A and 48% were assigned to hybrid classes (Table 4). Of the four hybrid classes, only F2 hybrids had individuals with qi > 0.5. In subpopulation B 39% of individuals were classified as F2, while in subpopulation A 4% of individuals were assigned as F2 hybrids (Table 4).

Table 4 Frequencies of purebred and hybrid individuals, as detected by NewHybrids[41] in each of the lineages determined by STRUCTURE analysis

Subpopulation diversity and demographic inference

Table 5 presents estimates of genetic diversity for subpopulations A and B (excluding admixed individuals). There were no significant differences between subpopulations in the mean allele richness (Wilcoxon signed-ranks test: P = 0.575), mean private allele richness (Wilcoxon signed-ranks test: P = 0.575) or expected heterozygosity (Wilcoxon signed-ranks test: P = 0.735).

Table 5 Estimates of genetic diversity for lineages A and B (as determined by STRUCTURE) of Anopheles darlingi

Since only two individuals were assigned to subpopulation A in early rainfall sampling period, heterozygosity tests implemented in BOTTLENECK were performed on three samples: subpopulation A (late rainfall), subpopulation B (early rainfall) and subpopulation B (late rainfall). Wilcoxon tests revealed a significant deficit of heterozygotes indicative of recent population expansion in subpopulation A while all other tests were non-significant (Additional file5).

Mosquito infection rates

Of the total mosquito DNA samples analysed (N = 440), 16 (3.6%) were positive for Plasmodium sp. DNA. Infected mosquitoes were found in nine out of the 12 collections analysed. The collections with no positive mosquitoes were Bate Estaca (early rainfall), Vila Candelária (early rainfall) and Teotônio (late rainfall). Infection rates were comparable between the two sampling periods: 3.2% (five out of 158) in the early and 3.9% (11 out of 282) in the late rainfall (χ2 = 0.157; d.f. 1; P = 0.692). Plasmodium falciparum was the predominant species, occurring in 14 (87.5%) of the positive mosquitoes. There was also one mosquito (subpopulation A, late rainfall) infected with P. vivax and one P. falciparum/P. vivax mixed infection (subpopulation A, late rainfall).

Concerning each A. darlingi genetic subpopulation, overall Plasmodium infection rate was higher in subpopulation A (4.6%) when compared to subpopulation B (3.0%) or admixed individuals (3.0%) but this difference was non-significant (Fisher’s exact test: P = 0.692). For subpopulation B and admixed individuals, it was possible to compare infection rates between seasons but no significant differences were found (Table 6). In subpopulation B, there was an increase of infection rates from the early (2.7%) to the late rainfall (4.2%) period (Fisher’s exact test, P = 0.550), while admixed individuals showed an opposite trend (early rain: 4.3%, late rain: 2.3%; Fisher’s exact test: P = 0.430).

Table 6 Number of Anopheles darlingi infected with Plasmodium sp. per lineage sampling period

Discussion

In Amazonian, suburban, riverine populations of A. darlingi near Porto Velho (RO), two genetically differentiated subpopulations associated with different seasonal sampling periods were detected. Subpopulation A is almost exclusively found during or at the end of the rainy season. However, subpopulation A was less abundant in the inland locality (Jaci Paraná) where it appeared to be replaced by subpopulation B earlier in the rainy season than in riverine localities. In this sense, subpopulation A seems to be ecologically closer to a typical forest/riverine A. darlingi population that displays the highest mosquito density after the peak of rainfall[18]. In contrast, subpopulation B was the most frequent at the beginning of the rainy season, before A. darlingi reaches its late rainy season peak. A previous mtDNA-based study conducted in the same region also reported a seasonal difference in the mosquito population of Portuchuelo, a riverine village ca. 30 km northeast of Porto Velho, but no seasonal mtDNA differences were found in the other localities analysed[9]. Other studies, using most of the same microsatellite markers here analysed, compared samples from broader geographical ranges (continental scale) and revealed patterns of population structure correlated with geographic distance[1517]. In these studies, the F ST estimates among samples at smaller distances (up to 80 km) were comparable with the F ST estimates within subpopulations obtained in this study. However, none of these microsatellite-based studies compared samples from the same locality at different seasons.

Two hypotheses can be proposed to explain the observed pattern of seasonal genetic variation. The genetic differences between sampling periods may reflect drastic fluctuations in the effective population size of a single genetic entity caused by the transition between rainy and dry seasons. This would agree with the marked seasonal variation in mosquito abundance reported for this region, where during the dry season A. darlingi densities decline to residual values[21]. If this was the case then one might expect differences in the level of genetic diversity between subpopulations. Subpopulation B, which prevails at the beginning of the dry season, would display lower levels of genetic diversity and a signal of a recent bottleneck. This was not the case in the present study. Both subpopulations presented similar levels of genetic diversity and a signal of population perturbation (expansion) was detected only in subpopulation A.

The other hypothesis assumes the occurrence of two genetically distinct and partially isolated subpopulations of A. darlingi. Subpopulation A would represent a subpopulation more dependent of the rainfall regime, probably adapted to exploit temporary larval habitats in forested, riverine, flooded areas. The bio-ecology of this subpopulation would thus be closer to the typical A. darlingi. This would agree with its predominance at the end of the rainy season, where a signal of population expansion was detected and the decline over the dry season as evidenced by its near absence from collections at the early rainfall period. On the other hand, subpopulation B would correspond to a subpopulation more adapted to the dry season, in agreement with its predominance in collections at the onset of the rains. This subpopulation could have derived from the typical subpopulation A through adaptation to more permanent anthropogenic larval habitats, in urban and semi-urban settings. This would result in the maintenance of a resident A. darlingi population in urban areas during the dry season with a relatively stable effective population size. In the urban and suburban localities sampled, this population should thus be the first to peak at the start of the rainy season while subpopulation A would peak later in the rainy season when the availability of temporal larval habitats is highest in the forest/riverine flooded areas surrounding these localities.

In the Amazonian region, the variety of larval habitats explored by A. darlingi is large, ranging from typical natural breeding sites such as the margins of slow flowing streams, river arms and lagoons, to more human-dependent sites such as irrigation canals, flooded crop fields and mine excavation pools. Such diversity in larval ecology may set a ground for genetic divergence driven by ecological adaptation at the larval stage. This was the case of the recent radiation of the malaria vector Anopheles coluzzii (formerly denoted as molecular form M[45]) within the Afrotropical Anopheles gambiae complex. This species is considered to be more adapted to permanent breeding sites due to superior predator avoidance[46]. In addition, there is evidence for a recent adaptation of central African populations of A. coluzzii to densely urbanized areas by occupation of permanent breeding sites of an anthropogenic nature and a high degree of pollution[4749]. These studies demonstrate that ecological heterogeneities derived from human-made environmental changes may promote rapid ecological niche adaptation associated with evolutionary divergence in mosquito vectors.

Further studies involving additional mosquito surveys in the region of Porto Velho are required to validate the proposed hypotheses. These studies will also allow verifying if the observed pattern of seasonal genetic partitioning within A. darlingi is a stable condition or if it was a sporadic/transient observation. The complete genome sequence of A. darlingi has been published recently[5]. This remarkable achievement will result in a much greater availability of polymorphic markers (e g, single nucleotide polymorphism) to study genetic variation of this vector. Combined with appropriate field surveys, it will be possible to obtain a much better clarification of the extent of genetic divergence between subpopulations and its association with local adaptation of A. darlingi in the region.

The vast majority of individuals within subpopulation A were classified as pure-bred while the genetic background of subpopulation B was more admixed, with 48% of individuals classified as hybrids. The different proportions of admixture within subpopulations suggest a pattern of asymmetric introgression, with more subpopulation A alleles being introgressed into subpopulation B. This pattern is consistent with the hypothesis of the presence of a resident subpopulation (subpopulation B) in the urban setting that is seasonally exposed to a genetically distinct subpopulation (subpopulation A) promoting local hybridization during the rainy season. It remains to be determined which factors may determine the unbalanced genetic exchange between subpopulations A and B.

Asymmetric introgression is not an uncommon event in mosquitoes. Similar findings have been reported between closely related taxa in Culex pipiens[50, 51] and A. gambiae[52]. These patterns have been related with ecological and behavioural differences between the interbreeding taxa (e g, mating strategies)[49]. While the detection of only F2 hybrids within each subpopulation may suggest that selection is acting against F1 hybrids[53], this result should be interpreted with caution. The Bayesian approach implemented in NewHybrids requires a large number of highly diagnostic markers between populations to accurately assign individuals to hybrid classes[40, 53].

From the locus specific F ST estimates, there were two loci (P2 and P9) that displayed the highest differentiation between subpopulations. These loci may be located near genomic regions under divergent selection that may harbour genes associated with phenotypic differences between the subpopulations that may promote sympatric divergence with gene flow[54, 55]. It is noteworthy, however, that these two loci differ from the remaining ones in two aspects: they were isolated from a salivary gland cDNA library[26] and they also differ in the repeat motif (trinucleotides), factors that can affect comparisons with the remainder loci.

The findings of this study suggest that A. darlingi appears to be adapting to human modified habitats by deviating its breeding preferences into more permanent larval habitats, which may lead to an increased persistence during the dry season. This adaptation may have important consequences in regional malaria transmission patterns. Both subpopulations were found infected with malaria parasites at comparable rates, with the highest values occurring at the end of the rainy season, in accordance with the seasonal patterns of malaria transmission intensity characteristic of the Amazonian region[56]. However, subpopulation B was also found to be infected by malaria parasites in the early rainfall period, and with considerable infection rates. Moreover, admixed individuals had a higher infection rate in the period after the dry season. Although the presence of parasite DNA in a mosquito does not fully incriminate it as infective (only presence of sporozoites in salivary glands), it is still an indicator of the potential for malaria transmission given the recognized vector competence of A. darlingi. The co-occurrence of both subpopulations in this region may thus promote a more perennial transmission of malaria throughout the year, with obvious implications for malaria epidemiology and control.

Several studies have demonstrated that deforestation and human colonization of deforested areas have a strong impact in mosquito populations, with an increase in larval breeding and human-biting activity[7, 57, 58] thereby increasing malaria risk[59, 60]. Moreover, the increased landscape disturbance resulting from the construction of two hydroelectric power plants in the Madeira River may promote greater changes in the genetic structure and bio-ecology of this malaria vector, with undetermined consequences for malaria transmission[61]. In this context, the results of this study may serve as a baseline for future bio-ecological and genetic monitoring of A. darlingi populations in this malaria endemic region.

Conclusions

The co-occurrence of two genetically distinct subpopulations of A. darlingi adapted to different periods of rainfall may promote a more perennial transmission of malaria throughout the year. These findings in a context of strong environmental impact due to deforestation and dam construction preclude serious implications for malaria epidemiology and control in the Amazonian region.

Abbreviations

mtDNA:

Mitochondrial DNA

INMET:

Instituto Nacional de Meteorologia

PCR:

Polimerase Chain Reaction

rRNA:

Ribosomal ribonucleic acid

dNTPs:

Deoxynucleoside triphosphate

SMM:

Stepwise Mutation Model

TPM:

Two Phase Model

He:

heterozygosity, heterozygosity (He) allele richness (Rs) and private allele richness (Ps).

References

  1. Hiwat H, Bretas G: Ecology of Anopheles darlingi Root with respect to vector importance: a review. Parasit Vectors. 2011, 4: 177-10.1186/1756-3305-4-177.

    Article  PubMed Central  PubMed  Google Scholar 

  2. WHO: World Malaria Report. 2013, Geneva: World Health Organization, 284-[http://www.who.int/malaria/publications/world_malaria_report_2013/en/]

    Google Scholar 

  3. Sinka ME, Rubio-Palis Y, Manguin S, Patil AP, Temperley WH, Gething PW, Van Boeckel T, Kabaria CW, Harbach RE, Hay SI: The dominant Anopheles vectors of human malaria in the Americas: occurrence data, distribution maps and bionomic précis. Parasit Vectors. 2010, 3: 72-10.1186/1756-3305-3-72.

    Article  PubMed Central  PubMed  Google Scholar 

  4. Loaiza JR, Bermingham E, Scott ME, Rovira JR, Conn JE: Species composition and distribution of adult Anopheles (Diptera; Culicidae) in Panama. J Med Entomol. 2008, 45: 841-851. 10.1603/0022-2585(2008)45[841:SCADOA]2.0.CO;2.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  5. Marinotti O, Cerqueira GC, de Almeida LG, Ferro MI, Loreto EL, Zaha A, Teixeira SM, Wespiser AR, Almeida E, Silva A, Schlindwein AD, Pacheco AC, Silva AL, Graveley BR, Walenz BP, Lima Bde A, Ribeiro CA, Nunes-Silva CG, de Carvalho CR, Soares CM, de Menezes CB, Matiolli C, Caffrey D, Araújo DA, de Oliveira DM, Golenbock D, Grisard EC, Fantinatti-Garboggini F, de Carvalho FM, Barcellos FG: The genome of Anopheles darlingi, the main neotropical malaria vector. Nucleic Acids Res. 2013, 41: 7387-7400. 10.1093/nar/gkt484.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  6. Moreno M, Marinotti O, Krzywinski J, Tadei WP, James AA, Achee NL, Conn JE: Complete mtDNA genomes of Anopheles darlingi and an approach to anopheline divergence time. Malar J. 2010, 9: 127-10.1186/1475-2875-9-127.

    Article  PubMed Central  PubMed  Google Scholar 

  7. Vittor AY, Gilman RH, Tielsch J, Glass G, Shields T, Lozano WS, Pinedo-Cancino V, Patz JA: The effect of deforestation on the human-biting rate of Anopheles darlingi, the primary vector of falciparum malaria in the Peruvian Amazon. Am J Trop Med Hyg. 2006, 74: 3-11.

    PubMed  Google Scholar 

  8. Charlwood J: Biological variation in Anopheles darlingi Root. Mem Inst Oswaldo Cruz. 1996, 91: 391-398.

    CAS  PubMed  Google Scholar 

  9. Angêlla AF, Gil LHS, Silva LHP, Ribolla PEM: Population structure of the malaria vector Anopheles darlingi in Rondônia, Brazilian Amazon, based on mitochondrial DNA. Mem Inst Oswaldo Cruz. 2007, 102: 953-958. 10.1590/S0074-02762007000800010.

    Article  PubMed  Google Scholar 

  10. Manguin S, Wilkerson RC, Conn JE, Rubio-Palis Y, Danoff-Burg JA, Roberts DR: Population structure of the primary malaria vector in South America, Anopheles darlingi, using isozyme, random amplified polymorphic DNA, internal transcribed spacer 2, and morphologic markers. Am J Trop Med Hyg. 1999, 60: 364-376.

    CAS  PubMed  Google Scholar 

  11. Santos JM, Lobo JA, Tadei WP, Contel EPB: Intrapopulational genetic differentiation in Anopheles (N.) darlingi Root, 1926 (Diptera: Culicidae) in the amazon region. Genet Mol Biol. 1999, 22: 325-331. 10.1590/S1415-47571999000300007.

    Article  Google Scholar 

  12. Conn JE: Systematics and population level analysis of Anopheles darlingi. Mem Inst Oswaldo Cruz. 1998, 93: 647-650. 10.1590/S0074-02761998000500015.

    Article  CAS  PubMed  Google Scholar 

  13. Mirabello L, Conn JE: Molecular population genetics of the malaria vector Anopheles darlingi in Central and South America. Heredity. 2006, 96: 311-321. 10.1038/sj.hdy.6800805.

    Article  CAS  PubMed  Google Scholar 

  14. Pedro PM, Sallum M: Spatial expansion and population structure of the neotropical malaria vector, Anopheles darlingi (Diptera: Culicidae). Biol J Linnean Soc. 2009, 97: 854-866. 10.1111/j.1095-8312.2009.01226.x.

    Article  Google Scholar 

  15. Conn JE, Vineis J, Bollback JP, Onyabe DY, Wilkerson C, Póvoa MM: Population structure of the malaria vector Anopheles darlingi in a malaria-endemic region of eastern Amazonian Brazil. Am J Trop Med Hyg. 2006, 74: 798-806.

    PubMed  Google Scholar 

  16. Mirabello L, Vineis JH, Yanoviak SP, Scarpassa VM, Póvoa MM, Padilla N, Achee NL, Conn JE: Microsatellite data suggest significant population structure and differentiation within the malaria vector Anopheles darlingi in Central and South America. BMC Ecol. 2008, 8: 3-10.1186/1472-6785-8-3.

    Article  PubMed Central  PubMed  Google Scholar 

  17. Scarpassa VM, Conn JE: Population genetic structure of the major malaria vector Anopheles darlingi (Diptera: Culicidae) from the Brazilian Amazon, using microsatellite markers. Mem Inst Oswaldo Cruz. 2007, 102: 319-327. 10.1590/S0074-02762007005000045.

    Article  CAS  PubMed  Google Scholar 

  18. Gil LHS, Alves FP, Zieler H, Salcedo JMV, Durlacher RR, Almeida Cunha RP, Tada MS, Camargo LMA, Camargo EP, Pereira Da Silva LH: Seasonal malaria transmission and variation of anopheline density in two distinct endemic areas in Brazilian Amazonia. J Med Entomol. 2003, 40: 636-641. 10.1603/0022-2585-40.5.636.

    Article  PubMed  Google Scholar 

  19. Gil LHS, Tada MS, Katsuragawa TH, Ribolla PEM, Pereira-Da-Silva LH: Urban and suburban malaria in Rondônia (Brazilian Western Amazon) II: Perennial transmissions with high anopheline densities are associated with human environmental changes. Mem Inst Oswaldo Cruz. 2007, 102: 271-276. 10.1590/S0074-02762007005000013.

    Article  PubMed  Google Scholar 

  20. Souza-Santos R: Seasonal distribution of malaria vectors in Machadinho d'Oeste, Rondônia State, Amazon Region, Brazil. Cad Saúde Pública. 2002, 18: 1813-1818. 10.1590/S0102-311X2002000600039.

    Article  PubMed  Google Scholar 

  21. Moutinho PR, Gil LHS, Cruz RB, Ribolla PEM: Population dynamics, structure and behavior of Anopheles darlingi in a rural settlement in the Amazon rainforest of Acre. Brazil Malar J. 2011, 10: 174-10.1186/1475-2875-10-174.

    Article  PubMed  Google Scholar 

  22. Lindsey R: NASA Earth Observatory. 2009, [http://earthobservatory.nasa.gov/Features/WorldOfChange/deforestation.php]

    Google Scholar 

  23. Silva MJ, Saraiva FA, Araujo ML: Aspectos Climáticos de Porto Velho - Rondônia. XIII Congresso Brasileiro de Meteorologia. 2004, Fortaleza – CE: Meteorologia e o Desenvolvimento Sustentável

    Google Scholar 

  24. Consoli RAGB, Lourenço de Oliveira R: Principais mosquitos de importância sanitária no Brasil. 1994, Rio de Janeiro: Editora FIOCRUZ, 228-

    Google Scholar 

  25. Conn JE, Bollback JP, Onyabe DY, Robinson TN, Wilkerson RC, Póvoa MM: Isolation of polymorphic microsatellite markers from the malaria vector Anopheles darlingi. Mol Ecol Not. 2001, 1: 223-225. 10.1046/j.1471-8278.2001.00078.x.

    Article  CAS  Google Scholar 

  26. Calvo E, Pham V, Marinotto O, Andersen JF, Ribeiro JM: The salivary gland transcriptome of the neotropical malaria vector Anopheles darlingi reveals accelerated evolution of genes relevant to hematophagy. BMC Genomics. 2009, 10: 57-10.1186/1471-2164-10-57.

    Article  PubMed Central  PubMed  Google Scholar 

  27. Kimura M, Kaneko O, Liu Q, Zhou M, Kawamoto F, Wataya Y, Otani S, Yamaguchi Y, Tanabe K: Identification of the four species of human malaria parasites by nested PCR that targets variant sequences in the small subunit rRNA gene. Parasitol Int. 1997, 46: 91-95. 10.1016/S1383-5769(97)00013-5.

    Article  CAS  Google Scholar 

  28. Goudet J: FSTAT (Version 1.2): A computer Program to calculate F-Statistics. J Hered. 1995, 86: 485-486.

    Google Scholar 

  29. Kalinowski ST: HP-RARE 1.0: a computer program for performing rarefaction on measures of allelic richness. Mol Ecol Not. 2004, 5: 187-189.

    Article  Google Scholar 

  30. Excoffier L, Laval G, Schneider S: ARLEQUIN ver. 3.0: An integrated software package for population genetics data analysis. Evol Bioinform Online. 2005, 1: 47-50.

    PubMed Central  CAS  Google Scholar 

  31. Raymond M, Rousset F: GENEPOP (version 1.2): population genetics software for exact tests and ecumenicism. J Hered. 1995, 86: 248-249.

    Google Scholar 

  32. Rousset F: Genepop'007: a complete re-implementation of the genepop software for Windows and Linux. Mol Ecol Resour. 2008, 8: 103-106. 10.1111/j.1471-8286.2007.01931.x.

    Article  PubMed  Google Scholar 

  33. Van Oosterhout C, Hutchinson W, Wills D, Shipley P: Micro-Checker: software for identifying and correcting genotyping errors in microsatellite data. Mol Ecol Not. 2004, 4: 535-538. 10.1111/j.1471-8286.2004.00684.x.

    Article  CAS  Google Scholar 

  34. Weir B, Cockerham C: Estimating F-statistics for the analysis of population structure. Evolution. 1984, 38: 1358-1370. 10.2307/2408641.

    Article  Google Scholar 

  35. Excoffier L, Smouse P, Quattro J: Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction data. Genetics. 1992, 131: 479-491.

    PubMed Central  CAS  PubMed  Google Scholar 

  36. Pritchard J, Stephens M, Donnelly P: Inference of population structure using multilocus genotype data. Genetics. 2000, 155: 945-959.

    PubMed Central  CAS  PubMed  Google Scholar 

  37. Evanno G, Regnaut S, Goudet J: Detecting the number of clusters of individuals using the software Structure: a simulation study. Mol Ecol. 2005, 14: 2611-2620. 10.1111/j.1365-294X.2005.02553.x.

    Article  CAS  PubMed  Google Scholar 

  38. Earl D, von Holdt B: STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv Genet Resour. 2011, 4: 359-361.

    Article  Google Scholar 

  39. Jakobsson M, Rosenberg NA: CLUMPP: a cluster matching and permutation program for dealing with label switching and multimodality in analysis of population structure. Bioinformatics. 2007, 23: 1801-1806. 10.1093/bioinformatics/btm233.

    Article  CAS  PubMed  Google Scholar 

  40. Vaha J-P, Primmer CR: Efficiency of model-based Bayesian methods for detecting hybrid individuals under different hybridization scenarios and with different numbers of loci. Mol Ecol. 2006, 15: 63-72.

    Article  CAS  PubMed  Google Scholar 

  41. Anderson EC, Thompson EA: A model-based method for identifying species hybrids using multilocus genetic data. Genetics. 2002, 160: 217-1229.

    Google Scholar 

  42. Piry S, Luikart G, Cornuet J: BOTTLENECK: a computer program for detecting recent reductions in the effective population size using allele frequency data. J Hered. 1999, 90: 502-503. 10.1093/jhered/90.4.502.

    Article  Google Scholar 

  43. Cornuet J-M, Luikart G: Description and power analysis of two tests for detecting recent population bottlenecks from allele frequency data. Genetics. 1996, 144: 2001-2014.

    PubMed Central  CAS  PubMed  Google Scholar 

  44. Holm S: A simple sequentially rejective multiple test procedure. Scand J Statist. 1979, 6: 65-70.

    Google Scholar 

  45. Coetzee M, Hunt RH, Wilkerson R, della Torre A, Coulibaly MB, Besansky NJ: Anopheles coluzzii and Anopheles amharicus, new members of the Anopheles gambiae complex. Zootaxa. 2013, 3619: 246-274.

    Article  PubMed  Google Scholar 

  46. Diabaté A, Dabiré RK, Heidenberger K, Crawford J, Lamp WO, Culler LE, Lehmann T: Evidence for divergent selection between the molecular forms of Anopheles gambiae: role of predation. BMC Evol Biol. 2008, 8: 5-10.1186/1471-2148-8-5.

    Article  PubMed Central  PubMed  Google Scholar 

  47. Simard F, Ayala D, Kamdem GC, Pombi M, Etouna J, Ose K, Fotsing J-M, Fontenille D, Besansky NJ, Costantini C: Ecological niche partitioning between Anopheles gambiae molecular forms in Cameroon: the ecological side of speciation. BMC Ecol. 2009, 9: 17-10.1186/1472-6785-9-17.

    Article  PubMed Central  PubMed  Google Scholar 

  48. Kamdem C, Tene Fossog B, Simard F, Etouna J, Ndo C, Kengne P, Boussès P, Etoa FX, Awono-Ambene P, Fontenille D, Antonio-Nkondjio C, Besansky NJ, Costantini C: Anthropogenic habitat disturbance and ecological divergence between incipient species of the malaria mosquito Anopheles gambiae. PLoS One. 2012, 7: e39453-10.1371/journal.pone.0039453.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  49. Fossog Tene B, Poupardin R, Costantini C, Awono-Ambene P, Wondji CS, Ranson H, Antonio-Nkondjio C: Resistance to DDT in an urban setting: common mechanisms implicated in both M and S forms of Anopheles gambiae in the city of Yaoundé, Cameroon. PLoS One. 2013, 8: e61408-10.1371/journal.pone.0061408.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  50. Gomes B, Sousa CA, Novo MT, Freitas FB, Alves R, Côrte-Real AR, Salgueiro P, Donnelly MJ, Almeida APJ, Pinto J: Asymmetric introgression between sympatric molestus and pipiens forms of Culex pipiens (Diptera: Culicidae) in the Comporta region, Portugal. BMC Evol Biol. 2009, 9: 262-10.1186/1471-2148-9-262.

    Article  PubMed Central  PubMed  Google Scholar 

  51. Gomes B, Kioulos E, Papa A, Almeida AP, Vontas J, Pinto J: Distribution and hybridization of Culex pipiens forms in Greece during the West Nile virus outbreak of 2010. Infect Genet Evol. 2013, 16: 218-225.

    Article  PubMed  Google Scholar 

  52. Marsden CD, Lee Y, Nieman CC, Sanford MR, Dinis J, Martins C, Rodrigues A, Cornel AJ, Lanzaro GC: Asymmetric introgression between the M and S forms of the malaria vector, Anopheles gambiae, maintains divergence despite extensive hybridization. Mol Ecol. 2011, 20: 4983-4994. 10.1111/j.1365-294X.2011.05339.x.

    Article  PubMed Central  PubMed  Google Scholar 

  53. Gow JL, Peichel CL, Taylor EB: Contrasting hybridization rates between sympatric three-spined sticklebacks highlight the fragility of reproductive barriers between evolutionarily young species. Mol Ecol. 2006, 15: 739-752. 10.1111/j.1365-294X.2006.02825.x.

    Article  PubMed  Google Scholar 

  54. Nosil P: Speciation with gene flow could be common. Mol Ecol. 2008, 17: 2103-2106. 10.1111/j.1365-294X.2008.03715.x.

    Article  PubMed  Google Scholar 

  55. Pinho C, Hey J: Divergence with gene flow: models and data. Annu Rev Ecol Evol Syst. 2010, 41: 215-230. 10.1146/annurev-ecolsys-102209-144644.

    Article  Google Scholar 

  56. Charlwood JD: Observations on the bionomics of Anopheles darlingi Root (Diptera; Culicidae) from Brazil. Bull Entomol Res. 1980, 70: 685-692. 10.1017/S0007485300007975.

    Article  Google Scholar 

  57. Vittor AY, Pan W, Gilman RH, Tielsch J, Sánchez-Lozano W, Salas-Cobos E, Flores S, Patz JA: Linking deforestation to malaria in the Amazon: characterization of the breeding habitat of the principal malaria vector, Anopheles darlingi. Am J Trop Med Hyg. 2009, 81: 5-12.

    PubMed Central  PubMed  Google Scholar 

  58. Hahn MB, Gangnon RE, Barcellos C, Asner GP, Patz JA: Influence of deforestation, logging, and fire on malaria in the Brazilian Amazon. Plos ONE. 2014, 9: e85725-10.1371/journal.pone.0085725.

    Article  PubMed Central  PubMed  Google Scholar 

  59. Povoa MM, Conn JE, Schlichting CD, Amaral JC, Sequra MN, Da Silva AN, Dos Santos CC, Lacerda RN, De Souza RT, Galiza D, Santa Rosa EP, Wirtz RA: Malaria vectors, epidemiology, and the re-emergence of Anopheles darlingi in Belem, Para, Brazil. J Med Entomol. 2003, 40: 379-386. 10.1603/0022-2585-40.4.379.

    Article  PubMed  Google Scholar 

  60. Harris AF, Matias-Arnez A, Hill N: Biting time of Anopheles darlingi in the Bolivian Amazon and implications for control of malaria. Trans R Soc Trop Med Hyg. 2006, 100: 45-47. 10.1016/j.trstmh.2005.07.001.

    Article  PubMed  Google Scholar 

  61. Keiser J, Castro MCD, Maltese MF, Bos R, Tanner M, Singer BH, Utzinger J: Effect of irrigation and large dams on the burden of malaria on a global and regional scale. Am J Trop Med Hyg. 2005, 72: 392-406.

    PubMed  Google Scholar 

Download references

Acknowledgements

This work received financial support from FAPESP and CAPES/FCT Bilateral Agreemen (Transnational Cooperation Project: capes128587311218103). PS was funded by a post-doc fellowship (SFRH/BPD/72532/2010) from Fundação para a Ciência e Tecnologia (FCT) /Ministério da Educação e Ciência (MEC) through the European Social Fund (ESF). AFA was funded by CAPES and PMR has a CNPq fellowship. We would like to thank Jan Conn for suggestions that improved the manuscript.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Paulo EM Ribolla.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

All the authors collaborated on the work presented in this study. AFA and PEMR defined the research theme; LHSG was responsible for the mosquitoes collection; AFA designed the methods, performed the experiments and interpreted the results; JLV performed some of the experiments; AFA, PS, JP, and PEMR analysed the data and drafted the manuscript. All authors read and approved the final version of the manuscript.

Aline F Angêlla, Patrícia Salgueiro contributed equally to this work.

Electronic supplementary material

12936_2014_3251_MOESM1_ESM.docx

Additional file 1: Title: Repeats and primer sequences of the microsatellite loci designed for A. darlingi. Description: Microssatelites used to genotype Anopheles darlingi, Locus position, numer of repeats, primer sequence , Tm of the reaction and fluorechrome of each loci used. (DOCX 64 KB)

12936_2014_3251_MOESM2_ESM.docx

Additional file 2: Title: Estimates of pairwise genetic differentiation among localities in A. darlingi pure subpopulation (without admixed): A) subpopulation A, B) subpopulation B. Description: Pairwise genetic differentiation of grouped subpopulations A and B excluding admixture individuals. (DOCX 72 KB)

12936_2014_3251_MOESM3_ESM.docx

Additional file 3: Title: Relative distribution of each genetic lineage revealed by STRUCTURE per locality/semester. Description: A table showing the distribution of subpopulation A, subpopulation B and admixture population at different collection sites and period. (DOCX 75 KB)

12936_2014_3251_MOESM4_ESM.docx

Additional file 4: Title: Apportionment of molecular variance measured among populations of Anopheles darlingi from all samples, or between populations from the two semesters. Description: Molecular variance of samples considering subpopulation A and subpopulation B. (DOCX 52 KB)

12936_2014_3251_MOESM5_ESM.docx

Additional file 5: Title: Summary of BOTTLENECK tests. Description: List of Bottleneck tests at Subpopulations A and B at each period of collection. (DOCX 42 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/.

The Creative Commons Public Domain Dedication waiver (https://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Angêlla, A.F., Salgueiro, P., Gil, L.H. et al. Seasonal genetic partitioning in the neotropical malaria vector, Anopheles darlingi. Malar J 13, 203 (2014). https://0-doi-org.brum.beds.ac.uk/10.1186/1475-2875-13-203

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/1475-2875-13-203

Keywords