A targeted amplicon sequencing panel to simultaneously identify mosquito species and Plasmodium presence across the entire Anopheles genus

Anopheles is a diverse genus of mosquitoes comprising over 500 described species, including all known human malaria vectors. While a limited number of key vector species have been studied in detail, the goal of malaria elimination calls for surveillance of all potential vector species. Here, we develop a multilocus amplicon sequencing approach that targets 62 highly variable loci in the Anopheles genome and two conserved loci in the Plasmodium mitochondrion, simultaneously revealing both the mosquito species and whether that mosquito carries malaria parasites. We also develop a cheap, nondestructive, and high-throughput DNA extraction workflow that provides template DNA from single mosquitoes for the multiplex PCR, which means specimens producing unexpected results can be returned to for morphological examination. Over 1000 individual mosquitoes can be sequenced in a single MiSeq run, and we demonstrate the panel’s power to assign species identity using sequencing data for 40 species from Africa, Southeast Asia, and South America. We also show that the approach can be used to resolve geographic population structure within An. gambiae and An. coluzzii populations, as the population structure determined based on these 62 loci from over 1000 mosquitoes closely mirrors that revealed through whole genome sequencing. The end-to-end approach is quick, inexpensive, robust, and accurate, which makes it a promising technique for very large-scale mosquito genetic surveillance and vector control.     

A Bayesian Hierarchical Model for Estimation of Abundance and Spatial Density of Aedes aegypti

Strategies to minimize dengue transmission commonly rely on vector control, which aims to maintain Ae. aegypti density below a theoretical threshold. Mosquito abundance is traditionally estimated from mark-release-recapture (MRR) experiments, which lack proper analysis regarding accurate vector spatial distribution and population density. Recently proposed strategies to control vector-borne diseases involve replacing the susceptible wild population by genetically modified individuals’ refractory to the infection by the pathogen. Accurate measurements of mosquito abundance in time and space are required to optimize the success of such interventions. In this paper, we present a hierarchical probabilistic model for the estimation of population abundance and spatial distribution from typical mosquito MRR experiments, with direct application to the planning of these new control strategies. We perform a Bayesian analysis using the model and data from two MRR experiments performed in a neighborhood of Rio de Janeiro, Brazil, during both low- and high-dengue transmission seasons. The hierarchical model indicates that mosquito spatial distribution is clustered during the winter (0.99 mosquitoes/premise 95% CI: 0.80–1.23) and more homogeneous during the high abundance period (5.2 mosquitoes/premise 95% CI: 4.3–5.9). The hierarchical model also performed better than the commonly used Fisher-Ford’s method, when using simulated data. The proposed model provides a formal treatment of the sources of uncertainty associated with the estimation of mosquito abundance imposed by the sampling design. Our approach is useful in strategies such as population suppression or the displacement of wild vector populations by refractory Wolbachia-infected mosquitoes, since the invasion dynamics have been shown to follow threshold conditions dictated by mosquito abundance. The presence of spatially distributed abundance hotspots is also formally addressed under this modeling framework and its knowledge deemed crucial to predict the fate of transmission control strategies based on the replacement of vector populations.     

Chromosome evolution in malaria mosquitoes

The genomic era offers excellent opportunities to improve our understanding of the genetic basis of mosquito adaptation, evolution, and competence to a pathogen. The availability of polytene chromosomes in anopheline mosquitoes makes them an excellent model system for studying genome organization, evolution, and function. Physical mapping facilitated the whole genome sequence assembly for the major malaria vector Anopheles gambiae and comparative genome mapping has determined types, patterns, and rates of chromosomal rearrangements in mosquito evolution. Together with sequencing projects, high-resolution physical mapping can shed light on mechanisms of chromosomal rearrangements and phylogenetic relation-ships among species.     

Evolutionary potential of the extrinsic incubation period of dengue virus in Aedes aegypti

Dengue fever is the most common arboviral disease worldwide. It is caused by dengue viruses (DENV) and the mosquito Aedes aegypti is its primary vector. One of the most powerful determinants of a mosquito's ability to transmit DENV is the length of the extrinsic incubation period (EIP), the time it takes for a virus to be transmitted by a mosquito after consuming an infected blood meal. Here, we repeatedly measured DENV load in the saliva of individual mosquitoes over their lifetime and used this in combination with a breeding design to determine the extent to which EIP might respond to the evolutionary forces of drift and selection. We demonstrated that genetic variation among mosquitoes contributes significantly to transmission potential and length of EIP. We reveal that shorter EIP is genetically correlated with reduced mosquito lifespan, highlighting negative life‐history consequences for virus‐infected mosquitoes. This work highlights the capacity for local genetic variation in mosquito populations to evolve and to dramatically affect the nature of human outbreaks. It also provides the impetus for isolating mosquito genes that determine EIP. More broadly, our dual experimental approach offers new opportunities for studying the evolutionary potential of transmission traits in other vector/pathogen systems.     

Wolbachia infection does not alter attraction of the mosquito Aedes (Stegomyia) aegypti to human odours

The insect endosymbiont Wolbachia pipientis (Rickettsiales: Rickettsiaceae) is undergoing field trials around the world to determine if it can reduce transmission of dengue virus from the mosquito Stegomyia aegypti to humans. Two different Wolbachia strains have been released to date. The primary effect of the wMel strain is pathogen protection whereby infection with the symbiont limits replication of dengue virus inside the mosquito. A second strain, wMelPop, induces pathogen protection, reduces the adult mosquito lifespan and decreases blood feeding success in mosquitoes after 15 days of age. Here we test whether Wolbachia infection affects mosquito attraction to host odours in adults aged 5 and 15 days. We found no evidence of reduced odour attraction of mosquitoes, even for those infected with the more virulent wMelPop. This bodes well for fitness and competitiveness in the field given that the mosquitoes must find hosts to reproduce for the biocontrol method to succeed.     

COSTS AND BENEFITS OF MOSQUITO REFRACTORINESS TO MALARIA PARASITES: IMPLICATIONS FOR GENETIC VARIABILITY OF MOSQUITOES AND GENETIC CONTROL OF MALARIA

The problem of fitness costs associated with host resistance to parasitism is related to the evolution of parasite virulence, population genetic diversity and the dynamics of host-parasite relationships, and proposed strategies for disease control through the genetic manipulation of mosquito vectors. Two Aedes aegypti populations, refractory and susceptible to Plasmodium gallinaceum, were previously selected from the Moyo-In-Dry strain (MOYO) through inbreeding (F = 0.5). Reproductive success and survivorship of the two populations were compared, and the influence of the parasite on mosquito fitness also was evaluated. Fitness components studied include fecundity, adult survivorship and egg-to-adult developmental time, blood-meal size, and adult body size. The refractory population has a significantly shorter egg-to-adult developmental time and a smaller body size, takes a smaller blood meal, and subsequently lays fewer eggs than the susceptible population. The mean longevity of the refractory population is significantly shorter than the susceptible population. Exposure to the parasite exhibited little effect on the survivorship and fecundity of either population. Several factors may contribute to the lower fitness of the refractory population, including founder effect, inbreeding depression, the effect of other uncharacterized genes linked to genes conferring refractoriness, and pleiotropic effects associated with these genes. The results are discussed in relation to the genetic diversity of natural mosquito populations and their implications for the genetic control of malaria.     

Quantitative trait loci controlling refractoriness to Plasmodium falciparum in natural Anopheles gambiae mosquitoes from a malaria-endemic region in western Kenya

Natural anopheline populations exhibit much variation in ability to support malaria parasite development, but the genetic mechanisms underlying this variation are not clear. Previous studies in Mali, West Africa, identified two quantitative trait loci (QTL) in Anopheles gambiae mosquitoes that confer refractoriness (failure of oocyst development in mosquito midguts) to natural Plasmodium falciparum parasites. We hypothesize that new QTL may be involved in mosquito refractoriness to malaria parasites and that the frequency of natural refractoriness genotypes may be higher in the basin region of Lake Victoria, East Africa, where malaria transmission intensity and parasite genetic diversity are among the highest in the world. Using field-derived F2 isofemale families and microsatellite marker genotyping, two loci significantly affecting oocyst density were identified: one on chromosome 2 between markers AG2H135 and AG2H603 and the second on chromosome 3 near marker AG3H93. The first locus was detected in three of the five isofemale families studied and colocalized to the same region as Pen3 and pfin1 described in other studies. The second locus was detected in two of the five isofemale families, and it appears to be a new QTL. QTL on chromosome 2 showed significant additive effects while those on chromosome 3 exhibited significant dominant effects. Identification of P. falciparum-refractoriness QTL in natural An. gambiae mosquitoes is critical to the identification of the genes involved in malaria parasite transmission in nature and for understanding the coevolution between malaria parasites and mosquito vectors.     

Genetic Mapping of Specific Interactions between Aedes aegypti Mosquitoes and Dengue Viruses

Specific interactions between host genotypes and pathogen genotypes (G×G interactions) are commonly observed in invertebrate systems. Such specificity challenges our current understanding of invertebrate defenses against pathogens because it contrasts the limited discriminatory power of known invertebrate immune responses. Lack of a mechanistic explanation, however, has questioned the nature of host factors underlying G×G interactions. In this study, we aimed to determine whether G×G interactions observed between dengue viruses and their Aedes aegypti vectors in nature can be mapped to discrete loci in the mosquito genome and to document their genetic architecture. We developed an innovative genetic mapping strategy to survey G×G interactions using outbred mosquito families that were experimentally exposed to genetically distinct isolates of two dengue virus serotypes derived from human patients. Genetic loci associated with vector competence indices were detected in multiple regions of the mosquito genome. Importantly, correlation between genotype and phenotype was virus isolate-specific at several of these loci, indicating G×G interactions. The relatively high percentage of phenotypic variation explained by the markers associated with G×G interactions (ranging from 7.8% to 16.5%) is consistent with large-effect host genetic factors. Our data demonstrate that G×G interactions between dengue viruses and mosquito vectors can be assigned to physical regions of the mosquito genome, some of which have a large effect on the phenotype. This finding establishes the existence of tangible host genetic factors underlying specific interactions between invertebrates and their pathogens in a natural system. Fine mapping of the uncovered genetic loci will elucidate the molecular mechanisms of mosquito-virus specificity.     

Isolation of microsatellite markers in the tiger mosquito Aedes albopictus (Skuse)

Among the numerous molecular markers available in population genetics, microsatellites are one of the most powerful tools developed in recent years. This paper describes the isolation of six polymorphic microsatellite loci in the tiger mosquito Aedes albopictus using an enriched genomic library technique. Such loci should be an efficient tool in population genetic studies for this mosquito species.     

A dengue receptor as possible genetic marker of vector competence in <Emphasis Type="Italic">Aedes aegypti</Emphasis>

Background  

Vector competence refers to the intrinsic permissiveness of an arthropod vector for infection, replication and transmission of a virus. Notwithstanding studies of Quantitative Trait Loci (QTL) that influence the ability of Aedes aegypti midgut (MG) to become infected with dengue virus (DENV), no study to date has been undertaken to identify genetic markers of vector competence. Furthermore, it is known that mosquito populations differ greatly in their susceptibility to flaviviruses. Differences in vector competence may, at least in part, be due to the presence of specific midgut epithelial receptors and their identification would be a significant step forward in understanding the interaction of the virus with the mosquito. The first interaction of DENV with the insect is through proteins in the apical membrane of the midgut epithelium resulting in binding and receptor-mediated endocytosis of the virus, and this determines cell permissiveness to infection. The susceptibility of mosquitoes to infection may therefore depend on their specific virus receptors. To study this interaction in Ae. aegypti strains that differ in their vector competence for DENV, we investigated the DS3 strain (susceptible to DENV), the IBO-11 strain (refractory to infection) and the membrane escape barrier strain, DMEB, which is infected exclusively in the midgut epithelial cells.     

Invasive mosquitoes, larval competition, and indirect effects on the vector competence of native mosquito species (Diptera: Culicidae)

Invasive arthropods that vector pathogens have the potential to influence pathogen transmission both directly, by becoming a novel pathogen vector, or indirectly, by interacting with native vectors. Adult mosquito size is influenced by food availability in the larval stage, and smaller, nutrient-deprived mosquitoes are, in some studies, more efficient viral vectors in the laboratory. This is the first study to examine the indirect impacts that larval competition between Aedes albopictus, an introduced mosquito species, and Ochlerotatus triseriatus, a native mosquito species and the primary vector for La Crosse virus (LACV) in the US, has on native mosquito larval survival, adult size, and vector competence. A. albopictus presence decreased Oc. triseriatus larval survival, but surviving Oc. triseriatus females were larger, potentially owing to a release from intraspecific competition. These larger, native females were more likely to develop both midgut and disseminated LACV infections than females emerging from monospecific treatments. Collectively, these results suggest a need to better understand the ecology of both native and invasive vector species, their interactions, and the potential for those interactions to alter vector-borne disease transmission.     

Effects of the removal or reduction in density of the malaria mosquito,Anopheles gambiae s.l., on interacting predators and competitors in local ecosystems

New genetic control methods for mosquitoes may reduce vector species without direct effects on other species or the physical environment common with insecticides or drainage. Effects on predators and competitors could, however, be a concern as Anopheles gambiae s.l. is preyed upon in all life stages. We overview the literature and assess the strength of the ecological interactions identified. Most predators identified consume many other insect species and there is no evidence that any species preys exclusively on any anopheline mosquito. There is one predatory species with a specialisation on blood‐fed mosquitoes including An. gambiae s.l.. Evarcha culicivora is a jumping spider, known as the vampire spider, found around Lake Victoria. There is no evidence that these salticids require Anopheles mosquitoes and will readily consume blood‐fed Culex. Interspecific competition studies focus on other mosquitoes of larval habitats. Many of these take place in artificial cosms and give contrasting results to semi‐field studies. This may limit their extrapolation regarding the potential impact of reduced An. gambiae numbers. Previous mosquito control interventions are informative and identify competitive release and niche opportunism; so while the identity and relative abundance of the species present may change, the biomass available to predators may not.     

Midgut Microbiota of the Malaria Mosquito Vector Anopheles gambiae and Interactions with Plasmodium falciparum Infection

The susceptibility of Anopheles mosquitoes to Plasmodium infections relies on complex interactions between the insect vector and the malaria parasite. A number of studies have shown that the mosquito innate immune responses play an important role in controlling the malaria infection and that the strength of parasite clearance is under genetic control, but little is known about the influence of environmental factors on the transmission success. We present here evidence that the composition of the vector gut microbiota is one of the major components that determine the outcome of mosquito infections. A. gambiae mosquitoes collected in natural breeding sites from Cameroon were experimentally challenged with a wild P. falciparum isolate, and their gut bacterial content was submitted for pyrosequencing analysis. The meta-taxogenomic approach revealed a broader richness of the midgut bacterial flora than previously described. Unexpectedly, the majority of bacterial species were found in only a small proportion of mosquitoes, and only 20 genera were shared by 80% of individuals. We show that observed differences in gut bacterial flora of adult mosquitoes is a result of breeding in distinct sites, suggesting that the native aquatic source where larvae were grown determines the composition of the midgut microbiota. Importantly, the abundance of Enterobacteriaceae in the mosquito midgut correlates significantly with the Plasmodium infection status. This striking relationship highlights the role of natural gut environment in parasite transmission. Deciphering microbe-pathogen interactions offers new perspectives to control disease transmission.     

A genetic module regulates the melanization response of Anopheles to Plasmodium

Two modes of refractoriness to Plasmodium, ookinete lysis and melanization, are known in the malaria vector, Anopheles gambiae. Melanization, a potent insect immune response, is manifested in a genetically selected refractory strain and in susceptible mosquitoes that are depleted of specific C-type lectins (CTLs). Here we use a systematic in vivo RNA interference-mediated reverse genetic screen and other recent results to define a melanization-regulating genetic module or network. It encompasses at least 14 genes, including those that encode five Easter-like clip domain serine proteases and four Masquerade-like serine protease homologues of the mosquito CLIPB and CLIPA subfamilies respectively. We show that several but not all CLIPB genes promote Plasmodium melanization, exhibiting partial functional overlap and synergy. We also report that several CLIPA genes have contrasting roles: CLIPA8 is essential for parasite melanization, while three other CLIPAs are novel synergistic inhibitors of this response. Importantly, the roles of certain CLIPAs and CLIPBs are strain specific, indicating that this network may differ between strains. Finally, we provide evidence that in susceptible mosquitoes melanization induced by knockdown of either CTL4 or CLIPA2/CLIPA5 directly kills ookinetes, in contrast to refractory mosquitoes where it merely disposes of dead parasites.     

Mosquito distribution in a saltmarsh: determinants of eggs in a variable environment

Two saltmarsh mosquitoes dominate the transmission of Ross River virus (RRV, Togoviridae: Alphavirus), one of Australia's most prominent mosquito‐borne diseases. Ecologically, saltmarshes vary in their structure, including habitat types, hydrological regimes, and diversity of aquatic fauna, all of which drive mosquito oviposition behavior. Understanding the distribution of vector mosquitoes within saltmarshes can inform early warning systems, surveillance, and management of vector populations. The aim of this study was to identify the distribution of Ae. camptorhynchus, a known vector for RRV, across a saltmarsh and investigate the influence that other invertebrate assemblage might have on Ae. camptorhynchus egg dispersal. We demonstrate that vegetation is a strong indicator for Ae. camptorhynchus egg distribution, and this was not correlated with elevation or other invertebrates located at this saltmarsh. Also, habitats within this marsh are less frequently inundated, resulting in dryer conditions. We conclude that this information can be applied in vector surveillance and monitoring of temperate saltmarsh environments and also provides a baseline for future investigations into understanding mosquito vector habitat requirements.     

Dinucleotide microsatellite markers from Anopheles funestus

A total of 39 dinucleotide microsatellite‐containing clones were sequenced from a principal African malaria vector, Anopheles funestus. Primers designed to amplify 20 loci were used to genotype A. funestus mosquitoes from Burkina Faso and Kenya. Of nine polymorphic loci that amplified reliably and could be scored unambiguously, the overall within‐sample gene diversity was similar between locales, 0.77 and 0.78, with an allelic richness per locus of five to 11. Both the high level of polymorphism and absence of significant heterozygote deficiency at any locus favour these markers for studies of population structure that are vital to controlling this medically important species.     

Polymorphic microsatellite markers for the human oomycete pathogen Pythium insidiosum

We describe the development and characterization of microsatellite loci from the human oomycete pathogen Pythium insidiosum. Nine of 15 microsatellite loci were shown to be appropriate for population genetic study. All loci were polymorphic with observed heterozygosities ranging from 0.241 to 0.912 and from five to 18 alleles per locus among 65 individuals in Thailand. These markers are being used to ascertain multilocus genotypes for molecular epidemiological and population genetic analyses of this little known human pathogen.     

CHROMOSOMAL REPATTERNING AND LINKAGE GROUP CONSERVATION IN MOSQUITO KARYOTYPIC EVOLUTION

Chromosome number and morphology in mosquitoes is remarkably uniform: virtually all mosquitoes have a diploid chromosome number of six (2N = 6), and their chromosomes are invariably metacentric or submetacentric. Numerical changes obviously have not been important in mosquito chromosomal evolution, and because of the morphological similarity of their chromosomes, it appears that structural changes have played little or no role in mosquito karyotypic evolution. The goal of the present study was to identify the types and relative numbers of chromosomal changes in mosquito evolution and to extend the comparison where possible to the higher diptera. To do this, we compared the enzyme linkage maps of six species of Aedes to each other and to enzyme maps of seven other mosquito species and to Drosophila melanogaster. Our results indicate that Aedes chromosomes have been modified by inversions, most which were paracentric, and by translocations, most which were Robertsonian. Intrageneric comparison of Aedes enzyme maps also revealed groups of linked enzyme loci whose integrity has been maintained throughout Aedes evolution (conserved linkages/syntenies). Intergeneric comparisons of Aedes enzyme maps with those of species in the genera Culex, Anopheles, and Toxorhynchites disclosed conserved associations of enzyme loci between mosquito genera. These findings lead us to postulate that the ancestral mosquito karyotype consisted of six chromosomal elements which, other than being combined in different ways in various mosquito groups, have remained essentially intact during mosquito evolution. Furthermore, the identification of groups of linked enzyme loci common to mosquitoes and to D. melanogaster indicates that linkage group conservation may characterize the karyotypic evolution of all dipteran insects.     

Vector competence and innate immune responses to dengue virus infection in selected laboratory and field‐collected Stegomyia aegypti (= Aedes aegypti)

Control of dengue virus (DenV) transmission, primarily based on strategies to reduce populations of the principle vector Stegomya aegypti (= Aedes aegypti) (Diptera: Culicidae), is difficult to sustain over time. Other potential strategies aim to manipulate characteristics such as vector competence (VC), the innate capacity of the vector to transmit the virus. Previous studies have identified genetic factors, including differential expression of apoptosis‐related genes, associated with the refractory and susceptible phenotypes in selected strains of S. aegypti from Cali, Colombia. The present study was designed to evaluate the variability of VC in selected strains against different DenV serotypes and to determine whether field‐collected mosquitoes respond similarly to selected laboratory strains in terms of enhanced or reduced expression of apoptosis‐related genes. Vector competence differed between strains, but did not differ in response to different DenV serotypes. Differences in VC were observed among mosquitoes collected from different localities in Cali. The overexpression of the pro‐apoptosis genes, caspase 16 and Aedronc, was conserved in field‐collected refractory mosquitoes and the selected laboratory refractory strain. The results suggest that the apoptosis response is conserved among all refractory mosquitoes to inhibit the development of all DenV serotypes.