Microevolution of Aedes aegypti

Scientific research into the epidemiology of dengue frequently focuses on the microevolution and dispersion of the mosquito Aedes aegypti. One of the world’s largest urban agglomerations infested by Ae. aegypti is the Brazilian megalopolis of Sao Paulo, where >26,900 cases of dengue were reported until June 2015. Unfortunately, the dynamics of the genetic variability of Ae. aegypti in the Sao Paulo area have not been well studied. To reduce this knowledge gap, we assessed the morphogenetic variability of a population of Ae. aegypti from a densely urbanised neighbourhood of Sao Paulo. We tested if allelic patterns could vary over a short term and if wing shape could be a predictor of the genetic variation. Over a period of 14 months, we examined the variation of genetic (microsatellites loci) and morphological (wing geometry) markers in Ae. aegypti. Polymorphisms were detected, as revealed by the variability of 20 microsatellite loci (115 alleles combined; overall F_st = 0.0358) and 18 wing landmarks (quantitative estimator Q_st = 0.4732). These levels of polymorphism are higher than typically expected to an exotic species. Allelic frequencies of the loci changed over time and temporal variation in the wing shape was even more pronounced, permitting high reclassification levels of chronological samples. In spite of the fact that both markers underwent temporal variation, no correlation was detected between their dynamics. We concluded that microevolution was detected despite the short observational period, but the intensities of change of the markers were discrepant. Wing shape failed from predicting allelic temporal variation. Possibly, natural selection (Q_st>F_st) or variance of expressivity of wing phenotype are involved in this discrepancy. Other possibly influential factors on microevolution of Ae. aegypti are worth searching. Additionally, the implications of the rapid evolution and high polymorphism of this mosquito vector on the efficacy of control methods have yet to be investigated.     

Wolbachia-Associated Bacterial Protection in the Mosquito Aedes aegypti

Background

Wolbachia infections confer protection for their insect hosts against a range of pathogens including bacteria, viruses, nematodes and the malaria parasite. A single mechanism that might explain this broad-based pathogen protection is immune priming, in which the presence of the symbiont upregulates the basal immune response, preparing the insect to defend against subsequent pathogen infection. A study that compared natural Wolbachia infections in Drosophila melanogaster with the mosquito vector Aedes aegypti artificially transinfected with the same strains has suggested that innate immune priming may only occur in recent host-Wolbachia associations. This same study also revealed that while immune priming may play a role in viral protection it cannot explain the entirety of the effect.

Methodology/Findings

Here we assess whether the level of innate immune priming induced by different Wolbachia strains in A. aegypti is correlated with the degree of protection conferred against bacterial pathogens. We show that Wolbachia strains wMel and wMelPop, currently being tested for field release for dengue biocontrol, differ in their protective abilities. The wMelPop strain provides stronger, more broad-based protection than wMel, and this is likely explained by both the higher induction of immune gene expression and the strain-specific activation of particular genes. We also show that Wolbachia densities themselves decline during pathogen infection, likely as a result of the immune induction.

Conclusions/Significance

This work shows a correlation between innate immune priming and bacterial protection phenotypes. The ability of the Toll pathway, melanisation and antimicrobial peptides to enhance viral protection or to provide the basis of malaria protection should be further explored in the context of this two-strain comparison. This work raises the questions of whether Wolbachia may improve the ability of wild mosquitoes to survive pathogen infection or alter the natural composition of gut flora, and thus have broader consequences for host fitness.     

TALEN-Based Gene Disruption in the Dengue Vector Aedes aegypti

In addition to its role as the primary vector for dengue viruses, Aedes aegypti has a long history as a genetic model organism for other bloodfeeding mosquitoes, due to its ease of colonization, maintenance and reproductive productivity. Though its genome has been sequenced, functional characterization of many Ae. aegypti genes, pathways and behaviors has been slow. TALE nucleases (TALENs) have been used with great success in a number of organisms to generate site-specific DNA lesions. We evaluated the ability of a TALEN pair to target the Ae. aegypti kmo gene, whose protein product is essential in the production of eye pigmentation. Following injection into pre-blastoderm embryos, 20–40% of fertile survivors produced kmo alleles that failed to complement an existing kh^w mutation. Most of these individuals produced more than 20% white-eyed progeny, with some producing up to 75%. Mutant alleles were associated with lesions of 1–7 bp specifically at the selected target site. White-eyed individuals could also be recovered following a blind intercross of G₁ progeny, yielding several new white-eyed strains in the genetic background of the sequenced Liverpool strain. We conclude that TALENs are highly active in the Ae. aegypti germline, and have the potential to transform how reverse genetic experiments are performed in this important disease vector.     

Wolbachia Reduces the Transmission Potential of Dengue-Infected Aedes aegypti

BackgroundDengue viruses (DENV) are the causative agents of dengue, the world’s most prevalent arthropod-borne disease with around 40% of the world’s population at risk of infection annually. Wolbachia pipientis, an obligate intracellular bacterium, is being developed as a biocontrol strategy against dengue because it limits replication of the virus in the mosquito. The Wolbachia strain wMel, which has been introduced into the mosquito vector, Aedes aegypti, has been shown to invade and spread to near fixation in field releases. Standard measures of Wolbachia’s efficacy for blocking virus replication focus on the detection and quantification of virus in mosquito tissues. Examining the saliva provides a more accurate measure of transmission potential and can reveal the extrinsic incubation period (EIP), that is, the time it takes virus to arrive in the saliva following the consumption of DENV viremic blood. EIP is a key determinant of a mosquito’s ability to transmit DENVs, as the earlier the virus appears in the saliva the more opportunities the mosquito will have to infect humans on subsequent bites.Conclusions/SignificanceThe saliva-based traits reported here offer more disease-relevant measures of Wolbachia’s effects on the vector and the virus. The lengthening of EIP highlights another means, in addition to the reduction of infection frequencies and DENV titers in mosquitoes, by which Wolbachia should operate to reduce DENV transmission in the field.     

Modeling the Non-Stationary Climate Dependent Temporal Dynamics of Aedes aegypti

Background

Temperature and humidity strongly affect the physiology, longevity, fecundity and dispersal behavior of Aedes aegypti, vector of dengue fever. Contrastingly, the statistical associations measured between time series of mosquito abundance and meteorological variables are often weak and contradictory. Here, we investigated the significance of these relationships at different time scales.

Methods and Findings

A time series of the adult mosquito abundance from a medium-sized city in Brazil, lasting 109 weeks was analyzed. Meteorological variables included temperature, precipitation, wind velocity and humidity. As analytical tools, generalized linear models (GLM) with time lags and interaction terms were used to identify average effects while the wavelet analysis was complementarily used to identify transient associations. The fitted GLM showed that mosquito abundance is significantly affected by the interaction between lagged temperature and humidity, and also by the mosquito abundance a week earlier. Extreme meteorological variables were the best predictors, and the mosquito population tended to increase at values above and 54% humidity. The wavelet analysis identified non-stationary local effects of these meteorological variables on abundance throughout the study period, with peaks in the spring-summer period. The wavelet detected weak but significant effects for precipitation and wind velocity.

Conclusion

Our results support the presence of transient relationships between meteorological variables and mosquito abundance. Such transient association may be explained by the ability of Ae. aegypti to buffer part of its response to climate, for example, by choosing sites with proper microclimate. We also observed enough coupling between the abundance and meteorological variables to develop a model with good predictive power. Extreme values of meteorological variables with time lags, interaction terms and previous mosquito abundance are strong predictors and should be considered when understanding the climate effect on mosquito abundance and population growth.     

The Role of piRNA-Mediated Epigenetic Silencing in the Population Dynamics of Transposable Elements in Drosophila melanogaster

The piwi-interacting RNAs (piRNA) are small RNAs that target selfish transposable elements (TEs) in many animal genomes. Until now, piRNAs’ role in TE population dynamics has only been discussed in the context of their suppression of TE transposition, which alone is not sufficient to account for the skewed frequency spectrum and stable containment of TEs. On the other hand, euchromatic TEs can be epigenetically silenced via piRNA-dependent heterochromatin formation and, similar to the widely known “Position-effect variegation”, heterochromatin induced by TEs can “spread” into nearby genes. We hypothesized that the piRNA-mediated spread of heterochromatin from TEs into adjacent genes has deleterious functional effects and leads to selection against individual TEs. Unlike previously identified deleterious effects of TEs due to the physical disruption of DNA, the functional effect we investigated here is mediated through the epigenetic influences of TEs. We found that the repressive chromatin mark, H3K9me, is elevated in sequences adjacent to euchromatic TEs at multiple developmental stages in Drosophila melanogaster. Furthermore, the heterochromatic states of genes depend not only on the number of and distance from adjacent TEs, but also on the likelihood that their nearest TEs are targeted by piRNAs. These variations in chromatin status probably have functional consequences, causing genes near TEs to have lower expression. Importantly, we found stronger selection against TEs that lead to higher H3K9me enrichment of adjacent genes, demonstrating the pervasive evolutionary consequences of TE-induced epigenetic silencing. Because of the intrinsic biological mechanism of piRNA amplification, spread of TE heterochromatin could result in the theoretically required synergistic deleterious effects of TE insertions for stable containment of TE copy number. The indirect deleterious impact of piRNA-mediated epigenetic silencing of TEs is a previously unexplored, yet important, element for the evolutionary dynamics of TEs.     

Functional and Genetic Characterization of Neuropeptide Y-Like Receptors in Aedes aegypti

Background

Female Aedes aegypti mosquitoes are the principal vector for dengue fever, causing 50–100 million infections per year, transmitted between human and mosquito by blood feeding. Ae. aegypti host-seeking behavior is known to be inhibited for three days following a blood meal by a hemolymph-borne humoral factor. Head Peptide-I is a candidate peptide mediating this suppression, but the mechanism by which this peptide alters mosquito behavior and the receptor through which it signals are unknown.

Methodology/Principal Findings

Head Peptide-I shows sequence similarity to short Neuropeptide-F peptides (sNPFs) that have been implicated in feeding behaviors and are known to signal through Neuropeptide Y (NPY)-Like Receptors (NPYLRs). We identified eight NPYLRs in the Ae. aegypti genome and screened each in a cell-based calcium imaging assay for sensitivity against a panel of peptides. Four of the Ae. aegypti NPYLRs responded to one or more peptide ligands, but only NYPLR1 responded to Head Peptide-I as well as sNPFs. Two NPYLR1 homologues identified in the genome of the Lyme disease vector, Ixodes scapularis, were also sensitive to Head Peptide-I. Injection of synthetic Head Peptide-I and sNPF-3 inhibited host-seeking behavior in non-blood-fed female mosquitoes, whereas control injections of buffer or inactive Head Peptide-I [Cys10] had no effect. To ask if NPYLR1 is necessary for blood-feeding-induced host-seeking inhibition, we used zinc-finger nucleases to generate five independent npylr1 null mutant strains and tested them for behavioral abnormalities. npylr1 mutants displayed normal behavior in locomotion, egg laying, sugar feeding, blood feeding, host seeking, and inhibition of host seeking after a blood meal.

Conclusions

In this work we deorphanized four Ae. aegypti NPYLRs and identified NPYLR1 as a candidate sNPF receptor that is also sensitive to Head Peptide-I. Yet npylr1 alone is not required for host-seeking inhibition and we conclude that other receptors, additional peptides, or both, regulate this important behavior.     

Repellents Inhibit P450 Enzymes in Stegomyia (Aedes) aegypti

The primary defence against mosquitoes and other disease vectors is often the application of a repellent. Despite their common use, the mechanism(s) underlying the activity of repellents is not fully understood, with even the mode of action of DEET having been reported to be via different mechanisms; e.g. interference with olfactory receptor neurones or actively detected by olfactory receptor neurones on the antennae or maxillary palps. In this study, we discuss a novel mechanism for repellence, one of P450 inhibition. Thirteen essential oil extracts from Colombian plants were assayed for potency as P450 inhibitors, using a kinetic fluorometric assay, and for repellency using a modified World Health Organisation Pesticide Evaluations Scheme (WHOPES) arm-in cage assay with Stegomyia (Aedes) aegypti mosquitoes. Bootstrap analysis on the inhibition analysis revealed a significant correlation between P450-inhibition and repellent activity of the oils.     

Influence of breeding site availability on the oviposition behaviour of Aedes aegypti

Despite the importance of the mosquito Aedes aegypti in thetransmission of arboviruses, such as yellow fever, Chikungunya fever and denguefever, some aspects of their behaviour remain unknown. In the present study, theoviposition behaviour of Ae. aegypti females that were exposed todifferent densities of breeding sites (2, 4, 8 and 16) was evaluated in laboratoryand semi-field conditions. The number of breeding sites that were used wasproportional to the number available, but tended towards stabilisation. Females usedfour-six breeding sites on average, with a maximum of 11. A high percentage of eggswas observed in the water, along with the presence of a breeding site termed“favourite”, which received at least 40% of the eggs. The results are discussed inecological, evolutionary and epidemiological approaches.     

Evolutionary Genomics of Transposable Elements in Saccharomyces cerevisiae

Saccharomyces cerevisiae is one of the premier model systems for studying the genomics and evolution of transposable elements. The availability of the S. cerevisiae genome led to unprecedented insights into its five known transposable element families (the LTR retrotransposons Ty1-Ty5) in the years shortly after its completion. However, subsequent advances in bioinformatics tools for analysing transposable elements and the recent availability of genome sequences for multiple strains and species of yeast motivates new investigations into Ty evolution in S. cerevisiae. Here we provide a comprehensive phylogenetic and population genetic analysis of all Ty families in S. cerevisiae based on a systematic re-annotation of Ty elements in the S288c reference genome. We show that previous annotation efforts have underestimated the total copy number of Ty elements for all known families. In addition, we identify a new family of Ty3-like elements related to the S. paradoxus Ty3p which is composed entirely of degenerate solo LTRs. Phylogenetic analyses of LTR sequences identified three families with short-branch, recently active clades nested among long branch, inactive insertions (Ty1, Ty3, Ty4), one family with essentially all recently active elements (Ty2) and two families with only inactive elements (Ty3p and Ty5). Population genomic data from 38 additional strains of S. cerevisiae show that the majority of Ty insertions in the S288c reference genome are fixed in the species, with insertions in active clades being predominantly polymorphic and insertions in inactive clades being predominantly fixed. Finally, we use comparative genomic data to provide evidence that the Ty2 and Ty3p families have arisen in the S. cerevisiae genome by horizontal transfer. Our results demonstrate that the genome of a single individual contains important information about the state of TE population dynamics within a species and suggest that horizontal transfer may play an important role in shaping the genomic diversity of transposable elements in unicellular eukaryotes.     

Splinkerette PCR for Mapping Transposable Elements in Drosophila

Transposable elements (such as the P-element and piggyBac) have been used to introduce thousands of transgenic constructs into the Drosophila genome. These transgenic constructs serve many roles, from assaying gene/cell function, to controlling chromosome arm rearrangement. Knowing the precise genomic insertion site for the transposable element is often desired. This enables identification of genomic enhancer regions trapped by an enhancer trap, identification of the gene mutated by a transposon insertion, or simplifying recombination experiments. The most commonly used transgene mapping method is inverse PCR (iPCR). Although usually effective, limitations with iPCR hinder its ability to isolate flanking genomic DNA in complex genomic loci, such as those that contain natural transposons. Here we report the adaptation of the splinkerette PCR (spPCR) method for the isolation of flanking genomic DNA of any P-element or piggyBac. We report a simple and detailed protocol for spPCR. We use spPCR to 1) map a GAL4 enhancer trap located inside a natural transposon, pinpointing a master regulatory region for olfactory neuron expression in the brain; and 2) map all commonly used centromeric FRT insertion sites. The ease, efficiency, and efficacy of spPCR could make it a favored choice for the mapping of transposable element in Drosophila.     

Application of wMelPop Wolbachia Strain to Crash Local Populations of Aedes aegypti

The endosymbiotic bacteria Wolbachia pipientis (wMel strain) has been successfully established in several populations of Aedes aegypti, the primary dengue vector. The virulent Wolbachia strain wMelPop is known to cause several pathological impacts (increased egg mortality, life shortening, etc.) reducing overall fitness in the mosquito Ae. aegypti. Increased egg mortality could substantially reduce egg banks in areas with a lengthy monsoonal dry season, and be employed to eliminate local populations. We tested this application under semi-field cage conditions. First, we determined that wMelPop infection significantly reduced the survival of desiccation-resistant eggs of the dengue vector Ae. aegypti, with shade and temperature having a significant impact; nearly all wMelPop-infected eggs failed to hatch after 6 and 10 weeks in summer and winter conditions, respectively. In laboratory selection experiments we found that egg desiccation resistance can be increased by selection, and that this effect of wMelPop infection is due to the nuclear background of the host rather than Wolbachia. We then conducted an invasion of wMelPop within a semi-field cage using sustained weekly releases of wMelPop infected mosquitoes, with fixation achieved after 9 weeks. The egg populations wMelPop infected and an uninfected control were then subjected to a simulated prolonged monsoonal dry season (2.5 months) before flooding to induce hatching. The wMelPop infected eggs suffered significantly greater mortality than the controls, with only 0.67% and 4.35% of respective infected and uninfected eggs held in 99% shade hatching after 80 days. These studies suggest that wMelPop could be used to locally eliminate populations of Ae. aegypti that are exposed to prolonged dry conditions, particularly if combined with vector control.     

Rapid Intraspecific Evolution of miRNA and siRNA Genes in the Mosquito Aedes aegypti

RNA silencing, or RNA interference (RNAi) in metazoans mediates development, reduces viral infection and limits transposon mobility. RNA silencing involves 21–30 nucleotide RNAs classified into microRNA (miRNA), exogenous and endogenous small interfering RNAs (siRNA), and Piwi-interacting RNA (piRNA). Knock-out, silencing and mutagenesis of genes in the exogenous siRNA (exo-siRNA) regulatory network demonstrate the importance of this RNAi pathway in antiviral immunity in Drosophila and mosquitoes. In Drosophila, genes encoding components for processing exo-siRNAs are among the fastest evolving 3% of all genes, suggesting that infection with pathogenic RNA viruses may drive diversifying selection in their host. In contrast, paralogous miRNA pathway genes do not evolve more rapidly than the genome average. Silencing of exo-siRNA pathway genes in mosquitoes orally infected with arboviruses leads to increased viral replication, but little is known about the comparative patterns of molecular evolution among the exo-siRNA and miRNA pathways genes in mosquitoes. We generated nearly complete sequences of all exons of major miRNA and siRNA pathway genes dicer-1 and dicer-2, argonaute-1 and argonaute-2, and r3d1 and r2d2 in 104 Aedes aegypti mosquitoes collected from six distinct geographic populations and analyzed their genetic diversity. The ratio of replacement to silent amino acid substitutions was 1.4 fold higher in dicer-2 than in dicer-1, 27.4 fold higher in argonaute-2 than in argonaute-1 and similar in r2d2 and r3d1. Positive selection was supported in 32% of non-synonymous sites in dicer-1, in 47% of sites in dicer-2, in 30% of sites in argonaute-1, in all sites in argonaute-2, in 22% of sites in r3d1 and in 55% of sites in r2d2. Unlike Drosophila, in Ae. aegypti, both exo-siRNA and miRNA pathway genes appear to be undergoing rapid, positive, diversifying selection. Furthermore, refractoriness of mosquitoes to infection with dengue virus was significantly positively correlated for nucleotide diversity indices in dicer-2.     

Characterization of an Enantioselective Odorant Receptor in the Yellow Fever Mosquito Aedes aegypti

Enantiomers differ only in the left or right handedness (chirality) of their orientations and exhibit identical chemical and physical properties. In chemical communication systems, enantiomers can be differentially active at the physiological and behavioral levels. Only recently were enantioselective odorant receptors demonstrated in mammals while their existence in insects has remained hypothetical. Using the two-microelectrode voltage clamp of Xenopus oocytes, we show that the yellow fever mosquito, Aedes aegypti, odorant receptor 8 (AaOR8) acts as a chiral selective receptor for the (R)-(—)-enantiomer of 1-octen-3-ol, which in the presence of other kairomones is an attractant used by blood-sucking insects to locate their hosts. In addition to steric constraints, chain length and degree of unsaturation play important roles in this recognition process. This is the first characterization of an enantioselective odorant receptor in insects and the results demonstrate that an OR alone, without helper proteins, can account for chiral specificity exhibited by olfactory sensory neurons (OSNs).     

QTL Mapping of Genome Regions Controlling Temephos Resistance in Larvae of the Mosquito Aedes aegypti

Introduction

The mosquito Aedes aegypti is the principal vector of dengue and yellow fever flaviviruses. Temephos is an organophosphate insecticide used globally to suppress Ae. aegypti larval populations but resistance has evolved in many locations.

Methodology/Principal Findings

Quantitative Trait Loci (QTL) controlling temephos survival in Ae. aegypti larvae were mapped in a pair of F₃ advanced intercross lines arising from temephos resistant parents from Solidaridad, México and temephos susceptible parents from Iquitos, Peru. Two sets of 200 F₃ larvae were exposed to a discriminating dose of temephos and then dead larvae were collected and preserved for DNA isolation every two hours up to 16 hours. Larvae surviving longer than 16 hours were considered resistant. For QTL mapping, single nucleotide polymorphisms (SNPs) were identified at 23 single copy genes and 26 microsatellite loci of known physical positions in the Ae. aegypti genome. In both reciprocal crosses, Multiple Interval Mapping identified eleven QTL associated with time until death. In the Solidaridad×Iquitos (SLD×Iq) cross twelve were associated with survival but in the reciprocal IqxSLD cross, only six QTL were survival associated. Polymorphisms at acetylcholine esterase (AchE) loci 1 and 2 were not associated with either resistance phenotype suggesting that target site insensitivity is not an organophosphate resistance mechanism in this region of México.

Conclusions/Significance

Temephos resistance is under the control of many metabolic genes of small effect and dispersed throughout the Ae. aegypti genome.     

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.