The Horticulture Journal 89 (4): 410–416. 2020. e Japanese Society for doi: 10.2503/hortj.UTD-175 JSHS Horticultural Science http://www.jshs.jp/

Molecular Characterization of Associated with Yellow Leaf Curl Disease in Solanaceae and Cucurbitaceae Crops from Northern Sumatra, Indonesia

Elly Kesumawati1,2*, Shoko Okabe3, Munawar Khalil2, Gian Alfan2, Putra Bahagia1, Nadya Pohan1, Sabaruddin Zakaria1,2 and Sota Koeda3

1Faculty of Agriculture, Syiah Kuala University, Darussalam, Banda Aceh 23111, Indonesia 2Graduate School of Agriculture, Syiah Kuala University, Darussalam, Banda Aceh 23111, Indonesia 3Graduate School of Agriculture, Kindai University, Nara 631-8505, Japan

Begomoviruses, transmitted by whiteflies (Bemisia tabaci), have emerged as serious constraints to the cultivation of a wide variety of vegetable crops worldwide. Leaf samples from Solanaceae (tomato, tobacco, and eggplant) and Cucurbitaceae (cucumber and squash) plants exhibiting typical begomoviral yellowing and/or curling symptoms were collected in Northern Sumatra, Aceh province, Indonesia. Rolling circle amplification was conducted using DNA isolated from cucumber, squash, eggplant, and tobacco, and the full- length sequences of the begomoviruses were evaluated. The following were isolated: bipartite begomoviruses Tomato leaf curl New Delhi (ToLCNDV), Squash leaf curl China virus (SLCCNV), Tomato yellow leaf curl Kanchanaburi virus (TYLCKaV), and a monopartite Ageratum yellow vein virus (AYVV). Begomovirus diagnosis was conducted by PCR using begomovirus species-specific primers for Pepper yellow leaf curl Indonesia virus (PepYLCIV), Pepper yellow leaf curl Aceh virus (PepYLCAV), ToLCNDV, SLCCNV, TYLCKaV, and AYVV, which are the predominant begomoviruses. The primary begomovirus species for each plant were as follows: PepYLCAV for tomato, AYVV for tobacco, TYLCKaV for eggplant, ToLCNDV for cucumber, and SLCCNV for squash. This study provides valuable information for breeding begomovirus-resistant cultivars as horticultural crops.

Key Words: diversity, geminivirus, horticultural crop, host plant, Southeast Asia.

sweet pepper ( spp.), 0.963 million tons of Introduction tomato (Solanum lycopersicum L.), 0.535 million tons Indonesia comprises more than seventeen thousand of eggplant (Solanum melongena L.), 0.152 million tons islands, and the total land area is 1,904,569 km2, which of tobacco (Nicotiana tabacum L.), 0.425 million tons is approximately five times larger than the Japanese of cucumber (Cucumis sativus L.), and 0.567 million land area. The current population is > 250 million peo‐ tons of squash (Cucurbita spp.) in 2017. Indonesia cur‐ ple, and it is the 4th most populous country in the world. rently has an expansive domestic market for horticultur‐ According to statistics from the Food and Agriculture al crops, and growth will continue along with economic Organization of the United Nations (FAO STAT), development. Indonesia produced 2.359 million tons of fresh chili/ Begomoviruses are transmitted by whiteflies in the Bemisia tabaci species complex. These viruses have emerged as serious constraints to the cultivation of a Received; December 9, 2019. Accepted; January 28, 2020. First Published Online in J-STAGE on April 23, 2020. wide variety of vegetable crops worldwide (Navas- This study was supported by the Directorate General of Resources for Castillo et al., 2011). The genus Begomovirus encom‐ Science Technology and Higher Education (DG-RSTHE) and the passes 409 viral species (ICTV, 2018); most Japan Society for the Promotion of Science (JSPS) under the DG- begomoviruses have bipartite genomes comprising two RSTHE-JSPS Joint Research Program to E. Kesumawati, S. Zakaria, circular DNA components of about 2,800 nucleotides and S. Koeda, and JSPS KAKENHI Grant Number 19H02950 to S. Koeda. (DNAs A and B), whereas monopartite begomoviruses * Corresponding author (E-mail: [email protected]). generally have only DNA-A like as a genome

© 2020 The Japanese Society for Horticultural Science (JSHS), All rights reserved. Hort. J. 89 (4): 410–416. 2020. 411

(Saunders et al., 2000). In some monopartite begomo‐ cumber and squash). However, the begomovirus species viruses, alphasatellite and betasatellite, formerly known infecting eggplant, cucumber, and squash have not yet as DNA-α and DNA-β, are also present (Hanley- been clarified. Although mixed infections of multiple Bowdoin et al., 2013). Tomato yellow leaf curl virus begomovirus species are frequently observed in chili (TYLCV), a monopartite begomovirus originating from peppers (Koeda et al., 2016), mixed infections in other the Eastern Mediterranean, is distributed around the plant species have not yet been examined. In this study, world and is seriously impacting tomato production molecular identification of the begomoviruses infecting (Mabvakure et al., 2016). TYLCV is currently ranked horticultural crops was conducted, and the diagnosis 3rd among the most economically and scientifically im‐ was made by PCR using newly developed begomovirus portant plant viruses worldwide (Scholthof et al., 2011). species-specific primers. However, in Southeast and East Asia, a wide variety of Materials and Methods distinct local begomovirus species have been identified from tomato plants, whereas TYLCV has only been de‐ Virus sources, cloning, and sequence determination tected in Japan, China, and Korea (Kenyon et al., 2014). Leaf samples of tomato (S. lycopersicum), tobacco In our previous study, geminivirus diagnosis was (N. tabacum), eggplant (S. melongena), cucumber made by PCR using DNA extracted from chili pepper (C. sativus), and squash (C. maxima) plants exhibiting leaf samples collected in the fields at five local farms yellowing and/or curling symptoms were collected at located in the suburbs of Banda Aceh, Aceh province, Aceh province, Indonesia from 2012 to 2017 (Table 1; Indonesia in 2012 (Koeda et al., 2016). Pepper yellow Fig. 1). Genomic DNA was extracted from the leaves leaf curl disease (PepYLCD) was observed in more using Nucleon PhytoPure (GE Healthcare, IL, USA). than 81% of plants from all of the fields studied, and One cucumber, one squash, two eggplant, and two to‐ symptoms reached 100% at four out of five fields. Most bacco samples were randomly selected for isolation of Table 1. Samples collected in this study. of the chili pepper plants were infected by a mixture of viruses such as Pepper yellow leaf curl Indonesia virus (PepYLCIV) with another bipartite begomovirus, Table 1. Samples collected in this study. Tomato yellow leaf curl Kanchanaburi virus Plant Collection Place of Number of (TYLCKaV) and/or a monopartite begomovirus species date sampling collected samples Ageratum yellow vein virus (AYVV). More recently, we Eggplant 2017 Feb Banda Aceh 2 identified a new bipartite begomovirus species, Pepper 2017 Feb Lambeugak 3 yellow leaf curl Aceh virus (PepYLCAV), that had in‐ 2017 Oct Saree 11 fected chili pepper, tomato, and tobacco from samples Tomato 2017 Feb Saree 14 collected in 2017 (Kesumawati et al., 2019). 2017 Jul Saree 29 From our preliminary study conducted in Aceh Tobacco 2016 Jul Lambeugak 3 province, yellow leaf curl disease was observed on 2017 Jul Lambeugak 24 many Solanaceae horticultural crops (chili pepper, Cucumber 2012 Oct Lambeugak 4 tomato, tobacco, and eggplant) and Cucurbitaceae (cu‐ Squash 2017 Jul Saree 5

a b c

d e f

Fig. 1. Symptoms of yellow leaf curl disease. (a) Chili pepper (C. annuum), (b) tobacco (N. tabacum), (c) tomato (S. lycopersicum), (d) eggplant (S. melongena), (e) squash (C. maxima), and (f) cucumber (C. sativus). Arrowheads indicate symptoms observed in each plant species. 412 E. Kesumawati, S. Okabe, M. Khalil, G. Alfan, P. Bahagia, N. Pohan, S. Zakaria and S. Koeda their viral-genome DNA. Full-length begomoviral Multiple sequences were aligned using MUSCLE genomes were amplified from the extracted DNA by (Edgar, 2004). Pairwise sequence identity comparisons rolling circle amplification (RCA) using a TempliPhi were performed using the Species Demarcation Tool DNA amplification kit (GE Healthcare). An RCA reac‐ v1.2 program (Muhire et al., 2014). Phylogenetic analy‐ tion was performed according to Koeda et al. (2015). ses were conducted on MEGA7 using the Neighbor- The concatemers produced in the reaction were joining method (Kumar et al., 2016). The bootstrap monomerized by restriction digestion with BamHI, consensus of the tree was inferred from 1000 replicates. HindIII, or PstI. Digested products were resolved by electrophoresis using 1% agarose gel, and the bands Begomovirus species-specific diagnosis by PCR that corresponded to 2.8 kbp were purified using a Gel DNA-A of PepYLCIV, PepYLCAV, and TYLCKaV; DNA Recovery Kit (ZYMO Research, CA, USA). The DNA-A like of AYVV; and DNA-B of Tomato leaf curl 2.8 kbp monomers were cloned into the appropriate site New Delhi virus (ToLCNDV) and Squash leaf curl of a pBlueScript II SK (+) vector (Agilent Technology, China virus (SLCCNV) were detected using species- CA, USA). Monomeric full-length clones were purified specific primers (Table 2). PCR was performed using using a FastGeneTM Plasmid Mini Kit (Nippon EmeraldAmp PCR Master Mix (Takara Bio). The PCR Genetics, Tokyo, Japan). Presence of an alphasatellite reaction mixtures were initially denatured at 94°C for and betasatellite was examined by PCR using universal 2 min, followed by 35 or 40 cycles at 94°C for 30 s, an‐ primers (Briddon et al., 2002; Bull et al., 2003). The nealing temperature for 30 s, as listed in Table 2, and full-length sequences for betasatellites were amplified 72°C for 1 min, terminating with 3 min of final exten‐ by specific primers (Beta full F, R) (Table 2). PCR was sion at 72°C. The amplified PCR products were sub‐ performed using EmeraldAmp PCR Master Mix jected to electrophoresis, using 1.0% (w/v) agarose gel. (Takara Bio, Shiga, Japan); the reaction mixtures were Results and Discussion initially denatured at 94°C for 2 min, followed by 35 cycles at 94°C for 30 s, 60°C for 30 s, and 72°C for Leaf samples of Solanaceae and Cucurbitaceae 2 min, terminating with 3 min of final extension at plants with typical begomoviral yellowing and/or curl‐ 72°C. The amplified PCR products were cloned into a ing symptoms were collected at Banda Aceh, pTAC-1 vector (BioDynamics Laboratory, Tokyo, Lambeugak, Laweung, and Saree in Northern Sumatra, Japan). Nucleotide sequencing was carried out in an Aceh province, Indonesia, between 2012 and 2017. Dis‐ ABI PRISM 3100 genetic analyzer with an ABI PRISM tinct yellowing and curling were observed in chili pep‐ BigDye Terminator v3.1 Cycle Sequencing Kit per and tobacco plants (Fig. 1a, b), whereas the (Applied Biosystems, CA, USA). symptoms observed in the tomato and eggplant were generally slight (Fig. 1c, d). In squash and cucumber, Sequence analysis yellowing was observed in the leaves (Fig. 1e, f). The The obtained sequences were assembled using ATGC same symptoms were observed in all of the other plants software (Genetyx, Tokyo, Japan) and analyzed by con‐ cultivated in these fields. ducting sequence similarity searches using BLASTn. In our previous studies, full-length sequences of Table 2. Begomovirus species-specific primers used in this study.

Table 2. Begomovirus species-specific primers used in this study.

Viurs species Fragment Annealing Cycle Primer name Primer sequence (5′–3′) Reference specificity size (bp) temperature number ToLCMYB Beta full F ACCACTACGCTACGCAGCAGCC 1345 60 35 This study Beta full R TTCTACCCTCCCAGGGGTACAC PepYLCIV Pep uni F GTGYWGTAYCTTCTGYGGAAYTKGA 666 58 35 Kesumawati et al., 2019 PI uni R ACGCCGTAAACGATGTTTAYGCG PepYLCAV Pep uni F GTGYWGTAYCTTCTGYGGAAYTKGA 466 58 35 Kesumawati et al., 2019 PA uni R GAACTTCAAGCGTGGGTGTC TYLCKaV TYLCKaV uni F TACATAATCAGCTGCGCGTATTACA 831 65 35 This study TYLCKaV uni R CATAGGAATKGTRATTGAAGGTGAATCT AYVV AYVV uni F GTGGTARTGGGCTTASMCATC 449 63 35 This study AYVV uni R GATYTTCARAGTAGCATACACWGGATTACT ToLCNDV ToLCNDV Uni F CKWTTKCTTCTMYTTATYCKAC 353 50 40 This study ToLCNDV Uni R GGCAATSGTYCCDGATAARKTC SLCCNV SLCCNV Uni F YACGRTCTTSMMTARTCGTTGC 396 60 35 This study SLCCNV Uni R AATGTTCCRGACATGGTAAGWC

K = G or T; M = A or C; R = A or G; S = C or G; W = A or T; Y = C or T Hort. J. 89 (4): 410–416. 2020. 413 begomovirus were isolated from chili pepper, tomato, Phylogenetic analysis was conducted using full- and tobacco (Koeda et al., 2016, 2017, 2018; length DNA-A sequences obtained in the present and Kesumawati et al., 2019). In the present study, full- previous studies (Koeda et al., 2016; Kesumawati et al., length sequences of the infecting begomoviruses affect‐ 2019) and with begomoviral sequences obtained in ing eggplant, cucumber, squash, and tobacco were other Asian countries (Fig. 2). AYVV is distributed isolated by an RCA method (Table 3). TYLCKaV se‐ widely from Southeast to East Asia. AYVV isolated in quences were isolated from two samples of eggplant our studies showed similarity with a Malaysian isolate and a single sample from tobacco, which showed high that was genetically distant from other isolates from similarity (above 99%) with TYLCKaV-[BA_B6] China, Taiwan, and Japan (Ishigaki island). Phylogenet‐ (DNA-A: LC051116, DNA-B: LC177332) isolated ic analysis implied that ToLCNDV and SLCCNV origi‐ from chili pepper showing PepYLCD in 2012 at Aceh. nated from India, and isolates from this study showed The AYVV sequence isolated from a tobacco sample the highest similarity with those found in Southeast showed high similarity (above 99%) with AYVV- Asian countries such as Thailand and Malaysia. [BA_D1-2] (DNA-A like: LC051119) isolated from TYLCKaV showed the highest similarity with these In‐ chili pepper. Furthermore, a betasatellite was isolated donesian isolates, and the phylogenetic relationship from the same tobacco sample and showed high simi‐ suggested that these isolates invaded from Thailand, larity with a Tomato leaf curl Malaysia betasatellite which is consistent with the report of Kenyon et al. (ToLCMYB) (LC511782) from a chili pepper that was (2014). TYLCKaV sequences from Vietnam and isolated in our previous study (unpublished data). How‐ Cambodia formed a different clade from those found in ever, infection with an alphasatellite was not observed Thailand, Laos, China (Yunnan province), and in the tobacco sample. These results suggested that Indonesia, indicating movement toward differentiation TYLCKaV and AYVV are capable of infecting various in Southeast Asia. species of Solanaceae crops. Wilisiani et al. (2019) detected PepYLCIV, ToLCNDV sequences were isolated from a cucumber TYLCKaV, and ToLCNDV from Solanaceae and sample. DNA-A exhibited the highest similarity with Cucurbitaceae crops with a loop-mediated isothermal sequence (AB330079) isolated from cucumber in amplification assay. In the present study conducted in Thailand in 1996 and DNA-B with sequence Aceh province, Indonesia, in addition to the ones men‐ (AB613826) isolated from cucumber in Indonesia in tioned before, PepYLCAV, SLCCNV, and AYVV were 2008. SLCCNV sequences were isolated from a squash the predominant begomoviruses for horticultural crops. sample, and DNA-A showed the highest similarity with Begomovirus species-specific primers were developed sequence (KY652743) isolated from a squash in East according to the sequences obtained in our studies and Timor in 2015, and DNA-B with sequence similar sequences registered in GenBank. As shown in (HM566113) isolated from melon in China in 2010. the phylogenetic tree (Fig. 2), similarity was high be‐ DNA-A and -B of ToLCNDV and SLCCNV showed tween PepYLCIV, PepYLCAV, and TYLCKaV, and 85% and 71% similarity, respectively, to each other, was also high between ToLCNDV and SLCCNV. For which indicates that genetically similar begomovirus ToLCNDV and SLCCNV, primer pairs to specifically species are also infective to Cucurbitaceae crops, even detect DNA-As of each begomovirus species were when genetically distant from the begomoviruses in‐ difficult to design because of their highly similar Tablefecting 3. Identitiesthe Solanaceae and sizes ofcrops. begomoviruses isolated in this study. sequences. Thus, primer pairs to detect DNA-Bs were

Table 3. Identities and sizes of begomoviruses isolated in this study.

Begomovirus DNA molecule Identity Isolate Accession no. Host species and size (nt) (%) Virus Accession no. TYLCKaV [BAEg-1] DNA A (2752) LC511771 S. melongena 99.7 TYLCKaV LC051116 DNA B (2752) LC511777 98.8 LC177332 [BAEg-49] DNA A (2752) LC511772 S. melongena 99.7 TYLCKaV LC051116 DNA B (2752) LC511778 99.0 LC177332 [BATa-20] DNA A (2752) LC511773 N. tabacum 99.5 TYLCKaV LC051116 DNA B (2752) LC511779 99.3 LC177332 AYVV [BATa-9] DNA A (2750) LC511774 N. tabacum 99.4 AYVV LC051119 ToLCMYB [BATa-9] betasatellite (1345) LC511783 N. tabacum 94.9 ToLCMYB LC511782 ToLCNDV [BACu-20] DNA A (2739) LC511775 C. sativus 97.5 ToLCNDV AB330079 DNA B (2680) LC511780 92.7 AB613826 SLCCNV [BASq-17] DNA A (2736) LC511776 C. maxima 95.0 SLCCNV KY652743 DNA B (2713) LC511781 87.6 HM566113 414 E. Kesumawati, S. Okabe, M. Khalil, G. Alfan, P. Bahagia, N. Pohan, S. Zakaria and S. Koeda

87 AYVV_LC051118_IDN * designed. AYVV had rather low similarity, less than AYVV_LC051119_IDN * 75%, to other bipartite begomoviruses. DNA samples AYVV_LC051120_IDN * AYVV-[BATa-9]_IDN *** of tobacco, eggplant, cucumber, and squash used for an 100 AYVV_LC051117_IDN * RCA reaction to clone the full-length begomoviral se‐ AYVV_KM051527_MYS quences in this study, were used as templates for PCR 100 AYVV_JN809821_THA AYVV_AB100305_IDN to check the specificities of each primer pair (Fig. 3). AYVV_EF544600_CHA For PepYLCIV, DNA extracted from agroinoculated 99 100 AYVV_AB306314_JPN 100 pepper No. 218 in our previous study (Koeda et al., 100 100 AYVV_DQ866134_TWN 2018) was used as template for PCR. All the primer MAYVCV_KP188831_KHM ToLCLV_AF195782_LAO pairs successfully detected specific begomovirus

100 PepLCV_AF414287_MYS species (Fig. 3). Diagnosis of begomoviruses infecting PepYLCTHV_KT322141_THA tomato, tobacco, eggplant, cucumber, and squash was TYLCTHV_AF206674_MMR made by PCR using these begomovirus species-specific 100 EYVThV_KY373213_THA TYLCVN_GQ246941_VNM primers (Table 4). TYLCKaV was detected in tomato 100 TLCCNV_KU980918_CHN and tobacco; however, the infection rate was highest in 99 LuYVV_DQ641708_VNM eggplant indicating that it is the primary host for CLCuMuV_KY373212_THA TYLCKaV. This result is consistent with other studies 89 LYMV_AF509739_VNM PuYMV_EF197941_MYS that have reported that most TYLCKaV have been iso‐ ToLCNDV_KU992383_THA 94 87 lated from eggplants in Southeast Asia (Green et al., 94 ToLCNDV_GU180095_TWN 2003; Ha et al., 2008; Tang et al., 2014; Bagewadi and 100 ToLCNDV-[BACu-20]_IDN *** 90 ToLCNDV_JN809814_THA Naidu, 2016; Zhang et al., 2018). Interestingly, 100 78 ToLCNDV_KM383737_BGD PepYLCIV, PepYLCAV, and AYVV were undetected in 100 ToLCNDV_LC431619_IDN eggplant, indicating that the primary hosts for these ToLCNDV_KP178729_IND viruses are chili pepper, tomato, and tobacco. AYVV 81 SLCCNV_KF188433_IND 100 SLCCNV_JN587811_IND was detected in tomato and tobacco, and the infection SLCCNV_EU487031_PHL rate was particularly high in tobacco, indicating that to‐ 89 100 SLCCNV_AB330078_THA bacco was the primary host. In contrast, ToLCNDV and 99 SLCCNV_KC857509_VNM 99 SLCCNV_AM260206_CHN SLCCNV, which are genetically distinct species in the SLCCNV_EF197940_MYS phylogenetic analysis (Fig. 2), were only detected from SLCCNV-[BASq-17]_IDN *** SLCCNV_KY652743_TLS

100 SLCuV_KY652743_TLS Fig. 2. Phylogenetic tree based on alignments of complete begomo‐ 100 PepYLCIV_LC051114_IDN * 100 PepYLCIV_LC051115_IDN * virus DNA-A sequences. The placements of the begomovirus 97 PepYLCIV_LC051113_IDN * isolates in Koeda et al. (2016), Kesumwati et al. (2019), and

93 PepYLCIV_KT809346_IDN this study are indicated by *, **, and ***, respectively. The se‐ PepYLCIV_DQ083765_DN quences of DNA-A were obtained from GenBank and used for 100 76 PepYLCIV_AB267838_IDN phylogenetic analysis. The bootstrap values are indicated at the PepYLCIV_LC051112_IDN * 100 nodes (based on 1000 replicates). The branch lengths are pro‐ PepYLCIV2_KT809345_IDN portional to the number of nucleotide changes, as indicated by 97 PepYLCAV_LC465989_IDN ** the scale bar (0.05 substitutions per site). The abbreviations PepYLCAV_LC465990_IDN ** 93 used are as follows: AYVV, Ageratum yellow vein virus; PepYLCAV_LC465992_IDN ** CLCuMuV, Cotton leaf curl Multan virus; EYVThV, Emilia 100 PepYLCAV_LC465991_IDN ** yellow vein Thailand virus; LYMV, Luffa yellow mosaic virus; PepYLCAV_LC387328_IDN ** LuYVV, Ludwigia yellow vein virus; MAYVCV, Malvastrum PepYLCAV_LC387327_IDN ** yellow vein Cambodia virus; PHYVV, Pepper huasteco yellow ToLCSuV_FJ237614_IDN vein virus; PepLCV, Pepper leaf curl virus; PepYLCAV, Pepper 92 TYLCKaV_DQ641702_VNM 90 yellow leaf curl Aceh virus; PepYLCIV, Pepper yellow leaf TYLCKaV_DQ169054_VNM 100 TYLCKaV_KR073086_KHM curl Indonesia virus; PepYLCIV2, Pepper yellow leaf curl TYLCKaV_KR073093_KHM Indonesia virus 2; PepYLCTHV, Pepper yellow leaf curl TYLCKaV_MF327262_CHN Thailand virus; PuYMV, Pumpkin yellow mosaic virus; 100 100 TYLCKaV_KU569606_VNM SLCCNV, Squash leaf curl China virus; SLCuV, Squash leaf

100 TYLCKaV_KF218820_LAO curl virus; TLCCNV, Tomato leaf curl China virus; ToLCLV, TYLCKaV_AF511530_THA Tomato leaf curl Laos virus; ToLCNDV, Tomato leaf curl New 90 TYLCKaV_AF511529_THA Delhi virus; ToLCSuV, Tomato leaf curl Sulawesi virus;

99 TYLCKaV_KF446661_IDN TYLCKaV, Tomato yellow leaf curl Kanchanaburi virus; TYLCKaV_KF446671_IDN TYLCTHV, Tomato yellow leaf curl Thailand virus; TYLCVN, 99 TYLCKaV_LC051116_IDN * Tomato leaf curl Vietnam virus. The abbreviations used for the TYLCKaV-[BAEg-1]_IDN *** 89 country of origin of the sequences are as follows: BGD, TYLCKaV-[BATa-20]_IDN *** 92 Bangladesh; CHN, China; IDN, Indonesia; IND, India; JPN, TYLCKaV-[BAEg-49]_IDN *** Japan; KHM, Cambodia; LAO, Laos; MEX, Mexico; MMR, PHYVV_MG582068_MEX Myanmar; MYS, Malaysia; PHL, Philippines; THA, Thailand;

0.050 TLS, East Timor; TWN, Taiwan; and VNM, Vietnam. Hort. J. 89 (4): 410–416. 2020. 415 the Cucurbitaceae crops. These results strongly suggest gation for PepYLCAV is required to clarify how that the infecting begomovirus species has differentiat‐ PepYLCAV emerged. ed from Solanaceae to Cucurbitaceae crops. In this study, isolation of the infecting begomovirus PepYLCAV is reported to be a recombinant between species and PCR-based diagnosis were conducted, and PepYLCIV infecting Solanaceae, the Pumpkin yellow this clarified the predominant begomovirus species for mosaic virus, and ToLCNDV infecting Cucurbitaceae each horticultural crop. Begomovirus species-specific (Kesumawati et al., 2019). Tobacco generally has a primers developed in our study will be useful for fur‐ broad spectrum of infecting viruses, but mix-infection ther phytopathological studies in other places in South‐ of Solanaceae-infecting and Cucurbitaceae-infecting east Asian countries, including Indonesia, where the begomoviruses was not detected in this study. Mix- same begomovirus species are present. This study pro‐ infection of these begomoviruses may be detected from vides valuable information for developing weeds in these fields. Therefore, further evaluation begomovirus-resistant cultivars in the future. Screening using larger numbers of samples or host-range investi‐ and breeding of begomovirus-resistant cultivars for chili pepper, tomato, eggplant, cucumber, and melon are ongoing research projects in our team. DNA templates used for PCR Acknowledgements We would like to thank Amirddin, Fitrizal, and Dharma Setiawan for supporting our fieldwork in Aceh. Author contribution statement (ToLCNDV) (TYLCKaV) (PepYLCAV, AYVV) (PepYLCIV) (SLCCNV) E. Kesumawati and S. Koeda designed experiment,

DNA ladder Pepper Tobacco Tobacco Cucumber Squash Water (Negative control) interpreted the results, and wrote the manuscript. S. Okabe conducted isolation of virus sequences and PCR. PepYLCIV M. Khalil, G. Alfan, P. Bahagia, N. Pohan, and S. Zakaria conducted field surveys and sampling. All au‐ PepYLCAV thors read and approved the final manuscript. Literature Cited TYLCKaV Bagewadi, B. and R. A. Naidu. 2016. First report of Tomato yel‐ low leaf curl Kanchanaburi virus in eggplant and tomato in AYVV Cambodia. Plant Dis. 100: 233. Briddon, R. W., S. E. Bull, S. Mansoor, I. Amin and P. G. Markham. 2002. Universal primers for the PCR-mediated ToLCNDV amplification of DNA beta: a molecule associated with some Begomovirus species specific primers monopartite begomoviruses. Mol. Biotechnol. 20: 315–318. SLCCNV Bull, S. E., R. W. Briddon and P. G. Markham. 2003. Universal primers for the PCR-mediated amplification of DNA 1: a satellite-like molecule associated with begomovirus-DNA Fig. 3. Detection of PepYLCIV, PepYLCAV, TYLCKaV, AYVV, ToLCNDV, and SLCCNV by PCR using begomovirus species- beta complexes. Mol. Biotechnol. 23: 83–86. specific primers. DNA of original samples used for isolating Edgar, R. C. 2004. MUSCLE: multiple sequence alignment with the full-length sequences by RCA (Table 3) were used as tem‐ high accuracy 5 and high throughput. Nucleic Acids Res. 32: plates for PCR; Tobacco [BATa-9] DNA for PepYLCAV and 1792–1797. AYVV, Tobacco [BATa-20] DNA for TYLCKaV, Cucumber Green, S. K., W. S. Tsai, S. L. Shih, M. A. Rezaian and U. [BACu-20] DNA for ToLCNDV, Squash [BASq-17] DNA for Duangsong. 2003. Molecular characterization of a new SLCCNV. For detection of PepYLCIV, DNA extracted from begomovirus associated with tomato yellow leaf curl and No. 218 (C. annuum) agroinoculated with PepYLCIV in our eggplant yellow mosaic diseases in Thailand. Plant Dis. 87: previous study (Koeda et al., 2018) was used as a template for 446. Table 4.PCR. Diagnosis of infecting begomovirus species.

Table 4. Diagnosis of infecting begomovirus species.

Number of Infected begomovirus species Host samples collected PepYLCIV PepYLCAV TYLCKaV AYVV ToLCNDV SLCCNV Eggplant 16 0 (0%) 0 (0%) 11 (69%) 0 (0%) 0 (0%) 0 (0%) Tomato 43 3 (7%) 30 (70%) 1 (2%) 7 (16%) 0 (0%) 0 (0%) Tobacco 27 0 (0%) 7 (26%) 5 (19%) 20 (74%) 0 (0%) 0 (0%) Cucumber 4 0 (0%) 0 (0%) 0 (0%) 0 (0%) 3 (75%) 1 (25%) Squash 5 0 (0%) 0 (0%) 0 (0%) 0 (0%) 0 (0%) 4 (80%) 416 E. Kesumawati, S. Okabe, M. Khalil, G. Alfan, P. Bahagia, N. Pohan, S. Zakaria and S. Koeda

Ha, C., S. Coombs, P. Revill, R. Harding, M. Vu and J. Dale. Kumar, S., G. Stecher and K. Tamura. 2016. MEGA7: Molecular 2008. Molecular characterization of begomoviruses and evolutionary genetics analysis version 10 7.0 for bigger DNA satellites from Vietnam: additional evidence that the datasets. Mol. Biol. Evol. 33: 1870–1874. New World geminiviruses were present in the Old World Mabvakure, B., D. P. Martin, S. Kraberger, L. Cloete, S. van prior to continental separation. J. Gen. Virol. 89: 312–326. Brunschot, A. D. Geering, J. E. Thomas, K. Bananej, J. M. Hanley-Bowdoin, L., E. R. Bejarano, D. Robertson and S. Lett, P. Lefeuvre, A. Varsani and G. W. Harkins. 2016. On‐ Mansoor. 2013. Geminiviruses: masters at redirecting and going geographical spread of Tomato yellow leaf curl virus. reprogramming plant processes. Nat. Rev. Microbiol. 11: Virology 498: 257–264. 777–788. Muhire, B. M., A. Varsani and D. P. Martin. 2014. SDT: a virus International Committee on Taxonomy of Viruses (ICTV). 2018. classifcation tool based on pairwise sequence alignment and Virus Taxonomy: 2018b Release. . journal.pone.0108277. Kenyon, L., W. S. Tsai, S. L. Shih and L. M. Lee. 2014. Emer‐ Navas-Castillo, J., E. Fiallo-Olivé and S. Sánchez-Campos. 2011. gence and diversity of begomoviruses infecting solanaceous Emerging virus diseases transmitted by whiteflies. Annu. crops in East and Southeast Asia. Virus Res. 186: 104–113. Rev. Phytopathol. 49: 219–248. Kesumawati, E., S. Okabe, K. Homma, I. Fujiwara, S. Zakaria, S. Saunders, K., I. D. Bedford, R. W. Briddon, P. G. Markham, Kanzaki and S. Koeda. 2019. Pepper yellow leaf curl Aceh S. M. Wong and J. Stanley. 2000. A unique virus complex virus: A novel bipartite begomovirus isolated from chili pep‐ causes Ageratum yellow vein disease. Proc. Natl. Acad. Sci. per, tomato, and tobacco plants in Indonesia. Arch. Virol. USA 97: 6890–6895. 164: 2379–2383. Scholthof, K. B. G., S. Adkins, H. Czosnek, P. Palukaitis, E. Koeda, S., K. Homma, Y. Tanaka, E. Kesumawati, S. Zakaria and Jacquot, T. Hohn, B. Hohn, K. Saunders, T. Candresse, P. S. Kanzaki. 2017. Highly efficient agroinoculation method Ahlquist, C. Hemenway and G. D. Foster. 2011. Top 10 for tomato plants with Tomato yellow leaf curl plant viruses in molecular plant pathology. Mol. Plant Kanchanaburi virus. Hort. J. 86: 479–486. Pathol. 12: 938–954. Koeda, S., K. Homma, Y. Tanaka, D. Onizaki, E. Kesumawati, S. Tang, Y. F., Z. F. He, Z. G. Du and L. H. Lu. 2014. First report of Zakaria and S. Kanzaki. 2018. Inoculation of with Tomato yellow leaf curl Kanchanaburi virus infecting egg‐ Pepper yellow leaf curl Indonesia virus by combining plant in Laos. Plant Dis. 98: 428. agroinoculation and grafting. Hort. J. 87: 364–371. Wilisiani, F., A. Tomiyama, H. Katoh, S. Hartono, Y. Neriya, H. Koeda, S., E. Kesumawati, Y. Tanaka, M. Hosokawa, M. Doi and Nishigawa and T. Natsuaki. 2019. Development of a LAMP A. Kitajima. 2016. Mixed infection of begomoviruses on assay with a portable device for real-time detection of pepper plants at northern Sumatra, Indonesia. Trop. Agric. begomoviruses under field conditions. J. Virol. Methods Dev. 60: 59–64. 265: 71–76. Koeda, S., R. Takisawa, T. Nabeshima, Y. Tanaka and A. Zhang, X. Y., L. L. Zhao, J. Zhong, T. T. Li, Z. K. Zhang and Kitajima. 2015. Production of Tomato yellow leaf curl virus M. Ding. 2018. First report of Tomato yellow leaf curl free parthenocarpic tomato plants by leaf primordia-free Kanchanaburi virus infecting eggplant in China. Plant Dis. shoot apical meristem culture combined with in vitro graft‐ 102: 257. ing. Hort. J. 84: 327–333.