Detection and analysis of non-retroviral RNA virus-like elements in plant, fungal, and insect genomes

Mycovirus Diversity and Evolution Revealed/Inferred from Recent Studies

High-throughput virome analyses with various fungi, from cultured or uncultured sources, have led to the discovery of diverse viruses with unique genome structures and even neo-lifestyles. Examples in the former category include splipalmiviruses and ambiviruses. Splipalmiviruses, related to yeast narnaviruses, have multiple positive-sense (+) single-stranded (ss) RNA genomic segments that separately encode the RNA-dependent RNA polymerase motifs, the hallmark of RNA viruses (members of the kingdom Orthornavirae). Ambiviruses appear to have an undivided ssRNA genome of 3∼5 kb with two large open reading frames (ORFs) separated by intergenic regions. Another narna-like virus group has two fully overlapping ORFs on both strands of a genomic segment that span more than 90% of the genome size. New virus lifestyles exhibited by mycoviruses include the yado-kari/yado-nushi nature characterized by the partnership between the (+)ssRNA yadokarivirus and an unrelated dsRNA virus (donor of the capsid for the former) and the hadaka nature of capsidless 10-11 segmented (+)ssRNA accessible by RNase in infected mycelial homogenates. Furthermore, dsRNA polymycoviruses with phylogenetic affinity to (+)ssRNA animal caliciviruses have been shown to be infectious as dsRNA-protein complexes or deproteinized naked dsRNA. Many previous phylogenetic gaps have been filled by recently discovered fungal and other viruses, which haveprovided interesting evolutionary insights. Phylogenetic analyses and the discovery of natural and experimental cross-kingdom infections suggest that horizontal virus transfer may have occurred and continue to occur between fungi and other kingdoms.

Endogenous reverse transcriptase and RNase H-mediated antiviral mechanism in embryonic stem cells

Nucleic acid-based systems play important roles in antiviral defense, including CRISPR/Cas that adopts RNA-guided DNA cleavage to prevent DNA phage infection and RNA interference (RNAi) that employs RNA-guided RNA cleavage to defend against RNA virus infection. Here, we report a novel type of nucleic acid-based antiviral system that exists in mouse embryonic stem cells (mESCs), which suppresses RNA virus infection by DNA-mediated RNA cleavage. We found that the viral RNA of encephalomyocarditis virus can be reverse transcribed into complementary DNA (vcDNA) by the reverse transcriptase (RTase) encoded by endogenous retrovirus-like elements in mESCs. The vcDNA is negative-sense single-stranded and forms DNA/RNA hybrid with viral RNA. The viral RNA in the heteroduplex is subsequently destroyed by cellular RNase H1, leading to robust suppression of viral growth. Furthermore, either inhibition of the RTase activity or depletion of endogenous RNase H1 results in the promotion of virus proliferation. Altogether, our results provide intriguing insights into the antiviral mechanism of mESCs and the antiviral function of endogenized retroviruses and cellular RNase H. Such a natural nucleic acid-based antiviral mechanism in mESCs is referred to as ERASE (endogenous RTase/RNase H-mediated antiviral system), which is an addition to the previously known nucleic acid-based antiviral mechanisms including CRISPR/Cas in bacteria and RNAi in plants and invertebrates.

Diverse Partitiviruses From the Phytopathogenic Fungus, Rosellinia necatrix

Partitiviruses (dsRNA viruses, family Partitiviridae) are ubiquitously detected in plants and fungi. Although previous surveys suggested their omnipresence in the white root rot fungus, Rosellinia necatrix, only a few of them have been molecularly and biologically characterized thus far. We report the characterization of a total of 20 partitiviruses from 16 R. necatrix strains belonging to 15 new species, for which “Rosellinia necatrix partitivirus 11–Rosellinia necatrix partitivirus 25” were proposed, and 5 previously reported species. The newly identified partitiviruses have been taxonomically placed in two genera, Alphapartitivirus, and Betapartitivirus. Some partitiviruses were transfected into reference strains of the natural host, R. necatrix, and an experimental host, Cryphonectria parasitica, using purified virions. A comparative analysis of resultant transfectants revealed interesting differences and similarities between the RNA accumulation and symptom induction patterns of R. necatrix and C. parasitica. Other interesting findings include the identification of a probable reassortment event and a quintuple partitivirus infection of a single fungal strain. These combined results provide a foundation for further studies aimed at elucidating mechanisms that underly the differences observed.

Virome Analysis of Aphid Populations That Infest the Barley Field: The Discovery of Two Novel Groups of Nege/Kita-Like Viruses and Other Novel RNA Viruses

Aphids (order Hemiptera) are important insect pests of crops and are also vectors of many plant viruses. However, little is known about aphid-infecting viruses, particularly their diversity and relationship to plant viruses. To investigate the aphid viromes, we performed deep sequencing analyses of the aphid transcriptomes from infested barley plants in a field in Japan. We discovered virus-like sequences related to nege/kita-, flavi-, tombus-, phenui-, mononega-, narna-, chryso-, partiti-, and luteoviruses. Using RT-PCR and sequence analyses, we determined almost complete sequences of seven nege/kitavirus-like virus genomes; one of which was a variant of the Wuhan house centipede virus (WHCV-1). The other six seem to belong to four novel viruses distantly related to Wuhan insect virus 9 (WhIV-9) or Hubei nege-like virus 4 (HVLV-4). We designated the four viruses as barley aphid RNA virus 1 to 4 (BARV-1 to -4). Moreover, some nege/kitavirus-like sequences were found by searches on the transcriptome shotgun assembly (TSA) libraries of arthropods and plants. Phylogenetic analyses showed that BARV-1 forms a clade with WHCV-1 and HVLV-4, whereas BARV-2 to -4 clustered with WhIV-9 and an aphid virus, Aphis glycines virus 3. Both virus groups (tentatively designated as Centivirus and Aphiglyvirus, respectively), together with arthropod virus-like TSAs, fill the phylogenetic gaps between the negeviruses and kitaviruses lineages. We also characterized the flavi/jingmen-like and tombus-like virus sequences as well as other RNA viruses, including six putative novel viruses, designated as barley aphid RNA viruses 5 to 10. Interestingly, we also discovered that some aphid-associated viruses, including nege/kita-like viruses, were present in different aphid species, raising a speculation that these viruses might be distributed across different aphid species with plants being the reservoirs. This study provides novel information on the diversity and spread of nege/kitavirus-related viruses and other RNA viruses that are associated with aphids.

Novel Victorivirus from a Pakistani Isolate of Alternaria alternata Lacking a Typical Translational Stop/Restart Sequence Signature

The family Totiviridae currently contains five genera Totivirus, Victorivirus, Leishmavirus, Trichomonasvirus, and Giardiavirus. Members in this family generally have a set of two-open reading frame (ORF) elements in their genome with the 5'-proximal ORF (ORF1) encoding a capsid protein (CP) and the 3'-proximal one (ORF2) for RNA-dependent RNA polymerase (RdRp). How the downstream open reading frames (ORFs) are expressed is genus-specific. All victoriviruses characterized thus far appear to use the stop/restart translation mechanism, allowing for the expression of two separate protein products from bicitronic genome-sized viral mRNA, while the totiviruses use a -1 ribosomal frame-shifting that leads to a fusion product of CP and RdRp. We report the biological and molecular characterization of a novel victorivirus termed Alternaria alternata victorivirus 1 (AalVV1) isolated from Alternaria alternata in Pakistan. The phylogenetic and molecular analyses showed AalVV1 to be distinct from previously reported victoriviruses. AalVV1 appears to have a sequence signature required for the -1 frame-shifting at the ORF1/2 junction region, rather than a stop/restart key mediator. By contrast, SDS-polyacrylamide gel electrophoresis and peptide mass fingerprinting analyses of purified virion preparations suggested the expression of two protein products, not a CP-RdRp fusion product. How these proteins are expressed is discussed in this study. Possible effects of infection by this virus were tested in two fungal species: A. alternata and RNA silencing proficient and deficient strains of Cryphonectria parasitica, a model filamentous fungus. AalVV1 showed symptomless infection in all of these fungal strains, even in the RNA silencing deficient C. parasitica strain.

Two novel fungal negative-strand RNA viruses related to mymonaviruses and phenuiviruses in the shiitake mushroom (Lentinula edodes)

There is still limited information on the diversity of (-)ssRNA viruses that infect fungi. Here, we have discovered two novel (-)ssRNA mycoviruses in the shiitake mushroom (Lentinula edodes). The first virus has a monopartite RNA genome and relates to that of mymonaviruses (Mononegavirales), especially to Hubei rhabdo-like virus 4 from arthropods and thus designated as Lentinula edodes negative-strand RNA virus 1. The second virus has a putative bipartite RNA genome and is related to the recently discovered bipartite or tripartite phenui-like viruses (Bunyavirales) associated with plants and ticks, and designated as Lentinula edodes negative-strand RNA virus 2 (LeNSRV2). LeNSRV2 is likely the first segmented (-)ssRNA virus known to infect fungi. Its smaller RNA segment encodes a putative nucleocapsid and a plant MP-like protein using a potential ambisense coding strategy. These findings enhance our understanding of the diversity, evolution and spread of (-)ssRNA viruses in fungi.

Identification of a Novel Hypovirulence-Inducing Hypovirus From Alternaria alternata

Mycoviruses are wide spread throughout almost all groups of fungi but only a small number of mycoviruses can attenuate the growth and virulence of their fungal hosts. Alternaria alternata is an ascomycete fungus that causes leaf spot diseases on various crop plants. In this study, we identified a novel ssRNA mycovirus infecting an A. alternata f. sp. mali strain isolated from an apple orchard in China. Sequence analyses revealed that this virus is related to hypoviruses, in particular to Wuhan insect virus 14, an unclassified hypovirus identified from insect meta-transcriptomics, as well as other hypoviruses belonging to the genus Hypovirus, and therefore this virus is designed as Alternaria alternata hypovirus 1 (AaHV1). The genome of AaHV1 contains a single large open-reading frame encoding a putative polyprotein (∼479 kDa) with a cysteine proteinase-like and replication-associated domains. Curing AaHV1 from the fungal host strain indicated that the virus is responsible for the slow growth and reduced virulence of the host. AaHV1 defective RNA (D-RNA) with internal deletions emerging during fungal subcultures but the presence of D-RNA does not affect AaHV1 accumulation and pathogenicities. Moreover, AaHV1 could replicate and confer hypovirulence in Botryosphaeria dothidea, a fungal pathogen of apple white rot disease. This finding could facilitate better understanding of A. alternata pathogenicity and is relevant for development of biocontrol methods of fungal diseases.

Protein Structure-Guided Hidden Markov Models (HMMs) as A Powerful Method in the Detection of Ancestral Endogenous Viral Elements

It has been believed for a long time that the transfer and fixation of genetic material from RNA viruses to eukaryote genomes is very unlikely. However, during the last decade, there have been several cases in which “virus-to-host” gene transfer from various viral families into various eukaryotic phyla have been described. These transfers have been identified by sequence similarity, which may disappear very quickly, especially in the case of RNA viruses. However, compared to sequences, protein structure is known to be more conserved. Applying protein structure-guided protein domain-specific Hidden Markov Models, we detected homologues of the Virgaviridae capsid protein in Schizophora flies. Further data analysis supported “virus-to-host” transfer into Schizophora ancestors as a single transfer event. This transfer was not identifiable by BLAST or by other methods we applied. Our data show that structure-guided Hidden Markov Models should be used to detect ancestral virus-to-host transfers.

Biological and Molecular Characterization of Chenopodium quinoa Mitovirus 1 Reveals a Distinct Small RNA Response Compared to Those of Cytoplasmic RNA Viruses

Indirect evidence of mitochondrial viruses in plants comes from discovery of genomic fragments integrated into the nuclear and mitochondrial DNA of a number of plant species. Here, we report the existence of replicating mitochondrial virus in plants: from transcriptome sequencing (RNA-seq) data of infected Chenopodium quinoa, a plant species commonly used as a test plant in virus host range experiments, among other virus contigs, we could assemble a 2.7-kb contig that had highest similarity to mitoviruses found in plant genomes. Northern blot analyses confirmed the existence of plus- and minus-strand RNA corresponding to the mitovirus genome. No DNA corresponding to the genomic RNA was detected, excluding the endogenization of such virus. We have tested a number of C. quinoa accessions, and the virus was present in a number of commercial varieties but absent from a large collection of Bolivian and Peruvian accessions. The virus could not be transmitted mechanically or by grafting, but it is transmitted vertically through seeds at a 100% rate. Small RNA analysis of a C. quinoa line carrying the mitovirus and infected by alfalfa mosaic virus showed that the typical antiviral silencing response active against cytoplasmic viruses (21- to 22-nucleotide [nt] vsRNA peaks) is not active against CqMV1, since in this specific case the longest accumulating vsRNA length is 16 nt, which is the same as that corresponding to RNA from mitochondrial genes. This is evidence of a distinct viral RNA degradation mechanism active inside mitochondria that also may have an antiviral effect.IMPORTANCE This paper reports the first biological characterization of a bona fide plant mitovirus in an important crop, Chenopodium quinoa, providing data supporting that mitoviruses have the typical features of cryptic (persistent) plant viruses. We, for the first time, demonstrate that plant mitoviruses are associated with mitochondria in plants. In contrast to fungal mitoviruses, plant mitoviruses are not substantially affected by the antiviral silencing pathway, and the most abundant mitovirus small RNA length is 16 nt.

Dicer functions transcriptionally and posttranscriptionally in a multilayer antiviral defense

In antiviral RNA interference (RNAi), Dicer plays a primary role in processing double-stranded RNA (dsRNA) molecules into small-interfering RNAs (siRNAs) that guide Argonaute effectors to posttranscriptional suppression of target viral genes. Here, we show a distinct role for Dicer in the siRNA-independent transcriptional induction of certain host genes upon viral infection in a filamentous fungus. Previous studies have shown that the two key players, dicer-like 2 (dcl2) and argonaute-like 2 (agl2), of antiviral RNAi in a phytopathogenic ascomycete, Cryphonectria parasitica, are highly transcriptionally induced upon infection with certain RNA mycoviruses, including the positive-stranded RNA hypovirus mutant lacking the RNAi suppressor (Cryphonectria hypovirus 1-Δp69, CHV1-Δp69). This induction is regulated by the Spt-Ada-Gcn5 acetyltransferase (SAGA) complex, a well-known transcriptional coactivator. The present study shows that diverse host genes, in addition to dcl2 and agl2, were up-regulated more than 10-fold by SAGA upon infection with CHV1-Δp69. Interestingly, DCL2, but not AGL2, was essential for SAGA-mediated global gene up-regulation. Moreover, deletion of certain virus-induced genes enhanced a CHV1-Δp69 symptom (growth rate) but not its accumulation. Constitutive, modest levels of dcl2 expression drastically reduced viral siRNA accumulation but were sufficient for full-scale up-regulation of host genes, suggesting that high induction of dcl2 and siRNA production are not essential for the transcriptional up-regulation function of DCL2. These data clearly demonstrate the dual functionality of DCL2: as a dsRNA-specific nuclease in posttranscriptional antiviral RNA silencing and as a key player in SAGA-mediated host gene induction, which independently represses viral replication and alleviates virus-induced symptom expression.

Cellular production of a counterfeit viral protein confers immunity to infection by a related virus

DNA copies of many non-retroviral RNA virus genes or portions thereof (NIRVs) are present in the nuclear genomes of many eukaryotes. These have often been preserved for millions of years of evolution, suggesting that they play an important cellular function. One possible function is resistance to infection by related viruses. In some cases, this appears to occur through the piRNA system, but in others by way of counterfeit viral proteins encoded by NIRVs. In the fungi, NIRVs may be as long as 1,400 uninterrupted codons. In one such case in the yeast Debaryomyces hansenii, one of these genes provides immunity to a related virus by virtue of expression of a counterfeit viral capsid protein, which interferes with assembly of viral capsids by negative complementation. The widespread occurrence of non-retroviral RNA virus genes in eukaryotes may reflect an underappreciated method of host resistance to infection. This work demonstrates for the first time that an endogenous host protein encoded by a gene that has been naturally acquired from a virus and fixed in a eukaryote can interfere with the replication of a related virus and do so by negative complementation.

Novel, diverse RNA viruses from Mediterranean isolates of the phytopathogenic fungus, Rosellinia necatrix: insights into evolutionary biology of fungal viruses

To reveal mycovirus diversity, we conducted a search of as-yet-unexplored Mediterranean isolates of the phytopathogenic ascomycete Rosellinia necatrix for virus infections. Of seventy-nine, eleven fungal isolates tested RNA virus-positive, with many showing coinfections, indicating a virus incidence of 14%, which is slightly lower than that (approximately 20%) previously reported for extensive surveys of over 1000 Japanese R. necatrix isolates. All viral sequences were fully or partially characterized by Sanger and next-generation sequencing. These sequences appear to represent isolates of various new species spanning at least 6 established or previously proposed families such as Partiti-, Hypo-, Megabirna-, Yado-kari-, Fusagra- and Fusarividae, as well as a newly proposed family, Megatotiviridae. This observation greatly expands the diversity of R. necatrix viruses, because no hypo-, fusagra- or megatotiviruses were previously reported from R. necatrix. The sequence analyses showed a rare horizontal gene transfer event of the 2A-like protease domain between a dsRNA (phlegivirus) and a positive-sense, single-stranded RNA virus (hypovirus). Moreover, many of the newly detected viruses showed the closest relation to viruses reported from fungi other than R. necatrix, such as Fusarium spp., which are sympatric to R. necatrix. These combined results imply horizontal virus transfer between these soil-inhabitant fungi.

The family Rhabdoviridae: mono- and bipartite negative-sense RNA viruses with diverse genome organization and common evolutionary origins

The family Rhabdoviridae consists of mostly enveloped, bullet-shaped or bacilliform viruses with a negative-sense, single-stranded RNA genome that infect vertebrates, invertebrates or plants. This ecological diversity is reflected by the diversity and complexity of their genomes. Five canonical structural protein genes are conserved in all rhabdoviruses, but may be overprinted, overlapped or interspersed with several novel and diverse accessory genes. This review gives an overview of the characteristics and diversity of rhabdoviruses, their taxonomic classification, replication mechanism, properties of classical rhabdoviruses such as rabies virus and rhabdoviruses with complex genomes, rhabdoviruses infecting aquatic species, and plant rhabdoviruses with both mono- and bipartite genomes.

Reprint of "Sequence and phylogenetic analyses of novel totivirus-like double-stranded RNAs from field-collected powdery mildew fungi"

The identification of mycoviruses contributes greatly to understanding of the diversity and evolutionary aspects of viruses. Powdery mildew fungi are important and widely studied obligate phytopathogenic agents, but there has been no report on mycoviruses infecting these fungi. In this study, we used a deep sequencing approach to analyze the double-stranded RNA (dsRNA) segments isolated from field-collected samples of powdery mildew fungus-infected red clover plants in Japan. Database searches identified the presence of at least ten totivirus (genus Totivirus)-like sequences, termed red clover powdery mildew-associated totiviruses (RPaTVs). The majority of these sequences shared moderate amino acid sequence identity with each other (<44%) and with other known totiviruses (<59%). Nine of these identified sequences (RPaTV1a, 1b and 2-8) resembled the genome of the prototype totivirus, Saccharomyces cerevisiae virus-L-A (ScV-L-A) in that they contained two overlapping open reading frames (ORFs) encoding a putative coat protein (CP) and an RNA dependent RNA polymerase (RdRp), while one sequence (RPaTV9) showed similarity to another totivirus, Ustilago maydis virus H1 (UmV-H1) that encodes a single polyprotein (CP-RdRp fusion). Similar to yeast totiviruses, each ScV-L-A-like RPaTV contains a -1 ribosomal frameshift site downstream of a predicted pseudoknot structure in the overlapping region of these ORFs, suggesting that the RdRp is translated as a CP-RdRp fusion. Moreover, several ScV-L-A-like sequences were also found by searches of the transcriptome shotgun assembly (TSA) libraries from rust fungi, plants and insects. Phylogenetic analyses show that nine ScV-L-A-like RPaTVs along with ScV-L-A-like sequences derived from TSA libraries are clustered with most established members of the genus Totivirus, while one RPaTV forms a new distinct clade with UmV-H1, possibly establishing an additional genus in the family. Taken together, our results indicate the presence of diverse, novel totiviruses in the powdery mildew fungus populations infecting red clover plants in the field.

Sequence and phylogenetic analyses of novel totivirus-like double-stranded RNAs from field-collected powdery mildew fungi

The identification of mycoviruses contributes greatly to understanding of the diversity and evolutionary aspects of viruses. Powdery mildew fungi are important and widely studied obligate phytopathogenic agents, but there has been no report on mycoviruses infecting these fungi. In this study, we used a deep sequencing approach to analyze the double-stranded RNA (dsRNA) segments isolated from field-collected samples of powdery mildew fungus-infected red clover plants in Japan. Database searches identified the presence of at least ten totivirus (genus Totivirus)-like sequences, termed red clover powdery mildew-associated totiviruses (RPaTVs). The majority of these sequences shared moderate amino acid sequence identity with each other (<44%) and with other known totiviruses (<59%). Nine of these identified sequences (RPaTV1a, 1b and 2-8) resembled the genome of the prototype totivirus, Saccharomyces cerevisiae virus-L-A (ScV-L-A) in that they contained two overlapping open reading frames (ORFs) encoding a putative coat protein (CP) and an RNA dependent RNA polymerase (RdRp), while one sequence (RPaTV9) showed similarity to another totivirus, Ustilago maydis virus H1 (UmV-H1) that encodes a single polyprotein (CP-RdRp fusion). Similar to yeast totiviruses, each ScV-L-A-like RPaTV contains a -1 ribosomal frameshift site downstream of a predicted pseudoknot structure in the overlapping region of these ORFs, suggesting that the RdRp is translated as a CP-RdRp fusion. Moreover, several ScV-L-A-like sequences were also found by searches of the transcriptome shotgun assembly (TSA) libraries from rust fungi, plants and insects. Phylogenetic analyses show that nine ScV-L-A-like RPaTVs along with ScV-L-A-like sequences derived from TSA libraries are clustered with most established members of the genus Totivirus, while one RPaTV forms a new distinct clade with UmV-H1, possibly establishing an additional genus in the family. Taken together, our results indicate the presence of diverse, novel totiviruses in the powdery mildew fungus populations infecting red clover plants in the field.

A novel betapartitivirus RnPV6 from Rosellinia necatrix tolerates host RNA silencing but is interfered by its defective RNAs

The family Partitiviridae comprises of five genera with bi-segmented dsRNA genomes that accommodate members infecting plants, fungi or protists. All partitiviruses with only a few exceptions cause asymptomatic infections. We report the characterization of a novel betapartitivirus termed Rosellinia necatrix partitivirus 6 (RnPV6) from a field isolate of a plant pathogenic fungus, white root rot fungus. RnPV6 has typical partitivirus features: dsRNA1 and dsRNA2 are 2462 and 2499bps in length encoding RNA-dependent RNA polymerase and capsid protein. Purified particles are spherical with a diameter of 30nm. Taking advantage of infectivity as virions, RnPV6 was introduced into a model filamentous fungal host, chestnut blight fungus to investigate virus/host interactions. Unlike other partitiviruses tested previously, RnPV6 induced profound phenotypic alterations with symptoms characterized by a reduced growth rate and enhanced pigmentation and was tolerant to host RNA silencing. In addition, a variety of defective RNAs derived from dsRNA1 appear after virion transfection. These sub-viral RNAs were shown to interfere with RnPV6 replication, at least for that of cognate segment dsRNA1. Presence of these sub-viral elements resulted in reduced symptom expression by RnPV6, suggesting their nature as defective-interfering RNAs. The features of RnPV6 are similar to but distinct from those of a previously reported alphapartitivirus, Rosellinia necatrix partitivirus 2 that is susceptible to RNA silencing.