Improved reference genome of the arboviral vector Aedes albopictus

CRISPR-editing of the virus vector Aedes albopictus cell line C6/36, illustrated by prohibitin 2 gene knockout

Aedes mosquitoes are important virus vectors. We provide a toolkit for CRISPR-Cas9-editing of difficult-to-knockdown gene previously shown to be refractory to siRNA silencing in mosquito cells, which is pivotal in understanding vector biology, vector competence, host-pathogen interactions and in gene annotations. Starting from database searches of Ae. albopictus and the C6/36 cell line whole genome shotgun sequences for the prohibitin 2 (PHB2) gene, primers were designed to confirm the gene sequence in our laboratory-passaged C6/36 cell line for the correct design and cloning of CRISPR RNA into an insect plasmid vector to create a single guide RNA for the PHB2 gene target. After transfection of this plasmid vector into the C6/36 cells, cell clones selected by puromycin and/or limiting dilution were analyzed for insertions and deletions (INDELs) using PCR, sequencing and computational sequence decomposition. From this, we have identified mono-allelic and bi-allelic knockout cell clones. Using a mono-allelic knockout cell clone as an example, we characterized its INDELs by molecular cloning and computational analysis. Importantly, mono-allelic knockout was sufficient to reduce >80 % of PHB2 expression, which led to phenotypic switching and the propensity to form foci but was insufficient to affect growth rate or to inhibit Zika virus infection.•We provide a toolkit for CRISPR-Cas9-editing of the virus vector, Aedes albopictus C6/36 cell line•We validate this using a difficult-to-knockdown gene prohibitin 2•This toolkit is pivotal in understanding vector biology, vector competence, host-pathogen interactions and in gene annotations.

Genome resequencing and genome-wide polymorphisms in mosquito vectors Aedes aegypti and Aedes albopictus from south India

Aedes aegypti and Aedes albopictus mosquitoes spread major vector-borne viral diseases in tropical and sub-tropical regions of the globe. In this study, we sequenced the genome of Indian Ae. aegypti and Ae. albopictus and mapped to their reference genomes. Comparative genomics were performed between our strain and the reference strains. A total of 14,416,484 single nucleotide polymorphisms (SNPs) and 156,487 insertions and deletions (InDels) were found in Ae. aegypti, and 28,940,433 SNPs and 188,987 InDels in Ae. albopictus. Particular emphasis was given to gene families involved in mosquito digestion, development, and innate immunity, which could be putative candidates for vector control. Serine protease cascades and their inhibitors called serpins, play a central role in these processes. We extracted high-impact variants in genes associated with serine proteases and serpins. This study reports for the first time a high coverage genome sequence data of an Indian Ae. albopictus mosquito. The results from this study will provide insights into Indian Aedes specific polymorphisms and the evolution of immune related genes in mosquitoes, which can serve as a resource for future comparative genomics and those pursuing the development of targeted biopesticides for effective mosquito control strategies.

Transposon and Transgene Tribulations in Mosquitoes: A Perspective of piRNA Proportions

Mosquitoes, like Drosophila, are dipterans, the order of "true flies" characterized by a single set of two wings. Drosophila are prime model organisms for biomedical research, while mosquito researchers struggle to establish robust molecular biology in these that are arguably the most dangerous vectors of human pathogens. Both insects utilize the RNA interference (RNAi) pathway to generate small RNAs to silence transposons and viruses, yet details are emerging that several RNAi features are unique to each insect family, such as how culicine mosquitoes have evolved extreme genomic feature differences connected to their unique RNAi features. A major technical difference in the molecular genetic studies of these insects is that generating stable transgenic animals are routine in Drosophila but still variable in stability in mosquitoes, despite genomic DNA-editing advances. By comparing and contrasting the differences in the RNAi pathways of Drosophila and mosquitoes, in this review we propose a hypothesis that transgene DNAs are possibly more intensely targeted by mosquito RNAi pathways and chromatin regulatory pathways than in Drosophila. We review the latest findings on mosquito RNAi pathways, which are still much less well understood than in Drosophila, and we speculate that deeper study into how mosquitoes modulate transposons and viruses with Piwi-interacting RNAs (piRNAs) will yield clues to improving transgene DNA expression stability in transgenic mosquitoes.

Unveiling candidate genes for metabolic resistance to malathion in Aedes albopictus through RNA sequencing-based transcriptome profiling

Aedes albopictus, also known as the Asian tiger mosquito, is indigenous to the tropical forests of Southeast Asia. Ae. albopictus is expanding across the globe at alarming rates, raising concern over the transmission of mosquito-borne diseases, such as dengue, West Nile fever, yellow fever, and chikungunya fever. Since Ae. albopictus was reported in Houston (Harris County, Texas) in 1985, this species has rapidly expanded to at least 32 states across the United States. Public health efforts aimed at controlling Ae. albopictus, including surveillance and adulticide spraying operations, occur regularly in Harris County. Despite rotation of insecticides to mitigate the development of resistance, multiple mosquito species including Culex quinquefasciatus and Aedes aegypti in Harris County show organophosphate and pyrethroid resistance. Aedes albopictus shows relatively low resistance levels as compared to Ae. aegypti, but kdr-mutation and the expression of detoxification genes have been reported in Ae. albopictus populations elsewhere. To identify potential candidate detoxification genes contributing to metabolic resistance, we used RNA sequencing of field-collected malathion-resistant and malathion-susceptible, and laboratory-maintained susceptible colonies of Ae. albopictus by comparing the relative expression of transcripts from three major detoxification superfamilies involved in malathion resistance due to metabolic detoxification. Between these groups, we identified 12 candidate malathion resistance genes and among these, most genes correlated with metabolic detoxification of malathion, including four P450 and one alpha esterase. Our results reveal the metabolic detoxification and potential cuticular-based resistance mechanisms associated with malathion resistance in Ae. albopictus in Harris County, Texas.

A high heterozygosity genome assembly of Aedes albopictus enables the discovery of the association of PGANT3 with blood-feeding behavior

The Asian tiger mosquito, Aedes albopictus, is a global invasive species, notorious for its role in transmitting dangerous human arboviruses such as dengue and Chikungunya. Although hematophagous behavior is repulsive, it is an effective strategy for mosquitoes like Aedes albopictus to transmit viruses, posing a significant risk to human health. However, the fragmented nature of the Ae. albopictus genome assembly has been a significant challenge, hindering in-depth biological and genetic studies of this mosquito. In this research, we have harnessed a variety of technologies and implemented a novel strategy to create a significantly improved genome assembly for Ae. albopictus, designated as AealbF3. This assembly boasts a completeness rate of up to 98.1%, and the duplication rate has been minimized to 1.2%. Furthermore, the fragmented contigs or scaffolds of AealbF3 have been organized into three distinct chromosomes, an arrangement corroborated through syntenic plot analysis, which compared the genetic structure of Ae. albopictus with that of Ae. aegypti. Additionally, the study has revealed a phylogenetic relationship suggesting that the PGANT3 gene is implicated in the hematophagous behavior of Ae. albopictus. This involvement was preliminarily substantiated through RNA interference (RNAi) techniques and behavioral experiment. In summary, the AealbF3 genome assembly will facilitate new biological insights and intervention strategies for combating this formidable vector of disease. The innovative assembly process employed in this study could also serve as a valuable template for the assembly of genomes in other insects characterized by high levels of heterozygosity.

The chromosome-scale genome assembly for the West Nile vector Culex quinquefasciatus uncovers patterns of genome evolution in mosquitoes

BACKGROUND

Understanding genome organization and evolution is important for species involved in transmission of human diseases, such as mosquitoes. Anophelinae and Culicinae subfamilies of mosquitoes show striking differences in genome sizes, sex chromosome arrangements, behavior, and ability to transmit pathogens. However, the genomic basis of these differences is not fully understood.

METHODS

In this study, we used a combination of advanced genome technologies such as Oxford Nanopore Technology sequencing, Hi-C scaffolding, Bionano, and cytogenetic mapping to develop an improved chromosome-scale genome assembly for the West Nile vector Culex quinquefasciatus.

RESULTS

We then used this assembly to annotate odorant receptors, odorant binding proteins, and transposable elements. A genomic region containing male-specific sequences on chromosome 1 and a polymorphic inversion on chromosome 3 were identified in the Cx. quinquefasciatus genome. In addition, the genome of Cx. quinquefasciatus was compared with the genomes of other mosquitoes such as malaria vectors An. coluzzi and An. albimanus, and the vector of arboviruses Ae. aegypti. Our work confirms significant expansion of the two chemosensory gene families in Cx. quinquefasciatus, as well as a significant increase and relocation of the transposable elements in both Cx. quinquefasciatus and Ae. aegypti relative to the Anophelines. Phylogenetic analysis clarifies the divergence time between the mosquito species. Our study provides new insights into chromosomal evolution in mosquitoes and finds that the X chromosome of Anophelinae and the sex-determining chromosome 1 of Culicinae have a significantly higher rate of evolution than autosomes.

CONCLUSION

The improved Cx. quinquefasciatus genome assembly uncovered new details of mosquito genome evolution and has the potential to speed up the development of novel vector control strategies.

A tale of caution: How endogenous viral elements affect virus discovery in transcriptomic data

Large-scale metagenomic and -transcriptomic studies have revolutionized our understanding of viral diversity and abundance. In contrast, endogenous viral elements (EVEs), remnants of viral sequences integrated into host genomes, have received limited attention in the context of virus discovery, especially in RNA-Seq data. EVEs resemble their original viruses, a challenge that makes distinguishing between active infections and integrated remnants difficult, affecting virus classification and biases downstream analyses. Here, we systematically assess the effects of EVEs on a prototypical virus discovery pipeline, evaluate their impact on data integrity and classification accuracy, and provide some recommendations for better practices. We examined EVEs and exogenous viral sequences linked to Orthomyxoviridae, a diverse family of negative-sense segmented RNA viruses, in 13 genomic and 538 transcriptomic datasets of Culicinae mosquitoes. Our analysis revealed a substantial number of viral sequences in transcriptomic datasets. However, a significant portion appeared not to be exogenous viruses but transcripts derived from EVEs. Distinguishing between transcribed EVEs and exogenous virus sequences was especially difficult in samples with low viral abundance. For example, three transcribed EVEs showed full-length segments, devoid of frameshift and nonsense mutations, exhibiting sufficient mean read depths that qualify them as exogenous virus hits. Mapping reads on a host genome containing EVEs before assembly somewhat alleviated the EVE burden, but it led to a drastic reduction of viral hits and reduced quality of assemblies, especially in regions of the viral genome relatively similar to EVEs. Our study highlights that our knowledge of the genetic diversity of viruses can be altered by the underestimated presence of EVEs in transcriptomic datasets, leading to false positives and altered or missing sequence information. Thus, recognizing and addressing the influence of EVEs in virus discovery pipelines will be key in enhancing our ability to capture the full spectrum of viral diversity.

De novo genome assembly of the invasive mosquito species Aedes japonicus and Aedes koreicus

BACKGROUND

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.

METHODS

A de novo genome assembly of Ae. japonicus (Ajap1) and Ae. koreicus (Akor1) has been produced based on a hybrid approach that combines Oxford Nanopore long-read and Illumina short-read data. Their quality was ascertained using various metrics. Masking of repetitive elements, gene prediction and functional annotation was performed.

RESULTS

Sequence analysis revealed a very high presence of repetitive DNA and, among others, thermal adaptation genes and insecticide-resistance genes. Through the RNA-seq analysis of larvae and adults of Ae. koreicus and Ae. japonicus exposed to different temperatures, we also identified genes showing a differential temperature-dependent activation.

CONCLUSIONS

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.

Dense residential areas promote gene flow in dengue vector mosquito Aedes albopictus

Aedes albopictus is a successful disease vector due to its ability to survive in a wide range of habitats. Despite its ubiquity and impact on public health, little is known about its differential gene flow capabilities across different city habitats. We obtained a comprehensive dataset of >27,000 genome-wide DNA markers across 105 wild-caught Ae. albopictus individuals from Singapore, a dengue-endemic tropical city with heterogeneous landscapes from densely populated urban areas to forests. Despite Singapore's challenging small-scale heterogeneity, our landscape-genomic approach indicated that dense urban areas are characterized by higher Aedes gene flow rates than managed parks and forests. We documented the incidence of Wolbachia infections of Ae. albopictus involving two strains (wAlbA and wAlbB). Our results dispel the misconception that substantial dispersal of Ae. albopictus is limited to urban greenery, with wide implications for vector management and critical insights into urban planning strategies to combat dengue transmission.

RNAi-mediated CHS-2 silencing affects the synthesis of chitin and the formation of the peritrophic membrane in the midgut of Aedes albopictus larvae

BACKGROUND

Mosquitoes are an important vector of viral transmission, and due to the complexity of the pathogens they transmit, vector control may be the most effective strategy to control mosquito-borne diseases. Chitin is required for insect growth and development and is absent in higher animals and plants, so regulating the chitin synthesis pathway can serve as a potentially effective means to control vector insects. Most of the current research on the chitin synthase (CHS) gene is focused on chitin synthase-1 (CHS-1), while relatively little is known about chitin synthase-2 (CHS-2).

RESULTS

The CHS-2 gene of Ae. albopictus is highly conserved and closely related to that of Aedes aegypti. The expression of CHS-2 in the third-instar larvae and pupal stage of Ae. albopictus was relatively high, and CHS-2 expression in adult mosquitoes reached the highest value 24 h after blood-feeding. In the fourth-instar larvae of Ae. albopictus, CHS-2 expression was significantly higher in the midgut than in the epidermis. Silencing CHS-2 in Ae. albopictus larvae had no effect on larval survival and emergence. The expression of four genes related to chitin synthesis enzymes was significantly upregulated, the expression level of three genes was unchanged, and only the expression level of GFAT was significantly downregulated. The expression of chitin metabolism-related genes was also upregulated after silencing. The level of chitin in the midgut of Ae. albopictus larvae was significantly decreased, while the chitinase activity was unchanged. The epithelium of the midgut showed vacuolization, cell invagination and partial cell rupture, and the structure of the peritrophic membrane was destroyed or even absent.

METHODS

The expression of CHS-2 in different developmental stages and tissues of Aedes albopictus was detected by real-time fluorescence quantitative PCR (qPCR). After silencing CHS-2 of the fourth-instar larvae of Ae. albopictus by RNA interference (RNAi), the expression levels of genes related to chitin metabolism, chitin content and chitinase activity in the larvae were detected. The structure of peritrophic membrane in the midgut of the fourth-instar larvae after silencing was observed by paraffin section and hematoxylin-eosin (HE) staining.

CONCLUSION

CHS-2 can affect midgut chitin synthesis and breakdown by regulating chitin metabolic pathway-related genes and is involved in the formation of the midgut peritrophic membrane in Ae. albopictus, playing an important role in growth and development. It may be a potential target for enhancing other control methods.

Optimization of Dieldrin Selection for the Genetic Sexing of Aedes albopictus

The mass production of mosquitoes at an industrial scale requires efficient sex separation, which can be achieved through mechanical, genetic or artificial intelligence means. Compared with other methods, the genetic sexing approach offers the advantage of limiting costs and space by removing females at the larval stage. We recently developed a Genetic Sexing Strain (GSS) in Aedes albopictus based on the sex linkage of the rdl^R allele, conferring resistance to dieldrin, to the male (M) locus. It has been previously reported that dieldrin ingested by larvae can be detected in adults and bioaccumulated in predators, raising the question of its use at a large scale. In this context, we performed several experiments aiming at optimizing dieldrin selection by decreasing both dieldrin concentration and exposure time while maintaining a stable percentage of contaminating females averaging 1%. We showed that the previously used dieldrin exposure induced an important toxicity as it killed 60% of resistant males at the larval stage. We lowered this toxicity by reducing the dose and/or the exposure time to recover nearly all resistant males. We then quantified the residues of dieldrin in resistant male adults and showed that dieldrin toxicity in larvae was positively correlated with dieldrin concentrations detected in adults. Interestingly, we showed that the use of reduced dieldrin exposure led to a dieldrin quantification in adult males that was below the quantity threshold of the Gas Chromatography-Mass Spectrometry detection method. Presented data show that dieldrin exposure can be adjusted to suppress toxicity in males while achieving efficient sexing and lowering the levels of dieldrin residues in adults to barely quantifiable levels.

Midgut transcriptomic responses to dengue and chikungunya viruses in the vectors Aedes albopictus and Aedes malayensis

Dengue (DENV) and chikungunya (CHIKV) viruses are among the most preponderant arboviruses. Although primarily transmitted through the bite of Aedes aegypti mosquitoes, Aedes albopictus and Aedes malayensis are competent vectors and have an impact on arbovirus epidemiology. Here, to fill the gap in our understanding of the molecular interactions between secondary vectors and arboviruses, we used transcriptomics to profile the whole-genome responses of A. albopictus to CHIKV and of A. malayensis to CHIKV and DENV at 1 and 4 days post-infection (dpi) in midguts. In A. albopictus, 1793 and 339 genes were significantly regulated by CHIKV at 1 and 4 dpi, respectively. In A. malayensis, 943 and 222 genes upon CHIKV infection, and 74 and 69 genes upon DENV infection were significantly regulated at 1 and 4 dpi, respectively. We reported 81 genes that were consistently differentially regulated in all the CHIKV-infected conditions, identifying a CHIKV-induced signature. We identified expressed immune genes in both mosquito species, using a de novo assembled midgut transcriptome for A. malayensis, and described the immune architectures. We found the JNK pathway activated in all conditions, generalizing its antiviral function to Aedines. Our comprehensive study provides insight into arbovirus transmission by multiple Aedes vectors.

Climate change and vector-borne diseases: a multi-omics approach of temperature-induced changes in the mosquito

BACKGROUND

Climate change and globalization contribute to the expansion of mosquito vectors and their associated pathogens. Long spared, temperate regions have had to deal with the emergence of arboviruses traditionally confined to tropical regions. Chikungunya virus (CHIKV) was reported for the first time in Europe in 2007, causing a localized outbreak in Italy, which then recurred repeatedly over the years in other European localities. This raises the question of climate effects, particularly temperature, on the dynamics of vector-borne viruses. The objective of this study is to improve the understanding of the molecular mechanisms set up in the vector in response to temperature.

METHODS

We combine three complementary approaches by examining Aedes albopictus mosquito gene expression (transcriptomics), bacterial flora (metagenomics) and CHIKV evolutionary dynamics (genomics) induced by viral infection and temperature changes.

RESULTS

We show that temperature alters profoundly mosquito gene expression, bacterial microbiome and viral population diversity. We observe that (i) CHIKV infection upregulated most genes (mainly in immune and stress-related pathways) at 20°C but not at 28°C, (ii) CHIKV infection significantly increased the abundance of Enterobacteriaceae Serratia marcescens at 28°C and (iii) CHIKV evolutionary dynamics were different according to temperature.

CONCLUSION

The substantial changes detected in the vectorial system (the vector and its bacterial microbiota, and the arbovirus) lead to temperature-specific adjustments to reach the ultimate goal of arbovirus transmission; at 20°C and 28°C, the Asian tiger mosquito Ae. albopictus was able to transmit CHIKV at the same efficiency. Therefore, CHIKV is likely to continue its expansion in the northern regions and could become a public health problem in more countries than those already affected in Europe.

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.

Ancient origin of Jingchuvirales derived glycoproteins integrated in arthropod genomes

Endogenous virus elements (EVEs) are viral-derived sequences integrated into their host genomes. EVEs of the Jingchuvirales order were detected in a wide range of insect genomes covering several distantly related families. Moreover, Jingchuvirales-derived glycoproteins were recently associated by our group with the origin of a putative new retrovirus based on a glycoprotein captured by a mosquito retrotransposon. But, except for mosquitoes, there is a lack of a more detailed understanding of the endogenization mechanism, timing, and frequency per Jingchuvirales viral lineages. Here we screened Jingchuvirales glycoprotein-derived EVEs (Jg-EVEs) in eukaryotic genomes. We found six distinct endogenization events of Jg-EVEs, that belong to two out of five known Jingchuvirales families (Chuviridae and Natareviridae). For seven arthropod families bearing Jg-EVEs there is no register of bona fide circulating chuvirus infection. Hence, our results show that Jingchuvirales viruses infected or still infect these host families. Although we found abundant evidence of LTR-Gypsy retrotransposons fragments associated with the glycoprotein in Hymenoptera and other insect orders, our results show that the widespread distribution of Jingchuvirales glycoproteins in extant Arhtropods is a result of multiple ancient endogenization events and that these virus fossils are being vertically inherited in Arthropods genomes for millions of years.

Molecular identification of the chitinase genes in Aedes albopictus and essential roles of AaCht10 in pupal-adult transition

BACKGROUND

Aedes albopictus is an increasingly serious threat in public health due to it is vector of multiple arboviruses that cause devastating human diseases, as well as its widening distribution in recent years. Insecticide resistance is a serious problem worldwide that limits the efficacy of chemical control strategies against Ae. albopictus. Chitinase genes have been widely recognized as attractive targets for the development of effective and environmentally safe insect management measures.

METHODS

Chitinase genes of Ae. albopictus were identified and characterized on the basis of bioinformatics search of the referenced genome. Gene characterizations and phylogenetic relationships of chitinase genes were investigated, and spatio-temporal expression pattern of each chitinase gene was evaluated using qRT-PCR. RNA interference (RNAi) was used to suppress the expression of AaCht10, and the roles of AaCht10 were verified based on phynotype observations, chitin content analysis and hematoxylin and eosin (H&E) stain of epidermis and midgut.

RESULTS

Altogether, 14 chitinase-related genes (12 chitinase genes and 2 IDGFs) encoding 17 proteins were identified. Phylogenetic analysis showed that all these AaChts were classified into seven groups, and most of them were gathered into group IX. Only AaCht5-1, AaCht10 and AaCht18 contained both catalytic and chitin-binding domains. Different AaChts displayed development- and tissue-specific expression profiling. Suppression of the expression of AaCht10 resulted in abnormal molting, increased mortality, decreased chitin content and thinning epicuticle, procuticle and midgut wall of pupa.

CONCLUSIONS

Findings of the present study will aid in determining the biological functions of AaChts and also contribute to using AaChts as potential target for mosquito management.

Chromatin profiling identifies transcriptional readthrough as a conserved mechanism for piRNA biogenesis in mosquitoes

The piRNA pathway in mosquitoes differs substantially from other model organisms, with an expanded PIWI gene family and functions in antiviral defense. Here, we define core piRNA clusters as genomic loci that show ubiquitous piRNA expression in both somatic and germline tissues. These core piRNA clusters are enriched for non-retroviral endogenous viral elements (nrEVEs) in antisense orientation and depend on key biogenesis factors, Veneno, Tejas, Yb, and Shutdown. Combined transcriptome and chromatin state analyses identify transcriptional readthrough as a conserved mechanism for cluster-derived piRNA biogenesis in the vector mosquitoes Aedes aegypti, Aedes albopictus, Culex quinquefasciatus, and Anopheles gambiae. Comparative analyses between the two Aedes species suggest that piRNA clusters function as traps for nrEVEs, allowing adaptation to environmental challenges such as virus infection. Our systematic transcriptome and chromatin state analyses lay the foundation for studies of gene regulation, genome evolution, and piRNA function in these important vector species.

What Are the Functional Roles of Piwi Proteins and piRNAs in Insects?

Research on Piwi proteins and piRNAs in insects has focused on three experimental models: oogenesis and spermatogenesis in Drosophila melanogaster, the antiviral response in Aedes mosquitoes and the molecular analysis of primary and secondary piRNA biogenesis in Bombyx mori-derived BmN4 cells. Significant unique and complementary information has been acquired and has led to a greater appreciation of the complexity of piRNA biogenesis and Piwi protein function. Studies performed in other insect species are emerging and promise to add to the current state of the art on the roles of piRNAs and Piwi proteins. Although the primary role of the piRNA pathway is genome defense against transposons, particularly in the germline, recent findings also indicate an expansion of its functions. In this review, an extensive overview is presented of the knowledge of the piRNA pathway that so far has accumulated in insects. Following a presentation of the three major models, data from other insects were also discussed. Finally, the mechanisms for the expansion of the function of the piRNA pathway from transposon control to gene regulation were considered.

Whole-transcriptome profiling across different developmental stages of Aedes albopictus (Diptera: Culicidae) provides insights into chitin-related non-coding RNA and competing endogenous RNA networks

BACKGROUND

The Asian tiger mosquito, Aedes albopictus, is one of the most invasive species and a vector of numerous arboviruses. The deleterious effects of long-term and inappropriate use of chemical pesticides have stimulated the exploration of new, environmentally friendly control strategies. Non-coding RNAs (ncRNAs) have been proven to participate in almost all biological processes of insects.

METHODS

In this study, circular RNAs (circRNAs) and microRNAs (miRNAs) covering five developmental stages [egg, early larvae, late larvae, pupae, adult (female and male)] of A. albopictus were obtained using whole-transcriptome sequencing technology. Combined with long non-coding RNAs (lncRNAs) from previous research, circRNA/lncRNA‒miRNA‒mitochondrial RNA (mRNA) networks were constructed.

RESULTS

A total of 1434 circRNAs and 208 miRNAs were identified. More differentially expressed circRNAs (DE circRNAs) and miRNAs (DE miRNAs) were found in the egg versus early larvae comparison group. Functional enrichment analysis demonstrated that most of the circRNA/lncRNA‒miRNA‒mRNA networks were involved in chitin metabolism. Hub genes of each circRNA/lncRNA‒miRNA‒mRNA network were screened out, which can be used as novel targets to disturb the molting process of A. albopictus.

CONCLUSIONS

Regulatory relationships obtained from competing endogenous RNA (ceRNA) networks provide more information to manipulate the metamorphosis process and are helpful for developing effective and sustainable methods to control mosquitoes.

Visualization of the Linear and Spatial Organization of Chromosomes in Mosquitoes

Mosquitoes are vectors of dangerous human diseases such as malaria, dengue, Zika, West Nile fever, and lymphatic filariasis. Visualization of the linear and spatial organization of mosquito chromosomes is important for understanding genome structure and function. Utilization of chromosomal inversions as markers for population genetics studies yields insights into mosquito adaptation and evolution. Cytogenetic approaches assist with the development of chromosome-scale genome assemblies that are useful tools for studying mosquito biology and for designing novel vector control strategies. Fluorescence in situ hybridization is a powerful technique for localizing specific DNA sequences within the linear chromosome structure and within the spatial organization of the cell nucleus. Here, we introduce protocols used in our laboratories for chromosome visualization and their application in mosquitoes.

Population genetics of an invasive mosquito vector, Aedes albopictus in the Northeastern USA

The Asian tiger mosquito (Aedes albopictus) arrived in the USA in the 1980’s and rapidly spread throughout eastern USA within a decade. The predicted northern edge of its overwintering distribution on the East Coast of the USA roughly falls across New York, Connecticut, and Massachusetts, where the species has been recorded as early as 2000. It is unclear whether Ae. albopictus populations have become established and survive the cold winters in these areas or are recolonized every year. We genotyped and analyzed populations of Ae. albopictus from the northeast USA using 15 microsatellite markers and compared them with other populations across the country and to representatives of the major global genetic clades to investigate their connectivity and stability. Founder effects or bottlenecks were rare at the northern range of the Ae. albopictus distribution in the northeastern USA, with populations displaying high levels of genetic diversity and connectivity along the East Coast. There is no evidence of population turnover in Connecticut during the course of three consecutive years, with consistent genetic structure throughout this period. Overall, these results support the presence of established populations of Ae. albopictus in New York, Connecticut, and Massachusetts, successfully overwintering and migrating in large numbers. Given the stability and interconnectedness of these populations, Ae. albopictus has the potential to continue to proliferate and expand its range northward under mean warming conditions of climate change. Efforts to control Ae. albopictus in these areas should thus focus on vector suppression rather than eradication strategies, as local populations have become firmly established and are expected to reemerge every summer.

Virome and nrEVEome diversity of Aedes albopictus mosquitoes from La Reunion Island and China

Background

Aedes albopictus is a public health threat for its worldwide spread and ability to transmit arboviruses. Understanding mechanisms of mosquito immunity can provide new tools to control arbovirus spread. The genomes of Aedes mosquitoes contain hundreds of nonretroviral endogenous viral elements (nrEVEs), which are enriched in piRNA clusters and produce piRNAs, with the potential to target cognate viruses. Recently, one nrEVE was shown to limit cognate viral infection through nrEVE-derived piRNAs. These findings suggest that nrEVEs constitute an archive of past viral infection and that the landscape of viral integrations may be variable across populations depending on their viral exposure.

Methods

We used bioinformatics and molecular approaches to identify known and novel (i.e. absent in the reference genome) viral integrations in the genome of wild collected Aedes albopictus mosquitoes and characterize their virome.

Results

We showed that the landscape of viral integrations is dynamic with seven novel viral integrations being characterized, but does not correlate with the virome, which includes both viral species known and unknown to infect mosquitoes. However, the small RNA coverage profile of nrEVEs and the viral genomic contigs we identified confirmed an interaction among these elements and the piRNA and siRNA pathways in mosquitoes.

Conclusions

Mosquitoes nrEVEs have been recently described as a new form of heritable, sequence-specific mechanism of antiviral immunity. Our results contribute to understanding the dynamic distribution of nrEVEs in the genomes of wild Ae. albopictus and their interaction with mosquito viruses.

Genome-wide SNPs reveal novel patterns of spatial genetic structure in Aedes albopictus (Diptera Culicidae) population in China

Introduction

Since the second half of the 20th century, Aedes albopictus, a vector for more than 20 arboviruses, has spread worldwide. Aedes albopictus is the main vector of infectious diseases transmitted by Aedes mosquitoes in China, and it has caused concerns regarding public health. A comprehensive understanding of the spatial genetic structure of this vector species at a genomic level is essential for effective vector control and the prevention of vector-borne diseases.

Methods

During 2016-2018, adult female Ae. albopictus mosquitoes were collected from eight different geographical locations across China. Restriction site-associated DNA sequencing (RAD-seq) was used for high-throughput identification of single nucleotide polymorphisms (SNPs) and genotyping of the Ae. albopictus population. The spatial genetic structure was analyzed and compared to those exhibited by mitochondrial cytochrome c oxidase subunit 1 (cox1) and microsatellites in the Ae. albopictus population.

Results

A total of 9,103 genome-wide SNP loci in 101 specimens and 32 haplotypes of cox1 in 231 specimens were identified in the samples from eight locations in China. Principal component analysis revealed that samples from Lingshui and Zhanjiang were more genetically different than those from the other locations. The SNPs provided a better resolution and stronger signals for novel spatial population genetic structures than those from the cox1 data and a set of previously genotyped microsatellites. The fixation indexes from the SNP dataset showed shallow but significant genetic differentiation in the population. The Mantel test indicated a positive correlation between genetic distance and geographical distance. However, the asymmetric gene flow was detected among the populations, and it was higher from south to north and west to east than in the opposite directions.

Conclusions

The genome-wide SNPs revealed seven gene pools and fine spatial genetic structure of the Ae. albopictus population in China. The RAD-seq approach has great potential to increase our understanding of the spatial dynamics of population spread and establishment, which will help us to design new strategies for controlling vectors and mosquito-borne diseases.

Impact of CHIKV Replication on the Global Proteome of Aedes albopictus Cells

Arboviruses are some of the important causative agents of mosquito-mediated viral diseases. These viruses are transmitted between vector and host during the blood meal. Upon viral entry, host replication machinery is hijacked, supporting new virus particle production and thereby allowing viral survival in the host. In this process, host proteins interact with viral proteins to either facilitate viral replication, or they may provide antiviral defense mechanisms. In this study, we analyzed the impact of chikungunya virus (CHIKV) infection on the global proteome of Dicer active Aedes albopictus cells during the early and late time points of infection. We utilized a bottom-up approach of global proteomics analysis, and we used label-free quantitative mass spectrometry to identify the global protein signatures of Ae. albopictus at two different time points upon CHIKV infection. The mass spectrometry data analysis of the early time point revealed that proteins belonging to pathways such as translation, RNA processing, and cellular metabolic processes were less in abundance, whereas those belonging to pathways such as cellular catabolic process and organic substance transport were significantly abundant. At later time points, proteins belonging to pathways such as cellular metabolic processes, primary metabolic process, organonitrogen compound metabolic process, and organic substance metabolic process were found to be decreased in their presence, whereas those belonging to pathways such as RNA processing, gene expression, macromolecule metabolic processing, and nitrogen compound metabolic processing were found to be abundant during CHIKV infection, indicating that modulation in gene expression favoring cell survival occurs at a later time point, suggesting a survival strategy of Aedes cells to counter prolonged CHIKV infection.

Revised Annotation and Characterization of Novel Aedes albopictus miRNAs and Their Potential Functions in Dengue Virus Infection

Simple Summary

Aedes albopictus (Ae. albopictus) is an important vector of the dengue virus. Genetics and molecular studies of virus infection in mosquito vectors are important to uncover the basic biology of the virus. It has been reported that miRNAs are important and possess functional roles in virus infection in Ae. albopictus. Here, we report a comprehensive catalog of miRNAs using the latest genome version of Ae. albopictus. We discovered a total of 72 novel mature miRNAs, 44 of which were differentially expressed in C6/36 cells infected with the dengue virus. Target prediction analysis revealed that the differentially expressed miRNAs were involved in lipid metabolism and protein processing in the endoplasmic reticulum. Results from this study provide a valuable resource for researchers to study miRNAs in this mosquito vector, especially in host–virus interactions.

Abstract

The Asian tiger mosquito, Ae. albopictus, is a highly invasive species that transmits several arboviruses including dengue (DENV), Zika (ZIKV), and chikungunya (CHIKV). Although several studies have identified microRNAs (miRNAs) in Ae. albopictus, it is crucial to extend and improve current annotations with both the newly improved genome assembly and the increased number of small RNA-sequencing data. We combined our high-depth sequence data and 26 public datasets to re-annotate Ae. albopictus miRNAs and found a total of 72 novel mature miRNAs. We discovered that the expression of novel miRNAs was lower than known miRNAs. Furthermore, compared to known miRNAs, novel miRNAs are prone to expression in a stage-specific manner. Upon DENV infection, a total of 44 novel miRNAs were differentially expressed, and target prediction analysis revealed that miRNA-target genes were involved in lipid metabolism and protein processing in endoplasmic reticulum. Taken together, the miRNA annotation profile provided here is the most comprehensive to date. We believed that this would facilitate future research in understanding virus–host interactions, particularly in the role of miRNAs.

Selective targeting of biting females to control mosquito-borne infectious diseases

Mosquitoes are vectors for a number of infectious diseases. Only females feed on blood to provision for their embryos and, in doing so, transmit pathogens to the associated vertebrate hosts. Therefore, sex is an important phenotype in the context of genetic control programs, both for sex separation in the rearing facilities to avoid releasing biting females and for ways to distort the sex ratio towards nonbiting males. We review recent progress in the fundamental knowledge of sex determination and sex chromosomes in mosquitoes and discuss new methods to achieve sex separation and sex ratio distortion to help control mosquito-borne infectious diseases. We conclude by suggesting a few critical areas for future research.

The worldwide spread of Aedes albopictus: New insights from mitogenomes

The tiger mosquito (Aedes albopictus) is one of the most invasive species in the world and a competent vector for numerous arboviruses, thus the study and monitoring of its fast worldwide spread is crucial for global public health. The small extra-nuclear and maternally-inherited mitochondrial DNA represents a key tool for reconstructing phylogenetic and phylogeographic relationships within a species, especially when analyzed at the mitogenome level. Here the mitogenome variation of 76 tiger mosquitoes, 37 of which new and collected from both wild adventive populations and laboratory strains, was investigated. This analysis significantly improved the global mtDNA phylogeny of Ae. albopictus, uncovering new branches and sub-branches within haplogroup A1, the one involved in its recent worldwide spread. Our phylogeographic approach shows that the current distribution of tiger mosquito mitogenome variation has been strongly affected by clonal and sub-clonal founder events, sometimes involving wide geographic areas, even across continents, thus shedding light on the Asian sources of worldwide adventive populations. In particular, different starting points for the two major clades within A1 are suggested, with A1a spreading mainly along temperate areas from Japanese and Chinese sources, and A1b arising and mainly diffusing in tropical areas from a South Asian source.

Multiple introductions and overwintering shape the progressive invasion of Aedes albopictus beyond the Alps

Aedes albopictus originates from Southeast Asia and is considered one of the most invasive species globally. This mosquito is a nuisance and a disease vector of significant public health relevance. In Europe, Ae. albopictus is firmly established and widespread south of the Alps, a mountain range that forms a formidable biogeographic barrier to many organisms. Recent reports of Ae. albopictus north of the Alps raise questions of (1) the origins of its recent invasion, and (2) if this mosquito has established overwintering populations north of the Alps. To answer these questions, we analyzed population genomic data from >4000 genome-wide SNPs obtained through double-digest restriction site-associated DNA sequencing. We collected SNP data from specimens from six sites in Switzerland, north and south of the Alps, and analyzed them together with specimens from other 33 European sites, five from the Americas, and five from its Asian native range. At a global level, we detected four genetic clusters with specimens from Indonesia, Brazil, and Japan as the most differentiated, whereas specimens from Europe, Hong Kong, and USA largely overlapped. Across the Alps, we detected a weak genetic structure and high levels of genetic admixture, supporting a scenario of rapid and human-aided dispersal along transportation routes. While the genetic pattern suggests frequent re-introductions into Switzerland from Italian sources, the recovery of a pair of full siblings in two consecutive years in Strasbourg, France, suggests the presence of an overwintering population north of the Alps. The suggestion of overwintering populations of Ae. albopictus north of the Alps and the expansion patterns identified points to an increased risk of further northward expansion and the need for increased surveillance of mosquito populations in Northern Europe.

Exploring new genomic territories with emerging model insects

Improvements in reference genome generation for insects and across the tree of life are extending the concept and utility of model organisms beyond traditional laboratory-tractable supermodels. Species or groups of species with comprehensive genome resources can be developed into model systems for studying a large variety of biological phenomena. Advances in sequencing and assembly technologies are supporting these emerging genome-enabled model systems by producing resources that are increasingly accurate and complete. Nevertheless, quality controls including assessing gene content completeness are required to ensure that these data can be included in expanding catalogues of high-quality references that will greatly advance understanding of insect biology and evolution.

The AalNix3&4 isoform is required and sufficient to convert Aedes albopictus females into males

Aedes albopictus is one of the most invasive insect species in the world and an effective vector for many important arboviruses. We reported previously that Ae. albopictus Nix (AalNix) is the male-determining factor of this species. However, whether AalNix alone is sufficient to initiate male development is unknown. Transgenic lines that express each of the three AalNix isoforms from the native promoter were obtained using piggyBac transformation. We verified the stable expression of AalNix isoforms in the transgenic lines and confirm that one isoform, AalNix3&4, is sufficient to convert females into fertile males (pseudo-males) that are indistinguishable from wild-type males. We also established a stable sex-converted female mosquito strain, AalNix3&4-♂4-pseudo-male. The pseudo-male mosquitoes can fly and mate normally with wild-type female, although their mating competitiveness is lower than wild-type. This work further clarifies the role of AalNix in the sex determination pathway and will facilitate the development of Ae. albopictus control strategies that rely on male-only releases such as SIT and sex-ratio distortion.

Maternally Instigated Diapause in Aedes albopictus: Coordinating Experience and Internal State for Survival in Variable Environments

The Asian tiger mosquito, Aedes albopictus, is one of the most dangerous invasive species in the world. Females bite mammalian hosts, including humans, to obtain blood for egg development. The ancestral range of Ae. albopictus likely spanned from India to Japan and this species has since invaded a substantial portion of the globe. Ae. albopictus can be broadly categorized into temperate and tropical populations. One key to their ability to invade diverse ecological spaces is the capacity of females to detect seasonal changes and produce stress-resistant eggs that survive harsh winters. Females living in temperate regions respond to cues that predict the onset of unfavorable environmental conditions by producing eggs that enter maternally instigated embryonic diapause, a developmentally arrested state, which allows species survival by protecting the embryos until favorable conditions return. To appropriately produce diapause eggs, the female must integrate environmental cues and internal physiological state (blood feeding and reproductive status) to allocate nutrients and regulate reproduction. There is variation in reproductive responses to environmental cues between interfertile tropical and temperate populations depending on whether females are actively producing diapause vs. non-diapause eggs and whether they originate from populations that are capable of diapause. Although diapause-inducing environmental cues and diapause eggs have been extensively characterized, little is known about how the female detects gradual environmental changes and coordinates her reproductive status with seasonal dynamics to lay diapause eggs in order to maximize offspring survival. Previous studies suggest that the circadian system is involved in detecting daylength as a critical cue. However, it is unknown which clock network components are important, how these connect to reproductive physiology, and how they may differ between behavioral states or across populations with variable diapause competence. In this review, we showcase Ae. albopictus as an emerging species for neurogenetics to study how the nervous system combines environmental conditions and internal state to optimize reproductive behavior. We review environmental cues for diapause induction, downstream pathways that control female metabolic changes and reproductive capacity, as well as diapause heterogeneity between populations with different evolutionary histories. We highlight genetic tools that can be implemented in Ae. albopictus to identify signaling molecules and cellular circuits that control diapause. The tools and discoveries made in this species could translate to a broader understanding of how environmental cues are interpreted to alter reproductive physiology in other species and how populations with similar genetic and circuit organizations diversify behavioral patterns. These approaches may yield new targets to interfere with mosquito reproductive capacity, which could be exploited to reduce mosquito populations and the burden of the pathogens they transmit.

Systematic identification and characterization of long noncoding RNAs (lncRNAs) during Aedes albopictus development

Background

Aedes albopictus originated in the tropical forests of Southeast Asia and can currently be found on all continents. As one of the main arboviral vectors, the control of Ae. albopictus requires novel strategies, informed by a deep knowledge of its biology. Little is known regarding mosquito long noncoding RNAs (lncRNAs), which are transcripts longer than 200 nucleotides that lack protein-coding potential and have roles in developmental regulation.

Results

Based on RNA-seq data from five developmental time points, eggs, early larvae, late larvae, pupae, and adults (female and male) of Ae. albopictus, 21,414 lncRNAs were characterized in this study. Differential expression analysis revealed that lncRNAs exhibited developmental stage specificity. The expression of most lncRNAs was upregulated at the onset of metamorphosis developmental stages. More differentially expressed lncRNAs were observed between eggs and early larvae. Weighted gene co-expression network analysis (WGCNA) further confirmed that the expression patterns of lncRNAs were obviously correlated with specific developmental time points. Functional annotation using co-expression analysis revealed that lncRNAs may be involved in the regulation of metamorphic developmental transitions of Ae. albopictus. The hub lncRNAs and hub gene clusters were identified for each module that were highly associated with specific developmental time points.

Conclusions

The results of this study will facilitate future researches to elucidate the regulatory mechanisms of lncRNAs in the development of Ae. albopictus and utilize lncRNAs to assist with mosquito control.

Transgenic expression of Nix converts genetic females into males and allows automated sex sorting in Aedes albopictus

Aedes albopictus is a major vector of arboviruses. Better understanding of its sex determination is crucial for developing mosquito control tools, especially genetic sexing strains. In Aedes aegypti, Nix is the primary gene responsible for masculinization and Nix-expressing genetic females develop into fertile, albeit flightless, males. In Ae. albopictus, Nix has also been implicated in masculinization but its role remains to be further characterized. In this work, we establish Ae. albopictus transgenic lines ectopically expressing Nix. Several are composed exclusively of genetic females, with transgenic individuals being phenotypic and functional males due to the expression of the Nix transgene. Their reproductive fitness is marginally impaired, while their flight performance is similar to controls. Overall, our results show that Nix is sufficient for full masculinization in Ae. albopictus. Moreover, the transgene construct contains a fluorescence marker allowing efficient automated sex sorting. Consequently, such strains constitute valuable sexing strains for genetic control.

Nix expression with a fluorescent marker in genetically female Ae. albopictus causes masculinization with minimal effects to reproductive fitness and flight performance.

The Peptide Hormone CNMa Influences Egg Production in the Mosquito Aedes aegypti

Mosquito reproduction is regulated by a suite of hormones, many acting through membrane-bound receptor proteins. The Aedes aegypti G protein-coupled receptors AAEL024199 (AeCNMaR-1a) and AAEL018316 (AeCNMaR-1b) were identified as orthologs of the Drosophila melanogaster CNMa receptor (DmCNMaR). The receptor was duplicated early in the evolution of insects, and subsequently in Culicidae, into what we refer to as CNMaR-1a and CNMaR-1b. AeCNMaR-1a is only detected in male mosquito antennae while AeCNMaR-1b is expressed at high levels in mosquito ovaries. Using a heterologous cell assay, we determined that AeCNMa activates AeCNMaR-1a with a ~10-fold lower concentration than it does AeCNMaR-1b, though both receptors displayed half maximal effective concentrations of AeCNMa in the low nanomolar range. Finally, we show that injections of AeCNMa into blood-fed mated female Ae. aegypti resulted in fewer eggs laid.

The diversity of endogenous viral elements in insects

We provide an overview of the currently known diversity of viral sequences integrated into insect genomes. Such endogenous viral elements (EVE) have so far been annotated in at least eight insect orders and can be assigned to at least three families of large double-stranded (ds) DNA viruses, at least 22 families of RNA viruses, and three families of single-stranded DNA viruses. The study of these EVE has already produced important insights into insect-virus interactions, including the discovery of a new form of adaptive antiviral immunity. Insect EVE diversity will continue to increase as new insect genomes and exogenous viruses are sequenced, which will continue to make paleovirology a vibrant research field in this group of animals in the years to come.

Endogenous viral elements in mosquito genomes: current knowledge and outstanding questions

Integrations from non-retroviral RNA viruses (nrEVEs) have been identified across several taxa, including mosquitoes. Amongst all Culicinae species, the viral vectors Aedes aegypti and Aedes albopictus stand out for their high number of nrEVEs. In addition, Aedes nrEVEs are enriched in piRNA clusters and generate piRNAs that can silence incoming viral genomes. As such, nrEVEs represent a new form of inherited antiviral immunity. To propel this discovery into novel transmission-blocking vector control strategies, a deeper understanding of nrEVE biology and evolution is essential because differences in the landscape of nrEVEs have been identified in wild-caught mosquitoes, the piRNA profile of nrEVEs is not homogeneous and nrEVEs outside piRNA clusters exist and are expressed at the mRNA level. Here we summarise current knowledge on nrEVEs in mosquitoes and we point out the many unanswered questions and potentials of these genomic elements.

Gosling: A Grammar-based Toolkit for Scalable and Interactive Genomics Data Visualization

The combination of diverse data types and analysis tasks in genomics has resulted in the development of a wide range of visualization techniques and tools. However, most existing tools are tailored to a specific problem or data type and offer limited customization, making it challenging to optimize visualizations for new analysis tasks or datasets. To address this challenge, we designed Gosling-a grammar for interactive and scalable genomics data visualization. Gosling balances expressiveness for comprehensive multi-scale genomics data visualizations with accessibility for domain scientists. Our accompanying JavaScript toolkit called Gosling.js provides scalable and interactive rendering. Gosling.js is built on top of an existing platform for web-based genomics data visualization to further simplify the visualization of common genomics data formats. We demonstrate the expressiveness of the grammar through a variety of real-world examples. Furthermore, we show how Gosling supports the design of novel genomics visualizations. An online editor and examples of Gosling.js, its source code, and documentation are available at https://gosling.js.org.

A Review of Discovery Profiling of PIWI-Interacting RNAs and Their Diverse Functions in Metazoans

PIWI-interacting RNAs (piRNAs) are a class of small non-coding RNAs (sncRNAs) that perform crucial biological functions in metazoans and defend against transposable elements (TEs) in germ lines. Recently, ubiquitously expressed piRNAs were discovered in soma and germ lines using small RNA sequencing (sRNA-seq) in humans and animals, providing new insights into the diverse functions of piRNAs. However, the role of piRNAs has not yet been fully elucidated, and sRNA-seq studies continue to reveal different piRNA activities in the genome. In this review, we summarize a set of simplified processes for piRNA analysis in order to provide a useful guide for researchers to perform piRNA research suitable for their study objectives. These processes can help expand the functional research on piRNAs from previously reported sRNA-seq results in metazoans. Ubiquitously expressed piRNAs have been discovered in the soma and germ lines in Annelida, Cnidaria, Echinodermata, Crustacea, Arthropoda, and Mollusca, but they are limited to germ lines in Chordata. The roles of piRNAs in TE silencing, gene expression regulation, epigenetic regulation, embryonic development, immune response, and associated diseases will continue to be discovered via sRNA-seq.

Uncovering the Worldwide Diversity and Evolution of the Virome of the Mosquitoes Aedes aegypti and Aedes albopictus

Aedes aegypti, the yellow fever mosquito, and Aedes albopictus, the Asian tiger mosquito, are the most significant vectors of dengue, Zika, and Chikungunya viruses globally. Studies examining host factors that control arbovirus transmission demonstrate that insect-specific viruses (ISVs) can modulate mosquitoes' susceptibility to arbovirus infection in both in vivo and in vitro co-infection models. While research is ongoing to implicate individual ISVs as proviral or antiviral factors, we have a limited understanding of the composition and diversity of the Aedes virome. To address this gap, we used a meta-analysis approach to uncover virome diversity by analysing ~3000 available RNA sequencing libraries representing a worldwide geographic range for both mosquitoes. We identified ten novel viruses and previously characterised viruses, including mononegaviruses, orthomyxoviruses, negeviruses, and a novel bi-segmented negev-like group. Phylogenetic analysis suggests close relatedness to mosquito viruses implying likely insect host range except for one arbovirus, the multi-segmented Jingmen tick virus (Flaviviridae) in an Italian colony of Ae. albopictus. Individual mosquito transcriptomes revealed remarkable inter-host variation of ISVs within individuals from the same colony and heterogeneity between different laboratory strains. Additionally, we identified striking virus diversity in Wolbachia infected Aedes cell lines. This study expands our understanding of the virome of these important vectors. It provides a resource for further assessing the ecology, evolution, and interaction of ISVs with their mosquito hosts and the arboviruses they transmit.

Competition between Usutu virus and West Nile virus during simultaneous and sequential infection of Culex pipiens mosquitoes

Usutu virus (USUV) and West Nile virus (WNV) are closely related mosquito-borne flaviviruses that are mainly transmitted between bird hosts by vector mosquitoes. Infections in humans are incidental but can cause severe disease. USUV is endemic in large parts of Europe, while WNV mainly circulates in Southern Europe. In recent years, WNV is also frequently detected in Northern Europe, thereby expanding the area where both viruses co-circulate. However, it remains unclear how USUV may affect the future spread of WNV and the likelihood of human co-infection. Here we investigated whether co-infections with both viruses in cell lines and their primary mosquito vector, Culex pipiens, affect virus replication and transmission dynamics. We show that USUV is outcompeted by WNV in mammalian, avian and mosquito cells during co-infection. Mosquitoes that were exposed to both viruses simultaneously via infectious blood meal displayed significantly reduced USUV transmission compared to mosquitoes that were only exposed to USUV (from 15% to 3%), while the infection and transmission of WNV was unaffected. In contrast, when mosquitoes were pre-infected with USUV via infectious blood meal, WNV transmission was significantly reduced (from 44% to 17%). Injection experiments established the involvement of the midgut in the observed USUV-mediated WNV inhibition. The competition between USUV and WNV during co-infection clearly indicates that the chance of concurrent USUV and WNV transmission via a single mosquito bite is low. The competitive relation between USUV and WNV may impact virus transmission dynamics in the field and affect the epidemiology of WNV in Europe.

Flavivirus integrations in Aedes aegypti are limited and highly conserved across samples from different geographic regions unlike integrations in Aedes albopictus

Mosquitoes of the genus Aedes are the main vectors of many viruses, e.g. dengue and Zika, which affect millions of people each year and for which there are limited treatment options. Understanding how Aedes mosquitoes tolerate high viral loads may lead to better disease control strategies. Elucidating endogenous viral elements (EVEs) within vector genomes may give exploitable biological insights. Previous studies have reported the presence of a large number of EVEs in Aedes genomes. Here we investigated if flavivirus EVEs are conserved across populations and different Aedes species by using ~ 500 whole genome sequence libraries from Aedes aegypti and Aedes albopictus, sourced from colonies and field mosquitoes across continents. We found that nearly all flavivirus EVEs in the Ae. aegypti reference genome originate from four separate putative viral integration events, and that they are highly conserved across geographically diverse samples. By contrast, flavivirus EVEs in the Ae. albopictus reference genome originate from up to nine distinct integration events and show low levels of conservation, even within samples from narrow geographical ranges. Our analysis suggests that flaviviruses integrated as long sequences and were subsequently fragmented and shuffled by transposable elements. Given that EVEs of Ae. aegypti and Ae. albopictus belong to different phylogenetic clades and have very differing levels of conservation, they may have different evolutionary origins and potentially different functional roles.

Genetic determinants of antiviral immunity in dipteran insects - Compiling the experimental evidence

The genetic basis of antiviral immunity in dipteran insects is extensively studied in Drosophila melanogaster and advanced technologies for genetic manipulation allow a better characterization of immune responses also in non-model insect species. Especially, immunity in vector mosquitoes is recently in the spotlight, due to the medical impact that these insects have by transmitting viruses and other pathogens. Here, we review the current state of experimental evidence that supports antiviral functions for immune genes acting in different cellular pathways. We discuss the well-characterized RNA interference mechanism along with the less well-defined JAK-STAT, Toll, and IMD signaling pathways. Furthermore, we highlight the initial evidence for antiviral activity observed for the autophagy pathway, transcriptional pausing, as well as piRNA production from endogenous viral elements. We focus our review on studies from Drosophila and mosquito species from the lineages Aedes, Culex, and Anopheles, which contain major vector species responsible for virus transmission.

Improving mosquito control strategies with population genomics

Mosquito control strategies increasingly apply knowledge from population genomics research. This review highlights recent applications to three research domains: mosquito invasions, insecticide resistance evolution, and rear and release programs. Current research trends follow developments in reference assemblies, either as improvements to existing assemblies (particularly Aedes) or assemblies for new taxa (particularly Anopheles). With improved assemblies, studies of invasive and rear and release target populations are better able to incorporate adaptive as well as demographic hypotheses. New reference assemblies are aiding comparisons of insecticide resistance across sister taxa while helping resolve taxon boundaries amidst frequent introgression. Anopheles gene drive deployments and improved Aedes genome assemblies should lead to a convergence in research aims for Anopheles and Aedes in the coming years.

Profile of Small RNAs, vDNA Forms and Viral Integrations in Late Chikungunya Virus Infection of Aedes albopictus Mosquitoes

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.

A mosquito small RNA genomics resource reveals dynamic evolution and host responses to viruses and transposons

Although mosquitoes are major transmission vectors for pathogenic arboviruses, viral infection has little impact on mosquito health. This immunity is caused in part by mosquito RNA interference (RNAi) pathways that generate antiviral small interfering RNAs (siRNAs) and Piwi-interacting RNAs (piRNAs). RNAi also maintains genome integrity by potently repressing mosquito transposon activity in the germline and soma. However, viral and transposon small RNA regulatory pathways have not been systematically examined together in mosquitoes. Therefore, we developed an integrated mosquito small RNA genomics (MSRG) resource that analyzes the transposon and virus small RNA profiles in mosquito cell cultures and somatic and gonadal tissues across four medically important mosquito species. Our resource captures both somatic and gonadal small RNA expression profiles within mosquito cell cultures, and we report the evolutionary dynamics of a novel Mosquito-Conserved piRNA Cluster Locus (MCpiRCL) made up of satellite DNA repeats. In the larger culicine mosquito genomes we detected highly regular periodicity in piRNA biogenesis patterns coinciding with the expansion of Piwi pathway genes. Finally, our resource enables detection of cross talk between piRNA and siRNA populations in mosquito cells during a response to virus infection. The MSRG resource will aid efforts to dissect and combat the capacity of mosquitoes to tolerate and spread arboviruses.

Chronological Incongruences between Mitochondrial and Nuclear Phylogenies of Aedes Mosquitoes

One-third of all mosquitoes belong to the Aedini, a tribe comprising common vectors of viral zoonoses such as Aedes aegypti and Aedes albopictus. To improve our understanding of their evolution, we present an updated multigene estimate of Aedini phylogeny and divergence, focusing on the disentanglement between nuclear and mitochondrial phylogenetic signals. We first show that there are some phylogenetic discrepancies between nuclear and mitochondrial markers which may be caused by wrong taxa assignment in samples collections or by some stochastic effect due to small gene samples. We indeed show that the concatenated dataset is model and framework dependent, indicating a general paucity of signal. Our Bayesian calibrated divergence estimates point toward a mosquito radiation in the mid-Jurassic and an Aedes radiation from the mid-Cretaceous on. We observe, however a strong chronological incongruence between mitochondrial and nuclear data, the latter providing divergence times within the Aedini significantly younger than the former. We show that this incongruence is consistent over different datasets and taxon sampling and that may be explained by either peculiar evolutionary event such as different levels of saturation in certain lineages or a past history of hybridization throughout the genus. Overall, our updated picture of Aedini phylogeny, reveal a strong nuclear-mitochondrial incongruence which may be of help in setting the research agenda for future phylogenomic studies of Aedini mosquitoes.

A Linkage-Based Genome Assembly for the Mosquito Aedes albopictus and Identification of Chromosomal Regions Affecting Diapause

The Asian tiger mosquito, Aedes albopictus, is an invasive vector mosquito of substantial public health concern. The large genome size (~1.19-1.28 Gb by cytofluorometric estimates), comprised of ~68% repetitive DNA sequences, has made it difficult to produce a high-quality genome assembly for this species. We constructed a high-density linkage map for Ae. albopictus based on 111,328 informative SNPs obtained by RNAseq. We then performed a linkage-map anchored reassembly of AalbF2, the genome assembly produced by Palatini et al. (2020). Our reassembled genome sequence, AalbF3, represents several improvements relative to AalbF2. First, the size of the AalbF3 assembly is 1.45 Gb, almost half the size of AalbF2. Furthermore, relative to AalbF2, AalbF3 contains a higher proportion of complete and single-copy BUSCO genes (84.3%) and a higher proportion of aligned RNAseq reads that map concordantly to a single location of the genome (46%). We demonstrate the utility of AalbF3 by using it as a reference for a bulk-segregant-based comparative genomics analysis that identifies chromosomal regions with clusters of candidate SNPs putatively associated with photoperiodic diapause, a crucial ecological adaptation underpinning the rapid range expansion and climatic adaptation of A. albopictus.

A Gene-Based Method for Cytogenetic Mapping of Repeat-Rich Mosquito Genomes

Long-read sequencing technologies have opened up new avenues of research on the mosquito genome biology, enabling scientists to better understand the remarkable abilities of vectors for transmitting pathogens. Although new genome mapping technologies such as Hi-C scaffolding and optical mapping may significantly improve the quality of genomes, only cytogenetic mapping, with the help of fluorescence in situ hybridization (FISH), connects genomic scaffolds to a particular chromosome and chromosome band. This mapping approach is important for creating and validating chromosome-scale genome assemblies for mosquitoes with repeat-rich genomes, which can potentially be misassembled. In this study, we describe a new gene-based physical mapping approach that was optimized using the newly assembled Aedes albopictus genome, which is enriched with transposable elements. To avoid amplification of the repetitive DNA, 15 protein-coding gene transcripts were used for the probe design. Instead of using genomic DNA, complementary DNA was utilized as a template for development of the PCR-amplified probes for FISH. All probes were successfully amplified and mapped to specific chromosome bands. The genome-unique probes allowed to perform unambiguous mapping of genomic scaffolds to chromosome regions. The method described in detail here can be used for physical genome mapping in other insects.

ViR: a tool to solve intrasample variability in the prediction of viral integration sites using whole genome sequencing data

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.