Evolution of plant virus movement proteins from the 30K superfamily and of their homologs integrated in plant genomes

Study of Endogenous Viruses in the Strawberry Plants

Endogenous viral elements (EVEs) have been reported to exist widely in the genomes of eukaryotic organisms, and they are closely associated with the growth, development, genetics, adaptation, and evolution of their hosts. In this study, two methods-homologous sequence search and genome alignment-were used to explore the endogenous viral sequences in the genomes of Fragaria species. Results revealed abundant endogenous pararetroviruses (EPRVs) in the genomes of Fragaria species, including 786 sequences belonging to five known taxa such as Caulimovirus and other unclassified taxa. Differences were observed in the detected EPRVs between the two methods, with the homologous sequence search having a greater number of EPRVs. On the contrary, genome alignment identified various types and sources of virus-like sequences. Furthermore, through genome alignment, a 267-bp sequence with 95% similarity to the gene encoding the aphid-transmitted protein of Strawberry vein banding virus (Caulimovirus venafragariae) was discovered in the F. chiloensis genome, which was likely a recent insertion. In addition, the statistical analysis of the genome alignment results indicated a remarkably higher abundance of virus-like sequences in the genomes of polyploid strawberries compared with diploid ones. Moreover, the differences in virus-like sequences were observed between the genomes of Fragaria species and those of their close relatives. This study enriched the diversity of viruses that infect strawberries, and laid a theoretical foundation for further research on the origin of endogenous viruses in the strawberry genome, host-virus interactions, adaptation, evolution, and their functions.

Emergence of two distinct spatial folds in a pair of plant virus proteins encoded by nested genes

Virus genomes may encode overlapping or nested open reading frames that increase their coding capacity. It is not known whether the constraints on spatial structures of the two encoded proteins limit the evolvability of nested genes. We examine the evolution of a pair of proteins, p22 and p19, encoded by nested genes in plant viruses from the genus Tombusvirus. The known structure of p19, a suppressor of RNA silencing, belongs to the RAGNYA fold from the alpha+beta class. The structure of p22, the cell-to-cell movement protein from the 30K family widespread in plant viruses, is predicted with the AlphaFold approach, suggesting a single jelly-roll fold core from the all-beta class, structurally similar to capsid proteins from plant and animal viruses. The nucleotide and codon preferences impose modest constraints on the types of secondary structures encoded in the alternative reading frames, nonetheless allowing for compact, well-ordered folds from different structural classes in two similarly-sized nested proteins. Tombusvirus p22 emerged through radiation of the widespread 30K family, which evolved by duplication of a virus capsid protein early in the evolution of plant viruses, whereas lineage-specific p19 may have emerged by a stepwise increase in the length of the overprinted gene and incremental acquisition of functionally active secondary structure elements by the protein product. This evolution of p19 toward the RAGNYA fold represents one of the first documented examples of protein structure convergence in naturally occurring proteins.

Deciphering the virome of Chunkung (Cnidium officinale) showing dwarfism-like symptoms via a high-throughput sequencing analysis

BACKGROUND

Viruses have notable effects on agroecosystems, wherein they can adversely affect plant health and cause problems (e.g., increased biosecurity risks and economic losses). However, our knowledge of their diversity and interactions with specific host plants in ecosystems remains limited. To enhance our understanding of the roles that viruses play in agroecosystems, comprehensive analyses of the viromes of a wide range of plants are essential. High-throughput sequencing (HTS) techniques are useful for conducting impartial and unbiased investigations of plant viromes, ultimately forming a basis for generating further biological and ecological insights. This study was conducted to thoroughly characterize the viral community dynamics in individual plants.

RESULTS

An HTS-based virome analysis in conjunction with proximity sampling and a tripartite network analysis were performed to investigate the viral diversity in chunkung (Cnidium officinale) plants. We identified 61 distinct chunkung plant-associated viruses (27 DNA and 34 RNA viruses) from 21 known genera and 6 unclassified genera in 14 known viral families. Notably, 12 persistent viruses (7 DNA and 5 RNA viruses) were exclusive to dwarfed chunkung plants. The detection of viruses from the families Partitiviridae, Picobirnaviridae, and Spinareoviridae only in the dwarfed plants suggested that they may contribute to the observed dwarfism. The co-infection of chunkung by multiple viruses is indicative of a dynamic and interactive viral ecosystem with significant sequence variability and evidence of recombination.

CONCLUSIONS

We revealed the viral community involved in chunkung. Our findings suggest that chunkung serves as a significant reservoir for a variety of plant viruses. Moreover, the co-infection rate of individual plants was unexpectedly high. Future research will need to elucidate the mechanisms enabling several dozen viruses to co-exist in chunkung. Nevertheless, the important insights into the chunkung virome generated in this study may be relevant to developing effective plant viral disease management and control strategies.

Identification of a negative-strand RNA virus with natural plant and fungal hosts

The presence of viruses that spread to both plant and fungal populations in nature has posed intriguingly scientific question. We found a negative-strand RNA virus related to members of the family Phenuiviridae, named Valsa mali negative-strand RNA virus 1 (VmNSRV1), which induced strong hypovirulence and was prevalent in a population of the phytopathogenic fungus of apple Valsa canker (Valsa mali) infecting apple orchards in the Shaanxi Province of China. Intriguingly, VmNSRV1 encodes a protein with a viral cell-to-cell movement function in plant tissue. Mechanical leaf inoculation showed that VmNSRV1 could systemically infect plants. Moreover, VmNSRV1 was detected in 24 out of 139 apple trees tested in orchards in Shaanxi Province. Fungal inoculation experiments showed that VmNSRV1 could be bidirectionally transmitted between apple plants and V. mali, and VmNSRV1 infection in plants reduced the development of fungal lesions on leaves. Additionally, the nucleocapsid protein encoded by VmNSRV1 is associated with and rearranged lipid droplets in both fungal and plant cells. VmNSRV1 represents a virus that has adapted and spread to both plant and fungal hosts and shuttles between these two organisms in nature (phyto-mycovirus) and is potential to be utilized for the biocontrol method against plant fungal diseases. This finding presents further insights into the virus evolution and adaptation encompassing both plant and fungal hosts.

Megataxonomy and global ecology of the virosphere

Nearly all organisms are hosts to multiple viruses that collectively appear to be the most abundant biological entities in the biosphere. With recent advances in metagenomics and metatranscriptomics, the known diversity of viruses substantially expanded. Comparative analysis of these viruses using advanced computational methods culminated in the reconstruction of the evolution of major groups of viruses and enabled the construction of a virus megataxonomy, which has been formally adopted by the International Committee on Taxonomy of Viruses. This comprehensive taxonomy consists of six virus realms, which are aspired to be monophyletic and assembled based on the conservation of hallmark proteins involved in capsid structure formation or genome replication. The viruses in different major taxa substantially differ in host range and accordingly in ecological niches. In this review article, we outline the latest developments in virus megataxonomy and the recent discoveries that will likely lead to reassessment of some major taxa, in particular, split of three of the current six realms into two or more independent realms. We then discuss the correspondence between virus taxonomy and the distribution of viruses among hosts and ecological niches, as well as the abundance of viruses versus cells in different habitats. The distribution of viruses across environments appears to be primarily determined by the host ranges, i.e. the virome is shaped by the composition of the biome in a given habitat, which itself is affected by abiotic factors.

WIV, a protein domain found in a wide number of arthropod viruses, which probably facilitates infection

The most powerful approach to detect distant homologues of a protein is based on structure prediction and comparison. Yet this approach is still inapplicable to many viral proteins. Therefore, we applied a powerful sequence-based procedure to identify distant homologues of viral proteins. It relies on three principles: (1) traces of sequence similarity can persist beyond the significance cutoff of homology detection programmes; (2) candidate homologues can be identified among proteins with weak sequence similarity to the query by using 'contextual' information, e.g. taxonomy or type of host infected; (3) these candidate homologues can be validated using highly sensitive profile-profile comparison. As a test case, this approach was applied to a protein without known homologues, encoded by ORF4 of Lake Sinai viruses (which infect bees). We discovered that the ORF4 protein contains a domain that has homologues in proteins from >20 taxa of viruses infecting arthropods. We called this domain 'widespread, intriguing, versatile' (WIV), because it is found in proteins with a wide variety of functions and within varied domain contexts. For example, WIV is found in the NSs protein of tospoviruses, a global threat to food security, which infect plants as well as their arthropod vectors; in the RNA2 ORF1-encoded protein of chronic bee paralysis virus, a widespread virus of bees; and in various proteins of cypoviruses, which infect the silkworm Bombyx mori. Structural modelling with AlphaFold indicated that the WIV domain has a previously unknown fold, and bibliographical evidence suggests that it facilitates infection of arthropods.

Endogenous Caulimovirids: Fossils, Zombies, and Living in Plant Genomes

The Caulimoviridae is a family of double-stranded DNA viruses that infect plants. The genomes of most vascular plants contain endogenous caulimovirids (ECVs), a class of repetitive DNA elements that is abundant in some plant genomes, resulting from the integration of viral DNA in the chromosomes of germline cells during episodes of infection that have sometimes occurred millions of years ago. In this review, we reflect on 25 years of research on ECVs that has shown that members of the Caulimoviridae have occupied an unprecedented range of ecological niches over time and shed light on their diversity and macroevolution. We highlight gaps in knowledge and prospects of future research fueled by increased access to plant genome sequence data and new tools for genome annotation for addressing the extent, impact, and role of ECVs on plant biology and the origin and evolutionary trajectories of the Caulimoviridae.

Plant virus movement proteins originated from jelly-roll capsid proteins

Numerous, diverse plant viruses encode movement proteins (MPs) that aid the virus movement through plasmodesmata, the plant intercellular channels. MPs are essential for virus spread and propagation in distal tissues, and several unrelated MPs have been identified. The 30K superfamily of MPs (named after the molecular mass of tobacco mosaic virus (TMV) MP, the classical model of plant virology) is the largest and most diverse MP variety, represented in 16 virus families, but its evolutionary origin remained obscure. Here, we show that the core structural domain of the 30K MPs is homologous to the jelly-roll domain of the capsid proteins (CPs) of small RNA and DNA viruses, in particular, those infecting plants. The closest similarity was observed between the 30K MPs and the CPs of the viruses in the families Bromoviridae and Geminiviridae. We hypothesize that the MPs evolved via duplication or horizontal acquisition of the CP gene in a virus that infected an ancestor of vascular plants, followed by neofunctionalization of one of the paralogous CPs, potentially through the acquisition of unique N- and C-terminal regions. During the subsequent coevolution of viruses with diversifying vascular plants, the 30K MP genes underwent explosive horizontal spread among emergent RNA and DNA viruses, likely permitting viruses of insects and fungi that coinfected plants to expand their host ranges, molding the contemporary plant virome.

Conserved Structure Associated with Different 3′CITEs Is Important for Translation of Umbraviruses

The cap-independent translation of plus-strand RNA plant viruses frequently depends on 3′ structures to attract translation initiation factors that bind ribosomal subunits or bind directly to ribosomes. Umbraviruses are excellent models for studying 3′ cap-independent translation enhancers (3′CITEs), as umbraviruses can have different 3′CITEs in the central region of their lengthy 3′UTRs, and most also have a particular 3′CITE (the T-shaped structure or 3′TSS) near their 3′ ends. We discovered a novel hairpin just upstream of the centrally located (known or putative) 3′CITEs in all 14 umbraviruses. These CITE-associated structures (CASs) have conserved sequences in their apical loops and at the stem base and adjacent positions. In 11 umbraviruses, CASs are preceded by two small hairpins joined by a putative kissing loop interaction (KL). Converting the conserved 6-nt apical loop to a GNRA tetraloop in opium poppy mosaic virus (OPMV) and pea enation mosaic virus 2 (PEMV2) enhanced translation of genomic (g)RNA, but not subgenomic (sg)RNA reporter constructs, and significantly repressed virus accumulation in Nicotiana benthamiana. Other alterations throughout OPMV CAS also repressed virus accumulation and only enhanced sgRNA reporter translation, while mutations in the lower stem repressed gRNA reporter translation. Similar mutations in the PEMV2 CAS also repressed accumulation but did not significantly affect gRNA or sgRNA reporter translation, with the exception of deletion of the entire hairpin, which only reduced translation of the gRNA reporter. OPMV CAS mutations had little effect on the downstream BTE 3′CITE or upstream KL element, while PEMV2 CAS mutations significantly altered KL structures. These results introduce an additional element associated with different 3′CITEs that differentially affect the structure and translation of different umbraviruses.

Transcriptome Mining Expands Knowledge of RNA Viruses across the Plant Kingdom

Our current understanding of plant viruses stems largely from those affecting economically important plants. Yet plant species in cultivation represent a small and biased subset of the plant kingdom. Here, we describe virus diversity and abundance in 1,079 transcriptomes from species across the breadth of the plant kingdom (Archaeplastida) by analyzing open-source data from the 1000 Plant Transcriptomes Initiative (1KP). We identified 104 potentially novel viruses, of which 40% were single-stranded positive-sense RNA viruses across eight orders, including members of the Hepelivirales, Tymovirales, Cryppavirales, Martellivirales, and Picornavirales. One-third of the newly described viruses were double-stranded RNA viruses from the orders Durnavirales and Ghabrivirales. The remaining were negative-sense RNA viruses from the Rhabdoviridae, Aspiviridae, Yueviridae, and Phenuiviridae and the newly proposed Viridisbunyaviridae. Our analysis considerably expands the known host range of 13 virus families to include lower plants (e.g., Benyviridae and Secoviridae) and 4 virus families to include alga hosts (e.g., Tymoviridae and Chrysoviridae). More broadly, however, a cophylogeny analysis revealed that the evolutionary history of these families is largely driven by cross-species transmission events. The discovery of the first 30-kDa movement protein in a nonvascular plant suggests that the acquisition of plant virus movement proteins occurred prior to the emergence of the plant vascular system. Together, these data highlight that numerous RNA virus families are associated with older evolutionary plant lineages than previously thought and that the apparent scarcity of RNA viruses found in lower plants likely reflects a lack of investigation rather than their absence. IMPORTANCE Our knowledge of plant viruses is mainly limited to those infecting economically important host species. In particular, we know little about those viruses infecting basal plant lineages such as the ferns, lycophytes, bryophytes, and charophytes. To expand this understanding, we conducted a broad-scale viral survey of species across the breadth of the plant kingdom. We found that basal plants harbor a wide diversity of RNA viruses, including some that are sufficiently divergent to likely compose a new virus family. The basal plant virome revealed offers key insights into the evolutionary history of core plant virus gene modules and genome segments. More broadly, this work emphasizes that the scarcity of viruses found in these species to date most likely reflects the limited research in this area.

An Evolved 5' Untranslated Region of Alfalfa Mosaic Virus Allows the RNA Transport of Movement-Defective Variants

Although the coat protein (CP) has a relevant role in the long-distance movement of alfalfa mosaic virus (AMV) and brome mosaic virus (BMV), its precise function is not fully understood. Previous results showed that a specific interaction between the C termini of the movement protein (MP) and the cognate CP is required for systemic transport. Thus, we have performed a compensatory evolution experiment using an AMV RNA3 derivative defective in long-distance transport that carries a BMV MP lacking the C-terminal 48 residues and unable to interact with the AMV CP. After several passages, five independent evolution lineages were able to move long distance. The analysis of the viral RNA of these lineages showed the presence of three different modifications located exclusively at the 5' untranslated region (5' UTR). The three evolved 5' UTR variants accumulated comparable levels of viral RNA and CP but reduced the accumulation of virus particles and the affinity between the 5' UTR and the AMV CP. In addition, the evolved 5' UTR increased cell-to-cell transport for both the AMV RNA3 carrying the BMV MP and that carrying the AMV MP. Finally, the evolved 5' UTRs allowed the systemic transport of an AMV RNA3 carrying a CP mutant defective in virus particles and increased the systemic transport of several AMV RNA3 derivatives carrying different viral MPs associated with the 30K superfamily. Altogether, our findings indicate that virus particles are not required for the systemic transport of AMV but also that BMV MP is competent for the short- and long-distance transport without the interaction with the CP. IMPORTANCE The results obtained in the present work could challenge the view of the role of the virus particle in the systemic transport of plant viruses. In this sense, we show that two different MPs are competent to systemically transport the AMV genome without the requirement of the virus particles, as reported for viruses lacking a CP (e.g., Umbravirus). The incapability of the viral MP to interact with the CP triggered virus variants that evolved to reduce the formation of virus particles, probably to increase the accessibility of the MP to the viral progeny. Our results point to the idea that virus particles would not be necessary for the viral systemic transport but would be necessary for vector virus transmission. This idea is reinforced by the observation that heterologous MPs also increased the systemic transport of the AMV constructs that have reduced encapsidation capabilities.

Viruses Infecting Trees and Herbs That Produce Edible Fleshy Fruits with a Prominent Value in the Global Market: An Evolutionary Perspective

Trees and herbs that produce fruits represent the most valuable agricultural food commodities in the world. However, the yield of these crops is not fully achieved due to biotic factors such as bacteria, fungi, and viruses. Viruses are capable of causing alterations in plant growth and development, thereby impacting the yield of their hosts significantly. In this work, we first compiled the world's most comprehensive list of known edible fruits that fits our definition. Then, plant viruses infecting those trees and herbs that produce fruits with commercial importance in the global market were identified. The identified plant viruses belong to 30 families, most of them containing single-stranded RNA genomes. Importantly, we show the overall picture of the host range for some virus families following an evolutionary approach. Further, the current knowledge about plant-virus interactions, focusing on the main disorders they cause, as well as yield losses, is summarized. Additionally, since accurate diagnosis methods are of pivotal importance for viral diseases control, the current and emerging technologies for the detection of these plant pathogens are described. Finally, the most promising strategies employed to control viral diseases in the field are presented, focusing on solutions that are long-lasting.

Spontaneous Mutation in the Movement Protein of Citrus Leprosis Virus C2, in a Heterologous Virus Infection Context, Increases Cell-to-Cell Transport and Generates Fitness Advantage

Previous results using a movement defective alfalfa mosaic virus (AMV) vector revealed that citrus leprosis virus C (CiLV-C) movement protein (MP) generates a more efficient local movement, but not more systemic transport, than citrus leprosis virus C2 (CiLV-C2) MP, MPs belonging to two important viruses for the citrus industry. Here, competition experiment assays in transgenic tobacco plants (P12) between transcripts of AMV constructs expressing the cilevirus MPs, followed by several biological passages, showed the prevalence of the AMV construct carrying the CiLV-C2 MP. The analysis of AMV RNA 3 progeny recovered from P12 plant at the second viral passage revealed the presence of a mix of progeny encompassing the CiLV-C2 MP wild type (MP^WT) and two variants carrying serines instead phenylalanines at positions 72 (MP^S72F) or 259 (MP^S259F), respectively. We evaluated the effects of each modified residue in virus replication, and cell-to-cell and long-distance movements. Results indicated that phenylalanine at position 259 favors viral cell-to-cell transport with an improvement in viral fitness, but has no effect on viral replication, whereas mutation at position 72 (MP^S72F) has a penalty in the viral fitness. Our findings indicate that the prevalence of a viral population may be correlated with its greater efficiency in cell-to-cell and systemic movements.

Frontiers in the Standardization of the Plant Platform for High Scale Production of Vaccines

The recent COVID-19 pandemic has highlighted the value of technologies that allow a fast setup and production of biopharmaceuticals in emergency situations. The plant factory system can provide a fast response to epidemics/pandemics. Thanks to their scalability and genome plasticity, plants represent advantageous platforms to produce vaccines. Plant systems imply less complicated production processes and quality controls with respect to mammalian and bacterial cells. The expression of vaccines in plants is based on transient or stable transformation systems and the recent progresses in genome editing techniques, based on the CRISPR/Cas method, allow the manipulation of DNA in an efficient, fast, and easy way by introducing specific modifications in specific sites of a genome. Nonetheless, CRISPR/Cas is far away from being fully exploited for vaccine expression in plants. In this review, an overview of the potential conjugation of the renewed vaccine technologies (i.e., virus-like particles-VLPs, and industrialization of the production process) with genome editing to produce vaccines in plants is reported, illustrating the potential advantages in the standardization of the plant platforms, with the overtaking of constancy of large-scale production challenges, facilitating regulatory requirements and expediting the release and commercialization of the vaccine products of genome edited plants.

Broken, silent, and in hiding: tamed endogenous pararetroviruses escape elimination from the genome of sugar beet (Beta vulgaris)

BACKGROUND AND AIMS

Endogenous pararetroviruses (EPRVs) are widespread components of plant genomes that originated from episomal DNA viruses of the Caulimoviridae family. Due to fragmentation and rearrangements, most EPRVs have lost their ability to replicate through reverse transcription and to initiate viral infection. Similar to the closely related retrotransposons, extant EPRVs were retained and often amplified in plant genomes for several million years. Here, we characterize the complete genomic EPRV fraction of the crop sugar beet (Beta vulgaris, Amaranthaceae) to understand how they shaped the beet genome and to suggest explanations for their absent virulence.

METHODS

Using next- and third-generation sequencing data and genome assembly, we reconstructed full-length in silico representatives for the three host-specific EPRVs (beetEPRVs) in the B. vulgaris genome. Focusing on the endogenous caulimovirid beetEPRV3, we investigated its chromosomal localization, abundance and distribution by fluorescent in situ and Southern hybridization.

KEY RESULTS

Full-length beetEPRVs range between 7.5 and 10.7 kb in size, are heterogeneous in structure and sequence, and occupy about 0.3 % of the beet genome. Although all three beetEPRVs were assigned to the florendoviruses, they showed variably arranged protein-coding domains, different fragmentation, and preferences for diverse sequence contexts. We observed small RNAs that specifically target the individual beetEPRVs, indicating stringent epigenetic suppression. BeetEPRV3 sequences occur along all sugar beet chromosomes, preferentially in the vicinity of each other and are associated with heterochromatic, centromeric and intercalary satellite DNAs. BeetEPRV3 members also exist in genomes of related wild species, indicating an initial beetEPRV3 integration 13.4-7.2 million years ago.

CONCLUSIONS

Our study in beet illustrates the variability of EPRV structure and sequence in a single host genome. Evidence of sequence fragmentation and epigenetic silencing implies possible plant strategies to cope with long-term persistence of EPRVs, including amplification, fixation in the heterochromatin, and containment of EPRV virulence.

The Bunyavirales: The Plant-Infecting Counterparts

Negative-strand (-) RNA viruses (NSVs) comprise a large and diverse group of viruses that are generally divided in those with non-segmented and those with segmented genomes. Whereas most NSVs infect animals and humans, the smaller group of the plant-infecting counterparts is expanding, with many causing devastating diseases worldwide, affecting a large number of major bulk and high-value food crops. In 2018, the taxonomy of segmented NSVs faced a major reorganization with the establishment of the order Bunyavirales. This article overviews the major plant viruses that are part of the order, i.e., orthospoviruses (Tospoviridae), tenuiviruses (Phenuiviridae), and emaraviruses (Fimoviridae), and provides updates on the more recent ongoing research. Features shared with the animal-infecting counterparts are mentioned, however, special attention is given to their adaptation to plant hosts and vector transmission, including intra/intercellular trafficking and viral counter defense to antiviral RNAi.

Molecular Epidemiology of Citrus Leprosis Virus C: A New Viral Lineage and Phylodynamic of the Main Viral Subpopulations in the Americas

Despite the importance of viral strains/variants as agents of emerging diseases, genetic and evolutionary processes affecting their ecology are not fully understood. To get insight into this topic, we assessed the population and spatial dynamic parameters of citrus leprosis virus C (CiLV-C, genus Cilevirus, family Kitaviridae). CiLV-C is the etiological agent of citrus leprosis disease, a non-systemic infection considered the main viral disorder affecting citrus orchards in Brazil. Overall, we obtained 18 complete or near-complete viral genomes, 123 complete nucleotide sequences of the open reading frame (ORF) encoding the putative coat protein, and 204 partial nucleotide sequences of the ORF encoding the movement protein, from 430 infected Citrus spp. samples collected between 1932 and 2020. A thorough examination of the collected dataset suggested that the CiLV-C population consists of the major lineages CRD and SJP, unevenly distributed, plus a third one called ASU identified in this work, which is represented by a single isolate found in an herbarium sample collected in Asuncion, Paraguay, in 1937. Viruses from the three lineages share about 85% nucleotide sequence identity and show signs of inter-clade recombination events. Members of the lineage CRD were identified both in commercial and non-commercial citrus orchards. However, those of the lineages SJP were exclusively detected in samples collected in the citrus belt of São Paulo and Minas Gerais, the leading Brazilian citrus production region, after 2015. The most recent common ancestor of viruses of the three lineages dates back to, at least, ∼1500 years ago. Since citrus plants were introduced in the Americas by the Portuguese around the 1520s, the Bayesian phylodynamic analysis suggested that the ancestors of the main CiLV-C lineages likely originated in contact with native vegetation of South America. The intensive expansion of CRD and SJP lineages in Brazil started probably linked to the beginning of the local citrus industry. The high prevalence of CiLV-C in the citrus belt of Brazil likely ensues from the intensive connectivity between orchards, which represents a potential risk toward pathogen saturation across the region.

Membrane Association and Topology of Citrus Leprosis Virus C2 Movement and Capsid Proteins

Although citrus leprosis disease has been known for more than a hundred years, one of its causal agents, citrus leprosis virus C2 (CiLV-C2), is poorly characterized. This study described the association of CiLV-C2 movement protein (MP) and capsid protein (p29) with biological membranes. Our findings obtained by computer predictions, chemical treatments after membrane fractionation, and biomolecular fluorescence complementation assays revealed that p29 is peripherally associated, while the MP is integrally bound to the cell membranes. Topological analyses revealed that both the p29 and MP expose their N- and C-termini to the cell cytoplasmic compartment. The implications of these results in the intracellular movement of the virus were discussed.

Unravelling the involvement of cilevirus p32 protein in the viral transport

Citrus leprosis (CL) is a severe disease that affects citrus orchards mainly in Latin America. It is caused by Brevipalpus-transmitted viruses from genera Cilevirus and Dichorhavirus. Currently, no reports have explored the movement machinery for the cilevirus. Here, we have performed a detailed functional study of the p32 movement protein (MP) of two cileviruses. Citrus leprosis-associated viruses are not able to move systemically in neither their natural nor experimental host plants. However, here we show that cilevirus MPs are able to allow the cell-to-cell and long-distance transport of movement-defective alfalfa mosaic virus (AMV). Several features related with the viral transport were explored, including: (i) the ability of cilevirus MPs to facilitate virus movement on a nucleocapsid assembly independent-manner; (ii) the generation of tubular structures from transient expression in protoplast; (iii) the capability of the N- and C- terminus of MP to interact with the cognate capsid protein (p29) and; (iv) the role of the C-terminus of p32 in the cell-to-cell and long-distance transport, tubule formation and the MP-plasmodesmata co-localization. The MP was able to direct the p29 to the plasmodesmata, whereby the C-terminus of MP is independently responsible to recruit the p29 to the cell periphery. Furthermore, we report that MP possess the capacity to enter the nucleolus and to bind to a major nucleolar protein, the fibrillarin. Based on our findings, we provide a model for the role of the p32 in the intra- and intercellular viral spread.

Insertion of Badnaviral DNA in the Late Blight Resistance Gene (R1a) of Brinjal Eggplant (Solanum melongena)

Endogenous viral elements (EVEs) are widespread in plant genomes. They result from the random integration of viral sequences into host plant genomes by horizontal DNA transfer and have the potential to alter host gene expression. We performed a large-scale search for co-transcripts including caulimovirid and plant sequences in 1,678 plant and 230 algal species and characterized 50 co-transcripts in 45 distinct plant species belonging to lycophytes, ferns, gymnosperms and angiosperms. We found that insertion of badnavirus EVEs along with Ty-1 copia mobile elements occurred into a late blight resistance gene (R1) of brinjal eggplant (Solanum melongena) and wild relatives in genus Solanum and disrupted R1 orthologs. EVEs of two previously unreported badnaviruses were identified in the genome of S. melongena, whereas EVEs from an additional novel badnavirus were identified in the genome of S. aethiopicum, the cultivated scarlet eggplant. Insertion of these viruses in the ancestral lineages of the direct wild relatives of the eggplant would have occurred during the last 3 Myr, further supporting the distinctiveness of the group of the eggplant within the giant genus Solanum.

Participation of Multifunctional RNA in Replication, Recombination and Regulation of Endogenous Plant Pararetroviruses (EPRVs)

Pararetroviruses, taxon Caulimoviridae, are typical of retroelements with reverse transcriptase and share a common origin with retroviruses and LTR retrotransposons, presumably dating back 1.6 billion years and illustrating the transition from an RNA to a DNA world. After transcription of the viral genome in the host nucleus, viral DNA synthesis occurs in the cytoplasm on the generated terminally redundant RNA including inter- and intra-molecule recombination steps rather than relying on nuclear DNA replication. RNA recombination events between an ancestral genomic retroelement with exogenous RNA viruses were seminal in pararetrovirus evolution resulting in horizontal transmission and episomal replication. Instead of active integration, pararetroviruses use the host DNA repair machinery to prevail in genomes of angiosperms, gymnosperms and ferns. Pararetrovirus integration - leading to Endogenous ParaRetroViruses, EPRVs - by illegitimate recombination can happen if their sequences instead of homologous host genomic sequences on the sister chromatid (during mitosis) or homologous chromosome (during meiosis) are used as template. Multiple layers of RNA interference exist regulating episomal and chromosomal forms of the pararetrovirus. Pararetroviruses have evolved suppressors against this plant defense in the arms race during co-evolution which can result in deregulation of plant genes. Small RNAs serve as signaling molecules for Transcriptional and Post-Transcriptional Gene Silencing (TGS, PTGS) pathways. Different populations of small RNAs comprising 21-24 nt and 18-30 nt in length have been reported for Citrus, Fritillaria, Musa, Petunia, Solanum and Beta. Recombination and RNA interference are driving forces for evolution and regulation of EPRVs.

Natural Defect of a Plant Rhabdovirus Glycoprotein Gene: A Case Study of Virus-Plant Coevolution

Seven isolates of a putative cytorhabdovirus (family Rhabdoviridae, order Mononegavirales) designated as citrus-associated rhabdovirus (CiaRV) were identified in citrus, passion fruit, and paper bush from the same geographical area in China. CiaRV, bean-associated cytorhabdovirus (Brazil), and papaya virus E (Ecuador) should be taxonomically classified in the species Papaya cytorhabdovirus. Due to natural mutations, the glycoprotein (G) and P4 genes were impaired in citrus-infecting isolates of CiaRV, resulting in an atypical rhabdovirus genome organization of 3' leader-N-P-P3-M-L-5' trailer. The P3 protein of CiaRV shared a common origin with begomoviral movement proteins (family Geminiviridae). Secondary structure analysis and trans-complementation of movement-deficient tomato mosaic virus and potato virus X mutants by CiaRV P3 supported its function in viral cell-to-cell trafficking. The wide geographical dispersal of CiaRV and related viruses suggests an efficient transmission mechanism, as well as an underlying risk to global agriculture. Both the natural phenomenon and experimental analyses demonstrated presence of the "degraded" type of CiaRV in citrus, in parallel to "undegraded" types in other host plant species. This case study shows a plant virus losing the function of an important but nonessential gene, likely due to host shift and adaption, which deepened our understanding of course of natural viral diversification.

Intercellular movement of plant RNA viruses: Targeting replication complexes to the plasmodesma for both accuracy and efficiency

Replication and movement are two critical steps in plant virus infection. Recent advances in the understanding of the architecture and subcellular localization of virus-induced inclusions and the interactions between viral replication complex (VRC) and movement proteins (MPs) allow for the dissection of the intrinsic relationship between replication and movement, which has revealed that recruitment of VRCs to the plasmodesma (PD) via direct or indirect MP-VRC interactions is a common strategy used for cell-to-cell movement by most plant RNA viruses. In this review, we summarize the recent advances in the understanding of virus-induced inclusions and their roles in virus replication and cell-to-cell movement, analyze the advantages of such coreplicational movement from a viral point of view and discuss the possible mechanical force by which MPs drive the movement of virions or viral RNAs through the PD. Finally, we highlight the missing pieces of the puzzle of viral movement that are especially worth investigating in the near future.

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.

Family Level Phylogenies Reveal Relationships of Plant Viruses within the Order Bunyavirales

Bunyavirales are negative-sense segmented RNA viruses infecting arthropods, protozoans, plants, and animals. This study examines the phylogenetic relationships of plant viruses within this order, many of which are recently classified species. Comprehensive phylogenetic analyses of the viral RNA dependent RNA polymerase (RdRp), precursor glycoprotein (preGP), the nucleocapsid (N) proteins point toward common progenitor viruses. The RdRp of Fimoviridae and Tospoviridae show a close evolutional relationship while the preGP of Fimoviridae and Phenuiviridae show a closed relationship. The N proteins of Fimoviridae were closer to the Phasmaviridae, the Tospoviridae were close to some Phenuiviridae members and the Peribunyaviridae. The plant viral movement proteins of species within the Tospoviridae and Phenuiviridae were more closely related to each other than to members of the Fimoviridae. Interestingly, distal ends of 3′ and 5′ untranslated regions of species within the Fimoviridae shared similarity to arthropod and vertebrate infecting members of the Cruliviridae and Peribunyaviridae compared to other plant virus families. Co-phylogeny analysis of the plant infecting viruses indicates that duplication and host switching were more common than co-divergence with a host species.

Deep Roots and Splendid Boughs of the Global Plant Virome

Land plants host a vast and diverse virome that is dominated by RNA viruses, with major additional contributions from reverse-transcribing and single-stranded (ss) DNA viruses. Here, we introduce the recently adopted comprehensive taxonomy of viruses based on phylogenomic analyses, as applied to the plant virome. We further trace the evolutionary ancestry of distinct plant virus lineages to primordial genetic mobile elements. We discuss the growing evidence of the pivotal role of horizontal virus transfer from invertebrates to plants during the terrestrialization of these organisms, which was enabled by the evolution of close ecological associations between these diverse organisms. It is our hope that the emerging big picture of the formation and global architecture of the plant virome will be of broad interest to plant biologists and virologists alike and will stimulate ever deeper inquiry into the fascinating field of virus-plant coevolution.

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.

Complete nucleotide sequence of loquat virus A, a member of the family Betaflexiviridae with a novel genome organization

Analysis of a loquat tree with leaf curl symptoms by deep sequencing revealed a novel virus with a single-stranded RNA genome, for which the name "loquat virus A" (LoVA) was proposed. The complete genome sequence comprised 7553 nucleotides (nt) and an additional poly(A) tail at the 3' terminus. Sequence comparisons of LoVA showed moderate similarity to cherry virus A (CVA), currant virus A (CuVA), and mume virus A (MuVA), which are members of the genus Capillovirus in the family Betaflexiviridae. Phylogenetic analysis of the full-genome nt sequence and replicase-like protein supported the placement of LoVA within the genus Capillovirus. However, it has a distinct genome organization, differing from recognized capilloviruses, as it contains three open reading frames (ORFs), with the coat protein (CP) expressed separately from the replication-associated protein (RP) rather than being encoded in the same ORF. This indicates that LoVA is a novel member of the genus Capillovirus in the family Betaflexiviridae with a distinct genomic organization.

Host-specific loss of sequences of an alfalfa mosaic virus isolate during systemic infection

Infectivity of an alfalfa mosaic virus (AMV) isolate from Leonotis nepetaefolia in different tomato cultivars was analyzed. Symptoms typical of AMV infection were observed in indicator plants, but not in Flora Dade and Rio Grande tomato cultivars; however, mild symptoms were observed in cv. Rutgers. Furthermore, at least 1 kb of the 3´ segment of RNA 2 and the coat protein gene were missing in systemic leaves of inoculated Rio Grande and Flora Dade plants, while in cv. Rutgers infected with this AMV strain all genomic components were detected. Northern blot analysis of plants infected with the aforementioned AMV isolate confirmed the absence of the CP gene, but suggested rearrangements in both RNA 2 and 3. Factors that may affect differential movement or systemic accumulation of genomic components in multipartite viruses in plants are discussed.

Specificity of Plant Rhabdovirus Cell-to-Cell Movement

Positive-stranded RNA virus movement proteins (MPs) generally lack sequence-specific nucleic acid-binding activities and display cross-family movement complementarity with related and unrelated viruses. Negative-stranded RNA plant rhabdoviruses encode MPs with limited structural and functional relatedness with other plant virus counterparts, but the precise mechanisms of intercellular transport are obscure. In this study, we first analyzed the abilities of MPs encoded by five distinct rhabdoviruses to support cell-to-cell movement of two positive-stranded RNA viruses by using trans-complementation assays. Each of the five rhabdovirus MPs complemented the movement of MP-defective mutants of tomato mosaic virus and potato X virus. In contrast, movement of recombinant MP deletion mutants of sonchus yellow net nucleorhabdovirus (SYNV) and tomato yellow mottle-associated cytorhabdovirus (TYMaV) was rescued only by their corresponding MPs, i.e., SYNV sc4 and TYMaV P3. Subcellular fractionation analyses revealed that SYNV sc4 and TYMaV P3 were peripherally associated with cell membranes. A split-ubiquitin membrane yeast two-hybrid assay demonstrated specific interactions of the membrane-associated rhabdovirus MPs only with their cognate nucleoproteins (N) and phosphoproteins (P). More importantly, SYNV sc4-N and sc4-P interactions directed a proportion of the N-P complexes from nuclear sites of replication to punctate loci at the cell periphery that partially colocalized with the plasmodesmata. Our data show that cell-to-cell movement of plant rhabdoviruses is highly specific and suggest that cognate MP-nucleocapsid core protein interactions are required for intra- and intercellular trafficking.IMPORTANCE Local transport of plant rhabdoviruses likely involves the passage of viral nucleocapsids through MP-gated plasmodesmata, but the molecular mechanisms are not fully understood. We have conducted complementation assays with MPs encoded by five distinct rhabdoviruses to assess their movement specificity. Each of the rhabdovirus MPs complemented the movement of MP-defective mutants of two positive-stranded RNA viruses that have different movement strategies. In marked contrast, cell-to-cell movement of two recombinant plant rhabdoviruses was highly specific in requiring their cognate MPs. We have shown that these rhabdovirus MPs are localized to the cell periphery and associate with cellular membranes, and that they interact only with their cognate nucleocapsid core proteins. These interactions are able to redirect viral nucleocapsid core proteins from their sites of replication to the cell periphery. Our study provides a model for the specific inter- and intracellular trafficking of plant rhabdoviruses that may be applicable to other negative-stranded RNA viruses.

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.

Discovery of highly divergent lineages of plant-associated astro-like viruses sheds light on the emergence of potyviruses

RNA viruses are believed to have originated from a common ancestor, but how this ancestral genome evolved into the large variety of genomic architectures and viral proteomes we see today remains largely unknown. Tackling this question is hindered by the lack of universally conserved proteins other than the RNA-dependent RNA polymerase (RdRp) as well as a limited RNA virus sampling. The latter is still heavily biased towards relatively few viral lineages from a non-representative collection of hosts, which complicates studies aiming to reveal possible trajectories during the evolution of RNA virus genomes that are favored over others. We report the discovery of 11 highly divergent lineages of viruses with genomic architectures that resemble those of the astroviruses. These genomes were initially identified through a sequence homology search in more than 6600 plant transcriptome projects from the Sequence Read Archive (SRA) using astrovirus representatives as query. Seed-based viral genome assembly of unprocessed SRA data for several dozens of the most promising hits resulted in two viral genome sequences with full-length coding regions, nine partial genomes and a much larger number of short sequence fragments. Genomic and phylogenetic characterization of the 11 discovered viruses, which we coined plastroviruses (plant-associated astro-like viruses), showed that they are related to both astro- and potyviruses and allowed us to identify divergent Serine protease, RdRp and viral capsid domains encoded in the plastrovirus genome. Interestingly, some of the plastroviruses shared different features with potyviruses including the replacement of the catalytic Ser by a Cys residue in the protease active site. These results suggest that plastroviruses may have reached different points on an evolutionary trajectory from astro-like to poty-like genomes. A model how potyviruses might have emerged from (pl)astro-like ancestors in a multi-step process is discussed.

A Negative-Stranded RNA Virus Infecting Citrus Trees: The Second Member of a New Genus Within the Order Bunyavirales

A new RNA virus has been identified from a sweet orange tree in southern Italy. This virus, tentatively named citrus virus A (CiVA), has a bipartite genome composed of (i) a negative-stranded (ns) RNA1, encoding the viral RNA-dependent RNA polymerase (RdRp), and (ii) an ambisense RNA2, coding for the putative movement protein (MP) and nucleocapsid protein (NP), with the two open reading frames separated by a long AU-rich intergenic region (IR) adopting a hairpin conformation. CiVA genomic RNAs and the encoded proteins resemble those of the recently discovered citrus concave gum-associated virus (CCGaV). This CCGaV, a nsRNA virus associated with the ancient citrus concave gum disease, has been proposed as the representative member of a new genus tentatively named Coguvirus. Molecular and phylogenetic analyses presented here support the classification of CiVA, and likely of other two recently described nsRNA viruses infecting plants, in this new genus. By showing that the evolutionary origin of the MP of all the putative coguviruses likely differs from that of their respective RdRp and NP, this study also provides evidence of a likely modular genome evolution for these viruses. Moreover, phylogenetic data support the proposal that, during the evolutionary history of nsRNA viruses, the plant-infecting viruses most likely emerged from an invertebrate-infecting ancestor several times as independent events. CiVA was identified in a field sweet orange tree not showing any obvious symptom and was graft-transmitted to sweet orange, grapefruit, rough lemon and Dweet tangor indicator plants that did not developed symptoms. The capacity of infecting citrus hosts of several species was also confirmed by a preliminary survey that identified orange, mandarin, clementine and lemon trees as natural hosts of CiVA in several fields of southern Italy, again without any obvious association with specific symptoms.

Discovery of Negative-Sense RNA Viruses in Trees Infected with Apple Rubbery Wood Disease by Next-Generation Sequencing

Apple rubbery wood is a disease of apple found around the world, often associated with Apple flat limb disease, and regulated in many countries. Despite its long history in apple cultivation, the disease's causal agent has remained elusive. In this study, next-generation sequencing (NGS) was used to identify and characterize several related novel viral agents from apple rubbery wood-infected plants, which have been named Apple rubbery wood virus (ARWV) 1 and 2. Additional specimens with apple rubbery wood disease tested positive by polymerase chain reaction with primers designed to ARWV 1 and 2 genomic RNA segments. In an NGS-based screening of over 100 Malus and 100 Prunus specimens from a collection of virus-infected trees, only one Malus specimen was found to be infected with ARWV not known to be infected with the disease, which strongly suggests that ARWV is not commonly found in Malus spp. or other fruit trees. The two viruses are most closely related to members of the order Bunyavirales. Three RNA segments (large, medium, and small) were characterized and the viruses likely represent a new genus under the family Phenuiviridae, with a suggested name of Rubodvirus (Rubbery wood virus).

The first phlebo-like virus infecting plants: a case study on the adaptation of negative-stranded RNA viruses to new hosts

A novel negative-stranded (ns) RNA virus associated with a severe citrus disease reported more than 80 years ago has been identified. Transmission electron microscopy showed that this novel virus, tentatively named citrus concave gum-associated virus, is flexuous and non-enveloped. Notwithstanding, its two genomic RNAs share structural features with members of the genus Phlebovirus, which are enveloped arthropod-transmitted viruses infecting mammals, and with a group of still unclassified phlebo-like viruses mainly infecting arthropods. CCGaV genomic RNAs code for an RNA-dependent RNA polymerase, a nucleocapsid protein and a putative movement protein showing structural and phylogenetic relationships with phlebo-like viruses, phleboviruses and the unrelated ophioviruses, respectively, thus providing intriguing evidence of a modular genome evolution. Phylogenetic reconstructions identified an invertebrate-restricted virus as the most likely ancestor of this virus, revealing that its adaptation to plants was independent from and possibly predated that of the other nsRNA plant viruses. These data are consistent with an evolutionary scenario in which trans-kingdom adaptation occurred several times during the history of nsRNA viruses and followed different evolutionary pathways, in which genomic RNA segments were gained or lost. The need to create a new genus for this bipartite nsRNA virus and the impact of the rapid and specific detection methods developed here on citrus sanitation and certification are also discussed.

Euphyllophyte Paleoviruses Illuminate Hidden Diversity and Macroevolutionary Mode of Caulimoviridae

Endogenous viral elements (paleoviruses) provide "molecular fossils" for studying the deep history and macroevolution of viruses. Endogenous plant pararetroviruses (EPRVs) are widespread in angiosperms, but little is known about EPRVs in earlier-branching plants. Here we use a large-scale phylogenomic approach to investigate the diversity and macroevolution of plant pararetroviruses (formally known as Caulimoviridae). We uncover an unprecedented and unappreciated diversity of EPRVs within the genomes of gymnosperms and ferns. The known angiosperm viruses constitute only a minor part of the Caulimoviridae diversity. By characterizing the distribution of EPRVs, we show that no major euphyllophyte lineages escape the activity of Caulimoviridae, raising the possibility that many exogenous Caulimoviridae remain to be discovered in euphyllophytes. We find that the copy numbers of EPRVs are generally high, suggesting that EPRVs might define a unique group of repetitive elements and represent important components of euphyllophyte genomes. Evolutionary analyses suggest an ancient origin of Caulimoviridae and at least three independent origins of Caulimoviridae in angiosperms. Our findings reveal the remarkable diversity of Caulimoviridae and have important implications for understanding the origin and macroevolution of plant pararetroviruses.IMPORTANCE Few viruses have been documented in plants outside angiosperms. Viruses can occasionally integrate into host genomes, forming endogenous viral elements (EVEs). Endogenous plant pararetroviruses (EPRVs) are widespread in angiosperms. In this study, we performed comprehensive comparative and phylogenetic analyses of EPRVs and found that EPRVs are present in the genomes of gymnosperms and ferns. We identified numerous EPRVs in gymnosperm and fern genomes, revealing an unprecedented depth in the diversity of plant pararetroviruses. Plant pararetroviruses mainly underwent cross-species transmission, and angiosperm pararetroviruses arose at least three times. Our study provides novel insights into the diversity and macroevolution of plant pararetroviruses.

Exploring the role of cellular homologous of the 30K-superfamily of plant virus movement proteins

Genes orthologous to the 30K-superfamily of movement proteins (MP) from plant viruses have been recently discovered by bioinformatics analyses as integrated elements in the genome of most vascular plants. However, their functional relevance for plants is still unclear. Here, we undertake some preliminary steps into the functional characterization of one of these putative MP genes found in Arabidopsis thaliana. We found that the AtMP gene is expressed at different stages of the plant development, with accumulation being highest in flowers but lowest in mature siliques. We also found down-regulation of the gene may result in a small delay in plant development and in an exacerbation of the negative effect of salinity in germination efficiency. We have also explored whether changes in expression of the endogenous AtMP have any effect on susceptibility to infection with several viruses, and found that the infectivity of tobacco rattle tobravirus was strongly dependent on the expression of the endogenous AtMP. Finally, we have cloned the endogenous MP from four different plant species into an expression vector that allows for specifically assessing their activity as cell-to-cell movement proteins and have shown that though some may still retain the ancestral activity, they do so in a quite inefficient manner, thus suggesting they have acquired a novel function during adaptation to the host genome.

Tracheophyte genomes keep track of the deep evolution of the Caulimoviridae

Endogenous viral elements (EVEs) are viral sequences that are integrated in the nuclear genomes of their hosts and are signatures of viral infections that may have occurred millions of years ago. The study of EVEs, coined paleovirology, provides important insights into virus evolution. The Caulimoviridae is the most common group of EVEs in plants, although their presence has often been overlooked in plant genome studies. We have refined methods for the identification of caulimovirid EVEs and interrogated the genomes of a broad diversity of plant taxa, from algae to advanced flowering plants. Evidence is provided that almost every vascular plant (tracheophyte), including the most primitive taxa (clubmosses, ferns and gymnosperms) contains caulimovirid EVEs, many of which represent previously unrecognized evolutionary branches. In angiosperms, EVEs from at least one and as many as five different caulimovirid genera were frequently detected, and florendoviruses were the most widely distributed, followed by petuviruses. From the analysis of the distribution of different caulimovirid genera within different plant species, we propose a working evolutionary scenario in which this family of viruses emerged at latest during Devonian era (approx. 320 million years ago) followed by vertical transmission and by several cross-division host swaps.

Identification and characterization of privet leaf blotch-associated virus, a novel idaeovirus

A novel virus has been identified by next-generation sequencing (NGS) in privet (Ligustrum japonicum L.) affected by a graft-transmissible disease characterized by leaf blotch symptoms resembling infectious variegation, a virus-like privet disease with an unclear aetiology. This virus, which has been tentatively named 'privet leaf blotch-associated virus' (PrLBaV), was absent in non-symptomatic privet plants, as revealed by NGS and reverse transcription-polymerase chain reaction (RT-PCR). Molecular characterization of PrLBaV showed that it has a segmented genome composed of two positive single-stranded RNAs, one of which (RNA1) is monocistronic and codes for the viral replicase, whereas the other (RNA2) contains two open reading frames (ORFs), ORF2a and ORF2b, coding for the putative movement (p38) and coat (p30) proteins, respectively. ORF2b is very probably expressed through a subgenomic RNA starting with six nucleotides (AUAUCU) that closely resemble those found in the 5'-terminal end of genomic RNA1 and RNA2 (AUAUUU and AUAUAU, respectively). The molecular signatures identified in the PrLBaV RNAs and proteins resemble those of Raspberry bushy dwarf virus (RBDV), currently the only member of the genus Idaeovirus. These data, together with phylogenetic analyses, are consistent with the proposal of considering PrLBaV as a representative of the second species in the genus Idaeovirus. Transient expression of a recombinant PrLBaV p38 fused to green fluorescent protein in leaves of Nicotiana benthamiana, coupled with confocal laser scanning microscopy assays, showed that it localizes at cell plasmodesmata, strongly supporting its involvement in viral movement/trafficking and providing the first functional characterization of an idaeovirus encoded protein.

A novel block of plant virus movement genes

Hibiscus green spot virus (HGSV) is a recently discovered and so far poorly characterized bacilliform plant virus with a positive-stranded RNA genome consisting of three RNA species. Here, we demonstrate that the proteins encoded by the ORF2 and ORF3 in HGSV RNA2 are necessary and sufficient to mediate cell-to-cell movement of transport-deficient Potato virus X in Nicotiana benthamiana. These two genes represent a specialized transport module called a 'binary movement block' (BMB), and ORF2 and ORF3 are termed BMB1 and BMB2 genes. In agroinfiltrated epidermal cells of N. benthamiana, green fluorescent protein (GFP)-BMB1 fusion protein was distributed diffusely in the cytoplasm and the nucleus. However, in the presence of BMB2, GFP-BMB1 was directed to cell wall-adjacent elongated bodies at the cell periphery, to cell wall-embedded punctate structures co-localizing with callose deposits at plasmodesmata, and to cells adjacent to the initially transformed cell. Thus, BMB2 can mediate the transport of BMB1 to and through plasmodesmata. In general, our observations support the idea that cell-to-cell trafficking of movement proteins involves an initial delivery to membrane compartments adjacent to plasmodesmata, subsequent entry of the plasmodesmata cavity and, finally, transport to adjacent cells. This process, as an alternative to tubule-based transport, has most likely evolved independently in triple gene block (TGB), double gene block (DGB), BMB and the single gene-coded transport system.

Movement protein of Apple chlorotic leaf spot virus is genetically unstable and negatively regulated by Ribonuclease E in E. coli

Movement protein (MP) of Apple chlorotic leaf spot virus (ACLSV) belongs to “30 K” superfamily of proteins and members of this family are known to show a wide array of functions. In the present study this gene was found to be genetically unstable in E. coli when transformed DH5α cells were grown at 28 °C and 37 °C. However, genetic instability was not encountered at 20 °C. Heterologous over expression failed despite the use of different transcriptional promoters and translational fusion constructs. Total cell lysate when subjected to western blotting using anti-ACLSV MP antibodies, showed degradation/cleavage of the expressed full-length protein. This degradation pointed at severe proteolysis or instability of the corresponding mRNA. Predicted secondary structure analysis of the transcript revealed a potential cleavage site for an endoribonuclease (RNase E) of E. coli. The negating effect of RNase E on transcript stability and expression was confirmed by northern blotting and quantitative RT-PCR of the RNA extracted from RNase E temperature sensitive mutant (strain N3431). The five fold accumulation of transcripts at non-permissive temperature (43 °C) suggests the direct role of RNase E in regulating the expression of ACLSV MP in E. coli.

The functional analysis of distinct tospovirus movement proteins (NSM) reveals different capabilities in tubule formation, cell-to-cell and systemic virus movement among the tospovirus species

The lack of infectious tospovirus clones to address reverse genetic experiments has compromised the functional analysis of viral proteins. In the present study we have performed a functional analysis of the movement proteins (NS_M) of four tospovirus species Bean necrotic mosaic virus (BeNMV), Chrysanthemum stem necrosis virus (CSNV), Tomato chlorotic spot virus (TCSV) and Tomato spotted wilt virus (TSWV), which differ biologically and molecularly, by using the Alfalfa mosaic virus (AMV) model system. All NS_M proteins were competent to: i) support the cell-to-cell and systemic transport of AMV, ii) generate tubular structures on infected protoplast and iii) transport only virus particles. However, the NS_M of BeNMV (one of the most phylogenetically distant species) was very inefficient to support the systemic transport. Deletion assays revealed that the C-terminal region of the BeNMV NS_M, but not that of the CSNV, TCSV and TSWV NS_M proteins, was dispensable for cell-to-cell transport, and that all the non-functional C-terminal NS_M mutants were unable to generate tubular structures. Bimolecular fluorescence complementation analysis revealed that the C-terminus of the BeNMV NS_M was not required for the interaction with the cognate nucleocapsid protein, showing a different protein organization when compared with other movement proteins of the '30K family'. Overall, our results revealed clearly differences in functional aspects among movement proteins from divergent tospovirus species that have a distinct biological behavior.

Bioinformatic and mutational analysis of ophiovirus movement proteins, belonging to the 30K superfamily

Ophioviridae is a family of segmented, negative-sense, single-stranded RNA plant viruses. We showed that their cell-to-cell movement protein (MP) is an isolated member of the 30K MP superfamily with a unique structural organization. All 30K MPs share a core domain that contains a nearly-invariant signature aspartate. We examined its role in the MP of Citrus psorosis virus (CPsV) and Mirafiori lettuce big-vein virus (MiLBVV). Alanine substitution of this aspartate prevented plasmodesmata accumulation of MP(MiLBVV), while MP(CPsV) was not affected. The capacity of ophiovirus MPs to increase the plasmodesmata size exclusion limit and non-cell autonomous protein feature was abolished in both mutants. To investigate the role of the signature aspartate in cell-to-cell movement, we constructed a new movement-deficient Tobacco mosaic virus vector used for trans-complementation assays. We showed that both ophiovirus MP mutants lack the cell-to-cell movement capacity, confirming that this signature aspartate is essential for viral cell-to-cell movement.

Identification of Ourmiavirus 30K movement protein amino acid residues involved in symptomatology, viral movement, subcellular localization and tubule formation

Several plant viruses encode movement proteins (MPs) classified in the 30K superfamily. Despite a great functional diversity, alignment analysis of MP sequences belonging to the 30K superfamily revealed the presence of a central core region, including amino acids potentially critical for MP structure and functionality. We performed alanine-scanning mutagenesis of the Ourmia melon virus (OuMV) MP, and studied the effects of amino acid substitutions on MP properties and virus infection. We identified five OuMV mutants that were impaired in systemic infection in Nicotiana benthamiana and Arabidopsis thaliana, and two mutants showing necrosis and pronounced mosaic symptoms, respectively, in N. benthamiana. Green fluorescent protein fusion constructs (GFP:MP) of movement-defective MP alleles failed to localize in distinct foci at the cell wall, whereas a GFP fusion with wild-type MP (GFP:MPwt) mainly co-localized with plasmodesmata and accumulated at the periphery of epidermal cells. The movement-defective mutants also failed to produce tubular protrusions in protoplasts isolated from infected leaves, suggesting a link between tubule formation and the ability of OuMV to move. In addition to providing data to support the importance of specific amino acids for OuMV MP functionality, we predict that these conserved residues might be critical for the correct folding and/or function of the MP of other viral species in the 30K superfamily.

Changes in the composition of the RNA virome mark evolutionary transitions in green plants

BACKGROUND

The known plant viruses mostly infect angiosperm hosts and have RNA or small DNA genomes. The only other lineage of green plants with a relatively well-studied virome, unicellular chlorophyte algae, is mostly infected by viruses with large DNA genomes. Thus RNA viruses and small DNA viruses seem to completely displace large DNA virus genomes in late branching angiosperms. To understand better the expansion of RNA viruses in the taxonomic span between algae and angiosperms, we analyzed the transcriptomes of 66 non-angiosperm plants characterized by the 1000 Plants Genomes Project.

RESULTS

We found homologs of virus RNA-dependent RNA polymerases in 28 non-angiosperm plant species, including algae, mosses, liverworts (Marchantiophyta), hornworts (Anthocerotophyta), lycophytes, a horsetail Equisetum, and gymnosperms. Polymerase genes in algae were most closely related to homologs from double-stranded RNA viruses leading latent or persistent lifestyles. Land plants, in addition, contained polymerases close to the homologs from single-stranded RNA viruses of angiosperms, capable of productive infection and systemic spread. For several polymerases, a cognate capsid protein was found in the same library. Another virus hallmark gene family, encoding the 30 K movement proteins, was found in lycophytes and monilophytes but not in mosses or algae.

CONCLUSIONS

The broadened repertoire of RNA viruses suggests that colonization of land and growth in anatomical complexity in land plants coincided with the acquisition of novel sets of viruses with different strategies of infection and reproduction.

Phylogenetic and Molecular Variability Studies Reveal a New Genetic Clade of Citrus leprosis virus C

Citrus leprosis virus C (CiLV-C) causes a severe disease affecting citrus orchards in the Western hemisphere. This study reveals the molecular variability of the virus by analyzing four genomic regions (p29, p15, MP and RNA2-intergenic region) distributed over its two RNAs. Nucleotide diversity (π) values were relatively low but statistically different over the analyzed genes and subpopulations, indicating their distinct evolutionary history. Values of πp29 and πMP were higher than those of πp15 and πRNA2-IR, whereas πMP was increased due to novel discovered isolates phylogenetically clustered in a divergent clade that we called SJP. Isolate BR_SP_SJP_01 RNA1 and RNA2 sequences, clade SJP, showed an identity of 85.6% and 88.4%, respectively, with those corresponding to CiLV-C, the type member of the genus Cilevirus, and its RNA2 5'-proximal region was revealed as a minor donor in a putative inter-clade recombination event. In addition to citrus, BR_SP_SJP_01 naturally infects the weed Commelina benghalensis and is efficiently transmitted by Brevipalpus yothersi mites. Our data demonstrated that negative selection was the major force operating in the evaluated viral coding regions and defined amino acids putatively relevant for the biological function of cilevirus proteins. This work provides molecular tools and sets up a framework for further epidemiological studies.

Cytorhabdovirus P3 genes encode 30K-like cell-to-cell movement proteins

Plant viruses encode movement proteins (MP) to facilitate cell-to-cell transport through plasmodesmata. In this study, using trans-complementation of a movement-defective turnip vein-clearing tobamovirus (TVCV) replicon, we show for the first time for cytorhabdoviruses (lettuce necrotic yellows virus (LNYV) and alfalfa dwarf virus (ADV)) that their P3 proteins function as MP similar to the TVCV P30 protein. All three MP localized to plasmodesmata when ectopically expressed. In addition, we show that these MP belong to the 30K superfamily since movement was inhibited by mutation of an aspartic acid residue in the critical 30K-specific LxD/N50-70G motif. We also report that Nicotiana benthamiana microtubule-associated VOZ1-like transcriptional activator interacts with LNYV P3 and TVCV P30 but not with ADV P3 or any of the MP point mutants. This host protein, which is known to interact with P3 of sonchus yellow net nucleorhabdovirus, may be involved in aiding the cell-to-cell movement of LNYV and TVCV.

Acquisition of Full-Length Viral Helicase Domains by Insect Retrotransposon-Encoded Polypeptides

Recent metagenomic studies in insects identified many sequences unexpectedly closely related to plant virus genes. Here we describe a new example of this kind, insect R1 LINEs with an additional C-terminal domain in their open reading frame 2. This domain is similar to NTPase/helicase (SF1H) domains, which are found in replicative proteins encoded by plant viruses of the genus Tobamovirus. We hypothesize that the SF1H domain could be acquired by LINEs, directly or indirectly, upon insect feeding on virus-infected plants. Possible functions of this domain in LINE transposition and involvement in LINEs counteraction the silencing-based cell defense against retrotransposons are discussed.

Sequence and Structure Analysis of Distantly-Related Viruses Reveals Extensive Gene Transfer between Viruses and Hosts and among Viruses

The origin and evolution of viruses is a subject of ongoing debate. In this study, we provide a full account of the evolutionary relationships between proteins of significant sequence and structural similarity found in viruses that belong to different classes according to the Baltimore classification. We show that such proteins can be found in viruses from all Baltimore classes. For protein families that include these proteins, we observe two patterns of the taxonomic spread. In the first pattern, they can be found in a large number of viruses from all implicated Baltimore classes. In the other pattern, the instances of the corresponding protein in species from each Baltimore class are restricted to a few compact clades. Proteins with the first pattern of distribution are products of so-called viral hallmark genes reported previously. Additionally, this pattern is displayed by the envelope glycoproteins from Flaviviridae and Bunyaviridae and helicases of superfamilies 1 and 2 that have homologs in cellular organisms. The second pattern can often be explained by horizontal gene transfer from the host or between viruses, an example being Orthomyxoviridae and Coronaviridae hemagglutinin esterases. Another facet of horizontal gene transfer comprises multiple independent introduction events of genes from cellular organisms into otherwise unrelated viruses.

In vitro properties of hordeivirus TGB1 protein forming ribonucleoprotein complexes

Hordeivirus movement protein encoded by the first gene of the triple gene block (TGB1 protein, TGBp1) interacts in vivo with viral genomic and subgenomic RNAs to form ribonucleoprotein (RNP) particles that are considered to be a form of viral genome (non-virion transport form) capable of cell-to-cell and long-distance transport in infected plants. The structures of these RNPs have not been elucidated. The poa semilatent virus (PSLV) TGBp1 contains a structured C-terminal NTPase/helicase domain and an N-terminal extension region consisting of two domains - a completely intrinsically disordered extreme N-terminal domain and an internal domain (ID) with structure resembling a partially disordered molten globule. Here, we characterized the structures assembled in vitro by the full-length PSLV TGBp1 alone or in the presence of viral RNA. The PSLV TGBp1 was capable of multimerization and self-assembly into extended high-molecular-mass complexes. These complexes disassembled to apparent monomers upon incubation with ATP. Upon incubation with viral RNA, the PSLV TGBp1 in vitro formed RNP structures that appeared as filamentous particles resembling virions of helical filamentous plant viruses in morphology and dimensions. By comparing the biophysical characteristics of PSLV TGBp1 and its domains in the presence and absence of RNA, we show that the ID plays the main structural role in the self-interactions and RNA interactions of TGBp1 leading to the assembly of virus-like RNP particles.