Experimental infections with Zika virus strains reveal high vector competence of Aedes albopictus and Aedes aegypti populations from Gabon (Central Africa) for the African virus lineage.

ABSTRACTThe two main Zika virus (ZIKV) vectors, Aedes albopictus and Aedes aegypti (invasive and native species, respectively), are present in Gabon (Central Africa). The aim of this study was to determine the entomological ZIKV risk associated with these mosquito species in Gabon by evaluating their vector competence for an African (i.e. representative of the endemic strains circulating in sub-Saharan Africa) and two Asian (i.e. representatives of exogenous epidemic strains that could be introduced) ZIKV strains. The transmission efficiency of one Ae. aegypti and two Ae. albopictus field-collected populations from Libreville and Franceville was assayed at day 7, 14 and 21 after experimental oral infection. The two mosquito species could transmit all three ZIKV strains already at day 7 post-infection, but transmission efficiency was higher for the African strain than the non-African strains (>60% versus <14%; incubation period of 14-21 days). The two mosquito species exhibited comparable vector competence for ZIKV, although the amount of viral particles (African strain) in saliva was significantly higher in Ae. albopictus than Ae. aegypti at day 14 post-infection. These findings suggest that overall, ZIKV risk in Gabon is mainly related to virus strains that circulate endemically across sub-Saharan Africa, although the transmission of non-African strains remain possible in case of introduction. Due to its high infestation indexes and ecological/geographical ranges, this risk appears mainly associated with Ae. albopictus. Vector surveillance and control methods against this invasive mosquito must be strengthened in the region to limit the risk of future outbreaks.

Non-retroviral Endogenous Viral Element Limits Cognate Virus Replication in Aedes aegypti Ovaries.

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.

Exome-wide association study reveals largely distinct gene sets underlying specific resistance to dengue virus types 1 and 3 in Aedes aegypti.

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.

Novel genome sequences of cell-fusing agent virus allow comparison of virus phylogeny with the genetic structure of Aedes aegypti populations.

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.

Tudor-SN Promotes Early Replication of Dengue Virus in the Aedes aegypti Midgut.

Diseases caused by mosquito-borne viruses have been on the rise for the last decades, and novel methods aiming to use laboratory-engineered mosquitoes that are incapable of carrying viruses have been developed to reduce pathogen transmission. This has stimulated efforts to identify optimal target genes that are naturally involved in mosquito antiviral defenses or required for viral replication. Here, we investigated the role of a member of the Tudor protein family, Tudor-SN, upon dengue virus infection in the mosquito Aedes aegypti. Tudor-SN knockdown reduced dengue virus replication in the midgut of Ae. aegypti females. In immunofluorescence assays, Tudor-SN localized to the nucleolus in both Ae. aegypti and Aedes albopictus cells. A reporter assay and small RNA profiling demonstrated that Tudor-SN was not required for RNA interference function in vivo. Collectively, these results defined a novel proviral role for Tudor-SN upon early dengue virus infection of the Ae. aegypti midgut.

Cell-Fusing Agent Virus Reduces Arbovirus Dissemination in Aedes aegypti Mosquitoes In Vivo.

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.

Larval Exposure to the Bacterial Insecticide Bti Enhances Dengue Virus Susceptibility of Adult Aedes aegypti Mosquitoes.

Understanding the interactions between pathogens sharing the same host can be complicated for holometabolous animals when larval and adult stages are exposed to distinct pathogens. In medically important insect vectors, the effect of pathogen exposure at the larval stage may influence susceptibility to human pathogens at the adult stage. We addressed this hypothesis in the mosquito Aedes aegypti, a major vector of arthropod-borne viruses (arboviruses), such as the dengue virus (DENV) and the chikungunya virus (CHIKV). We experimentally assessed the consequences of sub-lethal exposure to the bacterial pathogen Bacillus thuringiensis subsp. israelensis (Bti), during larval development, on arbovirus susceptibility at the adult stage in three Ae. aegypti strains that differ in their genetic resistance to Bti. We found that larval exposure to Bti significantly increased DENV susceptibility, but not CHIKV susceptibility, in the Bti-resistant strains. However, there was no major difference in the baseline arbovirus susceptibility between the Bti-resistant strains and their Bti-susceptible parental strain. Although the generality of our results remains to be tested with additional arbovirus strains, this study supports the idea that the outcome of an infection by a pathogen depends on other pathogens sharing the same host even when they do not affect the same life stage of the host. Our findings may also have implications for Bti as a mosquito biocontrol agent, indicating that the sub-optimal Bti efficacy may have counter-productive effects by increasing vector competence, at least for some combinations of arbovirus and mosquito strains.

Epidemiological significance of dengue virus genetic variation in mosquito infection dynamics.

The kinetics of arthropod-borne virus (arbovirus) transmission by their vectors have long been recognized as a powerful determinant of arbovirus epidemiology. The time interval between virus acquisition and transmission by the vector, termed extrinsic incubation period (EIP), combines with vector mortality rate and vector competence to determine the proportion of infected vectors that eventually become infectious. However, the dynamic nature of this process, and the amount of natural variation in transmission kinetics among arbovirus strains, are poorly documented empirically and are rarely considered in epidemiological models. Here, we combine newly generated empirical measurements in vivo and outbreak simulations in silico to assess the epidemiological significance of genetic variation in dengue virus (DENV) transmission kinetics by Aedes aegypti mosquitoes. We found significant variation in the dynamics of systemic mosquito infection, a proxy for EIP, among eight field-derived DENV isolates representing the worldwide diversity of recently circulating type 1 strains. Using a stochastic agent-based model to compute time-dependent individual transmission probabilities, we predict that the observed variation in systemic mosquito infection kinetics may drive significant differences in the probability of dengue outbreak and the number of human infections. Our results demonstrate that infection dynamics in mosquitoes vary among wild-type DENV isolates and that this variation potentially affects the risk and magnitude of dengue outbreaks. Our quantitative assessment of DENV genetic variation in transmission kinetics contributes to improve our understanding of heterogeneities in arbovirus epidemiological dynamics.

A Survey of Virus Recombination Uncovers Canonical Features of Artificial Chimeras Generated During Deep Sequencing Library Preparation.

Chimeric reads can be generated by in vitro recombination during the preparation of high-throughput sequencing libraries. Our attempt to detect biological recombination between the genomes of dengue virus (DENV; +ssRNA genome) and its mosquito host using the Illumina Nextera sequencing library preparation kit revealed that most, if not all, detected host-virus chimeras were artificial. Indeed, these chimeras were not more frequent than with control RNA from another species (a pillbug), which was never in contact with DENV RNA prior to the library preparation. The proportion of chimera types merely reflected those of the three species among sequencing reads. Chimeras were frequently characterized by the presence of 1-20 bp microhomology between recombining fragments. Within-species chimeras mostly involved fragments in opposite orientations and located less than 100 bp from each other in the parental genome. We found similar features in published datasets using two other viruses: Ebola virus (EBOV; -ssRNA genome) and a herpesvirus (dsDNA genome), both produced with the Illumina Nextera protocol. These canonical features suggest that artificial chimeras are generated by intra-molecular template switching of the DNA polymerase during the PCR step of the Nextera protocol. Finally, a published Illumina dataset using the Flock House virus (FHV; +ssRNA genome) generated with a protocol preventing artificial recombination revealed the presence of 1-10 bp microhomology motifs in FHV-FHV chimeras, but very few recombining fragments were in opposite orientations. Our analysis uncovered sequence features characterizing recombination breakpoints in short-read sequencing datasets, which can be helpful to evaluate the presence and extent of artificial recombination.

Individual co-variation between viral RNA load and gene expression reveals novel host factors during early dengue virus infection of the Aedes aegypti midgut.

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.

Carryover effects of larval exposure to different environmental bacteria drive adult trait variation in a mosquito vector.

Conditions experienced during larval development of holometabolous insects can affect adult traits, but whether differences in the bacterial communities of larval development sites contribute to variation in the ability of insect vectors to transmit human pathogens is unknown. We addressed this question in the mosquito Aedes aegypti, a major arbovirus vector breeding in both sylvatic and domestic habitats in Sub-Saharan Africa. Targeted metagenomics revealed differing bacterial communities in the water of natural breeding sites in Gabon. Experimental exposure to different native bacterial isolates during larval development resulted in significant differences in pupation rate and adult body size but not life span. Larval exposure to an Enterobacteriaceae isolate resulted in decreased antibacterial activity in adult hemolymph and reduced dengue virus dissemination titer. Together, these data provide the proof of concept that larval exposure to different bacteria can drive variation in adult traits underlying vectorial capacity. Our study establishes a functional link between larval ecology, environmental microbes, and adult phenotypic variation in a holometabolous insect vector.

Genetic Drift, Purifying Selection and Vector Genotype Shape Dengue Virus Intra-host Genetic Diversity in Mosquitoes.

Due to their error-prone replication, RNA viruses typically exist as a diverse population of closely related genomes, which is considered critical for their fitness and adaptive potential. Intra-host demographic fluctuations that stochastically reduce the effective size of viral populations are a challenge to maintaining genetic diversity during systemic host infection. Arthropod-borne viruses (arboviruses) traverse several anatomical barriers during infection of their arthropod vectors that are believed to impose population bottlenecks. These anatomical barriers have been associated with both maintenance of arboviral genetic diversity and alteration of the variant repertoire. Whether these patterns result from stochastic sampling (genetic drift) rather than natural selection, and/or from the influence of vector genetic heterogeneity has not been elucidated. Here, we used deep sequencing of full-length viral genomes to monitor the intra-host evolution of a wild-type dengue virus isolate during infection of several mosquito genetic backgrounds. We estimated a bottleneck size ranging from 5 to 42 founding viral genomes at initial midgut infection, irrespective of mosquito genotype, resulting in stochastic reshuffling of the variant repertoire. The observed level of genetic diversity increased following initial midgut infection but significantly differed between mosquito genetic backgrounds despite a similar initial bottleneck size. Natural selection was predominantly negative (purifying) during viral population expansion. Taken together, our results indicate that dengue virus intra-host genetic diversity in the mosquito vector is shaped by genetic drift and purifying selection, and point to a novel role for vector genetic factors in the genetic breadth of virus populations during infection. Identifying the evolutionary forces acting on arboviral populations within their arthropod vector provides novel insights into arbovirus evolution.

Excretion of dengue virus RNA by Aedes aegypti allows non-destructive monitoring of viral dissemination in individual mosquitoes.

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.