Non-replicative Integral Membrane Proteins Encoded by Plant Alpha-Like Viruses: Emergence of Diverse Orphan ORFs and Movement Protein Genes

Endogenous nege-like viral elements in arthropod genomes reveal virus-host coevolution and ancient history of two plant virus families

Negevirus is a recently proposed taxon of arthropod-infecting virus, which is associated with plant viruses of two families (Virgaviridae and Kitaviridae). Nevertheless, the evolutionary history of negevirus-host and its relationship with plant viruses remain poorly understood. Endogenous nege-like viral elements (ENVEs) are ancient nege-like viral sequences integrated into the arthropod genomes, which can serve as the molecular fossil records of previous viral infection. In this study, 292 ENVEs were identified in 150 published arthropod genomes, revealing the evolutionary history of nege-like viruses and two related plant virus families. We discovered three novel and eight strains of nege-like viruses in 11 aphid species. Further analysis indicated that 10 ENVEs were detected in six aphid genomes, and they were divided into four types (ENVE1-ENVE4). Orthologous integration and phylogenetic analyses revealed that nege-like viruses had a history of infection of over 60 My and coexisted with aphid ancestors throughout the Cenozoic Era. Moreover, two nege-like viral proteins (CP and SP24) were highly homologous to those of plant viruses in the families Virgaviridae and Kitaviridae. CP- and SP24-derived ENVEs were widely integrated into numerous arthropod genomes. These results demonstrate that nege-like viruses have a long-term coexistence with arthropod hosts and plant viruses of the two families, Virgaviridae and Kitaviridae, which may have evolved from the nege-like virus ancestor through horizontal virus transfer events. These findings broaden our perspective on the history of viral infection in arthropods and the origins of plant viruses.

IMPORTANCE

Although negevirus is phylogenetically related to plant virus, the evolutionary history of negevirus-host and its relationship with plant virus remain largely unknown. In this study, we used endogenous nege-like viral elements (ENVEs) as the molecular fossil records to investigate the history of nege-like viral infection in arthropod hosts and the evolution of two related plant virus families (Virgaviridae and Kitaviridae). Our results showed the infection of nege-like viruses for over 60 My during the arthropod evolution. ENVEs highly homologous to viral sequences in Virgaviridae and Kitaviridae were present in a wide range of arthropod genomes but were absent in plant genomes, indicating that plant viruses in these two families possibly evolved from the nege-like virus ancestor through cross-species horizontal virus transmission. Our findings provide a new perspective on the virus-host coevolution and the origins of plant viruses.

Properties of Plant Virus Protein Encoded by the 5'-Proximal Gene of Tetra-Cistron Movement Block

To move from cell to cell through plasmodesmata, many plant viruses require the concerted action of two or more movement proteins (MPs) encoded by transport gene modules of virus genomes. A tetra-cistron movement block (TCMB) is a newly discovered transport module comprising four genes. TCMB encodes three proteins, which are similar to MPs of the transport module known as the "triple gene block", and a protein unrelated to known viral MPs and containing a double-stranded RNA (dsRNA)-binding domain similar to that found in a family of cell proteins, including AtDRB4 and AtHYL1. Here, the latter TCMB protein, named vDRB for virus dsRNA-binding protein, is shown to bind both dsRNA and single-stranded RNA in vitro. In a turnip crinkle virus-based assay, vDRB exhibits the properties of a viral suppressor of RNA silencing (VSR). In the context of potato virus X infection, vDRB significantly decreases the number and size of "dark green islands", regions of local antiviral silencing, supporting the VSR function of vDRB. Nevertheless, vDRB does not exhibit the VSR properties in non-viral transient expression assays. Taken together, the data presented here indicate that vDRB is an RNA-binding protein exhibiting VSR functions in the context of viral infection.

Predicted Membrane-Associated Domains in Proteins Encoded by Novel Monopartite Plant RNA Viruses Related to Members of the Family Benyviridae

As a continuation of our previous work, in this paper, we examine in greater detail the genome organization and some protein properties of the members of a potential group named Reclovirids and belonging to Benyviridae-related viruses. It can be proposed that the single-component Reclovirid genomes encode previously undiscovered transport genes. Indeed, analysis of the coding potential of these novel viral genomes reveals one or more cistrons ranging in size from 40 to 80 to about 600 codons, located in the 3'-terminal region of the genomic RNA, encoding proteins with predicted hydrophobic segments that are structurally diverse among Reclovirids and have no analogues in other plant RNA viruses. Additionally, in many cases, the possible methyltransferase domain of Reclovirid replicases is preceded by membrane-embedded protein segments that are not present in annotated members of the Benyviridae family. These observations suggest a general association of most Reclovirid proteins with cell membranes.

Interaction between Movement Proteins of Hibiscus green spot virus

Movement proteins (MPs) of plant viruses enable the translocation of viral genomes from infected to healthy cells through plasmodesmata (PD). The MPs functions involve the increase of the PD permeability and routing of viral genome both to the PD entrance and through the modified PD. Hibiscus green spot virus encodes two MPs, termed BMB1 and BMB2, which act in concert to accomplish virus cell-to-cell transport. BMB1, representing an NTPase/helicase domain-containing RNA-binding protein, localizes to the cytoplasm and the nucleoplasm. BMB2 is a small hydrophobic protein that interacts with the endoplasmic reticulum (ER) membranes and induces local constrictions of the ER tubules. In plant cells, BMB2 localizes to PD-associated membrane bodies (PAMBs) consisting of modified ER tubules and directs BMB1 to PAMBs. Here, we demonstrate that BMB1 and BMB2 interact in vitro and in vivo, and that their specific interaction is essential for BMB2-directed targeting of BMB1 to PAMBs. Using mutagenesis, we show that the interaction involves the C-terminal BMB1 region and the N-terminal region of BMB2.

Uncovering Plant Virus Species Forming Novel Provisional Taxonomic Units Related to the Family Benyviridae

Based on analyses of recent open-source data, this paper describes novel horizons in the diversity and taxonomy of beny-like viruses infecting hosts of the plant kingdom (Plantae or Archaeplastida). First, our data expand the known host range of the family Benyviridae to include red algae. Second, our phylogenetic analysis suggests that the evolution of this virus family may have involved cross-kingdom host change events and gene recombination/exchanges between distant taxa. Third, the identification of gene blocks encoding known movement proteins in beny-like RNA viruses infecting non-vascular plants confirms other evidence that plant virus genomic RNAs may have acquired movement proteins simultaneously or even prior to the evolutionary emergence of the plant vascular system. Fourth, novel data on plant virus diversity highlight that molecular evolution gave rise to numerous provisional species of land-plant-infecting viruses, which encode no known potential movement genetic systems.

Distinct Mechanisms of Endomembrane Reorganization Determine Dissimilar Transport Pathways in Plant RNA Viruses

Plant viruses exploit the endomembrane system of infected cells for their replication and cell-to-cell transport. The replication of viral RNA genomes occurs in the cytoplasm in association with reorganized endomembrane compartments induced by virus-encoded proteins and is coupled with the virus intercellular transport via plasmodesmata that connect neighboring cells in plant tissues. The transport of virus genomes to and through plasmodesmata requires virus-encoded movement proteins (MPs). Distantly related plant viruses encode different MP sets, or virus transport systems, which vary in the number of MPs and their properties, suggesting their functional differences. Here, we discuss two distinct virus transport pathways based on either the modification of the endoplasmic reticulum tubules or the formation of motile vesicles detached from the endoplasmic reticulum and targeted to endosomes. The viruses with the movement proteins encoded by the triple gene block exemplify the first, and the potyviral system is the example of the second type. These transport systems use unrelated mechanisms of endomembrane reorganization. We emphasize that the mode of virus interaction with cell endomembranes determines the mechanism of plant virus cell-to-cell transport.

Bioinformatic Analysis Predicts a Novel Genetic Module Related to Triple Gene and Binary Movement Blocks of Plant Viruses: Tetra-Cistron Movement Block

Previous studies have shown that the RNA genomes of some plant viruses encode two related genetic modules required for virus movement over the host body, containing two or three genes and named the binary movement block (BMB) and triple gene block (TGB), respectively. In this paper, we predict a novel putative-related movement gene module, called the tetra-cistron movement block (TCMB), in the virus-like transcriptome assemblies of the moss Dicranum scoparium and the Antarctic flowering plant Colobanthus quitensis. These TCMBs are encoded by smaller RNA components of putative two-component viruses related to plant benyviruses. Similar to the RNA2 of benyviruses, TCMB-containing RNAs have the 5′-terminal coat protein gene and include the RNA helicase gene which is followed by two small overlapping cistrons encoding hydrophobic proteins with a distant sequence similarity to the TGB2 and TGB3 proteins. Unlike TGB, TCMB also includes a fourth 5′-terminal gene preceding the helicase gene and coding for a protein showing a similarity to the double-stranded RNA-binding proteins of the DSRM AtDRB-like superfamily. Additionally, based on phylogenetic analysis, we suggest the involvement of replicative beny-like helicases in the evolution of the BMB and TCMB movement genetic modules.

Poorly Conserved P15 Proteins of Cileviruses Retain Elements of Common Ancestry and Putative Functionality: A Theoretical Assessment on the Evolution of Cilevirus Genomes

The genus Cilevirus groups enveloped single-stranded (+) RNA virus members of the family Kitaviridae, order Martellivirales. Proteins P15, scarcely conserved polypeptides encoded by cileviruses, have no apparent homologs in public databases. Accordingly, the open reading frames (ORFs) p15, located at the 5'-end of the viral RNA2 molecules, are considered orphan genes (ORFans). In this study, we have delved into ORFs p15 and the relatively poorly understood biochemical properties of the proteins P15 to posit their importance for viruses across the genus and theorize on their origin. We detected that the ORFs p15 are under purifying selection and that, in some viral strains, the use of synonymous codons is biased, which might be a sign of adaptation to their plant hosts. Despite the high amino acid sequence divergence, proteins P15 show the conserved motif [FY]-L-x(3)-[FL]-H-x-x-[LIV]-S-C-x-C-x(2)-C-x-G-x-C, which occurs exclusively in members of this protein family. Proteins P15 also show a common predicted 3D structure that resembles the helical scaffold of the protein ORF49 encoded by radinoviruses and the phosphoprotein C-terminal domain of mononegavirids. Based on the 3D structural similarities of P15, we suggest elements of common ancestry, conserved functionality, and relevant amino acid residues. We conclude by postulating a plausible evolutionary trajectory of ORFans p15 and the 5'-end of the RNA2 of cileviruses considering both protein fold superpositions and comparative genomic analyses with the closest kitaviruses, negeviruses, nege/kita-like viruses, and unrelated viruses that share the ecological niches of cileviruses.

A Virus Infecting Hibiscus rosa-sinensis Represents an Evolutionary Link Between Cileviruses and Higreviruses

Hibiscus (Hibiscus spp.) are popular ornamental and landscape plants in Hawaii which are susceptible to foliar diseases caused by viruses belonging to the genera Cilevirus and Higrevirus (family Kitaviridae). In this study, a virus infecting H. rosa-sinensis plants displaying foliar symptoms consistent with infection by a kitavirus, including yellow chlorotic blotches with a green perimeter, was characterized. The genome consisted of two RNAs 8.4 and 4.4 kb in length, and was organized most similarly to cileviruses, but with important distinctions. These included the location of the p29 homolog as the 3′-terminal open reading frame (ORF) of RNA2 instead of its typical locus at the 3′-end of RNA1; the absence of a p15 homolog on RNA2 and the adjacent intergenic region which also harbors small putative ORFs of unknown function; and the presence of an ORF encoding a 10 kDa protein at the 3′-terminal end of RNA1 that was also found to be present in the hibiscus green spot virus 2 genome. Spherical particles approximately 55–65 nm in diameter were observed in infected leaf tissue, and viral RNA was detected by reverse-transcription PCR in individual mites collected from symptomatic plants tentatively identified as Brevipalpus yothersi. Although phylogenetic analyses placed this virus between the higrevirus and cilevirus clades, we propose the tentative taxonomic placement of this virus, designated hibiscus yellow blotch virus (HYBV), within the genus Cilevirus.

Reticulon-like properties of a plant virus-encoded movement protein

Plant viruses encode movement proteins (MPs) that ensure the transport of viral genomes through plasmodesmata (PD) and use cell endomembranes, mostly the endoplasmic reticulum (ER), for delivery of viral genomes to PD and formation of PD-anchored virus replication compartments. Here, we demonstrate that the Hibiscus green spot virus BMB2 MP, an integral ER protein, induces constrictions of ER tubules, decreases the mobility of ER luminal content, and exhibits an affinity to highly curved membranes. These properties are similar to those described for reticulons, cellular proteins that induce membrane curvature to shape the ER tubules. Similar to reticulons, BMB2 adopts a W-like topology within the ER membrane. BMB2 targets PD and increases their size exclusion limit, and these BMB2 activities correlate with the ability to induce constrictions of ER tubules. We propose that the induction of ER constrictions contributes to the BMB2-dependent increase in PD permeability and formation of the PD-associated replication compartments, therefore facilitating the virus intercellular spread. Furthermore, we show that the ER tubule constrictions also occur in cells expressing TGB2, one of the three MPs of Potato virus X (PVX), and in PVX-infected cells, suggesting that reticulon-like MPs are employed by diverse RNA viruses.

Small hydrophobic viral proteins involved in intercellular movement of diverse plant virus genomes

Most plant viruses code for movement proteins (MPs) targeting plasmodesmata to enable cell-to-cell and systemic spread in infected plants. Small membrane-embedded MPs have been first identified in two viral transport gene modules, triple gene block (TGB) coding for an RNA-binding helicase TGB1 and two small hydrophobic proteins TGB2 and TGB3 and double gene block (DGB) encoding two small polypeptides representing an RNA-binding protein and a membrane protein. These findings indicated that movement gene modules composed of two or more cistrons may encode the nucleic acid-binding protein and at least one membrane-bound movement protein. The same rule was revealed for small DNA-containing plant viruses, namely, viruses belonging to genus Mastrevirus (family Geminiviridae) and the family Nanoviridae. In multi-component transport modules the nucleic acid-binding MP can be viral capsid protein(s), as in RNA-containing viruses of the families Closteroviridae and Potyviridae. However, membrane proteins are always found among MPs of these multicomponent viral transport systems. Moreover, it was found that small membrane MPs encoded by many viruses can be involved in coupling viral replication and cell-to-cell movement. Currently, the studies of evolutionary origin and functioning of small membrane MPs is regarded as an important pre-requisite for understanding of the evolution of the existing plant virus transport systems. This paper represents the first comprehensive review which describes the whole diversity of small membrane MPs and presents the current views on their role in plant virus movement.

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.

Passion Fruit Green Spot Virus Genome Harbors a New Orphan ORF and Highlights the Flexibility of the 5'-End of the RNA2 Segment Across Cileviruses

Passion fruit green spot and passion fruit sudden death are two reportedly distinct viral diseases that recurrently affect passion fruit (Passiflora spp.) groves in Brazil. Here we used a systematic approach that interconnects symptoms, transmission electron microscopy, RT-PCR detection assays followed by Sanger sequencing, and high-throughput sequencing of the RNA of affected passion fruit plants to gain insights about these diseases. Our data confirmed not only the involvement of cileviruses in these two pathologies, as previously suggested, but also that these viruses belong to the same tentative species: passion fruit green spot virus (PfGSV). Results revealed that PfGSV has a positive-sense RNA genome split into two molecules of approximately 9 kb (RNA1) and 5 kb (RNA2), which share about 50–70% nucleotide sequence identity with other viruses in the genus Cilevirus. Genome sequences of five PfGSV isolates suggest that they have more conserved RNA1 (<5% of nucleotide sequence variability) compared to RNA2 (up to 7% of variability) molecules. The highest nucleotide sequence divergence among PfGSV isolates and other cileviruses is in the genomic segment covering from the 5′-end of the RNA2 until the 5′-end of the open reading frame (ORF) p61, which includes the ORF p15 and the intergenic region. This genomic stretch also harbors a novel orphan ORF encoding a 13 kDa protein presenting a cysteine-rich domain. High variability of 5′-end of the RNA2 in cileviruses is discussed in an evolutionary context assuming that they share putative common ancestors with unclassified arthropod-infecting single-strand positive RNA viruses, including mosquito-specific viruses of the group Negevirus (clades Nelorpivirus and Sandwavirus), and other viruses in the family Kitaviridae.