Genetic differentiation and bottleneck effects in the malaria vectors Anopheles farauti and Anopheles punctulatus after an LLIN-based vector control program in Papua New Guinea

Implementation of long-lasting insecticide-treated net (LLIN) programs to control human malaria transmission leads to substantial reductions in the abundance of Anopheles mosquitoes, but the impact on the population genetic structure of the malaria vectors is poorly known, nor has it been investigated in Papua New Guinea, where malaria is highly endemic and where several species of Anopheles have vector roles. Here, we applied Wright's F-statistic, analysis of molecular variance, Bayesian structure analysis, and discriminant analysis of principle components to microsatellite genotype data to analyze the population genetic structure of Anopheles farauti between and within the northern and southern lowland plains and of Anopheles punctulatus within the northern plain of Papua New Guinea after such a program. Bottleneck effects in the two malaria vectors were analyzed using Luikart and Cornuet's tests of heterozygosity. A large, panmictic population of An. punctulatus pre-LLIN program diverged into two subregional populations corresponding to Madang and East Sepik provinces post-LLIN distribution and experienced a genetic bottleneck during this process. By contrast, the An. farauti population existed as two regional populations isolated by mountain ranges pre-LLIN, a genetic structure that persisted after the distribution of LLINs with no further geographic differentiation nor evidence of a genetic bottleneck. These findings show the differential response of populations of different vector species to interventions, which has implications for program sustainability and gene flow.

Genotyping of Anopheles mosquito blood meals reveals nonrandom human host selection: implications for human-to-mosquito Plasmodium falciparum transmission

BACKGROUND

Control of malaria parasite transmission can be enhanced by understanding which human demographic groups serve as the infectious reservoirs. Because vector biting can be heterogeneous, some infected individuals may contribute more to human-to-mosquito transmission than others. Infection prevalence peaks in school-age children, but it is not known how often they are fed upon. Genotypic profiling of human blood permits identification of individual humans who were bitten. The present investigation used this method to estimate which human demographic groups were most responsible for transmitting malaria parasites to Anopheles mosquitoes. It was hypothesized that school-age children contribute more than other demographic groups to human-to-mosquito malaria transmission.

METHODS

In a region of moderate-to-high malaria incidence in southeastern Malawi, randomly selected households were surveyed to collect human demographic information and blood samples. Blood-fed, female Anopheles mosquitoes were sampled indoors from the same houses. Genomic DNA from human blood samples and mosquito blood meals of human origin was genotyped using 24 microsatellite loci. The resultant genotypes were matched to identify which individual humans were sources of blood meals. In addition, Plasmodium falciparum DNA in mosquito abdomens was detected with polymerase chain reaction. The combined results were used to identify which humans were most frequently bitten, and the P. falciparum infection prevalence in mosquitoes that resulted from these blood meals.

RESULTS

Anopheles females selected human hosts non-randomly and fed on more than one human in 9% of the blood meals. Few humans contributed most of the blood meals to the Anopheles vector population. Children ≤ 5 years old were under-represented in mosquito blood meals while older males (31-75 years old) were over-represented. However, the largest number of malaria-infected blood meals was from school age children (6-15 years old).

CONCLUSIONS

The results support the hypothesis that humans aged 6-15 years are the most important demographic group contributing to the transmission of P. falciparum to the Anopheles mosquito vectors. This conclusion suggests that malaria control and prevention programmes should enhance efforts targeting school-age children and males.

Blood-feeding patterns of Anopheles vectors of human malaria in Malawi: implications for malaria transmission and effectiveness of LLIN interventions

Background

Access to human hosts by Anopheles mosquitoes is a key determinant of vectorial capacity for malaria, but it can be limited by use of long-lasting insecticidal nets (LLINs). In Malawi, pyrethroid-treated LLINs with and without the synergist piperonyl butoxide (PBO) were distributed to control malaria. This study investigated the blood-feeding patterns of malaria vectors and whether LLINs containing pyrethroid and PBO led to a reduction of human blood feeding than those containing only pyrethroids.

Methods

Mosquitoes were sampled inside houses from May 2019 through April 2020 by aspiration, pyrethrum spray catch, and light trap methods in two sites. One site (Namanolo, Balaka district) had LLINs containing only pyrethroids whereas the other (Ntaja, Machinga district) had LLINs with both pyrethroids and PBO. Anopheles species, their blood-meal host, and infection with Plasmodium falciparum were determined using PCR methods.

Results

A total of 6585 female Anopheles were sampled in 203 houses. Of these, 633 (9.6%) were blood-fed mosquitoes comprising of 279 (44.1%) Anopheles arabiensis, 103 (16.3%) Anopheles gambiae 212 (33.5), Anopheles funestus, 2 (0.3%), Anopheles parensis and 37 (5.8%) were unidentified Anopheles spp. Blood meal hosts were successfully identified for 85.5% (n = 541) of the blood-fed mosquitoes, of which 436 (81.0%) were human blood meals, 28 (5.2%) were goats, 11 (2.0%) were dogs, 60 (11.1%) were mixed goat-human blood meals, 5 (0.9%) were dog–human, and 1 was a mixed dog-goat. Human blood index (fraction of blood meals that were humans) was significantly higher in Namanolo (0.96) than Ntaja (0.89). Even though human blood index was high, goats were over-selected than humans after accounting for relative abundance of both hosts. The number of infectious Anopheles bites per person-year was 44 in Namanolo and 22 in Ntaja.

Conclusion

Although LLINs with PBO PBO may have reduced human blood feeding, access to humans was extremely high despite high LLIN ownership and usage rates in both sites. This finding could explain persistently high rates of malaria infections in Malawi. However, this study had one village for each net type, thus the observed differences may have been a result of other factors present in each village.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12936-022-04089-7.

Vector composition, abundance, biting patterns and malaria transmission intensity in Madang, Papua New Guinea: assessment after 7 years of an LLIN-based malaria control programme

Background

A malaria control programme based on distribution of long-lasting insecticidal bed nets (LLINs) and artemisinin combination therapy began in Papua New Guinea in 2009. After implementation of the programme, substantial reductions in vector abundance and malaria transmission intensity occurred. The research reported here investigated whether these reductions remained after seven years of sustained effort.

Methods

All-night (18:00 to 06:00) mosquito collections were conducted using human landing catches and barrier screen methods in four villages of Madang Province between September 2016 and March 2017. Anopheles species identification and sporozoite infection with Plasmodium vivax and Plasmodium falciparum were determined with molecular methods. Vector composition was expressed as the relative proportion of different species in villages, and vector abundance was quantified as the number of mosquitoes per barrier screen-night and per person-night. Transmission intensity was quantified as the number of sporozoite-infective vector bites per person-night.

Results

Five Anopheles species were present, but vector composition varied greatly among villages. Anopheles koliensis, a strongly anthropophilic species was the most prevalent in Bulal, Matukar and Wasab villages, constituting 63.7–73.8% of all Anopheles, but in Megiar Anopheles farauti was the most prevalent species (97.6%). Vector abundance varied among villages (ranging from 2.8 to 72.3 Anopheles per screen-night and 2.2–31.1 Anopheles per person-night), and spatially within villages. Malaria transmission intensity varied among the villages, with values ranging from 0.03 to 0.5 infective Anopheles bites per person-night. Most (54.1–75.1%) of the Anopheles bites occurred outdoors, with a substantial proportion (25.5–50.8%) occurring before 22:00.

Conclusion

The estimates of vector abundance and transmission intensity in the current study were comparable to or higher than estimates in the same villages in 2010–2012, indicating impeded programme effectiveness. Outdoor and early biting behaviours of vectors are some of the likely explanatory factors. Heterogeneity in vector composition, abundance and distribution among and within villages challenge malaria control programmes and must be considered when planning them.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12936-021-04030-4.

Nonrandom Selection and Multiple Blood Feeding of Human Hosts by Anopheles Vectors: Implications for Malaria Transmission in Papua New Guinea

Nonrandom selection and multiple blood feeding of human hosts by Anopheles mosquitoes may exacerbate malaria transmission. Both patterns of blood feeding and their relationship to malaria epidemiology were investigated in Anopheles vectors in Papua New Guinea (PNG). Blood samples from humans and mosquito blood meals were collected in villages and human genetic profiles ("fingerprints") were analyzed by genotyping 23 microsatellites and a sex-specific marker. Frequency of blood meals acquired from different humans, identified by unique genetic profiles, was fitted to Poisson and negative binomial distributions to test for nonrandom patterns of host selection. Blood meals with more than one genetic profiles were classified as mosquitoes that fed on multiple humans. The age of a person bitten by a mosquito was determined by matching the blood-meal genetic profile to the villagers' genetic profiles. Malaria infection in humans was determined by PCR test of blood samples. The results show nonrandom distribution of blood feeding among humans, with biased selection toward males and individuals aged 15-30 years. Prevalence of Plasmodium falciparum infection was higher in this age group, suggesting males in this age range could be super-spreaders of malaria parasites. The proportion of mosquitoes that fed on multiple humans ranged from 6% to 13% among villages. The patterns of host utilization observed here can amplify transmission and contribute to the persistence of malaria in PNG despite efforts to suppress it with insecticidal bed nets. Excessive feeding on males aged 15-30 years underscores the importance of targeted interventions focusing on this demographic group.

Probe-based multiplex qPCR identifies blood-meal hosts in Anopheles mosquitoes from Papua New Guinea

BACKGROUND

Determination of blood-meal hosts in blood-fed female Anopheles mosquitoes is important for evaluating vectorial capacity of vector populations and assessing effectiveness of vector control measures. Sensitive molecular methods are needed to detect traces of host blood in mosquito samples, to differentiate hosts, and to detect mixed host blood meals. This paper describes a molecular probe-based quantitative PCR for identifying blood-meal hosts in Anopheles malaria vectors from Papua New Guinea.

METHODS

TaqMan oligonucleotide probes targeting specific regions of mitochondrial or nuclear DNA of the three primary Anopheles blood-meal hosts, humans, pigs and dogs, were incorporated into a multiplex, quantitative PCR which was optimized for sensitivity and specificity.

RESULTS

Amplification of serially diluted DNA showed that the quantitative PCR detected as low as 10^-5 ng/μl of host DNA. Application to field-collected, blood-fed Anopheles showed that the quantitative PCR identified the vertebrate hosts for 89% (335/375) of mosquitoes whereas only 55% (104/188) of blood-meal samples tested in a conventional PCR were identified. Of the 104 blood-fed Anopheles that were positive in both PCR methods, 16 (15.4%) were identified as mixed blood meals by the quantitative PCR whereas only 3 (2.9%) were mixed blood meals by the conventional PCR.

CONCLUSIONS

The multiplex quantitative PCR described here is sensitive at detecting low DNA concentration and mixed host DNA in samples and useful for blood-meal analysis of field mosquitoes, in particular mixed-host blood meals.

A Microsatellite Multiplex Assay for Profiling Pig DNA in Mosquito Bloodmeals

Genetic profiling has been used to link mosquito bloodmeals to the individual humans, but this analysis has not been done for other mammalian bloodmeals. In this study, we describe a microsatellite-based method for identifying individual pigs in mosquito bloodmeals based on their unique multilocus genotypes. Eleven tetranucleotide microsatellites and a sex-specific marker were selected based on Smith-Waterman DNA sequence alignment scores from the reference genome and primers were designed with features that reduce primer dimers, promote complete adenylation, and enable fluorescent labeling of amplicons. A multiplex polymerase chain reaction (PCR) assay was optimized and validated by analyzing DNA of individual pigs from several nuclear families and breeds before it was used to analyze genomic DNA of pig-derived mosquito bloodmeals from villages of Papua New Guinea. Population analysis of the nuclear families showed high expected and observed heterozygosity. The probability of observing two unrelated or sibling individuals sharing the same genotype at a single microsatellite locus or a combination of loci was vanishingly low. Samples had unique genotypes and gender was accurately predicted. Analysis of 129 pig bloodmeals identified 19 unique genotypes, which varied greatly in frequency in the mosquito bloodmeal samples. The high allelic diversity of the microsatellite loci and low probability of false attribution of identity show that this genotyping method reliably distinguishes distantly and closely related pigs and can be used to identify individual pigs from genotyped mosquito bloodmeals.

Species abundance, composition, and nocturnal activity of female Anopheles (Diptera: Culicidae) in malaria-endemic villages of Papua New Guinea: assessment with barrier screen sampling

Background

Community composition of Anopheles mosquitoes, and their host-seeking and peridomestic behaviour, are important factors affecting malaria transmission. In this study, barrier screen sampling was used to investigate species composition, abundance, and nocturnal activity of Anopheles populations in villages of Papua New Guinea.

Methods

Mosquitoes were sampled from 6 pm to 6 am in five villages from 2012 to 2016. The barrier screens were positioned between the village houses and the perimeter of villages where cultivated and wild vegetation (“the bush”) grew thickly. Female Anopheles that rested on either village or bush side of the barrier screens, as they commuted into and out of the villages, were captured. Similarity in species composition among villages was assessed. Mosquitoes captured on village and bush sides of the barrier screens were sorted by feeding status and by hour of collection, and their numbers were compared using negative binomial generalized linear models.

Results

Females of seven Anopheles species were present in the sample. Species richness ranged from four to six species per village, but relative abundance was highly uneven within and between villages, and community composition was similar for two pairs of villages and highly dissimilar in a fifth. For most Anopheles populations, more unfed than blood-fed mosquitoes were collected from the barrier screens. More blood-fed mosquitoes were found on the side of the barrier screens facing the village and relatively more unfed ones on the bush side, suggesting commuting behaviour of unfed host-seeking females into the villages from nearby bush and commuting of blood-fed females away from villages towards the bush. For most populations, the majority of host-seeking mosquitoes arrived in the village before midnight when people were active and unprotected from the mosquitoes by bed nets.

Conclusion

The uneven distribution of Anopheles species among villages, with each site dominated by different species, even among nearby villages, emphasizes the importance of vector heterogeneity in local malaria transmission and control. Yet, for most species, nocturnal activity patterns of village entry and host seeking predominantly occurred before midnight indicating common behaviours across species and populations relative to human risk of exposure to Anopheles bites.

Electronic supplementary material

The online version of this article (10.1186/s12936-019-2742-x) contains supplementary material, which is available to authorized users.

Mosquito Behavior Change After Distribution of Bednets Results in Decreased Protection Against Malaria Exposure

Background.

Behavioral resilience in mosquitoes poses a significant challenge to mosquito control. Although behavior changes in anopheline vectors have been reported over the last decade, there are no empirical data to suggest they compromise the efficacy of vector control in reducing malaria transmission.

Methods.

In this study, we quantified human exposure to both bites and infective bites of a major malaria vector in Papua New Guinea over the course of 4 years surrounding nationwide bednet distribution. We also quantified malaria infection prevalence in the human population during the same time period.

Results.

We observed a shift in mosquito biting to earlier hours of the evening, before individuals are indoors and protected by bednets, followed by a return to preintervention biting rates. As a result, net users and non–net users experienced higher levels of transmission than before the intervention. The personal protection provided by a bednet decreased over the study period and was lowest in the adult population, who may be an important reservoir for transmission. Malaria prevalence decreased in only 1 of 3 study villages after the distribution.

Discussion.

This study highlights the necessity of validating and deploying vector control measures targeting outdoor exposure to control and eliminate malaria.

Plasticity of host selection by malaria vectors of Papua New Guinea

Background

Host selection is an important determinant of vectorial capacity because malaria transmission increases when mosquitoes feed more on humans than non-humans. Host selection also affects the outcome of long-lasting insecticidal nets (LLIN). Despite the recent nationwide implementation of LLIN-based malaria control program in Papua New Guinea (PNG), little is known about the host selection of the local Anopheles vectors. This study investigated the host selection of Anopheles vectors in PNG.

Methods

Blood-engorged mosquitoes were sampled using the barrier screen method and blood meals analyzed for vertebrate host source with PCR-amplification of the mitochondrial cytochrome b gene. Abundance of common hosts was estimated in surveys. The test of homogeneity of proportions and the Manly resource selection ratio were used to determine if hosts were selected in proportion to their abundance.

Results

Two thousand four hundred and forty blood fed Anopheles females of seven species were sampled from five villages in Madang, PNG. Of 2,142 samples tested, 2,061 (96.2%) yielded a definitive host source; all were human, pig, or dog. Hosts were not selected in proportion to their abundance, but rather were under-selected or over-selected by the mosquitoes. Four species, Anopheles farauti (sensu stricto) (s.s.), Anopheles punctulatus (s.s.), Anopheles farauti no. 4 and Anopheles longirostris, over-selected humans in villages with low LLIN usage, but over-selected pigs in villages with high LLIN usage. Anopheles koliensis consistently over-selected humans despite high LLIN usage, and Anopheles bancroftii over-selected pigs.

Conclusions

The plasticity of host selection of an Anopheles species depends on its opportunistic, anthropophilic or zoophilic behavior, and on the extent of host availability and LLIN usage where the mosquitoes forage for hosts. The high anthropophily of An. koliensis increases the likelihood of contacting the LLIN inside houses. This allows its population size to be reduced to levels insufficient to support transmission. In contrast, by feeding on alternative hosts the likelihood of the opportunistic species to contact LLIN is lower, making them difficult to control. By maintaining high population size, the proportion that feed on humans outdoors can sustain residual transmission despite high LLIN usage in the village.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-017-2038-3) contains supplementary material, which is available to authorized users.

Malaria transmission dynamics surrounding the first nationwide long-lasting insecticidal net distribution in Papua New Guinea

BACKGROUND

The major malaria vectors of Papua New Guinea exhibit heterogeneities in distribution, biting behaviour and malaria infection levels. Long-lasting, insecticide-treated nets (LLINs), distributed as part of the National Malaria Control Programme, are the primary intervention targeting malaria transmission. This study evaluated the impact of LLINs on anopheline density, species composition, feeding behaviour, and malaria transmission.

METHODS

Mosquitoes were collected by human landing catch in 11 villages from East Sepik Province and Madang Province. Mosquitoes were collected for 3 years (1 year before distribution and 2 years after), and assayed to determine mosquito species and Plasmodium spp. infection prevalence. The influence of weather conditions and the presence of people and animals on biting density was determined. Determinants of biting density and sporozoite prevalence were analysed by generalized estimating equations (GEE).

RESULTS

Mosquito biting rates and entomological inoculation rates decreased significantly after the distribution. Plasmodium falciparum and P. vivax sporozoite prevalence decreased in year 2, but increased in year 3, suggesting the likelihood of resurgence in transmission if low biting rates are not maintained. An earlier shift in the median biting time of Anopheles punctulatus and An. farauti s.s. was observed. However, this was not accompanied by an increase in the proportion of infective bites occurring before 2200 hours. A change in species composition was observed, which resulted in dominance of An. punctulatus in Dreikikir region, but a decrease in An. punctulatus in the Madang region. When controlling for village and study year, An. farauti s.s., An. koliensis and An. punctulatus were equally likely to carry P. vivax sporozoites. However, An. punctulatus was significantly more likely than An. farauti s.s. (OR 0.14; p = 0.007) or An. koliensis (OR 0.27; p < 0.001) to carry P. falciparum sporozoites.

CONCLUSIONS

LLINs had a significant impact on malaria transmission, despite exophagic and crepuscular feeding behaviours of dominant vectors. Changes in species composition and feeding behaviour were observed, but their epidemiological significance will depend on their durability over time.

Entomological Monitoring and Evaluation: Diverse Transmission Settings of ICEMR Projects Will Require Local and Regional Malaria Elimination Strategies

The unprecedented global efforts for malaria elimination in the past decade have resulted in altered vectorial systems, vector behaviors, and bionomics. These changes combined with increasingly evident heterogeneities in malaria transmission require innovative vector control strategies in addition to the established practices of long-lasting insecticidal nets and indoor residual spraying. Integrated vector management will require focal and tailored vector control to achieve malaria elimination. This switch of emphasis from universal coverage to universal coverage plus additional interventions will be reliant on improved entomological monitoring and evaluation. In 2010, the National Institutes for Allergies and Infectious Diseases (NIAID) established a network of malaria research centers termed ICEMRs (International Centers for Excellence in Malaria Research) expressly to develop this evidence base in diverse malaria endemic settings. In this article, we contrast the differing ecology and transmission settings across the ICEMR study locations. In South America, Africa, and Asia, vector biologists are already dealing with many of the issues of pushing to elimination such as highly focal transmission, proportionate increase in the importance of outdoor and crepuscular biting, vector species complexity, and "sub patent" vector transmission.

Insecticide Resistance in Areas Under Investigation by the International Centers of Excellence for Malaria Research: A Challenge for Malaria Control and Elimination

Scale-up of the main vector control interventions, residual insecticides sprayed on walls or structures and/or impregnated in bed nets, together with prompt diagnosis and effective treatment, have led to a global reduction in malaria transmission. However, resistance in vectors to almost all classes of insecticides, particularly to the synthetic pyrethroids, is posing a challenge to the recent trend of declining malaria. Ten International Centers of Excellence for Malaria Research (ICEMR) located in the most malaria-endemic regions of the world are currently addressing insecticide resistance in the main vector populations, which not only threaten hope for elimination in malaria-endemic countries but also may lead to reversal where notable reductions in malaria have been documented. This communication illustrates the current status of insecticide resistance with a focus on the countries where activities are ongoing for 9 out of the 10 ICEMRs. Most of the primary malaria vectors in the ICEMR countries exhibit insecticide resistance, albeit of varying magnitude, and spanning all mechanisms of resistance. New alternatives to the insecticides currently available are still to be fully developed for deployment. Integrated vector management principles need to be better understood and encouraged, and viable insecticide resistance management strategies need to be developed and implemented.

High throughput multiplex assay for species identification of Papua New Guinea malaria vectors: members of the Anopheles punctulatus (Diptera: Culicidae) species group

Malaria and filariasis are transmitted in the Southwest Pacific region by Anopheles punctulatus sibling species including An. punctulatus, An. koliensis, the An. farauti complex 1-8 (includes An. hinesorum [An. farauti 2], An. torresiensis [An. farauti 3]). Distinguishing these species from each other requires molecular diagnostic methods. We developed a multiplex polymerase chain reaction (PCR)-based assay specific for known species-specific nucleotide differences in the internal transcribed spacer 2 region and identified the five species most frequently implicated in transmitting disease (An. punctulatus, An. koliensis, An. farauti 1, An. hinesorum, and An. farauti 4). A set of 340 individual mosquitoes obtained from seven Papua New Guinea provinces representing a variety of habitats were analyzed by using this multiplex assay. Concordance between molecular and morphological diagnosis was 56.4% for An. punctulatus, 85.3% for An. koliensis, and 88.9% for An. farauti. Among 158 mosquitoes morphologically designated as An. farauti, 33 were re-classified by PCR as An. punctulatus, 4 as An. koliensis, 26 as An. farauti 1, 49 as An. hinesorum, and 46 as An. farauti 4. Misclassification results from variable coloration of the proboscis and overlap of An. punctulatus, An. koliensis, the An. farauti 4. This multiplex technology enables further mosquito strain identification and simultaneous detection of microbial pathogens.

Pyrethroid susceptibility in natural populations of the Anopheles punctulatus group (Diptera: Culicidae) in Papua New Guinea

The development of insecticide resistance has compromised mosquito control efforts in many parts of the world. Papua New Guinea (PNG) has a long history of dichlorodiphenyltrichloroethane (DDT) use and currently distributes pyrethroid-treated nets for malaria control. This study is the first to investigate the status of pyrethroid resistance in the Anopheles punctulatus group, the major malaria and filariasis vectors of PNG. The study used World Health Organization standard susceptibility bioassays to detect knockdown phenotypes and a novel nested polymerase chain reaction to detect the knockdown resistant (kdr) allele in these vectors. Our results show 100% susceptibility to pyrethroids in all populations surveyed and an absence of the kdr allele.