Promoter-based identification of novel non-coding RNAs reveals the presence of dicistronic snoRNA-miRNA genes in Arabidopsis thaliana

Ge Qu, Katarzyna Kruszka, Patrycja Plewka, Shu-Yi Yang, Tzyy-Jen Chiou, Artur Jarmolowski, Zofia Szweykowska-Kulinska, Manuel Echeverria, Wojciech M. Karlowski

In the past few decades, non-coding RNAs (ncRNAs) have emerged as important regulators of gene expression in eukaryotes. Most studies of ncRNAs in plants have focused on the identification of silencing microRNAs (miRNAs) and small interfering RNAs (siRNAs). Another important family of ncRNAs that has been well characterized in plants is the small nucleolar RNAs (snoRNAs) and the related small Cajal body-specific RNAs (scaRNAs). Both target chemical modifications of ribosomal RNAs (rRNAs) and small nuclear RNAs (snRNAs). In plants, the snoRNA genes are organized in clusters, transcribed by RNA Pol II from a common promoter and subsequently processed into mature molecules. The promoter regions of snoRNA polycistronic genes in plants are highly enriched in two conserved cis-regulatory elements (CREs), Telo-box and Site II, which coordinate the expression of snoRNAs and ribosomal protein coding genes throughout the cell cycle. In order to identify novel ncRNA genes, we have used the snoRNA Telo-box/Site II motifs combination as a functional promoter indicator to screen the Arabidopsis genome. The predictions generated by this process were tested by detailed exploration of available RNA-Seq and expression data sets and experimental validation. As a result, we have identified several snoRNAs, scaRNAs and 'orphan' snoRNAs. We also show evidence for 16 novel ncRNAs that lack similarity to any reported RNA family. Finally, we have identified two dicistronic genes encoding precursors that are processed to mature snoRNA and miRNA molecules. We discuss the evolutionary consequences of this result in the context of a tight link between snoRNAs and miRNAs in eukaryotes. We present an alternative computational approach for non-coding RNA detection. Instead of depending on sequence or structure similarity in the whole genome screenings, we have explored the properties of promoter regions of well-characterized ncRNAs. Interestingly, besides expected ncRNAs predictions we were also able to recover single precursor arrangement for snoRNA-miRNA. Accompanied by analyses performed on rice sequences, we conclude that such arrangement might have interesting functional and evolutionary consequences and discuss this result in the context of a tight link between snoRNAs and miRNAs in eukaryotes.