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

Physiological variations in hypovirus-infected wild and model long-term laboratory strains of Cryphonectria parasitica

Introduction

Forest ecosystems are highly threatened by the simultaneous effects of climate change and invasive pathogens. Chestnut blight, caused by the invasive phytopathogenic fungus Cryphonectria parasitica, has caused severe damage to European chestnut groves and catastrophic dieback of American chestnut in North America. Within Europe, the impacts of the fungus are widely mitigated through biological control that utilizes the RNA mycovirus: Cryphonectria hypovirus 1 (CHV1). Viral infections, similarly to abiotic factors, can cause oxidative stress in their hosts leading to physiological attrition through stimulating ROS (reactive oxygen species) and NOx production.

Methods

To fully understand the interactions leading to the biocontrol of chestnut blight, it is vital to determine oxidative stress damage arising during CHV1 infection, especially considering that other abiotic factors, like long-term cultivation of model fungal strains, can also impact oxidative stress. Our study compared CHV1-infected C. parasitica isolates from two Croatian wild populations with CHV1-infected model strains (EP713, Euro7 and CR23) that have experienced long-term laboratory cultivation.

Results and Discussion

We determined the level of oxidative stress in the samples by measuring stress enzymes' activity and oxidative stress biomarkers. Furthermore, for the wild populations, we studied the activity of fungal laccases, expression of the laccase gene lac1, and a possible effect of CHV1 intra-host diversity on the observed biochemical responses. Relative to the wild isolates, the long-term model strains had lower enzymatic activities of superoxide dismutase (SOD) and glutathione S-transferase (GST), and higher content of malondialdehyde (MDA) and total non-protein thiols. This indicated generally higher oxidative stress, likely arising from their decades-long history of subculturing and freeze-thaw cycles. When comparing the two wild populations, differences between them in stress resilience and levels of oxidative stress were also observed, as evident from the different MDA content. The intra-host genetic diversity of the CHV1 had no discernible effect on the stress levels of the virus-infected fungal cultures. Our research indicated that an important determinant modulating both lac1 expression and laccase enzyme activity is intrinsic to the fungus itself, possibly related to the vc type of the fungus, i.e., vegetative incompatibility genotype.

Long-read sequencing reveals the evolutionary drivers of intra-host diversity across natural RNA mycovirus infections

Intra-host dynamics are a core component of virus evolution but most intra-host data come from a narrow range of hosts or experimental infections. Gaining broader information on the intra-host diversity and dynamics of naturally occurring virus infections is essential to our understanding of evolution across the virosphere. Here we used PacBio long-read HiFi sequencing to characterize the intra-host populations of natural infections of the RNA mycovirus Cryphonectria hypovirus 1 (CHV1). CHV1 is a biocontrol agent for the chestnut blight fungus (Cryphonectria parasitica), which co-invaded Europe alongside the fungus. We characterized the mutational and haplotypic intra-host virus diversity of thirty-eight natural CHV1 infections spread across four locations in Croatia and Switzerland. Intra-host CHV1 diversity values were shaped by purifying selection and accumulation of mutations over time as well as epistatic interactions within the host genome at defense loci. Geographical landscape features impacted CHV1 inter-host relationships through restricting dispersal and causing founder effects. Interestingly, a small number of intra-host viral haplotypes showed high sequence similarity across large geographical distances unlikely to be linked by dispersal.

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.

Transmission modes affect the population structure of potato virus Y in potato

Transmission is a crucial part of a viral life cycle and transmission mode can have an important impact on virus biology. It was demonstrated that transmission mode can influence the virulence and evolution of a virus; however, few empirical data are available to describe the direct underlying changes in virus population structure dynamics within the host. Potato virus Y (PVY) is an RNA virus and one of the most damaging pathogens of potato. It comprises several genetically variable strains that are transmitted between plants via different transmission modes. To investigate how transmission modes affect the within-plant viral population structure, we have used a deep sequencing approach to examine the changes in the genetic structure of populations (in leaves and tubers) of three PVY strains after successive passages by horizontal (aphid and mechanical) and vertical (via tubers) transmission modes. Nucleotide diversities of viral populations were significantly influenced by transmission modes; lineages transmitted by aphids were the least diverse, whereas lineages transmitted by tubers were the most diverse. Differences in nucleotide diversities of viral populations between leaves and tubers were transmission mode-dependent, with higher diversities in tubers than in leaves for aphid and mechanically transmitted lineages. Furthermore, aphid and tuber transmissions were shown to impose stronger genetic bottlenecks than mechanical transmission. To better understand the structure of virus populations within the host, transmission mode, movement of the virus within the host, and the number of replication cycles after transmission event need to be considered. Collectively, our results suggest a significant impact of virus transmission modes on the within-plant diversity of virus populations and provide quantitative fundamental data for understanding how transmission can shape virus diversity in the natural ecosystems, where different transmission modes are expected to affect virus population structure and consequently its evolution.

The Potyviruses: An Evolutionary Synthesis Is Emerging

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

Pepper veinal mottle virus in Japan is closely related to isolates from other Asian countries, but more distantly to most of those from Africa

Pepper veinal mottle virus (PVMV) is known to infect chilli pepper and belongs to the Chilli veinal mottle virus phylogroup of potyviruses. PVMV has recently appeared in Japan. In this study, we report six complete genomic sequences of PVMV isolates from chilli pepper (i.e. Capsicum annuum) in Okinawa Islands in Japan, and we determined the evolutionary relationships between Japanese isolates and the isolates reported earlier from African and Asian countries. Complete genomic sequences of the six Japanese PVMV isolates were 9760 nucleotides in length, excluding the nucleotide primer sequences used for amplifying 5' end of the genomes. The major findings of this study are as follows: (1) all the Japanese isolates of PVMV have similar biological and molecular characteristics, indicating the presence of only one population in Japan; (2) there are at least three major phylogenetic groups of PVMV worldwide; (3) PVMV probably originated in East Africa; and (4) all the Asian isolates are closely related to the Ghanaian isolate.

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.

Mature Luffa Leaves (Luffa cylindrica L.) as a Tool for Gene Expression Analysis by Agroinfiltration

We exploited the potential of cucurbits for ectopic gene expression. Agroinfiltration is a simple and commonly used method to obtain transient expression of foreign genes in plants. In contrast to in vitro transformation techniques, agroinfiltration can be used for genetic modification of mature plant tissues. Although the cucurbits are commonly used as model plants for molecular biology and biotechnology studies, to date there are no literature sources on the possibility of transient gene expression in mature cucurbit tissues. Our research has shown that mature leaves of Luffa cylindrica L. (luffa), in contrast to other cucurbit species, can be successfully transiently transformed with Agrobacterium tumefaciens. We efficiently transformed luffa leaves with a reporter gene encoding β-glucuronidase (GUS). The GUS activity in transiently transformed leaf tissues was detected within 24 h after the infiltration with bacteria. Additionally, we have shown that the activity of a transiently expressed the GUS gene can be monitored directly in the EDTA-exudates collected from the cut petioles of the agroinfiltrated leaves. The results suggest that luffa leaves can be useful as a plant expression system for studies of physiological and biochemical processes in cucurbits.

Deep sequencing for discovery and evolutionary analysis of plant viruses

The advent of next generation sequencing (NGS), or deep sequencing, has allowed great advances to be made in discovery, diagnostics, and evolutionary studies in plant viruses. Various methods have been used for enrichment for virus-specific nucleic acids, each of which have some drawbacks. Many novel viruses have been discovered in plants by NGS technologies, and there is a good deal of promise for more comprehensive studies in virus evolution. However, each aspect of using NGS has its caveats that need to be considered, and there is still a need for better tools of analysis, as well as method for validation of sequence variation.

Sequence polymorphism in an insect RNA virus field population: A snapshot from a single point in space and time reveals stochastic differences among and within individual hosts

Population structure of Homalodisca coagulata Virus-1 (HoCV-1) among and within field-collected insects sampled from a single point in space and time was examined. Polymorphism in complete consensus sequences among single-insect isolates was dominated by synonymous substitutions. The mutant spectrum of the C2 helicase region within each single-insect isolate was unique and dominated by nonsynonymous singletons. Bootstrapping was used to correct the within-isolate nonsynonymous:synonymous arithmetic ratio (N:S) for RT-PCR error, yielding an N:S value ~one log-unit greater than that of consensus sequences. Probability of all possible single-base substitutions for the C2 region predicted N:S values within 95% confidence limits of the corrected within-isolate N:S when the only constraint imposed was viral polymerase error bias for transitions over transversions. These results indicate that bottlenecks coupled with strong negative/purifying selection drive consensus sequences toward neutral sequence space, and that most polymorphism within single-insect isolates is composed of newly-minted mutations sampled prior to selection.

Next-Generation Sequencing and Genome Editing in Plant Virology

Next-generation sequencing (NGS) has been applied to plant virology since 2009. NGS provides highly efficient, rapid, low cost DNA, or RNA high-throughput sequencing of the genomes of plant viruses and viroids and of the specific small RNAs generated during the infection process. These small RNAs, which cover frequently the whole genome of the infectious agent, are 21-24 nt long and are known as vsRNAs for viruses and vd-sRNAs for viroids. NGS has been used in a number of studies in plant virology including, but not limited to, discovery of novel viruses and viroids as well as detection and identification of those pathogens already known, analysis of genome diversity and evolution, and study of pathogen epidemiology. The genome engineering editing method, clustered regularly interspaced short palindromic repeats (CRISPR)-Cas9 system has been successfully used recently to engineer resistance to DNA geminiviruses (family, Geminiviridae) by targeting different viral genome sequences in infected Nicotiana benthamiana or Arabidopsis plants. The DNA viruses targeted include tomato yellow leaf curl virus and merremia mosaic virus (begomovirus); beet curly top virus and beet severe curly top virus (curtovirus); and bean yellow dwarf virus (mastrevirus). The technique has also been used against the RNA viruses zucchini yellow mosaic virus, papaya ringspot virus and turnip mosaic virus (potyvirus) and cucumber vein yellowing virus (ipomovirus, family, Potyviridae) by targeting the translation initiation genes eIF4E in cucumber or Arabidopsis plants. From these recent advances of major importance, it is expected that NGS and CRISPR-Cas technologies will play a significant role in the very near future in advancing the field of plant virology and connecting it with other related fields of biology.

Triticum mosaic virus exhibits limited population variation yet shows evidence of parallel evolution after replicated serial passage in wheat

An infectious cDNA clone of Triticum mosaic virus (TriMV) (genus Poacevirus; family Potyviridae) was used to establish three independent lineages in wheat to examine intra-host population diversity levels within protein 1 (P1) and coat protein (CP) cistrons over time. Genetic variation was assessed at passages 9, 18 and 24 by single-strand conformation polymorphism, followed by nucleotide sequencing. The founding P1 region genotype was retained at high frequencies in most lineage/passage populations, while the founding CP genotype disappeared after passage 18 in two lineages. We found that rare TriMV genotypes were present only transiently and lineages followed independent evolutionary trajectories, suggesting that genetic drift dominates TriMV evolution. These results further suggest that experimental populations of TriMV exhibit lower mutant frequencies than that of Wheat streak mosaic virus (genus Tritimovirus; family Potyviridae) in wheat. Nevertheless, there was evidence for parallel evolution at a synonymous site in the TriMV CP cistron.

Transcriptional slippage in the positive-sense RNA virus family Potyviridae

The family Potyviridae encompasses ~30% of plant viruses and is responsible for significant economic losses worldwide. Recently, a small overlapping coding sequence, termed pipo, was found to be conserved in the genomes of all potyvirids. PIPO is expressed as part of a frameshift protein, P3N-PIPO, which is essential for virus cell-to-cell movement. However, the frameshift expression mechanism has hitherto remained unknown. Here, we demonstrate that transcriptional slippage, specific to the viral RNA polymerase, results in a population of transcripts with an additional "A" inserted within a highly conserved GAAAAAA sequence, thus enabling expression of P3N-PIPO. The slippage efficiency is ~2% in Turnip mosaic virus and slippage is inhibited by mutations in the GAAAAAA sequence. While utilization of transcriptional slippage is well known in negative-sense RNA viruses such as Ebola, mumps and measles, to our knowledge this is the first report of its widespread utilization for gene expression in positive-sense RNA viruses.

Temporal dynamics of intrahost molecular evolution for a plant RNA virus

Populations of plant RNA viruses are highly polymorphic in infected plants, which may allow rapid within-host evolution. To understand tobacco etch potyvirus (TEV) evolution, longitudinal samples from experimentally evolved populations in the natural host tobacco and from the alternative host pepper were phenotypically characterized and genetically analyzed. Temporal and compartmental variabilities of TEV populations were quantified using high throughput Illumina sequencing and population genetic approaches. Of the two viral phenotypic traits measured, virulence increased in the novel host but decreased in the original one, and viral load decreased in both hosts, though to a lesser extent in the novel one. Dynamics of population genetic diversity were also markedly different among hosts. Population heterozygosity increased in the ancestral host, with a dominance of synonymous mutations fixed, whereas it did not change or even decreased in the new host, with an excess of nonsynonymous mutations. All together, these observations suggest that directional selection is the dominant evolutionary force in TEV populations evolving in a novel host whereas either diversifying selection or random genetic drift may play a fundamental role in the natural host. To better understand these evolutionary dynamics, we developed a computer simulation model that incorporates the effects of mutation, selection, and drift. Upon parameterization with empirical data from previous studies, model predictions matched the observed patterns, thus reinforcing our idea that the empirical patterns of mutation accumulation represent adaptive evolution.