Genetic recombination among tomato yellow leaf curl virus isolates in commercial tomato crops in Kuwait drives emergence of virus diversity: a comparative genomic analysis
BMC Research Notes volume 16, Article number: 71 (2023)
Whitefly-transmitted tomato yellow leaf curl virus (TYLCV) continues to be a major constraint to tomato production in Kuwait. However, very limited information is available about the population structure and genetic diversity of TYLCV infecting tomato in Kuwait.
Whole genome sequences of 31 isolates of TYLCV, collected from commercial tomato crops grown in northern (Abdally) and southern (Al Wafra) parts of Kuwait, were deciphered. Eighteen isolates of TYLCV are identified as potential genetic recombinants. The isolates Abdally 6A and Abdally 3B reported in this study were identified to be potential recombinants. Compared to the 15 isolates from the Abdally area, and the three previously reported KISR isolates of Kuwait, six out of sixteen Al Wafra isolates showed an insertion of 19 extra nucleotides near the 5′-end. There are also four nucleotide variations before the 19-extra-nucleotides. The additional 19 nucleotides observed in nine isolates indicate that these isolates might have resulted from a single gene recombination/insertion event. Molecular phylogeny based on complete genome sequences of TYLCV isolates suggests transboundary movement of virus isolates due to geographic proximity. The information presented herein is quite useful for the comprehension of TYLCV biology, epidemiology and would aid in the management of disease in the long run.
Tomato yellow leaf curl virus (TYLCV) (Genus Begomovirus, Family Geminiviridae) , causes significant damages to tomato in both quantitative and qualitative terms. Symptoms associated with TYLCV infection in tomato include chlorotic leaf edges, upward leaf cupping, leaf mottling, reduced leaf size, and flower drop. Symptoms on the fruit make the agricultural product unmarketable. Thus, TYLCV can inflict a severe impact on tomato production. Tomato plants infected at an early stage won't bear fruit and plants are severely stunted. In early virus infections, 100% crop loss is not uncommon. The virus is generally transmitted by whitefly, Bemisia tabaci, and the vector has developed resistance to many of the insecticides, which limits the vector control option [2,3,4,5,6].
The virion of TYLCV is composed of twin icosahedral capsids. The viral genome consists of a single-stranded circular DNA genome of about 2.7 kb, which encodes six open reading frames (ORFs), two (V1 and V2) in viral sense and four (C1- C4) in complementary sense. All the proteins encoded by begomoviruses are multifunctional. Protein encoded by ORF V1 is involved in encapsidation of the genome into geminate particles, insect transmission and long distant movement. Protein derived from V2 ORF plays a role in cell-to-cell movement and suppression of host-mediated RNA silencing. C1 protein is involved in viral DNA replication, while ORF C2-derived protein activates viral DNA transcription and is also a host RNA silencing suppressor. Protein encoded by ORF C3 plays a role in enhancing viral DNA replication, and ORF C4 encoded protein is known to be a RNA silencing suppressor [7, 8]. Besides these ORFs, the TYLCV genome has a 300-nucleotide (nt) non-coding intergenic region (IR), which possesses a bidirectional promoter and a core conserved nona-nucleotide sequence (TAATATTAC). The nona-nucleotide sequence is recognized by the C1 protein, and it serves as the origin of replication. A great deal of information was generated over the years about the diversity, epidemiology and molecular biology of the disease and the virus [9,10,11].
Management of viral diseases, especially those caused by TYLCV, is challenging and expensive. The whitefly vector has a wide host range, and several weeds serve as reservoir green bridges for TYLCV to survive between crop seasons. Plant health management approaches include the use of resistant cultivars working on the principle of resistance (R) genes, RNA interference (RNAi), and vector management by insecticides, and cultural practices such as adjusting planting date and crop-free period [1, 12,13,14]. In tomato, host plant resistance to TYLCV was mainly based on the use of Ty genes . Nevertheless, the genomic features of TYLCV are altered owing to multiple factors such as virus recombination, genetic mutations, and inclusion of satellite components and invasion of exogenous whitefly species among others [16, 17]. Hence, TYLCV have the ability to overcome the endogenous and genetically engineered disease control measures, and there is a continuous need to develop novel virus management strategies that are effective and durable.
Viral diseases of plants are widespread in Kuwait and are causing significant economic losses reaching up to 95% in many vegetable crops . These crops comprise the majority of greenhouse and open field agriculture in Kuwait. Mosaic, abnormal leaf color, abnormal vein patterns of leaves, mottling in leaves, spotting patterns in leaves, and abnormal leaf shape, leaf curling, and yellowing were observed on plants grown in Al Wafra and Abdally during 2007, 2010, and 2014 seasons . Recent studies showed that TYLCV infection is causing major economic losses in tomato plants, even upto90%. TYLCV has been reported as a major tomato-infecting virus, but it has not been fully characterized at the molecular level and little is known about the genetic diversity of TYLCV populations in Kuwait. Whiteflies are the main vector for TYLCV [20, 21] and a considerable quantum of research was done regarding the management of whitefly vectors with a goal to reduce the impact of begomoviruses [20, 22, 23]. Considering the economic importance of the viral disease in tomato, caused due to TYLCV, and frequent disease outbreaks inflicting severe crop losses [18,19,20,21] this study was framed with an objective of characterizing the genome sequences of TYLCV infecting tomato in Kuwait and to perform a comprehensive genomic analysis. We report the molecular characterization of several TYLCV isolates collected from northern and southern Kuwait and the extent of genetic diversity among the isolates characterized herein and reported earlier from Kuwait, the role of virus recombination events, and the inferences from molecular phylogeny with a view to decipher the evolutionary genomics of the virus.
Materials and methods
Two commercial tomato farms (Abdally in northern Kuwait, and Al Wafra region in southern Kuwait) were surveyed and leaf samples from plants that displayed symptoms suggestive of TYLCV infection (640 sample) were collected and brought back to Kuwait Institute for Scientific Research, Environmental and Life Science Research Center, (KISR), Safat, 13109, Kuwait and stored at 4ºC until further use. The virus infected samples were collected during 2020 and 2021.
Total DNA extraction and TYLCV detection
Total genomic DNA was extracted using the CTAB method and tested for the presence of TYLCV by PCR using primer pair PTYc787 (GTTCGATAATGAGCCCAG) and PTYc1121(ATGTAACAGAAACTCATG) . The samples positive for TYLCV were used to sequence the full-length genome which was carried out in Department of Plant Pathology, Washington State University, Pullman, USA following rolling circle amplification (RCA) method. Briefly, the RCA reaction was performed by mixing the following components: genomic DNA (50 ng), 1 μL of exo-resistant random primers (500 μM, 35 optical density units (OU)/mL), 2 μL of 10X reaction buffer, to bring the volume to 10.2 μL using water. The RCA reaction mixture was heated at 95 °C for 5 min in water bath, followed by chilling on ice for 2 min. Following this, 2 μL of 10 mM dNTP, 1.6 μL of phi-29 DNA polymerase, 0.2 μL of pyrophosphate, inorganic (0.1U/μL) were added and incubated at 30 °C for 18 h. The RCA reaction was carried out in a thermal cycler machine. The reaction was stopped by heating the components at 65 °C for 10 min. The RCA derived reaction products were resolved in agarose gel electrophoresis and the amplicons were directly sequenced using the Sanger DNA sequencing protocol.
In silico analysis
Sequence comparisons of the fifteen Abdally TYLCV isolates, sixteen Al Wafra isolates and three previously reported KISR isolates were performed using the Sequence Demarcation Tool (SDT) . Sequence alignment was generated using ClustalW algorithm. A molecular phylogenetic tree was constructed using the maximum-likelihood method (default parameters with 1000 replicates in the bootstrap analysis) using these sequences along with TYLCV genome sequences available in GenBank utilizing MEGAX . Phylogeny was inferred from using, General Time Reversal (GTR) model with G + I (invariant sites and distributed range) . Potential recombination events were detected by Recombination Detection Program-4 (RDP 4 Beta 4.16) . All the complete TYLCV genome sequences reported from Kuwait and available in GenBank were used in the analysis and the sequence alignment was carried out using ClustalW in MEGA X  and the aligned sequences were used for recombination detection studies. For identifying recombination events, step-down correction with the highest acceptable p-value setting of 0.05 was used along with other default settings for all of nine methods (RDP, Chimaera, BootScan, 3Seq, GENECONV, MaxChi, SiScan and LARD, PhylPro) available in the RDP 4 .
Results and discussion
Molecular genomic analysis of geminiviruses has been the mainstay in deciphering the genomic diversity of the virus species infecting economically important crops and in devising suitable disease control measures [1, 2, 14]. Hence, in this study, the complete genomic sequences of thirty-one TYLCV isolates from Kuwait were sequenced and deposited in GenBank database (Table 1). The nucleotide identity analysis of various TYLCV genes from these isolates were shown along with the previously reported three isolates (KISR, KISR2 and KISR3) from Kuwait [18, 19] (Fig. 1). The nucleotide sequence identity of Kuwaiti isolates of TYLCV range from 91.20 to 100%. Compared to all fifteen Abdally TYLCV isolates, the three previously reported KISR isolates and six out of sixteen Al Wafra isolates have 19 extra nucleotides (TTCTTTCTAGGTGTGCCCC) in the intergenic region (Fig. 2). There exist 4-nucleotide variations before the 19 extra nucleotides. The three previously reported isolates and six Al Wafra isolates had “G/CCTT” before the 19 extra nucleotides, whereas the others possessed “AAA(A)”.
Phylogenetic analyses were conducted to compare complete nucleotide sequences with those from other parts of the world and to infer molecular evolution of current isolates (Fig. 3). A total of eighty-seven TYLCV isolates were analyzed, which formed two distinct clades. All 34 Kuwaiti isolates (31 from this study and three previously reported) grouped into a clade along with the isolates from China, Iran, Israel, Japan, Jordan, Mexico, Oman, Portugal, Turkey and the US (Fig. 3). TYLCV isolates reported from Cuba, Japan, Lebanon, Morocco, The Netherlands, and Spain formed a distinct clade (Fig. 3). Earlier study on molecular analysis of TYLCV in the Arabian Peninsula and Iran suggests these regions as the centre of diversity .
To identify putative recombination events, an analysis was performed using the RDP4 program. Fourteen recombination events in the 87 TYLCV full-length sequences from worldwide were detected by at least four methods in the RDP4 program (Table 2). The isolates Abdally 6A (OL890669) and Abdally 3B (OL890670) reported in this study were identified to be potential recombinants. In addition, TYLCV isolates Abdally 11B (OL890676), Abdally 13B (OL890678), and Abdally 3B (OL890670) served as major parents for generation of different recombinants. It suggests the occurrence of microevolution within the TYLCV populations in Kuwait (Table 2 and Fig. 4). Further, genetic recombinant TYLCV isolates were shown to have ecologically selective advantage over their parental viruses .
Comparative sequence analysis of Kuwaiti TYLCV isolates reveal > 90% sequence identity, suggesting that these isolates may have derived from a few parental strains originally imported into the country (Fig. 1). The additional 19 nucleotides observed in nine Kuwaiti isolates indicate that these isolates might have resulted from a single gene recombination/insertion event. The further integration of “CCTT” motif could represent a second mutation event occurred at a later time (Fig. 2). The phylogenetic analysis of all known TYLCV sequences reported worldwide suggests that TYLCV seems to undergo a relatively rapid evolution and exhibit significant genetic variation. The isolates Al Wafra 1, Al Wafra 7, Al Wafra2, AlWafra 16, Al Wafra 3, and Al Wafra 6 were genetically distinct and formed a separated cluster. However, some TYLCV isolates viz., Al Wafra 17, Al Wafra 23, Al Wafra 24 and Abdally 15 shared phylogenetic lineage with TYLCV isolates reported from Iran and Oman (Fig. 3). On the other hand, all the Abdally isolates, excluding Abdally 15, showed monophyletic origin along with TYLCV isolates reported from Oman (Fig. 3).The phylogenetic relatedness of Kuwait TYLCV isolates with that of Iran and Oman imply the cross border transfer or movement of virus isolates among the countries in Persian Gulf region (Fig. 3). If the observed diversity in sequences in selected isolates has any role in modulating the TYLCV-tomato interactions remains to be seen and the current study lays groundwork for further investigations.
In conclusion, we report here the complete genome sequences of 31 isolates of TYLCV circulating in the tomato fields of Northern and Southern Kuwait. A comprehensive comparative genomic analysis of the novel TYLCV isolates along with those reported earlier identified recombination events that could possibly involve in the evolution of TYLCV causing the generation of novel variants of concern. The information presented in this article will be quite useful for the comprehension of TYLCV biology, epidemiology and disease control measures.
No limitations were encountered during the study. Our findings form the basis for further studies on the functional significance of the extra 19 nucleotides found in the 5′ end of genomic sequences of KISR and six Al Wafra isolates.
Availability of data and materials
The viral genome sequence datasets generated and described in this research are deposited in the NCBI GenBank database (ncbi.nlm.nih.gov) and are available under accession numbers (Al-Wafra isolates: OM691678 to OM691692) and (Al-Abdally isolates: OL890666 to OL890680).The data described can be freely and openly accessed from the repository.
Rojas MR, Macedo MA, Maliano MR, Soto-Aguilar M, Souza JO, Briddon RW, Kenyon L, Rivera Bustamante RF, Zerbini FM, Adkins S, Legg JP. World management of geminiviruses. Annu Rev Phytopathol. 2018;56:637–77. https://doi.org/10.1146/annurev-phyto-080615-100327.
Moriones E, Navas-Castillo J. Tomato yellow leaf curl virus, an emerging virus complex causing epidemics worldwide. Virus Res. 2000;71:123–34. https://doi.org/10.1016/s0168-1702(00)00193-3.
Seal SE, Jeger MJ, Van den Bosch F. Begomovirus evolution and disease management. Adv Virus Res. 2006;67:297–316. https://doi.org/10.1016/S0065-3527(06)67008-5.
Ghanim M. A review of the mechanisms and components that determine the transmission efficiency of Tomato yellow leaf curl virus (Geminiviridae; Begomovirus) by its whitefly vector. Virus Res. 2014;186:47–54. https://doi.org/10.1016/j.virusres.2014.01.022.
Rosen R, Kanakala S, Kliot A, Pakkianathan BC, Farich BA, Santana-Magal N, Elimelech M, Kontsedalov S, Lebedev G, Cilia M, Ghanim M. (2015) Persistent, circulative transmission of begomoviruses by whitefly vectors. Curr Opin Virology. 2015;15:1–8. https://doi.org/10.1016/j.coviro.2015.06.008.
Diaz-Pendon JA, Cañizares MC, Moriones E, Bejarano ER, Czosnek H, Navas-Castillo J. Tomato yellow leaf curl viruses: menage a trois between the virus complex, the plant and the whitefly vector. Mol Plant Pathol. 2010;11:441–50. https://doi.org/10.1111/j.1364-3703.2010.00618.x.
Hanley-Bowdoin L, Settlage SB, Robertson D. Reprogramming plant gene expression: a prerequisite to geminivirus DNA replication. Mol Plant Pathol. 2004;5:149–56. https://doi.org/10.1111/j.1364-3703.2004.00214.x.
Fondong VN. Geminivirus protein structure and function. Mol Plant Pathol. 2013;14:635–49. https://doi.org/10.1111/mpp.12032.
Yang X, Wang B, Luan J, Xie Y, Liu S, Zhou X. Molecular variation of tomato yellow leaf curl virus in the insect vector Bemisia tabaci. Sci Reports. 2017;7:16427. https://doi.org/10.1038/s41598-017-16330-4.
Chen LF, Rojas M, Kon T, Gamby K, Xoconostle-Cazares B, Gilbertson RL. A severe symptom phenotype in tomato in Mali is caused by a reassortant between a novel recombinant begomovirus (Tomato yellow leaf curl Mali virus) and a betasatellite. Mol Plant Pathol. 2009;10:415–30. https://doi.org/10.1111/j.1364-3703.2009.00541.x.
Hosseinzadeh MR. Phylogenetic relationships, recombination analysis, and genetic variability among diverse variants of tomato yellow leaf curl virus in Iran and the Arabian Peninsula: further support for a TYLCV center of diversity. Arch Virol. 2014;159:485–97. https://doi.org/10.1007/s00705-013-1851-z.
Lapidot M, Legg JP, Wintermantel WM, Polston JE. Management of whitefly-transmitted viruses in open-field production systems. Adv Virus Res. 2014;90:147–206. https://doi.org/10.1016/B978-0-12-801246-8.00003-2.
Vu TV, Choudhury NR, Mukherjee SK. Transgenic tomato plants expressing artificial microRNAs for silencing the pre-coat and coat proteins of a begomovirus, Tomato leaf curl New Delhi virus, show tolerance to virus infection. Virus Res. 2013;172:35–45. https://doi.org/10.1016/j.virusres.2012.12.008.
Ramesh SV, Sahu PP, Prasad M, Praveen S, Pappu HR. Geminiviruses and plant hosts: a closer examination of the molecular arms race. Viruses. 2017;9:256. https://doi.org/10.3390/v9090256.
Gill U, Scott JW, Shekasteband R, Ogundiwin E, Schuit C, Francis DM, Sim SC, Smith H, Hutton SF. Ty-6, a major begomovirus resistance gene on chromosome 10, is effective against Tomato yellow leaf curl virus and Tomato mottle virus. TheorAppl Genet. 2019;132:1543–54. https://doi.org/10.1007/s00122-019-03298-0.
Czosnek H, Ghanim M. Bemisia tabaci–Tomato yellow leaf curl virus interaction causing worldwide epidemics. In: Thompson WMO, editor. The Whitefly, Bemisia tabaci (Homoptera: Aleyrodidae) interaction with geminivirus-infected host plants. Dordrecht: Springer; 2011. p. 51–67. https://doi.org/10.1007/978-94-007-1524-0_3.
Belabess Z, Peterschmitt M, Granier M, Tahiri A, Blenzar A, Urbino C. The non-canonical tomato yellow leaf curl virus recombinant that displaced its parental viruses in southern Morocco exhibits a high selective advantage in experimental conditions. J Gen Virol. 2016;97:3433–45. https://doi.org/10.1099/jgv.0.000633.
Al-Ali EM, Al-Hashash H, Al-Shayji N, Ben Hejji A, Al-Aqeel H.2014 Molecular characterization and identification of whitefly-transmitted TYLCV disease in tomato plant in Kuwait. Kuwait Institute for Scientific Research, Report No. KISR12261, Kuwait.
Al-Ali EM, Al-Hashash H, Al-Aqeel H, Ben HejjiA Application of detection technologies for viral diseases of important plant crops in Kuwait. Kuwait Institute for Scientific Research, Report No. KISR10251, 2010; Kuwait.
Czosnek H. Tomato yellow leaf curl virus disease: management, molecular biology, breeding for resistance. Dordrecht: Springer; 2007. https://doi.org/10.1007/978-1-4020-4769-5.
JiangYX, de Blas C, BarriosL, Fereres A,. Correlation between whitefly (homoptera: aleyrodidae) feeding behavior and transmission of tomato yellow leaf curl virus. Ann Entomol Soc. 2000;2000(93):573–9. https://doi.org/10.1603/00138746(2000)093[0573:CBWHAF]2.0.CO;2.
Czosnek H, Hariton-Shalev A, Sobol I, Gorovits R, Ghanim M. The incredible journey of begomoviruses in their whitefly vector. Viruses. 2017;9:273. https://doi.org/10.3390/v9100273.
Salati R, Nahkla MK, Rojas MR, Guzman P, Jaquez J, Maxwell DP, Gilbertson RL. Tomato yellow leaf curl virus in the Dominican Republic: characterization of an infectious clone, virus monitoring in whiteflies, and identification of reservoir hosts. Phytopathology. 2002;92:48C7-96. https://doi.org/10.1094/PHYTO.2002.92.5.487.
Muhire BM, Varsani A, Martin DP. a virus classification tool based on pairwise sequence alignment and identity calculation. PLoS ONE. 2014;9:e108277. https://doi.org/10.1371/journal.pone.0108277.
Kumar S, Stecher G, Li M, Knyaz C, Tamura KMEGAX. molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol. 2018;35(6):1547.
Martin DP, Lemey P, Lott M, Moulton V, Posada D, Lefeuvre P. RDP3: aflexible and fast computer program for analyzing recombination. Bioinformatics. 2010;26:2462–3.
This work was supported by Kuwait Foundation for the Advancement of Sciences under project code PR18-12SL-09.
The authors did not receive support from any organization for the submitted work.
Ethics approval and consent to participants
No human participants, human material, or human data were involved in this study.
Consent for publication
Authors declare no financial or non-financial interests that are directly or indirectly related to this work.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Al-Ali, E., Al-Hashash, H., Akbar, A. et al. Genetic recombination among tomato yellow leaf curl virus isolates in commercial tomato crops in Kuwait drives emergence of virus diversity: a comparative genomic analysis. BMC Res Notes 16, 71 (2023). https://doi.org/10.1186/s13104-023-06319-w