1. Introduction
Opuntia ficus-indica (L.) Miller is native to Mexico [
1] and is currently the most widely cultivated cactus in the world [
2]. Globally, it is primarily cultivated as forage for livestock feed, for its fruit production, and for erosion control in desert regions [
2,
3]. In Mexico, in addition to its fruit, the immature cladodes are consumed as vegetables known as “nopalitos” [
2,
4]. This species also has applications in the food industry and serves as a host for
Dactylopius coccus, which is used in the production of carminic acid, a dye employed in food, textile, and pharmaceutical products [
2,
5].
The State of Mexico is the leading producer of prickly pear cactus (nopal) in Mexico, with a cultivated area of approximately 17,500 ha. Of these, 15,800 ha are dedicated to the production of prickly pear fruit (tuna), 900 ha are dedicated to the production of nopal vegetable (nopalitos), and 800 ha are dedicated to the production of xoconostle nopal (
Opuntia joconostle) [
6]. The greatest production of nopal is concentrated in the eastern region of the State of Mexico, particularly in the municipalities of San Martín de las Pirámides, Axapusco, Otumba, Nopaltepec, Teotihuacan, Temascalapa, and Acolman [
6].
High-throughput sequencing (HTS) or next-generation sequencing has significantly increased virus detection [
7] and revolutionized the study of nucleic acids by allowing the sequencing of millions of nucleotides in a short period of time with very high redundancy (sequencing depth) [
8]. When combined with specific bioinformatics tools, HTS can be used for the detection of known viruses as well as the discovery of new viruses or viroids [
8,
9]. One of the major advantages over classical detection techniques such as enzyme-linked immunosorbent assay (ELISA) and various PCR methods is that HTS does not require prior knowledge of viral genomic information, and it allows for the identification of the plant’s entire virome in a single assay [
10].
Despite Mexico being the center of origin and domestication of
O. ficus-indica, as well as many other species of the genus
Opuntia [
1], there are limited studies regarding virus detection in these plants. In studies using high-throughput sequencing for virus detection in prickly pear cacti conducted in Mexico [
11,
12], opuntia virus 2 (OV2,
Tobamovirus genus), cactus carlavirus 1 (CCV-1,
Carlavirus genus), and opuntia potexvirus A (OPV-A,
Potexvirus genus) have been identified, along with a new species of viroid provisionally named
Opuntia viroid 1 (OVd-1,
Apscaviroid genus), the latter of which was found in nopal vegetable in Mexico City [
12].
In a study conducted on various wild species of the genus
Opuntia in Mexico to determine the presence of DNA genome viruses, squash leaf curl virus (SLCV,
Begomovirus genus) and watermelon chlorotic stunt virus (WCSV,
Begomovirus genus) were found [
13]. Additionally, a new species of geminivirus tentatively named Opuntia virus 1 (OpV1) was detected in
Opuntia santa-rita in the state of Sonora, Mexico [
14]. Other viruses detected in prickly pear cacti in Mexico through RT-PCR include Schlumbergera virus X (SVX,
Potexvirus genus) and rattail cactus necrosis-associated virus (RCNaV,
Tobamovirus genus) [
15,
16]. To date, no study has been conducted to determine the viruses present in the State of Mexico, which is the main region where prickly pear fruit, vegetables, and xoconostle are produced in the country.
In this study, we utilized HTS to investigate the virome of both nopal vegetable and prickly pear cacti. Previously, employing this approach, a new virus and a new viroid species were detected in nopal vegetable [
12]. In the present study, we leveraged the advantages of HTS technology to delve deeper and characterize the viruses and viroids in the primary production region of prickly pear cacti production regions in the State of Mexico. Our results identified various viruses and a previously reported viroid in nopal vegetable crops. Additionally, we detected two new viroid species, one provisionally named Mexican opuntia viroid (MOVd,
Pospiviroid genus) in nopal vegetable, and another provisionally named Opuntia viroid 2 (OVd-2,
Apscaviroid genus) in prickly pear cactus. This study sheds light on the diversity of viruses and viroids in
Opuntia plants in Mexico.
4. Discussion
In this study, we report on the viruses and viroids detected in the nopal vegetables and prickly pear cacti using HTS in the eastern region of the State of Mexico, Mexico. Two of these viroids are new species, one belonging to the Posviroid genus and the other to the Apscaviroid genus. In 2023, it was reported that OVd-1 (from the Apscaviroid genus) infected the nopal vegetable in Mexico, making it the first viroid reported in cacti worldwide. One of the criteria used to demarcate new viroid species is having less than 90% sequence identity. Therefore, OVd-2 should be considered a new species within the Apscaviroid genus, and MOVd should be considered a new species of the Pospiviroid genus.
Additionally, the conserved structural domains of MOVd are very similar to those of the reference genomes of Potato spindle tuber viroid (PSTVd) and Iresine viroid 1 (IrVd-1) (
Pospiviroid genus); the CCR of MOVd is identical to that of IrVd-1, while there are slight modifications in the TCR and CCR compared to those of PSTVd. In the case of OVd-2, the conserved structural domains were identical to those in the reference genome of the apple scar skin viroid (ASSVd) (
Apscaviroid genus).
Due to the presence of mixed infections of viruses and viroids in 82% of the samples analyzed individually by RT-PCR (Supplementary
Table 1), it was not possible to exclusively associate the symptoms with either of them. There is limited understanding of the role of virus/viroid coexistence, warranting a study on this interaction [
34]. In Mexico, at least three cases of mixed viroid infections have been reported: Mexican papita viroid (MPVd)/tomato chlorotic dwarf viroid and MPVd/tomato severe leaf curl virus in tomatoes under greenhouse conditions and hop stunt viroid/citrus exocortis viroid, citrus tristeza virus in oranges. The symptoms of plants with mixed infections were more severe than those plants with individual infections in all cases [
34]. Therefore, experiments on pathogenicity need to be performed, along with studies investigating the distribution, transmission modes, and effects on the production and quality of cladodes and prickly pears caused by these new viroids and viruses in nopal cultivation.
The detection of these new viroids in prickly pear cultivation reinforces the theory that Mexico, being the center of origin and domestication of
O. ficus-indica (the most cultivated cactus species in the world) and other species of the
Opuntia genus [
1], is a geographical region of viroid origin and possesses considerable biodiversity including endemic species that affect important crops such as tomatoes and high-value commercial products such as avocados [
35,
36]. The viroids detected in this study are among the most common viroids, such as the tomato bunchy top viroid, Mexican papita viroid, and avocado sunblotch viroid, which are considered endemic to Mexico [
34].
Moreover, of the detected viruses, CCV-1 and OPV-A were the most divergent (
Figure 3 and
Figure 4), while the OV2 isolates consistently clustered with previously reported isolates, forming a clade with viruses that naturally infect cacti (
Figure 6).
It is known that tobamoviruses, such as OV2, are easily transmitted mechanically [
37], so this is an important aspect to consider during crop management since many cultural practices in both prickly pear cacti and nopal vegetable involve the use of tools for pruning and harvesting that are not disinfected between plants.
OPV-A was previously detected in the nopal vegetable in Mexico City, but its range of hosts and economic impact on prickly pear cultivation are still unknown.
CCV-1 was detected alongside cactus carlavirus 2 (CCV-2) for the first time in the United States, asymptomatically, in a plant of the genus
Epiphyllum (hybrid ‘Professor Ebert’) [
38]. In Mexico, CCV-1 was detected in 93 out of 129 nopal vegetable samples collected in the states of Morelos, the State of Mexico, Hidalgo, and Mexico City [
12]. On the other hand, it is known that carlaviruses predominantly infect herbaceous plants; many cause latent or asymptomatic infections and are transmitted by aphids in a nonpersistent manner and by whiteflies [
38]. However, none of these insect groups are significant pests of prickly pear, hence transmission is not yet associated with any vector insect.
The vegetative and mechanical transmission of CCV-1, OV2 and OPV-A can contribute to their distribution in new prickly pear plantations not only in the study area but also throughout the country as the State of Mexico is the main provider of prickly pear propagative material. Therefore, phytosanitary measures should be implemented to help reduce the risk of the spread and movement of infected material carrying viruses or viroids.
Previously, a study was conducted to understand the virome of the nopal vegetable in the central region of Mexico, which included the states of Morelos, Hidalgo, Mexico City, and the State of Mexico. However, the two new viroids species reported in the present study were not detected in the State of Mexico [
12]. This suggests the need to consider more sampling sites and process a greater number of samples to generate more robust information. Additionally, in this study, we report for the first time the presence of OPV-A and OV2 in xoconostle prickly pear cacti and for the first time the presence of CCV-1 and OPV-A in prickly pear cacti.
Subsequent research on OVd-1, OVd-2, and MOVd should focus on understanding their host range, environmental conditions, or crop management practices that favor asymptomatic conditions, symptoms associated with individual and mixed infections, among other aspects. Although the main characteristic of viroids is their autonomous replication, their ability to be transmitted through grafting alone is not sufficient to demonstrate this feature [
9]; therefore, mechanical or biolistic transmission tests are necessary. Additionally, further investigations into other aspects of their transmission, such as the potential presence of an insect vector, will be needed to deepen our understanding of these viroids.
In this study, we confirmed the potential of HTS to provide an ideal methodology for determining the complete infection status of a plant by viruses and viroids. This technology does not depend on the availability of genomic sequence information such as RT-PCR and is limited only by the integrity of the reference database against which the sequences are compared [
9], as well as the correct use of different programs for bioinformatic analysis.
Figure 1.
Symptoms associated with virosis in different species of prickly pear cactus. (A-B), Chlorosis with green islands in nopal vegetable; (C), Irregular chlorotic patterns around the spines in nopal vegetable; (D), Mottling and chlorotic spots in fig prickly pear cactus. E, Ring spot symptoms in xoconostle cactus; (F) Symptomless nopal vegetable.
Figure 1.
Symptoms associated with virosis in different species of prickly pear cactus. (A-B), Chlorosis with green islands in nopal vegetable; (C), Irregular chlorotic patterns around the spines in nopal vegetable; (D), Mottling and chlorotic spots in fig prickly pear cactus. E, Ring spot symptoms in xoconostle cactus; (F) Symptomless nopal vegetable.
Figure 2.
A phylogenetic tree was constructed using the maximum likelihood method with complete genome sequences of the Tobamovirus genus obtained from GenBank (accession numbers shown). The tree was built with IQtree 2.3.1 using a color code to represent the different botanical families that are natural hosts of each virus. The identified OV2 sequences in this study are highlighted in bold. The sequence of soil-borne cereal mosaic virus (genus Furovirus) was used as an outgroup. Circles on branches indicate UFBoot support values > 70%.
Figure 2.
A phylogenetic tree was constructed using the maximum likelihood method with complete genome sequences of the Tobamovirus genus obtained from GenBank (accession numbers shown). The tree was built with IQtree 2.3.1 using a color code to represent the different botanical families that are natural hosts of each virus. The identified OV2 sequences in this study are highlighted in bold. The sequence of soil-borne cereal mosaic virus (genus Furovirus) was used as an outgroup. Circles on branches indicate UFBoot support values > 70%.
Figure 3.
A phylogenetic tree was constructed using the maximum likelihood method with complete genome sequences of the Carlavirus genus obtained from GenBank (accession numbers shown). The tree was built with IQtree 2.3.1. The CCV-1 sequence identified in this study is highlighted in bold. The sequence of Apricot vein clearing associated virus (genus Prunevirus) was used as an outgroup. Circles on branches indicate UFBoot support values > 70%.
Figure 3.
A phylogenetic tree was constructed using the maximum likelihood method with complete genome sequences of the Carlavirus genus obtained from GenBank (accession numbers shown). The tree was built with IQtree 2.3.1. The CCV-1 sequence identified in this study is highlighted in bold. The sequence of Apricot vein clearing associated virus (genus Prunevirus) was used as an outgroup. Circles on branches indicate UFBoot support values > 70%.
Figure 4.
A phylogenetic tree was constructed using the maximum likelihood method with complete genome sequences of the Potexvirus genus obtained from GenBank (accession numbers shown). The tree was built with IQtree 2.3.1. The identified OPV-A sequences in this study are highlighted in bold. The sequence of Senna severe yellow mosaic virus (genus Allexivirus) was used as an outgroup. Circles on branches indicate UFBoot support values > 70%.
Figure 4.
A phylogenetic tree was constructed using the maximum likelihood method with complete genome sequences of the Potexvirus genus obtained from GenBank (accession numbers shown). The tree was built with IQtree 2.3.1. The identified OPV-A sequences in this study are highlighted in bold. The sequence of Senna severe yellow mosaic virus (genus Allexivirus) was used as an outgroup. Circles on branches indicate UFBoot support values > 70%.
Figure 5.
Characterization of the predicted secondary structure and conserved structural elements of OVd-2. (A), Nucleotides in the upper and lower conserved central region (CCR) are indicated in red and blue, respectively. The nucleotides in green indicate the terminal conserved region (TCR) structural element. The two pairs of arrows pointing in opposite directions represent the primers used for amplifying the complete genome of the viroid; (B), The circular RNA structure of OVd-2 was determined using RT-PCR with the primer sets F90/89R and 208F/207R, each adjacent to the other in opposite directions (‘F’ = forward primers; ‘R’ = reverse primers). A healthy nopal vegetable sample was used as a negative control (NC). M, 100 bp molecular marker (Promega, Madison, USA).
Figure 5.
Characterization of the predicted secondary structure and conserved structural elements of OVd-2. (A), Nucleotides in the upper and lower conserved central region (CCR) are indicated in red and blue, respectively. The nucleotides in green indicate the terminal conserved region (TCR) structural element. The two pairs of arrows pointing in opposite directions represent the primers used for amplifying the complete genome of the viroid; (B), The circular RNA structure of OVd-2 was determined using RT-PCR with the primer sets F90/89R and 208F/207R, each adjacent to the other in opposite directions (‘F’ = forward primers; ‘R’ = reverse primers). A healthy nopal vegetable sample was used as a negative control (NC). M, 100 bp molecular marker (Promega, Madison, USA).
Figure 6.
Phylogenetic relationships of Opuntia viroid 2 (OPVd-2), Mexican opuntia viroid (MOVd) (highlighted in bold), and all current RefSeq viroid sequences available in GenBank for the family
Pospiviroidae. The phylogenetic tree was built with 1000 bootstrap replicates using maximum likelihood method and the Jukes-Cantor model in MEGA X [
9,
23]. A color code was used to represent the different genera within this family. Peach latent mosaic viroid (PLMVd) (family
Avsunviroidae) was included as an outgroup. Circles on branches indicate bootstrap values > 70% (generated from 1000 replicates).
Figure 6.
Phylogenetic relationships of Opuntia viroid 2 (OPVd-2), Mexican opuntia viroid (MOVd) (highlighted in bold), and all current RefSeq viroid sequences available in GenBank for the family
Pospiviroidae. The phylogenetic tree was built with 1000 bootstrap replicates using maximum likelihood method and the Jukes-Cantor model in MEGA X [
9,
23]. A color code was used to represent the different genera within this family. Peach latent mosaic viroid (PLMVd) (family
Avsunviroidae) was included as an outgroup. Circles on branches indicate bootstrap values > 70% (generated from 1000 replicates).
Figure 7.
The pairwise identity frequency distributions obtained with the Sequence Demarcation Tool (SDT) [
25] indicated that Opuntia viroid 2 (OVd-2) and the Mexican opuntia viroid (MOVd) are new viroid species. The sequences of OVd-2 and MOVd had less than 80% and 90% identity, respectively, with any other known viroid sequence across the entire genome. Different genera within the family
Pospiviroidae are indicated. Peach latent mosaic viroid (PLMVd, NC_003636) (family
Avsunviroidae) was included as an outgroup.
Figure 7.
The pairwise identity frequency distributions obtained with the Sequence Demarcation Tool (SDT) [
25] indicated that Opuntia viroid 2 (OVd-2) and the Mexican opuntia viroid (MOVd) are new viroid species. The sequences of OVd-2 and MOVd had less than 80% and 90% identity, respectively, with any other known viroid sequence across the entire genome. Different genera within the family
Pospiviroidae are indicated. Peach latent mosaic viroid (PLMVd, NC_003636) (family
Avsunviroidae) was included as an outgroup.
Figure 8.
Characterization of the predicted secondary structure and conserved structural elements of Mexican opuntia viroid (MOVd). (A), Nucleotides in the upper and lower conserved central region (CCR) are indicated in red and blue, respectively. The green nucleotides indicate the terminal conserved region (TCR) structural elements. The two pairs of arrows pointing in opposite directions represent the primers used for full genome amplification; (B) The circular nature of MOVd RNA was confirmed through RT-PCR with the primer sets F87/86R and 307F/306R, each adjacent to the other in opposite directions (‘F’ = forward primers; ‘R’ = reverse primers). A healthy nopal vegetable sample was used as a negative control (NC). M, 100 bp molecular marker (Promega, Madison, USA).
Figure 8.
Characterization of the predicted secondary structure and conserved structural elements of Mexican opuntia viroid (MOVd). (A), Nucleotides in the upper and lower conserved central region (CCR) are indicated in red and blue, respectively. The green nucleotides indicate the terminal conserved region (TCR) structural elements. The two pairs of arrows pointing in opposite directions represent the primers used for full genome amplification; (B) The circular nature of MOVd RNA was confirmed through RT-PCR with the primer sets F87/86R and 307F/306R, each adjacent to the other in opposite directions (‘F’ = forward primers; ‘R’ = reverse primers). A healthy nopal vegetable sample was used as a negative control (NC). M, 100 bp molecular marker (Promega, Madison, USA).
Table 1.
Primers used for the validation and detection of viruses/viroids in prickly pear cactus (fruit and nopal vegetable) from the eastern region of the State of Mexico.
Table 1.
Primers used for the validation and detection of viruses/viroids in prickly pear cactus (fruit and nopal vegetable) from the eastern region of the State of Mexico.
Virus/Viroid |
Primer |
Sequence (5′-3′) |
Amplicon (bp) |
Region |
Opuntia potexvirus A |
OPVA-RepF |
AAGCTCGCAGCATCCATCAA |
482 |
Viral replicase |
Opuntia potexvirus A |
OPVA-RepR |
GGGTGAAGGGACGGTAGTTG |
|
Cactus carlavirus 1 |
CCV-1F |
AATGGGCGCCTTTAGGTTCA |
559 |
Capsid protein |
Cactus carlavirus 1 |
CCV-1R |
AATTCCAAGCTCCCGTCAGG |
|
Opuntia virus 2 |
OV2-F |
CTTCCAAGAGTTCTAGCGCCT |
609 |
Capsid protein |
Opuntia virus 2 |
OV2-R |
ACCTGCAGGATTACCACCAC |
|
Opuntia viroid I |
OPVd_IF |
GACGGAGCGTCGAGAAGTAG |
412 |
Complete genome |
Opuntia viroid I |
OPVd_IR |
GCC GGC GCC GAA GCC CGA G |
|
Opuntia viroid II |
Opuntia_viroid_II_307F |
TCTGGCTACTACCCGGTGG |
407 |
|
Opuntia viroid II |
Opuntia_viroid_II_306R |
GCGACCAGCAGGGGAAG |
Complete genome |
Opuntia viroid II |
Opuntia_viroid_II_86R |
CCGGGGATCCCTGAAG |
407 |
|
Opuntia viroid II |
Opuntia_viroid_II_87F |
GGAAACCTGGAGCGAACTC |
Complete genome |
Opuntia viroid III |
Opuntia_viroid_III_F90 |
GAAGGCAGCTGAGTGGAG |
319 |
|
Opuntia viroid III |
Opuntia_viroid_III_R89 |
GTCGACGACGACAGGTGA |
Complete genome |
Opuntia viroid III |
Opuntia_viroid_III_208F |
AGGGCCACACTCGGTG |
319 |
Complete genome |
Opuntia viroid III |
Opuntia_viroid_III_207R |
CCGGAGGCAGAGGAGAG |
|
Table 2.
Locations sampled for the detection of viruses and viroids in different types of prickly pear cactus in the eastern region of the State of Mexico.
Table 2.
Locations sampled for the detection of viruses and viroids in different types of prickly pear cactus in the eastern region of the State of Mexico.
Municipality |
Location |
Type of Nopal |
Fruit |
Nopal Vegetable |
Xoconostle |
wild prickly pear cactus |
Axapusco |
Cuautlacingo |
3 |
38 |
0 |
0 |
Axapusco |
4 |
0 |
0 |
0 |
Temascalapa |
Santa Ana Tlachiahualpa |
11 |
1 |
0 |
0 |
Nopaltepec |
San Felipe Teotitlán |
20 |
0 |
4 |
2 |
Teotihuacán |
Teotihuacán |
0 |
4 |
0 |
0 |
Total |
|
37 |
43 |
4 |
2 |
Table 3.
Comparison of total or partial genome sequences of the viruses and viroids detected through HTS in nopal vegetable and prickly pear cactus samples with the most similar reference sequence available in NCBI’s GenBank.
Table 3.
Comparison of total or partial genome sequences of the viruses and viroids detected through HTS in nopal vegetable and prickly pear cactus samples with the most similar reference sequence available in NCBI’s GenBank.
HTS Sample |
Virus/Viroid Detected* |
GenBank Accession Number |
Isolation |
Genome Segment |
Reference Sequence Accession Number (NCBI GenBank) |
% Identity with Reference |
nopal vegetable |
CCV-1 |
PP579657 |
EM_C2 |
complete |
KU854930.4 |
92.14% |
OPV- A |
PP579658 |
EM_A2 |
complete |
OQ240443.1 |
85.50% |
OV2 |
PP579659 |
EM_T2 |
complete |
NC_040685.2 |
98.98% |
OVd 1 |
PP579662 |
EDMEX-V1 |
complete |
OQ240445.1 |
98.54% |
prickly pear cactus |
CCV-1 |
PP579654 |
EM_C1. |
**RdRp |
KU854930.4 |
88.83% |
OPV-A |
PP579655 |
EM_A1 |
complete |
OQ240443.1 |
90.16% |
OV2 |
PP579656 |
EM_T1 |
complete |
NC_040685.2 |
98.98% |
Table 4.
New viroid species detected by HTS in prickly pear cactus and nopal vegetable.
Table 4.
New viroid species detected by HTS in prickly pear cactus and nopal vegetable.
HTS Sample |
Viroid |
GenBank Accession Number |
Isolate |
Closest Viroid in BLASTn Analysis |
Reference Sequence, Accession Number (NCBI GenBank) |
% Identity with Reference (% Query Cover) |
nopal vegetable |
Mexican opuntia viroid |
PP579660 |
V1 |
Iresine viroid 1 |
OM108483.1 |
83.33 (99) |
prickly pear cactus |
Opuntia viroid 2 |
PP579661 |
EDMEX-T1 |
Grapevine latent viroid |
MG770884.1 |
82.81 (38) |
Table 5.
Location of conserved domains of viroids in the reference sequence of the Pospiviroid and Apscaviroid genera compared to the new viroid species detected in this study. Discrepancies in the sequences are underlined and in bold.
Table 5.
Location of conserved domains of viroids in the reference sequence of the Pospiviroid and Apscaviroid genera compared to the new viroid species detected in this study. Discrepancies in the sequences are underlined and in bold.
Viroid* |
TCR |
CCR Upper Strand |
CCR Lower Strand |
MOVd (PP579660.1) |
GGUUCCUGUGG |
CUUCAGGGAUCCCCGGGGAAACCUGGAG |
ACUACCCGGUGGAUACAACUGUAGCU |
PSTVd (NC_002030.1) |
GGUUCCUAUGG |
CUUCAGGGAUCCCCGGGGAAACCUGGAG |
ACUACCCGGUGGAAACAACUGAAGCU |
IrVd-1 (NC_003613.1) |
GGUUCCAAUGG |
CUUCAGGGAUCCCCGGGGAAACCUGGAG |
ACUACCCGGUGGAUACAACUGUAGCU |
OVd-2 (PP579661.1) |
GGUUCCUGUGG |
UCGUCGUCGACGAAGG |
CCGCUAGUCGAGCGGAC |
ASSVd (NC_001340) |
GGUUCCUGUGG |
UCGUCGUCGACGAAGG |
CCGCUAGUCGAGCGGAC |