Molecular variability and evolution of a natural population of tomato yellow leaf curl virus in Shanghai, China

Molecular variability and evolution of a natural population of tomato yellow leaf curl virus in Shanghai, China 番茄黄化曲叶病毒自然种群的分子变异和进化研究 first-author corresp ‡Corresponding Author ‡Corresponding Author †E-mail:xiey@zju.edu.cn †E-mail:xiey@zju.edu.cn xiey@zju.edu.cn State Key Laboratory of Rice Biology, Institute of Biotechnology, Zhejiang University, Hangzhou 310058, China State Key Laboratory of Rice Biology, Institute of Biotechnology, Zhejiang University, Hangzhou 310058, China Tomato yellow leaf curl virus (TYLCV), belonging to the genus Begomovirus of the family Geminiviridae, is emerging as the most destructive pathogen of tomato plants. Since the first report of TYLCV in Shanghai, China in 2006, TYLCV has spread rapidly to 13 provinces or autonomous regions of China. In this study, the molecular variability and evolution of TYLCV were monitored in Shanghai from its first upsurge in 2006 until 2010. Full-length genomic sequences of 26 isolates were obtained by rolling circle amplification. Sequence analysis showed that the intergenic region was the most variable, with a mean mutation rate of 4.81×10 −3 nucleotide substitutions per site per year. Genetic differentiation was found within isolates obtained from 2006, 2009, and 2010, though a linear increase in genetic diversity over time was not evident. Whilst significant parts of TYLCV genes were under negative selection, the C4 gene embedded entirely within the C1 gene had a tendency to undergo positive selection. Our results indicate that a mechanism of independent evolution of overlapping regions could apply to the natural population of TYLCV in Shanghai, China. 研究目的 2 2006年我国上海首次发现番茄黄化曲叶病毒（TYLCV），随后TYLCV迅速蔓延至全国13个省份和自治区。本研究分析了2006至2010年期间TYLCV在我国首发地上海市的分子变异规律。 创新要点 2 本研究持续五年追踪田间TYLCV，分析TYLCV的全长基因组序列、分子变异及种群遗传结构，为防控TYLCV提供理论依据。 研究方法 2 2006至2010年从上海采集了26个TYLCV的分离物，利用高保真性的滚环扩增技术获得TYLCV分离物的全长基因组。应用MEGA5等生物信息学软件分析TYLCV的分子变异。 重要结论 2 TYLCV自然种群具有与RNA病毒相似的突变率，以基因间隔区的分子变异最大，平均突变率为4.81×10 −3 （见图2和表2）。TYLCV的大部分基因都处于负向选择，但包含在C1开放阅读框内的C4，却与C1承受着不同的选择压而处于正向选择（见图3和表6）。 Molecular variability Tomato yellow leaf curl virus (TYLCV) Geminivirus Evolution 分子变异 番茄黄化曲叶病毒（TYLCV） 双生病毒 进化 Introduction 1 Geminiviruses are a group of circular, single-stranded DNA viruses that cause devastating diseases worldwide. Four genera ( Mastrevirus , Begomovirus , Curtovirus , and Topocuvirus ) have been taxonomically classified within the family Geminiviridae, with the genus Begomovirus being the largest, containing more than 200 species. Earlier studies have demonstrated that mutation, recombination, and pseudo-recombination are the major contributory factors driving the emergence of begomovirus variants (Czosnek and Ghanim, 2011 2011 ). Although geminiviruses depend on host DNA polymerase for replication, mutation occurs frequently in geminivirus genomes and progeny populations of high variability can be generated in a short period of time. Mean rates of 1.60×10 −3 and 1.33×10 −4 nucleotide substitutions per site per year (subs/site/year) have been estimated for the DNA-A and DNA-B components, respectively, of the East African cassava mosaic virus (EACMV) (Duffy and Holmes, 2009 2009 ). Similar rates of substitution have also been reported for artificially inoculated hosts (Ge et al., 2007 2007 ). It is now clear that geminiviruses evolve as rapidly as many RNA viruses, but only a few reports estimate the time required to generate such levels of genetic variability. Tomato yellow leaf curl virus (TYLCV) is a monopartite begomovirus belonging to the genus Begomovirus , and is transmitted by the whitefly Bemisia tabaci in a circulative and persistent manner. The genome of TYLCV encodes two partially overlapping genes in the viral strand: V1 encodes a coat protein, and V2 encodes a precoat protein and acts as a suppressor of RNA silencing (Glick et al., 2008 2008 ; Jiang et al., 2012 2012 ). The complementary strand of TYLCV encodes four proteins: a replication-associated protein (C1/Rep), a transcriptional activator protein (C2/TrAP), a replication enhancer protein (C3/REn), and protein C4, which is entirely embedded in the C1 gene, but has a different reading frame. Open reading frames (ORFs) are organized bi-directionally and are separated by an intergenic region (IR) which serves as the origin of viral replication and as a bi-directional promoter (Hanley-Bowdoin et al., 1999 1999 ). Tomato yellow leaf curl disease (TYLCD) caused by TYLCV is becoming the most destructive disease in tomato plants (Moriones and Navas-Castillo, 2000 2000 ). Plants infected by TYLCV during the early growth stages frequently suffer yield losses of up to 100%. Although tomato plants in Middle East countries have been severely affected by TYLCV only since the 1960s, the first report of damage caused by this disease dates back to 1931 in Israel (Cohen and Antignus, 1994 1994 ). Since then, epidemics have emerged and devastated tomato crops rapidly and widely. TYLCV was first reported in Shanghai, China in 2006 (Wu et al., 2006 2006 ), and it has since spread rapidly to 13 provinces or autonomous regions of China, including Guangdong, Zhejiang, Jiangsu, Anhui, Shandong, Hebei, Henan, Beijing, Tianjin, Shanxi, Shaanxi, and Inner Mongolia and Xinjiang Uygur Autonomous Regions (Zhang et al., 2009 2009 ). Since no other begomoviruses had been detected in Shanghai before 2006, the rapid spread of TYLCV provided a great opportunity to understand the evolution of a plant DNA virus soon after its first introduction to a new area. Twenty-six TYLCV isolates were collected from the same area of Shanghai during 2006–2010. Full-length genomic sequencing allows for a more comprehensive analysis of virus variability and of the different forces driving the evolution of TYLCV. Materials and methods 1 Collection of TYLCV-infected tomato plants 2 Twenty-six tomato plants (cultivar Pufen 1, formerly called 98-8) with typical TYLCV symptoms (leaf curling, yellowing, and stunted growth) were collected from Shanghai from 2006 to 2010. TYLCV isolates were named after the year in which they were collected (Table 1 1 ). Isolate Collection date SH2006-1 2006 SH2006-2 2006 SH2007-1 2007 SH2007-2 2007 SH2007-3 2007 SH2007-4 2007 SH2007-5 2007 SH2007-6 2007 SH2008-1 2008 SH2008-2 2008 SH2008-3 2008 SH2008-4 2008 SH2008-5 2008 SH2008-6 2008 SH2009-1 2009 SH2009-2 2009 SH2009-3 2009 SH2009-4 2009 SH2009-5 2009 SH2009-6 2009 SH2010-1 2010 SH2010-2 2010 SH2010-3 2010 SH2010-4 2010 SH2010-5 2010 SH2010-6 2010 Genomic DNA extraction and detection of TYLCV 2 Total genomic DNA was extracted from 0.5 g of leaf tissue using the cetyltrimethylammonium bromide (CTAB) method (Zhou et al., 2001 2001 ). After washing in 70% ethanol, the DNA pellets were dried at room temperature and re-suspended in about 200 μl of distilled water. Primers TYLCV/F (5'-GAGCTCTTAGCTGCCTGAATGTTC-3') and TYLCV/R (5'-GAGCTCAACAGATGTCAAGACCTAC-3'), which were expected to specifically amplify a 2.8-kb fragment of TYLCV, were used to detect TYLCV infection. Isolation, cloning, and sequencing of full-length genomes 2 Rolling circle amplification (RCA) based on the phi29 DNA polymerase available in the TempliPhi™ kit (GE Healthcare, UK) was used to amplify circular viral genomes from 1 μg of total plant DNA extractions (Guo et al., 2009 2009 ). The amplified DNA concatemers (5 μl) were cleaved with Sac I to generate full-length viral fragments, which were subsequently inserted into linearized pGEM-7Zf+ (Promega, Madison, WI). About five independent clones from each sample were randomly chosen for sequencing by an automated model 3730 DNA sequencer (Perkin-Elmer, USA) using M13 forward and reverse primers and a pair of walking primers, TYLCV-W/F (5'-TCTGCAATCCAGGACCTACC-3') and TYLCV-W/R (5'-AGTCTATCTTGCAATATGTG-3'). Sequences were assembled and edited using EditSeq (Lasergene Package of Programs from DNASTAR Inc.). The fidelity of phi29 DNA polymerase was pre-determined by sequencing 20 individual clones generated directly from the TYLCV infectious clone using this enzyme during the initial amplification trials, and the mutation frequency was estimated to be 1.8×10 −5 . Sequence analysis 2 The circular genomic sequences were arranged to begin at the nick-site in the invariant nonanucleotide sequence at the origin of replication (TAATATT-3'/5'-AC) (Padidam et al., 1995 1995 ). Multiple sequence alignment was generated with the aid of ClustalX 1.83 (Thompson et al., 2002 2002 ). The level of genetic variation was calculated by estimating Watterson’s estimator of θ and the average number of nucleotide differences per site between two sequences (Pi) using DnaSP v5 (Guo et al., 2008 2008 ; Librado and Rozas, 2009 2009 ). The ratio of pairwise nonsynonymous ( d n ) to synonymous ( d s ) nucleotide substitutions per site, was estimated using the Pamilo-Bianchi-Li method implemented in MEGA 5 (Tamura et al., 2011 2011 ). Confidence estimates for nonsynonymous and synonymous nucleotide substitutions were assessed by the bootstrap method based on 100 replicates. The extent of mutational bias and the number of unambiguous nucleotide changes of each type were calculated by MEGA 5 (Tamura et al., 2011 2011 ). Genetic distances within and between subpopulations were calculated by applying Kimura’s two-parameter method in MEGA 5 (Tamura et al., 2011 2011 ). Evolutionary constraints acting on the overlapping regions were assessed using entropy values and BioEdit software (Hall, 1999 1999 ). Cumulative synonymous and nonsynonymous mutations were determined using Synonymous Non-synonymous Analysis Program (SNAP) ( http://www.hiv.lanl.gov/content/sequence/SNAP/SNAP.html http://www.hiv.lanl.gov/content/sequence/SNAP/SNAP.html ) as described by Korber ( 2002 2002 ) and Melgarejo et al. ( 2013 2013 ). Results 1 Symptom observation and detection of TYLCV 2 From 2006 to 2010, a total of 26 samples were obtained from surveys conducted once a year in Shanghai, China. All the plants sampled showed typical TYLCV symptoms, such as leaf curling, yellowing, and stunted growth (Fig. 1a 1a ). The specific 2.8-kb DNA fragment was amplified by polymerase chain reaction (PCR) from all the samples using primers TYLCV/F and TYLCV/R (Fig. 1b 1b ), indicating that these symptomatic tomato plants were infected by TYLCV. To obtain the complete nucleotide sequence of each of the 26 isolates, RCA followed by Sac I digestion of the DNA concatemers was performed to obtain the full-length viral fragment (Fig. 1d 1d ). About five independent clones from each sample were chosen to be sequenced and were found to be 2781 nucleotides (nts) long with a genomic organization typical of begomoviruses (Fig. 1c 1c ). (a) TYLCV-infected tomato plants showing leaf curling, yellowing, and stunted growth symptoms. (b) PCR detection of TYLCV in six of the infected tomato plants (target size, 2.8 kb). (c) Genomic organization of TYLCV. The intergenic region and six ORFs encoded by the viral strand and complementary strand are shown as indicated. (d) Full-length TYLCV fragment obtained from Sac I digestion of rolling circle amplification products. IR: intergenic region; M: DNA marker Molecular variation of the TYLCV genome 2 To evaluate the extent of TYLCV variation, the mean pairwise diversity index π and the Watterson parameter θ were calculated at all sites along the TYLCV genome by DnaSP Version 5.0 (Fig. 2 2 ). Using a sliding window of 100 nts with a 25-step size, we found that while the overall variability for the TYLCV genome was below 3% at the nucleotide level, distribution patterns were uneven. The peak of variability was located within the intergenic region and 5'-terminal part of V2, whereas the most conserved regions were located in the overlapping regions of V1 and V2. Estimation of the average mutation rate showed that the mutation rate of the full-length TYLCV genome reached 1.69×10 −3 subs/site/year (Table 2 2 ), which is within the ranges for RNA viruses and EACMV (Jenkins et al., 2002 2002 ; Duffy and Holmes, 2009 2009 ). Further analysis showed that the substitution rate was even higher for the intergenic region, with a substitution rate of 4.81×10 −3 subs/site/year. Rates of nucleotide substitution were also higher in the genes ( C1 , C2 , C3 , and C4 ) encoded on the complementary strand than those ( V1 and V2 ) encoded on the virion strand (Table 2 2 ). Statistical tests revealed that the substitution rates estimated for the natural population of TYLCV were significantly higher than those generated from RCA error (1.8×10 −5 ). These results indicate that the intergenic region had evolved significantly faster than the other ORFs. The sliding window size is 100 sites wide, slide by 25-site intervals. The relative positions of the six ORFs and the intergenic region of the TYLCV genome are illustrated by lines above the plot. IR: intergenic region Genetic distance: the numbers of base substitutions per site from averaging over all sequence pairs; data are shown as mean±standard error ( n =100), obtained by a bootstrap procedure The analyses were conducted using the Kimura’s two-parameter model 2 Genomic region 2 Sequence length (nt) 2 Number of sequences 2 Genetic distance * * 2 Nucleotide substitution rate (substitutions per site per year) Mean value Corrected mean value Full genome 2781 122 0.007±0.001 1.69×10 −3 1.67×10 −3 IR 313 122 0.023±0.005 4.81×10 −3 4.79×10 −3 V1 777 122 0.005±0.002 9.52×10 −4 9.34×10 −4 V2 351 122 0.005±0.002 8.50×10 −4 8.32×10 −4 C1 1074 122 0.005±0.002 1.39×10 −3 1.37×10 −3 C2 408 122 0.009±0.002 2.64×10 −3 2.62×10 −3 C3 405 122 0.008±0.002 1.99×10 −3 1.97×10 −3 C4 294 122 0.005±0.002 1.38×10 −3 1.36×10 −3 To reveal further the process shaping TYLCV evolution, we tested whether specific kinds of mutational changes were overrepresented in the natural population of TYLCV. Our results showed that the shifts from C to T, T to C, G to A, and G to T were overrepresented relative to other kinds of mutational changes. These kinds of mutations constituted 62% of the mutations. However, the transversions of A to C, T to A, and C to G appeared to be underrepresented, constituting only about 10% of the mutations (Table 3 3 ). Therefore, these results indicate that transition biases exist in the evolution of TYLCV. These biases favor a shift from CG to AT genome content, especially with a bias in C to T and G to A transitions. 2 Original base 4 Percentage of nucleotide substitution (%) A T C G A 4.19 3.22 7.96 T 3.59 15.6 7.24 C 4.62 26.1 3.86 G 9.84 10.46 3.33 Genetic distances between subpopulations of TYLCV 2 Genetic distance refers to the average number of nucleotide substitutions between two randomly selected sequences in a population, and can be used to evaluate the genetic variation of a virus population. Using the Kimura’s two-parameter method, the mean genetic distances for the intergenic region and individual ORFs ranged from 0.005 to 0.023. Although the individual ORFs of the TYLCV genome presented similar genetic distances, the mean values for the intergenic region showed the highest genetic distances (Table 2 2 ). To understand better the genetic variation of the 26 isolates, the full-length genomic sequences were divided into several subgroups according to the year of collection. Analysis showed that the highest genetic distance was found for the intergenic region (IR) within subgroups (Table 4 4 ). Although we did not observe a quasi-linear increase in genetic diversity across years, calculation of the F st , the coefficient used to evaluate the extent of genetic differentiation or the gene flow between two populations, showed that the values of F st among subgroups 2006 and 2009, 2006 and 2010, as well as 2009 and 2010 were all < 0.33, suggesting frequent gene flow (Table 5 5 ). 2 Genomic region 5 Genetic distance within subgroups 2006 2007 2008 2009 2010 IR 0.007±0.004 0.003±0.002 0.009±0.002 0.018±0.004 0.017±0.004 V1 0.003±0.002 0.001±0.001 0.002±0.001 0.005±0.002 0.004±0.001 V2 0.002±0.002 0.001±0.000 0.002±0.001 0.004±0.002 0.005±0.002 C1 0.004±0.001 0.001±0.000 0.002±0.000 0.004±0.001 0.004±0.001 C2 0.008±0.003 0.003±0.002 0.004±0.001 0.005±0.002 0.007±0.003 C3 0.006±0.003 0.001±0.001 0.003±0.001 0.006±0.002 0.009±0.003 C4 0.003±0.002 0.002±0.001 0.003±0.001 0.004±0.002 0.004±0.002 F st is a coefficient of the extent of genetic differentiation and provides an estimate of the extent of gene flow. An absolute value of F st > 0.33 suggests infrequent gene flow 2 Year 4 F st 2006 2007 2008 2009 2007 0.701 18 2008 0.630 97 0.373 27 2009 0.167 67 0.606 62 0.548 62 2010 0.296 72 0.619 52 0.569 20 0.133 70 Selection constraints acting on TYLCV genes 2 Pairwise genetic differences at the synonymous ( d s ) and nonsynonymous ( d n ) nucleotide positions were estimated according to the Pamilo-Bianchi-Li method. The ratio between nucleotide diversity values in nonsynonymous and synonymous positions ( d n / d s ) (Table 6 6 ) provides an estimation of the degree and direction of the selective constraints acting on the coding regions of TYLCV. Overall, the values of the d n / d s ratio for the V1 , V2 , C1 , and C2 genes were markedly low (0.154–0.44), indicating that these genes are under negative selection. However, the ratios of d n / d s for the C3 and C4 genes were greater than 1. This was probably due to the overlapping regions of C3-C2 and C4-C1, respectively. In the overlapping, frame-shifted C1/C4 region of TYLCV, the first nucleotide position in the C1 codons corresponds to the third position in the C4 codons (C1-1/C4-3), the second position in the C1 codons corresponds to the first position in the C4 codons (C1-2/C4-1), and the third position in the C1 codons corresponds to the second position in the C4 codons (C1-3/C4-2) (Fig. 3a 3a ). Independent evolution of overlapping genes has been described for several viruses, with one of the two overlapping genes being subjected to negative selection, while the other is subjected to positive selection. To evaluate the variation between C1 and C4, entropy values for individual amino acids and nucleotides were determined using the BioEdit software. The amino acid changes in the Rep protein were largely determined by nucleotide mutations at the C1-1/C4-3 sites, whereas those in the C4 protein were primarily caused by mutations at the C1-3/C4-2 sites (Fig. 3b 3b ). Furthermore, mutations at the C1-2/C4-1 sites, which would affect the amino acid sequences of both the Rep and C4 proteins, were rare. Cumulative synonymous and nonsynonymous mutations plotted using SNAP showed that the C1 gene accumulated more synonymous mutations, while the C4 gene accumulated more nonsynonymous mutations (Fig. 3c 3c ). Taken together, these results revealed that the overlapping regions of C1 and C4 are under evolutionary constraints, allowing for more rapid evolution of the C4 protein while minimizing amino acid changes in the conserved Rep protein. d n : nonsynonymous substitution; d s : synonymous substitution. d n and d s are expressed as mean±standard error Gene d n d s d n / d s V1 0.002±0.001 0.013±0.004 0.154 V2 0.002±0.002 0.010±0.005 0.200 C1 0.004±0.001 0.009±0.004 0.440 C2 0.006±0.003 0.023±0.013 0.260 C3 0.009±0.004 0.004±0.004 2.250 C4 0.008±0.003 0.002±0.001 4.000 (a) Frame-shifted positions of the overlapping C1 and C4 genes. Mutations in the open reading frames of the overlapping C1 and C4 genes result in differential rates of amino acid changes. For example, a nucleotide mutation in the third position of the C1 codon is likely to be synonymous, allowing a nonsynonymous mutation in the C4 codon (arrows). (b) Levels of nucleotide and amino acid variation in the three sets of nucleotides in relation to the nucleotide positions of the codons in the overlapping Rep and C4 proteins. (c) Cumulative incidences of synonymous and nonsynonymous mutations in the overlapping region of C1/C4. The x axis represents the position of the codon in C4, and the y axis represents the cumulative values of synonymous or nonsynonymous mutations estimated at a specific codon position Discussion 1 Before 2006, no begomovirus or TYLCV had been detected in Shanghai. Due to the rapid introduction and severe outbreak of the Q biotype of whitefly in China, TYLCV has now moved to several provinces of China, causing unprecedented economic losses, especially in tomato plants. Prior to this field survey, no other competitive begomovirus had been described in Shanghai. This epidemiological situation provides a great opportunity to investigate the genetic structure and evolution of TYLCV following its introduction into a new region. This study analyzed the genetic variability of 26 TYLCV isolates collected from Shanghai over a five-year period. Analysis of the full-length genomic DNA sequences revealed that the level of genetic variation observed in the natural population of TYLCV (1.69×10 −3 subs/site/year) was similar to that found in RNA viruses (Duffy et al., 2008 2008 ; Simmons et al., 2008 2008 ; Pagan and Holmes, 2010 2010 ; Pagan et al., 2010 2010 ). These and other results suggest that a considerable amount of sequence variation occurs within and between populations of geminiviruses, challenging the view that RNA viruses are prone to higher mutation rates than DNA viruses (Domingo and Holland, 1997 1997 ; Roossinck, 1997 1997 ; Isnard et al., 1998 1998 ; Sanz et al., 1999 1999 ). It has been speculated that the high mutation rate of geminivirus genomes is in part due to improper methylation patterns occurring during their replication (Ge et al., 2007 2007 ; Duffy and Holmes, 2008 2008 ). However, recent studies have shown that geminivirus genomes are subjected to DNA methylation in infected host plants (Raja et al., 2008 2008 ; 2010 2010 ; Yang et al., 2011 2011 ; Zhang et al., 2011 2011 ). Thus, it is possible that base-excision repair may not act on geminivirus genomes because the double-stranded state lasts for a very short time during rolling circle replication. The high mutation rate is also probably in part due to base deamination occurring spontaneously or by the action of deaminating host enzymes (Duffy and Holmes, 2008 2008 ; van der Walt et al., 2008 2008 ). This may be supported by the observation of C to T transition biases in this study. Indeed, previous studies have also revealed the same transition biases as our genomic analysis (Ge et al., 2007 2007 ; Duffy and Holmes, 2008 2008 ; 2009 2009 ). Although we did not observe a significant quasi-linear increase in genetic diversity across years, we found frequent gene flow between isolates obtained from the years 2006, 2009, and 2010. As was reported for TYLCV surveyed in an insular environment, the amount of mutations analyzed in this study was greater than that found in tomato yellow leaf curl Sardinia virus (TYLCSV) over an eight-year period, using a single-strand conformation polymorphism technique (Sanchez-Campos et al., 2002 2002 ; Delatte et al., 2007 2007 ). This discrepancy might be explained by the different techniques used. The RCA technique used in this study has revolutionized the diagnosis and genomics of geminiviruses, allowing for a more efficient and reliable way to obtain full-length geminivirus DNA (Haible et al., 2006 2006 ). Our results also indicated that the mutation rate detected in the natural TYLCV population was greater than that caused by phi29 DNA polymerase and sequencing errors. Alternatively, this discrepancy might also come from the viruses used. It seems that TYLCV causes more severe damage to crop production, as surveys conducted between 1996 and 1998 in Málaga showed a progressive displacement of TYLCSV by TYLCV (Sanchez-Campos et al., 1999 1999 ). Indeed, since the first introduction of TYLCV into China, it has become the most prevalent begomovirus in China. Estimation of nucleotide diversities using the ratio between nonsynonymous and synonymous positions ( d n / d s ) suggested that V1, V2, C1, and C2 were more tolerant to amino acid changes than C3 and C4. Among all the coding regions analyzed, the C4 gene embedded in the C1 gene appeared to be the most flexible and might be subjected to positive selection. We demonstrated that the evolution of the overlapping C1 and C4 genes occurs primarily via C1-1/C4-3 and C1-3/C4-2 sites, respectively, and mutations at the C1-2/C4-1 sites, which were most likely to alter the amino acid sequences of both the Rep and C1 proteins, were rare. Further SNAP analysis showed that positions in the C-terminal of the C4 gene, which were under positive selection, were subjected to purifying selection in the overlapping region of the C1 gene. Indeed, a similar phenomenon has been demonstrated for the evolution of tomato leaf deformation virus, an indigenous new world (NW) monopartite begomovirus emerging from the DNA-A component of an NW bipartite progenitor (Melgarejo et al., 2013 2013 ). We speculated that this would be correlated with the function of the C4 protein. Compared to the conserved Rep protein, the functions of the C4 protein are more divergent among monopartite and bipartite begomoviruses. For monopartite begomoviruses, C4 is a pathogenicity determinant that is involved in movement and suppression of host defense responses, whereas the homologues of C4 encoded by most NW begomoviruses do not play a role in pathogenicity. This strategy, which allows viruses to condense a maximal amount of information into short genomes, has been widely used by a number of animal and plant-infecting viruses (Zaaijer et al., 2007 2007 ). Given the rapid spread of TYLCV in China, it will be interesting to monitor the evolutionary constraints in vivo in experimentally infected tomato plants. References S Cohen , Y Antignus . KF Harris . Tomato yellow leaf curl virus, a whitefly-borne geminivirus of tomatoes . Advances in Disease Vector Research , 1994 . : New York Springer-Verlag , 259 - 288 . DOI: 10.1007/978-1-4612-2590-4_10 http://doi.org/10.1007/978-1-4612-2590-4_10 . Tomato yellow leaf curl virus, a whitefly-borne geminivirus of tomatoes H Czosnek , M Ghanim . WMO Thompson . Bemisia tabaci —tomato yellow leaf curl virus interaction causing worldwide epidemics . The Whitefly, Bemisia Tabaci (Homoptera: Aleyrodidae) Interaction with Geminivirus-Infected Host Plants , 2011 . Springer Netherlands , 51 - 67 . DOI: 10.1007/978-94-007-1524-0_3 http://doi.org/10.1007/978-94-007-1524-0_3 . Bemisia tabaci—tomato yellow leaf curl virus interaction causing worldwide epidemics H Delatte , H Holota , B Moury , 等 . Evidence for a founder effect after introduction of tomato yellow leaf curl virus-mild in an insular environment . J Mol Evol , 2007 . 65 ( 1 ): 112 - 118 . 17609843 DOI: 10.1007/s00239-007-0005-x http://doi.org/10.1007/s00239-007-0005-x . Evidence for a founder effect after introduction of tomato yellow leaf curl virus-mild in an insular environment E Domingo , JJ Holland . RNA virus mutations and fitness for survival . Annu Rev Microbiol , 1997 . 51 ( 1 ): 151 - 178 . 9343347 DOI: 10.1146/annurev.micro.51.1.151 http://doi.org/10.1146/annurev.micro.51.1.151 . RNA virus mutations and fitness for survival S Duffy , EC Holmes . Phylogenetic evidence for rapid rates of molecular evolution in the single-stranded DNA begomovirus tomato yellow leaf curl virus . J Virol , 2008 . 82 ( 2 ): 957 - 965 . 17977971 DOI: 10.1128/JVI.01929-07 http://doi.org/10.1128/JVI.01929-07 . Phylogenetic evidence for rapid rates of molecular evolution in the single-stranded DNA begomovirus tomato yellow leaf curl virus S Duffy , EC Holmes . Validation of high rates of nucleotide substitution in geminiviruses: phylogenetic evidence from East African cassava mosaic viruses . J Gen Virol , 2009 . 90 ( 6 ): 1539 - 1547 . 19264617 DOI: 10.1099/vir.0.009266-0 http://doi.org/10.1099/vir.0.009266-0 . Validation of high rates of nucleotide substitution in geminiviruses: phylogenetic evidence from East African cassava mosaic viruses S Duffy , LA Shackelton , EC Holmes . Rates of evolutionary change in viruses: patterns and determinants . Nat Rev Genet , 2008 . 9 ( 4 ): 267 - 276 . 18319742 DOI: 10.1038/nrg2323 http://doi.org/10.1038/nrg2323 . Rates of evolutionary change in viruses: patterns and determinants L Ge , J Zhang , X Zhou , 等 . Genetic structure and population variability of tomato yellow leaf curl China virus . J Virol , 2007 . 81 ( 11 ): 5902 - 5907 . 17376922 DOI: 10.1128/JVI.02431-06 http://doi.org/10.1128/JVI.02431-06 . Genetic structure and population variability of tomato yellow leaf curl China virus E Glick , A Zrachya , Y Levy , 等 . Interaction with host SGS3 is required for suppression of RNA silencing by tomato yellow leaf curl virus V2 protein . PNAS , 2008 . 105 ( 1 ): 157 - 161 . 18165314 DOI: 10.1073/pnas.0709036105 http://doi.org/10.1073/pnas.0709036105 . Interaction with host SGS3 is required for suppression of RNA silencing by tomato yellow leaf curl virus V2 protein W Guo , T Jiang , X Zhang , 等 . Molecular variation of satellite DNAβ molecules associated with Malvastrum yellow vein virus and their role in pathogenicity . Appl Environ Microbiol , 2008 . 74 ( 6 ): 1909 - 1913 . 18245261 DOI: 10.1128/AEM.02461-07 http://doi.org/10.1128/AEM.02461-07 . Molecular variation of satellite DNAβ molecules associated with Malvastrum yellow vein virus and their role in pathogenicity W Guo , X Yang , Y Xie , 等 . Tomato yellow leaf curl Thailand virus-[Y72] from Yunnan is a monopartite begomovirus associated with DNAβ . Virus Genes , 2009 . 38 ( 2 ): 328 - 333 . 19165588 DOI: 10.1007/s11262-009-0327-4 http://doi.org/10.1007/s11262-009-0327-4 . Tomato yellow leaf curl Thailand virus-[Y72] from Yunnan is a monopartite begomovirus associated with DNAβ D Haible , S Kober , H Jeske . Rolling circle amplification revolutionizes diagnosis and genomics of geminiviruses . J Virol Methods , 2006 . 135 ( 1 ): 9 - 16 . 16513183 DOI: 10.1016/j.jviromet.2006.01.017 http://doi.org/10.1016/j.jviromet.2006.01.017 . Rolling circle amplification revolutionizes diagnosis and genomics of geminiviruses T Hall . BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT . Nucl Acids Symp Ser , 1999 . 41 95 - 98 . . BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT L Hanley-Bowdoin , SB Settlage , BM Orozco , 等 . Geminiviruses: models for plant DNA replication, transcription, and cell cycle regulation . Crit Rev Biochem Mol Biol , 1999 . 18 ( 1 ): 71 - 106 . DOI: 10.1080/07352689991309162 http://doi.org/10.1080/07352689991309162 . Geminiviruses: models for plant DNA replication, transcription, and cell cycle regulation M Isnard , M Granier , R Frutos , 等 . Quasispecies nature of three maize streak virus isolates obtained through different modes of selection from a population used to assess response to infection of maize cultivars . J Gen Virol , 1998 . 79 ( 12 ): 3091 - 3099 . 9880027 . Quasispecies nature of three maize streak virus isolates obtained through different modes of selection from a population used to assess response to infection of maize cultivars GM Jenkins , A Rambaut , OG Pybus , 等 . Rates of molecular evolution in RNA viruses: a quantitative phylogenetic analysis . J Mol Evol , 2002 . 54 ( 2 ): 156 - 165 . 11821909 DOI: 10.1007/s00239-001-0064-3 http://doi.org/10.1007/s00239-001-0064-3 . Rates of molecular evolution in RNA viruses: a quantitative phylogenetic analysis L Jiang , CH Wei , Y Li . Viral suppression of RNA silencing . Sci China C Life Sci , 2012 . 55 ( 2 ): 109 - 118 . DOI: 10.1007/s11427-012-4279-x http://doi.org/10.1007/s11427-012-4279-x . Viral suppression of RNA silencing B Korber . AG Rodrigo , JR Learn , H Gerald . HIV signature and sequence similarities . Computational and Evolutionary Analysis of HIV Molecular Sequences , 2002 . : Boston, Massachusetts Kluwer Academic Publishers , 55 - 72 . . HIV signature and sequence similarities P Librado , J Rozas . DnaSP v5: a software for comprehensive analysis of DNA polymorphism data . Bioinformatics , 2009 . 25 ( 11 ): 1451 - 1452 . 19346325 DOI: 10.1093/bioinformatics/btp187 http://doi.org/10.1093/bioinformatics/btp187 . DnaSP v5: a software for comprehensive analysis of DNA polymorphism data TA Melgarejo , T Kon , MR Rojas , 等 . Characterization of a new world monopartite begomovirus causing leaf curl disease of tomato in Ecuador and Peru reveals a new direction in geminivirus evolution . J Virol , 2013 . 87 ( 10 ): 5397 - 5413 . 23468482 DOI: 10.1128/JVI.00234-13 http://doi.org/10.1128/JVI.00234-13 . Characterization of a new world monopartite begomovirus causing leaf curl disease of tomato in Ecuador and Peru reveals a new direction in geminivirus evolution E Moriones , J Navas-Castillo . Tomato yellow leaf curl virus, an emerging virus complex causing epidemics worldwide . Virus Res , 2000 . 71 ( 1-2 ): 123 - 134 . 11137167 DOI: 10.1616/S0168-1702(00)00193-3 http://doi.org/10.1616/S0168-1702(00)00193-3 . Tomato yellow leaf curl virus, an emerging virus complex causing epidemics worldwide M Padidam , RN Beachy , CM Fauquet . Classification and identification of geminiviruses using sequence comparisons . J Gen Virol , 1995 . 76 ( Pt 2 ): 249 - 263 . 7844548 DOI: 10.1099/0022-1317-76-2-249 http://doi.org/10.1099/0022-1317-76-2-249 . Classification and identification of geminiviruses using sequence comparisons I Pagan , EC Holmes . Long-term evolution of the Luteoviridae : time scale and mode of virus speciation . J Virol , 2010 . 84 ( 12 ): 6177 - 6187 . 20375155 DOI: 10.1128/JVI.02160-09 http://doi.org/10.1128/JVI.02160-09 . Long-term evolution of the Luteoviridae: time scale and mode of virus speciation I Pagan , C Firth , EC Holmes . Phylogenetic analysis reveals rapid evolutionary dynamics in the plant RNA virus genus tobamovirus . J Mol Evol , 2010 . 71 ( 4 ): 298 - 307 . 20838783 DOI: 10.1007/s00239-010-9385-4 http://doi.org/10.1007/s00239-010-9385-4 . Phylogenetic analysis reveals rapid evolutionary dynamics in the plant RNA virus genus tobamovirus P Raja , BC Sanville , RC Buchmann , 等 . Viral genome methylation as an epigenetic defense against geminiviruses . J Virol , 2008 . 82 ( 18 ): 8997 - 9007 . 18596098 DOI: 10.1128/JVI.00719-08 http://doi.org/10.1128/JVI.00719-08 . Viral genome methylation as an epigenetic defense against geminiviruses P Raja , JN Wolf , DM Bisaro . RNA silencing directed against geminiviruses: post-transcriptional and epigenetic components . Biochim Biophys Acta , 2010 . 1799 ( 3 ): 337 - 351 . 20079472 DOI: 10.1016/j.bbagrm.2010.01.004 http://doi.org/10.1016/j.bbagrm.2010.01.004 . RNA silencing directed against geminiviruses: post-transcriptional and epigenetic components MJ Roossinck . Mechanisms of plant virus evolution . Annu Rev Phytopathol , 1997 . 35 ( 1 ): 191 - 209 . 15012521 DOI: 10.1146/annurev.phyto.35.1.191 http://doi.org/10.1146/annurev.phyto.35.1.191 . Mechanisms of plant virus evolution S Sanchez-Campos , J Navas-Castillo , R Camero , 等 . Displacement of tomato yellow leaf curl virus (TYLCV)-Sr by TYLCV-Is in tomato epidemics in Spain . Phytopathology , 1999 . 89 ( 11 ): 1038 - 1043 . 18944659 DOI: 10.1094/PHYTO.1999.89.11.1038 http://doi.org/10.1094/PHYTO.1999.89.11.1038 . Displacement of tomato yellow leaf curl virus (TYLCV)-Sr by TYLCV-Is in tomato epidemics in Spain S Sanchez-Campos , JA Diaz , F Monci , 等 . High genetic stability of the Begomovirus tomato yellow leaf curl Sardinia virus in southern Spain over an 8-year period . Phytopathology , 2002 . 92 ( 8 ): 842 - 849 . 18942962 DOI: 10.1094/PHYTO.2002.92.8.842 http://doi.org/10.1094/PHYTO.2002.92.8.842 . High genetic stability of the Begomovirus tomato yellow leaf curl Sardinia virus in southern Spain over an 8-year period AI Sanz , A Fraile , JM Gallego , 等 . Genetic variability of natural populations of cotton leaf curl geminivirus, a single-stranded DNA virus . J Mol Evol , 1999 . 49 ( 5 ): 672 - 681 . 10552048 DOI: 10.1007/PL00006588 http://doi.org/10.1007/PL00006588 . Genetic variability of natural populations of cotton leaf curl geminivirus, a single-stranded DNA virus HE Simmons , EC Holmes , AG Stephenson . Rapid evolutionary dynamics of zucchini yellow mosaic virus . J Gen Virol , 2008 . 89 ( 4 ): 1081 - 1085 . 18343852 DOI: 10.1099/vir.0.83543-0 http://doi.org/10.1099/vir.0.83543-0 . Rapid evolutionary dynamics of zucchini yellow mosaic virus K Tamura , D Peterson , N Peterson , 等 . MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods . Mol Biol Evol , 2011 . 28 ( 10 ): 2731 - 2739 . 21546353 DOI: 10.1093/molbev/msr121 http://doi.org/10.1093/molbev/msr121 . MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods JD Thompson , TJ Gibson , DG Higgins . Multiple sequence alignment using ClustalW and ClustalX . Curr Prot Bioinf , 2002 . 00 2.3.1 - 2.3.22 . DOI: 10.1002/0471250953.bi0203s00 http://doi.org/10.1002/0471250953.bi0203s00 . Multiple sequence alignment using ClustalW and ClustalX E van der Walt , DP Martin , A Varsani , 等 . Experimental observations of rapid Maize streak virus evolution reveal a strand-specific nucleotide substitution bias . Virol J , 2008 . 5 104 18816368 DOI: 10.1186/1743-422X-5-104 http://doi.org/10.1186/1743-422X-5-104 . Experimental observations of rapid Maize streak virus evolution reveal a strand-specific nucleotide substitution bias JB Wu , FM Dai , XP Zhou . First report of tomato yellow leaf curl virus in China . Plant Dis , 2006 . 90 ( 10 ): 1359 - 1359 . DOI: 10.1094/Pd-90-1359c http://doi.org/10.1094/Pd-90-1359c . First report of tomato yellow leaf curl virus in China X Yang , Y Xie , P Raja , 等 . Suppression of methylation-mediated transcriptional gene silencing by βC1-SAHH protein interaction during geminivirus-betasatellite infection . PLoS Pathog , 2011 . 7 ( 10 ): e1002329 22028660 DOI: 10.1371/journal.ppat.1002329 http://doi.org/10.1371/journal.ppat.1002329 . Suppression of methylation-mediated transcriptional gene silencing by βC1-SAHH protein interaction during geminivirus-betasatellite infection HL Zaaijer , FJ van Hemert , MH Koppelman , 等 . Independent evolution of overlapping polymerase and surface protein genes of hepatitis B virus . J Gen Virol , 2007 . 88 ( Pt 8 ): 2137 - 2143 . 17622615 DOI: 10.1099/vir.0.82906-0 http://doi.org/10.1099/vir.0.82906-0 . Independent evolution of overlapping polymerase and surface protein genes of hepatitis B virus H Zhang , H Gong , X Zhou . Molecular characterization and pathogenicity of tomato yellow leaf curl virus in China . Virus Genes , 2009 . 39 ( 2 ): 249 - 255 . 19590945 DOI: 10.1007/s11262-009-0384-8 http://doi.org/10.1007/s11262-009-0384-8 . Molecular characterization and pathogenicity of tomato yellow leaf curl virus in China Z Zhang , H Chen , X Huang , 等 . BSCTV C2 attenuates the degradation of SAMDC1 to suppress DNA methylation-mediated gene silencing in Arabidopsis . Plant Cell , 2011 . 23 ( 1 ): 273 - 288 . 21245466 DOI: 10.1105/tpc.110.081695 http://doi.org/10.1105/tpc.110.081695 . BSCTV C2 attenuates the degradation of SAMDC1 to suppress DNA methylation-mediated gene silencing in Arabidopsis XP Zhou , Y Xie , ZK Zhang . Molecular characterization of a distinct begomovirus infecting tobacco in Yunnan, China . Arch Virol , 2001 . 146 ( 8 ): 1599 - 1606 . 11676420 DOI: 10.1007/s007050170081 http://doi.org/10.1007/s007050170081 . Molecular characterization of a distinct begomovirus infecting tobacco in Yunnan, China Compliance with ethics guidelines Xiu-ling YANG, Meng-ning ZHOU, Ya-juan QIAN, Yan XIE, and Xue-ping ZHOU declare that they have no conflict of interest. This article does not contain any studies with human or animal subjects performed by any of the authors. 10.1631/jzus.B1300110 Project supported by the Special Fund for Agro-Scientific Research in the Public Interest from the Ministry of Agriculture, China (No. 201003065), the China Postdoctoral Science Foundation (No. 2012M510200), and the National Natural Science Foundation of China (No. 31200117) 2013-04-11 2014-01-19 2013-08-21 2022-01-05 Copyright © Zhejiang University and Springer-Verlag Berlin Heidelberg 2014 2014 Journal of Zhejiang University-SCIENCE B (Biomedicine and Biotechnology) 02 15