Background: The Asian tiger mosquito Aedes albopictus is globally expanding and has become the main vector for human arboviruses in Europe. With limited antiviral drugs and vaccines available, vector control is the primary approach to prevent mosquito-borne diseases. A reliable and accurate DNA sequence of the Ae. albopictus genome is essential to develop new approaches that involve genetic manipulation of mosquitoes. Results: We use long-read sequencing methods and modern scaffolding techniques (PacBio, 10X, and Hi-C) to produce AalbF2, a dramatically improved assembly of the Ae. albopictus genome. AalbF2 reveals widespread viral insertions, novel microRNAs and piRNA clusters, the sex-determining locus, and new immunity genes, and enables genome-wide studies of geographically diverse Ae. albopictus populations and analyses of the developmental and stage-dependent network of expression data. Additionally, we build the first physical map for this species with 75% of the assembled genome anchored to the chromosomes. Conclusion: The AalbF2 genome assembly represents the most up-to-date collective knowledge of the Ae. albopictus genome. These resources represent a foundation to improve understanding of the adaptation potential and the epidemiological relevance of this species and foster the development of innovative control measures.

The mosquitoAedes albopictusis one of the most dangerous invasive species. Its worldwide spread has created health concerns as it is a major vector of arboviruses of public health significance such as chikungunya (CHIKV). Dynamics of different genetic backgrounds and admixture events may have impacted competence for CHIKV in adventive populations. Using microsatellites, we infer the genetic structure of populations across the expansion areas that we then associate with their competence for different CHIKV genotypes. Here we show that the demographic history ofAe. albopictuspopulations is a consequence of rapid complex patterns of historical lineage diversification and divergence that influenced their competence for CHIKV. The history of adventive populations is associated with CHIKV genotypes in a genotype-by-genotype interaction that impacts their vector competence. Thus, knowledge of the demographic history and vector competence of invasive mosquitoes is pivotal for assessing the risk of arbovirus outbreaks in newly colonized areas.

Horizontal gene transfer from viruses to eukaryotic cells is a pervasive phenomenon. Somatic viral integrations are linked to persistent viral infection whereas integrations into germline cells are maintained in host genomes by vertical transmission and may be co‐opted for host functions. In the arboviral vector Aedes aegypti, an endogenous viral element from a nonretroviral RNA virus (nrEVE) was shown to produce PIWI‐interacting RNAs (piRNAs) to limit infection with a cognate virus. Thus, nrEVEs may constitute a heritable, sequence‐specific mechanism for antiviral immunity, analogous to piRNA‐mediated silencing of transposable elements. Here, we combine population genomics and evolutionary approaches to analyse the genomic architecture of nrEVEs in A. aegypti. We conducted a genome‐wide screen for adaptive nrEVEs and searched for novel population‐specific nrEVEs in the genomes of 80 individual wild‐caught mosquitoes from five geographical populations. We show a dynamic landscape of nrEVEs in mosquito genomes and identified five novel nrEVEs derived from two currently circulating viruses, providing evidence of the environmental‐dependent modification of a piRNA cluster. Overall, our results show that virus endogenization events are complex with only a few nrEVEs contributing to adaptive evolution in A. aegypti.

The arboviral vector Aedes albopictus became established on all continents except Antarctica in the past 50 years. A consequence of its rapid global invasion is the transmission of diseases previously confined to the tropics and subtropics occurring in temperate regions of the world, including the re-emergence of chikungunya and dengue in Europe. Application of pyrethroids is among the most widely-used interventions for vector control, especially in the presence of an arboviral outbreak. Studies are emerging that reveal phenotypic resistance and monitor mutations at the target site, the para sodium channel gene, primarily on a local scale. A total of 512 Ae. albopictus mosquitoes from twelve geographic sites, including those from the native home range and invaded areas, were sampled between 2011 and 2018, and were analyzed at five codons of the para sodium channel gene with mutations predictive of resistance phenotype. Additionally, to test for the origin of unique kdr mutations in Mexico, we analyzed the genetic connectivity of southern Mexico mosquitoes with mosquitoes from home range, the Reunion Island, America and Europe. We detected mutations at all tested positions of the para sodium channel gene, with heterozygotes predominating and rare instance of double mutants. We observed an increase in the distribution and frequency of F1534C/L/S mutations in the ancestral China population and populations in the Mediterranean Greece, the appearance of the V1016G/I mutations as early as 2011 in Italy and mutations at position 410 and 989 in Mexico. The analyses of the distribution pattern of kdr alleles and haplotype network analyses showed evidence for multiple origins of all kdr mutations. Here we provide the most-up-to-date survey on the geographic and temporal distribution of pyrethroid-predictive mutations in Ae. albopictus by combining kdr genotyping on current and historical samples with published data. While we confirm low levels of pyrethroid resistance in most analyzed samples, we find increasing frequencies of F1534C/S and V1016G in China and Greece or Italy, respectively. The observed patterns of kdr allele distribution support the hypothesis that on site emergence of resistance has contributed more than spread of resistance through mosquito migration/invasions to the current widespread of kdr alleles, emphasizing the importance of local surveillance programs and resistance management. Aedes albopictus is a highly invasive mosquito and a vector for a number of arboviruses. The arrival and establishment of Ae. albopictus in temperature regions of the world, such as Europe, has been accompanied by re-emergence of arboviral diseases. There are no effective therapeutic treatments for arboviruses meaning control of vector populations is the primary strategy to prevent arboviral disease transmission. Pyrethroids are frequently used for control of vectors based on their low mammalian toxicity and rapid knockdown effect on mosquitoes. The identification of mutations predictive of resistance phenotype in the para sodium channel gene, the target site of pyrethroids, has provided for molecular markers to test for resistance by genotyping wild-collected mosquitoes. Here we analysed all currently known predictive mutations for pyrethroid resistance in 512 geographic mosquitoes sampled in a span of seven years. Thus, we are able to add a temporal dimension to the analysis of geographic distribution of resistance alleles. Overall our data confirm a patchy distribution of mutations predictive of the resistance phenotype, but also reveal an alarming increase of resistance mutations in China, Greece and Italy. Resistance mutations appear to have arisen locally more than being spread through mosquito migration/invasions.

Current knowledge of the piRNA pathway is based mainly on studies on Drosophila melanogaster where three proteins of the Piwi subclade of the Argonaute family interact with PIWI-interacting RNAs to silence transposable elements in gonadal tissues. In mosquito species that transmit epidemic arboviruses such as dengue and chikungunya viruses, Piwi clade genes underwent expansion, are also expressed in the soma and cross-talk with proteins of recognized antiviral function cannot be excluded for some Piwi proteins. These observations underscore the importance of expanding our knowledge of the piRNA pathway beyond the model organism D. melanogaster. Here we focus on the emerging arboviral vector Aedes albopictus and we couple traditional approaches of expression and adaptive evolution analyses with most current computational predictions of protein structure to study evolutionary divergence among Piwi clade proteins. Superposition of protein homology models indicate possible high structure similarity among all Piwi proteins, with high levels of amino acid conservation in the inner regions devoted to RNA binding. On the contrary, solvent-exposed surfaces showed low conservation, with several sites under positive selection. Analysis of the expression profiles of Piwi transcripts during mosquito development and following infection with dengue serotype 1 or chikungunya viruses showed a concerted elicitation of all Piwi transcripts during viral dissemination of dengue viruses while maintenance of infection relied on expression of primarily Piwi5. Opposite, establishment of persistent infection by chikungunya virus is accompanied by increased expression of all Piwi genes, particularly Piwi4 and, again, Piwi5. Overall these results are consistent with functional specialization and a general antiviral role for Piwi5. Experimental evidences of sites under positive selection in Piwi1/3, Piwi4 and Piwi6, that have complex expression profiles, provide useful knowledge to design tailored functional experiments. Argonautes are ancient proteins involved in many cellular processes, including innate immunity. Early in eukaryote evolution, Argonautes separated into Ago and Piwi clades, which maintain a dynamic evolutionary history with frequent duplications and losses. The use of Drosophila melanogaster as a model organism proved fundamental to understand the function of Argonautes. However, recent studies showed that the patterns and observations made in D. melanogaster, including the number of Argonautes, their expression profile and their function, are a rarity among Dipterans. In vectors of epidemic arboviruses such as dengue and chikungunya viruses, Piwi genes underwent expansion, are expressed in the soma, and some of them appear to have antiviral functions. Besides being an important basic question, the identification of which (and how) Piwi genes have antiviral functions may be used for the development of novel genetic-based strategies of vector control. Here we coupled population genetics models with computational predictions of protein structure and expression analyses to investigate the evolution and function of Piwi genes of the emerging vector Aedes albopictus. Our data support a general antiviral role for Piwi5. Instead, the detection of complex expression profiles with the presence of sites under positive selection in Piwi1/3, Piwi4 and Piwi6 requires tailored functional experiments to clarify their antiviral role.

The Asian tiger mosquito Aedes albopictus is contributing to the (re)-emergence of Chikungunya virus (CHIKV). To gain insights into the molecular underpinning of viral persistence, which renders a mosquito a life-long vector, we coupled small RNA and whole genome sequencing approaches on carcasses and ovaries of mosquitoes sampled 14 days post CHIKV infection and investigated the profile of small RNAs and the presence of vDNA fragments. Since Aedes genomes harbor nonretroviral Endogenous Viral Elements (nrEVEs) which confers tolerance to cognate viral infections in ovaries, we also tested whether nrEVEs are formed after CHIKV infection. We show that while small interfering (si)RNAs are evenly distributed along the full viral genome, PIWI-interacting (pi)RNAs mostly arise from a ~1000 bp window, from which a unique vDNA fragment is identified. CHIKV infection does not result in the formation of new nrEVEs, but piRNAs derived from existing nrEVEs correlate with differential expression of an endogenous transcript. These results demonstrate that all three RNAi pathways contribute to the homeostasis during the late stage of CHIKV infection, but in different ways, ranging from directly targeting the viral sequence to regulating the expression of mosquito transcripts and expand the role of nrEVEs beyond immunity against cognate viruses.