The effect of ecosystem biodiversity on virus genetic diversity depends on virus species: A study of chiltepin-infecting begomoviruses in Mexico

Begomovirus populations in single plants are complex and may include both well-adapted and poorly-adapted viruses

Begomoviruses (single-stranded DNA plant viruses transmitted by whiteflies) are economically important pathogens causing epidemics worldwide. Tomato-infecting begomoviruses emerged in Brazil in the 1990's following the introduction of Bemisia tabaci Middle East-Asia Minor 1. It is believed that these viruses evolved from indigenous viruses infecting non-cultivated hosts. However, tomato-infecting viruses are rarely found in non-cultivated hosts, and vice-versa. It is possible that viral populations in a given host are composed primarily of viruses which are well adapted to this host, but also include a small proportion of poorly adapted viruses. Following transfer to a new host, the composition of the viral population would shift rapidly, with the viruses which are better adapted to the new host becoming predominant. To test this hypothesis, we collected tomato and Sida plants growing next to each other at two locations in 2014 and 2018. Total DNA was extracted from tomato and Sida samples from each location and year and used as a template for high-throughput sequencing. Reads were mapped following a highly stringent set of criteria. For the 2014 samples, >98% of the Sida reads mapped to Sida micrantha mosaic virus (SiMMV), but 0.1% of the reads mapped to tomato severe rugose virus (ToSRV). Conversely, >99% of the tomato reads mapped to ToSRV, with 0.18% mapping to SiMMV. For the 2018 samples, 41% of the Sida reads mapped to three Sida-adapted viruses and 0.1% of the reads mapped to ToSRV, while 99.9% of the tomato reads mapped to ToSRV. These results are consistent with the hypothesis that viral populations in a single plant are composed primarily of the virus that is better adapted to the host but also include a small proportion of viruses that are poorly adapted.

Determinants of Virus Variation, Evolution, and Host Adaptation

Virus evolution is the change in the genetic structure of a viral population over time and results in the emergence of new viral variants, strains, and species with novel biological properties, including adaptation to new hosts. There are host, vector, environmental, and viral factors that contribute to virus evolution. To achieve or fine tune compatibility and successfully establish infection, viruses adapt to a particular host species or to a group of species. However, some viruses are better able to adapt to diverse hosts, vectors, and environments. Viruses generate genetic diversity through mutation, reassortment, and recombination. Plant viruses are exposed to genetic drift and selection pressures by host and vector factors, and random variants or those with a competitive advantage are fixed in the population and mediate the emergence of new viral strains or species with novel biological properties. This process creates a footprint in the virus genome evident as the preferential accumulation of substitutions, insertions, or deletions in areas of the genome that function as determinants of host adaptation. Here, with respect to plant viruses, we review the current understanding of the sources of variation, the effect of selection, and its role in virus evolution and host adaptation.

Gene Overlapping as a Modulator of Begomovirus Evolution

In RNA viruses, which have high mutation—and fast evolutionary— rates, gene overlapping (i.e., genomic regions that encode more than one protein) is a major factor controlling mutational load and therefore the virus evolvability. Although DNA viruses use host high-fidelity polymerases for their replication, and therefore should have lower mutation rates, it has been shown that some of them have evolutionary rates comparable to those of RNA viruses. Notably, these viruses have large proportions of their genes with at least one overlapping instance. Hence, gene overlapping could be a modulator of virus evolution beyond the RNA world. To test this hypothesis, we use the genus Begomovirus of plant viruses as a model. Through comparative genomic approaches, we show that terminal gene overlapping decreases the rate of virus evolution, which is associated with lower frequency of both synonymous and nonsynonymous mutations. In contrast, terminal overlapping has little effect on the pace of virus evolution. Overall, our analyses support a role for gene overlapping in the evolution of begomoviruses and provide novel information on the factors that shape their genetic diversity.

Evolutionary dynamics of bipartite begomoviruses revealed by complete genome analysis

Several key evolutionary events marked the evolution of geminiviruses, culminating with the emergence of divided (bipartite) genomes represented by viruses classified in the genus Begomovirus. This genus represents the most abundant group of multipartite viruses, contributing significantly to the observed abundance of multipartite species in the virosphere. Although aspects related to virus-host interactions and evolutionary dynamics have been extensively studied, the bipartite nature of these viruses has been little explored in evolutionary studies. Here, we performed a parallel evolutionary analysis of the DNA-A and DNA-B segments of New World begomoviruses. A total of 239 full-length DNA-B sequences obtained in this study, combined with 292 DNA-A and 76 DNA-B sequences retrieved from GenBank, were analysed. The results indicate that the DNA-A and DNA-B respond differentially to evolutionary processes, with the DNA-B being more permissive to variation and more prone to recombination than the DNA-A. Although a clear geographic segregation was observed for both segments, differences in the genetic structure between DNA-A and DNA-B were also observed, with cognate segments belonging to distinct genetic clusters. DNA-B coding regions evolve under the same selection pressures than DNA-A coding regions. Together, our results indicate an interplay between reassortment and recombination acting at different levels across distinct subpopulations and segments.

Integrative Analysis of the Metabolome and Transcriptome of a Cultivated Pepper and Its Wild Progenitor Chiltepin (Capsicum annuum L. var. glabriusculum) Revealed the Loss of Pungency During Capsicum Domestication

Pungency is a unique characteristic of chili peppers (Capsicum spp.) caused by capsaicinoids. The evolutionary emergence of pungency is thought to be a derived trait within the genus Capsicum. However, it is not well-known how pungency has varied during Capsicum domestication and specialization. In this study, we applied a comparative metabolomics along with transcriptomics analysis to assess various changes between two peppers (a mildly pungent cultivated pepper BB3 and its hot progenitor chiltepin) at four stages of fruit development, focusing on pungency variation. A total of 558 metabolites were detected in two peppers. In comparison with chiltepin, capsaicinoid accumulation in BB3 was almost negligible at the early stage. Next, 412 DEGs associated with the capsaicinoid accumulation pathway were identified through coexpression analysis, of which 18 genes (14 TFs, 3 CBGs, and 1 UGT) were deemed key regulators due to their high coefficients. Based on these data, we speculated that downregulation of these hub genes during the early fruit developmental stage leads to a loss in pungency during Capsicum domestication (from chiltepin to BB3). Of note, a putative UDP-glycosyltransferase, GT86A1, is thought to affect the stabilization of capsaicinoids. Our results lay the foundation for further research on the genetic diversity of pungency traits during Capsicum domestication and specialization.

Diversity of tomato-infecting begomoviruses and spatiotemporal dynamics of an endemic viral species of the Brazilian Atlantic rain forest biome

Yield losses induced by a complex of begomoviruses are observed across all major tomato-producing areas in Brazil. Tomato severe rugose virus (ToSRV) is the most widespread begomovirus in the country. Conversely, tomato common mosaic virus (ToCmMV) displays a more restricted geographical distribution to areas associated with the Atlantic Rain Forest (ARF) biome, encompassing the States of Espírito Santo-ES, Minas Gerais-MG, and Rio de Janeiro-RJ. Here, we characterized 277 tomato-infecting isolates collected in fields located within the ARF biome from 2006 to 2018. ToSRV displayed the highest prevalence (n = 157), followed by ToCmMV (n = 95) and tomato interveinal chlorosis virus (n = 14). Four other begomoviruses were also detected, but with very low incidences. ToCmMV was the predominant begomovirus in the ARF biome up to 2014-2015 with very low ToSRV incidence. Subsequently, ToSRV became the most prevalent species in ES and RJ, but ToCmMV was still predominating in the "Zona da Mata" meso-region in MG. Due to the remarkable endemic distribution of ToCmMV, we carried out phylogeographical studies of this virus using information from all 28 available isolates with complete DNA-A sequences. The closest common ancestor of ToCmMV was more likely originated around Coimbra-MG area ≈ 25 years before the formal report of this viral species. So far, all surveys indicated tomatoes as the only natural hosts of ToCmMV with outbreaks occurring mainly (but not exclusively) in highland areas. ToSRV shows a more widespread incidence across both highland and lowland areas of the ARF biome.

Intra-host evolution of the ssDNA virus tomato severe rugose virus (ToSRV)

To evaluate and quantify the evolutionary dynamics of the bipartite begomovirus tomato severe rugose virus (ToSRV) in a cultivated and a non-cultivated host, plants of tomato and Nicandra physaloides were biolistically inoculated with an infectious clone and systemically infected leaves were sampled at 30, 75 and 120 days after inoculation. Total DNA was extracted and sequenced in the Illumina HiSeq 2000 platform. The datasets were trimmed with the quality score limit set to 0.01, and the assembly was performed using the infectious clone sequence as reference. SNPs were filtered using a minimum p-value of 0.001 and the sum frequencies were used to calculate the deviation from the original clone sequence. Nucleotide substitution rates were calculated for the two DNA components in both hosts: 1.73 × 10^-3 and 3.07 × 10^-4 sub/site/year for the DNA-A and DNA-B, respectively, in N. physaloides, and 8.05 × 10^-4 and 7.02 × 10^-5 sub/site/year the for DNA-A and DNA-B, respectively, in tomato. These values are in the same range of those estimated for viruses with single-stranded RNA genomes and for other begomoviruses. Strikingly, the number of substitutions decreased over time, suggesting the presence of bottlenecks during systemic infection. Determination of Shannon's entropy indicated different patterns of variation in the DNA-A and the DNA-B, suggesting distinct evolutionary forces acting upon each component.

Host Abundance and Identity Determine the Epidemiology and Evolution of a Generalist Plant Virus in a Wild Ecosystem

Increasing evidence indicates that in wild ecosystems plant viruses are important ecological agents, and with potential to jump into crops, but only recently have the diversity and population dynamics of wild plant viruses begun to be explored. Theory proposes that biotic factors (e.g., ecosystem biodiversity, host abundance, and host density) and climatic conditions would determine the epidemiology and evolution of wild plant viruses. However, these predictions seldom have been empirically tested. For 3 years, we analyzed the prevalence and genetic diversity of Potyvirus species in preserved riparian forests of Spain. Results indicated that potyviruses were always present in riparian forests, with a novel generalist potyvirus species provisionally named Iberian hop mosaic virus (IbHMV), explaining the largest fraction of infected plants. Focusing on this potyvirus, we analyzed the biotic and climatic factors affecting virus infection risk and population genetic diversity in its native ecosystem. The main predictors of IbHMV infection risk were host relative abundance and species richness. Virus prevalence and host relative abundance were the major factors determining the genetic diversity and selection pressures in the virus population. These observations support theoretical predictions assigning these ecological factors a key role in parasite epidemiology and evolution. Finally, our phylogenetic analysis indicated that the viral population was genetically structured according to host and location of origin, as expected if speciation is largely sympatric. Thus, this work contributes to characterizing viral diversity and provides novel information on the determinants of plant virus epidemiology and evolution in wild ecosystems.

Role of host genetic diversity for susceptibility-to-infection in the evolution of virulence of a plant virus

Predicting viral emergence is difficult due to the stochastic nature of the underlying processes and the many factors that govern pathogen evolution. Environmental factors affecting the host, the pathogen and the interaction between both are key in emergence. In particular, infectious disease dynamics are affected by spatiotemporal heterogeneity in their environments. A broad knowledge of these factors will allow better estimating where and when viral emergence is more likely to occur. Here, we investigate how the population structure for susceptibility-to-infection genes of the plant Arabidopsis thaliana shapes the evolution of Turnip mosaic virus (TuMV). For doing so we have evolved TuMV lineages in two radically different host population structures: (1) a metapopulation subdivided into six demes (subpopulations); each one being composed of individuals from only one of six possible A. thaliana ecotypes and (2) a well-mixed population constituted by equal number of plants from the same six A. thaliana ecotypes. These two populations were evolved for twelve serial passages. At the end of the experimental evolution, we found faster adaptation of TuMV to each ecotype in the metapopulation than in the well-mixed heterogeneous host populations. However, viruses evolved in well-mixed populations were more pathogenic and infectious than viruses evolved in the metapopulation. Furthermore, the viruses evolved in the demes showed stronger signatures of local specialization than viruses evolved in the well-mixed populations. These results illustrate how the genetic diversity of hosts in an experimental ecosystem favors the evolution of virulence of a pathogen.

High-Throughput Sequencing Reveals Differential Begomovirus Species Diversity in Non-Cultivated Plants in Northern-Pacific Mexico

Plant DNA viruses of the genus Begomovirus have been documented as the most genetically diverse in the family Geminiviridae and present a serious threat for global horticultural production, especially considering climate change. It is important to characterize naturally existing begomoviruses, since viral genetic diversity in non-cultivated plants could lead to future disease epidemics in crops. In this study, high-throughput sequencing (HTS) was employed to determine viral diversity of samples collected in a survey performed during 2012–2016 in seven states of Northern-Pacific Mexico, areas of diverse climatic conditions where different vegetable crops are subject to intensive farming. In total, 132 plant species, belonging to 34 families, were identified and sampled in the natural ecosystems surrounding cultivated areas (agro-ecological interface). HTS analysis and subsequent de novo assembly revealed a number of geminivirus-related DNA signatures with 80 to 100% DNA similarity with begomoviral sequences present in the genome databank. The analysis revealed DNA signatures corresponding to 52 crop-infecting and 35 non-cultivated-infecting geminiviruses that, interestingly, were present in different plant species. Such an analysis deepens our knowledge of geminiviral diversity and could help detecting emerging viruses affecting crops in different agro-climatic regions.

Life on the Edge: Geminiviruses at the Interface Between Crops and Wild Plant Hosts

Viruses constitute the largest group of emerging pathogens, and geminiviruses (plant viruses with circular, single-stranded DNA genomes) are the major group of emerging plant viruses. With their high potential for genetic variation due to mutation and recombination, their efficient spread by vectors, and their wide host range as a group, including both wild and cultivated hosts, geminiviruses are attractive models for the study of the evolutionary and ecological factors driving virus emergence. Studies on the epidemiological features of geminivirus diseases have traditionally focused primarily on crop plants. Nevertheless, knowledge of geminivirus infection in wild plants, and especially at the interface between wild and cultivated plants, is necessary to provide a complete view of their ecology, evolution, and emergence. In this review, we address the most relevant aspects of geminivirus variability and evolution in wild and crop plants and geminiviruses' potential to emerge in crops.

Natural Hosts and Genetic Diversity of the Emerging Tomato Leaf Curl New Delhi Virus in Spain

Knowledge about the host range and genetic structure of emerging plant viruses provides insights into fundamental ecological and evolutionary processes, and from an applied perspective, facilitates the design and implementation of sustainable disease control measures. Tomato leaf curl New Delhi virus (ToLCNDV) is an emerging whitefly transmitted begomovirus that is rapidly spreading and inciting economically important diseases in cucurbit crops of the Mediterranean basin. Genetic characterization of the ToLCNDV Mediterranean populations has shown that they are monophyletic in cucurbit plants. However, the extent to which other alternative (cultivated and wild) hosts may affect ToLCNDV genetic population structure and virus prevalence remains unknown. In this study a total of 683 samples from 13 cultivated species, and 203 samples from 24 wild species from three major cucurbit-producing areas of Spain (Murcia, Alicante and Castilla-La Mancha) from five cropping seasons (2012-2016) were analyzed for ToLCNDV infection. Except for watermelon, ToLCNDV was detected in all cultivated-cucurbit species as well as in tomato. Among weeds, Ecballium elaterium, Datura stramonium, Sonchus oleraceus, and Solanum nigrum were identified as alternative ToLCNDV plant hosts, which could act as new potential sources of virus inoculum. Furthermore, we performed full-genome deep-sequencing of 80 ToLCNDV isolates from different hosts, location and cropping year. Our phylogenetic analysis supports a Mediterranean virus population that is genetically very homogeneous, with no clustering pattern, and clearly different from Asian virus populations. Additionally, D. stramonium displayed higher levels of within-host genetic diversity than cultivated plants, and this variability appeared to increase with time. These results suggest that the potential ToLCNDV adaptive evolution occurring in wild plant hosts could serve as a source of virus genetic variability, thereby affecting the genetic structure and spatial-temporal dynamics of the viral population.

Addressing Research Needs in the Field of Plant Virus Ecology by Defining Knowledge Gaps and Developing Wild Dicot Study Systems

Viruses are ubiquitous within all habitats that support cellular life and represent the most important emerging infectious diseases of plants. Despite this, it is only recently that we have begun to describe the ecological roles of plant viruses in unmanaged systems and the influence of ecosystem properties on virus evolution. We now know that wild plants frequently harbor infections by diverse virus species, but much remains to be learned about how viruses influence host traits and how hosts influence virus evolution and vector interactions. To identify knowledge gaps and suggest avenues for alleviating research deficits, we performed a quantitative synthesis of a representative sample of virus ecology literature, developed criteria for expanding the suite of pathosystems serving as models, and applied these criteria through a case study. We found significant gaps in the types of ecological systems studied, which merit more attention. In particular, there is a strong need for a greater diversity of logistically tractable, wild dicot perennial study systems suitable for experimental manipulations of infection status. Based on criteria developed from our quantitative synthesis, we evaluated three California native dicot perennials typically found in Mediterranean-climate plant communities as candidate models: Cucurbita foetidissima (buffalo gourd), Cucurbita palmata (coyote gourd), and Datura wrightii (sacred thorn-apple). We used Illumina sequencing and network analyses to characterize viromes and viral links among species, using samples taken from multiple individuals at two different reserves. We also compared our Illumina workflow with targeted RT-PCR detection assays of varying costs. To make this process accessible to ecologists looking to incorporate virology into existing studies, we describe our approach in detail and discuss advantages and challenges of different protocols. We also provide a bioinformatics workflow based on open-access tools with graphical user interfaces. Our study provides evidence that dicot perennials in xeric habitats support multiple, asymptomatic infections by viruses known to be pathogenic in related crop hosts. Quantifying the impacts of these interactions on plant performance and virus epidemiology in our logistically tractable host systems will provide fundamental information about plant virus ecology outside of crop environments.

The diversity, evolution and epidemiology of plant viruses: A phylogenetic view

During the past four decades, the scientific community has seen an exponential advance in the number, sophistication, and quality of molecular techniques and bioinformatics tools for the genetic characterization of plant virus populations. Predating these advances, the field of Phylogenetics has significantly contributed to understand important aspects of plant virus evolution. This review aims at summarizing the impact of Phylogenetics in the current knowledge on three major aspects of plant virus evolution that have benefited from the development of phylogenetic inference: (1) The identification and classification of plant virus diversity. (2) The mechanisms and forces shaping the evolution of plant virus populations. (3) The understanding of the interaction between plant virus evolution, epidemiology and ecology. The work discussed here highlights the important role of phylogenetic approaches in the study of the dynamics of plant virus populations.

The impact of host genetic diversity on virus evolution and emergence

Accumulating evidence indicates that biodiversity has an important impact on parasite evolution and emergence. The vast majority of studies in this area have only considered the diversity of species within an environment as an overall measure of biodiversity, overlooking the role of genetic diversity within a particular host species. Although theoretical models propose that host genetic diversity in part shapes that of the infecting parasite population, and hence modulates the risk of parasite emergence, this effect has seldom been tested empirically. Using Rabies virus (RABV) as a model parasite, we provide evidence that greater host genetic diversity increases both parasite genetic diversity and the likelihood of a host being a donor in RABV cross-species transmission events. We conclude that host genetic diversity may be an important determinant of parasite evolution and emergence.

Ecological Factors Affecting Infection Risk and Population Genetic Diversity of a Novel Potyvirus in Its Native Wild Ecosystem

Increasing evidence indicates that there is ample diversity of plant virus species in wild ecosystems. The vast majority of this diversity, however, remains uncharacterized. Moreover, in these ecosystems the factors affecting plant virus infection risk and population genetic diversity, two traits intrinsically linked to virus emergence, are largely unknown. Along 3 years, we have analyzed the prevalence and diversity of plant virus species from the genus Potyvirus in evergreen oak forests of the Iberian Peninsula, the main wild ecosystem in this geographic region and in the entire Mediterranean basin. During this period, we have also measured plant species diversity, host density, plant biomass, temperature, relative humidity, and rainfall. Results indicated that potyviruses were always present in evergreen oak forests, with a novel virus species explaining the largest fraction of potyvirus-infected plants. We determined the genomic sequence of this novel virus and we explored its host range in natural and greenhouse conditions. Natural host range was limited to the perennial plant mountain rue (Ruta montana), commonly found in evergreen oak forests of the Iberian Peninsula. In this host, the virus was highly prevalent and was therefore provisionally named mediterranean ruda virus (MeRV). Focusing in this natural host-virus interaction, we analyzed the ecological factors affecting MeRV infection risk and population genetic diversity in its native wild ecosystem. The main predictor of virus infection risk was the host density. MeRV prevalence was the major factor determining genetic diversity and selection pressures in the virus populations. This observation supports theoretical predictions assigning these two traits a key role in parasite epidemiology and evolution. Thus, our analyses contribute both to characterize viral diversity and to understand the ecological determinants of virus population dynamics in wild ecosystems.

Environmental heterogeneity and the evolution of plant-virus interactions: Viruses in wild pepper populations

Understanding host-pathogen interactions requires analyses to address the multiplicity of scales in heterogeneous landscapes. Anthropogenic influence on plant communities, especially cultivation, is a major cause of environmental heterogeneity. We have approached the analysis of how environmental heterogeneity determines plant-virus interactions by studying virus infection in a wild plant currently undergoing incipient domestication, the wild pepper or chiltepin, across its geographical range in Mexico. We have shown previously that anthropogenic disturbance is associated with higher infection and disease risk, and with disrupted patterns of host and virus genetic spatial structure. We now show that anthropogenic factors, species richness, host genetic diversity and density in communities supporting chiltepin differentially affect infection risk according to the virus analysed. We also show that in addition to these factors, a broad range of abiotic and biotic variables meaningful to continental scales, have an important role on the risk of infection depending on the virus. Last, we show that natural virus infection of chiltepin plants in wild communities results in decreased survival and fecundity, hence negatively affecting fitness. This important finding paves the way for future studies on plant-virus co-evolution.