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 Asian tiger mosquito Aedes albopictus is globally expanding and has become the main vector for human arboviruses in Europe. Here we present AalbF2, a dramatically improved assembly of the Ae. albopictus genome that has revealed widespread viral insertions, novel microRNAs and piRNA clusters, the sex determining locus, new immunity genes, and has enabled genome-wide studies of geographically diverse Ae. albopictus populations and analyses of the developmental and stage-dependent network of expression data. Additionally, we built the first physical map for this species with 75% of the assembled genome anchored to the chromosomes. These up-to-date resources of the genome provide a foundation to improve understanding of the adaptation potential and the epidemiological relevance of this species and foster the development of innovative control measures.

This is an open access article under the terms of the Creative Commons Attribution-NonCommercial-NoDerivs License, which permits use and distribution in any medium, provided the original work is properly cited, the use is non-commercial and no modifications or adaptations are made.

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.

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.

Background Several bioinformatics pipelines have been developed to detect sequences from viruses that integrate into the human genome because of the health relevance of these integrations, such as in the persistence of viral infection and/or in generating genotoxic effects, often progressing into cancer. Recent genomics and metagenomics analyses have shown that viruses also integrate into the genome of non-model organisms (i.e., arthropods, fish, plants, vertebrates). However, rarely studies of endogenous viral elements (EVEs) in non-model organisms have gone beyond their characterization from reference genome assemblies. In non-model organisms, we lack a thorough understanding of the widespread occurrence of EVEs and their biological relevance, apart from sporadic cases which nevertheless point to significant roles of EVEs in immunity and regulation of expression. The concomitance of repetitive DNA, duplications and/or assembly fragmentations in a genome sequence and intrasample variability in whole-genome sequencing (WGS) data could determine misalignments when mapping data to a genome assembly. This phenomenon hinders our ability to properly identify integration sites. Results To fill this gap, we developed ViR, a pipeline which solves the dispersion of reads due to intrasample variability in sequencing data from both single and pooled DNA samples thus ameliorating the detection of integration sites. We tested ViR to work with both in silico and real sequencing data from a non-model organism, the arboviral vector Aedes albopictus. Potential viral integrations predicted by ViR were molecularly validated supporting the accuracy of ViR results. Conclusion ViR will open new venues to explore the biology of EVEs, especially in non-model organisms. Importantly, while we generated ViR with the identification of EVEs in mind, its application can be extended to detect any lateral transfer event providing an ad-hoc sequence to interrogate.

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.

In the model organism Drosophila melanogaster, the PIWI-interacing RNA pathway contributes in silencing transposable elements (TEs) through smallRNAs (piRNAs), which arise from genomic loci (piRNA clusters) that contain sequences of previously-acquired TEs. As such, they are a functionally-immune archive of previous TE invasions that is passed to the offspring. In the arboviral vector Aedes aegypti, piRNA clusters contain TEs and endogenous viral elements from nonretroviral RNA viruses (nrEVEs) which produce piRNAs, supporting the hypothesis that nrEVEs are heritable immunity effectors. However, direct evidence that nrEVEs mediate adaptive immunity is lacking. Here, by using an analytic approach intersecting population genomics with molecular biology we demonstrate that the composition of piRNA clusters is modular through acquisition and absence of nrEVEs. We show that the genomes of wild-caught mosquitoes have a different set of nrEVEs than those annotated in the reference genome, including population-specific integrations. nrEVEs are not distributed in mosquito genomes only by genetic drift, but some show signs of positive selection. Moreover, by comparing natural mosquito populations expressing or lacking two newly characterised nrEVEs with high sequence complementarity to cell fusing agent virus, we show that nrEVEs confer antiviral immunity in ovaries against the cognate virus. Our results confirm that some nrEVEs have been co-opted for adaptive immunity to viral infections. form an "archive" in the genome (Marraffini and Sontheimer, 2010; Koonin and Makarova 2017; Kofler 2019). CRISPR RNAs (Clustered Regularly Interspaced Palindromic Repeats, crRNAs) are produced from short fragments of foreign nucleic acids previously integrated into the repetitive CRISPR locus of the host genome. crRNAs assemble with specific CRISPR-associated proteins (Cas) proteins to form complexes that bind and cleave any phage or plasmid bearing sequence complementarity to the crRNA (Amitai and Sorek et al., 2016).piRNAs are a class of small RNAs of around 25-30 nt that guide PIWI proteins onto complementary target RNAs, resulting in gene silencing at the post-transcriptional or transcriptional level (Ozata et al., 2019). In D. melanogaster, piRNA precursors are generated from genomic loci, called piRNA clusters (Czech and Hannon, 2016). These regions are enriched for (remnants of) TE sequences and, as a consequence, cluster-derived piRNAs show sequence complementarity to TEs. It is thought that piRNA clusters act as traps for new TE invasions by horizontal transfer (Brennecke et al., 2007; Khurana et al., 2011;Parhad and Theurkauf, 2019). Thus, akin to CRISPR loci in prokaryotes for