Nanopore RNA Sequencing Revealed Long Non-Coding and LTR Retrotransposon-Related RNAs Expressed at Early Stages of Triticale SEED Development

• Published in: Plants (Basel) Vol. 9; no. 12; p. 1794
• Main Authors: Kirov, Ilya, Dudnikov, Maxim, Merkulov, Pavel, Shingaliev, Andrey, Omarov, Murad, Kolganova, Elizaveta, Sigaeva, Alexandra, Karlov, Gennady, Soloviev, Alexander
• Format: Journal Article
• Language: English
• Published: Switzerland MDPI AG 17-12-2020; MDPI
• Subjects: Annotations; Cell division; Deoxyribonucleic acid; DNA; DNA structure; Flowers & plants; Gag protein; Gene sequencing; Genomes; Interspecific; Loci; long non-coding RNAs; Nanopore sequencing; Non-coding RNA; Nucleotide sequence; Open reading frames; Porosity; Proteins; retrotransposons; Ribonucleic acid; RNA; Rye; seed development; Seeds; Sequence analysis; Transposition; Transposons; Triticale; Virus-like particles; Wheat; United States > US; Russia; seed development; triticale; long non-coding RNAs; retrotransposons; Nanopore sequencing
• ISSN: 2223-7747; 2223-7747
• DOI: 10.3390/plants9121794

The intergenic space of plant genomes encodes many functionally important yet unexplored RNAs. The genomic loci encoding these RNAs are often considered "junk", DNA as they are frequently associated with repeat-rich regions of the genome. The latter makes the annotations of these loci and the assembly of the corresponding transcripts using short RNAseq reads particularly challenging. Here, using long-read Nanopore direct RNA sequencing, we aimed to identify these "junk" RNA molecules, including long non-coding RNAs (lncRNAs) and transposon-derived transcripts expressed during early stages (10 days post anthesis) of seed development of triticale (AABBRR, 2
= 6
= 42), an interspecific hybrid between wheat and rye. Altogether, we found 796 lncRNAs and 20 LTR retrotransposon-related transcripts (RTE-RNAs) expressed at this stage, with most of them being previously unannotated and located in the intergenic as well as intronic regions. Sequence analysis of the lncRNAs provide evidence for the frequent exonization of Class I (retrotransposons) and class II (DNA transposons) transposon sequences and suggest direct influence of "junk" DNA on the structure and origin of lncRNAs. We show that the expression patterns of lncRNAs and RTE-related transcripts have high stage specificity. In turn, almost half of the lncRNAs located in Genomes A and B have the highest expression levels at 10-30 days post anthesis in wheat. Detailed analysis of the protein-coding potential of the RTE-RNAs showed that 75% of them carry open reading frames (ORFs) for a diverse set of GAG proteins, the main component of virus-like particles of LTR retrotransposons. We further experimentally demonstrated that some RTE-RNAs originate from autonomous LTR retrotransposons with ongoing transposition activity during early stages of triticale seed development. Overall, our results provide a framework for further exploration of the newly discovered lncRNAs and RTE-RNAs in functional and genome-wide association studies in triticale and wheat. Our study also demonstrates that Nanopore direct RNA sequencing is an indispensable tool for the elucidation of lncRNA and retrotransposon transcripts.