Forest cover percentage drives the peak biting time of Nyssorhynchus darlingi (Diptera: Culicidae) in the Brazilian Amazon.

BACKGROUND: Deforestation is an important driver of malaria dynamics, with a relevant impact on mosquito ecology, including larval habitat availability, blood-feeding behaviour, and peak biting time. The latter is one of several entomological metrics to evaluate vectorial capacity and effectiveness of disease control. This study aimed to test the effect of forest cover percentage on the peak biting time of Plasmodium-uninfected and infected Nyssorhynchus darlingi females. METHODS: Mosquitoes were captured utilizing human landing catch (HLC) in the peridomestic habitat in field collections carried out in the wet, wet-dry transition, and dry seasons from 2014 to 2017 in areas with active malaria transmission in Amazonian Brazil. The study locations were in rural settlements in areas with the mean annual malaria parasite incidence (Annual Parasite Incidence, API ≥ 30). All Ny. darlingi females were tested for Plasmodium spp. infection using real time PCR technique. Forest cover percentage was calculated for each collection site using QGIS v. 2.8 and was categorized in three distinct deforestation scenarios: (1) degraded, < 30% forest cover, (2) intermediate, 30-70% forest cover, and (3) preserved, > 70% forest cover. RESULTS: The highest number of uninfected female Ny. darlingi was found in degraded landscape-sites with forest cover < 30% in any peak biting time between 18:00 and 0:00. Partially degraded landscape-sites, with (30-70%) forest cover, showed the highest number of vivax-infected females, with a peak biting time of 21:00-23:00. The number of P. falciparum-infected mosquitoes was highest in preserved sites with > 70% forest cover, a peak biting at 19:00-20:00, and in sites with 30-70% forest cover at 22:00-23:00. CONCLUSIONS: Results of this study show empirically that degraded landscapes favour uninfected Ny. darlingi with a peak biting time at dusk (18:00-19:00), whereas partially degraded landscapes affect the behaviour of Plasmodium-infected Ny. darlingi by shifting its peak biting time towards hours after dark (21:00-23:00). In preserved sites, Plasmodium-infected Ny. darlingi bite around dusk (18:00-19:00) and shortly after (19:00-20:00).

Effect of spatiotemporal variables on abundance, biting activity and parity of Nyssorhynchus darlingi (Diptera: Culicidae) in peri-Iquitos, Peru.

BACKGROUND: In malaria endemic regions of the Peruvian Amazon, rainfall together with river level and breeding site availability drive fluctuating vector mosquito abundance and human malaria cases, leading to temporal heterogeneity. The main variables influencing spatial transmission include location of communities, mosquito behaviour, land use/land cover, and human ecology/behaviour. The main objective was to evaluate seasonal and microgeographic biting behaviour of the malaria vector Nyssorhynchus (or Anopheles) darlingi in Amazonian Peru and to investigate effects of seasonality on malaria transmission. METHODS: We captured mosquitoes from 18:00 to 06:00 h using Human Landing Catch in two riverine (Lupuna, Santa Emilia) and two highway (El Triunfo, Nuevo Horizonte) communities indoors and outdoors from 8 houses per community, during the dry and rainy seasons from February 2016 to January 2017. We then estimated parity rate, daily survival and age of a portion of each collection of Ny. darlingi. All collected specimens of Ny. darlingi were tested for the presence of Plasmodium vivax or Plasmodium falciparum sporozoites using real-time PCR targeting the small subunit of the 18S rRNA. RESULTS: Abundance of Ny. darlingi varied across village, season, and biting behaviour (indoor vs outdoor), and was highly significant between rainy and dry seasons (p < 0.0001). Biting patterns differed, although not significantly, and persisted regardless of season, with peaks in highway communities at ~ 20:00 h in contrast to biting throughout the night (i.e., 18:00-06:00) in riverine communities. Of 3721 Ny. darlingi tested for Plasmodium, 23 (0.62%) were infected. We detected Plasmodium-infected Ny. darlingi in both community types and most (20/23) were captured outdoors during the rainy season; 17/23 before midnight. Seventeen Ny. darlingi were infected with P. vivax, and 6 with P. falciparum. No infected Ny. darlingi were captured during the dry season. Significantly higher rates of parity were detected in Ny. darlingi during the rainy season (average 64.69%) versus the dry season (average 36.91%) and by community, Lupuna, a riverine village, had the highest proportion of parous to nulliparous females during the rainy season. CONCLUSIONS: These data add a seasonal dimension to malaria transmission in peri-Iquitos, providing more evidence that, at least locally, the greatest risk of malaria transmission is outdoors during the rainy season mainly before midnight, irrespective of whether the community was located adjacent to the highway or along the river.

Mutations Linked to Insecticide Resistance Not Detected in the Ace-1 or VGSC Genes in Nyssorhynchus darlingi from Multiple Localities in Amazonian Brazil and Peru.

Indoor residual spray (IRS), mainly employing pyrethroid insecticides, is the most common intervention for preventing malaria transmission in many regions of Latin America; the use of long-lasting insecticidal nets (LLINs) has been more limited. Knockdown resistance (kdr) is a well-characterized target-site resistance mechanism associated with pyrethroid and DDT resistance. Most mutations detected in acetylcholinesterase-1 (Ace-1) and voltage-gated sodium channel (VGSC) genes are non-synonymous, resulting in a change in amino acid, leading to the non-binding of the insecticide. In the present study, we analyzed target-site resistance in Nyssorhynchus darlingi, the primary malaria vector in the Amazon, in multiple malaria endemic localities. We screened 988 wild-caught specimens of Ny. darlingi from three localities in Amazonian Peru and four in Amazonian Brazil. Collections were conducted between 2014 and 2021. The criteria were Amazonian localities with a recent history as malaria hotspots, primary transmission by Ny. darlingi, and the use of both IRS and LLINs as interventions. Fragments of Ace-1 (456 bp) and VGSC (228 bp) were amplified, sequenced, and aligned with Ny. darlingi sequences available in GenBank. We detected only synonymous mutations in the frequently reported Ace-1 codon 280 known to confer resistance to organophosphates and carbamates, but detected three non-synonymous mutations in other regions of the gene. Similarly, no mutations linked to insecticide resistance were detected in the frequently reported codon (995) at the S6 segment of domain II of VGSC. The lack of genotypic detection of insecticide resistance mutations by sequencing the Ace-1 and VGSC genes from multiple Ny. darlingi populations in Brazil and Peru could be associated with low-intensity resistance, or possibly the main resistance mechanism is metabolic.

Natural Infection of Nyssorhynchus darlingi and Nyssorhynchus benarrochi B with Plasmodium during the Dry Season in the Understudied Low-Transmission Setting of Datem del Marañon Province, Amazonian Peru.

The persistence of malaria hotspots in Datem del Marañon Province, Peru, prompted vector control units at the Ministry of Health, Loreto Department, to collaborate with the Amazonian International Center of Excellence for Malaria Research to identify the main vectors in several riverine villages that had annual parasite indices > 15 in 2018-2019. Anophelinae were collected indoors and outdoors for two 12-hour nights/community during the dry season in 2019 using human landing catch. We identified four species: Nyssorhynchus benarrochi B, Nyssorhynchus darlingi, Nyssorhynchus triannulatus, and Anopheles mattogrossensis. The most abundant, Ny. benarrochi B, accounted for 96.3% of the total (7,550/7,844), of which 61.5% were captured outdoors (4,641/7,550). Six mosquitoes, one Ny. benarrochi B and five Ny. darlingi, were infected by Plasmodium falciparum or Plasmodium vivax. Human biting rates ranged from 0.5 to 592.8 bites per person per hour for Ny. benarrochi B and from 0.5 to 32.0 for Ny. darlingi, with entomological inoculation rates as high as 0.50 infective bites per night for Ny. darlingi and 0.25 for Ny. benarrochi B. These data demonstrate the risk of malaria transmission by both species even during the dry season in villages in multiple watersheds in Datem del Marañon province.

Vector role and human biting activity of Anophelinae mosquitoes in different landscapes in the Brazilian Amazon.

BACKGROUND: Environmental disturbance, deforestation and socioeconomic factors all affect malaria incidence in tropical and subtropical endemic areas. Deforestation is the major driver of habitat loss and fragmentation, which frequently leads to shifts in the composition, abundance and spatial distribution of vector species. The goals of the present study were to: (i) identify anophelines found naturally infected with Plasmodium; (ii) measure the effects of landscape on the number of Nyssorhynchus darlingi, presence of Plasmodium-infected Anophelinae, human biting rate (HBR) and malaria cases; and (iii) determine the frequency and peak biting time of Plasmodium-infected mosquitoes and Ny. darlingi. METHODS: Anopheline mosquitoes were collected in peridomestic and forest edge habitats in seven municipalities in four Amazon Brazilian states. Females were identified to species and tested for Plasmodium by real-time PCR. Negative binomial regression was used to measure any association between deforestation and number of Ny. darlingi, number of Plasmodium-infected Anophelinae, HBR and malaria. Peak biting time of Ny. darlingi and Plasmodium-infected Anophelinae were determined in the 12-h collections. Binomial logistic regression measured the association between presence of Plasmodium-infected Anophelinae and landscape metrics and malaria cases. RESULTS: Ninety-one females of Ny. darlingi, Ny. rangeli, Ny. benarrochi B and Ny. konderi B were found to be infected with Plasmodium. Analysis showed that the number of malaria cases and the number of Plasmodium-infected Anophelinae were more prevalent in sites with higher edge density and intermediate forest cover (30-70%). The distance of the drainage network to a dwelling was inversely correlated to malaria risk. The peak biting time of Plasmodium-infected Anophelinae was 00:00-03:00 h. The presence of Plasmodium-infected mosquitoes was higher in landscapes with > 13 malaria cases. CONCLUSIONS: Nyssorhynchus darlingi, Ny. rangeli, Ny. benarrochi B and Ny. konderi B can be involved in malaria transmission in rural settlements. The highest fraction of Plasmodium-infected Anophelinae was caught from midnight to 03:00 h. In some Amazonian localities, the highest exposure to infectious bites occurs when residents are sleeping, but transmission can occur throughout the night. Forest fragmentation favors increases in both malaria and the occurrence of Plasmodium-infected mosquitoes in peridomestic habitat. The use of insecticide-impregnated mosquito nets can decrease human exposure to infectious Anophelinae and malaria transmission.

Malaria transmission in landscapes with varying deforestation levels and timelines in the Amazon: a longitudinal spatiotemporal study.

The relationship between deforestation and malaria is a spatiotemporal process of variation in Plasmodium incidence in human-dominated Amazonian rural environments. The present study aimed to assess the underlying mechanisms of malarial exposure risk at a fine scale in 5-km2 sites across the Brazilian Amazon, using field-collected data with a longitudinal spatiotemporally structured approach. Anopheline mosquitoes were sampled from 80 sites to investigate the Plasmodium infection rate in mosquito communities and to estimate the malaria exposure risk in rural landscapes. The remaining amount of forest cover (accumulated deforestation) and the deforestation timeline were estimated in each site to represent the main parameters of both the frontier malaria hypothesis and an alternate scenario, the deforestation-malaria hypothesis, proposed herein. The maximum frequency of pathogenic sites occurred at the intermediate forest cover level (50% of accumulated deforestation) at two temporal deforestation peaks, e.g., 10 and 35 years after the beginning of the organization of a settlement. The incidence density of infected anophelines in sites where the original forest cover decreased by more than 50% in the first 25 years of settlement development was at least twice as high as the incidence density calculated for the other sites studied (adjusted incidence density ratio = 2.25; 95% CI, 1.38-3.68; p = 0.001). The results of this study support the frontier malaria as a unifying hypothesis for explaining malaria emergence and for designing specific control interventions in the Brazilian Amazon.

New Records of Anopheles benarrochi B (Diptera: Culicidae) in Malaria Hotspots in the Amazon Regions of Ecuador and Peru.

The increase in malaria transmission in the Amazon region motivated vector control units of the Ministry of Health of Ecuador and Peru to investigate Anopheles (Diptera: Culicidae) species present in transmission hotspots. Mosquitoes were collected using prokopack aspirators and CDC light traps (Ecuador) and human landing catch in Peru. In Ecuador, 84 Anopheles were captured from Pastaza, Morona Santiago, and Orellana provinces and identified morphologically [An. (An.) apicimacula Dyar and Knab, An. (Nys.) near benarrochi, An. (Nys.) near oswaldoi, An. (Nys.) near strodei, An. (An.) nimbus (Theobald, 1902), and An. (Nyssorhynchus) sp.]. In Peru, 1,150 Anopheles were collected in Andoas District. A subsample of 166 specimens was stored under silica and identified as An. near oswaldoi, An. darlingi, and An. (An.) mattogrossensis Lutz and Neiva. COI barcode region sequences were obtained for 137 adults (107 from Peru, 30 from Ecuador) identified by ITS2 PCR-RFLP as An. benarrochi Gabaldon, Cova Garcia, and Lopez and retained in the final analysis. Haplotypes from the present study plus An. benarrochi B GenBank sequences grouped separately from Brazilian An. benarrochi GenBank sequences by 44 mutation steps, indicating that the present study specimens were An. benarrochi B. Our findings confirm the presence of An. benarrochi B in Ecuador and reported here for the first time from the Amazonian provinces of Orellana and Morona Santiago. Furthermore, we confirm that the species collected in Andoas District in the Datem del Maranon Province, Peru, is An. benarrochi B, and we observed that it is highly anthropophilic. Overall, the known distribution of An. benarrochi B has been extended and includes southern Colombia, much of Peru and eastern Ecuador.

Minimal genetic differentiation of the malaria vector Nyssorhynchus darlingi associated with forest cover level in Amazonian Brazil.

The relationship between deforestation and malaria in Amazonian Brazil is complex, and a deeper understanding of this relationship is required to inform effective control measures in this region. Here, we are particularly interested in characterizing the impact of land use and land cover change on the genetics of the major regional vector of malaria, Nyssorhynchus darlingi (Root). We used nextera-tagmented, Reductively Amplified DNA (nextRAD) genotyping-by-sequencing to genotype 164 Ny. darlingi collected from 16 collection sites with divergent forest cover levels in seven municipalities in four municipality groups that span the state of Amazonas in northwestern Amazonian Brazil: São Gabriel da Cachoeira, Presidente Figueiredo, four municipalities in the area around Cruzeiro do Sul, and Lábrea. Using a dataset of 5,561 Single Nucleotide Polymorphisms (SNPs), we investigated the genetic structure of these Ny. darlingi populations with a combination of model- and non-model-based analyses. We identified weak to moderate genetic differentiation among the four municipality groups. There was no evidence for microgeographic genetic structure of Ny. darlingi among forest cover levels within the municipality groups, indicating that there may be gene flow across areas of these municipalities with different degrees of deforestation. Additionally, we conducted an environmental association analysis using two outlier detection methods to determine whether individual SNPs were associated with forest cover level without affecting overall population genetic structure. We identified 14 outlier SNPs, and investigated functions associated with their proximal genes, which could be further characterized in future studies.

Higher risk of malaria transmission outdoors than indoors by Nyssorhynchus darlingi in riverine communities in the Peruvian Amazon.

BACKGROUND: Malaria remains an important public health problem in Peru where incidence has been increasing since 2011. Of over 55,000 cases reported in 2017, Plasmodium vivax was the predominant species (76%), with P. falciparum responsible for the remaining 24%. Nyssorhynchus darlingi (previously Anopheles darlingi) is the main vector in Amazonian Peru, where hyperendemic Plasmodium transmission pockets have been found. Mazán district has pronounced spatial heterogeneity of P. vivax malaria. However, little is known about behavior, ecology or seasonal dynamics of Ny. darlingi in Mazán. This study aimed to gather baseline information about bionomics of malaria vectors and transmission risk factors in a hyperendemic malaria area of Amazonian Peru. METHODS: To assess vector biology metrics, five surveys (two in the dry and three in the rainy season), including collection of sociodemographic information, were conducted in four communities in 2016-2017 on the Napo (Urco Miraño, URC; Salvador, SAL) and Mazán Rivers (Visto Bueno, VIB; Libertad, LIB). Human-biting rate (HBR), entomological inoculation rate (EIR) and human blood index (HBI) were measured to test the hypothesis of differences in entomological indices of Ny. darlingi between watersheds. A generalized linear mixed effect model (GLMM) was constructed to model the relationship between household risk factors and the EIR. RESULTS: Nyssorhynchus darlingi comprised 95% of 7117 Anophelinae collected and its abundance was significantly higher along the Mazán River. The highest EIRs (3.03-4.54) were detected in March and June in URC, LIB and VIB, and significantly more Ny. darlingi were infected outdoors than indoors. Multivariate analysis indicated that the EIR was >12 times higher in URC compared with SAL. The HBI ranged from 0.42-0.75; humans were the most common blood source, followed by Galliformes and cows. There were dramatic differences in peak biting time and malaria incidence with similar bednet coverage in the villages. CONCLUSIONS: Nyssorhynchus darlingi is the predominant contributor to malaria transmission in the Mazán District, Peru. Malaria risk in these villages is higher in the peridomestic area, with pronounced heterogeneities between and within villages on the Mazán and the Napo Rivers. Spatiotemporal identification and quantification of the prevailing malaria transmission would provide new evidence to orient specific control measures for vulnerable or at high risk populations.

Malaria vector species in Amazonian Peru co-occur in larval habitats but have distinct larval microbial communities.

In Amazonian Peru, the primary malaria vector, Nyssorhynchus darlingi (formerly Anopheles darlingi), is difficult to target using standard vector control methods because it mainly feeds and rests outdoors. Larval source management could be a useful supplementary intervention, but to determine its feasibility, more detailed studies on the larval ecology of Ny. darlingi are essential. We conducted a multi-level study of the larval ecology of Anophelinae mosquitoes in the peri-Iquitos region of Amazonian Peru, examining the environmental characteristics of the larval habitats of four species, comparing the larval microbiota among species and habitats, and placing Ny. darlingi larval habitats in the context of spatial heterogeneity in human malaria transmission. We collected Ny. darlingi, Nyssorhynchus rangeli (formerly Anopheles rangeli), Nyssorhynchus triannulatus s.l. (formerly Anopheles triannulatus s.l.), and Nyssorhynchus sp. nr. konderi (formerly Anopheles sp. nr. konderi) from natural and artificial water bodies throughout the rainy and dry seasons. We found that, consistent with previous studies in this region and in Brazil, the presence of Ny. darlingi was significantly associated with water bodies in landscapes with more recent deforestation and lower light intensity. Nyssorhynchus darlingi presence was also significantly associated with a lower vegetation index, other Anophelinae species, and emergent vegetation. Though they were collected in the same water bodies, the microbial communities of Ny. darlingi larvae were distinct from those of Ny. rangeli and Ny. triannulatus s.l., providing evidence either for a species-specific larval microbiome or for segregation of these species in distinct microhabitats within each water body. We demonstrated that houses with more reported malaria cases were located closer to Ny. darlingi larval habitats; thus, targeted control of these sites could help ameliorate malaria risk. The co-occurrence of Ny. darlingi larvae in water bodies with other putative malaria vectors increases the potential impact of larval source management in this region.

Vector competence, vectorial capacity of Nyssorhynchus darlingi and the basic reproduction number of Plasmodium vivax in agricultural settlements in the Amazonian Region of Brazil.

BACKGROUND: Brazilian malaria control programmes successfully reduced the incidence and mortality rates from 2005 to 2016. Since 2017, increased malaria has been reported across the Amazon. Few field studies focus on the primary malaria vector in high to moderate endemic areas, Nyssorhynchus darlingi, as the key entomological component of malaria risk, and on the metrics of Plasmodium vivax propagation in Amazonian rural communities. METHODS: Human landing catch collections were carried out in 36 houses of 26 communities in five municipalities in the Brazilian states of Acre, Amazonas and Rondônia states, with API (> 30). In addition, data on the number of locally acquired symptomatic infections were employed in mathematical modelling analyses carried out to determine Ny. darlingi vector competence and vectorial capacity to P. vivax; and to calculate the basic reproduction number for P. vivax. RESULTS: Entomological indices and malaria metrics ranged among localities: prevalence of P. vivax infection in Ny. darlingi, from 0.243% in Mâncio Lima, Acre to 3.96% in Machadinho D'Oeste, Rondônia; daily human-biting rate per person from 23 ± 1.18 in Cruzeiro do Sul, Acre, to 66 ± 2.41 in Lábrea, Amazonas; vector competence from 0.00456 in São Gabriel da Cachoeira, Amazonas to 0.04764 in Mâncio Lima, Acre; vectorial capacity from 0.0836 in Mâncio Lima, to 1.5 in Machadinho D'Oeste. The estimated R0 for P. vivax (PvR0) was 3.3 in Mâncio Lima, 7.0 in Lábrea, 16.8 in Cruzeiro do Sul, 55.5 in São Gabriel da Cachoeira, and 58.7 in Machadinho D'Oeste. Correlation between P. vivax prevalence in Ny. darlingi and vector competence was non-linear whereas association between prevalence of P. vivax in mosquitoes, vectorial capacity and R0 was linear and positive. CONCLUSIONS: In spite of low vector competence of Ny. darlingi to P. vivax, parasite propagation in the human population is enhanced by the high human-biting rate, and relatively high vectorial capacity. The high PvR0 values suggest hyperendemicity in Machadinho D'Oeste and São Gabriel da Cachoeira at levels similar to those found for P. falciparum in sub-Saharan Africa regions. Mass screening for parasite reservoirs, effective anti-malarial drugs and vector control interventions will be necessary to shrinking transmission in Amazonian rural communities, Brazil.

Nyssorhynchus dunhami: bionomics and natural infection by Plasmodium falciparum and P. vivax in the Peruvian Amazon.

BACKGROUND Nyssorhynchus dunhami, a member of the Nuneztovari Complex, has been collected in Brazil, Colombia, and Peru and described as zoophilic. Although to date Ny. dunhami has not been documented to be naturally infected by Plasmodium, it is frequently misidentified as other Oswaldoi subgroup species that are local or regional malaria vectors. OBJECTIVES The current study seeks to verify the morphological identification of Nuneztovari Complex species collected in the peri-Iquitos region of Amazonian Peru, to determine their Plasmodium infection status, and to describe ecological characteristics of their larval habitats. METHODS We collected Ny. nuneztovari s.l. adults in 2011-2012, and Ny. nuneztovari s.l. larvae and adults in 2016-2017. When possible, samples were identified molecularly using cytochrome c oxidase subunit I (COI) barcode sequencing. Adult Ny. nuneztovari s.l. from 2011-2012 were tested for Plasmodium using real-time PCR. Environmental characteristics associated with Ny. nuneztovari s.l. larvae-positive water bodies were evaluated. FINDINGS We collected 590 Ny. nuneztovari s.l. adults and 116 larvae from eight villages in peri-Iquitos. Of these, 191 adults and 111 larvae were identified by COI sequencing; all were Ny. dunhami. Three Ny. dunhami were infected with P. falciparum, and one with P. vivax, all collected from one village on one night. Ny. dunhami larvae were collected from natural and artificial water bodies, and their presence was positively associated with other Anophelinae larvae and amphibians, and negatively associated with people living within 250m. MAIN CONCLUSIONS Of Nuneztovari Complex species, we identified only Ny. dunhami across multiple years in eight peri-Iquitos localities. This study is, to our knowledge, the first report of natural infection of molecularly identified Ny. dunhami with Plasmodium. We advocate the use of molecular identification methods in this region to monitor Ny. dunhami and other putative secondary malaria vectors to more precisely evaluate their importance in malaria transmission.

Decreasing proportion of Anopheles darlingi biting outdoors between long-lasting insecticidal net distributions in peri-Iquitos, Amazonian Peru.

BACKGROUND: In Loreto Department, Peru, a successful 2005-2010 malaria control programme (known as PAMAFRO) included massive distribution of long-lasting insecticidal nets (LLINs). Additional local distribution of LLINs occurred in individual villages, but not between 2012 and 2015. A 2011-2012 study of the primary regional malaria vector Anopheles darlingi detected a trend of increased exophagy compared with pre-PAMAFRO behaviour. For the present study, An. darlingi were collected in three villages in Loreto in 2013-2015 to test two hypotheses: (1) that between LLIN distributions, An. darlingi reverted to pre-intervention biting behaviour; and, (2) that there are separate sub-populations of An. darlingi in Loreto with distinct biting behaviour. RESULTS: In 2013-2015 An. darlingi were collected by human landing catch during the rainy and dry seasons in the villages of Lupuna and Cahuide. The abundance of An. darlingi varied substantially across years, villages and time periods, and there was a twofold decrease in the ratio of exophagic:endophagic An. darlingi over the study period. Unexpectedly, there was evidence of a rainy season population decline in An. darlingi. Plasmodium-infected An. darlingi were detected indoors and outdoors throughout the night, and the monthly An. darlingi human biting rate was correlated with the number of malaria cases. Using nextRAD genotyping-by-sequencing, 162 exophagic and endophagic An. darlingi collected at different times during the night were genotyped at 1021 loci. Based on model-based and non-model-based analyses, all genotyped An. darlingi belonged to a homogeneous population, with no evidence for genetic differentiation by biting location or time. CONCLUSIONS: This study identified a decreasing proportion of exophagic An. darlingi in two villages in the years between LLIN distributions. As there was no evidence for genetic differentiation between endophagic and exophagic An. darlingi, this shift in biting behaviour may be the result of behavioural plasticity in An. darlingi, which shifted towards increased exophagy due to repellence by insecticides used to impregnate LLINs and subsequently reverted to increased endophagy as the nets aged. This study highlights the need to target vector control interventions to the biting behaviour of local vectors, which, like malaria risk, shows high temporal and spatial heterogeneity.

Nyssorhynchus dunhami: bionomics and natural infection by Plasmodium falciparum and P. vivax in the Peruvian Amazon

BACKGROUND Nyssorhynchus dunhami, a member of the Nuneztovari Complex, has been collected in Brazil, Colombia, and Peru and described as zoophilic. Although to date Ny. dunhami has not been documented to be naturally infected by Plasmodium, it is frequently misidentified as other Oswaldoi subgroup species that are local or regional malaria vectors. OBJECTIVES The current study seeks to verify the morphological identification of Nuneztovari Complex species collected in the peri-Iquitos region of Amazonian Peru, to determine their Plasmodium infection status, and to describe ecological characteristics of their larval habitats. METHODS We collected Ny. nuneztovari s.l. adults in 2011-2012, and Ny. nuneztovari s.l. larvae and adults in 2016-2017. When possible, samples were identified molecularly using cytochrome c oxidase subunit I (COI) barcode sequencing. Adult Ny. nuneztovari s.l. from 2011-2012 were tested for Plasmodium using real-time PCR. Environmental characteristics associated with Ny. nuneztovari s.l. larvae-positive water bodies were evaluated. FINDINGS We collected 590 Ny. nuneztovari s.l. adults and 116 larvae from eight villages in peri-Iquitos. Of these, 191 adults and 111 larvae were identified by COI sequencing; all were Ny. dunhami. Three Ny. dunhami were infected with P. falciparum, and one with P. vivax, all collected from one village on one night. Ny. dunhami larvae were collected from natural and artificial water bodies, and their presence was positively associated with other Anophelinae larvae and amphibians, and negatively associated with people living within 250m. MAIN CONCLUSIONS Of Nuneztovari Complex species, we identified only Ny. dunhami across multiple years in eight peri-Iquitos localities. This study is, to our knowledge, the first report of natural infection of molecularly identified Ny. dunhami with Plasmodium. We advocate the use of molecular identification methods in this region to monitor Ny. dunhami and other putative secondary malaria vectors to more precisely evaluate their importance in malaria transmission.

Intensive trapping of blood-fed Anopheles darlingi in Amazonian Peru reveals unexpectedly high proportions of avian blood-meals.

Anopheles darlingi, the main malaria vector in the Neotropics, has been considered to be highly anthropophilic. However, many behavioral aspects of this species remain unknown, such as the range of blood-meal sources. Barrier screens were used to collect resting Anopheles darlingi mosquitoes from 2013 to 2015 in three riverine localities (Lupuna, Cahuide and Santa Emilia) in Amazonian Peru. Overall, the Human Blood Index (HBI) ranged from 0.58-0.87, with no significant variation among years or sites. Blood-meal analysis revealed that humans are the most common blood source, followed by avian hosts (Galliformes-chickens and turkeys), and human/Galliforme mixed-meals. The Forage Ratio and Selection Index both show a strong preference for Galliformes over humans in blood-fed mosquitoes. Our data show that 30% of An. darlingi fed on more than one host, including combinations of dogs, pigs, goats and rats. There appears to be a pattern of host choice in An. darlingi, with varying proportions of mosquitoes feeding only on humans, only on Galliformes and some taking mixed-meals of blood (human plus Galliforme), which was detected in the three sites in different years, indicating that there could be a structure to these populations based on blood-feeding preferences. Mosquito age, estimated in two localities, Lupuna and Cahuide, ranged widely between sites and years. This variation may reflect the range of local environmental factors that influence longevity or possibly potential changes in the ability of the mosquito to transmit the parasite. Of 6,204 resting An. darlingi tested for Plasmodium infection, 0.42% were infected with P. vivax. This study provides evidence for the first time of the usefulness of barrier screens for the collection of blood-fed resting mosquitoes to calculate the Human Blood Index (HBI) and other blood-meal sources in a neotropical malaria endemic setting.

Molecular operational taxonomic units of mosquitoes (Diptera: Culicidae) collected in high Andean mountain ecosystems of Antioquia, Colombia.

Accurate taxonomic identification of highland mosquito species may be complicated because of the lack of comprehensive regional morphological keys and taxonomic specialists, particularly for mosquitoes of medical or ecological importance. We applied a multi-locus approach to explore the diversity of genera/species collected, to define the Molecular Operational Taxonomic Units (MOTUs) and to perform phylogenetic clustering. Twenty MOTUs and three single sequences were revealed from 78 concatenated cox1 + ITS2 sequences, and the species name was allocated for five of these. This study provides molecular taxonomic information of culicid fauna present in high Andean mountain ecosystems in Antioquia, Colombia. However, future morphological and integrative taxonomic studies should be conducted to achieve the specific identity of all detected MOTUs.

Evidence for temporal population replacement and the signature of ecological adaptation in a major Neotropical malaria vector in Amazonian Peru.

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-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.

Changes in Genetic Diversity from Field to Laboratory During Colonization of Anopheles darlingi Root (Diptera: Culicidae).

The process of colonizing any arthropod species, including vector mosquitoes, necessarily involves adaptation to laboratory conditions. The adaptation and evolution of colonized mosquito populations needs consideration when such colonies are used as representative models for pathogen transmission dynamics. A recently established colony of Anopheles darlingi, the primary malaria vector in Amazonian South America, was tested for genetic diversity and bottleneck after 21 generations, using microsatellites. As expected, laboratory An. darlingi had fewer private and rare alleles (frequency < 0.05), decreased observed heterozygosity, and more common alleles (frequency > 0.50), but no significant evidence of a bottleneck, decrease in total alleles, or increase in inbreeding compared with field specimens (founder population). Low-moderate differentiation between field and laboratory populations was detected. With these findings, and the documented inherent differences between laboratory and field populations, results of pathogen transmission studies using this An. darlingi colony need to be interpreted cautiously.

Implications for changes in Anopheles darlingi biting behaviour in three communities in the peri-Iquitos region of Amazonian Peru.

BACKGROUND: Malaria transmission in the peri-Iquitos region of Amazonian Peru has been designated as seasonal and hypo-endemic with recently described hyper-endemic hotspots. Despite relatively recent distribution of long-lasting insecticidal bed nets (LLINs), malaria in Amazonian Peru persists and increased substantially in 2014 compared to previous years. Anopheles darlingi, identified as the main malaria vector, is known for its variable behaviour depending on locality and environment. METHODS: To evaluate vector biology metrics in relation to seasonality and malaria transmission, mosquito collections were carried out in three localities in the peri-Iquitos region, Loreto, Peru in 2011-2012. Human landing catch (HLC) collection method, Shannon (SHA) and CDC trap types were compared for effectiveness in a neotropical setting. Abundance, human biting rate and entomological inoculation rate (EIR) were measured to provide an updated view of transmission patterns post-LLIN distribution. RESULTS: HLC collected significantly more anopheline mosquitoes than SHA and CDC light traps. Anopheles darlingi was the most prevalent species in all three villages (84% overall). Biting patterns varied depending on trap type, season and village. EIR varied temporally (monthly) and spatially and the highest (2.52) occurred during the 2012 malaria outbreak in Cahuide. Unexpectedly there was a high infection rate (1.47 and 1.75) outside the normal malaria transmission season, coincident with a second local outbreak in Cahuide. The first identification of Anopheles dunhami and Anopheles oswaldoi C in Peru, using molecular markers, is also reported in this study. CONCLUSION: These data underscore the importance of HLC as the most meaningful collection method for measuring vector biology indices in this region. The highest monthly EIR provides additional evidence of seasonal transmission in riverine localities correlated with high river levels, and An. darlingi as the only contributor to transmission. The trend of an increase in outdoor-biting together with early-evening infected mosquitoes may undermine the effectiveness of LLINs as a primary malaria intervention.

A sensitive, specific and reproducible real-time polymerase chain reaction method for detection of Plasmodium vivax and Plasmodium falciparum infection in field-collected anophelines.

We describe a simple method for detection of Plasmodium vivax and Plasmodium falciparum infection in anophelines using a triplex TaqMan real-time polymerase chain reaction (PCR) assay (18S rRNA). We tested the assay on Anopheles darlingi and Anopheles stephensi colony mosquitoes fed with Plasmodium-infected blood meals and in duplicate on field collected An. darlingi. We compared the real-time PCR results of colony-infected and field collected An. darlingi, separately, to a conventional PCR method. We determined that a cytochrome b-PCR method was only 3.33% as sensitive and 93.38% as specific as our real-time PCR assay with field-collected samples. We demonstrate that this assay is sensitive, specific and reproducible.