Metagenomic sequencing suggests a diversity of RNA interference-like responses to viruses across multicellular eukaryotes

The protein structurome of Orthornavirae and its dark matter

Metatranscriptomics is uncovering more and more diverse families of viruses with RNA genomes comprising the viral kingdom Orthornavirae in the realm Riboviria. Thorough protein annotation and comparison are essential to get insights into the functions of viral proteins and virus evolution. In addition to sequence- and hmm profile‑based methods, protein structure comparison adds a powerful tool to uncover protein functions and relationships. We constructed an Orthornavirae "structurome" consisting of already annotated as well as unannotated ("dark matter") proteins and domains encoded in viral genomes. We used protein structure modeling and similarity searches to illuminate the remaining dark matter in hundreds of thousands of orthornavirus genomes. The vast majority of the dark matter domains showed either "generic" folds, such as single α-helices, or no high confidence structure predictions. Nevertheless, a variety of lineage-specific globular domains that were new either to orthornaviruses in general or to particular virus families were identified within the proteomic dark matter of orthornaviruses, including several predicted nucleic acid-binding domains and nucleases. In addition, we identified a case of exaptation of a cellular nucleoside monophosphate kinase as an RNA-binding protein in several virus families. Notwithstanding the continuing discovery of numerous orthornaviruses, it appears that all the protein domains conserved in large groups of viruses have already been identified. The rest of the viral proteome seems to be dominated by poorly structured domains including intrinsically disordered ones that likely mediate specific virus-host interactions.

IMPORTANCE

Advanced methods for protein structure prediction, such as AlphaFold2, greatly expand our capability to identify protein domains and infer their likely functions and evolutionary relationships. This is particularly pertinent for proteins encoded by viruses that are known to evolve rapidly and as a result often cannot be adequately characterized by analysis of the protein sequences. We performed an exhaustive structure prediction and comparative analysis for uncharacterized proteins and domains ("dark matter") encoded by viruses with RNA genomes. The results show the dark matter of RNA virus proteome consists mostly of disordered and all-α-helical domains that cannot be readily assigned a specific function and that likely mediate various interactions between viral proteins and between viral and host proteins. The great majority of globular proteins and domains of RNA viruses are already known although we identified several unexpected domains represented in individual viral families.

An evolutionary perspective to innate antiviral immunity in animals

Viruses pose a significant threat to cellular organisms. Innate antiviral immunity encompasses both RNA- and protein-based mechanisms designed to sense and respond to infections, a fundamental aspect present in all living organisms. A potent RNA-based antiviral mechanism is RNA interference, where small RNA-programmed nucleases target viral RNAs. Protein-based mechanisms often rely on the induction of transcriptional responses triggered by the recognition of viral infections through innate immune receptors. These responses involve the upregulation of antiviral genes aimed at countering viral infections. In this review, we delve into recent advances in understanding the diversification of innate antiviral immunity in animals. An evolutionary perspective on the gains and losses of mechanisms in diverse animals coupled to mechanistic studies in model organisms such as the fruit fly Drosophila melanogaster is essential to provide deep understanding of antiviral immunity that can be translated to new strategies in the treatment of viral diseases.

Insights into ADAR gene complement, expression patterns, and RNA editing landscape in Chlamys farreri

Adenosine Deaminases Acting on RNA (ADARs) are evolutionarily conserved enzymes known to convert adenosine to inosine in double-stranded RNAs and participate in host-virus interactions. Conducting a meta-analysis of available transcriptome data, we identified and characterised eight ADAR transcripts in Chlamys farreri, a farmed marine scallop susceptible to Acute viral necrosis virus (AVNV) infections and mortality outbreaks. Accordingly, we identified six ADAR genes in the Zhikong scallop genome, revised previous gene annotations, and traced alternative splicing variants. In detail, each ADAR gene encodes a unique combination of functional domains, always including the Adenosine deaminase domain, RNA binding domains and, in one case, two copies of a Z-DNA binding domain. After phylogenetic analysis, five C. farreri ADARs clustered in the ADAR1 clade along with sequences from diverse animal phyla. Gene expression analysis indicated CF051320 as the most expressed ADAR, especially in the eye and male gonad. The other four ADAR1 genes and one ADAR2 gene exhibited variable expression levels, with CF105370 and CF051320 significantly increasing during early scallop development. ADAR-mediated single-base editing, evaluated across adult C. farreri tissues and developmental stages, was mainly detectable in intergenic regions (83 % and 85 %, respectively). Overall, the expression patterns of the six ADAR genes together with the editing and hyper-editing values computed on scallops RNA-seq samples support the adaptive value of ADAR1-mediated editing, particularly in the pre-settling larval stages.

Discovery of Novel Viruses in Culicoides Biting Midges in Chihuahua, Mexico

Biting midges (Culicoides) are vectors of many pathogens of medical and veterinary importance, but their viromes are poorly characterized compared to certain other hematophagous arthropods, e.g., mosquitoes and ticks. The goal of this study was to use metagenomics to identify viruses in Culicoides from Mexico. A total of 457 adult midges were collected in Chihuahua, northern Mexico, in 2020 and 2021, and all were identified as female Culicoides reevesi. The midges were sorted into five pools and homogenized. An aliquot of each homogenate was subjected to polyethylene glycol precipitation to enrich for virions, then total RNA was extracted and analyzed by unbiased high-throughput sequencing. We identified six novel viruses that are characteristic of viruses from five families (Nodaviridae, Partitiviridae, Solemoviridae, Tombusviridae, and Totiviridae) and one novel virus that is too divergent from all classified viruses to be assigned to an established family. The newly discovered viruses are phylogenetically distinct from their closest known relatives, and their minimal infection rates in female C. reevesi range from 0.22 to 1.09. No previously known viruses were detected, presumably because viral metagenomics had never before been used to study Culicoides from the Western Hemisphere. To conclude, we discovered multiple novel viruses in C. reevesi from Mexico, expanding our knowledge of arthropod viral diversity and evolution.

Transcriptome Analysis Reveals a Diverse Range of Novel Viruses in Australian Sugarcane Soldier Fly (Inopus flavus) Larvae

In Australia, Soldier flies (Inopus spp.) are economically significant pests of sugarcane that currently lack a viable management strategy. Despite various research efforts, the mechanisms underlying the damage caused by soldier fly larvae remain poorly understood. Our study aims to explore whether this damage is associated with the transmission of plant viruses during larval feeding. We also explore the larval transcriptome to identify any entomopathogenic viruses with the potential to be used as biocontrol agents in future pest management programs. Seven novel virus sequences are identified and characterised using de novo assembly of RNA-Seq data obtained from salivary glands of larvae. The novel virus sequences belong to different virus families and are tentatively named SF-associated anphevirus (SFaAV), SF-associated orthomyxo-like virus (SFaOV), SF-associated narna-like virus (SFaNV), SF-associated partiti-like virus (SFaPV), SF-associated toti-like virus (SFaTV-1 and SFaTV-2) and SF-associated densovirus (SFaDV). These newly identified viruses are more likely insect-associated viruses, as phylogenetic analyses show that they cluster with other insect-specific viruses. Small RNA analysis indicates prominent peaks at both 21 nt and 26-29 nt, suggesting the activation of host siRNA and piwiRNA pathways. Our study helps to improve understanding of the virome of soldier flies and could identify insect viruses for deployment in novel pest management strategies.

Transcriptional landscape of small non-coding RNAs reveals diversity of categories and functions in molluscs

Small non-coding RNAs (sncRNAs) are non-coding RNA molecules that play various roles in metazoans. Among the sncRNAs, microRNAs (miRNAs) guide post-translational gene regulation during cellular development, proliferation, apoptosis, and differentiation, while PIWI-interacting RNAs (piRNAs) suppress transposon activity to safeguard the genome from detrimental insertion mutagenesis. While an increasing number of piRNAs are being identified in the soma and germlines of various organisms, they are scarcely reported in molluscs. To unravel the small RNA (sRNA) expression patterns and genomic function in molluscs, we generated a comprehensive sRNA dataset by sRNA sequencing (sRNA-seq) of eight mollusc species. Abundant miRNAs were identified and characterized in all investigated molluscs, and ubiquitous piRNAs were discovered in both somatic and gonadal tissues in six of the investigated molluscs, which are more closely associated with transposon silencing. Tens of piRNA clusters were also identified based on the genomic mapping results, which varied among different tissues and species. Our dataset serves as important reference data for future genomic and genetic studies on sRNAs in these molluscs and related species, especially in elucidating the ancestral state of piRNAs in bilaterians.

Virus replication in the honey bee parasite, Varroa destructor

IMPORTANCE

The parasitic mite Varroa destructor is a significant driver of worldwide colony losses of our most important commercial pollinator, the Western honey bee Apis mellifera. Declines in honey bee health are frequently attributed to the viruses that mites vector to honey bees, yet whether mites passively transmit viruses as a mechanical vector or actively participate in viral amplification and facilitate replication of honey bee viruses is debated. Our work investigating the antiviral RNA interference response in V. destructor demonstrates that key viruses associated with honey bee declines actively replicate in mites, indicating that they are biological vectors, and the host range of bee-associated viruses extends to their parasites, which could impact virus evolution, pathogenicity, and spread.

Viral Diversity in Samples of Freshwater Gastropods Benedictia baicalensis (Caenogastropoda: Benedictiidae) Revealed by Total RNA-Sequencing

Previously, the main studies were focused on viruses that cause disease in commercial and farmed shellfish and cause damage to food enterprises (for example, Ostreavirusostreidmalaco1, Aurivirus haliotidmalaco1 and Aquabirnavirus tellinae). Advances in high-throughput sequencing technologies have extended the studies to natural populations of mollusks (and other invertebrates) as unexplored niches of viral diversity and possible sources of emerging diseases. These studies have revealed a huge diversity of mostly previously unknown viruses and filled gaps in the evolutionary history of viruses. In the present study, we estimated the viral diversity in samples of the Baikal endemic gastropod Benedictia baicalensis using metatranscriptomic analysis (total RNA-sequencing); we were able to identify a wide variety of RNA-containing viruses in four samples (pools) of mollusks collected at three stations of Lake Baikal. Most of the identified viral genomes (scaffolds) had only distant similarities to known viruses or (in most cases) to metagenome-assembled viral genomes from various natural samples (mollusks, crustaceans, insects and others) mainly from freshwater ecosystems. We were able to identify viruses similar to those previously identified in mollusks (in particular to the picornaviruses Biomphalaria virus 1 and Biomphalaria virus 3 from the freshwater gastropods); it is possible that picorna-like viruses (as well as a number of other identified viruses) are pathogenic for Baikal gastropods. Our results also suggested that Baikal mollusks, like other species, may bioaccumulate or serve as a reservoir for numerous viruses that infect a variety of organisms (including vertebrates).

The structure and mechanism of action of a distinct class of dicistrovirus intergenic region IRESs

Internal ribosomal entry sites (IRESs) engage with the eukaryotic translation apparatus to promote end-independent initiation. We identified a conserved class of ∼150 nt long intergenic region (IGR) IRESs in dicistrovirus genomes derived from members of the phyla Arthropoda, Bryozoa, Cnidaria, Echinodermata, Entoprocta, Mollusca and Porifera. These IRESs, exemplified by Wenling picorna-like virus 2, resemble the canonical cricket paralysis virus (CrPV) IGR IRES in comprising two nested pseudoknots (PKII/PKIII) and a 3'-terminal pseudoknot (PKI) that mimics a tRNA anticodon stem-loop base-paired to mRNA. However, they are ∼50 nt shorter than CrPV-like IRESs, and PKIII is an H-type pseudoknot that lacks the SLIV and SLV stem-loops that are primarily responsible for the affinity of CrPV-like IRESs for the 40S ribosomal subunit and that restrict initial binding of PKI to its aminoacyl (A) site. Wenling-class IRESs bound strongly to 80S ribosomes but only weakly to 40S subunits. Whereas CrPV-like IRESs must be translocated from the A site to the peptidyl (P) site by elongation factor 2 for elongation to commence, Wenling-class IRESs bound directly to the P site of 80S ribosomes, and decoding begins without a prior translocation step. A chimeric CrPV clone containing a Wenling-class IRES was infectious, confirming that the IRES functioned in cells.

Initiation of translation on nedicistrovirus and related intergenic region IRESs by their factor-independent binding to the P site of 80S ribosomes

Initiation of translation on many viral mRNAs occurs by noncanonical mechanisms that involve 5' end-independent binding of ribosomes to an internal ribosome entry site (IRES). The ∼190-nt-long intergenic region (IGR) IRES of dicistroviruses such as cricket paralysis virus (CrPV) initiates translation without Met-tRNAi Met or initiation factors. Advances in metagenomics have revealed numerous dicistrovirus-like genomes with shorter, structurally distinct IGRs, such as nedicistrovirus (NediV) and Antarctic picorna-like virus 1 (APLV1). Like canonical IGR IRESs, the ∼165-nt-long NediV-like IGRs comprise three domains, but they lack key canonical motifs, including L1.1a/L1.1b loops (which bind to the L1 stalk of the ribosomal 60S subunit) and the apex of stem-loop V (SLV) (which binds to the head of the 40S subunit). Domain 2 consists of a compact, highly conserved pseudoknot (PKIII) that contains a UACUA loop motif and a protruding CrPV-like stem--loop SLIV. In vitro reconstitution experiments showed that NediV-like IRESs initiate translation from a non-AUG codon and form elongation-competent 80S ribosomal complexes in the absence of initiation factors and Met-tRNAi Met Unlike canonical IGR IRESs, NediV-like IRESs bind directly to the peptidyl (P) site of ribosomes leaving the aminoacyl (A) site accessible for decoding. The related structures of NediV-like IRESs and their common mechanism of action indicate that they exemplify a distinct class of IGR IRES.

The compact genome of the sponge Oopsacas minuta (Hexactinellida) is lacking key metazoan core genes

BACKGROUND

Explaining the emergence of the hallmarks of bilaterians is a central focus of evolutionary developmental biology-evodevo-and evolutionary genomics. For this purpose, we must both expand and also refine our knowledge of non-bilaterian genomes, especially by studying early branching animals, in particular those in the metazoan phylum Porifera.

RESULTS

We present a comprehensive analysis of the first whole genome of a glass sponge, Oopsacas minuta, a member of the Hexactinellida. Studying this class of sponge is evolutionary relevant because it differs from the three other Porifera classes in terms of development, tissue organization, ecology, and physiology. Although O. minuta does not exhibit drastic body simplifications, its genome is among the smallest of animal genomes sequenced so far, and surprisingly lacks several metazoan core genes (including Wnt and several key transcription factors). Our study also provides the complete genome of a symbiotic Archaea dominating the associated microbial community: a new Thaumarchaeota species.

CONCLUSIONS

The genome of the glass sponge O. minuta differs from all other available sponge genomes by its compactness and smaller number of encoded proteins. The unexpected loss of numerous genes previously considered ancestral and pivotal for metazoan morphogenetic processes most likely reflects the peculiar syncytial tissue organization in this group. Our work further documents the importance of convergence during animal evolution, with multiple convergent evolution of septate-like junctions, electrical-signaling and multiciliated cells in metazoans.

Plasticity, Paralogy, and Pseudogenization: Rhabdoviruses of Freshwater Mussels Elucidate Mechanisms of Viral Genome Diversification and the Evolution of the Finfish-Infecting Rhabdoviral Genera

Viruses in the family Rhabdoviridae display remarkable genomic variation and ecological diversity. This plasticity occurs despite the fact that, as negative sense RNA viruses, rhabdoviruses rarely if ever recombine. Here, we describe nonrecombinatorial evolutionary processes leading to genomic diversification in the Rhabdoviridae inferred from two novel rhabdoviruses of freshwater mussels (Mollusca: Bivalvia: Unionida). Killamcar virus 1 (KILLV-1) from a plain pocketbook (Lampsilis cardium) is closely related phylogenetically and transcriptionally to finfish-infecting viruses in the subfamily Alpharhabdovirinae. KILLV-1 offers a novel example of glycoprotein gene duplication, differing from previous examples in that the paralogs overlap. Evolutionary analyses reveal a clear pattern of relaxed selection due to subfunctionalization in rhabdoviral glycoprotein paralogs, which has not previously been described in RNA viruses. Chemarfal virus 1 (CHMFV-1) from a western pearlshell (Margaritifera falcata) is closely related phylogenetically and transcriptionally to viruses in the genus Novirhabdovirus, the sole recognized genus in the subfamily Gammarhabdovirinae, representing the first known gammarhabdovirus of a host other than finfish. The CHMFV-1 G-L noncoding region contains a nontranscribed remnant gene of precisely the same length as the NV gene of most novirhabdoviruses, offering a compelling example of pseudogenization. The unique reproductive strategy of freshwater mussels involves an obligate parasitic stage in which larvae encyst in the tissues of finfish, offering a plausible ecological mechanism for viral host-switching. IMPORTANCE Viruses in the family Rhabdoviridae infect a variety of hosts, including vertebrates, invertebrates, plants and fungi, with important consequences for health and agriculture. This study describes two newly discovered viruses of freshwater mussels from the United States. One virus from a plain pocketbook (Lampsilis cardium) is closely related to fish-infecting viruses in the subfamily Alpharhabdovirinae. The other virus from a western pearlshell (Margaritifera falcata) is closely related to viruses in the subfamily Gammarhabdovirinae, which until now were only known to infect finfish. Genome features of both viruses provide new evidence of how rhabdoviruses evolved their extraordinary variability. Freshwater mussel larvae attach to fish and feed on tissues and blood, which may explain how rhabdoviruses originally jumped between mussels and fish. The significance of this research is that it improves our understanding of rhabdovirus ecology and evolution, shedding new light on these important viruses and the diseases they cause.

Identification of RNA Virus-Derived RdRp Sequences in Publicly Available Transcriptomic Data Sets

RNA viruses are abundant and highly diverse and infect all or most eukaryotic organisms. However, only a tiny fraction of the number and diversity of RNA virus species have been catalogued. To cost-effectively expand the diversity of known RNA virus sequences, we mined publicly available transcriptomic data sets. We developed 77 family-level Hidden Markov Model profiles for the viral RNA-dependent RNA polymerase (RdRp)-the only universal "hallmark" gene of RNA viruses. By using these to search the National Center for Biotechnology Information Transcriptome Shotgun Assembly database, we identified 5,867 contigs encoding RNA virus RdRps or fragments thereof and analyzed their diversity, taxonomic classification, phylogeny, and host associations. Our study expands the known diversity of RNA viruses, and the 77 curated RdRp Profile Hidden Markov Models provide a useful resource for the virus discovery community.

Viruses of a key coral symbiont exhibit temperature-driven productivity across a reefscape

Viruses can affect coral health by infecting their symbiotic dinoflagellate partners (Symbiodiniaceae). Yet, viral dynamics in coral colonies exposed to environmental stress have not been studied at the reef scale, particularly within individual viral lineages. We sequenced the viral major capsid protein (mcp) gene of positive-sense single-stranded RNA viruses known to infect symbiotic dinoflagellates ('dinoRNAVs') to analyze their dynamics in the reef-building coral, Porites lobata. We repeatedly sampled 54 colonies harboring Cladocopium C15 dinoflagellates, across three environmentally distinct reef zones (fringing reef, back reef, and forereef) around the island of Moorea, French Polynesia over a 3-year period and spanning a reef-wide thermal stress event. By the end of the sampling period, 28% (5/18) of corals in the fringing reef experienced partial mortality versus 78% (14/18) of corals in the forereef. Over 90% (50/54) of colonies had detectable dinoRNAV infections. Reef zone influenced the composition and richness of viral mcp amino acid types ('aminotypes'), with the fringing reef containing the highest aminotype richness. The reef-wide thermal stress event significantly increased aminotype dispersion, and this pattern was strongest in the colonies that experienced partial mortality. These findings demonstrate that dinoRNAV infections respond to environmental fluctuations experienced in situ on reefs. Further, viral productivity will likely increase as ocean temperatures continue to rise, potentially impacting the foundational symbiosis underpinning coral reef ecosystems.

Origins and diversification of animal innate immune responses against viral infections

Immune systems are of pivotal importance to any living organism on Earth, as they protect the organism against deleterious effects of viral infections. Though the current knowledge about these systems is still biased towards the immune response in vertebrates, some studies have focused on the identification and characterization of components of invertebrate antiviral immune systems. Two classic model organisms, the insect Drosophila melanogaster and the nematode Caenorhabditis elegans, were instrumental in the discovery of several important components of the innate immune system, such as the Toll-like receptors and the RNA interference pathway. However, these two model organisms provide only a limited view of the evolutionary history of the immune system, as they both are ecdysozoan protostomes. Recent functional studies in non-classic models such as unicellular holozoans (for example, choanoflagellates), lophotrochozoans (for example, oysters) and cnidarians (for example, sea anemones) have added crucial information for understanding the evolution of antiviral systems, as they revealed unexpected ancestral complexity. This Review aims to summarize this information and present the ancestral nature of the antiviral immune response in animals. We also discuss lineage-specific adaptations and future perspectives for the comparative study of the innate immune system that are essential for understanding its evolution.

Small RNAs in Cnidaria: A review

As fundamental components of RNA silencing, small RNA (sRNA) molecules ranging from 20 to 32 nucleotides in length have been found as potent regulators of gene expression and genome stability in many biological processes of eukaryotes. Three major small RNAs are active in animals, including the microRNA (miRNA), short interfering RNA (siRNA), and PIWI-interacting RNA (piRNA). Cnidarians, the sister group to bilaterians, are at a critical phylogenetic node to better model eukaryotic small RNA pathway evolution. To date, most of our understanding of sRNA regulation and its potential contribution to evolution has been limited to a few triploblastic bilaterian and plant models. The diploblastic nonbilaterians, including the cnidarians, are understudied in this regard. Therefore, this review will present the current-known small RNA information in cnidarians to enhance our understanding of the development of the small RNA pathways in early branch animals.

A novel narnavirus is widespread in Saccharomyces cerevisiae and impacts multiple host phenotypes

RNA viruses are a widespread, biologically diverse group that includes the narnaviridiae, a family of unencapsidated RNA viruses containing a single ORF that encodes an RNA-dependent RNA polymerase. In the yeast Saccharomyces cerevisiae, the 20S and 23S RNA viruses are well-studied members of the narnaviridiae, which are present at low intracellular copy numbers, unless induced by stress or unfavorable growth conditions, and are not known to affect host fitness. In this study, we describe a new S. cerevisiae narnavirus that we designate as N1199. We show that N1199 is uniquely present as a double-stranded RNA at a high level relative to other known members of this family in 1 strain background, YJM1199, and is present as a single-stranded RNA at lower levels in 98 of the remaining 100-genomes strains. Furthermore, we see a strong association between the presence of high level N1199 and host phenotype defects, including greatly reduced sporulation efficiency and growth on multiple carbon sources. Finally, we describe associations between N1199 abundance and host phenotype defects, including autophagy.

Diverse RNA Viruses Associated with Diatom, Eustigmatophyte, Dinoflagellate, and Rhodophyte Microalgae Cultures

Unicellular microalgae are of immense ecological importance with growing commercial potential in industries such as renewable energy, food, and pharmacology. Viral infections can have a profound impact on the growth and evolution of their hosts. However, very little is known of the diversity within, and the effect of, unicellular microalgal RNA viruses. In addition, identifying RNA viruses in these organisms that could have originated more than a billion years ago constitutes a robust data set to dissect molecular events and address fundamental questions in virus evolution. We assessed the diversity of RNA viruses in eight microalgal cultures, including representatives from the diatom, eustigmatophyte, dinoflagellate, red algae, and euglenid groups. Using metatranscriptomic sequencing combined with bioinformatic approaches optimized to detect highly divergent RNA viruses, we identified 10 RNA virus sequences, with nine constituting new viral species. Most of the newly identified RNA viruses belonged to the double-stranded Totiviridae, Endornaviridae, and Partitiviridae, greatly expanding the reported host range for these families. Two new species belonging to the single-stranded RNA viral clade Marnaviridae, commonly associated with microalgal hosts, were also identified. This study highlights that a substantial diversity of RNA viruses likely exists undetected within the unicellular microalgae. It also highlights the necessity for RNA viral characterization and for investigation of the effects of viral infections on microalgal physiology, biology, and growth, considering their environmental and industrial roles. IMPORTANCE Our knowledge of the diversity of RNA viruses infecting microbial algae-the microalgae-is minimal. However, describing the RNA viruses infecting these organisms is of primary importance at both the ecological and economic scales because of the fundamental roles these organisms play in aquatic environments and their growing value across a range of industrial fields. Using metatranscriptomic sequencing, we aimed to reveal the RNA viruses present in cultures of eight microalgae species belonging to the diatom, dinoflagellate, eustigmatophyte, rhodophyte, and euglena major clades of algae. Accordingly, we identified 10 new divergent RNA virus species belonging to RNA virus families as diverse as the double-stranded Totiviridae, Endornaviridae, and Partitiviridae and the single-stranded Marnaviridae. By expanding the known diversity of RNA viruses infecting unicellular eukaryotes, this study contributes to a better understanding of the early evolution of the virosphere and will inform the use of microalgae in industrial applications.

Multiple waves of viral invasions in Symbiodiniaceae algal genomes

Dinoflagellates from the family Symbiodiniaceae are phototrophic marine protists that engage in symbiosis with diverse hosts. Their large and distinct genomes are characterized by pervasive gene duplication and large-scale retroposition events. However, little is known about the role and scale of horizontal gene transfer (HGT) in the evolution of this algal family. In other dinoflagellates, high levels of HGTs have been observed, linked to major genomic transitions, such as the appearance of a viral-acquired nucleoprotein that originated via HGT from a large DNA algal virus. Previous work showed that Symbiodiniaceae from different hosts are actively infected by viral groups, such as giant DNA viruses and ssRNA viruses, that may play an important role in coral health. Latent viral infections may also occur, whereby viruses could persist in the cytoplasm or integrate into the host genome as a provirus. This hypothesis received experimental support; however, the cellular localization of putative latent viruses and their taxonomic affiliation are still unknown. In addition, despite the finding of viral sequences in some genomes of Symbiodiniaceae, viral origin, taxonomic breadth, and metabolic potential have not been explored. To address these questions, we searched for putative viral-derived proteins in thirteen Symbiodiniaceae genomes. We found fifty-nine candidate viral-derived HGTs that gave rise to twelve phylogenies across ten genomes. We also describe the taxonomic affiliation of these virus-related sequences, their structure, and their genomic context. These results lead us to propose a model to explain the origin and fate of Symbiodiniaceae viral acquisitions.

Freshwater macrophytes harbor viruses representing all five major phyla of the RNA viral kingdom Orthornavirae

Research on aquatic plant viruses is lagging behind that of their terrestrial counterparts. To address this knowledge gap, here we identified viruses associated with freshwater macrophytes, a taxonomically diverse group of aquatic phototrophs that are visible with the naked eye. We surveyed pooled macrophyte samples collected at four spring sites in Florida, USA through next generation sequencing of RNA extracted from purified viral particles. Sequencing efforts resulted in the detection of 156 freshwater macrophyte associated (FMA) viral contigs, 37 of which approximate complete genomes or segments. FMA viral contigs represent putative members from all five major phyla of the RNA viral kingdom Orthornavirae. Similar to viral types found in land plants, viral sequences identified in macrophytes were dominated by positive-sense RNA viruses. Over half of the FMA viral contigs were most similar to viruses reported from diverse hosts in aquatic environments, including phototrophs, invertebrates, and fungi. The detection of FMA viruses from orders dominated by plant viruses, namely Patatavirales and Tymovirales, indicate that members of these orders may thrive in aquatic hosts. PCR assays confirmed the presence of putative FMA plant viruses in asymptomatic vascular plants, indicating that viruses with persistent lifestyles are widespread in macrophytes. The detection of potato virus Y and oat blue dwarf virus in submerged macrophytes suggests that terrestrial plant viruses infect underwater plants and highlights a potential terrestrial-freshwater plant virus continuum. Defining the virome of unexplored macrophytes will improve our understanding of virus evolution in terrestrial and aquatic primary producers and reveal the potential ecological impacts of viral infection in macrophytes.

The origin of RNA interference: Adaptive or neutral evolution?

The origin of RNA interference (RNAi) is usually explained by a defense-based hypothesis, in which RNAi evolved as a defense against transposable elements (TEs) and RNA viruses and was already present in the last eukaryotic common ancestor (LECA). However, since RNA antisense regulation and double-stranded RNAs (dsRNAs) are ancient and widespread phenomena, the origin of defensive RNAi should have occurred in parallel with its regulative functions to avoid imbalances in gene regulation. Thus, we propose a neutral evolutionary hypothesis for the origin of RNAi in which qualitative system drift from a prokaryotic antisense RNA gene regulation mechanism leads to the formation of RNAi through constructive neutral evolution (CNE). We argue that RNAi was already present in the ancestor of LECA before the need for a new defense system arose and that its presence helped to shape eukaryotic genomic architecture and stability.

Where does RNA interference come from? This Essay describes a new step-by-step evolutionary model of how RNA interference might have originated in early eukaryotes through neutral events from the molecular machinery present in prokaryotes.

Metatranscriptomic Comparison of Viromes in Endemic and Introduced Passerines in New Zealand

New Zealand/Aotearoa has many endemic passerine birds vulnerable to emerging infectious diseases. Yet little is known about viruses in passerines, and in some countries, including New Zealand, the virome of wild passerines has been only scarcely researched. Using metatranscriptomic sequencing we characterised the virome of New Zealand endemic and introduced species of passerine. Accordingly, we identified 34 possible avian viruses from cloacal swabs of 12 endemic and introduced bird species not showing signs of disease. These included a novel siadenovirus, iltovirus, and avastrovirus in the Eurasian blackbird (Turdus merula, an introduced species), song thrush (Turdus philomelos, introduced) and silvereye/tauhou (Zosterops lateralis, introduced), respectively. This is the first time novel viruses from these genera have been identified in New Zealand, likely reflecting prior undersampling. It also represents the first identification of an iltovirus and siadenovirus in blackbirds and thrushes globally. These three viruses were only found in introduced species and may pose a risk to endemic species if they were to jump species boundaries, particularly the iltoviruses and siadenoviruses that have a prior history of disease associations. Further virus study and surveillance are needed in New Zealand avifauna, particularly in Turdus populations and endemic species.

High Diversity of Novel Viruses in the Tree Pathogen Phytophthora castaneae Revealed by High-Throughput Sequencing of Total and Small RNA

Phytophthora castaneae, an oomycete pathogen causing root and trunk rot of different tree species in Asia, was shown to harbor a rich diversity of novel viruses from different families. Four P. castaneae isolates collected from Chamaecyparis hodginsii in a semi-natural montane forest site in Vietnam were investigated for viral presence by traditional and next-generation sequencing (NGS) techniques, i.e., double-stranded RNA (dsRNA) extraction and high-throughput sequencing (HTS) of small RNAs (sRNAs) and total RNA. Genome organization, sequence similarity, and phylogenetic analyses indicated that the viruses were related to members of the order Bunyavirales and families Endornaviridae, Megabirnaviridae, Narnaviridae, Totiviridae, and the proposed family "Fusagraviridae." The study describes six novel viruses: Phytophthora castaneae RNA virus 1-5 (PcaRV1-5) and Phytophthora castaneae negative-stranded RNA virus 1 (PcaNSRV1). All six viruses were detected by sRNA sequencing, which demonstrates an active RNA interference (RNAi) system targeting viruses in P. castaneae. To our knowledge, this is the first report of viruses in P. castaneae and the whole Phytophthora major Clade 5, as well as of the activity of an RNAi mechanism targeting viral genomes among Clade 5 species. PcaRV1 is the first megabirnavirus described in oomycetes and the genus Phytophthora.

Daphnia magna egg piRNA cluster expression profiles change as mothers age

BACKGROUND

PiRNAs prevent transposable elements wreaking havoc on the germline genome. Changes in piRNA expression over the lifetime of an individual may impact on ageing through continued suppression, or release, of transposable element expression. We identified piRNA producing clusters in the genome of Daphnia magna by a combination of bioinformatic methods, and then contrasted their expression between parthenogenetically produced eggs representing maternally-deposited germline piRNAs of young (having their 1st clutch) and old (having their 5th clutch) mothers. Results from eggs were compared to cluster expression in three generations of adults.

RESULTS

As for other arthropods, D. magna encodes long uni-directionally transcribed non-coding RNAs consisting of fragmented transposable elements which account for most piRNAs expressed. Egg tissues showed extensive differences between clutches from young mothers and those from old mothers, with 578 and 686 piRNA clusters upregulated, respectively. Most log fold-change differences for significant clusters were modest, however. When considering only highly expressed clusters, there was a bias towards 1st clutch eggs at 41 upregulated versus eight clusters in the eggs from older mothers. F0 generation differences between young and old mothers were fewer than eggs, as 179 clusters were up-regulated in young versus 170 old mothers. This dropped to 31 versus 22 piRNA clusters when comparing adults in the F1 generation, and no differences were detected in the F3 generation. Inter-generational losses of differential piRNA cluster were similar to that observed for D. magna micro-RNA expression.

CONCLUSIONS

Little overlap in differentially expressed clusters was found between adults containing mixed somatic and germline (ovary) tissues and germ-line representing eggs. A cluster encompassing a Tudor domain containing gene important in the piRNA pathway was upregulated in the eggs from old mothers. We hypothesise that regulation of this gene could form part of a feedback loop that reduces piRNA pathway activity explaining the reduced number of highly-expressed clusters in eggs from old mothers.

The RNA virome of echinoderms

Echinoderms are a phylum of marine invertebrates that include model organisms, keystone species, and animals commercially harvested for seafood. Despite their scientific, ecological, and economic importance, there is little known about the diversity of RNA viruses that infect echinoderms compared to other invertebrates. We screened over 900 transcriptomes and viral metagenomes to characterize the RNA virome of 38 echinoderm species from all five classes (Crinoidea, Holothuroidea, Asteroidea, Ophiuroidea and Echinoidea). We identified 347 viral genome fragments that were classified to genera and families within nine viral orders - Picornavirales, Durnavirales, Martellivirales, Nodamuvirales, Reovirales, Amarillovirales, Ghabrivirales, Mononegavirales, and Hepelivirales. We compared the relative viral representation across three life stages (embryo, larvae, adult) and characterized the gene content of contigs which encoded complete or near-complete genomes. The proportion of viral reads in a given transcriptome was not found to significantly differ between life stages though the majority of viral contigs were discovered from transcriptomes of adult tissue. This study illuminates the biodiversity of RNA viruses from echinoderms, revealing the occurrence of viral groups in natural populations.

Tracing RNA viruses associated with Nudibranchia gastropods

Background

Nudibranchia is an under-studied taxonomic group of gastropods, including more than 3,000 species with colourful and extravagant body shapes and peculiar predatory and defensive strategies. Although symbiosis with bacteria has been reported, no data are available for the nudibranch microbiome nor regarding viruses possibly associated with these geographically widespread species.

Methods

Based on 47 available RNA sequencing datasets including more than two billion reads of 35 nudibranch species, a meta-transcriptome assembly was constructed. Taxonomic searches with DIAMOND, RNA-dependent-RNA-polymerase identification with palmscan and viral hallmark genes identification by VirSorter2 in combination with CheckV were applied to identify genuine viral genomes, which were then annotated using CAT.

Results

A total of 20 viral genomes were identified as bona fide viruses, among 552 putative viral contigs resembling both RNA viruses of the Negarnaviricota, Pisuviricota, Kitrinoviricota phyla and actively transcribing DNA viruses of the Cossaviricota and Nucleocytoviricota phyla. The 20 commonly identified viruses showed similarity with RNA viruses identified in other RNA-seq experiments and can be putatively associated with bacteria, plant and arthropod hosts by co-occurence analysis. The RNA samples having the highest viral abundances showed a heterogenous and mostly sample-specific distribution of the identified viruses, suggesting that nudibranchs possess diversified and mostly unknown viral communities.

Human Land-Use Impacts Viral Diversity and Abundance in a New Zealand River

Although water-borne viruses have important implications for the health of humans and other animals, little is known about the impact of human land use on viral diversity and evolution in water systems such as rivers. We used metatranscriptomic sequencing to compare the diversity and abundance of viruses at sampling sites along a single river in New Zealand that differed in human land-use impacts, ranging from pristine to urban. From this, we identified 504 putative virus species, of which 97 per cent were novel. Many of the novel viruses were highly divergent and likely included a new subfamily within the Parvoviridae. We identified at least sixty-three virus species that may infect vertebrates—most likely fish and water birds—from the Astroviridae, Birnaviridae, Parvoviridae, and Picornaviridae. No putative human viruses were detected. Importantly, we observed differences in the composition of viral communities at sites impacted by human land use (farming and urban) compared to native forest sites (pristine). At the viral species level, the urban sites had higher diversity (327 virus species) than the farming (n = 150) and pristine sites (n = 119), and more viruses were shared between the urban and farming sites (n = 76) than between the pristine and farming or urban sites (n = 24). The two farming sites had a lower viral abundance across all host types, while the pristine sites had a higher abundance of viruses associated with animals, plants, and fungi. We also identified viruses linked to agriculture and human impact at the river sampling sites in farming and urban areas that were not present at the native forest sites. Although based on a small sample size, our study suggests that human land use can impact viral communities in rivers, such that further work is needed to reduce the impact of intensive farming and urbanisation on water systems.

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.

The Small RNA Universe of Capitella teleta

RNAi is an evolutionarily fluid mechanism with dramatically different activities across animal phyla. One major group where there has been little investigation is annelid worms. Here, the small RNAs of the polychaete developmental model Capitella teleta are profiled across development. As is seen with nearly all animals, nearly 200 microRNAs were found with 58 high-confidence novel species. Greater miRNA diversity was associated with later stages consistent with differentiation of tissues. Outside miRNA, a distinct composition of other small RNA pathways was found. Unlike many invertebrates, an endogenous siRNA pathway was not observed, indicating pathway loss relative to basal planarians. No processively generated siRNA-class RNAs could be found arising from dsRNA precursors. This has a significant impact on RNAi technology development for this group of animals. Unlike the apparent absence of siRNAs, a significant population of piRNAs was observed. For many piRNAs, phasing and ping-pong biogenesis pathways were identified. Interestingly, piRNAs were found to be highly expressed during early development, suggesting a potential role in regulation in metamorphosis. Critically, the configuration of RNAi factors in C. teleta is found in other annelids and mollusks, suggesting that similar biology is likely to be present in the wider clade. This study is the first in providing comprehensive analysis of small RNAs in annelids.

Metagenomic Assessment of DNA Viral Diversity in Freshwater Sponges, Baikalospongia bacillifera

Sponges (type Porifera) are multicellular organisms that give shelter to a variety of microorganisms: fungi, algae, archaea, bacteria, and viruses. The studies concerning the composition of viral communities in sponges have appeared rather recently, and the diversity and role of viruses in sponge holobionts remain largely undisclosed. In this study, we assessed the diversity of DNA viruses in the associated community of the Baikal endemic sponge, Baikalospongia bacillifera, using a metagenomic approach, and compared the virome data from samples of sponges and Baikal water (control sample). Significant differences in terms of taxonomy, putative host range of identified scaffolds, and functional annotation of predicted viral proteins were revealed in viromes of sponge B. bacillifera and the Baikal water. This is the evidence in favor of specificity of viral communities in sponges. The diversity shift of viral communities in a diseased specimen, in comparison with a visually healthy sponge, probably reflects the changes in the composition of microbial communities in affected sponges. We identified many viral genes encoding the proteins with metabolic functions; therefore, viruses in Baikal sponges regulate the number and diversity of their associated community, and also take a part in the vital activity of the holobiont, and this is especially significant in the case of damage (or disease) of these organisms in unfavorable conditions. When comparing the Baikal viromes with similar datasets of marine sponge (Ianthella basta), in addition to significant differences in the taxonomic and functional composition of viral communities, we revealed common scaffolds/virotypes in the cross-assembly of reads, which may indicate the presence of some closely related sponge-specific viruses in marine and freshwater sponges.

cGAS-like receptor-mediated immunity: the insect perspective

The cGAS-STING pathway plays a central role in the detection of DNA in the cytosol of mammalian cells and activation of immunity. Although the early evolutionary origin of this pathway in animals has been noted, its ancestral functions have remained elusive so far. We review here new findings in invertebrates establishing a role in sensing and signaling infection, triggering potent transcriptional responses, in addition to autophagy. Results from flies and moths/butterflies point to the importance of STING signaling in antiviral immunity in insects. The recent characterization of cGAS-like receptors in Drosophila reveals the plasticity of this family of pattern-recognition receptors, able to accommodate ligands different from DNA and to produce cyclic dinucleotides beyond 2'3'-cGAMP.

A glimpse on metazoan ZNFX1 helicases, ancient players of antiviral innate immunity

The human zinc finger NFX1-type containing 1 (ZNFX1) is an interferon-stimulated protein associated to the outer mitochondrial membrane, able to bind dsRNAs and interact with MAVS proteins, promoting type I IFN response in the early stage of viral infection. An N-terminal Armadillo (ARM)-type fold and a large helicase core (P-loop) and zinc fingers confer RNA-binding and ATPase activities to ZNFX1. We studied the phylogenetic distribution of metazoan ZNFX1s, ZNFX1 gene expression trends and genomic and protein signatures during viral infection of invertebrates. Based on 221 ZNFX1 sequences, we obtained a polyphyletic tree with a taxonomy-consistent branching at the phylum-level only. In metazoan genomes, ZNFX1 genes were found either in single copy, with up to some tens of exons in vertebrates, or in multiple copies, with one or a few exons and one of them sometimes encompassing most of the coding sequence, in invertebrates like sponges, sea urchins and mollusks. Structural analyses of selected ZNFX1 proteins showed high conservation of the helicase region (P-loop), an overall conserved region and domain architecture, an ARM-fold mostly traceable, and the presence of intrinsically disordered regions of varying length and position. The remarkable over-expression of ZNFX1 in bivalve and gastropod mollusks infected with dsDNA viruses underscores the antiviral role of ZNFX1, whereas nothing similar was found in virus-infected nematodes and corals. Whether the functional diversification reported in the C. elegans ZNFX1 occurs in other metazoan proteins remains to be established.

Exploring Virome Diversity in Public Data in South America as an Approach for Detecting Viral Sources From Potentially Emerging Viruses

The South American continent presents a great diversity of biomes, whose ecosystems are constantly threatened by the expansion of human activity. The emergence and re-emergence of viral populations with impact on the human population and ecosystem have shown increases in the last decades. In deference to the growing accumulation of genomic data, we explore the potential of South American-related public databases to detect signals that contribute to virosphere research. Therefore, our study aims to investigate public databases with emphasis on the surveillance of viruses with medical and ecological relevance. Herein, we profiled 120 "sequence read archives" metagenomes from 19 independent projects from the last decade. In a coarse view, our analyses identified only 0.38% of the total number of sequences from viruses, showing a higher proportion of RNA viruses. The metagenomes with the most important viral sequences in the analyzed environmental models were 1) aquatic samples from the Amazon River, 2) sewage from Brasilia, and 3) soil from the state of São Paulo, while the models of animal transmission were detected in mosquitoes from Rio Janeiro and Bats from Amazonia. Also, the classification of viral signals into operational taxonomic units (OTUs) (family) allowed us to infer from metadata a probable host range in the virome detected in each sample analyzed. Further, several motifs and viral sequences are related to specific viruses with emergence potential from Togaviridae, Arenaviridae, and Flaviviridae families. In this context, the exploration of public databases allowed us to evaluate the scope and informative capacity of sequences from third-party public databases and to detect signals related to viruses of clinical or environmental importance, which allowed us to infer traits associated with probable transmission routes or signals of ecological disequilibrium. The evaluation of our results showed that in most cases the size and type of the reference database, the percentage of guanine-cytosine (GC), and the length of the query sequences greatly influence the taxonomic classification of the sequences. In sum, our findings describe how the exploration of public genomic data can be exploited as an approach for epidemiological surveillance and the understanding of the virosphere.

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

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

Revealing RNA virus diversity and evolution in unicellular algae transcriptomes

Remarkably little is known about the diversity and evolution of RNA viruses in unicellular eukaryotes. We screened a total of 570 transcriptomes from the Marine Microbial Eukaryote Transcriptome Sequencing Project that encompasses a wide diversity of microbial eukaryotes, including most major photosynthetic lineages (i.e. the microalgae). From this, we identified thirty new and divergent RNA virus species, occupying a range of phylogenetic positions within the overall diversity of RNA viruses. Approximately one-third of the newly described viruses comprised single-stranded positive-sense RNA viruses from the order Lenarviricota associated with fungi, plants, and protists, while another third were related to the order Ghabrivirales, including members of the protist and fungi-associated Totiviridae. Other viral species showed sequence similarity to positive-sense RNA viruses from the algae-associated Marnaviridae, the double-stranded RNA (ds-RNA) Partitiviridae, as well as tentative evidence for one negative-sense RNA virus related to the Qinviridae. Importantly, we were able to identify divergent RNA viruses from distant host taxa, revealing the ancestry of these viral families and greatly extending our knowledge of the RNA viromes of microalgal cultures. Both the limited number of viruses detected per sample and the low sequence identity to known RNA viruses imply that additional microalgal viruses exist that could not be detected at the current sequencing depth or were too divergent to be identified using sequence similarity. Together, these results highlight the need for further investigation of algal-associated RNA viruses as well as the development of new tools to identify RNA viruses that exhibit very high levels of sequence divergence.

Functional characterization of the cnidarian antiviral immune response reveals ancestral complexity

Animals evolved a broad repertoire of innate immune sensors and downstream effector cascades for defense against RNA viruses. Yet, this system varies greatly among different bilaterian animals, masking its ancestral state. In this study, we aimed to characterize the antiviral immune response of the cnidarian Nematostella vectensis and decipher the function of the retinoic acid-inducible gene I (RIG-I)-like receptors (RLRs) known to detect viral double-stranded RNA (dsRNA) in bilaterians but activate different antiviral pathways in vertebrates and nematodes. We show that polyinosinic:polycytidylic acid (poly(I:C)), a mimic of long viral dsRNA and a primary ligand for the vertebrate RLR melanoma differentiation-associated protein 5 (MDA5), triggers a complex antiviral immune response bearing features distinctive for both vertebrate and invertebrate systems. Importantly, a well-characterized agonist of the vertebrate RIG-I receptor does not induce a significant transcriptomic response that bears signature of the antiviral immune response, which experimentally supports the results of a phylogenetic analysis indicating clustering of the two N. vectensis RLR paralogs (NveRLRa and NveRLRb) with MDA5. Furthermore, the results of affinity assays reveal that NveRLRb binds poly(I:C) and long dsRNA and its knockdown impairs the expression of putative downstream effector genes including RNA interference components. Our study provides for the first time the functional evidence for the conserved role of RLRs in initiating immune response to dsRNA that originated before the cnidarian–bilaterian split and lay a strong foundation for future research on the evolution of the immune responses to RNA viruses.

Digging into bivalve miRNAomes: between conservation and innovation

Bivalves are a diverse mollusc group of economic and ecological importance. An evident resilience to pollution, parasites and extreme environments makes some bivalve species important models for studying adaptation and immunity. Despite substantial progress in sequencing projects of bivalves, information on non-coding genes and gene-regulatory aspects is still lacking. Here, we review the current repertoire of bivalve microRNAs (miRNAs), important regulators of gene expression in Metazoa. We exploited available short non-coding RNA (sncRNA) data for Pinctada martensii, Crassostrea gigas, Corbicula fluminea, Tegillarca granosa and Ruditapes philippinarum, and we produced new sncRNA data for two additional bivalves, the Mediterranean mussel Mytilus galloprovincialis and the blood clam Scapharca broughtonii. We found substantial heterogeneity and incorrect annotations of miRNAs; hence, we reannotated conserved miRNA families using recently established criteria for bona fide microRNA annotation. We found 106 miRNA families missing in the previously published bivalve datasets and 89 and 87 miRNA complements were identified in the two additional species. The overall results provide a homogeneous and evolutionarily consistent picture of miRNAs in bivalves and enable future comparative studies. The identification of two bivalve-specific miRNA families sheds further light on the complexity of transcription and its regulation in bivalve molluscs. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

Single mosquito metatranscriptomics identifies vectors, emerging pathogens and reservoirs in one assay

Mosquitoes are major infectious disease-carrying vectors. Assessment of current and future risks associated with the mosquito population requires knowledge of the full repertoire of pathogens they carry, including novel viruses, as well as their blood meal sources. Unbiased metatranscriptomic sequencing of individual mosquitoes offers a straightforward, rapid, and quantitative means to acquire this information. Here, we profile 148 diverse wild-caught mosquitoes collected in California and detect sequences from eukaryotes, prokaryotes, 24 known and 46 novel viral species. Importantly, sequencing individuals greatly enhanced the value of the biological information obtained. It allowed us to (a) speciate host mosquito, (b) compute the prevalence of each microbe and recognize a high frequency of viral co-infections, (c) associate animal pathogens with specific blood meal sources, and (d) apply simple co-occurrence methods to recover previously undetected components of highly prevalent segmented viruses. In the context of emerging diseases, where knowledge about vectors, pathogens, and reservoirs is lacking, the approaches described here can provide actionable information for public health surveillance and intervention decisions.

Surveillance of densoviruses and mesomycetozoans inhabiting grossly normal tissues of three Aotearoa New Zealand asteroid species

Asteroid wasting events and mass mortality have occurred for over a century. We currently lack a fundamental understanding of the microbial ecology of asteroid disease, with disease investigations hindered by sparse information about the microorganisms associated with grossly normal specimens. We surveilled viruses and protists associated with grossly normal specimens of three asteroid species (Patiriella regularis, Stichaster australis, Coscinasterias muricata) on the North Island / Te Ika-a-Māui, Aotearoa New Zealand, using metagenomes prepared from virus and ribosome-sized material. We discovered several densovirus-like genome fragments in our RNA and DNA metagenomic libraries. Subsequent survey of their prevalence within populations by quantitative PCR (qPCR) demonstrated their occurrence in only a few (13%) specimens (n = 36). Survey of large and small subunit rRNAs in metagenomes revealed the presence of a mesomycete (most closely matching Ichthyosporea sp.). Survey of large subunit prevalence and load by qPCR revealed that it is widely detectable (80%) and present predominately in body wall tissues across all 3 species of asteroid. Our results raise interesting questions about the roles of these microbiome constituents in host ecology and pathogenesis under changing ocean conditions.

Sensing and signalling viral infection in drosophila

The fruitfly Drosophila melanogaster is a valuable model to unravel mechanisms of innate immunity, in particular in the context of viral infections. RNA interference, and more specifically the small interfering RNA pathway, is a major component of antiviral immunity in drosophila. In addition, the contribution of inducible transcriptional responses to the control of viruses in drosophila and other invertebrates is increasingly recognized. In particular, the recent discovery of a STING-IKKβ-Relish signalling cassette in drosophila has confirmed that NF-κB transcription factors play an important role in the control of viral infections, in addition to bacterial and fungal infections. Here, we review recent developments in the field, which begin to shed light on the mechanisms involved in sensing of viral infections and in signalling leading to production of antiviral effectors.

Viral RNA Genomes Identified from Marine Macroalgae and a Diatom

Protists provide insights into the diversity and function of RNA viruses in marine systems. Among them, marine macroalgae are good targets for RNA virome analyses because they have a sufficient biomass in nature. However, RNA viruses in macroalgae have not yet been examined in detail, and only partial genome sequences have been reported for the majority of RNA viruses. Therefore, to obtain further insights into the distribution and diversity of RNA viruses associated with marine protists, we herein examined RNA viruses in macroalgae and a diatom. We report the putative complete genome sequences of six novel RNA viruses from two marine macroalgae and one diatom holobiont. Four viruses were not classified into established viral genera or families. Furthermore, a virus classified into Totiviridae showed a genome structure that has not yet been reported in this family. These results suggest that a number of distinct RNA viruses are widespread in a broad range of protists.

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.

Polymorphism of genetic ambigrams

Double synonyms in the genetic code can be used as a tool to test competing hypotheses regarding ambigrammatic narnavirus genomes. Applying the analysis to recent observations of Culex narnavirus 1 and Zhejiang mosquito virus 3 ambigrammatic viruses indicates that the open reading frame on the complementary strand of the segment coding for RNA-dependent RNA polymerase does not code for a functional protein. Culex narnavirus 1 has been shown to possess a second segment, also ambigrammatic, termed 'Robin'. We find a comparable segment for Zhejiang mosquito virus 3, a moderately diverged relative of Culex narnavirus 1. Our analysis of Robin polymorphisms suggests that its reverse open reading frame also does not code for a functional protein. We make a hypothesis about its role.

Evolution of MicroRNA Biogenesis Genes in the Sterlet (Acipenser ruthenus) and Other Polyploid Vertebrates

MicroRNAs play a crucial role in eukaryotic gene regulation. For a long time, only little was known about microRNA-based gene regulatory mechanisms in polyploid animal genomes due to difficulties of polyploid genome assembly. However, in recent years, several polyploid genomes of fish, amphibian, and even invertebrate species have been sequenced and assembled. Here we investigated several key microRNA-associated genes in the recently sequenced sterlet (Acipenser ruthenus) genome, whose lineage has undergone a whole genome duplication around 180 MYA. We show that two paralogs of drosha, dgcr8, xpo1, and xpo5 as well as most ago genes have been retained after the acipenserid-specific whole genome duplication, while ago1 and ago3 genes have lost one paralog. While most diploid vertebrates possess only a single copy of dicer1, we strikingly found four paralogs of this gene in the sterlet genome, derived from a tandem segmental duplication that occurred prior to the last whole genome duplication. ago1,3,4 and exportins1,5 look to be prone to additional segment duplications producing up to four-five paralog copies in ray-finned fishes. We demonstrate for the first time exon microsatellite amplification in the acipenserid drosha2 gene, resulting in a highly variable protein product, which may indicate sub- or neofunctionalization. Paralogous copies of most microRNA metabolism genes exhibit different expression profiles in various tissues and remain functional despite the rediploidization process. Subfunctionalization of microRNA processing gene paralogs may be beneficial for different pathways of microRNA metabolism. Genetic variability of microRNA processing genes may represent a substrate for natural selection, and, by increasing genetic plasticity, could facilitate adaptations to changing environments.

Diversity of Sea Star-Associated Densoviruses and Transcribed Endogenous Viral Elements of Densovirus Origin

The primary interest in sea star densoviruses, specifically SSaDV, has been their association with sea star wasting syndrome (SSWS), a disease that has decimated sea star populations across the West Coast of the United States since 2013. The association of SSaDV with SSWS was originally drawn from metagenomic analysis, which was further studied through field surveys using quantitative PCR (qPCR), with the conclusion that it was the most likely viral candidate in the metagenomic data based on its representation in symptomatic sea stars compared to asymptomatic sea stars. We reexamined the original metagenomic data with additional genomic data sets and found that SSaDV was 1 of 10 densoviruses present in the original data set and was no more represented in symptomatic sea stars than in asymptomatic sea stars. Instead, SSaDV appears to be a widespread, generalist virus that exists among a large diversity of densoviruses present in sea star populations.

A viral etiology of sea star wasting syndrome (SSWS) was originally explored with virus-sized material challenge experiments, field surveys, and metagenomics, leading to the conclusion that a densovirus is the predominant DNA virus associated with this syndrome and, thus, the most promising viral candidate pathogen. Single-stranded DNA viruses are, however, highly diverse and pervasive among eukaryotic organisms, which we hypothesize may confound the association between densoviruses and SSWS. To test this hypothesis and assess the association of densoviruses with SSWS, we compiled past metagenomic data with new metagenomic-derived viral genomes from sea stars collected from Antarctica, California, Washington, and Alaska. We used 179 publicly available sea star transcriptomes to complement our approaches for densovirus discovery. Lastly, we focus the study on sea star-associated densovirus (SSaDV), the first sea star densovirus discovered, by documenting its biogeography and putative tissue tropism. Transcriptomes contained only endogenized densovirus elements similar to the NS1 gene, while numerous extant densoviral genomes were recovered from viral metagenomes. SSaDV was associated with nearly all tested species from southern California to Alaska, and in contrast to previous work, we show that SSaDV is one genotype among a high diversity of densoviruses present in sea stars across the West Coast of the United States and globally that are commonly associated with grossly normal (i.e., healthy or asymptomatic) animals. The diversity and ubiquity of these viruses in sea stars confound the original hypothesis that one densovirus is the etiological agent of SSWS.

IMPORTANCE The primary interest in sea star densoviruses, specifically SSaDV, has been their association with sea star wasting syndrome (SSWS), a disease that has decimated sea star populations across the West Coast of the United States since 2013. The association of SSaDV with SSWS was originally drawn from metagenomic analysis, which was further studied through field surveys using quantitative PCR (qPCR), with the conclusion that it was the most likely viral candidate in the metagenomic data based on its representation in symptomatic sea stars compared to asymptomatic sea stars. We reexamined the original metagenomic data with additional genomic data sets and found that SSaDV was 1 of 10 densoviruses present in the original data set and was no more represented in symptomatic sea stars than in asymptomatic sea stars. Instead, SSaDV appears to be a widespread, generalist virus that exists among a large diversity of densoviruses present in sea star populations.

The functional roles of the non-coding RNAs in molluscs

This review focus on the current knowledge of non-coding RNAs (ncRNAs) in molluscs. In this review, we provide an overview of long ncRNAs (lncRNA), microRNAs (miRNA) and piwi-interacting RNAs (piRNA), followed by evidence for the regulation of ncRNAs in variety of biological process in molluscs, including development, biomineralization and innate immune response. This review advances our understanding on the roles of ncRNAs in molluscs and suggest the future direction to fully understand the epigenetic regulatory network of molluscs.

Endophytic bacterial communities and spatiotemporal variations in cotton roots in Xinjiang, China

Endogenous bacteria are important for maintaining the health and other ecologically relevant functions of cotton plants. However, little is known about the community structures and diversity of endophytic bacteria in cotton plants. In our study, we used the Illumina amplicon sequencing technology to study the endophytic bacteria found in cotton root tissue in Xinjiang, China. A total of 60.84 × 10⁶ effective sequences of the 16S rRNA gene in the V5-V6 variable region revealed a large number of operational taxonomic units (OTUs), namely 81-338 OTUs, at a cut-off level of 3% and a sequencing depth of 50 000 sequences. Among the 23 classes identified, Gammaproteobacteria was the dominant group, followed by Alphaproteobacteria, Actinobacteria, and Bacillus. The diversity of endogenous bacteria differed at different growth periods, with the most OTUs detected in seedlings (654), followed by the budding stage (381), flowering stage (350), and flocking stage (351). A total of 217 OTUs were common to all four stages. Pantoea tags were more common to the Shihezi region, whereas Erwinia labels were more common to the Hami region. These results suggest that the dynamics of endophytic bacterial communities are affected by plant growth stage. This highlights the relevance of microbial diversity studies in improving our understanding of endophyte communities.

Metatranscriptomic Identification of Diverse and Divergent RNA Viruses in Green and Chlorarachniophyte Algae Cultures

Our knowledge of the diversity and evolution of the virosphere will likely increase dramatically with the study of microbial eukaryotes, including the microalgae within which few RNA viruses have been documented. By combining total RNA sequencing with sequence and structural-based homology detection, we identified 18 novel RNA viruses in cultured samples from two major groups of microbial algae: the chlorophytes and the chlorarachniophytes. Most of the RNA viruses identified in the green algae class Ulvophyceae were related to the Tombusviridae and Amalgaviridae viral families commonly associated with land plants. This suggests that the evolutionary history of these viruses extends to divergence events between algae and land plants. Seven Ostreobium sp-associated viruses exhibited sequence similarity to the mitoviruses most commonly found in fungi, compatible with horizontal virus transfer between algae and fungi. We also document, for the first time, RNA viruses associated with chlorarachniophytes, including the first negative-sense (bunya-like) RNA virus in microalgae, as well as a distant homolog of the plant virus Virgaviridae, potentially signifying viral inheritance from the secondary chloroplast endosymbiosis that marked the origin of the chlorarachniophytes. More broadly, these data suggest that the scarcity of RNA viruses in algae results from limited investigation rather than their absence.

Metatranscriptomic Identification of Diverse and Divergent RNA Viruses in Green and Chlorarachniophyte Algae Cultures

Our knowledge of the diversity and evolution of the virosphere will likely increase dramatically with the study of microbial eukaryotes, including the microalgae within which few RNA viruses have been documented. By combining total RNA sequencing with sequence and structural-based homology detection, we identified 18 novel RNA viruses in cultured samples from two major groups of microbial algae: the chlorophytes and the chlorarachniophytes. Most of the RNA viruses identified in the green algae class Ulvophyceae were related to the Tombusviridae and Amalgaviridae viral families commonly associated with land plants. This suggests that the evolutionary history of these viruses extends to divergence events between algae and land plants. Seven Ostreobium sp-associated viruses exhibited sequence similarity to the mitoviruses most commonly found in fungi, compatible with horizontal virus transfer between algae and fungi. We also document, for the first time, RNA viruses associated with chlorarachniophytes, including the first negative-sense (bunya-like) RNA virus in microalgae, as well as a distant homolog of the plant virus Virgaviridae, potentially signifying viral inheritance from the secondary chloroplast endosymbiosis that marked the origin of the chlorarachniophytes. More broadly, these data suggest that the scarcity of RNA viruses in algae results from limited investigation rather than their absence.