Reassessment of the Columnea latent viroid (CLVd) Taxonomic Classification

Columnea latent viroid (CLVd) is a member of the Pospiviroid family and its naked circular RNA genome typically forms native “rod-like” secondary structures. In this work, the CLVd taxonomy was reevaluated based on sequence similarity and phylogenetic analysis, as well as the evaluation of the symptom development and disease severity of four selected CLVd isolates in a range of host species. The phylogenetic analysis showed that all CLVd isolates were clustered into five distinct clades: (I) severe isolates originally found in tomato crops in Thailand, (II) ornamental isolates, (III) mild isolates originally found in tomato crops in Thailand, and two clades (IV and V) containing mild isolates originating mainly from tomato crops in European countries, with different virulence levels on several hosts. Our analysis demonstrated that some CLVd isolates have a sequence similarity of less than 90% within the species taxon, as well as distinct biological characteristics (symptom development and virulence), both of which are important ICTV criteria for viroid classification. For these reasons, we propose that CLVd should be re-classified into at least three main taxonomic lineages: a “CLVd-tomato Asian lineage” (I), a “CLVd-tomato European lineage” (IV) and a “CLVd-ornamental European lineage” (II), plus two minor lineages (III and V), fitting the ICTV criteria.

Viroid classification and nomenclature is managed by the International Committee on Virus Taxonomy (ICTV) and is based on molecular features (e.g., rod-like conformations, the type of central conserved region (CCR), the presence or absence of a terminal conserved region (TCR), etc.), as well as two additional mandatory criteria-an arbitrary level of less than 90% sequence similarity to other related viroid species, and distinct biological properties, particularly host range and symptoms [16][17][18]. The current taxonomy of CLVd is not completely consistent with these main taxonomic criteria. Firstly, CLVd most likely arose from an RNA recombination (Hammond et al. 1989) and even if it is classified within the genus Pospiviroid, it shares its CCR and TCR with the members of the Hostuviroid genus. In addition, within this taxonomic species, some CLVd isolates have a sequence similarity of less than 90% within the species taxon, especially when Asian and European isolates are compared, as we show in this study. Moreover, some biological properties within this species, such as symptom severity, are quite different between isolates. Ever since CLVd was found in South East Asia in tomato production crops, very severe symptoms and high crop losses have been observed and reported, particularly in Thailand. The Thai CLVd-infected tomatoes showed very severe stunting, strong leaf rugosity and strong vein necrosis [5,13,14,19,20], combined with very high yield losses (around 50% reduction in fruit size) in several commercial tomato cultivars [13]. On the other hand, mild symptoms and much less crop losses have been reported from most of the European CLVd isolates [12,21,22]. Since the ninth report, only one CLVd strain, designated as CLVd-col (Columnea latent viroid-Columnea isolate) has been recognized by ICTV [16,21], whereas previously three separate strains, namely, CLVd-bru (Brunfelsia undulata), CLVd-nem (Nematanthus wettsteinii) and CLVd-col, were assigned to the CLVd species [8,9,21,23].
In contrast, several other viroids have been recognized as distinct isolates by ICTV in both the Pospiviroidae and Avsunviroidae families due to differences in biological characteristics such as host range and symptoms [16]. However, they still have over 90% genome sequence similarity, and hence can still be considered the same species. For instance, Citrus exocortis viroid (CEVd) is assigned into two host-types; CEVd-citrus and CEVd-tomato (Indian tomato bunchy top viroid). Two host-types have been assigned for Potato spindle tuber viroid (PSTVd), namely, PSTVd-intermediate and PSTVd-tomato. Interestingly, Hop stunt viroid (HpSVd, currently HSVd) was split into five "species", linked to the host-types in the previous (ninth) ICTV report; HpSVd-hop, HpSVd-citrus HpSVdcucumber (Cucumber pale fruit viroid), HpSVd-peach (Peach dapple viroid) (HpSVdpch) and HpSVd-plum (Plum dapple viroid) [16].
In this study, the taxonomy of the currently known CLVd isolates is re-evaluated based on sequence similarity and phylogenic tree analysis, as well as taking into account the symptom development and disease severity across several Solanaceae hosts, such as tomatoes, bolo maka, Thai round eggplant, chili pepper, bell pepper (Capsicum annuum) and cucumber. Based on our findings, we propose re-designating the CLVd taxonomic classification into at least three distinct lineages (or strains), with most Asia-related CLVd isolates belonging to a "CLVd-tomato Asian lineage", the European tomato isolates as a "CLVd-tomato European lineage" and the European ornamental isolates as a "CLVdornamental European lineage". CLVd infection was confirmed using CLVd-specific (c-CLVd-infect: TGCAGGGTCAG-GTGTGAACCAC and h-CLVd-infect: GCCATGCAAAGRAAAAAGAAYGGG) [24] and two generic pospiviroid RT-PCR detection tests [14,25]. Total RNA was extracted from 100 mg fresh leaves from viroid-infected plants using the Spectrum™ Plant Total RNA Kit (Sigma-Aldrich, Saint Louis, MI, USA) following the manufacturer's instructions.

Materials and Methods
cDNA was synthesized by means of the iScript™ cDNA Synthesis kit (Bio-Rad, Hercules, CA, USA). PCR reactions were performed using the FastStart™ Taq DNA Polymerase (Sigma-Aldrich, Saint Louis, MI, USA) under conditions described in the respective publications. All steps were performed according to the respective manufacturer's protocols in a T100™ Thermal Cycler (Bio-Rad, Hercules, CA, USA). The RT-PCR amplification products were analyzed by means of gel electrophoresis using a 2% agarose gel, and Midori Green Direct as the staining dye (Nippon Gene Co., Ltd., Düren, Germany). The results were visualized under blue light using an Azure c150 Gel Imaging Workstation (Azure Biosystems, Inc., Dublin, CA, USA). Again, all steps were performed according to the manufacturer's protocols. Subsequent Sanger sequencing of the amplicon was performed, outsourced at Macrogen (Amsterdam, The Netherlands) and after contig assembly in Bionumerics (Bionumerics, version 7.6.3; Bioinformatics software; Applied Maths: Sint-Martens-Latem, Belgium, 2020), a BLASTn search against the non-redundant Nucleotide database (https://blast.ncbi.nlm.nih.gov/Blast.cgi, accessed on 29 January 2018) was used to confirm the viroid identity.

Symptoms and Severity Rating
CLVd-infected plants were photographed, disease severity was assessed and symptoms were observed and recorded using symptom rating scale tables (Tables 2 and 3). The symptom rating scale was calculated by scoring the main typical symptoms, multiplied by a weight score (based on damage and the severity of the symptoms). Five dominant symptoms-yellowing, necrosis, distortion and stunting of apical shoot and stem-were used for symptom rating on tomato plants ( Table 2). For bolo maka, Thai round eggplant, chili/bell pepper and cucumber, only three main symptoms, namely, leaf distortion, necrotic lesion and stunting, were scored and assessed ( Table 3). Each of the selected and rated symptoms were ranked into a scale containing six levels (from 0 (no symptoms) to 5 (most severe)). The symptom rating scores from each selected host were averaged for all plants (3 to 5 plants per host). The score x weight was calculated based on "the score" multiplied by the "weight of score" (of each symptom). The overall symptom rating scale was calculated using the equation presented below. The overall symptoms ranged from 0 to 5 and an average score of 0-1 was indicated as "mild", a score 2-3 was "intermediate", and a score of 4-5 was "severe".
Average overall symptom rating = ∑(average score for a specific symptom × average score weight for that symptom)

100
(1) In total, 83 sequences of the Columnea latent viroid (CLVd) from the NCBI GenBank database (https://www.ncbi.nlm.nih.gov/, accessed on 30 January 2020) (Table A1) were selected for multiple-sequence alignment and phylogenic tree analysis. For all these CLVd sequences, a multiple-sequence alignment was created with MUSCLE (MEGA7, version 7.0.21, Molecular Evolutionary Genetics Analysis software, Kumar, Stecher, Li, Knyaz, and Tamura: 2018) using default parameters [26]. The phylogenetic tree was created in MEGA7 using neighbor-joining (bootstrap 1000 replicates). The alignment was read into R v4.0.3 (R, version 4.0.3, A language and environment for statistical computing, R Core Team: Vienna, Austria, 2021) using the Biostrings package, and pairwise similarities between the sequences were calculated using the Jukes-Cantor 69 model (assuming equal base frequencies) in the phangorn package. Finally, the sequences were clustered based on their pairwise similarity using hierarchical clustering, and the similarities were visualized in a heatmap using the gplots package.

CLVd Symptom Severity and Sequence Similarity Analysis
The results of the symptom severity assessment for the four selected CLVd isolates (2.1) were evaluated with respect to the CLVd genome sequence similarity analysis from 83 CLVd isolates. A pairwise identity matrix (%) of four selected CLVd isolates was created using Clustal2.1 (https://www.ebi.ac.uk/Tools/msa/clustalo/, accessed 30 January 2020).

Phylogenetic Analysis of Four Selected CLVd Isolates
The phylogenetic analysis of the sequences of the four selected CLVd isolates (14A, Chaipayon-1, LP1-15c4 and Solanum-1) showed that the three Thai CLVd isolates clustered into the same clade, separately from the European isolate ( Figure 1). The percentage identity matrix showed the distinction between European and Thai isolates, close to the taxonomic species-threshold of 90% (Table 4). Phylogenetic relation between four CLVd isolates. The European CLVd isolate is indicated with a green dot, whereas the Thai CLVd isolates are indicated with a red dot. PCFVd LPng19-3c1 (GenBank Accession No. JF446920) was used as an outgroup and is indicated with a blue diamond shape. The tree was constructed using MEGA version 7.0.21 with 1000 bootstrap replicates and a 50% cut-off value. All four CLVd isolates induced typical viroid symptoms on tomato cv. Rutgers, such as strong vein necrosis, leaf rugosity and epinasty and severe stunting. CLVd 14A showed the least symptoms compared to the three Thai isolates ( Figure 2). A score of 3.2 on the symptom rating scale was recorded for CLVd 14A, inducing intermediate symptoms in tomato cv. Rutgers (Table 5). CLVd Chaipayon-1 showed more severe symptoms (stunting, epinasty, rugosity and necrosis) than CLVd 14A, and scored an average of 4.8 on the symptom rating scale ( Figure 2, Tables 5 and 6), whereas both CLVd LP1-15c4 and Solanum-1 resulted in the most severe symptoms in tomato cv. Rutgers. These two CLVd isolates induced extreme stunting of the whole plant, whereas leaves showed strong yellowing, epinasty and rugosity. Strong obvious necrotic lesions were found in several plant parts, such as the midrib, leaf vein, petiole and stem ( Figure 2). Both isolates scored an average of 5.0 on the symptom rating scale (Table 5).   Tables 2 and 3 for details about  the symptom rating scale used. For more detailed information on symptoms, disease severity and  symptom rating scales for each host plant, please see Tables A3-A7). For tomato cv. Seeda 50, the overall symptoms were quite similar, although generally less severe than in the cv. Rutgers. Notably, bunchy top symptoms were less pronounced in this tomato cultivar. Again, CLVd 14A (score 2.6) showed the mildest symptoms compared to all Thai isolates, mostly inducing moderate stunting and necrotic lesions, although it showed less yellowing, epinasty and rugosity (Figure 3, Tables 5 and 6). CLVd Chaipayon-1 (score 3.8) caused extensive yellowing and necrosis in the whole plant; however, less epinasty and rugosity was observed. (Figure 3, Tables 5 and 6). Furthermore, on cvs. Seeda 50, CLVd LP1-15c4 and Solanum-1 induced the strongest symptoms, with again yellowing of leaves, very severe necrosis and stunting, as well as slightly less epinasty and rugosity compared to Rutgers ( Figure 3 and Table 6), which had a score of 4.5 (severe) for both CLVd isolates (Table 5).

Bolo Maka
In bolo maka, a high variation in symptom development and severity was observed. Generally, symptoms on bolo maka were quite similar to those observed on the tomato cultivars. The plants showed stunting of the whole plant, leaf epinasty, rugosity, yellowing and apical stunting. Similarly as observed in tomatoes, CLVd 14A induced the least symptoms, compared to the three Thai isolates. The isolate only induced slight growth reduction and very mild leaf distortion (epinasty and rugosity). No necrosis was observed in any of the infected bolo maka plants ( Figure 5 and Table 6). CLVd Chaipayon-1 caused stronger symptoms, especially leaf epinasty and rugosity with strong vein necrosis on the leaf vein, midrib, petiole and stem. CLVd Chaipayon-1 caused more severe stunting than CLVd 14A (around a 30% reduction in plant height) ( Figure 5 and Table 6). CLVd LP1-15c4 and Solanum-1 induced the most severe symptoms; including very strong stunting (more than a 70% reduction in plant height), severe leaf distortion showing epinasty, rugosity and size reduction. Vein necrosis was more distinct than for CLVd Chaipayon-1 ( Figure 5 and Table 6). The obtained disease severity scores for these four CLVd isolates are presented in Table 5.

Thai Round Eggplant
In Thai round eggplant, generally the typical symptoms were much less severe than in tomatoes and bolo maka. The most typical symptoms that were observed were overall mild stunting and yellowing, with only slight leaf epinasty and rugosity. In addition, necrotic lesions were rarely found, and if so, then it was particularly for the more severe isolates CLVd LP1-15c4 and Solanum-1, and this occurred late in the infection period. CLVd 14A caused almost no symptoms on Thai round eggplant. Only a minor leaf distortion was observed. CLVd Chaipayon-1 also caused fewer symptoms than CLVd 14A. The symptoms largely remained limited to a minor reduction in plant height. CLVd LP1-15c4 and Solanum-1 showed slightly stronger symptoms than CLVd 14A and Chaipayon-1 and caused leaf distortion and slight vein necrosis on the midrib, petiole and stem ( Table 6). The obtained disease severity scores for these four CLVd isolates are presented in Table 5.

Chili Pepper, Bell Pepper and Cucumber
There was no apparent symptom development on chili or bell pepper for any of the CLVd isolates (Table 6) (Table 6).

Severity Evaluation of Four CLVd Isolates
Overall, CLVd 14A was the mildest of all four selected isolates in this study. It showed intermediate symptoms on tomato cvs. Rutgers and Seeda 50, whereas it caused mild symptoms on tomato cv. Insaf, bolo maka and Thai round eggplant (Tables 5 and 6 and Figure 6). CLVd Chaipayon-1 can be considered as an isolate inducing intermediate symptoms.
Although it did result in quite severe symptoms on tomato plants, more intermediate symptoms were observed on bolo maka, and its infection remained rather mild on Thai round eggplant (Table 6). Because the symptoms of CLVd Chaipayon-1 were generally less severe than those of the isolates CLVd LP1-15c4 and Solanum-1, this CLVd isolate was considered to be intermediately virulent ( Figure 6 and Table 6). Finally, the highly virulent CLVd LP1-15c4 and Solanum-1 isolates resulted in an almost identical severity. These two isolates caused extreme symptoms on tomato cvs. Rutgers and Seeda 50 as well as on bolo maka, whereas they induced intermediate symptoms in tomato cv. Insaf and remained rather mild in Thai round eggplant (Tables 5 and 6 and Figure 6). Phylogenetic analysis showed that the sequence similarity of the four selected CLVd isolates is correlated with the symptom severity that they induce. Of the three Thai CLVd isolates, the two severe isolates (LP1-15c4 and Solanum-1) were the most closely related, yet they also clustered with the intermediate isolate (Chaipayon-1), whereas the CLVd 14A mild isolate was completely segregated from the others (Figure 6).

CLVd Sequence Analysis
The phylogenetic tree (Figure 7), created based on the full genome sequences of all CLVd isolates from GenBank, shows that five distinct clades, or groups, can be differentiated-group I mainly contains tomato isolates from Asia (represented with a light blue triangle), except isolates Solanum-1, Solanum-4 and Solanum-16, which were obtained from bolo maka. Most isolates from this clade induce severe symptoms on tomatoes. All three selected CLVd isolates (Chaipayon-1, Solanum-1 and LP1-15c4) also group into this clade. Group II includes all ornamental isolates (represented with an orange triangle). This clade consists of isolates from European countries, plus three isolates from the United States and a Canadian isolate.
Group III contains a minor group of tomato isolates from Asia (Thailand) that cluster distantly from the major Asian tomato isolate group (I). The isolates in this group were all obtained during the very first CLVd outbreak in Thailand, which was eradicated, since related variants have never been found again in the country.
Group IV consists of tomato isolates from Europe (represented with a dark blue triangle), mainly originating from the Netherlands and the United Kingdom (GenBank Accession No. KY810771), plus one isolate from Mali. One of the selected isolates in this study, CLVd-14A, originally found in tomato crops in the Netherlands, clusters in this group.
Finally, group V contains some remaining deviating minor tomato isolates from Europe that were collected in France, Portugal and Italy.
A multiple-sequence similarity matrix was created based on the 83 publicly available complete sequences of CLVd isolates and visually presents the pairwise genome sequence similarity within the viroid species "Columnea latent viroid". The genome sequences of CLVd showed the variation in sequence similarity across this viroid species and also reflect the five groups that were discussed above ( Figure 8). The CLVd genome sequence similarity analysis (one of two mandatory viroid taxonomic criteria) also reveals that most CLVd isolates, particularly those belonging to CLVd group V, show a sequence similarity of less than 90% (indicated with a red color in Figure 8) with isolates of the other four groups (Figure 8). This is quite atypical for virus and viroid taxonomy, since the level of sequence similarity needs to be greater than or equal to 90% in order for isolates to be classified as members of the same species.

CLVd Symptom Severity and Sequence Similarity Analysis
All viroids from group I induced clearly different symptoms and severity, especially in several Solanaceae hosts, compared to, e.g., CLVd 14A, the representative isolate from group IV that was intensively studied. They are highly virulent and cause important yield losses in tomato crops and bolo maka [5,14,15,19]. The symptom severity of three selected CLVd isolates (Chaipayon-1, Solanum-1 and LP1-15c4) of this group was demonstrated and discussed in Section 3.1.2.
Isolates from group II, the ornamental European lineage, were previously reported not to cause any symptoms on their (original) natural ornamental host species; Columnea erythrophae, Nemathanthus wettsteini, Brunfelsia undulate, Gloxinia gymnostoma, G. nematanthodes and G. purpurascens [8][9][10]12]. Isolates from this CLVd group induced mild to intermediate symptoms on tomato plants. They caused symptoms on tomato and potato crops similar to those of other pospiviroids; however, they proved to be less severe than the highly virulent isolates of potato spindle viroid (PSTVd) [12].
Isolates from group III, the minor tomato Asia lineage, induce mild to intermediate symptoms on several Solanaceae hosts (Figure 7). They only induced mild to intermediate symptoms on tomato plants, and almost no symptoms on bolo maka and Thai round eggplant [14,19] This is again in contrast with the isolates from group I, the major tomato lineage from Asia, inducing severe symptoms on all Solanaceae vegetable hosts such as tomato, bolo maka and Thai round eggplant [5,15,19] and intermediate symptoms on Thai round eggplant and hot chili [5,14,19]. CLVd 14A perfectly clusters in CLVd group IV, the major tomato European lineage, and was originally found in a tomato plant in the Netherlands. Group IV isolates generally show mild to intermediate symptoms on several plant species and induce typical pospiviroid symptoms on tomato and potato plants [7,27]. Only one observation of a CLVd outbreak in UK has been recorded, which showed no symptoms on the fruit and no obvious negative effects on quality [28].
The minor tomato European lineage, group V, also induced viroid-like symptoms on tomato crops-stunting, vein necrosis and leaf yellowing and malformation [29,30].
The isolates from group I showed a sequence similarity lower than 90% with CLVd group V and some members in group IV. They also share a relatively low sequence similarity with isolates from CLVd group II and III, at 93.2% and 92.8% on average, respectively. In addition, isolates from this group clearly induced distinct symptoms on bolo maka and yield loss in tomatoes, when compared to CLVd group II, III and IV. However, it is impossible to discuss the biological properties of CLVd group V, because no data about symptom expression on bolo maka and yield loss in tomatoes of this group have been reported. Furthermore, the isolates from group II, the ornamental lineage, showed less than 90% sequence similarity to members of the CLVd group V, and some isolates from group III and IV. On average, they shared 93.2% sequence similarity with the isolates in CLVd group I. This CLVd lineage induced mild to intermediate symptom virulence in tomatoes and almost no symptoms on bolo maka, which is similar to CLVd group III and IV. Group III, the minor tomato Asia lineage, showed less than 90% sequence similarity to some members of CLVd group II, and shared low sequence similarity to CLVd groups I, IV and V (92.8%, 91.8% and 92.5% on average, respectively). Group IV, the major tomato European lineage, showed less than 90% sequence similarity to CLVd group V and some members of CLVd groups I and II. It shared on average only 91.8% sequence similarity to CLVd group III. Group V, the minor tomato European lineage, showed less than 90% sequence similarity to CLVd groups I, II and IV, whereas it shared on average 92.5% sequence similarity to CLVd group III (Figure 8).

Conserved Motifs along the CLVd Genome
The CLVd genome sequence multiple alignment analysis showed particular variations of conserved motifs in the variable (V) domain between the different CLVd groups. In detail, most CLVd groups showed a distinct sequence conserved motif in the upper V domain, except CLVd groups III and V (minor Asia and minor tomato European groups) which both contain an identical sequence motif (Figure 9). In the lower V domain, most CLVd groups show differences in the conserved motif as well. The main tomato European lineage (group IV) contains a distinct conserved motif, whereas in the main Asian (group I) and the ornamental European lineage (group II), a similar conserved motif can be observed. Furthermore, CLVd groups III and V contain an almost identical sequence in the lower V domain as well (Figure 9). This may suggest that the minor Asia lineage (group III) and minor tomato European lineage (group V) are more related to each other than to the other CLVd groups, and that they share a similar evolutionary route. On the other hand, the major Asian (group I), ornamental European (group II) and major tomato European (group IV) lineages contain a totally different conserved motif at the V domain, which suggests a dissimilar evolutionary route.

Discussion
In this work, we investigated the virulence of four Columnea latent viroid (CLVd) isolates (Chaipayon-1, Solanum-1, LP1-15c4 and 14A) on different host plants, revealing clear variations in symptom severity. The imperfect current taxonomic classification of Columnea latent viroid species was studied using genome sequence similarity analysis of publicly available CLVd sequences and reevaluated in light of some biological properties, such as symptom expression and virulence.
CLVd is a poorly studied member of the pospiviroids, and information about its taxonomy, host range, transmission methods and impact are limited [21]. To date, only one research article has proposed re-classifying CLVd's taxonomy (from Pospiviroid into the Hostuviroid genus) based on high-throughput selective 29-hydroxyl acylation analyzed by primer extension (SHAPE) and RNA structure prediction software. However, in this study, only one CLVd reference isolate was used (GenBank Accession No. AY367350, including in CLVd group V) [31].
Viroid pathogenicity and disease severity are directly related to the viroid genome sequence, in combination with the host plant species and age, as well as environmental factors (e.g., light intensity and ambient temperature) [32,33]. Four different CLVd isolates were selected to evaluate symptom development and disease severity across a range of plant species-tomato cvs. Rutgers, Seeda 50 and Insaf; bolo maka; Thai round eggplant; chili pepper; bell pepper and cucumber. CLVd 14A, isolated from tomatoes in Europe, proved to be a very mild isolate, overall causing much milder symptoms in all selected plant species than the three Thai isolates [7,28]. The intermediate isolate CLVd Chaipayon-1, found on tomato plants in the northeastern region of Thailand [5,19], induced strong symptoms only on tomato cvs. Rutgers and Seeda 50, while causing moderate symptoms on tomato cv. Insaf and bolo maka. Lastly, the severe isolates CLVd LP1-15c4 and Solanum-1 showed very strong symptoms in most of the plant species; tomato cvs. Rutgers and Seeda 50, and bolo maka, while causing intermediate symptoms on tomato cv. Insaf. CLVd LP1-15c4 was also found in tomato plants in the northern region of Thailand (personal communication, Dr. Kanungnit Reanwarakorn, Kasetsart University), whereas Solanum-1 was found in bolo maka in central Thailand [5,15,19]. In addition, several Asian CLVd isolates have been reported to cause serious yield losses in tomato crops [5,13]. This is the first report about the correlation between CLVd genome sequences and symptom severities, and three levels of virulence for CLVd were determined.
Undoubtedly, tomato cv. Rutgers was the most susceptible host for all selected CLVd isolates, showing typical strong pospiviroid symptoms. Tomato cv. Seeda 50 also proved to be susceptible to CLVd, and showed strong typical symptoms. Notably, the specific symptom expression on tomato cv. Seeda 50 was slightly different from that observed on Rutgers; Seeda 50 showed stronger leaf yellowing and necrotic lesions but less leaf epinasty and rugosity. Tomato cv. Insaf was the least susceptible tomato host showing mild to moderate symptoms for all four selected CLVd isolates. Interestingly, bolo maka was susceptible to all Thai CLVd isolates, with strong symptoms that varied according to the isolate and which were very similar to the typical symptoms on tomatoes, such as strong stunting, yellowing, epinasty, rugosity and necrosis ( Figure 6). In contrast, Thai round eggplant was a less susceptible host for CLVd, showing mild symptoms for all selected CLVd isolates, such as mild leaf distortion and little stunting. Conversely, chili pepper and bell pepper remained symptomless for all selected CLVd isolates, yet all isolates were able to infect and replicate strongly in pepper over the entire lifespan of the plants. This is similar to most natural ornamental hosts-Columnea erythrophae, Nematanthus wettsteinii, Brunfelsia undulata, Gloxinia gymnostoma, G. nematanthodes, G. purpurascenswhich showed no symptoms either [3,[8][9][10]12]. Interestingly, CLVd Solanum-1 was the only isolate that succeeded in infecting cucumber, however, at a 50% infection rate with no symptoms. The literature about CLVd infectivity in cucumber is somewhat contradictory (Table A2) [6,7,[12][13][14]34]. In hindsight, cucumber is probably not a main natural host of CLVd.
Both the multiple-sequence similarity matrix and phylogenetic analysis indicated that CLVd genome sequences cluster into five groups; a major tomato Asia group (I), an ornamental European group (II), a minor tomato Asia group (III), a major tomato European group (IV) and finally a minor tomato European group (V), which are all lineages with a low percentage sequence similarity, even less than 90% for several isolates. This observation corresponds to the results of the symptom severity assessment, in which CLVd 14A, a European isolate, showed different levels of symptom virulence when compared to the three Asia isolates (Chaipayon-1, Solanum-1 and LP1-15c4-group I) on several plant species such as tomatoes, bolo maka and Thai round eggplant. Based on the current ICTV taxonomic classification criteria, these groups would not belong to the same species [16]. In addition, to date, the taxonomic relationship between the ornamental isolates and the Asian and European tomato isolates are still uncertain [21].
Analysis of the conserved motifs on the V domain showed that most of the distinct CLVd lineages (groups) are characterized by different conserved sequences in the motifs (except the two minor groups, CLVd groups III and V) (Figure 9). This relates to recent findings in a CLVd quasi-species population study across several Solanaceae hosts, which demonstrated a high number of SNPs on P and TR domains among most CLVd progeny variants (mutant populations). In contrast, TL, CCR and V domains showed a limited number of mutations, confirming those regions to be very conserved for this viroid [24]. Therefore, it is almost impossible to create a new variant that surpasses its own lineage, due to the relatively wide genomic distance between the CLVd main lineages. This indicates that the three main CLVd groups observed, the main Asian, ornamental and major tomato European lineages (CLVd groups I, II and IV, respectively), possibly evolved from a very early evolutionary split. The two minor tomato Asian and European groups (CLVd groups III and V) possibly evolved from the same evolutionary route, since these two groups shared several similar features, such as genome sequence similarity and the conserved motif on the V domain. In hindsight, the CLVd group III (minor Asia lineage) might possibly be a bridge (or intermediate) group between group V and the main three CLVd groups. In addition, the main Asia lineage (group I) meets the two distinct mandatory criteria for taxonomic reclassification-low sequence similarity and especially different biological properties (symptom expression, disease virulence and yield loss).
For these reasons, we propose that CLVd should be reclassified into at least three main taxonomic lineages-a "CLVd-tomato Asian lineage" (I), a "CLVd-tomato European lineage" (IV) and a "CLVd-ornamental European lineage" (II), plus two minor lineages (III and V). This is similar to several viroid species that have been assigned into several strains by ICTV, such as PLMVd, CEVd, PSTVd, HpSVd, CCCVd, ASSVd and CDVd [16]. In addition, one could argue that these groups could possibly be designated as three different viroid species within the Pospiviroid genus. This proposed CLVd taxonomic reclassification fits with the two mandatory criteria required by ICTV; (1) an arbitrary level of less than 90% sequence similarity and (2) distinct biological properties, particularly symptom virulence and severity.

Conclusions
In this work, we propose the reclassification of the species Columnea latent viroid into at least three major lineages, namely, the "CLVd-tomato Asian" (severe isolates originally found in tomato crops in Thailand), "CLVd-ornamental European" (mild isolates originally found in ornamentals) and "CLVd-tomato European" (mild isolates originally found in tomato crops in European countries) lineages. In addition, we identified two minor lineages (corresponding to groups III and V in the phylogenetic study; Figures 7 and 8) containing mild isolates from Thailand and comprising mild isolates from European countries, respectively, of which the latter has been identified as the most distinct lineage, even meriting reclassification into a new viroid species, according to the ICTV criteria.