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

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

Phylogeography of the neotropical Anopheles triannulatus complex (Diptera: Culicidae) supports deep structure and complex patterns

Abstract Background The molecular phylogenetic relationships and population structure of the species of the Anopheles triannulatus complex: Anopheles triannulatus s.s., Anopheles halophylus and the putative species Anopheles triannulatus C were investigated. Methods\ud \ud The mitochondrial COI gene, the nuclear white gene and rDNA ITS2 of samples that include the known geographic distribution of these taxa were analyzed. Phylogenetic analyses were performed using Bayesian inference, Maximum parsimony and Maximum likelihood approaches. Results Each data set analyzed septely yielded a different topology but none provided evidence for the seption of An. halophylus and An. triannulatus C, consistent with the hypothesis that the two are undergoing incipient speciation. The phylogenetic analyses of the white gene found three main clades, whereas the statistical parsimony network detected only a single metapopulation of Anopheles triannulatus s.l. Seven COI lineages were detected by phylogenetic and network analysis. In contrast, the network, but not the phylogenetic analyses, strongly supported three ITS2 groups. Combined data analyses provided the best resolution of the trees, with two major clades, Amazonian (clade I) and trans-Andean + Amazon Delta (clade II). Clade I consists of multiple subclades: An. halophylus + An. triannulatus C; trans-Andean Venezuela; central Amazonia + central Bolivia; Atlantic coastal lowland; and Amazon delta. Clade II includes three subclades: Panama; cis-Andean Colombia; and cis-Venezuela. The Amazon delta specimens are in both clades, likely indicating local sympatry. Spatial and molecular variance analyses detected nine groups, corroborating some of subclades obtained in the combined data analysis. Conclusion Combination of the three molecular markers provided the best resolution for differentiation within An. triannulatus s.s. and An. halophylus and C. The latest two species seem to be very closely related and the analyses performed were not conclusive regarding species differentiation. Further studies including new molecular markers would be desirable to solve this species status question. Besides, results of the study indicate a trans-Andean origin for An. triannulatus s.l. The potential implications for malaria epidemiology remain to be investigated.This project has been financially supported by NIH grant AI ROI 54139–02 to JE Conn, PAPES IV Fiocruz-CNPq to RLO, and from the Fundação de Amparo à Pesquisa do Estado de São Paulo, FAPESP (Grant 2011/20397-7 to MAMS), CNPq (BPP no. 301666/2011-3 to MAMS).This research was performed under a Memorandum of Understanding between the Walter Reed Army Institute of Research and the Smithsonian Institution, with institutional support provided by both organizations. The material to be published reflects the views of the authors and should not be construed to represent those of the Department of the Army or the Department of Defense

Phylogeography of the neotropical Anopheles triannulatus complex (Diptera: Culicidae) supports deep structure and complex patterns

Abstract Background The molecular phylogenetic relationships and population structure of the species of the Anopheles triannulatus complex: Anopheles triannulatus s.s., Anopheles halophylus and the putative species Anopheles triannulatus C were investigated. Methods\ud \ud The mitochondrial COI gene, the nuclear white gene and rDNA ITS2 of samples that include the known geographic distribution of these taxa were analyzed. Phylogenetic analyses were performed using Bayesian inference, Maximum parsimony and Maximum likelihood approaches. Results Each data set analyzed septely yielded a different topology but none provided evidence for the seption of An. halophylus and An. triannulatus C, consistent with the hypothesis that the two are undergoing incipient speciation. The phylogenetic analyses of the white gene found three main clades, whereas the statistical parsimony network detected only a single metapopulation of Anopheles triannulatus s.l. Seven COI lineages were detected by phylogenetic and network analysis. In contrast, the network, but not the phylogenetic analyses, strongly supported three ITS2 groups. Combined data analyses provided the best resolution of the trees, with two major clades, Amazonian (clade I) and trans-Andean + Amazon Delta (clade II). Clade I consists of multiple subclades: An. halophylus + An. triannulatus C; trans-Andean Venezuela; central Amazonia + central Bolivia; Atlantic coastal lowland; and Amazon delta. Clade II includes three subclades: Panama; cis-Andean Colombia; and cis-Venezuela. The Amazon delta specimens are in both clades, likely indicating local sympatry. Spatial and molecular variance analyses detected nine groups, corroborating some of subclades obtained in the combined data analysis. Conclusion Combination of the three molecular markers provided the best resolution for differentiation within An. triannulatus s.s. and An. halophylus and C. The latest two species seem to be very closely related and the analyses performed were not conclusive regarding species differentiation. Further studies including new molecular markers would be desirable to solve this species status question. Besides, results of the study indicate a trans-Andean origin for An. triannulatus s.l. The potential implications for malaria epidemiology remain to be investigated.This project has been financially supported by NIH grant AI ROI 54139–02 to JE Conn, PAPES IV Fiocruz-CNPq to RLO, and from the Fundação de Amparo à Pesquisa do Estado de São Paulo, FAPESP (Grant 2011/20397-7 to MAMS), CNPq (BPP no. 301666/2011-3 to MAMS).This research was performed under a Memorandum of Understanding between the Walter Reed Army Institute of Research and the Smithsonian Institution, with institutional support provided by both organizations. The material to be published reflects the views of the authors and should not be construed to represent those of the Department of the Army or the Department of Defense

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

A sensitive, specific and reproducible real-time polymerase chain reaction method for detection of Plasmodium vivaxandPlasmodium 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

Molecular Taxonomy Provides New Insights into <i>Anopheles</i> Species of the Neotropical Arribalzagia Series

<div><p>Phylogenetic analysis of partial mitochondrial cytochrome oxidase c subunit I (<i>COI</i>) and nuclear internal transcribed spacer 2 (ITS2) sequences were used to evaluate initial identification and to investigate phylogenetic relationships of seven <i>Anopheles</i> morphospecies of the Arribalzagia Series from Colombia. Phylogenetic trees recovered highly supported clades for <i>An</i>. <i>punctimacula</i>s.s., <i>An</i>. <i>calderoni</i>, <i>An</i>. <i>malefactor</i> s.l., <i>An</i>. <i>neomaculipalpus</i>, <i>An</i>. <i>apicimacula</i> s.l., <i>An</i>. <i>mattogrossensis</i> and <i>An</i>. <i>peryassui</i>. This study provides the first molecular confirmation of <i>An</i>. <i>malefactor</i>from Colombia and discovered conflicting patterns of divergence for the molecular markers among specimens from northeast and northern Colombia suggesting the presence of two previously unrecognized Molecular Operational Taxonomic Units (MOTUs). Furthermore, two highly differentiated <i>An</i>. <i>apicimacula</i> MOTUs previously found in Panama were detected. Overall, the combined molecular dataset facilitated the detection of known and new Colombian evolutionary lineages, and constitutes the baseline for future research on their bionomics, ecology and potential role as malaria vectors.</p></div

Bayesian topology of ITS2 sequences for species of the Arribalzagia Series.

<p>Numbers on each branch represent posterior probabilities. <i>Anopheles pseudopunctipennis</i> was the outgroup.</p