Comparative analysis of flavivirus sfRNA dynamics and secondary structure

The accumulation of subgenomic flavivirus RNAs (sfRNAs) modulates viral fitness and pathogenicity in culture and in vivo. These noncoding RNAs are produced by incomplete digestion of the flavivirus genome by the cellular 5'-3' exoribonuclease (XRN1). Diverse flaviviruses have conserved RNA structural elements (RSEs) that map to their 3'-untranslated region (3'-UTR): Xrn-resistant RNA structures, dumbbell structures, and a 3'-stem loop (3'SL). Despite the importance of the 3'-UTR RSEs for flavivirus replication, the structural dynamics of sfRNA during flavivirus infection are understudied. Here, we use digital droplet PCR to quantify sfRNA levels during infection for a panel of mosquito-borne flaviviruses (MbFV) including dengue virus serotypes 1 (DENV1), 2 (DENV2), and 4 (DENV4), and Zika virus (ZIKV). We then used SHAPE-MaP on XRN1-digested, in vitro-transcribed sfRNAs from each virus to determine their secondary structures compared to the corresponding sfRNAs obtained from flavivirus-infected A549 cells. Results seen in-cell and in vitro were largely similar; however, motifs within the dumbbell, the small hairpin (sHP) directly upstream of the 3'-SL, and 3'-SL regions showed significant differences in the extent of nucleotide reactivity. These differences were consistent among the four flaviviruses examined and may indicate regions of sfRNA that are shielded by interaction with proteins or other nucleic acids during infection. However, strong protection indicative of sustained interaction was not apparent. Our findings suggest that sfRNA interactions with viral and host factors within the cell are few, occur via base-paired regions, or are highly transient.