Molecular and Biological Characterization of a New World Mono-/Bipartite Begomovirus/Deltasatellite Complex Infecting Corchorus siliquosus

The genus Begomovirus (family Geminiviridae) is the largest genus in the entire virosphere, with more than 400 species recognized. Begomoviruses are single-stranded DNA plant viruses transmitted by whiteflies of the Bemisia tabaci complex and are considered one of the most important groups of emerging plant viruses in tropical and subtropical regions. Several types of DNA satellites have been described to be associated with begomoviruses: betasatellites, alphasatellites, and deltasatellites. Recently, a family of single-stranded DNA satellites associated with begomoviruses has been created, Tolecusatellitidae, including the genera Betasatellite and Deltasatellite. In this work, we analyzed the population of begomoviruses and associated DNA satellites present in Corchorus siliquosus, a malvaceous plant growing wild in Central America, southeastern North America and the Caribbean, collected in Cuba. The genomes of isolates of two New World begomoviruses [(Desmodium leaf distortion virus (DesLDV) and Corchorus yellow vein Cuba virus (CoYVCUV)] and two deltasatellites [tomato yellow leaf distortion deltasatellite 2 (TYLDD2) and Desmodium leaf distortion deltasatellite (DesLDD)] have been cloned and sequenced from plants showing yellow vein symptoms. Isolates of one of the begomoviruses, CoYVCUV, and one of the deltasatellites, DesLDD, represent novel species. Experiments with infectious clones showed the monopartite nature of CoYVCUV and that DesLDD utilizes the bipartite DesLDV as helper virus, but not the monopartite CoYVCUV. Also, CoYVCUV was shown to infect common bean in addition to Nicotiana benthamiana. This is the first time that (i) a monopartite New World begomovirus is found in a host other than tomato and (ii) deltasatellites have been found in C. siliquosus, thus extending the host and helper virus ranges of this recently recognized class of DNA satellites.

The genus Begomovirus (family Geminiviridae) is the largest genus in the entire virosphere, with more than 400 species recognized. Begomoviruses are single-stranded DNA plant viruses transmitted by whiteflies of the Bemisia tabaci complex and are considered one of the most important groups of emerging plant viruses in tropical and subtropical regions. Several types of DNA satellites have been described to be associated with begomoviruses: betasatellites, alphasatellites, and deltasatellites. Recently, a family of single-stranded DNA satellites associated with begomoviruses has been created, Tolecusatellitidae, including the genera Betasatellite and Deltasatellite. In this work, we analyzed the population of begomoviruses and associated DNA satellites present in Corchorus siliquosus, a malvaceous plant growing wild in Central America, southeastern North America and the Caribbean, collected in Cuba. The genomes of isolates of two New World begomoviruses [(Desmodium leaf distortion virus (DesLDV) and Corchorus yellow vein Cuba virus (CoYVCUV)] and two deltasatellites [tomato yellow leaf distortion deltasatellite 2 (TYLDD2) and Desmodium leaf distortion deltasatellite (DesLDD)] have been cloned and sequenced from plants showing yellow vein symptoms. Isolates of one of the begomoviruses, CoYVCUV, and one of the deltasatellites, DesLDD, represent novel species. Experiments with infectious clones showed the monopartite nature of CoYVCUV and that DesLDD utilizes the bipartite DesLDV as helper virus, but not the monopartite CoYVCUV. Also, CoYVCUV was shown to infect common bean in addition to Nicotiana benthamiana. This is the first time that (i) a monopartite New World begomovirus is found in a host other than tomato and (ii) deltasatellites have been found in C. siliquosus, thus extending the host and helper virus ranges of this recently recognized class of DNA satellites.
Keywords: geminiviruses, begomoviruses, DNA satellites, deltasatellites, Corchorus siliquosus, Malvaceae, virus taxonomy, phylogenetic analysis INTRODUCTION Virus diseases that have emerged in the past three decades limit the production of important vegetable and fiber crops in tropical, subtropical, and temperate regions worldwide. Many of the causal viruses are transmitted by whiteflies (Hemiptera: Aleyrodidae), mainly by those belonging to the Bemisia tabaci cryptic species complex. Viruses known to be transmitted by whiteflies include members of the genera Begomovirus, Crinivirus, Ipomovirus, Torradovirus, and Carlavirus (Navas-Castillo et al., 2011), and two poleroviruses (genus Polerovirus) recently reported to be transmitted by B. tabaci (Ghosh et al., 2019;Costa et al., 2020).
The genus Begomovirus (family Geminiviridae) is the largest genus in the entire virosphere, with more than 400 species officially recognized. Begomoviruses are transmitted by whiteflies of the Bemisia tabaci complex Zerbini et al., 2017) and are considered among the most important emerging plant viruses. They infect a wide range of crops and wild plants in tropical and subtropical regions . Begomoviruses are circular singlestranded DNA plant viruses with twin (geminate) quasiicosahedral virions (Zerbini et al., 2017). Begomovirus genomes can be monopartite or bipartite depending on the presence of one or two (DNA-A and DNA-B) components, each of about 2.7 kb in size. Monopartite begomoviruses are widely spread in the Old World, but only a few examples have been reported from the New World. The genomes of monopartite begomoviruses resemble the DNA-A component. The DNA-A virion-sense strand encodes the coat and pre-coat proteins; the latter being present only in begomoviruses from the Old World. The DNA-A complementary-sense strand encodes the replication-associated protein (Rep), a transcriptional activator protein, a replication enhancer protein and C4 protein. DNA-B encodes a nuclear shuttle protein on the virion-sense strand and a movement protein on the complementary-sense strand. DNA-A and DNA-B share ∼200 nt in a common region (CR), located within the intergenic region that includes the replication origin. The Rep initiates viral DNA replication by binding to reiterated motifs (iterons) present in the CR and introducing a nick into the conserved non-anucleotide TAATATTAC. The CR shows a high level of nucleotide identity between both genome components of bipartite begomoviruses. Based on the phylogenetic analysis of full-length begomovirus genomic sequences (or DNA-A sequences), begomoviruses are classified into four lineages, Old World begomoviruses, New World begomoviruses, sweepoviruses and legumoviruses (Briddon et al., 2010).
Three types of DNA satellites have been described to be associated with begomoviruses, betasatellites (Briddon et al., 2003), alphasatellites (Briddon et al., 2004), and deltasatellites (Lozano et al., 2016). Recently, the International Committee on Taxonomy of Viruses (ICTV) modified the International Code of Virus Classification and Nomenclature to enable the classification of satellite nucleic acids (Adams et al., 2017). This change makes official the creation of a family of single-stranded DNA satellites associated with begomoviruses, Tolecusatellitidae Adams et al., 2017). The family Tolecusatellitidae includes two genera, Betasatellite and Deltasatellite. The genus Betasatellite includes DNA satellites about half the size of begomovirus genome components. Betasatellites encode the βC1 protein in the complementary-sense strand, which has important roles in symptom induction and suppression of transcriptional and post-transcriptional gene silencing (Zhou, 2013). Betasatellites have been found associated with monopartite begomoviruses in the Old World. Deltasatellites are non-coding DNA satellites associated with begomoviruses, of about a quarter size of begomoviral genome components (Fiallo-Olivé et al., 2012;Lozano et al., 2016). Tomato leaf curl deltasatellite, formerly ToLCV-sat, was the first DNA satellite identified in association with a plant virus, the monopartite Old World begomovirus tomato leaf curl virus (ToLCV) originating from Australia (Dry et al., 1997). Deltasatellites have been also found associated with bipartite New World begomoviruses infecting malvaceous weeds (Fiallo-Olivé et al., 2012 and sweepoviruses (Hassan et al., 2016;Lozano et al., 2016;Rosario et al., 2016). Also, the sequences of two uncharacterized deltasatellites from Philippines (KF433066) and India (AJ968684) have been deposited in GenBank. All deltasatellites contain common features; they share a small region with some sequence identity to a conserved region present in the betasatellites, an A-rich sequence, a predicted stemloop structure containing the nonanucleotide TAATATTAC, and a predicted secondary stem-loop (Fiallo-Olivé et al., 2012;Lozano et al., 2016).
The genus Deltasatellite comprises 11 species Adams et al., 2017), five of them reported associated with bipartite New World begomoviruses. Until now, New World deltasatellites have been found only in the malvaceous plants Malvastrum coromandelianum and Sidastrum micranthum and only associated with two helper begomoviruses, Sida golden yellow vein virus and tomato yellow leaf distortion virus (Fiallo-Olivé et al., 2012. Recently, biological evidence has been obtained that New World deltasatellites depend on a limited range of begomoviruses for maintenance in planta . In addition to the helper bipartite begomoviruses that the New World deltasatellites are associated with in nature, they can be transreplicated by the monopartite New World begomovirus tomato leaf deformation virus. In contrast, they cannot be maintained by the Old World begomoviruses tomato yellow leaf curl virus, tomato yellow leaf curl Sardinia virus and African cassava mosaic virus or the curtovirus beet curly top virus. New World deltasatellites do not affect the symptoms induced by the helper begomoviruses but in some instances they are able to reduce virus accumulation. Also, one New World deltasatellite has been shown to be transmitted by the whitefly B. tabaci, the natural vector of begomoviruses . In this study we analyzed the population of begomoviruses and associated DNA satellites present in wild Corchorus siliquosus plants collected in Cuba. The genomes of two New World begomoviruses (one bipartite and one monopartite) and two deltasatellites have been cloned and sequenced from plants showing yellow vein symptoms, the monopartite begomovirus and one of the deltasatellites representing novel species. Experiments with infectious clones showed that the monopartite begomovirus is also able to infect common bean and that the novel deltasatellite is maintained only by the bipartite begomovirus. This is the first time that a monopartite begomovirus native to the New World and deltasatellites have been found in C. siliquosus, extending the ranges of plant hosts and helper viruses for these recently recognized DNA satellites.

Plant Samples
Leaves of Corchorus siliquosus plants (family Malvaceae) showing yellow vein symptoms (Figure 1) were collected in the province of Matanzas, Cuba, on December 2013 ( Table 1). Each sample consisted of a few leaves which were dried, transported to the laboratory and held at 4 • C. Morphological identification of the plants was confirmed at the molecular level by DNA barcoding using chloroplast rbcL and matK genes (Hollingsworth et al., 2009).

DNA Extraction and Cloning
Total DNA was extracted from leaf samples using a cetyltrimethylammonium bromide-based method (Permingeat et al., 1998) and then used as a template in rolling-circle amplification (RCA) using φ29 DNA polymerase (TempliPhi kit, GE Healthcare). Amplified RCA products were first digested with HpaII, a 4-bp restriction enzyme, for preliminary restriction fragment length polymorphism (RFLP) analysis on 1.5% agarose gels. Then, RCA products were digested with a set of 6-bp restriction enzymes (EcoRI, BamHI, HindIII, NcoI, SacI). Putative full-length begomovirus genome components (∼2.6 kb) and deltasatellites (∼0.7 kb) obtained after digestion were sizeselected and cloned into the appropriate sites of pBluescriptII SK(+) (Stratagene) or pGEM-T Easy Vector (Promega). Recombinant plasmid DNAs were transformed into Escherichia coli DH5α by electroporation and selected clones were sequenced in both directions by the Sanger method using M13 forward and reverse [for pBluescriptII SK(+)] or T7 promoter and SP6 promoter (for pGEM-T Easy Vector) primers and then by primer-walking (Macrogen Inc., Seoul, South Korea).

Sequence Analysis
Sequences were assembled with SeqMan software included in the DNASTAR package (Lasergene). Open reading frames were identified using Open Reading Frame Finder (NCBI) and confirmed by a BLAST analysis 1 of their deduced amino acid sequences. Initial sequence identity comparison was performed using the BLAST program, sequences were aligned with MUSCLE (Edgar, 2004), and pairwise identity scores were calculated using SDT (Sequence demarcation tool) (Muhire et al., 2014). In silico digestion analysis of cloned begomovirus and deltasatellite sequences was performed with SnapGene (Insightful Science, available at snapgene.com). 1 https://blast.ncbi.nlm.nih.gov/Blast.cgi

Phylogenetic Analysis
For phylogenetic analysis of begomoviruses and deltasatellites, the maximum likelihood method was used after selecting the best-fit model of nucleotide substitution based on corrected Akaike Information Criterion and Bayesian Information Criterion as implemented in MEGA 7 (Kumar et al., 2016). For analysis of begomoviruses, the most closely related DNA-A and DNA-B components were selected as well as sequences belonging to different clades of New World begomoviruses. For phylogenetic analysis of deltasatelites, the deltasatellites isolates described in this work and one isolate belonging to each deltasatellite species Adams et al., 2017) were included.

Recombination Analysis
Recombination analysis of sequences reported in this study was performed using RDP4 (Martin et al., 2015) after alignment with MUSCLE (Edgar, 2004). Only recombination events detected with at least five methods were considered. To carry out the recombination analysis of the begomoviruses, the sequences with highest SWeBLAST scores (window size of 200 and step size of 200) were selected (Fourment et al., 2008). For deltasatellites, all the sequences available in GenBank were included in the analysis.

Construction of Infectious Clones
Infectious dimeric clones of the begomovirus genomic components and one of the deltasatellites characterized in this work were constructed as previously described . Inserts from monomeric clones of each viral component or deltasatellite were released from the plasmids, religated and subjected to RCA. RCA products were partially digested to produce dimeric molecules which were subsequently cloned in a plasmid vector. The inserts of the dimeric clones were excised and subcloned in the binary vector pCAMBIA0380. Details of the restriction enzymes and vectors used in each cloning step are given in Supplementary Table 1. Ligation reactions were transformed into Escherichia coli DH5α by electroporation (25 µF, 200 , 2500 V) in a Gene Pulser XCell (Bio-Rad). In each cloning step, clones were verified by digestion. Clones in the binary vector with inserts of the expected size were sequenced (Macrogen Inc., Seoul, South Korea) and transformed into Agrobacterium tumefaciens C58C1 by electroporation (25 µF, 200 , 2500 V).

Plant Inoculation and Begomovirus and Deltasatellite Detection
For agroinoculation assays, A. tumefaciens cultures harboring each dimeric construct were grown and inoculated as previously described . Plant species used for agroinoculation assays were Nicotiana benthamiana, Corchorus olitorius (seeds purchased from B&T World Seeds), common bean (Phaseolus vulgaris) cv. Donna, tomato (Solanum lycopersicum) cv. Moneymaker, pepper (Capsicum annuum) cv. California Wonder and zucchini (Cucurbita pepo) cv. Milenio. Plants were inoculated with A. tumefaciens cultures containing clones of viral DNA components and deltasatellites  For biolistic inoculation, the monomeric clone of CoYVCuV obtained in this work was used to inoculate N. benthamiana and C. olitorius plants. Viral DNA insert was released from the vector by digestion with SacI, circularized using T4 DNA ligase (Roche Diagnostics) and then subjected to RCA. The RCA product was precipitated onto 1.0 µm gold microcarriers (Bio-Rad) which were resuspended in ethanol. Plants were inoculated with a "Bim-Lab" instrument (Bio-Oz Biotechnologies). Both third and fourth leaves of each plant were shot twice with 5 µL of the RCA product-coated microcarriers. Biolistic inoculation experiments were repeated twice and each contained 12 (N. benthamiana) or 10 (C. olitorius) plants. Plants that served as negative controls were mock inoculated with the microcarriers in ethanol.
Apical leaves were used for tissue blots of petiole crosssections (tissue printing) performed on positively charged nylon membranes (Roche Diagnostics) at 33 days postinoculation. In the case of N. benthamiana, stem crosssections were employed instead. Hybridization was carried out using digoxigenin-labeled DNA probes specific to each viral genomic component and deltasatellites. The probes were prepared by PCR with primers specifically designed for each begomovirus genomic component and one of the deltasatellites (Supplementary Table 2) according to the DIG-labeling detection kit (Roche Diagnostics). Hybridization was carried out under high stringency conditions [washing steps at 65 • C in 0.1× SSC (15 mM NaCl, 1.5 mM sodium citrate) and 0.1% sodium dodecyl sulfate] following standard procedures. Hybridization signals were detected on X-ray film after treatment with CDP-Star (Roche Diagnostics).

A Bipartite and a Monopartite
Begomovirus, the Latter Representing a Novel Species, Infect Corchorus siliquosus in Cuba Digestion with HpaII of RCA-amplified products from the five C. siliquosus leaf samples yielded restriction patterns (Supplementary Figure 1) supporting the suspected begomovirus infections based on the symptomatology observed in the field, consisting of yellow veins (Figure 1). The complexity of the restriction patterns was suggestive of mixed infections. Sequencing of cloned DNA fragments following digestion with 6-bp enzymes confirmed the presence of DNA-A and DNA-B genomic components in the five samples. In total, 15 DNA-A and 21 DNA-B full-length clones were sequenced ( Table 2). All of them showed the organization of typical New World begomoviruses.
Pairwise comparison analysis carried out with SDT showed that 10 of the DNA-As (isolated from the five samples as BamHI or SacI fragments, GenBank acc. no. MF773881-MF773890) ( Table 2) showed a nucleotide identity of 97.3-100% between them and the highest identity (94.1-94.3%) with Desmodium leaf distortion virus (DesLDV) DNA-A (DQ875870). In accordance with the current taxonomic guidelines for the genus Begomovirus (a new DNA-A sequence with ≥91% pairwise identity to an available begomovirus DNA-A sequence will belong to the same species) (Brown et al., 2015), the isolates described here belong to the species Desmodium leaf distortion virus. DesLDV is a bipartite New World begomovirus previously found infecting only the wild fabaceous Desmodium glabrum in Yucatán, Mexico (Hernández-Zepeda et al., 2009). DNA-Bs (also isolated from the five samples as EcoRI, HindIII or SacI fragments, MF773891-MF773911) ( Table 2) showed nucleotide identity of 90.2-100% between them and the highest identity (87.9-88.9%) with DesLDV DNA-B from Mexico (DQ875871). A common region of 161 nt was identified in the cloned DNA-As and DNA-Bs of DesLDV which showed an identity of 93.8-96.9% between both components. Also, four copies of identical iterons, two of them inverted, were present in both genome components (Figure 2), thus supporting that they constitute cognate pairs. Phylogenetic analysis showed that both DNA- The other five full-length DNA-A components (isolated from samples 705 and 706 and cloned as EcoRI or SacI fragments, MF773912-MF773916) ( Table 2) showed an identity of 99.8-100% between them and only 76.6-77.4% with the cloned DesLDV isolates. The highest identity (81.9%) was with the DNA-A of two begomoviruses from the Caribbean, tobacco yellow crinkle virus (FJ213931) (Fiallo-Olivé et al., 2009) and Jatropha mosaic virus (KJ174331) . In accordance with the above-mentioned taxonomic guidelines for the genus Begomovirus, the isolates described here represent a novel species for which we propose the name Corchorus yellow vein Cuba virus (CoYVCUV). Phylogenetic analysis showed that the DNA-As of CoYVCUV isolates were highly divergent and grouped in a single cluster with a high bootstrap value, thus supporting that this begomovirus constitutes a novel species (Figure 3A). Recombination analysis of a representative CoYVCUV DNA-A (Cuba-Corchorus 706-1-2013), also showed its recombinant nature. The three recombination events detected involved begomoviruses infecting wild plants, Jatropha mosaic virus, Dalechampia chlorotic mosaic virus and DesLDV from Dominican Republic, Venezuela and Mexico, respectively (Hernández-Zepeda et al., 2009;Fiallo-Olivé et al., 2013;Melgarejo et al., 2015; Table 3).
The CR of the DesDLDV DNA-B and the equivalent region of the CoYVCUV DNA-A present in the same samples showed a low identity (65.8-68.0%) and dissimilar iterons (Figure 2). Based on this, the DNA-A of CoYVCUV and the DNA-B of DesLDV would not constitute a cognate pair. On the other hand, although the DNA-A of CoYVCUV showed the typical genome organization of a bipartite New World begomovirus DNA-A, no DNA-B could be found for this virus in the two samples analyzed despite repeated attempts to clone full-length or partial genome components different from the described above. This was done using the above-mentioned set of restriction enzymes as well as other 6-bp restriction enzymes whose recognition sequences were not present in the cloned genome components (data not shown). Also, in silico digestion analysis of cloned begomovirus and deltasatellite (see next section) sequences produced a RCA-RFLP pattern from samples 705 and 706 similar to those experimentally obtained (Figure 4). This analysis supports that no additional begomovirus genome components are present in an amplifiable amount in these samples. Taking together, these results suggested the possibility that CoYVCUV could be a monopartite New World begomovirus.
To assess the putative monopartite nature of CoYVCUV, a dimeric clone obtained for its DNA-A component was agroinoculated alone or in combination with a dimeric clone of DesLDV DNA-B in Nicotiana benthamiana and Corchorus olitorius plants. In four independent experiments, tissue printing hybridization of apical leaves showed that CoYVCUV was infectious in N. benthamiana, all the inoculated plants becoming infected (Table 4 and Supplementary Figure 2A) although they remained asymptomatic (data not shown). This included the two experiments in which CoYVCUV DNA-A was co-inoculated with DesLDV DNA-B, whose replication was not supported based on the tissue printing hybridization analysis (Table 4 and Supplementary Figure 2A). In contrast, C. olitorius plants were not infected in any of the four independent experiments (Table 4 and Supplementary Figure 2). Biolistic inoculation of C. olitorius plants was also assayed with CoYVCUV DNA-A for which RCA products obtained from circularized monomer genomes were employed. In two independent experiments, C. olitorius was not infected as determined by tissue printing hybridization analysis ( Table 4 and Supplementary Figure 2B). However, about 50% of the N. benthamiana plants inoculated in the same experiments were infected (Table 4 and Supplementary Figure 2B) but remained asymptomatic.  To further confirm the monopartite nature of CoYVCUV and to get insight on its possible host range, several plant species were inoculated with the CoYVCUV infectious clone. Two independent assays showed that CoYVCUV DNA-A was infectious in common bean, accumulating at detectable levels in apical leaves (Table 4 and Supplementary Figure 2A), but not in tomato, pepper and zucchini. Symptoms developed in common bean plants consisted of leaf crinkle and mild yellowing (Figure 5). These results confirm that CoYVCUV is a monopartite New World begomovirus.

Two Deltasatellites, One Representing a Novel Species, Are Associated With the Begomoviruses Infecting Corchorus siliquosus
Eighteen deltasatellite molecules were cloned and sequenced from three of the five C. siliquosus samples showing yellow vein symptoms analyzed in this work ( Table 2). Identification of such molecules as deltasatellites was based on their size (665-692 nt) and genome organization, including the features common to all deltasatellites: a small region with some sequence identity to the conserved region present in the betasatellites, an A-rich sequence, a predicted stem-loop structure containing the non-anucleotide TAATATTAC, and a predicted secondary stemloop. Pairwise comparison between these deltasatellites and all the deltasatellite sequences available in GenBank showed that the deltasatellites from C. siliquosus can be readily separated into two groups consisting of those molecules isolated from sample 680 (n = 3, 99.7-100% nucleotide identity between them, cloned as SacI fragments, MF773917-MF773919) and samples 704 and 705 (n = 15, 96.4-100% nucleotide identity between them, cloned as HindIII fragments, MF773920-MF773934), respectively. Nucleotide identity between isolates from both groups are only 73.0-74.9%. According to the proposed 91% deltasatellite species demarcation criteria , the deltasatellites isolated in this work belong to two different species. The three deltasatellites isolated from sample 680 belong to the species Tomato yellow leaf distortion deltasatellite 2, as having the highest nucleotide identity of 97.4% with the only previously known isolate of this species ( Phylogenetic analysis of all deltasatellites known to date, including those described in this work, showed two main clusters (Figure 6). One of these clusters contains deltasatellites associated with bipartite New World begomoviruses, Sida golden yellow vein deltasatellite 1 (SiGYVD1), SiGYVD2, SiGYVD3, tomato yellow leaf distortion deltasatellite 1 (ToYLDD1), ToYLDD2 and the novel deltasatellite described in this study, DesLDD. All DesLDD isolates clustered together and are separated from the rest of New World deltasatellites, supporting that this deltasatellite constitutes a novel species. No recombination events were detected in the genomes of the deltasatellites described here.

Desmodium Leaf Distortion Deltasatellite Is Replicated by the Bipartite Desmodium Leaf Distortion Virus but Not by the Monopartite Corchorus Yellow Vein Cuba Virus
Agroinoculation experiments were carried out to asses which begomovirus could act as a helper for DesLDD. For this, a dimeric clone obtained for DesLDD was co-inoculated on N. benthamiana and C. olitorius plants with dimeric clones of DesLDV DNA-A and DNA-B or of CoYVCUV. As it has been stated above, DesLDV and CoYVCUV did not infect C. olitorius in two different experiments. Thus, as expected, DesLDD did not infect this plant. In contrast, both begomoviruses were able to infect N. benthamiana (Table 5). In this host, tissue printing hybridization of petiole cross-sections showed that all    (Figure 7). The presence of DesLDD did not affect the symptoms caused by DesLDV (Figure 7). The fact that DesLDV is the helper virus of this deltasatellite supports the proposal to name the novel species Desmodium leaf distortion deltasatellite.

DISCUSSION
In this work, a complex of begomoviruses and associated deltasatellites has been characterized from C. siliquosus plants showing yellow vein symptoms in Cuba. Isolates of two New World begomovirus species, Desmodium leaf distortion virus and Corchorus yellow vein Cuba virus have been cloned and sequenced. DesLDV (DNA-A and DNA-B) was found in single infection in several samples, whereas CoYVCUV (DNA-A) was found only in mixed infection with DesLDV. DNA-A from both begomoviruses had the typical genome organization of bipartite New World begomoviruses, but despite repeated attempts, a DNA-B could not be detected for CoYVCUV. Thus, the question remained whether CoYVCUV could be a monopartite begomovirus or it would utilize the DNA-B of DesLDV. Although CoYVCUV DNA-A and DesLDV DNA-B have different iterons, the possibility of forming pseudorecombinants cannot be ruled out. A few cases are known of both New World and Old World begomovirus pairs that form viable pseudorecombinants, at least experimentally, even if their iterons are not identical (Andrade et al., 2006;Chakraborty et al., 2008). Agroinoculation assays carried out in this work showed that (i) CoYVCUV DNA-A alone was infectious in N. benthamiana and common bean plants and (ii) DesLDV DNA-B is not replicated in the presence of CoYVCUV DNA-A in N. benthamiana plants.
Although DNA-As of some bipartite begomoviruses have been shown to be able to infect N. benthamiana , the fact that CoYVCUV DNA-A was able to infect common bean definitely revealed the monopartite nature of this begomovirus. Although comparison of experimental RCA-RFLP data and in silico digestion analysis suggested that no additional begomovirus genome components are present in the studied samples, it does not escape our attention that a deep sequencing approach would have precisely defined if other viruses and/or subviral molecules were present at low concentration.
New World begomoviruses are classically bipartite (Brown et al., 2015;Zerbini et al., 2017). However, in 2011 a monopartite begomovirus native to the Americas was described infecting tomato plants in Peru and Ecuador (Márquez-Martín et al., 2011;Melgarejo et al., 2013;Sánchez-Campos et al., 2013). Later, four additional monopartite begomoviruses have been described from the Americas, all of them infecting tomato: tomato latent virus (TLV) in Cuba (Fuentes et al., 2016), tomato mottle leaf curl virus and tomato leaf curl purple vein virus in Brazil (Vu et al., 2015;Macedo et al., 2018), and tomato twisted leaf virus in Venezuela (Romay et al., 2019). TLV is in fact a recombinant consisting of a moiety from the monopartite Old World tomato yellow leaf curl virus (introduced into the Americas in the early 1990s) and a moiety from an unidentified bipartite New World begomovirus DNA-A (Fuentes et al., 2016). The results obtained in this work expands the number of New World begomoviruses FIGURE 6 | Phylogenetic tree illustrating the relationships of Desmodium leaf distortion deltasatellite (DesLDD) (bold in blue) and tomato yellow leaf distortion deltasatellite 2 (ToYLDD2) (bold in red) isolates obtained in this work with representative isolates of other deltasatellite species. The tree was constructed by the maximum-likelihood method (1000 replicates) with the MEGA 7 program using the best fit model, T92 + G. CrYVD, Croton yellow vein deltasatellite; MaLCuD, Malvastrum leaf curl deltasatellite; SiGYVD1, Sida golden yellow vein deltasatellite 1; SiGYVD2, Sida golden yellow vein deltasatellite 2; SiGYVD3, Sida golden yellow vein deltasatellite 3; SPLCD1, sweet potato leaf curl deltasatellite 1; SPLCD2, sweet potato leaf curl deltasatellite 2; SPLCD3, sweet potato leaf curl deltasatellite 3; ToLCD, tomato leaf curl deltasatellite; ToYLDD1, tomato yellow leaf distortion deltasatellite 1; ToYLDD2, tomato yellow leaf distortion deltasatellite 2. Major deltasatellite clusters associated with bipartite New World begomoviruses (NW) and monopartite Old World begomoviruses and sweepoviruses (OW and SWEEPO) are indicated. The bar below the tree indicates the number of nucleotide substitutions per site. Nodes with bootstrap values higher than 50% are marked with filled circles. whose monopartite nature has been shown experimentally to six, which is a very low number considering that there are about 170 New World begomovirus species recognized to date. It is generally believed that the first begomoviruses were monopartite (Briddon et al., 2010) and that the bipartite genome likely evolved before continental drift because bipartite begomoviruses occur in the Old and New World (Rojas et al., 2005). It has been suggested that the low number of indigenous monopartite begomoviruses in the New World could be explained by the predominance of bipartite progenitors in the areas of Pangea that gave rise to the Americas or by a later introduction (Melgarejo et al., 2013). The emergence of monopartite begomoviruses infecting tomato in the New World has been hypothesized to represent an evolutionary process which could have occurred due to a combination of factors including: (i) the extensive cultivation of this crop in which monopartite begomoviruses could evolve on their own, (ii) the existence of a local progenitor bipartite begomovirus in wild or cultivated plants, and (iii) the introduction of this local progenitor begomovirus into tomato crops by the invasive and polyphagous Middle East-Asia Minor 1 species (formerly B biotype) of the Bemisia tabaci complex (Melgarejo et al., 2013). In contrast to the rest of monopartite New World begomoviruses, CoYVCUV was not found infecting tomato but a wild malvaceous plant and tomato plants could not be experimentally infected by agroinoculation either. Thus, this could represent an additional evolutionary path in generating monopartite New World begomoviruses which could emerge in cultivated plants. Also, the facts that common bean became experimentally infected by CoYVCUV and showed symptoms of leaf crinkle and mild yellowing are relevant since this important crop could be threatened by this begomovirus in Cuba.
Since we were not successful in obtaining C. siliquosus seeds, C. olitorius, one of the two cultivated species in the genus Corchorus, was used in the agroinoculation and biolistics assays. C. olitorius was not infected by DesLDV or CoYVCUV, but it should be mentioned that this species is only distantly related to C. siliquosus (Benor, 2018).
Other partners in the viral complex found in C. siliquosus were subviral molecules, deltasatellites, which are the smallest DNA satellites associated with begomoviruses known to date. Deltasatellites have been previously found associated with two bipartite New World begomoviruses in Cuba, Sida golden yellow vein virus infecting M. coromandelianum and tomato yellow leaf distortion virus infecting S. micranthum (Fiallo-Olivé et al., 2012. Here, we have described the presence of two deltasatellites in a third malvaceous plant species, C. siliquosus, also from Cuba. One of these deltasatellites, Desmodium leaf distortion deltasatellite (DesLDD), constitutes a novel species, the twelfth of the genus Deltasatellite. This finding extends the ranges of natural hosts and helper viruses for these subviral molecules only recently recognized as a novel class of DNA satellites included in the family Tolecusatellitidae (Lozano et al., 2016;Adams et al., 2017).
Agroinoculation assays also showed that the novel deltasatellite described here, DesLDD, accumulated in N. benthamiana plants in the presence of DesLDV, but not CoYVCUV. Also, the presence of the deltasatellite did not modify the symptoms of leaf distortion and significant reduction in plant growth caused by DesLDV. This result is in agreement with those obtained with the two New World deltasatellites for which infectious clones were available to date, tomato yellow leaf distortion deltasatellite 2 and Sida golden yellow vein deltasatellite 1. Neither of the deltasatellites influenced the symptoms caused by their natural helper begomoviruses in N. benthamiana or in their natural hosts . In contrast, deltasatellites associated with sweet potato leaf curl virus, a sweepovirus, usually reduced symptom severity caused by its natural helper virus and by tomato leaf curl virus, a monopartite Old World begomovirus (Hassan et al., 2016).
Recombination analysis of representative isolates of the two begomoviruses described in this study showed the recombinant nature of both DNA-As. With the exception of cabbage leaf curl virus (Florida isolate) (Abouzid et al., 1992), all begomoviruses involved in the detected recombination events infect wild plants (C. siliquosus, Desmodium glabrum, Jatropha sp., Dalechampia sp., and Boerhavia diffusa) in countries from the Caribbean and surrounding areas, from Venezuela to Florida (Hernández-Zepeda et al., 2007Fiallo-Olivé et al., 2013;Polston et al., 2014;Melgarejo et al., 2015).
In this work, an additional begomovirus-deltasatellite complex has been characterized in a wild malvaceous plant, with known and novel partners which (i) have been shown to have interacted via recombination with other begomoviruses infecting both cultivated and wild plants and (ii) are able to infect new crops as common bean. In addition, the discovery of a deltasatellite that constitutes a novel species increases the complexity of a recently established, and still poorly understood, class of DNA satellites. This in an example which illustrates the phenomenon of emergence of plant viruses at the interface between natural vegetation and cultivated plants, which can cause outbreaks in nearby crops as well as in new areas driven by whitefly transmission, climate change and global trade.

DATA AVAILABILITY STATEMENT
The datasets generated for this study can be found in the online repositories. The names of the repository/repositories and accession number(s) can be found here: https://www.ncbi.nlm. nih.gov/genbank/, MF773881-MF773934.

AUTHOR CONTRIBUTIONS
EF-O and JN-C collected the Corchorus siliquosus leaf samples in the frame of the Ibero-American Network on Integrated Management of Vegetable Virus Diseases (CYTED 111RT0433), analyzed the data, and wrote the manuscript. EF-O and JN-C conceived the experiments that were performed by EF-O. Both authors contributed to the article and approved the submitted version.