Loss of genetic diversity in Culex quinquefasciatus targeted by a lymphatic filariasis vector control program in Recife, Brazil

Assessment of wing geometric morphometrics of urban Culex quinquefasciatus (Diptera: Culicidae) populations

Culex quinquefasciatus is a cosmopolitan species distributed throughout tropical and subtropical areas of the world. The species is of great epidemiological importance as it is responsible for vectoring the causative agent of lymphatic filariasis and several arboviruses, including West Nile virus. Wing geometric morphometrics has been widely used to assess phenotypic variations in mosquito species. Here, we hypothesize that Cx. quinquefasciatus populations in urban parks in the city of São Paulo, Brazil, have been subjected to anthropogenic selective pressures that are responsible for driving their ecology and behavior. Mosquitoes were collected by CDC traps in five municipal parks in the city of São Paulo. Eighteen anatomical landmark coordinates on each female right wing were digitized. Canonical variate analysis, wireframe graphs, cross-validated reclassification tests and the neighbor-joining method were used to assess phenotypical dissimilarity in wing shape between populations. Centroid size was calculated to assess differences in wing size between populations, which can result from different environmental conditions during immature mosquito development. Moderately heterogeneous wing shape and wing size patterns were found in the populations analyzed, indicating that selective pressures in the urban environment are affecting the wing patterns of Cx. quinquefasciatus populations in the city of São Paulo, Brazil.

Development of an urban molecular xenomonitoring system for lymphatic filariasis in the Recife Metropolitan Region, Brazil

Introduction

Molecular xenomonitoring (MX)—pathogen detection in the mosquito rather than human—is a promising tool for lymphatic filariasis (LF) surveillance. In the Recife Metropolitan Region (RMR), the last LF focus in Brazil, Culex quinquefasciatus mosquitoes have been implicated in transmitting Wuchereria bancrofti parasites. This paper presents findings on the ideal mosquito collection method, mosquito dispersion, W. bancrofti infection in mosquitoes and W. bancrofti antigen in humans to aid MX development.

Methods

Experiments occurred within two densely populated urban areas of Olinda, RMR, in July and August 2015. U.S. Centers for Disease Control and Prevention (CDC) light traps were compared to battery-powered aspirators as collection methods, and mosquito dispersion was measured by mosquito mark release recapture (MMRR). Female Cx. quinquefasciatus were tested by PCR for W. bancrofti infection, and study area residents were screened by rapid tests for W. bancrofti antigen.

Results

Aspirators caught 2.6 times more total Cx. quinquefasciatus, including 38 times more blood-fed and 5 times more gravid stages, than CDC light traps. They also collected 123 times more Aedes aegypti. Of the 9,644 marked mosquitoes released, only ten (0.01%) were recaptured, nine of which were < 50m (34.8m median, 85.4m maximum) from the release point. Of 9,169 unmarked mosquitoes captured in the MMR, 38.3% were unfed, 48.8% blood-fed, 5.5% semi-gravid, and 7.3% gravid. PCR on 182 pools (1,556 mosquitoes) found no evidence of W. bancrofti infection in Cx. quinquefasciatus. Rapid tests on 110 of 111 eligible residents were all negative for W. bancrofti antigen.

Conclusions

Aspirators were more effective than CDC light traps at capturing Ae. aegypti and all but unfed stages of Cx. quinquefasciatus. Female Cx. quinquefasciatus traveled short (< 86m) distances in this urban area. Lack of evidence for W. bancrofti infection in mosquitoes and antigen in humans in these fine-scale studies does not indicate that LF transmission has ceased in the RMR. A MX surveillance system should consider vector-specific collection methods, mosquito dispersion, and spatial scale but also local context, environmental factors such as sanitation, and host factors such as infection prevalence and treatment history.

Lymphatic filariasis (LF) is a parasitic disease transmitted by mosquitoes, and can cause elephantiasis. It is the world’s leading cause of disability due to infectious diseases, affects over 120 million people globally, and is scheduled for global elimination via mass drug administration (MDA) and mosquito control. Molecular xenomonitoring (MX) is a process of screening mosquitoes—not humans—for parasites to estimate whether they are circulating in human populations. MX is especially useful during and following MDA, when new case detection becomes difficult, but is challenging to design and conduct in cities. Using two study sites in the Recife Metropolitan Region, Brazil, we investigated two crucial questions for urban MX development—“What is the best operationally feasible tool to catch adult mosquitoes?” and “How far do mosquitoes disperse in cities?”—in order to determine placement of future surveillance sites. We also screened a proportion of mosquitoes and all eligible residents from the study sites for LF infection. We determined that handheld battery powered aspirators were the best mosquito collection tool; that mosquitoes flew no more than about 85m; and—in this small sample of mosquitoes and very small sample of humans—there was no evidence of LF infection in mosquitoes or study area residents.

Evolution of life in urban environments

Our planet is an increasingly urbanized landscape, with over half of the human population residing in cities. Despite advances in urban ecology, we do not adequately understand how urbanization affects the evolution of organisms, nor how this evolution may affect ecosystems and human health. Here, we review evidence for the effects of urbanization on the evolution of microbes, plants, and animals that inhabit cities. Urbanization affects adaptive and nonadaptive evolutionary processes that shape the genetic diversity within and between populations. Rapid adaptation has facilitated the success of some native species in urban areas, but it has also allowed human pests and disease to spread more rapidly. The nascent field of urban evolution brings together efforts to understand evolution in response to environmental change while developing new hypotheses concerning adaptation to urban infrastructure and human socioeconomic activity. The next generation of research on urban evolution will provide critical insight into the importance of evolution for sustainable interactions between humans and our city environments.

Local selection in the presence of high levels of gene flow: Evidence of heterogeneous insecticide selection pressure across Ugandan Culex quinquefasciatus populations

BACKGROUND

Culex quinquefasciatus collected in Uganda, where no vector control interventions directly targeting this species have been conducted, was used as a model to determine if it is possible to detect heterogeneities in selection pressure driven by insecticide application targeting other insect species.

METHODOLOGY/PRINCIPAL FINDINGS

Population genetic structure was assessed through microsatellite analysis, and the impact of insecticide pressure by genotyping two target-site mutations, Vgsc-1014F of the voltage-gated sodium channel target of pyrethroid and DDT insecticides, and Ace1-119S of the acetylcholinesterase gene, target of carbamate and organophosphate insecticides. No significant differences in genetic diversity were observed among populations by microsatellite markers with HE ranging from 0.597 to 0.612 and low, but significant, genetic differentiation among populations (FST = 0.019, P = 0.001). By contrast, the insecticide-resistance markers display heterogeneous allelic distributions with significant differences detected between Central Ugandan (urban) populations relative to Eastern and Southwestern (rural) populations. In the central region, a frequency of 62% for Vgsc-1014F, and 32% for the Ace1-119S resistant allele were observed. Conversely, in both Eastern and Southwestern regions the Vgsc-1014F alleles were close to fixation, whilst Ace1-119S allele frequency was 12% (although frequencies may be underestimated due to copy number variation at both loci).

CONCLUSIONS/SIGNIFICANCE

Taken together, the microsatellite and both insecticide resistance target-site markers provide evidence that in the face of intense gene flow among populations, disjunction in resistance frequencies arise due to intense local selection pressures despite an absence of insecticidal control interventions targeting Culex.

Morphometric Wing Characters as a Tool for Mosquito Identification

Mosquitoes are responsible for the transmission of important infectious diseases, causing millions of deaths every year and endangering approximately 3 billion people around the world. As such, precise identification of mosquito species is crucial for an understanding of epidemiological patterns of disease transmission. Currently, the most common method of mosquito identification relies on morphological taxonomic keys, which do not always distinguish cryptic species. However, wing geometric morphometrics is a promising tool for the identification of vector mosquitoes, sibling and cryptic species included. This study therefore sought to accurately identify mosquito species from the three most epidemiologically important mosquito genera using wing morphometrics. Twelve mosquito species from three epidemiologically important genera (Aedes, Anopheles and Culex) were collected and identified by taxonomic keys. Next, the right wing of each adult female mosquito was removed and photographed, and the coordinates of eighteen digitized landmarks at the intersections of wing veins were collected. The allometric influence was assessed, and canonical variate analysis and thin-plate splines were used for species identification. Cross-validated reclassification tests were performed for each individual, and a Neighbor Joining tree was constructed to illustrate species segregation patterns. The analyses were carried out and the graphs plotted with TpsUtil 1.29, TpsRelw 1.39, MorphoJ 1.02 and Past 2.17c. Canonical variate analysis for Aedes, Anopheles and Culex genera showed three clear clusters in morphospace, correctly distinguishing the three mosquito genera, and pairwise cross-validated reclassification resulted in at least 99% accuracy; subgenera were also identified correctly with a mean accuracy of 96%, and in 88 of the 132 possible comparisons, species were identified with 100% accuracy after the data was subjected to reclassification. Our results showed that Aedes, Culex and Anopheles were correctly distinguished by wing shape. For the lower hierarchical levels (subgenera and species), wing geometric morphometrics was also efficient, resulting in high reclassification scores.

Islands and stepping-stones: comparative population structure of Anopheles gambiae sensu stricto and Anopheles arabiensis in Tanzania and implications for the spread of insecticide resistance

Population genetic structures of the two major malaria vectors Anopheles gambiae s.s. and An. arabiensis, differ markedly across Sub-Saharan Africa, which could reflect differences in historical demographies or in contemporary gene flow. Elucidation of the degree and cause of population structure is important for predicting the spread of genetic traits such as insecticide resistance genes or artificially engineered genes. Here the population genetics of An. gambiae s.s. and An. arabiensis in the central, eastern and island regions of Tanzania were compared. Microsatellite markers were screened in 33 collections of female An. gambiae s.l., originating from 22 geographical locations, four of which were sampled in two or three years between 2008 and 2010. An. gambiae were sampled from six sites, An. arabiensis from 14 sites, and both species from two sites, with an additional colonised insectary sample of each species. Frequencies of the knock-down resistance (kdr) alleles 1014S and 1014F were also determined. An. gambiae exhibited relatively high genetic differentiation (average pairwise FST = 0.131), significant even between nearby samples, but without clear geographical patterning. In contrast, An. arabiensis exhibited limited differentiation (average FST = 0.015), but strong isolation-by-distance (Mantel test r = 0.46, p = 0.0008). Most time-series samples of An. arabiensis were homogeneous, suggesting general temporal stability of the genetic structure. An. gambiae populations from Dar es Salaam and Bagamoyo were found to have high frequencies of kdr 1014S (around 70%), with almost 50% homozygote but was at much lower frequency on Unguja Island, with no. An. gambiae population genetic differentiation was consistent with an island model of genetic structuring with highly restricted gene flow, contrary to An. arabiensis which was consistent with a stepping-stone model of extensive, but geographically-restricted gene flow.

Population genetics of neotropical Culex quinquefasciatus (Diptera: Culicidae)

Background

Culex quinquefasciatus mosquitoes can be found in almost every major city of Brazil and are vectors of filariasis and several arboviruses. Microsatellite markers have been widely used to uncover the genetic structure of various groups of insect populations. The aim of this study was to glimpse the genetic structure of Cx. quinquefasciatus in Brazil.

Methods

Nine populations were sampled across Brazil (one of them from a laboratory colony - COL) and another one from Argentina and process regarding the variability of six microsatellite loci.

Results

The analyzed loci revealed moderate population genetic structure (mean F_st = 0.12). Dendrograms of genetic distances evidenced two major population clusters, respectively corresponding to the northern and southern populations. The hybrid population Cx. pipiens/quinquefasciatus (from La Plata, Argentina) and the colony population fell outside the major clusters. Those clusters were substructured and there was a significant correlation between genetic and geographic distances and environmental variables (r = 0.51; p > 0.001 and r = 0.46; p > 0.004).

Conclusions

Multilocus cluster Bayesian analysis confirmed that populations are mutually distinct, and the set of results point to genetic differences among populations. The presumable low gene flow among them may be due to the large geographic distances (>1000 km) and to the environmental heterogeneity of the sampled areas. The genetic structure observed in this study may lead to the best understanding of Cx. quinquefasciatus demographical diversity as well as their genetic variations patterns in Brazil so far unknown.

Mitochondrial DNA analyses reveal low genetic diversity in Culex quinquefasciatus from residential areas in Malaysia

The present study explored the intraspecific genetic diversity, dispersal patterns and phylogeographic relationships of Culex quinquefasciatus Say (Diptera: Culicidae) in Malaysia using reference data available in GenBank in order to reveal this species' phylogenetic relationships. A statistical parsimony network of 70 taxa aligned as 624 characters of the cytochrome c oxidase subunit I (COI) gene and 685 characters of the cytochrome c oxidase subunit II (COII) gene revealed three haplotypes (A1-A3) and four haplotypes (B1-B4), respectively. The concatenated sequences of both COI and COII genes with a total of 1309 characters revealed seven haplotypes (AB1-AB7). Analysis using tcs indicated that haplotype AB1 was the common ancestor and the most widespread haplotype in Malaysia. The genetic distance based on concatenated sequences of both COI and COII genes ranged from 0.00076 to 0.00229. Sequence alignment of Cx. quinquefasciatus from Malaysia and other countries revealed four haplotypes (AA1-AA4) by the COI gene and nine haplotypes (BB1-BB9) by the COII gene. Phylogenetic analyses demonstrated that Malaysian Cx. quinquefasciatus share the same genetic lineage as East African and Asian Cx. quinquefasciatus. This study has inferred the genetic lineages, dispersal patterns and hypothetical ancestral genotypes of Cx. quinquefasciatus.

Susceptibility status of Culex quinquefasciatus (Diptera: Culicidae) populations to the chemical insecticide temephos in Pernambuco, Brazil

BACKGROUND

Culex quinquefasciatus is the vector of many agents of human diseases, including Wuchereria bancrofti, the parasite that causes bancroftian filariasis, an endemic disease in Pernambuco State, Brazil. Although temephos is not currently used to control C. quinquefasciatus, the species might be under a selection process from incidental exposure to this compound. This study aimed to evaluate the susceptibility status of C. quinquefasciatus to temephos, using bioassays, and to investigate its putative resistance mechanisms through biochemical assays and screening of the G119S mutation in the acetylcholinesterase gene, which is associated with organophosphate resistance, carried out by PCR and sequencing.

RESULTS

The results showed that only mosquitoes from Santa Cruz do Capibaribe (SC) had an alteration in their susceptibility status (RR = 7.2-fold), while the other populations were all susceptible to the insecticide. Biochemical assays showed increased activity for all esterases in SC, as well as evidence of acetylcholinesterase insensitivity. The G119S mutation was detected in this population with a frequency of 0.11, but it was not found in the remaining populations.

CONCLUSION

These data show that mechanisms of temephos resistance have been selected in natural C. quinquefasciatus populations from Pernambuco, which could undermine future control actions.

Large fluctuations in the effective population size of the malaria mosquito Anopheles gambiae s.s. during vector control cycle

On Bioko Island, Equatorial Guinea, indoor residual spraying (IRS) has been part of the Bioko Island Malaria Control Project since early 2004. Despite success in reducing childhood infections, areas of high transmission remain on the island. We therefore examined fluctuations in the effective population size (N e ) of the malaria vector Anopheles gambiae in an area of persistent high transmission over two spray rounds. We analyzed data for 13 microsatellite loci from 791 An. gambiae specimens collected at six time points in 2009 and 2010 and reconstructed the demographic history of the population during this period using approximate Bayesian computation (ABC). Our analysis shows that IRS rounds have a large impact on N e , reducing it by 65%-92% from prespray round N e . More importantly, our analysis shows that after 3-5 months, the An. gambiae population rebounded by 2818% compared shortly following the spray round. Our study underscores the importance of adequate spray round frequency to provide continuous suppression of mosquito populations and that increased spray round frequency should substantially improve the efficacy of IRS campaigns. It also demonstrates the ability of ABC to reconstruct a detailed demographic history across only a few tens of generations in a large population.

Sequential loss of genetic variation in flea beetle Agasicles hygrophila (Coleoptera: Chrysomelidae) following introduction into China

The flea beetle (Agasicles hygrophila) was imported to Florida, USA and then introduced from Florida into China in 1987 as a biological control agent for the invasive plant alligator weed (Alternanthera philoxeroides). The initial beetle population was subsequently used for sequential introductions in other areas of China, but little is known about the genetic consequences of the introductions. In this study, the genetic diversity and population structure of five beetle populations, the source Florida population, three intentionally introduced China populations and one accidentally dispersed China population, were examined using amplified fragment length polymorphisms. The results showed a clear pattern of decreasing genetic diversity with the sequential introductions. The diversity was highest in the Florida population followed by the first introduction to Chongqing and then in Kunming and Fuzhou. The lowest diversity was found in the accidentally dispersed Guangzhou population that was first recorded in 1996. Both loci parameters and Nei's genetic diversity showed a high variation among these populations. Genetic differentiation among populations was further verified by the GST statistic (0.136-0.432). Beetles in Kunming had the highest gene flow with those in Guangzhou, and therefore lowest differentiation and closest genetic distance. These data show that sequential introduction influenced the genetic diversity of populations in China. Genetic diversity should be considered in planning introduction and long-term maintenance of populations.

Population genetic and admixture analyses of Culex pipiens complex (Diptera: Culicidae) populations in California, United States

Microsatellite markers were used to genetically characterize 19 Culex pipiens complex populations from California. Two populations showed characteristics of earlier genetic bottlenecks. The overall F_ST value and a neighbor-joining tree suggested moderate amounts of genetic differentiation. Analyses using Structure indicated K = 4 genetic clusters: Cx. pipiens form pipiens L., Cx. quinquefasciatus Say, Cx. pipiens form molestus Forskäl, and a group of genetically similar individuals of hybrid origin. A Discriminant Analysis of Principal Components indicated that the latter group is a mixture of the other three taxa, with form pipiens and form molestus contributing somewhat more ancestry than Cx. quinquefasciatus. Characterization of 56 morphologically autogenous individuals classified most as Cx. pipiens form molestus, and none as Cx. pipiens form pipiens or Cx. quinquefasciatus. Comparison of California microsatellite data with those of Cx. pipiens pallens Coquillett from Japan indicated the latter does not contribute significantly to genotypes in California.

The effective population size of malaria mosquitoes: large impact of vector control

Malaria vectors in sub-Saharan Africa have proven themselves very difficult adversaries in the global struggle against malaria. Decades of anti-vector interventions have yielded mixed results--with successful reductions in transmission in some areas and limited impacts in others. These varying successes can be ascribed to a lack of universally effective vector control tools, as well as the development of insecticide resistance in mosquito populations. Understanding the impact of vector control on mosquito populations is crucial for planning new interventions and evaluating existing ones. However, estimates of population size changes in response to control efforts are often inaccurate because of limitations and biases in collection methods. Attempts to evaluate the impact of vector control on mosquito effective population size (N(e)) have produced inconclusive results thus far. Therefore, we obtained data for 13-15 microsatellite markers for more than 1,500 mosquitoes representing multiple time points for seven populations of three important vector species--Anopheles gambiae, An. melas, and An. moucheti--in Equatorial Guinea. These populations were exposed to indoor residual spraying or long-lasting insecticidal nets in recent years. For comparison, we also analyzed data from two populations that have no history of organized vector control. We used Approximate Bayesian Computation to reconstruct their demographic history, allowing us to evaluate the impact of these interventions on the effective population size. In six of the seven study populations, vector control had a dramatic impact on the effective population size, reducing N(e) between 55%-87%, the exception being a single An. melas population. In contrast, the two negative control populations did not experience a reduction in effective population size. This study is the first to conclusively link anti-vector intervention programs in Africa to sharply reduced effective population sizes of malaria vectors.