Detection of RNA viruses in potato seed-tubers from northern Antioquia (Colombia)¹

ABSTRACT

The most important viruses infecting potato in Colombia are PVY, PLRV, PVS, PYVV, PVV, PVX, PMTV, and PVB. Unfortunately, recent investigations on the health status of uncertified and certified seed-tubers of Solanum tuberosum cv. Diacol Capiro and Solanum phureja cv. Criolla Colombia used by farmers in eastern Antioquia (Colombia) have revealed alarming levels of infection by these viruses in all types of seeds distributed in the region. A similar study on northern Antioquia, testing all these eight viruses has not been performed. To close this gap, we present a complementary study on the prevalence of viruses in seed-tubers distributed and commercialized in this region using a combination of RT-qPCR and high-throughput sequencing. Our results further confirm the poor quality of seed-tubers used in Colombia and suggests infection by tobacco necrosis virus D (TNV-D), which is the first report for this virus in Colombia.

Keywords
plant viruses; RNAseq; RT-qPCR; seed certification; Solanaceae

INTRODUCTION

The province of Antioquia is the fourth most important potato producer in Colombia with fields comprising a total of 8619 hectares with annual yields estimated at 20.3 t/ha (Agronet, 2021) and dominated by the tetraploid cultivar Solanum tuberosum cv. Diacol Capiro, and the diploid Solanum phureja cv. Criolla Colombia (Porras & Herrera, 2015; MADR, 2019). Within Antioquia, potato production is concentrated in the eastern and northern regions, the latter of which comprises fields totaling 3585 ha across the municipalities of Santa Rosa de Osos, Entrerríos, San Pedro de los Milagros, and Belmira. Recent investigations on the viruses transmitted through seed-tubers in potato fields in eastern Antioquia have revealed a high prevalence of RNA viruses in formal- and informally- produced planting material, and it is believed that high viral loads are a major factor explaining the lower potato productivity of the region (Sierra et al., 2020; 2021; Gallo et al., 2019; 2021b). For example, in a study on seed-tubers from the Diacol Capiro cultivar it was found a prevalence for potato yellow vein virus (PYVV) of 73.3%, and up to 90% for potato virus Y (PVY) (Sierra et al., 2021). With respect to the cultivar Criolla Colombia, the prevalence of viruses was more dramatic as some viruses such as potato virus S (PVS) and PYVV reached 100%, potato leaf roll virus (PLRV) 47%, PVY 47%, potato virus X (PVX) 40%, and potato virus V (PVV) 40% (Sierra et al., 2020).

It is a well-stablished fact that the widespread use of low-quality seed-tubers in Colombia is a major factor in the high prevalence of viral diseases in potato fields throughout the country (Gallo et al., 2021b; Sierra et al., 2021). It is estimated that only 3-5% of farmers in Colombia use certified planting material (Fedepapa, 2019; MADR, 2019) which, unfortunately, is not of sufficient quality with respect to viral infections either (Sierra et al., 2021). To plant a new field, the average Colombian farmer uses seed-tubers produced by the same farmers from a preceding season or purchased in informal markets (MADR, 2019). Formal seed-tubers, on the other hand, are produced by specialized farmers in the highland plains of Cundinamarca, Boyacá and Nariño, where there is reduced insect pressure, from in vitro stock plants presumed to be free of viruses (ICA, 2015; Fedepapa, 2019). Unfortunately, the use of low sensitivity detection methods, and the circulation of local virus species and variants that are not detectable with available commercial serological kits has resulted in formal seeds with unacceptable levels of prevalence for some limiting viruses such as PLRV, PVY, PVX, PVS, and PYVV (Gallo et al., 2021b; Sierra et al., 2021).

PVY, PLRV, PVS, PYVV, PVV, PVX, potato mop-top virus (PMTV), and potato virus B (PVB) are the most prevalent viruses affecting potato in Colombia (Guzmán et al., 2010; Vallejo et al., 2016; Sierra et al., 2020; 2021; Gallo et al., 2019; 2021b; Giraldo et al., 2022). However, only one study on the prevalence of these eight viruses in seed-tubers has been performed to this day (García, 2022). On that work, the author performed a pilot seed-indexing study using the tuber sprouting test through RT-qPCR with an updated set of primers on potato producing regions in eastern Antioquia. The author found that viruses such as PYVV and PMTV were highly prevalent in seed lots from both cultivars; PLRV, PVY, and PVB were highly prevalent in all cv. Diacol Capiro lots, but mostly absent in cv. Criolla Colombia, and PVV was almost exclusively found in cv. Criolla Colombia. With respect to PVS, and PVX these viruses were only found at low prevalence in the Diacol Capiro cultivar. In agreement, with previous work, multiple infections were also common in both cultivars, and no uninfected tubers were identified in any of the tested lots (García, 2022). Little is known about the health status of potato seeds commercialized in the northern producing regions of the province of Antioquia. To close this gap, we have performed this study on the prevalence of PLRV, PVY, PVX, PVS, PVV, PYVV, PMTV, and PVB in potato seeds stocks from the municipalities of San Pedro de los Milagros, Santa Rosa de Osos and Entrerríos by real-time RT-PCR (RT-qPCR) and high-throughput sequencing (HTS).

MATERIALS AND METHODS

Plant material

This work involved sampling and analysis of seed-tuber samples from the municipalities of San Pedro de los Milagros (6° 27’ 0” N, 75° 33’ 0” W), Entrerríos (6° 33’ 55.44” N, 75° 31’ 0.84” W), and Santa Rosa de Osos (6° 38’ 50.57” N, 75° 27’ 37.12” W) in northern Antioquia. Samples were collected from eight uncertified S. tuberosum cv. Diacol Capiro seed lots, eight uncertified S. phureja cv. Criolla Colombia seed lots, and eight certified S. tuberosum cv. Diacol Capiro seed lots. Due to the local nature and short dormancy period of the S. phureja, certified seeds for this cultivar were not available in the region. Sprouting in dormant seed-tubers of cv. Diacol Capiro was induced by spraying a diluted solution of gibberellic acid at 50 ppm (Sigma-Aldrich) in distilled water, air drying, and storage in darkness at room temperature for 15-20 days. Finally, from each seed-tuber lot, three separate tubers and three sprouts per tuber were used for virus indexing.

RT-qPCR

Total RNA was extracted from 100 mg of the top two centimeters at the tip of each individual tuber sprout, following previous methods (García et al., 2016; Sierra et al., 2020). Prior to RNA extraction, plant tissue was ground in liquid nitrogen. RNA was purified using the GeneJET Plant RNA Purification Mini Kit (Thermo Fisher Scientific) and eluted in 40 µL of DEPC-treated water. RNA concentration and purity were determined using the A260/280 ratio measured in a Nanodrop 2000C spectrophotometer (Thermo Fisher Scientific). cDNA was synthesized using 80U of RevertAid Reverse Transcriptase (RT) (Thermo Fisher Scientific), RT buffer (1X), dNTPs at 1 mM, 20 U of RiboLock RNase Inhibitor, 1 µg of total RNA and 20 pmoles of the reverse primer in a reaction volume of 20 µL at 42 ^oC for 1 h. For PYVV, PLRV and PMTV, primers PYVV_R_CP, PLRV-R and 123-end were employed as reverse primers, respectively; for PVY, PVV, PVB, PVX and PVS, an Oligo-(dT)18 was employed (Table 1). cDNA synthesis was performed in a T3 thermal cycler (Biometra) at 65 °C for 5 min, 42 °C for 60 min, and a final step at 70 °C for 10 min. The real-time PCR (qPCR) step was performed in a final volume of 14 µL with 6.25 µL of the Maxima SYBR Green/ROX kit (Thermo Fisher Scientific) using 100 ng of cDNA and primers at 0.3 µM concentration (Table 1). The temperature profile consisted of 10 min at 95 ^oC, followed by 35 cycles at 95 ^oC for 15 s, and 52 ^oC for 60 s. Amplifications were carried out in a Rotor-Gene Q-5plex cycler (Qiagen). In addition to threshold cycle (Ct), the specificity of the RT-qPCR results was confirmed by High Resolution Melting (HRM) in the 50-99 ^oC range. Positive and negative controls were included in each test and consisted of samples previously isolated from infected and virus-free potato tissues. The identity of amplicons was confirmed for selected samples by Sanger sequencing.

Table 1
List of primers used in this study

Virus	Primers	Sequence	Source
PVY	PVY_Col_RV2	5’ GCA TAC TTG GAG ARA CAT CYT CG 3’	García et al. (2022)
	PVY_Col_qFV2	5’ TCG GAT TGG AYG GYG GCA T 3’
PVV	qPVV_phu_R	5’ CAT CCC GCT CCT CAA C 3’	Álvarez et al. (2016)
	PVV_phu_F	5’ ATG CTG GAA AAG ATC CAG C 3’
PVX	PVX_101-2_RP	5’ GCT TCA GAC GGT GGC CG 3’	Agindotan et al. (2007)
	PVX_101-2_FP	5’ AAG CCT GAG CAC AAA TTC GC 3’
PVS	PVS_gen_RV2	5’ ATT GRA ANC CCA TRG CTT GC 3’	García et al. (2022)
	qPVS_gen_F	5’ CGG GRY TGA GRA ARG TGT G 3’
PLRV	PLRV-Antisense	5’ TTG AAT GCC GGA CAG TCT GA 3’	Yang et al. (2014)
	PLRV-sense	5’ GCT CAA GCG AGA CAT TCG TG 3’
	PLRV-R	5’ GCA ATG GGG GTC CAA CTC CAA CTC AT 3’	Singh et al. (1995)
PYVV	qPYVV_R_CP	5’ AGG TCT CAG GAT CTG GAT CAA CT 3’	Álvarez et al. (2017)
	PYVV_F_CP	5’ TCA GGT TAG AGC AGA CAG AGG 3’
	PYVV_R_CP	5’ AGT TGC TGC ATT CTT GAA CAG G 3’
PMTV	PMTV-1948F	5’ GTG ATC AGA TCC GCG TCC TT 3’	Mumford et al. (2000)
	PMTV-2017R	5’ CCA CTG CAA AAG AAC CGA TTT C 3’
	123-end	5’ GTG AAC CAC GGT TTA RCC CTG KAA GC 3’	Savenkov et al. (1999)
PVB	Neponew_pol_F	5’ TCG GCA CAA ACC CCT ATT CC 3’	Giraldo et al. (2022)
	Neponew_qpol_R	5’ TCA ATA GAC CAT CGA ACC CAC T 3’

High-throughput sequencing

cDNA libraries for HTS were prepared with the TruSeq Stranded Total RNA LT Sample Prep kit using bulked total RNA from certified (BST19) and uncertified (BST20) S. tuberosum cv. Diacol Capiro seeds, and uncertified (BSP18) S. phureja cv. Criolla Colombia. Ribosomal RNA (rRNA) was removed with the Ribo-Zero Plant kit (Illumina). The quality and quantity of the nucleic acid preparations was assessed with a 2100 Bioanalyzer (Agilent Technologies). Sequencing was performed with the Illumina NovaSeq system service provided by Macrogen (South Korea). The quality of each dataset was verified with fastqc (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/), and the low-quality bases (Phred < 30) were removed using Seqtk (https://github.com/lh3/seqtk). Preliminary detection of viruses was done with the PVDP platform (Gutiérrez et al., 2021). Scaffolds from low coverage genomes were assembled with Magic-BLAST (Boratyn et al., 2019). Viruses with high coverage were assembled de novo with rnaSPAdes for each dataset (Bushmanova et al., 2019). Contigs and scaffolds were verified for assembly errors with the Integrative Genomics Viewer (IGV) (Robinson et al., 2011). Phylogenetic analyses were carried out in MEGA 11 (Tamura et al., 2021) using the Neighbor-Joining method (Saitou & Nei, 1987) with 1000 bootstrap replicates. Sequences were aligned with MUSCLE (Edgar, 2004). Ambiguous positions were removed for each sequence pair. Evolutionary distances were computed using the Tamura-Nei method with a gamma distribution (Tamura, 1992).

RESULTS

Detection of viruses by RT-qPCR

With exception of PVV and PVB, RT-qPCR revealed high levels of prevalence for most viruses in S. tuberosum cv. Diacol Capiro tubers. For uncertified tubers of this cultivar PVX, PYVV, and PMTV were detected in all samples (100%); PVY and PVS were detected in seven samples (87.5%); and PLRV in five samples (62.5%). All uncertified S. tuberosum cv. Diacol Capiro tubers tested negative for PVV and PVB (Figure 1A). Viruses detected in uncertified Diacol Capiro tubers were also detected in certified tubers at 100% prevalence, with the addition of one sample from San Pedro that tested positive for PVB (Figure 1A). On the other hand, the eight viruses were detected in S. phureja cv. Criolla Colombia; PVY, PVV, PVS, PYVV, and PLRV were found in all samples (100%), PVX in seven samples (87.5%), PMTV in six samples (75%), and PVB in two samples (25%). The most remarkable difference was the high prevalence of PVV in S. phureja which was not detected in S. tuberosum (Figure 1B). Coinfections were common in all tested seed-tuber lots. In S. tuberosum, samples were coinfected with at least four viruses with the combination PVY/PVS/PVX/PYVV/PLRV/PMTV being the most common as it was observed in seven certified samples, and five uncertified samples. One certified sample was infected with all viruses, except for PVV (Figure 1C). With respect to S. phureja, all samples were infected with at least six viruses, however, the large majority were infected with the seven-virus combination PVY/PVV/PVS/PVX/PYVV/PLRV/PMTV.

Figure 1
Virus prevalence in potato seed-tuber lots from northern Antioquia. A) Eight RNA viruses known to infect potato in Colombia were tested by RT-qPCR in certified, and uncertified seed-tubers from S. tuberosum cv. Diacol Capiro and B) S. phureja cv. Criolla Colombia. Numbers correspond to Ct values. C) Virus coinfections observed for each sample type. ER: Entrerríos; SP: San Pedro de los Milagros; SR: Santa Rosa de Osos. TC: S. tuberosum certified tubers; TU: S. tuberosum uncertified tubers; PU: S. phureja uncertified tubers.

High-throughput sequencing

RNAseq analysis revealed a higher proportion of viral reads, 6.21% and 0.95% in the S. tuberosum uncertified and certified samples, respectively. For S. phureja the proportion of viral reads was 1.01% (Figure 2). In uncertified S. tuberosum cv. Diacol Capiro seed lots, the dominant viruses were PLRV at 2165 reads per million (RPM), the PMTV segments at 105-287.9 RPM and PVY at 59334 RPM (Figure 2). Reads attributed to all RNA segments of PYVV were also detected in this dataset at 2.9 (RNA1), 18.9 (RNA2), and 8.6 (RNA3) RPM. Viral loads were lower in the certified tubers, in this case, the dominant viruses were also PVY at 9398.6 RPM and PLRV at 51 RPM (Figure 2). The set of viruses present in S. phureja seed lots was different from that of S. tuberosum. In this case, the dominant viruses were PVX (8889 RPM), and PVV (1158.6 RPM); PVS was detected a lower load (15.4 RPM). Interestingly, a significant number of reads with about 96% nucleotide sequence identities to tobacco necrosis virus D (TNV-D) (Betanecrovirus, Tombusviridae) were also identified in this dataset at 10.3 RPM, which would be the first report for this virus in Colombia (Figure 2).

Figure 2
HTS detection of viruses in potato seeds-tubers in northern Antioquia. Detection of viruses was performed using the PVDP package. The proportion of virus-like sequences is shown on the left in magenta and corresponds to the proportion of sequences sharing high nucleotide sequence identities with currently known plant viruses at NCBI. Boxplots illustrate the relative genome position of mapped reads against each virus genome; darker colors indicate higher sequence identity. RPM and total number of reads for each identified virus are shown on the right.

Phylogenetic analysis of scaffolds and/or contigs revealed the phylogenetic affinity of the viruses detected by RNAseq with previously known isolates in Colombia. For PLRV, sequences from both certified (BST19) and uncertified (BST20) S. tuberosum seed-tubers were part of the main clade found for this virus worldwide (Guyader & Ducray, 2002; Mesa et al., 2016), and closely related to isolate M1 obtained from S. quitoense in Antioquia by Gallo et al. (2021a). Analysis of the PMTV RNA3 scaffold obtained from the uncertified S. tuberosum seed-lots (BST20) formed an independent clade between isolates of worldwide distribution (Gil et al., 2011; 2016), and local isolates infecting S. tuberosum (CO2 and CO1) or P. peruviana (M7 and M2) recently reported by Gallo et al. (2020) (Figure 3).

Figure 3
Neighbor-joining trees of viruses detected by RNAseq in S. tuberosum cv. Diacol Capiro in seed-tubers from northern Antioquia. Each panel shows the phylogenetic tree obtained for partial genome sequences (bold) derived from the RNAseq data from certified (BST19), and uncertified (BST20) potato seed-tubers. Bootstrap values are shown above the tree branches. NCBI accession codes are indicated in parentheses, together with information about the host, when available.

PYVV sequence from the uncertified S. tuberosum seed lots (BST20) clustered within a clade containing isolates infecting different hosts such as S. lycopersicum, S. phureja, S. tuberosum, S. quitoense, and P. peruviana in Colombia (Guzmán-Barney et al., 2012; Álvarez et al., 2017; Gallo et al. 2020; 2021b) (Figure 3). Assembly of the PVY sequences resulted in three contigs clustering with the PVY^N, and the recombinant PVY^NTN strains previously detected by Muñoz et al. (2016a) in S. lycopersicum, Álvarez et al. (2018) in P. peruviana and S. tuberosum, including the isolate Yarumal from northern Antioquia (Muñoz et al., 2016b) (Figure 3). Surprisingly, uncertified seeds contained only the PVY^NTN strain while, certified S. tuberosum seeds were also carriers of the PVY^N strain.

Analysis of the S. phureja assemblies from the uncertified seed-tuber lots confirmed the presence of PVS, PVV, and PVX, in addition to TNV-D (Figure 4). The PVS scaffold (PVS-BSP18) clustered within the Andean PVS (PVS^A) group of strains infecting different solanaceous crops in Antioquia such as S. tuberosum, S. quitoense, and S. lycopersicum but formed a distinct clade with isolate RL5 reported to infect S. phureja in Antioquia (Gutiérrez et al., 2013; Vallejo et al., 2016) (Figure 4). The PVV assembly (PVV-BSP18), on the other hand, clustered within a distinct clade comprising other S. phureja isolates (Álvarez et al., 2016; Gutiérrez et al., 2016) in addition to isolate from P. peruviana identified in Antioquia (Álvarez et al., 2018). With respect to PVX, the assembled contig clustered within the Eurasian Clade at 100% bootstrap together with several worldwide isolates infecting S. tuberosum (García et al., 2016) (Figure 4). Finally, the partial genome sequence of TNV-D from uncertified (BSP18) S. phureja seeds was part of a clade comprising several TNV-D isolates infecting a wide diversity of hosts such as Valerianella locusta, Nicotiana tabacum and Glycine max; this TNV-D isolate was also clearly separated from the other two approved members of genus Betanecrovirus: leek white stripe virus (LWSV) and beet black scorch virus (BBSV) (ICTV, 2021) (Figure 4).

Figure 4
Neighbor-joining trees of viruses detected by RNAseq in S. phureja cv. Criolla Colombia in seed-tubers from northern Antioquia. Each panel shows the phylogenetic tree obtained for partial genome sequences (bold) derived from the RNAseq data from uncertified (BSP18) potato seed-tubers. Bootstrap values are shown above the tree branches. NCBI accession codes are indicated in parentheses, together with information about the host, if available.

DISCUSSION

In this study, we have evaluated the prevalence of RNA viruses in seed-tubers from S. tuberosum cv. Diacol Capiro and S. phureja cv. Criolla Colombia, used by farmers in northern Antioquia. Our results suggest a high prevalence of viruses in seeds from all sources. With respect to cv. Criolla Colombia, our analyses confirmed infection with all the viruses tested and the high prevalence of PVV in contrast to S. tuberosum (Álvarez et al., 2016). A comparison with similar studies performed in the eastern region of Antioquia (Sierra et al., 2020; 2021; Gallo et al., 2021b; García, 2022), confirms that seed-tuber degeneration resulting from the accumulation of viruses is a serious problem affecting the potato production in Antioquia that must be addressed urgently by the corresponding phytosanitary agencies. In perspective, the maximum productivity of cv. Diacol Capiro is estimated at 40 t/ha (Porras & Herrera, 2015) which is almost twice the productivity of northern Antioquia of 24 t/ha; with respect to S. phureja the maximum productivity is approximately 25 t/ha (Rodríguez et al., 2009; Seminario et al., 2018) but current productivity in fields from Antioquia rarely surpass 15 t/ha (Agronet, 2021). Of particular concern are the levels of infection of certified cv. Diacol Capiro tubers, which are expected to have a lower prevalence of viruses but tested at 100% for PVY, PVS, PVX, PYVV, and PLRV as current legislation limits their prevalence at 5% for PLRV, PVY, PVS, and PVX; and 1% for PYVV (ICA, 2015); it is important to clarify that these values correspond to prevalence as measured by visual inspection and/or serological tests which are less sensitive than the methods used in this work. Another worrying fact is the high number of different viruses detected in certified seed lots, which were infected with at least six viruses. We strongly believe that the current potato seed certification program in Colombia should be reformulated to include more sensitive detection methods, particularly during the screening phase of the seed stock plants used to produce the material for further certification steps. The plantation of virus-free tubers can substantially improve yields in potato, particularly when combined with methods that reduce potential sources of inoculum such as early-season crop hygiene, roguing, and positive selection (Gildemacher et al., 2011; Schulte-Geldermann et al., 2012; Thomas-Sharma et al., 2016). Some reports indicate that planting reliable certificated seed-tubers can result in about an initial 30% increment in yields (Halterman et al., 2012; Frost et al., 2013; Thomas-Sharma et al., 2016; 2017), but can be increased further to more than 70% with the continued used of high-quality seed, as evidenced by data from western Europe, the United States, and Canada (Halterman et al., 2012; Frost et al., 2013; Thomas-Sharma et al., 2016; 2017).

The production of healthy seed-tubers should be a concerted effort involving farmers, phytosanitary surveillance institutions, and basic research groups that results in the routine implementation of highly sensitive diagnostic techniques such as RT-qPCR or isothermal amplification tests such as RPA and LAMP for monitoring the prevalence of viruses in the field and the validation of virus clean-up protocols (Raigond et al., 2020; Kumar et al., 2021; Schumpp et al., 2021). In addition, HTS methods should be used as a complementary technique for detecting the emergence of new virus strains and species and for monitoring the movement of viruses across regional or country borders. The usefulness of HTS as a complementary diagnostic method has been demonstrated in this work, allowing the detection of PLRV, PMTV, PYVV, and PVY in seed-tubers of S. tuberosum cv. Diacol Capiro, and of PVX, PVV, and PVS in S. phureja tubers. As a proof of concept, HTS also allowed the identification of TNV-D, a virus that was not initially included as part of this study and therefore should be investigated further to clarify its impact on S. phureja. TNV-D is a Betanecrovirus (Tombusviridae) first described in N. tabacum that can infect a wide spectrum of hosts since biological studies have shown that this virus can be transmitted to at least 88 dicotyledonous and monocotyledon species in 37 families (Price, 1938). TNV-D is naturally transmitted by zoospores of the root-infecting chytrid fungus Olpidium brassicae (Teakle & Gold, 1963) and despite its early discovery in 1935 (Smith & Bald, 1935), only a few isolates have been sequenced: TNV-D infecting french bean (Phaseolus vulgaris) from England, TNV-DH infecting Nicotiana clevelandii from Hungary, TNV-DP infecting olive trees from Portugal and the Manitoban isolate (TNV-D^MB) from soybean (Díaz-Cruz et al., 2017). Genome analysis of TNV-D^MB revealed a genome of 3,742 nt encoding six proteins: 22-kDa and 82-kDa proteins products required for RNA replication (p22 and p82), three 7-kDa proteins (71, 7a, and 7b) believed to be involved in cell-to-cell movement, and a coat protein (CP) gene (Díaz-Cruz et al., 2017). Infection of potato with TNV-D is uncommon, however, it has been included as part of the EPPO Certification Schemes for potato seeds (PM4/028-1) and all nuclear stock microplants commercially produced in Europe should test negative for this virus (EPPO, 1999). Future work must investigate the prevalence of TNV-D, and its associated vector O. brassicae, in potato fields in northern Antioquia.

In more general terms, we believe that similar studies like this one should be extended to other potato cultivating regions in Colombia and to include additional cultivars such as Parda Pastusa, Superior or Pastusa suprema, widely grown in Cundinamarca, Boyacá, and Nariño (MADR, 2019), to better understand the current phytosanitary status of this crop and have a more complete picture of the virome of potato in Colombia. Potato is affected by more than fifty viruses across the world (Kreuze et al., 2020) and it is likely that more viruses remain to be discovered in Colombia and the rest of the South American Andes, the center of diversity of potato.

CONCLUSIONS

A high level of prevalence of viruses: PVX, PYVV, PMTV, PVY, PVS and PLRV was found in potato seed-tubers from the northern region of Antioquia (Colombia) in both certified and uncertified material of cultivar Diacol Capiro. Additionally to these viruses, PVV and PVB were also found in uncertified tubers of S. phureja cv. Criolla Colombia. Using high-throughput sequencing, the partial genome sequence of tobacco necrosis virus D (TNV-D) was assembled from uncertified S. phureja. This is the first report for TNV-D in Colombia.

ACKNOWLEDGEMENTS, FINANCIAL SUPPORT AND FULL DISCLOSURE

This work was funded by Fondo de Ciencia Tecnología e Innovación del Sistema General de Regalías del Departamento de Antioquia, under grant #1101-805-62787 (Agreement #4600007658-779) and supervision of Ministerio de Ciencia, Tecnología e Innovación (Colombia). Plant samples were collected using permission RGE152-27 from Ministerio del Medio Ambiente y Desarrollo Sostenible (Resolution 0208, 9/03/2020). The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

• 1
  This article is part of the first author’s Master’s Dissertation.

REFERENCES

• Agindotan BO, Shiel PJ & Berger PH (2007) Simultaneous detection of potato viruses, PLRV, PVA, PVX and PVY from dormant potato tubers by TaqMan real-time RT-PCR. Journal of Virological Methods, 142:01-09.
• Agronet (2021) Red de información y comunicación del sector Agropecuario Colombiano. Available at: <https//www.agronet.gov.co/estadistica/Paginas/home.aspx>. Accessed on: October 15^th, 2021.
  » https//www.agronet.gov.co/estadistica/Paginas/home.aspx
• Álvarez D, Gutiérrez P & Marín M (2016) Caracterización molecular del Potato virus V (PVV) infectando Solanum phureja mediante secuenciación de nueva generación. Acta Biológica Colombiana, 21:521-531.
• Álvarez D, Gutiérrez-Sánchez P & Marín-Montoya M (2017) Genome sequencing of Potato yellow vein virus (PYVV) and development of a molecular test for its detection. Bioagro, 29:03-14.
• Álvarez N, Jaramillo H, Gallo Y, Gutiérrez P & Marín M (2018) Molecular characterization of Potato virus Y (PVY) and Potato virus V (PVV) isolates naturally infecting Cape gooseberry (Physalis peruviana) in Antioquia, Colombia. Agronomía Colombiana, 36:13-23.
• Boratyn GM, Thierry-Mieg D, Busby B & Madden TL (2019) Magic-BLAST, an accurate DNA and RNA-seq aligner for long and short reads. BMC Bioinformatics, 20:405.
• Bushmanova E, Antipov D, Lapidus A & Prjibelski AD (2019) rnaSPAdes: a de novo transcriptome assembler and its application to RNA-Seq data. GigaScience, 8:01-13.
• Díaz-Cruz GA, Smith CM, Wiebe KF & Cassone BJ (2017) First complete genome sequence of Tobacco necrosis virus D isolated from soybean and from North America. Genome Announcements, 5:e00781.
• Edgar RC (2004) MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Research, 32:1792-1797.
• EPPO Global Database (1999) Tobacco necrosis virus (TNV000). Available at: <https://gd.eppo.int/taxon/TNV000/documents>. Accessed on: April 25^th, 2021.
  » https://gd.eppo.int/taxon/TNV000/documents
• Fedepapa (2019) Informe de gestión. Vigencia 2019. Available at: <https://fedepapa.com/wp-content/uploads/2020/05/INFORME-DE-GESTIO%CC%81N-VIGENCIA-2019.pdf>. Accessed on: October 20^th, 2021.
  » https://fedepapa.com/wp-content/uploads/2020/05/INFORME-DE-GESTIO%CC%81N-VIGENCIA-2019.pdf
• Frost KE, Groves RL & Charkowski AO (2013) Integrated control of potato pathogens through seed potato certification and provision of clean seed potatoes. Plant Disease, 97:1268-1280.
• Gallo Y, Sierra A, Donaire L, Aranda MA, Gutiérrez PA & Marín MA (2019) Coinfección natural de virus de ARN en cultivos de papa (Solanum tuberosum subsp. Andigena) en Antioquia (Colombia). Acta Biológica Colombiana, 24:546-560.
• Gallo Y, Marín M & Gutiérrez PA (2020) Detection of RNA viruses in Cape gooseberry (Physalis peruviana L.) by RNAseq using total RNA and dsRNA inputs. Archives of Phytopathology and Plant Protection, 53:395-413.
• Gallo Y, Marín M & Gutiérrez PA (2021a) Detection of RNA viruses in Solanum quitoense by high-throughput sequencing (HTS) using total and double stranded RNA inputs. Physiological and Molecular Plant Pathology, 113:101570.
• Gallo Y, Sierra A, Marín M & Gutiérrez PA (2021b) Prevalencia de cinco virus de ARN en tubérculos-semilla de papa cultivados en Antioquia (Colombia). Biotecnología en el Sector Agropecuario y Agroindustrial, 19:66-78.
• García D, Olarte MA, Gutiérrez P & Marín MA (2016) Detección serológica y molecular del Potato virus X (PVX) en tubérculos-semilla de papa (Solanum tuberosum L. y Solanum phureja Juz. y Bukasov) en Antioquia. Revista Colombiana de Biotecnología, 18:104-111.
• García A (2022) Detección molecular y limpieza de virus en material de siembra de Solanum tuberosum y S. phureja Master Dissertation. Universidad Nacional de Colombia, Medellín. 358p.
• Gil JF, Cotes JM, González EP & Marín M (2011) Caracterización genotípica de aislamientos colombianos del potato mop-top virus (PMTV, Pomovirus). Actualidades Biológicas, 33:69-84.
• Gil JF, Adams I, Boonham N, Nielsen SL & Nicolaisen M (2016) Molecular and biological characterization of Potato mop-top virus (PMTV, Pomovirus) isolates from the potato-growing regions of Colombia. Plant Pathology, 65:1210-1220.
• Gildemacher PR, Schulte-Geldermann E, Borus D, Demo P, Kinyae P, Mundia P & Struik PC (2011) Seed potato quality improvement through positive selection by smallholder farmers in Kenya. Potato Research, 54:253-266.
• Giraldo S, Sierra A, Ospina M, Higuita M, Gallo Y, Gutiérrez P & Marín M (2022) Detección y caracterización molecular del potato virus B (PVB) en papa criolla (Solanum phureja) en Antioquia. Acta Biológica Colombiana, 27:258-268.
• Gutiérrez PA, Alzate JF & Marín-Montoya MA (2013) Complete genome sequence of a novel potato virus S strain infecting Solanum phureja in Colombia. Archives of Virology, 158:2205-2208.
• Gutiérrez P, Mesa HJ & Marín M (2016) Genome sequence of a divergent Colombian isolate of potato virus V (PVV) infecting Solanum phureja. Acta Virologica, 60:49-54.
• Gutiérrez P, Rivillas A, Tejada D, Giraldo S, Restrepo A, Ospina M, Cadavid S, Gallo Y & Marín M (2021) PVDP: A portable open source pipeline for detection of plant viruses in RNAseq data. A case study on potato viruses in Antioquia (Colombia). Physiological and Molecular Plant Pathology, 113:101604.
• Guyader S & Ducray DG (2002) Sequence analysis of Potato leafroll virus isolates reveals genetic stability, major evolutionary events and differential selection pressure between overlapping reading frame products. Journal of General Virology, 83:1799-1807.
• Guzmán M, Román V, Franco L & Rodríguez P (2010) Presencia de cuatro virus en algunas accesiones de la Colección Central Colombiana de papa mantenida en campo. Agronomía Colombiana, 28:225-233.
• Guzmán-Barney M, Franco-Lara L, Rodríguez D, Vargas L & Fierro JE (2012) Yield losses in Solanum tuberosum Group Phureja cultivar criolla Colombia in plants with symptoms of PYVV in field trials. American Journal of Potato Research, 89:438-447.
• Halterman D, Charkowski A & Verchot J (2012) Potato viruses and seed certification in the USA to provide healthy propagated tubers. Pest Technology, 6:01-14.
• ICA (2015) Resolución 3168 de 2015. Available at: <https://www.ica.gov.co/getattachment/4e8c3698-8fcb-4e42-80e7-a6c7acde9bf8/2015R3168.aspx>. Accessed on: September 20^th, 2022.
  » https://www.ica.gov.co/getattachment/4e8c3698-8fcb-4e42-80e7-a6c7acde9bf8/2015R3168.aspx
• ICTV (2021) International Committee on Taxonomy of Viruses. Available at: <https://ictv.global/msl>. Accessed on: October 20^th, 2021.
  » https://ictv.global/msl
• Kreuze JF, Souza-Dias JAC, Jeevalatha A, Figueira AR, Valkonen JPT & Jones RAC (2020) Viral diseases in potato. In: Campos H & Ortiz O (Eds.) The potato crop. Cham, Springer. pp.389-430.
• Kumar R, Kaundal P, Kumar R, Siddappa S, Kumari H, Chandra K, Sharma S & Kumar M (2021) Rapid and sensitive detection of potato virus X by one-step reverse transcription-recombinase polymerase amplification method in potato leaves and dormant tubers. Molecular and Cellular Probes, 58:101743.
• MADR (2019) Estrategia de ordenamiento de la producción. Cadena productiva de la papa y su industria. Available at: <https://sioc.minagricultura.gov.co/Papa/Normatividad/Plan%20de%20Ordenamiento%20papa%202019-2023.pdf>. Accessed on: April 25^th, 2021.
  » https://sioc.minagricultura.gov.co/Papa/Normatividad/Plan%20de%20Ordenamiento%20papa%202019-2023.pdf
• Mesa ME, González MI, Gutiérrez PA & Marín MA (2016) Diagnóstico serológico y molecular del Potato leafroll virus (PLRV) en tubérculos-semilla de papa en Antioquia. Acta Agronómica, 65:204-210.
• Mumford RA, Walsh K, Barker I & Boonham N (2000) Detection of Potato mop top virus and Tobacco rattle virus Using a Multiplex Real-Time Fluorescent Reverse-Transcription Polymerase Chain Reaction Assay. Phytopathology, 90:448-453.
• Muñoz L, Gutiérrez-Sánchez PA & Marín-Montoya M (2016a) Detección y secuenciación del genoma del Potato Virus Y (PVY) que infecta plantas de tomate en Antioquia, Colombia. Bioagro, 28:69-80.
• Muñoz D, Gutiérrez P & Marín M (2016b) Detección y caracterización molecular del Potato virus Y (PVY) en cultivos de papa (Solanum tuberosum L.) del norte de Antioquia, Colombia. Protección Vegetal, 31:09-19.
• Porras P & Herrera C (2015) Modelo productivo de la papa variedad Diacol Capiro para el departamento de Antioquia. Mosquera, Corpoica. 92p.
• Price WC (1938) Studies on the virus of tobacco necrosis. American Journal of Botany, 25:603.
• Raigond B, Verma A, Pathania S, Sridhar J, Kochhar T & Chakrabarti SK (2020) Development of a reverse transcription loop-mediated isothermal amplification for detection of potato virus a in potato and in insect vector aphids. Crop Protection, 137:105296.
• Robinson JT, Thorvaldsdóttir H, Winckler W, Guttman M, Lander ES, Getz G & Mesirov JP (2011) Integrative Genomics Viewer. Nature Biotechnology, 29:24-26.
• Rodríguez LE, Ñustez CE & Estrada N (2009) Criolla Latina, Criolla Paisa y Criolla Colombia, nuevos cultivares de papa criolla para el departamento de Antioquia (Colombia). Agronomía Colombiana, 27:289-303.
• Saitou N & Nei M (1987) The neighbor-joining method: a new method for reconstructing phylogenetic trees. Molecular Biology and Evolution, 4:406-425.
• Savenkov EI, Sandgren MY & Valkonen JPT (1999) Complete sequence of RNA 1 and the presence of tRNA-like structures in all RNAs of Potato mop-top virus, genus Pomovirus Journal of General Virology, 80:2779-2784.
• Schulte-Geldermann E, Gildemacher PR & Struik PC (2012) Improving seed health and seed performance by positive selection in three Kenyan potato varieties. American Journal of Potato Research, 89:429-437.
• Schumpp O, Bréchon A, Brodard J, Dupuis B, Farinelli L, Frei P, Otten P & Pellet D (2021) Large-Scale RT-qPCR diagnostics for seed potato certification. Potato Research, 64:553-569.
• Seminario JF, Villanueva-Guevara R & Valdez-Yopla MH (2018) Rendimiento de cultivares de papa (Solanum tuberosum L.) amarillos precoces del grupo Phureja. Agronomía Mesoamericana, 29:639-653.
• Sierra A, Gallo Y, Estrada M, Gutiérrez PA & Marín M (2020) Detección molecular de seis virus de ARN en brotes de tubérculos de papa criolla (Solanum phureja) en Antioquia, Colombia. Bioagro, 32:03-14.
• Sierra A, Gallo Y, Estrada M, Gutiérrez P & Marín M (2021) Detection of four RNA viruses in commercial and informal potato seed tubers in Antioquia (Colombia). Archives of Phytopathology and Plant Protection, 54:273-294.
• Singh RP, Kurz J, Boiteau G & Bernard G (1995) Detection of potato leafroll virus in single aphids by the reverse transcription polymerase chain reaction and its potential epidemiological application. Journal of Virological Methods, 1:133-143.
• Smith KM & Bald JG (1935) A description of a necrotic virus disease affecting tobacco and other plants. Parasitology, 27:231-245.
• Tamura K (1992) Estimation of the number of nucleotide substitutions when there are strong transition-transversion and G + C-content biases. Molecular Biology and Evolution, 9:678-687.
• Tamura K, Stecher G & Kumar S (2021) MEGA11: Molecular Evolutionary Genetics Analysis Version 11. Molecular Biology and Evolution, 38:3022-3027.
• Teakle DS & Gold AH (1963) Further studies of Olpidium as a vector of tobacco necrosis virus. Virology, 19:310-315.
• Thomas-Sharma S, Abdurahman A, Ali S, Andrade-Piedra JL, Bao S, Charkowski AO, Crook D, Kadian M, Kromann P, Struik PC, Torrance L, Garrett KA & Forbes GA (2016) Seed degeneration in potato: the need for an integrated seed health strategy to mitigate the problem in developing countries. Plant Pathology, 65:03-16.
• Thomas-Sharma S, Andrade-Piedra J, Carvajal M, Hernandez JF, Jeger MJ, Jones RAC, Kromann P, Legg JP, Yuen J, Forbes GA & Garrett KA (2017) A risk assessment framework for seed degeneration: Informing an integrated seed Health strategy for vegetatively propagated crops. Phytopathology, 107:1123-1135.
• Vallejo D, Gutiérrez P & Marín M (2016) Genome characterization of a Potato virus S (PVS) variant from tuber sprouts of Solanum phureja Juz. et Buk. Agronomía Colombiana, 34:51-60.
• Yang L, Nie B, Liu J & Song B (2014) A reexamination of the effectiveness of ribavirin on eradication of viruses in potato plantlets in vitro using ELISA and quantitative RT–PCR. American Journal of Potato Research, 91:304-311.

Publication Dates

• Publication in this collection
  27 Oct 2023
• Date of issue
  2023

History

• Received
  04 May 2022
• Accepted
  28 Mar 2023