Deep-Sequencing of the Peach Latent Mosaic Viroid Reveals New Aspects of Population Heterogeneity

Viroids are small circular single-stranded infectious RNAs characterized by a relatively high mutation level. Knowledge of their sequence heterogeneity remains largely elusive and previous studies, using Sanger sequencing, were based on a limited number of sequences. In an attempt to address sequence heterogeneity from a population dynamics perspective, a GF305-indicator peach tree was infected with a single variant of the Avsunviroidae family member Peach latent mosaic viroid (PLMVd). Six months post-inoculation, full-length circular conformers of PLMVd were isolated and deep-sequenced. We devised an original approach to the bioinformatics refinement of our sequence libraries involving important phenotypic data, based on the systematic analysis of hammerhead self-cleavage activity. Two distinct libraries yielded a total of 3,939 different PLMVd variants. Sequence variants exhibiting up to ∼17% of mutations relative to the inoculated viroid were retrieved, clearly illustrating the high level of divergence dynamics within a unique population. While we initially assumed that most positions of the viroid sequence would mutate, we were surprised to discover that ∼50% of positions remained perfectly conserved, including several small stretches as well as a small motif reminiscent of a GNRA tetraloop which are the result of various selective pressures. Using a hierarchical clustering algorithm, the different variants harvested were subdivided into 7 clusters. We found that most sequences contained an average of 4.6 to 6.4 mutations compared to the variant used to initially inoculate the plant. Interestingly, it was possible to reconstitute and compare the sequence evolution of each of these clusters. In doing so, we identified several key mutations. This study provides a reliable pipeline for the treatment of viroid deep-sequencing. It also sheds new light on the extent of sequence variation that a viroid population can sustain, and which may give rise to a quasispecies.

[1]  C. Hammann,et al.  Idiosyncratic cleavage and ligation activity of individual hammerhead ribozymes and core sequence variants thereof , 2007, Biological chemistry.

[2]  E. Domingo,et al.  The quasispecies (extremely heterogeneous) nature of viral RNA genome populations: biological relevance--a review. , 1985, Gene.

[3]  Mapping studies of the Peach latent mosaic viroid reveal novel structural features , 2011 .

[4]  Tengke Xiong,et al.  A New MCA-Based Divisive Hierarchical Algorithm for Clustering Categorical Data , 2009, 2009 Ninth IEEE International Conference on Data Mining.

[5]  G. Steger,et al.  Analysis of thermal stress-mediated PSTVd variation and biolistic inoculation of progeny of viroid "thermomutants" to tomato and Brassica species. , 2004, Virology.

[6]  J. Motard,et al.  Molecular features of new Peach Latent Mosaic Viroid variants suggest that recombination may have contributed to the evolution of this infectious RNA. , 2007, Virology.

[7]  S. B. Needleman,et al.  A general method applicable to the search for similarities in the amino acid sequence of two proteins. , 1970, Journal of molecular biology.

[8]  Pavel Skums,et al.  Efficient error correction for next-generation sequencing of viral amplicons , 2012, BMC Bioinformatics.

[9]  J. Daròs,et al.  Arabidopsis thaliana has the enzymatic machinery for replicating representative viroid species of the family Pospiviroidae. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[10]  P. Lemey,et al.  Covering all bases in HIV research: unveiling a hidden world of viral evolution. , 2010, AIDS reviews.

[11]  C. Hernández,et al.  Genomic Structure of Three Phenotypically Different Isolates of Peach Latent Mosaic Viroid: Implications of the Existence of Constraints Limiting the Heterogeneity of Viroid Quasispecies , 1998, Journal of Virology.

[12]  Jiong Yang,et al.  CLUSEQ: efficient and effective sequence clustering , 2003, Proceedings 19th International Conference on Data Engineering (Cat. No.03CH37405).

[13]  C. Hernández,et al.  Rapid generation of genetic heterogeneity in progenies from individual cDNA clones of peach latent mosaic viroid in its natural host. , 1999, The Journal of general virology.

[14]  C. Masuta,et al.  A Viral Satellite RNA Induces Yellow Symptoms on Tobacco by Targeting a Gene Involved in Chlorophyll Biosynthesis using the RNA Silencing Machinery , 2011, PLoS pathogens.

[15]  M. Nei,et al.  Mathematical model for studying genetic variation in terms of restriction endonucleases. , 1979, Proceedings of the National Academy of Sciences of the United States of America.

[16]  Mário A. T. Figueiredo,et al.  Similarity-Based Clustering of Sequences Using Hidden Markov Models , 2003, MLDM.

[17]  Zasha Weinberg,et al.  Identification of Hammerhead Ribozymes in All Domains of Life Reveals Novel Structural Variations , 2011, PLoS Comput. Biol..

[18]  Ricard V Solé,et al.  The Fittest versus the Flattest: Experimental Confirmation of the Quasispecies Effect with Subviral Pathogens , 2006, PLoS pathogens.

[19]  Emese Meglécz,et al.  Accuracy and quality assessment of 454 GS-FLX Titanium pyrosequencing , 2011, BMC Genomics.

[20]  Rafael Sanjuán,et al.  Extremely High Mutation Rate of a Hammerhead Viroid , 2009, Science.

[21]  Martin Pelchat,et al.  The Subviral RNA Database: a toolbox for viroids, the hepatitis delta virus and satellite RNAs research , 2006, BMC Microbiology.

[22]  A. Gisel,et al.  Small RNAs containing the pathogenic determinant of a chloroplast-replicating viroid guide the degradation of a host mRNA as predicted by RNA silencing. , 2012, The Plant journal : for cell and molecular biology.

[23]  B. Ding.  The biology of viroid-host interactions. , 2009, Annual review of phytopathology.

[24]  M. Greenacre,et al.  Multiple Correspondence Analysis and Related Methods , 2006 .

[25]  P. Waterhouse,et al.  On the role of RNA silencing in the pathogenicity and evolution of viroids and viral satellites. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[26]  T. O. Diener.  Potato Spindle Tuber , 1987 .

[27]  S Ambrós,et al.  In vitro and in vivo self-cleavage of a viroid RNA with a mutation in the hammerhead catalytic pocket. , 1998, Nucleic acids research.

[28]  Xuehua Zhong,et al.  A Genomic Map of Viroid RNA Motifs Critical for Replication and Systemic Trafficking[W] , 2008, The Plant Cell Online.

[29]  R Cedergren,et al.  Numbering system for the hammerhead. , 1992, Nucleic acids research.

[30]  Joshua Zhexue Huang,et al.  A New Markov Model for Clustering Categorical Sequences , 2011, 2011 IEEE 11th International Conference on Data Mining.

[31]  J. Perreault,et al.  Characterization of a viroid-derived RNA promoter for the DNA-dependent RNA polymerase from Escherichia coli. , 2002, Biochemistry.

[32]  T. O. Diener.  Potato spindle tuber "virus". IV. A replicating, low molecular weight RNA. , 1971, Virology.

[33]  Jiong Yang,et al.  Towards automatic clustering of protein sequences , 2002, Proceedings. IEEE Computer Society Bioinformatics Conference.

[34]  J. Ouellet,et al.  Mapping in Solution Shows the Peach Latent Mosaic Viroid To Possess a New Pseudoknot in a Complex, Branched Secondary Structure , 2000, Journal of Virology.

[35]  Anton I. Petrov,et al.  A Three-Dimensional RNA Motif in Potato spindle tuber viroid Mediates Trafficking from Palisade Mesophyll to Spongy Mesophyll in Nicotiana benthamiana[W] , 2011, Plant Cell.