The Asian tiger mosquito Aedes albopictus is an invasive mosquito and a competent vector for public-health relevant arboviruses such as Chikungunya ( Alphavirus ), Dengue and Zika ( Flavivirus ) viruses. Unexpectedly, the sequencing of the genome of this mosquito revealed an unusually high number of integrated sequences with similarities to non-retroviral RNA viruses of the Flavivirus and Rhabdovirus genera. These Non-retroviral Integrated RNA Virus Sequences (NIRVS) are enriched in piRNA clusters and coding sequences and have been proposed to constitute novel mosquito immune factors. However, given the abundance of NIRVS and their variable viral origin, their relative biological roles remain unexplored. Here we used an analytical approach that intersects computational, evolutionary and molecular methods to study the genomic landscape of mosquito NIRVS. We demonstrate that NIRVS are differentially distributed across mosquito genomes, with a core set of seemingly the oldest integrations with similarity to Rhabdoviruses . Additionally, we compare the polymorphisms of NIRVS with respect to that of fast and slow-evolving genes within the Ae. albopictus genome. Overall, NIRVS appear to be less polymorphic than slow-evolving genes, with differences depending on whether they occur in intergenic regions or in piRNA clusters. Finally, two NIRVS that map within the coding sequences of genes annotated as Rhabdovirus RNA-dependent RNA polymerase and the nucleocapsid-encoding gene, respectively, are highly polymorphic and are expressed, suggesting exaptation possibly to enhance the mosquito’s antiviral responses. These results greatly advance our understanding of the complexity of the mosquito repeatome and the biology of viral integrations in mosquito genomes.

The sequenced genome of the arboviral vector mosquito, Aedes albopictus, is replete with repetitive DNA and it harbors an unusually large number of endogenous viral sequences, collectively called Nonretroviral Integrated RNA Virus Sequences (NIRVS). NIRVS are enriched both within protein-coding gene exons and PIWI-interacting RNA (piRNA) clusters, where they encode piRNAs. Based on these features, NIRVS have been proposed to function as novel mosquito antiviral immune factors. However, the relative importance and contributions of different NIRVS as functional antiviral elements and their mechanisms of action remain open questions.We apply an analytical approach that intersects computational, evolutionary and molecular methods to identify NIRVS most likely affecting mosquito immunity. Using this strategy, we show that NIRVS are a highly dynamic component of the Ae. albopictus repeatome, which nevertheless maintains a core set of seemingly the oldest NIRVS with similarity to Rhabdoviruses. Population-level polymorphism of NIRVS varies depending on whether they occur in intergenic regions, piRNA clusters or are part of gene exons. NIRVS from piRNA clusters are differentially widespread in diverse populations but conserved at the sequence level. This is consistent with the hypothesis that they act analogously to fragments of transposable elements in piRNA clusters and contribute to piRNA-based immunity. Among NIRVS from gene exons, AlbRha52 and AlbRha12 have the hallmarks of domestication as they are fixed across populations, stably expressed, and as polymorphic at the sequence level as fast-evolving genes. Overall these results support the hypothesis that NIRVS contribute to mosquito immunity, potentially through diverse modes of action.

No abstract

Recently, two invasive Aedes mosquito species, Ae. japonicus and Ae. koreicus, are circulating in several European countries posing potential health risks to humans and animals. Vector control is the main option to prevent mosquito-borne diseases, and an accurate genome sequence of these mosquitoes is essential to better understand their biology and to develop effective control strategies. Here, we present a de novo genome assembly of the Ae. japonicus (Ajap1) and Ae. koreicus (Akor1) based on a hybrid approach that combines PacBio long reads and Illumina short reads data. Their quality was ascertained using various metrics. Masking of repetitive elements, gene prediction and functional annotation was performed. Sequence analysis revealed a very high presence of repetitive DNA and, among others, thermal adaptation genes and insecticide-resistance genes. Through the RNAseq analysis of larvae and adults of Ae. koreicus and Ae. japonicus exposed to different temperatures (15 and 18 C) we also identified genes showing a differential temperature-dependent activation. The assembly of Akor1 and Ajap1 genomes constitutes the first updated collective knowledge of the genomes of both mosquito species, providing the possibility of understanding key mechanisms of their biology such as the ability to adapt to harsh climates and to develop insecticide-resistance mechanisms.