Background: Non-coding RNAs (ncRNAs), including microRNAs (miRNAs), long ncRNAs (lncRNAs) and circular RNAs (circRNAs), accomplish remarkable variety of biological functions. However, the composition of ncRNAs and their interactions with coding RNAs in modulating and controlling of cellular process in plants is largely unknown. Using a diverse group of high-throughput sequencing strategies, the mRNA, miRNA, lncRNA and circRNA compositions of tobacco (Nicotiana tabacum) roots determined and their alteration and potential biological functions in response to topping treatment analyzed.Results: A total of 688 miRNAs, 7423 non-redundant lncRNAs and 12,414 circRNAs were identified, among which, some selected differentially expressed RNAs were verified by quantitative real-time PCR. Using the differentially expressed RNAs, a co-expression network was established that included all four types of RNAs. The number of circRNAs identified were higher than that of miRNAs and lncRNAs, but only two circRNAs were present in the co-expression network. LncRNAs appear to be the most active ncRNAs based on their numbers presented in the co-expression network, but none of them seems to be an eTM (endogenous Target Mimicry) of miRNAs. Integrated with analyses of sequence interaction, several mRNA-circRNA-miRNA interaction networks with a potential role in the regulation of nicotine biosynthesis were uncovered, including a QS-circQS-miR6024 interaction network. In this network miR6024 was significantly down-regulated, while the expression levels of its two targets, circQS and its host geneQS, were sharply increased following the topping treatment.Conclusions: These results illustrated the transcriptomic profiles of tobacco roots, the organ responsible for nicotine biosynthesis. mRNAs always play the most important roles, while ncRNAs are also expressed extensively for topping treatment response, especially circRNAs are the most activated in the ncRNA pool. These studies also provided insights on the coordinated regulation module of coding and non-coding RNAs in a single plant biological sample. The findings reported here indicate that ncRNAs appear to form interaction complex for the regulation of stress response forming regulation networks with transcripts involved in nicotine biosynthesis in tobacco.

Extracellular microRNAs (miRNAs) have been proposed to function in cross-kingdom gene regulation. Among these, plant-derived miRNAs of dietary origin have been reported to survive the harsh conditions of the human digestive system, enter the circulatory system, and regulate gene expression and metabolic function. However, definitive evidence supporting the presence of plant-derived miRNAs of dietary origin in mammals has been difficult to obtain due to limited sample sizes. We have developed a bioinformatics pipeline (ePmiRNA_finder) that provides strident miRNA classification and applied it to analyze 421 small RNA sequencing data sets from 10 types of human body fluids and tissues and comparative samples from carnivores and herbivores. A total of 35 miRNAs were identified that map to plants typically found in the human diet and these miRNAs were found in at least one human blood sample and their abundance was significantly different when compared to samples from human microbiome or cow. The plant-derived miRNA profiles were body fluid/tissue-specific and highly abundant in the brain and the breast milk samples, indicating selective absorption and/or the ability to be transported across tissue/organ barriers. Our data provide conclusive evidence for the presence of plant-derived miRNAs as a consequence of dietary intake and their cross-kingdom regulatory function within human circulating system.