Molecular Evolutionary Analysis of Potato Virus Y Infecting Potato Based on the VPg Gene

Molecular Footprints of Potato Virus Y Isolate Infecting Potatoes (Solanum tuberosum) in Kenya

Potato virus Y (PVY) is a highly diverse and genetically variable virus with various strains. Differential evolutionary routes have been reported in the genus Potyvirus, caused by natural selection pressure, mutation, and recombination, with their virulence being dependent on different environmental conditions. Despite its significance and economic impact on Solanaceous species, the understanding of PVY's phylogeography in Kenya remains limited and inadequately documented. The study centers on the molecular characterization of a Kenyan PVY isolate, GenBank accession number PP069009. In-depth phylogenetic analysis unveiled a strong evolutionary association between the Kenyan isolate and isolate [JQ924287] from the United States of America, supported by a robust 92% probability. Recombinant analyses exposed a mosaic-like genetic architecture within the Kenyan isolate, indicating multiple gene recombination events. Selection pressure scrutiny identified specific sites under selective pressure, with evidence of positive/diversifying and negative/purifying selection. Population genetics analysis revealed a calculated nucleotide diversity (π) of 0.00354881, while analysis of molecular variance (AMOVA) unveiled a structured genetic landscape with an øST value of 0.45224. The extensive haplotype network depicted the possibility of diverse PVY strains occurring across continents. This analysis provides valuable insights into the genetic diversity and distribution of PVY globally, highlighting the importance of understanding evolutionary dynamics for effective management and control strategies of PVY on a global scale.

Interaction of Solanum tuberosum L. translation initiation factors eIF4E with potato virus Y VPg: Apprehend and avoid

Potato virus Y (PVY) is one of the most dangerous agricultural pathogens that causes substantial harm to vegetative propagated crops, such as potatoes (Solanum tuberosum L.). A necessary condition for PVY infection is an interaction between the plant cap-binding translation initiation factors eIF4E and a viral protein VPg, which mimics the cap-structure. In this study, we identified the point mutations in potato eIF4E1 and eIF4E2 that disrupt VPg binding while preserving the functional activity. For the structural interpretation of the obtained results, molecular models of all the studied forms of eIF4E1 and eIF4E2 were constructed and analyzed via molecular dynamics. The results of molecular dynamics simulations corresponds to the biochemical results and suggests that the β1β2 loop plays a key role in the stabilization of both eIF4E-cap and eIF4E-VPg complexes.

Molecular evolution and phylogeographic analysis of wheat dwarf virus

Wheat dwarf virus (WDV) has caused considerable economic loss in the global production of grain crops. Knowledge of the evolutionary biology and population history of the pathogen remain poorly understood. We performed molecular evolution and worldwide phylodynamic analyses of the virus based on the genes in the protein-coding region of the entire viral genome. Our results showed that host-driven and geography-driven adaptation are major factors that affects the evolution of WDV. Bayesian phylogenetic analysis estimates that the average WDV substitution rate was 4.240 × 10-4 substitutions/site/year (95% credibility interval, 2.828 × 10-4-5.723 × 10-4), and the evolutionary rates of genes encoding proteins with virion-sense transcripts and genes encoding proteins with complementary-sense transcripts were different. The positively selected sites were detected in only two genes encoding proteins with complementary-sense, and WDV-barley are subject to stronger purifying selection than WDV-wheat. The time since the most recent common WDV ancestor was 1746 (95% credibility interval, 1517-1893) CE. Further analyses identified that the WDV-barley population and WDV-wheat population experienced dramatic expansion-decline episodes, and the expansion time of the WDV-barley population was earlier than that of the WDV-wheat population. Our phylogeographic analysis showed that the WDV population originating in Iran was subsequently introduced to Europe, and then spread from Eastern Europe to China.

Global molecular evolution and phylogeographic analysis of barley yellow dwarf virus based on the cp and mp genes

Barley yellow dwarf virus (BYDV) has caused considerable losses in the global production of grain crops such as wheat, barley and maize. We investigated the phylodynamics of the virus by analysing 379 and 485 nucleotide sequences of the genes encoding the coat protein and movement protein, respectively. The maximum clade credibility tree indicated that BYDV-GAV and BYDV-MAV, BYDV-PAV and BYDV-PAS share the same evolutionary lineage, respectively. The diversification of BYDV arises from its adaptability to vector insects and geography. Bayesian phylogenetic analyses showed that the mean substitution rates of the coat and movement proteins of BYDV ranged from 8.327 × 10^- 4 (4.700 × 10^- 4-1.228 × 10^- 3) and 8.671 × 10^- 4 (6.143 × 10^- 4-1.130 × 10^- 3) substitutions/site/year, respectively. The time since the most recent common BYDV ancestor was 1434 (1040-1766) CE (Common Era). The Bayesian skyline plot (BSP) showed that the BYDV population experienced dramatic expansions approximately 8 years into the 21st century, followed by a dramatic decline in less than 15 years. Our phylogeographic analysis showed that the BYDV population originating in the United States was subsequently introduced to Europe, South America, Australia and Asia. The migration pathways of BYDV suggest that the global spread of BYDV is associated with human activities.

Insight into aphid mediated Potato Virus Y transmission: A molecular to bioinformatics prospective

Potato, the world's most popular crop is reported to provide a food source for nearly a billion people. It is prone to a number of biotic stressors that affect yield and quality, out of which Potato Virus Y (PVY) occupies the top position. PVY can be transmitted mechanically and by sap-feeding aphid vectors. The application of insecticide causes an increase in the resistant vector population along with detrimental effects on the environment; genetic resistance and vector-virus control are the two core components for controlling the deadly PVY. Using transcriptomic tools together with differential gene expression and gene discovery, several loci and genes associated with PVY resistance have been widely identified. To combat this virus we must increase our understanding on the molecular response of the PVY-potato plant-aphid interaction and knowledge of genome organization, as well as the function of PVY encoded proteins, genetic diversity, the molecular aspects of PVY transmission by aphids, and transcriptome profiling of PVY infected potato cultivars. Techniques such as molecular and bioinformatics tools can identify and monitor virus transmission. Several studies have been conducted to understand the molecular basis of PVY resistance/susceptibility interactions and their impact on PVY epidemiology by studying the interrelationship between the virus, its vector, and the host plant. This review presents current knowledge of PVY transmission, epidemiology, genome organization, molecular to bioinformatics responses, and its effective management.

Bioinformatics in Plant Breeding and Research on Disease Resistance

In the context of plant breeding, bioinformatics can empower genetic and genomic selection to determine the optimal combination of genotypes that will produce a desired phenotype and help expedite the isolation of these new varieties. Bioinformatics is also instrumental in collecting and processing plant phenotypes, which facilitates plant breeding. Robots that use automated and digital technologies to collect and analyze different types of information to monitor the environment in which plants grow, analyze the environmental stresses they face, and promptly optimize suboptimal and adverse growth conditions accordingly, have helped plant research and saved human resources. In this paper, we describe the use of various bioinformatics databases and algorithms and explore their potential applications in plant breeding and for research on plant disease resistance.

Molecular Detection of Southern Tomato Amalgavirus Prevalent in Tomatoes and Its Genomic Characterization with Global Evolutionary Dynamics

Southern tomato amalgavirus (STV) is a cryptic pathogen that is abundant in tomato production fields and intensifies the resurgence of tomato yellow stunt disease (ToYSD), together with other phytoviruses. Here, we mapped the geographical and genomic diversity, phylogenetics, and evolutionary dynamics of STV. We found that STV prevailed across China and Pakistan, with a maximum average rate of infection of 43.19% in Beijing, China, and 40.08% in Punjab, Pakistan. Subsequently, we amplified, cloned, and annotated the complete genome sequences of STV isolates from Solanum lycopersicum L. in China (OP548653 and OP548652) and Pakistan (MT066231) using Sanger and next-generation sequencing (NGS). These STV isolates displayed close evolutionary relationships with others from Asia, America, and Europe. Whole-genome-based molecular diversity analysis showed that STV populations had 33 haplotypes with a gene diversity (Hd) of 0.977 and a nucleotide diversity (π) of 0.00404. The genetic variability of RNA-dependent RNA-polymerase (RdRp) was higher than that of the putative coat protein (CP) p42. Further analysis revealed that STV isolates were likely to be recombinant but with a lower-to-moderate level of confidence. With a variable distribution pattern of positively and negatively selected sites, negative selection pressure predominantly acted on p42 and RdRp. These findings elaborated on the molecular variability and evolutionary trends among STV populations across major tomato-producing regions of the world.

Characterization of yam mosaic viruses from Brazil reveals a new phylogenetic group and possible incursion from the African continent

Yam (Dioscorea spp.) is an important crop for smallholder farmers in the Northeast region of Brazil. Wherever yam is grown, diseases caused by yam mosaic virus (YMV) are prevalent. In the present study, the diversity of YMV infecting Dioscorea cayennensis-rotundata was analyzed. In addition, five species of Dioscorea (D. alata, D. altissima, D. bulbifera, D. subhastata, and D. trifida) commonly found in Brazil were analyzed using ELISA and high-throughput sequencing (HTS). YMV was detected only in D. cayennensis-rotundata, of which 66.7% of the samples tested positive in ELISA. Three YMV genome sequences were assembled from HTS and one by Sanger sequencing to group the sequences in a clade phylogenetically distinct from YMV from other origins. Temporal phylogenetic analyses estimated the mean evolutionary rate for the CP gene of YMV as 1.76 × 10^-3 substitutions per site per year, and the time to the most recent common ancestor as 168.68 years (95% Highest Posterior Density, HPD: 48.56-363.28 years), with a most likely geographic origin in the African continent. The data presented in this study contribute to reveal key aspects of the probable epidemiological history of YMV in Brazil.

Geography-Driven Evolution of Potato Virus A Revealed by Genetic Diversity Analysis of the Complete Genome

Potato virus A (PVA), a member of the genus Potyvirus, is an important potato pathogen that causes 30%-40% yield reduction to global potato production. Knowledge on the genetic structure and the evolutionary forces shaping the structure of this pathogen is limited but vital in developing effective management strategies. In this study, we investigated the population structure and molecular evolution of PVA by analyzing novel complete genomic sequences from Chinese isolates combined with available sequences from Europe, South America, Oceania, and North America. High nucleotide diversity was discovered among the populations studied. Pairwise F _ST values between geographical populations of PVA ranged from 0.22 to 0.46, indicating a significant spatial structure for this pathogen. Although purifying selection was detected at the majority of polymorphic sites, significant positive selection was identified in the P1, NIa, and NIb proteins, pointing to adaptive evolution of PVA. Further phylogeny-trait association analysis showed that the clustering of PVA isolates was significantly correlated with geographic regions, suggesting that geography-driven adaptation may be an important determinant of PVA diversification.

Viruses Without Borders: Global Analysis of the Population Structure, Haplotype Distribution, and Evolutionary Pattern of Iris Yellow Spot Orthotospovirus (Family Tospoviridae, Genus Orthotospovirus)

Iris yellow spot, caused by Iris yellow spot orthotospovirus (IYSV) (Genus: Orthotospovirus, Family: Tospoviridae), is an important disease of Allium spp. The complete N gene sequences of 142 IYSV isolates of curated sequence data from GenBank were used to determine the genetic diversity and evolutionary pattern. In silico restriction fragment length polymorphism (RFLP) analysis, codon-based maximum likelihood studies, genetic differentiation and gene flow within the populations of IYSV genotypes were investigated. Bayesian phylogenetic analysis was carried out to estimate the evolutionary rate. In silico RFLP analysis of N gene sequences categorized IYSV isolates into two major genotypes viz., IYSV Netherlands (IYSV_NL; 55.63%), IYSV Brazil (IYSV_BR; 38.73%) and the rest fell in neither group [IYSV other (IYSV_other; 5.63%)]. Phylogenetic tree largely corroborated the results of RFLP analysis and the IYSV genotypes clustered into IYSV_NL and IYSV_BR genotypes. Genetic diversity test revealed IYSV_other to be more diverse than IYSV_NL and IYSV_BR. IYSV_NL and IYSV_BR genotypes are under purifying selection and population expansion, whereas IYSV_other showed decreasing population size and hence appear to be under balancing selection. IYSV_BR is least differentiated from IYSV_other compared to IYSV_NL genotype based on nucleotide diversity. Three putative recombinant events were found in the N gene of IYSV isolates based on RDP analysis, however, RAT substantiated two among them. The marginal likelihood mean substitution rate was 5.08 × 10^–5 subs/site/year and 95% highest posterior density (HPD) substitution rate between 5.11 × 10^–5 and 5.06 × 10^–5. Findings suggest that IYSV continues to evolve using population expansion strategies. The substitution rates identified are similar to other plant RNA viruses.

Genetic Variability and Molecular Evolution of Maize Yellow Mosaic Virus Populations from Different Geographic Origins

Maize yellow mosaic virus (MaYMV) hosted in various gramineous plants was assigned to the genus Polerovirus (family Luteoviridae) in 2018. However, little is known about its genetic diversity and population structure. In this study, 509 sugarcane leaf samples with mosaic symptoms were collected in 2017 to 2019 from eight sugarcane-growing provinces in China. Reverse-transcription PCR results revealed that four positive-sense RNA viruses were found to infect sugarcane, and the incidence of MaYMV among samples from Fujian, Sichuan, and Guangxi Provinces was 52.1, 9.8, and 2.5%, respectively. Based on 82 partial MaYMV sequences and 46 whole-genome sequences from different host plants, phylogenetic analysis revealed that MaYMV populations are very closely associated with their source geographical regions (China, Africa, and South America). Pairwise identity analysis showed significant variability in genome sequences among MaYMV isolates with genomic nucleotide identities of 91.1 to 99.9%. In addition to codon mutations, insertions or deletions also contributed to genetic variability in individual coding regions, especially in the readthrough protein (P3-P5 fusion protein). Low gene flow and significant genetic differentiation of MaYMV were observed among the three geographical populations, suggesting that environmental adaptation is an important evolutionary force that shapes the genetic structure of MaYMV. Genes in the MaYMV genome were subject to strong negative or purification selection during evolution, except for the movement protein (MP), which was under positive selection pressure. This finding suggests that the MP may play an important role in MaYMV evolution. Taken together, our findings provide basic information for the development of an integrated disease management strategy against MaYMV.

Phylodynamics and Codon Usage Pattern Analysis of Broad Bean Wilt Virus 2

Broad bean wilt virus 2 (BBWV-2), which belongs to the genus Fabavirus of the family Secoviridae, is an important pathogen that causes damage to broad bean, pepper, yam, spinach and other economically important ornamental and horticultural crops worldwide. Previously, only limited reports have shown the genetic variation of BBWV2. Meanwhile, the detailed evolutionary changes, synonymous codon usage bias and host adaptation of this virus are largely unclear. Here, we performed comprehensive analyses of the phylodynamics, reassortment, composition bias and codon usage pattern of BBWV2 using forty-two complete genome sequences of BBWV-2 isolates together with two other full-length RNA1 sequences and six full-length RNA2 sequences. Both recombination and reassortment had a significant influence on the genomic evolution of BBWV2. Through phylogenetic analysis we detected three and four lineages based on the ORF1 and ORF2 nonrecombinant sequences, respectively. The evolutionary rates of the two BBWV2 ORF coding sequences were 8.895 × 10⁻⁴ and 4.560 × 10⁻⁴ subs/site/year, respectively. We found a relatively conserved and stable genomic composition with a lower codon usage choice in the two BBWV2 protein coding sequences. ENC-plot and neutrality plot analyses showed that natural selection is the key factor shaping the codon usage pattern of BBWV2. Strong correlations between BBWV2 and broad bean and pepper were observed from similarity index (SiD), codon adaptation index (CAI) and relative codon deoptimization index (RCDI) analyses. Our study is the first to evaluate the phylodynamics, codon usage patterns and adaptive evolution of a fabavirus, and our results may be useful for the understanding of the origin of this virus.

Analysis of selection in protein-coding sequences accounting for common biases

The evolution of protein-coding genes is usually driven by selective processes, which favor some evolutionary trajectories over others, optimizing the subsequent protein stability and activity. The analysis of selection in this type of genetic data is broadly performed with the metric nonsynonymous/synonymous substitution rate ratio (dN/dS). However, most of the well-established methodologies to estimate this metric make crucial assumptions, such as lack of recombination or invariable codon frequencies along genes, which can bias the estimation. Here, we review the most relevant biases in the dN/dS estimation and provide a detailed guide to estimate this metric using state-of-the-art procedures that account for such biases, along with illustrative practical examples and recommendations. We also discuss the traditional interpretation of the estimated dN/dS emphasizing the importance of considering complementary biological information such as the role of the observed substitutions on the stability and function of proteins. This review is oriented to help evolutionary biologists that aim to accurately estimate selection in protein-coding sequences.

The evolutionary history and global spatio-temporal dynamics of potato virus Y

Potato virus Y (PVY) is a destructive plant pathogen that causes considerable losses to global potato and tobacco production. Although the molecular structure of PVY is well characterized, the evolutionary and global transmission dynamics of this virus remain poorly understood. We investigated the phylodynamics of the virus by analysing 253 nucleotide sequences of the genes encoding the third protein (P3), cylindrical inclusion protein (CI), and the nuclear inclusion protein (NIb). Our Bayesian phylogenetic analyses showed that the mean substitution rates of different regions of the genome ranged from 8.50 × 10^-5 to 1.34 × 10^-4 substitutions/site/year, whereas the time to the most recent common ancestor of PVY varied with the length of the genomic regions and with the number of viral isolates being analysed. Our phylogeographic analysis showed that the PVY population originated in South America and was introduced into Europe in the 19th century, from where it spread around the globe. The migration pathways of PVY correlate well with the trade routes of potato tubers, suggesting that the global spread of PVY is associated with human activities.

Genetic changes and host adaptability in sugarcane mosaic virus based on complete genome sequences

Sugarcane mosaic virus (SCMV), a member of the genus Potyvirus in the family Potyviridae, is an important pathogen that causes mosaic diseases in maize, sugarcane, canna and other graminaceous species worldwide. Previously, several reports have showed the genetic variation and population structure of SCMV. However, the evolutionary dynamics, synonymous codon usage pattern and adaptive evolution of the virus is unclear. In this study, we performed comprehensive analyses of phylodynamics, composition bias and codon usage of SCMV using 108 complete genomic sequences. Our phylogenetic analysis found six host- and geographically confined phylogenetic lineages within the SCMV non-recombinant isolates. We found a relatively stable and conserved genomic composition with a lower codon usage choice in the SCMV protein coding sequences. Mutation pressure and natural selection have shaped the codon usage patterns of the SCMV protein coding sequences with natural selection being the dominant factor. The codon adaptation index (CAI), relative codon deoptimization index (RCDI) and similarity index (SiD) analyses revealed a stronger correlation between SCMV and maize than between SCMV and sugarcane or canna. Our study is the first to evaluate the codon usage pattern of SCMV based on complete sequences and may provide a better understanding of the origin of SCMV and its evolutionary patterns for future research.

The Potyviruses: An Evolutionary Synthesis Is Emerging

In this review, encouraged by the dictum of Theodosius Dobzhansky that "Nothing in biology makes sense except in the light of evolution", we outline the likely evolutionary pathways that have resulted in the observed similarities and differences of the extant molecules, biology, distribution, etc. of the potyvirids and, especially, its largest genus, the potyviruses. The potyvirids are a family of plant-infecting RNA-genome viruses. They had a single polyphyletic origin, and all share at least three of their genes (i.e., the helicase region of their CI protein, the RdRp region of their NIb protein and their coat protein) with other viruses which are otherwise unrelated. Potyvirids fall into 11 genera of which the potyviruses, the largest, include more than 150 distinct viruses found worldwide. The first potyvirus probably originated 15,000-30,000 years ago, in a Eurasian grass host, by acquiring crucial changes to its coat protein and HC-Pro protein, which enabled it to be transmitted by migrating host-seeking aphids. All potyviruses are aphid-borne and, in nature, infect discreet sets of monocotyledonous or eudicotyledonous angiosperms. All potyvirus genomes are under negative selection; the HC-Pro, CP, Nia, and NIb genes are most strongly selected, and the PIPO gene least, but there are overriding virus specific differences; for example, all turnip mosaic virus genes are more strongly conserved than those of potato virus Y. Estimates of dN/dS (ω) indicate whether potyvirus populations have been evolving as one or more subpopulations and could be used to help define species boundaries. Recombinants are common in many potyvirus populations (20%-64% in five examined), but recombination seems to be an uncommon speciation mechanism as, of 149 distinct potyviruses, only two were clear recombinants. Human activities, especially trade and farming, have fostered and spread both potyviruses and their aphid vectors throughout the world, especially over the past five centuries. The world distribution of potyviruses, especially those found on islands, indicates that potyviruses may be more frequently or effectively transmitted by seed than experimental tests suggest. Only two meta-genomic potyviruses have been recorded from animal samples, and both are probably contaminants.