Genetic diversity of Nyssorhynchus (Anopheles) darlingi related to biting behavior in western Amazon

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

Heritability of biting time behaviours in the major African malaria vector Anopheles arabiensis

BACKGROUND

The use of insecticide-treated nets for malaria control has been associated with shifts in mosquito vector feeding behaviour including earlier and outdoor biting on humans. The relative contribution of phenotypic plasticity and heritability to these behavioural shifts is unknown. Elucidation of the mechanisms behind these shifts is crucial for anticipating impacts on vector control.

METHODS

A novel portable semi-field system (PSFS) was used to experimentally measure heritability of biting time in the malaria vector Anopheles arabiensis in Tanzania. Wild An. arabiensis from hourly collections using the human landing catch (HLC) method were grouped into one of 3 categories based on their time of capture: early (18:00-21:00), mid (22:00-04:00), and late (05:00-07:00) biting, and placed in separate holding cages. Mosquitoes were then provided with a blood meal for egg production and formation of first filial generation (F1). The F1 generation of each biting time phenotype category was reared separately, and blood fed at the same time as their mothers were captured host-seeking. The resultant eggs were used to generate the F2 generation for use in heritability assays. Heritability was assessed by releasing F2 An. arabiensis into the PSFS, recording their biting time during a human landing catch and comparing it to that of their F0 grandmothers.

RESULTS

In PSFS assays, the biting time of F2 offspring (early: 18:00-21:00, mid: 22:00-04:00 or late: 05:00-07:00) was significantly positively associated with that of their wild-caught F0 grandmothers, corresponding to an estimated heritability of 0.110 (95% CI 0.003, 0.208). F2 from early-biting F0 were more likely to bite early than F2 from mid or late-biting F0. Similarly, the probability of biting late was higher in F2 derived from mid and late-biting F0 than from early-biting F0.

CONCLUSIONS

Despite modest heritability, our results suggest that some of the variation in biting time is attributable to additive genetic variation. Selection can, therefore, act efficiently on mosquito biting times, highlighting the need for control methods that target early and outdoor biting mosquitoes.

The impact of sustained malaria control in the Loreto region of Peru: a retrospective, observational, spatially-varying interrupted time series analysis of the PAMAFRO program

Background

Although malaria control investments worldwide have resulted in dramatic declines in transmission since 2000, progress has stalled. In the Amazon, malaria resurgence has followed withdrawal of Global Fund support of the Project for Malaria Control in Andean Border Areas (PAMAFRO). We estimate intervention-specific and spatially-explicit effects of the PAMAFRO program on malaria incidence across the Loreto region of Peru, and consider the influence of the environmental risk factors in the presence of interventions.

Methods

We conducted a retrospective, observational, spatial interrupted time series analysis of malaria incidence rates among people reporting to health posts across Loreto, Peru between the first epidemiological week of January 2001 and the last epidemiological week of December 2016. Model inference is at the smallest administrative unit (district), where the weekly number of diagnosed cases of Plasmodium vivax and Plasmodium falciparum were determined by microscopy. Census data provided population at risk. We include as covariates weekly estimates of minimum temperature and cumulative precipitation in each district, as well as spatially- and temporally-lagged malaria incidence rates. Environmental data were derived from a hydrometeorological model designed for the Amazon. We used Bayesian spatiotemporal modeling techniques to estimate the impact of the PAMAFRO program, variability in environmental effects, and the role of climate anomalies on transmission after PAMAFRO withdrawal.

Findings

During the PAMAFRO program, incidence of P. vivax declined from 42.8 to 10.1 cases/1000 people/year. Incidence for P. falciparum declined from 14.3 to 2.5 cases/1000 people/year over this same period. The effects of PAMAFRO-supported interventions varied both by geography and species of malaria. Interventions were only effective in districts where interventions were also deployed in surrounding districts. Further, interventions diminished the effects of other prevailing demographic and environmental risk factors. Withdrawal of the program led to a resurgence in transmission. Increasing minimum temperatures and variability and intensity of rainfall events from 2011 onward and accompanying population displacements contributed to this resurgence.

Interpretation

Malaria control programs must consider the climate and environmental scope of interventions to maximize effectiveness. They must also ensure financial sustainability to maintain local progress and commitment to malaria prevention and elimination efforts, as well as to offset the effects of environmental change that increase transmission risk.

Funding

National Aeronautics and Space Administration, National Institutes of Health, Bill and Melinda Gates Foundation.

Host Feeding Patterns of Mansonia (Diptera, Culicidae) in Rural Settlements near Porto Velho, State of Rondonia, Brazil

Mosquito females of the genus Mansonia (Blanchard) can be a nuisance to humans and animals since they are voraciously hematophagous and feed on the blood of a variety of vertebrates. Despite their relevance, there is a lack of investigation into the blood-feeding patterns of the Mansonia species. Knowledge of the host preference is crucial in establishing the public health importance of a mosquito species and its potential to be involved in the transmission dynamics of pathogens. Species that are primarily anthropophilic can be more effective in spreading vector-borne pathogens to humans. In this study, we used an Illumina Nextera sequencing protocol and the QIIME2 workflow to assess the diversity of DNA sequences extracted in the ingested blood of mosquito species to evaluate the overall and local host choices for three species: Ma. titillans, Ma. Amazonensis, and Ma. humeralis, in rural areas alongside the Madeira River in the vicinities of the Santo Antonio Energia (SAE) reservoir in the municipality of Porto Velho, Rondônia, Western Brazil. By performing our analysis pipeline, we have found that host diversity per collection site showed a significant heterogeneity across the sample sites. In addition, in rural areas, Ma. amazonensis present a high affinity for B. taurus, Ma. humeralis shows an overall preference for C. familiaris and B. taurus, but also H. sapiens and E. caballus in urban areas, and Ma. titillans showed more opportunistic behavior in rural areas, feeding on wild animals and G. gallus, though with an overall preference for H. sapiens.

Variation in Pyrethroid Resistance Phenotypes in Anopheles darlingi in an Area with Residual Malaria Transmission: Warning of Suspected Resistance in French Guiana

Anopheles darlingi is the main vector of malaria in South America. In French Guiana, malaria transmission occurs inland and along the rivers with a regular reemergence in the lower Oyapock area. Control against malaria vectors includes indoor residual spraying of deltamethrin and the distribution of long-lasting impregnated bednets. In this context, the level of resistance to pyrethroids was monitored for 4 years using CDC bottle tests in An. darlingi populations. A loss of susceptibility to pyrethroids was recorded with 30-minute knock-down measured as low as 81%. However, no pyrethroid molecular resistance was found by sequencing a 170 base pair fragment of the S6 segment of domain II of the voltage-gated sodium channel gene. Fluctuation of resistance phenotypes may be influenced by the reintroduction of susceptible alleles from sylvatic populations or by other mechanisms of metabolic resistance.

Entomological parameters and population structure at a microgeographic scale of the main Colombian malaria vectors Anopheles albimanus and Anopheles nuneztovari

Population subdivision among several neotropical malaria vectors has been widely evaluated; however, few studies have analyzed population variation at a microgeographic scale, wherein local environmental variables may lead to population differentiation. The aim of the present study was to evaluate the genetic and geometric morphometric structure of Anopheles nuneztovari and Anopheles albimanus in endemic localities of northwestern Colombia. Genetic and phenetic structures were evaluated using microsatellites markers and wing geometric morphometrics, respectively. In addition, entomological indices of importance in transmission were calculated. Results showed that the main biting peaks of Anopheles nuneztovari were between 20:00 and 22:00, whereas Anopheles albimanus exhibited more variation in biting times among localities. Infection in An. nuneztovari by Plasmodium spp. (IR: 4.35%) and the annual entomological inoculation rate (30.31), indicated high vector exposure and local transmission risk. We did not detect Plasmodium-infected An. albimanus in this study. In general, low genetic and phenetic subdivision among the populations of both vectors was detected using a combination of phenotypic, genetic and environmental data. The results indicated high regional gene flow, although local environmental characteristics may be influencing the wing conformation differentiation and behavioral variation observed in An. albimanus. Furthermore, the population subdivision detected by microsatellite markers for both species by Bayesian genetic analysis provides a more accurate picture of the current genetic structure in comparison to previous studies. Finally, the biting behavior variation observed for both vectors among localities suggests the need for continuous malaria vector surveys covering the endemic region to implement the most effective integrated local control interventions.

Malaria Resilience in South America: Epidemiology, Vector Biology, and Immunology Insights from the Amazonian International Center of Excellence in Malaria Research Network in Peru and Brazil

The 1990s saw the rapid reemergence of malaria in Amazonia, where it remains an important public health priority in South America. The Amazonian International Center of Excellence in Malaria Research (ICEMR) was designed to take a multidisciplinary approach toward identifying novel malaria control and elimination strategies. Based on geographically and epidemiologically distinct sites in the Northeastern Peruvian and Western Brazilian Amazon regions, synergistic projects integrate malaria epidemiology, vector biology, and immunology. The Amazonian ICEMR's overarching goal is to understand how human behavior and other sociodemographic features of human reservoirs of transmission-predominantly asymptomatically parasitemic people-interact with the major Amazonian malaria vector, Nyssorhynchus (formerly Anopheles) darlingi, and with human immune responses to maintain malaria resilience and continued endemicity in a hypoendemic setting. Here, we will review Amazonian ICEMR's achievements on the synergies among malaria epidemiology, Plasmodium-vector interactions, and immune response, and how those provide a roadmap for further research, and, most importantly, point toward how to achieve malaria control and elimination in the Americas.

Anopheline diversity in urban and peri-urban malaria foci: comparison between alternative traps and seasonal effects in a city in the Western Brazilian Amazon

BACKGROUND

Continuous vector surveillance and sustainable interventions are mandatory in order to prevent anopheline proliferation (or spread to new areas) and interrupt malaria transmission. Anopheline abundance and richness were evaluated in urban and peri-urban malaria foci at a medium-sized city in the Brazilian Amazon, comparing the protected human landing catch technique (PHLC) and alternative sampling methods over different seasonal periods. Additional information was assessed for female feeding behaviour and faunal composition.

METHODS

Anophelines were sampled bimonthly in four urban and peri-urban sites in the city of Porto Velho, state of Rondônia, Brazil. The average number of captured mosquitoes was compared between an PHLC (gold standard), a tent trap (Gazetrap), and a barrier screen by means of generalized linear mixed models (GLMM), which also included season and environment (peri-urban/urban) as predictors.

RESULTS

Overall, 2962 Anopheles individuals belonging to 12 species and one complex were caught; Anopheles darlingi represented 86% of the individuals. More mosquitoes were captured in the peri-urban setting, and the urban setting was more diverse. The model estimates that significantly more anophelines were collected by PHLC than by the Screen method, and Gazetrap captured fewer individuals. However, the Screen technique yielded more blood-engorged females. The peak hours of biting activity were from 6 to 7 p.m. in urban areas and from 7 to 8 p.m. in peri-urban areas.

CONCLUSIONS

Although peri-urban settings presented a greater abundance of anophelines, Shannon and Simpson diversities were higher in urban sites. Each technique proved to be useful, depending on the purpose: PHLC was more effective in capturing the highest anopheline densities, Gazetrap caught the greatest number of species, and the barrier screen technique captured more engorged individuals. There was no seasonal effect on Anopheles assemblage structure; however, a more diverse fauna was caught in the transitional season. Biting activity was more intense from 6 to 8 p.m., with a predominance of An. darlingi.

Nyssorhynchus darlingi genome-wide studies related to microgeographic dispersion and blood-seeking behavior

BACKGROUND

In Brazil, malaria is concentrated in the Amazon Basin, where more than 99% of the annual cases are reported. The main goal of this study was to investigate the population structure and genetic association of the biting behavior of Nyssorhynchus (also known as Anopheles) darlingi, the major malaria vector in the Amazon region of Brazil, using low-coverage genomic sequencing data.

METHODS

Samples were collected in the municipality of Mâncio Lima, Acre state, Brazil between 2016 and 2017. Different approaches using genotype imputation and no gene imputation for data treatment and low-coverage sequencing genotyping were performed. After the samples were genotyped, population stratification analysis was performed.

RESULTS

Weak but statistically significant stratification signatures were identified between subpopulations separated by distances of approximately 2-3 km. Genome-wide association studies (GWAS) were performed to compare indoor/outdoor biting behavior and blood-seeking at dusk/dawn. A statistically significant association was observed between biting behavior and single nucleotide polymorphism (SNP) markers adjacent to the gene associated with cytochrome P450 (CYP) 4H14, which is associated with insecticide resistance. A statistically significant association between blood-seeking periodicity and SNP markers adjacent to genes associated with the circadian cycle was also observed.

CONCLUSION

The data presented here suggest that low-coverage whole-genome sequencing with adequate processing is a powerful tool to genetically characterize vector populations at a microgeographic scale in malaria transmission areas, as well as for use in GWAS. Female mosquitoes entering houses to take a blood meal may be related to a specific CYP4H14 allele, and female timing of blood-seeking is related to circadian rhythm genes.

Susceptibility of Field-Collected Nyssorhynchus darlingi to Plasmodium spp. in Western Amazonian Brazil

Mosquito susceptibility to Plasmodium&nbsp;spp. infection is of paramount importance for malaria occurrence and sustainable transmission. Therefore, understanding the genetic features underlying the mechanisms of susceptibility traits is pivotal to assessing malaria transmission dynamics in endemic areas. The aim of this study was to investigate the susceptibility of Nyssorhynchus darlingi-the dominant malaria vector in Brazil-to Plasmodium spp. using a reduced representation genome-sequencing protocol. The investigation was performed using a genome-wide association study (GWAS) to identify mosquito genes that are predicted to modulate the susceptibility of natural populations of the mosquito to Plasmodium infection. After applying the sequence alignment protocol, we generated the variant panel and filtered variants; leading to the detection of 202,837 SNPs in all specimens analyzed. The resulting panel was used to perform GWAS by comparing the pool of SNP variants present in Ny.&nbsp;darlingi infected with Plasmodium spp. with the pool obtained in field-collected mosquitoes with no evidence of infection by the parasite (all mosquitoes were tested separately using RT-PCR). The GWAS results for infection status showed two statistically significant variants adjacent to important genes that can be associated with susceptibility to Plasmodium infection: Cytochrome P450 (cyp450) and chitinase. This study provides relevant knowledge on malaria transmission dynamics by using a genomic approach to identify mosquito genes associated with susceptibility to Plasmodium infection in Ny. darlingi in western Amazonian Brazil.

Natural malaria infection in anophelines vectors and their incrimination in local malaria transmission in Darién, Panama

Background

More than 85% of the malaria cases in Panama occur in poor, rural and indigenous regions like Darien Province. Vector diversity, infection rate and spatial distribution are important entomological parameters of malaria transmission dynamics. Their understanding is crucial for the development of effective disease control strategies. The objective of this study was to determine the composition of Anopheles species, their natural infection rate and their geographic distribution to better understand the malaria transmission dynamics in Darién, Panama.

Methods

Anophelines mosquitoes were captured during the rainy and dry season of 2016. We selected five communities where adult anophelines were collected using CDC light-traps, and through protective human-baited traps. Detection of natural infection and Plasmodium genotype were detected via nested PCR through the amplification of ssrRNA and the circumsporozoite protein gene (csp), respectively.

Results

A total of 1,063 mosquitoes were collected mosquitoes were collected for the detection of natural infection with Plasmodium spp. Nine Anophelines species were identified, with the predominant species being: An. (Nys.) darlingi (45.0%) and An. (Nys.) albimanus (42.6%). Natural infection in An. (Nys.) albimanus with P. vivax was detected in one mosquito pool from the community Pueblo Tortuga (0.6%), three from Marraganti (1.7%), two from Bajo Chiquito (1.1%) and three pools from Alto Playona 3 (1.7%). For An. (Nys.) darlingi mosquitoes, we detected seven positive pools from the community Bajo Chiquito (4.0%), two pools from Marraganti (1.1%) and two pools from Alto Playona (1.1%). The P. vivax allelic variant VK210 was detected in infected mosquitoes.

Conclusion

The results from this study provide new information on the transmission dynamics associated with anophelines vectors in the Darién region. This is the first report of natural P. vivax infection in An. (Nys.) darlingi and its incrimination as a potential malaria vector in this region of Panama. Additional studies are necessary to expand our knowledge and determine crucial parameters in malaria transmission in Darién, which in turn will aid the National Malaria Program in attaining an adequate malaria control strategy towards malaria elimination.

Variability in human attractiveness to mosquitoes

Blood-feeding mosquitoes locate humans spatially by detecting a combination of human-derived chemical signals, including carbon dioxide, lactic acid, and other volatile organic compounds. Mosquitoes use these signals to differentiate humans from other animals. Spatial abiotic factors (e.g. humidity, heat) are also used by mosquitoes to find a host. Mosquitoes cause discomfort and harm to humans, being vectors of many pathogens. However, not all humans suffer from mosquito bites with the same frequency or intensity. Some individuals are more attractive to mosquitoes than others, and this has an important impact on the risk of infection by pathogens transmitted by these vectors, such as arboviruses and malaria parasites. Variability in human attractiveness to mosquitoes is partially due to individual characteristics in the composition and intensity in the release of mosquito attractants. The factors that determine these particularities are diverse, modestly understood and still quite controversial. Thus, this review discusses the role of pregnancy, infection with malaria parasites (Plasmodium spp.), skin microbiota, diet, and genetics in human attractiveness to mosquitoes. In brief, pregnancy and Plasmodium infection increase the host attractiveness to mosquitoes. Skin microbiota and human genetics (especially HLA alleles) modulate the production of mosquito attractants and therefore influence individual susceptibility to these insects. There is evidence pointing to a role of diet on human susceptibility to mosquitoes, with some dietary components having a bigger influence than others. In the last part of the review, other factors affecting human-mosquito interactions are debated, with a special focus on the role of mosquito genetics, pathogens and environmental factors (e.g. wind, environmental disturbances). This work highlights that individual susceptibility to mosquitoes is composed of interactions of different human-associated components, environmental factors, and mosquito characteristics. Understanding the importance of these factors, and how they interact with each other, is essential for the development of better mosquito control strategies and studies focused on infectious disease dynamics.

•

Individual human attractiveness to mosquitoes is highly variable.

•

Mosquito attractants released into the air vary from person to person.

•

Variation in attractiveness to mosquitoes alters the risk of mosquito-borne infections.

•

Pregnancy, malaria infection, skin microbiota and genetic factors alter the release of mosquito attractants.

•

Environment and mosquito-related factors affect human–mosquito interactions.

The influence of anthropogenic habitat fragmentation on the genetic structure and diversity of the malaria vector Anopheles cruzii (Diptera: Culicidae)

Fragmentation of natural environments as a result of human interference has been associated with a decrease in species richness and increase in abundance of a few species that have adapted to these environments. The Brazilian Atlantic Forest, which has been undergoing an intense process of fragmentation and deforestation caused by human-made changes to the environment, is an important hotspot for malaria transmission. The main vector of simian and human malaria in this biome is the mosquito Anopheles cruzii. Anthropogenic processes reduce the availability of natural resources at the tree canopies, An. cruzii primary habitat. As a consequence, An. cruzii moves to the border of the Atlantic Forest nearing urban areas seeking resources, increasing their contact with humans in the process. We hypothesized that different levels of anthropogenic changes to the environment can be an important factor in driving the genetic structure and diversity in An. cruzii populations. Five different hypotheses using a cross-sectional and a longitudinal design were tested to assess genetic structure in sympatric An. cruzii populations and microevolutionary processes driving these populations. Single nucleotide polymorphisms were used to assess microgeographic genetic structure in An. cruzii populations in a low-endemicity area in the city of São Paulo, Brazil. Our results show an overall weak genetic structure among the populations, indicating a high gene flow system. However, our results also pointed to the presence of significant genetic structure between sympatric An. cruzii populations collected at ground and tree-canopy habitats in the urban environment and higher genetic variation in the ground-level population. This indicates that anthropogenic modifications leading to habitat fragmentation and a higher genetic diversity and structure in ground-level populations could be driving the behavior of An. cruzii, ultimately increasing its contact with humans. Understanding how anthropogenic changes in natural areas affect An. cruzii is essential for the development of more effective mosquito control strategies and, on a broader scale, for malaria-elimination efforts in the Brazilian Atlantic Forest.

Evidence for family-level variation of phenotypic traits in response to temperature of Brazilian Nyssorhynchus darlingi

BACKGROUND

Nyssorhynchus darlingi (also known as Anopheles darlingi) is the primary malaria vector in the Amazon River Basin. In Brazil, analysis of single nucleotide polymorphisms (SNPs) previously detected three major population clusters, and a common garden experiment in a laboratory setting revealed significant population variation in life history traits. Increasing temperatures and local level variation can affect life history traits, i.e. adult longevity, that alter vectorial capacity with implications for malaria transmission in Ny. darlingi.

METHODS

We investigated the population structure of Ny. darlingi from 7 localities across Brazil utilizing SNPs and compared them to a comprehensive Ny. darlingi catalog. To test the effects of local level variation on life history traits, we reared F1 progeny from the 7 localities at three constant temperatures (20, 24 and 28 °C), measuring key life history traits (larval development, food-starved adult lifespan, adult size and daily survival).

RESULTS

Using nextRAD genotyping-by-sequencing, 93 of the field-collected Ny. darlingi were genotyped at 33,759 loci. Results revealed three populations (K = 3), congruent with major biomes (Amazonia, Cerrado and Mata Atlântica), with greater FST values between biomes than within. In the life history experiments, increasing temperature reduced larval development time, adult lifespan, and wing length in all localities. The variation of family responses for all traits within four localities of the Amazonia biome was significant (ANOVA, P < 0.05). Individual families within localities revealed a range of responses as temperature increased, for larval development, adult lifespan, wing length and survival time.

CONCLUSIONS

SNP analysis of several Brazilian localities provided results in support of a previous study wherein populations of Ny. darlingi were clustered by three major Brazilian biomes. Our laboratory results of temperature effects demonstrated that population variation in life history traits of Ny. darlingi exists at the local level, supporting previous research demonstrating the high plasticity of this species. Understanding this plasticity and inherent variation between families of Ny. darlingi at the local level should be considered when deploying intervention strategies and may improve the likelihood of successful malaria elimination in South America.