SciELO - Scientific Electronic Library Online

 
vol.72 issue1Ochenta años de la revista Facultad Nacional de Agronomía, honores y reflexionesMixture of glufosinate and atrazine for ryegrass ( Lolium multiflorum Lam.) control and its effect on seeds’ quality author indexsubject indexarticles search
Home Pagealphabetic serial listing  

Services on Demand

Journal

Article

Indicators

Related links

  • On index processCited by Google
  • Have no similar articlesSimilars in SciELO
  • On index processSimilars in Google

Share


Revista Facultad Nacional de Agronomía Medellín

Print version ISSN 0304-2847On-line version ISSN 2248-7026

Rev. Fac. Nac. Agron. Medellín vol.72 no.1 Medellín Jan./Apr. 2019

https://doi.org/10.15446/rfnam.v72n1.69438 

Artículos

Complete genome sequence of a Passion fruit yellow mosaic virus (PFYMV) isolate infecting purple passion fruit ( Passiflora edulis f. edulis )

Secuenciación del genoma completo de un aislamiento del Passion fruit yellow mosaic virus (PFYMV) que infecta gulupa ( Passiflora edulis f. edulis )

Helena Jaramillo Mesa1 

Mauricio Marín Montoya1 

Pablo Andrés Gutiérrez Sánchez1 

1 Facultad de Ciencias. Universidad Nacional de Colombia. AA 3840, Medelllín, Colombia. <mamarinm@unal.edu.co>


ABSTRACT

Purple passion fruit (Passiflora edulis f. edulis), also known as gulupa, is a vine plant of the family Passifloraceae, which in recent years has gained importance in the world fruit market due to its exotic nature and excellent organoleptic properties. Although the demand for gulupa in Colombia has increased significantly to become one of the most important fruit exports, the cultivated area has been in decline since 2009 due to the impact of plant diseases. Cucumber mosaic virus (CMV), Soybean mosaic virus (SMV) and Cowpea aphid-borne mosaic virus (CABMV) are amongst the main viruses found infecting gulupa in Colombia. To further characterize the virome of gulupa, a deep sequencing transcriptome study was performed from a producing region in eastern Antioquia. Based on the results of next-generation sequencing (NGS), we report the genome sequence of a tymovirus infecting this plant. Phylogenetic analysis revealed this virus to be a close relative of Passion fruit yellow mosaic virus (PFYMV), Cassia yellow mosaic-associated virus (CYMaV) and Calopogonium yellow vein virus (CYVV). To date, only a 1,115 nt segment comprising the RdRp-CP region of PFYMV has been reported; this sequence shares 84.79% and 95.24% identities at the nucleotide and amino acid levels with isolate PFYMV_Antioquia. Finally, RT-qPCR and Sanger sequencing using specific primers confirmed the presence of PFYMV in different purple passion fruit crops in Antioquia. The present work is the first report of a complete genome sequence of PFYMV.

Keywords: Next-generation sequencing; Passifloraceae; Tymovirus; Viral genomes

RESUMEN

El maracuyá morado (Passiflora edulis f. edulis), también conocido como gulupa, es una planta enredadera de la familia Passifloraceae, que recientemente ha tomado importancia en el mercado mundial de frutas exóticas, gracias a sus excelentes características organolépticas. Aunque la demanda de gulupa en Colombia se ha incrementado significativamente en los últimos años al ser una de las frutas que más exporta el país, su área cultivada ha venido reduciéndose desde 2009 debido al impacto negativo de diferentes enfermedades. Hasta el momento se conocía que el Cucumber mosaic virus (CMV), Soybean mosaic virus (SMV) y Cowpea aphid-borne mosaic virus (CABMV) eran los principales virus que afectaban este cultivo en Colombia. Con el fin de continuar caracterizando el viroma de la gulupa, en el presente estudio se realizó una secuenciación del transcriptoma de plantas cultivadas en el oriente de Antioquia. Con base en los resultados obtenidos para la secuenciación de nueva generación (NGS), se detectó un tymovirus infectando estas plantas. Los análisis filogenéticos indicaron que dicho virus está relacionado con: Passion fruit yellow mosaic virus (PFYMV), Cassia yellow mosaic-associated virus (CYMaV) y Calopogonium yellow vein virus (CYVV). La secuencia aquí reportada compartió niveles de identidad de 84,79% y 95,24% en nucleótidos y aminoácidos, con una secuencia de 1115 nt del segmento que comprende la region RdRP-CP del aislamiento PFYMV_Antioquia. Finalmente, mediante reacciones de RT-qPCR y secuenciación Sanger utilizando cebadores específicos, se confirmó la presencia del PFYMV en diferentes cultivos de gulupa de Antioquia. Este es el primer reporte de un genoma completo del PFYMV.

Palabras clave: Secuenciación de nueva generación; Passifloraceae; Tymovirus; Genomas virales

Purple passion fruit (Passiflora edulis f. edulis), also known as gulupa, is a vine plant of the family Passifloraceae, which in recent years has gained importance in the world fruit market due to its exotic nature, excellent organoleptic properties and high content of bioactive compounds (Dhawan et al., 2004; Fischer and Rezende, 2008). Purple passion fruit is closely related to yellow passion fruit (P. edulis f. flavicarpa) but has a sweeter taste and is less acidic due to a higher content of sugars and a lower concentration of citric acid (Morton, 1987; Manicom et al., 2003). The gulupa fruit is also a good source of iron, vitamins A, B and C, amino acids, carotenoids and flavonoids and can be eaten fresh or used in the preparation of beverages, jams, chutneys and desserts in general (Morton, 1987; Dhawan et al., 2004; Rodríguez-Amaya, 2012). The leaves of purple passion fruit are frequently used in infusions due to their sedative and anti-inflammatory properties (Morton 1987; Ferreres et al., 2007, da Silva et al., 2013); the flowers, on the other hand, are often used as ornamentals (Manicom et al., 2003).

The origin of purple passion fruit has been located to an area comprising southern Brazil, Paraguay and northern Argentina (Morton, 1987) and its cultivation is restricted to tropical and subtropical regions at altitudes between 1400-2500 MASL with temperatures in the 15 to 30 °C range (Jiménez et al., 2012). Gulupa is mainly cultivated in Colombia, Ecuador, Peru, South Africa, Kenya, Zimbabwe, India, New Zealand, Sri Lanka, Uganda, Australia and some Pacific islands (Fischer and Rezende, 2008). Total production of passiflora fruits is estimated at 805,000 t year-1, with Brazil (56%), Ecuador (31.0%), and Colombia (3.7%) as the biggest producers (Manicom et al., 2003; Fischer and Rezende, 2008; Rodríguez-Amaya, 2012).

During the past decade, the demand of gulupa in Colombia has increased significantly to become one of the most important exotic fruit exports in the country, just behind banana and Cape gooseberry (Jiménez et al., 2012). The total cultivated area of gulupa in Colombia comprises approximately 500 ha mostly concentrated in the provinces of Antioquia (36.87%), Cundinamarca (28.8%) and Boyacá (13.51%) and smaller amounts produced in Valle del Cauca, Magdalena, Huila and Santander (Guerrero et al., 2011; Agronet, 2015). In spite of the increasing demand, the cultivated area of gulupa has been in decline since 2009 due to the impact of fungal, bacterial and viral diseases. This effect has been most dramatic in the province of Cundinamarca where the total cultivated area went down from 630 ha in 2009 to 204 ha in 2013, and it has had one of the poorest yields in the country (Agronet, 2015). In 2010, a study on plants showing symptoms of viral infection in the Sumapaz region of Cundinamarca confirmed the infection of gulupa by the potyvirus Soybean mosaic virus (SMV) (Camelo, 2010). Further work, reported SMV incidences of 23.7 and 27.8% in two gulupa producing regions of Cundinamarca (Gordillo, 2011). Besides SMV, there are reports of gulupa infections by another potyvirus, Cowpea aphid-borne mosaic virus (CABMV), which was detected using ELISA in leaves showing mosaic symptoms in Cundinamarca and Valle del Cauca (Gordillo, 2011).

In spite of ranking third in the cultivated area (74 ha), Antioquia has the highest yield of purple passion fruit (31.2 t ha-1) and is also the largest producer with 36.8% of the national production (Agronet, 2015). Unfortunately, commercial gulupa cultivation plots in this region have started to be affected by an unknown disease that results in significant yield reductions. The symptoms of this pathology have been previously observed in gulupa fields in Cundinamarca and are suggestive of an emergent disease of viral origin. A Next-generation sequencing (NGS) study on gulupa fields in the municipality of Jardin demonstrated the presence of a SMV strain infecting this plant in Antioquia and reported the first complete SMV genome infecting a host different to soybean (Jaramillo, 2017). As different viruses have been found to infect P. edulis in the world (i.e., Cucumber mosaic virus, Passionfruit mottle virus, Bean yellow mosaic virus), it is likely that other viruses could be found infecting this crop in Antioquia as a result of the free movement of infected plants and the lack of seed certification programs in Colombia. To further characterize the virome infecting gulupa in Antioquia, in this study a NGS analysis was performed from a producing region in the eastern municipalities of El Retiro and La Ceja. The molecular properties of this virus are discussed and its widespread presence in this region demonstrated by RT-PCR and RT-qPCR tests.

MATERIALS AND METHODS

Next-generation Sequencing

High-throughput sequencing of the gulupa transcriptome was performed on a bulked sample of two leaves from approximately 10 plants with symptoms of viral infection such as mosaic, mottling along the veins, enation and severe leaflet deformation (Figure 1) collected at the municipalities of El Retiro (6°03’15.60” N75°30’5.39”W) and La Ceja (6°01’39.4”N, 75°25’53.8”W ) (Antioquia, Colombia). After grinding the leaf tissue with liquid nitrogen, total RNA was extracted with Trizol (Ambion, Thermo Fisher Scientific, EEUU) and rRNA depleted with the Ribo-Zero Plant kit (Illumina, EEUU). Libraries were constructed using the TruSeq RNA Sample Preparation Kit (Illumina, EEUU) and sequencing performed with the Illumina HiSeq 2500 system provided by Macrogen (South Korea) (Muñoz-Baena et al., 2017). Adapter sequences and low-quality bases were removed with Trimommatic (Bolger et al., 2014). Sequence assembly was performed with Trinity (Grabherr et al., 2011) and viral contigs identified by a local BLASTN search using a database of viral reference genomes. The database comprised a collection of 10,884 completely sequenced plant viruses representative of all the taxa and genome types currently known for this group of pathogens and were downloaded from NCBI. The PFYMV genome assembly was confirmed by mapping with BWA (Li and Durbin, 2009) using the Trinity contig as the reference. Sequencing of the purple passion fruit transcriptome resulted in a paired-end library of 21,480,984 reads (101 nt/read) for a total of 2,169,579,384 sequenced nucleotides. The consensus genome sequence was deposited in GenBank under accession KY823429 with PFYMV_Antioquia as isolate name.

Figure 1 Leaf symptoms of virus infection in P. edulis f. edulis; A. Asymptomatic leaf; B. Mottling along the veins; C. Mosaic and enation; D. Severe leaflet deformation. 

RT-PCR and RT-qPCR

To confirm the presence of PFYMV_Antioquia in eastern Antioquia, leaf samples were collected in one gulupa plot (~0.4 ha) in the municipality of La Ceja (LCS= La Ceja symptomatic sample; LCA= La Ceja asymptomatic sample) and three more (<1 ha) in El Retiro (R#S=Retiro symptomatic sample; R#A= Retiro asymptomatic sample) (Antioquia, Colombia); two bulked samples consisting of nearly 20 leaves with and without viral symptoms were randomly collected at each plot. RNA was extracted with Trizol (Ambion, Thermo Fisher Scientific, EEUU) and eluted in 40 μL of DEPC treated water; the purity and concentration were determined by absorbance readings at 260 and 280 nm using a Nanodrop 2000C (Thermo Fisher Scientific, EEUU). Retrotranscription was performed for 30 min at 50 °C in 20 μL containing 200 U of Maxima Reverse Transcriptase (Thermo Fisher Scientific, EEUU), 1X RT Buffer, 0.5 mM dNTP Mix, 100 pmol of the reverse primer, 20 U of RiboLock RNase Inhibitor and 100-500 ng of total RNA (Muñoz-Baena et al., 2017). All primers were designed in this study using the PFYMV_Antioquia genome sequence reported in this work. RT-qPCR amplification of the RNA-dependent RNA polymerase (RdRp) region was carried out with primers Tymo_F_RdRp1 (5’-CGG TAC CCT TCC ACA ACT CC-3’) and qTymo_R_RdRp1 (5’-AGG GAG CAG AAG AAT TTC G-3’), targeting a 137 bp fragment of the RdRp region; primers Tymo_F_CP (5’-AGT CTC CGA CTC AAT CAG CG-3’) and qTymo_R_CP (5’-GCG ACG AGA GAG GTG AGT CG-3’) were used for the capsid protein (CP) segment, which resulted in a fragment of 90 bp. For RT-PCR, primer Tymo_R_RdRp1 (5’-CCT GGT GTA AGT GAA CAG GGC-3’) was used in combination with Tymo_F_RdRp1 resulting in an amplicon of 576 bp and primer Tymo_R_CP (5’-GAA TGG AGT TCG ATC GTG C-3’) was used in combination with Tymo_F_CP, resulting in an amplicon of 405 bp.

For the qPCR, the Maxima SYBR Green/ROX qPCR Master Mix (2X) kit (Thermo Fisher Scientific, EEUU) was used in 25 μL of the reaction containing 12.5 μL mix, 10 μL DEPC water, primers at 0.3 μM each and 50-100 ng cDNA. Samples were considered positive if they exhibited fluorescence values higher than the threshold before the 35th cycle. Primer specificity was verified by High-Resolution Melting (HRM) in the 50 and 99 °C range. For RT-PCR, the amplification mix included 17.8 μL water, 1X enzyme buffer, 1.8 mM MgCl2, 0.2 mM dNTPs, 0.2 μM primers, 1 U Taq DNA polymerase (Thermo Fisher Scientific, EEUU) and 50-100 ng de cDNA in a final volume of 25 μL. The amplification consisted of an initial 3 min incubation at 95 °C, followed by 35 cycles that included denaturation at 94 °C (30 s), annealing at 52 °C (1 min) and extension at 72 °C (1 min); the amplification ended with an extension cycle at 72 °C for 5 min. Amplicon size was determined by 1.8% agarose gel electrophoresis stained with GelRed 1X (Biotium, EEUU) in a Bio Doc Analyze transilluminator (Biometra, Alemania). Sequences of the RT-qPCR and RT-PCR amplification products were confirmed by the Sanger method using an ABI Prism 3730xl sequencer (PE Applied Biosystems, EEUU) at Macrogen (South Korea). Sanger sequences were deposited in GenBank under accession codes KY823430-KY823445.

Bioinformatic analysis

Open reading frames (ORF) were identified by comparison against known tymoviral sequences using BLASTX with default parameters (scoring matrix: BLOSUM62, word size: 6, gap opening: 11, gap extension: 1) with the exception of E-value, which was set up with a cut-off of 1e-10 (Gish and States, 1993). Functional domains were identified by a search against the Pfam database using the Conserved Domains search at NCBI using default parameters with an expect value threshold of 0.01 (Marchler-Bauer et al., 2017). Protease cleavage sites were identified by finding matches to the sequences reported by Jakubiec et al. (2007) in Turnip yellow mosaic virus (TYMV) using an amino acid sequence alignment with the PFYMV_Antioquia RdRp protein. The 5’ hairpin loops 1 and 2 (HP1-2) and the tRNAval-like motif at the 3’ end of this tymovirus genome were predicted using the minimum free energy algorithm (Zuker and Stiegler, 1981) available in the ViennaRNA Package (Gruber et al., 2008). Evolutionary analyses were conducted in MEGA7 (Kumar et al., 2016) using the Maximum Likelihood (ML) method based on the General Reverse Time Reversible model with 1,000 bootstrap replicates (Nei and Kumar, 2000) using a discrete Gamma distribution with five categories. The evolutionary model was selected with ModelTest (Posada and Crandall, 1998). All Sequence alignments were performed with MUSCLE using gap opening and extension penalties of 1 and 0, respectively, a maximum of 16 iterations, and UPGMA as clustering method (Edgar, 2004).

RESULTS AND DISCUSSION

Next-generation Sequencing

Sequencing of the purple passion fruit transcriptome resulted in a total of 2,169,579,384 nucleotides. After sequence assembly, a contig of 6,088 nt with high sequence similarity to members of the genus Tymovirus (family Tymoviridae) such as PFYMV (84.79%, AF5467107), TYMV (73%, KJ690173), Grapevine fleck virus (GFkV) (73%, JF927942) and Watercress white vein virus (WWVV) (71%, JQ001816) (Table 1). Due to high sequence similarity with the partial RdRp (96.3%) and CP (97.0%) protein sequences of PFYMV, the viral contig from P. edulis f. edulis was named as PFYMV_Antioquia. The PFYMV_Antioquia sequence was assembled from a total of 3,823,738 paired-end reads for an average sequence depth of 123.58x and an abundance of 65,922.13 fragments per kilobase per million reads (FPKM) (Figure 2).

Table 1 Percent nucleotide (nt) and amino acid (a.a.) sequence identities of the PFYMV_Antioquia genome and proteins against selected tymovirus species representing the genus variability. Nucleotide sequence identities are shown for the complete genome and the predicted ORFs encoding the Movement protein (MP), RdRp and CP proteins. Protein comparisons are with respect to predicted amino acid sequences of MP, RdRp, and CP in PFYMV_Antioquia. Dashes represent missing sequences. 

Figure 2 Sequence coverage of the PFYMV_Antioquia genome. The contig was assembled from a total of 3,823,738 paired-end reads for an average sequence depth of 123.58x and an abundance of 65,922.13 fragments per kilobase per million reads (FPKM). The relative position of the ORFs encoding for the MP, RdRp and CP proteins are shown on top and include the location of predicted methyltransferase (MT), protease (PRO), helicase (HEL) and polymerase (POL) domains. 

Three ORFs were identified in the PFYMV_Antioquia genome (Figure 3A). ORF1 is located at nucleotide positions 112 to 5,418 and encodes for a putative replication polyprotein of 1,768 a.a. (197.6 kDa) containing methyltransferase (Mtr, pfam01660), helicase (HEL, pfam01443), RNA-dependent RNA polymerase (POL, pfam00978) and protease (PRO, pfam05381) domains. In TYMV, the protease domain mediates cleavage of the 206K polyprotein precursor between the HEL/POL and PRO/HEL domains (Jakubiec et al., 2007); in the PFYMV_Antioquia genome, the polyprotein is predicted to be cleaved at amino acid positions 802 and 1,186 with recognition sequences PKLLGS/S and PRLLGS/H, respectively. ORF2 is located at positions 105-2,006, 7 nt upstream of ORF1, and encodes a movement protein (MP) of 68.6 kDa. ORF3 was identified at positions 5,422-6,003 and encodes the capsid protein (CP, 19.8 kDa). This protein product is predicted to be translated from a subgenomic mRNA transcribed from a viral promoter, or tymobox, with sequence GAGUCUGAAUUGCUUC and identical to that found in TYMV (Schirawski et al., 2000). The tymobox is located at positions 5,383-5,404 within ORF1 and 39 nt upstream of the ORF3 initiation codon (Figure 3A).

Figure 3 Genome annotation and phylogenetic relationships of the PFYMV isolate from Purple passion fruit in Antioquia. A. Molecular features of the PFYMV_Antioquia genome. The first bar diagram represents the complete genome; nucleotide positions of HP1, HP2 and the tymobox are indicated. The diagrams below are drawn to scale and indicate the position of predicted ORFs and functional domains. Putative protease cleavage sites within the RdRp polyprotein sites are indicated as triangles; B. Maximum-likelihood (ML) phylogenetic tree using available completely sequence tymovirus genomes; C. ML tree using CP sequences. GenBank accession numbers are shown in parentheses. 

In agreement with other tymovirus genomes (Giegé et al., 1993; Bink et al., 2003) RNA secondary structure predictions suggests the formation of two hairpin sequences at positions 33-63 (HP1) and 75-98 (HP2) in the 5’ UTR and a tRNA-like structure with a valine anticodon sequence (CAC) in 3’ UTR (Figure 4).

Figure 4 Predicted RNA secondary structure elements in PFYMV_Antioquia. A. Hairpin loops HP1 and HP1. B. tRNAval-like structure at the 3&apos;UTR. Secondary structure predictions were performed using the minimum free energy algorithm package available in the ViennaRNA Package. The color scheme represents the probability for a given base to adopt the secondary structure presented in the diagram. 

Maximum likelihood analysis of complete tymoviral genomes, using the marafivirus Maize rayado fino virus (MRFV) as outgroup, confirms PFYMV_Antioquia to be a member of the genus Tymovirus with phylogenetic affinity to a clade comprising Asclepias asymptomatic virus (AsAV), Kennedya yellow mottle virus (KYMV), Okra mosaic virus (OMV) and TYMV (Figure 3B). Analysis of CP sequences confirms that PFYMV_Antioquia is more closely related to PFYMV and forms part of a subclade that also includes Cassia yellow mosaic-associated virus (CYMaV) and Calopogonium yellow vein virus (CAYVV) (Figure 3C). The only PFYMV sequence available corresponds to a 1,115 nt segment of the 3’ end obtained from a passion fruit plant with conspicuous yellow mosaic leaf symptoms in Colombia (Morales et al., 2002). It is important to note that the currently available PFYMV sequence (AF5467107) seem to contain errors in the 3’-end of the CP ORF which could mislead future work on PFYMV (Figure 5).

Figure 5 Multiple sequence alignment of the CP region of tymovirus closely related to PFYMV_Antioquia. A. Despite high amino acid sequence identity (97.0%) at the N-terminus (1-143) with PFYMV isolate (AF467107), premature codons are observed at the C-terminus suggesting errors in the deposited sequence (in red); B. Removal of the previously sequenced PFYMV segment (AF467107) reveals that the PFYMV_Antioquia is in good agreement with other previously sequenced tymovirus. 

RT-PCR and RT-qPCR

The presence of PFYMV_Antioquia was confirmed in composite samples consisting of either symptomatic (R1S, R2S, R3S, and LCS) or asymptomatic (R1A, R2A, R3A, and LCA) leaves from four purple passion fruit plots in Antioquia using RT-qPCR and RT-PCR with primers targeting RdRp and CP regions (Figure 6). All samples tested positive for this tymovirus. In RT-qPCR, Tymo_F_RdRp1/qTymo_R_RdRp1 reactions exhibited Ct values between 4.16-18.57 and a narrow Tm range (83.6±0.15 °C), while reactions for Tymo_F_CP/qTymo_R_CP had Ct values between 2.33 and 17.72 and Tm of 86.6±0.15 °C, suggesting low sequence variability (Figure 7). Sanger Sequencing of the RT-PCR amplicons and phylogenetic analysis confirmed the presence of the PFYMV (Figure 7).

Figure 6 Maximum-likelihood phylogenetic trees of RT-PCR amplicons from PFYMV isolates infecting commercial Purple passion fruit plots. A. ML tree derived from RT-PCR amplification and sequencing of a RdRp segment using primers Tymo_F_RdRp1/Tymo_R_RdRp1. B. ML tree derived from RT-PCR amplification and sequencing of a CP segment using primers Tymo_F_CP/Tymo_R_CP. 

Figure 7 Real-time PCR detection of PFYMV_Antioquia in field samples. Amplifications were performed on symptomatic (R1S, R2S, R3S, and LCS) and asymptomatic (R1A, R2A, R3A, and LCA) bulked leaf samples from commercial P. edulis f. edulis plots in Antioquia with primers targeting either A. the RdRp or B. CP regions. Primers were designed based on the PFYMV genome presented in this work. Bar diagrams on top of each panel represent the relative position of each primer within the PFYMV genome. Panels to the right and left correspond to amplification curves and melting temperature profiles, respectively 

This work confirms the presence of a tymovirus infecting purple passion fruit and phylogenetically related to the PFYMV isolate infecting the yellow passion fruit (P. edulis f. flavicarpa). Before this work, the only available PFYMV sequence was obtained from passion fruit plants with conspicuous yellow mosaic leaf symptoms in Colombia. Partial genome sequencing showed this isolate to be distinct from other tymoviruses but with similar pathogenic and antigenic properties to PFYMV infecting Passiflora spp. (Morales et al., 2002). The Brazilian PFYMV isolate was identified in yellow passion fruit plants with symptoms of yellow net, yellow mosaic, and leaf crinkle at Papucaia, in the State of Rio de Janeiro (Crestani et al., 1986). The Brazilian PFYMV isolate has a host range limited to species of the genus Passiflora and can be transmitted by graft or mechanical inoculation and by the chrysomelid beetle Diabrotica speciosa (Crestani et al., 1986). The Colombian PFYMV isolate, on the other hand, was first detected in the passion fruit-producing regions of Antioquia, Caldas, Santander and Valle del Cauca as an isometric virus with serological relationship to Desmodium yellow mottle virus with particle morphology, which induced symptoms similar to tymoviruses (Varón de Agudelo et al., 1992). In contrast to the Brazilian PFYMV, inoculation of the Colombian PFYMV isolate can induce symptoms in other Passiflora species such as Passiflora adenopoda, P. edulis, P. pinnatistipula, P. cuadrangularis as well as solanaceous plants such as Physalis angulata, P. floridiana and P. peruviana (Morales et al., 2002). However, the Colombian PFYMV isolate failed to be transmitted by chrysomelid bettles of the genus Cerotoma, Colaspis and Diabrotica (Morales et al., 2002).

This report, on the complete genome sequence of a PFYMV isolate infecting purple passion fruit, is a new addition to the growing number of NGS sequencing reports on viruses affecting Passiflora species in South America. For example, a recent work by Vidal et al. (2018) showed that blisters, mosaics and leaf deformation symptoms observed in passion fruit crops in Bahia (Brazil) was due to a mixed infection with CABMV and Cucurbit aphid-borne yellows virus (CABYV) (Polerovirus, Luteoviridae). This work is also the first record of CABYV in passion fruit in the world (Vidal et al., 2018). Fontenele et al. (2018), on the other hand, discovered a divergent geminivirus co-infecting a diseased passion fruit plant together with CABMV in the state of Mato Grosso do Sul (Brazil) using NGS sequencing of circular DNA. This new geminivirus has been tentatively named as Passion fruit chlorotic mottle virus (PCMoV) and is the third ssDNA virus affecting Passiflora species together with Passion fruit severe leaf distortion virus (PSLDV), identified infecting this plant in Brazil and Passion fruit leaf distortion virus (PLDV), a proposed species found on Passiflora edulis in the province of Valle del Cauca (Colombia) (Vaca-Vaca et al., 2017). Finally, Jaramillo et al. (2018) reported the complete genome sequence of a SMV strain infecting purple passion fruit in south eastern Antioquia (Colombia) and associated with variegation, rugose mosaics and dark green islands symptoms using RNAseq of symptomatic leaves.

We have given evidence for the existence of a new PFYMV isolate within the genus Tymovirus. The species demarcation criteria within this genus are: overall genome sequence identity of less than 80%, capsid protein sequences less than 90% identical, differences in the 3’-terminal structure, differential host range and serological specificity (King et al., 2012). At the genome level, the PFYMV_Antioquia shares between 54.8% to 63.8% nucleotide sequence identities with available complete genome sequences of tymovirus being most similar to AsAV (62.5%), KYMV (63.8%) and OkMV (63.2%). A comparison of sequences of the capsid protein reveals that the more closely related species are CYMaV and CYVV with 66.4% and 66.3% amino acid sequence identity, respectively. Tymoviruses have a very restrictive host range and can be transmitted mechanically or by beetles of the families Chrysomelidae and Curculionidae in a semi-persistent manner (King et al., 2012).

CONCLUSIONS

Using Next-generation sequencing of the gulupa transcriptome from a bulk sample of leaves from eastern Antioquia, a new isolate of Passion fruit yellow mosaic virus (PFYMV) was identified and its complete genome sequence of 6,088 nt was obtained. To our knowledge, this represents the first complete PFYMV genome in the world. The three ORFs of tymoviruses were identified: ORF1 encodes a replication polyprotein of 1,768 a.a. (197.6 kDa), ORF2 encodes the movement protein (MP), it is located 7 nt upstream of ORF1 and has a predicted molecular weight of 68.6 kDa (633 a.a.) and finally, ORF3 encodes a coat protein (CP) of 19.8 kDa.

Two primer sets targeting the RdRp polyprotein and CP were designed using the PFYMV_Antioquia genome as the reference to confirm the presence of PFYMV by RT-qPCR. Primers were validated using the NGS sample as a positive control and tested in composite samples of both symptomatic and asymptomatic leaves from gulupa plots. Except for the negative control, all samples tested positive for PFYMV. Further work should investigate the host range and vector transmission of PFYMV as well as its serological specificity.

ACKNOWLEDGEMENTS

This work was funded by Universidad Nacional de Colombia through Grant VRI: 34579.

REFERENCES

Agronet. 2015. Área, producción y rendimiento de Gulupa en Colombia. In: Agronet, In: Agronet, http://www.agronet.gov.co/estadistica/Paginas/default.aspx ; accessed: November 2016. [ Links ]

Bink HH, Schirawski J, Haenni AL and Pleij CW. 2003. The 5’-proximal hairpin of turnip yellow mosaic virus RNA: its role in translation and encapsidation. Journal of Virology 77(13): 7452-7458. doi: 10.1128/JVI.77.13.7452-7458.2003 [ Links ]

Bolger AM, Lohse M and Usadel B. 2014. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30(15): 2114-2120. doi: 10.1093/bioinformatics/btu170 [ Links ]

Camelo VM. 2010. Detección e identificación de los virus patógenos de cultivos de gulupa (Passiflora edulis Sims) en la región del Sumapaz (Cundinamarca). Master’s Thesis in Phytopathology. Faculty of Agronomy. Universidad Nacional de Colombia. Bogotá. 60 p. [ Links ]

Crestani OA, Kitajima EW, Lin MT and Marinho VLA. 1986. Passion fruit yellow mosaic virus, a new tymovirus found in Brazil. Phytopathology 76: 951-955. doi: 10.1094/Phyto-76-951 [ Links ]

da Silva JK, Cazarin CBB, Colomeu TC, Batista AG, Meletti LMM, Paschoal JAR, Bogusz Júnior S, Furlan MF, Reyes FGR, Augusto F, Maróstica M and Zollnerb R. 2013. Antioxidant activity of aqueous extract of passion fruit (Passiflora edulis) leaves: In vitro and in vivo study. Food Research International 53: 882-890. doi: 10.1016/j.foodres.2012.12.043 [ Links ]

Dhawan K, Dhawan S and Sharma A. 2004. Passiflora: a review update. Journal of Ethnopharmacology 94(1): 1-23. doi: 10.1016/j.jep.2004.02.023 [ Links ]

Edgar RC. 2004. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Research 32(5): 1792-1797. doi: 10.1093/nar/gkh340 [ Links ]

Ferreres F, Sousa C, Valentão P, Andrade PB, Seabra RM and Gil-Izquierdo A. 2007. New C-deoxyhexosyl flavones and antioxidant properties of Passiflora edulis leaf extract. Journal of Agriculture and Food Chemistry 55(25): 10187-10193. doi: 10.1021/jf072119y [ Links ]

Fischer IH and Rezende JAM. 2008. Diseases of Passion Flower (Passiflora spp.). Pest Technology 2(1): 1-19. [ Links ]

Fontenele RS, Abreu RA, Lamas NS, Alves-Freitas DMT, Vidal AH, Poppiel RR, Melo FL, Lacorte C, Martin DP, Campos MA, Varsani A and Ribeiro SG. 2018. Passion fruit chlorotic mottle virus: molecular characterization of a new divergent geminivirus in Brazil. Viruses 10(4): E169. doi: 10.3390/v10040169. [ Links ]

Giegé R, Florentz C and Dreher TW. 1993. The TYMV tRNA-like structure. Biochimie 75(7): 569-582. doi: 10.1016/0300-9084(93)90063-X [ Links ]

Gish W and States DJ. 1993. Identification of protein coding regions by database similarity search. Nature Genetics 3: 266-272. doi: 10.1038/ng0393-266 [ Links ]

Gordillo LA. 2011. Incidencia del Soybean mosaic virus en cultivos de Gulupa (Passiflora edulis Sims) en Cundinamarca y estudio de su diversidad en Colombia. Master’s Thesis in Microbiology. Faculty of Sciences. Universidad Nacional de Colombia. Bogotá. 106 p. [ Links ]

Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng Q, Chen Z, Mauceli E, Hacohen N, Gnirke A, Rhind N, di Palma F, Birren BW, Nusbaum C, Lindblad-Toh K, Friedman N and Regev A. 2011. Full-length transcriptome assembly from RNA-seq data without a reference genome. Nature Biotechnology 29: 644-652. doi: 10.1038/nbt.1883 [ Links ]

Guerrero E, Potosí Guampe C, Melgarejo LM and Hoyos-Carvajal L. 2011. Chapter 7: Manejo agronómico de gulupa (Passiflora edulis Sims) en el marco de las Buenas Prácticas Agrícolas (BPA). pp 123-144. In: Melgarejo LM (ed). Ecofisiología del cultivo de la gulupa (Passiflora Edulis Sims). Produmedios, Bogotá. 144p. [ Links ]

Gruber AR, Lorenz R, Bernhart SH, Neuböck R and Hofacker IL. 2008. The Vienna RNA Websuite. Nucleic Acids Research 36 Suppl. 2: W70-W74. doi: 10.1093/nar/gkn188 [ Links ]

King AMQ. 2012. Virus Taxonomy. Ninth Report of the International Committee on Taxonomy of Viruses . Academic Press, San Diego. pp. 946-948. [ Links ]

Jakubiec A, Drugeon G, Camborde L and Jupin I. 2007. Proteolytic processing of turnip yellow mosaic virus replication proteins and functional impact on infectivity. Journal of Virology 81(20): 11402-11412. doi: 10.1128/JVId.01428-07 [ Links ]

Jaramillo H. 2017. Análisis del transcriptoma y viroma de Passiflora edulis f. edulis en cultivos de Antioquia utilizando métodos de secuenciación de nueva generación. Master’s Thesis in Science. Faculty of Sciences. Universidad Nacional de Colombia. Medellín. 255 p. [ Links ]

Jaramillo H, Marín M and Gutiérrez P. 2018. Molecular characterization of Soybean mosaic virus (SMV) infecting Purple passion fruit (Passiflora edulis f. edulis) in Antioquia, Colombia. Archives of Phytopathology and Plant Protection. In press.Links ]

Jiménez Y, Carranza C and Rodríguez M. 2012. Gulupa (Passiflora edulis Sims). pp 579-599. In: Fischer G (ed). Manual para el cultivo de frutales en el trópico. Produmedios, Bogotá . [ Links ]

Kumar S, Stecher G and Tamura K. 2016. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Molecular Biology and Evolution 33(7): 1870-1874. doi: 10.1093/molbev/msw054 [ Links ]

Li H and Durbin R. 2009. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25(14): 1754-1760. doi: 10.1093/bioinformatics/btp324 [ Links ]

Manicom B, Ruggiero C, Ploetz RC and Goes A. 2003. Diseases of passion fruit. pp: 413-441. In: Ploetz RC (Ed). Diseases of tropical fruit crops. CAB international, London. [ Links ]

Marchler-Bauer A, Bo Y, Han L, He J, Lanczycki CJ, Lu S, Chitsaz F, Derbyshire MK, Geer RC, Gonzales NR, Gwadz M, Hurwitz DI, Lu F, Marchler GH, Song JS, Thanki N, Wang Z, Yamashita RA, Zhang D, Zheng C, Geer LY and Bryant SH. 2017. CDD/SPARCLE: functional classification of proteins via subfamily domain architectures. Nucleic Acids Research 45(D1): D200-D203. doi: 10.1093/nar/gkw1129 [ Links ]

Morton, J. 1987. Passionfruit. pp. 320-328. In: Morton JF (ed). Fruits of warm climates. Miami. [ Links ]

Morales FJ, Lozano I, Castaño M, Arroyave J, Velasco AC and Varon F. 2002. Partial characterization of a tymovirus infecting passion fruit in Colombia, South America. Journal of Phytopathology 150(4-5): 292-296. doi: 10.1046/j.1439-0434.2002.00740.x [ Links ]

Muñoz-Baena L, Marín-Montoya M and Gutiérrez P. 2017. Genome sequencing of two Bell pepper endornavirus (BPEV) variants infecting Capsicum annuum in Colombia. Agronomía Colombiana 35(1): 44-53. doi: 10.15446/agron.colomb.v35n1.60626. [ Links ]

Nei M and Kumar S. 2000. Molecular evolution and phylogenetics. Oxford University Press, Oxford. 333 p. [ Links ]

Posada D and Crandall KA. 1998. MODELTEST: testing the model of DNA substitution. Bioinformatics 14(9): 817-818. doi: 10.1093/bioinformatics/14.9.817 [ Links ]

Rodríguez-Amaya DB. 2012. Passionfruit. pp 321-332. In: Siddiq M (ed). Tropical and Subtropical Fruits. Postharvest, Physiology, Processing and Packaging. Wiley-Blackwell, New Jersey. 648p. [ Links ]

Schirawski J, Voyatzakis A, Zaccomer B, Bernardi F and Haenni AL. 2000. Identification and functional analysis of the turnip yellow mosaic tymovirus subgenomic promoter. Journal of Virology 74(23): 11073-11080. doi: 10.1128/JVI.74.23.11073-11080.2000 [ Links ]

Vaca-Vaca JC, Carrasco-Lozano EC and López-López K. 2017. Molecular identification of a new begomovirus infecting yellow passion fruit (Passiflora edulis) in Colombia. Archives of Virology 162(2): 573-576. doi: 10.1007/s00705-016-3098-y [ Links ]

Varón de Agudelo F, Castaño M, Arroyave JA, Velasco AC, Vuillaume C and Morales FJ. 1992. Complejo viral que afecta plantaciones de maracuyá (Passiflora edulis Sims) en el Valle del Cauca. Fruits 47: 321-328. [ Links ]

Vidal AH, Sanches MM, Alves-Freitas DMT, Abreu EFM, Lacorte C, Pinheiro-Lima B, Rosa RCC, Jesus ON, Campos MA, Varsani A and Ribeiro SG. 2018. First World Report of Cucurbit aphid-borne yellows virus infecting Passionfruit. Plant Disease, In press. doi: 10.1094/PDIS-04-18-0694-PDN. [ Links ]

Zuker M and Stiegler P. 1981. Optimal computer folding of large RNA sequences using thermodynamics and auxiliary information. Nucleic Acids Research 9(1): 133-148. [ Links ]

Received: December 14, 2017; Accepted: September 28, 2018

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License