The global population at risk from mosquito-borne diseases-including dengue, yellow fever, chikungunya and Zika-is expanding in concert with changes in the distribution of two key vectors: Aedes aegypti and Aedes albopictus. The distribution of these species is largely driven by both human movement and the presence of suitable climate. Using statistical mapping techniques, we show that human movement patterns explain the spread of both species in Europe and the United States following their introduction. We find that the spread of Ae. aegypti is characterized by long distance importations, while Ae. albopictus has expanded more along the fringes of its distribution. We describe these processes and predict the future distributions of both species in response to accelerating urbanization, connectivity and climate change. Global surveillance and control efforts that aim to mitigate the spread of chikungunya, dengue, yellow fever and Zika viruses must consider the so far unabated spread of these mosquitos. Our maps and predictions offer an opportunity to strategically target surveillance and control programmes and thereby augment efforts to reduce arbovirus burden in human populations globally.

Female Aedes aegypti mosquitoes infect more than 400 million people each year with dangerous viral pathogens including dengue, yellow fever, Zika and chikungunya. Progress in understanding the biology of mosquitoes and developing the tools to fight them has been slowed by the lack of a high-quality genome assembly. Here we combine diverse technologies to produce the markedly improved, fully re-annotated AaegL5 genome assembly, and demonstrate how it accelerates mosquito science. We anchored physical and cytogenetic maps, doubled the number of known chemosensory ionotropic receptors that guide mosquitoes to human hosts and egg-laying sites, provided further insight into the size and composition of the sex-determining M locus, and revealed copy-number variation among glutathione S-transferase genes that are important for insecticide resistance. Using high-resolution quantitative trait locus and population genomic analyses, we mapped new candidates for dengue vector competence and insecticide resistance. AaegL5 will catalyse new biological insights and intervention strategies to fight this deadly disease vector.

Dengue virus (DENV) is an RNA virus transmitted among humans by mosquito vectors, mainly Aedes aegypti. DENV transmission requires viral dissemination from the mosquito midgut to the salivary glands. During this process the virus undergoes several population bottlenecks, which are stochastic reductions in population size that restrict intra-host viral genetic diversity and limit the efficiency of natural selection. Despite the implications for virus transmission and evolution, DENV replication in salivary glands has not been directly demonstrated. Here, we used a strand-specific quantitative RT-PCR assay to demonstrate that negative-strand DENV RNA is produced in Ae. aegypti salivary glands, providing conclusive evidence that viral replication occurs in this tissue. Furthermore, we showed that the concentration of DENV genomic RNA in salivary glands increases significantly over time, indicating that active replication likely replenishes DENV genetic diversity prior to transmission. These findings improve our understanding of the biological determinants of DENV fitness and evolution.

Arthropod-borne viruses (arboviruses) transmitted by Aedes aegypti mosquitoes are an increasing threat to global health. The small interfering RNA (siRNA) pathway is considered the main antiviral immune pathway of insects, but its effective impact on arbovirus transmission is surprisingly poorly understood. Here, we use CRISPR-Cas9-mediated gene editing in vivo to mutate Dicer2, a gene encoding the RNA sensor and key component of the siRNA pathway. The loss of Dicer2 enhances early viral replication and systemic viral dissemination of four medically significant arboviruses (chikungunya, Mayaro, dengue, and Zika viruses) representing two viral families. However, Dicer2 mutants and wild-type mosquitoes display overall similar levels of vector competence. In addition, Dicer2 mutants undergo significant virus-induced mortality during infection with chikungunya virus. Together, our results define a multifaceted role for Dicer2 in the transmission of arboviruses by Ae. aegypti mosquitoes and pave the way for further mechanistic investigations.

Endogenous viral elements derived from nonretroviral RNA viruses have been described in various animal genomes. Whether they have a biological function, such as host immune protection against related viruses, is a field of intense study. Here, we investigated the repertoire of endogenous flaviviral elements (EFVEs) in Aedes mosquitoes, the vectors of arboviruses such as dengue and chikungunya viruses. Previous studies identified three EFVEs from Aedes albopictus cell lines and one from Aedes aegypti cell lines. However, an in-depth characterization of EFVEs in wild-type mosquito populations and individual mosquitoes in vivo has not been performed. We detected the full-length DNA sequence of the previously described EFVEs and their respective transcripts in several A. albopictus and A. aegypti populations from geographically distinct areas. However, EFVE-derived proteins were not detected by mass spectrometry. Using deep sequencing, we detected the production of PIWI-interacting RNA-like small RNAs, in an antisense orientation, targeting the EFVEs and their flanking regions in vivo The EFVEs were integrated in repetitive regions of the mosquito genomes, and their flanking sequences varied among mosquito populations. We bioinformatically predicted several new EFVEs from a Vietnamese A. albopictus population and observed variation in the occurrence of those elements among mosquitoes. Phylogenetic analysis of an A. aegypti EFVE suggested that it integrated prior to the global expansion of the species and subsequently diverged among and within populations. The findings of this study together reveal the substantial structural and nucleotide diversity of flaviviral integrations in Aedes genomes. Unraveling this diversity will help to elucidate the potential biological function of these EFVEs.IMPORTANCE Endogenous viral elements (EVEs) are whole or partial viral sequences integrated in host genomes. Interestingly, some EVEs have important functions for host fitness and antiviral defense. Because mosquitoes also have EVEs in their genomes, characterizing these EVEs is a prerequisite for their potential use to manipulate the mosquito antiviral response. In the study described here, we focused on EVEs related to the Flavivirus genus, to which dengue and Zika viruses belong, in individual Aedes mosquitoes from geographically distinct areas. We show the existence in vivo of flaviviral EVEs previously identified in mosquito cell lines, and we detected new ones. We show that EVEs have evolved differently in each mosquito population. They produce transcripts and small RNAs but not proteins, suggesting a function at the RNA level. Our study uncovers the diverse repertoire of flaviviral EVEs in Aedes mosquito populations and contributes to an understanding of their role in the host antiviral system.

Endogenous viral elements (EVEs) are viral sequences integrated in host genomes. A large number of non-retroviral EVEs was recently detected in Aedes mosquito genomes, leading to the hypothesis that mosquito EVEs may control exogenous infections by closely related viruses. Here, we experimentally investigated the role of an EVE naturally found in Aedes aegypti populations and derived from the widespread insect-specific virus, cell-fusing agent virus (CFAV). Using CRISPR-Cas9 genome editing, we created an Ae. aegypti line lacking the CFAV EVE. Absence of the EVE resulted in increased CFAV replication in ovaries, possibly modulating vertical transmission of the virus. Viral replication was controlled by targeting of viral RNA by EVE-derived P-element-induced wimpy testis-interacting RNAs (piRNAs). Our results provide evidence that antiviral piRNAs are produced in the presence of a naturally occurring EVE and its cognate virus, demonstrating a functional link between non-retroviral EVEs and antiviral immunity in a natural insect-virus interaction.

Environmental factors such as temperature can alter mosquito vector competence for arboviruses. Results from recent studies indicate that daily fluctuations around an intermediate mean temperature (26°C) reduce vector competence of Aedes aeygpti for dengue viruses (DENV). Theoretical predictions suggest that the mean temperature in combination with the magnitude of the diurnal temperature range (DTR) mediate the direction of these effects.

Mechanisms and evolutionary dynamics of sex-determination systems are of particular interest in insect vectors of human pathogens like mosquitoes because novel control strategies aim to convert pathogen-transmitting females into nonbiting males, or rely on accurate sexing for the release of sterile males. In Aedes aegypti, the main vector of dengue and Zika viruses, sex determination is governed by a dominant male-determining locus, previously thought to reside within a small, nonrecombining, sex-determining region (SDR) of an otherwise homomorphic sex chromosome. Here, we provide evidence that sex chromosomes in Ae. aegypti are genetically differentiated between males and females over a region much larger than the SDR. Our linkage mapping intercrosses failed to detect recombination between X and Y chromosomes over a 123-Mbp region (40% of their physical length) containing the SDR. This region of reduced male recombination overlapped with a smaller 63-Mbp region (20% of the physical length of the sex chromosomes) displaying high male-female genetic differentiation in unrelated wild populations from Brazil and Australia and in a reference laboratory strain originating from Africa. In addition, the sex-differentiated genomic region was associated with a significant excess of male-to-female heterozygosity and contained a small cluster of loci consistent with Y-specific null alleles. We demonstrate that genetic differentiation between sex chromosomes is sufficient to assign individuals to their correct sex with high accuracy. We also show how data on allele frequency differences between sexes can be used to estimate linkage disequilibrium between loci and the sex-determining locus. Our discovery of large-scale genetic differentiation between sex chromosomes in Ae. aegypti lays a new foundation for mapping and population genomic studies, as well as for mosquito control strategies targeting the sex-determination pathway.

In most of the world, Dengue virus (DENV) is mainly transmitted by the mosquito Aedes aegypti while in Europe, Aedes albopictus is responsible for human DENV cases since 2010. Identifying mutations that make DENV more competent for transmission by Ae. albopictus will help to predict emergence of epidemic strains. Ten serial passages in vivo in Ae. albopictus led to select DENV-1 strains with greater infectivity for this vector in vivo and in cultured mosquito cells. These changes were mediated by multiple adaptive mutations in the virus genome, including a mutation at position 10,418 in the DENV 3'UTR within an RNA stem-loop structure involved in subgenomic flavivirus RNA production. Using reverse genetics, we showed that the 10,418 mutation alone does not confer a detectable increase in transmission efficiency in vivo. These results reveal the complex adaptive landscape of DENV transmission by mosquitoes and emphasize the role of epistasis in shaping evolutionary trajectories of DENV variants.

Host-associated microbes, collectively known as the microbiota, play an important role in the biology of multicellular organisms. In mosquito vectors of human pathogens, the gut bacterial microbiota influences vectorial capacity and has become the subject of intense study. In laboratory studies of vector biology, genetic effects are often inferred from differences between geographically and genetically diverse colonies of mosquitoes that are reared in the same insectary. It is unclear, however, to what extent genetic effects can be confounded by uncontrolled differences in the microbiota composition among mosquito colonies. To address this question, we used 16S metagenomics to compare the midgut bacterial microbiome of six laboratory colonies of Aedes aegypti recently derived from wild populations representing the geographical range and genetic diversity of the species.

Aedes aegypti mosquitoes are the main vectors of arthropod-borne viruses (arboviruses) of public health significance, such as the flaviviruses dengue virus (DENV) and Zika virus (ZIKV). Mosquitoes are also the natural hosts of a wide range of viruses that are insect specific, raising the question of their influence on arbovirus transmission in nature. Cell-fusing agent virus (CFAV) was the first described insect-specific flavivirus, initially discovered in an A. aegypti cell line and subsequently detected in natural A. aegypti populations. It was recently shown that DENV and the CFAV strain isolated from the A. aegypti cell line have mutually beneficial interactions in mosquito cells in culture. However, whether natural strains of CFAV and DENV interact in live mosquitoes is unknown. Using a wild-type CFAV isolate recently derived from Thai A. aegypti mosquitoes, we found that CFAV negatively interferes with both DENV type 1 and ZIKV in vitro and in vivo For both arboviruses, prior infection by CFAV reduced the dissemination titer in mosquito head tissues. Our results indicate that the interactions observed between arboviruses and the CFAV strain derived from the cell line might not be a relevant model of the viral interference that we observed in vivo Overall, our study supports the hypothesis that insect-specific flaviviruses may contribute to reduce the transmission of human-pathogenic flaviviruses.IMPORTANCE The mosquito Aedes aegypti carries several arthropod-borne viruses (arboviruses) that are pathogenic to humans, including dengue and Zika viruses. Interestingly, A. aegypti is also naturally infected with insect-only viruses, such as cell-fusing agent virus. Although interactions between cell-fusing agent virus and dengue virus have been documented in mosquito cells in culture, whether wild strains of cell-fusing agent virus interfere with arbovirus transmission by live mosquitoes was unknown. We used an experimental approach to demonstrate that cell-fusing agent virus infection reduces the propagation of dengue and Zika viruses in A. aegypti mosquitoes. These results support the idea that insect-only viruses in nature can modulate the ability of mosquitoes to carry arboviruses of medical significance and that they could possibly be manipulated to reduce arbovirus transmission.

African populations of the mosquito Aedes aegypti are usually considered less susceptible to infection by human-pathogenic flaviviruses than globally invasive populations found outside Africa. Although this contrast has been well documented for Zika virus (ZIKV), it is unclear to what extent it is true for dengue virus (DENV), the most prevalent flavivirus of humans. Addressing this question is complicated by substantial genetic diversity among DENV strains, most notably in the form of four genetic types (DENV1 to DENV4), that can lead to genetically specific interactions with mosquito populations. Here, we carried out a continent-wide survey of DENV susceptibility using a panel of field-derived Ae. aegypti colonies from across the African range of the species and a colony from Guadeloupe, French West Indies as non-African reference. We found considerable variation in the ability of African Ae. aegypti populations to acquire and replicate a panel of six DENV strains spanning the four DENV types. Although African Ae. aegypti populations were generally less susceptible than the reference non-African population from Guadeloupe, in several instances some African populations were equally or more susceptible than the Guadeloupe population. Moreover, the relative level of susceptibility between African mosquito populations depended on the DENV strain, indicating genetically specific interactions. We conclude that unlike ZIKV susceptibility, there is no clear-cut dichotomy in DENV susceptibility between African and non-African Ae. aegypti. DENV susceptibility of African Ae. aegypti populations is highly heterogeneous and largely governed by the specific pairing of mosquito population and DENV strain.

Successful transmission of a vector-borne pathogen relies on a complex life cycle in the arthropod vector that requires initial infection of the digestive tract followed by systemic viral dissemination. The time interval between acquisition and subsequent transmission of the pathogen, called the extrinsic incubation period, is one of the most influential parameters of vector-borne pathogen transmission. However, the dynamic nature of this process is often ignored because vector competence assays are sacrificial and rely on end-point measurements. Here, we report that individual Aedes aegypti mosquitoes release large amounts of dengue virus (DENV) RNA in their excreta that can be non-sacrificially detected over time following oral virus exposure. Further, we demonstrate that detection of DENV RNA in excreta from individual mosquitoes is correlated to systemic viral dissemination with high specificity (0.9-1) albeit moderate sensitivity (0.64-0.89). Finally, we illustrate the potential of our finding to detect biological differences in the dynamics of DENV dissemination in a proof-of-concept experiment. Individual measurements of the time required for systemic viral dissemination, a prerequisite for transmission, will be valuable to monitor the dynamics of DENV vector competence, to carry out quantitative genetics studies, and to evaluate the risk of DENV transmission in field settings.

The recent Zika virus (ZIKV) epidemic in the Americas highlights its potential public health threat. While the Asian ZIKV lineage has been identified as the main cause of the epidemic, the African lineage, which has been primarily confined to Africa, has shown evidence of higher transmissibility in Aedes mosquitoes. To gain a deeper understanding of this differential transmissibility, our study employed a combination of tissue-level infection kinetics and single-cell-level infection kinetics using in situ immunofluorescent staining. We discovered that the African ZIKV lineage propagates more rapidly and spreads more efficiently within mosquito cells and tissues than its Asian counterpart. This information lays the groundwork for future exploration of the viral and host determinants driving these variations in propagation efficiency.

The exogenous small interfering RNA (exo-siRNA) pathway is a key antiviral mechanism in the Aedes aegypti mosquito, a widely distributed vector of human-pathogenic arboviruses. This pathway is induced by virus-derived double-stranded RNAs (dsRNA) that are cleaved by the ribonuclease Dicer 2 (Dcr2) into predominantly 21 nucleotide (nt) virus-derived small interfering RNAs (vsiRNAs). These vsiRNAs are used by the effector protein Argonaute 2 within the RNA-induced silencing complex to cleave target viral RNA. Dcr2 contains several domains crucial for its activities, including helicase and RNase III domains. In Drosophila melanogaster Dcr2, the helicase domain has been associated with binding to dsRNA with blunt-ended termini and a processive siRNA production mechanism, while the platform-PAZ domains bind dsRNA with 3' overhangs and subsequent distributive siRNA production. Here we analyzed the contributions of the helicase and RNase III domains in Ae. aegypti Dcr2 to antiviral activity and to the exo-siRNA pathway. Conserved amino acids in the helicase and RNase III domains were identified to investigate Dcr2 antiviral activity in an Ae. aegypti-derived Dcr2 knockout cell line by reporter assays and infection with mosquito-borne Semliki Forest virus (Togaviridae, Alphavirus). Functionally relevant amino acids were found to be conserved in haplotype Dcr2 sequences from field-derived Ae. aegypti across different continents. The helicase and RNase III domains were critical for silencing activity and 21 nt vsiRNA production, with RNase III domain activity alone determined to be insufficient for antiviral activity. Analysis of 21 nt vsiRNA sequences (produced by functional Dcr2) to assess the distribution and phasing along the viral genome revealed diverse yet highly consistent vsiRNA pools, with predominantly short or long sequence overlaps including 19 nt overlaps (the latter representing most likely true Dcr2 cleavage products). Combined with the importance of the Dcr2 helicase domain, this suggests that the majority of 21 nt vsiRNAs originate by processive cleavage. This study sheds new light on Ae. aegypti Dcr2 functions and properties in this important arbovirus vector species.

Like other pathogens with high mutation and replication rates, within-host dengue virus (DENV) populations evolve during infection of their main mosquito vector, Aedes aegypti. Within-host DENV evolution during transmission provides opportunities for adaptation and emergence of novel virus variants. Recent studies of DENV genetic diversity failed to detect convergent evolution of adaptive mutations in mosquito tissues such as midgut and salivary glands, suggesting that convergent positive selection is not a major driver of within-host DENV evolution in the vector. However, it is unknown whether this conclusion extends to the transmitted viral subpopulation because it is technically difficult to sequence DENV genomes in mosquito saliva. Here, we achieved DENV full-genome sequencing by pooling saliva samples collected non-sacrificially from 49 to 163 individual Ae. aegypti mosquitoes previously infected with one of two DENV-1 genotypes. We compared the transmitted viral subpopulations found in the pooled saliva samples collected in time series with the input viral population present in the infectious blood meal. In all pooled saliva samples examined, the full-genome consensus sequence of the input viral population was unchanged. Although the pooling strategy prevents analysis of individual saliva samples, our results demonstrate the lack of strong convergent positive selection during a single round of DENV transmission by Ae. aegypti. This finding reinforces the idea that genetic drift and purifying selection are the dominant evolutionary forces shaping within-host DENV genetic diversity during transmission by mosquitoes.

Flaviviruses encompass not only medically relevant arthropod-borne viruses (arboviruses) but also insect-specific flaviviruses (ISFs) that are presumably maintained primarily through vertical transmission in the insect host. Interestingly, ISFs are commonly found infecting important arbovirus vectors such as the mosquito Aedes aegypti. Cell-fusing agent virus (CFAV) was the first described ISF of mosquitoes more than four decades ago. Despite evidence for widespread CFAV infections in A.aegypti populations and for CFAV potential to interfere with arbovirus transmission, little is known about CFAV evolutionary history. Here, we generated six novel CFAV genome sequences by sequencing three new virus isolates and subjecting three mosquito samples to untargeted viral metagenomics. We used these new genome sequences together with published ones to perform a global phylogenetic analysis of CFAV genetic diversity. Although there was some degree of geographical clustering among CFAV sequences, there were also notable discrepancies between geography and phylogeny. In particular, CFAV sequences from Cambodia and Thailand diverged significantly, despite confirmation that A.aegypti populations from both locations are genetically close. The apparent phylogenetic discrepancy between CFAV and its A.aegypti host in Southeast Asia indicates that other factors than host population structure shape CFAV genetic diversity.

Although specific interactions between host and pathogen genotypes have been well documented in invertebrates, the identification of host genes involved in discriminating pathogen genotypes remains a challenge. In the mosquito Aedes aegypti, the main dengue virus (DENV) vector worldwide, statistical associations between host genetic markers and DENV types or strains were previously detected, but the host genes underlying this genetic specificity have not been identified. In particular, it is unknown whether DENV type- or strain-specific resistance relies on allelic variants of the same genes or on distinct gene sets. Here, we investigated the genetic architecture of DENV resistance in a population of Ae. aegypti from Bakoumba, Gabon, which displays a stronger resistance phenotype to DENV type 1 (DENV-1) than to DENV type 3 (DENV-3) infection. Following experimental exposure to either DENV-1 or DENV-3, we sequenced the exomes of large phenotypic pools of mosquitoes that are either resistant or susceptible to each DENV type. Using variation in single-nucleotide polymorphism (SNP) frequencies among the pools, we computed empirical p values based on average gene scores adjusted for the differences in SNP counts, to identify genes associated with infection in a DENV type-specific manner. Among the top 5% most significant genes, 263 genes were significantly associated with resistance to both DENV-1 and DENV-3, 287 genes were only associated with DENV-1 resistance and 290 were only associated with DENV-3 resistance. The shared significant genes were enriched in genes with ATP binding activity and sulfur compound transmembrane transporter activity, whereas the genes uniquely associated with DENV-3 resistance were enriched in genes with zinc ion binding activity. Together, these results indicate that specific resistance to different DENV types relies on largely non-overlapping sets of genes in this Ae. aegypti population and pave the way for further mechanistic studies.

Dengue virus (DENV) causes more human infections than any other mosquito-borne virus. The current lack of antiviral strategies has prompted genome-wide screens for host genes that are required for DENV infectivity. Earlier transcriptomic studies that identified DENV host factors in the primary vector Aedes aegypti used inbred laboratory colonies and/or pools of mosquitoes that erase individual variation. Here, we performed transcriptome sequencing on individual midguts in a field-derived Ae. aegypti population to identify new candidate host factors modulating DENV replication. We analyzed the transcriptomic data using an approach that accounts for individual co-variation between viral RNA load and gene expression. This approach generates a prediction about the agonist or antagonist effect of candidate genes on DENV replication based on the sign of the correlation between gene expression and viral RNA load. Using this method, we identified 39 candidate genes that went undetected by conventional pairwise comparison of gene expression levels between DENV-infected midguts and uninfected controls. Only four candidate genes were detected by both methods, emphasizing their complementarity. We demonstrated the value of our approach by functional validation of a candidate agonist gene encoding a sterol regulatory element-binding protein (SREBP), which was identified by correlation analysis but not by pairwise comparison. We confirmed that SREBP promotes DENV infection in the midgut by RNAi-mediated gene knockdown in vivo. We suggest that our approach for transcriptomic analysis can empower genome-wide screens for potential agonist or antagonist factors by leveraging inter-individual variation in gene expression. More generally, this method is applicable to a wide range of phenotypic traits displaying inter-individual variation.

Female mosquitoes need a blood meal to reproduce, and in obtaining this essential nutrient they transmit deadly pathogens. Although crucial for the spread of mosquito-borne diseases, blood feeding remains poorly understood due to technological limitations. Indeed, studies often expose human subjects to assess biting behavior. Here, we present the biteOscope, a device that attracts mosquitoes to a host mimic which they bite to obtain an artificial blood meal. The host mimic is transparent, allowing high-resolution imaging of the feeding mosquito. Using machine learning, we extract detailed behavioral statistics describing the locomotion, pose, biting, and feeding dynamics of Aedes aegypti, Aedes albopictus, Anopheles stephensi, and Anopheles coluzzii. In addition to characterizing behavioral patterns, we discover that the common insect repellent DEET repels Anopheles coluzzii upon contact with their legs. The biteOscope provides a new perspective on mosquito blood feeding, enabling the high-throughput quantitative characterization of this lethal behavior.