- Open Access
Characterization of Cronartium ribicola dsRNAs reveals novel members of the family Totiviridae and viral association with fungal virulence
Virology Journal volume 16, Article number: 118 (2019)
Mycoviruses were recently discovered in the white pine blister rust (WPBR) fungus Cronartium ribicola (J.C. Fisch.). Detection and characterization of their double stranded RNA (dsRNA) would facilitate understanding of pathogen virulence and disease pathogenesis in WPBR systems.
Full-length cDNAs were cloned from the dsRNAs purified from viral-infected C. ribicola, and their cDNA sequences were determined by DNA sequencing. Evolutionary relationships of the dsRNAs with related mycoviruses were determined by phylogenetic analysis. Dynamic distributions of the viral RNAs within samples of their fungal host C. ribicola were investigated by measurement of viral genome prevalence and viral gene expression.
In this study we identified and characterized five novel dsRNAs from C. ribicola, designated as Cronartium ribicola totivirus 1–5 (CrTV1 to CrTV5). These dsRNA sequences encode capsid protein and RNA-dependent RNA polymerase with significant homologies to dsRNA viruses of the family Totiviridae. Phylogenetic analysis showed that the CrTVs were grouped into two distinct clades. CrTV2 through CrTV5 clustered within the genus Totivirus. CrTV1 along with a few un-assigned dsRNAs constituted a distinct phyletic clade that is genetically distant from presently known genera in the Totiviridae family, indicating that CrTV1 represents a novel genus in the Totiviridae family. The CrTVs were prevalent in fungal samples obtained from infected western white pine, whitebark pine, and limber pines. Viral RNAs were generally expressed at higher levels during in planta mycelium growth than in aeciospores and urediniospores. CrTV4 was significantly associated with C. ribicola virulent pathotype and specific C. ribicola host tree species, suggesting dsRNAs as potential tools for dissection of pathogenic mechanisms of C. ribicola and diagnosis of C. ribicola pathotypes.
Phylogenetic and expression analyses of viruses in the WPBR pathogen, C. ribicola, have enchanced our understanding of virus diversity in the family Totiviridae, and provided a potential strategy to utilize pathotype-associated mycoviruses to control fungal forest diseases.
The investigation of viruses that infect pathogenic fungi is important to understand interactions of plant-fungal pathogen-mycovirus. Infection of fungi by mycoviruses may result in no symptoms or in phenotypic alteration to the fungi, such as hypovirulence or hypervirulence. Elucidation of the mechanisms underlying virus-mediated hypovirulence or hypervirulence facilitates development of biological control tools for effective management of pathogenic fungi and parasites in plant and animal health [1,2,3].
The Totiviridae family is characterized by a double-stranded RNA (dsRNA) genome of 4–7 kilobases encapsidated in an icosahedral virus particle, about 40 nm in diameter. The virus genome often contains two overlapping open reading frames (ORF), encoding the viral capsid protein (CP) and RNA-dependent RNA polymerase (RdRp) . To date, the virus family includes five approved genera, namely Giardiavirus, Leishmaniavirus, Trichomonasvirus, Totivirus and Victorivirus [4, 5]. Viruses in the genera Totivirus and Victorivirus mainly infect fungi while Giardiavirus, Leishmaniavirus and Trichomonasvirus are commonly present in protozoa. Two additional genera were recently proposed, including Artivirus with arthropod and fish hosts , and Insevirus with insect hosts . Several totiviruses, mainly the members of genus Totivirus, have been reported in obligate biotrophs such as powdery mildew (ascomycetes) and rust fungi (basidiomycetes) [8, 9]. Further, some totiviruses (genus Totivirus) have also been detected in higher plants, including black currant (Ribes cultivars) .
Cronartium ribicola (J.C. Fisch) infects five-needle pine species of all ages, causing white pine blister rust (WPBR) disease in forests world-wide. The fungus was accidently introduced into western North America in 1910, killing 95% or more of white pine trees in the most seriously affected stands . C. ribicola requires two host plant species to complete its life cycle, including development of pycniospore and aeciospore stages in five-needle pines, and development of urediniospore, teliospore, and basidiospore stages on an alternate host, principally Ribes. There are several virulent C. ribicola pathotypes detected in North America, and evidence suggests cytoplasmic inheritance of virulent traits , but molecular and genetic mechanisms underlying fungal pathogenicity and virulence are poorly understood in WPBR pathosystems. The factors that contribute to C. ribicola virulence and the severity of WPBR disease on susceptible five-needle pines are also not well understood [12, 13]. RNA viruses are found in various fungal pathogen species and are considered one of the important factors influencing disease development in plants following infection by a fungal pathogen . Like other fungi, C. ribicola is susceptible to virus infection; we recently detected five novel mitoviruses in this species . An investigation into mycovirus prevalence and phenotypic impact on virulence of C. ribicola may provide an effective approach to dissect white pine-blister rust interactions.
In the present work we report detection and characterization of novel viral dsRNAs in C. ribicola. Five new dsRNA viruses were discovered, and based on their genomic sequences and organization, four were designated as new members of the genus Totivirus. The fifth dsRNA virus belongs to a novel and previously unclassified genus of the family Totiviridae. The prevalence of the mycoviruses and their RNA expression were examined in various C. ribicola samples. One of the viruses was found to be significantly associated with C. ribicola virulence.
Collection of C. ribicola samples
Fungal samples were collected at aeciospore, urediniospore, or in planta mycelium growth stages. Aeciospores were collected from intact aecial blisters occurring on western white pine (Pinus monticola, Douglas ex D. Don) stems in the Vancouver Island (British Columbia-BC, Canada) and the Cottage Grove regions (Oregon-OR, USA). Black currant plants (Ribes nigrum, cultivar Ben Nevis) were inoculated using aeciospores of BC samples in a greenhouse at the Canadian Forest Service (CFS), Victoria, BC. Urediniospores (containing some mycelia) were collected from the undersides of infected Ribes leaves when the spores began to develop on currant plants. Infected tree tissues of P. monticola and limber pine (P. flexilis) inoculated at Dorena GRC, and whitebark pine (P. albicaulis) inoculated at CFS-Victoria were collected as fungal samples during the in planta growth stage of the rust lifecycle . There was no restriction to sampling aeciospores in Canada and USA. A permit was issued by the Canadian Food Inspection Agency (CFIA) for CFS to import C. ribicola samples from OR, USA. Two fungal pathotypes (avirulent avcr2 and virulent vcr2), isolated from P. monticola trees with the susceptible genotype (cr2/cr2) and the resistant genotype (Cr2/Cr2 or Cr2/cr2), respectively; were included in the present study. A total of 33 C. ribicola samples (Additional file 5: Table S1), including six vcr2 samples, were harvested, frozen in liquid nitrogen and immediately stored at − 80 °C before RNA extraction.
dsRNA detection and purification
Total RNA was extracted from aeciospores and cankered stems of white pines as previously described , then treated with Turbo DNase I (Thermo Fisher Scientific), and nuclease S1 at 37 °C. The effect of enzymatic digestion was checked by agarose gel electrophoresis in TAE buffer (40 mM Tris base, 20 mM acetic acid, 1 mM EDTA, pH 8.0) and the gel was stained by ethidium bromide (0.5 μg/mL). Enzyme-treated dsRNA was purified again using an RNeasy Mini kit (Qiagen) before cDNA synthesis.
cDNA cloning and sanger sequencing analysis
Totivirus-like contigs were discovered in the de novo assembly of a C. ribicola RNA-seq dataset  by Blastx searching against the Viral RefSeq database (viral.1.protein.faa.gz) downloaded from the National Center for Biotechnology Information (NCBI, Bethesda, Maryland, USA) (Additional file 5: Table S2). The virus-related contigs were used to design primers for cDNA cloning of the full virus genomes (Additional file 5: Table S3).
The dsRNA molecules purified from C. ribicola urediniospores and cankered western white pine stem samples (Additional file 5: Table S1) was used as templates for cDNA synthesis using a SuperScript III First-Strand Synthesis System (Life Technologies; Burlington, ON). The 5′- and 3′-ends of the viral genomes were determined by the method of rapid amplification of complementary DNA ends (RACE) using a 5′-RACE System (Invitrogen, Grand Island, NY). Following separation of RT-PCR products on a 1% (w/v) agarose gel, DNA fragments were purified from gel using a MinElute gel extraction kit (Qiagen) and cloned into the pGEM-T Easy vector (Promega). Recombinant plasmids were purified using a QIAprep Spin Miniprep kit (Qiagen) and the sequences of the inserts were determined by Sanger sequencing using T7 and SP6 universal primers, or internal primers as needed, on an ABI 3730XL sequencer.
DNA sequences were assembled into overlapping contigs for each individual viral genome using the software Sequencer (Gene Codes, Ann Arbor, MI). Complete nucleotide sequences of the viral assemblies were deposited in GenBank under accession numbers MK967418-MK967422. Open reading frames (ORFs) were determined using the NCBI ORF Finder. Protein domains and motifs were revealed by searching against the NCBI conserved domain database (CDD)  and the Pfam protein domain database at EMBL-EBI . Multiple sequence alignment of protein sequences was conducted and identities between the sequences were calculated using Clustal-Omega at EMBL-EBI website. On the basis of the aligned sequences, phylogenetic trees were constructed using the statistical method of maximum likelihood (ML) with the MEGA version 6.0 program . A pseudoknot structure in the viral genomes was predicted using the HPknotter program .
Prevalence of viral genomes across fungal samples assessed by quantitative reverse transcription PCR (qRT-PCR)
Virus prevalence was evaluated in 33 rust samples, including three host plants infected by C. ribicola. Two micrograms of total RNA, including viral dsRNA genomes and ssRNA transcripts or replication intermediates, were used to synthesise the first strand cDNA with random primers using a Superscript VILO master mix kit (Life Technologies). Viral genome-specific primers were designed at 3′-terminal regions showing no homology to the other four viral dsRNA genomes using the software PrimerExp (ABI) (Additional file 5: Table S3), and qRT-PCR was run with three replicates per sample on an Applied Biosystems 7500 Fast Real-time PCR System (Life Technologies). Melt curve analysis showed a single sharp peak for each reaction, indicating target-specific amplification of a unique amplicon for each virus. No cross-reactions against the rest of the viruses were observed for each pair of PCR primers. C. ribicola and P. monticola genomic DNA, total RNA with no reverse transcriptase, and water were included as negative controls for each viral genome on each plate. The C. ribicola α-tubulin gene was included as an internal control for normalization of fungal RNA levels across the fungal samples.
The 2−ΔΔCT method was used to calculate relative RNA level using Expression Suite Software v1.0.3 (Life Technologies) . RNA levels were compared between fungal samples, and significance of the differences were determined using one-tailed independent t-tests (P < 0.05).
Assessment of viral RNA levels by RNA-seq analysis
Relative levels of viral transcript-related RNAs were compared across different viral genomes as well as across different stages of the fungal life-cycle using RNA-seq data from previous transcriptome analyses [13, 20]. The RNA-seq cDNA libraries were constructed from total RNA with depletion of rRNAs using a TruSeq RNA-seq Sample Preparation Kit (Illumina, San Diego, CA, US). Transcriptomic data from RNA-seq analysis were thus presumably generated mainly from transcripts or ssRNA. Nine avcr2 samples were available for viral data analysis (Additional file 5: Table S1), including three aeciospore samples (SRR1583540, SRR1583545, and SRR1583552), three urediniospore samples (SRR1583557-SRR1583559) and three western white pine cankered stem samples (SRR3273235-SRR3273237). In addition, three healthy western white pine stem samples (SRR1574690–1574692) were included as negative controls, and targeted viral transcripts were not detected in the RNA-seq data of these samples which were not infected by C. ribicola. CLC Genomics Workbench v5.5 was used for read mapping with parameter setting at a mismatch cost of 1, indel cost of 3, length fraction of 0.95, and similarity fraction of 0.95. Trimmed reads were mapped to viral genomes and only paired reads (fragments) were counted in mapping. Fragments per kilobase of exon per million reads mapped (FPKM) were used to evaluate the relative level of ssRNA(+) or transcripts. FPKM values were compared using one-tailed independent t-tests (P < 0.05) to assess differences in viral transcript levels among fungal samples or among different viruses inside the same fungal isolate.
Detection of dsRNAs in C. ribicola
Double-strand RNA (dsRNA) was detected in C. ribicola aeciospores, urediniospores, and P. monticola stems infected by C. ribicola with cankered disease symptoms (Additional file 5: Table S1). After treatment of total RNAs with DNase and nuclease S1, clear dsRNA fragments were estimated with a size of about 5 Kb by agarose gel electrophoresis (Fig. 1a). Purified dsRNA samples were then used for cDNA cloning and Sanger DNA sequencing.
Full-length genome sequences of the dsRNA viruses infecting C. ribicola
Our previously assembled C. ribicola transcriptome  was used to search for candidate dsRNA sequences by a Blast analysis. A Blastx search revealed several C. ribicola contigs with significant homology hits to the viral database. The top-hit viral species included Black raspberry virus F (BRV-F, 6.23E-05~ 3.36E-39), Saccharomyces virus L-A (SCV-LA, 5.08E-53), Xanthophyllomyces dendrorhous virus L1A (XDV-L1A, 8.21E-23~ 2.31E-120), and Persimmon latent virus (PeLV, 1.99E-12) (Additional file 5: Table S2). Because BRV-F, SCV-LA, and XDV-L1A are totiviruses and the C. ribicola dsRNA was in the size range of totivirus genomes, these C. ribicola contig sequences were used to design oligonucleotide primers for cDNA cloning of the viral genomes (Additional file 5: Table S3).
Following RT-PCR, RACE, and cDNA cloning of the dsRNA (Fig. 1a), five independent sequences were determined by Sanger DNA sequencing. The full-length dsRNAs, from large to small, were 5112, 5045, 5008, 4947 and 4706 bp in length (Additional file 5: Table S4). All sequences contain two open reading frames (ORF1 and ORF2) as displayed in Fig. 1b and Additional file 5: Table S4; their full-length cDNA sequences are shown in Additional file 5: Table S5.
A domain search against the NCBI conserved domain database (CDD) found that ORF1 of dsRNA 2, 4, and 5 encoded putative proteins with significant homologies (E-values ranging from 3.44E-03 to 3.32E-87) to the major coat protein of SCV-LA, a typical totivirus, while ORF1 of dsRNA1 encoded a putative protein with limited homology to Marek’s disease glycoprotein A (E-value at 2.85E-03) (Additional file 5: Table S4). The ORF2 of all five dsRNAs had significant homologies to the reverse transcriptase (RT, RNA-dependent DNA polymerase)-like family (1.81E-10 ~ 8.75E-60) (Additional file 5: Table S4). A domain search against the Pfam at EMBL-EBI database revealed the ORF2 of all five dsRNAs encoded putative proteins with significant homologies (E-values ranging from 3.8E-14 to 5.4E-56) to viral RNA-directed RNA-polymerases of the RdRP_4 family (Additional file 5: Table S5). Blastx search against NBCI nr database indicated that the five dsRNAs respectively had the highest identities to Pythium polare RNA virus 2 (PpRV2, LC376045), Puccinia striiformis totivirus 2 (PsV2, KY207362), Phakopsora pachyrhizi mycovirus (PaMV, KT222801), PsV3 (KY207363), and PsV5 (KY207365) (Additional file 5: Table S6 and S7). Our phylogenetic and sequence homology data indicated that the five dsRNAs originated from totiviruses that infected the fungus C. ribicola. Therefore, we tentatively named the new fungal viruses Cronartium ribicola totivirus 1 to 5 (CrTV1 to CrTV5) (Fig. 1b).
Genome organization of CrTVs
ORF1 and ORF2 overlapped in the genomes of CrTV1 to CrTV4, but were separated by a 349-bp internal sequence in the CrTV5 genome. CrTV2, CrTV3, and CrTV4 were predicted to be translated with ORF1 and ORF2 as a fusion protein through a − 1 ribosomal frameshift, supported by the presence of a GGGUUUU-like slippery site in the overlapping regions of ORF1 and ORF2 (Fig. 1b). This slippery site resembles the consensus XXXYYYZ motif presumed in the red-yeast totiviruses (XdV-LA1 and XdV-L2) and other viruses of genera Totivirus and Giardiavirus [21, 22]. A pseudoknot structure was predicted to form downstream of the putative slippery sites in the transcripts of CrTV2, CrTV3 and CrTV4 (Fig. 1c). The pseudoknot structure was proposed to pause the ribosomes, thus favoring frameshifting-translation of a fusion protein . CrTV1 was predicted to be translated as a fusion protein with a rare + 1 ribosomal frameshift, like Leishmania RNA virus 1 (LRV1) of the genus Leishmaniavirus, in the Totiviridae family . A putative pseudoknot structure was predicted in CrTV1, but no consensus slippery site was characterized for a + 1 ribosomal frameshift (Fig. 1c).
Alignments of putative protein sequences and phylogeny of CrTVs
Both assigned and unassigned representative members of the family Totiviridae were included in an alignment analysis of protein sequences using Clustal Omega. CrTV1 to CrTV5 showed the highest identities to PpRV2 (CP at 22% and RdRp at 32%), PaMV (CP at 41% and RdRP at 52%), PsV2 (CP at 35% and RdRp at 47%), PsV3 (CP at 46% and RdRp at 56%), and PsV5 (CP at 46% and RdRP at 46%), respectively (Additional file 5: Table S6 and S7). Because our samples included urediniospores collected from Ribes, CrTVs were compared to an available RdRp fragment (AJ133821) of an unidentified Totiviridae genomic RNA from black currant , and we found that CrTV5 had the highest identity at 41.03% while CrTV1 showed the lowest identity at 13.91%. Protein sequence identities among the five CrTVs themselves were much lower, ranging from 8 to 34% for CPs, and from 15 to 47% for RdRps (Additional file 5: Table S6 and S7). This suggests that the CrTVs probably evolved from different ancestors.
Based on the alignments of the putative CP and RdRp sequences, phylogenetic trees were constructed using the ML method in MEGA (Fig. 2). Five approved genera of the Totiviridae family displayed well-supported clades in both CP- and RdRp-based phylogenetic trees (Fig. 2). The viruses in the genus Totivirus were grouped into four subclades (I-A to I-D), as previously denoted [8, 9]. CrTV2, 3, and 4 were clustered in pairs with PsV2, red clover powdery mildew-associated totivirus 5 (RPaTV5), and PsV3, respectively, inside the subclade I-D of the genus Totivirus, while CrTV5 was clustered with RPaTV1a and PsV5 inside the Totivirus subclade I-A. The phylogenetic relationships indicated that CrTV2 to CrTV5 are likely new species of the Totivirus genus.
In contrast, CrTV1 was distant from the other four CrTVs and did not associate with any previously reported genera. In the RdRp-based phylogenetic tree, CrTV1 clustered together with several other unclassified members, including Hubei toti-like virus 6 (HTLV6, KX882940), PpRV2 (LC376045), Beihai barnacle virus 15 (BBV15, NC_032413), Diatom colony associated dsRNA virus (DcaRV-17A and B, BAU79517 and BAU79519), RPaTV9 (NC_028490), Ustilago maydis virus H1 (UmV-H1, NC_003823), forming a separated phylogenetic branch clade. This clade was parallel to the other three major branches displayed in the RdRp-based phylogenetic tree (Fig. 2a), including one branch that contained only the genus Totivirus, another branch that was formed by the genera Tricomonasvirus, Leishmaniavirus, and Victorivirus (including eimeriaviruses as a potential subgenus of the genus) , and a third branch that was composed of the genera Artivirus, Giardivirus and Insevirus (Fig. 2a).
In the CP phylogenetic tree (Fig. 2b), CrTV1 clustered together with PpRV2 (LC376045) and BBLV15 (NC_032413) in a separate clade. The CP phylogenetic tree displayed four major branches that were parallel with each other. The CrTV1-containing clade was closely grouped with the clade of Trichomonasvirus in one major branch. The other three major branches included a clade of Artivirus, Giardiavirus, Leishmaniavirus and Victorivirus, a clade with the genus Insevirus only, and a clade with the genus Totivirus only. In general, CP genes between two groups of toti- and toti-like viruses are not well conserved. Thus, CP-based phylogenetic trees were separately created using two datasets, one for the members of genus Totivirus (except for UmV-HI) (Additional file 1: Figure S1-a) and DcaRV17/PpRV2 and UmV-HI related viruses, respectively (Additional file 1: Figure S1-b). Similar phylogenetic relationships were observed when the two groups were separated. Phylogenetic distances of CrTV1 and several other unclassified members to the characterized genera suggests the existence of another genus in the Totiviridae family which we tentatively named as Cronartivirus in this report.
Prevalence of CrTVs affected by C. ribicola life cycle, genotypes, and host plants
The relative abundance of CrTV genomes was evaluated by qRT-PCR across 33 C. ribicola samples (Fig. 3). PCR primers were designed at the highly variable regions to distinguish different CrTV genomes (Additional file 5: Table S3). We used the C. ribicola α-tubulin gene as an internal control for normalization of total RNA extracted from different C. ribicola samples. No PCR product was detected in negative controls (including C. ribicola and P. monticola genomic DNA, total RNAs without reverse-transcriptase, and water). RNA-seq analysis also revealed the presence of CrTVs in cankered western white pines and absence in healthy pine stems. All of the available evidence indicated no integration of CrTVs into the genomes of C. ribicola and five-needle pines. qRT-PCR detected these totiviruses repeatedly in spore samples with dsRNA levels undetectable by agarose gel electrophoresis (Fig. 1a). Significant changes in genome abundance of CrTV1 through 4 during the life-cycle of C. ribicola were discovered. CrTV1, CrTV2, and CrTV3 showed higher levels during in-planta mycelium growth than during the aeciospore and urediniospore stages. CrTV4 showed an opposite pattern, with 7.5 times more abundance in fungal spores than at the in planta mycelium growth stage (15.423 ± 11.157 vs. 2.069 ± 3.661, t-test P = 0.03).
When viral abundance was compared between C. ribicola avirulent samples (avcr2) and virulent samples (vcr2), CrTV3 and CrTV5 showed no significant difference. However, CrTV1 and CrTV2 accumulated significantly more in avcr2 than in vcr2 (t-test, P < 0.05). In contrast, CrTV4 was on average 20 times more abundant in vcr2 than in avcr2 samples (t-test, P < 0.01) (Fig. 3). Pearson correlation analysis further demonstrated that abundance of CrTV4 was positively correlated to the virulence of vcr2 samples(R = 0.5127, P = 0.002, N = 33 (Additional file 2: Figure S2).
We further examined the effects of C. ribicola host plants on the prevalence of CrTVs. No significant change was detected for the relative amounts of CrTV3 and CrTV5 genomes across the three host plants (P. monticola, P. albicaulis, and P. flexilis) infected by C. ribicola. In contrast, CrTV1 showed higher genome abundance in P. flexilis than in P. monticola while CrTV2 displayed higher genome abundance in P. albicaulis than in P. flexilis (Fig. 3). Of all five CrTVs, only CrTV4 showed a host range restriction. CrTV1, 2, 3 and 5 were detected in all 33 C. ribicola samples, but the CrTV4 genome was not detectable in all nine C. ribicola-infected P. albicaulis samples and three of nine C. ribicola-infected P. flexilis samples. CrTV4 was detected in all six cankered P. monticola samples with significantly higher levels than in infected P. flexilis (5.997 ± 2.276 vs. 1.290 ± 0.535, t-test, P = 0.014). CrTV4 accumulation across all 33 C. ribicola samples are shown in Additional file 3: Figure S3. In particular, vcr2 aeciospores had a far higher CrTV4 abundance than all other sample types (one-way ANOVA, P = 0.014). Our qRT-PCR analysis results suggest that C. ribicola genotypes, developmental stages of the fungal life-cycle, and plant hosts are important factors affecting the prevalence of CrTVs.
Differential expression of CrTV RNAs assessed by RNA-seq
RNA-seq-based FPKM was used to evaluate relative viral transcript-related RNA levels across three stages of the C. ribicola lifecycle as well as across five CrTVs (Fig. 4). All of the CrTVs except CrTV5 showed relative transcript-related RNA levels significantly higher at the in planta mycelium growth stage (fold changes ranging from 4 to 36) than that at the aeciospore stage (t-test, P < 0.05). Of the five CrTVs, only CrTV3 showed a significantly higher level of RNA expression in urediniospores than in aeciospores (Fig. 4b), suggesting a possible interaction between Ribes and C. ribicola that may affect the transcription and replication activity of CrTV3.
Because data were normalized across the transcriptome, FPKM can be used to compare transcript-related RNA levels among the five viruses in the same sample (Fig. 4b). At the aeciospore stage, all five CrTVs showed similar transcript-related RNA levels. In urediniospores, CrTV1 had the highest average transcript-related RNA level, which was significantly higher than that in CrTV4 and CrTV5, while CrTV2 and CrTV3 had transcript-related RNA at intermediate levels. At the in planta growth stage in cankered P. monticola stems with both Cr2/− and cr2/cr2 genotypes, CrTV1 showed an average transcript-related RNA level 10 to 40 times higher than the other four viruses (t-test, P < 0.01), followed by CrTV2 with 4 times higher levels than CrTV5 (t-test, P < 0.05).
The accumulation profiles of CrTVs showed a significant correlation between the qRT-PCR and RNA-seq datasets for all five CrTVs across the nine avcr2 samples assessed with both methods (Pearson Correlation Coefficient R = 0.47, P = 0.001, Additional file 4: Figure S4). CrTV1 showed a Pearson correlation at the highest level (R = 0.845, P = 0.004) and CrTV4 showed a correlation at an intermediate level (R = 0.631, P = 0.068) probably due to the lower abundance of CrTV4 RNA relative to CrTV1 and the limited sample size (N = 9).
The five CrTV dsRNA sequences discovered in C. ribicola shared similar genome sizes ranging from 4.7 to 5.1 kb, but the sequence homologies between them were limited. All five CrTV genomes contained two ORFs encoding unique CP and RdRp, putatively assigned to the family Totiviridea. We identified three types of genome organizations in the CrTVs with respect to ORF translation. CrTV2, CrTV3, and CrTV4 were predicted to have two overlapping ORFs; the RdRp is likely translated as a fusion protein with the CP as a consequence of a ribosomal − 1 frameshift event. A putative canonical slippery site with GGGUUUU-like sequences was further predicted in CrTV2, CrTV3, and CrTV4, which resembles the slippery site presumed in the red-yeast totivirus XdV-LA1 and XdV-L2 . In contrast, CrTV1 RdRp was predicted to be translated as a fusion protein by a + 1 ribosomal frameshift, similar to Leishmania RNA virus 1 and a number of other members of the family Totiviridae [26,27,28,29,30]. A canonical slippery site is required for a + 1 ribosomal frameshift, but none has been confirmed in the family Totiviridae yet. CrTV5 represented a different type of genome organization, with its CP and RdRp separated by an internal non-translated sequence. This pattern was also observed in the Helminthosporium victorivirus , and the myonecrosis virus (IMNV, AY570982) that infects penaeid shrimp .
With advancements in next generation sequencing (NGS) technology and wide application of metagenomics, totiviruses have been sequenced from fungi, protozoa, arthropods, insects, plants, and fish. Small RNAs and ribosomal RNA depleted total RNA are commonly used for NGS-based viral RNA discovery . NGS does not require any prior knowledge for identification of novel viruses in field samples , especially those previously overlooked by traditional virus diagnosis tools . More recently, transcriptome sequencing was used for detection of viral agents in fungi and higher plants, such as citrus . Our previous work detected several mitoviruses in C. ribicola using RNA-seq-based NGS . In the present study, the same transcriptomic assemblies were used to search the NCBI Blast database, facilitating discovery of the CrTVs. These studies showed that sequencing of three transcriptomes in one sample allows simultaneous identification of novel RNA sequences originating from plant, fungi, and mycoviruses in fungal infected plant tissues. Co-existence of RNA viruses with different genome organizations may prevent competition among mycoviruses in the same host.
To date five genera are approved by ICTV, and two additional genera were recently proposed. The ICTV sets species demarcation threshold at 50% identity for both CP and RdRp in the family Totiviridae. All five CrTV CPs shared amino acid identities below 50% to previously reported viruses in this family. Both phylogenetic analyses of the putative RdRp and CP revealed that CrTV1 was only distantly related to the members of the above seven genera (five approved and two proposed). CrTV1 RdRp was clustered with PpRV2, DcaRV17, UmV-H1, and RPaTV9 in the RdRp phylogenetic tree, but only CrTV1 and RpRV2 were clustered together in the CP phylogenetic tree. This is probably due to the latter two sharing fungi as their hosts . CPs of DcaRV17, UmV-H1, and RPaTV9 were clustered together in one clade, much closer to the clade with Giardiavirus and Insevirus but distantly separated from the clade with CrTV1 and PpRV2. UmV-H1 and RPaTV9 likely established an additional genus in the family . DcaRV17 was detected in diatom hosts . Long-term adaptation to different host groups (fungi vs. diatom) might have caused their CPs to be phylogenetically more distant than their RdRps. CP is a structural protein with a function more related to host selection and host adaptation, while RdRp is a functional protein with replication domains more conserved across all different viruses. This may help explain differences in the phylogenetic relationships between CP- and RdRp-based analyses. In contrast, evolutionary relationships of four CrTVs with their closest homologs (CrTV2 vs. PsV2, CrTV3 vs. RPaTV5, CrTV4 vs. PsV3, and CrTV5 vs. RPaTV1a/PsV5) were very consistent between CP- and RdRp-based phylogenetic analyses. The hosts of CrTV, PsV, and RPaTV (C. ribicola, Puccinia striiformis, Erysiphe trifoliorum) are obligate biotrophs [8, 9]. The genus Totivirus infects a wide range of hosts, and totiviruses have been detected in other obligate biotrophs belonging to both the Basidiomycota and the Ascomycota, as well as higher plants, such as Ribes cultivars . Although our study detected CrTVs in C. ribicola urediniospores collected from Ribes, low identities of RdRp sequences between CrTVs and the totivirus-like sequence from Ribes black currant indicted they were different members of Totivirus. Infection of closely related hosts may lead to viral CP and RdRp having similar amino acid changes and evolutionary rates. Based on available data regarding genome organization, amino acid identity, phylogeny, and host range and host speciality, we propose CrTV1 as a representative of a new genus Cronartivirus in the family Totiviridae. As more and more virus-like sequences have been identified in recent years, classification of the family Totiviridae may need to be revised [39,40,41,42]. Discovery of new dsRNA viruses with ancient origins may provide insight for understanding the evolution and distribution of RNA viruses .
Virus infection in fungi and parasites may be symptom-less or have phenotypic effects by association with either hypovirulence or hypervirulence. To look for potential relationships between CrTV infection and C. ribicola virulence, we surveyed the prevalence of five CrTVs across 33 rust samples. We found all of the fungal samples from the infected P. albicaulis trees were devoid of CrTV4 but abundantly infected by the other four CrTVs. Absence of CrTV4 in C. ribicola-infected P. albicaulis, as well as low frequency of CrTV4 in C. ribicola-infected P. flexilis, suggests that there may be a significant disadvantage for CrTV4-infected C. ribicola to infect and propagate in the two five-needle pine species well adapted at high elevations. The pine hosts may affect the success of a mycovirus in a fungus even after infection and colonization.
Furthermore, a significant association was detected between presence of CrTV4 and C. ribicola vcr2 virulence, suggesting that there may be some epigenetic factors that caused changes in the pathogenicity of C. ribicola. Several case studies demonstrated that totiviruses play important roles in altering the fitness and virulence of the fungi and parasites. Expression of viral genome-encoded toxin in yeast and fungal species killed uninfected peers, providing the viral hosts with a survival advantage [44, 45]. LRV is associated with some benefits in survival advantages for its host parasite, such as reduced sensitivity to oxidative stress, and avoidance of immune destruction, increasing host-host parasite transmission . On the other hand, the virulence level was decreased in Beauveria bassiana and Bostrytis cinera following infection by viruses, suggesting potential application of mycoviruses as bio-control agents against agricultural pests [46, 47]. The beneficial effects of totiviral infection on virulence or survival of their fungal hosts have been reported [48,49,50], and our data suggest that mycoviruses in C. ribicola may enhance virulence on host trees. We recognize that the sampling size in the present study was limited, and a fungal inoculation test to P. albicaulis or P. flexilis as well as a large survey of field infected trees are required to monitor CrTV4 accumulation in a future study. A future comparison of virus-free and virus-infected isogenic lines would conclusively prove any link between the presence of a specific virus and fungal virulence.
It is challenging to assess the effects of the virus infections on C. ribicola in the WPBR pathosystem because invitro culture of biotrophic rust fungi is very difficult, making it unfeasible to produce virus-free and virus-infected isogenic C. ribicola lines. To date, virulence in C. ribicola is determined by observation of canker development on pine trees with major gene resistance post rust infection . Although genetic diversity of the fungal pathogen was likely limited when C. ribicola was introduced, 100 years of interaction with new hosts (especially those trees with major gene resistance) and environments has apparently begun to produce distinct fungal populations . Given the economic and ecological significance of C. ribicola infections in five-needle pines, it would be worthwhile to investigate the presence of distinct mycoviruses in a large population of field rust-infected tree samples, exploring viral association with severity of canker development and tree survival. In addition to screening of genetic resistance in different five-needle pine populations [52, 53] and characterization of C. ribicola virulent factors , confirmation of hypo- or hyper-virulent effects of the viruses on C. ribicola may provide alternative tools for WPBR management.
In the present work five novel dsRNA sequences were characterized in C. ribicola. Four (CrTV2-CrTV5) were identified as new members of the genus Totivirus, and the other (CrTV1) was proposed to be a representative of previously unclassified genus of the family Totiviridae. Both genome prevalence and transcript-related RNA expression analyses revealed dynamic changes of totiviruses over stages of the rust fungal life cycle across 33 C. ribicola samples. Significant association of CrTV4 with C. ribicola virulence suggests that mycoviruses are potential diagnostic markers for assessment of WPBR pathogenic diversity.
Availability of data and materials
Sequence data were deposited in GenBank and materials were available with reasonable request.
Cronartium ribicola totivirus
Double stranded RNA
Open reading frame
Quantitative reverse transcription-polymerase chain reaction
Rapid amplification of cDNA ends
RNA-dependent RNA polymerase
White pine blister rust
Nuss DL. Hypovirulence: mycoviruses at the fungal–plant interface. Nat Rev Microbiol. 2005;3:632–42.
Xie J, Jiang D. New insights into mycoviruses and exploration for the biological control of crop fungal diseases. Annu Rev Phytopathol. 2014;52:45–68.
Ghabrial SA, Castón JR, Jiang D, Nibert ML, Suzuki N. 50-plus years of fungal viruses. Virol. 2015;479–480:356–68.
King AM, Adams MJ, Lefkowitz EJ, editors. Virus Taxonomy: Ninth Report of the International Committee on Taxonomy of Viruses, vol. 9. San Diego: Elsevier; 2012.
ICTV—International Committee on Taxonomy of Viruses. Virus Taxonomy: 2014 Release. <http://www.ictvonline.org/virusTaxonomy.asp>. Accessed 13 Jul 2015).
Zhai Y, Attoui H, Mohd Jaafar F, Wang HQ, Cao YX, Fan SP, Sun YX, Liu LD, Mertens PP, Meng WS, Wang D, Liang G. Isolation and full-length sequence analysis of Armigeres subalbatus totivirus, the first totivirus isolate from mosquitoes representing a proposed novel genus (Artivirus) of the family Totiviridae. J Gen Virol. 2010;91:2836–45.
Zhang P, Liu W, Cao M, Massart S, Wang X. Two novel totiviruses in the white-backed planthopper, Sogatella furcifera. J Gen Virol. 2018;99:710–6.
Kondo H, Hisano S, Chiba S, Maruyama K, Andika IB, Toyoda K, Fujimori F, Suzuki N. Sequence and phylogenetic analyses of novel totivirus-like double-stranded RNAs from field-collected powdery mildew fungi. Virus Res. 2016;213:353–64.
Zheng L, Lu X, Liang X, Jiang S, Zhao J, Zhan G, Liu P, Wu J, Kang Z. Molecular characterization of novel totivirus-like double-stranded RNAs from Puccinia striiformis f. sp. tritici, the causal agent of wheat stripe rust. Front Microbiol. 2017;8:1960.
Cox S, Mayo MA, Jones AT. The occurrence of dsRNA species in apparently healthy and virus-infected ribes cultivars, and evidence that one such species originates from a member of the virus family Totiviridae. Eur J Plant Pathol. 2000;106:353–64.
Fins L, Byler JW, Ferguson D, Harvey AE, Mahalovich MF, McDonald GI, Miller D, et al. Return of the giants: restoring western white pine to the inland northwest. J Forest. 2002;100:20–6.
Kinloch BB, Sniezko RA, Dupper GE. Virulence gene distribution and dynamics of the white pine blister rust pathogen in western North America. Phytopathol. 2004;94:751–8.
Liu J-J, Sturrock RN, Sniezko RA, Williams H, Benton R, Zamany A. Transcriptome analysis of the white pine blister rust pathogen Cronartium ribicola: de novo assembly, expression profiling, and identification of candidate effectors. BMC Genomics. 2015;16:678.
Liu J-J, Chan D, Xiang Y, Williams H, Li XR, Sniezko RA, Sturrock RN. Characterization of five novel mitoviruses in the white pine blister rust fungus Cronartium ribicola. PLoS ONE. 2016;11:e0154267.
Marchler-Bauer A, Bo Y, Han L, He J, Lanczycki CJ, Lu S, Chitsaz F, Derbyshire MK, Geer RC, Gonzales NR, Gwadz M, Hurwitz DI, Lu F, Marchler GH, Song JS, Thanki N, Wang Z, Yamashita RA, Zhang D, Zheng C, Geer LY, Bryant SH. CDD/SPARCLE: functional classification of proteins via subfamily domain architectures. Nucleic Acids Res. 2017;45(D):200–3.
El-Gebali S, Mistry J, Bateman A, Eddy SR, Luciani A, Potter SC, Qureshi M, Richardson LJ, Salazar GA, Smart A, Sonnhammer ELL, Hirsh L, Paladin L, Piovesan D, Tosatto SCE, Finn RD. The Pfam protein families database in 2019. Nucleic Acids Res. 2019. https://doi.org/10.1093/nar/gky995.
Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30:2725–9.
Huang C-H, Lu CL, Chiu HT. HPknotter: a heuristic approach for detecting RNA H-type Pseudoknots - offers a variety of tools including pknotsRG, PNOTS and NUPACK. Bioinformatics. 2005;21:3501–8.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) method. Methods. 2001;25:402–8.
Liu J-J, Sniezko RA, Sturrock RN, Chen H. Western white pine SNP discovery and high-throughput genotyping for breeding and conservation applications. BMC Plant Biol. 2014;4:380.
Baeza M, Bravo N, Sanhueza M, Flores O, Villarreal P, Cifuentes V. Molecular characterization of totiviruses in Xanthophyllomyces dendrorhous. Virol J. 2012;9:140.
Wickner RB, Fujimura T, Esteban R. Viruses and prions of Saccharomyces cerevisiae. Adv Virus Res. 2013;86:1–36.
Plant EP, Jacobs KL, Harger JW, Meskauskas A, Jacobs AL, Baxter JL, Petrov AN, Dinman JD. The 9-a solution: how mRNA pseudoknots promote efficient programmed −1 ribosomal frameshifting. RNA. 2003;9:168–74.
Kim SN, Choi JH, Park MW, Jeong SJ, Han KS, Kim HJ. Identification of the +1 ribosomal frameshifting site of LRV1-4 by mutational analysis. Arch Pharm Res. 2005;28:956–62.
Wu B, Zhang X, Gong P, Li M, Ding H, Xin C, Zhao N, Li J. Eimeria tenella: a novel dsRNA virus in E. tenella and its complete genome sequence analysis. Virus Genes. 2016;52:244–52.
Maga JA, Widmer G, Lebowitz JH. Leishmania RNA virus1-mediated cap-independent translation. Mol Cell Biol. 1995;15:4884–9.
Lee SE, Suh JM, Scheffter S, Patterson JL, Chung IK. Identification of a ribosomal frame shift in Leishmania RNA virus1–4. J Biochem. 1996;120:22–5.
Ro YT, Scheffter SM, Patterson JL. Specific in vitro cleavage of a Leishmania virus capsid-RNA-dependent RNA polymerase polyprotein by a host cysteine-like protease. J Virol. 1997;71:8983–90.
Firth AE, Brierley I. Non-canonical translation in RNA viruses. J Gen Virol. 2012;93:1385–409.
Hartley M-A, Ronet C, Zangger H, Beverley SM, Fasel N. Leishmania RNA virus: when the host pays the toll. Front Cell Infect Microbiol. 2012;2:99.
Huang S, Ghabrial SA. Organization and expression of the double-stranded RNA genome of Helminthosporium victoriae 190S virus, a totivirus infecting a plant pathogenic filamentous fungus. Proc Natl Acad Sci U S A. 1996;93:12541–6.
Poulos BT, Tang KF, Pantoja CR, Bonami JR, Lightner DV. Purification and characterization of infectious myonecrosis virus of penaeid shrimp. J Gen Virol. 2006;87:987–96.
Pecman A, Kutnjak D, Gutiérrez-Aguirre I, Adams I, Fox A, Boonham N, Ravnikar M. Next generation sequencing for detection and discovery of Plant viruses and Viroids: comparison of two approaches. Front Microbiol. 2017;8:1998.
Adams I, Fox A. Diagnosis of plant viruses using next-generation sequencing and metagenomic analysis. In: Wang A, Zhou X, editors. Current Research Topics in Plant Virology. Cham: Springer; 2016. p. 323–35.
Candresse T, Filloux D, Muhire B, Julian C, Galzi S, Fort G, et al. Appearances can be deceptive: revealing a hidden viral infection with deep sequencing in a plant quarantine context. PLoS One. 2014;9:e102945.
Wang Y, Atta S, Wang X, Yang F, Zhou C, Cao M. Transcriptome sequencing reveals novel Citrus bark cracking viroid (CBCVd) variants from citrus and their molecular characterization. PLoS One. 2018;13:e0198022.
Sasai S, Tamura K, Tojo M, Herrero M-L, Hoshino T, Ohki ST, Mochizuki T. A novel non-segmented double stranded virus from an Arctic isolate of Pythium polare. Virol. 2018;522:234–43.
Urayama S-I, Takaki Y, Nunoura T. FLDS: a comprehensive dsRNA sequencing method for intracellular RNA virus surveillance. Microbes Environ. 2016;31:33–40.
Shi M, Lin X-D, Tian J-H, Chen L-J, Chen X, Li C-X, Qin X-C, Li J, Cao J-P, Eden J-S, Buchmann J, Wang W, Xu J, Holmes EC, Zhang Y-Z. Redefining the invertebrate RNA virosphere. Nature. 2016;540:539–43.
Colmant AMG, Etebari K, Webb CE, Ritchie SA, Jansen CC, van den Hurk AF, Bielefeldt-Ohmann H, Hobson-Peters J, Asgari S, Hall RA. Discovery of new orbiviruses and totivirus from Anopheles mosquitoes in eastern Australia. Arch Virol. 2017;162:3529–34.
Akinyemi IA, Wang F, Chang ZX, Wu Q. Genome characterization of the newly identified maize-associated totivirus Anhui. Arch Virol. 2018;163:2929–31.
Huang Y, Guo X, Zhang S, Zhao Q, Sun Q, Zhou H, Zhang J, Tong Y. Discovery of two novel totiviruses from Culex tritaeniorhynchus classifiable in a distinct clade with arthropod-infecting viruses within the family Totiviridae. Arch Virol. 2018;163:2899–902.
de Lima JGS, Teixeira DG, Freitas TT, Lima JPMS, Lanza DCF. Evolutionary origin of 2A-like sequences in Totiviridae genomes. Virus Res. 2019;259:1–9.
Wickner RB. Double-stranded RNA viruses of Saccharomyces cerevisiae. Microbiol Rev. 1996;60:250–65.
Schmitt MJ, Breinig F. Yeast viral killer toxins: lethality and self-protection. Nat Rev Microbiol. 2006;4:212–21.
Castro M, Kramer K, Valdivia L, Ortiz S, Castillo A. A double-stranded RNA mycovirus confers hypovirulence-associated traits to Botrytis cinerea. FEMS Microbiol Lett. 2003;228:87–91.
Dalzoto PR, Glienke-Blanco C, Kava-Cordeiro V, Ribeiro JZ, Kitajima EW, Azevedo JL. Horizontal transfer and hypovirulence associated with double-stranded RNA in Beauveria bassiana. Mycol Res. 2006;110:1475–81.
Koltin Y, Day PR. Inheritance of killer phenotypes and double-stranded RNA in Ustilago maydis. Proc Natl Acad Sci U S A. 1976;73:594–8.
Bostian KA, Elliott Q, Bussey H, Burn V, Smith A, Tipper DJ. Sequence of the preprotoxin dsRNA gene of type I killer yeast: multiple processing events produce a two-component toxin. Cell. 1984;36:741–51.
Ives A, Ronet C, Prevel F, Ruzzante G, Fuertes-Marraco S, Schutz F, Zangger H, Revaz-Breton M, Lye L-F, Hickerson SM, Beverley SM, Acha-Orbea H, Launois P, Fasel N, Masina S. Leishmania RNA virus controls the severity of Mucocutaneous Leishmaniasis. Science. 2011;331:775–8.
Richardson BA, Klopfenstein NB, Zambino PJ, McDonald GI, Geils BW, Carris LM. Influence of host resistance on the genetic structure of the white pine blister rust fungus in the western United States. Phytopathol. 2008;98:413–20.
Geils BW, Hummer KE, Hunt RS. White pines, Ribes, and blister rust: a review and synthesis. Forest Pathol. 2010;40:147–85.
Sniezko R, Smith J, Liu J-J, Hamelin R. Genetic resistance to fusiform rust in southern pines and white pine blister rust in white pines—a contrasting tale of two rust pathosystems—current status and future prospects. Forests. 2014;5:2050–83.
Ma Z, Liu JJ, Zamany A. Identification and functional characterization of an effector secreted by Cronartium ribicola. Phytopathol. 2019;109:942–51.
Authors are grateful to Dr. M. Murray at BC Ministry of Forests, Lands, Natural Resource Operations and Rural Development for providing whitebark pine seedlings; D. Chan and X.-R. Li at CFS for their technical support; and G. Roke at CFS and colleagues at Dorena Genetic Resource Center for sample collection.
This work was partly supported by the CFS-GRDI (Genomics Research and Development Initiative) fund and PFC DG’s (Director General’s) fund.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Phylogenetic trees of the family Totiviridae based on the deduced amino acid sequences of capsid proteins. Sequences were aligned using Clustal Omega, and trees were generated using the maximum likelihood method with a bootstrap test of 100 iterations in MEGA 6 software. (a) Totiviruses grouping in the genus Totivirus. (b) Toti-like sequences potentially belonging to other genera of the family Totiviridae. Sequence identifications are the same as shown in Fig. 2.
Pearson correlation analysis of fungal virulence and viral RNA levels of CrTV4 as detected by qRT-PCR. Cronartium ribicola vcr2 samples were set at a virulence level of 1 (N = 6) and avcr2 samples were set at a virulence level of 0 (N = 27).
CrTV4 RNA levels measured by qRT-PCR in 33 Cronartium ribicola samples. Bars show average RNA levels of each virus with standard error of the mean (SEM). One-way Analysis of Variance (ANOVA) showed significant difference among seven types of samples (P = 0.014).
Pearson correlation analysis of viral RNA levels as measured by RNA-seq analysis and qRT-PCR.
Cronartium ribicola samples collected in the present study. Table S2. Contigs of the Cronartium ribicola transcriptome with significant homologies to Totiviruses revealed by Blastx search. Table S3. Oligonucleotide sequences designed for PCR. Table S4. Open reading frames (ORFs) identified in the C. ribicola viral full-length genomes. Table S5. CDD search of putative C. ribicola Totivirus dsRNA genomes. Table S6. Pfam search of putative C. ribicola Totivirus dsRNA genomes. Table S7. Amino acid sequence identities among putative CrTV proteins and other members of the family Totiviridae.
About this article
Cite this article
Liu, JJ., Xiang, Y., Sniezko, R.A. et al. Characterization of Cronartium ribicola dsRNAs reveals novel members of the family Totiviridae and viral association with fungal virulence. Virol J 16, 118 (2019). https://doi.org/10.1186/s12985-019-1226-5
- Cronartium ribicola
- Double-stranded RNA
- Five-needle pines
- White pine blister rust