Next Article in Journal
Antimicrobial Efficiency of Some Essential Oils in Antibiotic-Resistant Pseudomonas aeruginosa Isolates
Next Article in Special Issue
High-Throughput Sequencing Identified Distinct Bipartite and Monopartite Begomovirus Variants Associated with DNA-Satellites from Tomato and Muskmelon Plants in Saudi Arabia
Previous Article in Journal
Analysis of Tissue-Specific Defense Responses to Sclerotinia sclerotiorum in Brassica napus
Previous Article in Special Issue
Reproducibility and Sensitivity of High-Throughput Sequencing (HTS)-Based Detection of Citrus Tristeza Virus and Three Citrus Viroids
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Discovery and Diagnosis of a New Sobemovirus Infecting Cyperus esculentus Showing Leaf Yellow Mosaic and Dwarfism Using Small-RNA High Throughput Sequencing

by
Luis Rubio
1,
Francisco J. J. Guinot-Moreno
1,2,
Carmen Sanz-López
1,3 and
Luis Galipienso
1,*
1
Instituto Valenciano de Investigaciones Agrarias (IVIA), 46113 Moncada, Valencia, Spain
2
Universidad Católica de Valencia “San Vicente Mártir”, 46002 Valencia, Valencia, Spain
3
Universitat Politècnica de València, 46022 Valencia, Valencia, Spain
*
Author to whom correspondence should be addressed.
Plants 2022, 11(15), 2002; https://doi.org/10.3390/plants11152002
Submission received: 27 June 2022 / Revised: 19 July 2022 / Accepted: 27 July 2022 / Published: 31 July 2022

Abstract

:
C. esculentus is a profitable crop in Valencia, Spain, but the emergence of a disease causing of leaf yellow mosaic, dwarfism, and a drastic decrease in tuber production has become a problem. The small-RNA high-throughput sequencing (HTS) of a diseased C. esculentus plant identified only one virus, which could be the causal agent of this disease. The amino-acid comparison with viral sequences from GenBank and phylogenetic analyses indicated that this was a new species of genus Sobemovirus, and the name Xufa yellow dwarf virus was proposed. Completion with Sanger sequencing yielded a contig of 3072 nt corresponding to about 75% of the typical genome of sobemoviruses, including ORFs 2a (polyprotein-containing protease, VPG, and other proteins), 2b (RNA-dependent RNA polymerase), and 3 (coat protein). The nucleotide sequence was used to develop fast and accurate methods for the detection and quantification of xufa yellow dwarf virus (XYDV) based on reverse transcription (RT) and DNA amplification. XYDV was detected in leaves and tubers and showed a high incidence in the field in both symptomatic (almost 100%) and asymptomatic (70%) plants, but its accumulation was much higher in symptomatic plants. The relevance of these results for disease control was discussed.

1. Introduction

Cyperus esculentus, known as tiger nut, chufa or xufa, yellow nutsedge, iron water chestnut, underground chestnut, underground walnut, earth almond, ginseng fruit, ginseng bean, etc., is considered an invasive weed in most tropical, subtropical, and Mediterranean regions. However, C. esculentus is cultivated in a few places, since it produces edible sweet nut-like tubers, which are used mainly as a side dish in some Western African countries and as fodder for animals in some American countries [1]. In Europe, C. esculentus is cultivated mainly in a small region of about 400 ha in Valencia, Spain, and is used to elaborate a sweet milky beverage named horchata that is served cold in summer as a natural refreshment. The elaboration and commercialization of horchata represent a very profitable industry in Valencia, producing about 50 million liters per year that are prepared from about 6000 tons of dried tubers, where half of them is produced in Valencia and the other half is imported from Western Africa.
Since 2015, farmers have found diseased plants showing leaves with a yellowed mosaic pattern of parallel discontinuous streaks, dwarfism in the aerial parts and roots, and a severe reduction in tuber production (Figure 1).
The first step to eradicate or mitigate the disease is the identification of the causal agent. The symptoms of leaf yellowing and plant stunting present in these C. esculentus plants are typical of viral infections. The high-throughput sequencing (HTS) of viral small RNAs (produced by the plant defense mechanism based on RNA silencing) is a powerful tool generating in a short time thousands of unbiased sequences corresponding to all viruses present in the sample [2,3]. HTS does not require previous knowledge of the targeted organisms, and it is especially suited to discovering new viruses. The obtained nucleotide sequences can be ascribed to different taxonomic levels, such as virus species, genera, o families, according to the nucleotide or amino-acid identity with known viral sequences from the GenBank database or the presence of sequence motifs [4].
However, HTS is too expensive for routine diagnosis, and it is necessary to develop rapid and accurate detection techniques for each virus to study virus incidence and epidemiology and apply disease-control measures such as certification and eradication. The techniques based on reverse transcription (RT) and DNA amplification are the fastest to be developed and can be designed from the sequences obtained with HTS. Polymerase chain reaction (PCR) is the most widely used, enabling fast, highly specific, and sensitive diagnosis, but real-time quantitative PCR (qPCR) with SYBR Green or TaqMan probes is more sensitive and can quantify the virus titer, which is very useful to evaluate disease-control measures such as the plant breeding of resistant cultivars and cross protection [4]. Loop-mediated isothermal amplification (LAMP) is a specific and sensitive technique with the advantage of not requiring expensive thermal-cycling instruments, being less sensitive to inhibitors, unlike PCR and qPCR, and having a great potential to be used directly on field [5].
In this work, a C. esculentus plant with leaf mosaic yellowing and dwarfism was analyzed using small-RNA HTS, showing that only one virus was present. The amino-acid sequence comparison and phylogenetic analyses showed that it was a new species of genus Sobemovirus that was proposed to be named Xufa yellow dwarf virus. The obtained nucleotide sequence was used to develop RT-PCR, RT-PCR with SYBR and a TaqMan probe, and RT-LAMP for the fast and specific detection and quantification of xufa yellow dwarf virus (XYDV). These detection methods were applied to detect XYDV in leaves and tubers, evaluate its incidence on field and estimate XYDV accumulation in plants.

2. Results

2.1. Identification of a New Sobemovirus

The HTS of RNA extracts from a C. esculentus diseased plant yielded 622 contigs, but only 15 corresponded to viruses or viroids, and the rest were from the host. The Blastn analyses showed a 59 nt contig with a 92% identity with Hop stunt viroid (HSVd). However, the real-time RT-qPCR analysis with specific HSVd primers [6] was negative for this plant and other plants with the same symptoms. In addition, the 59 nt contig was formed by small RNAs with low reads and no similarity with HSVd, indicating that HSVd was not present in these plants and that this contig was an assembly artifact. The Blastx analyses showed 14 contigs with amino-acid similarity with viruses of genus Sobemovirus, such as rottboelia yellow mottle virus (RYMoV) and ryegrass mottle virus (RGMoV). To fill the gaps between contigs and confirm the sequences obtained with HTS, primers were designed (Table 1) and used for RT-PCR and Sanger sequencing. The sequence was extended by 105 and 91 nt towards the 5′ and 3′ ends, respectively, with a 5′/3′ RACE system. When all sequences were assembled, a single contig of 3072 nt was obtained; this corresponded to the genome of a virus with amino-acid similarity with members of genus Sobemovirus (Figure 2a). This virus was named Xufa yellow dwarf virus based on C. esculentus’ name in the local language, “xufa” (Valencian/Catalan), and the symptoms. No other viruses or viroids were identified.
About 75% of the genome of xufa yellow dwarf virus (XYDV) was sequenced and showed the typical genomic organization of the genus Sobemovirus [7]. The XYDV sequence included three ORFs (Figure 2a). ORF 2a (about 90% sequenced) encodes polyprotein P2a containing a transmembrane-anchoring domain, a serine protease, a VPg, and other products called P10 and P8. ORF 2b (about 94% sequenced) encodes an RNA polymerase RNA dependent (RdRp) expressed as a fusion protein via a −1 programmed ribosomal frameshift (−1PRF) with respect to ORF 2a. The XYDV sequence contains the −1 PRF signal consisting of a slippery sequence, 5′-UUUAAAC-3′, followed by a stem-loop structure. ORF 3 (about 75% sequenced) encodes the coat protein (CP) expressed by a subgenomic RNA. This nucleotide sequence was deposited in GenBank under accession number ON828429.
The comparison of the XYDV sequence with those of the members of genus Sobemovirus showed amino-acid identities ranging from 10.8 to 63.2% (Table 2). The most similar viruses to XYDV were RYMoV, RGMoV, and artemisia virus A (ArtVA). The three ORFs showed different levels of variability, with ORF 2b (RdRp) being the most conserved and ORF 3 (CP) the most variable. The phylogenetic relationships were similar for the three ORFs, so only ORF 3 (CP) is shown in Figure 2b. XYDV formed a statistically supported clade with RGMoV, RyMoV, and ArtVA, confirming that XYDV is a member of genus Sobemovirus, and the low amino-acid identity between these viruses indicates that XYDV is a different species of this genus.

2.2. Development of Methods for Fast Detection and Quantification of XYDV

Based on the obtained nucleotide sequence, primers and a TaqMan probe (Table 1). were designed for techniques based on reverse transcription (RT) and DNA amplification: polymerase chain reaction (PCR), real-time quantitative PCR (qPCR) with SYBR Green, qPCR with a TaqMan probe, and loop-mediated isothermal amplification (LAMP). These techniques detected XYDV in the leaves and tubers of C. esculentus plants and were negative for the negative controls, which were water, RNA and crude extracts from healthy wild C. esculentus plants and Nicotiana benthamiana plants and serial dilutions from these extracts (data not shown).
The sensitivity of these techniques was compared using serial ten-fold dilutions of 10 ng/µL total-RNA extracts from an XYDV-infected plant (Figure 3a). RT-PCR and RT-LAMP detected XYDV until the fourth dilution, corresponding to 1 pg, whereas RT-qPCR with SYBR or the TaqMan probe detected it in the fifth dilution, corresponding to 100 fg. The Ct values were similar for both RT-qPCR techniques, being slightly lower for RT-qPCR with the TaqMan probe (Figure 3a). No amplifications were obtained from the negative controls.
To find out if the RNA-extraction step could be avoided, these amplification methods were assayed with concentrated and ten-fold serial dilutions of crude plant extracts (only ground with buffer) from two XYDV isolates: CM1, obtained from a plant with severe symptoms, and CMN1, from an asymptomatic plant (hypothesized having lower viral titer; see below). RT-PCR detected XYDV in the CM1 crude-extract concentrate and the two subsequent 10-fold dilutions but not in any of the CMN1 crude extracts, despite being positive for the RNA extracts (Figure 3b). The intensity of the RT-PCR bands for CM1 was similar for the crude extract and the three serial dilutions, suggesting that the viral RNA titer at the lowest dilution was enough to attain maximum amplification. In addition, the concentrated crude extracts did not contain enough inhibitors to interfere with RT-PCR.
RT-LAMP detected XYDV in the CM1 concentrate and the two dilutions, as well as the CMN1 concentrate. Finally, RT-qPCR with SYBR Green and the TaqMan probe detected XYDV in CM1 crude extracts until dilution 1/1000 and CMN1 until dilution 1/10 (Figure 3b). Curiously, the Cts of the CM1 concentrate were higher than those of dilutions 1/10 and 1/100, suggesting the presence of inhibitors.

2.3. Incidence of XYDV in the Field and Viral Accumulation

Leaf samples were collected from 40 symptomatic and 20 asymptomatic C. esculentus plants from eight plots of the five most productive municipalities of Horta Nord comarca (county), Valencia province, Spain, in 2019. Crude extracts were analyzed with RT-LAMP, and RNA extracts were analyzed with RT-qPCR using the TaqMan probe and with RT-LAMP (Figure 4a). The RT-qPCR analysis of RNA extracts detected XYDV in 39 symptomatic plants (97.5 %), whereas RT-LAMP was a little less sensitive and detected XYDV in 37 or 38 (92.5 or 95.5%) symptomatic-plant RNA or crude extracts, respectively. Curiously, XYDV was also detected in a high proportion of asymptomatic plants (70.0% with RT-qPCR). The difference in XYDV incidence between symptomatic and asymptomatic plants was not statistically significant. XYDV was also detected in the tubers from five symptomatic plants with RT-LAMP and RT-qPCR.
The RT-qPCR analyses showed a correlation between XYDV accumulation in leaves and symptom manifestation, since the XYDV titer was about 100 times higher in symptomatic (4.8 ± 0.2) than in asymptomatic plants (2.2 ± 0.4), and this difference was statistically significant (Figure 4b). The accumulation of XYDV in tubers from the five symptomatic plants was high (4.2 ± 0.3), although it was significantly higher in the leaves of these plants (5.4 ± 0.1).

3. Discussion

The cultivation of C. esculentus and the elaboration of horchata have become substantial revenue sources for the county of Horta Nord in Valencia, Spain. However, some diseases have risen in the last decades, such as tuber rot [8] and leaf necrosis [9] produced by fungi and the black spot of unknown etiology [10]. These diseases have caused severe economic losses and threatened tuber production. This work described the first viral disease detected in C. esculentus in Spain. Small-RNA HTS enabled the identification of a new virus of genus Sobemovirus, and the name Xufa yellow dwarf virus was proposed. To our knowledge, only four viruses infecting C. esculentus have been reported: turnip mosaic virus (TuMV) in Zimbabwe [11]; impatiens necrotic spot virus (INSV) in Georgia, USA [12]; brome streak mosaic virus (BrSMV) in Hungary [13]; and rice yellow mottle virus (RYMV) in Nigeria [14]. The low number of viruses described worldwide is probably due to having been poorly investigated, since C. esculentus is not widely used in agriculture. The emergence of viral diseases in different crops in Valencia and other Mediterranean areas is frequent, mainly caused by intensive agriculture, climate change, and the global movement of plant material [15].
The detection of xufa yellow dwarf virus (XYDV) in tubers could be important, since tubers are used as seeds and could be a means of transmission, as reported for some sobemoviruses [16]. Thus, XYDV could have been introduced in Valencia by the importation of C. esculentus tubers from Western Africa. Tubers are imported without a phytosanitary passport, as they can be only used for human or animal consumption. However, some farmers cultivate African tubers despite this being forbidden, with a high risk of dispersing pathogens and pests. Curiously, another sobemovirus, RYMV, has been also detected in C. esculentus in Western Africa [14]. Knowing the transmission means of viruses is key to applying prophylactic measures for disease control, e.g., the certification and control of vectors. Transmission assays and/or epidemiological studies are necessary to know if XYDV is transmitted by tubers, mechanical wounding, and/or insects as other sobemoviruses [16]. In case XYDV is tuber-transmitted, in addition to the sanitary certification of imported tubers, other control measures should be tested, such as the selection of XYDV-free tubers, chemical and/or thermal disinfection, or in vitro culture [4,10]. The molecular techniques developed here are very sensitive and can be used to test the effectiveness of removing XYDV from tubers. RT-LAMP could be used by farmers in situ; since it only requires a thermal bath, it can be directly used with plant crude extracts, and the amplification can also be visualized by means of naked-eye inspection using turbidity-changing ion indicators or color-changing fluorescent molecules [5].
These techniques were applied in the survey of a field, and XYDV was detected in almost all C. esculentus plants with leaf yellow mosaic and dwarfism and in 70% of the asymptomatic plants. The presence of XYDV in asymptomatic plants and the high overall incidence of this virus do not make the use of eradication strategies such as roguing feasible [4]. In this scenario, the most plausible disease-control strategy would be using virus-free tubers as seeds. Farmers could test the presence of XYDV in these tubers with RT-LAMP in each sowing season and evaluate the incidence in the field. However, as mentioned above, it is necessary to know if XYDV can be transmitted by other means, which would require additional measures. The RT-qPCR analyses showed that the XYDV titer was significantly much higher in symptomatic than in asymptomatic plants. This suggests that XYDV could interfere with root development and plant growth in seedlings when XYCV concentration surpasses a certain threshold. A positive correlation between symptoms and viral accumulation has been also reported for several plant viruses [17,18,19], although this does not apply to tolerant hosts [20]. The correlation among XYDV accumulation and symptom manifestation, as well as the fact that XYDV was the only viral sequence found using HTS strongly supports that XYDV is the causal agent of the disease of leaf yellow mosaic and dwarfism in C. esculentus.

4. Materials and Methods

4.1. Plant Material

Leaf samples and tubers were collected from 40 symptomatic and 20 asymptomatic C. esculentus plants from eight plots in the municipalities of Alboraya, Burjassot, Godella, Massalfasar, Moncada, and Vinalesa, which belong to the comarca (county) of Horta Nord, Valencia province, Spain.

4.2. Extract Preparation

About 200 mg of leaf tissue or tuber was ground with liquid nitrogen in microtubes containing glass beads in a TissueLyser power homogenizer (Qiagen, Hilden, Germany). Total-RNA extracts were obtained with Spectrum Plant Total RNA Kit (Sigma-Aldrich, San Luis, MO, USA). Crude extracts were prepared by means of the resuspension of ground plant tissue in 500 µL of STE buffer (100 mM NaCl, 10 mM Tris-HCl,1 mM EDTA (pH 7.5)) and homogenization with a vortex mixer for 2 min. RNA and crude extracts were kept at −80 °C until use.

4.3. Sequencing

For the HTS of small RNAs, total-RNA concentration and purity were determined with the Qubit® RNA assay kit in a Qubit® 3.0 fluorometer (Thermo Fisher Scientific, Waltham, MA, USA) and a NanoPhotometer® spectrophotometer (Implen, Westlake Village, CA, USA), respectively. RNA integrity was determined in an Agilent Bioanalyzer 2100 system with the RNA Nano 6000 assay kit (Agilent Technologies, Santa Clara, CA, USA). cDNA was obtained from 1 µg of total-RNA with NEBNext® Multiplex Small RNA Library Prep Set for Illumina® (Sigma-Aldrich, San Luis, MO, USA) and sequenced in the Illumina NextSeq550 platform (Illumina, San Diego, CA, USA). Sequencing adapters were trimmed from the reads, and low-quality reads were filtered in November 2019 using SeqTrimNext V2.0.67 software (https://github.com/dariogf/SeqtrimNext, accessed on 30 November 2019) by applying the standard parameters for Illumina short reads [21]. VirusDetect V 1.7 [22] and the virus nucleotide database (http://bioinfo.bti.cornell.edu/ftp/program/VirusDetect/virus_database/v239, accessed on 30 November 2019) were used for virus identification. Reads were aligned with the virus nucleotide database with the Burrows–Wheler Aligner program [23] and assembled with Velvet V.1.2.09 software [24]. Reads not aligning with viral sequences were recovered, assembled, and compared to the GenBank database with Blastn and Blastx [25] to find similar viral nucleotide and/or protein sequences. Sanger sequencing was performed on RT-PCR products purified using QIAquick PCR Purification Kit (Qiagen) using a Big Dye Terminator V. 3.0 Cycle Sequencing kit (Thermo Fisher Scientific, Waltham, MA, USA) in an ABI 3130 XL capillary sequencer (Applied Biosystems, Thermo Fisher Scientific, Waltham, MA, USA). Afterwards, 5′/3′ RACE Kit, 2nd Generation (Roche, Basilea, Switzerland), was used to determine the 5′ and 3′ terminal sequences.

4.4. Amino-Acid Identity and Phylogenetic Analyses

The nucleotide sequences of ORFs 2a, 2b, and 3 of XYDV were translated to amino acids (https://web.expasy.org/translate, accessed on 30 November 2019), and the equivalent sequences of members of genus Sobemovirus were retrieved from GenBank (Table 2). These sequences were aligned with the CLUSTALW algorithm [26] implemented in MEGA X software [27]. Amino-acid identities were estimated with formula 100 × (1 − pdistance) in MEGA X. The phylogenetic relationships of the three ORFs were inferred with the Maximum Likelihood method [28] using the LG + F + G model, which is based on the LG matrix of amino-acid substitution [29] with empirical amino-acid frequencies and five discrete gamma rates. The statistical significance of the internal nodes was estimated with a bootstrap analysis with 100 replicates [30]. All these analyses were carried out in MEGA X.

4.5. Primer and Probe Design

The nucleotide sequence obtained for XYDV was used to design primers for RT-PCR and RT-qPCR and a TaqMan probe with Primer3 [31] and Primer Express (Thermo Fisher Scientific, Waltham, MA, USA). Primers for RT-LAMP were designed with LAMP Designer 1.16 (Premier Biosoft, Palo Alto, CA, USA). The list of the designed primers is shown in Table 1.

4.6. RT and DNA Amplification

RT was performed by denaturing total RNA or plant crude extracts mixed with 5 µM of random hexamers by heating at 65 °C for 5 min and incubating on ice for 1 min; then, they were added to a 20 µL reaction mixture containing 5 mM DTT, 1 mM of each dNTP, first-strand buffer, 20 units of RNaseOUT™ RNase Inhibitor (Invitrogen), and 100 units of SuperScript IV Reverse Transcriptase (Invitrogen, Thermo Fisher Scientific, Waltham, MA, USA). The reaction mixture was incubated at 23 °C for 10 min, at 50 °C for 10 min, and at 80 °C for 10 min to inactivate the reaction.
PCR was performed in a 20 µL reaction mixture containing PCR buffer, 1.5 mM MgCl2, 1 mM of each dNTP, 1 µM of primers X1F and X1R (Table 1), and 0.5 U of Taq DNA polymerase (Invitrogen, Thermo Fisher Scientific, Waltham, MA, USA). Eppendorf® Mastercycler EP Gradient was used under the following thermocycling conditions: denaturation at 94 °C for 2 min; 40 cycles at 94 °C for 30 s, at 52 °C for 1 min, and at 72 °C for 45 s; and an extension step at 72 °C for 5 min. PCR products were analyzed by means of electrophoresis in 2% agarose gels and visualized with Gel Red (Biotium. Fremont, CA, USA) under UV light.
For qPCR, total-RNA concentration was measured with a NanoDrop™ 1000 spectrophotometer and adjusted to 10 ng/µL to normalize the different extractions. Real-time qPCR with SYBR Green was carried out in LightCycler®480 (Roche Molecular Systems Inc.) in a 20 µL mixture containing TB Green Premix Ex Taq II (Tli RNase H Plus) (Takara), 0.2 µM of primers qX1F and qX1R (Table 1), and 1/10 of cDNA from the RT reaction. The thermocycling conditions included denaturation at 95 °C for 30 s and 35 cycles of 95 °C for 5 s and at 60 °C for 30 s. The melting-curve analysis of the PCR products showed only one peak, indicating the absence of nonspecific products and/or primer-dimers. One-step real-time RT-qPCR with a TaqMan probe was performed in LightCycler®480 (Roche Molecular Systems, Pleasanton, CA, USA) in a 20 µL mixture prepared with One Step PrimeScript RT-PCR Kit (TaKaRa, Shiga, Japan ) and containing 0.2 µM of primers qX2F and X2R and 0.4 µM TaqMan probe (Table 1). The thermocycling conditions consisted of RT at 42 °C for 25 min, followed by denaturation at 95 °C for 15 s and 35 cycles of 95 °C for 5 s and at 60 °C for 20 s. For the relative quantification of XYDV, a standard curve was obtained with the RT-qPCR of the sample with the lowest Ct and 10-fold serial dilutions assigning arbitrary values starting from seven. The standard curve showed a high determination coefficient (R2 = 0.9994), indicating a strong linear relationship, and the linear regression formula was y = −3.642x + 27.02.
One-step RT-LAMP was performed in a 25 µL mixture containing a sample (RNA of crude extracts) previously denatured at 95 °C for 5 min, Isotermal Amplification Buffer, 1.4 mM dNTPS, 5 U WarmStar RT (New England Biolabs, Ipswich, MA, USA), 8 U Bst DNA polymerase (New England Biolabs, Ipswich, MA, USA), and the three XYDV-specific primer pairs: 0.2 µ XF3 and XB3, 1.6 µM XFIP and XBIP, and 0.4 µM XLoopF and XLoopR (Table 1). The mixture was incubated at 68 °C for 1 h in a water bath and heated to 80 °C for 10 min to stop the reaction. The amplification products were analyzed by means of electrophoresis in 2% agarose gels.

4.7. Statistical Analysis

Symptomatic and asymptomatic plants were compared for XYDV incidence and accumulation with chi-squared and unpaired two-tailed Student’s t-tests, respectively. The analyses were performed with Statgraphics plus Software Version 5.1 using an alpha level of 0.05.

Author Contributions

Experimental design, L.R. and L.G.; Collection of samples, L.R., L.G. and F.J.J.G.-M.; Nucleotide sequence analyses and primer design, L.R. and L.G.; Experimental procedures—sample processing and analysis using RT-PCR, RT-qPCR, and RT-LAMP, F.J.J.G.-M. and C.S.-L.; Writing, L.R.; Funding acquisition, L.G.; Supervision, L.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research study was funded by Instituto Valenciano de Investigaciones Agrarias (IVIA) grants 51903, 51912, and 52202C.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The nucleotide sequence of the virus with the proposed name xufa yellow dwarf virus generated in this work was deposited in GenBank under accession number ON828429. All the sequencing-output datasets generated in the study are freely available upon request to the corresponding author.

Acknowledgments

We thank Consejo Regulador D.O. Chufa de Valencia for assistance in the field surveys.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Sánchez-Zapata, E.; Fernández-López, J.; Pérez-Álvarez, J.A. Tiger nut (Cyperus esculentus) commercialization: Health aspects, composition, properties, and food applications. Compr. Rev. Food Sci. Food Saf. 2012, 11, 366–377. [Google Scholar] [CrossRef]
  2. Adams, I.P.; Glover, R.H.; Monger, W.A.; Mumford, R.; Jackeviciene, E.; Navalinskiene, M.; Samuitiene, M.; Boonham, N. Next-generation sequencing and metagenomic analysis: A universal diagnostic tool in plant virology. Mol. Plant Pathol. 2009, 10, 537–545. [Google Scholar] [CrossRef]
  3. Kreuze, J.F.; Pérez, A.; Untiveros, M.; Quispe, D.; Fuentes, S.; Barker, I.; Simon, R. Complete viral genome sequence and discovery of novel viruses by deep sequencing of small RNAs: A generic method for diagnosis, discovery and sequencing of viruses. Virology 2009, 388, 1–7. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Rubio, L.; Galipienso, L.; Ferriol, I. Detection of plant viruses and disease management: Relevance of genetic diversity and evolution. Front. Plant Sci. 2020, 11, 1092. [Google Scholar] [CrossRef]
  5. Panno, S.; Matić, S.; Tiberini, A.; Caruso, A.G.; Bella, P.; Torta, L.; Stassi, R.; Davino, S. Loop mediated isothermal amplification: Principles and applications in plant virology. Plants 2020, 9, 461. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Papayiannis, L.C. Diagnostic real-time RT-PCR for the simultaneous detection of Citrus exocortis viroid and Hop stunt viroid. J. Virol. Methods 2014, 196, 93–99. [Google Scholar] [CrossRef] [PubMed]
  7. Sõmera, M.; Fargette, D.; Hébrard, E.; Sarmiento, C. ICTV Virus Taxonomy Profile: Solemoviridae 2021. J. Gen. Virol. 2021, 102, 001707. [Google Scholar] [CrossRef] [PubMed]
  8. García-Jiménez, J.; Busto, J.; Vicent, A.; Sales, R.; Armengol, J. A tuber rot of Cyperus esculentus caused by Rosellinia necatrix. Plant Dis. 1998, 82, 1281. [Google Scholar] [CrossRef] [PubMed]
  9. Montano-Mata, N.J. Etiologia, Epidemiologia y Control de la necrosis foliar de la chufa (Ciperus esculentus L.). Ph.D. Thesis, Universitat Politècnica de València, Valencia, Spain, 2008. [Google Scholar]
  10. Alvares, D.; Armero, C.; Forte, A.; Serra, J.; Galipienso, L.; Rubio, L. Incidence and control of black spot syndrome of tiger nut. Ann. Appl. Biol. 2017, 171, 417–423. [Google Scholar] [CrossRef]
  11. Chivasa, S.; Ekpo, E.J.A.; Hicks, R.G.T. New hosts of Turnip mosaic virus in Zimbabwe. Plant Pathol. 2002, 51, 386. [Google Scholar] [CrossRef]
  12. Martínez-Ochoa, N.; Mullis, S.W.; Csinos, A.S. First report of yellow nutsedge (Cyperus esculentus) and purple nutsedge (C. rotundus) in Georgia naturally infected with Impatiens necrotic spot virus. Plant Dis. 2004, 88, 771. [Google Scholar] [CrossRef] [PubMed]
  13. Takács, A.P.; Kazinczi, G.; Horvath, J.; Gaborjanyi, R. Cyperus esculentus L. a new host of Brome streak mosaic virus (BrSMV). J. Plant Dis. Prot. 2008, 21, 527–528. [Google Scholar]
  14. Salaudeen, M.T.; Banwo, O.O.; Kashina, B.D.; Alegbejo, M.D. Possible wild hosts of rice yellow mottle Sobemovirus in northern Nigeria. Nigerian J. Weed Sci. 2008, 21, 67–71. [Google Scholar]
  15. Velasco, L.; Ruiz, L.; Galipienso, L.; Rubio, L.; Janssen, D. A historical account of viruses in intensive horticultural crops in the Spanish Mediterranean Arc: New challenges for a sustainable agriculture. Agronomy 2020, 10, 860. [Google Scholar] [CrossRef]
  16. Sõmera, M.; Sarmiento, C.; Truve, E. Overview on sobemoviruses and a proposal for the creation of the family Sobemoviridae. Viruses 2015, 7, 3076–3115. [Google Scholar] [CrossRef] [Green Version]
  17. Zhu, Y.J.; Lim, S.T.S.; Schenck, S.; Arcinas, A.; Komor, E. RT-PCR and quantitative real-time RT-PCR detection of Sugarcane yellow leaf virus (SCYLV) in symptomatic and asymptomatic plants of Hawaiian sugarcane cultivars and the correlation of SCYLV titre to yield. Eur. J. Plant Pathol. 2010, 127, 263–273. [Google Scholar] [CrossRef]
  18. Galipienso, L.; Jassen, D.; Rubio, L.; Aramburu, J.M.; Velasco, L. Cucumber vein yellowing virus isolate-specific expression of symptoms and viral RNA accumulation in susceptible and resistant cucumber cultivars. Crop Prot. 2013, 43, 141–145. [Google Scholar] [CrossRef]
  19. Elvira-González, L.; Peiró, R.; Rubio, L.; Galipienso, L. Persistent Southern tomato virus (STV) interacts with Cucumber mosaic and/or Pepino mosaic virus in mixed- infections modifying plant symptoms, viral titer and small RNA accumulation. Microorganisms 2021, 9, 689. [Google Scholar] [CrossRef] [PubMed]
  20. Soler, S.; Debreczeni, D.E.; Vidal, E.; Aramburu, J.; López, C.; Galipienso, L.; Rubio, L. A new Capsicum baccatum accession shows tolerance to wild-type and resistance-breaking isolates of Tomato spotted wilt virus. Ann. Appl. Biol. 2015, 167, 343–353. [Google Scholar] [CrossRef] [Green Version]
  21. Falgueras, J.; Lara, A.J.; Fernández-Pozo, N.; Cantón, F.R.; Pérez-Trabado, G.; Claros, M.G. SeqTrim: A high-throughput pipeline for pre-processing any type of sequence read. BMC Bioinform. 2010, 11, 38. [Google Scholar] [CrossRef] [Green Version]
  22. Zheng, Y.; Gao, S.; Padmanabhan, C.; Li, R.; Galvez, M.; Gutierrez, D.; Fuentes, S.; Ling, K.S.; Kreuze, J.; Fei, Z. VirusDetect: An automated pipeline for efficient virus discovery using deep sequencing of small RNAs. Virology 2017, 500, 130–138. [Google Scholar] [CrossRef] [PubMed]
  23. Li, H.; Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 2009, 25, 1754–1760. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Zerbino, D.R.; Birney, E. Velvet: Algorithms for de novo short read assembly using de Bruijn graphs. Genome Res. 2008, 18, 821–829. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic local alignment search tool. J. Mol. Biol. 1990, 215, 403–410. [Google Scholar] [CrossRef]
  26. Higgins, D.; Thompson, J.; Gibson, T. CLUSTAL W: Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22, 4673–4680. [Google Scholar]
  27. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 2018, 35, 1547–1549. [Google Scholar] [CrossRef] [PubMed]
  28. Nei, M.; Kumar, S. Molecular Evolution and Phylogenetics; Oxford University Press: Oxford, UK, 2000. [Google Scholar]
  29. Le, S.Q.; Gascuel, O. An improved general amino acid replacement matrix. Mol. Biol. Evol. 2008, 25, 1307–1320. [Google Scholar] [CrossRef] [Green Version]
  30. Efron, B.; Halloran, E.; Holmes, S. Bootstrap confidence levels for phylogenetic trees. Proc. Natl. Acad. Sci. USA 1996, 93, 13429–13434. [Google Scholar] [CrossRef] [Green Version]
  31. Untergasser, A.; Cutcutache, I.; Koressaar, T.; Ye, J.; Faircloth, B.C.; Remm, M.; Rozen, S.G. Primer3—New capabilities and interfaces. Nucleic Acids Res. 2012, 40, e115. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Symptoms found in some C. esculentus plants: leaves with yellowing or chlorosis following a mosaic pattern in streaks and stunting in the aerial parts and roots with lower tuber production.
Figure 1. Symptoms found in some C. esculentus plants: leaves with yellowing or chlorosis following a mosaic pattern in streaks and stunting in the aerial parts and roots with lower tuber production.
Plants 11 02002 g001
Figure 2. (a) Genome organization of sobemoviruses [7] and the sequenced genome of xufa yellow dwarf virus (XYDV) with the contigs obtained with HTS shown below. (b) Unrooted Maximum Likelihood phylogenetic tree of coat-protein amino-acid sequences of members of genus Sobemovirus (Table 2) and XYDV. GenBank accession numbers and significant bootstrap values are indicated.
Figure 2. (a) Genome organization of sobemoviruses [7] and the sequenced genome of xufa yellow dwarf virus (XYDV) with the contigs obtained with HTS shown below. (b) Unrooted Maximum Likelihood phylogenetic tree of coat-protein amino-acid sequences of members of genus Sobemovirus (Table 2) and XYDV. GenBank accession numbers and significant bootstrap values are indicated.
Plants 11 02002 g002
Figure 3. Evaluation of XYDV detection techniques: (a) sensitivity assessed using serial dilutions of 10 ng/µL RNA extracts; (b) serial dilutions of plant crude extracts from two XYDV isolates, CM1 and CNM1. CMN1 RNA was used as an RT-PCR positive control. M = GeneRuler 1 kb plus ladder (Thermo Fisher Scientific, Waltham, MA, USA). Negative controls: H = RNA or crude extracts from C. esculentus healthy plants and W = water. No amplifications were obtained with RT-qPCR of the negative controls.
Figure 3. Evaluation of XYDV detection techniques: (a) sensitivity assessed using serial dilutions of 10 ng/µL RNA extracts; (b) serial dilutions of plant crude extracts from two XYDV isolates, CM1 and CNM1. CMN1 RNA was used as an RT-PCR positive control. M = GeneRuler 1 kb plus ladder (Thermo Fisher Scientific, Waltham, MA, USA). Negative controls: H = RNA or crude extracts from C. esculentus healthy plants and W = water. No amplifications were obtained with RT-qPCR of the negative controls.
Plants 11 02002 g003
Figure 4. Analysis of leaves from symptomatic and asymptomatic C. esculentus plants collected in the field: (a) percentage of positive samples; (b) relative quantification expressed in a logarithm scale.
Figure 4. Analysis of leaves from symptomatic and asymptomatic C. esculentus plants collected in the field: (a) percentage of positive samples; (b) relative quantification expressed in a logarithm scale.
Plants 11 02002 g004
Table 1. Oligonucleotides (primers and TaqMan probe) designed and used in this work.
Table 1. Oligonucleotides (primers and TaqMan probe) designed and used in this work.
UseOligoSequence
RT-PCR and Sanger479F5′-AAGATGTGATCCTCCAGCC-3′
sequencing402R5′-CAGCTTGGACCAGACAGAA-3′
541R5′-GGGTATATCTAGCGAAGT-3′
253F5′-TCAAATTTAGAGAGTCTGGTCG-3′
253R5′-CAGACTCTCTAAATTTGAC-3′
171F5′-AGGACGATTCCGCTTGATATC-3′
171R5′-ATATCAAGCGGAATCGTCCTTC-3′
148R5′-TGCAGTACGATCCAGATTTC-3′
220F5′-GTTAAACTTTAACGCTAGGAATG-3′
RT-PCRX1F5′-ACGACTTAGTCGTTGAAGC-3′
X1R5′-TCCGCGTATTCCCAGATAGC-3′
RT-qPCRqX1FACGTGCTTGATGCCGCTAAG
(SYBR Green)qX1RGGAACCTGTACCGCGGAGAT
RT-qPCRqX2F5′-GTGCAATGCGGGAAATCC-3′
(TaqMan probe)qX2R5′-AGCTTAGCGGCATCAAGCA-3′
Xprobe5′Fam-CCGTGTTGCTCACAGCATGGCA-Tamra3′
RT-LAMPXF35′-TTCCACCGTCTCCTACAG-3′
XB35′-TCCACACCTGCGTATGTA-3′
XFIP5′-GCGGAGAAGAATCTCACTCGGTGAACTTCAGTGGCTTGC-3′
XBIP5′-AGTGGCGACATTGCGATAGGCTGGTATGGTATGGTGACTGTAGTTG-3′
XLoopF5′-CTCGGAACTTCTGGTATCTCTG-3′
XLoopR5′-GTTGTGTATGACTCCGCTGA-3′
Table 2. Amino-acid identities between xufa yellow dwarf virus (XYDV) and viruses of genus Sobemovirus for ORFs 2a, 2b, and 3.
Table 2. Amino-acid identities between xufa yellow dwarf virus (XYDV) and viruses of genus Sobemovirus for ORFs 2a, 2b, and 3.
VirusAcronymGenBank2a2b3
Artemisia virus AArtVAJN62080232.558.324.4
Cocksfoot mottle virus CoMVAB04044725.151.311.0
Cymbidium chlorotic mosaic virus CyCMVLC01976428.651.420.5
Imperata yellow mottle virus IYMVAM99092821.151.310.8
Lucerne transient streak virus LTSVJQ78221326.754.517.6
Papaya lethal yellowing virus PLYVJX12331826.753.316.3
Physalis rugose mosaic virusPhyRMVMN78230027.751.316.1
Rice yellow mottle virus RYMVAJ60821026.153.617.6
Rottboellia yellow mottle virus RYMoVKC57746940.663.233.7
Rubus chlorotic mottle virus RuCMVAM94043726.052.318.1
Ryegrass mottle virus RGMoVEF09171439.661.135.1
Sesbania mosaic virus SeMVAY00429127.050.820.7
Solanum nodiflorum mottle virus SNMoVKC57747026.048.720.7
Southern bean mosaic virus SBMVDQ87559425.051.821.6
Southern cowpea mosaic virus SCPMVNC_00162527.752.819.6
Sowbane mosaic virus SoMVGQ84500226.051.919.4
Soybean yellow common mosaic virus SYCMVJF49512726.550.422.0
Subterranean clover mottle virus SCMoVAF20800127.754.918.0
Turnip rosette virus TRoVKC77872029.450.016.6
Velvet tobacco mottle virus VTMoVHM75426327.248.721.0
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Rubio, L.; Guinot-Moreno, F.J.J.; Sanz-López, C.; Galipienso, L. Discovery and Diagnosis of a New Sobemovirus Infecting Cyperus esculentus Showing Leaf Yellow Mosaic and Dwarfism Using Small-RNA High Throughput Sequencing. Plants 2022, 11, 2002. https://doi.org/10.3390/plants11152002

AMA Style

Rubio L, Guinot-Moreno FJJ, Sanz-López C, Galipienso L. Discovery and Diagnosis of a New Sobemovirus Infecting Cyperus esculentus Showing Leaf Yellow Mosaic and Dwarfism Using Small-RNA High Throughput Sequencing. Plants. 2022; 11(15):2002. https://doi.org/10.3390/plants11152002

Chicago/Turabian Style

Rubio, Luis, Francisco J. J. Guinot-Moreno, Carmen Sanz-López, and Luis Galipienso. 2022. "Discovery and Diagnosis of a New Sobemovirus Infecting Cyperus esculentus Showing Leaf Yellow Mosaic and Dwarfism Using Small-RNA High Throughput Sequencing" Plants 11, no. 15: 2002. https://doi.org/10.3390/plants11152002

APA Style

Rubio, L., Guinot-Moreno, F. J. J., Sanz-López, C., & Galipienso, L. (2022). Discovery and Diagnosis of a New Sobemovirus Infecting Cyperus esculentus Showing Leaf Yellow Mosaic and Dwarfism Using Small-RNA High Throughput Sequencing. Plants, 11(15), 2002. https://doi.org/10.3390/plants11152002

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop