Figures
Abstract
Begomovirus (genus Begomovirus, family Geminiviridae) infection is devastating to a wide variety of agricultural crops including tomato, squash, and cassava. Thus, understanding the replication and adaptation of begomoviruses has important translational value in alleviating substantial economic loss, particularly in developing countries. The bipartite genome of begomoviruses prevalent in the New World and their counterparts in the Old World share a high degree of genome homology except for a partially overlapping reading frame encoding the pre-coat protein (PCP, or AV2). PCP contributes to the essential functions of intercellular movement and suppression of host RNA silencing, but it is only present in the Old World viruses. In this study, we analyzed a set of non-redundant bipartite begomovirus genomes originating from the Old World (N = 28) and the New World (N = 65). Our bioinformatic analysis suggests ∼120 nucleotides were deleted from PCP’s proximal promoter region that may have contributed to its loss in the New World viruses. Consequently, genomes of the New World viruses are smaller than the Old World counterparts, possibly compensating for the loss of the intercellular movement functions of PCP. Additionally, we detected substantial purifying selection on a portion of the New World DNA-B movement protein (MP, or BC1). Further analysis of the New World MP gene revealed the emergence of a putative tyrosine phosphorylation site, which likely explains the increased purifying selection in that region. These findings provide important information about the strategies adopted by bipartite begomoviruses in adapting to new environment and suggest future in planta experiments.
Citation: Ho ES, Kuchie J, Duffy S (2014) Bioinformatic Analysis Reveals Genome Size Reduction and the Emergence of Tyrosine Phosphorylation Site in the Movement Protein of New World Bipartite Begomoviruses. PLoS ONE 9(11): e111957. https://doi.org/10.1371/journal.pone.0111957
Editor: Darren P. Martin, Institute of Infectious Disease and Molecular Medicine, South Africa
Received: August 18, 2014; Accepted: October 9, 2014; Published: November 10, 2014
Copyright: © 2014 Ho et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: The authors confirm that all data underlying the findings are fully available without restriction. Genomes and protein sequences used in this study are included in the Supporting information section or through the web link: (http://sites.lafayette.edu/hoe/files/2014/01/bipartite_seqs_eh_jk_sd.tar_.gz).
Funding: ESH is supported by NIH K12 GM093854-01. JK is supported by New Jersey City University and the RiSE program at Rutgers. SD is supported by NSF DEB 1026095 and BMGF/DFID OPP1052391. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Begomoviruses (genus Begomovirus, family Geminiviridae) are single-stranded DNA viruses of dicots with small genomes - one or two circular segments of ∼2.5–2.9 K nucleotides (nts). Begomoviruses are transmitted by the whitefly Bemisia tabaci [1], [2] and their damaging infections pose a severe threat to commercial and subsistence production of key crops worldwide, including tomato, squash, cassava and bean [3]. Understanding the molecular biology and adaptation of begomoviruses to novel hosts has an important socioeconomic impact as they are emerging problems in developing countries [3]. The vast majority of begomovirus sequences also exhibit a classic biogeographic pattern: they fall into clades of New World (the Americas, and Caribbean) and Old World (rest of the world) viruses, with New World viruses thought to be derived from those in the Old World [4]–[6]. Bipartite begomoviruses, which have two similarly-sized, ambisense genomic segments termed DNA-A and DNA-B, are found worldwide, with monopartite begomoviruses largely restricted to the Old World [7]. The DNA-A segment contains five or six genes, including the capsid protein (CP, also known as AV1), the replication-associated protein (REP, also known as AC1), a transcriptional activator (TrAP, also known as AC2), a replication enhancer (REn, also known as AC3) that overlaps with both the REP and TrAP genes and a virulence factor (AC4) that overlaps the reading frame within REP. The DNA-B segment contains two non-overlapping genes: the nuclear shuttle protein (NSP, also known as BV1), and the movement protein (MP, also known as BC1).
Old and New World bipartite begomoviruses share a high degree of homology, with the largest exception being the gene for the pre-coat protein (PCP, also known as AV2), which partially overlaps the CP gene and is only present in Old World viruses [8]. PCP and the monopartite V2 has been shown to localize at the cell periphery and is thought to act as a “movement protein” by increasing the size exclusion limit of the plasmodesmata [9], [10]. They also suppresses RNA silencing by binding to the host’s SGS3 protein [11]. V2 is thought to be the key movement protein in monopartite Old World viruses, but two genes on the DNA-B segment (NSP and MP) also contribute to systemic infection of plants by bipartite begomoviruses [12]. Virulent New World begomoviruses must rely on their other seven proteins to cope with the loss of PCP, and this is frequently invoked as the reason the DNA-B segment is required for infectivity of the overwhelming majority of New World begomoviruses [13]. Despite this assumption, the selective pressures imposed by the loss of PCP on the remaining New World viral genes have not been examined.
In this report we have compared the genome size, degree of variability and purifying selection of the viral genes between the Old and New World. Results indicate a loss of 100 nts in PCP’s promoter region, stronger purifying selection on the two DNA-B genes in the New World, and the emergence of a putative tyrosine phosphorylation site in the New World MP. Studies with RNA plant viruses have shown that phosphorylation of MP regulates their localization and may account for cell-to-cell movement [14]. We speculate that the reduction in viral cell-to-cell movement caused by the loss of the PCP in the New World begomoviruses may be compensated by systematic genome size reduction and/or the gain of additional phosphorylation activity in the MP.
Materials and Methods
Compilation of bipartite begomovirus genomes
Genomes of begomovirus were downloaded from the June-2012 release of the viral genome database hosted in NCBI (ftp://ftp.ncbi.nih.gov/refseq/release/viral/). Only genomes containing the distinct, invariant nonamer “TAATATT|AC” were included in this study (the vertical bar represents the cleavage site). The pairing of DNA-A and DNA-B genomes, and the classification of genomes into Old and New Worlds were done semi-automatically according to the information stated in NCBI’s RefSeq records [15] and ICTV report [16]. To ease sequence comparison, the beginning of the cleavage site “AC” was adopted as reference position 1 and the original genomic coordinates stated in NCBI’s RefSeq records of the begomoviruses were adjusted accordingly. 33 and 83 Old and New World bipartite begomoviruses were collected, respectively, before further redundancy checking.
Identification of common regions
The DNA-A and DNA-B genomes of a bipartite begomovirus share a 200- to 250-nt long highly identical segment (>85%), namely the common region (CR), in which the invariant nonamer “TAATATT|AC” resides near to the middle of it. To determine the 5′ and 3′ termini of the CR, a pair of segments consisting of 250 nts upstream and downstream flanking regions of the invariant nonamer from DNA-A and DNA-B was aligned. Based on the alignment, the longest stretch of highly identical (at least 20 nts long and 80% identity) segment flanking the invariant nonamer was taken as the CR.
Identification of non-redundant genomes and ORFs
We clustered DNA-As together if their CRs shared >80% similarity. If more than one species was found in a cluster, only one species was retained arbitrarily for further analysis. A Peruvian begomovirus, Tomato leaf deformation virus (ToLDeV), was confirmed to be the first New World monopartite begomovirus in 2013 [7], but this was after our dataset had been finalized. ToLDeV does not appear to have a PCP gene. As a result, 26 out of 33 (85%) and 65 out of 83 (78%) non-redundant Old World and New World bipartite begomovirus genomes were included in this analysis (Table 1). The full list of bipartite begomovirus genomes used and their sizes can be found in the Table S1 in File S1. Additionally, ORFs specified in RefSeq records were verified. We required the stated coding sequences or ORFs be translated exactly to the protein sequences specified in the RefSeq records. Genes failed to meet this requirement were excluded from this study (Table 1). Genomes and viral protein sequences used in this study can be downloaded as Data S1 or through this web link: http://sites.lafayette.edu/hoe/files/2014/01/bipartite_seqs_eh_jk_sd.tar_.gz.
dN/dS calculation
dN/dS represents the log ratio of the rate of non-synonymous substitutions to the rate of synonymous (silent) substitutions. A negative, zero, or positive dN/dS value indicates purifying (negative), neutral, or positive selection, respectively. Protein sequences were aligned by T_COFFEE [17] using default parameters. Protein alignments were converted to codon alignments using pal2nal v14 [18]. The codon alignments were submitted to the tool SLAC [19] hosted in the Datamonkey web server http://www.datamonkey.org/ [20] for site dN/dS calculation. Substitution models were selected by iterating the likelihood ratio tests between nested and non-nested models. This procedure is implemented in Datamonkey web server and detailed discussion of the procedure can be found in [21]. The results calculated by SLAC were downloaded in CSV format for analysis.
Pairwise protein sequence alignment
As dynamic programming approach to local pairwise sequence alignment produces the optimal alignment for a given scoring scheme. We used the percentage of identity calculated by an implementation of such approach i.e. Smith-Waterman water program [22], to determine the diversity of each viral protein for either Old or New World regions. BLOSUM62 score matrix was used and gap opening penalty and gap extension penalty were 10 and 0.5, respectively.
D-statistic of the Kolmogorov-Smirnov test
In order to ascertain the statistical significance of the difference between two non-Gaussian, cumulative distributions of protein sequence similarities and dN/dS values, we quantified the difference using the D-statistic of the two-sample Kolmogorov-Smirnov (KS) test. In both worlds, the viral protein AC4 exhibited the highest diversity. Thus, AC4 was chosen as the reference for two-sample D-statistic calculation. D-statistics were computed using the R function ks. test () [23]. All the values of D-statistic calculated showed significant differences between the two worlds with p-value in the range of 10−16.
Scanning of functional sites in the movement protein
We developed a Python script (available upon request) to scan for functional sites in protein sequences using the BioPython scanProsite package [24], where the option for skipping of high probability of occurrence was turned off. In addition, our script used the bootstrap approach to compute the p-value of hits through these steps: 1. Obtain the list of functional sites detected in the input sequences through scanProsite, 2. Scramble input sequences, 3. Scan for functional sites in scrambled sequences, 4. Register the list of functional sites found in scrambled sequences, 5. Repeat steps 2 to 4 100 times (a user-defined parameter), 6. Estimate the p-value of a functional site by dividing the occurrence of the functional site in scrambled sequences by the occurrence of the same site in the original input sequences.
Results and Discussion
New World begomoviruses have smaller segments
We discovered that the genome size of DNA-A in the New World is on average 121 nts shorter than their counterparts in the Old World (Figure 1C). Intriguingly, though no apparent gene loss event was reported previously in the New World DNA-B, their genomes (mean size is 2,589 nts, standard deviation or s.d. 43) are also on average 113 nts smaller than the Old World DNA-B (mean is 2,702 nts and s.d. 55) as shown in Figure 1C. This commensurate genome size reduction does not seem to be coincidental as bipartite genome segments (DNA-A and DNA-B) in the New World begomoviruses show a higher correlation in size (R = 0.91, p-value <2.2×10−16) than those in the Old World (R = 0.74, p-value <8.1×10−6). Besides, the genome size differences between DNA-A and DNA-B concurred this point as we found smaller and less variable differences between DNA-A and DNA-B in the New World (mean 37 nts, s.d. 18) than the Old World (mean 45 nts, s.d. 37). Regardless of the geographical factor, this result may suggest size codependency of the bipartite genomes, which is still largely unknown. Our findings are unlikely confounded by biased samples as viral genes AC4 (pink) and REP (blue) exhibit similar spectra of sequence diversity between the two worlds (Figure 2). We further investigated whether or not deletions are localized at a particular region and how it may explain the loss of PCP in the New World begomoviruses.
A) C+G profiles of the homologous regions upstream from the CP from Old and New World bipartite begomoviruses. All positions labeled in the gene structure diagram are average values. Each plot represents the average number of CC, CG, GC, or GG in a 60-nt window. B) Distributions of the distance between the cleavage site “AC” of the invariant nonamer and the beginning of CP gene. C) Distribution of genome size.
A) Cumulative distributions of percentage of identity (%id) of viral proteins from the Old World bipartite begomoviruses. The D values printed beside the protein name in the legend represent the magnitude of deviation of the plot from the AC4’s curve and it was determined by two-sample Kolmogorov test. Larger the D value, the great is the deviation from AC4. B) Proteins from the New World.
Deletions are localized at PCP’s promoter region
We compared dinucleotide profiles in the 400-nt upstream, homologous region of all viral genes using a 60-nt sliding window between the two worlds (see Materials and Methods, and Figure S1A–G in File S1). We found that dinucleotides were better than single nucleotides in insulating the profiles from random nucleotide fluctuation. If short (<5 nts) insertions or deletions are scattered, dinucleotide profiles between the two worlds should exhibit similar patterns; otherwise we should see a direct shift between the two profiles. Among dinucleotide profiles of all genes, only C+G profiles, i.e. CC, CG, GC, and CG, of the CP gene were found to differ between the two worlds in which the region with high concentration of C+G in the New World was shifted ∼100 nts closer to the start of the ORF (Figure 1A). The elevated C+G content is chiefly due to the stem of the highly conserved hairpin structure found in all begomoviruses in which the loop region contains the invariant “TAATATT|AC” nonamer (“|” represents the cleavage site during complementary strand synthesis). Corroborating results were found when we examined the distance between the cleavage site “AC” and the start of the CP genes in both worlds (Figure 1B) where the New World’s CP gene is on average 100 nts closer to the cleavage site “AC” than those in the Old World. This accounts for much of the 121 nts shorter average genome size of DNA-A of the New World viruses compared to those of the Old World (Figure 1C). In New World begomoviruses the distance from the cleavage site to CP is highly correlated with DNA-A size (R = 0.93), but to a lesser extent in the Old World (R = 0.72). Additionally, the C+G content in the non-overlapping region of PCP (from −164 to 0 in Figure 1A) remains at similar level between the two worlds. Therefore our analysis indicates one or more deletions totaling more than 100 nts were mainly localized in the proximal promoter region of PCP, not in other genomic regions, and that these deletions may have led to PCP inactivation in the New World begomoviruses.
Currently, little is known about the effect of genome size on cell-to-cell transport through plasmodesma but studies have shown that the plasmodesmata impose a size limit [25], [26]. Effective shuttling of viral genomes between cells without passing through the cell wall is critical for maintaining infectivity of plant viruses as small viruses do not encode enzymes to breakdown the cell wall, which other phytopathogens such as fungi employ [27]. This finding suggests genome size reduction may be one of the evolutionary paths selected for in the New World begomoviruses in order to maintain virulence despite the loss of cell-to-cell movement conferred by PCP.
New World NSP and MP are under enhanced purifying selection
The compact begomovirus genome encodes only a small number of highly overlapping genes in ambisense, most known to have multiple functions during infection. The lost functions of the PCP gene are likely compensated by remaining genes in New World begomoviruses. Therefore, we took a comparative approach to identify the presence of purifying selection in the New World viral proteins. We measured the within-world diversity of each gene by pairwise protein sequence alignment. A high sequence similarity indicates strong conservation pressures on the genes. Figure 2 shows the cumulative distributions of pairwise identity (%id) of seven or eight viral proteins from the two worlds. In the Old World (Figure 2A), AC4 exhibits the highest variability followed by the two DNA-B proteins NSP and MP, then REn, TrAP, REP and finally CP. The CP is known to be the most conserved of all begomovirus proteins, and under the greatest amount of purifying selection [12], [28]. All distributions were tested for statistical significance (p-value ∼10−26) according to the two-sample Kolmogorov-Smirnov test with AC4 as the reference protein. Results from the New World viral proteins show much lower levels of diversity except for AC4 and REP. Such results are consistent with the presumed more recent origin of New World begomoviruses [4], but show a different pattern in protein variability (Figure 2B). The seven proteins are still bounded by AC4 as the most variable protein, and CP being the most conserved. The most striking difference is in the reduced variability of the New World MP (black plot in Figure 2B), which has become the most conserved protein after CP. Additionally, the New World NSP (gray plot in Figure 2B) is also found to show significant reduction in variability, comparable to the essential replication protein REP (blue plot in Figure 2B). It appears that both DNA-B genes are under stronger selective pressure in the absence of PCP. These results suggest that the less genomically compressed DNA-B genomic segment was more able to accommodate new or enhanced functions than the more constrained DNA-A segment, which already has several overlapping open reading frames.
To further confirm this point, we sought evidence for adaptive evolution at the nucleotide level to corroborate these protein sequence analyses. Adaptive evolution is measured by the log ratio of the rate of non-synonymous substitutions versus synonymous substitutions (dN/dS). If the rate of non-synonymous substitution is lower than synonymous substitution, dN/dS will yield a negative value, indicating amino acid substitution is unfavorable. Conversely, a positive dN/dS value indicates amino acid substitution is permissible, suggesting the protein is under positive or adaptive selection. Aligned protein sequences were converted to corresponding codon alignments before dN/dS calculation using the Single Likelihood Ancestor Counting (SLAC) method from the Datamonkey website [20], [29]. Figure 3A–B show the cumulative dN/dS ratios for MP and NSP. The New World MP shows the biggest deviation (D = 0.54, p-value = 0) from the Old World counterpart and nearly all dN/dS values fall in the negative region, reconfirming elevated purifying selection in the New World MP. But we did not see this in other viral genes (Figure S2 in File S1). We further explored whether or not purified residues in the New World MP constitute to any functional motif(s).
D-value, with p-value, represents the deviation of the New World curve from the Old World curve. D-value was calculated by the two-sample Kolmogorov test. Averaged site dN/dS is displayed on the top left, which reflects the overall selection pressure on the protein. A) MP. B) NSP. C) The emergence of putative tyrosine phosphorylation site in the New World MP. Consensus sequence pictures were created using Weblogo [34]. Searching is based on PROSITE database [35]. PROSITE ID of the tyrosine phosphorylation site is PS00007 where its consensus is [RK]-x(2,3)-[DE]-x(2,3)-Y. D) Codon alignment of the homologous region of the site.
The emergence of tyrosine phosphorylation site in the New World MP
In order to uncover the specific nucleotides subjected to elevated purifying selection in the New World MP, we compared the functional sites of MP in both worlds. Among all functional sites discovered, a putative tyrosine phosphorylation site [RK]-x(2,3)-[DE]-x(2,3)-Y (PROSITE ID: PS00007; the notation means the site starts with either R or K, followed by any 2 to 3 residues, and then a D or E residue, followed by any 2 to 3 residues, and ends with Y) shows the greatest difference between the two worlds: 62/65 New World MP sequences were found to carry the tyrosine phosphorylation site compared to only 1/28 from the Old World.
We also evaluated the likelihood for the emergence of this site by comparing the putative tyrosine phosphorylation site in the New World MP with the homologous region in the Old World MP. We aligned the eight-residue site and the corresponding codons. The amino acid consensus of the Old World (Figure 3C) shows only two residue substitutions are needed to transform the functionally indeterminate eight-residue site in the Old World MP to the tyrosine phosphorylation site found in the New World MP. From the codon perspective, four nucleotide substitutions from the first and fifth codons are sufficient to transform the site (Figure 3D).
Conclusions
Our analysis strongly suggests one or more deletions of 100 nts in the promoter region of the New World PCP, which may be linked to the inactivation of PCP. The resultant shrunken genome may have had an advantage in cell-to-cell movement through plasmodesmata. Furthermore, our genome size analysis unraveled putative size codependency of the bipartite genomes. As the New World begomoviruses are presumably originated from the Old World counterparts recently [4], the conspicuous correlation between the New World DNA-A and DNA-B could be alluded to bottleneck effect. However, the more diverse Old World begomoviruses still maintain a high level of correlation (R = 0.74, p-value <8.1×10−6) between segment size of their bipartite genomes. DNA-B’s functional sequences – the common region (∼200 nts), ORFs of NSP (∼800 nts) and MP (∼900 nts) and their promoter regions (∼200 nts) – occupy ∼2,100 nts of the 2,700-nt genome on average, leaving ∼600 nts (22% of the genome) available for size reduction. According to our data (Table S1 in File S1), the mean, median and maximum difference between the bipartite genomes of the Old World viruses are only 45, 40, and 170 nts, respectively, which are far smaller than the 600 nts permissible range without interrupting the genomic structure of the viruses. This result is surprising as begomoviruses are fast mutating [30] and recombining [31], [32] ssDNA viruses, indicating the presence of unknown constraints that limit the variance in size between segments in bipartite genomes.
Our prediction aligns with findings in closely related monopartite begomoviruses. PCP has been reported previously to perform some MP functions, such as intracellular movement and cell periphery localization [9], [25], [26]. Additionally, in-vitro phosphorylation activity was reported in MP of Abutilon mosaic virus [33]. Our thorough bioinformatic comparison of geographically separated begomovirus species has produced a candidate region for detailed wet lab analysis. If the tyrosine phosphorylation site is critical to infectivity of New World begomoviruses, it will be a novel target for sequence-specific, anti-viral strategies.
Supporting Information
File S1.
Contains the following files: Table S1: Selected bipartite begomoviruses and their genome size. Figure S1: Dinucleotide profiles in 400-nt upstream regions. Window size is 60 nts. Y-axis denotes the average occurrences of the specified dinucleotide in the 60-nt window. Plots of dinucleotides AA, AC, …, TG, TT are arranged from top left to bottom right. Figure S2: Cumulative dN/dS values by gene.
https://doi.org/10.1371/journal.pone.0111957.s001
(DOC)
Data S1.
Genomes and protein sequences used in this study.
https://doi.org/10.1371/journal.pone.0111957.s002
(TAR.GZ)
Author Contributions
Conceived and designed the experiments: ESH JK SD. Performed the experiments: ESH JK. Analyzed the data: ESH SD. Contributed reagents/materials/analysis tools: ESH JK. Contributed to the writing of the manuscript: ESH SD.
References
- 1. Nault LR (1997) Arthropod transmission of plant viruses: A new synthesis. Annals of the Entomological Society of America 90: 521–541.
- 2. Zhang WM, Fu HB, Wang WH, Piao CS, Tao YL, et al. (2014) Rapid Spread of a Recently Introduced Virus (Tomato Yellow Leaf Curl Virus) and Its Vector Bemisia tabaci (Hemiptera: Aleyrodidae) in Liaoning Province, China. Journal of Economic Entomology 107: 98–104.
- 3.
Seal SE, Jeger MJ, Van den Bosch F (2006) Begomovirus evolution and disease management. Plant Virus Epidemiology. 297–316.
- 4. Rybicki EP (1994) A phylogenetic and evolutionary justification for 3 genera of geminiviridae. Archives of Virology 139: 49–77.
- 5. Xu XZ, Liu QP, Fan LJ, Cui XF, Zhou XP (2008) Analysis of synonymous codon usage and evolution of begomoviruses. Journal of Zhejiang University Science B 9: 667–674.
- 6. Rojas MR, Hagen C, Lucas WJ, Gilbertson RL (2005) Exploiting chinks in the plant's armor: evolution and emergence of geminiviruses. Annu Rev Phytopathol 43: 361–394.
- 7. Melgarejo TA, Kon T, Rojas MR, Paz-Carrasco L, Zerbini FM, et al. (2013) 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 87: 5397–5413.
- 8. Ha C, Coombs S, Revill P, Harding R, Vu M, et al. (2008) Molecular characterization of begomoviruses and DNA satellites from Vietnam: additional evidence that the New World geminiviruses were present in the Old World prior to continental separation. The Journal of general virology 89: 312–326.
- 9. Rothenstein D, Krenz B, Selchow O, Jeske H (2007) Tissue and cell tropism of Indian cassava mosaic virus (ICMV) and its AV2 (precoat) gene product. Virology 359: 137–145.
- 10. Poornima Priyadarshini CG, Ambika MV, Tippeswamy R, Savithri HS (2011) Functional characterization of coat protein and V2 involved in cell to cell movement of Cotton leaf curl Kokhran virus-Dabawali. PLoS ONE 6: e26929.
- 11. Glick E, Zrachya A, Levy Y, Mett A, Gidoni D, et al. (2008) Interaction with host SGS3 is required for suppression of RNA silencing by tomato yellow leaf curl virus V2 protein. Proc Natl Acad Sci U S A 105: 157–161.
- 12. Padidam M, Beachy RN, Fauquet CM (1995) Classification and identification of geminiviruses using sequence comparisons. Journal of General Virology 76: 249–263.
- 13. Briddon RW, Patil BL, Bagewadi B, Nawaz-ul-Rehman MS, Fauquet CM (2010) Distinct evolutionary histories of the DNA-A and DNA-B components of bipartite begomoviruses. BMC Evol Biol. 10: 97.
- 14. Modena NA, Zelada AM, Conte F, Mentaberry A (2008) Phosphorylation of the TGBp1 movement protein of Potato virus X by a Nicotiana tabacum CK2-like activity. Virus Res 137: 16–23.
- 15. Pruitt KD, Brown GR, Hiatt SM, Thibaud-Nissen F, Astashyn A, et al. (2014) RefSeq: an update on mammalian reference sequences. Nucleic Acids Res 42: D756–763.
- 16.
King AMQ, Adams MJ, Carstens EB, Lefkowitz EJ, editors (2012) Virus taxonomy: classification and nomenclature of viruses: Ninth Report of the International Committee on Taxonomy of Viruses. San Diego: Elsevier.
- 17. Notredame C, Higgins DG, Heringa J (2000) T-Coffee: A novel method for fast and accurate multiple sequence alignment. J Mol Biol 302: 205–217.
- 18. Suyama M, Torrents D, Bork P (2006) PAL2NAL: robust conversion of protein sequence alignments into the corresponding codon alignments. Nucleic Acids Res 34: W609–612.
- 19. Kosakovsky Pond SL, Frost SDW (2005) Datamonkey: rapid detection of selective pressure on individual sites of codon alignments. Bioinformatics 21: 2531–2533.
- 20. Delport W, Poon AF, Frost SD, Kosakovsky Pond SL (2010) Datamonkey 2010: a suite of phylogenetic analysis tools for evolutionary biology. Bioinformatics 26: 2455–2457.
- 21. Posada D, Buckley TR (2004) Model selection and model averaging in phylogenetics: advantages of akaike information criterion and bayesian approaches over likelihood ratio tests. Syst Biol 53: 793–808.
- 22. Rice P, Longden I, Bleasby A (2000) EMBOSS: the European Molecular Biology Open Software Suite. Trends Genet 16: 276–277.
- 23.
Team RDC (2012) R: a languate and environment for statistical computing. Vienna: R Foundation for Statistical Computing.
- 24. Cock PJ, Antao T, Chang JT, Chapman BA, Cox CJ, et al. (2009) Biopython: freely available Python tools for computational molecular biology and bioinformatics. Bioinformatics 25: 1422–1423.
- 25. Gilbertson RL, Sudarshana M, Jiang H, Rojas MR, Lucas WJ (2003) Limitations on geminivirus genome size imposed by plasmodesmata and virus-encoded movement protein: insights into DNA trafficking. Plant Cell 15: 2578–2591.
- 26. Rojas MR, Jiang H, Salati R, Xoconostle-Cazares B, Sudarshana MR, et al. (2001) Functional analysis of proteins involved in movement of the monopartite begomovirus, Tomato yellow leaf curl virus. Virology 291: 110–125.
- 27. Tonukari NJ, Scott-Craig JS, Walton JD (2000) The Cochliobolus carbonum SNF1 gene is required for cell wall-degrading enzyme expression and virulence on maize. Plant Cell 12: 237–248.
- 28. Duffy S, Holmes EC (2009) Validation of high rates of nucleotide substitution in geminiviruses: phylogenetic evidence from East African cassava mosaic viruses. J Gen Virol 90: 1539–1547.
- 29. Kosakovsky Pond SL, Frost SDW (2005) Not so different after all: a comparison of methods for detecting amino acid sites under selection. Mol Biol Evol 22: 1208–1222.
- 30. Duffy S, Shackelton LA, Holmes EC (2008) Rates of evolutionary change in viruses: patterns and determinants. Nat Rev Genet 9: 267–276.
- 31. Monjane AL, Pande D, Lakay F, Shepherd DN, van der Walt E, et al. (2012) Adaptive evolution by recombination is not associated with increased mutation rates in Maize streak virus. BMC Evol Biol 12: 252.
- 32. Rocha CS, Castillo-Urquiza GP, Lima AT, Silva FN, Xavier CA, et al. (2013) Brazilian begomovirus populations are highly recombinant, rapidly evolving, and segregated based on geographical location. J Virol 87: 5784–5799.
- 33. Kleinow T, Holeiter G, Nischang M, Stein M, Karayavuz M, et al. (2008) Post-translational modifications of Abutilon mosaic virus movement protein (BC1) in fission yeast. Virus Res 131: 86–94.
- 34. Crooks GE, Hon G, Chandonia JM, Brenner SE (2004) WebLogo: a sequence logo generator. Genome research 14: 1188–1190.
- 35. Sigrist CJ, de Castro E, Cerutti L, Cuche BA, Hulo N, et al. (2013) New and continuing developments at PROSITE. Nucleic Acids Res 41: D344–347.