Modelling the evolutionary dynamics of viruses within their hosts: a case study using high-throughput sequencing

Plant eIF4E isoforms as factors of susceptibility and resistance to potyviruses

Potyviruses are the largest group of plant-infecting RNA viruses that affect a wide range of crop plants. Plant resistance genes against potyviruses are often recessive and encode translation initiation factors eIF4E. The inability of potyviruses to use plant eIF4E factors leads to the development of resistance through a loss-of-susceptibility mechanism. Plants have a small family of eIF4E genes that encode several isoforms with distinct but overlapping functions in cell metabolism. Potyviruses use distinct eIF4E isoforms as susceptibility factors in different plants. The role of different members of the plant eIF4E family in the interaction with a given potyvirus could differ drastically. An interplay exists between different members of the eIF4E family in the context of plant-potyvirus interactions, allowing different eIF4E isoforms to modulate each other's availability as susceptibility factors for the virus. In this review, possible molecular mechanisms underlying this interaction are discussed, and approaches to identify the eIF4E isoform that plays a major role in the plant-potyvirus interaction are suggested. The final section of the review discusses how knowledge about the interaction between different eIF4E isoforms can be used to develop plants with durable resistance to potyviruses.

Reversion of a resistance-breaking mutation shows reversion costs and high virus diversity at necrotic local lesions

An instance of host range evolution relevant to plant virus disease control is resistance breaking. Resistance breaking can be hindered by across-host fitness trade-offs generated by negative effects of resistance-breaking mutations on the virus fitness in susceptible hosts. Different mutations in pepper mild mottle virus (PMMoV) coat protein result in the breaking in pepper plants of the resistance determined by the L³ resistance allele. Of these, mutation M138N is widespread in PMMoV populations, despite associated fitness penalties in within-host multiplication and survival. The stability of mutation M138N was analysed by serial passaging in L³ resistant plants. Appearance on passaging of necrotic local lesions (NLL), indicating an effective L³ resistance, showed reversion to nonresistance-breaking phenotypes was common. Most revertant genotypes had the mutation N138K, which affects the properties of the virus particle, introducing a penalty of reversion. Hence, the costs of reversion may determine the evolution of resistance-breaking in addition to resistance-breaking costs. The genetic diversity of the virus population in NLL was much higher than in systemically infected tissues, and included mutations reported to break L³ resistance other than M138N. Infectivity assays on pepper genotypes with different L alleles showed high phenotypic diversity in respect to L alleles in NLL, including phenotypes not reported in nature. Thus, high diversity at NLL may potentiate the appearance of genotypes that enable the colonization of new host genotypes or species. Collectively, the results of this study contribute to better understanding the evolutionary dynamics of resistance breaking and host-range expansions.

Plant Virus Adaptation to New Hosts: A Multi-scale Approach

Viruses are studied at each level of biological complexity: from within-cells to ecosystems. The same basic evolutionary forces and principles operate at each level: mutation and recombination, selection, genetic drift, migration, and adaptive trade-offs. Great efforts have been put into understanding each level in great detail, hoping to predict the dynamics of viral population, prevent virus emergence, and manage their spread and virulence. Unfortunately, we are still far from this. To achieve these ambitious goals, we advocate for an integrative perspective of virus evolution. Focusing in plant viruses, we illustrate the pervasiveness of the above-mentioned principles. Beginning at the within-cell level, we describe replication modes, infection bottlenecks, and cellular contagion rates. Next, we move up to the colonization of distal tissues, discussing the fundamental role of random events. Then, we jump beyond the individual host and discuss the link between transmission mode and virulence. Finally, at the community level, we discuss properties of virus-plant infection networks. To close this review we propose the multilayer network theory, in which elements at different layers are connected and submit to their own dynamics that feed across layers, resulting in new emerging properties, as a way to integrate information from the different levels.

Phylodynamics of a regional SARS-CoV-2 rapid spreading event in Colorado in late 2020

Since the initial reported discovery of SARS-CoV-2 in late 2019, genomic surveillance has been an important tool to understand its transmission and evolution. Here, we sought to describe the underlying regional phylodynamics before and during a rapid spreading event that was documented by surveillance protocols of the United States Air Force Academy (USAFA) in late October-November of 2020. We used replicate long-read sequencing on Colorado SARS-CoV-2 genomes collected July through November 2020 at the University of Colorado Anschutz Medical campus in Aurora and the United States Air Force Academy in Colorado Springs. Replicate sequencing allowed rigorous validation of variation and placement in a phylogenetic relatedness network. We focus on describing the phylodynamics of a lineage that likely originated in the local Colorado Springs community and expanded rapidly over the course of two months in an outbreak within the well-controlled environment of the United States Air Force Academy. Divergence estimates from sampling dates indicate that the SARS-CoV-2 lineage associated with this rapid expansion event originated in late October 2020. These results are in agreement with transmission pathways inferred by the United States Air Force Academy, and provide a window into the evolutionary process and transmission dynamics of a potentially dangerous but ultimately contained variant.

A point mutation and large deletion at the candidate avirulence locus AvrMlp7 in the poplar rust fungus correlate with poplar RMlp7 resistance breakdown

The deployment of plant varieties carrying resistance genes (R) exerts strong selection pressure on pathogen populations. Rapidly evolving avirulence genes (Avr) allow pathogens to escape R-mediated plant immunity through a variety of mechanisms, leading to virulence. The poplar rust fungus Melampsora larici-populina is a damaging pathogen of poplars in Europe. It underwent a major adaptive event in 1994, with the breakdown of the poplar RMlp7 resistance gene. Population genomics studies identified a locus in the genome of M. larici-populina that probably corresponds to the candidate avirulence gene AvrMlp7. Here, to further characterize this effector, we used a population genetics approach on a comprehensive set of 281 individuals recovered throughout a 28-year period encompassing the resistance breakdown event. Using two dedicated molecular tools, genotyping at the candidate locus highlighted two different alterations of a predominant allele found mainly before the resistance breakdown: a nonsynonymous mutation and a complete deletion of this locus. This results in six diploid genotypes: three genotypes related to the avirulent phenotype and three related to the virulent phenotype. The temporal survey of the candidate locus revealed that both alterations were found in association during the resistance breakdown event. They pre-existed before the breakdown in a heterozygous state with the predominant allele cited above. Altogether, these results suggest that the association of both alterations at the candidate locus AvrMlp7 drove the poplar rust adaptation to RMlp7-mediated immunity. This study demonstrates for the first time a case of adaptation from standing genetic variation in rust fungi during a qualitative resistance breakdown.

Integrating CRISPR-Cas and Next Generation Sequencing in Plant Virology

Plant pathology has been revolutionized by the emergence and intervention of next-generation sequencing technologies (NGS) which provide a fast, cost-effective, and reliable diagnostic for any class of pathogens. NGS has made tremendous advancements in the area of research and diagnostics of plant infecting viromes and has bridged plant virology with other advanced research fields like genome editing technologies. NGS in a broader perspective holds the potential for plant health improvement by diagnosing and mitigating the new or unusual symptoms caused by novel/unidentified viruses. CRISPR-based genome editing technologies can enable rapid engineering of efficient viral/viroid resistance by directly targeting specific nucleotide sites of plant viruses and viroids. Critical genes such as eIf (iso) 4E or eIF4E have been targeted via the CRISPR platform to produce plants resistant to single-stranded RNA (ssRNA) viruses. CRISPR/Cas-based multi-target DNA or RNA tests can be used for rapid and accurate diagnostic assays for plant viruses and viroids. Integrating NGS with CRISPR-based genome editing technologies may lead to a paradigm shift in combating deadly disease-causing plant viruses/viroids at the genomic level. Furthermore, the newly discovered CRISPR/Cas13 system has unprecedented potential in plant viroid diagnostics and interference. In this review, we have highlighted the application and importance of sequencing technologies on covering the viral genomes for precise modulations. This review also provides a snapshot vision of emerging developments in NGS technologies for the characterization of plant viruses and their potential utilities, advantages, and limitations in plant viral diagnostics. Furthermore, some of the notable advances like novel virus-inducible CRISPR/Cas9 system that confers virus resistance with no off-target effects have been discussed.

The Epidemiology of Plant Virus Disease: Towards a New Synthesis

Epidemiology is the science of how disease develops in populations, with applications in human, animal and plant diseases. For plant diseases, epidemiology has developed as a quantitative science with the aims of describing, understanding and predicting epidemics, and intervening to mitigate their consequences in plant populations. Although the central focus of epidemiology is at the population level, it is often necessary to recognise the system hierarchies present by scaling down to the individual plant/cellular level and scaling up to the community/landscape level. This is particularly important for diseases caused by plant viruses, which in most cases are transmitted by arthropod vectors. This leads to range of virus-plant, virus-vector and vector-plant interactions giving a distinctive character to plant virus epidemiology (whilst recognising that some fungal, oomycete and bacterial pathogens are also vector-borne). These interactions have epidemiological, ecological and evolutionary consequences with implications for agronomic practices, pest and disease management, host resistance deployment, and the health of wild plant communities. Over the last two decades, there have been attempts to bring together these differing standpoints into a new synthesis, although this is more apparent for evolutionary and ecological approaches, perhaps reflecting the greater emphasis on shorter often annual time scales in epidemiological studies. It is argued here that incorporating an epidemiological perspective, specifically quantitative, into this developing synthesis will lead to new directions in plant virus research and disease management. This synthesis can serve to further consolidate and transform epidemiology as a key element in plant virus research.

Integrating mathematical models with experimental data to investigate the within-host dynamics of bacterial infections

Bacterial infections still constitute a major cause of mortality and morbidity worldwide. The unavailability of therapeutics, antimicrobial resistance and the chronicity of infections due to incomplete clearance contribute to this phenomenon. Despite the progress in antimicrobial and vaccine development, knowledge about the effect that therapeutics have on the host–bacteria interactions remains incomplete. Insights into the characteristics of bacterial colonization and migration between tissues and the relationship between replication and host- or therapeutically induced killing can enable efficient design of treatment approaches. Recently, innovative experimental techniques have generated data enabling the qualitative characterization of aspects of bacterial dynamics. Here, we argue that mathematical modeling as an adjunct to experimental data can enrich the biological insight that these data provide. However, due to limited interdisciplinary training, efforts to combine the two remain limited. To promote this dialogue, we provide a categorization of modeling approaches highlighting their relationship to data generated by a range of experimental techniques in the area of in vivo bacterial dynamics. We outline common biological themes explored using mathematical models with case studies across all pathogen classes. Finally, this review advocates multidisciplinary integration to improve our mechanistic understanding of bacterial infections and guide the use of existing or new therapies.

Genetic diversity of porcine reproductive and respiratory syndrome virus 1 in the United States of America from 2010 to 2018

Porcine reproductive and respiratory syndrome virus 1 (PRRSV-1) was first detected in the United States of America (USA) in 1999, several strains were also recognized soon later, and these isolates are typically called North American (NA) PRRSV-1. However, few reports have characterized PRRSV-1 viruses in the USA. We explored the genetic characteristics and diversity of PRRSV-1 viruses circulating in the USA. PRRSV-1 PCR-positive samples collected from seven states in 2010-2018 (n = 27) were subjected to next-generation sequencing. The 27 PRRSV-1 viruses had 88.4-91.3% nucleotide identity to the PRRSV-1 Lelystad-virus strain (the type 1 prototype strain) and 87.4-89.8% to the previously reported NA PRRSV-1 viruses. Individual proteins had several unique genetic characteristics and only one of the 27 tested samples had the characteristic 17-amino acid (aa) deletion in Nsp2, a genetic marker of NA PRRSV-1 viruses described previously. Fourteen isolates displayed a 3-aa C-terminal truncation in the highly conserved Nsp12 gene; 16 samples had a 21- or 18-aa C-terminal truncation in GP3 gene; and one was observed with a 1-aa deletion at the overlapping region of GP3 and GP4. In addition, the GP5 protein in most isolates, excluding one exception, demonstrated similar genetic variation as other reported NA PRRSV-1 isolates. All tested isolates clustered within subtype 1 together with other available NA PRRSV-1 viruses. Collectively, our results provide up-to-date information on PRRSV-1 viruses circulating in the USA in the past 9 years although the number of PRRSV-1 isolates included in this study is limited. These PRRSV-1 viruses have undergone gradual genetic variation and exhibited some previously undescribed genetic characteristics and diversity, which complicates the diagnosis and control of NA PRRSV-1.

Exploring the causes of small effective population sizes in cyst nematodes using artificial Globodera pallida populations

The effective size of a population is the size of an ideal population which would undergo genetic drift at the same rate as the real population. The balance between selection and genetic drift depends on the effective population size ( Ne), rather than the real numbers of individuals in the population ( N). The objectives of the present study were to estimate Ne in the potato cyst nematode Globodera pallida and to explore the causes of a low Ne/ N ratio in cyst nematodes using artificial populations. Using a temporal analysis of 24 independent populations, the median Ne was 58 individuals (min Ne = 25 and max Ne = 228). Ne is commonly lower than N but in the case of cyst nematodes, the Ne/ N ratio was extremely low. Using artificial populations showed that this low ratio did not result from migration, selection and overlapping generations, but could be explain by the fact that G. pallida populations deviate in structure from the assumptions of the ideal population by having unequal sex ratios, high levels of inbreeding and a high variance in family sizes. The consequences of a low Ne, resulting in a strong intensity of genetic drift, could be important for their control because G. pallida populations will have a low capacity to adapt to changing environments.

Impact of genetic drift, selection and accumulation level on virus adaptation to its host plants

The efficiency of plant major resistance genes is limited by the emergence and spread of resistance-breaking mutants. Modulation of the evolutionary forces acting on pathogen populations constitutes a promising way to increase the durability of these genes. We studied the effect of four plant traits affecting these evolutionary forces on the rate of resistance breakdown (RB) by a virus. Two of these traits correspond to virus effective population sizes (N_e ) at either plant inoculation or during infection. The third trait corresponds to differential selection exerted by the plant on the virus population. Finally, the fourth trait corresponds to within-plant virus accumulation (VA). These traits were measured experimentally on Potato virus Y (PVY) inoculated to a set of 84 pepper doubled-haploid lines, all carrying the same pvr2³ resistance gene, but having contrasting genetic backgrounds. The lines showed extensive variation for the rate of pvr2³ RB by PVY and for the four other traits of interest. A generalized linear model showed that three of these four traits, with the exception of N_e at inoculation, and several pairwise interactions between them had significant effects on RB. RB increased with increasing values of N_e during plant infection or VA. The effect of differential selection was more complex because of a strong interaction with VA. When VA was high, RB increased as the differential selection increased. An opposite relationship between RB and differential selection was observed when VA was low. This study provides a framework to select plants with appropriate virus evolution-related traits to avoid or delay RB.

Estimating virus effective population size and selection without neutral markers

By combining high-throughput sequencing (HTS) with experimental evolution, we can observe the within-host dynamics of pathogen variants of biomedical or ecological interest. We studied the evolutionary dynamics of five variants of Potato virus Y (PVY) in 15 doubled-haploid lines of pepper. All plants were inoculated with the same mixture of virus variants and variant frequencies were determined by HTS in eight plants of each pepper line at each of six sampling dates. We developed a method for estimating the intensities of selection and genetic drift in a multi-allelic Wright-Fisher model, applicable whether these forces are strong or weak, and in the absence of neutral markers. This method requires variant frequency determination at several time points, in independent hosts. The parameters are the selection coefficients for each PVY variant and four effective population sizes N_e at different time-points of the experiment. Numerical simulations of asexual haploid Wright-Fisher populations subjected to contrasting genetic drift (N_e ∈ [10, 2000]) and selection (|s| ∈ [0, 0.15]) regimes were used to validate the method proposed. The experiment in closely related pepper host genotypes revealed that viruses experienced a considerable diversity of selection and genetic drift regimes. The resulting variant dynamics were accurately described by Wright-Fisher models. The fitness ranks of the variants were almost identical between host genotypes. By contrast, the dynamics of N_e were highly variable, although a bottleneck was often identified during the systemic movement of the virus. We demonstrated that, for a fixed initial PVY population, virus effective population size is a heritable trait in plants. These findings pave the way for the breeding of plant varieties exposing viruses to stronger genetic drift, thereby slowing virus adaptation.

A growing number of experimental evolution studies are using an “evolve-and-resequence” approach to observe the within-host dynamics of pathogen variants of biomedical or ecological interest. The resulting data are particularly appropriate for studying the effects of evolutionary forces, such as selection and genetic drift, on the emergence of new pathogen variants. However, it remains challenging to unravel the effects of selection and genetic drift in the absence of neutral markers, a situation frequently encountered for microbes, such as viruses, due to their small constrained genomes. Using such an approach on a plant virus, we observed that the same set of virus variants displayed highly diverse dynamics in closely related plant genotypes. We developed and validated a method that does not require neutral markers, for estimating selection coefficients and effective population sizes from these experimental evolution data. We found that the viruses experienced considerable diversity in genetic drift regimes, depending on host genotype. Importantly, genetic drift experienced by virus populations was shown to be a heritable plant trait. These findings pave the way for the breeding of plant varieties exposing viruses to strong genetic drift, thereby slowing virus adaptation.

Analysis of intra-host genetic diversity of Prunus necrotic ringspot virus (PNRSV) using amplicon next generation sequencing

PCR amplicon next generation sequencing (NGS) analysis offers a broadly applicable and targeted approach to detect populations of both high- or low-frequency virus variants in one or more plant samples. In this study, amplicon NGS was used to explore the diversity of the tripartite genome virus, Prunus necrotic ringspot virus (PNRSV) from 53 PNRSV-infected trees using amplicons from conserved gene regions of each of PNRSV RNA1, RNA2 and RNA3. Sequencing of the amplicons from 53 PNRSV-infected trees revealed differing levels of polymorphism across the three different components of the PNRSV genome with a total number of 5040, 2083 and 5486 sequence variants observed for RNA1, RNA2 and RNA3 respectively. The RNA2 had the lowest diversity of sequences compared to RNA1 and RNA3, reflecting the lack of flexibility tolerated by the replicase gene that is encoded by this RNA component. Distinct PNRSV phylo-groups, consisting of closely related clusters of sequence variants, were observed in each of PNRSV RNA1, RNA2 and RNA3. Most plant samples had a single phylo-group for each RNA component. Haplotype network analysis showed that smaller clusters of PNRSV sequence variants were genetically connected to the largest sequence variant cluster within a phylo-group of each RNA component. Some plant samples had sequence variants occurring in multiple PNRSV phylo-groups in at least one of each RNA and these phylo-groups formed distinct clades that represent PNRSV genetic strains. Variants within the same phylo-group of each Prunus plant sample had ≥97% similarity and phylo-groups within a Prunus plant sample and between samples had less ≤97% similarity. Based on the analysis of diversity, a definition of a PNRSV genetic strain was proposed. The proposed definition was applied to determine the number of PNRSV genetic strains in each of the plant samples and the complexity in defining genetic strains in multipartite genome viruses was explored.

High-throughput whole genome sequencing of Porcine reproductive and respiratory syndrome virus from cell culture materials and clinical specimens using next-generation sequencing technology

Next-generation sequencing (NGS) technologies have increasingly played crucial roles in biological and medical research, but are not yet in routine use in veterinary diagnostic laboratories. We developed and applied a procedure for high-throughput RNA sequencing of Porcine reproductive and respiratory syndrome virus (PRRSV) from cell culture-derived isolates and clinical specimens. Ten PRRSV isolates with known sequence information, 2 mixtures each with 2 different PRRSV isolates, and 51 clinical specimens (19 sera, 16 lungs, and 16 oral fluids) with various PCR threshold cycle (Ct) values were subjected to nucleic acid extraction, cDNA library preparation (24-plexed), and sequencing. Whole genome sequences were obtained from 10 reference isolates with expected sequences and from sera with a PRRSV real-time reverse transcription PCR Ct ≤ 23.6, lung tissues with Ct ≤ 21, and oral fluids with Ct ≤ 20.6. For mixtures with PRRSV-1 and -2 isolates (57.8% nucleotide identity), NGS was able to distinguish them as well as obtain their respective genome sequences. For mixtures with 2 PRRSV-2 isolates (92.4% nucleotide identity), sequence reads with nucleotide ambiguity at numerous sites were observed, indicating mixed infection; however, individual virus sequences could only be separated when 1 isolate identity and sequence in the mixture is known. The NGS approach described herein offers the prospect of high-throughput sequencing and could be adapted to routine workflows in veterinary diagnostic laboratories, although further improvement of sequencing outcomes from clinical specimens with higher Ct values remains to be investigated.

Mutations in Rice yellow mottle virus Polyprotein P2a Involved in RYMV2 Gene Resistance Breakdown

Rice yellow mottle virus (RYMV) is one of the major diseases of rice in Africa. The high resistance of the Oryza glaberrima Tog7291 accession involves a null allele of the RYMV2 gene, whose ortholog in Arabidopsis, CPR5, is a transmembrane nucleoporin involved in effector-triggered immunity. To optimize field deployment of the RYMV2 gene and improve its durability, which is often a weak point in varietal resistance, we analyzed its efficiency toward RYMV isolates representing the genetic diversity of the virus and the molecular basis of resistance breakdown. Tog7291 resistance efficiency was highly variable depending on the isolate used, with infection rates ranging from 0 to 98% of plants. Back-inoculation experiments indicated that infection cases were not due to an incomplete resistance phenotype but to the emergence of resistance-breaking (RB) variants. Interestingly, the capacity of the virus to overcome Tog7291 resistance is associated with a polymorphism at amino-acid 49 of the VPg protein which also affects capacity to overcome the previously studied RYMV1 resistance gene. This polymorphism appeared to be a main determinant of the emergence of RB variants. It acts independently of the resistance gene and rather reflects inter-species adaptation with potential consequences for the durability of resistance. RB mutations were identified by full-length or partial sequencing of the RYMV genome in infected Tog7291 plants and were validated by directed mutagenesis of an infectious viral clone. We found that Tog7291 resistance breakdown involved mutations in the putative membrane anchor domain of the polyprotein P2a. Although the precise effect of these mutations on rice/RYMV interaction is still unknown, our results offer a new perspective for the understanding of RYMV2 mediated resistance mechanisms. Interestingly, in the susceptible IR64 variety, RB variants showed low infectivity and frequent reversion to the wild-type genotype, suggesting that Tog7291 resistance breakdown is associated with a major loss of viral fitness in normally susceptible O. sativa varieties. Despite the high frequency of resistance breakdown in controlled conditions, this loss of fitness is an encouraging element with regards to RYMV2 resistance durability.

Human Management of a Wild Plant Modulates the Evolutionary Dynamics of a Gene Determining Recessive Resistance to Virus Infection

This work analyses the genetic variation and evolutionary patterns of recessive resistance loci involved in matching-allele (MA) host-pathogen interactions, focusing on the pvr2 resistance gene to potyviruses of the wild pepper Capsicum annuum glabriusculum (chiltepin). Chiltepin grows in a variety of wild habitats in Mexico, and its cultivation in home gardens started about 25 years ago. Potyvirus infection of Capsicum plants requires the physical interaction of the viral VPg with the pvr2 product, the translation initiation factor eIF4E1. Mutations impairing this interaction result in resistance, according to the MA model. The diversity of pvr2/eIF4E1 in wild and cultivated chiltepin populations from six biogeographical provinces in Mexico was analysed in 109 full-length coding sequences from 97 plants. Eleven alleles were found, and their interaction with potyvirus VPg in yeast-two-hybrid assays, plus infection assays of plants, identified six resistance alleles. Mapping resistance mutations on a pvr2/eIF4E1 model structure showed that most were around the cap-binding pocket and strongly altered its surface electrostatic potential, suggesting resistance-associated costs due to functional constraints. The pvr2/eIF4E1 phylogeny established that susceptibility was ancestral and resistance was derived. The spatial structure of pvr2/eIF4E1 diversity differed from that of neutral markers, but no evidence of selection for resistance was found in wild populations. In contrast, the resistance alleles were much more frequent, and positive selection stronger, in cultivated chiltepin populations, where diversification of pvr2/eIF4E1 was higher. This analysis of the genetic variation of a recessive resistance gene involved in MA host-pathogen interactions in populations of a wild plant show that evolutionary patterns differ according to the plant habitat, wild or cultivated. It also demonstrates that human management of the plant population has profound effects on the diversity and the evolution of the resistance gene, resulting in the selection of resistance alleles.

Temporal sampling helps unravel the genetic structure of naturally occurring populations of a phytoparasitic nematode. 1. Insights from the estimation of effective population sizes

The sustainability of modern agriculture relies on strategies that can control the ability of pathogens to overcome chemicals or genetic resistances through natural selection. This evolutionary potential, which depends partly on effective population size (N e ), is greatly influenced by human activities. In this context, wild pathogen populations can provide valuable information for assessing the long-term risk associated with crop pests. In this study, we estimated the effective population size of the beet cyst nematode, Heterodera schachtii, by sampling 34 populations infecting the sea beet Beta vulgaris spp. maritima twice within a one-year period. Only 20 populations produced enough generations to analyze the variation in allele frequencies, with the remaining populations showing a high mortality rate of the host plant after only 1 year. The 20 analyzed populations showed surprisingly low effective population sizes, with most having N e close to 85 individuals. We attribute these low values to the variation in population size through time, systematic inbreeding, and unbalanced sex-ratios. Our results suggest that H. schachtii has low evolutionary potential in natural environments. Pest control strategies in which populations on crops mimic wild populations may help prevent parasite adaptation to host resistance.

Evolutionary and Epidemiological Implications of Multiple Infection in Plants

Recent methodological advances have uncovered tremendous microbial diversity cohabiting in the same host plant, and many of these microbes cause disease. In this review we highlight how the presence of other pathogen species, or other pathogen genotypes, within a plant can affect key components of host-pathogen interactions: (i) within-plant virulence and pathogen accumulation, through direct and host-mediated mechanisms; (ii) evolutionary trajectories of pathogen populations, through virulence evolution, generation of novel genetic combinations, and maintenance of genetic diversity; and (iii) disease dynamics, with multiple infection likely to render epidemics more devastating. The major future challenges are to couple a community ecology approach with a molecular investigation of the mechanisms operating under coinfection and to evaluate the evolution and effectiveness of resistance within a coinfection framework.

Epidemiological and evolutionary management of plant resistance: optimizing the deployment of cultivar mixtures in time and space in agricultural landscapes

The management of genes conferring resistance to plant–pathogens should make it possible to control epidemics (epidemiological perspective) and preserve resistance durability (evolutionary perspective). Resistant and susceptible cultivars must be strategically associated according to the principles of cultivar mixture (within a season) and rotation (between seasons). We explored these questions by modeling the evolutionary and epidemiological processes shaping the dynamics of a pathogen population in a landscape composed of a seasonal cultivated compartment and a reservoir compartment hosting pathogen year-round. Optimal deployment strategies depended mostly on the molecular basis of plant–pathogen interactions and on the agro-ecological context before resistance deployment, particularly epidemic intensity and landscape connectivity. Mixtures were much more efficient in landscapes in which between-field infections and infections originating from the reservoir were more prevalent than within-field infections. Resistance genes requiring two mutations of the pathogen avirulence gene to be broken down, rather than one, were particularly useful when infections from the reservoir predominated. Combining mixture and rotation principles were better than the use of the same mixture each season as (i) they controlled epidemics more effectively in situations in which within-field infections or infections from the reservoir were frequent and (ii) they fulfilled the epidemiological and evolutionary perspectives.

Deep sequencing of virus-derived small interfering RNAs and RNA from viral particles shows highly similar mutational landscapes of a plant virus population

RNA viruses exist within a host as a population of mutant sequences, often referred to as quasispecies. Within a host, sequences of RNA viruses constitute several distinct but interconnected pools, such as RNA packed in viral particles, double-stranded RNA, and virus-derived small interfering RNAs. We aimed to test if the same representation of within-host viral population structure could be obtained by sequencing different viral sequence pools. Using ultradeep Illumina sequencing, the diversity of two coexisting Potato virus Y sequence pools present within a plant was investigated: RNA isolated from viral particles and virus-derived small interfering RNAs (the derivatives of a plant RNA silencing mechanism). The mutational landscape of the within-host virus population was highly similar between both pools, with no notable hotspots across the viral genome. Notably, all of the single-nucleotide polymorphisms with a frequency of higher than 1.6% were found in both pools. Some unique single-nucleotide polymorphisms (SNPs) with very low frequencies were found in each of the pools, with more of them occurring in the small RNA (sRNA) pool, possibly arising through genetic drift in localized virus populations within a plant and the errors introduced during the amplification of silencing signal. Sequencing of the viral particle pool enhanced the efficiency of consensus viral genome sequence reconstruction. Nonhomologous recombinations were commonly detected in the viral particle pool, with a hot spot in the 3' untranslated and coat protein regions of the genome. We stress that they present an important but often overlooked aspect of virus population diversity.

IMPORTANCE

This study is the most comprehensive whole-genome characterization of a within-plant virus population to date and the first study comparing diversity of different pools of viral sequences within a host. We show that both virus-derived small RNAs and RNA from viral particles could be used for diversity assessment of within-plant virus population, since they show a highly congruent portrayal of the virus mutational landscape within a plant. The study is an important baseline for future studies of virus population dynamics, for example, during the adaptation to a new host. The comparison of the two virus sequence enrichment techniques, sequencing of virus-derived small interfering RNAs and RNA from purified viral particles, shows the strength of the latter for the detection of recombinant viral genomes and reconstruction of complete consensus viral genome sequence.

Analysis of viral (zucchini yellow mosaic virus) genetic diversity during systemic movement through a Cucurbita pepo vine

Determining the extent and structure of intra-host genetic diversity and the magnitude and impact of population bottlenecks is central to understanding the mechanisms of viral evolution. To determine the nature of viral evolution following systemic movement through a plant, we performed deep sequencing of 23 leaves that grew sequentially along a single Cucurbita pepo vine that was infected with zucchini yellow mosaic virus (ZYMV), and on a leaf that grew in on a side branch. Strikingly, of 112 genetic (i.e. sub-consensus) variants observed in the data set as a whole, only 22 were found in multiple leaves. Similarly, only three of the 13 variants present in the inoculating population were found in the subsequent leaves on the vine. Hence, it appears that systemic movement is characterized by sequential population bottlenecks, although not sufficient to reduce the population to a single virion as multiple variants were consistently transmitted between leaves. In addition, the number of variants within a leaf increases as a function of distance from the inoculated (source) leaf, suggesting that the circulating sap may serve as a continual source of virus. Notably, multiple mutational variants were observed in the cylindrical inclusion (CI) protein (known to be involved in both cell-to-cell and systemic movement of the virus) that were present in multiple (19/24) leaf samples. These mutations resulted in a conformational change, suggesting that they might confer a selective advantage in systemic movement within the vine. Overall, these data reveal that bottlenecks occur during systemic movement, that variants circulate in the phloem sap throughout the infection process, and that important conformational changes in CI protein may arise during individual infections.

Interaction patterns between potato virus Y and eIF4E-mediated recessive resistance in the Solanaceae

The structural pattern of infectivity matrices, which contains infection data resulting from inoculations of a set of hosts by a set of parasites, is a key parameter for our understanding of biological interactions and their evolution. This pattern determines the evolution of parasite pathogenicity and host resistance, the spatiotemporal distribution of host and parasite genotypes, and the efficiency of disease control strategies. Two major patterns have been proposed for plant-virus genotype infectivity matrices. In the gene-for-gene model, infectivity matrices show a nested pattern, where the host ranges of specialist virus genotypes are subsets of the host ranges of less specialized viruses. In contrast, in the matching-allele (MA) model, each virus genotype is specialized to infect one (or a small set of) host genotype(s). The corresponding infectivity matrix shows a modular pattern where infection is frequent for plants and viruses belonging to the same module but rare for those belonging to different modules. We analyzed the structure of infectivity matrices between Potato virus Y (PVY) and plant genotypes in the family Solanaceae carrying different eukaryotic initiation factor 4E (eIF4E)-coding alleles conferring recessive resistance. Whereas this system corresponds mechanistically to an MA model, the expected modular pattern was rejected based on our experimental data. This was mostly because PVY mutations involved in adaptation to a particular plant genotype displayed frequent pleiotropic effects, conferring simultaneously an adaptation to additional plant genotypes with different eIF4E alleles. Such effects should be taken into account for the design of strategies of sustainable control of PVY through plant varietal mixtures or rotations.

IMPORTANCE

The interaction pattern between host and virus genotypes has important consequences on their respective evolution and on issues regarding the application of disease control strategies. We found that the structure of the interaction between Potato virus Y (PVY) variants and host plants in the family Solanaceae departs significantly from the current model of interaction considered for these organisms because of frequent pleiotropic effects of virus mutations. These mutational effects allow the virus to expand rapidly its range of host plant genotypes, make it very difficult to predict the effects of mutations in PVY infectivity factors, and raise concerns about strategies of sustainable management of plant genetic resistance to viruses.

Current impact and future directions of high throughput sequencing in plant virus diagnostics

The ability to provide a fast, inexpensive and reliable diagnostic for any given viral infection is a key parameter in efforts to fight and control these ubiquitous pathogens. The recent developments of high-throughput sequencing (also called Next Generation Sequencing - NGS) technologies and bioinformatics have drastically changed the research on viral pathogens. It is now raising a growing interest for virus diagnostics. This review provides a snapshot vision on the current use and impact of high throughput sequencing approaches in plant virus characterization. More specifically, this review highlights the potential of these new technologies and their interplay with current protocols in the future of molecular diagnostic of plant viruses. The current limitations that will need to be addressed for a wider adoption of high-throughput sequencing in plant virus diagnostics are thoroughly discussed.

Historical perspective, development and applications of next-generation sequencing in plant virology

Next-generation high throughput sequencing technologies became available at the onset of the 21st century. They provide a highly efficient, rapid, and low cost DNA sequencing platform beyond the reach of the standard and traditional DNA sequencing technologies developed in the late 1970s. They are continually improved to become faster, more efficient and cheaper. They have been used in many fields of biology since 2004. In 2009, next-generation sequencing (NGS) technologies began to be applied to several areas of plant virology including virus/viroid genome sequencing, discovery and detection, ecology and epidemiology, replication and transcription. Identification and characterization of known and unknown viruses and/or viroids in infected plants are currently among the most successful applications of these technologies. It is expected that NGS will play very significant roles in many research and non-research areas of plant virology.

Narrow bottlenecks affect Pea seedborne mosaic virus populations during vertical seed transmission but not during leaf colonization

The effective size of populations (Ne) determines whether selection or genetic drift is the predominant force shaping their genetic structure and evolution. Populations having high Ne adapt faster, as selection acts more intensely, than populations having low Ne, where random effects of genetic drift dominate. Estimating Ne for various steps of plant virus life cycle has been the focus of several studies in the last decade, but no estimates are available for the vertical transmission of plant viruses, although virus seed transmission is economically significant in at least 18% of plant viruses in at least one plant species. Here we study the co-dynamics of two variants of Pea seedborne mosaic virus (PSbMV) colonizing leaves of pea plants (Pisum sativum L.) during the whole flowering period, and their subsequent transmission to plant progeny through seeds. Whereas classical estimators of Ne could be used for leaf infection at the systemic level, as virus variants were equally competitive, dedicated stochastic models were needed to estimate Ne during vertical transmission. Very little genetic drift was observed during the infection of apical leaves, with Ne values ranging from 59 to 216. In contrast, a very drastic genetic drift was observed during vertical transmission, with an average number of infectious virus particles contributing to the infection of a seedling from an infected mother plant close to one. A simple model of vertical transmission, assuming a cumulative action of virus infectious particles and a virus density threshold required for vertical transmission to occur fitted the experimental data very satisfactorily. This study reveals that vertically-transmitted viruses endure bottlenecks as narrow as those imposed by horizontal transmission. These bottlenecks are likely to slow down virus adaptation and could decrease virus fitness and virulence.

Applications of next generation high throughput sequencing technologies in characterization, discovery and molecular interaction of plant viruses

Present era of molecular biology is witnessing revolutionary developments in sequencing technology. This advancement has considerably influenced plant virology in the field of diagnostics and host virus interaction. Next generation high-throughput sequencing technology has made it possible to directly detect, identify and discover novel viruses in several plants in an unbiased manner without antibodies or prior knowledge of the virus sequences. Entire viral genome could be sequenced from symptomatic or asymptomatic plants through next generation sequencing of total nucleic acids including small RNAs. It provides census of viral population in a particular ecosystem or cropping system. Viral genome variability, evolution within the host and virus defence mechanism in plants can also be easily understood by massive parallel sequencing. In this article, we provide an overview of the applications of next generation sequencing technology in characterization, discovery and molecular interaction of plant viruses.

Potato virus Y: a major crop pathogen that has provided major insights into the evolution of viral pathogenicity

TAXONOMY

Potato virus Y (PVY) is the type member of the genus Potyvirus in the family Potyviridae. VIRION AND GENOME PROPERTIES: PVY virions have a filamentous, flexuous form, with a length of 730 nm and a diameter of 12 nm. The genomic RNA is single stranded, messenger sense, with a length of 9.7 kb, covalently linked to a viral-encoded protein (VPg) at the 5' end and to a 3' polyadenylated tail. The genome is expressed as a polyprotein of approximately 3062 amino acid residues, processed by three virus-specific proteases into 11 mature proteins.

HOSTS

PVY is distributed worldwide and has a broad host range, consisting of cultivated solanaceous species and many solanaceous and nonsolanaceous weeds. It is one of the most economically important plant pathogens and causes severe diseases in cultivated hosts, such as potato, tobacco, tomato and pepper, as well as in ornamental plants.

TRANSMISSION

PVY is transmitted from plant to plant by more than 40 aphid species in a nonpersistent manner and, in potato, by planting contaminated seed tubers. DIVERSITY: Five major clades, named C1, C2, Chile, N and O, have been described within the PVY species. In recent decades, a strong increase in prevalence of N × O recombinant isolates has been observed worldwide. A correlation has been observed between PVY phylogeny and certain pathogenicity traits. GENETIC CONTROL OF PVY: Resistance genes against PVY have been used widely in breeding programmes and deployed in the field. These resistance genes show a large diversity of spectrum of action, durability and genetic determinism. Notably, recessive and dominant major resistance genes show highly contrasting patterns of interaction with PVY populations, displaying rapid co-evolution or stable relationships, respectively.

Within-host dynamics of the emergence of Tomato yellow leaf curl virus recombinants

Tomato yellow leaf curl virus (TYLCV) is a highly damaging begomovirus native to the Middle East. TYLCV has recently spread worldwide, recombining with other begomoviruses. Recent analysis of mixed infections between TYLCV and Tomato leaf curl Comoros begomovirus (ToLCKMV) has shown that, although natural selection preserves certain co-evolved intra-genomic interactions, numerous and diverse recombinants are produced at 120 days post-inoculation (dpi), and recombinant populations from different tomato plants are very divergent. Here, we investigate the population dynamics that lead to such patterns in tomato plants co-infected with TYLCV and ToLCKMV either by agro-inoculation or using the natural whitefly vector Bemisia tabaci. We monitored the frequency of parental and recombinant genotypes independently in 35 plants between 18 and 330 dpi and identified 177 recombinants isolated at different times. Recombinants were detected from 18 dpi and their frequency increased over time to reach about 50% at 150 dpi regardless of the inoculation method. The distribution of breakpoints detected on 96 fully sequenced recombinants was consistent with a continuous generation of new recombinants as well as random and deterministic effects in their maintenance. A severe population bottleneck of around 10 genomes was estimated during early systemic infection–a phenomenon that could account partially for the heterogeneity in recombinant patterns observed among plants. The detection of the same recombinant genome in six of the thirteen plants analysed beyond 30 dpi supported the influence of selection on observed recombination patterns. Moreover, a highly virulent recombinant genotype dominating virus populations within one plant has, apparently, the potential to be maintained in the natural population according to its infectivity, within-host accumulation, and transmission efficiency - all of which were similar or intermediate to those of the parent genotypes. Our results anticipate the outcomes of natural encounters between TYLCV and ToLCKMV.

Farther, slower, stronger: how the plant genetic background protects a major resistance gene from breakdown

Genetic resistance provides efficient control of crop diseases, but is limited by pathogen evolution capacities which often result in resistance breakdown. It has been demonstrated recently, in three different pathosystems, that polygenic resistances combining a major-effect gene and quantitative resistance controlled by the genetic background are more durable than monogenic resistances (with the same major gene in a susceptible genetic background), but the underlying mechanisms are unknown. Using the pepper-Potato virus Y system, we examined three mechanisms that could account for the greater durability of the polygenic resistances: (i) the additional quantitative resistance conferred by the genetic background; (ii) the increase in the number of mutations required for resistance breakdown; and (iii) the slower selection of adapted resistance-breaking mutants within the viral population. The three mechanisms were experimentally validated. The first explained a large part of the variation in resistance breakdown frequency and is therefore expected to be a major determinant of resistance durability. Quantitative resistance factors also had an influence on the second mechanism by modifying the virus mutational pathways towards resistance breakdown and could also have an influence on the third mechanism by increasing genetic drift effects on the viral population. The relevance of these results for other plant-pathogen systems and their importance in plant breeding are discussed.

Experimental measures of pathogen competition and relative fitness

Competition among pathogen strains for limited host resources can have a profound effect on pathogen evolution. A better understanding of the principles and consequences of competition can be useful in designing more sustainable disease management strategies. The competitive ability and relative fitness of a pathogen strain are determined by its intrinsic biological properties, the resistance and heterogeneity of the corresponding host population, the population density and genetic relatedness of the competing strains, and the physical environment. Competitive ability can be inferred indirectly from fitness components, such as basic reproduction rate or transmission rate. However, pathogen strains that exhibit higher fitness components when they infect a host alone may not exhibit a competitive advantage when they co-infect the same host. The most comprehensive measures of competitive ability and relative fitness come from calculating selection coefficients in a mixed infection in a field setting. Mark-release-recapture experiments can be used to estimate fitness costs associated with unnecessary virulence and fungicide resistance.

Dynamics of the establishment of systemic Potyvirus infection: independent yet cumulative action of primary infection sites

In the clinic, farm, or field, for many viruses there is a high prevalence of mixed-genotype infections, indicating that multiple virions have initiated infection and that there can be multiple sites of primary infection within the same host. The dynamic process by which multiple primary infection sites interact with each other and the host is poorly understood, undoubtedly due to its high complexity. In this study, we attempted to unravel the basic interactions underlying this process using a plant RNA virus, as removing the inoculated leaf can instantly and rigorously eliminate all primary infection sites. Effective population size in the inoculated leaf and time of removal of the inoculated leaf were varied in experiments, and it was found that both factors positively influenced if the plant became systemically infected and what proportion of cells in the systemic tissue were infected, as measured by flow cytometry. Fitting of probabilistic models of infection to our data demonstrated that a null model in which the action of each focus is independent of the presence of other foci was better supported than a dependent-action model. The cumulative effect of independently acting foci therefore determined when plants became infected and how many individual cells were infected. There was no evidence for interference between primary infection sites, which is surprising given the planar structure of leaves. By showing that a simple null model is supported, we experimentally confirmed--to our knowledge for the first time--the minimal components that dictate interactions of a conspecific virus population establishing systemic infection.