Genomic sequence analysis of four new chrysanthemum virus B isolates: evidence of RNA recombination

Evaluation of Virus-Free Chrysanthemum 'Hangju' Productivity and Response to Virus Reinfection in the Field: Molecular Insights into Virus-Host Interactions

The shoot apical meristem culture has been used widely to produce virus-free plantlets which have the advantages of strong disease resistance, high yield, and prosperous growth potential. However, this virus-free plant will be naturally reinfected in the field. The physiological and metabolic responses in the reinfected plant are still unknown. The flower of chrysanthemum 'Hangju' is a traditional medicine which is unique to China. In this study, we found that the virus-free 'Hangju' (VFH) was reinfected with chrysanthemum virus B/R in the field. However, the reinfected VFH (RVFH) exhibited an increased yield and medicinal components compared with virus-infected 'Hangju' (VIH). Comparative analysis of transcriptomes was performed to explore the molecular response mechanisms of the RVFH to CVB infection. A total of 6223 differentially expressed genes (DEGs) were identified in the RVFH vs. the VIH. KEGG enrichment and physiological analyses indicated that treatment with the virus-free technology significantly mitigated the plants' lipid and galactose metabolic stress responses in the RVFH. Furthermore, GO enrichment showed that plant viral diseases affected salicylic acid (SA)-related processes in the RVFH. Specifically, we found that phenylalanine ammonia-lyase (PAL) genes played a major role in defense-related SA biosynthesis in 'Hangju'. These findings provided new insights into the molecular mechanisms underlying plant virus-host interactions and have implications for developing strategies to improve plant resistance against viruses.

Identification and Molecular Characterization of a Novel Carlavirus Infecting Chrysanthemum morifolium in China

Chrysanthemum (Chrysanthemum morifolium) is an important ornamental and medicinal plant suffering from many viruses and viroids worldwide. In this study, a new carlavirus, tentatively named Chinese isolate of Carya illinoinensis carlavirus 1 (CiCV1-CN), was identified from chrysanthemum plants in Zhejiang Province, China. The genome sequence of CiCV1-CN was 8795 nucleotides (nt) in length, with a 68-nt 5'-untranslated region (UTR) and a 76-nt 3'-UTR, which contained six predicted open reading frames (ORFs) that encode six corresponding proteins of various sizes. Phylogenetic analyses based on full-length genome and coat protein sequences revealed that CiCV1-CN is in an evolutionary branch with chrysanthemum virus R (CVR) in the Carlavirus genus. Pairwise sequence identity analysis showed that, except for CiCV1, CiCV1-CN has the highest whole-genome sequence identity of 71.3% to CVR-X6. At the amino acid level, the highest identities of predicted proteins encoded by the ORF1, ORF2, ORF3, ORF4, ORF5, and ORF6 of CiCV1-CN were 77.1% in the CVR-X21 ORF1, 80.3% in the CVR-X13 ORF2, 74.8% in the CVR-X21 ORF3, 60.9% in the CVR-BJ ORF4, 90.2% in the CVR-X6 and CVR-TX ORF5s, and 79.4% in the CVR-X21 ORF6. Furthermore, we also found a transient expression of the cysteine-rich protein (CRP) encoded by the ORF6 of CiCV1-CN in Nicotiana benthamiana plants using a potato virus X-based vector, which can result in a downward leaf curl and hypersensitive cell death over the time course. These results demonstrated that CiCV1-CN is a pathogenic virus and C. morifolium is a natural host of CiCV1.

Virus classification based on in-depth sequence analyses and development of demarcation criteria using the Betaflexiviridae as a case study

Currently, many viruses are classified based on their genome organization and nucleotide/amino acid sequence identities of their capsid and replication-associated proteins. Although biological traits such as vector specificities and host range are also considered, this later information is scarce for the majority of recently identified viruses, characterized only from genomic sequences. Accordingly, genomic sequences and derived information are being frequently used as the major, if not only, criteria for virus classification and this calls for a full review of the process. Herein, we critically addressed current issues concerning classification of viruses in the family Betaflexiviridae in the era of high-throughput sequencing and propose an updated set of demarcation criteria based on a process involving pairwise identity analyses and phylogenetics. The proposed framework has been designed to solve the majority of current conundrums in taxonomy and to facilitate future virus classification. Finally, the analyses performed herein, alongside the proposed approaches, could be used as a blueprint for virus classification at-large.

Integrated next-generation sequencing and comparative transcriptomic analysis of leaves provides novel insights into the ethylene pathway of Chrysanthemum morifolium in response to a Chinese isolate of chrysanthemum virus B

Background

Chrysanthemum virus B (CVB), a key member of the genus Carlavirus, family Betaflexiviridae, causes severe viral diseases in chrysanthemum (Chrysanthemum morifolium) plants worldwide. However, information on the mechanisms underlying the response of chrysanthemum plants to CVB is scant.

Methods

Here, an integrated next-generation sequencing and comparative transcriptomic analysis of chrysanthemum leaves was conducted to explore the molecular response mechanisms of plants to a Chinese isolate of CVB (CVB-CN) at the molecular level.

Results

In total, 4934 significant differentially expressed genes (SDEGs) were identified to respond to CVB-CN, of which 4097 were upregulated and 837 were downregulated. Gene ontology and functional classification showed that the majority of upregulated SDEGs were categorized into gene cohorts involved in plant hormone signal transduction, phenylpropanoid and flavonoid biosynthesis, and ribosome metabolism. Enrichment analysis demonstrated that ethylene pathway-related genes were significantly upregulated following CVB-CN infection, indicating a strong promotion of ethylene biosynthesis and signaling. Furthermore, disruption of the ethylene pathway in Nicotiana benthamiana, a model plant, using virus-induced gene silencing technology rendered them more susceptible to cysteine-rich protein of CVB-CN induced hypersensitive response, suggesting a crucial role of this pathway in response to CVB-CN infection.

Conclusion

This study provides evidence that ethylene pathway has an essential role of plant in response to CVB and offers valuable insights into the defense mechanisms of chrysanthemum against Carlavirus.

Highly divergent isolates of chrysanthemum virus B and chrysanthemum virus R infecting chrysanthemum in Russia

BACKGROUND

Chrysanthemum is a popular ornamental and medicinal plant that suffers from many viruses and viroids. Among them, chrysanthemum virus B (CVB, genus Carlavirus, family Betaflexiviridae) is widespread in all chrysanthemum-growing regions. Another carlavirus, chrysanthemum virus R (CVR), has been recently discovered in China. Information about chrysanthemum viruses in Russia is very scarce. The objective of this work was to study the prevalence and genetic diversity of CVB and CVR in Russia.

METHODS

We surveyed the chrysanthemum (Chrysanthemum morifolium Ramat.) germplasm collection in the Nikita Botanical Gardens, Yalta, Russia. To detect CVB and CVR, we used RT-PCR with virus-specific primers. To reveal the complete genome sequences of CVB and CVR isolates, metatransciptomic analysis of the cultivars Ribonette, Fiji Yellow, and Golden Standard plants, naturally co-infected with CVB and CVR, was performed using Illumina high-throughput sequencing. The recombination detection tool (RDP4) was employed to search for recombination in assembled genomes.

RESULTS

A total of 90 plants of 23 local and introduced chrysanthemum cultivars were surveyed. From these, 58 and 43% plants tested positive for CVB and CVR, respectively. RNA-Seq analysis confirmed the presence of CVB and CVR, and revealed tomato aspermy virus in each of the three transcriptomes. Six near complete genomes of CVB and CVR were assembled from the RNA-Seq reads. The CVR isolate X21 from the cultivar Golden Standard was 92% identical to the Chinese isolate BJ. In contrast, genomes of the CVR isolates X6 and X13 (from the cultivars Ribonette and Fiji Yellow, respectively), were only 76% to 77% identical to the X21 and BJ, and shared 95% identity to one another and appear to represent a divergent group of the CVR. Two distantly related CVB isolates, GS1 and GS2, were found in a plant of the cultivar Golden Standard. Their genomes shared from 82% to 87% identity to each other and the CVB genome from the cultivar Fiji Yellow (isolate FY), as well as to CVB isolates from Japan and China. A recombination event of 3,720 nucleotides long was predicted in the replicase gene of the FY genome. It was supported by seven algorithms implemented in RDP4 with statistically significant P-values. The inferred major parent was the Indian isolate Uttar Pradesh (AM765837), and minor parent was unknown.

CONCLUSION

We found a wide distribution of CVB and CVR in the chrysanthemum germplasm collection of the Nikita Botanical Gardens, which is the largest in Russia. Six near complete genomes of CVR and CVB isolates from Russia were assembled and characterized for the first time. This is the first report of CVR in Russia and outside of China thus expanding the information on the geographical distribution of the virus. Highly divergent CVB and CVR isolates have been identified that contributes the better understanding the genetic diversity of these viruses.

Gynura japonica: A new host of Apple stem grooving virus and Chrysanthemum virus B in China

Gynura japonica (Thunb.) Juel [Asteraceae; syn: G. segetum (Lour.) Merr] is an important perennial medicinal herb used in China for topical treatment of trauma injuries (Lin et al. 2003). It grows naturally in the southern provinces of China and is also sometimes cultivated. During 2018-2020, wild G. japonica plants exhibiting chlorotic spots and mosaic symptoms were observed in Zhejiang province, China. To identify the possible causal agents of the disease, a single symptomatic leaf sample was collected in August 2019 and sent to Zhejiang Academy of Agricultural Sciences (Hangzhou, China) for next generation sequencing (NGS). Total RNAs extracted with TRIzol (Invitrogen, Carlsbad, USA) were subjected to high throughput sequencing on the Illumina NovaSeq 6000 platform with PE150bp and data analysis was performed by CLC Genomic Workbench 11 with default parameters (QIAGEN, Hilden, Germany). A total of 37,314,080 paired-end reads were obtained, and 11,785 contigs (961 to 10,964 bp) were generated and compared with sequences in GenBank using BLASTn or BLASTx. Of the total of 12 viral-related contigs obtained, one with a length of 6,442 nt mapped to the genomic RNA of ASGV (MN495979), seven contigs with lengths ranging from 1,034 to 2,901 nt mapped to Chrysanthemum virus B (CVB), and four mapped to broad bean wilt virus 2 (BBWV2), a virus which is known to infect G. procumbens (Kwak et al. 2017). To further confirm the presence of ASGV and CVB, primers were designed and the complete nucleotide sequences of both viruses were amplified from the original NGS sample using reverse transcription polymerase chain reaction (RT-PCR) and rapid amplification of cDNA ends (RACE) according to the manufacturer's instructions (Tiosbio, Beijing, China). BLASTn analysis revealed that the complete 6,451 nt sequence of ASGV (GenBank accession No. MW259059) shared the highest identity (81.2%) with a Chinese isolate of ASGV from citrus (MN495979). The two isolates grouped with another Chinese isolate (from pear) in phylogenetic analysis. The predicted coat protein of the virus had the highest nt identity of 93.7% (96.2% amino acid sequence identity) with that of the Chinese ASGV isolate XY from apple (KX686100). The complete genomes of two distinct molecular variants of CVB (both 8,987 nt in length) were also obtained from this sample (GenBank accession Nos. MW269552, MW269553). They shared 86.8% nt identity with each other and had 81.1% and 82.1% identity to the only known complete sequence of CVB from chrysanthemum (AB245142). Ten additional wild G. japonica plants with mosaic symptoms were collected randomly during 2019-2020 from Hangzhou (n=6) and Ningbo (n=4) in Zhejiang province and tested by RT-PCR with specific primer pairs to detect BBWV2, ASGV and CVB. RT-PCR and subsequent sequencing revealed that these three viruses were present in all the samples tested, indicating that co-infection of G. japonica by ASGV, CVB and BBWV2 is common. CVB mainly infects chrysanthemum (Singh et al. 2012), while ASGV is known as a pathogen of various fruit trees especially in the family Rosaceae, although there are recent reports that it can also infect some plants in Gramineae, Asparagaceae and Nelumbonaceae (Bhardwaj et al. 2017; Chen et al. 2019; He et al. 2019). Our results provide the first report that Gynura is a natural host of CVB and ASGV. Further surveys and biological studies are underway to evaluate the importance of Gynura as a virus reservoir for epidemics among the various hosts.

Identification and molecular characterization of a novel foveavirus from Rubus spp. in Turkey

A novel plant virus was identified by high-throughput sequencing analysis from a raspberry plant showing slight mottling symptom. The complete genome sequence of this virus is 8645 nucleotides long, including the 5' and 3' UTRs. Its genome contains five ORFs and is very close to members of the genus Foveavirus (Quinvirinae, Betaflexiviridae) in terms of genome organization, TGB presence and the sizes of the RdRp and CP proteins. The novel virus shares 33.5-51.3 % and 23.3-41.3 % nucleotide identity to other genera of the Betaflexifiviridae family based on polymerase (RdRp) and CP genes, respectively. Compared to other foveavirus species, the RdRp protein showed the highest sequence identity (45.3 %) to the RdRp of peach chlorotic mottle virus (PCMV) while the maximal sequence identity for the CP protein was 33.9 % with grapevine rupestris stem pitting-associated virus (GRSPaV). The low nucleotide and amino acid sequence identity with known foveaviruses indicated that it was a novel virus, for which the provisional name "rubus virus 1 (RuV1)" is proposed. The phylogenetic analysis supports the assignment of this virus as a new species of the genus Foveavirus. A survey of 537 Rubus spp. samples grown in six provinces of Turkey, including some symptomatic samples, showed a RuV1 prevalence of 2.2 %, confirming its presence in both raspberry and blackberry plants in a single province, although no obvious association between virus infection and specific symptoms was found.

Evolving by deleting: patterns of molecular evolution of Apple stem pitting virus isolates from Poland

In this study, 267 coat protein gene (CP) sequences from 48 Polish isolates of Apple stem pitting virus (ASPV) were determined. The genetic structure of the virus population was analysed and possible mechanisms of molecular evolution explored. We found evidence of recombination within the ASPV population and the presence of 17 ASPV molecular variants that differ in the length, number and arrangement of deletions in the CP. Population genetic analyses showed significant variation among isolates from pear and apple trees, between isolates from the same host species and, more interestingly, within isolates, supporting the existence of significant levels of variability within individual hosts, as expected by a quasispecies population structure. In addition, different tests support that selection might have been an important force driving diversification within isolates: positive selection was found acting upon certain amino acids. Phylogenetic analyses also showed that isolates did not classify according to the host species (pear or apple trees) but according to the pattern of deletions, suggesting a possible role for deletions during clade diversification.

A Virus in American Blackcurrant (Ribes americanum) with Distinct Genome Features Reshapes Classification in the Tymovirales

A novel virus with distinct genome features was discovered by high throughput sequencing in a symptomatic blackcurrant plant. The virus, tentatively named Ribes americanum virus A (RAVA), has distinct genome organization and molecular features bridging genera in the order Tymovirales. The genome consists of 7106 nucleotides excluding the poly(A) tail. Five open reading frames were identified, with the first encoding a putative viral replicase with methyl transferase (MTR), AlkB, helicase, and RNA dependent RNA polymerase (RdRp) domains. The genome organization downstream of the replicase resembles that of members of the order Tymovirales with an unconventional triple gene block (TGB) movement protein arrangement with none of the other four putative proteins exhibiting significant homology to viral proteins. Phylogenetic analysis using replicase conserved motifs loosely placed RAVA within the Betaflexiviridae. Data strongly suggest that RAVA is a novel virus that should be classified as a species in a new genus in the Betaflexiviridae or a new family within the order Tymovirales.

Complete nucleotide sequence of a new carlavirus in chrysanthemums in China

A new virus causing a serious stunt disease of chrysanthemum was identified in China by high-throughput sequencing (HTS) and named chrysanthemum virus R (CVR). The complete sequence of CVR was determined by reverse transcription polymerase chain reaction (RT-PCR) and rapid amplification of cDNA ends (RACE). The genomic RNA of CVR consists of 8,874 nucleotides (nt), excluding the poly(A) tail, contains six putative open reading frames (ORFs), and has a genomic organization typical of members of the genus Carlavirus. BLAST analysis of the full genome sequence showed low similarity (38%-56% sequence identity) to other members of the genus Carlavirus. BLAST analysis and phylogenetic analysis based on the amino acid (aa) sequences of the CVR replicase and coat protein (CP) confirmed that CVR is a distinct member of the genus Carlavirus.

New Insights into Asian Prunus Viruses in the Light of NGS-Based Full Genome Sequencing

Double stranded RNAs were purified from five Prunus sources of Asian origin and submitted to 454 pyrosequencing after a random, whole genome amplification. Four complete genomes of Asian prunus virus 1 (APV1), APV2 and APV3 were reconstructed from the sequencing reads, as well as four additional, near-complete genome sequences. Phylogenetic analyses confirmed the close relationships of these three viruses and the taxonomical position previously proposed for APV1, the only APV so far completely sequenced. The genetic distances in the respective polymerase and coat protein genes as well as their gene products suggest that APV2 should be considered as a distinct viral species in the genus Foveavirus, even if the amino acid identity levels in the polymerase are very close to the species demarcation criteria for the family Betaflexiviridae. However, the situation is more complex for APV1 and APV3, for which opposite conclusions are obtained depending on the gene (polymerase or coat protein) analyzed. Phylogenetic and recombination analyses suggest that recombination events may have been involved in the evolution of APV. Moreover, genome comparisons show that the unusually long 3’ non-coding region (3' NCR) is highly variable and a hot spot for indel polymorphisms. In particular, two APV3 variants differing only in their 3’ NCR were identified in a single Prunus source, with 3' NCRs of 214–312 nt, a size similar to that observed in other foveaviruses, but 567–850 nt smaller than in other APV3 isolates. Overall, this study provides critical genome information of these viruses, frequently associated with Prunus materials, even though their precise role as pathogens remains to be elucidated.

Characterization of New Isolates of Apricot vein clearing-associated virus and of a New Prunus-Infecting Virus: Evidence for Recombination as a Driving Force in Betaflexiviridae Evolution

Double stranded RNAs from Prunus samples gathered from various surveys were analyzed by a deep-sequencing approach. Contig annotations revealed the presence of a potential new viral species in an Azerbaijani almond tree (Prunus amygdalus) and its genome sequence was completed. Its genomic organization is similar to that of the recently described Apricot vein clearing associated virus (AVCaV) for which two new isolates were also characterized, in a similar fashion, from two Japanese plums (Prunus salicina) from a French germplasm collection. The amino acid identity values between the four proteins encoded by the genome of the new virus have identity levels with those of AVCaV which fall clearly outside the species demarcation criteria. The new virus should therefore be considered as a new species for which the name of Caucasus prunus virus (CPrV) has been proposed. Phylogenetic relationships and nucleotide comparisons suggested that together with AVCaV, CPrV could define a new genus (proposed name: Prunevirus) in the family Betaflexiviridae. A molecular test targeting both members of the new genus was developed, allowing the detection of additional AVCaV isolates, and therefore extending the known geographical distribution and the host range of AVCaV. Moreover, the phylogenetic trees reconstructed with the amino acid sequences of replicase, movement and coat proteins of representative Betaflexiviridae members suggest that Citrus leaf blotch virus (CLBV, type member of the genus Citrivirus) may have evolved from a recombination event involving a Prunevirus, further highlighting the importance of recombination as a driving force in Betaflexiviridae evolution. The sequences reported in the present manuscript have been deposited in the GenBank database under accession numbers KM507061-KM504070.

Molecular variability of cowpea mild mottle virus infecting soybean in Brazil

Molecular variability was assessed for 18 isolates of cowpea mild mottle virus (CPMMV, genus Carlavirus, family Betaflexiviridae) found infecting soybean in various Brazilian states (Bahia, Goiás, Maranhão, Mato Grosso, Minas Gerais, Pará) in 2001 and 2010. A variety of symptoms was expressed in soybean cv. CD206, ranging from mild (crinkle/blistering leaves, mosaic and vein clearing) to severe (bud blight, dwarfing, leaf and stem necrosis). Recombination analysis revealed only one CPMMV isolate to be recombinant. Pairwise comparisons and phylogenetic analysis were performed for partial genomes (ORF 2 to the 3' terminus) and for each ORF individually (ORFs 2 to 6), showing the isolates to be distinct. The topology of the phylogenetic tree could be related to symptoms, but not to the year of collection or geographical origin. Additionally, the phylogenetic analysis supported the existence of two distinct strains of the virus, designated CPMMV-BR1 and CPMMV-BR2, with molecular variations between them.

Molecular characterization of coat protein gene of Garlic common latent virus isolates from India: an evidence for distinct phylogeny and recombination

The coat protein (CP) gene of five Indian Garlic common latent virus (GarCLV) isolates was sequenced and it was 960 bp long in all the five isolates, encoding a protein of 319 amino acids. Comparative nucleotide sequence analysis revealed diversity of 4.3% among the Indian isolates and of 11.9% among all isolates worldwide. Amino acid sequence comparison showed a significant variability in the N-terminal of CP of GarCLV. Various protein analysis tools identified thirteen conserved domains and motifs including Carlavirus and Potexvirus-specific Flexi CP and Flexi N CP. Phylogenetic analysis clustered GarCLV isolates in the subgroup II with isolates from Australia, Brazil, Japan, and South Korea. Intraspecies recombination study revealed that only one of the Indian isolates was a recombinant. Interspecies recombination study suggested the absence of genetic exchange from Carlavirus species to GarCLV; conversely, GarCLV was identified as a putative donor for at least two other Carlavirus species. This is the first report of molecular variability and recombination in GarCLV isolates.

An isoform of eukaryotic initiation factor 4E from Chrysanthemum morifolium interacts with Chrysanthemum virus B coat protein

Background

Eukaryotic translation initiation factor 4E (eIF4E) plays an important role in plant virus infection as well as the regulation of gene translation.

Methodology/Principal Findings

Here, we describe the isolation of a cDNA encoding CmeIF(iso)4E (GenBank accession no. JQ904592), an isoform of eIF4E from chrysanthemum, using RACE PCR. We used the CmeIF(iso)4E cDNA for expression profiling and to analyze the interaction between CmeIF(iso)4E and the Chrysanthemum virus B coat protein (CVBCP). Multiple sequence alignment and phylogenetic tree analysis showed that the sequence similarity of CmeIF(iso)4E with other reported plant eIF(iso)4E sequences varied between 69.12% and 89.18%, indicating that CmeIF(iso)4E belongs to the eIF(iso)4E subfamily of the eIF4E family. CmeIF(iso)4E was present in all chrysanthemum organs, but was particularly abundant in the roots and flowers. Confocal microscopy showed that a transiently transfected CmeIF(iso)4E-GFP fusion protein distributed throughout the whole cell in onion epidermis cells. A yeast two hybrid assay showed CVBCP interacted with CmeIF(iso)4E but not with CmeIF4E. BiFC assay further demonstrated the interaction between CmeIF(iso)4E and CVBCP. Luminescence assay showed that CVBCP increased the RLU of Luc-CVB, suggesting CVBCP might participate in the translation of viral proteins.

Conclusions/Significance

These results inferred that CmeIF(iso)4E as the cap-binding subunit eIF(iso)4F may be involved in Chrysanthemum Virus B infection in chrysanthemum through its interaction with CVBCP in spatial.