Improved reference genome of Aedes aegypti informs arbovirus vector control.

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.

Peri-domestic entomological surveillance using private traps allows detection of dengue virus in Aedes albopictus during an autochthonous transmission event in mainland France, late summer 2023.

While locally-acquired dengue virus (DENV) human infections occur in mainland France since 2010, data to identify the mosquito species involved and to trace the virus are frequently lacking. Supported by a local network gathering public health agencies and research laboratories, we analysed, in late summer 2023, mosquitoes from privately-owned traps within a French urban neighbourhood affected by a dengue cluster. The cluster, in Auvergne-Rhône-Alpes, comprised three cases, including two autochthonous ones. Upon return from a recent visit to the French Caribbean Islands, the third case had consulted healthcare because of dengue-compatible symptoms, but dengue had not been recognised. For the two autochthonous cases, DENV-specific antibodies in serum or a positive quantitative PCR for DENV confirmed DENV infection. The third case had anti-flavivirus IgMs. No DENV genetic sequences were obtained from affected individuals but Aedes albopictus mosquitoes trapped less than 200 m from the autochthonous cases' residence contained DENV. Genetic data from the mosquito-derived DENV linked the cluster to the 2023-2024 dengue outbreak in the French Caribbean Islands. This study highlights the importance of raising mosquito-borne disease awareness among healthcare professionals. It demonstrates Ae. albopictus as a DENV vector in mainland France and the value of private mosquito traps for entomo-virological surveillance.

Molecular epidemiology identifies the expansion of the DENV2 epidemic lineage from the French Caribbean Islands to French Guiana and mainland France, 2023 to 2024.

In 2023, dengue virus serotype 2 (DENV2) affected most French overseas territories. In the French Caribbean Islands, viral circulation continues with > 30,000 suspected infections by March 2024. Genome sequence analysis reveals that the epidemic lineage in the French Caribbean islands has also become established in French Guiana but not Réunion. It has moreover seeded autochthonous circulation events in mainland France. To guide prevention of further inter-territorial spread and DENV introduction in non-endemic settings, continued molecular surveillance and mosquito control are essential.

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.

Modeling intra-mosquito dynamics of Zika virus and its dose-dependence confirms the low epidemic potential of Aedes albopictus.

Originating from African forests, Zika virus (ZIKV) has now emerged worldwide in urbanized areas, mainly transmitted by Aedes aegypti mosquitoes. Although Aedes albopictus can transmit ZIKV experimentally and was suspected to be a ZIKV vector in Central Africa, the potential of this species to sustain virus transmission was yet to be uncovered until the end of 2019, when several autochthonous transmissions of the virus vectored by Ae. albopictus occurred in France. Aside from these few locally acquired ZIKV infections, most territories colonized by Ae. albopictus have been spared so far. The risk level of ZIKV emergence in these areas remains however an open question. To assess Ae. albopictus' vector potential for ZIKV and identify key virus outbreak predictors, we built a complete framework using the complementary combination of (i) dose-dependent experimental Ae. albopictus exposure to ZIKV followed by time-dependent assessment of infection and systemic infection rates, (ii) modeling of intra-human ZIKV viremia dynamics, and (iii) in silico epidemiological simulations using an Agent-Based Model. The highest risk of transmission occurred during the pre-symptomatic stage of the disease, at the peak of viremia. At this dose, mosquito infection probability was estimated to be 20%, and 21 days were required to reach the median systemic infection rates. Mosquito population origin, either temperate or tropical, had no impact on infection rates or intra-host virus dynamic. Despite these unfavorable characteristics for transmission, Ae. albopictus was still able to trigger and yield large outbreaks in a simulated environment in the presence of sufficiently high mosquito biting rates. Our results reveal a low but existing epidemic potential of Ae. albopictus for ZIKV, that might explain the absence of large scale ZIKV epidemics so far in territories occupied only by Ae. albopictus. They nevertheless support active surveillance and eradication programs in these territories to maintain the risk of emergence to a low level.

Interactions between vector competence to chikungunya virus and resistance to deltamethrin in Aedes aegypti laboratory lines?

The urban mosquito species Aedes aegypti is the main vector of arboviruses worldwide. Mosquito control with insecticides is the most prevalent method for preventing transmission in the absence of effective vaccines and available treatments; however, the extensive use of insecticides has led to the development of resistance in mosquito populations throughout the world, and the number of epidemics caused by arboviruses has increased. Three mosquito lines with different resistance profiles to deltamethrin were isolated in French Guiana, including one with the I1016 knock-down resistant allele. Significant differences were observed in the cumulative proportion of mosquitoes with a disseminated chikungunya virus infection over time across these lines. In addition, some genes related to resistance (CYP6BB2, CYP6N12, GST2, trypsin) were variably overexpressed in the midgut at 7 days after an infectious bloodmeal in these three lines. Our work shows that vector competence for chikungunya virus varied between Ae. aegypti laboratory lines with different deltamethrin resistance profiles. More accurate verification of the functional association between insecticide resistance and vector competence remains to be demonstrated.

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.

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.

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.

Three-way interactions between mosquito population, viral strain and temperature underlying chikungunya virus transmission potential.

Interactions between pathogens and their insect vectors in nature are under the control of both genetic and non-genetic factors, yet most studies on mosquito vector competence for human pathogens are conducted in laboratory systems that do not consider genetic and/or environmental variability. Evaluating the risk of emergence of arthropod-borne viruses (arboviruses) of public health importance such as chikungunya virus (CHIKV) requires a more realistic appraisal of genetic and environmental contributions to vector competence. In particular, sources of variation do not necessarily act independently and may combine in the form of interactions. Here, we measured CHIKV transmission potential by the mosquito Aedes albopictus in all combinations of six worldwide vector populations, two virus strains and two ambient temperatures (20°C and 28°C). Overall, CHIKV transmission potential by Ae. albopictus strongly depended on the three-way combination of mosquito population, virus strain and temperature. Such genotype-by-genotype-by-environment (G × G × E) interactions question the relevance of vector competence studies conducted with a simpler set of conditions. Our results highlight the need to account for the complex interplay between vectors, pathogens and environmental factors to accurately assess the potential of vector-borne diseases to emerge.

Impact on Aedes aegypti Mosquitoes Exposed to Honey-Impregnated Flinders Technology Associates (FTA®) Cards.

Programs to control viruses transmitted by mosquitoes requires the implementation of surveillance tools. Over the past decade, Flinders Technology Associates (FTA®) cards, which preserve nucleic acids, have emerged as an innovating surveillance system for collecting arboviruses expectorated during mosquito sugar feeding. In this study, we evaluate the survival rates of two strain of Aedes aegypti (New Orleans (NO) and Cayenne (CAY)) in the laboratory after exposing to either honey-impregnated FTA® cards or untreated filter paper (UFP) card. Experimental exposure of mosquitoes to FTA® cards during sugar feeding significantly negatively impacted their survival, as compared to untreated filter paper. The median survival time was 2 days (95% confidence interval [CI] 1 day, 3 days) for mosquitoes exposed to FTA cards from strain NO and 3 days (95% CI 2 days, 3 days) for mosquitoes exposed to FTA cards from strain CAY. Mosquitoes exposed to UFP did not survive until the end of the experiment (4 days for strain NO and 7 days for strain CAY). Although this finding does not preclude the use of FTA® cards in surveillance, it is crucial to consider and incorporate this factor into study designs.

Sand flies and Toscana virus: Intra-vector infection dynamics and impact on Phlebotomus perniciosus life-history traits.

Toscana virus (TOSV) is a leading cause of summer viral meningitis in Southern Europe (Central Italy, south of France, Spain and Portugal) and can cause severe neurological cases. Within the Mediterranean basin, it is transmitted by hematophagous sand flies belonging to the Phlebotomus genus. Despite the identification of the primary TOSV vectors, the viral developmental cycle in vector species remains largely unknown. Limited research has been conducted on transmission dynamics and the vector competence and vectorial capacity of the principal TOSV vector, Phlebotomus perniciosus. In this context, we investigated the intra-vector TOSV infection dynamics in Ph. perniciosus, as well as its impact on the vector life history traits. Female sand flies were experimentally infected with TOSV through an artificial blood meal. Systemic dissemination of the virus was observed approximately three days post-infection, potentially resulting in a short extrinsic incubation period. Moreover, the study revealed a longer hatching time for eggs laid by infected females. This research brought additional experimental insights regarding the vector competence of Ph. perniciosus but also provided the first insight into TOSV developmental cycle and its impact on the vector. These findings prompt further exploration of TOSV transmission dynamics, raise new hypotheses on the virus transmission and highlight the importance of follow-up studies.

ß-D-N4-hydroxycytidine (NHC, EIDD-1931) inhibits chikungunya virus replication in mosquito cells and ex vivo Aedes aegypti guts, but not when ingested during blood-feeding.

Chikungunya virus (CHIKV) is a mosquito-borne virus transmitted by Aedes mosquitoes. While there are no antiviral therapies currently available to treat CHIKV infections, several licensed oral drugs have shown significant anti-CHIKV activity in cells and in mouse models. However, the efficacy in mosquitoes has not yet been assessed. Such cross-species antiviral activity could be favorable, since virus inhibition in the mosquito vector might prevent further transmission to vertebrate hosts. Here, we explored the antiviral effect of ß-d-N4-hydroxycytidine (NHC, EIDD-1931), the active metabolite of molnupiravir, on CHIKV replication in Aedes aegypti mosquitoes. Antiviral assays in mosquito cells and in ex vivo cultured mosquito guts showed that NHC had significant antiviral activity against CHIKV. Exposure to a clinically relevant concentration of NHC did not affect Ae. aegypti lifespan when delivered via a bloodmeal, but it slightly reduced the number of eggs developed in the ovaries. When mosquitoes were exposed to a blood meal containing both CHIKV and NHC, the compound did not significantly reduce virus infection and dissemination in the mosquitoes. This was confirmed by modelling and could be explained by pharmacokinetic analysis, which revealed that by 6 h post-blood-feeding, 90% of NHC had been cleared from the mosquito bodies. Our data show that NHC inhibited CHIKV replication in mosquito cells and gut tissue, but not in vivo when mosquitoes were provided with a CHIKV-infectious bloodmeal spiked with NHC. The pipeline presented in this study offers a suitable approach to identify anti-arboviral drugs that may impede replication in mosquitoes.

Genomic surveillance reveals a dengue 2 virus epidemic lineage with a marked decrease in sensitivity to Mosnodenvir.

Dengue fever is the most important arbovirosis for public health, with more than 5 million cases worldwide in 2023. Mosnodenvir is the first anti-dengue compound with very high preclinical pan-serotype activity, currently undergoing phase 2 clinical evaluation. Here, by analyzing dengue virus (DENV) genomes from the 2023-2024 epidemic in the French Caribbean Islands, we show that they all exhibit mutation NS4B:V91A, previously associated with a marked decrease in sensitivity to mosnodenvir in vitro. Using antiviral activity tests on four clinical and reverse-genetic strains, we confirm a marked decrease in mosnodenvir sensitivity for DENV-2 ( > 1000 fold). Finally, combining phylogenetic analysis and experimental testing for resistance, we find that virus lineages with low sensitivity to mosnodenvir due to the V91A mutation likely emerged multiple times over the last 30 years in DENV-2 and DENV-3. These results call for increased genomic surveillance, in particular to track lineages with resistance mutations. These efforts should allow to better assess the activity profile of DENV treatments in development against circulating strains.

Noninvasive detection of Zika virus in mosquito excreta sampled from wild mosquito populations in French Guiana.

Arboviruses can be difficult to detect in the field due to relatively low prevalence in mosquito populations. The discovery that infected mosquitoes can release viruses in both their saliva and excreta gave rise to low-cost methods for the detection of arboviruses during entomological surveillance. We implemented both saliva and excreta-based entomological surveillance during the emergence of Zika virus (ZIKV) in French Guiana in 2016 by trapping mosquitoes around households of symptomatic cases with confirmed ZIKV infection. ZIKV was detected in mosquito excreta and not in mosquito saliva in 1 trap collection out of 85 (1.2%). One female Ae. aegypti L. (Diptera: Culicidae) was found with a ZIKV systemic infection in the corresponding trap. The lag time between symptom onset in a ZIKV-infected individual living near the trap site and ZIKV detection in this mosquito was 1 wk. These results highlight the potential of detection in excreta from trapped mosquitoes as a sensitive and cost-effective method to non invasively detect arbovirus circulation.

Updated estimation of cutaneous leishmaniasis incubation period in French Guiana.

BACKGROUND: The cutaneous leishmaniasis (CL) incubation period (IP) is defined as the time between parasite inoculation by sandfly bite and the onset of the first CL lesion. IP distribution is difficult to assess for CL because the date of exposure to an infectious bite cannot be accurately determined in endemic areas. IP current estimates for CL range from 14 days to several months with a median around 30-60 days, as established by a few previous studies in both New and Old Worlds. METHODOLOGY: We estimated CL incubation period distribution using time-to-event models adapted to interval-censored data based on declared date of travels from symptomatic military personnel living in non-endemic areas that were exposed during their short stays in French Guiana (FG) between January 2001 and December 2021. PRINCIPAL FINDINGS: A total of 180 patients were included, of which 176 were men (97.8%), with a median age of 26 years. When recorded, the parasite species was always Leishmania guyanensis (31/180, 17.2%). The main periods of CL diagnosis spread from November to January (84/180, 46.7%) and over March-April (54/180, 30.0%). The median IP was estimated at 26.2 days (95% Credible Level, 23.8-28.7 days) using a Bayesian accelerated failure-time regression model. Estimated IP did not exceed 62.1 days (95% CI, 56-69.8 days) in 95% of cases (95th percentile). Age, gender, lesion number, lesion evolution and infection date did not significantly modify the IP. However, disseminated CL was significantly associated with a 2.8-fold shortening of IP. CONCLUSIONS: This work suggests that the CL IP distribution in French Guiana is shorter and more restricted than anticipated. As the incidence of CL in FG usually peaks in January and March, these findings suggest that patients are contaminated at the start of the rainy season.

Combining two genetic sexing strains allows sorting of non-transgenic males for Aedes genetic control.

Chemical control of disease vectoring mosquitoes Aedes albopictus and Aedes aegypti is costly, unsustainable, and increasingly ineffective due to the spread of insecticide resistance. The Sterile Insect Technique is a valuable alternative but is limited by slow, error-prone, and wasteful sex-separation methods. Here, we present four Genetic Sexing Strains (two for each Aedes species) based on fluorescence markers linked to the m and M sex loci, allowing for the isolation of transgenic males. Furthermore, we demonstrate how combining these sexing strains enables the production of non-transgenic males. In a mass-rearing facility, 100,000 first instar male larvae could be sorted in under 1.5 h with an estimated 0.01-0.1% female contamination on a single machine. Cost-efficiency analyses revealed that using these strains could result in important savings while setting up and running a mass-rearing facility. Altogether, these Genetic Sexing Strains should enable a major upscaling in control programmes against these important vectors.

Analysis of trapped mosquito excreta as a noninvasive method to reveal biodiversity and arbovirus circulation.

Emerging and endemic mosquito-borne viruses can be difficult to detect and monitor because they often cause asymptomatic infections in human or vertebrate animals or cause nonspecific febrile illness with a short recovery waiting period. Some of these pathogens circulate into complex cryptic cycles involving several animal species as reservoir or amplifying hosts. Detection of cases in vertebrate hosts can be complemented by entomological surveillance, but this method is not adapted to low infection rates in mosquito populations that typically occur in low or nonendemic areas. We identified West Nile virus circulation in Camargue, a wetland area in South of France, using a cost-effective xenomonitoring method based on the molecular detection of virus in excreta from trapped mosquitoes. We also succeeded at identifying the mosquito species community on several sampling sites, together with the vertebrate hosts on which they fed prior to being captured using amplicon-based metabarcoding on mosquito excreta without processing any mosquitoes. Mosquito excreta-based virus surveillance can complement standard surveillance methods because it is cost-effective and does not require personnel with a strong background in entomology. This strategy can also be used to noninvasively explore the ecological network underlying arbovirus circulation.

Anopheles salivary gland proteomes from major malaria vectors.

BACKGROUND: Antibody responses against Anopheles salivary proteins can indicate individual exposure to bites of malaria vectors. The extent to which these salivary proteins are species-specific is not entirely resolved. Thus, a better knowledge of the diversity among salivary protein repertoires from various malaria vector species is necessary to select relevant genus-, subgenus- and/or species-specific salivary antigens. Such antigens could be used for quantitative (mosquito density) and qualitative (mosquito species) immunological evaluation of malaria vectors/host contact. In this study, salivary gland protein repertoires (sialomes) from several Anopheles species were compared using in silico analysis and proteomics. The antigenic diversity of salivary gland proteins among different Anopheles species was also examined. RESULTS: In silico analysis of secreted salivary gland protein sequences retrieved from an NCBInr database of six Anopheles species belonging to the Cellia subgenus (An. gambiae, An. arabiensis, An. stephensi and An. funestus) and Nyssorhynchus subgenus (An. albimanus and An. darlingi) displayed a higher degree of similarity compared to salivary proteins from closely related Anopheles species. Additionally, computational hierarchical clustering allowed identification of genus-, subgenus- and species-specific salivary proteins. Proteomic and immunoblot analyses performed on salivary gland extracts from four Anopheles species (An. gambiae, An. arabiensis, An. stephensi and An. albimanus) indicated that heterogeneity of the salivary proteome and antigenic proteins was lower among closely related anopheline species and increased with phylogenetic distance. CONCLUSION: This is the first report on the diversity of the salivary protein repertoire among species from the Anopheles genus at the protein level. This work demonstrates that a molecular diversity is exhibited among salivary proteins from closely related species despite their common pharmacological activities. The involvement of these proteins as antigenic candidates for genus-, subgenus- or species-specific immunological evaluation of individual exposure to Anopheles bites is discussed.