Analysis of the virome associated to grapevine downy mildew lesions reveals new mycovirus lineages.

The obligatory biotrophic oomycetes Plasmopara viticola is the causal agent of downy mildew, a destructive disease of grapevine worldwide. So far, chemical fungicides are widely employed to limit this pathogen, but their adverse effects are stimulating the quest for environmentally friendly alternative approaches. Here, we report on the search for mycoviruses that might be later developed as biocontrol agents for this pathogen. Symptomatic leaves were collected from various regions in Spain and Italy and mycelia associated to leaf lesions was harvested. Total RNA extractions were depleted of rRNA and metatranscriptomes were generated using a high-throughput sequencing approach. The virome associated to leaf lesions was then characterized through a bioinformatic pipeline relying on blast searches against current viral databases. Here, we present an inventory of 283 new RNA viruses: 222 positive strand RNA viruses, 29 negative strand RNA viruses, 27 double-stranded RNA viruses and 5 ORFan virus RdRP segments, which could not be reliably assigned to any existing group in the Riboviria. In addition to ORFan viruses, we found other surprising new evolutionary trajectories in this wide inventory of viruses. The most represented viruses in our collection are those in phylum Lenarviricota, and, among them, a group of mycovirus segments distantly related to narnaviruses, but characterized by a polymerase palm domain lacking subdomain C, with the putative GDD catalytic triad. We also provided evidence of a strict association between two RNA segments that form a new mycovirus clade of positive strand RNA in the phylum Kitrinoviricota, order Martellivirales. In the phylum Negarnaviricota, we report for the first time in the order Mononegavirales a clade of viruses that is ambisense, a feature that so far was present only in the order Bunyavirales. Furthermore, in the same phylum we detected the widespread occurrence and abundant accumulation in our libraries of a distinct mycovirus clade distantly related to the Muvirales and Goujanvirales orders, which so far include only viruses infecting invertebrates. Possible new oomycetes-specific virus clades are also described in the phylum Duplornaviricota. These data greatly expand the evolutionary history of mycoviruses adding new layers of diversity to the realm Riboviria.

Landscape heterogeneity shapes host-parasite interactions and results in apparent plant-virus codivergence.

Knowledge on how landscape heterogeneity shapes host-parasite interactions is central to understand the emergence, dynamics and evolution of infectious diseases. However, this is an underexplored subject, particularly for plant-virus systems. Here, we analyse how landscape heterogeneity influences the prevalence, spatial genetic structure, and temporal dynamics of Pepper golden mosaic and Pepper huasteco yellow vein begomoviruses infecting populations of the wild pepper Capsicum annuum glabriusculum (chiltepin) in Mexico. Environmental heterogeneity occurred at different nested spatial scales (host populations within biogeographical provinces), with levels of human management varying among host population within a province. Results indicate that landscape heterogeneity affects the epidemiology and genetic structure of chiltepin-infecting begomoviruses in a scale-specific manner, probably related to conditions favouring the viruses' whitefly vector and its dispersion. Increased levels of human management of the host populations were associated with higher virus prevalence and erased the spatial genetic structure of the virus populations. Also, environmental heterogeneity similarly shaped the spatial genetic structures of host and viruses. This resulted in the congruence between host and virus phylogenies, which does not seem to be due to host-virus co-evolution. Thus, results provide evidence of the key role of landscape heterogeneity in determining plant-virus interactions.

Evolutionary analysis of genetic variation observed in citrus tristeza virus (CTV) after host passage.

We have studied the genetic variability in two genes (p18 and p20) from two groups of Citrus tristeza virus (CTV) isolates. One group (isolates T385, T317, T318, and T305) was derived from a Spanish source by successive host passages while the other (isolates T388 and T390) was obtained after aphid transmission from a Japanese source. A total of 274 sequences were obtained for gene p18 and 451 for p20. In the corresponding phylogenetic trees, sequences derived from the severe isolates (T318, T305, and T388) clustered together and separately from those derived from mild or moderate isolates (T385, T317, and T390), regardless of their geographic origin. Hierarchical analyses of molecular variance showed that up to 53% of the total genetic variability in p18 and up to 87% of the variation in p20 could be explained by differences in the pathogenicity features of the isolates. Neutrality tests revealed that different selection forces had been acting between isolates and between genes, with purifying selection being suggested for p18 from isolates T385 and T390 and for p20 from isolates T385, T317, and T388, and balancing selection for p18 from isolates T318, T305, and T388 and for p20 from isolates T318 and T390. Furthermore, several models of codon selection were observed, with purifying selection being the most notable one, compatible with low effective population size of the virus populations resulting from transmission bottlenecks. We found no evidence of recombination playing a significant role during p18 and p20 evolution in these isolates. These results suggest that hosts can be an important evolutionary factor for CTV isolates.

Genetic variation of Citrus tristeza virus isolates from California and Spain: evidence for mixed infections and recombination.

We examined the population structure and genetic variation of four genomic regions within and between 30 Citrus tristeza virus (CTV) isolates from Spain and California. Our analyses showed that most isolates contained a population of sequence variants, with one being predominant. Four isolates showed two major sequence variants in some genomic regions. The two major variants of three of these isolates showed very low nucleotide identity to each other but were very similar to those of other isolates, suggesting the possibility of mixed infections with two divergent isolates. Incongruencies of phylogenetic relationships in the different genomic regions and statistical analyses suggested that the genomes of some CTV sequence variants originated by recombination events between diverged sequence variants. No correlation was observed between geographic origin and nucleotide distance, and thus from a genetic view, the Spanish and Californian isolates analyzed here could be considered members of the same population.

Characterization of the cis-acting elements controlling subgenomic mRNAs of citrus tristeza virus: production of positive- and negative-stranded 3'-terminal and positive-stranded 5'-terminal RNAs.

Citrus tristeza virus (CTV), a member of the Closteroviridae, has an approximately 20-kb positive-sense RNA genome with two 5' ORFs translated from the genomic RNA and 10 3' genes expressed via nine or ten 3'-terminal subgenomic (sg) RNAs. The expression of the 3' genes appears to have properties intermediate between the smaller viruses of the "alphavirus supergroup" and the larger viruses of the Coronaviridae. The sgRNAs are contiguous with the genome, without a common 5' leader, and are associated with large amounts of complementary sgRNAs. Production of the different sgRNAs is regulated temporally and quantitatively, with the highly expressed genes having noncoding regions (NCR) 5' of the ORFs. The cis-acting elements that control the highly expressed major coat protein (CP) gene and the intermediately expressed minor coat protein (CPm) gene were mapped and compared. Mutational analysis showed that the CP sgRNA controller element mapped within nts -47 to -5 upstream of the transcription start site, entirely within the NCR, while the CPm control region mapped within a 57 nt sequence within the upstream ORF. Although both regions were predicted to fold into two stem-loop structures, mutagenesis suggested that primary structure might be more important than the secondary structure. Because each controller element produced large amounts of 3'-terminal positive- and negative-stranded sgRNAs, we could not differentiate whether the cis-acting element functioned as a promoter or terminator, or both. Reversal of the control element unexpectedly produced large amounts of a negative-stranded sgRNA apparently by termination of negative-stranded genomic RNA synthesis. Further examination of controller elements in their native orientation showed normal production of abundant amounts of positive-stranded sgRNAs extending to near the 5'-terminus, corresponding to termination at each controller element. Thus, each controller element produced three sgRNAs, a 5'-terminal positive strand and both positive- and negative-stranded 3'-terminal RNAs. Therefore, theoretically CTV could produce 30-33 species of RNAs in infected cells.

Polymorphism of the 5' terminal region of Citrus tristeza virus (CTV) RNA: incidence of three sequence types in isolates of different origin and pathogenicity.

Sequences of the 5' terminal region of the genomic RNA from eight isolates of Citrus tristeza virus (CTV) were previously classified into three types (I, II and III), with intragroup sequence identity higher than 88% and intergroup sequence identity as low as 44%. Sequencing of an additional 58 cDNA clones from 15 virus isolates showed that all sequences could be unequivocally assigned to one of the three types previously established. The relative frequency of each sequence type was assessed in 57 CTV isolates of different geographic origin and pathogenic characteristics by RT-PCR with sets of type-specific primers using CTV dsRNA as template. None of the isolates yielded amplification of the type I or II sequences alone, but in 19 of them type III sequences were the only amplification product detected. Within isolates containing more than one sequence type, eight had type II and III sequences, 11 had type I and III sequences, and 19 had sequences of the three types. Isolates containing only type III sequences caused only mild to moderate symptoms in Mexican lime, an indicator species for most CTV isolates, whereas isolates causing stem pitting in sweet orange an/or grapefruit, generally contained sequences type II. None of the sequence types could be traced to a precise geographic area, as all types were detected in isolates from at least nine of the 12 countries from which samples were taken.

Amplification of Citrus tristeza virus from a cDNA clone and infection of citrus trees.

Isolates of the Closterovirus, Citrus tristeza virus (CTV), are populations of disparate genotypes and defective RNAs developed during long periods of vegetative propagation of citrus trees. Because it has not been possible to obtain pure cultures of the virus, it is not known what components of the population are primarily responsible for induction of diseases. We previously developed an infectious cDNA clone from which in vitro-produced RNA transcripts could infect protoplasts (Satyanarayana et al., 1999, Proc. Natl. Acad. Sci. USA 96, 7433-7438). However, neither the RNA transcripts nor virions from transcript-infected protoplasts were competent for infection of citrus trees. Using a green fluorescent protein-marked virus as inoculum, we found that the approximately 20-kb RNA from virions or transcripts of cDNA infected only a small percentage of protoplasts ( approximately 0.01%), but virions could infect more than 80% of the protoplasts. Based on this information, we amplified the virus from the cDNA clone (recombinant virus) by successive passages in protoplasts using virions in crude sap as inoculum. By the third to seventh passages in protoplasts, maximal amounts of recombinant progeny virus were produced, which were used for inoculation of small citrus trees by slashing stems in the presence of virion preparations. A relatively high percentage of plants became infected with the recombinant virus from protoplasts, resulting in the first defined pure culture of CTV in plants. The comparative biology of the pure culture of recombinant CTV with that of the parental population in planta demonstrated that the recombinant virus retained through all of the recombinant DNA manipulations the normal functions of replication, movement, and aphid transmissibility, and had a symptom phenotype indistinguishable from that of the parental population. Additionally, fulfilling Koch's postulates of the first pure culture of CTV in plants suggested that the major genotype of the CTV T36 population is the primary determinant of the symptom phenotype. We could distinguish no biological contributions resulting from the minor genotypes and defective RNAs of the parental population.

Closterovirus encoded HSP70 homolog and p61 in addition to both coat proteins function in efficient virion assembly.

Assembly of the viral genome into virions is a critical process of the virus life cycle often defining the ability of the virus to move within the plant and to be transmitted horizontally to other plants. Closteroviridae virions are polar helical rods assembled primarily by a major coat protein, but with a related minor coat protein at one end. The Closteroviridae is the only virus family that encodes a protein with similarity to cellular chaperones, a 70-kDa heat-shock protein homolog (HSP70h). We examined the involvement of gene products of Citrus tristeza virus (CTV) in virion formation and found that the chaperone-like protein plus the p61 and both coat proteins were required for efficient virion assembly. Competency of virion assembly of different CTV mutants was assayed by their ability to be serially passaged in Nicotiana benthamiana protoplasts using crude sap as inoculum, and complete and partial virus particles were analyzed by serologically specific electron microscopy. Deletion mutagenesis revealed that p33, p6, p18, p13, p20, and p23 genes were not needed for virion formation. However, deletion of either minor- or major-coat protein resulted in formation of short particles which failed to be serially transferred in protoplasts, suggesting that both coat proteins are required for efficient virion assembly. Deletion or mutation of HSP70h and/or p61 dramatically reduced passage and formation of full-length virions. Frameshift mutations suggested that the HSP70h and p61 proteins, not the RNA sequences, were needed for virion assembly. Substitution of the key amino acid residues in the ATPase domain of HSP70h, Asp(7) to Lys or Glu(180) to Arg, reduced assembly, suggesting that the chaperone-like ATPase activity is involved in assembly. Both HSP70h and p61 proteins appeared to contribute equally to assembly, consistent with coordinate functions of these proteins in closterovirus virion formation. The requirement of two accessory proteins in addition to both coat proteins for efficient assembly is uniquely complex for helical virions.

The fitness of citrus tristeza virus defective RNAs is affected by the lengths of their 5'- and 3'-termini and by the coding capacity.

Populations of the Closterovirus Citrus tristeza virus (CTV) generally contain defective RNAs (dRNAs) that vary in size, abundance, and sequence. The variation in abundance of the different dRNAs in a population suggests selection for those of higher fitness. To examine factors affecting fitness of dRNAs, we investigated a series of in vitro constructed dRNAs for their ability to be amplified in protoplasts by an efficiently replicated CTV deletion mutant. The minimal sequences required for accumulation of the dRNAs were within the genomic 5' proximal approximately 1 kb and the 3' 270 nucleotides. However, other factors were involved, because a dRNA with only the minimal sequences failed to be replicated. Rescue of a nonviable dRNA by insertion of nonviral sequences between the termini suggested that "spacing" between terminal cis-acting signals influenced fitness. A continuous open reading frame (ORF) through most of the sequences derived from the 5' of the genome was a requirement for dRNA amplification. In general, insertions, deletions, or nucleotide substitutions were tolerated in the dRNAs as long as an ORF was retained, whereas dRNAs with mutations that prematurely terminated the ORF were not viable. To discriminate between a requirement for an essential protein and ribosomal travel, perhaps to present replication signals to the replicase complex, mutations were made to modify the potential protein but still maintain an ORF. Deletions, insertions of nonviral sequences, or switching of reading frames that altered the amino acid sequence of the protein, except the N-terminal 161 amino acids, did not destroy the fitness of the dRNAs. Yet termination of the ORF in the middle of nonviral sequences did destroy the ability of the dRNAs to be amplified. These results suggest that even though a continuous ORF was needed for fitness, its protein product did not affect the amplification of the dRNAs.

A new procedure to differentiate citrus tristeza virus isolates by hybridisation with digoxigenin-labelled cDNA probes.

A non-isotopic hybridisation procedure was developed to differentiate isolates of citrus tristeza virus (CTV) using digoxigenin (DIG)-labelled cDNA probes and different kinds of target RNA. Hybridisation of DIG-probes with purified double-stranded RNA (dsRNA) or concentrated total RNA extracts spotted on nylon membranes allowed detection of CTV nucleic acid equivalent to as little as 0.1-1 mg infected tissue, when the reaction was developed with a chemiluminiscent substrate. This sensitivity was similar to or slightly better than that obtained by hybridisation with a 32P-labelled probe. CTV was also detectable by hybridisation of DIG-probes with tissue prints from freshly cut young citrus shoots. Hybridisation of tissue prints with DIG-probes under stringent conditions (60 degrees C and 50% formamide) could differentiate CTV isolates in citrus, whether grown in the greenhouse or in the field. This rapid and sensitive procedure can easily be applied to many samples, even under field conditions, and opens the way to monitoring spatio-temporal movement of specific CTV strains (or groups of strains) in epidemiological studies.

The haplotype distribution of two genes of citrus tristeza virus is altered after host change or aphid transmission.

Genetic variability of citrus tristeza virus (CTV) was studied using the haplotypes detected by single-strand conformation polymorphism (SSCP) analysis of genes p18 and p20 in six virus populations of two origins. The Spanish group included a CTV isolate and subisolates obtained by graft-transmission to different host species. The other included two subisolates aphid-transmitted from a single Japanese isolate. The homozygosity observed for gene p20 was always significantly higher than that expected under neutral evolution, whereas only three populations showed high homozygosity for p18, suggesting stronger host constraints for p20 than for p18. Sequential transmissions of a Spanish isolate to new host species increased the difference between its population and that of the successive subisolates for gene p18, as estimated by the F statistic. Analysis of molecular variance indicated that variation between both groups of populations was not statistically significant, whereas variations between populations of the same group or within populations were significant for both genes studied. Our data indicate that selection affects the haplotype distribution and that adaptation to a new host can be as important or more as the geographical origin. Variation of the CTV populations after host change or aphid transmission may explain in part the wide biological variability observed among CTV isolates.

Molecular variability of the 5'- and 3'-terminal regions of citrus tristeza virus RNA.

ABSTRACT Isolates of citrus tristeza virus (CTV) differ widely in their biological properties. These properties may depend on the structure of viral RNA populations comprising the different isolates. As a first approach to study the molecular basis of the biological variability, we have compared the sequences of multiple cDNA clones of the two terminal regions of the RNA from different CTV isolates. The polymorphism of the 5' untranslated region (UTR) allowed the classification of the sequences into three groups, with intragroup sequence identity higher than 88% and intergroup sequence identity as low as 44%. The variability of an open reading frame (ORF) 1a segment adjacent to the 5' UTR supports the same grouping. Some CTV isolates contained sequences of more than one group. Most sequences from Spanish isolates belonged to group III, whereas a Japanese isolate was composed mostly of sequences of groups I and II. The mildest isolates contained only sequences of group III, whereas the most severe isolates also contained sequences of groups I, II, or both. The most stable secondary structure predicted for the 5' UTR was composed of two stem-loops and remained essentially unchanged as a result of compensatory mutations in the stems and accommodation of most of the variability in the loops. In contrast to the 5'-terminal region, the variability of the 3'-terminal region of CTV RNA was very much restricted, with nucleotide identity values higher than 90%. The presence of a conserved putative "zinc-finger" domain adjacent to a basic region in p23, the predicted product of ORF 11, suggests that this protein might act as a regulatory factor during virus replication.