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

Analysis of somatic piRNAs in the malaria mosquito Anopheles coluzzii reveals atypical classes of genic small RNAs

Piwi-interacting small RNAs (piRNA) play a key role in controlling the activity of transposable elements (TEs) in the animal germline. In diverse arthropod species, including the pathogen vectors mosquitoes, the piRNA pathway is also active in nongonadal somatic tissues, where its targets and functions are less clear. Here, we studied the features of small RNA production in head and thorax tissues of an uninfected laboratory strain of Anopheles coluzzii focusing on the 24-32-nt-long RNAs. Small RNAs derived from repetitive elements constitute a minor fraction while most small RNAs process from long noncoding RNAs (lncRNAs) and protein-coding gene mRNAs. The majority of small RNAs derived from repetitive elements and lncRNAs exhibited typical piRNAs features. By contrast, majority of protein-coding gene-derived 24-32 nt small RNAs lack the hallmarks of piRNAs and have signatures of nontemplated 3' end tailing. Most of the atypical small RNAs exhibit female-biased expression and originate from mitochondrial and nuclear genes involved in energy metabolism. We also identified atypical genic small RNAs in Anopheles gambiae somatic tissues, which further validates the noncanonical mechanism of their production. We discuss a novel mechanism of small RNA production in mosquito somatic tissues and the possible functional significance of genic small RNAs.

Traffic Jam activates the Flamenco piRNA cluster locus and the Piwi pathway to ensure transposon silencing and Drosophila fertility

Flamenco (Flam) is a prominent Piwi-interacting RNA (piRNA) locus expressed in Drosophila ovarian follicle cells that silences gypsy/mdg4 transposons to ensure female fertility. Promoter-bashing reporter assays in ovarian somatic sheet (OSS) cells uncover compact enhancer sequences within Flam. We confirm the enhancer sequence relevance in vivo with Drosophila Flam deletion mutants that compromise Flam piRNA levels and female fertility. Proteomic analysis of proteins associated with Flam enhancer sequences discover the transcription factor Traffic Jam (TJ). Tj knockdown in OSS cells causes a decrease in Flam transcripts, Flam piRNAs, and multiple Piwi pathway genes. TJ chromatin immunoprecipitation sequencing (ChIP-seq) analysis confirms TJ binding at enhancer sequences deleted in our distinct Flam mutants. TJ also binds multiple Piwi pathway gene enhancers and long terminal repeats of transposons that decrease in expression after Tj knockdown. TJ plays an integral role in the ongoing arms race between selfish transposons and their suppression by the host Piwi pathway and Flam piRNA locus.

Somatic Miwi2 modulates mitochondrial function in airway multiciliated cells and exacerbates influenza pathogenesis

MIWI2, a P element-induced wimpy testes (PIWI) argonaute protein known for suppressing retrotransposons during male gonadogenesis, has an unexplored role in mammalian somatic cells. We identify MIWI2 multiciliated (M2MC) cells as a rare subset of airway multiciliated cells and investigate MIWI2's function in antiviral host defense. We analyzed transcriptomes from Miwi2 heterozygous (Miwi2 +/tom) and deficient (Miwi2 tom/tom) mice following influenza A infection. During infection, Miwi2 deficiency was associated with reduced mitochondrial and ribosomal gene expression in M2MC cells, increased mitochondrial reactive oxygen species (ROS) production and ADP/ATP ratios in multiciliated cells, and enhanced viral clearance and recovery. Additionally, Miwi2-expressing cells exhibited reduced levels of small RNAs derived from nuclear mitochondrial DNA. These findings reveal a previously unrecognized role for Miwi2 in regulating small non-coding RNAs and mitochondrial oxidant production in somatic cells, indicating a function beyond its established germline activities. Our study identifies Miwi2/Piwil4 as a potential factor influencing susceptibility to severe respiratory infections.

Lipid Peroxidation and Type I Interferon Coupling Fuels Pathogenic Macrophage Activation Causing Tuberculosis Susceptibility

A quarter of human population is infected with Mycobacterium tuberculosis, but less than 10% of those infected develop pulmonary TB. We developed a genetically defined sst1-susceptible mouse model that uniquely reproduces a defining feature of human TB: the development of necrotic lung granulomas and determined that the sst1-susceptible phenotype was driven by the aberrant macrophage activation. This study demonstrates that the aberrant response of the sst1-susceptible macrophages to prolonged stimulation with TNF is primarily driven by conflicting Myc and antioxidant response pathways leading to a coordinated failure 1) to properly sequester intracellular iron and 2) to activate ferroptosis inhibitor enzymes. Consequently, iron-mediated lipid peroxidation fueled Ifnβ superinduction and sustained the Type I Interferon (IFN-I) pathway hyperactivity that locked the sst1-susceptible macrophages in a state of unresolving stress and compromised their resistance to Mtb. The accumulation of the aberrantly activated, stressed, macrophages within granuloma microenvironment led to the local failure of anti-tuberculosis immunity and tissue necrosis. The upregulation of Myc pathway in peripheral blood cells of human TB patients was significantly associated with poor outcomes of TB treatment. Thus, Myc dysregulation in activated macrophages results in an aberrant macrophage activation and represents a novel target for host-directed TB therapies.

Prolonged exposure to heat enhances mosquito tolerance to viral infection

How and to what extent mosquito-virus interaction is influenced by climate change is a complex question of ecological and epidemiological relevance. We worked at the intersection between thermal biology and vector immunology and studied shifts in tolerance and resistance to the cell fusing agent virus (CFAV), a prominent component of the mosquito virome known to contribute to shaping mosquito vector competence, in warm-acclimated and warm-evolved Aedes albopictus mosquitoes. We show that the length of the thermal challenge influences the outcome of the infection with warm-evolved mosquitoes being more tolerant to CFAV infection, while warm-acclimated mosquitoes being more resistant and suffering from extensive fitness costs. These results highlight the importance of considering fluctuations in vector immunity in relation to the length of a thermal challenge to understand natural variation in vector response to viruses and frame realistic transmission models.

Traffic Jam activates the Flamenco piRNA cluster locus and the Piwi pathway to ensure transposon silencing and Drosophila fertility

Flamenco (Flam) is the most prominent piRNA cluster locus expressed in Drosophila ovarian follicle cells, and it is required for female fertility to silence gypsy/mdg4 transposons. To determine how Flam is regulated, we used promoter-bashing reporter assays in OSS cells to uncover novel enhancer sequences within the first exons of Flam . We confirmed the enhancer sequence relevance in vivo with new Drosophila Flam deletion mutants of these regions that compromised Flam piRNA expression and lowered female fertility from activated transposons. Our proteomic analysis of proteins associated with these enhancer sequences discovered the transcription factor Traffic Jam (TJ). Tj knockdowns in OSS cells caused a decrease in Flam transcripts, Flam piRNAs, and multiple Piwi pathway genes. A TJ ChIP-seq analysis from whole flies and OSS cells confirmed TJ binding exactly at the enhancer that was deleted in the new Flam mutant as well as at multiple Piwi pathway gene enhancers. Interestingly, TJ also bound the Long Terminal Repeats of transposons that had decreased expression after Tj knockdowns in OSS cells. Our study reveals the integral role TJ plays in the on-going arms race between selfish transposons and their suppression by the host Piwi pathway and the Flam piRNA cluster locus.

Rapid evolution of piRNA clusters in the Drosophila melanogaster ovary

The piRNA pathway is a highly conserved mechanism to repress transposable element (TE) activity in the animal germline via a specialized class of small RNAs called piwi-interacting RNAs (piRNAs). piRNAs are produced from discrete genomic regions called piRNA clusters (piCs). Although the molecular processes by which piCs function are relatively well understood in Drosophila melanogaster, much less is known about the origin and evolution of piCs in this or any other species. To investigate piC origin and evolution, we use a population genomic approach to compare piC activity and sequence composition across eight geographically distant strains of D. melanogaster with high-quality long-read genome assemblies. We perform annotations of ovary piCs and genome-wide TE content in each strain. Our analysis uncovers extensive variation in piC activity across strains and signatures of rapid birth and death of piCs. Most TEs inferred to be recently active show an enrichment of insertions into old and large piCs, consistent with the previously proposed "trap" model of piC evolution. In contrast, a small subset of active LTR families is enriched for the formation of new piCs, suggesting that these TEs have higher proclivity to form piCs. Thus, our findings uncover processes leading to the origin of piCs. We propose that piC evolution begins with the emergence of piRNAs from individual insertions of a few select TE families prone to seed new piCs that subsequently expand by accretion of insertions from most other TE families during evolution to form larger "trap" clusters. Our study shows that TEs themselves are the major force driving the rapid evolution of piCs.

Global genomics of Aedes aegypti unveils widespread and novel infectious viruses capable of triggering a small RNA response

The mosquito Aedes aegypti is a prominent vector for arboviruses, but the breadth of mosquito viruses that infects this specie is not fully understood. In the broadest global survey to date of over 200 Ae. aegypti small RNA samples, we detected viral small interfering RNAs (siRNAs) and Piwi interacting RNAs (piRNAs) arising from mosquito viruses. We confirmed that most academic laboratory colonies of Ae. aegypti lack persisting viruses, yet two commercial strains were infected by a novel tombus-like virus. Ae. aegypti from North to South American locations were also teeming with multiple insect viruses, with Anphevirus and a bunyavirus displaying geographical boundaries from the viral small RNA patterns. Asian Ae. aegypti small RNA patterns indicate infections by similar mosquito viruses from the Americas and reveal the first wild example of dengue virus infection generating viral small RNAs. African Ae. aegypti also contained various viral small RNAs including novel viruses only found in these African substrains. Intriguingly, viral long RNA patterns can differ from small RNA patterns, indicative of viral transcripts evading the mosquitoes' RNA interference (RNAi) machinery. To determine whether the viruses we discovered via small RNA sequencing were replicating and transmissible, we infected C6/36 and Aag2 cells with Ae. aegypti homogenates. Through blind passaging, we generated cell lines stably infected by these mosquito viruses which then generated abundant viral siRNAs and piRNAs that resemble the native mosquito viral small RNA patterns. This mosquito small RNA genomics approach augments surveillance approaches for emerging infectious diseases.

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.

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.

Aedes aegypti Argonaute 2 controls arbovirus infection and host mortality

Ae. aegypti mosquitoes transmit some of the most important human viral diseases that are responsible for a significant public health burden worldwide. The small interfering RNA (siRNA) pathway is considered the major antiviral defense system in insects. Here we show that siRNA pathway disruption by CRISPR/Cas9-based Ago2 knockout impaired the mosquitoes' ability to degrade arbovirus RNA leading to hyper-infection accompanied by cell lysis and tissue damage. Ago2 disruption impaired DNA repair mechanisms and the autophagy pathway by altering histone abundance. This compromised DNA repair and removal of damaged cellular organelles and dysfunctional aggregates promoted mosquito death. We also report that hyper-infection of Ago2 knockout mosquitoes stimulated a broad-spectrum antiviral immunity, including apoptosis, which may counteract infection. Taken together, our studies reveal novel roles for Ago2 in protecting mosquitoes from arbovirus infection and associated death.

Rapid evolution of piRNA clusters in the Drosophila melanogaster ovary

Animal genomes are parasitized by a horde of transposable elements (TEs) whose mutagenic activity can have catastrophic consequences. The piRNA pathway is a conserved mechanism to repress TE activity in the germline via a specialized class of small RNAs associated with effector Piwi proteins called piwi-associated RNAs (piRNAs). piRNAs are produced from discrete genomic regions called piRNA clusters (piCs). While piCs are generally enriched for TE sequences and the molecular processes by which they are transcribed and regulated are relatively well understood in Drosophila melanogaster, much less is known about the origin and evolution of piCs in this or any other species. To investigate piC evolution, we use a population genomics approach to compare piC activity and sequence composition across 8 geographically distant strains of D. melanogaster with high quality long-read genome assemblies. We perform extensive annotations of ovary piCs and TE content in each strain and test predictions of two proposed models of piC evolution. The 'de novo' model posits that individual TE insertions can spontaneously attain the status of a small piC to generate piRNAs silencing the entire TE family. The 'trap' model envisions large and evolutionary stable genomic clusters where TEs tend to accumulate and serves as a long-term "memory" of ancient TE invasions and produce a great variety of piRNAs protecting against related TEs entering the genome. It remains unclear which model best describes the evolution of piCs. Our analysis uncovers extensive variation in piC activity across strains and signatures of rapid birth and death of piCs in natural populations. Most TE families inferred to be recently or currently active show an enrichment of strain-specific insertions into large piCs, consistent with the trap model. By contrast, only a small subset of active LTR retrotransposon families is enriched for the formation of strain-specific piCs, suggesting that these families have an inherent proclivity to form de novo piCs. Thus, our findings support aspects of both 'de novo' and 'trap' models of piC evolution. We propose that these two models represent two extreme stages along an evolutionary continuum, which begins with the emergence of piCs de novo from a few specific LTR retrotransposon insertions that subsequently expand by accretion of other TE insertions during evolution to form larger 'trap' clusters. Our study shows that piCs are evolutionarily labile and that TEs themselves are the major force driving the formation and evolution of piCs.

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.

Culex Mosquito Piwi4 Is Antiviral against Two Negative-Sense RNA Viruses

Culex spp. mosquitoes transmit several pathogens concerning public health, including West Nile virus and Saint Louis encephalitis virus. Understanding the antiviral immune system of Culex spp. mosquitoes is important for reducing the transmission of these viruses. Mosquitoes rely on RNA interference (RNAi) to control viral replication. While the siRNA pathway in mosquitoes is heavily studied, less is known about the piRNA pathway. The piRNA pathway in mosquitoes has recently been connected to mosquito antiviral immunity. In Aedes aegypti, Piwi4 has been implicated in antiviral responses. The antiviral role of the piRNA pathway in Culex spp. mosquitoes is understudied compared to Ae. aegypti. Here, we aimed to identify the role of PIWI genes and piRNAs in Culex quinquefasciatus and Culex tarsalis cells during virus infection. We examined the effect of PIWI gene silencing on virus replication of two arboviruses and three insect-specific viruses in Cx. quinquefasciatus derived cells (Hsu) and Cx. tarsalis derived (CT) cells. We show that Piwi4 is antiviral against the La Crosse orthobunyavirus (LACV) in Hsu and CT cells, and the insect-specific rhabdovirus Merida virus (MERDV) in Hsu cells. None of the silenced PIWI genes impacted replication of the two flaviviruses Usutu virus (USUV) and Calbertado virus, or the phasivirus Phasi-Charoen-like virus. We further used small RNA sequencing to determine that LACV-derived piRNAs, but not USUV-derived piRNAs were generated in Hsu cells and that PIWI gene silencing resulted in a small reduction in vpiRNAs. Finally, we determined that LACV-derived DNA was produced in Hsu cells during infection, but whether this viral DNA is required for vpiRNA production remains unclear. Overall, we expanded our knowledge on the piRNA pathway and how it relates to the antiviral response in Culex spp mosquitoes.

Vertical and Horizontal Transmission of Cell Fusing Agent Virus in Aedes aegypti

Cell fusing agent virus (CFAV) is an insect-specific flavivirus (ISF) found in Aedes aegypti mosquitoes. ISFs have demonstrated the ability to modulate the infection or transmission of arboviruses such as dengue, West Nile, and Zika viruses. It is thought that vertical transmission is the main route for ISF maintenance in nature. This has been observed with CFAV, but there is evidence of horizontal and venereal transmission in other ISFs. Understanding the route of transmission can inform strategies to spread ISFs to vector populations as a method of controlling pathogenic arboviruses. We crossed individually reared male and female mosquitoes from both a naturally occurring CFAV-positive Ae. aegypti colony and its negative counterpart to provide information on maternal, paternal, and horizontal transmission. RT-PCR was used to detect CFAV in individual female pupal exuviae and was 89% sensitive, but only 42% in male pupal exuviae. This is a possible way to screen individuals for infection without destroying the adults. Female-to-male horizontal transmission was not observed during this study. However, there was a 31% transmission rate from mating pairs of CFAV-positive males to negative female mosquitoes. Maternal vertical transmission was observed with a filial infection rate of 93%. The rate of paternal transmission was 85% when the female remained negative, 61% when the female acquired CFAV horizontally, and 76% overall. Maternal and paternal transmission of CFAV could allow the introduction of this virus into wild Ae. aegypti populations through male or female mosquito releases, and thus provides a potential strategy for ISF-derived arbovirus control. IMPORTANCE Insect-specific flaviviruses (ISFs), are a group of nonpathogenic flaviviruses that only infect insects. ISFs can have a high prevalence in mosquito populations, but their transmission routes are not well understood. The results of this study confirm maternal transmission of cell fusing agent virus (CFAV) and demonstrate that paternal transmission is also highly efficient. Horizontal transmission of CFAV was also observed, aided by evaluation of the pupal infection status before mating with an infected individual. This technique of detecting infection in discarded pupae exuviae has not been evaluated previously and will be a useful tool for others in the field of studying viral transmission in mosquitoes. Identifying these routes of transmission provides information about how CFAV could be maintained in wild populations of mosquitoes and can aid future studies focusing on interactions of CFAV with their hosts and other viruses that infect mosquitoes.

MeShClust v3.0: high-quality clustering of DNA sequences using the mean shift algorithm and alignment-free identity scores

Background

Tools for accurately clustering biological sequences are among the most important tools in computational biology. Two pioneering tools for clustering sequences are CD-HIT and UCLUST, both of which are fast and consume reasonable amounts of memory; however, there is a big room for improvement in terms of cluster quality. Motivated by this opportunity for improving cluster quality, we applied the mean shift algorithm in MeShClust v1.0. The mean shift algorithm is an instance of unsupervised learning. Its strong theoretical foundation guarantees the convergence to the true cluster centers. Our implementation of the mean shift algorithm in MeShClust v1.0 was a step forward. In this work, we scale up the algorithm by adapting an out-of-core strategy while utilizing alignment-free identity scores in a new tool: MeShClust v3.0.

Results

We evaluated CD-HIT, MeShClust v1.0, MeShClust v3.0, and UCLUST on 22 synthetic sets and five real sets. These data sets were designed or selected for testing the tools in terms of scalability and different similarity levels among sequences comprising clusters. On the synthetic data sets, MeShClust v3.0 outperformed the related tools on all sets in terms of cluster quality. On two real data sets obtained from human microbiome and maize transposons, MeShClust v3.0 outperformed the related tools by wide margins, achieving 55%–300% improvement in cluster quality. On another set that includes degenerate viral sequences, MeShClust v3.0 came third. On two bacterial sets, MeShClust v3.0 was the only applicable tool because of the long sequences in these sets. MeShClust v3.0 requires more time and memory than the related tools; almost all personal computers at the time of this writing can accommodate such requirements. MeShClust v3.0 can estimate an important parameter that controls cluster membership with high accuracy.

Conclusions

These results demonstrate the high quality of clusters produced by MeShClust v3.0 and its ability to apply the mean shift algorithm to large data sets and long sequences. Because clustering tools are utilized in many studies, providing high-quality clusters will help with deriving accurate biological knowledge.

Supplementary Information

The online version contains supplementary material available at (10.1186/s12864-022-08619-0).

VectorBase.org updates: bioinformatic resources for invertebrate vectors of human pathogens and related organisms

VectorBase (VectorBase.org) is part of the VEuPathDB Bioinformatics Resource Center, providing free online access to multi-omics and population biology data, focusing on arthropod vectors and invertebrates of importance to human health. VectorBase includes genomics and functional genomics data from bed bugs, biting midges, body lice, kissing bugs, mites, mosquitoes, sand flies, ticks, tsetse flies, stable flies, house flies, fruit flies, and a snail intermediate host. Tools include the Search Strategy system and MapVEu, enabling users to interrogate and visualize diverse 'omics and population-level data using a graphical interface (no programming experience required). Users can also analyze their own private data, such as transcriptomic sequences, exploring their results in the context of other publicly-available information in the database. Help Desk: help@vectorbase.org.

Mucus piRNAs profiles of Vibrio harveyi-infected Cynoglossus semilaevis: A hint for fish disease monitoring

The half-smooth tongue sole, Cynoglossus semilaevis, is an important cultured flatfish species. Vibrio harveyi is a common pathogen to this fish, which may result in great economic loss to C. semilaevis culture industry. piRNAs, a non-coding RNAs with 26-32 nt, have been regarded as promising biomarkers for cancer diagnosis and fish diseases. Here, we extracted the RNA from mucus of C. semilaevis and constructed the differential expression profiles of piRNAs between the sick fish (MS) and healthy fish (MC). We identified 45,696 differentially expressed piRNAs including 22,735 up-regulated piRNAs and 22,961 down-regulated piRNAs in MS group compared with MC group. The GO enrichment and KEGG pathway enrichment analyses of the differential piRNAs were carried out. The result showed immunity-related target genes mainly involved in immune system process, response to stimulus, cell killing, immune system, infectious diseases and cell growth and death. The 10 most differentially expressed piRNAs were chosen to perform the qRT-PCR, while only seven piRNAs were consistent with the sequence result. Compared with MC group, the expression levels of piR-mmu-72173>piR-rno-62831>piR-xtr-704880, piR-dme-15546979, piR-mmu-49941660, piR-mmu-29283297 and piR-mmu-1758399 were significantly lower, and piR-gga-10574 and piR-gga-134812 were significantly higher in MS group. These piRNAs may be potential biomarkers during the V. harveyi infection of C. semilaevis. This study could provide a new method to identify the infection status of C. semilaevis and understand better about the innate and adaptive immune system in C. semilaevis during bacterial infection.

Extending and Running the Mosquito Small RNA Genomics Resource Pipeline

The Mosquito Small RNA Genomics (MSRG) resource is a repository of analyses on the small RNA transcriptomes of mosquito cell cultures and somatic and gonadal tissues. This resource allows for comparing the regulation dynamics of small RNAs generated from transposons and viruses across mosquito species. This chapter covers the procedures to set up the MSRG resource pipeline as a new installation by detailing the necessary collection of genome reference and annotation files and lists of microRNAs (miRNAs) hairpin sequences, transposon repeats consensus sequences, and virus genome sequences. Proper execution of the MSRG resource pipeline yields outputs amenable to biologists to further analyze with desktop and spreadsheet software to gain insights into the balance between arthropod endogenous small RNA populations and the proportions of virus-derived small RNAs that include Piwi-interacting RNAs (piRNAs) and endogenous small interfering RNAs (siRNAs).

Aedes aegypti Piwi4 Structural Features Are Necessary for RNA Binding and Nuclear Localization

The PIWI-interacting RNA (piRNA) pathway provides an RNA interference (RNAi) mechanism known from Drosophila studies to maintain the integrity of the germline genome by silencing transposable elements (TE). Aedes aegypti mosquitoes, which are the key vectors of several arthropod-borne viruses, exhibit an expanded repertoire of Piwi proteins involved in the piRNA pathway, suggesting functional divergence. Here, we investigate RNA-binding dynamics and subcellular localization of A. aegypti Piwi4 (AePiwi4), a Piwi protein involved in antiviral immunity and embryonic development, to better understand its function. We found that AePiwi4 PAZ (Piwi/Argonaute/Zwille), the domain that binds the 3′ ends of piRNAs, bound to mature (3′ 2′ O-methylated) and unmethylated RNAs with similar micromolar affinities (K_D = 1.7 ± 0.8 μM and K_D of 5.0 ± 2.2 μM, respectively; p = 0.05) in a sequence independent manner. Through site-directed mutagenesis studies, we identified highly conserved residues involved in RNA binding and found that subtle changes in the amino acids flanking the binding pocket across PAZ proteins have significant impacts on binding behaviors, likely by impacting the protein secondary structure. We also analyzed AePiwi4 subcellular localization in mosquito tissues. We found that the protein is both cytoplasmic and nuclear, and we identified an AePiwi4 nuclear localization signal (NLS) in the N-terminal region of the protein. Taken together, these studies provide insights on the dynamic role of AePiwi4 in RNAi and pave the way for future studies aimed at understanding Piwi interactions with diverse RNA populations.

Bioinformatic and cell-based tools for pooled CRISPR knockout screening in mosquitos

Mosquito-borne diseases present a worldwide public health burden. Current efforts to understand and counteract them have been aided by the use of cultured mosquito cells. Moreover, application in mammalian cells of forward genetic approaches such as CRISPR screens have identified essential genes and genes required for host-pathogen interactions, and in general, aided in functional annotation of genes. An equivalent approach for genetic screening of mosquito cell lines has been lacking. To develop such an approach, we design a new bioinformatic portal for sgRNA library design in several mosquito genomes, engineer mosquito cell lines to express Cas9 and accept sgRNA at scale, and identify optimal promoters for sgRNA expression in several mosquito species. We then optimize a recombination-mediated cassette exchange system to deliver CRISPR sgRNA and perform pooled CRISPR screens in an Anopheles cell line. Altogether, we provide a platform for high-throughput genome-scale screening in cell lines from disease vector species.

Small RNA Response to Infection of the Insect-Specific Lammi Virus and Hanko Virus in an Aedes albopictus Cell Line

RNA interference (RNAi)-mediated antiviral immunity is believed to be the primary defense against viral infection in mosquitoes. The production of virus-specific small RNA has been demonstrated in mosquitoes and mosquito-derived cell lines for viruses in all of the major arbovirus families. However, many if not all mosquitoes are infected with a group of viruses known as insect-specific viruses (ISVs), and little is known about the mosquito immune response to this group of viruses. Therefore, in this study, we sequenced small RNA from an Aedes albopictus-derived cell line infected with either Lammi virus (LamV) or Hanko virus (HakV). These viruses belong to two distinct phylogenetic groups of insect-specific flaviviruses (ISFVs). The results revealed that both viruses elicited a strong virus-derived small interfering RNA (vsiRNA) response that increased over time and that targeted the whole viral genome, with a few predominant hotspots observed. Furthermore, only the LamV-infected cells produced virus-derived Piwi-like RNAs (vpiRNAs); however, they were mainly derived from the antisense genome and did not show the typical ping-pong signatures. HakV, which is more distantly related to the dual-host flaviviruses than LamV, may lack certain unknown sequence elements or structures required for vpiRNA production. Our findings increase the understanding of mosquito innate immunity and ISFVs' effects on their host.

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.

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.

Argonaute-CLIP delineates versatile, functional RNAi networks in Aedes aegypti, a major vector of human viruses

Argonaute (AGO) proteins bind small RNAs to silence complementary RNA transcripts, and they are central to RNA interference (RNAi). RNAi is critical for regulation of gene expression and antiviral defense in Aedes aegypti mosquitoes, which transmit Zika, chikungunya, dengue, and yellow fever viruses. In mosquitoes, AGO1 mediates miRNA interactions, while AGO2 mediates siRNA interactions. We applied AGO-crosslinking immunoprecipitation (AGO-CLIP) for both AGO1 and AGO2, and we developed a universal software package for CLIP analysis (CLIPflexR), identifying 230 small RNAs and 5,447 small RNA targets that comprise a comprehensive RNAi network map in mosquitoes. RNAi network maps predicted expression levels of small RNA targets in specific tissues. Additionally, this resource identified unexpected, context-dependent AGO2 target preferences, including endogenous viral elements and 3'UTRs. Finally, contrary to current thinking, mosquito AGO2 repressed imperfect targets. These findings expand our understanding of small RNA networks and have broad implications for the study of antiviral RNAi.

Leaning Into the Bite: The piRNA Pathway as an Exemplar for the Genetic Engineering Need in Mosquitoes

The piRNA pathway is a specialized small RNA interference that in mosquitoes is mechanistically distant from analogous biology in the Drosophila model. Current genetic engineering methods, such as targeted genome manipulation, have a high potential to tease out the functional complexity of this intricate molecular pathway. However, progress in utilizing these methods in arthropod vectors has been geared mostly toward the development of new vector control strategies rather than to study cellular functions. Herein we propose that genetic engineering methods will be essential to uncover the full functionality of PIWI/piRNA biology in mosquitoes and that extending the applications of genetic engineering on other aspects of mosquito biology will grant access to a much larger pool of knowledge in disease vectors that is just out of reach. We discuss motivations for and impediments to expanding the utility of genetic engineering to study the underlying biology and disease transmission and describe specific areas where efforts can be placed to achieve the full potential for genetic engineering in basic biology in mosquito vectors. Such efforts will generate a refreshed intellectual source of novel approaches to disease control and strong support for the effective use of approaches currently in development.