Virus-derived small RNAs in the penaeid shrimp Fenneropenaeus chinensis during acute infection of the DNA virus WSSV

vsRNAfinder: a novel method for identifying high-confidence viral small RNAs from small RNA-Seq data

Virus-encoded small RNAs (vsRNA) have been reported to play an important role in viral infection. Unfortunately, there is still a lack of an effective method for vsRNA identification. Herein, we presented vsRNAfinder, a de novo method for identifying high-confidence vsRNAs from small RNA-Seq (sRNA-Seq) data based on peak calling and Poisson distribution and is publicly available at https://github.com/ZenaCai/vsRNAfinder. vsRNAfinder outperformed two widely used methods namely miRDeep2 and ShortStack in identifying viral miRNAs with a significantly improved sensitivity. It can also be used to identify sRNAs in animals and plants with similar performance to miRDeep2 and ShortStack. vsRNAfinder would greatly facilitate effective identification of vsRNAs from sRNA-Seq data.

A Method to Produce vsiRNAs in Plants with Cross-Kingdom Gene Silencing Capacity

Plants have evolved defense mechanisms to suppress viral transcription and replication by transcriptional and post-transcriptional gene silencing mediated by virus-derived small interfering RNAs (vsiRNAs). Based on this response, virus-induced gene silencing (VIGS)-based technology has been developed to silence target genes on either host plants or insect pests. This mechanism could also be used for the silencing of genes of interest in the medical field. We used the VIGS vector pEuMV-YP:Krt18, which was obtained by inserting the Mus musculus (M. musculus) Krt18 sequence into pEuMV-YP:ΔAV1. The objective was to evaluate the capacity of pEuMV-YP:Krt18 to induce Nicotiana benthamiana (N. benthamiana) production of vsiRNAs of a specific sequence that belongs to neither the plant genome nor the wild virus genome, which were used to induce cross-kingdom gene silencing between plants and mammals. The percentage of vsiRNA for each viral gene was calculated from an sRNA library of N. benthamiana plants infected by pEuMV-YP: Krt18. When the vsiRNAs were characterized, it was found that they corresponded to all the genes of the pEuMV-YP:Krt18 vector. These vsiRNAs induced the silencing of the Krt18 gene in M. musculus macrophages, supporting the ability to use VIGS vectors in plants as biofactories for the production of sRNAs that induce gene silencing in mammals.

Genome-wide transcriptomic analysis of highly virulent African swine fever virus infection reveals complex and unique virus host interaction

African swine fever virus (ASFV), one of the most devastating emerging swine pathogens in China, causes nearly 100 % mortality in naive herds. Here, whole-transcriptome RNA-seq analysis was conducted in porcine alveolar macrophages (PAMs) infected with Pig/Heilongjiang/2018 (Pig/HLJ/18) ASFV at different time points. Our data suggested that ASFV genes expression demonstrated a time-depended pattern and ASFV early genes were involved in antagonizing host innate immunity. Moreover, viral small RNA (vsRNA) was generated as well. Meanwhile, transcriptome analysis of host genes suggested a strong inhibition host immunity-related genes by ASFV infection in PAMs, while enhanced chemokine-mediated signaling pathways and neutrophil chemotaxis were observed in ASFV infected PAMs. Furthermore, ASFV infection also down-regulated host microRNAs (miRNAs) that putatively targeted viral genes, while also triggering dysregulation of host metabolism that promoted virus replication at transcription level. Most importantly, infection of PAMs with ASFV induced a different transcriptome pattern from that of highly pathogenic porcine reproductive and respiratory syndrome virus (HP-PRRSV), which is known to trigger a host cytokine storm. In conclusion, our transcriptome data implied that ASFV infection in PAMs appeared to be associated with strong inhibition of host immune responses, dysregulation of host chemokine axis and metabolic pathways.

Global mRNA and miRNA Analysis Reveal Key Processes in the Initial Response to Infection with WSSV in the Pacific Whiteleg Shrimp

White Spot Disease (WSD) presents a major barrier to penaeid shrimp production. Mechanisms underlying White Spot Syndrome Virus (WSSV) susceptibility in penaeids are poorly understood due to limited information related to early infection. We investigated mRNA and miRNA transcription in Penaeus vannamei over 36 h following infection. Over this time course, 6192 transcripts and 27 miRNAs were differentially expressed—with limited differential expression from 3–12 h post injection (hpi) and a more significant transcriptional response associated with the onset of disease symptoms (24 hpi). During early infection, regulated processes included cytoskeletal remodelling and alterations in phagocytic activity that may assist WSSV entry and translocation, novel miRNA-induced metabolic shifts, and the downregulation of ATP-dependent proton transporter subunits that may impair cellular recycling. During later infection, uncoupling of the electron transport chain may drive cellular dysfunction and lead to high mortalities in infected penaeids. We propose that post-transcriptional silencing of the immune priming gene Dscam (downregulated following infections) by a novel shrimp miRNA (Pva-pmiR-78; upregulated) as a potential mechanism preventing future recognition of WSSV that may be suppressed in surviving shrimp. Our findings improve our understanding of WSD pathogenesis in P. vannamei and provide potential avenues for future development of prophylactics and treatments.

A genome-wide screening for RNAi pathway proteins in Acari

BACKGROUND

RNA interference (RNAi) is a highly conserved, sequence-specific gene silencing mechanism present in Eukaryotes. Three RNAi pathways are known, namely micro-RNA (miRNA), piwi-interacting RNA (piRNA) and short interfering RNA (siRNA). However, little knowledge exists about the proteins involved in these pathways in Acari. Moreover, variable successes has been obtained in gene knockdown via siRNA pathway in their functional genomics and management. We hypothesized that the clue may be in the variability of the composition and the efficacy of siRNA machinery among Acari.

RESULTS

Both comparative genomic analyses and domain annotation suggest that all the analyzed species have homologs of putative core proteins that mediate cleaving of targeted genes via the three RNAi pathways. We identified putative homologs of Caenorhabditis elegans RNA-dependent RNA polymerase (RdRP) protein in all species though no secondary Argonaute homologs that operate with this protein in siRNA amplification mechanism were found, suggesting that the siRNA amplification mechanism present in Acari may be distinct from that described in C. elegans. Moreover, the genomes of these species do not encode homologs of C. elegans systemic RNAi defective-1 (Sid-1) protein that mediate silencing of the mRNA target throughout the treated organisms suggesting that the phenomena of systemic RNAi that has been reported in some Acari species probably occur through a different mechanism. However, homologs of putative RNAi spreading defective-3 (Rsd-3) protein and scavenger receptors namely Eater and SR-CI that mediate endocytosis cellular update of dsRNA in C. elegans and Drosophila melanogaster were found in Acari genomes. This result suggests that cellular dsRNA uptake in Acari is endocytosis-dependent. Detailed phylogenetic analyses of core RNAi pathway proteins in the studied species revealed that their evolution is compatible with the proposed monophyletic evolution of this group.

CONCLUSIONS

Our analyses have revealed the potential activity of all three pathways in Acari. Still, much experimental work remains to be done to confirm the mechanisms behind these pathways in particular those that govern systemic/parental RNAi and siRNA amplification in Acari. Disclosure of these mechanisms will facilitate the development of new and specific management tools for the harmful species and enrichment of the beneficial species.

Comparative analysis of shrimp (Penaeus vannamei) miRNAs expression profiles during WSSV infection under experimental conditions and in pond culture

In recent years, the importance of viral and host microRNAs (miRNAs) in mediating viral replication and control of host cellular machinery, has been realised and increasing efforts have been taken in order to understand the interactions of miRNAs from host and pathogen during infection. However, all existing studies has thus far been conducted in controlled experimental conditions and the veracity of these data for field conditions are yet to be established. In this framework, small RNA sequencing was performed to identify the miRNAs involved in shrimp (Penaeus vannamei) immune responses under two different WSSV infection conditions of natural infection and experimentally challenged conditions. The expression profiles of miRNAs of shrimp infected with WSSV under two contrasting conditions were compared and as a result, 23365 known miRNAs and 481 novel miRNAs were identified. Amongst the most abundantly expressed miRNAs, the hypoxia related miR-210 and immune pathway related miR-29b were expressed only in infected shrimps of both conditions. miR-8-5p, having a functional role in modulation of chitin biosynthesis was exclusively represented in higher numbers in the WSSV -infected shrimps under natural conditions whilst four of the miRNAs (mja-miR-6493-5p, mja-miR-6492, mmu-miR-3968, tcf-miR-9b-5p) identified from shrimps collected from pond culture targeted chitinase, an important enzyme involved in growth and moulting in shrimps, indicating an interaction between WSSV infection and moult cycle under culture conditions. Some of the miRNAs (tca-miR-87b-3p, cte-miR-277a) and miRNAs belonging to class miR-9, miR-981 that were identified only in WSSV infected shrimps under experimental conditions, are known to respond against WSSV infection in shrimps. Moreover, the miRNA target prediction revealed several immune-related gene targets such as cathepsin, c-type lectin, haemocyanin and ubiquitin protein ligase were commonly identified under both the conditions. However, the miRNAs identified from challenge experiment had wide number of gene targets as compared to the miRNAs of natural infection. The shrimp miRNA mja-miR-6489-3p, was also found to target early virus gene wsv001 of WSSV. Our study, therefore, provides the comparative analysis of miRNA expression from shrimp during WSSV infection in two different conditions.

Discovery of Novel Crustacean and Cephalopod Flaviviruses: Insights into the Evolution and Circulation of Flaviviruses between Marine Invertebrate and Vertebrate Hosts

Most described flaviviruses (family Flaviviridae) are disease-causing pathogens of vertebrates maintained in zoonotic cycles between mosquitoes or ticks and vertebrate hosts. Poor sampling of flaviviruses outside vector-borne flaviviruses such as Zika virus and dengue virus has presented a narrow understanding of flavivirus diversity and evolution. In this study, we discovered three crustacean flaviviruses (Gammarus chevreuxi flavivirus, Gammarus pulex flavivirus, and Crangon crangon flavivirus) and two cephalopod flaviviruses (Southern Pygmy squid flavivirus and Firefly squid flavivirus). Bayesian and maximum likelihood phylogenetic methods demonstrate that crustacean flaviviruses form a well-supported clade and share a more closely related ancestor with terrestrial vector-borne flaviviruses than with classical insect-specific flaviviruses. In addition, we identify variants of Wenzhou shark flavivirus in multiple gazami crab (Portunus trituberculatus) populations, with active replication supported by evidence of an active RNA interference response. This suggests that Wenzhou shark flavivirus moves horizontally between sharks and gazami crabs in ocean ecosystems. Analyses of the mono- and dinucleotide composition of marine flaviviruses compared to that of flaviviruses with known host status suggest that some marine flaviviruses share a nucleotide bias similar to that of vector-borne flaviviruses. Furthermore, we identify crustacean flavivirus endogenous viral elements that are closely related to elements of terrestrial vector-borne flaviviruses. Taken together, these data provide evidence of flaviviruses circulating between marine vertebrates and invertebrates, expand our understanding of flavivirus host range, and offer potential insights into the evolution and emergence of terrestrial vector-borne flaviviruses.IMPORTANCE Some flaviviruses are known to cause disease in vertebrates and are typically transmitted by blood-feeding arthropods such as ticks and mosquitoes. While an ever-increasing number of insect-specific flaviviruses have been described, we have a narrow understanding of flavivirus incidence and evolution. To expand this understanding, we discovered a number of novel flaviviruses that infect a range of crustaceans and cephalopod hosts. Phylogenetic analyses of these novel marine flaviviruses suggest that crustacean flaviviruses share a close ancestor to all terrestrial vector-borne flaviviruses, and squid flaviviruses are the most divergent of all known flaviviruses to date. Additionally, our results indicate horizontal transmission of a marine flavivirus between crabs and sharks. Taken together, these data suggest that flaviviruses move horizontally between invertebrates and vertebrates in ocean ecosystems. This study demonstrates that flavivirus invertebrate-vertebrate host associations have arisen in flaviviruses at least twice and may potentially provide insights into the emergence or origin of terrestrial vector-borne flaviviruses.

Dual Analysis of Virus-Host Interactions: The Case of Ostreid herpesvirus 1 and the Cupped Oyster Crassostrea gigas

Dual analyses of the interactions between Ostreid herpesvirus 1 (OsHV-1) and the bivalve Crassostrea gigas during infection can unveil events critical to the onset and progression of this viral disease and can provide novel strategies for mitigating and preventing oyster mortality. Among the currently used “omics” technologies, dual transcriptomics (dual RNA-seq) coupled with the analysis of viral DNA in the host tissues has greatly advanced the knowledge of genes and pathways mostly contributing to host defense responses, expression profiles of annotated and unknown OsHV-1 open reading frames (ORFs), and viral genome variability. In addition to dual RNA-seq, proteomics and metabolomics analyses have the potential to add complementary information, needed to understand how a malacoherpesvirus can redirect and exploit the vital processes of its host. This review explores our current knowledge of “omics” technologies in the study of host-pathogen interactions and highlights relevant applications of these fields of expertise to the complex case of C gigas infections by OsHV-1, which currently threaten the mollusk production sector worldwide.

Aquatic animal viruses mediated immune evasion in their host

Viruses are important and lethal pathogens that hamper aquatic animals. The result of the battle between host and virus would determine the occurrence of diseases. The host will fight against virus infection with various responses such as innate immunity, adaptive immunity, apoptosis, and so on. On the other hand, the virus also develops numerous strategies such as immune evasion to antagonize host antiviral responses. Here, We review the research advances on virus mediated immune evasions to host responses containing interferon response, NF-κB signaling, apoptosis, and adaptive response, which are executed by viral genes, proteins, and miRNAs from different aquatic animal viruses including Alloherpesviridae, Iridoviridae, Nimaviridae, Birnaviridae, Reoviridae, and Rhabdoviridae. Thus, it will facilitate the understanding of aquatic animal virus mediated immune evasion and potentially benefit the development of novel antiviral applications.

Changes of microRNAs expression profiles from red swamp crayfish (Procambarus clarkia) hemolymph exosomes in response to WSSV infection

MicroRNAs (miRNAs) as short noncoding RNAs play important regulatory roles in diverse biological processes by degrading the target mRNAs, and could be delivered by exosomes. WSSV is a highly pathogenic and prevalent virus, and has brought high mortality of P. clarkia. Till present, no studies focus on the miRNAs changes in exosomes during WSSV infection. To understand the different virulence of WSSV on miRNAs expression in P. clarkia hemolymph exosome, the deep sequencing was performed to compare the small RNA libraries from the hemolymph exosome of P. clarkia individuals with or without WSSV infections. From the TEM observations, NTA and Western Blot analysis, the extracted exosomes were well identified with classic characteristics. The 209 conserved miRNAs and 250 novel miRNAs were identified from the small RNA libraries. In response to WSSV infection, there were about 98 miRNAs significantly up-regulated and 59 miRNAs significantly down-regulated. The target genes prediction, GO and KEGG enrichment analysis revealed that some target genes of P. clarkia miRNAs were grouped mainly into the categories of biological regulation, immune system process, signal pathway and other more functions. This is the first report of comprehensive identification of P. clarkia hemolymph exosome miRNAs being differentially regulated in response to WSSV infection. These results will help to understand the hemolymph exosome miRNAs response to different virulence WSSV infection.

High-throughput sequencing analysis of microRNAs in gills of red swamp crayfish, Procambarus clarkii infected with white spot syndrome virus

MicroRNAs (miRNAs) are important posttranscriptional regulators. They play an important role in the antiviral innate immunity of invertebrates. In the present study, high-throughput small RNAs Illumina sequencing systems were carried out to identify differentially expressed miRNAs (DEMs) in the gills of Procambarus clarkii, which was challenged with white spot syndrome virus (WSSV). Our results identified 11,617 known and 6 novel miRNAs in normal group (NG) and WSSV-challenged group (WG) small RNA libraries. Additionally, 27 DEMs were shown to participate in the antiviral innate immunity of P. clarkii and were significantly upregulated or downregulated. In addition, the results of the KEGG pathway prediction of the DEMs target genes showed that putative target genes of these 27 DEMs were related mainly to the RNA transport pathway, tight junction pathway, mRNA surveillance pathway, regulation actin cytoskeleton pathway, focal adhesion pathway, and MAPK signaling pathway. These results provide important information for future studies about the antiviral innate immunity of crustaceans.

Identification of MicroRNAs and Their Target Genes Associated with Ovarian Development in Black Tiger Shrimp (Penaeus monodon) Using High-Throughput Sequencing

Plenty of evidence showing that microRNAs (miRNAs) post-transcriptionally regulate gene expression and are involved in a wide range of biological processes. However, the roles of miRNAs in ovarian development process remain largely unknown in shrimp. In the present study, high-throughput sequencing of small RNAs was performed to find specific miRNAs that are involved in ovarian development process in Penaeus monodon. Two small RNA libraries were constructed from undeveloped (UNDEV group) and developed (DEV group) ovarian tissues in P. monodon. In total, 43 differentially expressed miRNAs were identified between the two groups (P ≤ 0.05, |log2 ratio| ≥1), and their expression profiles were validated by qRT-PCR. In order to further clarify the functional roles of these differentially expressed miRNAs during ovarian development process, target gene prediction was performed. In total, 4,102 target genes of 43 miRNAs were predicted, then clustered by the Kyoto Encyclopedia of Genes and Genomes (KEGG) database; only four specific pathways related to ovarian development were obtained (P < 0.05). Dual-luciferase reporter assays and integrated expression analysis were also conducted to further clarify the interaction between the miRNAs and their target mRNAs. This study provides important information about the function of miRNAs involved in ovarian developmental stages in P. monodon.

Identification and characterization of intestine microRNAs and targets in red swamp crayfish, Procambarus clarkii infected with white spot syndrome virus

MicroRNAs (miRNAs) are small non-coding endogenous RNA molecules that play important roles in the innate immunity system of invertebrates, especially in the aspect of antivirus. In the present study, high-throughput small RNA Illumina sequencing systems were used to identify differentially expressed miRNAs (DEMs) from the intestines of Procambarus clarkii that were infected with white spot syndrome virus (WSSV). As a result, 39 known and 12 novel miRNAs were identified in both NG and WG small RNA libraries. Seven DEMs were determined to be involved in the antiviral innate immunity in the intestines of P. clarkii. The results of the target gene predictions of the DEMs showed that the putative target genes of these 7 DEMs are related to tight junctions, vascular smooth muscle contraction regulation of the actin cytoskeleton, focal adhesion, RNA transport, mRNA surveillance, viral carcinogenesis, and Salmonella infection. These results provide theoretical insights for future studies on the antiviral immunity of crustaceans.

Identification and characterization of lymph organ microRNAs in red swamp crayfish, Procambarus clarkii infected with white spot syndrome virus

MicroRNAs (miRNAs) were important post-transcriptional regulators and played vital roles in innate immunity system of invertebrates, especially in the aspect of antivirus. In this study, using high-throughput small RNAs Illumina sequencing system, differentially expressed miRNAs (DEMs) from lymph organs in red swamp crayfish, Procambarus clarkii, infected with white spot syndrome virus, were identified. As a result, 32 known miRNAs and 7 novel miRNAs were identified in crayfish lymph organ small RNAs library of NG and WG. Among them, 7 differentially expressed miRNAs (DEMs) were predicted to be involved in the lymph organ antiviral innate immunity of P. clarkii. Besides, the results showed that putative target genes of these DEMs were related with tight junction, RNA transport, regulation of actin cytoskeleton, focal adhesion, vascular smooth muscle contraction, mRNA surveillance pathway, NOD-like receptor signaling pathway, leukocyte transendothelial migration, and protein processing in endoplasmic reticulum. These results might provide the guiding theoretical foundation for future studies about crustaceans' antiviral innate immunity.

White spot syndrome virus VP28 specific double-stranded RNA provides protection through a highly focused siRNA population

Several studies have demonstrated that injection of double-stranded RNAs (dsRNA) homologous to mRNA for the white spot syndrome virus (WSSV) viral protein 28 (VP28) can induce protection in shrimp against WSSV through RNA interference (RNAi). In comparison to shrimp injected with either PBS or a green fluorescent protein (GFP) nonspecific dsRNA, we obtained nearly complete protection against WSSV infection in shrimp injected with VP28 dsRNA. Upregulation of host genes associated with small RNA silencing was measured 48 hours post treatment in groups injected with dsRNA, and although the VP28-treated group remained moderately upregulated after challenge with WSSV, many-fold higher induction was observed in both control groups reflecting the ongoing viral infection. RNA sequencing of VP28-treated shrimp demonstrated a siRNA population dominated by high levels of 22 nt long molecules narrowly targeting the VP28 mRNA both before and after challenge with WSSV. Conversely, while no siRNAs targeting WSSV were detected before challenge, a broad response of 22 nt siRNAs mapping across the entire WSSV genome were found in both control groups after challenge. These results give detailed insight to how dsRNA targeting VP28 function to induce protection against WSSV, by generating a highly focused population of 22 nt long siRNA molecules.

Advances in the application of high-throughput sequencing in invertebrate virology

Over the last decade, advances in high-throughput sequencing technologies have revolutionised biological research, making it possible for DNA/RNA sequencing of any organism of interest to be undertaken. Sequencing approaches are now routinely used in the detection and characterisation of (novel) viruses, investigation of host-pathogen interactions, and effective development of disease treatment strategies. For the sequencing and identification of viruses of interest, metagenomics approaches using infected host tissue are frequently used, as it is not always possible to culture and isolate these pathogens. High-throughput sequencing can also be used to investigate host-pathogen interactions by investigating (temporal) transcriptomic responses of both the host and virus, potentially leading to the discovery of novel opportunities for treatment and drug targets. In addition, viruses in environmental samples (e.g. water or soil samples) can be identified using eDNA/metagenomics approaches. The promise that recent developments in sequencing brings to the field of invertebrate virology are not devoid of technical challenges, including the need for better laboratory and bioinformatics strategies to sequence and assemble virus genomes within complex tissue or environmental samples, and the difficulties associated with the annotation of the large number of novel viruses being discovered.

Differential expression of microRNAs of Litopenaeus vannamei in response to different virulence WSSV infection

WSSV is one of the most harmful pathogeny in the pacific white shrimp, and genetic variations caused the strains of different virulence. MicroRNAs (miRNAs) involved in the regulation of virus defense. To understand the different virulence of WSSV on miRNA expression in Litopeneaus vannamei, the deep sequencing was performed to compare two small RNA libraries prepared from hepatopancreas of Litopeneaus vannamei infected with normal-virulence or low-virulence WSSV. Approximately 29,398,623 raw reads from normal-virulence library and 35,291,803 raw reads from low-virulence library were obtained. There were about 37 miRNAs homologs identified. Sixteen miRNAs were significantly up-regulated and twenty-one miRNAs were significantly down-regulated in normal-virulence infection library compared with low-virulence infection library. Of these, Igi-miR-1175-3p was the most significant different miRNA, followed by bmo-miR-1175-3p and ipu-miR-26b, respectively. The putative target genes for differentially expressed miRNAs were concerned with biological processes, signal meditated, cell differentiation and apoptosis, immune recognition and other more functions. The results will help to understand the miRNAs response to different virulence WSSV infection.