Diverse patterns of secondary structure across genes and transposable elements are associated with siRNA production and epigenetic fate

Abstract

RNA molecules carry information in their primary sequence and also their secondary structure. Secondary structure can confer important functional information, but it is also a potential signal for an RNAi-like host epigenetic response mediated by small interfering RNAs (siRNAs). In this study, we predicted local secondary structures in features of the maize genome, focusing on small regions that had folding energies similar to pre-miRNA loci. We found secondary structures to be common in retrotransposons, in Helitrons, and in genes. These structured regions mapped higher diversities of siRNAs than regions without structure, explaining up to 24% of variation of the siRNA distribution across some TE types. Among genes, those with secondary structure were 1.5-fold more highly expressed, on average, than genes without secondary structure. However, these genes were also more variably expressed across the 26 NAM lines, and this variability correlated with the number of mapping siRNAs. We conclude that local stem-loop structures are a nearly ubiquitous feature of expressed regions of the maize genome, that they correlate with higher siRNA mapping, and that they can represent a trade-off between functional need and the potentially negative consequences of siRNA production.