The transcriptomics of an experimentally evolved plant-virus interaction

Transcriptomic insights into the epigenetic modulation of turnip mosaic virus evolution in Arabidopsis thaliana

BACKGROUND

Plant-virus interaction models propose that a virus's ability to infect a host genotype depends on the compatibility between virulence and resistance genes. Recently, we conducted an evolution experiment in which lineages of turnip mosaic virus (TuMV) were passaged in Arabidopsis thaliana genotypes carrying mutations in components of the DNA methylation and the histone demethylation epigenetic pathways. All evolved lineages increased infectivity, virulence and viral load in a host genotype-dependent manner.

RESULTS

To better understand the underlying reasons for these evolved relationships, we delved into the transcriptomic responses of mutant and WT plant genotypes in mock conditions and infected with either the ancestral or evolved viruses. Such a comparison allowed us to classify every gene into nine basic expression profiles. Regarding the targets of viral adaptation, our analyses allowed the identification of common viral targets as well as host genotype-specific genes and categories of biological processes. As expected, immune response-related genes were found to be altered upon infection. However, we also noticed the pervasive over-representation of other functional groups, suggesting that viral adaptation was not solely driven by the level of expression of plant resistance genes. In addition, a significant association between the presence of transposable elements within or upstream the differentially expressed genes was observed. Finally, integration of transcriptomic data into a virus-host protein-protein interaction network highlighted the most impactful interactions.

CONCLUSIONS

These findings shed extra light on the complex dynamics between plants and viruses, indicating that viral infectivity depends on various factors beyond just the plant's resistance genes.

Transcriptional and hormonal profiling uncovers the interactions between plant developmental stages and RNA virus infection

Arabidopsis thaliana is more susceptible to certain viruses during its later developmental stages. The differential responses and the mechanisms behind this development-dependent susceptibility to infection are still not fully understood. Here we explored the outcome of a viral infection at different host developmental stages by studying the response of A. thaliana to infection with turnip mosaic virus at three developmental stages: juvenile vegetative, bolting, and mature flowering plants. We found that infected plants at later stages downregulate cell wall biosynthetic genes and that this downregulation may be one factor facilitating viral spread and systemic infection. We also found that, despite being more susceptible to infection, infected mature flowering plants were more fertile (i.e. produce more viable seeds) than juvenile vegetative and bolting infected plants; that is, plants infected at the reproductive stage have greater fitness than plants infected at earlier developmental stages. Moreover, treatment of mature plants with salicylic acid increased resistance to infection at the cost of significantly reducing fertility. Together, these observations support a negative trade-off between viral susceptibility and plant fertility. Our findings point towards a development-dependent tolerance to infection.

Historical Perspective on the Discovery of the Quasispecies Concept

Viral quasispecies are dynamic distributions of nonidentical but closely related mutant and recombinant viral genomes subjected to a continuous process of genetic variation, competition, and selection that may act as a unit of selection. The quasispecies concept owes its theoretical origins to a model for the origin of life as a collection of mutant RNA replicators. Independently, experimental evidence for the quasispecies concept was obtained from sampling of bacteriophage clones, which revealed that the viral populations consisted of many mutant genomes whose frequency varied with time of replication. Similar findings were made in animal and plant RNA viruses. Quasispecies became a theoretical framework to understand viral population dynamics and adaptability. The evidence came at a time when mutations were considered rare events in genetics, a perception that was to change dramatically in subsequent decades. Indeed, viral quasispecies was the conceptual forefront of a remarkable degree of biological diversity, now evident for cell populations and organisms, not only for viruses. Quasispecies dynamics unveiled complexities in the behavior of viral populations,with consequences for disease mechanisms and control strategies. This review addresses the origin of the quasispecies concept, its major implications on both viral evolution and antiviral strategies, and current and future prospects.

A genome-wide association study identifies Arabidopsis thaliana genes that contribute to differences in the outcome of infection with two Turnip mosaic potyvirus strains that differ in their evolutionary history and degree of host specialization

Viruses lie in a continuum between generalism and specialism depending on their ability to infect more or less hosts. While generalists are able to successfully infect a wide variety of hosts, specialists are limited to one or a few. Even though generalists seem to gain an advantage due to their wide host range, they usually pay a pleiotropic fitness cost within each host. On the contrary, a specialist has maximal fitness within its own host. A relevant yet poorly explored question is whether viruses differ in the way they interact with their hosts' gene expression depending on their degree of specialization. Using a genome-wide association study approach, we have identified host genes whose expression depends on whether hosts were infected with more or less specialized viral strains. Four hundred fifty natural accessions of Arabidopsis thaliana were inoculated with Turnip mosaic potyvirus strains with different past evolutionary histories and that shown different degrees of specialization. Three disease-related traits were measured and associated with different sets of host genes for each strain. The genetic architectures of these traits differed among viral strains and, in the case of the more specialized virus, also varied along the duration of infection. While most of the mapped loci were strain specific, one shared locus was mapped for both strains, a disease-resistance TIR-NBS-LRR class protein. Likewise, only putative cysteine-rich receptor-like protein kinases were involved in all three traits. The impact on disease progress of 10 selected genes was validated by studying the infection phenotypes of loss-of-function mutant plants. Nine of these mutants have altered the disease progress and/or symptoms intensity between both strains. Compared to wild-type plants six had an effect on both viral strains, three had an effect only on the more specialized, and two were significant during infection with the less specialized.

The Oxylipin Signaling Pathway Is Required for Increased Aphid Attraction and Retention on Virus-Infected Plants

Many studies have shown that virus infection alters phytohormone signaling and insect vector contact with hosts. Increased vector contact and movement among plants should increase virus survival and host range. In this study we examine the role of virus-induced changes in phytohormone signaling in plant-aphid interactions, using Pea enation mosaic virus (PEMV), pea aphids (Acyrthosiphon pisum), and pea (Pisum sativum) as a model. We observed that feeding by aphids carrying PEMV increases salicylic acid and jasmonic acid accumulation in pea plants compared to feeding by virus-free aphids. To determine if induction of the oxylipin jasmonic acid is critical for aphid settling, attraction, and retention on PEMV-infected plants, we conducted insect bioassays using virus-induced gene silencing (VIGS), an oxylipin signaling inducer, methyl jasmonate (MeJA), and a chemical inhibitor of oxylipin signaling, phenidone. Surprisingly, there was no impact of phenidone treatment on jasmonic acid or salicylic acid levels in virus-infected plants, though aphid attraction and retention were altered. These results suggest that the observed impacts of phenidone on aphid attraction to and retention on PEMV-infected plants are independent of the jasmonic acid and salicylic acid pathway but may be mediated by another component of the oxylipin signaling pathway. These results shed light on the complexity of viral manipulation of phytohormone signaling and vector-plant interactions.

Viral Fitness Determines the Magnitude of Transcriptomic and Epigenomic Reprograming of Defense Responses in Plants

Although epigenetic factors may influence the expression of defense genes in plants, their role in antiviral responses and the impact of viral adaptation and evolution in shaping these interactions are still poorly explored. We used two isolates of turnip mosaic potyvirus with varying degrees of adaptation to Arabidopsis thaliana to address these issues. One of the isolates was experimentally evolved in the plant and presented increased load and virulence relative to the ancestral isolate. The magnitude of the transcriptomic responses was larger for the evolved isolate and indicated a role of innate immunity systems triggered by molecular patterns and effectors in the infection process. Several transposable elements located in different chromatin contexts and epigenetic-related genes were also affected. Correspondingly, mutant plants having loss or gain of repressive marks were, respectively, more tolerant and susceptible to turnip mosaic potyvirus, with a more efficient response against the ancestral isolate. In wild-type plants, both isolates induced similar levels of cytosine methylation changes, including in and around transposable elements and stress-related genes. Results collectively suggested that apart from RNA silencing and basal immunity systems, DNA methylation and histone modification pathways may also be required for mounting proper antiviral defenses and that the effectiveness of this type of regulation strongly depends on the degree of viral adaptation to the host.

Impact of Mutations in Arabidopsis thaliana Metabolic Pathways on Polerovirus Accumulation, Aphid Performance, and Feeding Behavior

During the process of virus acquisition by aphids, plants respond to both the virus and the aphids by mobilizing different metabolic pathways. It is conceivable that the plant metabolic responses to both aggressors may be conducive to virus acquisition. To address this question, we analyze the accumulation of the phloem-limited polerovirus Turnip yellows virus (TuYV), which is strictly transmitted by aphids, and aphid's life traits in six Arabidopsis thaliana mutants (xth33, ss3-2, nata1, myc234, quad, atr1D, and pad4-1). We observed that mutations affecting the carbohydrate metabolism, the synthesis of a non-protein amino acid and the glucosinolate pathway had an effect on TuYV accumulation. However, the virus titer did not correlate with the virus transmission efficiency. Some mutations in A. thaliana affect the aphid feeding behavior but often only in infected plants. The duration of the phloem sap ingestion phase, together with the time preceding the first sap ingestion, affect the virus transmission rate more than the virus titer did. Our results also show that the aphids reared on infected mutant plants had a reduced biomass regardless of the mutation and the duration of the sap ingestion phase.

Next generation sequencing and proteomics in plant virology: how is Colombia doing?

Crop production and trade are two of the most economically important activities in Colombia, and viral diseases cause a high negative impact to agricultural sector. Therefore, the detection, diagnosis, control, and management of viral diseases are crucial. Currently, Next-Generation Sequencing (NGS) and ‘Omic’ technologies constitute a right-hand tool for the discovery of novel viruses and for studying virus-plant interactions. This knowledge allows the development of new viral diagnostic methods and the discovery of key components of infectious processes, which could be used to generate plants resistant to viral infections. Globally, crop sciences are advancing in this direction. In this review, advancements in ‘omic’ technologies and their different applications in plant virology in Colombia are discussed. In addition, bioinformatics pipelines and resources for omics data analyses are presented. Due to their decreasing prices, NGS technologies are becoming an affordable and promising means to explore many phytopathologies affecting a wide variety of Colombian crops so as to improve their trade potential.

Evolution and ecology of plant viruses

The discovery of the first non-cellular infectious agent, later determined to be tobacco mosaic virus, paved the way for the field of virology. In the ensuing decades, research focused on discovering and eliminating viral threats to plant and animal health. However, recent conceptual and methodological revolutions have made it clear that viruses are not merely agents of destruction but essential components of global ecosystems. As plants make up over 80% of the biomass on Earth, plant viruses likely have a larger impact on ecosystem stability and function than viruses of other kingdoms. Besides preventing overgrowth of genetically homogeneous plant populations such as crop plants, some plant viruses might also promote the adaptation of their hosts to changing environments. However, estimates of the extent and frequencies of such mutualistic interactions remain controversial. In this Review, we focus on the origins of plant viruses and the evolution of interactions between these viruses and both their hosts and transmission vectors. We also identify currently unknown aspects of plant virus ecology and evolution that are of practical importance and that should be resolvable in the near future through viral metagenomics.

The Respiratory Burst Oxidase Homolog D (RbohD) Cell and Tissue Distribution in Potato-Potato Virus Y (PVYNTN) Hypersensitive and Susceptible Reactions

The respiratory burst oxidase homolog D (RbohD) acts as a central driving force of reactive oxygen species signaling in plant cells by integrating many different signal transduction pathways in plants, including incompatible interactions with pathogens. This study demonstrated the localization and distribution of RbohD in two types of potato–potato virus Y (PVY) interactions: Compatible and incompatible (resistant). The results indicated a statistically significant induction of the RbohD antigen signal in both interaction types. In the hypersensitive response (resistant reaction) of potato with a high level of resistance to the potato tuber necrotic strain of PVY (PVY^NTN), RbohD localization followed by hydrogen peroxide (H₂O₂) detection was concentrated in the apoplast. In contrast, in the hypersensitive response of potato with a low resistance level to PVY^NTN, the distribution of RbohD was concentrated more in the plant cell organelles than in the apoplast, resulting in the virus particles being present outside the inoculation area. Moreover, when compared to mock-inoculated plants and to the hypersensitive response, the PVY^NTN-compatible potato interaction triggered high induction in the RbohD distribution, which was associated with necrotization. Our findings indicated that RbohD and hydrogen peroxide deposition was associated with the hypersensitive response, and both were detected in the vascular tissues and chloroplasts. These results suggest that the RbohD distribution is actively dependent on different types of PVY ^NTN-potato plant interactions. Additionally, the RbohD may be involved in the PVY^NTN tissue limitation during the hypersensitive response, and it could be an active component of the systemic signal transduction in the susceptible host reaction.

Spatiotemporal Changes in Xylan-1/Xyloglucan and Xyloglucan Xyloglucosyl Transferase (XTH-Xet5) as a Step-In of Ultrastructural Cell Wall Remodelling in Potato⁻Potato Virus Y (PVYNTN) Hypersensitive and Susceptible Reaction

One type of monitoring system in a plant cell is the cell wall, which intensively changes its structure during interaction with pathogen-stress factors. The wall plays a role as a dynamic and controlled structure, although it is not fully understood how relevant these modifications are to the molecular mechanisms during plant–virus interactions. In this work we localise the non-cellulosic polysaccharides such as xyloglucan, xylan (xylan-1) and xyloglucosyl transferase (XTH-Xet5), the enzyme that participates in the metabolism of xyloglucan. This provided us with information about the in situ distribution of the components of the hemicellulotic cell wall matrix in hypersensitive and susceptible potato–PVY^NTN interactions. The loosening of the cell wall was accompanied by an increase in xylan depositions during susceptible interactions, whereas, during the hypersensitive response, when the cell wall was reinforced, the xylan content decreased. Moreover, the PVY inoculation significantly redirected XTH-Xet5 depositions, regardless of types of interactions, compared to mock-inoculated tissues. Furthermore, the immunogold localisation clearly revealed the domination of Xet5 in the cell wall and in vesicles in the susceptible host. In contrast, in the resistant host increased levels of Xet5 were observed in cytoplasm, in the cell wall and in the trans-Golgi network. These findings show that the hypersensitive reaction activated XTH-Xet5 in the areas of xyloglucan endo-transglycosylase (XET) synthesis, which was then actively transported to cytoplasm, cell wall and to vacuoles. Our results provide novel insight into cell wall reorganisation during PVY^NTN infection as a response to biotic stress factors. These novel findings help us to understand the mechanisms of defence responses that are incorporated into the cell wall signalling network.

Viral Fitness Correlates with the Magnitude and Direction of the Perturbation Induced in the Host's Transcriptome: The Tobacco Etch Potyvirus-Tobacco Case Study

Determining the fitness of viral genotypes has become a standard practice in virology as it is essential to evaluate their evolutionary potential. Darwinian fitness, defined as the advantage of a given genotype with respect to a reference one, is a complex property that captures, in a single figure, differences in performance at every stage of viral infection. To what extent does viral fitness result from specific molecular interactions with host factors and regulatory networks during infection? Can we identify host genes in functional classes whose expression depends on viral fitness? Here, we compared the transcriptomes of tobacco plants infected with seven genotypes of tobacco etch potyvirus that differ in fitness. We found that the larger the fitness differences among genotypes, the more dissimilar the transcriptomic profiles are. Consistently, two different mutations, one in the viral RNA polymerase and another in the viral suppressor of RNA silencing, resulted in significantly similar gene expression profiles. Moreover, we identified host genes whose expression showed a significant correlation, positive or negative, with the virus' fitness. Differentially expressed genes which were positively correlated with viral fitness activate hormone- and RNA silencing-mediated pathways of plant defense. In contrast, those that were negatively correlated with fitness affect metabolism, reducing growth, and development. Overall, these results reveal the high information content of viral fitness and suggest its potential use to predict differences in genomic profiles of infected hosts.

Different Metabolic Pathways Are Involved in Response of Saccharomyces cerevisiae to L-A and M Viruses

Competitive and naturally occurring yeast killer phenotype is governed by coinfection with dsRNA viruses. Long-term relationship between the host cell and viruses appear to be beneficial and co-adaptive; however, the impact of viral dsRNA on the host gene expression has barely been investigated. Here, we determined the transcriptomic profiles of the host Saccharomyces cerevisiae upon the loss of the M-2 dsRNA alone and the M-2 along with the L-A-lus dsRNAs. We provide a comprehensive study based on the high-throughput RNA-Seq data, Gene Ontology and the analysis of the interaction networks. We identified 486 genes differentially expressed after curing yeast cells of the M-2 dsRNA and 715 genes affected by the elimination of both M-2 and L-A-lus dsRNAs. We report that most of the transcriptional responses induced by viral dsRNAs are moderate. Differently expressed genes are related to ribosome biogenesis, mitochondrial functions, stress response, biosynthesis of lipids and amino acids. Our study also provided insight into the virus–host and virus–virus interplays.

Virus-host interactome: Putting the accent on how it changes

Viral infections are extremely complex processes that could only be well understood by precisely characterizing the interaction networks between the virus and the host components. In recent years, much effort has gone in this direction with the aim of unveiling the molecular basis of viral pathology. These networks are mostly formed by viral and host proteins, and are expected to be dynamic both with time and space (i.e., with the progression of infection, as well as with the virus and host genotypes; what we call plastodynamic). This largely overlooked spatio-temporal evolution urgently calls for a change both in the conceptual paradigms and experimental techniques used so far to characterize virus-host interactions. More generally, molecular plasticity and temporal dynamics are unavoidable components of the mechanisms that underlie any complex disease; components whose understanding will eventually enhance our ability to modulate those networks with the aim of improving disease treatments.

Effect of Host Species on Topography of the Fitness Landscape for a Plant RNA Virus

Adaptive fitness landscapes are a fundamental concept in evolutionary biology that relate the genotypes of individuals to their fitness. In the end, the evolutionary fate of evolving populations depends on the topography of the landscape, that is, the numbers of accessible mutational pathways and possible fitness peaks (i.e., adaptive solutions). For a long time, fitness landscapes were only theoretical constructions due to a lack of precise information on the mapping between genotypes and phenotypes. In recent years, however, efforts have been devoted to characterizing the properties of empirical fitness landscapes for individual proteins or for microbes adapting to artificial environments. In a previous study, we characterized the properties of the empirical fitness landscape defined by the first five mutations fixed during adaptation of tobacco etch potyvirus (TEV) to a new experimental host, Arabidopsis thaliana Here we evaluate the topography of this landscape in the ancestral host Nicotiana tabacum By comparing the topographies of the landscapes for the two hosts, we found that some features remained similar, such as the existence of fitness holes and the prevalence of epistasis, including cases of sign and reciprocal sign epistasis that created rugged, uncorrelated, and highly random topographies. However, we also observed significant differences in the fine-grained details between the two landscapes due to changes in the fitness and epistatic interactions of some genotypes. Our results support the idea that not only fitness tradeoffs between hosts but also topographical incongruences among fitness landscapes in alternative hosts may contribute to virus specialization. IMPORTANCE Despite its importance for understanding virus evolutionary dynamics, very little is known about the topography of virus adaptive fitness landscapes, and even less is known about the effects that different host species and environmental conditions may have on this topography. To bridge this gap, we evaluated the topography of a small fitness landscape formed by all genotypes that result from every possible combination of the first five mutations fixed during adaptation of TEV to the novel host A. thaliana To assess the effect that host species may have on this topography, we evaluated the fitness of every genotype in both the ancestral and novel hosts. We found that both landscapes share some macroscopic properties, such as the existence of holes and being highly rugged and uncorrelated, yet they differ in microscopic details due to changes in the magnitude and sign of fitness and epistatic effects.

Local adaptation of plant viruses: lessons from experimental evolution

For multihost pathogens, adaptation to multiple hosts has important implications for both applied and basic research. At the applied level, it is one of the main factors determining the probability and severity of emerging disease outbreaks. At the basic level, it is thought to be a key mechanism for the maintenance of genetic diversity both in host and pathogen species. In recent years, a number of evolution experiments have assessed the fate of plant virus populations replicating within and adapting to one single or to multiple hosts species. A first group of these experiments tackled the existence of trade-offs in fitness and virulence for viruses evolving either within a single hosts species or alternating between two different host species. A second set of experiments explored the role of genetic variability in susceptibility and resistance to infection among individuals from the same host species in the extent of virus local adaptation and of virulence. In general, when a single host species or genotype is available, these experiments show that local adaptation takes place, often but not always associated with a fitness trade-off. However, alternating between different host species or infecting resistant host genotypes may select for generalist viruses that experience no fitness cost. Therefore, the expected cost of generalism, arising from antagonistic pleiotropy and other genetic mechanisms generating fitness trade-offs between hosts, could not be generalized and strongly depend on the characteristics of each particular pathosystem. At the genomic level, these studies show pervasive convergent molecular evolution, suggesting that the number of accessible molecular pathways leading to adaptation to novel hosts is limited.