Mutation rates vary by orders of magnitude across biological systems, being higher for simpler genomes. The simplest known genomes correspond to viroids, subviral plant replicons constituted by circular non-coding RNAs of few hundred bases. Previous work has revealed an extremely high mutation rate for chrysanthemum chlorotic mottle viroid, a chloroplast-replicating viroid. However, whether this is a general feature of viroids remains unclear. Here, we have used high-fidelity ultra-deep sequencing to determine the mutation rate in a common host (eggplant) of two viroids, each representative of one family: the chloroplastic eggplant latent viroid (ELVd, Avsunviroidae) and the nuclear potato spindle tuber viroid (PSTVd, Pospiviroidae). This revealed higher mutation frequencies in ELVd than in PSTVd, as well as marked differences in the types of mutations produced. Rates of spontaneous mutation, quantified in vivo using the lethal mutation method, ranged from 1/1000 to 1/800 for ELVd and from 1/7000 to 1/3800 for PSTVd depending on sequencing run. These results suggest that extremely high mutability is a common feature of chloroplastic viroids, whereas the mutation rates of PSTVd and potentially other nuclear viroids appear significantly lower and closer to those of some RNA viruses. Spontaneous mutations are the ultimate source of genetic variation and their characterization provides fundamental information about evolutionary processes. The highest mutation rate so far described corresponds to a hammerhead viroid infecting plant chloroplasts. Viroids are plant-exclusive parasites constituted by 250–400 nt-long, non-protein-coding RNAs, and are divided into two families with distinct mechanisms of replication and localization: chloroplastic (Avsunviroidae), and nuclear (Pospiviroidae). Here, we have used high-fidelity ultra-deep sequencing to compare side by side the mutation rates of one representative member of each viroid family in the same host. We found that the mutation rate of the nuclear viroid was several fold lower than that of the chloroplastic viroid.

With a minimal (250–400 nt), non-protein-coding, circular RNA genome, viroids rely on sequence/structural motifs for replication and colonization of their host plants. These motifs are embedded in a compact secondary structure whose elucidation is crucial to understand how they function. Viroid RNA structure has been tackled in silico with algorithms searching for the conformation of minimal free energy, and in vitro by probing in solution with RNases, dimethyl sulphate and bisulphite, and with selective 2′-hydroxyl acylation analyzed by primer extension (SHAPE), which interrogates the RNA backbone at single-nucleotide resolution. However, in vivo approaches at that resolution have not been assayed. Here, after confirming by 3 termodynamics-based predictions and by in vitro SHAPE that the secondary structure adopted by the infectious monomeric circular (+) RNA of potato spindle tuber viroid (PSTVd) is a rod-like conformation with double-stranded segments flanked by loops, we have probed it in vivo with a SHAPE modification. We provide direct evidence that a similar, but not identical, rod-like conformation exists in PSTVd-infected leaves of Nicotiana benthamiana, verifying the long-standing view that this RNA accumulates in planta as a “naked” form rather than tightly associated with protecting host proteins. However, certain nucleotides of the central conserved region, including some of the loop E involved in key functions such as replication, are more SHAPE-reactive in vitro than in vivo. This difference is most likely due to interactions with proteins mediating some of these functions, or to structural changes promoted by other factors of the in vivo habitat.

Eggplant latent viroid (ELVd), like other members of family Avsunviroidae, replicates in plastids through a symmetric rolling-circle mechanism in which elongation of RNA strands is most likely catalyzed by a nuclear-encoded polymerase (NEP) translocated to plastids. Here we have addressed where NEP initiates transcription of viroid strands. Because this step is presumably directed by sequence/structural motifs, we have previously determined the conformation of the monomeric linear (+) and (−) RNAs of ELVd resulting from hammerhead-mediated self-cleavage. In silico predictions with 3 softwares led to similar bifurcated conformations for both ELVd strands. In vitro examination by non-denaturing PAGE showed that they migrate as prominent single bands, with the ELVd (+) RNA displaying a more compact conformation as revealed by its faster electrophoretic mobility. In vitro SHAPE analysis corroborated the ELVd conformations derived from thermodynamics-based predictions in silico. Moreover, sequence analysis of 94 full-length natural ELVd variants disclosed co-variations, and mutations converting canonical into wobble pairs or vice versa, which confirmed in vivo most of the stems predicted in silico and in vitro, and additionally helped to introduce minor structural refinements. Therefore, results from the 3 experimental approaches were essentially consistent among themselves. Application to RNA preparations from ELVd-infected tissue of RNA ligase-mediated rapid amplification of cDNA ends, combined with pretreatments to modify the 5′ ends of viroid strands, mapped the transcription initiation sites of ELVd (+) and (−) strands in vivo at different sequence/structural motifs, in contrast with the situation previously observed in 2 other members of the family Avsunviroidae.

Spatial ITH is defined by genomic and biological variations within a tumour acquired by tumour cell evolution under diverse microenvironments, and its role in NB patient prognosis is understudied. In this work, we applied pangenomic techniques to detect chromosomal aberrations in at least two different areas of each tumour and/or in simultaneously obtained solid and liquid biopsies, detecting ITH in the genomic profile of almost 40% of HR-NB. ITH was better detected when comparing one or more tumour pieces and liquid biopsy (50%) than between different tumour pieces (21%). Interestingly, we found that patients with ITH analysed by pangenomic techniques had a significantly better survival rate that those with non-heterogeneous tumours, especially in cases without MYCN amplification. Moreover, all patients in the studied cohort with high ITH (defined as 50% or more genomic aberration differences between areas of a tumour or simultaneously obtained samples) survived after 48 months. These results clearly support analysing at least two solid tumour areas (separately or mixed) and liquid samples to provide more accurate genomic diagnosis, prognosis and therapy options in HR-NB.

Viroids are small infectious, non-protein-coding circular RNAs that replicate independently and, in some cases, incite diseases in plants. They are classified into two families: Pospiviroidae, composed of species that have a central conserved region (CCR) and replicate in the cell nucleus, and Avsunviroidae, containing species that lack a CCR and whose multimeric replicative intermediates of either polarity generated in plastids self-cleave through hammerhead ribozymes. The compact, rod-like or branched, secondary structures of viroid RNAs have been predicted by RNA folding algorithms and further examined using different in vitro and in vivo experimental techniques. However, direct data about their native tertiary structure remain scarce. Here we have applied atomic force microscopy (AFM) to image at single-molecule resolution different variant RNAs of three representative viroids: potato spindle tuber viroid (PSTVd, family Pospiviroidae), peach latent mosaic viroid and eggplant latent viroid (PLMVd and ELVd, family Avsunviroidae). Our results provide a direct visualization of their native, three-dimensional conformations at 0 and 4 mM Mg²⁺ and highlight the role that some elements of tertiary structure play in their stabilization. The AFM images show that addition of 4 mM Mg²⁺ to the folding buffer results in a size contraction in PSTVd and ELVd, as well as in PLMVd when the kissing-loop interaction that stabilizes its 3D structure is preserved.