Background:The major Neotropical malaria vector, Anopheles darlingi, was reintroduced into the Iquitos, Loreto, Peru area during the early 1990s, where it displaced other anophelines and caused a major malaria epidemic. Since then, case numbers in Loreto have fluctuated, but annual increases have been reported since 2012. Methods:The population genetic structure of An. darlingi sampled before and after the introduction of long-lasting insecticidal nets (LLINs) was investigated to test the hypothesis of temporal population change (2006 vs. 2012). Current samples of An. darlingi were used to test the hypothesis of ecological adaptation to human modified (highway) compared with wild (riverine) habitat, linked to forest cover. In total, 693 An. darlingi from nine localities in Loreto, Peru area were genotyped using 13 microsatellite loci. To test the hypothesis of habitat differentiation in An. darlingi biting time patterns, HBR and EIR, four collections of An. darlingi from five localities (two riverine and three highway) were analysed.Results: Analyses of microsatellite loci from seven (2006) and nine settlements (2012)(2013)(2014) in the Iquitos area detected two distinctive populations with little overlap, although it is unclear whether this population replacement event is associated with LLIN distribution or climate. Within the 2012-2014 population two admixed subpopulations, A and B, were differentiated by habitat, with B significantly overrepresented in highway, and both in near-equal proportions in riverine. Both subpopulations had a signature of expansion and there was moderate genetic differentiation between them. Habitat and forest cover level had significant effects on HBR, such that Plasmodium transmission risk, as measured by EIR, in peridomestic riverine settlements was threefold higher than in peridomestic highway settlements. HBR was directly associated with available host biomass rather than forest cover. Conclusions: A population replacement event occurred between 2006 and 2012-2014, concurrently with LLIN distribution and a moderate El Niño event, and prior to an increase in malaria incidence. The likely drivers of this replacement cannot be determined with current data. The present-day An. darlingi population is composed of two highly admixed subpopulations, which appear to be in an early stage of differentiation, triggered by anthropogenic alterations to local habitat.

The primary Brazilian malaria vector, Nyssorhynchus darlingi (formerly Anopheles darlingi ), ranges from 0°S–23°S across three biomes (Amazonia, Cerrado, Mata Atlântica). Rising temperatures will increase mosquito developmental rates, and models predict future malaria transmission by Ny . darlingi in Brazil will shift southward. We reared F 1 Ny . darlingi (progeny of field-collected females from 4 state populations across Brazil) at three temperatures (20, 24, 28 °C) and measured key life-history traits. Our results reveal geographic variation due to both genetic differences among localities and plastic responses to temperature differences. Temperature significantly altered all traits: faster larval development, shorter adult life and overall lifespan, and smaller body sizes were seen at 28 °C versus 20 °C. Low-latitude Amazonia mosquitoes had the fastest larval development at all temperatures, but at 28 °C, average development rate of high-latitude Mata Atlântica mosquitoes was accelerated and equivalent to low-latitude Amazonia. Body size of adult mosquitoes from the Mata Atlântica remained larger at all temperatures. We detected genetic variation in the plastic responses among mosquitoes from different localities, with implications for malaria transmission under climate change. Faster development combined with larger body size, without a tradeoff in adult longevity, suggests vectorial capacities of some Mata Atlântica populations may significantly increase under warming climates.