Detection of single nucleotide polymorphisms in virus genomes assembled from high-throughput sequencing data: large-scale performance testing of sequence analysis strategies.

Recent developments in high-throughput sequencing (HTS) technologies and bioinformatics have drastically changed research in virology, especially for virus discovery. Indeed, proper monitoring of the viral population requires information on the different isolates circulating in the studied area. For this purpose, HTS has greatly facilitated the sequencing of new genomes of detected viruses and their comparison. However, bioinformatics analyses allowing reconstruction of genome sequences and detection of single nucleotide polymorphisms (SNPs) can potentially create bias and has not been widely addressed so far. Therefore, more knowledge is required on the limitations of predicting SNPs based on HTS-generated sequence samples. To address this issue, we compared the ability of 14 plant virology laboratories, each employing a different bioinformatics pipeline, to detect 21 variants of pepino mosaic virus (PepMV) in three samples through large-scale performance testing (PT) using three artificially designed datasets. To evaluate the impact of bioinformatics analyses, they were divided into three key steps: reads pre-processing, virus-isolate identification, and variant calling. Each step was evaluated independently through an original, PT design including discussion and validation between participants at each step. Overall, this work underlines key parameters influencing SNPs detection and proposes recommendations for reliable variant calling for plant viruses. The identification of the closest reference, mapping parameters and manual validation of the detection were recognized as the most impactful analysis steps for the success of the SNPs detections. Strategies to improve the prediction of SNPs are also discussed.

The Expanding Menagerie of Prunus-Infecting Luteoviruses.

Members of the genus Luteovirus are responsible for economically destructive plant diseases worldwide. Over the past few years, three luteoviruses infecting Prunus trees have been characterized. However, the biological properties, prevalence, and genetic diversity of those viruses have not yet been studied. High-throughput sequencing of samples of various wild, cultivated, and ornamental Prunus species enabled the identification of four novel species in the genus Luteovirus for which we obtained complete or nearly complete genomes. Additionally, we identified another new putative species recovered from Sequence Read Archive data. Furthermore, we conducted a survey on peach-infecting luteoviruses in eight European countries. Analyses of 350 leaf samples collected from germplasm, production orchards, and private gardens showed that peach-associated luteovirus (PaLV), nectarine stem pitting-associated virus (NSPaV), and a novel luteovirus, peach-associated luteovirus 2 (PaLV2), are present in all countries; the most prevalent virus was NSPaV, followed by PaLV. The genetic diversity of these viruses was also analyzed. Moreover, the biological indexing on GF305 peach indicator plants demonstrated that PaLV and PaLV2, like NSPaV, are transmitted by graft at relatively low rates. No clear viral symptoms have been observed in either graft-inoculated GF305 indicators or different peach tree varieties observed in an orchard. The data generated during this study provide a broader overview of the genetic diversity, geographical distribution, and prevalence of peach-infecting luteoviruses and suggest that these viruses are likely asymptomatic in peach under most circumstances.

Molecular and biological characterization of a new mulberry idaeovirus.

A virus-like disease with symptoms including mosaic structure, deformation, vein clearing and necrosis on the leaves and deformation, crumbling, and scab on the fruits was detected in black mulberry trees (Morus nigra L.) in Kayseri province of Turkey. A novel positive single-stranded RNA virus with a bipartite genome and the mulberry badnavirus 1 (MBV-1) were detected in the black mulberry trees by high throughput sequencing and bioinformatic analyses. The novel virus RNA1 (5,796/7 nt) encodes a polyprotein (1,808 aa, 204.31 kDa) with three conserved domains, [MTR (aa 294-705), Hel (aa 971-1,226) and RdRp (aa 1,348-1,788)], whereas RNA2 (2,243 nt) encodes two putative proteins, MP (374 aa, 40.98 kDa), and CP (272 aa, 30.59 kDa), separated by an intergenic region of 97 nt. The highest amino acids identities were 70, 57 and 70 % with raspberry bushy dwarf virus (RBDV) for ORF1, MP and CP genes, respectively. The genome organization and phylogenetic analyses suggested that the novel virus is likely a putative new member of the genus Idaeovirus and it has been tentatively named black mulberry idaeovirus (BMIV). Virus survey showed both the BMIV and MBV-1 are likely prevalent in the region. Seven complete (six Turkish and one Iranian) and 41 partial genome sequences of the BMIV isolates revealed moderate genetic diversity (0.033 ± 0.001 %, 0.020 ± 0.002 % and 0.016 ± 0.002 % for RNA1, RNA2, and partial genomes, respectively). Both the BMIV and MBV-1 were detected in all tested pollens (n = 24, 100 %), in seed-borne balck mulberry saplings (n = 96, 100 %).This situation clearly revealed the potential spread risk of both viruses in black mulberry plantations and the necessity of taking precautions.

Identification and molecular characterization of a novel foveavirus from Rubus spp. in Turkey.

A novel plant virus was identified by high-throughput sequencing analysis from a raspberry plant showing slight mottling symptom. The complete genome sequence of this virus is 8645 nucleotides long, including the 5' and 3' UTRs. Its genome contains five ORFs and is very close to members of the genus Foveavirus (Quinvirinae, Betaflexiviridae) in terms of genome organization, TGB presence and the sizes of the RdRp and CP proteins. The novel virus shares 33.5-51.3 % and 23.3-41.3 % nucleotide identity to other genera of the Betaflexifiviridae family based on polymerase (RdRp) and CP genes, respectively. Compared to other foveavirus species, the RdRp protein showed the highest sequence identity (45.3 %) to the RdRp of peach chlorotic mottle virus (PCMV) while the maximal sequence identity for the CP protein was 33.9 % with grapevine rupestris stem pitting-associated virus (GRSPaV). The low nucleotide and amino acid sequence identity with known foveaviruses indicated that it was a novel virus, for which the provisional name "rubus virus 1 (RuV1)" is proposed. The phylogenetic analysis supports the assignment of this virus as a new species of the genus Foveavirus. A survey of 537 Rubus spp. samples grown in six provinces of Turkey, including some symptomatic samples, showed a RuV1 prevalence of 2.2 %, confirming its presence in both raspberry and blackberry plants in a single province, although no obvious association between virus infection and specific symptoms was found.

Identification and Characterization of a Novel Robigovirus Species from Sweet Cherry in Turkey.

High throughput sequencing of total RNA isolated from symptomatic leaves of a sweet cherry tree (Prunus avium cv. 0900 Ziraat) from Turkey identified a new member of the genus Robigovirus designated cherry virus Turkey (CVTR). The presence of the virus was confirmed by electron microscopy and overlapping RT-PCR for sequencing its whole-genome. The virus has a ssRNA genome of 8464 nucleotides which encodes five open reading frames (ORFs) and comprises two non-coding regions, 5' UTR and 3' UTR of 97 and 296 nt, respectively. Compared to the five most closely related robigoviruses, RdRp, TGB1, TGB2, TGB3 and CP share amino acid identities ranging from 43-53%, 44-60%, 39-43%, 38-44% and 45-50%, respectively. Unlike the four cherry robigoviruses, CVTR lacks ORFs 2a and 5a. Its genome organization is therefore more similar to African oil palm ringspot virus (AOPRV). Using specific primers, the presence of CVTR was confirmed in 15 sweet cherries and two sour cherries out of 156 tested samples collected from three regions in Turkey. Among them, five samples were showing slight chlorotic symptoms on the leaves. It seems that CVTR infects cherry trees with or without eliciting obvious symptoms, but these data should be confirmed by bioassays in woody and possible herbaceous hosts in future studies.

Genetic variation and possible mechanisms driving the evolution of worldwide fig mosaic virus isolates.

Fig mosaic virus (FMV) is a multipartite negative-sense RNA virus infecting fig trees worldwide. FMV is transmitted by vegetative propagation and grafting of plant materials, and by the eriophyid mite Aceria ficus. In this work, the genetic variation and evolutionary mechanisms shaping FMV populations were characterized. Nucleotide sequences from four genomic regions (each within the genomic RNAs 1, 2, 3, and 4) from FMV isolates from different countries were determined and analyzed. FMV genetic variation was low, as is seen for many other plant viruses. Phylogenetic analysis showed some geographically distant FMV isolates which clustered together, suggesting long-distance migration. The extent of migration was limited, although varied, between countries, such that FMV populations of different countries were genetically differentiated. Analysis using several recombination algorithms suggests that genomes of some FMV isolates originated by reassortment of genomic RNAs from different genetically similar isolates. Comparison between nonsynonymous and synonymous substitutions showed selection acting on some amino acids; however, most evolved neutrally. This and neutrality tests together with the limited gene flow suggest that genetic drift plays an important role in shaping FMV populations.

Endophytic bacterial community living in roots of healthy and 'Candidatus Phytoplasma mali'-infected apple (Malus domestica, Borkh.) trees.

'Candidatus Phytoplasma mali', the causal agent of apple proliferation (AP) disease, is a quarantine pathogen controlled by chemical treatments against insect vectors and eradication of diseased plants. In accordance with the European Community guidelines, novel strategies should be developed for sustainable management of plant diseases by using resistance inducers (e.g. endophytes). A basic point for the success of this approach is the study of endophytic bacteria associated with plants. In the present work, endophytic bacteria living in healthy and 'Ca. Phytoplasma mali'-infected apple trees were described by cultivation-dependent and independent methods. 16S rDNA sequence analysis showed the presence of the groups Proteobacteria, Acidobacteria, Bacteroidetes, Actinobacteria, Chlamydiae, and Firmicutes. In detail, library analyses underscored 24 and 17 operational taxonomic units (OTUs) in healthy and infected roots, respectively, with a dominance of Betaproteobacteria. Moreover, differences in OTUs number and in CFU/g suggested that phytoplasmas could modify the composition of endophytic bacterial communities associated with infected plants. Intriguingly, the combination of culturing methods and cloning analysis allowed the identification of endophytic bacteria (e.g. Bacillus, Pseudomonas, and Burkholderia) that have been reported as biocontrol agents. Future research will investigate the capability of these bacteria to control 'Ca. Phytoplasma mali' in order to develop sustainable approaches for managing AP.

Multilocus sequence analysis reveals the genetic diversity of European fruit tree phytoplasmas and supports the existence of inter-species recombination.

The genetic diversity of three temperate fruit tree phytoplasmas 'Candidatus Phytoplasma prunorum', 'Ca. P. mali' and 'Ca. P. pyri' has been established by multilocus sequence analysis. Among the four genetic loci used, the genes imp and aceF distinguished 30 and 24 genotypes, respectively, and showed the highest variability. Percentage of substitution for imp ranged from 50 to 68 % according to species. Percentage of substitution varied between 9 and 12 % for aceF, whereas it was between 5 and 6 % for pnp and secY. In the case of 'Ca P. prunorum' the three most prevalent aceF genotypes were detected in both plants and insect vectors, confirming that the prevalent isolates are propagated by insects. The four isolates known to be hypo-virulent had the same aceF sequence, indicating a possible monophyletic origin. Haplotype network reconstructed by eBURST revealed that among the 34 haplotypes of 'Ca. P. prunorum', the four hypo-virulent isolates also grouped together in the same clade. Genotyping of some Spanish and Azerbaijanese 'Ca. P. pyri' isolates showed that they shared some alleles with 'Ca. P. prunorum', supporting for the first time to our knowledge, the existence of inter-species recombination between these two species.

Comparison by Sequence-Based and Electron Microscopic Analyses of Fig mosaic virus Isolates Obtained from Field and Experimentally Inoculated Fig Plants.

Fig mosaic disease (FMD) and the fig mite, Aceria ficus, are widespread in different fig growing provinces of Turkey. Fig trees (Ficus carica) cv. Bursa siyahi (D1) and an unknown seedling (D2) that showed typical FMD symptoms and was heavily infested by fig mites were used as donor plants for attempted mite transmissions to healthy fig seedlings. Transmission electron microscopy observations of donor plant samples prior to the transmission tests were performed and showed the presence of double membrane bodies (DMBs) in the palisade mesophyll cells. Electron microscopy of all experimentally inoculated fig seedlings showed the same bodies. This result reinforced the suggestion that an agent that elicits the production of DMBs in infected cells is involved in the etiology of FMD. Double-stranded (ds)RNA analyses were also performed from experimentally inoculated plants, and dsRNAs with sizes approximately 1.30 and 1.96 kb were obtained. Reverse transcription-polymerase chain reaction (RT-PCR) products of 468 and 298 bp specific to Fig mosaic virus (FMV) were amplified from both donor and experimentally inoculated plants. BLAST analyses of nucleotide sequences of these fragments showed 90% identity with FMV for the donor plant and 94 to 96% for experimentally inoculated plants. According to these results, FMV is present in both donor and experimentally inoculated plants in Turkey, and this virus is transmissible by A. ficus from fig plant to fig plant.

Further characterization of a new recombinant group of Plum pox virus isolates, PPV-T, found in orchards in the Ankara province of Turkey.

Sixteen Plum pox virus (PPV) isolates collected in the Ankara region of Turkey were analyzed using available serological and molecular typing assays. Surprisingly, despite the fact that all isolates except one, which was a mix infection, were typed as belonging to the PPV-M strain in four independent molecular assays, nine of them (60%) reacted with both PPV-M specific and PPV-D specific monoclonal antibodies. Partial 5' and 3' genomic sequence analysis on four isolates demonstrated that irrespective of their reactivity towards the PPV-D specific monoclonal antibody, they were all closely related to a recombinant PPV isolate from Turkey, Ab-Tk. All three isolates for which the relevant genomic sequence was obtained showed the same recombination event as Ab-Tk in the HC-Pro gene, around position 1566 of the genome. Complete genomic sequencing of Ab-Tk did not provide evidence for additional recombination events in its evolutionary history. Taken together, these results indicate that a group of closely related PPV isolates characterized by a unique recombination in the HC-Pro gene is prevalent under field conditions in the Ankara region of Turkey. Similar to the situation with the PPV-Rec strain, we propose that these isolates represent a novel strain of PPV, for which the name PPV-T (Turkey) is proposed. Given that PPV-T isolates cannot be identified by currently available typing techniques, it is possible that their presence has been overlooked in other situations. Further efforts should allow a precise description of their prevalence and of their geographical distribution in Turkey and, possibly, in other countries.

Oligonucleotide microarray-based detection and genotyping of Plum pox virus.

Plum pox virus (PPV) is the most damaging viral pathogen of stone fruits. The detection and identification of its strains are therefore of critical importance to plant quarantine and certification programs. Existing methods to screen strains of PPV suffer from significant limitations such as the simultaneous detection and genotyping of several strains of PPV in samples infected with different isolates of the virus. A genomic strategy for PPV screening based on the viral nucleotide sequence was developed to enable the detection and genotyping of the virus from infected plant tissue or biological samples. The basis of this approach is a long 70-mer oligonucleotide DNA microarray capable of simultaneously detecting and genotyping PPV strains. Several 70-mer oligonucleotide probes were specific for the detection and genotyping of individual PPV isolates to their strains. Other probes were specific for the detection and identification of two or three PPV strains. One probe (universal), derived from the genome highly conserved 3' non-translated region, detected all individual strains of PPV. This universal PPV probe, combined with probes specific for each known strain, could be used for new PPV strain discovery. Finally, indirect fluorescent labeling of cDNA with cyanine after cDNA synthesis enhanced the sensitivity of the virus detection without the use of the PCR amplification step. The PPV microarray detected and identified efficiently the PPV strains in PPV-infected peach, apricot and Nicotiana benthamiana leaves. This PPV detection method is versatile, and enables the simultaneous detection of plant pathogens.