Abstract
Background: Long non-coding RNAs (lncRNAs) are involved in variable cleavage, transcriptional interference, regulation of DNA methylation and protein modification. However, the regulation of lncRNAs in plant somatic embryos remains unclear. The longan (Dimocarpus longan) somatic embryogenesis (SE) system is a good system for research on longan embryo development. Results: In this study, 7643 lncRNAs obtained during early SE in D. longan were identified by high-throughput sequencing, among which 6005 lncRNAs were expressed. Of the expressed lncRNAs, 4790 were found in all samples and 160 were specifically expressed in embryogenic callus (EC), 154 in incomplete embryogenic compact structures (ICpECs), and 376 in globular embryos (GEs). We annotated the 6005 expressed lncRNAs, and 1404 lncRNAs belonged to 506 noncoding RNA (ncRNA) families and 4682 lncRNAs were predicted to target protein-coding genes. The target genes included 5051 cis-regulated target genes (5712 pairs) and 1605 trans-regulated target genes (3618 pairs). KEGG analysis revealed that most of the differentially expressed target genes (mRNAs) of the lncRNAs were enriched in the "plant-pathogen interaction" and "plant hormone signaling" pathways during early longan SE. Real-time quantitative PCR confirmed that 20 selected lncRNAs showed significant differences in expression and that five lncRNAs were related to auxin response factors. Compared with the FPKM expression trends, 16 lncRNA expression trends were the same in qPCR. In lncRNA-miRNA-mRNA relationship prediction, 40 lncRNAs were predicted to function as eTMs for 15 miRNAs and 7 lncRNAs were identified as potential miRNA precursors. In addition, we verified the lncRNA-miRNA-mRNA regulatory relationships by transient expression of miRNAs (miR172a, miR159a.1 and miR398a). Conclusion: Analyses of lncRNAs during early longan SE showed that differentially expressed lncRNAs were involved in expression regulation at each SE stage, and may form a regulatory network with miRNAs and mRNAs. These findings provide new insights into lncRNAs and lay a foundation for future functional analysis of lncRNAs during early longan SE.