# New World Cactaceae Plants Harbor Diverse Geminiviruses

Rafaela Fontenele, Andrew Salywon, Lucas Majure, Ilaria Cobb, Amulya Bhaskara, Jesús Avalos-Calleros, Gerardo Argüello-Astorga, Kara Schmidlin, Anthony Khalifeh, Kendal Smith, et al.

► To cite this version:

Rafaela Fontenele, Andrew Salywon, Lucas Majure, Ilaria Cobb, Amulya Bhaskara, et al. New World Cactaceae Plants Harbor Diverse Geminiviruses. Viruses, 2021, 13 (4), pp.694. 10.3390/v13040694. hal-03230622

HAL Id: hal-03230622

https://hal.inrae.fr/hal-03230622

Submitted on 25 May 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# New World Cactaceae Plants Harbor Diverse Geminiviruses

Rafaela S. Fontenele<sup>1,2 ID</sup>, Andrew M. Salywon<sup>3</sup>, Lucas C. Majure<sup>4 ID</sup>, Ilaria N. Cobb<sup>1,5</sup>, Amulya Bhaskara<sup>1,6</sup>, Jesús A. Avalos-Calleros<sup>7</sup>, Gerardo R. Argüello-Astorga<sup>7</sup>, Kara Schmidlin<sup>1,2</sup>, Anthony Khalifeh<sup>1,2</sup>, Kendal Smith<sup>1,2</sup>, Joshua Schreck<sup>1,2</sup>, Michael C. Lund<sup>1,2 ID</sup>, Matias Köhler<sup>8</sup>, Martin F. Wojciechowski<sup>2</sup>, Wendy C. Hodgson<sup>3 ID</sup>, Raul Puente-Martinez<sup>3 ID</sup>, Koenraad Van Doorslaer<sup>9</sup>, Safaa Kumari<sup>10 ID</sup>, Kehinde A. Oyeniran<sup>11 ID</sup>, Christian Vernière<sup>12,13</sup>, Denis Filloux<sup>12,13 ID</sup>, Philippe Roumagnac<sup>12,13</sup>, Pierre Lefeuvre<sup>14 ID</sup>, Simone G. Ribeiro<sup>15</sup>, Simona P. Kraberger<sup>1</sup>, Darren P. Martin<sup>11 ID</sup> and Arvind Varsani<sup>1,2,16,17,*</sup>

<figure><img src="image_2.png" /></figure>

check for
updates

**Citation:** Fontenele, R.S.; Salywon, A.M.; Majure, L.C.; Cobb, I.N.; Bhaskara, A.; Avalos-Calleros, J.A.; Argüello-Astorga, G.R.; Schmidlin, K.; Khalifeh, A.; Smith, K.; et al. New World Cactaceae Plants Harbor Diverse Geminiviruses. *Viruses* **2021**, *13*, 694. https://doi.org/10.3390/v13040694

Academic Editor: Jan Frederik Kreuze

Received: 15 March 2021

Accepted: 13 April 2021

Published: 16 April 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2021 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** The family Cactaceae comprises a diverse group of typically succulent plants that are native to the American continent but have been introduced to nearly all other continents, predominantly for ornamental purposes. Despite their economic, cultural, and ecological importance, very little research has been conducted on the viral community that infects them. We previously identified a highly divergent geminivirus that is the first known to infect cacti. Recent research efforts in non-cultivated and asymptomatic plants have shown that the diversity of this viral family has been under-sampled. As a consequence, little is known about the effects and interactions of geminiviruses in many plants, such as cacti. With the objective to expand knowledge on the diversity of geminiviruses infecting cacti, we used previously acquired high-throughput sequencing results to search for viral sequences using BLASTx against a viral RefSeq protein database. We identified two additional sequences with similarity to geminiviruses, for which we designed abutting primers and recovered full-length

<footnote n="1">¹ The Biodesign Center for Fundamental and Applied Microbiomics, Arizona State University, Tempe, AZ 85287, USA; rafasfontenele@asu.edu (R.S.F.); ilaria.cobb@gmail.com (I.N.C.); amulyabhaskara@gmail.com (A.B.); kara.schmidlin@asu.edu (K.S.); akhalif5@asu.edu (A.K.); krsmit39@asu.edu (K.S.); jrschrec@asu.edu (J.S.); mclund2@asu.edu (M.C.L.); simona.kraberger@asu.edu (S.P.K.)</footnote>

<footnote n="2">² School of Life Sciences, Arizona State University, Tempe, AZ 85287, USA; Martin.Wojciechowski@asu.edu</footnote>

<footnote n="3">³ Desert Botanical Garden, Phoenix, AZ 85008, USA; asalywon@dbg.org (A.M.S.); whodgson@dbg.org (W.C.H.); rpuente@dbg.org (R.P.-M.)</footnote>

<footnote n="4">⁴ Florida Museum of Natural History, University of Florida, Gainesville, FL 32611, USA; lmajure@floridamuseum.ufl.edu</footnote>

<footnote n="5">⁵ The University of British Columbia, Vancouver, BC V6T 1Z4, Canada</footnote>

<footnote n="6">⁶ School of Behavioral and Brain Sciences, University of Texas at Dallas, 800 W Campbell Rd, Richardson, TX 75080, USA</footnote>

<footnote n="7">⁷ División de Biología Molecular, Instituto Potosino de Investigación Científica y Tecnológica, A.C., Camino a la Presa de San José 2055, Lomas 4ta Secc, San Luis Potosí, San Luis Potosí 78216, Mexico; jesus.avalos@ipicyt.edu.mx (J.A.A.-C.); grarguel@ipicyt.edu.mx (G.R.A.-A.)</footnote>

<footnote n="8">⁸ Programa de Pós-Graduação em Botânica, Universidade Federal do Rio Grande do Sul, Porto Alegre-RS 90040-060, Brazil; matias.k@ufrgs.br</footnote>

<footnote n="9">⁹ UA Cancer Center, Department of Immunobiology, School of Animal and Comparative Biomedical Sciences, University of Arizona, Tucson, AZ 85721, USA; vandoorslaer@email.arizona.edu</footnote>

<footnote n="10">¹⁰ International Center for Agricultural Research in the Dry Areas (ICARDA), Terbol Station, Beqa'a, Zahle, Lebanon; s.kumari@cgiar.org</footnote>

<footnote n="11">¹¹ Department of Integrative Biomedical Sciences, Computational Biology Division, Institute of Infectious Diseases and Molecular Medicine, University of Cape Town, Observatory, Cape Town 7925, South Africa; oyenirankehinde@gmail.com (K.A.O.); darrenpatrickmartin@gmail.com (D.P.M.)</footnote>

<footnote n="12">¹² CIRAD, UMR PHIM, 34090 Montpellier, France; christian.verniere@cirad.fr (C.V.); denis.filloux@cirad.fr (D.F.); philippe.roumagnac@cirad.fr (P.R.)</footnote>

<footnote n="13">¹³ Plant Health Institute, University Montpellier, INRAE, CIRAD, 34090 Montpellier, France</footnote>

<footnote n="14">¹⁴ CIRAD, UMR PVBMT, F-97410 St. Pierre, France; pierre.lefeuvre@cirad.fr</footnote>

<footnote n="15">¹⁵ Embrapa Recursos Genéticos e Biotecnologia, Brasília 70770-917, Brazil; simone.ribeiro@embrapa.br</footnote>

<footnote n="16">¹⁶ Center for Evolution and Medicine, Arizona State University, Tempe, AZ 85287, USA</footnote>

<footnote n="17">¹⁷ Structural Biology Research Unit, Department of Clinical Laboratory Sciences, University of Cape Town, Cape Town 7925, South Africa</footnote>

<footnote>* Correspondence: Arvind varsani@asu.edu</footnote>
genomes. From 42 cacti and five scale insects, we derived 42 complete genome sequences of a novel geminivirus species that we have tentatively named Opuntia virus 2 (OpV2) and 32 genomes of an Opuntia-infecting becurtovirus (which is a new strain of the spinach curly top Arizona virus species). Interspecies recombination analysis of the OpV2 group revealed several recombinant regions, in some cases spanning half of the genome. Phylogenetic analysis demonstrated that OpV2 is a novel geminivirus more closely related to viruses of the genus *Curtovirus*, which was further supported by the detection of three recombination events between curtoviruses and OpV2. Both OpV2 and Opuntia becurtoviruses were identified in mixed infections, which also included the previously characterized Opuntia virus 1. Viral quantification of the co-infected cactus plants compared with single infections did not show any clear trend in viral dynamics that might be associated with the mixed infections. Using experimental *Rhizobium*-mediated inoculations, we found that the initial accumulation of OpV2 is facilitated by co-infection with OpV1. This study shows that the diversity of geminiviruses that infect cacti is under-sampled and that cacti harbor diverse geminiviruses. The detection of the Opuntia becurtoviruses suggests spill-over events between viruses of cultivated species and native vegetation. The threat this poses to cacti needs to be further investigated.

**Keywords:** Cactaceae; *Geminiviridae*; *Becurtovirus*; recombination; diversity

# 1. Introduction

Plants in the family Cactaceae have adapted to a wide-range of environmental conditions, which is reflected in a remarkable diversity of growth forms [1,2]. They are characterized by having adapted to diverse edaphically arid and semi-arid environments and, with the exception of the species *Rhipsalis baccifera* (Sol.) Stearn, which naturally occurs in Madagascar and South Africa [3,4], are native to the Americas. One of the regions with the highest diversity and abundance of cacti is North American deserts, such as the Sonoran, Chihuahuan and Mojave deserts [5]. However, several species have been introduced to all continents except Antarctica [3], especially as ornamental plants. Cacti are important for the ecosystems in which they exist by providing shade, water and sources of food for a variety of animals. The cactus genus *Opuntia* (commonly referred to as prickly pears and nopales) is also economically and culturally relevant for certain communities since their fruits, and sometimes stem segments, are also used as food sources [6]. In addition, cacti are farmed for the breeding of cochineal insects (genus *Dactylopius*) that produce the carminic acid dye used as a colorant for food, textiles, and in the pharmaceutical industry [7].

The known diversity of viruses associated with cacti has mostly been restricted to RNA viruses belonging to the single-stranded RNA virus families *Alphaflexiviridae*, *Betaflexiviridae*, *Puribunyaviridae*, *Tombusviridae* and *Virgaviridae* [8–14], but through recent research efforts, DNA viruses that either infect or are associated with cacti have been identified [15–18].

The known plant-infecting DNA viruses are classified into three families: *Geminiviridae* [19], *Nanoviridae* [20] and *Caulimoviridae* [21]. Two caulimoviruses have been described as infecting the cactus *Epiphyllum* [17,18] and our group has identified a divergent geminivirus named Opuntia virus 1 (OpV1) found to infect several cactus species in the USA and Mexico [15]. These recent findings suggest that the diversity of DNA viruses that infect cacti have been seriously under-sampled and therefore highlight the importance of future research efforts in this area.

Members of the family Geminiviridae [19] are plant-infecting viruses that contain a circular single-stranded DNA genome encapsidated in twinned semi-icosahedral particles [22,23]. Currently, the family Geminiviridae is divided into nine genera: *Becurtovirus*, *Begomovirus*, *Capulavirus*, *Curtovirus*, *Eragrovirus*, *Grablovirus*, *Mastrevirus*, *Topocuvirus*, and *Turncurtovirus* [19,24]. There are also several highly divergent unclassified geminivirus species (n = 14) [15,25–35]. The International Committee on Taxonomy of Viruses (ICTV) is reviewing a proposal to accommodate some of these divergent species and create five new
genera named *Citlodavirus*, *Maldovirus*, *Mulcrilevirus*, *Opunvirus* (to accommodate OpV1), and *Topilevirus* [36].

Geminiviruses are well known as causative agents of viral diseases in various crop plants and are associated with high yield losses [37,38]. However, many of the divergent members of the family (including most of the unclassified species) have been identified in non-cultivated native plants. In some instances, plants showed no observable visual symptoms, or only very mild symptoms caused by the virus [15,27,33,39–41]. The native and non-cultivated plants within which these divergent geminiviruses occur can act as reservoirs for viral species with the potential to emerge as crop pathogens [42,43].

Another area of concern is the far less well studied possibility of geminiviruses and other plant-infecting viruses spilling-over from agricultural crops to native vegetation. Studying such spillovers requires the study of viruses within an ecological framework, which some studies are beginning to address [44–46], and these are reviewed in [47,48]. In the case of plants such as cacti, which are often vegetatively propagated, the movement of infected plant material can mediate long-distance viral spreading [49]. This is especially likely to occur when the viruses concerned have no, or only very mild, associated symptoms. Although long-lived and/or vegetatively propagated plants are expected to select for viral variants with low levels of virulence, such plants can also be more susceptible to multiple or mixed infections, which can potentially select for increased virulence [50]. Indeed, geminiviruses that infect long-lived perennial plants are commonly found in mixed infections [51], which provides such viruses with more opportunities to recombine with members of other virus species.

Following the initial identification of OpV1, we decided to broaden the search for other DNA viruses infecting cacti in order to understand more about their viral community structure. Upon further investigation using high-throughput sequencing (HTS) data we identified fragments of sequences with similarity to other geminiviruses. Here we describe and characterize a novel divergent group of geminiviruses that we name Opuntia virus 2 (OpV2; 42 genomes), and a group of viruses (32 genomes) that fall within the genus *Becurtovirus* (sharing 82% identity with the member of *Spinach curly top Arizona virus* species). Although most of these viruses were found within single infections, in a few cases they were within detectably mixed infections, in some cases with OpV1 [15]. Infectivity assays were conducted with OpV2, the Opuntia becurtovirus, and OpV1 to evaluate viral load during co-infections, as well as to compare viral load between laboratory-infected *Nicotiana benthamiana* plants and naturally infected cactus samples.

## 2. Materials and Methods

### 2.1. *Sample Collecting and Processing*

Tissue samples were obtained from 577 plants in the Cactoideae and Opuntioideae clades [2]. These samples were collected from 19 countries: Argentina (n = 14), Bolivia (n = 8), Brazil (n = 8), Cuba (n = 1), Curaçao (n = 1), Dominican Republic (n = 2), France (n = 20), Haiti (n = 2), Lebanon (n = 1), Morocco (n = 1), Mexico (n = 31), Nigeria (n = 50), Paraguay (n = 3), Réunion (19), Spain (n = 6), Tunisia (n = 10), Uruguay (n = 5), the United States (n = 394), and Venezuela (n = 1) [15]. Samples were collected using either a 3-mm biopsy punch or scalpel blades and total DNA was extracted and processed as described previously [15]. In cases where plants were infested with cochineal insects, a cohort of about 5–10 insects were collected and used for total DNA extraction, as described previously [15].

### 2.2. *Geminivirus Genome Identification and Recovery*

Using the previously acquired high-throughput data from cacti samples, the contigs were mined for those with similarities to geminiviruses [15]. Two contigs with detectable homology to known geminiviruses were identified. Pairs of abutting primers were designed (OpV2_F 5'-CAT GTA TTT CAT CAT TTA CAA AAA GCA GAC TTA-3' and OpV2_R 5'-ATT TAG ATA TGG AGC AGA TTT GTT CCT CTT TT-3'; Bec_F 5'-TTG ATT TCG TTA GGC AAC CTA TTG AAT TCT-3' and Bec_R 5'-AGA GTG GGC AGA ACA TAA
TAT TTA TTT CGT-3') to recover potentially full-length virus genomes from plant and insect samples. The primers were used to amplify the geminivirus genomes using KAPA HiFi HotStart DNA polymerase (KAPA Biosystems, Wilmington, MA, USA), following the thermal cycling protocol: initial denaturation at 95 °C for 3 min, followed by 25 cycles at 98 °C for 20 s, at 60 °C for 15 s, at 72 °C for 3 min, and a final elongation at 72 °C for 3 min and a final renaturation at 4 °C for 10 min. Amplicons with a size between 2.5 and 3.5 kb (the size range of geminiviruses genomes) were resolved in 0.7% agarose gels, excised, purified, and cloned in the pJET1.2 cloning vector (Thermo Fisher Scientific, Waltham, MA, USA). Cloned amplicons were Sanger sequenced by means of primer walking at Macrogen Inc. (South Korea). Genome assemblies and annotations were performed using Geneious 11.1.5 (Biomatters Ltd., Auckland, New Zealand).

### 2.3. *Infectivity Assays*

One isolate of each virus was used for the infectivity assays. The OpV1 isolate infectious clone construction has been previously described [15]. The OpV2 isolateDBG_56 (GenBank accession MT840871) was used for the infectious clone, and this genome was recovered from an *Opuntia basilaris* sample from the Desert Botanical Garden (Phoenix, AZ, USA) collected in 2018. The Opuntia-derived becurtovirus isolate S18_40 (GenBank accession MT840851) that was recovered from an *Opuntia aciculata* sample growing in an urban area in Tempe (Arizona, USA) in 2018 was used to construct the infectious clone. Specific primers were designed to amplify two copies of each Opuntia-derived geminivirus that were cloned in tandem to the binary vector pJL-89 [52], excluding the 35S promoter region of the vector, using the Gibson assembly [53] (New England Biolabs, Ipswich, MA, USA), as previously described by Ferro et al. [54]. Each clone was transformed into competent *Escherichia coli* XL1 Blue cells, and to confirm that the ligation occurred correctly, clones were analyzed by digestion with EcoRV (OpV2) and NdeI/SalI (becurtovirus). The clone containing the two tandemly cloned copies of the Opuntia-derived geminiviruses were then used to transform *Rhizobium radiobacter* (synonymous species name for *Agrobacterium tumefaciens*) GV3101.

Initial infection assays were performed in *N. benthamiana* plants that were inoculated with OpV2 and the Opuntia becurtovirus as single infections. Furthermore, to evaluate mixed infection dynamics, *N. benthamiana* plants were inoculated with both OpV1/OpV2, OpV1/Opuntia becurtovirus, or OpV2/Opuntia becurtovirus to mirror mixed infections found in cactus plants. In all *Rhizobium*-inoculations, *R. radiobacter* was grown for 20 h in Luria broth with kanamycin (50 g/mL) and rifampicin (50 g/mL). The culture was then centrifuged for 10 min at 4600 rpm to pellet the cells before resuspension in MES buffer (10 mM MES hydrate and 10 mM MgSO4•7H2O) with acetosyringone 150 µM to an OD of 1.0. In the mixed infection, equal volumes of both viruses (OD 1.0) were mixed together prior to inoculation.

### 2.4. *Southern Blot Analysis*

Southern blots were used to detect the replicative forms of OpV2 during infection. Total nucleic acid was extracted from inoculated *N. benthamiana* plants as described above. Five µg total DNA from each plant, along with a positive control (5 ng of OpV2 PCR amplified genome), and one negative control of an uninoculated *N. benthamiana* plant, were resolved in a 1% agarose gel. The nucleic acid from the gel was transferred to a positively charged Hybond-N+ nylon membrane (GE Healthcare, Chicago, IL, USA) and UV-crosslinked using the default setting of the UV Stratalinker 1800 (Stratagene, La Jolla, CA, USA). The membrane was hybridized with a digoxigenin-labeled probe of the full-length OpV2 genome. The hybridization and probe preparation were carried out using a DIG High Prime DNA Labeling and Detection Starter Kit I (Roche, Basel, Switzerland), as directed by the manufacturer.
## 2.5. Viral Particle Purification and Transmission Electron Microscopy

Approximately 100 g of *N. benthamiana* leaf material 45 days post-inoculation with the OpV2 infectious clone (GenBank accession MT840871) was homogenized in 100 mL of extraction buffer (1 × PBS, 10 mg/mL sodium ascorbate, 2 mM PMSF, 1 mM EDTA). The homogenate was filtered through two layers of miracloth and two layers of cheesecloth, and subsequently centrifuged three times for 30 min at 14,800 rpm until clear. The clear homogenate was centrifuged for 4 h at 32,000 rpm using a Beckman 32 Ti rotor (Beckman, USA) onto a 10% sucrose cushion, and the pellet was resuspended in 1 mL of 1× PBS. Ten microliters of the resuspended pellet was absorbed onto carbon-coated copper grids for 2 min, washed three times, and stained with 2% uranyl acetate. The grids were viewed using a Philips CM 200 transmission electron microscope (Philips, Amsterdam, The Netherlands).

## 2.6. Viral Load Quantification

Absolute viral quantification was performed by means of real-time quantitative PCR (qPCR). Reactions were performed with either 40 ng of total *N. benthamiana* DNA or 10 ng of infected cactus plant DNA in a 10-µL qPCR reaction with 2x iTaq Universal SYBR Green Supermix (Bio-rad, USA), and 500 nM of each primer pair (Supplementary Table S1). Standard curves (StC) for each virus were prepared using a tenfold serial dilution of the plasmid containing the viral full-length genome diluted in 2.5 ng/µL of total DNA extracted from uninoculated *N. benthamiana* plants. The dilutions ranged from 10⁸ to 10 copies of viral DNA (genomic units) per microliter. All qPCR experiments contained a negative control and a non-template control, and were performed in triplicate. Reactions were performed in a Bio-rad CFX96 Real-time PCR System with the following conditions: 95 °C for 3 min, followed by 40 cycles of 95 °C for 30 s and 60 °C (varied by primer, Supplementary Table S1) for 30 s. A final melting curve analysis was performed ranging from 65 °C to 95 °C, with 0.5 °C increments every 5 s.

## 2.7. Pairwise Identity and Phylogenetic Analyses

Genome-wide pairwise identity comparison of the 32 cactus-derived becurtoviruses and the 42 OpV2 sequences were performed using SDT v.1.2 [55].

The full-length nucleotide sequences of representatives from each genotype (OpV2 and becurtovirus), together with those encoded by representatives of each geminivirus genera, including currently unclassified geminiviruses (n = 49), were aligned using MAFFT v.7 [56]. The resulting alignment was used to infer a neighbor-joining phylogenetic tree using the Jukes-Cantor nucleotide substitution model with 1000 bootstrap replicates to test branch support. Branches with <60% bootstrap support were collapsed using TreeGraph2 [57] and the tree was midpoint rooted.

In addition, the full-length nucleotide sequences of the OpV2 group (n = 42) and becurtovirus (n = 32) were each aligned using MUSCLE [58]. Each alignment was used to infer a neighbor-joining tree using the Jukes-Cantor substitution model with 1000 bootstrap replicates to test branch support. Branches with <60% bootstrap branch support were collapsed using TreeGraph2 [57]. The 42 and 32 genomes sequences, respectively, of the group OpV2 and cactus-derived becurtovirus, with fragments derived through interspecies recombination removed (see below), were used to infer a maximum-likelihood (ML) phylogenetic tree using IQ-TREE [59]. The model of nucleotide substitution used for the OpV2 tree was GTR+F+I+G4 and for the cactus-derived becurtovirus was TIM2+F+I+G4, as determined by ModelFinder [60] with 1000 bootstrap replicates to test for branch support.

Datasets were assembled for the inferred replication-associated protein (Rep) and capsid protein (CP) amino acid sequences from representatives of the OpV2 and becurtovirus genotype groups, along with representative sequences from the geminivirus genera and unclassified geminiviruses (n = 49). The Rep and CP amino acid sequence datasets were aligned using MAFFT v.7 [56], and the alignments were used to infer ML phylogenetic trees with models rtRev+G+I for the CP data and rtRev+G+I+F for the Rep data, as deter-
mined to be best-fitting by ProtTest [61] and with the approximate maximum likelihood
ratio test (aLRT) of branch support. Branches with <80% aLRT support were collapsed
using TreeGraph2 [57], and each ML tree was rooted with inferred Rep or CP amino acid
sequences from two members of the viral family Genomoviridae [62].

## 2.8. Recombination Analyses

The 42 OpV2 genomes and 32 becurtovirus genome sequences were aligned using MUSCLE [58], and the resulting alignments were used for intra-species recombination analysis using RDP v5.5 [63] with default settings using the detection methods RDP [64], GENECONV [65], BOOTSCAN [66], MAXCHI [67], CHIMERA [68], SISCAN [69], and 3SEQ [70]. Only recombination events that were detected by more than three methods with a *p*-value <0.05 were accepted.

Based on the ML phylogeny of the encoded Rep and CP amino acid sequences, we also decided to detect recombination events at an inter-generic level using representative sequences of the *Becurtovirus* and *Curtovirus* genera, together with one representative genome from OpV2, and the cactus-derived becurtovirus sequences were aligned using MUSCLE [58]. The nucleotide sequences on either side of the Rep and CP coding regions were removed from the alignment (these regions are non-homologous between the groups and can therefore not be aligned in any meaningful way), and this was then used for recombination analysis using RDP v5.5 [63] with the same standards used for the inter-species analysis.

## 2.9. Capsid Protein Cluster Analysis

The CP amino acid sequences of all geminiviruses available in GenBank were extracted and a representative dataset was generated through cluster analysis, using CD-HIT [71] with a 90% identity threshold. One representative from each cluster, together with representatives of the 11 genotypes of OpV2, 15 genotypes of OpV1 and one representative from the Opuntia becurtoviruses were used to generate a sequence similarity network using the enzyme function initiative-enzyme similarity tool (EFI-EST) [72] with a similarity score of 65 and an E-value threshold of $1 \times 10^{-5}$. The network was visualized in Cytoscape v3.8.2 [73] with the organic layout.

# 3. Results and Discussion

As a result of the use of HTS technologies, there has been a steady increase in the num-
ber of divergent geminiviruses that have been discovered. These newly discovered viruses
have further illuminated the vast breadth of geminivirus diversity and have revealed the
complex evolutionary history of this family. The combination of HTS with the rolling circle
amplification technique has helped to facilitate the identification of geminiviruses without
prior knowledge of the circulating viral population. Additionally, research efforts, which
have been expanded to include uncultivated native plant species instead of exclusively
crop plants displaying visible disease, have contributed tremendously to revealing the
true evolutionary and ecological contexts of the known geminivirus crop pathogens. As a
consequence, in the past five years, four new genera within the family Geminiviridae have
been recognized [24,74] to accommodate the novel divergent viruses identified through
HTS technologies, whereas an additional five new genera are currently under consideration
by the ICTV [36]. These findings demonstrate that the diversity of plant-infecting viruses
associated with crop diseases only represent a small proportion of the global plant-virus
diversity [75]. The study of geminivirus diversity will be of the utmost importance in the
context of agro-ecological interfaces [42,44,76,77], where spill-over between agricultural
and native vegetation is most likely to occur, and where these spill-overs are likely to
have the most significant ecological and/or agricultural impact [78,79]. Particularly in
areas where the environment has been anthropogenically modified, any sudden changes in
host population structure and viral transmission can lead to viral population changes and
possible disease emergence [80]. For long-lived and vegetatively propagated plants with
little industry, such as cacti, no sanitation measures exist and the management of infected plant material can be difficult, especially when pathogens have low virulence and plants do not display apparent symptoms.

### 3.1. *Geminiviruses Infecting Cacti*

In an initial analysis of cactus samples collected from 18 countries, we identified a novel geminivirus named OpV1 [15]. That initial finding led us to take a deeper look into the HTS data obtained from the cactus plants in this study, and consequently two other sequence fragments with similarity to geminiviruses were identified. By using two pairs of abutting primers, amplicons of ~3 kb were amplified from the cactus samples in question.

Sequence analyses of the putative full-length genomes revealed a divergent geminivirus, which we have tentatively named Opuntia virus 2 (OpV2), and a new becurtovirus strain of the spinach curly top Arizona virus species that we refer to as Opuntia becurtovirus (Table 1).

<table><caption>**Table 1.** Summary of the Opuntia virus 2 and Opuntia becurtoviruses identified in this study with isolate names, accession numbers, genotypes, host species, collection dates, and locations. The instances where the cochineal insect was associated with any cactus host species are highlighted.</caption><thead><tr><th>Virus</th><th>Host Species</th><th>ID</th><th>Accession Number</th><th>Genotype</th><th>Collection Year</th><th>Associated Scale Insect</th><th>Region of Collection</th></tr></thead><tbody><tr><td rowspan="24">Opuntia becurtovirus</td><td>Opuntia spp.</td><td>ASU_PP13</td><td>MT840840</td><td>2</td><td>2018</td><td>SI_47</td><td>Arizona, USA</td></tr><tr><td>Opuntia spp.</td><td>ASU_PP7</td><td>MT840839</td><td>2</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>Opuntia martiniana</td><td>DBG_38</td><td>MT840843</td><td>2</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>Cylindropuntia echinocarpa</td><td>DBG_80</td><td>MT840841</td><td>2</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>Cylindropuntia spinosior</td><td>DBG_86</td><td>MT840842</td><td>2</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>Opuntia phaeacantha</td><td>LCM_23</td><td>MT840844</td><td>2</td><td>2006</td><td></td><td>Texas, USA</td></tr><tr><td>Opuntia stenopetala</td><td>2014</td><td>MT840845</td><td>2</td><td>2015</td><td></td><td>Arizona, USA</td></tr><tr><td>Opuntia spp.</td><td>S18_100</td><td>MT840861</td><td>2</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>Opuntia spp.</td><td>S18_101</td><td>MT840862</td><td>2</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>Opuntia microdasys</td><td>S18_12</td><td>MT840846</td><td>2</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>Opuntia microdasys</td><td>S18_14</td><td>MT840847</td><td>2</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>Opuntia santa-rita</td><td>S18_17</td><td>MT840870</td><td>2</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>Opuntia basilaris</td><td>S18_24</td><td>MT840848</td><td>2</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>Opuntia engelmannii var. linguiformis</td><td>S18_30</td><td>MT840849</td><td>2</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>Opuntia spp.</td><td>S18_34</td><td>MT840850</td><td>2</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>Opuntia aciculata</td><td>S18_40</td><td>MT840851</td><td>2</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>Opuntia aciculata</td><td>S18_54</td><td>MT840863</td><td>2</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>Opuntia microdasys</td><td>S18_56</td><td>MT840864</td><td>2</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>Opuntia spp.</td><td>S18_59_1</td><td>MT840852</td><td>2</td><td>2018</td><td>SI_68</td><td>Arizona, USA</td></tr><tr><td>Opuntia spp.</td><td>S18_59_2</td><td>MT840865</td><td>2</td><td>2018</td><td>SI_68</td><td>Arizona, USA</td></tr><tr><td>Opuntia engelmannii var. lindheimeri</td><td>S18_69</td><td>MT840853</td><td>2</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>Opuntia santa-rita</td><td>S18_71_1</td><td>MT840854</td><td>2</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>Opuntia santa-rita</td><td>S18_71_2</td><td>MT840855</td><td>2</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>Opuntia phaeacantha</td><td>S18_77</td><td>MT840856</td><td>2</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>Opuntia microdasys</td><td>S18_84_1</td><td>MT840857</td><td>2</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>Opuntia microdasys</td><td>S18_84_2</td><td>MT840858</td><td>2</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>Opuntia engelmannii</td><td>S18_89_1</td><td>MT840859</td><td>2</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>Opuntia engelmannii</td><td>S18_89_2</td><td>MT840860</td><td>2</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>Scale insect</td><td>SI_47</td><td>MT840866</td><td>2</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>Scale insect</td><td>SI_68</td><td>MT840867</td><td>2</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>Opuntia sp.</td><td>TM3_2</td><td>MT840868</td><td>2</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>Opuntia engelmannii</td><td>UTH_RH6</td><td>MT840869</td><td>2</td><td>2018</td><td></td><td>Utah, USA</td></tr></tbody></table>
<table><caption>Table 1. *Cont.*</caption><thead><tr><th>Virus</th><th>Host Species</th><th>ID</th><th>Accession Number</th><th>Genotype</th><th>Collection Year</th><th>Associated Scale Insect</th><th>Region of Collection</th></tr></thead><tbody><tr><td rowspan="20">Opuntia virus 2</td><td>Opuntia basilaris</td><td>DBG_56</td><td>MT840871</td><td>8</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>Opuntia basilaris</td><td>DBG_57</td><td>MT840872</td><td>8</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>Opuntia santa-rita</td><td>DBG_61</td><td>MT840873</td><td>8</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>Opuntia santa-rita</td><td>DBG_62</td><td>MT840874</td><td>8</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>Opuntia santa-rita</td><td>DBG_63</td><td>MT840875</td><td>8</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>Opuntia microdasys</td><td>S18_12_1</td><td>MT840888</td><td>3</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>S18_12_2</td><td>MT840889</td><td>3</td><td>2018</td><td></td><td></td><td>Arizona, USA</td></tr><tr><td>Opuntia microdasys</td><td>S18_14</td><td>MT840890</td><td>3</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td rowspan="2">Opuntia phaeacantha</td><td>S18_25_1</td><td>MT840891</td><td>5</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>S18_25_2</td><td>MT840892</td><td>5</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td rowspan="3">Opuntia phaeacantha</td><td>S18_26_1</td><td>MT840893</td><td>5</td><td>2018</td><td>SI_63</td><td>Arizona, USA</td></tr><tr><td>S18_26_2</td><td>MT840894</td><td>5</td><td>2018</td><td>SI_63</td><td>Arizona, USA</td></tr><tr><td>S18_26_3</td><td>MT840895</td><td>5</td><td>2018</td><td>SI_63</td><td>Arizona, USA</td></tr><tr><td>Opuntia phaeacantha</td><td>S18_27</td><td>MT840896</td><td>3</td><td>2018</td><td>SI_64</td><td>Arizona, USA</td></tr><tr><td rowspan="3">Opuntia engelmannii</td><td>S18_3_1</td><td>MT840876</td><td>6</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>S18_3_2</td><td>MT840877</td><td>6</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>S18_3_3</td><td>MT840878</td><td>6</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td rowspan="3">Opuntia engelmannii</td><td>S18_4_1</td><td>MT840879</td><td>7</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>S18_4_2</td><td>MT840880</td><td>7</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>S18_4_3</td><td>MT840881</td><td>7</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td rowspan="3">Opuntia aciculata</td><td>S18_40_1</td><td>MT840897</td><td>9</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>S18_40_2</td><td>MT840898</td><td>4</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>S18_40_3</td><td>MT840899</td><td>8</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td rowspan="2">Opuntia microdasys</td><td>S18_41_1</td><td>MT840900</td><td>10</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>S18_41_2</td><td>MT840901</td><td>9</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td rowspan="2">Opuntia basilaris</td><td>S18_5_1</td><td>MT840882</td><td>7</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>S18_5_2</td><td>MT840883</td><td>7</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td rowspan="4">Opuntia santa-rita</td><td>S18_75_1</td><td>MT840902</td><td>1</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>S18_75_2</td><td>MT840903</td><td>1</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>S18_75_3</td><td>MT840904</td><td>1</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>S18_75_4</td><td>MT840905</td><td>11</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td rowspan="4">Opuntia santa-rita</td><td>S18_8_1</td><td>MT840884</td><td>7</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>S18_8_2</td><td>MT840885</td><td>7</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>S18_8_3</td><td>MT840886</td><td>7</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>S18_8_4</td><td>MT840887</td><td>7</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td rowspan="2">Scale insect</td><td>SI_63_1</td><td>MT840906</td><td>5</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>SI_63_2</td><td>MT840907</td><td>5</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>Scale insect</td><td>SI_64</td><td>MT840908</td><td>3</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>Scale insect</td><td>SI_70</td><td>MT840909</td><td>4</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td rowspan="2">Opuntia sp.</td><td>TM_3_1</td><td>MT840910</td><td>2</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>TM_3_2</td><td>MT840911</td><td>2</td><td>2018</td><td></td><td>Arizona, USA</td></tr><tr><td>Opuntia santa-rita</td><td>UTH_RH4</td><td>MT840912</td><td>8</td><td>2018</td><td></td><td>Utah, USA</td></tr></tbody></table>

## 3.2. *Opuntia* Virus 2

A total of 42 OpV2 genome sequences were recovered from 16 cactus plants and three cochineal insects collected in 2018 from the USA. With the exception of one plant sample from Utah, the remaining samples were collected in the state of Arizona (Table 1). Additionally, two of the groups of cochineal insects from which OpV2 genomes were recovered were collected from OpV2-positive plants (Table 1). None of the plants presented
any apparent symptoms.

### 3.2.1. Genome Organization and Diversity

The OpV2 sequences recovered ranged from 3194–3247 nt in size and have a genomic organization similar to other geminiviruses, containing at least six recognizable open-reading frames (ORFs) encoding proteins >84 amino acids. Interestingly, the OpV2 genome size is about 10% larger than most other known geminiviruses; however, other divergent geminiviruses recently identified through HTS approaches also present genomes > 3000 nt [26,28,32]. Based on sequence similarity with geminivirus-expressed proteins, the OpV2 ORF's in the complementary strand likely encode a replication-associated protein (Rep), a symptom determinant protein (C4), and a hypothetical protein. Additionally, the OpV2 virion strand likely encodes a capsid protein (CP), a regulatory protein (V2), and a possible movement protein (V3) (Figure 1). The 42 isolated genomes from OpV2 share 90.3%–100% identity to each other and were grouped into 11 genotypes, based on a 95% nucleotide identity cut-off (Table 1; Supplementary Data S1). Interestingly, two of the 16 plants analyzed were infected with more than two genotypes (Table 1).

All the OpV2 genomes contained the conserved geminivirus virion-strand origin of replication nonanucleotide sequence "TAATATTAC", which is located within a sequence that is capable of forming a hairpin-loop structure with the conserved nonanucleotide in its loop. In the long intergenic region, we also identified iterative potential replication-associated protein recognition sequences, commonly referred to as "iterons," and a potential TATA box (Figure 1). Among the OpV2 genotypes two "iterons" were identified, one adjacent to the TATA box and the other immediately upstream of the Rep gene. It is interesting to note that in eight of the eleven OpV2 genotypes, the "iteron" adjacent to the TATA box occurs as two tandem TATA repeats. Sequences from genotypes 3 and 6 present two different "iteron" sequences.

The OpV2 representative genotypes all share <74.1% nucleotide pairwise identity with other geminiviruses (Supplementary Data S1). Their encoded Rep and CP proteins respectively share <76.8% and <85.4% amino acid sequence identity with other geminiviruses (Supplementary Data S2). Hence, the OpV2 sequences are diverse compared to the currently classified geminiviruses and likely represent a new genus. The Rep amino acid sequences encoded by OpV2 all have a rolling-circle replication, geminivirus Rep sequence (GRS), a helicase SF3, and Walker motifs [81], which are other important conserved features shared with geminiviruses and which are essential for the rolling circle replication functions of Rep (Supplementary Table S2).

Phylogenetic analysis of representative OpV2 sequences, together with representatives of other geminiviruses (Figure 1), show that OpV2 clusters evolutionarily together with species in the genus *Curtovirus*. Not surprisingly, the highest genome-wide pairwise identities of the OpV2 sequences are with the members of the curtovirus species, such as that of beet curly top virus (AF379637) (Supplementary Data S2). In contrast, ML phylogenetic analysis of representative amino acid sequences of the encoded OpV2 Rep and CP proteins, together with those of other geminiviruses, demonstrate a more complex scenario. The OpV2 encoded CPs cluster with those of the curtoviruses, becurtoviruses, and two presently unclassified divergent geminiviruses, *Limeum africanum*-associated virus [31] and parsley yellow leaf curl virus [35] (Figure 2). The OpV2 encoded Rep sequences group together with those of two curtoviruses, horseradish curly top virus and spinach severe curly top virus (Figure 2).
<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /><figcaption>**Figure 1.** (A) Neighbor-joining phylogenetic tree of the genome sequences of the eleven genotypes of OpV2 and the Opuntia becurtovirus, together with representative sequences from various genera in the family *Geminiviridae*. Branches with <60% bootstrap support have been collapsed and the phylogenetic tree is midpoint-rooted. The genomic organization of OpV2 and the Opuntia becurtovirus are illustrated on the right side of the phylogenetic tree near their respective groups. (B) Organization and orientation of the replication-associated interactive sequences "iterons" in the intergenic region of both eleven genotypes of OpV2 and the Opuntia becurtovirus. The arrows indicate the orientation of the iteron sequences relative to the nonanucleotide and lower-case letters indicate that the nucleotide is variable among sequences in that genotype. Some genotypes presented more than one type of "iteron" sequence.</figcaption></figure>
<figure><img src="image_6.png" /><figcaption>**Figure 2.** Maximum-likelihood phylogenetic tree of the Rep and CP amino acid sequences of the eleven representative genotypes of OpV2 and Opuntia becurtovirus, and representative sequences from various genera in the family *Geminiviridae*. Branches with <0.8 aLRT support have been collapsed and the trees are rooted with sequences of genomoviruses [62].</figcaption></figure>

### 3.2.2. Identification of Recombination in OpV2

Recombination plays an important role in the evolution of geminiviruses [82–84]. We analyzed the 42 OpV2 sequences for evidence of recombination using RDP5 [63] and identified 21 well-supported recombination events (Figure 3, Table 2). Interestingly, genotypes 5 and 7 had no detectable recombinant regions. In genotype 8, only one of the seven sequences that comprise this group had detectable evidence of recombination. Recombinant region sizes ranged from 142 to 1790 nt. One sequence from genotype 3 had the largest recombination transferred genome fragment, which corresponds to nearly half of the genome, spanning the *cp* gene, small intergenic region, and a gene encoding a hypothetical protein (Figure 3).
<figure><img src="image_7.png" /><figcaption>**Figure 3.** (A) Neighbor-joining phylogenetic tree of the 42 OpV2 sequences with genotypes denoted by symbols. (B) Maximum-likelihood phylogenetic tree of the 42 OpV2 sequences with recombination regions removed. The eleven genotypes of OpV2 are represented by symbols. Graphical representation of each genome representing the recombination event with the breakpoint location within the genome.</figcaption></figure>

<table><caption>**Table 2.** Summary of the 21 recombination events detected in the OpV2 sequences by RDP5 v5.5. The methods used to detect recombination were RDP (R), GENCONV (G), BOOTSCAN (B), MAXCHI (M), CHIMERA (C), SISCAN (S) and 3SEQ (T). The method with the highest *p*-value for each recombination event is bolded. Sites where the actual breakpoint is undetermined are marked with *. (T) denotes traces of recombination signals and (P) denotes partial evidence. Recombinant sequences marked with ^ indicate that the recombinant sequence may have been misidentified (one of the identified parents might be the recombinant). Please refer to the table for accession # of the genotypes.</caption><thead><tr><th>Recombination EVENT</th><th>Region</th><th>Recombinant Sequence(s)</th><th>Minor Parental Sequence(s)</th><th>Major Parental Sequence(s)</th><th>Detection Methods</th><th>p-Value</th></tr></thead><tbody><tr><td>1</td><td>2242–3220</td><td>Genotype 10</td><td>Genotype 9</td><td>Genotype 4</td><td>GBMCST</td><td>3.15 × 10<sup>-42</sup></td></tr><tr><td>2~</td><td>2294–3203</td><td>Genotype 4</td><td>Genotype 3</td><td>Genotype 8</td><td>RGBMCST</td><td>1.91 × 10<sup>-34</sup></td></tr><tr><td>3</td><td>2318–3210</td><td>Genotype 8</td><td>Genotype 8<br>Genotype 7</td><td>Genotype 7<br>Genotype 8</td><td>RGBMCT</td><td>1.12 × 10<sup>-17</sup></td></tr><tr><td>4</td><td>531–2241</td><td>Genotype 3</td><td>Genotype 3</td><td>Genotype 5</td><td>RGBMCST</td><td>7.72 × 10<sup>-22</sup></td></tr><tr><td>5</td><td>2325–3200</td><td>^Genotype 4</td><td>Genotype 3<br>Genotype 8</td><td>Genotype 8<br>Genotype 3</td><td>RBMCT</td><td>1.52 × 10<sup>-15</sup></td></tr><tr><td>6</td><td>1051–2144</td><td>^Genotype 3</td><td>Unknown(Genotype 7)</td><td>Genotype 3</td><td>RBMCT</td><td>3.23 × 10<sup>-8</sup></td></tr><tr><td>7</td><td>3014–3186</td><td>^Genotype 3<br>Genotype 4</td><td>Genotype 5</td><td>Genotype 3</td><td>GMCST</td><td>7.31 × 10<sup>-14</sup></td></tr><tr><td>8</td><td>2147–173</td><td>^Genotype 6</td><td>Genotype 7</td><td>Genotype 6</td><td>RGBMCST</td><td>7.09 × 10<sup>-18</sup></td></tr><tr><td>9</td><td>2585–2950</td><td>Genotype 6</td><td>Genotype 6</td><td>Genotype 3</td><td>RGBMCST</td><td>4.51 × 10<sup>-13</sup></td></tr><tr><td>10</td><td>3009–3178</td><td>^Genotype 6</td><td>Genotype 1</td><td>Unknown<br>(Genotype 7)</td><td>RGBMCS</td><td>9.77 × 10<sup>-16</sup></td></tr></tbody></table>
<table><caption>Table 2. *Cont.*</caption><thead><tr><th>Recombination EVENT</th><th>Region</th><th>Recombinant Sequence(s)</th><th>Minor Parental Sequence(s)</th><th>Major Parental Sequence(s)</th><th>Detection Methods</th><th>p-Value</th></tr></thead><tbody><tr><td>11</td><td>1466–3037</td><td>Genotype 3<br>Genotype 6[T]<br>Genotype 3[P]</td><td>Genotype 5</td><td>Genotype 3</td><td>MCT</td><td>3.69 × 10<sup>-6</sup></td></tr><tr><td>12</td><td>2611–3142</td><td>Genotype 11</td><td>Genotype 9</td><td>Genotype 1</td><td>RBMCST</td><td>7.37 × 10<sup>-15</sup></td></tr><tr><td>13</td><td>3204–192</td><td>Genotype 3</td><td>Genotype 1</td><td>Genotype 5</td><td>RBT</td><td>1.93 × 10<sup>-3</sup></td></tr><tr><td>14</td><td>1740–2322</td><td>^Genotype 7</td><td>Unknown (Genotype 5)</td><td>Genotype 7</td><td>RGMCST</td><td>1.78 × 10<sup>-11</sup></td></tr><tr><td>15</td><td>2493–3161</td><td>^Genotype 9<br>Genotype 10</td><td>Genotype 8</td><td>Unknown<br>(Genotype 7)</td><td>GBMCS</td><td>2.69 × 10<sup>-18</sup></td></tr><tr><td>16</td><td>3169 *–1027</td><td>^Genotype 3</td><td>Genotype 3</td><td>Genotype 6</td><td>RGBMCST</td><td>2.95 × 10<sup>-9</sup></td></tr><tr><td>17</td><td>2008–2991</td><td>^Genotype 1</td><td>Genotype 5</td><td>Genotype 8</td><td>RGBMCS</td><td>5.19 × 10<sup>-15</sup></td></tr><tr><td>18</td><td>1034–1650 *</td><td>^Genotype 8</td><td>Genotype 6</td><td>Genotype 6</td><td>RGMCS</td><td>9.12 × 10<sup>-11</sup></td></tr><tr><td>19</td><td>1862–2995 *</td><td>^Genotype 6</td><td>Genotype 5</td><td>Genotype 8</td><td>RGMCST</td><td>2.68 × 10<sup>-12</sup></td></tr><tr><td>20</td><td>1907–2049</td><td>^Genotype 3</td><td>Genotype 4</td><td>Genotype 4<br>Genotype 10</td><td>RGMCS</td><td>1.11 × 10<sup>-6</sup></td></tr><tr><td>21</td><td>349–606</td><td>^Genotype 3</td><td>Genotype 2</td><td>Genotype 3</td><td>RBT</td><td>5.23 × 10<sup>-4</sup></td></tr></tbody></table>

The majority of the recombination events occur in the virion-strand gene coding regions. Additionally, the breakpoints for several recombination events are located in the long intergenic region where the origin of replication is located, a feature that has been extensively reported for other geminiviruses [15,43,84–86]. Collectively, the same recombination event seemed to occur in only one or up to three sequences and always within the same genotype, which indicates that the detected events likely occurred more recently than when the last common ancestor of the 42 OpV2 sequences existed. This demonstrates that coinfection between these viruses are common.

### 3.2.3. Infection Assays

To assess the infectivity of OpV2, an infectious clone with two tandem copies of OpV2 DBG_56 (GenBank accession MT840871) was constructed. The infectivity assays in *N. benthamiana* plants showed that based on three independent assays performed with 10 inoculated plants in each, OpV2 had a 40–50% rate of systemic infection. OpV2 infectivity was further confirmed with Southern blot analysis of DNA from *N. benthamiana*-infected plants that revealed evidence of circular covalently closed DNA, super coiled dsDNA and single-strand DNA. No replicative forms of OpV2 were observed for the uninoculated control plants (Supplementary Figure S1A).

Transmission electron microscopy analysis of OpV2 particles recovered from inoculated *N. benthamiana* plant material revealed the presence of twinned icosahedral particles (Supplementary Figure S1B). This suggests that despite the fact that the OpV2 genome is approximately 10% larger than the most other characterized geminiviruses, the ~3.2 kb genome is still packaged into geminate particles.

### 3.3. *Opuntia Becurtovirus*

The *Becurtovirus* genus is currently composed of three species, beet curly top Iran virus, spinach curly top Arizona virus, and *Exomis microphylla* latent virus. Their genomic organization consists of a capsid protein gene (*cp*), a regulatory protein gene (*v2*) and a movement protein gene (*v3*) on the virion-sense strand and on the complementary-sense strand, a replication-associated protein gene (*rep*), potentially from alternatively spliced transcripts. The nonanucleotide “TAATATTAC” sequence that is found in the virion-strand origin of replication is highly conserved among the geminiviruses; however, viruses in the genus *Becurtovirus* have a different nonanucleotide sequence, i.e., “TAAGATTCC”.
Thus far, becurtoviruses have been found to infect eudicotyledonous plants of the species *Beta vulgaris* (beet), *Beta vulgaris* subsp. *maritima* (beet), *Capsicum frutescens* (pepper), *Exomis microphylla*, *Phaseolus vulgaris* (bean), *Spinacia oleracea* (spinach), *Solanum lycopersicum* (tomato), *Solanum melongena* (eggplant) and *Vigna unguiculata* (black-eyed pea) [31,87–94].

Beet curly top Iran virus is transmitted by the leafhopper *Circulifer haematoceps*, a species commonly found in Iran [92]. The only becurtovirus previously described in the Americas is spinach curly top Arizona virus, a virus found in symptomatic spinach plants from Arizona, USA, in 2009 [93]. Since then, to our knowledge this becurtovirus has not been reported, likely due to under-sampling. The identification of the Opuntia becurtovirus could suggest a spill-over event in agro-ecological interface areas.

### 3.3.1. Diversity of Opuntia Becurtoviruses

The becurtoviruses identified in the cactus samples share ~82% genome-wide pairwise identity with spinach curly top Arizona virus (Supplementary Data S2). According to the current species and strain demarcation for the genus *Becurtovirus* [74], the cactus-derived becurtovirus is a new strain from the species spinach curly top Arizona virus, which we have named Opuntia becurtovirus.

A total of 32 genome sequences of the Opuntia becurtovirus were obtained from 26 cactus plants and two cochineal insects collected in 2018 in the USA in the states of Arizona (n = 26), Texas (n = 1) and Utah (n = 1) (Table 1). The cactus-derived becurtovirus genomes range from 2899 to 2982 nt and have a similar genome organization as compared to other becurtoviruses (Figure 1). Based on the pairwise sequence identity results, the 32 genomes of cactus-infecting becurtovirus range from 94.8%–100% sequence identity among themselves (Supplementary Data S3) and based on a 95% identity cut-off they are all members of the same genotype. SCTAV is therefore tentatively assigned to genotype 1 and the Opuntia becurtovirus to genotype 2.

Phylogenetic analysis of the genome sequences from the Opuntia becurtoviruses with representatives of the geminivirus family shows, as expected, that they group with other becurtoviruses, being more closely related to SCTAV (Figure 1). The same is observed in the ML phylogenetic tree of the Opuntia becurtovirus encoded Rep and CP amino acid sequences (Figure 2). The rolling circle amplification GRS, helicase SF3 and Walker motifs [81] were also present in the predicted Rep amino acid sequences of the Opuntia becurtoviruses and all were very similar to their counterparts in the SCTAV isolates (Supplementary Table S2).

### 3.3.2. Identification of Recombination in Becurtoviruses

In the interspecies recombination analysis of the Opuntia becurtoviruses, three well supported events were detected (Figure 4; Table 3). The recombinant regions ranged from 522 to 1640 nt. Again, the largest of the detected events (event 1 in Table 1) implied the recombination transfer of nearly half of the genome, spanning the intergenic region and the complementary sense-encoded proteins. Another event (number 2 in Table 1) was evident in two sequences, both derived from the same cactus sample (Table 1).

### 3.3.3. Infectivity Assays

To evaluate the infectivity of the Opuntia becurtovirus, an infectious clone was constructed and used to inoculate *N. benthamiana* plants. In three independent experiments with 10 plants each, only two plants were detectably infected through conventional PCR at 45 days post-inoculation. We were not able to detect the replicative forms of the Opuntia becurtovirus after Southern blot analysis of the infected plants (data not shown). However, qPCR quantification of both *N. benthamiana* plants inoculated with Opuntia becurtovirus showed their presence at very low levels with approximately 100 viral copies in 40 ng of plant material or 2 to 3 copies/ng of plant material (Supplementary Figure S2). It is likely that Opuntia becurtovirus does not replicate well in *N. benthamiana* plants or a longer period of viral infection might be required to establish infection.
<figure><img src="image_8.png" /><figcaption>**Figure 4.** (A) Neighbor-joining phylogenetic tree of the 32 Opuntia becurtovirus sequences. Branches with <60% bootstrap branch support have been collapsed and the tree is midpoint-rooted. (B) Maximum-likelihood phylogenetic tree of the 32 Opuntia becurtovirus sequences with recombination regions removed. Graphical representation of each genome representing the recombination event with the breakpoint location within the genome.</figcaption></figure>

<table><caption>**Table 3.** Summary of the three recombination events detected in the Opuntia becurtovirus by RDP5 v.5.5. The methods used to detect recombination were RDP (R) GENCONV (G), BOOTSCAN (B), MAXCHI (M), CHIMERA (C), SISCAN (S), and 3SEQ (T). The method with the highest *p*-value for each recombination event is bolded. Sites where the actual breakpoint is undetermined are marked with *. Recombinant sequences marked with ^ indicate that the recombinant sequence may have been misidentified (one of the identified parents might be the recombinant).</caption><thead><tr><th>Recombination Event</th><th>Region</th><th>Recombinant Sequence(s)</th><th>Minor Parental Sequence(s)</th><th>Major Parental Sequence(s)</th><th>Detection Methods</th><th>p-Value</th></tr></thead><tbody><tr><td>1</td><td>2790–1165</td><td>^MT840851</td><td>MT840856</td><td>Unknown</td><td>RBMCST</td><td>1.29 × 10<sup>-7</sup></td></tr><tr><td rowspan="2">2</td><td rowspan="2">1546–2590 *</td><td>^MT840860</td><td>MT840870</td><td rowspan="2">MT840848</td><td rowspan="2">MCST</td><td rowspan="2">2.86 × 10<sup>-11</sup></td></tr><tr><td>MT840859</td><td>MT840856</td></tr><tr><td>3~</td><td>1525–2036</td><td>^MT840839</td><td>MT840841</td><td>MT840868</td><td>RBCS</td><td>4.86 × 10<sup>-7</sup></td></tr></tbody></table>

### 3.4. Inter-Genus Recombination

Viruses in several geminivirus genera show clear evidence of inter-genus recombination [39,77,93,95,96]. The geminivirus genomic organization, which involves bi-directionally transcribed genes and the geminivirus mode of rolling circle replication, are both factors that seem to influence the rates and patterns of recombination in these viruses [43,84,97].

Recombination is clearly evident in the incongruence observed in the phylogenetic trees of the encoded Rep and CP amino acid sequences (Figure 2). Since the OpV2 encoded Rep and CP protein sequences cluster together with those of viruses in the genera *Becurtovirus* and *Curtovirus*, we undertook a recombination analysis to detect any events of recombination between members of these genera. We took representatives of each
genus, including one sequence from OpV2 and the Opuntia becurtovirus to create a dataset. The dataset sequences were aligned and this alignment was trimmed to retain only the nucleotide sequences of the Rep and CP coding regions (the only parts of the genomes that were clearly homologous).

In the recombination analysis, seven well supported events were identified, all of them in the sequences of curtoviruses (Table 4, Figure 5). Three events were identified in the CP sequence of curtoviruses, two of which were interspecies events. In one event (event 3 in Figure 5), the recombinant region appears to be derived from curtoviruses and becurtoviruses. Recombination events 1 and 3 span almost the entire CP nucleotide sequence. In the Rep coding region, four events were detected. One displays contributions from both curtoviruses and becurtoviruses and the other three are recombinants between curtoviruses and OpV2 (Figure 5).

<table><caption>**Table 4.** Summary of the seven recombination events from a dataset of representative sequences of OpV2 and Opuntia becurtovirus, with representatives from the genera *Becurtovirus* and *Curtovirus*, as detected by RDP5 5.5 [63]. The methods used to detect recombination were RDP (R) GENCONV (G), BOOTSCAN (B), MAXCHI (M), CHIMERA (C), SISCAN (S), and 3SEQ (T). The method with the highest *p*-value for each recombination event is bolded. Sites where the actual breakpoint is undetermined are marked with *. Recombinant sequences marked with ^ indicate that the recombinant sequence may have been misidentified (one of the identified parents might be the recombinant).</caption><thead><tr><th>Recombination Event</th><th>Region</th><th>Recombinant Sequence(s)</th><th>Minor Parental Sequence(s)</th><th>Major Parental Sequence(s)</th><th>Detection Methods</th><th>p-Value</th></tr></thead><tbody><tr><td>1</td><td>726–1474</td><td>*KX529650 curtovirus<br>EU921828 curtovirus</td><td>Unknown<br>(AF379637 curtovirus)</td><td>AY134867 curtovirus<br>JN817383 curtovirus</td><td>RGBMCST</td><td>8.79 × 10<sup>-33</sup></td></tr><tr><td>2~</td><td>2719–3418 *</td><td>*AF379637 curtovirus</td><td>MT840871 OpV2</td><td>AY134867 curtovirus<br>JN817383 curtovirus<br>U02311 curtovirus</td><td>RGBMCS</td><td>2.56 × 10<sup>-24</sup></td></tr><tr><td>3</td><td>729 *–1498 *</td><td>*GU734126 curtovirus</td><td>HQ443515 becurtovirus</td><td>U49907 curtovirus</td><td>RMCST</td><td>1.22 × 10<sup>-14</sup></td></tr><tr><td>4</td><td>2978–3374</td><td>*U02311 curtovirus</td><td>MT840871 OpV2</td><td>JN817383 curtovirus</td><td>RGBMCST</td><td>5.13 × 10<sup>-46</sup></td></tr><tr><td>5</td><td>2843 *–3375</td><td>*U49907 curtovirus</td><td>Unknown<br>(GU734126 curtovirus)<br>Unknown<br>(MT840871 OpV2)</td><td>MT840871 OpV2<br>GU734126 curtovirus</td><td>RMCT</td><td>1.08 × 10<sup>-11</sup></td></tr><tr><td>6</td><td>727 *–928</td><td>EU921828 curtovirus</td><td>JN817383 curtovirus<br>AY134867 curtovirus</td><td>KX529650 curtovirus</td><td>GBMCST</td><td>4.57 × 10<sup>-12</sup></td></tr><tr><td>7</td><td>2073–2502</td><td>*JN817383 curtovirus</td><td>U02311 curtovirus</td><td>AY134867<br>Becurtovirus</td><td>RGBMST</td><td>4.11 × 10<sup>-12</sup></td></tr></tbody></table>

<figure><img src="image_9.png" /><figcaption>**Figure 5.** (A) Neighbor-joining phylogenetic tree of the OpV2 and becurtovirus representative sequences, together with representatives from the genera *Becurtovirus* and *Curtovirus*. Branches with <60% bootstrap branch support have been collapsed and the tree is midpoint-rooted. (B) Maximum-likelihood (ML) phylogenetic tree of the OpV2 and becurtovirus representative sequences, together with representatives from the genera *Becurtovirus* and *Curtovirus* with recombination regions removed. Branches with <60% bootstrap support have been collapsed and the tree is midpoint-rooted. Graphical representation of each genome, representing the recombination event with the breakpoint location within the OpV2 genome for reference, which is color coded according to the genera of the major and minor parent sequences of each event. In bold are the two sequences that represent OpV2 and Opuntia becurtovirus.</figcaption></figure>
### 3.5. Mixed Infections of Geminiviruses in Cacti

Geminivirus mixed infections are frequently found in nature. There are several records of mixed infection from species in the same genera, especially species in the genus *Begomovirus* [85,98–100]. In some cases, co-infecting viruses can have synergistic interactions, leading to increases in pathogenicity [101–105]. Mixed infections also allow for the exchange of genetic material through recombination [65,86], which can lead to the emergence of new variants, species, and genera.

Naturally occurring mixed infections were detected in a few cactus samples in our study (Table 1). Two plants were infected with both OpV1 and OpV2, whereas Opuntia becurtovirus was also found in co-infections with OpV1 (n = 6) and OpV2 (n = 2). We used qPCR to quantify the viral loads of these co-infected plants compared to single infections in other cactus samples. Unsurprisingly, Opuntia becurtovirus was not detected by qPCR in either the co-infected or single infected *N. benthamiana* plants, which means the virus in these plants was at very low titers, since we used 10 ng of DNA from a total plant DNA extract for the quantification reactions. Nevertheless, both OpV1 and OpV2 were detected by means of qPCR, and cactus samples had varying viral loads (genomic units/ng) (Figure 6). There is no clear evidence in the cactus samples that co-infection seems to be favoring one virus over the other (Figure 6). It should be noted although the cacti analyzed are different species (Table 1), which might help explain the variable qPCR results. However, there were also variations in viral load within the same cactus plants at different locations within the plants. In the *Opuntia basilaris* plants co-infected with OpV1 and OpV2, two pads from the same plant were collected (DBG_56 andDBG_57) and showed different co-infection dynamics. In one pad (DBG_56), OpV1 had a higher viral load than OpV2, whereas in the second pad (DBG_57), OpV2 had a higher viral load than OpV1 (Figure 6). It is unknown if geminiviruses can establish a systemic infection in cacti, and further research is certainly needed to address this.

A

<figure><img src="image_10.png" /></figure>

Cactus plants with single infection

<figure><img src="image_11.png" /></figure>

<figure><img src="image_12.png" /></figure>

B

<figure><img src="image_13.png" /></figure>

Cactus plants with mixed infection

<figure><img src="image_14.png" /></figure>

<figure><img src="image_15.png" /></figure>

**Figure 6.** Bar graph of the viral loads (genomic units/ng of total DNA) determined by real-time quantitative PCR for the Opuntia virus 1 (OpV1), Opuntia virus 2 (OpV2), and Opuntia becurtovirus in a subset of cactus samples collected in this study. The graphs show (A) the group of single-infected cacti and (B) the cactus samples that presented mixed infection with standard deviation for the triplicate reactions. The Opuntia becurtovirus was the only virus that did not present detectable levels of genomic units.
To further investigate the potential dynamics of those co-infections, infectivity assays were performed with OpV1 and OpV2 in *N. benthamiana*. Of 10 co-inoculated *N. benthamiana* plants, four were found to be co-infected, one plant was found to be infected with only OpV2, and five plants only with OpV1. The four plants with mixed infections and four of the plants with single infections were used for viral load quantification using qPCR. All the samples were analyzed at two time points, 17- and 45-days post-inoculation (dpi). OpV1 infected plants showed little difference in viral load between the time points (Figure 7). However, OpV2 could only be detected at 45 dpi but with similar viral loads to plants with OpV1 infections (Figure 7). The mixed-infection plants showed some interesting dynamics. At 17 dpi OpV2 was detectable, and its viral load significantly increased by 45 dpi (Figure 7). On the other hand, OpV1 maintained similar viral loads over the two time points. These results suggests that OpV2 favors mixed infections with OpV1, but in this case it did not seem to affect the viral load of OpV1.

In attempted co-infections with OpV1 and Opuntia becurtovirus and OpV2 and Opuntia becurtovirus, the fact that Opuntia becurtovirus was not detected was not surprising based on the initial single infection experiments. Therefore, we could not draw any conclusions on coinfections involving this species (data not shown).

### 3.6. Cluster Analysis of Geminivirus Capsid Protein

Geminivirus CP amino acid sequences are possibly co-diverging with their specific insect vectors, as highlighted in the phylogenetic analysis presented in [48]. Hence, we undertook a sequence similarity network analysis of the CP amino acid sequences of all geminiviruses (with a >90% identity cut off) and further associated the identified clusters with the known insect vectors (Figure 8). Sequence similarity network analyses are useful for clustering groups of similar protein sequences in large datasets relatively quickly (see [48,106–108] for examples). Unlike OpV1, which forms its own distinct cluster [15], indicating that it is likely transmitted by an insect vector that has not yet been associated with geminiviruses, OpV2 CPs cluster with those of curtoviruses, becurtoviruses (including Opuntia becurtovirus), and two unclassified viruses (parsley yellow leaf curl virus and *Limeum africanum*-associated virus). Becurtoviruses and curtoviruses are known to be transmitted by leafhoppers (*Circulifer* sp.), and based on the cluster analysis, it is highly likely that OpV2 is also transmitted by this insect vector. Controlled insect transmission experiments will be necessary to properly test this hypothesis.

The identification of OpV2 and Opuntia becurtoviruses expands the diversity of known geminiviruses that are associated with plants in the family Cactaceae. OpV2 is a novel geminivirus that will likely be assigned to a new genus. Nevertheless, the molecular characterization of OpV2 shows that it shares several similarities with other geminiviruses including its genomic organization, the presence of conserved motifs at the origin of replication, intergenic region "iterons", and characteristic Rep motifs. OpV2 has been identified in several species of *Opuntia*. Based on phylogenetic analysis, it is clear that OpV2 is very closely related to viruses in the genus *Curtovirus*. This is further supported by the presence of several recombination events in the Rep coding region between curtoviruses and OpV2 sequences.

It is plausible that the presence of Opuntia becurtovirus in cacti could be an example of spill-over from agricultural areas into native vegetation, given that spinach curly top Arizona virus was initially identified in spinach fields in Arizona, where it caused severe symptoms [93]. However, Opuntia becurtovirus was identified in three different states in the USA, suggesting it might have been circulating in the natural ecosystem for a longer time and that the spinach infection might have been a spillover or emergence event that originated in native uncultivated plants. More viral surveys with an ecological focus are needed to address this question. Nonetheless, our viral load analysis suggests that the viral load of Opuntia becurtovirus in cactus samples and experimentally inoculated *N. benthamiana* plants is low compared to those of OpV1 and OpV2.
<figure><img src="image_16.png" /><figcaption>**Figure 7.** Bar graph of the viral loads (genomic units/ng total DNA) observed through real-time quantitative PCR for the Opuntia virus 1 (OpV1) and Opuntia virus 2 (OpV2) in *N. benthamiana* plants at 17 and 45 days post-inoculation. (A) *N. benthamiana* inoculated with OpV1; (B) *N. benthamiana* inoculated with OpV2 and (C) *N. benthamiana* plants with mixed infection of OpV1 and OpV2 with standard deviation for the triplicate reactions.</figcaption></figure>
<figure><img src="image_17.png" /><figcaption>**Figure 8.** Sequence similarity network analysis of the CP amino acid sequences of geminiviruses present in GenBank (dataset was created with an amino acid identity cut-off of 90%), together with representatives of the 11 genotypes from OpV2, 15 genotypes from OpV1, and one Opuntia becurtovirus. Each dot represents a sequence and is colored based on the genus or group. The genera that have known insect vectors are highlighted in a light grey box with the insect vector name displayed at the top. Clusters or singletons marked with a brown halo have no known insect vector associated with them. ApGV, apple geminivirus; ACSV, *Axonopus compressus* streak virus; CaCDaV, Camellia chlorotic dwarf-associated virus; CCDaV, citrus chlorotic dwarf-associated virus; DfasMV, dragonfly-associated mastrevirus; CBCSV, common bean curly stunt virus; ECSV, *Eragrostis curvula* streak virus; EMSV, *Eragrostis minor* streak virus; GraGV, grapevine geminivirus; JmaV, *Juncus maritimus*-associated virus; LaaV, *Limeum africanum*-associated virus; MCaV, mulberry crinkle-associated virus; MiSV, *Miscanthus* streak virus; MSMV, maize streak Reunion virus; PYLCV, parsley yellow leaf curl virus; PMLCV-1, paper mulberry leaf curl virus 1; PMLCV-1, paper mulberry leaf curl virus 2; PCMoV, passion fruit chlorotic mottle virus; PgaV, *Polygala garcinii*-associated virus; RLV1, rice latent virus 1; RLV2, rice latent virus 2; SMAv, switchgrass mosaic-associated virus; SpSMV1, sweet potato symptomless mastrevirus 1; SSMV1, *Sporobolus striate* mosaic virus 1; SSMV2, *Sporobolus striate* mosaic virus 2; SStV-A, sugarcane striate virus A; SStV-D, sugarcane striate virus D; SWSV, sugarcane white streak virus; TaGV, tomato-associated geminivirus; ToALCV, tomato apical leaf curl virus.</figcaption></figure>

Mixed infections are commonly reported for geminiviruses, and here we identified a few cactus samples infected with either OpV1 and OpV2 (n = 2) or Opuntia becurtovirus and either OpV2 (n = 2) or OpV1 (n = 6). Viral load quantification of the co-infected plants in comparison with cacti infected with only one of those viruses did not show any clear trend that would indicate that mixed infections favor any of the viruses. Further analysis of the OpV1 and OpV2 mixed infections in *N. benthamiana* plants showed that OpV2 is favored in mixed infections with OpV1, in that it displayed higher titers earlier in infections in the presence of OpV1 than it did in its absence.

Until very recently, cactus plants were mainly known to be infected by RNA viruses. Based on this study and our previous one [15], there are now three geminiviruses known to infect cacti. OpV1, OpV2 and Opuntia becurtovirus were found to infect cacti in the USA and OpV1 was found in cacti from Mexico. The fact that two novel geminiviruses and a distinct becurtovirus have been identified in an initial survey of cactus plants indicates that these plants may be hosts to other viruses. Further research efforts are needed to broaden our knowledge on viral diversity in cactus and the ways that these cactus viruses interact. The identification of Opuntia becurtovirus, which is most closely related to spinach curly top Arizona virus, may hint towards possible spill-overs between natural and agricultural areas and thus raises questions as to the extent to which viruses from cropping systems spill over into endemic uncultivated plants and vice-versa.
**Supplementary Materials:** The following are available online at https://www.mdpi.com/article/10.3390/v13040694/s1, Supplementary Table S1: Primer pairs used for real-time quantitative PCR analysis of the OpV1, OpV2 and Opuntia becurtovirus, including the sequence and melting temperature. Supplementary Table S2: Summary of the rolling circle, SF3, and Walker motifs identified in the Rep amino acid sequences encoded by the Opuntia Virus 2, Opuntia becurtovius, and spinach curly top Arizona virus sequences. Supplementary Figure S1: (A) Southern blot showing replicative forms of OpV2 infectious clone from *N. benthamiana* plants inoculated with OpV2 and Opuntia becurtovirus, including a non-inoculated control and positive control with an amplicon of linearized full-length OpV2 genome (5 ng). The gel images on the top show the DNA (5 µg) loaded for each sample. CCC—circular covalently closed; SC—super coiled dsDNA; SS—single-stranded. (B) Transmission electron micrographs (negatively stained with 2% uranyl acetate) of geminate particles (shown with white arrows) from OpV2-inoculated plants. Bottom image scale bar is 100 nm and the top zoomed section scale bar is 50 nm. Supplementary Figure S2: Bar graph represents the viral load quantification of two *N. benthamiana* plants infected with Opuntia becurtovirus at 45 days post-inoculation, acquired by real-time quantitative PCR. Standard deviation represents triplicate reactions. Supplementary Data S1: Pairwise sequence identity of the full-length genome sequences and encoded amino acid sequences of the Rep and CP from the 42 full-length OpV2 genomes. Supplementary Data S2: Pairwise sequence identity of the full-length genome sequences and encoded amino acid sequences of the Rep and CP of OpV2 genotypes and Opuntia becurtovirus, together with representatives of the geminiviruses genera and unclassified geminiviruses. Supplementary Data S3: Pairwise sequence identity of the full-length genome sequences and encoded amino acid sequences of the Rep and CP from the 32 cactus-derived becurtovirus genomes.

**Author Contributions:** Conceptualization, R.S.F. and A.V.; methodology, R.S.F. and A.V.; software, R.S.F., D.P.M., S.P.K. and A.V.; validation, R.S.F. and A.V.; formal analysis, R.S.F., J.A.A.-C., G.R.A.-A., D.P.M., S.P.K. and A.V.; investigation, R.S.F., A.M.S., L.C.M., I.N.C., A.B., J.A.A.-C., G.R.A.-A., K.S. (Kara Schmidlin), A.K., K.S. (Kendal Smith), J.S., M.C.L., M.K., M.F.W., W.C.H., R.P.-M., K.V.D., S.K., C.V., D.F., P.R., P.L., S.G.R., S.P.K., D.P.M., K.A.O. and A.V.; resources, R.S.F., A.M.S., L.C.M., I.N.C., A.B., K.S. (Kara Schmidlin), M.C.L., M.K., M.F.W., W.C.H., R.P.-M., K.V.D., S.K., C.V., D.F., P.R., P.L., S.G.R., S.P.K., D.P.M., K.A.O. and A.V.; data curation, R.S.F., S.P.K. and A.V.; writing—original draft preparation, R.S.F., D.P.M., S.P.K. and A.V.; writing—review and editing, R.S.F., A.M.S., L.C.M., I.N.C., A.B., J.A.A.-C., G.R.A.-A., K.S., A.K., K.S., J.S., M.C.L., M.K., M.F.W., W.C.H., R.P.-M., K.V.D., S.K., C.V., D.F., P.R., P.L., S.G.R., S.P.K., D.P.M., K.A.O. and A.V.; visualization, R.S.F. and A.V.; supervision, A.M.S. and A.V.; project administration, A.V.; funding acquisition, A.V. All authors have read and agreed to the published version of the manuscript.

**Funding:** The molecular work was supported by a start-up grant awarded to A.V. from Arizona State University, USA.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** Sequence determined as part of this study have been deposited in GenBank under accession #s MT840839—MT840912.

**Acknowledgments:** We thank The School of Life Sciences Electron Microscopy Lab at Arizona State University for the electron microscopy support.

**Conflicts of Interest:** The authors declare no conflict of interest.

## References

1. Hernandez-Hernandez, T.; Brown, J.W.; Schlumpberger, B.O.; Eguiarte, L.E.; Magallon, S. Beyond aridification: Multiple explanations for the elevated diversification of cacti in the New World Succulent Biome. *New Phytol.* **2014**, *202*, 1382–1397.
[CrossRef] [PubMed]

2. Hernandez-Hernandez, T.; Hernandez, H.M.; De-Nova, J.A.; Puente, R.; Eguiarte, L.E.; Magallon, S. Phylogenetic relationships and evolution of growth form in Cactaceae (Caryophyllales, Eudicotyledoneae). *Am. J. Bot.* **2011**, *98*, 44–61.
[CrossRef]

3. Anderson, E.F. *The Cactus Family*; Timber Press: Portland, OR, USA, 2001; p. 776.

4. Guerrero, P.C.; Majure, L.C.; Cornejo-Romero, A.; Hernandez-Hernandez, T. Phylogenetic Relationships and Evolutionary Trends in the Cactus Family. *J. Hered.* **2019**, *110*, 4–21.
[CrossRef] [PubMed]
5. Ortega-Baes, P.; Godínez-Alvarez, H. Global diversity and conservation priorities in the Cactaceae. *Biodivers. Conserv.* **2006**, *15*, 817–827. [CrossRef]

6. Nobel, P.S. *Cacti: Biology and Uses*; University of California Press: Los Angeles, CA, USA, 2002; p. 280.

7. Greenfield, A.B. *A Perfect Red: Empire, Espionage, and the Quest for the Color of Desire*; Harper Collins: New York, NY, USA, 2009; p. 352.

8. Milbrath, G.M. Isolation and Characterization of a Virus from Saguaro Cactus. *Phytopathology* **1972**, *62*, 739. [CrossRef]

9. Weng, Z.; Xiong, Z. Genome organization and gene expression of saguaro cactus carmovirus. *J. Gen. Virol.* **1997**, *78*, 525–534. [CrossRef] [PubMed]

10. Koenig, R.; Pleij, C.W.; Loss, S.; Burgermeister, W.; Aust, H.; Schiemann, J. Molecular characterisation of potexviruses isolated from three different genera in the family Cactaceae. *Arch. Virol.* **2004**, *149*, 903–914. [CrossRef] [PubMed]

11. Min, B.E.; Song, Y.S.; Ryu, K.H. Complete sequence and genome structure of cactus mild mottle virus. *Arch. Virol.* **2009**, *154*, 1371–1374. [CrossRef]

12. Min, B.E.; Chung, B.N.; Kim, M.J.; Ha, J.H.; Lee, B.Y.; Ryu, K.H. Cactus mild mottle virus is a new cactus-infecting tobamovirus. *Arch. Virol.* **2006**, *151*, 13–21. [CrossRef] [PubMed]

13. Liou, M.R.; Chen, Y.R.; Liou, R.F. Complete nucleotide sequence and genome organization of a Cactus virus X strain from Hylocereus undatus (Cactaceae). *Arch. Virol.* **2004**, *149*, 1037–1043. [CrossRef]

14. Sanches, M.M.; Lamas, N.S.; Reis, M.B.; Arieta-Sosa, J.G.; Romano, E.; Melo, F.L.; Ribeiro, S.G. Genome Assembly of Schlumbergera Virus X Infecting Prickly Pear (Opuntia cochenillifera) in Brazil. *Genome Announc.* **2015**, *3*. [CrossRef]

15. Fontenele, R.S.; Salywon, A.M.; Majure, L.C.; Cobb, I.N.; Bhaskara, A.; Avalos-Calleros, J.A.; Arguello-Astorga, G.R.; Schmidlin, K.; Khalifeh, A.; Smith, K.; et al. A Novel Divergent Geminivirus Identified in Asymptomatic New World Cactaceae Plants. *Viruses* **2020**, *12*, 398. [CrossRef] [PubMed]

16. Fontenele, R.S.; Roumagnac, P.; Richet, C.; Kraberger, S.; Stainton, D.; Aleamotu'a, M.; Filloux, D.; Bernardo, P.; Harkins, G.W.; McCarthy, J.; et al. Diverse genomoviruses representing twenty-nine species identified associated with plants. *Arch. Virol.* **2020**, *165*, 2891–2901. [CrossRef] [PubMed]

17. Zheng, L.; Cao, M.; Wu, L.; Liu, H.; Chen, M.; Li, R. First identification and molecular characterization of a novel cavemovirus infecting Epiphyllum spp. *Arch. Virol.* **2020**, *165*, 2083–2086. [CrossRef] [PubMed]

18. Lan, P.; Tian, T.; Pu, J.; Rao, W.; Li, F.; Li, R. Characterization and detection of a new badnavirus infecting Epiphyllum spp. *Arch. Virol.* **2019**, *164*, 1837–1841. [CrossRef] [PubMed]

19. Zerbini, F.M.; Briddon, R.W.; Idris, A.; Martin, D.P.; Moriones, E.; Navas-Castillo, J.; Rivera-Bustamante, R.; Roumagnac, P.; Varsani, A.; Ictv Report, C. ICTV Virus Taxonomy Profile: Geminiviridae. *J. Gen. Virol.* **2017**, *98*, 131–133. [CrossRef] [PubMed]

20. Thomas, J.E.; Gronenborn, B.; Harding, R.M.; Mandal, B.; Grigoras, I.; Randles, J.W.; Sano, Y.; Timchenko, T.; Vetten, H.J.; Yeh, H.H.; et al. ICTV Virus Taxonomy Profile: Nanoviridae. *J. Gen. Virol.* **2021**, *102*. [CrossRef]

21. Teycheney, P.Y.; Geering, A.D.W.; Dasgupta, I.; Hull, R.; Kreuze, J.F.; Lockhart, B.; Muller, E.; Olszewski, N.; Pappu, H.; Pooggin, M.M.; et al. ICTV Virus Taxonomy Profile: Caulimoviridae. *J. Gen. Virol.* **2020**, *101*, 1025–1026. [CrossRef]

22. Zhang, W.; Olson, N.H.; Baker, T.S.; Faulkner, L.; Agbandje-McKenna, M.; Boulton, M.I.; Davies, J.W.; McKenna, R. Structure of the Maize streak virus geminate particle. *Virology* **2001**, *279*, 471–477. [CrossRef] [PubMed]

23. Hesketh, E.L.; Saunders, K.; Fisher, C.; Potze, J.; Stanley, J.; Lomonossoff, G.P.; Ranson, N.A. The 3.3 A structure of a plant geminivirus using cryo-EM. *Nat. Commun.* **2018**, *9*, 2369. [CrossRef]

24. Varsani, A.; Roumagnac, P.; Fuchs, M.; Navas-Castillo, J.; Moriones, E.; Idris, A.; Briddon, R.W.; Rivera-Bustamante, R.; Murilo Zerbini, F.; Martin, D.P. Capulavirus and Grablovirus: Two new genera in the family Geminiviridae. *Arch. Virol.* **2017**, *162*, 1819–1831. [CrossRef]

25. Ma, Y.; Navarro, B.; Zhang, Z.; Lu, M.; Zhou, X.; Chi, S.; Di Serio, F.; Li, S. Identification and molecular characterization of a novel monopartite geminivirus associated with mulberry mosaic dwarf disease. *J. Gen. Virol.* **2015**, *96*, 2421–2434. [CrossRef] [PubMed]

26. Loconsole, G.; Saldarelli, P.; Doddapaneni, H.; Savino, V.; Martelli, G.P.; Saponari, M. Identification of a single-stranded DNA virus associated with citrus chlorotic dwarf disease, a new member in the family Geminiviridae. *Virology* **2012**, *432*, 162–172. [CrossRef] [PubMed]

27. Liang, P.; Navarro, B.; Zhang, Z.; Wang, H.; Lu, M.; Xiao, H.; Wu, Q.; Zhou, X.; Di Serio, F.; Li, S. Identification and characterization of a novel geminivirus with a monopartite genome infecting apple trees. *J. Gen. Virol.* **2015**, *96*, 2411–2420. [CrossRef]

28. Zhang, S.; Shen, P.; Li, M.; Tian, X.; Zhou, C.; Cao, M. Discovery of a novel geminivirus associated with camellia chlorotic dwarf disease. *Arch. Virol.* **2018**, *163*, 1709–1712. [CrossRef] [PubMed]

29. Al Rwahnih, M.; Alabi, O.J.; Westrick, N.M.; Golino, D.; Rowhani, A. Description of a Novel Monopartite Geminivirus and Its Defective Subviral Genome in Grapevine. *Phytopathology* **2017**, *107*, 240–251. [CrossRef] [PubMed]

30. Vaghi Medina, C.G.; Teppa, E.; Bornancini, V.A.; Flores, C.R.; Marino-Buslje, C.; Lopez Lambertini, P.M. Tomato Apical Leaf Curl Virus: A Novel, Monopartite Geminivirus Detected in Tomatoes in Argentina. *Front. Microbiol.* **2017**, *8*, 2665. [CrossRef] [PubMed]

31. Claverie, S.; Bernardo, P.; Kraberger, S.; Hartnady, P.; Lefeuvre, P.; Lett, J.M.; Galzi, S.; Filloux, D.; Harkins, G.W.; Varsani, A.; et al. From Spatial Metagenomics to Molecular Characterization of Plant Viruses: A Geminivirus Case Study. *Adv. Virus Res.* **2018**, *101*, 55–83. [CrossRef] [PubMed]
32. Fontenele, R.S.; Abreu, R.A.; Lamas, N.S.; Alves-Freitas, D.M.T.; Vidal, A.H.; Poppiel, R.R.; Melo, F.L.; Lacorte, C.; Martin, D.P.; Campos, M.A.; et al. Passion Fruit Chlorotic Mottle Virus: Molecular Characterization of a New Divergent Geminivirus in Brazil. *Viruses* **2018**, *10*, 169. [CrossRef] [PubMed]

33. Fontenele, R.S.; Lamas, N.S.; Lacorte, C.; Lacerda, A.L.M.; Varsani, A.; Ribeiro, S.G. A novel geminivirus identified in tomato and cloome plants sampled in Brazil. *Virus Res.* **2017**, *240*, 175–179. [CrossRef] [PubMed]

34. Zhang, R.; Wu, X.; Jiang, X.; Wu, X.; Luan, X.; Cheng, X. Molecular characterization of common bean curly stunt virus: A novel recombinant geminivirus in China. *Arch. Virol.* **2020**, *165*, 257–260. [CrossRef]

35. Hasanvand, V.; Heydanejad, J.; Massumi, H.; Kleinow, T.; Jeske, H.; Varsani, A. Isolation and characterization of a novel geminivirus from parsley. *Virus Res.* **2020**, *286*, 198056. [CrossRef] [PubMed]

36. Roumagnac, P.; Varsani, A.; Martin, D.P.; Lett, J.-M. 2020.008P.A.v1.Geminiviridae_5ng_11nsp. 2020. Available online: https://talk.ictvonline.org/files/proposals/taxonomy-proposals_plant1/m/plant04/10509 (accessed on 15 March 2021).

37. Rojas, M.R.; Macedo, M.A.; Maliano, M.R.; Soto-Aguilar, M.; Souza, J.O.; Briddon, R.W.; Kenyon, L.; Rivera Bustamante, R.F.; Zerbini, F.M.; Adkins, S.; et al. World Management of Geminiviruses. *Annu. Rev. Phytopathol.* **2018**, *56*, 637–677. [CrossRef]

38. Moffat, A.S. PLANT PATHOLOGY: Geminiviruses Emerge as Serious Crop Threat. *Science* **1999**, *286*, 1835. [CrossRef]

39. Bernardo, P.; Golden, M.; Akram, M.; Nadarajan, N.; Fernandez, E.; Granier, M.; Rebelo, A.G.; Peterschmitt, M.; Martin, D.P.; Roumagnac, P.; et al. Identification and characterisation of a highly divergent geminivirus: Evolutionary and taxonomic implications. *Virus Res.* **2013**, *177*, 35–45. [CrossRef] [PubMed]

40. Kraberger, S.; Geering, A.D.W.; Walters, M.; Martin, D.P.; Varsani, A. Novel mastreviruses identified in Australian wild rice. *Virus Res.* **2017**, *238*, 193–197. [CrossRef]

41. Susi, H.; Laine, A.L.; Filloux, D.; Kraberger, S.; Farkas, K.; Bernardo, P.; Frilander, M.J.; Martin, D.P.; Varsani, A.; Roumagnac, P. Genome sequences of a capulavirus infecting Plantago lanceolata in the Aland archipelago of Finland. *Arch. Virol.* **2017**, *162*, 2041–2045. [CrossRef]

42. Rocha, C.S.; Castillo-Urquiza, G.P.; Lima, A.T.; Silva, F.N.; Xavier, C.A.; Hora-Junior, B.T.; Beserra-Junior, J.E.; Malta, A.W.; Martin, D.P.; Varsani, A.; et al. Brazilian begomovirus populations are highly recombinant, rapidly evolving, and segregated based on geographical location. *J. Virol.* **2013**, *87*, 5784–5799. [CrossRef] [PubMed]

43. Varsani, A.; Shepherd, D.N.; Monjane, A.L.; Owor, B.E.; Erdmann, J.B.; Rybicki, E.P.; Peterschmitt, M.; Briddon, R.W.; Markham, P.G.; Oluwafemi, S.; et al. Recombination, decreased host specificity and increased mobility may have driven the emergence of maize streak virus as an agricultural pathogen. *J. Gen. Virol.* **2008**, *89*, 2063–2074. [CrossRef] [PubMed]

44. Bernardo, P.; Charles-Dominique, T.; Barakat, M.; Ortet, P.; Fernandez, E.; Filloux, D.; Hartnady, P.; Rebelo, T.A.; Cousins, S.R.; Mesleard, F.; et al. Geometagenomics illuminates the impact of agriculture on the distribution and prevalence of plant viruses at the ecosystem scale. *ISME J.* **2018**, *12*, 173–184. [CrossRef]

45. Sallinen, S.; Norberg, A.; Susi, H.; Laine, A.L. Intraspecific host variation plays a key role in virus community assembly. *Nat. Commun.* **2020**, *11*, 5610. [CrossRef] [PubMed]

46. Susi, H.; Laine, A.L. Agricultural land use disrupts biodiversity mediation of virus infections in wild plant populations. *New Phytol.* **2020**. [CrossRef]

47. McLeish, M.J.; Fraile, A.; Garcia-Arenal, F. Population Genomics of Plant Viruses: The Ecology and Evolution of Virus Emergence. *Phytopathology* **2021**, *111*, 32–39. [CrossRef]

48. Lefeuvre, P.; Martin, D.P.; Elena, S.F.; Shepherd, D.N.; Roumagnac, P.; Varsani, A. Evolution and ecology of plant viruses. *Nat. Rev. Microbiol.* **2019**, *17*, 632–644. [CrossRef] [PubMed]

49. Gaur, R.K.; Petrov, N.M.; Patil, B.L.; Stoyanova, M.I. *Plant Viruses: Evolution and Management*; Springer: Berlin/Heidelberg, Germany, 2016; p. 312.

50. Van Mölken, T.; Stuefer, J.F. Virulence in clonal plants: Conflicting selection pressures at work? *Evol. Ecol.* **2007**, *22*, 467–470. [CrossRef]

51. Garcia-Arenal, F.; Zerbini, F.M. Life on the Edge: Geminiviruses at the Interface between Crops and Wild Plant Hosts. *Annu. Rev. Virol.* **2019**, *6*, 411–433. [CrossRef] [PubMed]

52. Lindbo, J.A. TRBO: A high-efficiency tobacco mosaic virus RNA-based overexpression vector. *Plant. Physiol.* **2007**, *145*, 1232–1240. [CrossRef] [PubMed]

53. Gibson, D.G.; Young, L.; Chuang, R.Y.; Venter, J.C.; Hutchison, C.A., 3rd; Smith, H.O. Enzymatic assembly of DNA molecules up to several hundred kilobases. *Nat. Methods* **2009**, *6*, 343–345. [CrossRef] [PubMed]

54. Ferro, M.M.M.; Ramos-Sobrinho, R.; Xavier, C.A.D.; Zerbini, F.M.; Lima, G.S.A.; Nagata, T.; Assuncao, I.P. New approach for the construction of infectious clones of a circular DNA plant virus using Gibson Assembly. *J. Virol. Methods* **2019**, *263*, 20–23. [CrossRef] [PubMed]

55. Muhire, B.M.; Varsani, A.; Martin, D.P. SDT: A virus classification tool based on pairwise sequence alignment and identity calculation. *PLoS ONE* **2014**, *9*, e108277. [CrossRef] [PubMed]

56. Katoh, K.; Standley, D.M. MAFFT multiple sequence alignment software version 7: Improvements in performance and usability. *Mol. Biol. Evol.* **2013**, *30*, 772–780. [CrossRef] [PubMed]

57. Stover, B.C.; Muller, K.F. TreeGraph 2: Combining and visualizing evidence from different phylogenetic analyses. *BMC Bioinform.* **2010**, *11*, 7. [CrossRef] [PubMed]
58. Edgar, R.C. MUSCLE: Multiple sequence alignment with high accuracy and high throughput. *Nucleic Acids Res.* **2004**, *32*, 1792–1797. [CrossRef] [PubMed]

59. Nguyen, L.T.; Schmidt, H.A.; von Haeseler, A.; Minh, B.Q. IQ-TREE: A fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. *Mol. Biol. Evol.* **2015**, *32*, 268–274. [CrossRef] [PubMed]

60. Kalyaanamoorthy, S.; Minh, B.Q.; Wong, T.K.F.; von Haeseler, A.; Jermiin, L.S. ModelFinder: Fast model selection for accurate phylogenetic estimates. *Nat. Methods* **2017**, *14*, 587–589. [CrossRef]

61. Darriba, D.; Taboada, G.L.; Doallo, R.; Posada, D. ProtTest 3: Fast selection of best-fit models of protein evolution. *Bioinformatics* **2011**, *27*, 1164–1165. [CrossRef]

62. Varsani, A.; Krupovic, M. Sequence-based taxonomic framework for the classification of uncultured single-stranded DNA viruses of the family Genomoviridae. *Virus Evol.* **2017**, *3*, vew037. [CrossRef] [PubMed]

63. Martin, D.P.; Varsani, A.; Roumagnac, P.; Botha, G.; Maslamoney, S.; Schwab, T.; Kelz, Z.; Kumar, V.; Murrell, B. RDP5: A computer program for analysing recombination in, and removing signals of recombination from, nucleotide sequence datasets. *Virus Evol.* **2020**. [CrossRef]

64. Martin, D.; Rybicki, E. RDP: Detection of recombination amongst aligned sequences. *Bioinformatics* **2000**, *16*, 562–563. [CrossRef] [PubMed]

65. Padidam, M.; Sawyer, S.; Fauquet, C.M. Possible emergence of new geminiviruses by frequent recombination. *Virology* **1999**, *265*, 218–225. [CrossRef] [PubMed]

66. Martin, D.P.; Posada, D.; Crandall, K.A.; Williamson, C. A modified bootscan algorithm for automated identification of recombinant sequences and recombination breakpoints. *AIDS Res. Hum. Retrovir.* **2005**, *21*, 98–102. [CrossRef] [PubMed]

67. Smith, J.M. Analyzing the mosaic structure of genes. *J. Mol. Evol.* **1992**, *34*, 126–129. [CrossRef] [PubMed]

68. Posada, D.; Crandall, K.A. Evaluation of methods for detecting recombination from DNA sequences: Computer simulations. *Proc. Natl. Acad. Sci. USA* **2001**, *98*, 13757–13762. [CrossRef] [PubMed]

69. Gibbs, M.J.; Armstrong, J.S.; Gibbs, A.J. Sister-scanning: A Monte Carlo procedure for assessing signals in recombinant sequences. *Bioinformatics* **2000**, *16*, 573–582. [CrossRef] [PubMed]

70. Boni, M.F.; Posada, D.; Feldman, M.W. An exact nonparametric method for inferring mosaic structure in sequence triplets. *Genetics* **2007**, *176*, 1035–1047. [CrossRef] [PubMed]

71. Huang, Y.; Niu, B.; Gao, Y.; Fu, L.; Li, W. CD-HIT Suite: A web server for clustering and comparing biological sequences. *Bioinformatics* **2010**, *26*, 680–682. [CrossRef]

72. Gerlt, J.A.; Bouvier, J.T.; Davidson, D.B.; Imker, H.J.; Sadkhin, B.; Slater, D.R.; Whalen, K.L. Enzyme Function Initiative-Enzyme Similarity Tool (EFI-EST): A web tool for generating protein sequence similarity networks. *Biochim. Biophys. Acta* **2015**, *1854*, 1019–1037. [CrossRef] [PubMed]

73. Shannon, P.; Markiel, A.; Ozier, O.; Baliga, N.S.; Wang, J.T.; Ramage, D.; Amin, N.; Schwikowski, B.; Ideker, T. Cytoscape: A software environment for integrated models of biomolecular interaction networks. *Genome Res.* **2003**, *13*, 2498–2504. [CrossRef] [PubMed]

74. Varsani, A.; Navas-Castillo, J.; Moriones, E.; Hernandez-Zepeda, C.; Idris, A.; Brown, J.K.; Murilo Zerbini, F.; Martin, D.P. Establishment of three new genera in the family Geminiviridae: Becurtovirus, Eragrovirus and Turncurtovirus. *Arch. Virol.* **2014**, *159*, 2193–2203. [CrossRef] [PubMed]

75. Wren, J.D.; Roossinck, M.J.; Nelson, R.S.; Scheets, K.; Palmer, M.W.; Melcher, U. Plant virus biodiversity and ecology. *PLoS Biol.* **2006**, *4*, e80. [CrossRef] [PubMed]

76. Rodriguez-Negrete, E.A.; Morales-Aguilar, J.J.; Dominguez-Duran, G.; Torres-Devora, G.; Camacho-Beltran, E.; Leyva-Lopez, N.E.; Voloudakis, A.E.; Bejarano, E.R.; Mendez-Lozano, J. High-Throughput Sequencing Reveals Differential Begomovirus Species Diversity in Non-Cultivated Plants in Northern-Pacific Mexico. *Viruses* **2019**, *11*, 594. [CrossRef] [PubMed]

77. Varsani, A.; Shepherd, D.N.; Dent, K.; Monjane, A.L.; Rybicki, E.P.; Martin, D.P. A highly divergent South African geminivirus species illuminates the ancient evolutionary history of this family. *Virol. J.* **2009**, *6*, 36. [CrossRef]

78. Jones, R.A.; Coutts, B.A. Spread of introduced viruses to new plants in natural ecosystems and the threat this poses to plant biodiversity. *Mol. Plant. Pathol.* **2015**, *16*, 541–545. [CrossRef] [PubMed]

79. Alexander, H.M.; Mauck, K.E.; Whitfield, A.E.; Garrett, K.A.; Malmstrom, C.M. Plant-virus interactions and the agro-ecological interface. *Eur. J. Plant. Pathol.* **2013**, *138*, 529–547. [CrossRef]

80. Jones, R.A. Plant virus emergence and evolution: Origins, new encounter scenarios, factors driving emergence, effects of changing world conditions, and prospects for control. *Virus Res.* **2009**, *141*, 113–130. [CrossRef]

81. Nash, T.E.; Dallas, M.B.; Reyes, M.I.; Buhrman, G.K.; Ascencio-Ibanez, J.T.; Hanley-Bowdoin, L. Functional analysis of a novel motif conserved across geminivirus Rep proteins. *J. Virol.* **2011**, *85*, 1182–1192. [CrossRef]

82. Lefeuvre, P.; Moriones, E. Recombination as a motor of host switches and virus emergence: Geminiviruses as case studies. *Curr. Opin. Virol.* **2015**, *10*, 14–19. [CrossRef] [PubMed]

83. Martin, D.P.; Lefeuvre, P.; Varsani, A.; Hoareau, M.; Semegni, J.Y.; Dijoux, B.; Vincent, C.; Reynaud, B.; Lett, J.M. Complex recombination patterns arising during geminivirus coinfections preserve and demarcate biologically important intra-genome interaction networks. *PLoS Pathog.* **2011**, *7*, e1002203. [CrossRef] [PubMed]
84. Lefeuvre, P.; Martin, D.P.; Hoareau, M.; Naze, F.; Delatte, H.; Thierry, M.; Varsani, A.; Becker, N.; Reynaud, B.; Lett, J.M. Begomovirus 'melting pot' in the south-west Indian Ocean islands: Molecular diversity and evolution through recombination. *J. Gen. Virol.* **2007**, *88*, 3458–3468. [CrossRef]

85. Garcia-Andres, S.; Tomas, D.M.; Sanchez-Campos, S.; Navas-Castillo, J.; Moriones, E. Frequent occurrence of recombinants in mixed infections of tomato yellow leaf curl disease-associated begomoviruses. *Virology* **2007**, *365*, 210–219. [CrossRef] [PubMed]

86. Lefeuvre, P.; Lett, J.M.; Varsani, A.; Martin, D.P. Widely conserved recombination patterns among single-stranded DNA viruses. *J. Virol.* **2009**, *83*, 2697–2707. [CrossRef] [PubMed]

87. Gharouni Kardani, S.; Heydarnejad, J.; Zakiaghl, M.; Mehrvar, M.; Kraberger, S.; Varsani, A. Diversity of beet curly top Iran virus isolated from different hosts in Iran. *Virus Genes* **2013**, *46*, 571–575. [CrossRef] [PubMed]

88. Mandal, B.; Rao, G.P.; Baranwal, V.K.; Jain, R.K. *A Century of Plant. Virology in India*; Springer: Berlin/Heidelberg, Germany, 2017; p. 805.

89. Eini, O.; Ebadzad Sahraei, G.; Behjatnia, S.A. Molecular characterization and construction of an infectious clone of a pepper isolate of Beet curly top Iran virus. *Mol. Biol. Res. Commun.* **2016**, *5*, 101–113.

90. Heydarnejad, J.; Keyvani, N.; Razavinejad, S.; Massumi, H.; Varsani, A. Fulfilling Koch's postulates for beet curly top Iran virus and proposal for consideration of new genus in the family Geminiviridae. *Arch. Virol.* **2013**, *158*, 435–443. [CrossRef]

91. Yazdi, H.R.; Heydarnejad, J.; Massumi, H. Genome characterization and genetic diversity of beet curly top Iran virus: A geminivirus with a novel nonanucleotide. *Virus Genes* **2008**, *36*, 539–545. [CrossRef] [PubMed]

92. Soleimani, R.; Matic, S.; Taheri, H.; Behjatnia, S.A.A.; Vecchiati, M.; Izadpanah, K.; Accotto, G.P. The unconventional geminivirusBeet curly top Iran virus: Satisfying Koch's postulates and determining vector and host range. *Ann. Appl. Biol.* **2013**, *162*, 174–181. [CrossRef]

93. Hernández-Zepeda, C.; Varsani, A.; Brown, J.K. Intergeneric recombination between a new, spinach-infecting cur-
tovirus and a new geminivirus belonging to the genus Becurtovirus: First New World exemplar. *Arch. Virol.* **2013**, *158*,
2245–2254. [CrossRef] [PubMed]

94. Tahan, V.; Heydarnejad, J.; Jafarpour, B. Characterization of Beet curly top Iran virus infecting eggplant and pepper in north-eastern Iran. *Indian Phytopathol.* **2020**, *73*, 577–581. [CrossRef]

95. Briddon, R.W.; Bedford, I.D.; Tsai, J.H.; Markham, P.G. Analysis of the nucleotide sequence of the treehopper-transmitted geminivirus, tomato pseudo-curly top virus, suggests a recombinant origin. *Virology* **1996**, *219*, 387–394. [CrossRef] [PubMed]

96. Rybicki, E.P. A phylogenetic and evolutionary justification for three genera of Geminiviridae. *Arch. Virol.* **1994**, *139*, 49–77. [CrossRef] [PubMed]

97. Jeske, H.; Lutgemeier, M.; Preiss, W. DNA forms indicate rolling circle and recombination-dependent replication of Abutilon mosaic virus. *EMBO J.* **2001**, *20*, 6158–6167. [CrossRef] [PubMed]

98. Barboza, N.; Blanco-Meneses, M.; Esker, P.; Moriones, E.; Inoue-Nagata, A.K. Distribution and diversity of begomoviruses in tomato and sweet pepper plants in Costa Rica. *Ann. Appl. Biol.* **2018**, *172*, 20–32. [CrossRef]

99. Alabi, O.J.; Ogbe, F.O.; Bandyopadhyay, R.; Lava Kumar, P.; Dixon, A.G.; Hughes, J.; Naidu, R.A. Alternate hosts of African cassava mosaic virus and East African cassava mosaic Cameroon virus in Nigeria. *Arch. Virol.* **2008**, *153*, 1743–1747. [CrossRef]

100. Ferro, C.G.; Silva, J.P.; Xavier, C.A.D.; Godinho, M.T.; Lima, A.T.M.; Mar, T.B.; Lau, D.; Zerbini, F.M. The ever increasing diversity of begomoviruses infecting non-cultivated hosts: New species from Sida spp. and Leonurus sibiricus, plus two New World alphasatellites. *Ann. Appl. Biol.* **2017**, *170*, 204–218. [CrossRef]

101. Fondong, V.N.; Pita, J.S.; Rey, M.E.; de Kochko, A.; Beachy, R.N.; Fauquet, C.M. Evidence of synergism between African cassava mosaic virus and a new double-recombinant geminivirus infecting cassava in Cameroon. *J. Gen. Virol.* **2000**, *81*, 287–297. [CrossRef] [PubMed]

102. Macedo, M.A.; Albuquerque, L.C.; Maliano, M.R.; Souza, J.O.; Rojas, M.R.; Inoue-Nagata, A.K.; Gilbertson, R.L. Characterization of tomato leaf curl purple vein virus, a new monopartite New World begomovirus infecting tomato in Northeast Brazil. *Arch. Virol.* **2018**, *163*, 737–743. [CrossRef]

103. Mendez-Lozano, J.; Torres-Pacheco, I.; Fauquet, C.M.; Rivera-Bustamante, R.F. Interactions Between Geminiviruses in a Naturally Occurring Mixture: Pepper huasteco virus and Pepper golden mosaic virus. *Phytopathology* **2003**, *93*, 270–277. [CrossRef] [PubMed]

104. Pita, J.S.; Fondong, V.N.; Sangare, A.; Otim-Nape, G.W.; Ogwal, S.; Fauquet, C.M. Recombination, pseudorecombination and synergism of geminiviruses are determinant keys to the epidemic of severe cassava mosaic disease in Uganda. *J. Gen. Virol.* **2001**, *82*, 655–665. [CrossRef] [PubMed]

105. Oyeniran, K.A.; Hartnady, P.; Claverie, S.; Lefeuvre, P.; Monjane, A.L.; Donaldson, L.; Lett, J.M.; Varsani, A.; Martin, D.P. How virulent are emerging maize-infecting mastreviruses? *Arch. Virol.* **2021**, *166*, 955–959. [CrossRef]

106. De la Higuera, I.; Kasun, G.W.; Torrance, E.L.; Pratt, A.A.; Maluenda, A.; Colombet, J.; Bisseux, M.; Ravet, V.; Dayaram, A.; Stainton, D.; et al. Unveiling Crucivirus Diversity by Mining Metagenomic Data. *mBio* **2020**, *11*. [CrossRef] [PubMed]

107. Levy, H.; Fontenele, R.S.; Harding, C.; Suazo, C.; Kraberger, S.; Schmidlin, K.; Djurhuus, A.; Black, C.E.; Hart, T.; Smith, A.L.; et al. Identification and Distribution of Novel Cressdnaviruses and Circular molecules in Four Penguin Species in South Georgia and the Antarctic Peninsula. *Viruses* **2020**, *12*, 1029. [CrossRef] [PubMed]

108. Tisza, M.J.; Pastrana, D.V.; Welch, N.L.; Stewart, B.; Peretti, A.; Starrett, G.J.; Pang, Y.S.; Krishnamurthy, S.R.; Pesavento, P.A.; McDermott, D.H.; et al. Discovery of several thousand highly diverse circular DNA viruses. *Elife* **2020**, *9*. [CrossRef] [PubMed]