The Diversity, Structure, and Function of Heritable Adaptive Immunity Sequences in the Aedes aegypti Genome

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.

Diversity and evolution analysis of RNA viruses in three wheat aphid species

BACKGROUND

Although advances in metagenomics, viral diversity and non-retroviral endogenous viral elements (EVEs) in wheat aphids remain underexplored. By analyzing 470 publicly available datasets and one laboratory-generated transcriptome, the RNA virome and EVEs in the genomes of Sitobion avenae, Schizaphis graminum, and Rhopalosiphum padi were systematically investigated.

RESULTS

We identified 43 RNA viruses, including 12 novel and 31 known RNA viruses. These viruses were widely distributed and abundant in different geographic populations of three wheat aphid species. +ssRNA viruses were the dominant type of aphid viruses. Besides, 90 EVEs were discovered in the genomes of three aphid species. In addition, the EVEs exhibit potential domestication and novel functional roles within aphid genomes.

CONCLUSIONS

This study expands the understanding of RNA virus diversity in aphids and provides valuable insights into the potential functions of EVEs in virus-host coevolution.

piRNA/PIWI pathways and epigenetic crosstalk in human diseases: Molecular insights into HIV-1 infection and drugs of abuse

P-element-induced wimpy (PIWI)-interacting RNAs (piRNAs) and PIWI proteins have long been studied in insects and germline cells for their roles in regulating transposable elements (TEs). However, emerging evidence suggests that piRNAs and PIWI proteins also play crucial roles in human diseases beyond gametocyte protection, and these molecules are implicated in the onset and progression of various human diseases, particularly those arising in somatic cells. Notably, piRNAs and PIWI proteins are increasingly recognized for their involvement in cancers, cardiovascular diseases, neurodegenerative disorders, and viral infections, including HIV. This review first provides an overview of piRNAs/PIWIs and their interactions with TEs and primary targets. We then explore the molecular mechanisms and signaling pathways through which piRNAs and PIWIs modulate human disease processes, focusing on neurodegeneration, cancers, and HIV. Special attention is given to the role of piRNA/PIWI complexes in regulating gene transcription, translation, and post-translational modifications in the context of disease. Additionally, we address emerging research into the role of piRNAs/PIWIs in HIV- and drug abuse or substance abuse-associated neurodegenerative diseases, highlighting existing knowledge gaps. Finally, we discuss future research directions to understand better the functions of piRNAs/PIWI proteins in human health and disease.

From macro to micro: De novo genomes of Aedes mosquitoes enable comparative genomics among close and distant relatives

The yellow fever mosquito (Aedes aegypti) is an organism of high medical importance because it is the primary vector for diseases such as yellow fever, Zika, dengue, and chikungunya. Its medical importance has made it a subject of numerous efforts to understand their biology. One such effort, was the development of a high-quality reference genome (AaegL5). However, this reference genome was sourced from a highly inbred laboratory strain with unknown geographic origin. Thus, the reference is not representative of a wild mosquito, let alone one from its native range in sub-Saharan Africa. To better understand the genetic architecture of Ae. aegypti and their sister species, we developed two de novo chromosome-scale genomes with sequences sourced from single individuals: one of Ae. aegypti formosus (Aaf) from Burkina Faso and one of Ae. mascarensis (Am) from Mauritius. Both genomes exhibit high contiguity and gene completeness, comparable to AaegL5. While Aaf exhibits high degree of synteny to AaegL5, it also exhibits several large inversions. We further conducted comparative genomic analyses using our genomes and other publicly available culicid reference genomes to find extensive chromosomal rearrangements between major lineages. Overrepresentation analysis of expanded genes in Aaf, AaegL5, and Am revealed that while the overarching category of genes that have expanded are similar, the specific genes that have expanded differ. Our findings elucidate novel insights into chromosome evolution at both microevolutionary and macroevolutionary scales. The genomic resources we present are additions to the arsenal of biologists in understanding mosquito biology and genome evolution.

EEfinder, a general purpose tool for identification of bacterial and viral endogenized elements in eukaryotic genomes

Horizontal gene transfer is a phenomenon of genetic material transmission between species with no parental relationship. It has been characterized among several major branches of life, including among prokaryotes, viruses and eukaryotes. The characterization of endogenous elements derived from viruses or bacteria provides a snapshot of past host-pathogen interactions and coevolution as well as reference information to remove false positive results from metagenomic studies. Currently there is a lack of general purpose standardized tools for endogenous elements screening which limits reproducibility and hinder comparative analysis between studies. Here we describe EEfinder, a new general purpose tool for identification and classification of endogenous elements derived from viruses or bacteria found in eukaryotic genomes. The tool was developed to include six common steps performed in this type of analysis: data cleaning, similarity search through sequence alignment, filtering candidate elements, taxonomy assignment, merging of truncated elements and flanks extraction. We evaluated the sensitivity of EEfinder to identify endogenous elements through comparative analysis using data from the literature and showed that EEfinder automatically detected 97 % of the EVEs compared to published results obtained by manual curation and detected an almost exact full integration of a Wolbachia genome described using wet-lab experiments. Therefore, EEfinder can effectively and systematically identify endogenous elements with bacterial/viral origin integrated in eukaryotic genomes. EEfinder is publicly available on https://github.com/WallauBioinfo/EEfinder.

Evidence of endogenous non-retroviral RNA virus sequences into the genome and transcriptome of the malaria vector Anopheles darlingi

The characterization of non-retroviral integrated RNA virus sequences (NIRVS) in mosquitoes has emerged as a significant area of research that could yield insight into virus-host interactions. This study aimed to characterize NIRVS in the Anopheles darlingi reference genome and identify putative transcribed NIRVS in field-collected mosquitoes from Colombia. The An. darlingi reference genome was analyzed to identify and characterize NIRVS by conducting a BLAST query with all the virus sequences previously identified in arthropods available in the NCBI-virus repository. In addition, An. darlingi field-collected mosquitoes were examined for NIRVS using a metatranscriptomic approach. As a result, 44 NIRVS were identified in the An. darlingi genome, constituting integrations of negative single-stranded RNA viruses (ssRNA-) from the families Rhabdoviridae, Chuviridae and Phasmaviridae, and integrations of double-stranded RNA viruses (dsRNA) from the families Partitiviridae and Sedoreoviridae. These NIRVS were not randomly distributed but clustered in specific regions of the genome enriched with BEL/Pao and Ty3/Gypsy long terminal repeat elements. Furthermore, putative NIRVS-like sequences were present in the transcriptomic data from all the Colombian An. darlingi natural populations. This study is significant as it represents the first identification of NIRVS in the most important malaria vector of the Neotropics. The findings help in understanding the intricate relationship between the mosquito and its virome, and the regulation of viruses' mechanisms in the Anopheles genus.

Endogenous viral elements: insights into data availability and accessibility

Endogenous viral elements (EVEs) are remnants of viral genetic material endogenized into the host genome. They have, in the last decades, attracted attention for their role as potential contributors to pathogenesis, drivers of selective advantage for the host, and genomic remnants of ancient viruses. EVEs have a nuanced and complex influence on both host health and evolution, and can offer insights on the deep evolutionary history of viruses. As an emerging field of research, several factors limit a comprehensive understanding of EVEs: they are currently underestimated and periodically overlooked in studies of the host genome, transcriptome, and virome. The absence of standardized guidelines for ensuring EVE-related data availability and accessibility following the FAIR ('findable, accessible, interoperable, and reusable') principles obstructs our ability to gather and connect information. Here, we discuss challenges to the availability and accessibility of EVE-related data and propose potential solutions. We identified the biological and research focus imbalance between different types of EVEs, and their overall biological complexity as genomic loci with viral ancestry, as potential challenges that can be addressed with the development of a user-oriented identification tool. In addition, reports of EVE identification are scattered between different subfields under different keywords, and EVE sequences and associated data are not properly gathered in databases. While developing an open and dedicated database might be ideal, targeted improvements of generalist databases might provide a pragmatic solution to EVE data and metadata accessibility. The implementation of these solutions, as well as the collective effort by the EVE scientific community in discussing and setting guidelines, is now drastically needed to lead the development of EVE research and offer insights into host-virus interactions and their evolutionary history.

Rich diversity of RNA viruses in the biological control agent, Orius laevigatus

Orius laevigatus (Hemiptera, Anthocoridae) is a generalist predator extensively used for the biocontrol of diverse agricultural pests. Previous studies on O. laevigatus have focused on the improvement of insect genetic traits, but little is known about its association with microbes, especially viruses that may influence its production and efficacy. More than 280 RNA viruses have been described in other Hemiptera insects, in line with the continuous discovery of insect-specific viruses (ISVs) boosted by next-generation sequencing. In this study, we characterized the repertoire of RNA viruses associated with O. laevigatus. Its virome comprises 27 RNA viruses, classified within fourteen viral families, of which twenty-three viruses are specific to O. laevigatus and four are likely associated with fungal microbiota. The analysis of viral abundance in five O. laevigatus populations confirmed the presence of simultaneous viral infections and highlighted the ubiquitous presence and high abundance of one solinvivirus and three totiviruses. Moreover, we identified 24 non-retroviral endogenous viral elements (nrEVEs) in the genome of O. laevigatus, suggesting a long-term relationship between the host and its virome. Although no symptoms were described in the insect populations under study, the high diversity of viral species and the high abundance of certain RNA viruses identified indicate that RNA viruses may be significant for the applicability and efficacy of O. laevigatus in biocontrol programs.

Innate and germline immune memory: specificity and heritability of the ancient immune mechanisms for adaptation and survival

The immune memory is one of the defensive strategies developed by both unicellular and multicellular organisms for ensuring their integrity and functionality. While the immune memory of the vertebrate adaptive immune system (based on somatic recombination) is antigen-specific, encompassing the generation of memory T and B cells that only recognize/react to a specific antigen epitope, the capacity of vertebrate innate cells to remember past events is a mostly non-specific mechanism of adaptation. This "innate memory" can be considered as germline-encoded because its effector tools (such as innate receptors) do not need somatic recombination for being active. Also, in several organisms the memory-related information is integrated in the genome of germline cells and can be transmitted to the progeny for several generations, but it can also be erased depending on the environmental conditions. Overall, depending on the organism, its environment and its living habits, innate immune memory appears to be a mechanism for achieving better protection and survival against repeated exposure to microbes/stressful agents present in the same environment or occurring in the same anatomical district, able to adapt to changes in the environmental cues. The anatomical and functional complexity of the organism and its lifespan drive the generation of different immune memory mechanisms, for optimal adaptation to changes in the living/environmental conditions. The concept of innate immunity being non-specific needs to be revisited, as a wealth of evidence suggests a significant degree of specificity both in the primary immune reaction and in the ensuing memory-like responses. This is clearly evident in invertebrate metazoans, in which distinct scenarios can be observed, with both non-specific (immune enhancement) or specific (immune priming) memory-like responses. In the case of mammals, there is evidence that some degree of specificity can be attained in different situations, for instance as organ-specific protection rather than microorganism-specific reaction. Thus, depending on the challenges and conditions, innate memory can be non-specific or specific, can be integrated in the germline and transmitted to the progeny or be short-lived, thereby representing an exceptionally plastic mechanism of defensive adaptation for ensuring individual and species survival.

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.

Effects of Arboviral Infections on Transposable Element Transcript Levels in Aedes aegypti

Transposable elements are mobile repeated sequences found in all genomes. Transposable elements are controlled by RNA interference pathways in most organisms, and this control involves the PIWI-interacting RNA pathway and the small interfering RNA pathway, which is also known to be the first line of antiviral defense in invertebrates. Using Drosophila, we recently showed that viral infections result in the modulation of transposable element transcript levels through modulation of the small RNA repertoire. The Aedes aegypti mosquito is of particular interest because almost half of its genome is made of transposable elements, and it is described as a major vector of viruses (such as the dengue [DENV], Zika [ZIKV], and chikungunya [CHIKV] arboviruses). Moreover, Aedes mosquitoes are unique among insects in that the PIWI-interacting RNA pathway is also involved in the somatic antiviral response, in addition to the transposable element control and PIWI-interacting RNA pathway genes expanded in the mosquito genome. For these reasons, we studied the impacts of viral infections on transposable element transcript levels in A. aegypti samples. We retrieved public datasets corresponding to RNA-seq data obtained from viral infections by DENV, ZIKV, and CHIKV in various tissues. We found that transposable element transcripts are moderately modulated following viral infection and that the direction of the modulation varies greatly across tissues and viruses. These results highlight the need for an in-depth investigation of the tightly intertwined interactions between transposable elements and viruses.

In vitro characterization of cell-fusing agent virus DNA forms in Aedes aegypti mosquitoes

How non-retroviral endogenous viral elements (EVEs) are established is a long-standing question. Viral DNA (vDNA) forms of RNA viruses are likely to be EVE precursors. Cell-fusing agent virus (CFAV) is a major insect-specific virus (ISV) in the Aedes aegypti mosquitoes and one of the few existing non-retroviral RNA viruses found as EVEs. We characterized CFAV-derived vDNA in the cell line to understand the mechanism of why current viruses are rarely endogenized. vDNA production was affected by cell culture media independent of CFAV replication. vDNAs that correspond to different regions covering the entire viral genome were detected, implying multiple initiation sites exist. A considerable proportion of vDNAs corresponded to ssDNA. Higher vDNA copies were detected in the cytoplasm than the nucleus. Our findings provide valuable insights into the intracellular characteristics of ISV-derived vDNAs, which will aid in understanding the underlying mechanisms of non-retroviral EVE formation.

Mosquito transgenerational antiviral immunity is mediated by vertical transfer of virus DNA sequences and RNAi

Mosquitoes are important vectors for transmission of many viruses of public and veterinary health concern. These viruses most commonly have an RNA genome and infect mosquitoes for life. The principal mosquito antiviral response is the RNAi system which destroys virus RNA. Here, we confirm an earlier study that Aedes aegypti mosquitoes infected with positive-stranded RNA arboviruses can transmit specific immunity to their offspring. We show that this trans-generational immunity requires replication of virus RNA and reverse transcription of vRNA to vDNA in the infected parents and intergenerational transfer of vDNA. This vDNA is both genome-integrated and episomal. The episomal vDNA sequences are flanked by retrotransposon long-terminal repeats, predominantly Copia-like. Integrated vDNA sequences are propagated along several generations but specific immunity is effective only for a few generations and correlates with the presence of vRNA and episomal vDNA. This understanding raises new possibilities for the control of important mosquito-borne virus diseases.

Invertebrate Immunity, Natural Transplantation Immunity, Somatic and Germ Cell Parasitism, and Transposon Defense

While the vertebrate immune system consists of innate and adaptive branches, invertebrates only have innate immunity. This feature makes them an ideal model system for studying the cellular and molecular mechanisms of innate immunity sensu stricto without reciprocal interferences from adaptive immunity. Although invertebrate immunity is evolutionarily older and a precursor of vertebrate immunity, it is far from simple. Despite lacking lymphocytes and functional immunoglobulin, the invertebrate immune system has many sophisticated mechanisms and features, such as long-term immune memory, which, for decades, have been exclusively attributed to adaptive immunity. In this review, we describe the cellular and molecular aspects of invertebrate immunity, including the epigenetic foundation of innate memory, the transgenerational inheritance of immunity, genetic immunity against invading transposons, the mechanisms of self-recognition, natural transplantation, and germ/somatic cell parasitism.

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.

Intrinsic factors driving mosquito vector competence and viral evolution: a review

Mosquitoes are responsible for the transmission of numerous viruses of global health significance. The term "vector competence" describes the intrinsic ability of an arthropod vector to transmit an infectious agent. Prior to transmission, the mosquito itself presents a complex and hostile environment through which a virus must transit to ensure propagation and transmission to the next host. Viruses imbibed in an infectious blood meal must pass in and out of the mosquito midgut, traffic through the body cavity or hemocoel, invade the salivary glands, and be expelled with the saliva when the vector takes a subsequent blood meal. Viruses encounter physical, cellular, microbial, and immunological barriers, which are influenced by the genetic background of the mosquito vector as well as environmental conditions. Collectively, these factors place significant selective pressure on the virus that impact its evolution and transmission. Here, we provide an overview of the current state of the field in understanding the mosquito-specific factors that underpin vector competence and how each of these mechanisms may influence virus evolution.

Non-retroviral Endogenous Viral Elements in Tephritid Fruit Flies Reveal Former Viral Infections Not Related to Known Circulating Viruses

A wide variety of insect-specific non-retroviral RNA viruses specifically infect insects. During viral infection, fragments of viral sequences can integrate into the host genomes creating non-retroviral endogenous viral elements (nrEVEs). Although the exact function of nrEVEs is so far unknown, some studies suggest that nrEVEs may interfere with virus replication by producing PIWI-interacting RNAs (piRNAs) that recognize and degrade viral RNAs through sequence complementarity. In this article, we identified the nrEVEs repertoire of ten species within the dipteran family Tephritidae (true fruit flies), which are considered a major threat to agriculture worldwide. Our results suggest that each of these species contains nrEVEs, although in limited numbers, and that nrEVE integration may have occurred both before and after speciation. Furthermore, the majority of nrEVEs originated from viruses with negative single-stranded RNA genomes and represent structural viral functions. Notably, these nrEVEs exhibit low similarity to currently known circulating viruses. To explore the potential role of nrEVEs, we investigated their transcription pattern and the production of piRNAs in different tissues of Ceratitis capitata. We successfully identified piRNAs that are complementary to the sequence of one nrEVE in C. capitata, thereby highlighting a potential link between nrEVEs and the piRNA pathway. Overall, our results provide valuable insights into the comparative landscape of nrEVEs in true fruit flies, contributing to the understanding of the intimate relation between fruit flies and their past and present viral pathogens.

Co-option of a non-retroviral endogenous viral element in planthoppers

Non-retroviral endogenous viral elements (nrEVEs) are widely dispersed throughout the genomes of eukaryotes. Although nrEVEs are known to be involved in host antiviral immunity, it remains an open question whether they can be domesticated as functional proteins to serve cellular innovations in arthropods. In this study, we found that endogenous toti-like viral elements (ToEVEs) are ubiquitously integrated into the genomes of three planthopper species, with highly variable distributions and polymorphism levels in planthopper populations. Three ToEVEs display exon‒intron structures and active transcription, suggesting that they might have been domesticated by planthoppers. CRISPR/Cas9 experiments revealed that one ToEVE in Nilaparvata lugens, NlToEVE14, has been co-opted by its host and plays essential roles in planthopper development and fecundity. Large-scale analysis of ToEVEs in arthropod genomes indicated that the number of arthropod nrEVEs is currently underestimated and that they may contribute to the functional diversity of arthropod genes.

Accumulated metagenomic studies reveal recent migration, whole genome evolution, and undiscovered diversity of orthomyxoviruses

IMPORTANCE

The number of known virus species has increased dramatically through metagenomic studies, which search genetic material sampled from a host for non-host genes. Here, we focus on an important viral family that includes influenza viruses, the Orthomyxoviridae, with over 100 recently discovered viruses infecting hosts from humans to fish. We find that one virus called Wǔhàn mosquito virus 6, discovered in mosquitoes in China, has spread across the globe very recently. Surface proteins used to enter cells show signs of rapid evolution in Wǔhàn mosquito virus 6 and its relatives which suggests an ability to infect vertebrate animals. We compute the rate at which new orthomyxovirus species discovered add evolutionary history to the tree of life, predict that many viruses remain to be discovered, and discuss what appropriately designed future studies can teach us about how diseases cross between continents and species.

Mosquito community composition shapes virus prevalence patterns along anthropogenic disturbance gradients

Previously unknown pathogens often emerge from primary ecosystems, but there is little knowledge on the mechanisms of emergence. Most studies analyzing the influence of land-use change on pathogen emergence focus on a single host-pathogen system and often observe contradictory effects. Here, we studied virus diversity and prevalence patterns in natural and disturbed ecosystems using a multi-host and multi-taxa approach. Mosquitoes sampled along a disturbance gradient in Côte d'Ivoire were tested by generic RT-PCR assays established for all major arbovirus and insect-specific virus taxa including novel viruses previously discovered in these samples based on cell culture isolates enabling an unbiased and comprehensive approach. The taxonomic composition of detected viruses was characterized and viral infection rates according to habitat and host were analyzed. We detected 331 viral sequences pertaining to 34 novel and 15 previously identified viruses of the families Flavi-, Rhabdo-, Reo-, Toga-, Mesoni- and Iflaviridae and the order Bunyavirales. Highest host and virus diversity was observed in pristine and intermediately disturbed habitats. The majority of the 49 viruses was detected with low prevalence. However, nine viruses were found frequently across different habitats of which five viruses increased in prevalence towards disturbed habitats, in congruence with the dilution effect hypothesis. These viruses were mainly associated with one specific mosquito species (Culex nebulosus), which increased in relative abundance from pristine (3%) to disturbed habitats (38%). Interestingly, the observed increased prevalence of these five viruses in disturbed habitats was not caused by higher host infection rates but by increased host abundance, an effect tentatively named abundance effect. Our data show that host species composition is critical for virus abundance. Environmental changes that lead to an uneven host community composition and to more individuals of a single species are a key driver of virus emergence.

Characterization of Viral Interference in Aedes albopictus C6/36 Cells Persistently Infected with Dengue Virus 2

Arboviruses are an important group of pathogens that cause diseases of medical and veterinary concern worldwide. The interactions of these viruses with their host cells are complex, and frequently, the coexistence of two different viruses in the same cell results in the inhibition of replication in one of the viruses, which is a phenomenon called viral interference. This phenomenon can be exploited to develop antiviral strategies. Insect cell lines persistently infected with arboviruses are useful models with which to study viral interference. In this work, a model of C6/36-HT cells (from Aedes albopictus mosquitoes) persistently infected with Dengue virus, serotype 2, was used. Viral interference was evaluated via plaque and flow cytometry assays. The presence of heterotypic interference against the other serotypes of the same virus and homologous interference against yellow fever virus was determined; however, this cell line did not display heterologous viral interference against Sindbis virus. The mechanisms responsible for viral interference have not been fully elucidated, but small RNAs could be involved. However, the silencing of Ago3, a key protein in the genome-derived P-element-induced wimpy testis pathway, did not alter the viral interference process, suggesting that viral interference occurs independent of this pathway.

Large-scale reference-free analysis of flavivirus sequences in Aedes aegypti whole genome DNA sequencing data

Flaviviruses are a diverse group of RNA viruses, which include the etiological agents of Zika, dengue and yellow fever that are transmitted by mosquitoes. Flaviviruses do not encode reverse transcriptase and cannot reverse transcribe into DNA, yet DNA sequences of flaviviruses are found both integrated in the chromosomes of Aedes aegypti mosquitoes and as extrachromosomal sequences. We have previously examined the Ae. aegypti reference genome to identify flavivirus integrations and analyzed conservation of these sequences among whole-genome data of 464 Ae. aegypti collected across 10 countries globally. Here, we extended this analysis by identifying flavivirus sequences in these samples independently of the Ae. aegypti reference assembly. Our aim was to identify the complete set of viral sequences, including those absent in the reference genome, and their geographical distribution. We compared the identified sequences using BLASTn and applied machine learning methods to identify clusters of similar sequences. Apart from clusters of sequences that correspond to the four viral integration events that we had previously described, we identified 19 smaller clusters. The only cluster with a strong geographic association consisted of Cell-fusing agent virus-like sequences specific to Thailand. The remaining clusters did not have a geographic association and mostly consisted of near identical short sequences without strong similarity to any known flaviviral genomes. The short read sequencing data did not permit us to determine whether identified sequences were extrachromosomal or integrated into Ae. aegypti chromosomes. Our results suggest that Liverpool strain and field Ae. aegypti mosquitoes have a similar variety of conserved flaviviral DNA, whose functional role should be investigated in follow-up studies.

Trans-generational immune priming is not mediated by the sex of the parent primed: a meta-analysis of invertebrate data

Traditionally, only vertebrates were thought capable of acquired immune responses, such as the ability to transfer immunological experience vertically to their offspring (known as trans-generational immune priming, TGIP). Increasing evidence challenges this belief and it is now clear that invertebrates also have the ability to exhibit functionally equivalent TGIP. This has led to a surge in papers exploring invertebrate TGIP, with most focusing on the costs, benefits or factors that affect the evolution of this trait. Whilst many studies have found support for the phenomenon, not all studies do, and there is considerable variation in the strength of positive results. To address this, we conducted a meta-analysis to answer the question: what is the overall effect of TGIP in invertebrates? Then, to understand the specific factors that affect its presence and intensity, we conducted a moderator analysis. Our results corroborate that TGIP occurs in invertebrates (demonstrated by a large, positive effect size). The strength of the positive effect was related to if and how offspring were immune challenged (i.e. whether they were challenged with the same or different insult as their parents or not challenged at all). Interestingly, there was no effect of the ecology or life history of the species or the sex of the parent or the offspring primed, and responses were comparable across different immune elicitors. Our publication bias testing suggests that the literature may suffer from some level of positive-result bias. However, even after accounting for potential bias, our effect size remains positive. Publication bias testing can be influenced by diversity in the data set, which was considerable in our data, even after moderator analysis. It is therefore conceivable that differences among studies could be caused by other moderators that were unable to be included in our meta-analysis. Nonetheless, our results suggest that TGIP does occur in invertebrates, whilst providing some potential avenues to examine the factors that account for variation in effect sizes.

Comparative evolutionary analyses of eight whitefly Bemisia tabaci sensu lato genomes: cryptic species, agricultural pests and plant-virus vectors

BACKGROUND

The group of > 40 cryptic whitefly species called Bemisia tabaci sensu lato are amongst the world's worst agricultural pests and plant-virus vectors. Outbreaks of B. tabaci s.l. and the associated plant-virus diseases continue to contribute to global food insecurity and social instability, particularly in sub-Saharan Africa and Asia. Published B. tabaci s.l. genomes have limited use for studying African cassava B. tabaci SSA1 species, due to the high genetic divergences between them. Genomic annotations presented here were performed using the 'Ensembl gene annotation system', to ensure that comparative analyses and conclusions reflect biological differences, as opposed to arising from different methodologies underpinning transcript model identification.

RESULTS

We present here six new B. tabaci s.l. genomes from Africa and Asia, and two re-annotated previously published genomes, to provide evolutionary insights into these globally distributed pests. Genome sizes ranged between 616-658 Mb and exhibited some of the highest coverage of transposable elements reported within Arthropoda. Many fewer total protein coding genes (PCG) were recovered compared to the previously published B. tabaci s.l. genomes and structural annotations generated via the uniform methodology strongly supported a repertoire of between 12.8-13.2 × 10³ PCG. An integrative systematics approach incorporating phylogenomic analysis of nuclear and mitochondrial markers supported a monophyletic Aleyrodidae and the basal positioning of B. tabaci Uganda-1 to the sub-Saharan group of species. Reciprocal cross-mating data and the co-cladogenesis pattern of the primary obligate endosymbiont 'Candidatus Portiera aleyrodidarum' from 11 Bemisia genomes further supported the phylogenetic reconstruction to show that African cassava B. tabaci populations consist of just three biological species. We include comparative analyses of gene families related to detoxification, sugar metabolism, vector competency and evaluate the presence and function of horizontally transferred genes, essential for understanding the evolution and unique biology of constituent B. tabaci. s.l species.

CONCLUSIONS

These genomic resources have provided new and critical insights into the genetics underlying B. tabaci s.l. biology. They also provide a rich foundation for post-genomic research, including the selection of candidate gene-targets for innovative whitefly and virus-control strategies.

Endoparasitoid lifestyle promotes endogenization and domestication of dsDNA viruses

The accidental endogenization of viral elements within eukaryotic genomes can occasionally provide significant evolutionary benefits, giving rise to their long-term retention, that is, to viral domestication. For instance, in some endoparasitoid wasps (whose immature stages develop inside their hosts), the membrane-fusion property of double-stranded DNA viruses have been repeatedly domesticated following ancestral endogenizations. The endogenized genes provide female wasps with a delivery tool to inject virulence factors that are essential to the developmental success of their offspring. Because all known cases of viral domestication involve endoparasitic wasps, we hypothesized that this lifestyle, relying on a close interaction between individuals, may have promoted the endogenization and domestication of viruses. By analyzing the composition of 124 Hymenoptera genomes, spread over the diversity of this clade and including free-living, ecto, and endoparasitoid species, we tested this hypothesis. Our analysis first revealed that double-stranded DNA viruses, in comparison with other viral genomic structures (ssDNA, dsRNA, ssRNA), are more often endogenized and domesticated (that is, retained by selection) than expected from their estimated abundance in insect viral communities. Second, our analysis indicates that the rate at which dsDNA viruses are endogenized is higher in endoparasitoids than in ectoparasitoids or free-living hymenopterans, which also translates into more frequent events of domestication. Hence, these results are consistent with the hypothesis that the endoparasitoid lifestyle has facilitated the endogenization of dsDNA viruses, in turn, increasing the opportunities of domestications that now play a central role in the biology of many endoparasitoid lineages.

Endogenous viral elements reveal associations between a non-retroviral RNA virus and symbiotic dinoflagellate genomes

Endogenous viral elements (EVEs) offer insight into the evolutionary histories and hosts of contemporary viruses. This study leveraged DNA metagenomics and genomics to detect and infer the host of a non-retroviral dinoflagellate-infecting +ssRNA virus (dinoRNAV) common in coral reefs. As part of the Tara Pacific Expedition, this study surveyed 269 newly sequenced cnidarians and their resident symbiotic dinoflagellates (Symbiodiniaceae), associated metabarcodes, and publicly available metagenomes, revealing 178 dinoRNAV EVEs, predominantly among hydrocoral-dinoflagellate metagenomes. Putative associations between Symbiodiniaceae and dinoRNAV EVEs were corroborated by the characterization of dinoRNAV-like sequences in 17 of 18 scaffold-scale and one chromosome-scale dinoflagellate genome assembly, flanked by characteristically cellular sequences and in proximity to retroelements, suggesting potential mechanisms of integration. EVEs were not detected in dinoflagellate-free (aposymbiotic) cnidarian genome assemblies, including stony corals, hydrocorals, jellyfish, or seawater. The pervasive nature of dinoRNAV EVEs within dinoflagellate genomes (especially Symbiodinium), as well as their inconsistent within-genome distribution and fragmented nature, suggest ancestral or recurrent integration of this virus with variable conservation. Broadly, these findings illustrate how +ssRNA viruses may obscure their genomes as members of nested symbioses, with implications for host evolution, exaptation, and immunity in the context of reef health and disease.

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.

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.

Composition and global distribution of the mosquito virome - A comprehensive database of insect-specific viruses

Mosquitoes are vectors for emerging and re-emerging infectious viral diseases of humans, livestock and other animals. In addition to these arthropod-borne (arbo)viruses, mosquitoes are host to an array of insect-specific viruses, collectively referred to as the mosquito virome. Mapping the mosquito virome and understanding if and how its composition modulates arbovirus transmission is critical to understand arboviral disease emergence and outbreak dynamics. In recent years, next-generation sequencing as well as PCR and culture-based methods have been extensively used to identify mosquito-associated viruses, providing insights into virus ecology and evolution. Until now, the large amount of mosquito virome data, specifically those acquired by metagenomic sequencing, has not been comprehensively integrated. We have constructed a searchable database of insect-specific viruses associated with vector mosquitoes from 175 studies, published between October 2000 and February 2022. We identify the most frequently detected and widespread viruses of the Culex, Aedes and Anopheles mosquito genera and report their global distribution. In addition, we highlight the challenges of extracting and integrating published virome data and we propose that a standardized reporting format will facilitate data interpretation and re-use by other scientists. We expect our comprehensive database, summarizing mosquito virome data collected over 20 years, to be a useful resource for future studies.

Mosquito vector competence for dengue is modulated by insect-specific viruses

Aedes aegypti and A. albopictus mosquitoes are the main vectors for dengue virus (DENV) and other arboviruses, including Zika virus (ZIKV). Understanding the factors that affect transmission of arboviruses from mosquitoes to humans is a priority because it could inform public health and targeted interventions. Reasoning that interactions among viruses in the vector insect might affect transmission, we analysed the viromes of 815 urban Aedes mosquitoes collected from 12 countries worldwide. Two mosquito-specific viruses, Phasi Charoen-like virus (PCLV) and Humaita Tubiacanga virus (HTV), were the most abundant in A. aegypti worldwide. Spatiotemporal analyses of virus circulation in an endemic urban area revealed a 200% increase in chances of having DENV in wild A. aegypti mosquitoes when both HTV and PCLV were present. Using a mouse model in the laboratory, we showed that the presence of HTV and PCLV increased the ability of mosquitoes to transmit DENV and ZIKV to a vertebrate host. By transcriptomic analysis, we found that in DENV-infected mosquitoes, HTV and PCLV block the downregulation of histone H4, which we identify as an important proviral host factor in vivo.

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.

Understanding pathogen survival and transmission by arthropod vectors to prevent human disease

Many endemic poverty-associated diseases, such as malaria and leishmaniasis, are transmitted by arthropod vectors. Pathogens must interact with specific molecules in the vector gut, the microbiota, and the vector immune system to survive and be transmitted. The vertebrate host, in turn, is infected when the pathogen and vector-derived factors, such as salivary proteins, are delivered into the skin by a vector bite. Here, we review recent progress in our understanding of the biology of pathogen transmission from the human to the vector and back, from the vector to the host. We also highlight recent advances in the biology of vector-borne disease transmission, which have translated into additional strategies to prevent human disease by either reducing vector populations or by disrupting their ability to transmit pathogens.

Comparative analysis of genome-encoded viral sequences reveals the evolutionary history of flavivirids (family Flaviviridae)

Flavivirids (family Flaviviridae) are a group of positive-strand ribonucleic acid (RNA) viruses that pose serious risks to human and animal health on a global scale. Here, we use flavivirid-derived deoxyribonucleic acid (DNA) sequences, identified in animal genomes, to reconstruct the long-term evolutionary history of family Flaviviridae. We demonstrate that flavivirids are >100 million years old and show that this timing can be combined with dates inferred from co-phyletic analysis to produce a cohesive overview of their evolution, distribution, and diversity wherein the main flavivirid subgroups originate in early animals and broadly co-diverge with major animal phyla. In addition, we reveal evidence that the 'classical flaviviruses' of vertebrates, most of which are transmitted via blood-feeding arthropod vectors, originally evolved in haematophagous arachnids and later acquired the capacity to be transmitted by insects. Our findings imply that the biological properties of flavivirids have been acquired gradually over the course of animal evolution. Thus, broad-scale comparative analysis will likely reveal fundamental insights into their biology. We therefore published our results via an open, extensible, database (Flavivirid-GLUE), which we constructed to facilitate the wider utilisation of genomic data and evolution-related domain knowledge in flavivirid research.

Integrated Jingmenvirus Polymerase Gene in Ixodes ricinus Genome

Members of the jingmenviruses group have been found in arthropods and mammals on all continents except Australia and Antarctica. Two viruses of this group were isolated from patients with fever after a tick bite. Using a nested RT-PCR assay targeting a jingmenvirus polymerase gene fragment, we screened ticks collected in seven regions of Russia and found that the abundant jingmenvirus-positive were of Ixodes ricinus species, with the prevalence ranging from 19.8% to 34.3%. In all cases, DNase/RNase treatment suggested that the detected molecule was DNA and subsequent next generation sequencing (NGS) proved that the viral polymerase gene was integrated in the I. ricinus genome. The copy number of the integrated polymerase gene was quantified by qPCR relative to the ITS2 gene and estimated as 1.32 copies per cell. At least three different genetic variants of the integrated polymerase gene were found in the territory of Russia. Phylogenetic analysis of the integrated jingmenvirus polymerase gene showed the highest similarity with the sequence of the correspondent gene obtained in Serbia from I. ricinus.

Shrimp genome sequence contains independent clusters of ancient and current Endogenous Viral Elements (EVE) of the parvovirus IHHNV

Background

Shrimp have the ability to accommodate viruses in long term, persistent infections without signs of disease. Endogenous viral elements (EVE) play a role in this process probably via production of negative-sense Piwi-interacting RNA (piRNA)-like fragments. These bind with Piwi proteins to dampen viral replication via the RNA interference (RNAi) pathway. We searched a genome sequence (GenBank record JABERT000000000) of the giant tiger shrimp (Penaeus monodon for the presence of EVE related to a shrimp parvovirus originally named infectious hypodermal and hematopoietic necrosis virus (IHHNV).

Results

The shrimp genome sequence contained three piRNA-like gene clusters containing scrambled IHHNV EVE. Two clusters were located distant from one another in pseudochromosome 35 (PC35). Both PC35 clusters contained multiple sequences with high homology (99%) to GenBank records DQ228358 and EU675312 that were both called “non-infectious IHHNV Type A” (IHHNV-A) when originally discovered. However, our results and those from a recent Australian P. monodon genome assembly indicate that the relevant GenBank records for IHHNV-A are sequence-assembly artifacts derived from scrambled and fragmental IHHNV-EVE. Although the EVE in the two PC35 clusters showed high homology only to IHHNV-A, the clusters were separate and distinct with respect to the arrangement (i.e., order and reading direction) and proportional content of the IHHNV-A GenBank records. We conjecture that these 2 clusters may constitute independent allele-like clusters on a pair of homologous chromosomes. The third EVE cluster was found in pseudochromosome 7 (PC7). It contained EVE with high homology (99%) only to GenBank record AF218266 with the potential to protect shrimp against current types of infectious IHHNV. One disadvantage was that some EVE in PC7 can give false positive PCR test results for infectious IHHNV.

Conclusions

Our results suggested the possibility of viral-type specificity in EVE clusters. Specificity is important because whole EVE clusters for one viral type would be transmitted to offspring as collective hereditary units. This would be advantageous if one or more of the EVE within the cluster were protective against the disease caused by the cognate virus. It would also facilitate gene editing for removal of non-protective EVE clusters or for transfer of protective EVE clusters to genetically improve existing shrimp breeding stocks that might lack them.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08802-3.

A gene-based capture assay for surveying patterns of genetic diversity and insecticide resistance in a worldwide group of invasive mosquitoes

Understanding patterns of diversification, genetic exchange, and pesticide resistance in arthropod disease vectors is necessary for effective population management. With the availability of next-generation sequencing technologies, one of the best approaches for surveying such patterns involves the simultaneous genotyping of many samples for a large number of genetic markers. To this end, the targeting of gene sequences of known function can be a cost-effective strategy. One insect group of substantial health concern are the mosquito taxa that make up the Culex pipiens complex. Members of this complex transmit damaging arboviruses and filariae worms to humans, as well as other pathogens such as avian malaria parasites that are detrimental to birds. Here we describe the development of a targeted, gene-based assay for surveying genetic diversity and population structure in this mosquito complex. To test the utility of this assay, we sequenced samples from several members of the complex, as well as from distinct populations of the relatively under-studied Culex quinquefasciatus. The data generated was then used to examine taxonomic divergence and population clustering between and within these mosquitoes. We also used this data to investigate genetic variants present in our samples that had previously been shown to correlate with insecticide-resistance. Broadly, our gene capture approach successfully enriched the genomic regions of interest, and proved effective for facilitating examinations of taxonomic divergence and geographic clustering within the Cx. pipiens complex. It also allowed us to successfully survey genetic variation associated with insecticide resistance in Culex mosquitoes. This enrichment protocol will be useful for future studies that aim to understand the genetic mechanisms underlying the evolution of these ubiquitous and increasingly damaging disease vectors.

The mosquito taxa that make up the Culex pipiens complex are important vectors of the agents of several human diseases such as West Nile and St. Louis encephalitides, and lymphatic filariasis. They are also important vectors of avian malaria, which impacts livestock and wildlife. The development of effective strategies for the control of these mosquitoes requires knowledge of their origins, distribution, dispersal patterns, and the extent to which discreet taxonomic entities within the complex interbreed. To achieve these objectives, it is necessary to compare patterns of genetic diversity across many mosquito samples, which can be cost-prohibitive. To address this limitation, we developed a targeted, gene-based assay that allowed us to cost-effectively genotype a large number of genetic variants from a representative global sampling of individual Cx. pipiens complex mosquitoes. We show that this assay is a powerful tool for examining genetic structure and hybridization among populations. We also explore its utility for surveying alleles previously shown to be associated with insecticide resistance. Future use of this enrichment assay and the bioinformatics methods described here will allow researchers to study evolutionary patterns across the Cx. pipiens complex as well as monitor the presence of genetic variation that could affect control efforts.

Transposable element repression using piRNAs, and its relevance to endogenous viral elements (EVEs) and immunity in insects

The piRNA system controls transposable element (TE) mobility by transcriptional gene silencing and post-transcriptional gene silencing. Dispersed in insect genomes, piRNA clusters contain TE copies, from which they produce piRNAs (specific small RNAs). These piRNAs can both target the nascent transcripts produced by active TE copies and directly repress them by heterochromatinization. They can also target mature transcripts and cleave them following amplification by the so-called 'ping-pong' loop mechanism. Moreover, piRNA clusters contain endogenous viral elements (EVEs), from which they produce piRNAs. The current idea is that these piRNAs could participate in the antiviral response against exogenous viral infection. In this review, we show that among insects, to date, this antiviral response by the piRNA system appears mainly restricted to mosquitoes, but this could be due to the focus of most studies on arboviruses.

Jingchuvirales: a New Taxonomical Framework for a Rapidly Expanding Order of Unusual Monjiviricete Viruses Broadly Distributed among Arthropod Subphyla

Technical advances in metagenomics and metatranscriptomics have dramatically accelerated virus discovery in recent years. "Chuviruses" were first described in 2015 as obscure negative-sense RNA viruses of diverse arthropods. Although "chuviruses" first appeared to be members of the negarnaviricot order Mononegavirales in phylogenetic analyses using RNA-directed RNA polymerase sequences, further characterization revealed unusual gene orders in genomes that are nonsegmented, segmented, and/or possibly circular. Consequently, a separate order, Jingchuvirales, was established to include a monospecific family, Chuviridae. Recently, it has become apparent that jingchuvirals are broadly distributed and are therefore likely of ecological and economic importance. Here, we describe recent and ongoing efforts to create the necessary taxonomic framework to accommodate the expected flood of novel viruses belonging to the order.

Intrinsic variation in the vertically transmitted core virome of the mosquito Aedes aegypti

Virome studies among metazoans have revealed the ubiquity of RNA viruses in animals, contributing to a fundamental re-thinking of the relationships between organisms and their microbiota. Mosquito viromes, often scrutinized due to their public health relevance, may also provide insight into broadly applicable concepts, such as a "core virome," a set of viruses consistently associated with a host species or population that may fundamentally impact its basic biology. A subset of mosquito-associated viruses (MAVs) could comprise such a core, and MAVs can be categorized as (i) arboviruses, which alternate between mosquito and vertebrate hosts, (ii) insect-specific viruses, which cannot replicate in vertebrate cells, and (iii) viruses with unknown specificity. MAVs have been widely characterized in the disease vector Aedes aegypti, and the occurrence of a core virome in this species has been proposed but remains unclear. Using a wild population previously surveyed for MAVs and a common laboratory strain, we investigated viromes in reproductive tissue via metagenomic RNA sequencing. Virome composition varied across samples, but four groups comprised >97% of virus sequences: a novel partiti-like virus (Partitiviridae), a toti-like virus (Totiviridae), unclassified Riboviria, and four orthomyxo-like viruses (Orthormyxoviridae). Whole or partial genomes for the partiti-like virus, toti-like virus, and one orthomyxo-like virus were assembled and analyzed phylogenetically. Multigenerational maintenance of these MAVs was confirmed by RT-PCR, indicating vertical transmission as a mechanism for persistence. This study provides fundamental information regarding MAV ecology and variability in A. aegypti and the potential for vertically maintained core viromes at the population level.

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.

Innate immune memory in invertebrates: Concept and potential mechanisms

Invertebrates are the protagonists of a recent paradigm shift because they now show that vertebrates are not the only group with immune memory. This review discusses the concept of immune priming, its characteristics, and differences with trained immunity and immune enhancement. We include an update of the current status of immune priming within generations in different groups of invertebrates which now include work in 5 Phyla: Ctenophora, Cnidaria, Mollusca, Nematoda, and Arthropoda. Clearly, few Phyla have been studied. We also resume and discuss the effector mechanism related to immune memory, including integrating viral elements into the genome, endoreplication, and epigenetics. The roles of other elements are incorporated, such as hemocytes, immune pathways, and metabolisms. We conclude that taking care of the experimental procedure will discern if results provide or do not support the invertebrates' immune memory and that regarding mechanisms, indeed, there are no studies on the immune memory mechanisms, this is how specificity is reached, and how and where the immune memory is stored and how is recall upon subsequent encounters. Finally, we discuss the possibility of having more than one mechanism working in different groups of invertebrates depending on the environmental conditions.

RNA virus EVEs in insect genomes

Insects are infected by a diverse set of RNA viruses that are more broadly distinguished by their ability to infect single or multiple host species. During replication into the host cell, partial or complete double strand DNA derived from the viral genome may be integrated into their host genomes giving origin to endogenous viral elements (EVEs). EVEs from RNA viruses have been identified in a variety of insect genomes showing different evolutionary trajectories: from highly degraded viral genomic remains to partial and complete viral coding regions. Limited functional knowledge exists about RNA EVEs impact on hosts and circulating viruses, but exciting results are emerging showing a complex arms race interplay that influences the evolutionary trajectory of these interacting entities.

Virus and endogenous viral element-derived small non-coding RNAs and their roles in insect-virus interaction

RNA interference pathways mediated by different types of small non-coding RNAs (siRNAs, miRNAs and piRNAs) are conserved biological responses to exotic stresses, including viral infection. Aside from the well-established siRNA pathway, the miRNA pathway and the piRNA pathway process viral sequences, exogenously or endogenously, into miRNAs and piRNAs, respectively. During the host-virus interaction, viral sequences, including both coding and non-coding sequences, can be integrated as endogenous viral elements (EVEs) and thereby become present within the germline of a non-viral organism. In recent years, significant progress has been made in characterizing the biogenesis and function of viruses and EVEs associated with snRNAs. Overall, the siRNA pathway acts as the primarily antiviral defense against a wide range of exogenous viruses; the miRNA pathways associated with viruses or EVEs function in antiviral response and host gene regulation; EVE derived piRNAs with a ping-pong signature have the potential to limit cognate viral infection.

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.

Somatic piRNAs and Transposons are Differentially Expressed Coincident with Skeletal Muscle Atrophy and Programmed Cell Death

PIWI-interacting RNAs (piRNAs) are small single-stranded RNAs that can repress transposon expression via epigenetic silencing and transcript degradation. They have been identified predominantly in the ovary and testis, where they serve essential roles in transposon silencing in order to protect the integrity of the genome in the germline. The potential expression of piRNAs in somatic cells has been controversial. In the present study we demonstrate the expression of piRNAs derived from both genic and transposon RNAs in the intersegmental muscles (ISMs) from the tobacco hawkmoth Manduca sexta. These piRNAs are abundantly expressed, ∼27 nt long, map antisense to transposons, are oxidation resistant, exhibit a 5’ uridine bias, and amplify via the canonical ping-pong pathway. An RNA-seq analysis demonstrated that 19 piRNA pathway genes are expressed in the ISMs and are developmentally regulated. The abundance of piRNAs does not change when the muscles initiate developmentally-regulated atrophy, but are repressed coincident with the commitment of the muscles undergo programmed cell death at the end of metamorphosis. This change in piRNA expression is correlated with the repression of several retrotransposons and the induction of specific DNA transposons. The developmentally-regulated changes in the expression of piRNAs, piRNA pathway genes, and transposons are all regulated by 20-hydroxyecdysone, the steroid hormone that controls the timing of ISM death. Taken together, these data provide compelling evidence for the existence of piRNA in somatic tissues and suggest that they may play roles in developmental processes such as programmed cell death.

Computational Methods for the Discovery and Annotation of Viral Integrations

The transfer of genetic material between viruses and eukaryotic cells is pervasive. Somatic integrations of DNA viruses and retroviruses have been linked to persistent viral infection and genotoxic effects. Integrations into germline cells, referred to as Endogenous Viral Elements (EVEs), can be co-opted for host functions. Besides DNA viruses and retroviruses, EVEs can also derive from nonretroviral RNA viruses, which have often been observed in piRNA clusters. Here, we describe a bioinformatic framework to annotate EVEs in a genome assembly, study their widespread occurrence and polymorphism and identify sample-specific viral integrations using whole genome sequencing data.

Detection of Insect-Specific Flaviviruses in Mosquitoes (Diptera: Culicidae) in Northeastern Regions of South Africa

Mosquitoes in the Aedes and Culex genera are considered the main vectors of pathogenic flaviviruses worldwide. Entomological surveillance using universal flavivirus sets of primers in mosquitoes can detect not only pathogenic viruses but also insect-specific ones. It is hypothesized that insect-specific flaviviruses, which naturally infect these mosquitoes, may influence their vector competence for zoonotic arboviruses. Here, entomological surveillance was performed between January 2014 and May 2018 in five different provinces in the northeastern parts of South Africa, with the aim of identifying circulating flaviviruses. Mosquitoes were sampled using different carbon dioxide trap types. Overall, 64,603 adult mosquitoes were collected, which were screened by RT-PCR and sequencing. In total, 17 pools were found positive for insect-specific Flaviviruses in the mosquito genera Aedes (12/17, 70.59%) and Anopheles (5/17, 29.41%). No insect-specific viruses were detected in Culex species. Cell-fusing agent viruses were detected in Aedes aegypti and Aedes caballus. A range of anopheline mosquitoes, including Anopheles coustani, An. squamosus and An. maculipalpis, were positive for Culex flavivirus-like and Anopheles flaviviruses. These results confirm the presence of insect-specific flaviviruses in mosquito populations in South Africa, expands their geographical range and indicates potential mosquito species as vector species.