RNA-RNA recombination in plant virus replication and evolution

The human microbiota is a beneficial reservoir for SARS-CoV-2 mutations

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) mutations are rapidly emerging. In particular, beneficial mutations in the spike (S) protein, which can either make a person more infectious or enable immunological escape, are providing a significant obstacle to the prevention and treatment of pandemics. However, how the virus acquires a high number of beneficial mutations in a short time remains a mystery. We demonstrate here that variations of concern may be mutated due in part to the influence of the human microbiome. We searched the National Center for Biotechnology Information database for homologous fragments (HFs) after finding a mutation and the six neighboring amino acids in a viral mutation fragment. Among the approximate 8,000 HFs obtained, 61 mutations in S and other outer membrane proteins were found in bacteria, accounting for 62% of all mutation sources, which is 12-fold higher than the natural variable proportion. A significant proportion of these bacterial species-roughly 70%-come from the human microbiota, are mainly found in the lung or gut, and share a composition pattern with COVID-19 patients. Importantly, SARS-CoV-2 RNA-dependent RNA polymerase replicates corresponding bacterial mRNAs harboring mutations, producing chimeric RNAs. SARS-CoV-2 may collectively pick up mutations from the human microbiota that change the original virus's binding sites or antigenic determinants. Our study clarifies the evolving mutational mechanisms of SARS-CoV-2.

IMPORTANCE

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) mutations are rapidly emerging, in particular advantageous mutations in the spike (S) protein, which either increase transmissibility or lead to immune escape and are posing a major challenge to pandemic prevention and treatment. However, how the virus acquires a high number of advantageous mutations in a short time remains a mystery. Here, we provide evidence that the human microbiota is a reservoir of advantageous mutations and aids mutational evolution and host adaptation of SARS-CoV-2. Our findings demonstrate a conceptual breakthrough on the mutational evolution mechanisms of SARS-CoV-2 for human adaptation. SARS-CoV-2 may grab advantageous mutations from the widely existing microorganisms in the host, which is undoubtedly an "efficient" manner. Our study might open a new perspective to understand the evolution of virus mutation, which has enormous implications for comprehending the trajectory of the COVID-19 pandemic.

Characterization of Variant RNAs Encapsidated during Bromovirus Infection by High-Throughput Sequencing

Previously, we described the RNA recombinants accumulating in tissues infected with the bromoviruses BMV (Brome mosaic virus) and CCMV (Cowpea chlorotic mottle virus). In this work, we characterize the recombinants encapsidated inside the purified virion particles of BMV and CCMV. By using a tool called the Viral Recombination Mapper (ViReMa) that detects recombination junctions, we analyzed a high number of high-throughput sequencing (HTS) short RNA sequence reads. Over 28% of BMV or CCMV RNA reads did not perfectly map to the viral genomes. ViReMa identified 1.40% and 1.83% of these unmapped reads as the RNA recombinants, respectively, in BMV and CCMV. Intra-segmental crosses were more frequent than the inter-segmental ones. Most intra-segmental junctions carried short insertions/deletions (indels) and caused frameshift mutations. The mutation hotspots clustered mainly within the open reading frames. Substitutions of various lengths were also identified, whereas a small fraction of crosses occurred between viral and their host RNAs. Our data reveal that the virions can package detectable amounts of multivariate recombinant RNAs, contributing to the flexible nature of the viral genomes.

A novel chimeric RNA originating from BmCPV S4 and Bombyx mori HDAC11 transcripts regulates virus proliferation

Polymerases encoded by segmented negative-strand RNA viruses cleave 5'-m7G-capped host transcripts to prime viral mRNA synthesis ("cap-snatching") to generate chimeric RNA, and trans-splicing occurs between viral and cellular transcripts. Bombyx mori cytoplasmic polyhedrosis virus (BmCPV), an RNA virus belonging to Reoviridae, is a major pathogen of silkworm (B. mori). The genome of BmCPV consists of 10 segmented double-stranded RNAs (S1-S10) from which viral RNAs encoding a protein are transcribed. In this study, chimeric silkworm-BmCPV RNAs, in which the sequence derived from the silkworm transcript could fuse with both the 5' end and the 3' end of viral RNA, were identified in the midgut of BmCPV-infected silkworms by RNA_seq and further confirmed by RT-PCR and Sanger sequencing. A novel chimeric RNA, HDAC11-S4 RNA 4, derived from silkworm histone deacetylase 11 (HDAC11) and the BmCPV S4 transcript encoding viral structural protein 4 (VP4), was selected for validation by in situ hybridization and Northern blotting. Interestingly, our results indicated that HDAC11-S4 RNA 4 was generated in a BmCPV RNA-dependent RNA polymerase (RdRp)-independent manner and could be translated into a truncated BmCPV VP4 with a silkworm HDAC11-derived N-terminal extension. Moreover, it was confirmed that HDAC11-S4 RNA 4 inhibited BmCPV proliferation, decreased the level of H3K9me3 and increased the level of H3K9ac. These results indicated that during infection with BmCPV, a novel mechanism, different from that described in previous reports, allows the genesis of chimeric silkworm-BmCPV RNAs with biological functions.

Phylogeny and Genetic Divergence among Sorghum Mosaic Virus Isolates Infecting Sugarcane

Sorghum mosaic virus (SrMV, the genus Potyvirus of the family Potyviridae) is a causal agent of common mosaic in sugarcane and poses a threat to the global sugar industry. In this study, a total of 901 sugarcane leaf samples with mosaic symptom were collected from eight provinces in China and were detected via RT-PCR using a primer pair specific to the SrMV coat protein (CP). These leaf samples included 839 samples from modern cultivars (Saccharum spp. hybrids) and 62 samples from chewing cane (S. officinarum). Among these, 632 out of 901 (70.1%) samples were tested positive for SrMV. The incidences of SrMV infection were 72.3% and 40.3% in modern cultivars and chewing cane, respectively. Phylogenetic analysis showed that all tested SrMV isolates were clustered into three clades consisting of six phylogenetic groups based on 306 CP sequences (this study = 265 and GenBank database = 41). A total of 10 SrMV isolates from South America (the United States and Argentina) along with 106 isolates from China were clustered in group D, while the remaining 190 SrMV isolates from Asia (China and Vietnam) were dispersed in five groups. The SrMV isolates in group F were limited to Yunnan province in China, and those in group A were spread over eight provinces. A significant genetic heterogeneity was elucidated in the nucleotide sequence identities of all SrMV CPs, ranging from 69.0% to 100%. A potential recombination event was postulated among SrMV isolates based on CP sequences. All tested SrMV CPs underwent dominant negative selection. Geographical isolation (South America vs. Asia) and host types (modern cultivars vs. chewing cane) are important factors promoting the genetic differentiation of SrMV populations. Overall, this study contributes to the global understanding of the genetic evolution of SrMV and provides a valuable resource for the epidemiology and management of the mosaic in sugarcane.

Mixed infections with new emerging viruses associated with jujube mosaic disease

BACKGROUND

Jujube is an economically important fruit tree and native to China. Viral disease is a new threat to jujube production, and several new viruses have been identified infecting jujube plants. During our field survey, jujube mosaic disease was widely distributed in Beijing, but the associated causal agents are still unknown.

METHODS

Small RNA deep sequencing was conducted to identify the candidate viruses associated with jujube mosaic. Further complete genome sequences of the viruses were cloned, and the genomic characterization of each virus was analyzed. The field distribution of these viruses was further explored with PCR/RT-PCR detection of field samples.

RESULTS

Mixed infection of four viruses was identified in a plant sample with the symptom of mosaic and leaf twisting, including the previously reported jujube yellow mottle-associated virus (JYMaV), persimmon ampelovirus (PAmpV), a new badnavirus tentatively named jujube-associated badnavirus (JaBV), and a new secovirus tentatively named jujube-associated secovirus (JaSV). PAmpV-jujube was 14,093 nt in length with seven putative open reading frames (ORFs) and shared highest (79.4%) nucleotide (nt) sequence identity with PAmpV PBs3. Recombination analysis showed that PAmpV-jujube was a recombinant originating from plum bark necrosis stem pitting-associated virus isolates nanjing (KC590347) and bark (EF546442). JaBV was 6449 bp in length with conserved genomic organization typical of badnaviruses. The conserved RT and RNAse H region shared highest 67.6% nt sequence identity with jujube mosaic-associated virus, which was below the 80% nt sequence identity value used as the species demarcation threshold in Badnavirus. The genome of JaSV composed of two RNA molecules of 5878 and 3337 nts in length, excluding the polyA tails. Each genome segment contained one large ORF that shared homology and phylogenetic identity with members of the family Secoviridae. Field survey showed JYMaV and JaBV were widely distributed in jujube trees in Beijing.

CONCLUSION

Two new viruses were identified from jujube plants, and mixed infections of JYMaV and JaBV were common in jujube in Beijing.

Zucchini Yellow Mosaic Virus (ZYMV) as a Serious Biotic Stress to Cucurbits: Prevalence, Diversity, and Its Implications for Crop Sustainability

Zucchini yellow mosaic virus (ZYMV) is a severe threat to cucurbit crops worldwide, including Pakistan. This study was pursued to evaluate the prevalence, geographic distribution, and molecular diversity of ZYMV isolates infecting cucurbits in Pakistan's Pothwar region. Almost all the plant viruses act as a biotic stress on the host plants, which results in a yield loss. These viruses cause losses in single-infection or in mixed-infection cucurbit crops, and we have found a number of mixed-infected samples belonging to the Curubitaceae family. Serological detection of the tested potyviruses in the collected cucurbit samples revealed that ZYMV was the most prevalent virus, with a disease incidence (DI) at 35.2%, followed by Papaya ringspot virus (PRSV) with an incidence of 2.2%, and Watermelon mosaic virus (WMV) having an incidence as little as 0.5% in 2016. In the year 2017, a relatively higher disease incidence of 39.7%, 2.4%, and 0.3% for ZYMV, WMV, and PRSV, respectively, was recorded. ZYMV was the most prevalent virus with the highest incidence in Attock, Rawalpindi, and Islamabad, while PRSV was observed to be the highest in Islamabad and Jhelum. WMV infection was observed only in Rawalpindi and Chakwal. Newly detected Pakistani ZYMV isolates shared 95.8-97.0% nucleotide identities among themselves and 77.1-97.8% with other isolates retrieved from GenBank. Phylogenetic relationships obtained using different ZYMV isolates retrieved from GenBank and validated by in silico restriction analysis revealed that four Pakistani isolates clustered with other ZYMV isolates in group IIb with Chinese, Italian, Polish, and French isolates, while another isolate (MK848239) formed a separate minor clade within IIb. The isolate MK8482490, reported to infect bitter gourd in Pakistan, shared a minor clade with a Chinese isolate (KX884570). Recombination analysis revealed that the recently found ZYMV isolate (MK848239) is most likely a recombinant of Pakistani (MK848237) and Italian (MK956829) isolates, with a recombinant breakpoint between 266 and 814 nucleotide positions. Local isolate comparison and recombination detection may aid in the development of a breeding program that identifies resistant sources against recombinant isolates because the ZYMV is prevalent in a few cucurbit species grown in the surveyed areas and causes heavy losses and economic damage to the agricultural community.

Virus-virus interactions alter the mechanical transmissibility and host range of begomoviruses

Introduction

Begomoviruses are mainly transmitted by whiteflies. However, a few begomoviruses can be transmitted mechanically. Mechanical transmissibility affects begomoviral distribution in the field.

Materials and methods

In this study, two mechanically transmissible begomoviruses, tomato leaf curl New Delhi virus-oriental melon isolate (ToLCNDV-OM) and tomato yellow leaf curl Thailand virus (TYLCTHV), and two nonmechanically transmissible begomoviruses, ToLCNDV-cucumber isolate (ToLCNDV-CB) and tomato leaf curl Taiwan virus (ToLCTV), were used to study the effects of virus-virus interactions on mechanical transmissibility.

Results

Nicotiana benthamiana and host plants were coinoculated through mechanical transmission with inoculants derived from plants that were mix-infected or inoculants derived from individually infected plants, and the inoculants were mixed immediately before inoculation. Our results showed that ToLCNDV-CB was mechanically transmitted with ToLCNDV-OM to N. benthamiana, cucumber, and oriental melon, whereas ToLCTV was mechanically transmitted with TYLCTHV to N. benthamiana and tomato. For crossing host range inoculation, ToLCNDV-CB was mechanically transmitted with TYLCTHV to N. benthamiana and its nonhost tomato, while ToLCTV with ToLCNDV-OM was transmitted to N. benthamiana and its nonhost oriental melon. For sequential inoculation, ToLCNDV-CB and ToLCTV were mechanically transmitted to N. benthamiana plants that were either preinfected with ToLCNDV-OM or TYLCTHV. The results of fluorescence resonance energy transfer analyses showed that the nuclear shuttle protein of ToLCNDV-CB (CBNSP) and the coat protein of ToLCTV (TWCP) localized alone to the nucleus. When coexpressed with movement proteins of ToLCNDV-OM or TYLCTHV, CBNSP and TWCP relocalized to both the nucleus and the cellular periphery and interacted with movement proteins.

Discussion

Our findings indicated that virus-virus interactions in mixed infection circumstances could complement the mechanical transmissibility of nonmechanically transmissible begomoviruses and alter their host range. These findings provide new insight into complex virus-virus interactions and will help us to understand the begomoviral distribution and to reevaluate disease management strategies in the field.

Molecular Variation and Genomic Function of Citrus Vein Enation Virus

In this study, we identified a new citrus vein enation virus (CVEV) isolate (named CVEV-DT1) through sRNA high-throughput sequencing and traditional sequencing. Phylogenetic analysis based on whole genome sequences of all known CVEV isolates revealed that CVEV-DT1 was in an evolutionary branch with other isolates from China. Molecular variation analysis showed that the single nucleotide variability along CVEV full-length sequences was less than 8%, with more transitions (60.55%) than transversions (39.43%), indicating a genetically homogeneous CVEV population. In addition, non-synonymous nucleotide mutations mainly occurred in ORF1 and ORF2. Based on disorder analysis of all encoded ORF by CVEV-DT1, we identified that the CVEV-DT1 coat protein (CP) formed spherical granules, mainly in the cell nucleus and partly throughout the cytoplasm, with liquid properties through subcellular localization and photobleaching assay. Furthermore, we also confirmed that the CVEV P0 protein has weak post-transcriptional RNA-silencing suppressor activity and could elicit a strong hypersensitive response (HR) in tobacco plants. Collectively, to the best of our knowledge, our study was the first to profile the genomic variation in all the reported CVEV isolates and reveal the functions of CVEV-DT1-encoded proteins.

Multiple Levels of Triggered Factors and the Obligated Requirement of Cell-to-Cell Movement in the Mutation Repair of Cucumber Mosaic Virus with Defects in the tRNA-like Structure

Simple Summary

Based on analysis of the tRNA-like structure (TLS) mutation in cucumber mosaic virus (CMV), mutation repair is correlated with several levels of triggered factors, including the dose of inoculation of virus mutants, the quantity effect on corresponding viral RNA, and the quality effect on corresponding viral RNA. All types of TLS mutation in different RNAs of CMV can be repaired at a low dose around the dilution end-point. At a high dose of inoculation, TLS mutations in RNA2 and RNA3, but not RNA1, can be repaired, which correlates with the relative quantity defect of RNA2 or the genome size defect of RNA3. In addition, all the above types of mutation repair necessarily require cell-to-cell movement, which presents the obligated effect of cell-to-cell movement on mutation repair.

Abstract

Some debilitating mutations in RNA viruses are repairable; however, the triggering factors of mutation repair remain largely unknown. In this study, multiple triggering factors of mutation repair are identified based on genetic damage to the TLS in CMV. TLS mutations in different RNAs distinctively impact viral pathogenicity and present different types of mutation repair. RNA2 relative reduction level or RNA3 sequence change resulting from TLS mutation is correlated with a high rate of mutation repair, and the TLS mutation of RNA1 fails to be repaired at the high inoculum dose. However, the TLS mutation of RNA1 can be repaired at a low dose of inoculation, particularly around the dilution end-point or in the mixed inoculation with RNA2 having a pre-termination mutation of the 2b gene, an RNAi suppressor. Taken together, TLS mutations resulting in quality or quantity defects of the viral genome or TLS mutations at low doses around the dilution end-point are likely to be repaired. Different levels of TLS mutation repair necessarily require cell-to-cell movement, therefore implying its obligated effect on the evolution of low-fitness viruses and providing a new insight into Muller’s ratchet. This study provides important information on virus evolution and the application of mild viral vaccines.

Recombination in Positive-Strand RNA Viruses

RNA recombination is a major driver of genetic shifts tightly linked to the evolution of RNA viruses. Genomic recombination contributes substantially to the emergence of new viral lineages, expansion in host tropism, adaptations to new environments, and virulence and pathogenesis. Here, we review some of the recent progress that has advanced our understanding of recombination in positive-strand RNA viruses, including recombination triggers and the mechanisms behind them. The study of RNA recombination aids in predicting the probability and outcome of viral recombination events, and in the design of viruses with reduced recombination frequency as candidates for the development of live attenuated vaccines. Surveillance of viral recombination should remain a priority in the detection of emergent viral strains, a goal that can only be accomplished by expanding our understanding of how these events are triggered and regulated.

Genomic Characterisation of an Isolate of Brassica Yellows Virus Associated with Brassica Weed in Tasmania

Brassica yellows virus (BrYV), a tentative species in the genus Polerovirus, of the Solemoviridae family, is a phloem-restricted and aphid-transmitted virus with at least three genotypes (A, B, and C). It has been found across mainland China, South Korea, and Japan. BrYV was previously undescribed in Tasmania, and its genetic variability in the state remains unknown. Here, we describe a near-complete genome sequence of BrYV (genotype A) isolated from Raphanus raphanistrum in Tasmania using next-generation sequencing and sanger sequencing of RT-PCR products. BrYV-Tas (GenBank Accession no. OM469309) possesses a genome of 5516 nucleotides (nt) and shares higher sequence identity (about 90%) with other BrYV isolates. Phylogenetic analyses showed variability in the clustering patterns of the individual genes of BrYV-Tas. Recombination analysis revealed beginning and ending breakpoints at nucleotide positions 1922 to 5234 nt, with the BrYV isolate LC428359 and BrYV isolate KY310572 identified as major and minor parents, respectively. Results of the evolutionary analysis showed that the majority of the codons for each gene are evolving under purifying selection, though a few codons were also detected to have positive selection pressure. Taken together, our findings will facilitate an understanding of the evolutionary dynamics and genetic diversity of BrYV.

Transcriptome and Small RNA Profiling of Potato Virus Y Infected Potato Cultivars, Including Systemically Infected Russet Burbank

Potatoes are the world’s most produced non-grain crops and an important food source for billions of people. Potatoes are susceptible to numerous pathogens that reduce yield, including Potato virus Y (PVY). Genetic resistance to PVY is a sustainable way to limit yield and quality losses due to PVY infection. Potato cultivars vary in their susceptibility to PVY and include susceptible varieties such as Russet Burbank, and resistant varieties such as Payette Russet. Although the loci and genes associated with PVY-resistance have been identified, the genes and mechanisms involved in limiting PVY during the development of systemic infections have yet to be fully elucidated. To increase our understanding of PVY infection, potato antiviral responses, and resistance, we utilized RNA sequencing to characterize the transcriptomes of two potato cultivars. Since transcriptional responses associated with the extreme resistance response occur soon after PVY contact, we analyzed the transcriptome and small RNA profile of both the PVY-resistant Payette Russet cultivar and PVY-susceptible Russet Burbank cultivar 24 h post-inoculation. While hundreds of genes, including terpene synthase and protein kinase encoding genes, exhibited increased expression, the majority, including numerous genes involved in plant pathogen interactions, were downregulated. To gain greater understanding of the transcriptional changes that occur during the development of systemic PVY-infection, we analyzed Russet Burbank leaf samples one week and four weeks post-inoculation and identified similarities and differences, including higher expression of genes involved in chloroplast function, photosynthesis, and secondary metabolite production, and lower expression of defense response genes at those time points. Small RNA sequencing identified different populations of 21- and 24-nucleotide RNAs and revealed that the miRNA profiles in PVY-infected Russet Burbank plants were similar to those observed in other PVY-tolerant cultivars and that during systemic infection ~32% of the NLR-type disease resistance genes were targeted by 21-nt small RNAs. Analysis of alternative splicing in PVY-infected potato plants identified splice variants of several hundred genes, including isoforms that were more dominant in PVY-infected plants. The description of the PVYN-Wi-associated transcriptome and small RNA profiles of two potato cultivars described herein adds to the body of knowledge regarding differential outcomes of infection for specific PVY strain and host cultivar pairs, which will help further understanding of the mechanisms governing genetic resistance and/or virus-limiting responses in potato plants.

New Insights into the Genome Organization of Yeast Double-Stranded RNA LBC Viruses

The yeasts Torulaspora delbrueckii (Td) and Saccharomyces cerevisiae (Sc) may show a killer phenotype that is encoded in dsRNA M viruses (V-M), which require the helper activity of another dsRNA virus (V-LA or V-LBC) for replication. Recently, two TdV-LBCbarr genomes, which share sequence identity with ScV-LBC counterparts, were characterized by high-throughput sequencing (HTS). They also share some similar characteristics with Sc-LA viruses. This may explain why TdV-LBCbarr has helper capability to maintain M viruses, whereas ScV-LBC does not. We here analyze two stretches with low sequence identity (LIS I and LIS II) that were found in TdV-LBCbarr Gag-Pol proteins when comparing with the homologous regions of ScV-LBC. These stretches may result from successive nucleotide insertions or deletions (indels) that allow compensatory frameshift events required to maintain specific functions of the RNA-polymerase, while modifying other functions such as the ability to bind V-M (+)RNA for packaging. The presence of an additional frameshifting site in LIS I may ensure the synthesis of a certain amount of RNA-polymerase until the new compensatory indel appears. Additional 5'- and 3'-extra sequences were found beyond V-LBC canonical genomes. Most extra sequences showed high identity to some stretches of the canonical genomes and can form stem-loop structures. Further, the 3'-extra sequence of two ScV-LBC genomes contains rRNA stretches. The origin and possible functions of these extra sequences are here discussed.

OUP accepted manuscript

Potyviridae, the largest family of known RNA viruses (realm Riboviria), belongs to the picorna-like supergroup and has important agricultural and ecological impacts. Potyvirid genomes are translated into polyproteins, which are in turn hydrolyzed to release mature products. Recent sequencing efforts revealed an unprecedented number of potyvirids with a rich variability in gene content and genomic layouts. Here, we review the heterogeneity of non-core modules that expand the structural and functional diversity of the potyvirid proteomes. We provide a family-wide classification of P1 proteinases into the functional Types A and B, and discuss pretty interesting sweet potato potyviral ORF (PISPO), putative zinc fingers, and alkylation B (AlkB)—non-core modules found within P1 cistrons. The atypical inosine triphosphate pyrophosphatase (ITPase/HAM1), as well as the pseudo tobacco mosaic virus-like coat protein (TMV-like CP) are discussed alongside homologs of unrelated virus taxa. Family-wide abundance of the multitasking helper component proteinase (HC-pro) is revised. Functional connections between non-core modules are highlighted to support host niche adaptation and immune evasion as main drivers of the Potyviridae evolutionary radiation. Potential biotechnological and synthetic biology applications of potyvirid leader proteinases and non-core modules are finally explored.

Viruses Without Borders: Global Analysis of the Population Structure, Haplotype Distribution, and Evolutionary Pattern of Iris Yellow Spot Orthotospovirus (Family Tospoviridae, Genus Orthotospovirus)

Iris yellow spot, caused by Iris yellow spot orthotospovirus (IYSV) (Genus: Orthotospovirus, Family: Tospoviridae), is an important disease of Allium spp. The complete N gene sequences of 142 IYSV isolates of curated sequence data from GenBank were used to determine the genetic diversity and evolutionary pattern. In silico restriction fragment length polymorphism (RFLP) analysis, codon-based maximum likelihood studies, genetic differentiation and gene flow within the populations of IYSV genotypes were investigated. Bayesian phylogenetic analysis was carried out to estimate the evolutionary rate. In silico RFLP analysis of N gene sequences categorized IYSV isolates into two major genotypes viz., IYSV Netherlands (IYSV_NL; 55.63%), IYSV Brazil (IYSV_BR; 38.73%) and the rest fell in neither group [IYSV other (IYSV_other; 5.63%)]. Phylogenetic tree largely corroborated the results of RFLP analysis and the IYSV genotypes clustered into IYSV_NL and IYSV_BR genotypes. Genetic diversity test revealed IYSV_other to be more diverse than IYSV_NL and IYSV_BR. IYSV_NL and IYSV_BR genotypes are under purifying selection and population expansion, whereas IYSV_other showed decreasing population size and hence appear to be under balancing selection. IYSV_BR is least differentiated from IYSV_other compared to IYSV_NL genotype based on nucleotide diversity. Three putative recombinant events were found in the N gene of IYSV isolates based on RDP analysis, however, RAT substantiated two among them. The marginal likelihood mean substitution rate was 5.08 × 10^–5 subs/site/year and 95% highest posterior density (HPD) substitution rate between 5.11 × 10^–5 and 5.06 × 10^–5. Findings suggest that IYSV continues to evolve using population expansion strategies. The substitution rates identified are similar to other plant RNA viruses.

Genetic Variability and Evidence of a New Subgroup in Watermelon Mosaic Virus Isolates

Watermelon mosaic virus (WMV) is one of the important Potyviruses that infect cucurbits worldwide. To better understand the population structure of WMV in the United States (U.S.), 57 isolates were collected from cucurbit fields located in nine southern states. The complete coat protein gene of all WMV isolates was cloned, sequenced and compared with 89 reported WMV isolates. The nucleotide and amino acid sequence identities among the U.S. WMV isolates ranged from 88.9 to 99.7% and from 91.5 to 100%, respectively. Phylogenetic analysis revealed that all the U.S. WMV isolates irrespective of their geographic origin or hosts belonged to Group 3. However, the fifty-seven isolates made three clusters in G3, where two clusters were similar to previously reported subgroups EM1 and EM2, and the third cluster, containing nine WMV isolates, formed a distinct subgroup named EM5 in this study. The ratio of non-synonymous to synonymous nucleotide substitution was low indicating the occurrence of negative purifying selection in the CP gene of WMV. Phylogenetic analysis of selected 37 complete genome sequences of WMV isolates also supported the above major grouping. Recombination analysis in the CP genes confirmed various recombinant events, indicating that purifying selection and recombination are the two dominant forces for the evolution of WMV isolates in the U.S.

P1 of Sweet Potato Feathery Mottle Virus Shows Strong Adaptation Capacity, Replacing P1-HCPro in a Chimeric Plum Pox Virus

Potyviridae is the largest family of plant RNA viruses. Their genomes are expressed through long polyproteins that are usually headed by the leader endopeptidase P1. This protein can be classified as type A or type B based on host proteolytic requirements and RNA silencing suppression (RSS) capacity. The main Potyviridae genus is Potyvirus, and a group of potyviruses infecting sweet potato presents an enlarged P1 protein with a polymerase slippage motif that produces an extra product termed P1N-PISPO. These two proteins display some RSS activity and are expressed followed by HCPro, which appears to be the main RNA silencing suppressor in these viruses. Here, we studied the behavior of the P1 protein of Sweet potato feathery mottle virus (SPFMV) using a viral system based on a canonical potyvirus, Plum pox virus (PPV), and discovered that this protein is able to replace both PPV P1 and HCPro. We also found that P1N-PISPO, produced after polymerase slippage, provides extra RNA silencing suppression capacity to SPFMV P1 in this viral context. In addition, the results showed that presence of two type A P1 proteins was detrimental for viral viability. The ample recombination spectrum that we found in the recovered viruses supports the strong adaptation capacity of P1 proteins and signals the N-terminal part of SPFMV P1 as essential for RSS activity. Further analyses provided data to add extra layers to the evolutionary history of sweet potato-infecting potyvirids. IMPORTANCE Plant viruses represent a major challenge for agriculture worldwide and Potyviridae, being the largest family of plant RNA viruses, is one of the primary players. P1, the leader endopeptidase, is a multifunctional protein that contributes to the successful spread of these viruses over a wide host range. Understanding how P1 proteins work, their dynamic interplay during viral infection, and their evolutionary path is critical for the development of strategic tools to fight the multiple diseases these viruses cause. We focused our efforts on the P1 protein of Sweet potato feathery mottle virus, which is coresponsible for the most devastating disease in sweet potato. The significance of our research is in understanding the capacity of this protein to perform several independent functions, using this knowledge to learn more about P1 proteins in general and the potyvirids infecting this host.

From RNA World to SARS-CoV-2: The Edited Story of RNA Viral Evolution

The current SARS-CoV-2 pandemic underscores the importance of understanding the evolution of RNA genomes. While RNA is subject to the formation of similar lesions as DNA, the evolutionary and physiological impacts RNA lesions have on viral genomes are yet to be characterized. Lesions that may drive the evolution of RNA genomes can induce breaks that are repaired by recombination or can cause base substitution mutagenesis, also known as base editing. Over the past decade or so, base editing mutagenesis of DNA genomes has been subject to many studies, revealing that exposure of ssDNA is subject to hypermutation that is involved in the etiology of cancer. However, base editing of RNA genomes has not been studied to the same extent. Recently hypermutation of single-stranded RNA viral genomes have also been documented though its role in evolution and population dynamics. Here, we will summarize the current knowledge of key mechanisms and causes of RNA genome instability covering areas from the RNA world theory to the SARS-CoV-2 pandemic of today. We will also highlight the key questions that remain as it pertains to RNA genome instability, mutations accumulation, and experimental strategies for addressing these questions.

Molecular Characterization of Potato Virus Y (PVY) Using High-Throughput Sequencing: Constraints on Full Genome Reconstructions Imposed by Mixed Infection Involving Recombinant PVY Strains

In recent years, high throughput sequencing (HTS) has brought new possibilities to the study of the diversity and complexity of plant viromes. Mixed infection of a single plant with several viruses is frequently observed in such studies. We analyzed the virome of 10 tomato and sweet pepper samples from Slovakia, all showing the presence of potato virus Y (PVY) infection. Most datasets allow the determination of the nearly complete sequence of a single-variant PVY genome, belonging to one of the PVY recombinant strains (N-Wi, NTNa, or NTNb). However, in three to-mato samples (T1, T40, and T62) the presence of N-type and O-type sequences spanning the same genome region was documented, indicative of mixed infections involving different PVY strains variants, hampering the automated assembly of PVY genomes present in the sample. The N- and O-type in silico data were further confirmed by specific RT-PCR assays targeting UTR-P1 and NIa genomic parts. Although full genomes could not be de novo assembled directly in this situation, their deep coverage by relatively long paired reads allowed their manual re-assembly using very stringent mapping parameters. These results highlight the complexity of PVY infection of some host plants and the challenges that can be met when trying to precisely identify the PVY isolates involved in mixed infection.

Genetic Variability and Molecular Evolution of Maize Yellow Mosaic Virus Populations from Different Geographic Origins

Maize yellow mosaic virus (MaYMV) hosted in various gramineous plants was assigned to the genus Polerovirus (family Luteoviridae) in 2018. However, little is known about its genetic diversity and population structure. In this study, 509 sugarcane leaf samples with mosaic symptoms were collected in 2017 to 2019 from eight sugarcane-growing provinces in China. Reverse-transcription PCR results revealed that four positive-sense RNA viruses were found to infect sugarcane, and the incidence of MaYMV among samples from Fujian, Sichuan, and Guangxi Provinces was 52.1, 9.8, and 2.5%, respectively. Based on 82 partial MaYMV sequences and 46 whole-genome sequences from different host plants, phylogenetic analysis revealed that MaYMV populations are very closely associated with their source geographical regions (China, Africa, and South America). Pairwise identity analysis showed significant variability in genome sequences among MaYMV isolates with genomic nucleotide identities of 91.1 to 99.9%. In addition to codon mutations, insertions or deletions also contributed to genetic variability in individual coding regions, especially in the readthrough protein (P3-P5 fusion protein). Low gene flow and significant genetic differentiation of MaYMV were observed among the three geographical populations, suggesting that environmental adaptation is an important evolutionary force that shapes the genetic structure of MaYMV. Genes in the MaYMV genome were subject to strong negative or purification selection during evolution, except for the movement protein (MP), which was under positive selection pressure. This finding suggests that the MP may play an important role in MaYMV evolution. Taken together, our findings provide basic information for the development of an integrated disease management strategy against MaYMV.

Different potato virus Y strains frequently co-localize in single epidermal leaf cells and in the aphid stylet

Single aphids can simultaneously or sequentially acquire and transmit multiple potato virus Y (PVY) strains. Multiple PVY strains are often found in the same field and occasionally within the same plant, but little is known about how PVY strains interact in plants or in aphid stylets. Immuno-staining and confocal microscopy were used to examine the spatial and temporal dynamics of PVY strain mixtures (PVYO and PVYNTN or PVYO and PVYN) in epidermal leaf cells of 'Samsun NN' tobacco and 'Goldrush' potato. Virus binding and localization was also examined in aphid stylets following acquisition. Both strains systemically infected tobacco and co-localized in cells of all leaves examined; however, the relative amounts of each virus changed over time. Early in the tobacco infection, when mosaic symptoms were observed, PVYO dominated the infection although PVYNTN was detected in some cells. As the infection progressed and vein necrosis developed, PVYNTN was prevalent. Co-localization of PVYO and PVYN was also observed in epidermal cells of potato leaves with most cells infected with both viruses. Furthermore, two strains could be detected binding to the distal end of aphid stylets following virus acquisition from a plant infected with a strain mixture. These data are in contrast with the traditional belief of spatial separation of two closely related potyviruses and suggest apparent non-antagonistic interaction between PVY strains that could help explain the multitude of emerging recombinant PVY strains discovered in potato in recent years.

Enhancing Capsid Proteins Capacity in Plant Virus-Vector Interactions and Virus Transmission

Vector transmission of plant viruses is basically of two types that depend on the virus helper component proteins or the capsid proteins. A number of plant viruses belonging to disparate groups have developed unusual capsid proteins providing for interactions with the vector. Thus, cauliflower mosaic virus, a plant pararetrovirus, employs a virion associated p3 protein, the major capsid protein, and a helper component for the semi-persistent transmission by aphids. Benyviruses encode a capsid protein readthrough domain (CP-RTD) located at one end of the rod-like helical particle, which serves for the virus transmission by soil fungal zoospores. Likewise, the CP-RTD, being a minor component of the luteovirus icosahedral virions, provides for persistent, circulative aphid transmission. Closteroviruses encode several CPs and virion-associated proteins that form the filamentous helical particles and mediate transmission by aphid, whitefly, or mealybug vectors. The variable strategies of transmission and evolutionary ‘inventions’ of the unusual capsid proteins of plant RNA viruses are discussed.

Bromoviruses (Bromoviridae)

The members of the family Bromoviridae have spherical or bacilliform virions with tri-segmented, single-stranded genomic RNAs, packaged in separate particles. Six genera including Alfamovirus, Anulavirus, Bromovirus, Cucumovirus, Ilarvirus, and Oleavirus are part of the family. RNA1 and RNA2 code for the replicase whereas RNA3 codes for movement and coat proteins. Genomic RNAs are infectious, but some species also require CP for infectivity. Members can encapsidate/accumulate sub-genomic RNAs, and/or defective or satellite RNAs. RNA replication occurs inside membranous spherules, and the role of host factors in RNA replication have been documented. Frequent RNA-RNA recombination and segment reassortment processes were observed among bromovirids. Transmission occurs mechanically, via pollen, seeds or insects. Host range varies from narrow to wide, infecting herbaceous plants, shrubs and trees, with some members causing major epidemics.

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.

Next-Generation Sequencing Combined With Conventional Sanger Sequencing Reveals High Molecular Diversity in Actinidia Virus 1 Populations From Kiwifruit Grown in China

Kiwifruit (Actinidia spp.) is native to China. Viral disease–like symptoms are common on kiwifruit plants. In this study, six libraries prepared from total RNA of leaf samples from 69 kiwifruit plants were subjected to next-generation sequencing (NGS). Actinidia virus 1 (AcV-1), a tentative species in the family Closteroviridae, was discovered in the six libraries. Two full-length and two near-full genome sequences of AcV-1 variants were determined by Sanger sequencing. The genome structure of these Chinese AcV-1 variants was identical to that of isolate K75 and consisted of 12 open reading frames (ORFs). Analyses of these sequences together with the NGS-derived contig sequences revealed high molecular diversity in AcV-1 populations, with the highest sequence variation occurring at ORF1a, ORF2, and ORF3, and the available variants clustered into three phylogenetic clades. For the first time, our study revealed different domain compositions in the viral ORF1a and molecular recombination events among AcV-1 variants. Specific reverse transcriptase–polymerase chain reaction assays disclosed the presence of AcV-1 in plants of four kiwifruit species and unknown Actinidia spp. in seven provinces and one city.

Genome Organization of a New Double-Stranded RNA LA Helper Virus From Wine Torulaspora delbrueckii Killer Yeast as Compared With Its Saccharomyces Counterparts

Wine killer yeasts such as killer strains of Torulaspora delbrueckii and Saccharomyces cerevisiae contain helper large-size (4.6 kb) dsRNA viruses (V-LA) required for the stable maintenance and replication of killer medium-size dsRNA viruses (V-M) which bear the genes that encode for the killer toxin. The genome of the new V-LA dsRNA from the T. delbrueckii Kbarr1 killer yeast (TdV-LAbarr1) was characterized by high-throughput sequencing (HTS). The canonical genome of TdV-LAbarr1 shares a high sequence identity and similar genome organization with its Saccharomyces counterparts. It contains all the known conserved motifs predicted to be necessary for virus translation, packaging, and replication. Similarly, the Gag-Pol amino-acid sequence of this virus contains all the features required for cap-snatching and RNA polymerase activity, as well as the expected regional variables previously found in other LA viruses. Sequence comparison showed that two main clusters (99.2-100% and 96.3-98.8% identity) include most LA viruses from Saccharomyces, with TdV-LAbarr1 being the most distant from all these viruses (61.5-62.5% identity). Viral co-evolution and cross transmission between different yeast species are discussed based on this sequence comparison. Additional 5' and 3' sequences were found in the TdV-LAbarr1 genome as well as in some newly sequenced V-LA genomes from S. cerevisiae. A stretch involving the 5' extra sequence of TdV-LAbarr1 is identical to a homologous stretch close to the 5' end of the canonical sequence of the same virus (self-identity). Our modeling suggests that these stretches can form single-strand stem loops, whose unpaired nucleotides could anneal to create an intramolecular kissing complex. Similar stem loops are also found in the 3' extra sequence of the same virus as well as in the extra sequences of some LA viruses from S. cerevisiae. A possible origin of these extra sequences as well as their function in obviating ssRNA degradation and allowing RNA transcription and replication are discussed.

Modulation of disease severity by plant positive-strand RNA viruses: The complex interplay of multifunctional viral proteins, subviral RNAs and virus-associated RNAs with plant signaling pathways and defense responses

Plant viruses induce a range of symptoms of varying intensity, ranging from severe systemic necrosis to mild or asymptomatic infection. Several evolutionary constraints drive virus virulence, including the dependence of viruses on host factors to complete their infection cycle, the requirement to counteract or evade plant antiviral defense responses and the mode of virus transmission. Viruses have developed an array of strategies to modulate disease severity. Accumulating evidence has highlighted not only the multifunctional role that viral proteins play in disrupting or highjacking plant factors, hormone signaling pathways and intracellular organelles, but also the interaction networks between viral proteins, subviral RNAs and/or other viral-associated RNAs that regulate disease severity. This review focusses on positive-strand RNA viruses, which constitute the majority of characterized plant viruses. Using well-characterized viruses with different genome types as examples, recent advances are discussed as well as knowledge gaps and opportunities for further research.

Genetic diversity of strawberry mild yellow edge virus from eastern Canada

Strawberry mild yellow edge virus (SMYEV) is a member of the genus Potexvirus, family Alphaflexiviridae. It is one of the most common pathogenic viruses infecting cultivated strawberries worldwide. In this study, we investigated the genetic diversity of SMYEV in strawberry fields that were severely affected by strawberry decline disease in the eastern Canadian provinces of New Brunswick, Nova Scotia, Prince Edward Island and Quebec. A total of 134 SMYEV coat protein (CP) gene sequences, representing 85 nucleic acid haplotypes, were identified in 56 field samples. A highly divergent SMYEV population was found in all four provinces, but there was little genetic differentiation among the populations, and moreover, the Canadian SMYEV isolates formed a unique dissimilar, genetically divergent population group when compared to those reported in other countries. Phylogenetic analysis revealed three new SMYEV subclades that consisted mainly of Canadian variants and were composed of 76 sequence haplotypes (76/85, 88%). Mixed infections by different SMYEV variants were observed in 38 samples (38/56, 68%). Evolutionary analysis suggested that the SMYEV strains in eastern Canada possibly originated outside of Canada but adapted to conditions in the region through genetic mutations.

Molecular basis of genetic variation of viruses

Genetic variation is a necessity of all biological systems. Viruses use all known mechanisms of variation; mutation, several forms of recombination, and segment reassortment in the case of viruses with a segmented genome. These processes are intimately connected with the replicative machineries of viruses, as well as with fundamental physical-chemical properties of nucleotides when acting as template or substrate residues. Recombination has been viewed as a means to rescue viable genomes from unfit parents or to produce large modifications for the exploration of phenotypic novelty. All types of genetic variation can act conjointly as blind processes to provide the raw materials for adaptation to the changing environments in which viruses must replicate. A distinction is made between mechanistically unavoidable and evolutionarily relevant mutation and recombination.

Detection of Plant Viruses and Disease Management: Relevance of Genetic Diversity and Evolution

Plant viruses cause considerable economic losses and are a threat for sustainable agriculture. The frequent emergence of new viral diseases is mainly due to international trade, climate change, and the ability of viruses for rapid evolution. Disease control is based on two strategies: i) immunization (genetic resistance obtained by plant breeding, plant transformation, cross-protection, or others), and ii) prophylaxis to restrain virus dispersion (using quarantine, certification, removal of infected plants, control of natural vectors, or other procedures). Disease management relies strongly on a fast and accurate identification of the causal agent. For known viruses, diagnosis consists in assigning a virus infecting a plant sample to a group of viruses sharing common characteristics, which is usually referred to as species. However, the specificity of diagnosis can also reach higher taxonomic levels, as genus or family, or lower levels, as strain or variant. Diagnostic procedures must be optimized for accuracy by detecting the maximum number of members within the group (sensitivity as the true positive rate) and distinguishing them from outgroup viruses (specificity as the true negative rate). This requires information on the genetic relationships within-group and with members of other groups. The influence of the genetic diversity of virus populations in diagnosis and disease management is well documented, but information on how to integrate the genetic diversity in the detection methods is still scarce. Here we review the techniques used for plant virus diagnosis and disease control, including characteristics such as accuracy, detection level, multiplexing, quantification, portability, and designability. The effect of genetic diversity and evolution of plant viruses in the design and performance of some detection and disease control techniques are also discussed. High-throughput or next-generation sequencing provides broad-spectrum and accurate identification of viruses enabling multiplex detection, quantification, and the discovery of new viruses. Likely, this technique will be the future standard in diagnostics as its cost will be dropping and becoming more affordable.

Extensive recombination challenges the utility of Sugarcane mosaic virus phylogeny and strain typing

Sugarcane mosaic virus (SCMV) is distributed worldwide and infects three major crops: sugarcane, maize, and sorghum. The impact of SCMV is increased by its interaction with Maize chlorotic mottle virus which causes the synergistic maize disease maize lethal necrosis. Here, we characterised maize lethal necrosis-infected maize from multiple sites in East Africa, and found that SCMV was present in all thirty samples. This distribution pattern indicates that SCMV is a major partner virus in the East African maize lethal necrosis outbreak. Consistent with previous studies, our SCMV isolates were highly variable with several statistically supported recombination hot- and cold-spots across the SCMV genome. The recombination events generate conflicting phylogenetic signals from different fragments of the SCMV genome, so it is not appropriate to group SCMV genomes by simple similarity.

Interviral Recombination between Plant, Insect, and Fungal RNA Viruses: Role of the Intracellular Ca2+/Mn2+ Pump

Recombination is one of the driving forces of viral evolution. RNA recombination events among similar RNA viruses are frequent, although RNA recombination could also take place among unrelated viruses. In this paper, we have established efficient interviral recombination systems based on yeast and plants. We show that diverse RNA viruses, including the plant viruses tomato bushy stunt virus, carnation Italian ringspot virus, and turnip crinkle virus-associated RNA; the insect plus-strand RNA [(+)RNA] viruses Flock House virus and Nodamura virus; and the double-stranded L-A virus of yeast, are involved in interviral recombination events. Most interviral recombinants are minus-strand recombinant RNAs, and the junction sites are not randomly distributed, but there are certain hot spot regions. Formation of interviral recombinants in yeast and plants is accelerated by depletion of the cellular SERCA-like Pmr1 ATPase-driven Ca2+/Mn2+ pump, regulating intracellular Ca2+ and Mn2+ influx into the Golgi apparatus from the cytosol. The interviral recombinants are generated by a template-switching mechanism during RNA replication by the viral replicase. Replication studies revealed that a group of interviral recombinants is replication competent in cell-free extracts, in yeast, and in the plant Nicotiana benthamiana We propose that there are major differences among the viral replicases to generate and maintain interviral recombinants. Altogether, the obtained data promote the model that host factors greatly contribute to the formation of recombinants among related and unrelated viruses. This is the first time that a host factor's role in affecting interviral recombination is established.IMPORTANCE Viruses with RNA genomes are abundant, and their genomic sequences show astonishing variation. Genetic recombination in RNA viruses is a major force behind their rapid evolution, enhanced pathogenesis, and adaptation to their hosts. We utilized a previously identified intracellular Ca2+/Mn2+ pump-deficient yeast to search for interviral recombinants. Noninfectious viral replication systems were used to avoid generating unwanted infectious interviral recombinants. Altogether, interviral RNA recombinants were observed between plant and insect viruses, and between a fungal double-stranded RNA (dsRNA) virus and an insect virus, in the yeast host. In addition, interviral recombinants between two plant virus replicon RNAs were identified in N. benthamiana plants, in which the intracellular Ca2+/Mn2+ pump was depleted. These findings underline the crucial role of the host in promoting RNA recombination among unrelated viruses.

Characterization of a highly divergent Sugarcane mosaic virus from Canna indica L. by deep sequencing

Background

Cannas are popular ornamental plants and widely planted for the beautiful foliage and flower. Viral disease is a major threaten to canna horticulture industry. In the city of Beijing, mosaic disease in canna was frequently observed, but the associated causal agent and its biological characterization is still unknown.

Results

After small RNA deep sequencing, 36,776 contigs were assembled and 16 of them shared high sequence identities with the different proteins of Sugarcane mosaic virus (SCMV) of the size ranging from 86 to 1911 nt. The complete genome of SCMV isolate (canna) was reconstructed by sequencing all cDNA clones obtained from RT-PCR and 5′\3′ RACE amplifications. SCMV-canna isolate showed to have a full RNA genome of 9579 nt in length and to share 78% nt and 85% aa sequence identities with SCMV isolates from other hosts. The phylogenetic tree constructed based on the full genome sequence of SCMV isolates allocated separately the canna-isolate in a distinct clade, indicating a new strain. Recombination analyses demonstrated that SCMV-canna isolate was a recombinant originating from a sugarcane-infecting isolate (major parent, acc. no. AJ310103) and a maize-infecting isolate (minor parent, acc. no. AJ297628). Pathogenicity test showed SCMV-canna could cause typical symptoms of mosaic and necrosis in some tested plants with varying levels of severity but was less virulent than the isolate SCMV-BJ. Field survey showed that the virus was widely distributed.

Conclusions

This study identified SCMV as the major agent causing the prevalent mosaic symptom in canna plants in Beijing and its genomic and biological characterizations were further explored. All these data enriched the knowledge of the viruses infecting canna and would be helpful in effective disease management in canna.

The Role of Noncoding RNAs in Double-Strand Break Repair

Genome stability is constantly threatened by DNA lesions generated by different environmental factors as well as endogenous processes. If not properly and timely repaired, damaged DNA can lead to mutations or chromosomal rearrangements, well-known reasons for genetic diseases or cancer in mammals, or growth abnormalities and/or sterility in plants. To prevent deleterious consequences of DNA damage, a sophisticated system termed DNA damage response (DDR) detects DNA lesions and initiates DNA repair processes. In addition to many well-studied canonical proteins involved in this process, noncoding RNA (ncRNA) molecules have recently been discovered as important regulators of the DDR pathway, extending the broad functional repertoire of ncRNAs to the maintenance of genome stability. These ncRNAs are mainly connected with double-strand breaks (DSBs), the most dangerous type of DNA lesions. The possibility to intentionally generate site-specific DSBs in the genome with endonucleases constitutes a powerful tool to study, in vivo, how DSBs are processed and how ncRNAs participate in this crucial event. In this review, we will summarize studies reporting the different roles of ncRNAs in DSB repair and discuss how genome editing approaches, especially CRISPR/Cas systems, can assist DNA repair studies. We will summarize knowledge concerning the functional significance of ncRNAs in DNA repair and their contribution to genome stability and integrity, with a focus on plants.

Recombination in Enteroviruses, a Multi-Step Modular Evolutionary Process

RNA recombination is a major driving force in the evolution and genetic architecture shaping of enteroviruses. In particular, intertypic recombination is implicated in the emergence of most pathogenic circulating vaccine-derived polioviruses, which have caused numerous outbreaks of paralytic poliomyelitis worldwide. Recent experimental studies that relied on recombination cellular systems mimicking natural genetic exchanges between enteroviruses provided new insights into the molecular mechanisms of enterovirus recombination and enabled to define a new model of genetic plasticity for enteroviruses. Homologous intertypic recombinant enteroviruses that were observed in nature would be the final products of a multi-step process, during which precursor nonhomologous recombinant genomes are generated through an initial inter-genomic RNA recombination event and can then evolve into a diversity of fitter homologous recombinant genomes over subsequent intra-genomic rearrangements. Moreover, these experimental studies demonstrated that the enterovirus genome could be defined as a combination of genomic modules that can be preferentially exchanged through recombination, and enabled defining the boundaries of these recombination modules. These results provided the first experimental evidence supporting the theoretical model of enterovirus modular evolution previously elaborated from phylogenetic studies of circulating enterovirus strains. This review summarizes our current knowledge regarding the mechanisms of recombination in enteroviruses and presents a new evolutionary process that may apply to other RNA viruses.

Sequence analysis of the medium and small RNAs of impatiens necrotic spot virus reveals segment reassortment but not recombination

The complete sequence of the medium (M) and small (S) RNA genome segments were determined for twelve isolates of impatiens necrotic spot virus from eight plant species. The M- and S-RNAs of these isolates shared 97-99% and 93-98% nucleotide sequence identity, respectively, with the corresponding full-length sequences available in public databases. Phylogenetic analysis based on the M- or S-RNA sequences showed incongruence in the phylogenetic position of some isolates, suggesting intraspecies segment reassortment. The lack of phylogenetic discordance in individual and concatenated sequences of individual genes encoded by M- or S-RNAs suggests that segment reassortment rather than recombination is driving evolution of these INSV isolates.

Genome characterization of citrus yellow vein-clearing virus: limited heterogeneity of viral genomes in Mandarivirus-infecting different citrus species

Citrus yellow vein-clearing virus (CYVCV) is a mandarivirus infecting citrus producing yellow vein-clearing symptoms. The leaf samples collected during surveys of different citrus-growing areas in India exhibited diverse symptoms and 40% of the plants were positive for CYVCV in RT-PCR, indicating the wide distribution of the virus in India. It was reported for first time that CYVCV infects kinnow mandarin and sweet oranges and produces chlorotic ringspots symptoms identical to Indian citrus ringspot virus (ICRSV). The complete genome sequences of CYVCV infecting four citrus cultivars have been deciphered through overlapping primers. All the four genomes comprise of 7531 nucleotides excluding the 3' poly (A) tail. The sequence identity of genomes of four CYVCV isolates in the present study ranged from 95.2 to 99.8% with genome sequences of 31 CYVCV isolates available in public domain and the mean genomic diversity was 0.017, indicating low level of heterogeneity. The phylogenetic analysis revealed that CYVCV isolates from India, Pakistan, and Turkey were clustered in the same clad apart from China isolates. The least normalized dN/dS mean value (0.092) indicated that RdRP region evolved under relatively stronger selection constraints than the other five coding regions of CYVCV. The four intragenic putative recombination events detected in RDP4 program occurred naturally in CYVCV genome, indicating the evolutionary progress of the virus. Tajima's and Fu and Li's D parameters were performed using genomic sequences in DnaSP v5 program and the retrieved negative values indicated the presence of limited genetic variability in CYVCV genomes. To the best of our knowledge, this is the first comprehensive report on molecular characterization of CYVCV from India. It will be helpful in understanding the evolutionary relationship of CYVCV and ICRSV.

Molecular characterization of a new recombinant brassica yellows virus infecting tobacco in China

Brassica yellows virus (BrYV), prevalently distributed throughout mainland China and South Korea while triggering serious diseases in cruciferous crops, is proposed to be a new species in the genus Polerovirus within the family Luteoviridae. There are three distinct genotypes (BrYV-A, BrYV-B and BrYV-C) reported in cabbage and radish. Here, we describe a new BrYV isolate infecting tobacco plants in the field, which was named BrYV-NtabQJ. The complete genome sequence of BrYV-NtabQJ is 5741 nt in length, and 89% of the sequence shares higher sequence identities (about 90%) with different BrYV isolates. However, it possesses a quite divergent region within ORF5, which is more close to Beet western yellows virus (BWYV), Beet mild yellowing virus (BMYV) and Beet chlorosis virus (BChV). A significant recombination event was then detected among BrYV-NtabQJ, BrYV-B Beijng isolate (BrYV-BBJ) and BWYV Leonurus sibiricus isolate (BWYV-LS). It is proposed that BrYV-NtabQJ might be an interspecific recombinant between BrYV-BBJ and BWYV-LS, and the recombination might result in the successful aphid transmission of BrYV from cruciferous crops to tobacco. And it also poses new challenges for BrYV diagnosis and the vegetable production.

Impact of Cultivated Hosts on the Recombination of Cucumber Mosaic Virus

Cucumber mosaic virus (CMV) is one of the most successful viruses known, infecting over 1,200 species of plants. Like other single-stranded RNA viruses, CMV is known to have a high potential for population diversity due to error-prone replication and short generation times. Recombination is also a mechanism that allows viruses to adapt to new hosts. Host genes have been identified that impact the recombination of RNA viruses by using single-cell yeast systems. To determine the impact that the natural plant host has on virus recombination, we used a high-recombination-frequency strain of CMV, LS-CMV, which belongs to subgroup II, in three different cultivated hosts: Capsicum annuum cv. Marengo (pepper), Nicotiana tabacum cv. Xanthi nc (tobacco), and Cucurbita pepo cv. Black Beauty (zucchini). The recombination frequency was calculated by using an RNA 3 reporter carrying restriction enzyme sites created by introducing silent mutations. Our results show that the recombination frequency of LS-CMV is correlated with the infected host. The recombination events in pepper were 1.8-fold higher than those in tobacco and 5-fold higher than those in zucchini. Furthermore, we observed the generation of defective RNAs in inoculated pepper plants, but not in tobacco or zucchini. These results indicate that the host is involved in both intra- and intermolecular recombination events and that hosts like pepper could foster more rapid evolution of the virus. In addition, we report for the first time the production of defective RNAs in a CMV subgroup II isolate.IMPORTANCE Recombination is an important mechanism used by viruses for their diversification and to adapt to diverse hosts. Understanding the host role in the mechanisms of evolution is important for virus disease management and controlling the emergence of new strains. This study shows the impact that cultivated hosts are playing in the evolution of CMV. Furthermore, our results and previous studies show how some specific hosts could be an ideal environment for the emergence of new viral strains.

Cowpea chlorotic mottle bromovirus replication proteins support template-selective RNA replication in Saccharomyces cerevisiae

Positive-strand RNA viruses generally assemble RNA replication complexes on rearranged host membranes. Alphaviruses, other members of the alpha-like virus superfamily, and many other positive-strand RNA viruses invaginate host membrane into vesicular RNA replication compartments, known as spherules, whose interior is connected to the cytoplasm. Brome mosaic virus (BMV) and its close relative, cowpea chlorotic mottle virus (CCMV), form spherules along the endoplasmic reticulum. BMV spherule formation and RNA replication can be fully reconstituted in S. cerevisiae, enabling many studies identifying host factors and viral interactions essential for these processes. To better define and understand the conserved, core pathways of bromovirus RNA replication, we tested the ability of CCMV to similarly support spherule formation and RNA replication in yeast. Paralleling BMV, we found that CCMV RNA replication protein 1a was the only viral factor necessary to induce spherule membrane rearrangements and to recruit the viral 2a polymerase (2a^pol) to the endoplasmic reticulum. CCMV 1a and 2a^pol also replicated CCMV and BMV genomic RNA2, demonstrating core functionality of CCMV 1a and 2a^pol in yeast. However, while BMV and CCMV 1a/2a^pol strongly replicate each others’ genomic RNA3 in plants, neither supported detectable CCMV RNA3 replication in yeast. Moreover, in contrast to plant cells, in yeast CCMV 1a/2a^pol supported only limited replication of BMV RNA3 (<5% of that by BMV 1a/2a^pol). In keeping with this, we found that in yeast CCMV 1a was significantly impaired in recruiting BMV or CCMV RNA3 to the replication complex. Overall, we show that many 1a and 2a^pol functions essential for replication complex assembly, and their ability to be reconstituted in yeast, are conserved between BMV and CCMV. However, restrictions of CCMV RNA replication in yeast reveal previously unknown 1a-linked, RNA-selective host contributions to the essential early process of recruiting viral RNA templates to the replication complex.

First natural crossover recombination between two distinct species of the family Closteroviridae leads to the emergence of a new disease

Lettuce chlorosis virus-SP (LCV-SP) (family Closteroviridae, genus Crinivirus), is a new strain of LCV which is able to infect green bean plants but not lettuce. In the present study, high-throughput and Sanger sequencing of RNA was used to obtain the LCV-SP full-length sequence. The LCV-SP genome comprises 8825 nt and 8672 nt long RNA1 and RNA2 respectively. RNA1 of LCV-SP contains four ORFs, the proteins encoded by the ORF1a and ORF1b are closely related to LCV RNA1 from California (FJ380118) whereas the 3´ end encodes proteins which share high amino acid sequence identity with RNA1 of Bean yellow disorder virus (BnYDV; EU191904). The genomic sequence of RNA2 consists of 8 ORFs, instead of 10 ORFs contained in LCV-California isolate. The distribution of vsiRNA (virus-derived small interfering RNA) along the LCV-SP genome suggested the presence of subgenomic RNAs corresponding with HSP70, P6.4 and P60. Results of the analysis using RDP4 and Simplot programs are the proof of the evidence that LCV-SP is the first recombinant of the family Closteroviridae by crossover recombination of intact ORFs, being the LCV RNA1 (FJ380118) and BnYDV RNA1 (EU191904) the origin of the new LCV strain. Genetic diversity values of virus isolates in the recombinant region obtained after sampling LCV-SP infected green bean between 2011 and 2017 might suggest that the recombinant virus event occurred in the area before this period. The presence of LCV-SP shows the role of recombination as a driving force of evolution within the genus Crinivirus, a globally distributed, emergent genus.

Identification of a putative novel genotype 3/rabbit hepatitis E virus (HEV) recombinant

Hepatitis E virus (HEV) is a viral pathogen transmitted by the fecal-oral route and is a major cause of waterborne acute hepatitis in many developing countries. In addition to infecting humans, HEV has been identified in swine, wild boars, rabbits and other mammals; with swine and wild boars being main reservoirs for zoonotic transmission of HEV. There are four major HEV genotypes known to infect humans; genotypes 1 (HEV-1) and 2 (HEV-2) are restricted to humans, and genotypes 3 (HEV-3) and 4 (HEV-4) are zoonotic. Herein, three human HEV strains originating in France were sequenced and near full-length genomes were characterized. Phylogenetic analysis showed that two strains were genotype 3 and closely grouped (a 100% bootstrap value) with subtype 3i reference strains. In percent nucleotide identities, these two strains were 94% identical to each other, 90–93% identical to subtype 3i strains, 82–86% identical to other HEV-3, and 77–79% identical to rabbit HEV strains excluding the two divergent strains KJ013414 and KJ013415 (74%); these two strains were less than 77% identical to strains of HEV genotypes 1, 2 and 4. The third strain was found distinct from any known HEV strains in the database, and located between the clusters of HEV-3 and rabbit HEV strains. This unique strain was 74–75% identical to HEV-1, 73% to HEV-2, 81–82% to HEV-3, 77–79% to rabbit HEV again excluding the two divergent strains KJ013414 and KJ013415 (74%), and 74–75% to HEV-4, suggesting a novel unclassified strain associated with HEV-3 and rabbit HEV. SimPlot and BootScan analyses revealed a putative recombination of HEV-3 and rabbit HEV sequences at four breakpoints. Phylogenetic trees of the five fragments of the genome confirmed the presence of two HEV-3 derived and three unclassified sequences. Analyses of the amino acid sequences of the three open reading frames (ORF1-3) encoded proteins of these three novel strains showed that some amino acid residues specific to rabbit HEV strains were found solely in this unclassified strain but not in the two newly identified genotype 3i strains. The results obtained by SimPlots, BootScans, phylogenetic analyses, and amino acid sequence comparisons in this study all together appear to suggest that this novel unclassified strain is likely carrying a mosaic genome derived from HEV-3 and rabbit HEV sequences, and is thus designated as a putative genotype 3/rabbit HEV recombinant.

Genetic diversity and biological characterization of sugarcane streak mosaic virus isolates from Iran

Sugarcane streak mosaic virus (SCSMV; genus Poacevirus, family Potyviridae) is a major causal agent of sugarcane mosaic disease in Asia. A survey of SCSMV was conducted in cultivated fields in Khuzestan province, southwestern Iran. Sixty-five sugarcane leaf samples showing mosaic symptoms were collected and investigated by RT-PCR. Almost one-fourth of the samples were found to be infected by SCSMV. To verify molecular variability, 12 SCSMV isolates were sequenced and analyzed by comparing partial NIb-CP gene sequences. The nucleotide identity among Iranian isolates was 83.1-99.8%, indicating high nucleotide variability, while amino acid identity was 95.2-100%, which suggesting selection for amino acid conservation. They shared nucleotide identities of 76.2-99.1% with those of other SCSMV isolates available in GenBank, the highest with isolates from Pakistan (PAK), India (IND671) and China (M117, KT257289). Further analysis was conducted based on complete CI coding region to gain more insight into the phylogenetic relationships of Iranian SCSMV compared to those from other Asian countries. Iranian isolates shared identities of 79.8-89.0% (nucleotide) and 94.8-98.6% (amino acid) with those from other geographical regions in the CI gene. The highest nucleotide identity of Iranian isolates was with isolates PAK (Pakistan), M121 (JQ975096, China) and IND671 (India), respectively. The phylogenetic trees (based on CI and NIb-CP) revealed the segregation of SCSMV isolates into two major divergent evolutionary lineages that reflect geographical origin of the isolates (with minor exception). Phylogenetic analyses grouped Iranian SCSMV isolates together with isolates from Pakistan, India and just one Chinese isolate in group II. Biological results showed that Iranian SCSMV isolates infect sugarcane, sorghum, maize and some wild grasses, causing mosaic symptoms on the leaves.

Intermolecular RNA Recombination Occurs at Different Frequencies in Alternate Forms of Brome Mosaic Virus RNA Replication Compartments

Positive-strand RNA viruses replicate their genomes in membrane-bound replication compartments. Brome mosaic virus (BMV) replicates in vesicular invaginations of the endoplasmic reticulum membrane. BMV has served as a productive model system to study processes like virus-host interactions, RNA replication and recombination. Here we present multiple lines of evidence showing that the structure of the viral RNA replication compartments plays a fundamental role and that recruitment of parental RNAs to a common replication compartment is a limiting step in intermolecular RNA recombination. We show that a previously defined requirement for an RNA recruitment element on both parental RNAs is not to function as a preferred crossover site, but in order for individual RNAs to be recruited into the replication compartments. Moreover, modulating the form of the replication compartments from spherular vesicles (spherules) to more expansive membrane layers increased intermolecular RNA recombination frequency by 200- to 1000-fold. We propose that intermolecular RNA recombination requires parental RNAs to be recruited into replication compartments as monomers, and that recruitment of multiple RNAs into a contiguous space is much more common for layers than for spherules. These results could explain differences in recombination frequencies between viruses that replicate in association with smaller spherules versus larger double-membrane vesicles and convoluted membranes.

Emergency Services of Viral RNAs: Repair and Remodeling

Reproduction of RNA viruses is typically error-prone due to the infidelity of their replicative machinery and the usual lack of proofreading mechanisms. The error rates may be close to those that kill the virus. Consequently, populations of RNA viruses are represented by heterogeneous sets of genomes with various levels of fitness. This is especially consequential when viruses encounter various bottlenecks and new infections are initiated by a single or few deviating genomes. Nevertheless, RNA viruses are able to maintain their identity by conservation of major functional elements. This conservatism stems from genetic robustness or mutational tolerance, which is largely due to the functional degeneracy of many protein and RNA elements as well as to negative selection. Another relevant mechanism is the capacity to restore fitness after genetic damages, also based on replicative infidelity. Conversely, error-prone replication is a major tool that ensures viral evolvability. The potential for changes in debilitated genomes is much higher in small populations, because in the absence of stronger competitors low-fit genomes have a choice of various trajectories to wander along fitness landscapes. Thus, low-fit populations are inherently unstable, and it may be said that to run ahead it is useful to stumble. In this report, focusing on picornaviruses and also considering data from other RNA viruses, we review the biological relevance and mechanisms of various alterations of viral RNA genomes as well as pathways and mechanisms of rehabilitation after loss of fitness. The relationships among mutational robustness, resilience, and evolvability of viral RNA genomes are discussed.

Molecular Characterization and Distribution of Two Strains of Dasheen mosaic virus on Taro in Hawaii

Dasheen mosaic virus (DsMV) is one of the major viruses affecting taro (Colocasia esculenta) production worldwide. Whole genome sequences were determined for two DsMV strains, Hawaii Strain I (KY242358) and Hawaii Strain II (KY242359), from taro in Hawaii. They represent the first full-length coding sequences of DsMV reported from the United States. Hawaii Strains I and II were 77 and 85% identical, respectively, with other completely sequenced DsMV isolates. Hawaii Strain I was most closely related to vanilla mosaic virus (VanMV) (KX505964.1), a strain of DsMV infecting vanilla in the southern Pacific Islands. Hawaii Strain II was most closely related to a taro DsMV isolate CTCRI-II-14 (KT026108.1) from India. Phylogenetic analysis of all available DsMV isolates based on amino acid sequences of their coat protein showed some correlation between host plant and genetic diversity. Analyses of DsMV genome sequences detected three recombinants from China and India among the six isolates with known complete genome sequences. The DsMV strain NC003537.1 from China is a recombinant of KJ786965.1 from India and Hawaii Strain II. Another DsMV strain KT026108.1 is a recombinant of Hawaii Strain II and NC003537.1 from China. The third DsMV strain KJ786965.1 from India is a recombinant of Hawaii Strain II and NC003537.1 from China. To our knowledge, this is the first report of recombination events in DsMV. Both Hawaii Strains I and II of DsMV were found widespread throughout the Hawaiian islands.

Non-enzymatic recombination of RNA: Ligation in loops

BACKGROUND

While the RNA world hypothesis is widely accepted, it is still far from complete: the existence of self-replicating ribozyme, consisting of potentially hundreds of nucleotides, is a core assumption for the majority of RNA world models. The appearance of such long RNA molecules under prebiotic conditions is not self-evident. Recombination seems to be a plausible way of creating RNA diversity, resulting in the appearance of functional RNAs, capable of self-replicating.

METHODS

We report here on the study of recombination process modelled with two 96 nts RNA fragments. Detection of recombination products was performed with RT-PCR followed by TA-cloning and Sanger sequencing.

RESULTS

A wide range of recombinant products was detected. We found that (i) the most efficient ligation was observed for RNA species forming bulges or internal loops, with ligation partners located within the loop; (ii) a strong preference was observed for formation of a few types of major products with a large variety of minor products; (iii) ligation could occur with participation of either 2',3'-cyclophosphate or 5'-ppp; (iv) the presence of key reaction components, i.e. 5'ppp-RNAs, enabled the formation of additional types of product; (v) molecular dynamics simulations of one of the most abundant products suggests that the ligation results in a preferable formation of 2'-5'- rather than 3'-5'-linkages.

CONCLUSIONS

The study demonstrates regularities of new RNA molecules formation with non-enzymatic recombination process.

GENERAL SIGNIFICANCE

Our findings provide new data supporting the RNA World hypothesis and show the way of new RNA sequences emergence under prebiotic conditions.

Sequence analysis of six full-length bean yellow mosaic virus genomes reveals phylogenetic diversity in India strains, suggesting subdivision of phylogenetic group-IV

We report the complete genome sequence of five bean yellow mosaic virus (BYMV) isolates (CK-GL1, CK-GL3, CK-GL4, CK-GL5 and Vfaba2) that share 74.6-98.9% (nucleotide) and 81.5-99.1% (amino acid) identity with globally available BYMV sequences. Phylogenetic analysis clustered them specifically in BYMV phylogenetic group-IV within the existing nine groups. The CK-GL1, CK-GL2, CK-GL4 and CK-GL5 isolates formed a discrete cluster within group-IV. The present study suggests subdivision of group-IV into subgroup-IVa and IVb. Moreover, infectivity assays using in vitro RNA transcripts from subgroup-IVa (CK-GL3 isolate) and IVb (CK-GL1 isolate) showed distinct biological differences between the isolates supporting subdivision.

New Insights into the Genome Organization of Yeast Killer Viruses Based on "Atypical" Killer Strains Characterized by High-Throughput Sequencing

Viral M-dsRNAs encoding yeast killer toxins share similar genomic organization, but no overall sequence identity. The dsRNA full-length sequences of several known M-viruses either have yet to be completed, or they were shorter than estimated by agarose gel electrophoresis. High-throughput sequencing was used to analyze some M-dsRNAs previously sequenced by traditional techniques, and new dsRNAs from atypical killer strains of Saccharomyces cerevisiae and Torulaspora delbrueckii. All dsRNAs expected to be present in a given yeast strain were reliably detected and sequenced, and the previously-known sequences were confirmed. The few discrepancies between viral variants were mostly located around the central poly(A) region. A continuous sequence of the ScV-M2 genome was obtained for the first time. M1 virus was found for the first time in wine yeasts, coexisting with Mbarr-1 virus in T. delbrueckii. Extra 5'- and 3'-sequences were found in all M-genomes. The presence of repeated short sequences in the non-coding 3'-region of most M-genomes indicates that they have a common phylogenetic origin. High identity between amino acid sequences of killer toxins and some unclassified proteins of yeast, bacteria, and wine grapes suggests that killer viruses recruited some sequences from the genome of these organisms, or vice versa, during evolution.

Exploiting Genetic Information to Trace Plant Virus Dispersal in Landscapes

During the past decade, knowledge of pathogen life history has greatly benefited from the advent and development of molecular epidemiology. This branch of epidemiology uses information on pathogen variation at the molecular level to gain insights into a pathogen's niche and evolution and to characterize pathogen dispersal within and between host populations. Here, we review molecular epidemiology approaches that have been developed to trace plant virus dispersal in landscapes. In particular, we highlight how virus molecular epidemiology, nourished with powerful sequencing technologies, can provide novel insights at the crossroads between the blooming fields of landscape genetics, phylogeography, and evolutionary epidemiology. We present existing approaches and their limitations and contributions to the understanding of plant virus epidemiology.