Abstract
The 5'-proximal region of the citrus tristeza virus (CTV) RNA genome is a hub where several elements involved in different facets of the virus cycle reside, including the sequences driving the production of the viral long non-coding RNA (lncRNA) LMT1. The sequence of this region is one of the most divergent genome areas, allowing for strain differentiation. Beyond its use in assessing viral population diversity, the region provides a valuable model for studying the conservation of RNA structure and function despite sequence variation. Here, we integrated comparative in silico analysis of the LMT1 region from variants of eight CTV strains with selective 2'-hydroxyl acylation, analyzed by primer extension and mutational profiling (SHAPE-MaP) probing of in vitro-generated LMT1 RNAs from two divergent strains, T36 and T68. The predicted consensus structures revealed 19 putative, conserved stem-loops. The SHAPE-MaP reactivity data supported and substantiated the thermodynamics-based predictions for the 15 previously uncharacterized stem-loops and two functional elements identified earlier. The strong structural conservation across strains highlights that the LMT1 RNA structure contributes to its function during CTV infection. These results provide the first experimentally supported structure of this viral lncRNA and lay the foundation for defining how individual RNA motifs influence CTV biology.