Sara A. Bickersmith

Review of genetic diversity in malaria vectors (Culicidae: Anophelinae)

We review previous studies on the genetic diversity of malaria vectors to highlight the major trends in population structure and demographic history. In doing so, we outline key information about molecular markers, sampling strategies and approaches to investigate the causes of genetic structure in Anopheles mosquitoes. Restricted gene flow due to isolation by distance and physical barriers to dispersal may explain the spatial pattern of current genetic diversity in some Anopheles species. Nonetheless, there is noteworthy disagreement among studies, perhaps due to variation in sampling methodologies, choice of molecular markers, and/or analytical approaches. More refined genealogical methods of population analysis allowing for the inclusion of the temporal component of genetic diversity facilitated the evaluation of the contribution of historical demographic processes to genetic structure. A common pattern of past unstable demography (i.e., historical fluctuation in the effective population size) by several Anopheles species, regardless of methodology (DNA markers), mosquito ecology (anthropophilic vs zoophilic), vector status (primary vs secondary) and geographical distribution, suggests that Pleistocene environmental changes were major drivers of divergence at population and species levels worldwide.

Late Pleistocene environmental changes lead to unstable demography and population divergence of Anopheles albimanus in the northern Neotropics

We investigated the historical demography of Anopheles albimanus using mosquitoes from five countries and three different DNA regions, the mitochondrial cytochrome oxidase subunit I gene (COI), the single copy nuclear white gene and the ribosomal internal transcribed spacer two (ITS2). All the molecular markers supported the taxonomic status of a single species of An. albimanus. Furthermore, agreement between the COI and the white genes suggested a scenario of Pleistocene geographic fragmentation (i.e., population contraction) and subsequent range expansion across southern Central America.

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

AbstractBackgroundMalaria 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.MethodsTo 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.ResultsHLC 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.ConclusionThese 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.

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

BackgroundThe 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.MethodsThe 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.ResultsEach 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.ConclusionCombination 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.

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

BackgroundThe 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.MethodsThe 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.ResultsAnalyses 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.ConclusionsA 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.

Molecular Taxonomy Provides New Insights into Anopheles Species of the Neotropical Arribalzagia Series

Phylogenetic analysis of partial mitochondrial cytochrome oxidase c subunit I (COI) and nuclear internal transcribed spacer 2 (ITS2) sequences were used to evaluate initial identification and to investigate phylogenetic relationships of seven Anopheles morphospecies of the Arribalzagia Series from Colombia. Phylogenetic trees recovered highly supported clades for An. punctimaculas.s., An. calderoni, An. malefactor s.l., An. neomaculipalpus, An. apicimacula s.l., An. mattogrossensis and An. peryassui. This study provides the first molecular confirmation of An. malefactorfrom 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 An. apicimacula 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.

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.

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.

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.

Molecular Taxonomy of Anopheles (Nyssorhynchus) benarrochi (Diptera: Culicidae) and Malaria Epidemiology in Southern Amazonian Peru

Anopheline specimens were collected in 2011 by human landing catch, Shannon and CDC traps from the malaria endemic localities of Santa Rosa and San Pedro in Madre de Dios Department, Peru. Most specimens were either Anopheles (Nyssorhynchus) benarrochi B or An. (Nys.) rangeli, confirmed by polymerase chain reaction-restriction fragment length polymorphism-internal transcribed spacer 2 (PCR-RFLP-ITS2) and, for selected individuals, ITS2 sequences. A few specimens from Lupuna, Loreto Department, northern Amazonian Peru, were also identified as An. benarrochi B. A statistical parsimony network using ITS2 sequences confirmed that all Peruvian An. benarrochi B analyzed were identical to those in GenBank from Putumayo, southern Colombia. Sequences of the mtDNA COI BOLD region of specimens from all three Peruvian localities were connected using a statistical parsimony network, although there were multiple mutation steps between northern and southern Peruvian sequences. A Bayesian inference of concatenated Peruvian sequences of ITS2 + COI detected a single clade with very high support for all An. benarrochi B except one individual from Lupuna that was excluded. No samples were positive for Plasmodium by CytB-PCR.